Linkage QTL Mapping and Genome-Wide Association Study on Resistance in Chickpea to Pythium ultimum

Abstract

The soilborne oomycete plant pathogen Pythium ultimum causes seed rot and pre-emergence damping-off of chickpea (Cicer arietinum L.). The pathogen has been controlled for several decades using the fungicide metalaxyl as seed treatment but has re-emerged as a severe problem with the detection of metalaxyl-resistant isolates of the pathogen from infested fields in the United States Pacific Northwest. The objective of this study was to identify genetic markers and candidate genes associated with resistance to P. ultimum in an interspecific recombinant inbred line population (CRIL-7) derived from a cross between C. reticulatum (PI 599072) x C. arietinum (FLIP 84-92C) and conduct genome-wide association studies (GWAS) for disease resistance using a chickpea diversity panel consisting of 184 accessions. CRIL-7 was examined using 1029 SNP markers spanning eight linkage groups. A major QTL, “qpsd4-1,” was detected on LG 4 that explained 41.8% of phenotypic variance, and a minor QTL, “qpsd8-1,” was detected on LG8 that explained 4.5% of phenotypic variance. Seven candidate genes were also detected using composite interval mapping including several genes previously associated with disease resistance in other crop species. A total of 302,902 single nucleotide polymorphic (SNP) markers were used to determine population structure and kinship of the diversity panel. Marker–trait associations were established by employing different combinations of principal components (PC) and kinships (K) in the FarmCPU model. Genome-wide association studies detected 11 significant SNPs and seven candidate genes associated with disease resistance. SNP Ca4_1765418, detected by GWAS on chromosome 4, was located within QTL qpsd4-1 that was revealed in the interspecific CRIL-7 population. The present study provides tools to enable MAS for resistance to P. ultimum and identified genomic domains and candidate genes involved in the resistance of chickpea to soilborne diseases.

Introduction

Chickpea (Cicer arietinum L) is one of the most historically significant field crops, being among the eight “founder crops” domesticated by Neolithic societies 8,000–12,000 years ago in the “Fertile Crescent” of present-day Iraq, Syria, and Turkey (Zohary et al., 2012). Its importance continues to this day, with more than 14.2 million tonnes of chickpea produced globally in 2019, making it the third most important pulse crop in terms of global production, after dry beans (Phaseolus vulgaris L.) and peas (Pisum sativum L.) (FAOSTAT, 2022). Chickpea is produced in more than 50 nations, with India being the largest producer, accounting for approximately 68% of global production (FAOSTAT, 2022). There are two major market classes of chickpea based on seed traits; “Desi” chickpea, which have a “teardrop”-shaped seed and a pigmented seed coat, and “Kabuli” chickpea, which has an “owl head” shape, a light beige or cream-colored seed coat, and is typically larger than the Desi chickpea (Toker, 2009).

In the United States, chickpea is primarily grown in dryland production systems in rotation with wheat or barley in the Pacific Northwest (Idaho and Washington) and Northern Plains (Montana and North Dakota). In 2019, more than 194,000 tonnes of chickpea were produced in the United States, with an estimated value greater than $116 million (NASS, 2022). However, several diseases challenge farmers in the United States, including Ascochyta blight caused by Ascochyta rabiei (Pass.) Lab. (Bayaa et al., 2011). Recently, seed rot and pre-emergent damping-off of chickpea caused by metalaxyl-resistant P. ultimum Trow has re-emerged as a significant disease in the Pacific Northwest.

The genus Pythium includes several soilborne species that cause seed and seedling diseases across a wide range of crops (Martin and Loper, 1999). P. ultimum causes seed rot, damping-off, and root rot in several legumes, including soybean (Glycine max L.), common bean (Phaseolus vulgaris L.), and pea (Pisum sativum L.) (Plaats-Niterink, 1981). Seed rot and pre-emergence damping-off of chickpea caused by P. ultimum were first detected in the United States in 1979 in Washington State (Kaiser and Hannan, 1983). Historically, these diseases have been controlled using pre-plant seed treatments containing metalaxyl or its stereoisomer mefenoxam (Casas et al., 1990). However, in 2014 isolates of P. ultimum var. ultimum with metalaxyl resistance were collected from unsprouted and rotten chickpea seeds obtained from a field in Washington exhibiting poor sprouting (Chen and Van Vleet, 2016). Subsequently, metalaxyl-resistant isolates of P. ultimum have been collected from several chickpea production fields in Idaho and Washington, where poor seedling sprouting was observed (Wang et al., 2020). Greenhouse and field tests showed that ethaboxam effectively manages metalaxyl-resistant P. ultimum, and commercial chickpea farmers now commonly apply a seed treatment containing both metalaxyl and ethaboxam for disease control (Wang et al., 2020). However, this increases production costs for farmers, and the effective use of ethaboxam depends on pathogen populations not developing resistance to the fungicide. Although resistance to ethaboxam in P. ultimum has not been detected, resistance has been detected in P. aphanidermatum (Edson) Fitzp. and P. deliense Meurs isolated from soybean (Glycine max L.) and dry bean (Phaseolus vulgaris L.), respectively. Increased production costs and concerns about the development of metalaxyl-resistant P. ultimum suggest that effective disease control in the future may require other approaches, including the use of disease-resistant chickpea cultivars.

Initial studies indicated that only small, dark-seeded Desi chickpeas were resistant to P. ultimum and all Kabuli chickpeas tested were susceptible (Kaiser and Hannan, 1983). We recently evaluated a collection of commercial chickpea cultivars and accessions from the United States National Plant Germplasm System (NPGS) and the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) for resistance to metalaxyl-resistant P. ultimum (Agarwal et al., 2020). The great majority of resistant accessions were Desi types and accessions with pigmented seed coats. Although the popular Kabuli cultivars ‘Sierra’ (Muehlbauer et al., 2004) and ‘Nash’ (Vandemark et al., 2015) were susceptible, three Kabuli accessions W625864, W625882, and W625884 were identified that were significantly more resistant than Sierra to two different isolates of metalaxyl-resistant P. ultimum (Agarwal et al., 2020). These results were promising because chickpea production in the United States is almost entirely composed of Kabuli types (Vandemark et al., 2014). Although these partially resistant accessions may be useful as parents for developing Kabuli cultivars with improved resistance to metalaxyl-resistant P. ultimum, it may also be possible to accelerate the development of resistant cultivars through the use of marker-assisted breeding approaches.

