Genome-wide association analyses of agronomic traits and Striga hermonthica resistance in pearl millet

Pearl millet (Pennisetum glaucum [L.] R. Br.) is a nutrient-dense, relatively drought-tolerant cereal crop cultivated in dry regions worldwide. The crop is under-researched, and its grain yield is low (< 0.8 tons ha−1) and stagnant in the major production regions, including Burkina Faso. The low productivity of pearl millet is mainly attributable to a lack of improved varieties, Striga hermonthica [Sh] infestation, downy mildew infection, and recurrent heat and drought stress. Developing high-yielding and Striga-resistant pearl millet varieties that satisfy the farmers’ and market needs requires the identification of yield-promoting genes linked to economic traits to facilitate marker-assisted selection and gene pyramiding. The objective of this study was to undertake genome-wide association analyses of agronomic traits and Sh resistance among 150 pearl millet genotypes to identify genetic markers for marker-assisted breeding and trait introgression. The pearl millet genotypes were phenotyped in Sh hotspot fields and screen house conditions. Twenty-nine million single nucleotide polymorphisms (SNPs) initially generated from 345 pearl millet genotypes were filtered, and 256 K SNPs were selected and used in the present study. Phenotypic data were collected on days to flowering, plant height, number of tillers, panicle length, panicle weight, thousand-grain weight, grain weight, number of emerged Striga and area under the Striga number progress curve (ASNPC). Agronomic and Sh parameters were subjected to combined analysis of variance, while genome-wide association analysis was performed on phenotypic and SNPs data. Significant differences (P < 0.001) were detected among the assessed pearl millet genotypes for Sh parameters and agronomic traits. Further, there were significant genotype by Sh interaction for the number of Sh and ASNPC. Twenty-eight SNPs were significantly associated with a low number of emerged Sh located on chromosomes 1, 2, 3, 4, 6, and 7. Four SNPs were associated with days-to-50%-flowering on chromosomes 3, 5, 6, and 7, while five were associated with panicle length on chromosomes 2, 3, and 4. Seven SNPs were linked to thousand-grain weight on chromosomes 2, 3, and 6. The putative SNP markers associated with a low number of emerged Sh and agronomic traits in the assessed genotypes are valuable genomic resources for accelerated breeding and variety deployment of pearl millet with Sh resistance and farmer- and market-preferred agronomic traits.

