
==== Front
BMC Plant Biol
BMC Plant Biol
BMC Plant Biology
1471-2229
BioMed Central London

39294608
5590
10.1186/s12870-024-05590-8
Research
Genome-wide association analysis of grain yield and Striga hermonthica and S. asiatica resistance in tropical and sub-tropical maize populations
Dossa Emeline N. dossaemeline@gmail.com

Shimelis Hussein
Shayanowako Admire I. T
https://ror.org/04qzfn040 grid.16463.36 0000 0001 0723 4123 School of Agricultural, Earth and Environmental Sciences, University of KwaZulu-Natal, Private Bag X01, Scottsville, Pietermaritzburg, 3209 South Africa
19 9 2024
19 9 2024
2024
24 87112 6 2024
12 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Background

Genetic improvement for Striga hermonthica (Sh) and S. asiatica (Sa) resistance is the most economical and effective control method to enhance the productivity of maize and other major cereal crops. Hence, identification of quantitative trait loci (QTL) associated with Striga resistance and economic traits will guide the pace and precision of resistance breeding in maize. The objective of this study was to undertake a genome-wide association analysis of grain yield and Sh and Sa resistance among tropical and sub-tropical maize populations to identify putative genetic markers and genes for resistance breeding. 126 maize genotypes were evaluated under controlled environment conditions using artificial infestation of Sh and Sa. The test genotypes were profiled for grain yield (GY), Striga emergence counts at 8 (SEC8) and 10 (SEC10) weeks after planting, and Striga damage rate scores at 8 (SDR8) and 10 (SDR10) weeks after planting. Population structure analysis and genome-wide association mapping were undertaken based on 16,000 single nucleotide polymorphism (SNP) markers.

Results

A linkage disequilibrium (LD) analysis in 798,675 marker pairs revealed that 21.52% of pairs were in significant linkage (P < 0.001). Across the chromosomes, the LD between SNPs decayed below a critical level (r2 = 0.1) at a map distance of 0.19 Mbp. The genome-wide association study identified 50 significant loci associated with Sh resistance and 22 significant loci linked to Sa resistance, corresponding to 39 and 19 candidate genes, respectively.

Conclusion

The study found non-significant QTL associated with dual resistance to the two examined Striga species Some of the detected genes reportedly conditioned insect and pathogen resistance, plant cell development, variable senescence, and pollen fertility. The markers detected in the present study for Sa resistance were reported for the first time. The gene Zm00001eb219710 was pleiotropic, and conditioned GY and SEC10, while Zm00001eb165170 affected SDR8 and SDR10, and Zm00001eb112030 conditioned SDR8 and SDR10 associated with Sh resistance. The candidate genes may facilitate simultaneous selection for Sh and Sa resistance and grain yield in maize after further validation and introgression in breeding pipelines. Overall, we recommend breeding maize specifically for resistance to each Striga species using germplasm adapted to the endemic region of each parasite.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12870-024-05590-8.

Keywords

Grain yield
GWAS
Maize
Resistance breeding
Striga resistance
issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Maize (Zea mays L., 2n = 2x = 20) is the global food, feed and industrial cereal. In sub-Saharan Africa (SSA), the per capita consumption of maize is > 330 g/person/day [1]. The productivity of maize in SSA has remained stagnant and low (< 3 t/ha) compared to the potential productivity reaching 5 to 10 t/ha [2]. The parasitic weeds, Striga hermonthica (Del.) Benth and S. asiatica (L.) Kuntze are among the leading causes of low productivity of maize, especially in communal and small-scale maize production systems with low agricultural inputs [3, 4].

Striga hermonthica (Sh) and S. asiatica (Sa) are cosmopolitan parasitic weeds and destructive species affecting major cereals crops, including maize in SSA [3–6]. The parasites damages the host species by extracting its photoassimilates, diminishing its growth and productivity, especially under moisture stress and sub-optimal soil nutrient conditions, which are common in marginal maize production areas of SSA [7]. Under severe Striga infestation of both species, yield loss reaching 100% has been recorded in susceptible maize varieties [8]. Hence, there is a need for genetic improvement for Striga resistance and grain yield to enhance the productivity of maize and other major cereal crops in SSA.

Various Striga control strategies have been recommended, such as cultural practices (e.g., crop rotation, intercropping, catch-cropping, trap-cropping, and the application of organic and inorganic fertilizers), chemical and bio-herbicides. The control practices are used solo or in combination, which is referred to as integrated Striga management [9, 10]. Most of the described control options have limited application, and some are inaccessible and unaffordable to many smallholder farmers [8, 11, 12]. Striga-resistant varieties are the most economical and environmentally friendly option to control the parasitic weed under low-input farming systems [13, 14].

Several Striga-resistant varieties have been developed by the International Institute of Tropical Agriculture (IITA) in West Africa, where Sh is prevalent. However, only partial to moderate Striga resistance has been attained so far [15, 16], due to the low heritability of Striga resistance traits [8]. Low heritable traits are subject to genotype, environment, and genotype x environment interaction effects [17]. Molecular marker-assisted breeding tools could improve selection gains for Striga resistance and low heritable traits. Molecular markers can identify, locate, and map genes conditioning economic traits and Striga resistance. Linkage mapping has successfully identified quantitative trait loci (QTL) for complex traits using bi-parental populations [18]. However, detecting recombination events within pedigree and families is minimal in linkage mapping, leading to low mapping resolution of genetic markers and genes. Genome-wide association study (GWAS) is a cost-effective tool for discerning marker-trait association and dissecting the genetic architecture of complex quantitative traits, which can provide a high-resolution and high allelic richness. GWAS saves time in genetic analysis compared to linkage mapping for identifying QTL [19]. Various populations can be used in association mapping, including elite cultivars, landraces, wild relatives, and introductions [20, 21]. Limited studies documented QTL mapping associated with Striga resistance, yield, and yield component traits.

The first QTL associated with Striga resistance in maize were reported by Amusan [22]. The authors reported two QTL mapped on chromosome 6 using an F2 mapping population involving a cross between a susceptible (5057) and a resistant (ZD05) maize Inbred lines. The two QTL accounted for 55% of observed phenotypic variation for incompatible responses to Sh infestation on host roots. Adewale, et al. [13] used GWAS and reported three markers located close to the putative genes GRMZM2G164743 (bin 10.05), GRMZM2G060216 (bin 3.06), and GRMZM2G103085 (bin 5.07), linked to grain yield, the number of ears per plant, and Striga damage under Striga infestation. The three QTL explained 9 to 42% of the phenotypic variance for the incompatible response of Sh on the agronomic traits. Further, using QTL mapping, Badu-Apraku et al. (2020a) identified 12 QTL associated with Sh resistance traits from an F2:3 population. The authors further identified 14 other QTL and 154 candidate genes associated with Striga resistance/tolerance traits using QTL mapping [23]. Gowda, et al. [24] identified 57 SNPs significantly associated with Sh resistance indicator traits and grain yield, pin-pointing 32 candidate genes near the significant SNPs using GWAS. Recently, 22 SNP markers associated with grain yield, Striga damage, number of emerged Striga plants after planting, and ear aspect were reported by Okunlola, et al. [25] using GWAS. So far, only one study [8] reported on QTL associated with Sa resistance. The authors identified 3 SNP markers on chromosome numbers 5, 6, and 7 for total Sa plants that emerged. There is a need to identify QTL associated with Striga resistance and economic traits to guide the pace and precision of resistance breeding in maize aiming Sh and Sa.

Thus far, QTL associated with Striga-resistance have been detected using a bi-parental population, which may have low mapping resolution due to the small number of accumulated recombination events [26]. The quantitative trait loci intervals found are extended over several centimorgans (cM). Compared with Sh, there is a dearth of information on major QTL for Sa resistance. Hence, this study will be the first attempt to report QTL associated with both Sa and Sh resistance for resistance breeding in SSA. The two species are the most devastating parasitic weed populations occurring in tandem over space and time in SSA, especially in East Africa. Therefore, the objective of this study was to undertake a genome-wide association analysis of grain yield and Striga resistance among 126 tropical and subtropical maize genotypes to identify genetic markers linked to resistance to Sh and Sa. The QTL associated with resistance to both parasites will be useful in developing Striga-resistant maize with wide adaptability across the SSA region.

Materials and methods

Plant material and phenotyping

The current study used a panel of 126 maize germplasm. Of these, 70 were acquired from the Striga-resistant pool developed by IITA/Nigeria, 45 were locally adapted varieties from the International Maize and Wheat Improvement Centre (CIMMYT)/Zimbabwe, and 11 were from the National Plant Genetic Resources Centre, South Africa (NPGRC)/South Africa (Supplemental Table 1). The germplasm included tropical and sub-tropical inbred lines, open-pollinated varieties (OPVs), and top-cross hybrids, as described in [27]. The 126 genotypes were rigorously phenotyped under Sa and Sh-infested conditions separately at the University of KwaZulu-Natal Controlled Environment Facilities (UKZN-CEF)/South Africa in two cropping seasons from December 2021 to April 2022 and from August 2022 to December 2022. The average maximum temperature during those seasons is between 26 and 28 °C, while the minimum is 10 °C. Each of the experiments were arranged in a 13 × 10 alpha lattice design with two replications. The UKZN-CEF is situated at the UKZN College of Agriculture, Engineering, and Science (29.62° S, 30.40° E). Two weeks before planting, each pot was infested with a scoop of sand mixed with 0.03 g of 2-year-old Sa or Sh seed containing approximately 3000 Striga seeds.