A range of genomic resources is available for chickpeas, including a draft sequence of the cultivated chickpea genome, which has an estimated size of approximately 738 Mb and contains 28,269 genes (Varshney et al., 2013). By using different marker genotyping platforms and molecular mapping approaches, significant associations have been identified between molecular markers and several diseases of chickpea, including Ascochyta blight (Tekeoglu et al., 2004; Sabbavarapu et al., 2013; Jendoubi et al., 2016; Garg et al., 2018; Mannur et al., 2019), Fusarium wilt (Winter et al., 2000; Sharma and Muehlbauer, 2007; Jingade and Ravikumar, 2015; Li et al., 2015; Meng et al., 2015; Mannur et al., 2019; Karadi et al., 2021), and dry root rot (Karadi et al., 2021). The objective of this study was to detect significant marker–trait associations and identify candidate genes for resistance in chickpea to metalaxyl-resistant P. ultimum.

Materials and Methods

Plant Material

A chickpea mapping population (CRIL-7) that included 177 chickpea recombinant inbred lines (RIL) derived from an interspecific cross C. reticulatum (PI 599072) x C. arietinum (FLIP 84–92C) was used to conduct QTL analysis. PI 599072 is a Desi type and is resistant to P. ultimum, whereas FLIP 84-92C is a disease-susceptible Kabuli type. The RILs were increased under greenhouse conditions by single seed descent to F₇ (Tekeoglu et al., 2000). Genome-wide association study was carried out on 184 taxonomically, morphologically, and geographically diverse accessions obtained from ICRISAT, Patancheru, India. This collection included 34 Kabuli, 144 Desi, and 6 pea-shaped accessions (Upadhyaya et al., 2001).

Disease Screening Assay and Resistance Scoring

P. ultimum strain PT410 (Agarwal et al., 2020) was used for all disease screening. The isolate was originally obtained from decaying chickpea seeds collected from a field in Patterson, WA, and its resistance to metalaxyl was confirmed through serial subculturing on media containing 50 ppm metalaxyl (Wang et al., 2020). The isolate was cultured and maintained on sucrose yeast extract agar in Petri plates at room temperature. CRIL-7 recombinant inbred lines along with the parental genotypes were screened for disease reaction to P. ultimum isolate (PT410) under controlled growth chamber conditions (12°C night–14°C days, 12 h day length) at Washington State University, Pullman, United States. Five seeds were planted in 10-cm pots containing 70 g of soil mix infested with 25,000 CFU of P. ultimum oospores for each entry. Myles and Sierra chickpea cultivars were used as resistant and susceptible checks, respectively (Agarwal et al., 2020). Pots were arranged in a completely randomized design. The number of seedlings that emerged from each pot was counted 14 days after planting. The experiment was repeated once.

Similarly, 184 accessions of mini-core collection from ICRISAT were evaluated for resistance to P. ultimum PT410 using the aforementioned methods. Again, the cultivars Sierra and Myles were included as susceptible and resistant controls. The experiment was repeated once. Results of these evaluations were previously reported (Agarwal et al., 2020).

Statistical Analysis

