A genome-wide association and fine-mapping study of white rust resistance in hexaploid chrysanthemum cultivars with a wild diploid reference genome

Abstract

White rust caused by Puccinia horiana is one of the most serious diseases of chrysanthemum (Chrysanthemum × morifolium). In this study, we report the DNA markers associated with resistance against P. horiana via a simple approach using the genome of a wild diploid relative, Chrysanthemum seticuspe. First, we identified the important region of the genome in the resistant cultivar “Ariesu” via a genome-wide association study. Simplex single nucleotide polymorphism (SNP) markers mined from ddRAD-Seq were used in a biparental population originating from crosses between resistant “Ariesu” and susceptible “Yellow Queen”. The C. seticuspe genome was used as a reference. For the fine mapping of P. horiana resistance locus 2 (Phr2), a comparative whole genome sequencing study was conducted. Although the genome sequences of chrysanthemum cultivars assembled via the short-read approach were fragmented, reliable genome alignments were reconstructed by mapping onto the chromosome level of the C. seticuspe pseudomolecule. Base variants were then identified by comparing the assembled genome sequences of resistant “Ariesu” and susceptible “Yellow Queen”. Consequently, SNP markers that were closer to Phr2 compared with ddRAD-Seq markers were obtained. These SNP markers co-segregated with resistance in F₁ progenies originating from resistant “Ariesu” and showed robust transferability for detecting Phr2-conferring resistance among chrysanthemum genetic resources. The wild C. seticuspe pseudomolecule, a de facto monoploid genome used for ddRAD-Seq analysis and assembled genome sequence comparison, demonstrated this method’s utility as a model for developing DNA markers in hexaploid chrysanthemum cultivars.

References

[1] Alaei H, De Backer M, Nuytinck J et al. Phylogenetic relationships of Puccinia horiana and other rust pathogens of Chrysanthemum × morifolium based on rDNA ITS sequence analysis . Mycol Res. 2009; 113: 668-83.

[2] Ashraf BH, Jensen J, Asp T et al. Association studies using family pools of outcrossing crops based on allele-frequency estimates from DNA sequencing. Theor Appl Genet. 2014; 127: 1331-41.

[3] Baker JJ . Chrysanthemum white rust in England and Wales 1963-66. Plant Pathol. 1967; 16: 162-6.

[4] Benjamini Y, Hochberg Y . Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995; 57: 289-300.

[5] Bradbury PJ, Zhang Z, Kroon DE et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633-5.

[6] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[7] De Backer M, Alaei H, van Bockstaele E et al. Identification and characterization of pathotypes in Puccinia horiana, a rust pathogen of Chrysanthemum x morifolium . Eur J Plant Pathol. 2011; 130: 325-38.

[8] De Jong J, Rademaker W . The reaction of Chrysanthemum cultivars to Puccinia horiana and the inheritance of resistance . Euphytica. 1986; 35: 945-52.

[9] Dickens JSW . The resistance of various cultivars and species of chrysanthemum to white rust (Puccinia horiana Henn.) . Plant Pathol. 1968; 17: 19-22.

[10] Edger PP, Poorten TJ, Van Buren R et al. Origin and evolution of the octoploid strawberry genome. Nat Genet. 2019; 51: 541-7.

[11] Eitas TK, Dangl JL . NB-LRR proteins: pairs, pieces, perception, partners, and pathways. Curr Opin Plant Biol. 2010; 13: 472-7.

[12] Ellis JG, Lagudah ES, Spielmeyer W et al. The past, present and future of breeding rust resistant wheat. Front Plant Sci. 2014; 5: 641.

[13] Elshire RJ, Glaubitz JC, Sun Q et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011; 6: e19379.

[14] Flor HH . The complementary genic systems in flax and flax rust. Adv Genet. 1956; 8: 29-54.

[15] Flor HH . Current status of the gene-for-gene concept. Annu Rev Phytopathol. 1971; 9: 275-96.

[16] van Geest G, Bourke PM, Voorrips RE et al. An ultra-dense integrated linkage map for hexaploid chrysanthemum enables multi-allelic QTL analysis. Theor Appl Genet. 2017a; 130: 2527-41.

[17] van Geest G, Voorrips RE, Esselink D et al. Conclusive evidence for hexasomic inheritance in chrysanthemum based on analysis of a 183 k SNP array. BMC Genomics. 2017b; 18: 585.

[18] Hagan WL, Hooker AL . Genetics of reaction to Puccinia sorghi in eleven corn inbred lines from Central and South America. Phytopathology. 1965; 55: 193-7.

[19] Herrera-Foessel SA, Huerta-Espino J, Calvo-Salazar V et al. Lr72 confers resistance to leaf rust in durum wheat cultivar Atil C2000 . Plant Dis. 2014; 98: 631-5.

[20] Hirakawa H, Sumitomo K, Hisamatsu T et al. De novo whole-genome assembly in Chrysanthemum seticuspe, a model species of Chrysanthemums, and its application to genetic and gene discovery analysis . DNA Res. 2019; 26: 195-203.

[21] Hiratsuka N. Three species of chrysanthemum rust in Japan and its neighboring districts. Sydowia. 1957; 2: 34-44.

[22] Hooker AL, Russell WA . Inheritance of resistance to Puccinia sorghi in six corn inbred lines. Phytopathology. 1962; 52: 122-8.

[23] Iwai T, Nakamura S, Itabashi T et al. Resistance to Puccinia horiana in chrysanthemum cultivars (in Japanese) . Plant Prot. 2009; 63: 26-30.