in the farmers' fields is attributable to various biotic and abiotic constraints, including the use of low-yielding landraces, Striga hermonthica (Sh) infestation, bird damage, insect pests, diseases, heat and drought stresses 3,4 .
Striga hermonthica (Del.)Benth. is the most significant biotic constraint to pearl millet production and productivity in Burkina Faso and yield losses vary between 40 and 55% 5,6 .The parasite infests several other major cereal crops, including rice (Oryza glaberrima Steudel and O. Sativa L.), maize (Zea mays L.), sorghum (Sorghum bicolour [L.] Moench), and fonio (Digitaria exilis [Kippist] Stapf) 7,8 .Farmers often abandon Sh-infested fields and switch from pearl millet to other non-host crops, reducing the crop's overall production and economic value 9 .Striga is a highly prolific parasite in Burkina Faso attributed to the host crop being mostly grown in semi-arid parts of the Sahelian and Sudano-Sahelian zones, which are dominated by poor soil fertility, low and erratic rainfall, and high temperatures that favour germination, growth and spread of the weed 3 .
Striga control is difficult because each parasitic plant can quickly disperse and deposit thousands of seeds into the soil seedbank.Furthermore, Striga seeds can remain viable in the soils for more than 14 years 10 .Striga hermonthica is a major threat to food security, exacerbating hunger and poverty in many African countries 11,12 .Monetary losses ranging from 117 to 200 billion US$ is incurred annually due to crop damage, and increases by US $30 million annually 13 .
Several Striga control strategies are recommended, including hand weeding, mulching crop fields with biomass of the shea tree (Vitellaria paradoxa C.F. Gaertn.) as a bio-control agent, optimal fertilizer application, and soil moisture management 7,14 .These strategies improve the soil properties, promote crop growth and development, and retard germination and growth of Striga 15 .Herbicides are less effective in controlling the effect of the parasite after emergence, and they are unaffordable for smallholder farmers.
The use of Striga-resistant pearl millet varieties is the most sustainable and environment-friendly management option for smallholder farmers in semi-arid regions.Resistant cultivars support fewer Striga plants and yield higher 16,17 .However, with the paucity of locally adapted and Sh-resistant donor sources, breeding for Striga resistance in pearl millet is still challenging compared to other cereals [18][19][20][21] .In the past decade, intensive research on the interaction of Striga with the host at the molecular level has opened opportunities to develop new management strategies 22 .For instance, 154 candidate genes associated with Sh resistance traits were identified in maize 23 .Adewale et al. 24 reported 13 associated markers with the Sh resistance trait in early maturing tropical white maize inbred lines.
Genome-wide association studies (GWAS) has been used in pearl millet for the identification of putative genes related to flowering time 25 , iron, zinc and protein content 26 , downy mildew resistance 27 and Sh resistance 28 .Also, GWAS has been used in finger millet for the identification of genes associated with Striga resistance 29 and grain nutritional contents 30 and for genetic diversity analysis 31 .GWAS is a valuable genomic tool to identify quantitative trait loci (QTLs) linked to Striga resistance for marker-assisted selection.GWAS results depend on the genetic marker and its density, genetic composition and diversity of the test populations.
Genetic markers are landmarks on chromosomes that help pinpoint the location of genes of interest.They can be detected through morphological and molecular markers 9 .Genetic markers such as GRMZM2G077208, GRMZM2G164502, GRMZM2G018508, and GRMZM2G171986, located on chromosomes 3, 5, 7, and 9 were reportedly significantly associated with Sh count in tropical maize germplasm 32 .SNP markers are instrumental in the dissection of complex traits such as Striga resistance, and their association with the trait can be revealed through GWAS.Identification of genomic regions linked to Striga resistance in pearl millet breeding would speed up the development of Striga-resistant varieties.Genetic markers improve the efficiency of novel Striga-resistant genes introgression and pyramiding into high-yielding elite varieties through backcross method.Therefore, the objective of this study was to undertake genome-wide association analyses of agronomic traits and Sh resistance among 150 pearl millet genotypes to identify genetic markers for marker-assisted breeding and trait introgression.
Table 2 presented the best linear unbiased prediction means for the response of pearl millet genotypes evaluated under Sh infestation.The DTF and PH ranged from 60.13 to 64.95 days and 132.57 to 159.74 cm in the naturally-infested field and the screen house conditions, respectively.The TGW ranged from 7.60 to 9.03 g under Striga infestation in field and in screen house.Several emerged Striga were recorded during the second counting, particularly in the plastic pots (Fig. 1A) and the hotspot field (Fig. 1B) conditions compared to Strigafree field condition (Fig. 1C).Thousand-grain weight markedly reduced due to high Striga infestations compared to Striga-free field (Table 2).
A low broad-sense heritability value was computed for the number of emerged Sh SN1 (18.17%) and SN2 (28.61%), while high heritability values ranging from 75.07 to 92.42% were recorded for the DTF, PCL, TGW, and PH (

Genome-wide association mapping
The number of emerged Sh count on the Manhattan plot is presented in Fig. 2. BLINK model analysis for Sh traits led to identifying candidate genetic regions associated with Striga resistance.Twenty-eight SNP markers were significantly (P < 0,001) associated with Sh resistance located on chromosomes 1, 2, 3, 4, 6, and 7. Two SNPs, S3_113184999 and S3_113184999, were located on the same position and were associated with the number of Sh counted (Table 3).Three significant (P < 0,001) SNP markers, S1_75620319, S3_1159738, and S6_231436300, were associated with Striga resistance in the naturally Striga-infested field on chromosomes 1, 3, and 6.In the greenhouse, 10 significant SNPs were associated with Sh-resistance, while 15 SNPs were associated with Sh-resistance on chromosomes 1, 2, 3, 4, 6, and 7 in a pooled analysis.Table 1.Mean squares and significant tests for pearl millet and Striga parameters when evaluating 137 genotypes with and without Striga infestation under the field and screen house environments in Burkina Faso.*, and *** = denote significant differences at 0.05, and 0.001 probability levels, respectively; ns = not significant; DF = degree of freedom; DTF = days to 50% flowering; PH = plant height at maturity (cm); NT = number of tillers; PCL = panicle length (cm); PWT = panicle weight (g); TGW = thousand-grain weight (g); GW = grain weight per plant (g); ASNPC = area under the Striga number progress curve.

Agronomic traits
BLINK model analysis of pearl millet agronomic traits under Sh infestation identified candidate genetic regions associated with DTF, PCL, and TGW (Fig. 3).Eleven SNP markers were associated with the assessed pearl millet phenotypic traits.Four SNPs were associated with DTF on chromosomes 3, 5, 6, and 7 in the naturally  www.nature.com/scientificreports/Striga-infested field; five with PCL on chromosomes 2, 3, and 4 with which two in the naturally Striga-infested field and three in the screen house; and two with TGW on chromosome 6 in the screen house (Table 4).

Discussion
The assessed genotypes exhibited significant differences in agronomic traits and Sh parameters (Table 1), suggesting substantial genetic variation for selection.The results allowed marker-trait association analysis to be valuable for current and future selection and new variety design and commercialisation.A significant genotype by the Sh interaction effect existed (Table 1), revealing the potential existence of genes controlling Sh resistance among the populations.This concurs with the deductions made by Mrema et al. 8 and Shayanowako et al. 33 , who reported significant variation and differential genotypic responses to Sh infestation among sorghum and maize genotypes, in that order.Thus, the population set evaluated in the current study was suitable for marker-trait association analyses for Sh resistance and economic traits.Reduced DTF and PH were recorded on pearl millet genotypes under Striga infestation (Table 2), indicating the negative impact of the parasite on the measured traits.These findings concur with reports by Ransom and Odhiambo 34 .Wilson et al. 18 reported a negative correlation between the number of emerged Sh and DTF in pearl millet.Similarly, a reduction of PH by 28% under Sh infestation was reported in pearl millet by Graves et al. 35 .Badu-Apraku 36 reported a negative correlation between Striga damage rating and DTF in maize in Sub-Saharan Africa.Poor crop growth and subsequent low productivity result from the Striga plant attachment and siphoning of nutrients from the host plant's roots.
The recorded TGW of 7.60 to 9.03 g in the present study (Table 2) aligns with the 6.9 to 12.9 g reported by Kanatti et al. 37 .The relatively high ranges of values in the current results are probably attributed to the large grain size, which increased the yield of pearl millet 37,38 .The broad-sense heritability (H 2 ) for agronomic traits ranged from 75.07 to 92.42% and 18.17 to 28.61% for emerged Sh count (Table 2).The high broad-sense heritability estimates of 75.07 to 92.42% computed for DTF, PH, PCL, and TGW indicated that the traits are mainly governed by genes with limited influence by the test environment 39 .Traits with high heritability are easy  www.nature.com/scientificreports/ to select and improve using marker-assisted selections and pyramiding in a desirable genetic background.The lower broad-sense heritability estimates of 18.17 and 28.61% for the number of emerged Striga (SN1 and SN2) (Table 2) suggests that the genetic variation was small and genetic gain for those traits will be slow because both genetic and phenotypic constituents of the genotypes are affected by Striga infestation stress 39 .Robert 40 and Kaewchumnong and Price 41 reported a low heritability estimate for Striga resistance-related traits in sorghum and rice, respectively.The seven pearl millet chromosomes harbour several genes (Table 4) conditioning Striga resistance and agronomic traits.Each chromosome had at least two significant marker-trait associations in the present study.After successful validation, the 28 significant SNP markers associated with Sh emergence on chromosomes 1, 2, 3, 4, 6, and 7 (Table 3) can be used for marker-assisted selection and trait introgression to improve Sh resistance in pearl millet.Dawud 28 identified 16 SNP markers associated with the area under the Striga number progress curve on chromosomes 2, 3, 4, 5, and 7 in pearl millet.The findings confirm that the chromosomes harbour some beneficial alleles influencing Sh resistance.Dawud 28 reported significant gene markers related to Striga resistance on chromosomes 1, 2, 3, and 5 in pearl millet.Markers associated with Striga resistance traits have also been reported in sorghum and maize 44,45 .The two SNP makers located adjacent to each other and associated with the low number of Sh count could be tightly linked and co-segregating.Hence, the respective candidate genes can be selected and introgressed simultaneously 24 .Identifying genetic markers associated with agronomic traits will facilitate marker-assisted breeding in pearl millet (Fig. 3).Using SSR markers, Kannan et al. 46 detected significant markers associated with pearl millet panicle length and thousand-grain weight in Striga-free conditions.In maize, some quantitative trait loci associated with grain yield and ear aspect have been reported by Stanley et al. 47 and Badu-Apraku et al. 23 .Dawud 28 also reported Significant SNP markers related to the number of tillers in pearl millet.This study identified 28 SNP markers associated with low Sh emergence on chromosomes 1, 2, 3, 4, 6, and 7, involving genetic analysis of 150 genetically diverse pearl millet genotypes in Burkina Faso (Fig. 3).The candidate markers and genotypes are novel genomic and genetic resources for Striga resistance breeding programs in the country and elsewhere.

Study sites
A field experiment was conducted in the 2019/20 main growing season in a naturally infested Striga hotspot field at the Didri site in Burkina Faso, and a screen house evaluation was conducted at the main station of the Institute of Environment and Agricultural Research (INERA) in the offseason of 2020/21.The Didri site is located at 12° 12′ 15" N and 1° 14′ 13" W and is a hotspot site for Sh affecting pearl millet, maize and sorghum crops.The site received an annual rainfall of 748.5 mm for 46 days during the 2019/20 rainy season and has sandy soils.The INERA site is located at 12°28/27 N and 1°33/31W.

Plant materials
The study used 148 pearl millet genotypes collected from the International Crop Research Institute for the Semiarid Tropics (ICRISAT) in Niger and two elite breeding lines from INERA/Burkina Faso.The descriptions of the test genotypes are summarised in Table 5.The pearl millet genotypes acquired from ICRISAT are part of the pearl millet germplasm association panel (PMiGAP) comprising 250 inbred lines representing cultivated germplasm from Africa and Asia.They are included in the present study to identify unique genetic resources with unique agronomic and farmers' preferred traits, and because of their wide genetic diversity.

Experimental design and trial management
The field and greenhouse experiments were laid out using a 10 × 15 alpha lattice design with two replicates.In the greenhouse, 5L plastic pots were used and filled with a soil medium composed of clay, sand, and organic manure at a ratio of 2:1:1 respectively.Two weeks before planting, each pot was infested with a scoop of sand mixed with 0.05 g of 1-year-old Sh seed collected from farmers' fields in Burkina Faso.Pearl millet seeds in the naturally Striga-infested field (hereafter designated as GenInf), were sown during the main crop growing season from June to October 2019.Genotypes were established in 4.2 m long rows spaced at an inter-row spacing of 160 cm and intra-row spacing of 60 cm, providing a total plot size of 6.72 m 2 per genotype.Four seeds were initially sown per hill and later thinned to one plant two weeks after planting.A total of three plants was selected randomly from the middle of the experimental unit and tagged for agronomic data collection.In the greenhouse, one healthy and vigorous plant was grown per pot for the test genotypes with Sh (hereafter denoted as GenStr), and genotypes without Sh (hereafter referred to as control) treatments.Standard agronomic practices recommended for pearl millet production were followed.Experimental units were fertilized using nitrogen, phosphorus and potassium (NPK: 14:23:14) and applied as a microdose of 3 g per hill 15 days after planting.Hand weeding was routinely done after the first hoeing to remove all other weeds except Striga.

Phenotypic data
The following agronomic parameters were collected from pearl millet: days-to-50%-flowering (DTF) were recorded as the days when 50% of the plants in each plot had intruded stigma.Plant height (PH) was measured in cm from the base to the top of the panicle of the main tiller.The number of tillers per plant (NT) was recorded by counting the number of tillers with panicles for the tagged plants.Panicle length (PCL) was measured in cm from the base to the top of the main tiller panicle.Panicle weight (PWT) was recorded in grams by weighing the harvested panicles for each entry after 14 days of sun-dry, and thousand-grain weight (TGW) was determined www.nature.com/scientificreports/ in grams by weighing one thousand-grain for each of the entries.Grain weight per plant (GW) was determined in grams by weighing the grain after threshing and dividing it by the number of harvested plants for each plot.
For Striga parameters, the number of emerged Sh plants per plot were recorded at 70 and 96 days after sowing in the naturally infested field for each row, excluding the borders.The number of emerged Sh plants were counted per pot 116 and 144 days after sowing in the greenhouse.The area under the Striga Number Progress Curve (ASNPC) was computed using the successive Striga counts according to Haussmann et al. 48

as follows:
where n is the number of Striga assessment dates, Y i is the Striga count at the ith assessment date, Y (i+1) is the Striga count at the ith plus 1 assessment date, t i is the number of days after planting (DAP) at the ith assessment date, t (i+1) is DAP at the ith plus 1 assessment date.

Phenotypic data analysis
Both the crop and Striga data were subjected to analysis of variance using GenStat 19th Edition (http:// www.genst at.co.uk).Homogeneity of variance test was done for each site using the Bartlett 49 procedure before combined analyses.The treatment, genotype, and genotype × treatment interaction significance tests were computed using GenStat.The Best Linear Unbiased Prediction (BLUP) was calculated according to Haslett and Puntanen 50 to predict the accuracy and to aid selection.The area under the Striga number progress curve was drawn using R. ASReml-R Version 4 was used to fit the linear mixed models using Residual Maximum Likelihood (REML) in R 51 .

Genotyping and GWAS analysis
To assemble the pearl millet genome, whole genome shotgun (WGS) and bacterial artificial chromosome (BAC) sequencing were used.Ten small inserts (of ~ 170, 250, 500 and 800 bp) and 13 large inserts (of ~ 2, 5, 10, 20 and 40 kb) WGS libraries were constructed using Tift 23D2B1-P1-P510 genotype.These libraries were sequenced on the Illumina HiSeq 2000, and 520 Gb of sequence data, representing 296 × genome coverage.Two BAC libraries, with an average insert size of ~ 120 kb, were constructed from Tift 23D2B1-P1-P5 using EcoRI and HindIII.Nine hundred seventy-two Gb of sequence data were generated from 100,608 BAC clones at ~ 80 × genome coverage.In brief, 1.49 Tb of sequence data, after stringent filtering and correction steps, were assembled into 1.58 Gb of contigs and 1.82 Gb of scaffolds (https:// doi.org/ 10. 1038/ nbt.3943).A raw marker set consisting of 29 million SNPs generated from 345 pearl millet genotypes from which 148 genotypes used in the study was sourced from the Nature Paper Pearl Millet (https:// doi.org/ 10. 1038/ nbt.3943) and filtered using Tassel v4.2 for site coverage of 90%, minor allele frequency of 0.1, taxa coverage of 30% and maximum heterozygosity of 50%.The resulting final set of variants contained 256 K SNPs was used for the current analysis.Phenotypic data collected from 150 genotypes were available for marker-trait association analysis.Principal component analysis (PCA) was calculated, and the resulting eigenvalue (7) was used for genome-wide association analysis (GWAS) following multiple methods procedures generated from GAPIT v3.0 52 .The Bayesian-information and linkagedisequilibrium iteratively nested keyway (BLINK) software was used to determine the significant variations among pearl millet and Sh traits owing to its ability to produce fewer false positives and more true positives than the GWAS method, FarmCPU 42 .Liu et al. 43 reported the power of BLINK to outperform FarmCPU relative to statistical capabilities vs False Discovery Rate (FDR) and statistical power vs type I error.

Conclusion
The current study detected significant genetic variability and markers for Sh resistance and agronomic traits through GWAS involving 150 pearl millet genotypes in Burkina Faso.There were significant genotypes by Sh interaction for assessed agronomic traits, the number of Sh and ASNPC.Twenty-eight SNPs were associated with Sh traits on chromosomes 1, 2, 3, 4, 6, and 7. SNPs markers associated with DTF, PCL, and TGW were located on chromosomes 2, 5, 6, and 7; chromosomes 2, 3, and 4; and chromosome 6, respectively.After successful validation, the new markers would be deployed for marker-assisted breeding emphasising the above agronomic traits and Sh resistance in pearl millet in Burkina Faso or related agro-ecologies.
naturally Sh infested field and their genetic parameters.Trt = treatment; DTF = days to 50% flowering; PH = plant height at maturity (cm); PCL = panicle length (cm); TGW = thousand-grain weight (g); SN1 = the number of Striga counted 70 days after sowing in the Striga-infested field and 116 days after planting in the screen house; SN2 = Striga number counted 96 days after sowing in the Striga-infested field and 144 days after sowing in the screen house; GenInf = genotypes in naturally Striga infested field; GenStr = genotypes with Striga infestation in the screen house; SED = standard error of the mean difference; H 2 (bs) = broad sense heritability; LSD = least significant difference; CV = coefficient of variation; ns = non-significant.

Figure 1 .
Figure 1.Pearl millet crop and Striga infestation in Burkina Faso.Note: Striga hermonthica infestation under screen house (A) and field (B) conditions.Photo C denotes Striga free field at the Didri site in Burkina Faso.(Photos supplied by Armel Rouamba).

Figure 2 .
Figure 2. Manhattan plots showing associations between single nucleotide polymorphisms and SN1 and SN2 in naturally Striga-infested fields and greenhouse conditions.Single nucleotide polymorphisms were plotted on the x-axis according to their positions on each chromosome against association with Sh-related traits on the y-axis (−log 10 p-value).The top line indicates genome-wide significant threshold.Note: Blink.Field_SN2 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 96 days after sowing in the Striga-infested field; Blink.Pooled_SN1 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 70 days after sowing in the Strigainfested field and 116 days after planting in the greenhouse; Blink.Pooled_SN2 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 96 days after sowing in the Strigainfested field and 144 days after sowing in the greenhouse; Blink.Greenhouse_SN1 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 116 days after planting in the greenhouse, Blink.VarStr_SN1 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 116 days after planting in the greenhouse for genotypes with Striga infestation; Blink.VarStr_SN2 = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of Striga number counted 144 days after sowing in the greenhouse for genotypes with Striga infestation.

Figure 3 .
Figure 3. Manhattan plots showing associations between single nucleotide polymorphisms and panicle length, flowering time and thousand-grain weight under Sh conditions in naturally Striga infested field and greenhouse.Single nucleotide polymorphisms were plotted on the x-axis according to their positions on each chromosome against association with pearl millet-related traits on the y-axis (−log 10 p-value).The top line indicates genomewide significant threshold.Note: Blink.Greenhouse_PclL = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of panicle length in screen house; Blink.Field_PclL = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of panicle length in the field; Blink.VarStr_PclL = the Bayesianinformation and linkage-disequilibrium iteratively nested keyway of panicle length for genotypes with Striga infestation; Blink.Greenhouse_TGW = the Bayesian-information and linkage-disequilibrium iteratively nested keyway of thousand-grain weight in the screen house, Blink.Var_TGW = the Bayesian-information and linkagedisequilibrium iteratively nested keyway of thousand-grain weight for genotypes in the screen house. https://doi.org/10.1038/s41598-023-44046-1

Table 2 .
Best linear unbiased prediction means and standard error for 137 pearl millet genotypes evaluated in

Table 3 .
Markers significantly associated with Striga resistance traits in 150 pearl millet genotypes assessed in naturally infested fields and greenhouse conditions in Burkina Faso.SNP = single nucleotide polymorphism, Chr.= chromosome, MAF = minor allele frequency, AdjP = false discovery rate adjusted P-values, SN1 = the number of Striga counted 70 days after sowing in the Striga-infested field and 116 days after planting in the screen house; SN2 = Striga number counted 96 days after sowing in the Striga-infested field and 144 days after sowing in the screen house.

Table 4 .
Markers significantly associated with agronomic traits in 150 pearl millet genotypes assessed in naturally Striga-infested fields and screen house conditions in Burkina Faso.SNP = single nucleotide polymorphism, Chr.= chromosome, MAF = minor allele frequency, AdjP = false discovery rate adjusted P.values.

Table 5 .
Description of the pearl millet genotypes used in the study.
E.N. = entry number; ICRISAT = International Crop Research Institute for the Semi-Arid Tropics; INERA = Institute of Environment and Agricultural Research/Burkina Faso; C = check; R = resistance; T = tolerance; S = susceptible,-= denote data not available.