Development of novel InDel markers by whole-genome sequence comparison and genetic diversity assessment of Thailand rice blast fungus populations

Abstract

InDel markers are commonly used to assess genetic relationships among populations. In this study, we employed a whole-genome sequence comparison method to identify and develop InDel markers for the rice blast fungus Pyricularia oryzae. We analyzed 152 whole-genome sequences of P. oryzae isolates from diverse global regions, including Brazil, Burundi, China, Colombia, Côte d'Ivoire, France, Ghana, Hungary, India, Japan, Korea, Laos, Madagascar, Mali, Morocco, Nepal, the Philippines, Portugal, Spain, Suriname, Thailand, the UK, the USA, and Zambia. Our analysis identified a total of 233,595 InDel loci distributed across the seven chromosomes of P. oryzae. From these, 82 loci were selected based on their high polymorphism across the 152 genome sequences. The effectiveness of these 82 loci was assessed by analyzing the genetic diversity of 47 Thai rice blast isolates alongside two reference isolates, GUY11 (France) and KJ201 (Korea). Of the 82 InDel loci, 33 exhibited polymorphisms, with 2–4 alleles per locus and polymorphic information content (PIC) scores ranging from 0.04 to 0.67. Principal coordinate and structure analyses revealed two genetic subgroups among the Thai rice blast isolates, categorized according to host specificity. Genetic relationships highlighted disparities among rice blast populations based on their respective hosts: rice and grassy weeds. This finding suggests a correlation between genetic relatedness and the plant hosts susceptible to rice blast disease. The newly developed InDel markers provide a valuable resource for future research in this field.

Keywords

Genetic diversity / Rice blast fungus / Rice blast / Polymorphism / Biological Sciences / Genetics

References

[1] AdedzeYMN, LuX, XiaY, et al. . Agarose-resolvable InDel markers based on whole genome re-sequencing in cucumber. Sci Rep, 2021, 11: 3872.

[2] AgarwalM, ShrivastavaN, PadhH. Advances in molecular marker techniques and their applications in plant sciences. Plant Cell Rep, 2008, 274617-631.

[3] AntunesM, Sá-CorreiaI. The NPR/Hal family of protein kinases in yeasts: biological role, phylogeny and regulation under environmental challenges. CSBJ, 2022, 20: 5698-5712.

[4] BaudinM, Le Naour-VernetM, GladieuxP, et al. . Pyriculariaoryzae: Lab star and field scourge. Mol Plant Pathol, 2024, 254e13449.

[5] BhattaraiG, MehlenbacherSA. In silico development and characterization of tri-nucleotide simple sequence repeat markers in hazelnut (Corylus avellana L.). PLoS One, 2017, 125e0178061.

[6] BolgerAM, LohseM, UsadelB. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30152114-2120.

[7] ChungH, GohJ, HanSS, RohJH, KimY, HeuS, ShimHK, JeongDG, KangIJ, YangJW. Comparative pathogenicity and host ranges of Magnaportheoryzae and related species. Plant Pathol J, 2020, 373195-205.

[8] CouchBC, FudalI, LebrunMH, TharreauD, ValentB, van KimP, NottéghemJL, KohnLM. Origins of host-specific populations of the blast pathogen Magnaportheoryzae in crop domestication with subsequent expansion of pandemic clones on rice and weeds of rice. Genetics, 2005, 1702613-630.

[9] DanecekP, BonfieldJK, LiddleJ, MarshallJ, et al. . Twelve years of SAMtools and BCFtools. Gigascience, 2021, 102giab008.

[10] DoyleJJ, DoyleJL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull, 1987, 19: 11-15

[11] Gladieux P, Condon B, Ravel S, Soanes D, Maciel JLN, Nhani A Jr, Chen L, Terauchi R, Lebrun MH, Tharreau D, Mitchell T, Pedley KF, Valent B, Talbot NJ, Farman M, Fournier E (2018a) Gene flow between divergent cereal- and grass-specific lineages of the rice blast fungus Magnaportheoryzae. mBio 9(1):e01219-17. https://doi.org/10.1128/mbio.01219-17

[12] GladieuxP, RavelS, RieuxA, Cros-ArteilS, et al. . Coexistence of multiple endemic and pandemic lineages of the rice blast pathogen. mBio, 2018, 92e01806-17.