[24] Kavanagh PJ, Singh D, Bansal UK et al. Inheritance and characterization of the new and rare gene Rph25 conferring seedling resistance in Hordeum vulgare against Puccinia hordei . Plant Breed. 2017; 136: 908-12.

[25] Koboldt DC, Zhang Q, Larson DE et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012; 22: 568-76.

[26] Langmead B, Salzberg SL . Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9: 357-9.

[27] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34: 3094-100.

[28] Li H, Handsaker B, Wysoker A et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[29] van Lieshout N, van Kaauwen M, Kodde L et al. De novo whole-genome assembly of Chrysanthemum makinoi, a key wild chrysanthemum . G3 Genes|Genomes|Genetics. 2022; 23: jkab358.

[30] Luo R, Liu B, Xie Y et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012; 1: 18.

[31] Martin PH, Firman ID . Resistance of chrysanthemum cultivars to white rust (Puccinia horiana) . Plant Pathol. 1970; 19: 180-4.

[32] Miler N, Kulus D, Sliwinska E . Nuclear DNA content as an indicator of inflorescence colour stability of in vitro propagated solid and chimera mutants of chrysanthemum. Plant Cell Tissue Organ Cult. 2020; 143: 421-30.

[33] Nakano M, Hirakawa H, Fukai E et al. A chromosome-level genome sequence of chrysanthemum seticuspe, a model species for hexaploid cultivated Chrysanthemum. Communications Biology. 2021; 4: 1167.

[34] Nakano M, Taniguchi K, Masuda Y et al. A pure line derived from a self-compatible Chrysanthemum seticuspe mutant as a model strain in the genus Chrysanthemum . Plant Sci. 2019; 287: 110174.

[35] O’Keefe G, Davis DD . Morphology of Puccinia horiana, causal agent of chrysanthemum white rust, sampled from naturally infected plants. Plant Dis. 2015; 99: 1738-43.

[36] Park RF, Golegaonkar PG, Derevnina L et al. Leaf rust of cultivated barley: pathology and control. Annu Rev Phytopathol. 2015; 53: 565-89.

[37] Park RF, Mohler V, Nazari K et al. Characterisation and mapping of gene Lr73 conferring seedling resistance to Puccinia triticina in common wheat . Theor Appl Genet. 2014a; 127: 2041-9.

[38] Park SK, Lim JH, Shin HK et al. Identification of chrysanthemum genetic resources resistant to white rust caused by Puccinia horiana . Plant Breeding and Biotechnology. 2014b; 2: 184-93.

[39] Qureshi N, Bariana H, Kumran V et al. A new leaf rust resistance gene Lr79 mapped in chromosome 3BL from the durum wheat landrace Aus26582. Theor Appl Genet. 2018; 131: 1091-8.

[40] Robinson JT, Thorvaldsdóttir H, Winckler W et al. Integrative genomics viewer. Nat Biotechnol. 2011; 29: 24-6.

[41] Saxena KMS, Hooker AL . On the structure of a gene for disease resistance in maize. Proc Natl Acad Sci U S A. 1968; 61: 1300-5.

[42] Shirasawa K, Tanaka M, Takahata Y et al. A high-density SNP genetic map consisting of a complete set of homologous groups in autohexaploid sweetpotato (Ipomoea batatas) . Sci Rep. 2017; 7: 44207.

[43] Singh D, Mohler V, Park R . Discovery, characterisation and mapping of wheat leaf rust resistance gene Lr71 . Euphytica. 2013; 190: 131-6.

[44] Song C, Liu Y, Song A et al. The Chrysanthemum nankingense genome provides insights into the evolution and diversification of chrysanthemum flowers and medicinal traits. Mol Plant. 2018; 11: 1482-91.

[45] Sumitomo K, Shirasawa K, Isobe S et al. Genome-wide association study overcomes the genome complexity in autohexaploid chrysanthemum and tags SNP markers onto the flower color genes. Sci Rep. 2019; 9: 13947.

[46] Sumitomo K, Shirasawa K, Isobe S et al. DNA marker for resistance to Puccinia horiana in chrysanthemum (Chrysanthemum morifolium Ramat.) “Southern Pegasus”. Breed Sci. 2021; 71: 261-7.

[47] Van der Vossen E, Sikkema A, Hekkert B 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. 2003; 36: 867-82.

[48] Wang J, Roe B, Macmil S et al. Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes. BMC Genomics. 2010; 11: 261.

[49] Wen X, Li J, Wang L et al. The Chrysanthemum lavandulifolium genome and the molecular mechanism underlying diverse capitulum types. Horticulture Research. 2022; 9: unab022.

[50] Wilkinson DR, Hooker AL . Genetics of reaction to Puccinia sorghi in ten corn inbred lines from Africa and Europe. Phytopathology. 1968; 58: 605-8.

[51] Yamaguchi T . Chrysanthemum breeding for resistance to white rust. Japanese Journal of Breeding. 1981; 31: 121-32.

[52] Yang L, Wang D, Xu Y et al. A new resistance gene against potato late blight originating from Solanum pinnatisectum located on potato chromosome 7. Front Plant Sci. 2017; 8: 1729.

[53] Yu X, Kong HY, Meiyalaghan V et al. Genetic mapping of a barley leaf rust resistance gene Rph26 introgressed from Hordeum bulbosum . Theor Appl Genet. 2018; 131: 2567-80.

[54] Ziems LA, Hickey LT, Platz GJ et al. Characterization of Rph24: a gene conferring adult plant resistance to Puccinia hordei in barley . Phytopathology. 2017; 107: 834-41.