Sudhir Navathe1
Xinyao He2
Umesh Kamble3
Manjeet Kumar4
Madhu Patial4
Gyanendra Singh3
Gyanendra Pratap Singh5
Arun Kumar Joshi6
Pawan Kumar Singh2*- 1Genetics and Plant Breeding Group, Agharkar Research Institute, Pune, India
- 2Global Wheat Program, International Maize and Wheat Improvement Centre (CIMMYT), Texcoco, Mexico
- 3Division of Crop Improvement, ICAR-Indian Institute of Wheat and Barley Research, Karnal, India
- 4Division of Genetics, ICAR-Indian Agricultural Research Institute, New Delhi, India
- 5Indian Council of Agricultural Research (ICAR)-National Bureau of Plant Genetic Resources, New Delhi, India
- 6International Maize and Wheat Improvement Centre (CIMMYT) & Borlaug Institute for South Asia (BISA), New Delhi, India
The leaf blight diseases, Septoria nodorum blotch (SNB), and tan spot (TS) are emerging due to changing climatic conditions in the northern parts of India. We screened 296 bread wheat cultivars released in India over the past 20 years for seedling resistance against SNB (three experiments) and TS (two experiments). According to a genome-wide association study, six QTLs on chromosome arms 1BL, 2AS, 5BL, and 6BL were particularly significant for SNB across all three years, of which Q.CIM.snb.1BL, Q.CIM.snb.2AS1, Q.CIM.snb.2AS.2, and Q.CIM.snb.6BL appeared novel. In contrast, those on 5BS and 5BL may correspond to Snn3 and Tsn1, respectively. The allelic combination of tsn1/snn3 conferred resistance to SNB, whereas that of Tsn1/Snn3 conferred high susceptibility. As for TS, Tsn1 was the only stably significant locus identified in this panel. Several varieties like PBW 771, DBW 277, and HD 3319, were identified as highly resistant to both diseases that can be used in future wheat improvement programs as resistant donors.
1 Introduction
Wheat (Triticum spp.) is a staple food and a crucial element of global food security. However, of several factors, fungal diseases are the most important that limit wheat production. Septoria nodorum blotch (SNB) is caused by the necrotrophic fungal pathogen Parastagonospora nodorum (syn. Phaeosphaeria nodorum [E. Müll.], syn. Leptosphaeria nodorum [E. Müll.], syn. Stagonospora nodorum [Berk.], syn. Septoria nodorum [Berk.]). SNB frequently co-occurs with other necrotrophic fungal diseases like tan spot (TS, caused by Pyrenophora tritici-repentis) and Septoria tritici blotch (STB, caused by Zymoseptoria tritici). This disease is common in areas that experience frequent or high rainfall during the wheat growing season, such as Australia, Canada, Scandinavia, Central and Eastern Europe, eastern USA, and South America (Bearchell et al., 2005; Solomon et al., 2006; Shaw et al., 2008). In India, the disease was first recorded in the Nilgiri hills of south India (Chona and Munjal, 1952) and later from the Kumaon hills in Northern India (Joshi et al., 1971). In the last few decades, it has been encountered frequently in the northwestern plains zone (NWPZ) of the country, especially during cool and wet Rabi seasons (Rana et al., 2000). Due to the changing climatic conditions, SNB has been expanding into new niches. For example, it was first observed on emmer wheat (T. dicoccoides) in Turkey in 2017 and has been reported recently in Himachal Pradesh, India (Cat et al., 2018; Katoch et al., 2019).
Tan spot (TS), or yellow leaf spot, is a serious wheat disease affecting temperate and tropical wheat-growing regions. The fungal pathogen, a necrotroph that causes minor to severe spotting in wheat, was first described in 1823 (Hosford, 1982). The disease was subsequently reported in Europe, the USA, and Japan in the early 1900s (Wegulo, 2011). In India, reports on TS infection from northern plains and central regions were documented between 1934 and 1972 (Mitra, 1934; Misra and Singh, 1972), followed by a more recent one in 2007 (Singh, 2007). The development of the disease is encouraged by the favourable climate in South Asia, particularly in the Himalayan and eastern Gangetic plains, where low temperature and humidity favour prolonged leaf wetness. It has been reported that TS co-occurs with spot blotch in some parts of Nepal, causing yield losses of up to 20-30% (Duveiller et al., 2005). In Nepal, there is a higher incidence of foliar blight in the last week of January to mid-February when the temperature is still cool, and a few genotypes were identified as tolerant to the disease, like NL750, Milan, and Shanghai-7 (Dubin and Bimb, 1994; Duveiller et al., 2005; Gurung et al., 2012). So far, no systematic work has been done on SNB and TS in India except for reports on its occurrence. Moreover, no isolates have been deposited in the national fungal culture collection and type culture collections.
It is thought that P. nodorum obtains nutrients from dying plant tissue caused by secreted effectors. These effectors cause host hypersensitivity and result in programmed cell death (Friesen et al., 2007; Oliver et al., 2012). So far, eight effectors (SnToxA, SnTox1, SnTox2, SnTox3, SnTox4, SnTox5, SnTox6, and SnTox7) have been identified to date, along with nine major wheat sensitivity loci that correspond to them: Tsn1, Snn1, Snn2, Snn3-B1/Snn3-D1, Snn4, Snn5, Snn6, and Snn7, respectively (Friesen et al., 2007; Friesen and Faris, 2012; Tan et al., 2012; Phan et al., 2016).
Marker-trait associations (MTAs) of polygenetic traits in plants have been widely identified using genome-wide association studies (GWAS). Adhikari et al. (2011); Gurung et al. (2014); Liu et al. (2015), and Phan et al. (2018) have all discovered MTAs for SNB resistance at seedling resistance. In seedling experiments for resistance against a P. nodorum isolate lacking SnTox3, QTL was found on 2D, 3A, and 5B (Adhikari et al., 2011; Gurung et al., 2014). The 5B QTL’s association with the Tsn1 gene, which has been extensively studied (Liu et al., 2009; Gurung et al., 2014), supports the ability of GWAS to identify QTL for SNB resistance. TS resistance is both quantitatively and qualitatively inherited, where toxicity resistance genes and QTL have been identified. Tan spot resistance (Tsr) refers to the quality genes discovered through conidial inoculations, and “Tsc” and “Tsn” refer to genes for chlorosis and necrosis reactions against HST-containing cultures, respectively (McIntosh et al., 2008). In addition to the host susceptibility gene Tsn1, eight significant Tsr genes have so far been discovered, i.e., Tsrl, Tsr2, Tsr3, Tsr4, Tsr5, Tsr6, TsrHar, and TsrAri, being located on the chromosomes 2BS, 3AS, 3BL, 3DS, and 5BL (Kokhmetova et al., 2021; Lozano-Ramírez et al., 2022).
Though spot blotch remained a major concern in the Indian subcontinent among the foliar blotch diseases, it is important to do pre-emptive screening work on SNB and TS, considering the past reports and future outbreak risk under climate changing scenarios. In the present study, a panel of 296 bread wheat genotypes released in India over the past 20 years was screened at seedling stages for SNB and TS resistance. Further, the population was studied for the presence of the major toxin sensitivity gene Tsn1, and genome-wide association studies were performed for the genetic basis of the resistance.
2 Materials and methods
2.1 Plant material and genotyping
This study used a panel of 296 bread wheat genotypes from released varieties and advanced breeding lines developed primarily over the last twenty years from 25 research centres across India. DNA extraction and DArT sequencing of the genotypes were done per the protocol demonstrated by Li et al. (2016). The wheat accessions were sequenced with the DArTseq® technology at the Genetic Analysis Service for Agriculture (SAGA) at CIMMYT, Mexico. Markers with a minor allele frequency of less than 10% (2804 markers) or more than 30% missing data points (96 markers) were excluded from further analysis. A total of 9668 SNPs were finally used for the GWAS analysis, and their physical positions on the reference whole genome sequence (IWGSC: Chinese Spring RefSeq v1.0, International Wheat Genome Sequencing Consortium et al., 2018) were acquired from the database https://wheat-urgi.versailles.inra.fr/). The SNP markers were given names based on their chromosome location followed by clone ID, e.g., 5BL:3955588.
2.2 Disease screening for SNB and TS
The panel was evaluated for 3 years (2019, 2021, 2022) for SNB and 2 years for TS (2019, 2021) in greenhouse at the seedling stage. Mexican P. tritici-repentis (Ptr) isolate MexPtr1 and P. nodorum isolate MexSn4 were used for resistance screening against TS and SNB, respectively. Both isolates are ToxA producers based on inoculation experiments using differential genotypes, infiltration experiments, and PCR with the ToxA-specific marker (data not shown). The isolates were grown on V8-PDA media (Lamari and Bernier, 1989), and conidiospore concentrations for inoculation were adjusted to 4 × 103 spores mL−1 (MexPtr1) and 1 × 107 spores mL−1 (MexSn4) (Singh et al., 2007; Singh et al., 2016). The response of TS and SNB was tested on seedlings in a greenhouse at 22°C day and 18°C night temperatures with a 16-h photoperiod. Experiments were set up in a randomized complete block design with two replicates, with four plants grown in plastic containers as experimental units to derive mean values for further analysis. Erik and Glenlea were used as resistant and susceptible controls, respectively. Inoculations were performed when the second leaf was fully expanded at around 14 days after planting. The inoculum was applied to seedlings with a hand sprayer until runoff (about 0.5 mL inoculum per plant). The trays were transferred to a humid chamber (RH 100%, 20°C) once the leaves were dry to promote infection, and the plants were returned to the greenhouse bench after 24 hours. Both diseases were rated on a linear scale of 1–5 at seven days after inoculation (dpi) (Feng et al., 2004; Hu et al., 2019).
2.3 Linkage disequilibrium and population structure
All 9668 SNP markers were used to calculate a kinship matrix, clusters among individual genotypes, and a heat map using the traditional Van Raden (2008) equation using TASSEL v5 (http://www.maizegenetics.net, accessed on 25 Oct 2022). Linkage disequilibrium (LD) was analysed with R2 among SNP markers plotted against the physical distances in mega base pairs (Mb) across the 21 wheat chromosomes. The genotypic data was numerically transformed for population structure analysis utilizing XLSTAT (v. 2022.1). Structure 2.3.4 software was used to obtain the population structure (Pritchard et al., 2000). The admixture model was adjusted with a 100,000 burn-in period followed by 500,000 marker chain Monte Carlo (MCMC) iterations. The subpopulation test range was maintained at K1 - K5, with five iterations. The actual subpopulations were assessed using the DK approach (Earl and vonHoldt, 2012), which was verified using the Structure Harvester program (Web v0.6.94, Earl and vonHoldt, 2012) as per the method described by Evanno et al. (2005). The output summary calculated the standard deviation and average logarithm of the probability of the observed likelihood [LnP(D)]. The log-likelihood of the data was computed for each class (K = 1 to 5) to determine LnP(D) for each MCMC step. A neighbour-joining tree was created using TASSEL 5.0 (Bradbury et al., 2007) and visualized using the iTOL website (Letunic and Bork, 2021).
2.4 Genome-wide association analysis
For GWAS, three models (MLM, MLMM and FarmCPU) implemented in the R package GAPIT3 (Wang and Zhang, 2021) were used. The first was the mixed linear model (MLM=K+Q), which was based on the kinship matrix (K) and the principal component (PC) (Yu et al., 2006). The second method was the Multiple Loci Mixed Linear Model (MLMM), which employs forward-backward stepwise linear mixed-model regression to include associated markers as covariates. In contrast to MLM and MLMM, the third model FarmCPU (fixed and random, circulating probability unification); used all linked markers in a fixed-effect model and optimised the linked markers in a separate random-effect model, reducing false positives and false negatives and enabling quick computation. Population structure, principal components, and kinship were used as covariates in MLM and MLMM, and their QQ plots were compared. Our phenotypic data fit the FarmCPU model better, and QTL showed higher significance; therefore, only data analyzed by the FarmCPU model was chosen for additional examination. The Bonferroni correction (α=0.1) was used for an exploratory significance threshold to uncover putative QTL (Chan et al., 2010). QTL was recognized as robust when linked markers reached the strict -log10(p) criteria of 3.0 in at least being detected across-year by either threshold adopted in this investigation and consistent across the models. The Quantile-Quantile (QQ) plots were also examined to determine the point at which the observed p-values diverge from those predicted by the null hypothesis. SNP marker annotations were obtained using the databases http://www.cerealsdb.uk.net and https://triticeaetoolbox.org. The physical positions of markers on the reference whole genome sequence (IWGSC: Chinese spring RefSeq v1.0, International Wheat Genome Sequencing Consortium et al., 2018) were acquired from the database https://wheat-urgi.versailles.inra.fr/ and https://plants.ensembl.org/Triticum_aestivum/Info/Index. If two significant markers shared a 10-Mbp interval or had substantial LD (R2>0.8) with one another, they were regarded as belonging to the same QTL region.
2.5 KASP genotyping for Tsn1
Kompetitive Allele-Specific PCR (KASP) marker Ktsn1 specific for Tsn1 (Dreisigacker et al., 2016) was used for genotyping the panel. Genomic DNA was isolated from 21-day-old seedling leaves, DNA concentration was measured with a Nanodrop 8000 spectrophotometer (Thermo Scientific), and the final concentration of 10 ng/ml was obtained by diluting the DNA with sterile PCR-grade water. KASP genotyping was undertaken using the PACE master mix per the manufacturer’s guidelines (3CR biosciences). In brief, the assay mix for 25 samples was prepared as 3µl VIC primer, 7.5 µl common primer, and 14.5 µl of PCR grade H2O. Finally, for reaction assembly, individual reactions contained 5 µl 2X PACE, 0.138 µl assay mix, and 5 µl genomic DNA. The samples were run on a CFX96 Real-time PCR system (Bio-Rad). The PCR steps included 94°C for 15 min, 10 cycles of touch down program (94°C-0.20 min, 65 to 57°C-1 min), followed by 30 cycles of (94°C-0.20 min, 57°C-1 min). The results were analyzed in Bio-Rad CFX manager 3.1, and output was obtained as.csv files. The results obtained were compared against the linked SNP calls from the SNP array.
2.6 Haplotype analysis for the stable markers
For haplotype analysis, stable MTAs on chromosomes 1B, 2A, 5B, and 6B across experiments were selected. Eight markers that fitted the requirement that they were over the -log10(p) 3.0 threshold across the environments were chosen for the haplotype analysis: 5BS:1102120, 2AS:7487614, 2AS:1094287, 5BL:5324846, 5BL:3955588, 6B:1085698, 1B:1129298, and 5BL: Ktsn1. For TS, three markers, 5BL:5324846, 5BL:3955588, and 5BL : Ktsn1 were selected for the analysis. The corrected disease severities between haplotypes were compared using the Wilcoxon test implemented in R package ggpubr (Kassambara, 2020). The phenotypic data available for the seedling resistance to SNB in 2019, 2021, 2022, and TS in 2019 and 2021 was used for the analysis.
2.7 Stacking resistance allele and assessing interaction effect of Tsn1 and Snn3
Data for three years (2019, 2021, and 2022) were used for this analysis. The corrected mean disease index from BLUP was used for haplotype analysis. Marker trait associations for SNB and TS were chosen from various models (MLM, MLMM, FarmCPU) to investigate the effect of stacking resistance alleles (Table S1). Resistant alleles were determined by mean comparison of corrected disease severity between alleles based on the Wilcoxon test, using the R package ‘ggpubr’ (Kassambara, 2020). Wheat lines were classified according to the number of resistant alleles present. The t-test (p < 0.05) was implemented in the R package ‘multcomView’ (Graves et al., 2015) to differentiate the significant groups.
To assess the gene effect of Tsn1 and Snn3, two makers, 5BL:5324846 and 5BL:3955588, linked to Tsn1 and data obtained from KASP assay Ktsn1 was compared with SNP 5BS:1102120 linked with Snn3. The interaction was discovered by comparing the mean between alleles using the t-test in the R package ‘ggpubr’ (Kassambara, 2020).
2.8 Statistical analysis
The corrected disease index, genotypic and phenotypic variance, and heritability estimates were obtained from META-R software (Alvarado et al., 2020). Additionally, mean comparisons with the Wilcoxon test, t-test, and data visualization were executed using packages `ggpubr`, ‘dplyr’, and ‘multcomView’ using R software (Version 4.2.1, R Core Team, 2022).
3 Results
3.1 Disease evaluations
The population displayed significant phenotypic variation for SNB and TS resistance with a skewed distribution towards the lower disease in all experiments (Figure S1). A high proportion of genotypes exhibited high resistance to both diseases (See Table S1 for the disease scores) and the top 20 entries resistant to both diseases are presented in Table 1, including PBW 277, PBW 771, and HD 3319.
Table 1 List of wheat genotypes showing good resistance to Septoria nodorum blotch (SNB) and tan spot (TS), and their Tsn1/tsn1 allele status.
Analysis of variance revealed significant effects on genotype and genotype-by-environment for both diseases (Table 2). SNB showed higher overall heritability (0.87) than TS (0.75), with heritability in individual experiments ranging from 0.90 to 0.96.
Table 2 Basic statistics for Septoria nodorum blotch (SNB) and tan spot (TS) and analysis of variance across the environments.
3.2 Population structure and linkage disequilibrium analysis
Among the 9668 SNP markers selected for GWAS in the 296 bread wheat genotypes, 21.78% were from the A genome, 28.63% from the B genome, 12.03% from the D genome and 21.33% from unknown chromosomes. Population structure analysis based on the K means cluster approach divided the population into 4 major clusters that were depicted as a neighbour-joining tree (Figure 1). The kinship analysis also distinguished the population into 4 major groups and presented it as a kinship matrix-based heat map (Figure 2). The average extent of LD considered physical distance taken for the decay of R2 to reach a critical value of 0.10 across the genome, was approximately 10 Mb (Figures 2, S2, S3).
Figure 1 Population structures of 296 genotypes revealed by STRUCTURE 2.3.4 and neighbour-joining tree.
Figure 2 (A) 3D plots of principal components (PC). (B) The kinship matrix used in GWAS is visualized through a heat map. (C) LD matrix obtained from TASSEL 5.0.
3.3 Genome-wide association study for SNB and TS
The exploratory -log10(p) threshold for the panel ranged from 2.92 to 11.53 for SNB and 2.96 to 15.86 for TS (Tables S1, S2). In total, 49 marker-trait associations (MTAs) were detected on various chromosomes for SNB (Table S1; Figures 3, S4). Six QTL on chromosomes 1BL, 2AS, 5BL, and 6BL were particularly significant for SNB over all three years and were considered robust QTLs. The QTL Q.CIM.snb.1BL on chromosome 1BL was identified with FarmCPU, and its peak marker 1BL:1129298 (p = 2.86E-04 to 8.97E-06) was located at 450.50Mbp on 1BL (Table 3). Two QTLs on 2AS, Q.CIM.snb.2AS1 and Q.CIM.snb.2AS.2, were detected in all models with their peak markers being 2AS:7487614 (59.43Mbp) and 2AS:1094287 (88.18Mbp), respectively. Another robust QTL, Q.CIM.snb.5BS, was identified at 5BS:1102120 (p=2.31E-05 to 5.99E-12) and was associated with Snn3-B1. Two Tsn1-associated markers 5BL:5324846 and 5BL:3955588 were detected across years and models on chromosome 5B at 566.04 and 588.4 Mbp, respectively (Table 3). The QTL Q.CIM.snb.6BL (p= 4.16E-04 to 1.61E-09) was detected using farmCPU, on the long arm of the chromosome 6B with peak marker 6BL:1085698 (669.13 Mbp).
Figure 3 Manhattan plots of marker-trait associations detected by the FarmCPU model using 9668 single nucleotide polymorphisms in 296 Indian spring wheat genotypes evaluated for the seedling resistance to Septoria nodorum blotch (SNB) during years 2019, 2021, 2022 and tan spot (TS) during years 2019, 2021. The horizontal line represents the significant threshold –log10(P) = 5.00.
Table 3 Markers significantly associated with Septoria nodorum blotch (SNB) and tan spot (TS) resistance through genome-wide association mapping using a mixed linear model (MLM), multi-locus mixed model (MLMM) and fixed and random model circulating probability unification (FarmCPU).
For TS, 79 MTA (-log10(p) >3) were detected on various chromosomes (Table S2; Figure 3). Two Tsn1-associated markers 5BL:5324846 and 5BL:3955588 were significantly detected across years. The SNP 5BL:5324846 (p=3.51E-05 to 5.41E-06) was detected across the three models (MLM, MLMM, FarmCPU). While the SNP 5BL:3955588 (p=1.33E-05 to 3.63E-13) was detected with a negative effect across the years (2019, 2021) from models MLM and MLMM only (Table 3).
Out of the 49 MTA identified for SNB resistance, 15 were selected to analyse their association with predicted genes and the corresponding biological functions. The criteria for shortlisting MTAs was that they were detected in at least two models or two experiments (See Table 4). Two Tsn1-associated markers, 5BL:5324846 and 5BL:3955588, were found to be associated with the transcripts TraesCS5B02G409000, TraesCS5B02G387500, and TraesCS5B02G387000. These transcripts are known to possess biological function pathogen stress response (PADRE) domain, Serine/threonine-protein kinase and Ubiquitin-like proteins respectively. Another marker, 2AS:7487614, has been found to be associated with two genes, namely TraesCS2A02G107100 and TraesCS2A02G107200. These genes are known to have biological roles as a potassium transporter and a plant peroxidase, respectively. As for other MTAs, several noteworthy biological functions have been identified, including cysteine-type endopeptidase inhibitor activity (6B:1085698), NAD(P)-binding domain (4A:1082366), Cytochrome c oxidase subunit 5c, and Serine/threonine-protein kinase, active site (both for 3B:2251334). These functions are believed to play an active role in disease resistance (Table 4).
Table 4 SNPs associated with SNB and TS resistance and possible function elucidated based on the gene annotation using wheat reference sequence (RefSeq V1.0).
3.4 The significant role of Tsn1 and Snn3 on SNB resistance
The KASP assay revealed that 83 genotypes (25.85%) carry the susceptible allele Tsn1, whose SNB indices ranged mostly between 2.4 to 4.2 (Table S3). The allelic combinations at Tsn1 and Snn3 were tested using the results obtained from the Ktsn1 marker and two Tsn1 linked markers, 5B:5324846 and 5B:3955588, along with Snn3-B1 linked marker 5BS:1102120. The comparison revealed that the allelic combination of recessive alleles tsn1/snn3 imparts best resistance to SNB, followed by the tsn1/Snn3 combination, whereas the allele combinations Tsn1/snn3 and Tsn1/Snn3 showed similar levels of resistance, indicating the higher phenotypic effects of tsn1(Figure 4).
Figure 4 Boxplots showing gene interaction effect for Tsn1 and Snn3 in response to Septoria nodorum blotch (SNB). *** on boxplots indicates a significant difference in the mean disease index between groups by Wilcoxon test (p < 0.05).
3.5 Haplotype analysis and staking of R alleles
There were significant differences in the corrected SNB index between resistant and susceptible haplotypes. When compared to the resistance haplotypes “AA/TT” for tsn1 (mean SNB index 1.69), the susceptibility haplotype “GG” for Tsn1 had a higher disease index (3.4) (Figure 5). Furthermore, the haplotypes were tested for the four significant SNB QTLs, Q.CIM.snb.1BL, Q.CIM.snb.2AS.1, Q.CIM.snb.2AS.2, and Q.CIM.snb.6BL (Figure S5). Similarly, there was a significant difference between susceptible and resistant haplotypes for TS index (3.2 for the susceptibility haplotypes and 1.57 for the resistance haplotypes; Figure 5).
Figure 5 Haplotype analysis of the markers (A) 5BL:5324846 (B) 5BL:3955588 (C) KASP marker Ktsn1 indicated their significant association with Septoria nodorum blotch (SNB) and tans spot (TS) resistance determined by the Wilcoxon test.
The stacking of resistance alleles for SNB and TS exhibited a clear trend that the number of resistance alleles is positively correlated with disease resistance (Figure 6). For example, cultivars DBW 242 and DBW 246, each possessing 35 resistant alleles, exhibited very high SNB resistance, whereas HD 2932 and Kharchiya 65 having the least number of resistance alleles were highly susceptible. A similar trend was observed for TS resistance.
Figure 6 Boxplots showing effects of stacking resistant alleles for (A) Septoria nodorum blotch (SNB) and (B) Tan spot (TS). *** on boxplots indicates a significant difference in the mean disease index between groups by Wilcoxon test (p < 0.05).
4 Discussion
Utilizing resistant cultivars coupled with management approaches contribute to environmentally and economically sustainable wheat production. The discovery of host-specific gene-for-gene interactions determining the P. nodorum-wheat pathosystem and residual genetic resistance (Liu et al., 2004) offers immediate opportunities to further investigate host-genetic resistance in wheat breeding (Ruud and Lillemo, 2018; Peters Haugrud et al., 2022). Evaluations of the Indian wheat cultivars reveal at least 113 genotypes resistant for the TS, 174 genotypes for SNB, and 118 for both the SNB and TS. Out of the top 20 highly resistant genotypes, WH1256, GW 519 and DBW 285 carry the susceptibility allele Tsn1, implying that they must have other resistance genes/QTL to counteract the negative effects of Tsn1. By checking the pedigree information, it was found that WH1256 and DBW 285 have Chinese progenitors and might have inherited unknown resistance genes.
We found seedling resistance QTL Q.CIM.snb.2AS for SNB close to previously reported QTL QSnb.nmbu-2AS (4-24 Mbp), which may be associated with two adult-plant leaf blotch QTLs at 2-20 and 15-16 Mbp on the reference genome (Lin et al., 2022). Additionally, this region was also mapped near the seedling resistance QTLs Snb.niab-2A.1 (0.8-2.4 Mbp) and Qsnb.cur-2AS2 (2.3-3.8 Mbp) reported by Lin et al. (2020) and Phan et al. (2016), respectively. Snn3-B1 has been associated with the short arm of chromosome 5B in previous research (Phan et al., 2016; Ruud et al., 2017; Ruud et al., 2019; Lin et al., 2020). Here, genetic analysis in the Indian wheat cultivars identified one marker 5BS:1102120 tightly linked to Snn3-B1 as representing the most significant genetic determinant of SnTox3 sensitivity in Indian germplasm.
Numerous shared QTLs between SNB and TS have been found in an increasing number of publications on the QTL mapping of both diseases (Phuke et al., 2020). A well-known example is Tsn1 which confers sensitivity to both SnToxA and PtrToxA (Friesen et al., 2006). The current study has proven that Tsn1 plays a major role in the susceptibility of the Indian panel to SNB and TS. TS resistance in seedlings and adult plants has also been discovered to be significantly influenced by the P. nodorum resistance/sensitivity QTL Qsnb.cur-2AS.1 (Phan et al., 2016), discovered at the seedling and adult plant stage (Phan et al., 2016; Shankar et al., 2017). This phenomenon might suggest similar susceptibility/resistance mechanisms between the two diseases (Lin and Lillemo, 2021; Lin et al., 2022). It would be intriguing to learn if they share any additional effectors. To breed wheat varieties resistant to both diseases, wheat breeders can focus on mutual interactions, particularly for those QTL with relatively large effects and at both the seedling and adult stages.
The primary factor causing leaf necrosis is the HST ToxA, produced by the pathogens P. nodorum and P. tritici-repentis (produced by races 1, 2, 7, and 8). These pathogens are highly susceptible due to the sensitive gene Tsn1. The most prevalent Ptr race, Race 1, can also be found in South Asia and Mexico (Duveiller and Singh, 2009). However, P. nodorum from South Asia lacks information on the racial composition, variability, and toxins produced. In the last 25 years, South Asian breeding programs have incorporated new knowledge about resistance to leaf blight disease complex. Researchers and breeders worldwide could ascertain the relationship between effector sensitivity and cultivar susceptibility using expressed ToxA since 2005 and Tox3 since 2011. The ToxA sequence in various B. sorokiniana isolates from India has been analyzed recently, and the results show that the gene is under positive selection (Navathe et al., 2020). The investigations on the pathogenicity factors, variability, and toxin profiling of P. nodorum and P. tritici-repentis from South Asia is lacking, irrespective of their numerous reports from India and Nepal. The horizontal transfer of ToxA between these three pathogens underlines that unless Tsn1/Snn3 is selectively bred out of widely planted wheat germplasms, it is likely that ToxA will continuously evolve into forms that are more effective in inducing host cell death. Other cutting-edge technologies will accelerate the discovery and functional characterization of effector resistance genes in the coming years and offer effective methods for utilizing these in breeding programs.
India and Nepal have reported SNB, TS, and spot blotch occurrences individually or in the complex. Moreover, climate change has warranted cross-continent jump in wheat diseases and we know very little about the interactions between these three diseases. Further, this subject is especially fascinating given that P. nodorum, P. tritici-repentis and B. sorokiniana share effectors. Because a new cultivar must be fully assessed at adult plant stage in field conditions, breeding for SNB/TS/SB resistance has always been difficult. The difficulties are made worse by inoculation using a representative group of isolates. Making significant annual isolate collections, especially from the cultivar that is currently most resistant, is one clear recommendation that has come out of recent studies. These novel isolates can be examined phenotypically for novel effectors and virulence traits and genotypically to follow specific chromosomal regions. Any new effector can be expressed, and their contribution to virulence can be evaluated. The main benefit of isolate collections is that they make it possible to rationally choose the smallest set that accurately captures all the pathogen’s phenotypic diversity, against which resistance should be sought.
At present, complete genetic resistance to SNB and TS has not been discovered. Therefore, a multifaceted strategy based on agronomic practices, disease surveillance, and genetic resistance will be required for the effective management of SNB and TS. Nevertheless, genetic resistance will continue to play a major role in managing these diseases. It was interesting to see that some varieties like PBW 771, DBW 277, and HD 3319 displayed high levels of resistance to both diseases. Hence, they must be deployed in the ongoing breeding programs for further enhancement of resistance for SNB and TS. Using field resistance and knowledge about new QTLs will certainly help breeders find much better resistance for these diseases facilitating higher production in the farmer’s fields.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. The genotypic data is available at Dataverse CIMMYT data repository and can be accessed at https://hdl.handle.net/11529/10548934. Further inquiries can be directed to the corresponding author/s.
Author contributions
SN contributed to manuscript writing and data analysis, XH contributed to conceptualization, disease scoring and data analysis, and UK, MK, and MP helped in data analysis and manuscript writing. PKS and AJ contributed to the conceptualization, project monitoring, fund acquisition and revising. GPS and GS contributed to germplasm resources and fund acquisition. All authors contributed to the article and approved the submitted version.
Funding
Financial support from the Indian Council of Agriculture Research (ICAR), India, the Bill and Melinda Gates Foundation through project Accelerating Genetic Gain (AGG) in Maize and Wheat Project Grant INV-003439, USAID, and One CGIAR Initiatives for conducting this research is gratefully acknowledged.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1223959/full#supplementary-material
References
Adhikari, T. B., Jackson, E. W., Gurung, S., Hansen, J. M., Bonman, J. M. (2011). Association mapping of quantitative resistance to Phaeosphaeria nodorum in spring wheat landraces from the USDA National Small Grains Collection. Phytopathology 101, 1301–1310. doi: 10.1094/phyto-03-11-0076
Alvarado, G., Rodríguez, F. M., Pacheco, A., Burgueño, J., Crossa, J., Vargas, M., et al. (2020). META-R: a software to analyze data from multi-environment plant breeding trials. Crop J. 8, 745–756. doi: 10.1016/j.cj.2020.03.010
Bearchell, S. J., Fraaije, B. A., Shaw, M. W., Fitt, B. D. L. (2005). Wheat archive links long-term fungal pathogen population dynamics to air pollution. Proc. Natl. Acad. Sci. U.S.A. 102, 5438–5442. doi: 10.1073/pnas.0501596102
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Cat, A., Tekin, M., Akar, T., Catal, M. (2018). First report of Stagonospora nodorum blotch caused by Parastagonospora nodorum on emmer wheat (Triticum dicoccum Schrank) in Turkey. J. Plant Pathol. 17, 42161. doi: 10.1007/s42161-018-00205-2
Chan, E. K., Rowe, H. C., Kliebenstein, D. J. (2010). Understanding the evolution of defense metabolites in Arabidopsis thaliana using genome-wide association mapping. Genetics 185, 991–1007. doi: 10.1534/genetics.109.108522
Dreisigacker, S., Sehgal, D., Reyes Jaimez, A. E., Luna Garrido, B., Muñoz Zavala, S., Núñez Ríos, C., et al. (2016). CIMMYT wheat molecular genetics: laboratory protocols and applications to wheat breeding (Mexico, D.F: CIMMYT).
Dubin, H. J., Bimb, H. P. (1994). Studies of soilborne diseases and foliar blights of wheat at the national wheat research experiment station, Bhairahawa, Nepal. Wheat Special Report No. 36 (Mexico DF, Mexico: CIMMYT).
Duveiller, E., Kandel, Y. R., Sharma, R. C., Shrestha, S. M. (2005). Epidemiology of foliar blights (spot blotch and tan spot) of wheat in the plains bordering the Himalayas. Phytopathology 95, 248–256. doi: 10.1094/phyto-95-0248
Duveiller, E., Singh, P. K. (2009). Improving Wheat Resistance to Tan Spot caused by Pyrenophora tritici-repentis, Newsletter (Mexico: CIMMYT). Available at: http://hdl.handle.net/10883/1293.
Earl, D. A., vonHoldt, B. M. (2012). Structure Harvester: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7
Evanno, G., Regnaut, S., Goudet, J. (2005). Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294x.2005.02553.x
Feng, J., Ma, H., Hughes, G. R. (2004). Genetics of resistance to Stagonospora nodorum blotch of hexaploid wheat. Crop Sci. 44, 2043–2048. doi: 10.2135/cropsci2004.2043
Friesen, T. L., Faris, J. D. (2012). Characterization of plant-fungal interactions involving necrotrophic effector–producing plant pathogens. Methods Mol. Biol. 835, 191–207. doi: 10.1007/978-1-61779-501-5_12
Friesen, T. L., Meinhardt, S. W., Faris, J. D. (2007). The Stagonospora nodorum–wheat pathosystem involves multiple proteinaceous host-selective toxins and corresponding host sensitivity genes that interact in an inverse gene-for-gene manner. Plant J. 51, 681–692. doi: 10.1111/j.1365-313x.2007.03166.x
Friesen, T. L., Stukenbrock, E. H., Liu, Z. H., Meinhardt, S., Ling, H., Faris, J. D. (2006). Emergence of a new disease as a result of interspecific virulence gene transfer. Nat. Genet. 38, 953–956. doi: 10.1038/ng1839
Graves, S., Piepho, H., Selzer, L. (2015) Dorai-Raj, S. multcompView: visualizations of paired comparisons version 01–7. Available at: https://cranr-project.org/web/packages/multcompView/index.html.
Gurung, S., Mamidi, S., Bonman, J. M., Xiong, M., Brown-Guedira, G., Adhikari, T. B. (2014). Genome-wide association study reveals novel quantitative trait loci associated with resistance to multiple leaf spot diseases of spring wheat. PloS One 9, e108179. doi: 10.1371/journal.pone.0108179
Gurung, S., Sharma, R. C., Duveiller, E., Shrestha, S. M. (2012). Comparative analyses of spot blotch and tan spot epidemics on wheat under optimum and late sowing period in South Asia. Eur. J. Plant Pathol. 134, 257–266. doi: 10.1007/s10658-012-9984-6
Hosford, R. M., Jr. (1982). Tan spot of wheat and related diseases workshop. Ed. Hosford, R. M., Jr (North Dakota State Univ., Fargo USA), 1–24.
Hu, W., He, X., Dreisigacker, S., Sansaloni, C. P., Juliana, P., Singh, P. K. (2019). A wheat chromosome 5AL region confers seedling resistance to both tan spot and Septoria nodorum blotch in two mapping populations. Crop J. 7, 809–818. doi: 10.1016/j.cj.2019.05.004
International Wheat Genome Sequencing Consortium, Appels, R., Eversole, K., Stein, N., Feuillet, C., Keller, B., Rogers, J., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361, eaar7191. doi: 10.1126/science.aar7191
Joshi, L. M., Saari, E. E., Wilcoxson, R. D. (1971). Epidemic of glume blotch on wheat and reaction of wheat varieties to. Septoria nodorum. Indian Phytopathol. 24, 413–415.
Kassambara, A. (2020). ggpubr: ‘ggplot2’ based publication ready plots. R package version 0.4.0. Available at: https://CRAN.R-project.org/package=ggpubr.
Katoch, S., Rana, S. K., Sharma, P. N. (2019). Application of PCR based diagnostics in the exploration of Parastagonospora nodorum prevalence in wheat growing regions of Himachal Pradesh. J. Plant Biochem. Biotechnol. 28, 169–175. doi: 10.1007/s13562-018-0481-7
Kokhmetova, A., Sehgal, D., Ali, S., Atishova, M., Kumarbayeva, M., Leonova, I., et al. (2021). Genome-wide association study of tan spot resistance in a hexaploid wheat collection from Kazakhstan. Front. Genet. 11. doi: 10.3389/fgene.2020.581214
Lamari, L., Bernier, C. C. (1989). Evaluation of wheat lines and cultivars to tan spot (Pyrenophora tritici-repentis) based on lesion type. Can. J. Plant Pathol. 11, 49–56. doi: 10.1080/07060668909501146
Letunic, I., Bork, P. (2021). Interactive Tree of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49 (W1), W293–W296. doi: 10.1093/nar/gkab301
Li, H., Vikram, P., Singh, R. P., Kilian, A., Carling, J., Song, J., et al. (2016). A high-density GBS map of bread wheat and its application for dissecting complex disease resistance traits. BMC Genom. 16, 216. doi: 10.1186/s12864-015-1424-5
Lin, M., Corsi, B., Ficke, A., Tan, K. C., Cockram, J., Lillemo, M. (2020). Genetic mapping using a wheat multi-founder population reveals a locus on chromosome 2A controlling resistance to both leaf and glume blotch caused by the necrotrophic fungal pathogen Parastagonospora nodorum. Theor. Appl. Genet. 133, 785–808. doi: 10.1007/s00122-019-03507-w
Lin, M., Ficke, A., Dieseth, J. A., Lillemo, M. (2022). Genome-wide association mapping of Septoria nodorum blotch resistance in Nordic winter and spring wheat collections. Theor. Appl. Genet. 135, 4169–4182. doi: 10.1007/s00122-022-04210-z
Lin, M., Lillemo, M. (2021). “Advances in genetic mapping of Septoria nodorum blotch resistance in wheat and applications in resistance breeding,” in Achieving durable disease resistance in cereals. Ed. Oliver, R. P. (Burleigh Dodds Science Publishing, Cambridge, UK), 393–433. doi: 10.1201/9781003180715-15
Liu, Z., El-Basyoni, I., Kariyawasam, G., Zhang, G., Fritz, A., Hansen, J. (2015). Evaluation and association mapping of resistance to tan spot and Stagonospora nodorum blotch in adapted winter wheat germplasm. Plant Dis. 99, 1333–1341. doi: 10.1094/pdis-11-14-1131-re
Liu, Z. H., Faris, J. D., Meinhardt, S. W., Ali, S., Rasmussen, J. B., Friesen, T. L. (2004). Genetic and physical mapping of a gene conditioning sensitivity in wheat to a partially purified host-selective toxin produced by Stagonospora nodorum. Phytopathology 94, 1056–1060. doi: 10.1094/phyto.2004.94.10.1056
Liu, Z. H., Faris, J. D., Oliver, R. P., Tan, K. C., Solomon, P. S., McDonald, M. C. (2009). SnTox3 acts in effector triggered susceptibility to induce disease on wheat carrying the Snn3 gene. PloS Pathog. 5, e1000581. doi: 10.1371/journal.ppat.1000581
Lozano-Ramírez, N., Dreisigacker, S., Sansaloni, C. P., He, X., Islas, S. S., Pérez-Rodríguez, P., et al. (2022). Genome-wide association study for resistance to tan spot in synthetic hexaploid wheat. Plants 11, 433. doi: 10.3390/plants11030433
McIntosh, R. A., Yamazaki, Y., Dubcovsky, J., Rogers, J., Morris, C., Somers, D. J., et al. (2008). Catalogue of gene symbols for wheat (MacGene). Available at: http://wheat.pw.usda.gov/GG2/Triticum/wgc/2008/.
Misra, A. P., Singh, R. A. (1972). Pathogenic differences among three isolates of Helminthosporium tritici - repentis and the performance of wheat varieties against them. Indian Phytopathol. 25, 350–353.
Mitra, M. (1934). A leaf spot disease of wheat caused by Helminthosporium tritici-repentis Died. Indian J. Agric. Sci. 4, 692–700.
Navathe, S., Yadav, P. S., Chand, R., Mishra, V. K., Vasistha, N. K., Meher, P. K., et al. (2020). ToxA–Tsn1 interaction for spot blotch susceptibility in Indian wheat: An example of inverse gene-for-gene relationship. Plant Dis. 104, 71–81. doi: 10.1094/pdis-05-19-1066-re
Oliver, R. P., Friesen, T. L., Faris, J. D., Solomon, P. S. (2012). Stagonospora nodorum: from pathology to genomics and host resistance. Annu. Rev. Phytopathol. 50, 23–43. doi: 10.1146/annurev-phyto-081211-173019
Peters Haugrud, A. R., Zhang, Z., Friesen, T. L., Faris, J. D. (2022). Genetics of resistance to Septoria nodorum blotch in wheat. Theor. Appl. Genet. 135, 3685–3707. doi: 10.1007/s00122-022-04036-9
Phan, H. T., Rybak, K., Bertazzoni, S., Furuki, E., Dinglasan, E., Hickey, L. T., et al. (2018). Novel sources of resistance to Septoria nodorum blotch in the Vavilov wheat collection identified by genome-wide association studies. Theor. Appl. Genet. 131, 1223–1238. doi: 10.1007/s00122-018-3073-y
Phan, H. T., Rybak, K., Furuki, E., Breen, S., Solomon, P. S., Oliver, R. P. (2016). Differential effector gene expression underpins epistasis in a plant fungal disease. Plant J. 87, 343–354. doi: 10.1111/tpj.13203
Phuke, R. M., He, X., Juliana, P., Bishnoi, S. K., Singh, G. P., Kabir, M. R., et al. (2020). Association mapping of seedling resistance to tan spot (Pyrenophora tritici-repentis race 1) in CIMMYT and South Asian wheat germplasm. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.01309
Pritchard, J. K., Stephens, M., Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945
Rana, S. K., Paul, Y. S., Singh, D. (2000). Occurrence of glume blotch of wheat in Himachal Pradesh. J. Mycol Plant Pathol. 30, 427–428.
R Core Team. (2022). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing). Available at: https://www.R-project.org/.
Ruud, A. K., Dieseth, J. A., Ficke, A., Furuki, E., Phan, H. T., Oliver, R. P., et al. (2019). Genome-wide association mapping of resistance to Septoria nodorum leaf blotch in a Nordic spring wheat collection. Plant Genome 12, 180105. doi: 10.3835/plantgenome2018.12.0105
Ruud, A. K., Lillemo, M. (2018). “Diseases affecting wheat: Septoria nodorum blotch,” in Integrated disease management of wheat and barley. Ed. Oliver, R. P. (Cambridge, UK: Burleigh Dodds Science Publishing Ltd.). doi: 10.19103/as.2018.0039.06
Ruud, A. K., Windju, S., Belova, T., Friesen, T. L., Lillemo, M. (2017). Mapping of SnTox3–Snn3 as a major determinant of field susceptibility to Septoria nodorum leaf blotch in the SHA3/CBRD Naxos population. Theor. Appl. Genet. 130, 1361–1374. doi: 10.1007/s00122-017-2893-5
Shankar, M., Jorgensen, D., Taylor, J., Chalmers, K. J., Fox, R., Hollaway, G. J., et al. (2017). Loci on chromosomes 1A and 2A affect resistance to tan (yellow) spot in wheat populations not segregating for tsn1. Theor. Appl. Genet. 130, 2637–2654. doi: 10.1007/s00122-017-2981-6
Shaw, M. W., Bearchell, S. J., Fitt, B. D. L., Fraaije, B. A. (2008). Long term relationships between environment and abundance in wheat of Phaeosphaeria nodorum and Mycosphaerella graminicola. New Phytol. 177, 229–238. doi: 10.1111/j.1469-8137.2007.02236.x
Singh, D. P. (2007). First report of tan spot of wheat caused by Pyrenophora tritici-repentis in the Northern Hills and Northwestern Plains Zones of India. Plant Dis. 91, 460. doi: 10.1094/pdis-91-4-0460c
Singh, P. K., Crossa, J., Duveiller, E., Singh, R. P., Djurle, A. (2016). Association mapping for resistance to tan spot induced by Pyrenophora tritici-repentis race 1 in CIMMYTs historical bread wheat set. Euphytica 207, 515–525. doi: 10.1007/s10681-015-1528-7
Singh, P. K., Mergoum, M., Adhikari, T. B., Kianian, S. F., Elias, E. M. (2007). Chromosomal location of genes for seedling resistance to tan spot and Stagonospora nodorum blotch in tetraploid wheat. Euphytica 155, 27–34. doi: 10.1007/s10681-006-9297-y
Solomon, P., Lowe, R. G. T., Tan, K.-C., Waters, O. D. C., Oliver, R. P. (2006). Stagonospora nodorum: cause of Stagonospora nodorum blotch of wheat. Mol. Plant Pathol. 7, 147–156. doi: 10.1111/j.1364-3703.2006.00326.x
Tan, K. C., Ferguson-Hunt, M., Rybak, K., Waters, O. D., Stanley, W. A., Bond, C. S., et al. (2012). Quantitative variation in effector activity of ToxA isoforms from Stagonospora nodorum and Pyrenophora tritici-repentis. MPMI 25, 515–522. doi: 10.1094/mpmi-10-11-0273
Van Raden, P. M. (2008). Efficient methods to compute genomic predictions. J. Dairy Sci. 91, 4414–4423. doi: 10.3168/jds.2007-0980
Wang, J., Zhang, Z. (2021). GAPIT Version 3: Boosting power and accuracy for genomic association and prediction. Genomics Proteomics. Bioinformatics 19, 629–640. doi: 10.1016/j.gpb.2021.08.005
Wegulo, S. N. (2011). Tan spot of cereals (The Plant Health Instructor, The American Phytopathological Society, St. Paul, MN, USA). doi: 10.1094/PHI-I-2011-0426-01
Keywords: Parastagonospora nodorum, Pyrenophora tritici-repentis, quantitative inheritance, Triticum aestivum, Tsn1
Citation: Navathe S, He X, Kamble U, Kumar M, Patial M, Singh G, Singh GP, Joshi AK and Singh PK (2023) Assessment of Indian wheat germplasm for Septoria nodorum blotch and tan spot reveals new QTLs conferring resistance along with recessive alleles of Tsn1 and Snn3. Front. Plant Sci. 14:1223959. doi: 10.3389/fpls.2023.1223959
Received: 16 May 2023; Accepted: 20 September 2023;
Published: 10 October 2023.
Edited by:
Runsheng Ren, Jiangsu Academy of Agricultural Sciences (JAAS), ChinaReviewed by:
Jaspal Kaur, Punjab Agricultural University, IndiaAamir W. Khan, University of Missouri, United States
Satinder Kaur, Punjab Agricultural University, India
Copyright © 2023 Navathe, He, Kamble, Kumar, Patial, Singh, Singh, Joshi and Singh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pawan Kumar Singh, pk.singh@cgiar.org