[13] Gupta P, Roy J, Prasad M (2001) Single nucleotide polymorphisms: a new paradigm for molecular marker technology and DNA polymorphism detection with emphasis on their use in plants. Curr Sci 80(4):524–535

[14] JainA, RoorkiwalM, KaleS, GargV, YadalaR, VarshneyRK. InDel markers: an extended marker resource for molecular breeding in chickpea. PLoS One, 2019, 143e0213999.

[15] KhanMAI, AliMA, MonsurMA, et al. . Diversity and distribution of rice blast (Pyriculariaoryzae Cavara) races in Bangladesh. Plant Dis, 2016, 100: 2025-2033.

[16] KhodaeiaminjanM, KafkasS, MotalebipourEZ, CobanN. Insilico polymorphic novel SSR marker development and the first SSR-based genetic linkage map in pistachio. Tree Genet Genomes, 2018, 14: 45.

[17] KorinsakS, TangphatsornruangS, PootakhamW. Genome-wide association mapping of virulence gene in rice blast fungus Magnaportheoryzae using a genotyping by sequencing approach. Genomics, 2019, 1114661-668.

[18] LabbéJ, ZhangX, YinT, SchmutzJ, GrimwoodJ, MartinF, TuskanGA, Le TaconF. A genetic linkage map for the ectomycorrhizal fungus Laccariabicolor and its alignment to the whole-genome sequence assemblies. New Phytol, 2008, 1802316-328.

[19] LevyM, Correa-VictoriaFJ, ZeiglerRS, XuS, HamerJE. Genetic diversity of the rice blast fungus in a disease nursery in Colombia. Phytopathol, 1993, 83121427-1433.

[20] LiH, DurbinR. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics, 2010, 265589-595.

[21] LiYL, LiuJX. StructureSelector: a web based software to select and visualize the optimal number of clusters using multiple methods. Mol Ecol Resour, 2018, 18: 176-177.

[22] LiY, LuoX, PengX, JinY, TanH, WuL, LiJ, PeiY, XuX, ZhangW. Development of InDel markers and identify candidate flowering time genes in radish (Raphanus sativus L.). BMC Genomics, 2023, 24: 445.

[23] LvHH, YangLM, KangJG, WangQB, et al. . Development of InDel markers linked to Fusarium wilt resistance in cabbage. Mol Breeding, 2013, 32: 961-967.

[24] MacielJLN. Magnaportheoryzae, the blast pathogen: current status and options for its control. Plant Sci Rev, 2011, 264: 233-240.

[25] MacLeodA, TaitA, TurnerCM. The population genetics of Trypanosoma brucei and the origin of human infectivity. Philos Trans R Soc Lond B Biol Sci, 2001, 35614111035-1044.

[26] MiahG, RafiiMY, IsmailMR, PutehAB, RahimHA, IslamK, LatifMA. A review of microsatellite markers and their applications in rice breeding programs to improve blast disease resistance. Int J Mol Sci, 2013, 141122499-22528.

[27] NgernmuenA, SuktrakulW, DamchuayK, LongyaA, Kate-NgamS, JantasuriyaratC. Substantial enhancement of high polymorphic SSR marker development using in silico method from 18 available rice blast fungus genome sequences and its application in genetic diversity assessment. Biol, 2019, 74: 1181-1189.

[28] NiihamaM, MochizukiM, KurataN, NonomuraK. PCR-based INDEL markers co-dominant between Oryzasativa, japonica cultivars and closely-related wild Oryza species. Breed Sci, 2015, 654357-361.

[29] OliveiraM, AzevedoL. Molecular markers: an overview of data published for fungi over the last ten years. J Fungi (Basel), 2022, 88803.

[30] OsorioCE, UdallJA, SalvoGH, MaureiraBIJ. Development and characterization of InDel markers for Lupinusluteus L. (Fabaceae) and cross-species amplification in other Lupin species. Electron J Biotechnol, 2018, 31: 44-47.

[31] PaolaCF, EdwinRP, VictorSS, BruceWH, HankTS, MarshallCL, et al. . Study of the genetic diversity of the aflatoxin biosynthesis cluster in Aspergillus section Flavi using insertion/deletion markers in peanut seeds from Georgia, USA. Mycol, 2017, 1092200-209.