Statistical analysis was conducted using JMP 14 software (MP^®, Version <14>. SAS Institute Inc. Cary, NC, 1989–2019). Broad-sense heritability (H²) was calculated based on the average mean seedling emergence values from repeated experiments with R software (http://www.R-project.org/) using MME-based algorithms (Alexander and Lange, 2011).

Linkage Group Construction and Linkage QTL Mapping

Genomic DNA was extracted from young leaves of both parents and 177 RILs. Single nucleotide polymorphic (SNP) markers were detected using the Genotyping by Sequencing (GBS) approach. The libraries from the parental lines and RILs were prepared using MslI restriction enzymes and sequenced using the Illumina NextSeq 500 V² to generate 150 bp paired-end reads by LGC Company (https://www.lgcgroup.com/). For processing reads, demultiplexing of all library groups was done using Illumina bcl2fastq 2.17.1.14 software, followed by demultiplexing of library groups into samples according to their barcodes. Quality trimming was done by discarding low-quality reads with a final length <20 bases, and filtered data were used to call SNPs. The filtered, high-quality reads from each sample were aligned to the chickpea reference genome (Cicer arietinum CDC Frontier whole-genome assembly v1.0) (Varshney et al., 2013). The variant discovery was made using Freebayes v1.0.2-16 with the parameters min-base-quality: 10; min-supportingallele-qsum: 10; read-mismatch-limit: 3; min-coverage: 5; min-alternate count: 4, excluding unobserved genotypes; and mismatch-base-quality-threshold: 10. Variant filtering was done by removing markers with missing allele calls and minor allele frequency (MAF) < 0.05. Redundant markers were excluded from the analysis by implementing the BIN function in QTL IciMaping 4.1 (Meng et al., 2015). SNP markers with highly distorted segregation ratios at probability level (p ≤ 0.0001) were excluded. However, SNP that were slightly distorted (0.0001 ≤ p ≤ 0.05) from the Mendelian ratio were included in the linkage map.

These filtered markers were used to construct linkage groups (LGs) using the “Map” function in QTL IciMaping 4.1 and were assigned numbers (LG1–LG8) based on the genomic position of SNP markers. LGs with unlinked markers were removed from further construction. The remaining SNP were grouped with a logarithm-of-odds (LOD) threshold of 9.0. Recombination counting and ordering (RECORD) and “COUNT” (number of recombination events) algorithms were used in ordering and rippling. The linkage map and the best linear unbiased predictions (BLUP) value of phenotypic data of the CRIL-7 population from repeated experiments were used for QTL analysis. QTL was detected with composite inclusive composite interval mapping of additive (ICIM-ADD) function in QTL IciMaping 4.1 (Li et al., 2007). The threshold used to declare significant QTL was the permutation test with 1,000 permutations at the 0.05 significance level. Mapping parameters to detect additive QTL were set as step = 1.0 cM and PIN = 0.001 (PIN: the largest value for entering variables in stepwise regression of residual phenotype on marker variables). Parents for trait enhancing alleles were detected using the sign of the additive effects; the positive sign denotes that trait enhancing allele is from parent PI599072, whereas the negative sign indicates that the trait-enhancing allele is from FLIP 84-92C.

SNP Panel for GWAS

An SNP dataset from the Center of Excellence in Genomics & Systems Biology, ICRISAT that contained approximately 900,000 SNPs across 184 chickpea accessions was used for GWAS. First, SNP data were filtered by removing markers with more than 80% missing data and minor allele frequency smaller than 0.05. Duplicate markers and duplicate genotype samples in the dataset were then removed along with contigs and scaffolds using VCFtools (Hyun et al., 2008). Next, pairwise r² was calculated for all SNPs across each chromosome of the chickpea genome. SNPs with significant r² values (p < 0.001) were considered informative and were pruned using a linkage disequilibrium (LD) pruning method implemented in PLINK software v1.09 using “—indep-pairwise 50 5 0.5” command line in Linux (Purcell et al., 2007). The pruned set of SNP markers was used for association analysis.

Population Structure and Relatedness

Population structure and kinship were estimated using ADMIXTURE (v1.23) software. The ADMIXTURE tool uses a model-based algorithm to estimate the ancestry of unrelated individuals (Alexander and Lange, 2011). The number of underlying population groups (k) was estimated from 1 to 10 using the maximum likelihood estimation approach with a fast numerical optimization algorithm. A Q-matrix file representing the least number of population groups (k) was used for GWAS. Population structure was further estimated by principal component analysis (PCA) using the PLINK function. The EMMA algorithm embedded in the GAPIT package of R software was used to account for kinship (Hyun et al., 2008; Lipka et al., 2012). Finally, a dendrogram was generated using a neighbor-joining (NJ) algorithm to assess the relationship between mini-core accessions.

Marker–Trait Association Analysis

Association analysis was done using disease scores of the 184 chickpea accessions and 302,902 filtered SNPs. Kinship relatedness (K) was considered a random effect, and population structure based on the number of principal components (PC) that explained 25%–50% of the total phenotypic variance was considered a fixed effect. For marker–trait associations, models with different combinations of the population (PC)/admixture(Q) and family (K) structures were applied using the GAPIT package of R software: FarmCPU with Kinship (K), FarmCPU with Kinship (K) + PC(2), FarmCPU with Kinship (K) + PC(3), FarmCPU with Kinship (K) + PC(4), FarmCPU with Kinship (K) + PC(5), and FarmCPU with Kinship (K) + Admixture (Q) (Lipka et al., 2012). The best model was selected based on the mean squared difference (MSD) value between observed and expected p-values of all SNPs (Mamidi et al., 2011). The final Q–Q plots and Manhattan plots were created, and significant SNPs were calculated based on a p-value < 10^–5 and Bonferroni cut-off (p-value of 0.05/(the total number of SNP markers) (i.e., 0.05/302902 = 1.65 × 10^–7). Genes located within a 100-Kb region centered on a significant SNP were selected as candidate genes. The SNPeff tool (http://snpeff.sourceforge.net/) was used to detect genome coordinates of candidate genes and amino acid changes due to the SNPs (Cingolani et al., 2012).

Results

Phenotypic Evaluation for Resistance to Pythium ultimum.

A total of 177 CRIL-7 lines and both parents were screened for resistance to P. ultimum in a repeated experiment. A summary of descriptive statistics of disease reactions is presented in Supplementary Table S1. The mean of the susceptible check Sierra was 0.2 and the mean of the resistant check Myles was 4.1 in both experiments. The mean of PI 599072 was >2.5 in both experiments, indicating partial resistance to P. ultimum. The mean of FLIP 84-92C was ≤0.3 in both experiments, indicating susceptibility to P. ultimum. The means of all RILs in experiments 1 and 2 were 1.7 and 1.8, respectively. Transgressive segregants were observed among RILs in both experiments. RIL effects, experiment effects, and their interaction effects were all significant (Supplementary Table S1), with the greatest magnitude for the RIL effect and the least for the interaction effect. A broad-sense heritability estimate of 0.78 suggests that disease resistance is highly heritable (Supplementary Table S2).

SNP-Based Interspecific Genetic Map

A total of 65,112 SNP markers were obtained by GBS on CRIL-7 population [C. reticulatum (PI 599072) x C. arietinum (FLIP 84–92C)]. After filtering, 1,029 SNP markers were used to construct a linkage map and were assigned to the eight chickpea linkage groups (LG1–LG8) (Supplementary Figure S3). The eight linkage groups covered 1,186.30 cM (Supplementary Table S4). LG 8 was the smallest linkage group with 65 markers and a length of 50.59 cM, while the largest was LG 5 with 112 markers and a total length of 281.55 cM (Supplementary Table S3). Gaps in marker coverage of LGs are due to a high proportion of distorted segregation markers among the interspecific mapping population.

QTLs and Candidate Genes Detected by Linkage Mapping

Two QTLs were detected that were significantly associated with resistance to P. ultimum, QTL qpsd4-1 on LG 4, and QTL qpsd8-1 on LG 8. qpsd4-1 explained 46.5% of total variance with LOD = 25.2, whereas qpsd8-1 explained 4.5% of the total phenotypic variance with LOD = 3.3 (Table 1). qpsd4-1 had a positive additive effect value, indicating that positive alleles came from PI 599072, while the additive effect value was negative for qpsd8-1, indicating that negative alleles came from FLIP 84-92C.

The physical positions of flanking markers were used to identify putative genes associated with resistance to P. ultimum located within two QTLs, qpsd4-1 and qpsd8-1. Cicer arietinum cv. CDC Frontier (Kabuli) reference genome on the Pulse Crop database was used for identifying candidate genes (Varshney et al., 2013). Three candidate genes related to disease resistance were identified within the region flanking qpsd4-1 on LG 4, and four additional candidate genes were detected within the flanking qpsd8-1 region of LG 8 (Table 2).

Marker–Trait Association and Candidate Genes Through GWAS Analysis

For GWAS analysis, the SNP dataset on the ICRISAT accessions was filtered based on the minor allele frequencies (MAF <5%), missing data, and duplicate markers, and 229,965 SNPs were removed. Further linkage disequilibrium (LD) pruning removed 367,133 SNPs. After filtering and pruning, a total of 302,902 polymorphic SNPs remained for GWAS (Supplementary Figure S4).

PCA was performed to estimate population structure. The first two principal components (PC) explained 25% of the total variance, and the first five principal components (PC) explained 50% of the total variance. Graphs plotted using the first two PC explained the distribution of genotypes within different subpopulations (Figure 1A). In ADMIXTURE analysis with a total of 10 numbers of ancestral populations, the lowest cross-validation error was observed at K = 4, followed by k = 2 (with a minimum difference) were used as cofactors for GWAS (Figure 1B). Using K = 2, the panel was split into two subpopulations corresponding to Desi and Kabuli classes (Figure 1C), while with k = 4, genotypes were grouped according to seed shapes; angular, owl, and pea-shaped.

FIGURE 1

The greatest number of significant SNPs was observed on chromosome 4 (Ca4) (Figure 2). Marker Ca4_1765418 on chromosome 4 with a p-value of 1.65 × 10⁻⁷ was also detected within QTL qpsd4-1 that was revealed in the interspecific CRIL-7 population. Among the 35 significant SNPs, 29 were found in intergenic regions and 6 within genic regions (Table 3). Seven candidate genes that could be related to disease resistance were found on chromosomes 2, 4, 6, 7, and 8 within 100 kb flanking regions of significant SNPs (Table 4).

FIGURE 2

Discussion

GBS can be used to generate large-scale SNPs that are abundant in the genome allowing for the construction of high-density genetic maps and higher statistical power in association studies (Spencer et al., 2009; Yang et al., 2017). The CRIL-7 population has been used previously to develop linkage maps. These maps covered 981.6 cM and 1,174.4 cM across 9 linkage groups based on isozyme, inter simple sequence repeat (ISSR) loci, and RAPD markers (Santra et al., 2000b; Tekeoglu et al., 2002). In the present study, SNPs were used to resolve eight linkage groups spanning 1,186.30 cM.

QTLs and Candidate Genes Detected by Linkage Mapping

A major QTL for resistance to metalaxyl-resistant P. ultimum, qpsd4-1, was detected on LG 4 (Table 1). Previously, CRIL-7 has been evaluated for Ascochyta blight resistance, and QTLs associated with disease resistance were detected on LG 1, LG 4, and LG 8 (Santra et al., 2000a; Tekeoglu et al., 2004; Kumar et al., 2018). Two significant clusters of QTLs (QTLAR1 and QTLAR2) associated with resistance to Ascochyta blight have also been detected on LG 4 based on an analysis of RILs from a cross of cultivars Amit and ICCV 96029 (Deokar et al., 2019). Additionally, two QTLs (AB-Q-SR-4-1 and AB-Q-SR-4-2) on LG 4 associated with Ascochyta blight disease resistance were detected using mapping population (C 214’ × ‘ILC 3279) (Gisi and Sierotzki, 2008). These results and our detection of a major QTL for disease resistance on LG 4 suggest that genes for resistance to diverse chickpea pathogens are located on LG 4.

A total of seven genes associated with disease resistance were found in the QTL region on LG 4 and LG 8. Three candidate genes, Ca_07798, Ca_07797, and Ca_07799, were identified in QTL qpsd4-1. Ca_07798 is a JmjC domain-containing protein D, which is a group of histone lysine demethylases. This protein positively regulates rice defense against bacterial blight pathogen Xanthomonas oryzae pv. oryzae by epigenetically suppressing negative defense regulator H3K4me2/3 (Hou et al., 2015). The Ca_07797 gene belongs to a WD-repeat family that is involved in plant innate immune signaling pathway. Studies on maize showed that WDR-containing TTG1 protein–induced resistance against leaf blights (Ibraheem et al., 2010; Ibraheem et al., 2015). In tobacco, the interaction of TTG1-WDR with an elicitin protein (ParA1) from a pathogenic oomycete Phytophthora parasitica var. nicotianae-activated plant immune responses, including the generation of reactive oxygen species and programmed cell death (Wang et al., 2009). The Ca_07799 gene belongs to a zinc finger protein family that has been shown to play a critical role in disease resistance across many plant species (Gupta et al., 2012).

The Ca_02390 gene found on LG 8 encodes for 1-aminocyclopropane-1-carboxylate synthase. This enzyme initiates the conversion of S-adenosyl-l-methionine (SAM) into 1-aminocyclopropane-1-carboxylate (ACC), which is the precursor of ethylene and acts as a signaling molecule to regulate plant growth and reduce stress response (Polko and Kieber, 2019). The role of ACC and ethylene biosynthesis in plant defense against bacterial pathogens, including Erwinia carotovora subsp. carotovora and Pseudomonas syringae have been studied in Arabidopsis thaliana (L.) Heynh mutants where plants with reduced ACC production showed greater disease susceptibility (Norman-Setterblad et al., 2000; Guan et al., 2015). Recently, upregulation of genes involved in ethylene biosynthesis was detected in resistant apple seedling reactions to infection by P. ultimum (Zhu and Saltzgiver, 2020). Ca_02384 on LG 8 encodes an AT-hook DNA-binding protein that binds to minor groove DNA and alters gene expression. Genes such as AHL19 encoding an AT-hook DNA-binding protein are associated with enhanced disease resistance in A. thaliana to Verticillium wilt caused by V. dahliae, V. albo-atrum, and V. longisporum (Yadeta et al., 2011). The Ca_02383 gene on LG 8 belongs to the family of multidrug and toxic compound extrusion (MATE) transporters associated with plant disease resistance during pathogen interaction. The expression of MATE genes in plants is induced by pathogen attack (Sun et al., 2011). Members of the MATE family such as enhanced disease susceptibility 5 (EDS5) and activated disease susceptibility 1 (ADS1) function as negative regulators of plant immune systems by reducing basal resistance during pathogen interaction or by negatively regulating the accumulation of salicylic acid and pathogenesis-related 1 (PR1) gene expression (Nawrath et al., 2002; Sun et al., 2011). Another candidate gene on LG 8, Ca_02389, belongs to the protein kinase family. Members of this family have also been shown to be upregulated in response to pathogens, for example, Xanthomonas oryzae pv. oryzicola (Xoc) in Oryza sativa and Pseudomonas syringae pv tomato (pto) in Solanum lycopersicum (Martin et al., 1993; Xu et al., 2013). Recently, GWAS of common bean identified several candidate genes associated with resistance to P. ultimum, including genes for protein kinase superfamily proteins and MAPK/ERK kinase 1 (Dramadri et al., 2020).

GWAS Analysis of Resistance in Chickpea to P. ultimum

Different statistical models were deployed using FarmCPU to assess population structure and kinship in the chickpea diversity panel. A combination of these models using FarmCPU separates a mixed linear model (MLM) into a random effect and a fixed-effect model, which reduces false positives and false negatives caused by kinship and population structure and gives highly significant SNP markers (Lipka et al., 2012). The best model was selected based on MSD value (Supplementary Table S5), with a low MSD value indicating less deviation from the expected distribution of p-values, signifying lower type I error of the selected model. This study identified a total of seven candidate gene mapping to 11 loci associated with resistance in chickpea to P. ultimum. PCA and ADMIXTURE analyses revealed two major groups within the core collections (Figure 1) corresponding to Desi and Kabuli market classes.

GWAS identified many SNPs associated with disease-resistance–related traits. Gene Ca_19996 encodes cellulose synthase–like protein, which inhibits the progress of the fungal penetration peg during powdery mildew infection caused by Blumeria graminis f. sp. hordei (Bgh) in barley (Hordeum vulgare L.) (Douchkov et al., 2016). Gene Ca_09957 encodes a histidine kinase protein involved in seed maturation and disease resistance against fungal and bacterial pathogens (Pham et al., 2012). Ca_17277 belongs to the O-methyltransferase family of enzymes that play a significant role in plant stress and disease resistance. Studies on corn (Zea mays L.) and wheat (Triticum aestivum L.) demonstrated that caffeoyl-CoA O-methyltransferase conferred resistance against southern leaf blight, gray leaf spot, and sharp eyespot disease, respectively (Yang et al., 2017; Wang et al., 2018). Ca_04625 encodes a thiamine thiazole synthase that increases resistance to fungal pathogens by enhancing anti-oxidative capacity and inducing systemic acquired resistance (SAR) in diverse plant species, including Oryza sativa L. A. thaliana (L.) Heynh., Nicotiana sp., and Cucumis sativus L. (Goyer, 2010). Ca_22742 encodes an ethylene-responsive transcription factor 1–like protein. Ethylene-responsive transcription factors play a critical role in the plant defense system by regulating pathogenesis-related (PR) gene expression, including effectors GmERF5 and GmERF113, and contribute to resistance against root and seed rot caused by Oomycete pathogens, including Phytopthora nicotianae and Py. sojae (Goyer, 2010; Yang et al., 2017). The Ca_17276 gene encodes LUPR1 protein that is upregulated in response to the Hyaloperonospora parasitica (LURP) cluster in part of the A. thaliana. The LUPR1 gene has been associated with resistance to oomycetes Hyaloperonospora parasitica and Py. infestans (Dong et al., 2015). The Ca_10436 gene encodes a calmodulin-binding protein, which has activated and enhanced resistance to a broad spectrum of pathogens in Nicotiana tabacum and A. thaliana (Zhao et al., 2017). Additionally, calmodulin-binding proteins in A. thaliana and Hordeum vulgare have been shown to confer resistance to powdery mildew by interacting with MLO (powdery mildew-resistance gene o) protein (Heo et al., 1999).

Conclusion

In this study, we used an interspecific chickpea population to identify one major and one minor QTL associated with resistance to P. ultimum. We also identified 35 SNPs and 14 candidate genes associated with disease resistance based on the GWAS of a chickpea diversity panel. SNP Ca4_1765418, detected by GWAS on chromosome 4, was located within QTL qpsd4-1 that was revealed in the interspecific CRIL-7 population. These findings suggest this region of the genome should be examined more closely to identify genes conditioning disease resistance. Significant QTLs must be validated in different chickpea populations before the markers can be widely used for breeding. The present study provides tools to enable MAS for resistance to P. ultimum and identified genomic domains and candidate genes involved in the resistance of chickpea to soilborne diseases.

References

• 1

  AgarwalC.ChenW.CoyneC.VandemarkG. (2020). Identifying Sources of Resistance in Chickpea to Seed Rot and Seedling Damping‐off Caused by Metalaxyl‐resistant Pythium Ultimum, [Internet]. Crop Sci.61, 1739. [cited 2021 Jan 20];csc2.20424. 10.1002/csc2.20424

  • CrossRef
  • Google Scholar

• 2

  AlexanderD. H.LangeK. (2011). Enhancements to the ADMIXTURE Algorithm for Individual Ancestry Estimation. BMC Bioinforma.12 (1), 246. [cited 2020 Apr 20]Available from: https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-12-246. 10.1186/1471-2105-12-246

  • CrossRef
  • Google Scholar

• 3

  BayaaB.ChenW. (2011). “Ascochyta Blight of Chickpea,” in Compendium of Chickpea and Lentil Diseases and Pests. Editors ChenW.SharmaH. C.MuehlbauerF. J. (St. Paul, MN: APS Press).

  • Google Scholar

• 4

  CasasA. T.KaiserW. J.IngramD. M. (1990). Control of Pythium Seed Rot and Preemergence Damping-Off of Chickpea in the U.S. Pacific Northwest and Spain. Plant Dis.74, 563–569. 10.1094/pd-74-0563

  • CrossRef
  • Google Scholar

• 5

  ChenW.Van VleetS. (2016). WSU - Chickpea Damping-Off Due to Metalaxyl-Resistant Pythium: An Emerging Disease in the Palouse. [Internet]. Washington: Washingtom State University Extension. Pullman, Washington : Washington State University. Extension., 1–4. [cited 2020 Apr 20]Available from: https://pubs.wsu.edu/ItemDetail.aspx?ProductID=15853&SeriesCode=&CategoryID=150&Keyword=.

  • Google Scholar

• 6

  CingolaniP.PlattsA.Wangle L.CoonM.NguyenT.WangL.et al (2012). A Program for Annotating and Predicting the Effects of Single Nucleotide Polymorphisms, SnpEff: SNPs in the Genome of Drosophila melanogaster Strain W1118; Iso-2; Iso-3. Fly. (Austin)6 (2), 80–92. 10.4161/fly.19695

  • CrossRef
  • Google Scholar

• 7

  DeokarA.SagiM.Tar’anB. (2019). Genome-wide SNP Discovery for Development of High-Density Genetic Map and QTL Mapping of Ascochyta Blight Resistance in Chickpea (Cicer Arietinum L.). Theor. Appl. Genet.132, 1861–1872. 10.1007/s00122-019-03322-3

  • CrossRef
  • Google Scholar

• 8

  DongL.ChengY.WuJ.ChengQ.LiW.FanS.et al (2015). Overexpression of GmERF5, a New Member of the Soybean EAR Motif-Containing ERF Transcription Factor, Enhances Resistance to Phytophthora Sojae in Soybean. J. Exp. Bot.66, 2635–2647. 10.1093/jxb/erv078

  • CrossRef
  • Google Scholar

• 9

  DouchkovD.LueckS.HenselG.KumlehnJ.RajaramanJ.JohrdeA.et al (2016). The Barley (Hordeum Vulgare) Cellulose Synthase-like D2 Gene (HvCslD2) Mediates Penetration Resistance to Host-Adapted and Nonhost Isolates of the Powdery Mildew Fungus, [Internet]. New Phytol.212 (2), 421–433. [cited 2021 Oct 24]. 10.1111/nph.14065

  • CrossRef
  • Google Scholar

• 10

  DramadriI. O.AmongiW.KellyJ. D.Mugisha MukankusiC. (2020). “Genome‐wide Association Analysis of Resistance to Pythium Ultimum in Common Bean ( Phaseolus vulgaris ),” [Internet] in Plant Breed. Editor LinkW. (New Jersey: Wiley). [cited 2020 Nov 17];pbr.12855. 10.1111/pbr.12855

  • CrossRef
  • Google Scholar

• 11

  GargT.MallikarjunaB. P.ThudiM.SamineniS.SinghS.SandhuJ. S.et al (2018). Identification of QTLs for Resistance to Fusarium Wilt and Ascochyta Blight in a Recombinant Inbred Population of Chickpea (Cicer Arietinum L.). Euphytica214 (3), 45. [cited 2020 Nov 16]. 10.1007/s10681-018-2125-3

  • CrossRef
  • Google Scholar

• 12

  GisiU.SierotzkiH. (2008). Fungicide Modes of Action and Resistance in Downy Mildews. Eur. J. Plant Pathol.122 (1), 157–167. 10.1007/s10658-008-9290-5

  • CrossRef
  • Google Scholar

• 13

  GoyerA. (2010). Thiamine in Plants: Aspects of its Metabolism and Functions. Phytochemistry71, 1615–1624. 10.1016/j.phytochem.2010.06.022

  • CrossRef
  • Google Scholar

• 14

  GuanR.SuJ.MengX.LiS.LiuY.XuJ.et al (2015). Multilayered Regulation of Ethylene Induction Plays a Positive Role in Arabidopsis Resistance against Pseudomonas syringae. Plant Physiol.169 (1), 299–312. 10.1104/pp.15.00659

  • CrossRef
  • Google Scholar

• 15

  GuptaS. K.RaiA. K.KanwarS. S.SharmaT. R. (2012). Comparative Analysis of Zinc Finger Proteins Involved in Plant Disease Resistance, [Internet]. PLoS One7 (8), e42578. [cited 2020 Apr 20]Available from: http://www.ncbi.nlm.nih.gov/pubmed/22916136. 10.1371/journal.pone.0042578

  • CrossRef
  • Google Scholar

• 16

  HeoW. D.LeeS. H.KimM. C.KimJ. C.ChungW. S.ChunH. J.et al (1999). Involvement of Specific Calmodulin Isoforms in Salicylic Acid-independent Activation of Plant Disease Resistance Responses. Proc. Natl. Acad. Sci. U. S. A.96 (2), 766–771. 10.1073/pnas.96.2.766

  • CrossRef
  • Google Scholar

• 17

  HouY.WangL.WangL.LiuL.LiL.SunL.et al (2015). JMJ704 Positively Regulates Rice Defense Response against Xanthomonas Oryzae Pv. Oryzae Infection via Reducing H3K4me2/3 Associated with Negative Disease Resistance Regulators. BMC Plant Biol.15 (1), 286. 10.1186/s12870-015-0674-3

  • CrossRef
  • Google Scholar

• 18

  HyunM. K.ZaitlenN. A.WadeC. M.KirbyA.HeckermanD.DalyM. J.et al (2008). Efficient Control of Population Structure in Model Organism Association Mapping, [Internet]. Genetics178 (3), 1709–1723. [cited 2020 Apr 20]Available from: http://www.ncbi.nlm.nih.gov/pubmed/18385116. 10.1534/genetics.107.080101

  • CrossRef
  • Google Scholar

• 19

  IbraheemF.GaffoorI.ChopraS. (2010). Flavonoid Phytoalexin-dependent Resistance to Anthracnose Leaf Blight Requires a Functional Yellow Seed1 in Sorghum Bicolor. Genetics184 (4), 915–926. 10.1534/genetics.109.111831

  • CrossRef
  • Google Scholar

• 20

  IbraheemF.GaffoorI.TanQ.ShyuC.-R.ChopraS. (2015). A Sorghum MYB Transcription Factor Induces 3-deoxyanthocyanidins and Enhances Resistance against Leaf Blights in Maize. Molecules20 (2), 2388–2404. 10.3390/molecules20022388

  • CrossRef
  • Google Scholar

• 21

  JendoubiW.BouhadidaM.MillanT.KharratM.GilJ.RubioJ.et al (2016). Identification of the Target Region Including the Foc0 1/foc0 1 Gene and Development of Near Isogenic Lines for Resistance to Fusarium Wilt Race 0 in Chickpea, [Internet]. Euphytica210, 119–133. [cited 2020 Nov 17]. 10.1007/s10681-016-1712-4

  • CrossRef
  • Google Scholar

• 22

  JingadeP.RavikumarR. L. (2015). Development of Molecular Map and Identification of QTLs Linked to Fusarium Wilt Resistance in Chickpea, [Internet]. J. Genet.94, 723–729. [cited 2020 Nov 17]Available from: http://www.ias.ac.in/. 10.1007/s12041-015-0589-7

  • CrossRef
  • Google Scholar

• 23

  KaiserW. J.HannanR. M. (1983). Etiology and Control of Seed Decay and Preemergence Damping-Off of Chickpea byPythium Ultimum. Plant Dis.67 (1), 77. [cited 2020 Apr 19]Available from: http://agris.fao.org/agris-search/search.do?recordID=US19900085897. 10.1094/pd-67-77

  • CrossRef
  • Google Scholar

• 24

  KaradiA.SamineniS.SajjaS.SharmaM.ThudiM.MallikarjunaB. P.et al (2021). Molecular Mapping of Dry Root Rot Resistance Genes in Chickpea (Cicer Arietinum L.), [Internet]. Euphytica217, 123. [cited 2021 Oct 24]Available from: https://research-repository.uwa.edu.au/en/publications/molecular-mapping-of-dry-root-rot-resistance-genes-in-chickpea-ci. 10.1007/s10681-021-02854-4

  • CrossRef
  • Google Scholar

• 25

  KumarK.PurayannurS.KaladharV. C.ParidaS. K.VermaP. K. (2018). mQTL-Seq and Classical Mapping Implicates the Role of an AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) Family Gene in Ascochyta Blight Resistance of Chickpea. Plant Cell Environ.41 (9), 2128–2140. 10.1111/pce.13177

  • CrossRef
  • Google Scholar

• 26

  LiH.YeG.WangJ. (2007). A Modified Algorithm for the Improvement of Composite Interval Mapping, [Internet]. Genetics175 (1), 361–374. [cited 2022 Apr 26]Available from: https://academic.oup.com/genetics/article/175/1/361/6062043. 10.1534/genetics.106.066811

  • CrossRef
  • Google Scholar

• 27

  LiH.RoddaM.GnanasambandamA.AftabM.ReddenR.HobsonK.et al (2015). Breeding for Biotic Stress Resistance in Chickpea: Progress and Prospects. Euphytica204 (2), 257–288. 10.1007/s10681-015-1462-8

  • CrossRef
  • Google Scholar

• 28

  LipkaA. E.TianF.WangQ.PeifferJ.LiM.BradburyP. J.et al (2012). GAPIT: Genome Association and Prediction Integrated Tool. Bioinforma. [Internet]28 (18), 2397–2399. 10.1093/bioinformatics/bts444

  • CrossRef
  • Google Scholar

• 29

  MamidiS.ChikaraS.GoosR. J.HytenD. L.AnnamD.MoghaddamS. M.et al (2011). Genome-Wide Association Analysis Identifies Candidate Genes Associated with Iron Deficiency Chlorosis in Soybean, [Internet]. Plant Genome4 (3), 154–164. 10.3835/plantgenome2011.04.0011

  • CrossRef
  • Google Scholar

• 30

  MannurD. M.BabbarA.ThudiM.SabbavarapuM. M.RoorkiwalM.YeriS. B.et al (2019). Super Annigeri 1 and Improved JG 74: Two Fusarium Wilt-Resistant Introgression Lines Developed Using Marker-Assisted Backcrossing Approach in Chickpea (Cicer Arietinum L.), [Internet]. Mol. Breed.39 (1), 2. [cited 2020 Nov 17]. 10.1007/s11032-018-0908-9

  • CrossRef
  • Google Scholar

• 31

  MartinF. N.LoperJ. E. (1999). Soilborne Plant Diseases Caused by Pythium spp.: Ecology, Epidemiology, and Prospects for Biological Control, [Internet]. CRC Crit. Rev. Plant Sci.18, 111–181. [cited 2021 Jan 20]Available from: https://www.tandfonline.com/action/journalInformation?journalCode=bpts20.

  • Google Scholar

• 32

  MartinG. B.BrommonschenkelS. H.ChunwongseJ.FraryA.GanalM. W.SpiveyR.et al (1993). Map-based Cloning of a Protein Kinase Gene Conferring Disease Resistance in Tomato. Science262 (5138), 1432–1436. 10.1126/science.7902614

  • CrossRef
  • Google Scholar

• 33

  MengL.LiH.ZhangL.WangJ. (2015). QTL IciMapping: Integrated Software for Genetic Linkage Map Construction and Quantitative Trait Locus Mapping in Biparental Populations. Crop J.3 (3), 269–283. 10.1016/j.cj.2015.01.001

  • CrossRef
  • Google Scholar

• 34

  MuehlbauerF. J.TempleS. R.ChenW. (2004). Registration of ‘Sierra’ Chickpea, [Internet]. Crop Sci.44, 1864. [cited 2020 Apr 20]. 10.2135/cropsci2004.1864

  • CrossRef
  • Google Scholar

• 35

  NawrathC.HeckS.ParinthawongN.MétrauxJ.-P. (2002). EDS5, an Essential Component of Salicylic Acid-dependent Signaling for Disease Resistance in Arabidopsis, Is a Member of the MATE Transporter Family. Plant Cell14 (1), 275–286. 10.1105/tpc.010376

  • CrossRef
  • Google Scholar

• 36

  Norman-SetterbladC.VidalS.PalvaE. T. (2000). Interacting Signal Pathways Control Defense Gene Expression in Arabidopsis in Response to Cell Wall-Degrading Enzymes from Erwinia Carotovora, [Internet]. Mol. Plant-Microbe Interact.13, 430–438. [cited 2020 Apr 20]. 10.1094/MPMI.2000.13.4.430

  • CrossRef
  • Google Scholar

• 37

  PhamJ.LiuJ.BennettM. H.MansfieldJ. W.DesikanR. (2012). Arabidopsis Histidine Kinase 5 Regulates Salt Sensitivity and Resistance against Bacterial and Fungal Infection. New Phytol.194 (1), 168–180. 10.1111/j.1469-8137.2011.04033.x

  • CrossRef
  • Google Scholar

• 38

  Plaats-NiterinkA. van der. (1981). Monograph of the Genus Pythium[Internet]. [cited 2020 Apr 18] Available from: http://www.wi.knaw.nl/publications/1021/content_files/content.htm

  • Google Scholar

• 39

  PolkoJ. K.KieberJ. J. (2019). 1-Aminocyclopropane 1-Carboxylic Acid and its Emerging Role as an Ethylene-independent Growth Regulator. Front. Plant Sci.10, 1602. 10.3389/fpls.2019.01602

  • CrossRef
  • Google Scholar

• 40

  PurcellS.NealeB.Todd-BrownK.ThomasL.FerreiraM. A. R.BenderD.et al (2007). PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet.81 (3), 559–575. 10.1086/519795

  • CrossRef
  • Google Scholar

• 41

  SabbavarapuM. M.SharmaM.ChamarthiS. K.SwapnaN.RathoreA.ThudiM.et al (2013). Molecular Mapping of QTLs for Resistance to Fusarium Wilt (Race 1) and Ascochyta Blight in Chickpea (Cicer Arietinum L.). Euphytica193 (1), 121–133. 10.1007/s10681-013-0959-2

  • CrossRef
  • Google Scholar

• 42

  SantraD. K.A4 - TekeogluM.A4 - RatnaparkheM.A4 - KaiserM.KaiserW. J.MuehlbauerF. J.et al (2000). Identification and Mapping of QTLs Conferring Resistance to Ascochyta Blight in Chickpea. Crop Sci.40 (6), 1606–1612. 10.2135/cropsci2000.4061606x

  • CrossRef
  • Google Scholar

• 43

  SantraD. K.TekeogluM.RatnaparkheM.KaiserW. J.MuehlbauerF. J. (2000). Identification and Mapping of QTLs Conferring Resistance to Ascochyta Blight in Chickpea. Crop Sci.40 (6), 1606–1612. 10.2135/cropsci2000.4061606x

  • CrossRef
  • Google Scholar

• 44

  SharmaK. D.MuehlbauerF. J. (2007). Fusarium Wilt of Chickpea: Physiological Specialization, Genetics of Resistance and Resistance Gene Tagging. Euphytica157, 1–14. 10.1007/s10681-007-9401-y

  • CrossRef
  • Google Scholar

• 45

  SpencerC. C.SuZ.DonnellyP.MarchiniJ. (2009). Designing Genome-wide Association Studies: Sample Size, Power, Imputation, and the Choice of Genotyping Chip. PLoS Genet.5 (5), e1000477. 10.1371/journal.pgen.1000477

  • CrossRef
  • Google Scholar

• 46

  SunX.GilroyE. M.ChiniA.NurmbergP. L.HeinI.LacommeC.et al (2011). ADS1 Encodes a MATE-Transporter that Negatively Regulates Plant Disease Resistance. New Phytol.192 (2), 471–482. 10.1111/j.1469-8137.2011.03820.x

  • CrossRef
  • Google Scholar

• 47

  TekeogluM.IşıkM.MuehlbauerF. (2004). QTL Analysis of Ascochyta Blight Resistance in Chickpea. Turkish J Agric28, 183–187.

  • Google Scholar

• 48

  TekeogluM.RajeshP. N.MuehlbauerF. J. (2002). Integration of Sequence Tagged Microsatellite Sites to the Chickpea Genetic Map, [Internet]. Theor. Appl. Genet.105 (6–7), 847–854. [cited 2020 Apr 19]Available from: http://www.ncbi.nlm.nih.gov/pubmed/12582909. 10.1007/s00122-002-0993-2

  • CrossRef
  • Google Scholar

• 49

  TekeogluM.SantraD. K.KaiserW. J.MuehlbauerF. J. (2000). Ascochyta Blight Resistance Inheritance in Three Chickpea Recombinant Inbred Line Populations. Crop Sci.40, 1251–1256. 10.2135/cropsci2000.4051251x

  • CrossRef
  • Google Scholar

• 50

  TokerC. (2009). A Note on the Evolution of Kabuli Chickpeas as Shown by Induced Mutations in Cicer Reticulatum Ladizinsky. Genet. Resour. Crop Evol.56, 7–12. 10.1007/s10722-008-933

  • CrossRef
  • Google Scholar

• 51

  UpadhyayaH. D.BramelP. J.SinghS. (2001). Development of a Chickpea Core Subset Using Geographic Distribution and Quantitative Traits. Crop Sci.41 (1), 206–210. 10.2135/cropsci2001.411206x

  • CrossRef
  • Google Scholar

• 52

  VandemarkG.GuyS. O.ChenW.McPheeK.PfaffJ.LauverM.et al (2015). Registration of 'Nash' Chickpea. J. Plant Registrations9, 275–278. 10.3198/jpr2014.07.0047crc

  • CrossRef
  • Google Scholar

• 53

  VandemarkG. J.BrickM. A.OsornoJ. M.KellyJ. D.UrreaC. A. (2014). “Edible Grain Legumes,” in Madison, WI: Yield Gains in Major U.S. Field Crops. Editors SmithS.DiersB.Specht BCJ. (US: ASA, CSSA, SSSA), 87–124.

  • Google Scholar

• 54

  VarshneyR. K.SongC.SaxenaR. K.AzamS.YuS.SharpeA. G.et al (2013). Draft Genome Sequence of Chickpea (Cicer Arietinum) Provides a Resource for Trait Improvement. Nat. Biotechnol.31 (3), 240–246. 10.1038/nbt.2491

  • CrossRef
  • Google Scholar

• 55

  WangM.Van VleetS.McGeeR.PaulitzT. C.PorterL. D.VandemarkG.et al (2020). Chickpea Seed Rot and Damping-Off Caused by Metalaxyl-Resistant Pythium Ultimum and its Management with Ethaboxam, [Internet]. Plant Dis. 105, 1728. [cited 2021 Jan 20]; Available from: https://apsjournals.apsnet.org/doi/abs/10.1094/PDIS-08-20-1659-RE. 10.1094/PDIS-08-20-1659-RE

  • CrossRef
  • Google Scholar

• 56

  WangM.ZhuX.WangK.LuC.LuoM.ShanT.et al (2018). A Wheat Caffeic Acid 3-O-Methyltransferase TaCOMT-3D Positively Contributes to Both Resistance to Sharp Eyespot Disease and Stem Mechanical Strength. Sci. Rep.8 (1), 6543–6614. 10.1038/s41598-018-24884-0

  • CrossRef
  • Google Scholar

• 57

  WangY.LiuR.ChenL.WangY.LiangY.WuX.et al (2009). Nicotiana tabacumTTG1 Contributes to ParA1-Induced Signalling and Cell Death in Leaf Trichomes. J. Cell Sci.122 (15), 2673–2685. 10.1242/jcs.049023

  • CrossRef
  • Google Scholar

• 58

  WinterP.Benko-IsepponA.-M.HüttelB.RatnaparkheM.TulluA.SonnanteG.et al (2000). A Linkage Map of the Chickpea (Cicer Arietinum L.) Genome Based on Recombinant Inbred Lines from a C. arietinum×C. reticulatum Cross: Localization of Resistance Genes for fusarium Wilt Races 4 and 5. Theor. Appl. Genet.101 (7), 1155–1163. 10.1007/s001220051592

  • CrossRef
  • Google Scholar

• 59

  XuM.-R.HuangL.-Y.ZhangF.ZhuL.-H.ZhouY.-L.LiZ.-K. (2013). Genome-Wide Phylogenetic Analysis of Stress-Activated Protein Kinase Genes in Rice (OsSAPKs) and Expression Profiling in Response to Xanthomonas Oryzae Pv. Oryzicola Infection. Plant Mol. Biol. Rep.31 (4), 877–885. 10.1007/s11105-013-0559-2

  • CrossRef
  • Google Scholar

• 60

  YadetaK. A.HanemianM.SmitP.HiemstraJ. A.PereiraA.MarcoY.et al (2011). The Arabidopsis thaliana DNA-Binding Protein AHL19 Mediates Verticillium Wilt Resistance. MPMI24 (12), 1582–1591. 10.1094/mpmi-04-11-0090

  • CrossRef
  • Google Scholar

• 61

  YangQ.HeY.KabahumaM.ChayaT.KellyA.BorregoE.et al (2017). A Gene Encoding Maize Caffeoyl-CoA O-Methyltransferase Confers Quantitative Resistance to Multiple Pathogens. Nat. Genet.49 (9), 1364–1372. 10.1038/ng.3919

  • CrossRef
  • Google Scholar

• 62

  ZhaoY.ChangX.QiD.DongL.WangG.FanS.et al (2017). A Novel Soybean ERF Transcription Factor, GmERF113, Increases Resistance to Phytophthora Sojae Infection in Soybean. Front. Plant Sci.8, 299. Available from: https://www.frontiersin.org/article/10.3389/fpls.2017.00299. 10.3389/fpls.2017.00299

  • CrossRef
  • Google Scholar

• 63

  ZhuY.SaltzgiverM. (2020). A Systematic Analysis of Apple Root Resistance Traits to Pythium Ultimum Infection and the Underpinned Molecular Regulations of Defense Activation, [Internet]. Hortic. Res.7, 62. [cited 2021 Oct 25]Available from: https://www.nature.com/articles/s41438-020-0286-4. 10.1038/s41438-020-0286-4

  • CrossRef
  • Google Scholar

• 64

  ZoharyD.HopfM.WeissE. (2012). Domestication of Plants in the Old World: The Origin and Spread of Domesticated Plants in Southwest Asia, Europe, and the Mediterranean Basin. [Internet]. 4th ed. Oxford: Oxford University Press, 264. Available from: https://www.oxfordscholarship.com/10.1093/acprof:osobl/9780199549061.001.0001/acprof-9780199549061.

  • Google Scholar