Grain yield and Striga resistance parameters were collected. The GY (expressed in g/plant) was determined as the weight of the grain from the ears of an individual plant after shelling, adjusted to 12.5% moisture content. Striga parameters comprised the number of emerged Sa and Sh plants 8 and 10 weeks after planting, denoted as Sect. 8 and Sect. 10. A rating of host plant damage 8 and 10 weeks after planting, designated as SDR8 and SDR10, was done using a visual score of 1 to 9, where 1 = no damage, indicating normal plant growth and a high level of tolerance, and 9 = complete collapse or death of the maize plant, i.e., highly susceptible [28].

Genotyping

Following the plant DNA extraction protocol for DArT, genomic DNA was extracted from fresh leaves of 126 genotypes at the three-leaf stage [29]. The quality and purity of the extracted DNA were rigorously assessed using the NanoDrop 2000 spectrophotometer (ND-2000 V3.5, NanoDrop Technologies Inc) [30]. Each genotype provided a 20 µl DNA sample with concentrations between 50 and 100 ng and absorbance values ranging from 1.75 to 2.05. These samples were sent to Sequential Art (SEQAT) in Kenya (https://www.seqart.net/) for high-throughput genotyping using the DArTseq protocol. The genotyping yielded 70,197 SNPs, of which 16,000 informative SNPs were selected for further analysis after filtering out DArT loci with unknown chromosomal positions and markers with more than 20% missing data.

Data analysis for grain yield and Striga parameters

Data collected on GY and Striga parameters were subjected to Bartlet’s homogeneity of variance test prior to analysis of variance (ANOVA) using a lattice procedure, using the Deltagen platform. Genotypes mean comparison were done at the 5% significance level using Fisher’s least significance difference (LSD). Broad sense heritability (H2) was computed using DeltaGen with the following formula:H2=σ2gσ2g+σ2sns+σ2rnr+σ2bnb+σε2ns+nr+nb

where σ2g,σ2s,σ2r,σ2b,and σε2 are the variance components for genotype, season, replication, block, and the pooled error, respectively, while ns, nr, and nb are the number of seasons, replications, and blocks, in that order. The linear mixed model was best-fit by the restricted maximum likelihood and used for phenotypic data analysis.

Genome-wide analysis to decipher marker-trait associations

The initial 70,197 SNPs were imputed by removing SNPs with > 20% missing data and < 5% minor allele frequency (MAF) on the KDCompute server (https://kdcompute.seqart.net/). A filtered set of 16,000 highly informative SNP markers was used for GWAS. The 126 maize genotypes were assessed for kinship using the software Structure [31]. The parameters of the analysis were set at 10,000 burn-in periods, with 10,000 Markov chain–Monte Carlo (MCMC) repetitions after burn-in. Five iterations were run for population number (K) values of 1 to 10 to allow the selection of the replication with the highest mean value of ln likelihood. The optimum K value was determined by the ad-hoc delta K method [32]. The linkage disequilibrium (LD) decay was estimated using the LD function in TASSEL software version 5.2.92 [33], and plotted with R version 4.3.0 (R Core Team, 2022) as described by Remington, et al. [34]. To identify the possible number of loci that are associated with Striga resistance in maize, pairwise LD of the populations was estimated using squared allele frequency correlations (r2) for the DArTseq SNPs markers [30, 35]. This was achieved using the LD measure, r2 program (Version 2.1) within the KDCompute plugin system. The r2 was estimated by an LD sliding window size of 50 and a threshold of r2 set at 0.1, such that any SNP that was below 0.1 was considered to have a weak LD [36]. The SNPs significantly associated with Striga resistance were compared with the rest of the SNPs within the linkage group. The distribution pattern for the whole genome LD was visualized using graphs generated as LD heatmap from TASSEL v5.2.5. The inter-chromosomal LD among the SNPs was achieved by examining the LD data to determine whether some SNPs had an r2 equal to or above 0.1. SNPs significantly associated with the traits of interest were recorded for each chromosome pair. Analysis for the genomic regions associated with Striga resistance and grain yield, probability values, and percentage of the effect of the markers were computed using the genome association and prediction integrated tool (GAPIT) package via the KDCompute interface (https://kdcompute.seqart.net/).

The associations between Single Nucleotide Polymorphisms (SNPs) and phenotypic traits were analyzed using the Compressed Mixed Linear Model (CMLM) implemented in the GAPIT software, as outlined by Yu et al. (2006). In this model, SNP markers were treated as fixed effects and were assessed individually. The model can be expressed as: Y = Xβ + Wα + Qv + Zu + ε, where Y is the observed vector for the phenotypic records of the traits; β is the fixed-effect vector (p × 1) other than the molecular marker effects and the population structure; α is the fixed-effect vector of the molecular markers; v is the fixed-effect vector from the population structure; u is the random-effect vector from the polygenic background effect; X, W, Q, and Z are the incidence matrixes from the associated β, α, v, and u parameters; and ε is the residual effect vector, in that order [37]. The P values and false discovery rate thresholds were set at P < 1 × 10−4 and – log(p) = 3 respectively. The kinship matrix was calculated from the 16,000 markers, and the quantile-quantile (QQ) and Manhattan plots were generated using the log10 (p) value distributions of the GAPIT function in KDCompute.

The SNP(s) above the threshold value were used to indicate the genomic region associated with Striga resistance and grain yield traits. Further, the QQ plot was used to evaluate how well the model used in GWAS for this study accounted for population structure. The significant markers identified for each trait were blasted on Ensemble based on the maize genome version B73_V4.0. assembly to identify candidate genes associated with the markers using the maizeGDB (https://www.maizegdb.org/) platform.

Results

Phenotyping for grain yield response and Striga parameters

Analysis of variance indicated significant genotypic variation (P < 0.05) for grain yield and Striga parameters in both Sa and Sh-infested environments (Table 1). Testing seasons had significant effects (P < 0.001) on GY and Striga parameters under both Striga-infested conditions except for SEC10 under Sa-infested conditions. Also, significant effects of replications in seasons were noted for the assessed traits under the two conditions, except for SDR8 and SDR10 under Sa conditions. Significant differences were also attributed to the block nested into replication-by-season interaction effect under both Sa and Sh-infested environments except for GY under Sa-infested conditions and SDR10 under Sh-infested conditions.

Table 1 Analysis of variance and significance tests for grain yield and Striga parameters of 126 maize genotypes evaluated under Striga asiatica and S. Hermonthica infested conditions

S. asiatica	
Source of variation	Df	GY (g/plant)	SEC8	SEC10	SDR8	SDR10	
Genotypes (G)	125	1205.52***	0.32***	52.94***	0.78***	0.55***	
Seasons (S)	1	268.64***	6.55**	0.70	53.50***	0.26*	
G x S	125	0.00	0.45	0.00	0.00	0.00	
Replications in seasons	1	58.37***	2.22***	34.96***	0.40	0.04	
Block/(replication x season)	13	0.00	0.79***	2.69**	0.32***	0.73***	
Residual	238	924.67	0.95	100.66	1.54	1.47	
S. hermonthica	
Source of variation	Df	GY	SEC8	SEC10	SDR8	SDR10	
Genotypes (G)	125	579.51***	4.29***	1.75***	0.88***	0.93***	
Seasons (S)	1	460.89***	8.43***	1.73***	1.88***	0.44***	
G x S	125	152.51	0.00	0.46	0.00	0.00	
Replications in seasons	1	270.00***	2.27***	1.47***	1.97***	0.92***	
Block/(replication x season)	13	41.39**	0.61*	0.33**	2.44**	0.03	
Residual	238	1126.63	5.84	4.83	2.21	1.28	
*,**, and *** denote significance at P < 0.05, P < 0.01 and P < 0.001, respectively

Df degrees of freedom, GY grain yield, SEC8 Striga emergence counts eight weeks after sowing, SEC10 Striga emergence counts ten weeks after sowing, SDR8 Striga damage rating eight weeks after sowing, and SDR10 Striga damage rating 10 weeks after sowing

Low broad-sense heritability values were recorded for SDR8 (< 0.10) and SDR10 (0.11) in Sa-infested conditions. In contrast, a high heritability value was recorded for Sect. 8, Sect. 10, SDR8, and SDR10 (H2 > 0.50) but low for GY (H2 = 0.02) under Sh-infested conditions (Table 2).

Table 2 Broad sense heritability estimates for grain yield and Striga resistance traits of 126 maize genotypes evaluated over two seasons, under two Striga-infested conditions and two replications

	S. asiatica	S. hermonthica	
Parameter	GY (g/plant)	SEC8	SEC10	SDR8	SDR10	GY	SEC8	SEC10	SDR8	SDR10	
Broad sense heritability	0.96	0.88	0.34	0.01	0.11	0.02	0.88	0.87	0.92	0.82	
LSD (5%)	1.93	37.73	1.87	12.44	1.64	1.39	34.87	34.87	2.37	1.62	
GY grain yield, SEC8 Striga emergence counts eight weeks after sowing, SEC10 Striga emergence counts ten weeks after sowing, SDR8 Striga damage rating eight weeks after sowing, and SDR10 Striga damage rating 10 weeks after sowing, LSD least significant difference

Population structure and linkage disequilibrium

After removing markers with > 20% missing values and a marker frequency (MAF) < 5%, 16,000 SNP markers were found polymorphic and were used for GWAS. The 16,000 markers were distributed on the 10 chromosomes of the maize genome, with the highest and lowest marker densities observed on chromosome 8 (1272 SNPs) and chromosome 4 (835 SNPs), respectively. A population structure was constructed to reveal genetic relationships and aid genotype selection using the software Structure. The highest value for ΔK occurred at K = 8, showing that the genotypes could be clustered into eight sub-populations (Fig. 1A and B). The minimum coefficient for membership to a particular sub-population was 0.70. A similar trend was also observed using the Kinship analysis, where eight clusters were identified (Fig. 2).

Fig. 1 The population structure ΔK (A) and the eight populations (B) of 126 maize genotypes resolved by the Evanno method based on 16,000 SNP markers

The LD was estimated by calculating the squared correlation coefficient (r2) for all the 16,000 SNPs. Pairwise LD analysis between 16,000 SNP markers generated 798,675 comparisons within a physical distance extending up to 40,000,000 bp and was found to decay rapidly with the genetic distance. It was observed that LD varied along the chromosomes, with regions of high LD interspersed with regions of low LD (Fig. 3A). About 353,719 (21.52%) loci pairs were in significant LD (P < 0.001). Further, 78,157 (9.78%) were in complete LD (r2= 1), while 0.19% were at completely no LD (r2 = 0). A critical value of r2 was calculated from inter-chromosomal LD analysis and is estimated to be 0.1, beyond which LD is assumed to be caused by genetic linkage. The point at which the locally estimated scatterplot smoothing (LOESS) curve intercepts the critical r2 is determined as the average LD decay of the population. Based on these criteria, the intra-chromosomal LD decayed at 0.19 Mb for the whole genome (Fig. 3B). The significant intra-chromosomal LD (P < 0.001) ranged from 0.00 to 1.00 with an average of 0.11 Mb for the whole panel.Fig. 2 Kinship heatmap showing the phylogenetic relationship among 126 maize genotypes using 16,000 SNP markers. Red small squares indicate genetic relatedness, dark orange colours represent high kinship relations, while lighter colours (yellow) represent weaker relations. The dendrogram out of the matrix indicated the resulting clusters

Fig. 3 Plots showing linkage disequilibrium (LD) heatmap (A) and intra-chromosomal LD decay (B) for 126 maize genotypes using 16,000 SNP markers. Note: in Figure A the triangle on the right diagonal revealed the squared correlation coefficients (r2), while the P values are shown on the lower diagonal of the triangle for each marker pair. In Figure B, the horizontal blue line indicates the 95th percentile distribution of unlinked r2, while the red cover shows the locally estimated scatterplot smoothing (LOESS) illustrating the LD decay

Marker-trait association (MTAs) under S. Asiatica and S. hermonthica-infested conditions

Grain yield and Striga parameters (SEC8, SEC10, SDR8, and SDR10) were subjected to GWAS using 16,000 SNP markers. The traits SDR8 under Sa-infested conditions and GY under Sh-infested conditions were not included in the GWAS due to their low heritability (H2 < 0.1). The GWAS results across environments are shown in Manhattan plots and Q-Q plots of P values comparing the expected − log10 p values to the observed − log10 p values in Figs. 4 and 5 for Sa and Sh environments, respectively. A total of 72 MTAs were identified at P ≤ 1 × 10−4 (Tables 3 and 4). Two MTAs were detected for GY under Sa (Table 3; Fig. 4), while seven significant MTAs were associated with the same trait under Sh-infested conditions (Table 4; Fig. 5). There were two, seven, and eleven MTAs for SEC8, SEC10, and SDR10, respectively, under Sa-infested conditions (Table 3; Fig. 4), while twelve, fourteen, seventeen, and seven MTAs were detected for the same traits, respectively under Sh conditions (Table 4; Fig. 5). Table 3 Striga asiatica resistance-associated SNP markers and candidate genes identified in the association mapping panel of 126 maize genotypes

Traits	Marker codes	Chr	Position	p value	MAF	R2	Allelic effect	Gene ID	Annotation	
GY	2,459,912	chr1	2.87E + 08	9.67 × 10−6	0.15	0.25	2.86			
	100,065,851	chr5	19,528,687	5.10 × 10−5	0.10	0.22	0.47			
SEC8	100,034,001	chr10	1.49E + 08	2.50 × 10−5	0.12	0.16	0.59	Zm00001eb432870	Putative oligouridylate binding protein	
	100,096,602	chr4	2.15E + 08	8.60 × 10−5	0.13	0.14	-0.54			
SEC10	2,567,877	chr3	2.23E + 08	1.60 × 10−6	0.11	0.20	5.89	Zm00001eb159960	Transcription factor bHLH87	
	100,070,348	chr7	9,877,800	4.50 × 10−6	0.35	0.19	4.03	Zm00001eb301250	Ethylene-overproduction protein 1	
	4,771,900	chr10	1.41E + 08	5.90 × 10−6	0.39	0.18	3.19	Zm00001eb430230	UBC core domain-containing protein	
	4,590,609	chr1	1.94E + 08	6.30 × 10−6	0.16	0.18	-4.88	Zm00001eb036050	Ras-related protein ARA-4	
	4,774,001	chr3	2.23E + 08	8.60 × 10−6	0.26	0.18	3.54	Zm00001eb160080	HSF transcription factor	
	100,047,118	chr5	74,829,144	9.20 × 10−6	0.42	0.17	-2.70	Zm00001eb230770	tRNA (guanine(26)-N(2))-dimethyltransferase	
	4,584,680	chr1	4,146,000	9.70 × 10−6	0.38	0.17	-3.54	Zm00001eb001400	NADH ubiquinone oxidoreductase B14 subunit	
SDR10	4,579,420	chr7	1.48E + 08	8.70 × 10−5	0.22	0.21	0.87	Zm00001eb319510	Oxidoreductase-like domain-containing protein	
	4,582,405	chr3	2.08E + 08	1.42 × 10−6	0.25	0.20	0.054			
	4,770,550	chr8	16,558,612	3.20 × 10−6	0.22	0.19	0.63			
	100,092,893	chr2	2.33E + 08	4.90 × 10−6	0.13	0.19	0.66	Zm00001eb336310	Myosin-binding protein 7; GTD-binding domain-containing protein	
	100,048,286	chr7	1.76E + 08	5.05 × 10−6	0.19	0.18	0.85	Zm00001eb114280	PCI domain-containing protein; Uncharacterized protein	
	2,464,091	chr6	1.25E + 08	6.09 × 10−6	0.27	0.18	0.52			
	4,773,089	chr3	2.13E + 08	6.76 × 10−6	0.28	0.18	0.56	Zm00001eb156390	Uncharacterized protein	
	7,060,947	chr10	16,779,790	6.97 × 10−6	0.25	0.18	0.62	Zm00001eb409270	Ureide permease 5	
	2,431,593	chr1	1.98E + 08	7.50 × 10−6	0.21	0.18	0.36			
	100,056,350	chr5	15,102,289	9.20 × 10−6	0.30	0.18	0.54	Zm00001eb217880	Protein CDC73-like protein	
GY grain yield, SEC8 Striga emergence counts 8 weeks after sowing, SEC10 Striga emergence counts 10 weeks after sowing, SDR10 Striga damage rating 10 weeks after sowing, MAF minor allele frequency, and R2 squared correlation coefficient, Gene ID gene identity

Table 4 Striga hermonthica resistance-associated SNP markers and their candidate genes identified in the association mapping panel of 126 maize genotypes

Traits	Marker codes	Chr	Position	p value	MAF	R2	Allelic effect	Gene ID	Annotation	
SEC8	100,083,158	chr2	2.27E + 08	9.5 × 10−6	0.25	0.41	0.42			
	4,773,476	chr7	1.22E + 08	1.75 × 10−6	0.18	0.41	0.42			
	100,095,995	chr6	1.46E + 08	1.84 × 10−6	0.29	0.41	0.45	Zm00001eb313080	protein_coding	
	4,582,578	chr5	2.06E + 08	3.77 × 10−6	0.21	0.40	0.54			
	100,085,103	chr9	1.28E + 08	4.92 × 10−6	0.16	0.39	0.42			
	2,375,758	chr5	2.13E + 08	5.51 × 10−6	0.19	0.39	-0.05			
	100,077,767	chr5	4,211,185	6.16 × 10−4	0.20	0.39	0.14	Zm00001eb253620	Sadenosylmethionine carrier 1	
	4,775,335	chr6	88,822,670	7.25 × 10−4	0.15	0.39	0.59			
	100,072,527	Chr10	1.39E + 08	7.99 × 10−4	0.13	0.39	0.41			
	100,081,160	chr3	2.33E + 08	8.01 × 10−4	0.35	0.39	0.50	Zm00001eb163250	Clathrin heavy chain	
	2,443,196	chr1	2.91E + 08	8.98 × 10−4	0.19	0.39	0.49	Zm00001eb060090	Nuclear pore complex protein NUP1	
	4,772,592	chr10	1.33E + 08	9.96 × 10−4	0.14	0.39	0.42	Zm00001eb427110	Serine/arginine repetitive matrix protein 2	
SEC10	4,771,611	chr3	27,066,610	2.7 × 10−6	0.12	0.32	1.65	Zm00001eb126130	protein_coding	
	2,378,993	chr3	27,133,911	4.9 × 10−5	0.15	0.31	1.59	Zm00001eb428830	Dof zinc finger protein DOF4.6	
	5,583,989	chr10	1.28E + 08	6.92 × 10−5	0.09	0.30	1.79			
	2,447,143	chr5	66,769,956	1.1 × 10−4	0.25	0.30	1.06			
	4,773,342	chr10	20,619,497	1.35 × 10−4	0.18	0.29	1.18	Zm00001eb409850	Phosphoglucan	
	2,429,588	chr6	1.33E + 08	3.29 × 10−4	0.36	0.28	0.84	Zm00001eb281590	proline-rich receptor-like protein kinase PERK8	
	2,448,259	chr5	76,338,383	3.30 × 10−4	0.24	0.28	1.17			
	2,404,041	chr10	1.38E + 08	3.37 × 10−6	0.22	0.28	1.20	Zm00001eb428830	Dof zinc finger protein DOF4.6	
	2,543,649	chr8	1.68E + 08	4.06 × 10−4	0.15	0.28	1.30			
	4,770,876	chr5	20,571,879	4.78 × 10−6	0.31	0.28	0.93	Zm00001eb219710	Terpene cyclase/mutase family member	
	4,586,050	chr7	20,584,385	6.16 × 10−6	0.30	0.27	0.95	Zm00001eb303380	Eukaryotic translation initiation factor 3 subunit E	
	2,485,365	chr1	4,403,200	6.2 × 10−6	0.17	0.27	-0.13			
	2,468,459	chr10	1.22E + 08	6.51 × 10−4	0.30	0.27	0.80	Zm00001eb001470	Transcription factor TFIIIC	
	2,463,944	chr1	2.37E + 08	7.57 × 10−4	0.22	0.27	1.13	Zm00001eb046220	Protein kinase domain-containing protein	
SDR8	4,764,587	chr8	11,275,425	1.99 × 10−6	0.18	0.20	1.01	Zm00001eb334960	Actin-interacting protein 1–2	
	4,772,065	chr6	1.68E + 08	1.46 × 10−6	0.46	0.16	-0.46	Zm00001eb295600	Golgin candidate 5	
	9,681,009	chr2	2.23E + 08	2.30 × 10−6	0.45	0.16	0.42	Zm00001eb112030	Sister chromatid cohesion protein DCC1	
	100,068,489	chr10	81,558,128	2.55 × 10−5	0.18	0.15	0.70			
	2,512,743	chr4	2,964,998	2.79 × 10−5	0.37	0.15	0.97	Zm00001eb165170	protein_coding	
	2,554,932	chr8	822,625	3.20 × 10−5	0.28	0.15	0.47	Zm00001eb332250	protein_coding	
	2,447,891	chr3	1.77E + 08	3.7 × 10−5	0.23	0.15	-0.08	Zm00001eb145870	LigB domain-containing protein	
	4,778,435	chr8	5,857,145	3.90 × 10−6	0.33	0.15	0.52			
	4,774,228	chr1	60,226,359	5.04 × 10−4	0.21	0.14	0.54	Zm00001eb016900	protein_coding	
	2,434,422	chr10	1.46E + 08	6.63 × 10−4	0.35	0.14	0.41			
	4,774,905	chr2	1.13E + 08	6.85 × 10−4	0.18	0.14	-0.01			
	2,469,044	chr2	27,774,525	7.5 × 10−4	0.19	0.14	0.16	Zm00001eb076570	Tubulin-folding cofactor B	
	100,020,060	chr5	63,352,656	7.77 × 10−4	0.13	0.14	-0.63			
	2,539,015	chr3	1.33E + 08	7.89 × 10−4	0.31	0.14	0.49	Zm00001eb137350	ATP10 protein expressed	
	4,775,714	chr7	4,773,701	8.87 × 10−4	0.15	0.13	0.59	Zm00001eb299670	GYF domain-containing protein	
	5,585,130	chr2	27,502,516	9.44 × 10−4	0.21	0.13	0.47	Zm00001eb112030	Sister chromatid cohesion protein DCC1	
SDR10	2,512,743	chr4	2,964,998	2.01 × 10−6	0.37	0.14	0.94	Zm00001eb165170	protein_coding	
	2,620,744	chr8	88,538,381	4.47 × 10−6	0.37	0.14	0.52	Zm00001eb346540	protein_coding	
	7,049,807	chr9	328,756	5.92 × 10−5	0.19	0.14	0.79	Zm00001eb371500	microtubule-associated protein 70 − 2 isoform X2	
	2,520,511	chr10	24,263,900	6.58 × 10−5	0.20	0.13	0.50	Zm00001eb410250	Kinesin-like protein KIN-14Q	
	4,582,343	chr10	96,965,850	7.56 × 10−5	0.25	0.13	-0.53	Zm00001eb419180	Receptor kinase	
	4,770,996	chr8	40,631,480	8.81 × 10−5	0.46	0.13	0.41	Zm00001eb340640	Polygalacturonase	
	4,590,802	chr7	54,060,484	9.64 × 10−5	0.25	0.12	0.43	Zm00001eb307490	Phytocyanin domain-containing protein	
SEC8 Striga emergence counts 8 weeks after sowing, SEC10 Striga emergence counts 10 weeks after sowing, SDR8 Striga damage rating 8 weeks after sowing, SDR10 Striga damage rating 10 weeks after sowing, MAF minor allele frequency, and R2 squared correlation coefficient, Gene ID gene identity

Fig. 4 Manhattan (right) and quantile-quantile (Q-Q) (left) plots showing SNPs associated with grain yield and Striga parameters in assessing 126 maize genotypes under Striga asiatica-infested conditions and 16,000 SNP markers. The plots were drawn using the compressed mixed linear model at P < 0.001. GY = grain yield, SEC8 = Striga emergence counts 8 weeks after sowing, SEC10 = Striga emergence counts 10 weeks after sowing, and SDR10 = Striga damage rating 10 weeks after sowing

Fig. 5 Manhattan (right) and quantile-quantile (Q-Q) (left) plots showing SNPs associated with grain yield and Striga parameters in assessing 126 maize genotypes under Striga hermonthica-infested conditions and 16,000 SNP markers. The plots were drawn using the compressed mixed linear model at P < 0.001. SEC8 = Striga emergence counts 8 weeks after sowing, SEC10 = Striga emergence counts 10 weeks after sowing, SDR8 = Striga damage rating 8 weeks after sowing, and SDR10 = Striga damage rating 10 weeks after sowing

The two MTAs recorded for GY under Sa-infested conditions were located on chromosomes 1 and 5 (Table 3). For SEC8, the 14 MTAs were identified on chromosomes 1 (one marker), 2 (one marker), 3 (one marker), 5 (three markers), 6 (two markers), 7 (one marker), 9 (one marker), and 10 (two markers) for Sh-infested conditions and chromosomes 10 (one marker) and 4 (one marker) for Sa-infested conditions. The 21 MTAs observed for SEC10 were on chromosomes 1 (two markers), 3 (two markers), 5 (three markers), 6 (one marker), 7 (one marker), 8 (one marker), and 10 (four markers) for Sh-infested conditions and on chromosomes 1 (two markers), 3 (two markers), 5 (one marker), seven (one marker), and 10 (ten markers) for Sa-infested conditions. SDR8 had 15 MTAs from chromosomes 1 (one marker), 2 (five markers), 3 (one marker), 4 (one marker), 5 (one marker), 6 (one marker), 7 (one marker), 8 (two markers), and 10 (two markers) in Sh-infested environment. MTAS for SDR10 were identified on chromosomes 4 (one marker), 8 (two markers), 9 (one marker), 10 (two markers), and 7 (one marker) under Sh-infested conditions and on chromosomes 1 (one marker), 2 (one marker), 3 (two markers), 5 (one marker), 6 (one marker), 7 (three markers), 8 (one marker), and 10 (one marker) for Sa-infested environment. Ten of the significant MTAs showed a negative allelic effect under Sa-infested conditions (Table 3), and six had a negative allelic effect under Sh-infested conditions (Table 4). One pleiotropic MTA was identified for SDR8 and SDR10 on chromosome 4 under Sh-infested conditions (Table 4).

Candidate genes associated with grain yield and Striga resistance traits

Based on SNP genome-wide association mapping, 33 and 14 candidate genes were identified for Sh and Sa resistance, respectively. The significant SNPs associated with GY under Sh were linked to six candidate genes Zm00001eb288770, Zm00001eb014030, Zm00001eb014020, Zm00001eb381210, Zm00001eb412710, and Zm00001eb219710. However, no candidate genes were associated with GY under Sa-infested conditions. Five of the markers linked to SEC8 under Sh-infested environment flanked the regions overlapping the candidate genes Zm00001eb313080, Zm00001eb253620, Zm00001eb272670, Zm00001eb163250, and Zm00001eb060090 on chromosome 1, 5, 6, and 10. In Sa-infested environment, only one candidate gene, Zm00001eb432870, was associated with SEC8. Eight and seven putative candidate genes were detected for SEC10 under Sh and Sa-infested conditions, respectively. These genes were located on chromosomes 5, 6, 7, 8, and 10, under Sa conditions, and 1, 5, 6, 7, 8, and 10 under Sh-conditions. Under Sa-conditions, the genes associated with SEC10 were located on chromosomes 1, 3, 5, 7, and 10. Twelve genes were identified for SDR8 under Sh-infested conditions and were located on chromosomes 1, 2, 3, 4, 5, 6, 7, and 8. The significant SNPs involved in SDR10 under Sh-infested conditions were localized within the following candidate genes: Zm00001eb165170, Zm00001eb346540, Zm00001eb371500, Zm00001eb410250, Zm00001eb419180, Zm00001eb340640, and Zm00001eb307490. Out of the 12 MTAs for SDR10 in Sa-infested environment, nine were found in regions covering the following genes: Zm00001eb319510, Zm00001eb336310, Zm00001eb114280, Zm00001eb329320, Zm00001eb156390, Zm00001eb409270, and Zm00001eb217880 on chromosomes 2, 3, 5, 7, and 10.

Discussion

Phenotyping

The parasitic weeds Sh and Sa have been recognized as major constraints to maize production in SSA. Striga resistance is a polygenic trait and is subject to genotype-environment interaction effects needing accurate selection strategies. Agronomic management options (e.g. crop rotation, intercropping, catch-cropping, and trap-cropping) are often insufficient. There is a need for integrated Striga management spearheaded by resistance breeding programmes. Breeding for Striga resistance is the most economical approach to minimise yield loss of major cereal crops caused by Striga infestation in SSA. Hence, genetic markers associated with Striga resistance need to be identified to guide the pace and precision of resistance breeding in maize [38, 39]. Integrating marker-assisted and conventional breeding techniques will accelerate breeding and genetic gain of maize for Striga resistance [40]. In this study, a GWAS was conducted based on a panel of genetically diverse maize populations comprising 98 inbred lines, 21 OPVs, and 11 top cross hybrids to discern genetic loci associated with resistance to the two dominant Striga species (Supplemental Table 1).

The phenotypic evaluation revealed significant variability for both Sa and Sh resistance among the tested maize genotypes, thus confirming the availability of relevant alleles for future breeding and genetic improvement (Table 1). The broad sense heritability observed for Striga damage under Sa-infested conditions was low (SDR8 = 0.01 and SDR10 = 0.11) compared to Sh resistance (SDR8 = 0.92 and SDR10 = 0.82) (Table 2). This suggests that the genotypic variance of the Sa resistance traits was low compared to the phenotypic variance, making direct selection difficult for Sa resistance. This confirms that Sh resistance has been the focus of more research that yielded higher genetic variation than Sa in maize. The high heritability of Sh resistance suggests that the genotypic variance was relatively higher supporting phenotypic selection. Hence, Sh traits can be relatively quickly selected (Table 2). The present results agree with Olakojo and Olaoye [41], who reported low heritability of some Striga parameters and high heritability on other associated traits when assessing maize agronomic traits for yield improvement and Sa tolerance. However, the high heritability values obtained could suggest the involvement of high nonadditive genetic variance and the confounding effects of the test environments, as Striga resistance traits are known as low heritable. This implies the need to test the population across multiple growing environments to reduce the confounding effect of the test environments and the genotype-by-environment interaction effect.

Population structure and linkage disequilibrium

Estimating population structure and within-group familial relatedness in maize genomic association studies is important to reduce the risks of false positives [42]. The Structure and Kinship analyses revealed population stratifications and admixtures, suggesting the need to use a robust statistical model in association analysis to control spurious MTAs [43]. The present analyses suggested eight sub-groups in the assessed tropical and sub-tropical maize population (Fig. 1), and the clustering was mainly based on genotype and geographic origin. This was expected, given that the genetic material used in the study was from diverse sources. The kinship heatmap showed different genetic relatedness in each group (Fig. 2). Hybrids and OPVs from IITA/Nigeria, NPGRC/South Africa, and CIMMYT/Zimbabwe form a distinct cluster from inbred lines from IITA/Nigeria and CIMMYT/Zimbabwe. Some clusters consisted of inbred lines from IITA (Clusters VI, VII, and VIII) or inbred lines from CIMMYT (Cluster V) (Figs. 1 and 2). The genetic distinctness of tropical inbred lines from CIMMYT lines suggests that crosses could be explored to achieve high heterosis and Striga resistance with broader adaptation.

Based on pairwise LD analysis involving 16,000 DArTSeq SNP markers, a low number of markers (9.78%) were in perfect LD (r2 = 1), which are more likely to be inherited across the whole genome. This reveals that this marker population was ideal for association mapping, as SNP markers in a strong LD provide redundant genotyping information during association analysis [44]. Furthermore, the rapid LD decay points to the high genetic diversity of the present panel and its usefulness for GWAS with a genetic mapping resolution (Fig. 3). The high genetic diversity is also associated with the genetic material that is composed of tropical and sub-tropical populations. The rapid LD decay is in agreement with previous studies in maize [45, 46] and other crops [47–49].

Marker-trait association for Striga resistance

For GWAS analysis, the study used CMLM model in GAPIT. The CMLM model of GAPIT was demonstrated to be more powerful and effective in association studies and can provide accurate predictions with less computing time [50]. Also, the method significantly reduces false positives [51]. The QQ plots revealed that the population structure and the phenotypic data were well distributed in both Striga-infested conditions, as the expected P-values were close to the observed P-values (Figs. 4 and 5). The plots revealed a good fit of the model. The kinship deviation at the top of the null hypothesis diagonal points to a great association between the markers and the traits (Figs. 4 and 5). The present study identified 22 MTAs for grain yield and Striga resistance traits under Sa infestations. The MTAs comprised of 2 MTAs conditioning GY, 2 MTAs conditioning SEC8, 7 MTAs conditioning SEC10, and 11 MTAs conditioning SDR10 (Table 3; Fig. 4). Fifty-seven MTAs were found for GY (7), SEC8 (12), SEC10 (14), SDR8 (17), and SDR10 (7) under Sh infestation (Table 4; Fig. 5). The marker 2,512,743 was pleiotropic for SDR8 and SDR10 for Sh resistance (Table 4) and could be useful for simultaneous selection. The SNPs were distributed across almost all chromosomes for all the assessed traits in the present study. The variable chromosome suggests that Striga resistance is a complex trait governed by polygenes. The significant markers identified in the current populations were located on chromosomes previously reported to harbour genes linked to Striga resistance traits including Striga emergence counts and Striga damage rating scores [13, 24, 52]. The significant markers for GY were located on chromosomes 1, 5, 6, 8, 9, and 10 and were associated with Sh resistance, and on chromosomes 1 and 5, they were associated with Sa resistance. These results corroborate previous studies that highlighted the importance of chromosome 10 [13], chromosome 6 [24], and chromosome 9 [52] for GY response under Sh conditions. In the present study, MTAs for Striga emergence counts were recorded on chromosome 5, previously reported to have a significant association with the same trait using the marker GRMZM2G018508 by [24], and GRMZM2G129543 and GRMZM5G823157 [52] in tropical maize lines. The same authors reported the ten chromosomes harbouring genes for Sh host damage in agreement with the present study. The markers reported in the present study for Sa resistance are novel except for the three SNPs reported by Pfunye, et al. [8] on chromosomes 5, 6, and 7 for total Striga plants emerged. There is a knowledge gap in QTL analysis associated with Sa resistance. The results revealed that different genomic regions are involved in Striga resistance (Tables 3 and 4). None of the genetic markers were associated with dual resistance of both Sa and Sh, suggesting that different genomics regions govern resistance for the two Striga species. Further association studies should be done using different populations to search for markers and genes controlling GY and Striga resistance traits under Sa and Sh infestation, as the two Striga species occur in tandem, especially in East Africa.

Seven of the significant markers had a negative allelic effect under Sh-infested conditions and four significant MTAs showed a negative allelic effect under Sa-infested conditions (Tables 3 and 4). This implies that these genetic markers are linked with Striga resistance proper for introgression into desirable parents. A negative allelic effect indicates that the allele is associated with a more resistant phenotype in a desirable direction. On the other hand, a positive effect indicates the opposite [53]. Stably expressing SNPs are useful for maize breeders to introduce target QTL in pipeline breeding materials. In this study, markers 4,772,065 and 2,447,891 are more stable for SDR8 under Sh-infested conditions.

Forty-one and 23 candidate functional genes were identified for Striga resistance under Sh and Sa-infested environments, respectively (Tables 3 and 4). Most detected genes are involved in various cellular and metabolic processes, plant defense, and cell development. The candidate gene Zm00001eb319510 is significantly associated with SDR10 under Sa-infested conditions and encodes Oxidoreductase-like domain-containing protein. Other candidate genes, such as Zm00001eb163250, Zm00001eb295600, Zm00001eb165170, and Zm00001eb165170 associated with SEC8, SDR8, and SDR10 under Sh-infested conditions, encode Clathrin heavy chain, Golgin candidate 5, and AAA + ATPase domain-containing protein respectively. The candidate gene Zm00001eb288770 associated with GY under Sh-conditions on chromosome 6 encodes for Coronatine-insensitive protein 1. Zm00001eb288770 is linked with jasmonate production, which regulates defense against insects and pathogens, wound healing, and pollen fertility [54]. The results show the involvement of Zm00001eb288770 in maize’s tolerance response to Striga. Another putative gene Zm00001eb253620, encodes for S-adenosylmethionine carrier 1 chloroplastic/mitochondrial that acts as a precursor of ethylene (a senescence inducer) and polyamines (antisenescence molecules) and controls the regulation of senescence in plants [55]. One of the main symptoms of Striga infestation is rapid senescence. This suggests that the gene Zm00001eb253620 is vital to unlocking Striga resistance breeding. Zm00001eb409850 encodes CBM20 domain-containing protein is implicated in both glycogen metabolism and autophagy [56]. The gene model Zm00001eb219710 associated with GY under Sh-infested conditions on chromosome 8 is associated with SEC10 under the same conditions on chromosome 7. Zm00001eb165170 and Zm00001eb112030 are linked with SDR8 and SDR10 under Sh-infested conditions.

Conclusion and recommendation

The present study is the first attempt to identify genomic regions associated with the dual resistance of maize to Sa and Sh. A population of 130 maize germplasm was used, comprising tropical and sub-tropical inbred lines, OPVs, and hybrids. The study identified 72 SNPs associated with grain yield and Striga resistance parameters under Sa and Sh-infested conditions. However, markers 4,772,065 and 2,447,891 are considered more stable for SDR8 under Sh-infested conditions. No marker was linked with dual resistance to Sh and Sa. Fifty significant markers were associated with Sh resistance and 22 were linked to Sa resistance. Significant SNPs were found flanking 47 protein-coding putative genes, of which one was associated with GY, while six, sixteen, eleven, and thirteen genes linked to SEC8, SEC10, SDR8, and SDR10 were associated with Sh and Sa resistance. Some of the detected genes reportedly conditioned insect and pathogen resistance, plant cell development, variable senescence, and pollen fertility. The gene Zm00001eb219710 was pleiotropic, and conditioned GY and SEC10, while Zm00001eb165170 affected SDR8 and SDR10, and Zm00001eb112030 conditioned SDR8 and SDR10 associated with Sh resistance. The candidate markers may facilitate simultaneous selection for Sh and Sa resistance and grain yield in maize after further evaluation and introgression in breeding pipelines. Based on the current results, the two Striga species have different genetic compositions. The tested maize genotypes were resistant to Sh but not necessarily resistant to Sa. There is a need to collect separate data sets from multiple studies on the two species. This will enable to document the economic importance and genetic basis of the two Striga species to guide future research and development. We recommend resistance breeding of maize on separate genetic backgrounds for Sa and Sh resistance. This will allow for a much clearer perspective on the economic importance of each of the two major Striga species. However, gene pyramiding on a desirable genetic background of maize could harness Striga resistance against the two species.

Supplementary Information

Supplementary Material 1.

Acknowledgements

The authors acknowledge the African Centre for Crop Improvement (ACCI) and the Organisation for Women in Science for the Developing World (OWSD) for funding this research work.

Additional information

The Striga asiatica and S. hermonthica seeds used in this study were separately packed and stored in the Soil Science Laboratory of TARI Tumbi in Tanzania for further use, and are available on demand.

Authors’ contributions

Author contributions: E.N.D: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Resources; Software; Writing – original draft; Writing – review & editing. H.S: Conceptualization; Funding acquisition; Investigation; Methodology; Project administration; Resources; Supervision; Validation; Visualization; Writing – review & editing. A.S: Conceptualization; Funding acquisition; Investigation; Methodology; Resources; Validation; Visualization; Writing – review & editing.

Funding

This work was conducted through the funding of African Centre for Crop Improvement (ACCI) and the Organisation for Women in Science for the Developing World (OWSD).

Availability of data and materials

The datasets generated during/or analyzed during the current study are publicly available and deposited in the Mendeley data archive. These data can be found here: https://data.mendeley.com/datasets/fd2mf3zbhw/1.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Prasanna BM, Palacios-Rojas N, Hossain F, Muthusamy V, Menkir A, Dhliwayo T, Ndhlela T, San Vicente F, Nair SK, Vivek BS et al. Molecular breeding for nutritionally enriched maize: status and prospects. Front Genet. 2020;10:1–16. 10.3389/fgene.2019.01392.
2. Lee D, Davenport F, Shukla S, Husak G, Funk C, Harrison L, McNally A, Rowland J, Budde M, Verdin J. Maize yield forecasts for sub-Saharan Africa using earth observation data and machine learning. Global Food Secur. 2022;33:1–11. 10.1016/j.gfs.2022.100643.
3. Lobulu J, Shimelis H, Laing M, Mushongi AA. Maize production constraints, traits preference and current Striga control options in western Tanzania: farmers’ consultation and implications for breeding. Acta Agriculturae Scand Sect B-Soil Plant Sci. 2019;69:734–46. 10.1080/09064710.2019.1652680.
4. David OG, Ayangbenro AS, Odhiambo JJ, Babalola OO. Striga hermonthica: a highly destructive pathogen in maize production. Environ Challenges. 2022;8:1–9.
5. Badu-Apraku B Akinwale R Cultivar evaluation and trait analysis of tropical early maturing maize under Striga-infested and Striga-free environments Field Crops Res 2011 121 186 194 10.1016/j.fcr.2010.12.011
Badu-Apraku B, Akinwale R. Cultivar evaluation and trait analysis of tropical early maturing maize under Striga-infested and Striga-free environments. Field Crops Res. 2011;121:186–94.
6. Menkir A Makumbi D Franco J Assessment of reaction patterns of hybrids to Striga hermonthica (Del.) Benth. Under artificial infestation in Kenya and Nigeria Crop Sci 2012 52 2528 2537 10.2135/cropsci2012.05.0307
Menkir A, Makumbi D, Franco J. Assessment of reaction patterns of hybrids to Striga hermonthica (Del.) Benth. Under artificial infestation in Kenya and Nigeria. Crop Sci. 2012;52:2528–37.
7. Menkir A Kling JG Response to recurrent selection for resistance to Striga hermonthica (Del.) Benth in a tropical maize population Crop Sci 2006 47 674 84 10.2135/cropsci2006.07.0494
Menkir A, Kling JG. Response to recurrent selection for resistance to Striga hermonthica (Del.) Benth in a tropical maize population. Crop Sci. 2006;47:674–84. 10.2135/cropsci2006.07.0494.
8. Pfunye A, Rwafa R, Mabasa S, Gasura E. Genome-wide association studies for Striga asiatica resistance in tropical maize. Int J Genomics. 2021:1–8. 10.1155/2021/9979146.
9. Kavuluko J Kibe M Sugut I Kibet W Masanga J Mutinda S Wamalwa M Magomere T Odeny D Runo S GWAS provides biological insights into mechanisms of the parasitic plant (Striga) resistance in sorghum BMC Plant Biol 2021 21 1 15 10.1186/s12870-021-03155-7 33386069
Kavuluko J, Kibe M, Sugut I, Kibet W, Masanga J, Mutinda S, Wamalwa M, Magomere T, Odeny D, Runo S. GWAS provides biological insights into mechanisms of the parasitic plant (Striga) resistance in sorghum. BMC Plant Biol. 2021;21:1–15. 10.1186/s12870-021-03155-7.33386069
10. Mrema E, Shimelis H, Laing M, Mwadzingeni L. Integrated management of’ Striga hermonthica and S. Asiatica in sorghum: a review. Aust J Crop Sci. 2020;14:36–45.
11. Gasura E Setimela P Mabasa S Rwafa R Kageler S Nyakurwa C Response of IITA maize inbred lines bred for Striga hermonthica resistance to Striga asiatica and associated resistance mechanisms in southern Africa Euphytica 2019 215 1 15 10.1007/s10681-019-2467-5
Gasura E, Setimela P, Mabasa S, Rwafa R, Kageler S, Nyakurwa C. Response of IITA maize inbred lines bred for Striga hermonthica resistance to Striga asiatica and associated resistance mechanisms in southern Africa. Euphytica. 2019;215:1–15. 10.1007/s10681-019-2467-5.
12. Rich PJ Ejeta G Towards effective resistance to Striga in African maize Plant Signal Behav 2008 3 618 621 10.4161/psb.3.9.5750 19513251
Rich PJ, Ejeta G. Towards effective resistance to Striga in African maize. Plant Signal Behav. 2008;3:618–21. 10.4161/psb.3.9.5750.19513251
13. Adewale SA Badu-Apraku B Akinwale RO Paterne AA Gedil M Garcia-Oliveira AL Genome-wide association study of Striga resistance in early maturing white tropical maize inbred lines BMC Plant Biol 2020 20 1 20 10.1186/s12870-020-02360-0 31898482
Adewale SA, Badu-Apraku B, Akinwale RO, Paterne AA, Gedil M, Garcia-Oliveira AL. Genome-wide association study of Striga resistance in early maturing white tropical maize inbred lines. BMC Plant Biol. 2020;20:1–20. 10.1186/s12870-020-02360-0.31898482
14. Dossa EN, Shimelis H, Shayanowako AIT, Laing MD. A meta-analysis of the effects of Striga control methods on maize, sorghum, and major millets production in sub‐Saharan Africa. Crop Sci. 2023:1–20. 10.1002/csc2.20889.
15. Shayanowako AIT Shimelis H Laing MD Mwadzingeni L Striga resistance and compatibility of maize genotypes to a biocontrol agent, Fusarium oxysporum f.sp.strigea J Crop Improv 2020 34 437 454 10.1080/15427528.2020.1728599
Shayanowako AIT, Shimelis H, Laing MD, Mwadzingeni L. Striga resistance and compatibility of maize genotypes to a biocontrol agent, Fusarium oxysporum f.sp.strigea. J Crop Improv. 2020;34:437–54. 10.1080/15427528.2020.1728599.
16. De Groote H Wangare L Kanampiu F Odendo M Diallo A Karaya H Friesen D The potential of a herbicide resistant maize technology for Striga control in Africa Agric Syst 2008 97 83 94 10.1016/j.agsy.2007.12.003
De Groote H, Wangare L, Kanampiu F, Odendo M, Diallo A, Karaya H, Friesen D. The potential of a herbicide resistant maize technology for Striga control in Africa. Agric Syst. 2008;97:83–94. 10.1016/j.agsy.2007.12.003.
17. Bhat J Salgotra R Dar M Phenomics: a challenge for crop improvement in genomic era Mol Plant Breed 2015 6 1 11
Bhat J, Salgotra R, Dar M. Phenomics: a challenge for crop improvement in genomic era. Mol Plant Breed. 2015;6:1–11.
18. Zhang X Guan Z Li Z Liu P Ma L Zhang Y Pan L He S Zhang Y Li P A combination of linkage mapping and GWAS brings new elements on the genetic basis of yield-related traits in maize across multiple environments Theor Appl Genet 2020 133 2881 2895 10.1007/s00122-020-03639-4 32594266
Zhang X, Guan Z, Li Z, Liu P, Ma L, Zhang Y, Pan L, He S, Zhang Y, Li P, et al. A combination of linkage mapping and GWAS brings new elements on the genetic basis of yield-related traits in maize across multiple environments. Theor Appl Genet. 2020;133:2881–95. 10.1007/s00122-020-03639-4.32594266
19. Sukumaran S Dreisigacker S Lopes M Chavez P Reynolds MP Genome-wide association study for grain yield and related traits in an elite spring wheat population grown in temperate irrigated environments Theor Appl Genet 2015 128 353 363 10.1007/s00122-014-2435-3 25490985
Sukumaran S, Dreisigacker S, Lopes M, Chavez P, Reynolds MP. Genome-wide association study for grain yield and related traits in an elite spring wheat population grown in temperate irrigated environments. Theor Appl Genet. 2015;128:353–63.25490985
20. Garcia M Eckermann P Haefele S Satija S Sznajder B Timmins A Baumann U Wolters P Mather DE Fleury D Genome-wide association mapping of grain yield in a diverse collection of spring wheat (Triticum aestivum L.) evaluated in southern Australia PLoS ONE 2019 14 1 18 10.1371/journal.pone.0211730
Garcia M, Eckermann P, Haefele S, Satija S, Sznajder B, Timmins A, Baumann U, Wolters P, Mather DE, Fleury D. Genome-wide association mapping of grain yield in a diverse collection of spring wheat (Triticum aestivum L.) evaluated in southern Australia. PLoS ONE. 2019;14:1–18.
21. Xu Y Li P Yang Z Xu C Genetic mapping of quantitative trait loci in crops Crop J 2017 5 175 184 10.1016/j.cj.2016.06.003
Xu Y, Li P, Yang Z, Xu C. Genetic mapping of quantitative trait loci in crops. Crop J. 2017;5:175–84.
22. Amusan IO. Mechanisms and quantitative trait loci for Striga hermonthica resistance in maize (Zea mays L.) inbred line. Purdue University. 2010: pp. 1–188.
23. Badu-Apraku B Adewale S Paterne AA Gedil M Toyinbo J Asiedu R Identification of QTLs for grain yield and other traits in tropical maize under Striga infestation PLoS ONE 2020 15 1 20 10.1371/journal.pone.0239205
Badu-Apraku B, Adewale S, Paterne AA, Gedil M, Toyinbo J, Asiedu R. Identification of QTLs for grain yield and other traits in tropical maize under Striga infestation. PLoS ONE. 2020;15:1–20. 10.1371/journal.pone.0239205.
24. Gowda M, Makumbi D, Das B, Nyaga C, Kosgei T, Crossa J, Beyene Y, Montesinos-Lopez OA, Olsen MS, Prasanna BM. Genetic dissection of Striga hermonthica (Del.) Benth. Resistance via genome-wide association and genomic prediction in tropical maize germplasm. Theor Appl Genet. 2021;134:941–58. 10.1007/s00122-020-03744-4.
25. Okunlola G, Badu-Apraku B, Ariyo O, Agre P, Offernedo Q, Ayo-Vaughan M. Genome-wide association studies of Striga resistance in extra-early maturing quality protein maize inbred lines. G3 (Bethesda). 2023;G3(2):1–11. 10.1093/g3journal/jkac237.
26. Arrones A Vilanova S Plazas M Mangino G Pascual L Díez MJ Prohens J Gramazio P The dawn of the age of multi-parent MAGIC populations in plant breeding: novel powerful next-generation resources for genetic analysis and selection of recombinant elite material Biology 2020 9 229 10.3390/biology9080229 32824319
Arrones A, Vilanova S, Plazas M, Mangino G, Pascual L, Díez MJ, Prohens J, Gramazio P. The dawn of the age of multi-parent MAGIC populations in plant breeding: novel powerful next-generation resources for genetic analysis and selection of recombinant elite material. Biology. 2020;9: 229.32824319
27. Dossa EN, Shimelis H, Shayanowako AIT, Laing MD. Screening tropical and sub-tropical maize germplasm for resistance to Striga hermonthica and S. asiatica and yield-related traits. Euphytica. 2024;220:1–23. 10.1007/s10681-024-03309-2.
28. Kim SK Akintunde AY Walker P Responses of maize, sorghum and millet host plants to infestation by Striga hermonthica Crop Prot 1994 13 582 90 10.1016/0261-2194(94)90003-5
Kim SK, Akintunde AY, Walker P. Responses of maize, sorghum and millet host plants to infestation by Striga hermonthica. Crop Prot. 1994;13:582–90. 10.1016/0261-2194(94)90003-5.
29. Kilian A, Wenzl P, Huttner E, Carling J, Xia L, Blois H, Caig V, Heller-Uszynska K, Jaccoud D, Hopper C. Diversity arrays technology: a generic genome profiling technology on open platforms. Data production and analysis in population genomics: methods and protocols 2012;67–89.
30. Pritchard JK Stephens M Donnelly P Inference of population structure using multilocus genotype data Genetics 2000 1552 945 59 10.1093/genetics/155.2.945
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;1552:945–59.
31. Evanno G egnaut S Goudet J Detecting the number of clusters of individuals using the software structure: a simulation study Mol Ecol 2005 14 2611 2620 10.1111/j.1365-294X.2005.02553.x 15969739
Evanno G, egnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14:2611–20. 10.1111/j.1365-294X.2005.02553.x.15969739
32. Bradbury PJ Zhang Z Kroon DE Casstevens TM Ramdoss Y Buckler ES TASSEL: software for association mapping of complex traits in diverse samples Bioinformatics 2007 23 2633 2635 10.1093/bioinformatics/btm308 17586829
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5. 10.1093/bioinformatics/btm308.17586829
33. Remington DL Thornsberry JM Matsuoka Y Wilson LM Whitt SR Doebley J Kresovich S Goodman MM Buckler ES IV Structure of linkage disequilibrium and phenotypic associations in the maize genome Proc Natl Acad Sci 2001 98 11479 11484 10.1073/pnas.201394398 11562485
Remington DL, Thornsberry JM, Matsuoka Y, Wilson LM, Whitt SR, Doebley J, Kresovich S, Goodman MM, Buckler ES IV. Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc Natl Acad Sci. 2001;98:11479–84.11562485
34. Gimase JM, Thagana WM, Omondi CO, Cheserek JJ, Gichimu BM, Gichuru EK, Ziyomo C, Sneller CH. Genome-wide Association Study identify the genetic loci conferring resistance to Coffee Berry Disease (Colletotrichum Kahawae) in Coffea arabica var. Rume Sudan. Euphytica. 2020;216:216. 10.1007/s10681-020-02621-x.
35. Nkhata W Shimelis H Melis R Chirwa R Mzengeza T Mathew I Shayanowako A Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean PLoS ONE 2021 16 1 24 10.1371/journal.pone.0250729
Nkhata W, Shimelis H, Melis R, Chirwa R, Mzengeza T, Mathew I, Shayanowako A. Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean. PLoS ONE. 2021;16:1–24. 10.1371/journal.pone.0250729.
36. Nyine M Uwimana B Akech V Brown A Ortiz R Doležel J Lorenzen J Swennen R Association genetics of bunch weight and its component traits in east African highland banana (Musa Spp. AAA group) Theor Appl Genet 2019 132 3295 3308 10.1007/s00122-019-03425-x 31529270
Nyine M, Uwimana B, Akech V, Brown A, Ortiz R, Doležel J, Lorenzen J, Swennen R. Association genetics of bunch weight and its component traits in east African highland banana (Musa Spp. AAA group). Theor Appl Genet. 2019;132:3295–308.31529270
37. Bhadmus OA Badu-Apraku B Adeyemo OA Agre PA Queen ON Ogunkanmi AL Genome-wide association analysis reveals genetic architecture and candidate genes associated with grain yield and other traits under low soil nitrogen in early-maturing white quality protein maize inbred lines Genes (Basel) 2022 13 1 20 10.3390/genes13050826
Bhadmus OA, Badu-Apraku B, Adeyemo OA, Agre PA, Queen ON, Ogunkanmi AL. Genome-wide association analysis reveals genetic architecture and candidate genes associated with grain yield and other traits under low soil nitrogen in early-maturing white quality protein maize inbred lines. Genes (Basel). 2022;13:1–20. 10.3390/genes13050826.
38. Longmei N Gill GK Zaidi PH Kumar R Nair SK Hindu V Vinayan MT Vikal Y Genome-wide association mapping for heat tolerance in sub-tropical maize BMC Genomics 2021 22 1 14 10.1186/s12864-021-07463-y 33388042
Longmei N, Gill GK, Zaidi PH, Kumar R, Nair SK, Hindu V, Vinayan MT, Vikal Y. Genome-wide association mapping for heat tolerance in sub-tropical maize. BMC Genomics. 2021;22:1–14. 10.1186/s12864-021-07463-y.33388042
39. Chen J Xu W Velten J Xin Z Stout J Characterization of maize inbred lines for drought and heat tolerance J Soil Water Conserv 2012 67 354 364 10.2489/jswc.67.5.354
Chen J, Xu W, Velten J, Xin Z, Stout J. Characterization of maize inbred lines for drought and heat tolerance. J Soil Water Conserv. 2012;67:354–64.
40. Badu-Apraku B Adewale S Paterne A Gedil M Asiedu R Identification of QTLs controlling resistance/tolerance to Striga hermonthica in an extra-early maturing yellow maize population Agronomy-Basel 2020 10 1 18 10.3390/agronomy10081168
Badu-Apraku B, Adewale S, Paterne A, Gedil M, Asiedu R. Identification of QTLs controlling resistance/tolerance to Striga hermonthica in an extra-early maturing yellow maize population. Agronomy-Basel. 2020;10:1–18. 10.3390/agronomy10081168.
41. Olakojo SA, Olaoye G. Correlation and heritability estimates of maize agronomic traits for yield improvement and Striga asiatica (L.) kuntze tolerance. Afr J Plant Sci. 2011;5(6):365–9.
42. Yu J Buckler ES Genetic association mapping and genome organization of maize Curr Opin Biotechnol 2006 17 155 160 10.1016/j.copbio.2006.02.003 16504497
Yu J, Buckler ES. Genetic association mapping and genome organization of maize. Curr Opin Biotechnol. 2006;17:155–60. 10.1016/j.copbio.2006.02.003.16504497
43. Mekonnen T, Sneller CH, Haileselassie T, Ziyomo C, Abeyo BG, Goodwin SB, Lule D, Tesfaye K. Genome-wide association study reveals novel genetic loci for quantitative resistance to Septoria Tritici Blotch in Wheat (Triticum aestivum L). Front Plant Sci. 2021;12(1–20). 10.3389/fpls.2021.671323.
44. Mohammadi M Xavier A Beckett T Beyer S Chen L Chikssa H Cross V Moreira FF French E Gaire R Identification, deployment, and transferability of quantitative trait loci from genome-wide association studies in plants Curr Plant Biology 2020 24 1 12 10.1016/j.cpb.2020.100145
Mohammadi M, Xavier A, Beckett T, Beyer S, Chen L, Chikssa H, Cross V, Moreira FF, French E, Gaire R. Identification, deployment, and transferability of quantitative trait loci from genome-wide association studies in plants. Curr Plant Biology. 2020;24:1–12.
45. de Faria SV Zuffo LT Rezende WM Caixeta DG Pereira HD Azevedo CF DeLima RO Phenotypic and molecular characterization of a set of tropical maize inbred lines from a public breeding program in Brazil BMC Genomics 2022 23 1 17 10.1186/s12864-021-08127-7 34979896
de Faria SV, Zuffo LT, Rezende WM, Caixeta DG, Pereira HD, Azevedo CF, DeLima RO. Phenotypic and molecular characterization of a set of tropical maize inbred lines from a public breeding program in Brazil. BMC Genomics. 2022;23:1–17. 10.1186/s12864-021-08127-7.34979896
46. Nyaga C Gowda M Beyene Y Muriithi WT Makumbi D Olsen MS Suresh LM Bright JM Das B Prasanna BM Genome-wide analyses and prediction of resistance to MLN in large tropical maize germplasm Genes (Basel) 2019 11 16 10.3390/genes11010016 31877962
Nyaga C, Gowda M, Beyene Y, Muriithi WT, Makumbi D, Olsen MS, Suresh LM, Bright JM, Das B, Prasanna BM. Genome-wide analyses and prediction of resistance to MLN in large tropical maize germplasm. Genes (Basel). 2019;11:16. 10.3390/genes11010016.31877962
47. Leiser WL Rattunde HF Weltzien E Cisse N Abdou M Diallo A Toure AO Magalhaes JV Haussmann BI Two in one sweep: aluminum tolerance and grain yield in P-limited soils are associated to the same genomic region in west African sorghum BMC Plant Biol. 2014 14 14 10.1186/s12870-014-0206-6 24405939
Leiser WL, Rattunde HF, Weltzien E, Cisse N, Abdou M, Diallo A, Toure AO, Magalhaes JV, Haussmann BI. Two in one sweep: aluminum tolerance and grain yield in P-limited soils are associated to the same genomic region in west African sorghum. BMC Plant Biol. 2014;14:14. 10.1186/s12870-014-0206-6.24405939
48. Odilbekov F Armoniene R Koc A Svensson J Chawade A GWAS-assisted genomic prediction to predict resistance to Septoria Tritici Blotch in nordic winter wheat at seedling stage Front Genet 2019 10 1 10 10.3389/fgene.2019.01224 30804975
Odilbekov F, Armoniene R, Koc A, Svensson J, Chawade A. GWAS-assisted genomic prediction to predict resistance to Septoria Tritici Blotch in nordic winter wheat at seedling stage. Front Genet. 2019;10:1–10. 10.3389/fgene.2019.01224.30804975
49. Pasam RK Sharma R Malosetti M Eeuwijk FAv Haseneyer G Kilian B Graner A Genome-wide association studies for agronomical traits in a world wide spring barley collection BMC Plant Biol 2012 121 1 22
Pasam RK, Sharma R, Malosetti M, Eeuwijk FAv, Haseneyer G, Kilian B, Graner A. Genome-wide association studies for agronomical traits in a world wide spring barley collection. BMC Plant Biol. 2012;121:1–22.
50. Chang F Guo C Sun F Zhang J Wang Z Kong J He Q Sharmin RA Zhao T Genome-wide association studies for dynamic plant height and number of nodes on the main stem in summer sowing soybeans Front Plant Sci 2018 9 1 13 10.3389/fpls.2018.01184 29410674
Chang F, Guo C, Sun F, Zhang J, Wang Z, Kong J, He Q, Sharmin RA, Zhao T. Genome-wide association studies for dynamic plant height and number of nodes on the main stem in summer sowing soybeans. Front Plant Sci. 2018;9:1–13. 10.3389/fpls.2018.01184.29410674
51. Zhang J Song Q Cregan PB Nelson RL Wang X Wu J Jiang G-L Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm BMC Genomics 2015 16 1 11 25553907
Zhang J, Song Q, Cregan PB, Nelson RL, Wang X, Wu J, Jiang G-L. Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm. BMC Genomics. 2015;16:1–11.25553907
52. Stanley AE, Menkir A, Ifie B, Paterne AA, Unachukwu NN, Meseka S, Mengesha WA, Bossey B, Kwadwo O, Tongoona PB, et al. Association analysis for resistance to Striga hermonthica in diverse tropical maize inbred lines. Sci Rep. 2021;11:24193. 10.1038/s41598-021-03566-4.
53. Bossa-Castro AM Tekete C Raghavan C Delorean EE Dereeper A Dagno K Koita O Mosquera G Leung H Verdier V Allelic variation for broad-spectrum resistance and susceptibility to bacterial pathogens identified in a rice MAGIC population Plant Biotechnol J 2018 16 1559 68 10.1111/pbi.12895 29406604
Bossa-Castro AM, Tekete C, Raghavan C, Delorean EE, Dereeper A, Dagno K, Koita O, Mosquera G, Leung H, Verdier V, et al. Allelic variation for broad-spectrum resistance and susceptibility to bacterial pathogens identified in a rice MAGIC population. Plant Biotechnol J. 2018;16:1559–68. 10.1111/pbi.12895.29406604
54. Xie D-X Feys BF James S Nieto-Rostro M Turner JG COI1: an Arabidopsis gene required for jasmonate-regulated defense and fertility Science 1998 280 1091 1094 10.1126/science.280.5366.1091 9582125
Xie D-X, Feys BF, James S, Nieto-Rostro M, Turner JG. COI1: an Arabidopsis gene required for jasmonate-regulated defense and fertility. Science. 1998;280:1091–4.9582125
55. Kozbial PZ Mushegian AR Natural history of S-adenosylmethionine-binding proteins BMC Struct Biol 2005 5 1 19 10.1186/1472-6807-5-19 15663787
Kozbial PZ, Mushegian AR. Natural history of S-adenosylmethionine-binding proteins. BMC Struct Biol. 2005;5:1–19. 10.1186/1472-6807-5-19.15663787
56. Itakura AK Chan KX Atkinson N Pallesen L Wang L Reeves G Patena W Caspari O Roth R Goodenough U A rubisco-binding protein is required for normal pyrenoid number and starch sheath morphology in Chlamydomonas reinhardtii Proc Natl Acad Sci United States Am 2019 116 18445 18454 10.1073/pnas.1904587116
Itakura AK, Chan KX, Atkinson N, Pallesen L, Wang L, Reeves G, Patena W, Caspari O, Roth R, Goodenough U, et al. A rubisco-binding protein is required for normal pyrenoid number and starch sheath morphology in Chlamydomonas reinhardtii. Proc Natl Acad Sci United States Am. 2019;116:18445–54. 10.1073/pnas.1904587116.