[32] PeakallR, SmousePE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes, 2006, 6: 288-295.

[33] PeakallR, SmousePE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics, 2012, 28: 2537-2539.

[34] QiH, YangJ, YinC, ZhaoJ, RenX, JiaS, ZhangG. Analysis of Pyricularia oryzae and P. grisea from different hosts based on multilocus phylogeny and pathogenicity associated with host preference in China. Phytopathol, 2019, 10981433-1440.

[35] SerroteCML, ReinigerLRS, SilvaKB, RabaiolliSMDS, StefanelCM. Determining the polymorphism information content of a molecular marker. Gene, 2020, 5726144175.

[36] SirisathawornT, SriratT, LongyaA, JantasuriyaratC. Evaluation of mating type distribution and genetic diversity of three Magnaportheoryzae avirulence genes, PWL-2, AVR-Pii and Avr-Piz-t, in Thailand rice blast isolates. Agr Nat Resour, 2017, 5117-14.

[37] SirithunyaP, SreewongchaiT, SriprakhonS, ToojindaT, PimpisithavornS, KosawangC. Assessment of genetic diversity in Thai isolates of Pyricularia grisea by random amplification of polymorphic DNA. J Phytopathol, 2008, 1564196-204.

[38] SorkhehK, PrudencioAS, GhebinejadA, DehkordiMK, ErogulD, RubioM, Martínez-GómezP. Insilico search, characterization and validation of new EST-SSR markers in the genus Prunus. BMC Res Notes, 2016, 9: 336.

[39] Suzuki R, Terada Y, Shimodaira H (2017) Pvclust: hierarchical clustering with p values via multiscale bootstrap resampling. R package version 1.3–2. http://CRAN.R-project.org/package=pvclust. Accessed 15 Jan. 2023

[40] TalbotNJ. On the trail of a cereal killer: exploring the biology of Magnaporthegrisea. Annu Rev Microbiol, 2003, 57: 177-202.

[41] ThierryM, GladieuxP, FournierE, TharreauD, IoosR. A genomic approach to develop a new qPCR test enabling detection of the Pyriculariaoryzae lineage causing wheat blast. Plant Dis, 2020, 104160-70.

[42] TsukiboshiT. Blast disease of ryegrass in Japan: The characteristics, geographical distribution, and control by fungicide-coated seeds. JARQ-Jpn Agric Res Q, 2022, 564313-320.

[43] TuanT, DuongVX, KhanhTD. Evaluation of genetic diversity of rice blast fungus ( Magnaporthe oryzae Barr) isolates collected from South central coast areas of Vietnam. Development, 2020, 2: 3

[44] Xu P, Lu C, Sun Z, Kuang Y, Cao D, Huo T, Li C, Jin H, Zheng X (2022) InSilico screening and development of microsatellite markers for genetic analysis in Percafluviatilis. Animals (Basel) 12(14):1809. https://doi.org/10.3390/ani12141809

[45] YadavMK, AravindanS, RaghuS, PrabhukarthikeyanSR, KeerthanaU, NgangkhamU, et al. . Assessment of genetic diversity and population structure of Magnaporthe oryzae causing rice blast disease using SSR markers. Physiol Mol Plant Pathol, 2019, 106: 157-165.

[46] YangJ, WangY, ShenH, YangW. Insilico identification and experimental validation of insertion–deletion polymorphisms in tomato genome. DNA Res, 2014, 214429-438.

[47] YangJ, HeJ, WangD, ShiE, YangW, GengQ, WangZ. Progress in research and application of InDel markers. Biodiv Sci., 2016, 242237.

[48] YounasMU, WangG, DuH, ZhangY, AhmadI, RajputN, et al. . Approaches to reduce rice blast disease using knowledge from host resistance and pathogen pathogenicity. Int J Mol Sci, 2023, 2454985.

[49] ZhongZ, ChenM, LinL, HanY, BaoJ, et al. . Population genomic analysis of the rice blast fungus reveals specific events associated with expansion of three main clades. ISME J, 2018, 1281867-1878.

RIGHTS & PERMISSIONS

The Author(s)