Abstract

Sorghum (Sorghum bicolor (L.) Moench) is an important food crop in semi-arid tropics. The crop grain yield ranges from 0.5 t/ha to 0.8 t/ha compared to potential yields of 10 t/ha. The African stem borer Busseola fusca Fuller (Noctuidae) and the spotted stem borer Chilo partellus Swinhoe (Crambidae), are among the most economically important insect pests of sorghum. The two borers can cause 15% - 80% grain yield loss in sorghum. Mapping of QTLs associated with resistance traits to the two stem borers is important towards marker-assisted breeding. The objective of this study was to map QTLs associated with resistance traits to B. fusca and C. partellus in sorghum. 243 F9:10 sorghum RILs derived from ICSV 745 (S) and PB 15520-1 (R) were selected for the study with 4,955 SNP markers. The RILs were evaluated in three sites. Data was collected on leaf feeding, deadheart, exit holes, stem tunnels, leaf toughness, seedling vigour, bloom waxiness, and leaf glossiness. ANOVA for all the traits was done using Genstat statistical software. Insect damage traits and morphological traits were correlated using Pearson’s correlation coefficients. Genetic mapping was done using JoinMap 4 software, while QTL analysis was done using PLABQTL software. A likelihood odds ratio (LOD) score of 3.0 was used to declare linkage. Joint analyses across borer species and sites revealed 4 QTLs controlling deadheart formation; 6 controlling leaf feeding damage; 5 controlling exit holes and stem tunneling damages; 2 controlling bloom waxiness, leaf glossiness, and seedling vigour; 4 conditioning trichome density; and 6 conditioning leaf toughness. Joint analyses for B. fusca and C. partellus further revealed that marker CS132-2 colocalised for leaf toughness and stem tunneling traits on QTLs 1 and 2, respectively; thus, the two traits can be improved using the same linked marker. This study recommended further studies to identify gene(s) underlying the mapped QTLs.

1. Introduction

Sorghum (Sorghum bicolor (L.) Moench) is an important food crop in drought prone areas in the tropics. The crop is cultivated by more than 500 million resource-challenged smallholder farmers, mostly women [1]. The crop is an important industrial crop in East Africa and has untapped potential in bioenergy production [2]. The African stem borer, Busseola fusca Fuller (Noctuidae) and the spotted stem borer Chilo partellus Swinhoe (Crambidae), are among the most economically important insect pests of sorghum and maize in Eastern and South Africa [3]. The two stem borer species cause leaf feeding, deadheart formation, exit holes, and stem tunneling damages in cereals [3]. Stem borers are associated with grain yield loss of 15% - 80% depending on borer species population and variety phenological stage at the time of attack [4]. Management approaches such as cultural practices and use of synthetic chemical pesticides have yielded little success in the management of the stem borers. Host plant resistance can reduce sorghum grain yield losses to stem borers and thus improve food security among rural households [5]. Breeding for host plant resistance is part of integrated pest management that can contribute in an economic way in the management of insect pests in indigenous cereals cultivated by resource-challenged smallholder farmers [5]. Breeding for stem borer resistance has been slow partly due to the inadequate understanding of inheritance of traits conditioning resistance to the pests. Stem borer resistance is a polygenic trait controlled by many genes of small effects [6]. Grain sorghum with dual resistance to B. fusca and C. partellus has been reported [7]. There is no information about genomic regions associated with dual resistance to B. fusca and C. partellus in grain sorghum. Traits such as leaf toughness, trichomes, bloom waxiness, leaf glossiness, and seedling vigour have been reported to condition resistance to borer and foliar insect resistance in cereals [5, 811].

Identification and mapping of quantitative trait loci (QTLs) associated with sorghum resistance to stem borers could enhance the efficiency and effectiveness in marker-assisted breeding [12]. QTL mapping enhances the biological understanding of inheritance of quantitative traits, and the markers identified can be used to select for a complex trait [13]. QTL mapping using recombinant inbred lines (RILs) increases the power of QTL detection compared to the F2:3 population because of complete homozygosity at QTLs and marker loci [14]. Marker-assisted selection for stem borer resistance can enhance breeding for B. fusca and C. partellus since traditional breeding has been unsuccessful. Genetic linkage maps are essential for localization of genes conferring resistance/tolerance to stem borer damage in sorghum [15]. Single nucleotide polymorphism (SNP) markers have become widely accepted as a tool for understanding complex genetic traits and evolution [16]. SNPs have been found to be the most efficient genetic marker for gene identification since they are codominant, highly polymorphic, and have good reproducibility [16]. SNPs represent the finest resolution of a DNA sequence [17]. The genotyping-by-sequencing approach has been employed in whole-genome sequencing to discover SNPs for mapping studies in crop plants [18, 19]. The objective of this study was to identify and characterize QTLs associated with resistance traits to B. fusca and C. partellus in sorghum.

2. Materials and Methods

2.1. Development of Mapping Population

243 F9:10 recombinant inbred lines (RILs) derived from a cross between the stem borer-susceptible cultivar ICSV 745 and stem borer-resistant PB 15520-1 were selected for the mapping study. The progenies and their parents had been developed following a single seed descent approach in ICRISAT, Patancheru, India [20]. The two parents were crossed, and the resulting F1 seeds were advanced to F2 by selfing of single F1 plants. The F2 seeds were selfed and the resulting F3 population grown in progeny rows [20]. Single sorghum plants in each of the progeny row were selfed and the process continued up to the F9 generation [20]. Seeds from the F9 sorghum plants of each row were bulked to produce the 243 F9:10 RILs used in this study.

2.2. Experimental Design, Site, and Source of Stem Borer Larvae

The 243 RILs along with their parents were imported from the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru, India to Kenya. The sorghum plants were phenotyped at Embu and Kabete for resistance against B. fusca for one season in each environment and at Kiboko for resistance against C. partellus for two seasons thus totaling to four environments. Each experiment in each site was laid out in a 0 alpha-lattice design consisting of twenty-five plots in ten blocks, replicated twice. Each plot consisted of 2 m rows with plants spaced at inter- and intra-rows, respectively. First-instar neonates of the two borer species were obtained from the International Centre of Insect Physiology and Ecology (ICIPE), Nairobi, Kenya. At 30 days after sowing, five plants in each row were tagged and artificially infested with five larvae of the respective stem borer species using a camel hair brush as described by Singh et al. [21] in all the three sites.

2.3. Phenotypic Data Collection and Data Analysis

Data on stem borer leaf feeding damage, deadheart incidence, number of exit holes, and stem tunnel length were scored as described by Muturi et al. [7]. Morphological traits measured included leaf toughness, seedling vigour, bloom waxiness, leaf glossiness, and total grain yield as described by Kumar et al. [22]. Data were subjected to analysis of variance using the residual maximum likelihood model (ReML) in Genstat Version 14 statistical package. The predicted means for each genotype were estimated with genotypes as fixed and reps as random effects in individual and across environment analyses. Phenotypic and genetic correlation coefficients were calculated from adjusted entry means across environments for each parameter. Estimates of variance components, including genotypic variance (), interaction (), and residual () were calculated by equating the mean squares to their expected values as described by Shimelis and Shiringani [23]. Heritability () was estimated using Meta-R software version 6 using , where = genotypic variance, = genotype by environment variance, = error variance, = number of replications, and env = number of environments as described by Hallauer and Miranda, [24]. Direct and indirect effect analyses were conducted using path analysis to study interrelationship among resistance parameters and their relationship to grain yield reduction [25]. Genotypic and phenotypic correlation coefficients were calculated following the method described by Holland [26].

2.4. DNA Extraction and Quantification

Leaf tissue from the parental lines and the 243 RILs were harvested from 10-day-old sorghum seedlings. DNA was extracted according to [27]. The quality of DNA in each sample was checked using 0.8% agarose gel stained with ethidium bromide. Each well of the agarose gel was loaded with 5 μl of the sample, and the gel was allowed to run at 100 V for 5 minutes. After electrophoresis, DNA banding patterns on the gel were visualized under UV light. A smear of DNA indicated poor quality, and the DNA was re-extracted, whereas a clear band indicated good quality DNA. The quantity of DNA in each sample was assessed using a fluorescence spectrophotometer (Switzerland) by staining the DNA with PicoGreen™ (1/200 dilution). The DNA concentration of each sample was calculated and then normalized to 2.5 ng/μl for the PCR.

2.5. Scoring of Sequenced Products and Construction of Linkage Map

The SNP markers screened on the RILs were scored as follows: A = homozygote carrying allele from female parent (ICSV 745); B = homozygote carrying allele from male parent (PB 15520-1); and - = missing data for an individual at a locus. The genotypic data was used to construct a genetic linkage map, which spanned 4,692.4 cM, with a total of 4,955 SNP markers distributed into the 10 sorghum linkage groups (Figure 1).

Segregation of each marker was tested with a chi-square goodness of fit test to the expected Mendelian segregation ratio (1 : 1) of parental configuration. The names of the markers were coded for ease of analysis (Supplementary Material Table 14). The markers that did not conform to the expected segregation ratio were excluded from the analysis. Three-point linkage analysis was performed for each linkage group, and the most saturated linkage was adopted as described by Ooijen [28]. After the addition of each loci, a ripple was used to verify local locus orders [28].

2.6. QTL Analysis

QTL analyses were performed on morphological data and on leaf feeding, deadheart, exit holes, and stem tunneling damage under infestation with B. fusca at the Kabete and Embu sites and the same traits under infestation with C. partellus at the Kiboko site. Joint analysis was also done using averaged means across environments and insect species. Composite interval mapping (CIM) was performed on the data using PLABQTL software, version 1.2 [29]. Whole-genome scan with CIM was conducted using an automatic cofactor selection model to determine additive effects at individual QTL, and the value of ten options were selected. Cofactors were chosen by step-wise regression and Akaike’s information criterion (AIC) of 3.0 with the “cov” statement in PLABQTL. Dependence of QTL estimation on sampling effects was estimated by a five-fold cross validation by dividing the genotypes into five subsets. The LOD threshold for declaring a putative QTL was set to 3.0 after performing 1,000 permutation tests (). All QTLs identified in this study explained ≥10% of the total phenotypic variation and were thus classified as major QTLs [30].

3. Results

3.1. Phenotypic Analysis
3.1.1. Results for Combined Analyses of Variance across Environments

The combined results for the two borer species suggested that plant damage (deadheart formation, leaf feeding damage, exit holes, and stem tunnels) and agro-morphological traits (seedling vigour and bloom waxiness) evaluated in this study were influenced significantly () by borer species tested, genotype, and the environment where the experiment was conducted (Table 1).

Results of mean squares for leaf feeding, deadheart, exit holes and morphological traits at Kiboko are presented in Table 2. There were significant differences observed for all the traits measured except for leaf feeding damage. Exit holes and stem tunneling damages were the only plant damage traits that were significantly () influenced by genotype by season interaction (Table 2). Results of mean squares for leaf feeding, deadheart, exit holes, and morphological traits at Embu are presented in Table 3. There were significant () genotypic differences on all the traits evaluated at Embu site for B. fusca (Table 3).

Results of mean squares for leaf feeding, deadheart, exit holes, and morphological traits at Kabete are presented in Table 4. Genotype significantly () influenced all the traits evaluated except leaf feeding damage (Table 4). Heritability estimates and their standard errors for the sorghum damage and morphological traits across the three environments are presented in Table 5. It was noted that heritability estimates for C. partellus at Kiboko were generally low for most of the traits compared to the two other sites where B. fusca was tested.

The contribution of B. fusca damage parameters to total grain yield through leaf feeding, deadheart, stem tunnels, and number of exit hole damages at Embu are presented in Table 6. Damage effects were partitioned into direct and indirect associations through path coefficient analysis, and grain yield was used as the resultant variable (Table 6). Deadheart, exit holes, leaf feeding, and stem tunneling damages had a negative correlation with grain yield. Deadheart had a significant direct positive effect on grain yield supported by an indirect positive effect through stem tunneling and indirect negative effects through the number of exit holes and leaf feeding damages. Exit holes had a direct positive effect on grain yield supported by indirect positive effects through deadheart and leaf feeding damages and an indirect negative effect through stem tunneling. Leaf feeding had a direct positive effect on grain yield supported by indirect positive effects through deadheart and the number of exit holes and an indirect negative effect through stem tunneling. Stem tunneling had a direct positive effect on grain yield supported by indirect positive effects through deadheart, number of exit holes, and leaf feeding damages. The contribution of B. fusca damage parameters to total grain yield through leaf feeding, deadheart formation, stem tunnels, and number of exit hole damages at Kabete is presented in Table 7.

Deadheart, exit holes, leaf feeding, and stem tunneling damages had a negative correlation with grain yield. Deadheart had a direct positive effect on grain yield supported by indirect positive effects on exit holes, leaf feeding, and stem tunneling damages. The number of exit holes had a direct positive effect on grain yield supported by indirect positive effects through deadheart, leaf feeding, and stem tunneling damages. Leaf feeding damage had a direct positive effect on grain yield supported by an indirect positive effect through deadheart damage and indirect negative effects through exit hole and stem tunneling damages. Stem tunneling had direct positive effects on grain yield supported by indirect negative effects through deadheart, the number of exit holes, and leaf feeding damages. The contribution of C. partellus damage traits to total grain yield through leaf feeding, deadheart, stem tunneling, and the number of exit hole damages at Kiboko are presented in Table 8.

Deadheart, exit holes, leaf feeding, and stem tunneling damages had a negative correlation with grain yield (Table 8). Deadheart had a direct positive effect on grain yield supported by an indirect positive effect through stem tunneling and an indirect negative effect through exit holes and leaf feeding damages. Exit holes had direct positive effects on grain yield supported by an indirect positive effect through leaf feeding and stem tunneling damages and an indirect positive effect through deadheart damage. Leaf feeding had a direct positive effect on grain yield supported by indirect negative effects through deadheart, exit holes, and stem tunneling damages. Stem tunneling had a direct positive effect on grain yield supported by indirect positive effects through deadheart, exit holes, and leaf feeding damages.

3.1.2. Correlations of Means among Stem Borer Damage and Agro-morphological Traits

The results for correlations of mean for B. fusca damage and agro-morphological traits in Embu, Kenya are presented in Table 9. A significant positive correlation was observed between panicle length and plant height (, ), implying that taller plants had a long panicle length. Stem tunneling had a significant positive association with plant height (, ) and with exit holes (, ). Seedling vigour had a positive and significant relationship with plant height (, ). Exit holes had a significant positive relationship with deadheart damage (, ). Leaf damage had a significant positive relationship with deadheart damage (, ). A significant positive correlation was observed between trichome density and leaf glossiness (, ), with seedling vigour (, ), and with bloom waxiness (, ). A significant positive correlation was recorded between seedling vigour and panicle length (, ). Bloom waxiness had a significant positive correlation with seedling vigour (, ). A significant negative relationship was recorded between seedling vigour and deadheart formation (, ) and with exit holes (, ). Bloom waxiness had a significant negative correlation with deadheart damage (, ), with exit holes (, ), and with stem tunneling (, ), respectively.

The results for correlations of mean for B. fusca damage and agro-morphological traits at Kabete, Kenya are presented in Table 10.

Leaf feeding had a significant positive correlation with deadheart damage (, ). Stem tunneling had a significant positive relationship with exit holes (, ). Bloom waxiness had a significant positive correlation with leaf glossiness (, ) and with seedling vigour (, ). Yield had a significant positive correlation with plant height (, ) and with panicle length (, ). Seedling vigour had a significant negative correlation with deadheart damage (, ), with exit holes (, ), and with stem tunneling (, ). Days to 50% flowering had a significant negative correlation with the number of exit holes (, ).

The results for correlations of mean for C. partellus damage and agromorphological traits at Kiboko, Kenya are presented in Table 11. Stem tunneling had a significant positive correlation with exit holes (, ). Days to 50% flowering had a positive correlation with leaf damage (, ) and with plant height (, ). Panicle length had a significant positive relationship with plant height (, ), with bloom waxiness (, ), and with total grain yield (, ). Seedling vigour had a significant positive correlation with bloom waxiness (, ). Days to 50% flowering had a significant negative correlation with stem tunneling (, ) and with seedling vigour (, ). Stem tunneling had a significant negative relationship with bloom waxiness (, ). Leaf feeding damage had a significant negative correlation with seedling vigour (, ) and with total grain yield (, ).

3.1.3. QTL Mapping

QTL analysis with composite interval mapping was performed for the damage traits under artificial infestation with B. fusca at the Embu and Kabete sites, and the same traits were evaluated for C. partellus at Kiboko under artificial infestation as described in Materials and Methods. The results for damage and morphological traits are shown in Table 12. Three QTLs for deadheart damage due to B. fusca were detected on chromosomes 1, 3, and 5 (flanked by markers SB261_1, DB103_3, and AD323_5) at Embu (Table 12) that explained 3.9% of the phenotypic variation, with individual QTLs accounting for 7.5% - 11.2%. A simultaneous fit with all the three QTLs based on cross validation explained 28.3% of the adjusted genetic variance. The additive gene action ranged from -7.79 to 35.13. Two of the QTLs on chromosomes 1 and 5 (flanked by markers SB261_1 and AD323_5) were from the susceptible parent, while the remaining QTLs on chromosome 3 (marker DB103_3) came from the resistant parent (PB 15520-1). The resistant parent contributed to less deadheart damage, which is a putative trait associated with resistance to stem borers in sorghum.

Five QTLs for leaf feeding damage due to B. fusca at Embu site were detected on chromosomes 2 and 10 (flanked by markers CS115_2, CS190_2, CS389_2, CS414_2, and CPS162_10, respectively) (Table 12). These explained 10.5% of the phenotypic variation, with individual QTLs accounting for 7.1% - 9.1%. A simultaneous fit with all the five QTLs based on cross validation explained 40.6% of the adjusted genetic variance. The additive gene action ranged from -6.48 to 8.13. Three of the QTLs on chromosomes 2 and 10 (flanked by markers CS190_2, CS389_2, and CPS162_10) came from the susceptible parent. The rest of the QTLs on chromosome 2 (flanked by markers CS115_2 and CS414_2) came from the resistant parent (PB 15520-1), which contributed to less leaf feeding damage, which is a trait associated with resistance to stem borers in sorghum.

One QTL for exit hole damage due to B. fusca at Embu site was detected on chromosome 7 flanked by marker GH89_7 (Table 12). This explained 6.1% of the phenotypic variation, with individual QTL accounting for 9%. A simultaneous fit with the single QTL based on cross validation explained 9% of the adjusted genetic variance. The additive gene action was -1.42. This QTL came from the resistant parent (PB 15520-1) which contributed to less exit hole damage, a trait associated with resistance to stem borers in sorghum.

Two QTLs for stem tunneling damage by B. fusca at Embu site were detected on chromosomes 7 and 9, respectively (flanked by markers GH87_7 and BF123_9) (Table 12). These explained 5.9% of the phenotypic variation, with individual QTLs accounting for 8.1% - 8.2%. A simultaneous fit with the two QTLs based on cross validation explained 16.3% of the adjusted genetic variance. The additive gene action ranged from -2.39 to -7.33. The two QTLs came from the resistant parent (PB 15520-1), an indication of the transfer of genes associated with resistance to stem tunneling by borers in sorghum.

Ten QTLs for deadheart damage due to B. fusca were detected on chromosomes 1, 2, 3, 4, 8, and 5, respectively, (flanked by markers SB691_1, CS402_2, CS259_2, CS350_2, DB172_3, DB169_3, BC149_4, JK399_8, BF97_9, and BF106_9) (Table 12). These explained 2.8% of the phenotypic variation, with individual QTLs accounting for 7.6% - 17.5%. A simultaneous fit with all the ten QTLs based on cross validation explained 80.6% of the adjusted genetic variance. The additive gene action ranged from -20.68 to 12.95. Four of the QTLs on chromosomes 1, 2, 4, and 9 (flanked by markers on markers SB691_1, CS259_2, BC149_4, and BF97_9) were from the susceptible parent, while the QTLs on chromosomes 2, 3, 8, and 9 (flanked by markers CS402_2, CS350_2, DB172_3, DB169_3, JK399_8, and BF106_9) came from the resistant parent (PB 15520) which contributed to less deadheart damage.

Seven QTLs for leaf feeding damage due to B. fusca at Kabete site were detected on markers on chromosomes 2, 3, 6, 8, and 10 (flanked by markers CS403_2, CS111_2, CS397_2, DB164_3, EF334_6, m08/014.9, and CPS158_10) (Table 12). These explained 1.3% of the phenotypic variation, with individual QTLs accounting for 7% - 13%. A simultaneous fit with all the seven QTLs based on cross validation explained 66.8% of the adjusted genetic variance. The additive gene action ranged from -14.42 to 6.56. Three of the QTLs on chromosomes 2 and 10 (flanked by markers CS111_2, CS397_2, and CPS158_10, respectively) came from the susceptible parent (ICSV 745). QTLs on markers on chromosomes 2, 3, 6, and 8 (flanked by markers CS403_2, DB164_3, EF334_6, and m08/014.9, respectively) came from the resistant parent (PB 15520-1) which contributed to less leaf feeding damage.

Eight QTLs for exit holes due to B. fusca at Kabete site were detected on chromosomes 2, 5, 6, 8, and 9 (flanked by markers CS342_2, CS190_2, AD25_5, AD185_5, EF255_6, JK208_8, BF138_9, and BF137_9) (Table 12). These explained 6.1% of the phenotypic variation, with individual QTLs accounting for 7.5%–16.5%. A simultaneous fit with all the eight QTLs based on cross validation explained 84.7% of the adjusted genetic variance. The additive gene action ranged from -7.28 to 1.17. Two of the QTLs on chromosomes 2 and 5 (flanked by markers CS342_2 and AD185_5, respectively) came from the susceptible parent (ICSV 745). The rest of the six QTLs on chromosomes 2, 5, 6, 8, and 9 (flanked by markers CS190_2, AD25_5, EF255_6, JK208_8, BF138_9, and BF137_9) came from the resistant parent (PB 15520-1) which contributed to less exit hole damage.

Nine QTLs for stem tunneling damage due to B. fusca at Kabete site were on chromosomes 1, 2, 3, 4, and 7 (flanked by markers SB466_1, CS132_2, CS133_2, CS150_2, DB140_3, DB143_3, BC140_4, GH87_7, and GH108_7) (Table 12). These QTLs explained 0.16% of the phenotypic variation, with individual QTLs accounting for 7.1% – 13.9%. A simultaneous fit with all the nine QTLs based on cross validation explained 74.29% of the adjusted genetic variance. The additive gene action ranged from -4.53 to 3.98. Four of the QTLs on chromosomes 1, 3, and 4 (flanked by markers SB466_1, DB140_3, DB143_3, and BC140_4) came from the susceptible parent (ICSV 745). The other five QTLs on chromosomes 2 and 7 (flanked by markers CS132_2, CS133_2, CS150_2, GH87_7, and GH108_7) came from the resistant parent (PB 15520-1) which contributed to less stem tunneling damage.

Four QTLs for deadheart damage due to C. partellus at Kiboko site were detected on chromosomes 2, 6, and 9 (flanked by markers CS369_2, CS389_2, EF322_6, and BF152_9) (Table 12). These QTLs explained 2.7% of the phenotypic variation with individual QTLs accounting for 8.3% – 12.6%. A simultaneous fit with all the four QTLs based on cross validation explained 41.2% of the adjusted genetic variance. The additive gene action ranged from -8.46 to 3.32. One of the QTLs on chromosome 2 (flanked by marker CS369_2) came from the susceptible parent (ICSV 745), while the other three QTLs came from the resistant parent (PB 15520-1) which contributed to less deadheart damage.

Three QTLs for leaf feeding damage due to C. partellus at Kiboko site were detected on chromosomes 2 and 6 (flanked by markers CS133_2, CS397_2, and EF184_6) (Table 12). The markers explained 4.1% of the phenotypic variation with individual QTLs accounting for 7.2% – 8.7%. A simultaneous fit with all the three QTLs based on cross validation explained 24.5% of the adjusted genetic variance. The additive gene action ranged from -3.07 to 2.76. All the three QTLs came from the resistant parent (PB 15520-1) which contributed to reduce leaf feeding damage.

Seven QTLs for exit hole damage due to C. partellus at Kiboko site were detected on chromosomes 2, 3, 4, 5, 6, and 7 (flanked by markers CS414_2, DB153_3, DB208_3, BC222_4, m05/015.6, EF255_6, and GH66_7) (Table 12). The QTLs explained 1.4% of the phenotypic variation with individual QTLs accounting for 7.2% – 10.3%. A simultaneous fit with all the seven QTLs based on cross validation explained 58.5% of the adjusted genetic variance. The additive gene action ranged from -8.04 to 0.61. Two of the QTLs on chromosomes 2 and 3 (flanked by markers CS414_2 and DB153_3) came from the susceptible parent (ICSV 745). The other five QTLs came from the resistant parent (PB 15520-1), an indication of the transfer of the resistance genes to the progenies.

Three QTLs for stem tunneling damage due to C. partellus at Kiboko site were detected on chromosomes 3 and 7 (flanked by markers DB152_3, GH70_7, and GH118_7) (Table 12). These explained 0.29% of the phenotypic variation, with individual QTLs accounting for 7.8% –11.4%. A simultaneous fit with all the three QTLs based on cross validation explained 28.7% of the adjusted genetic variance. The additive gene action ranged from -3.56 to 2.64. One of the QTLs on chromosome 3 flanked by marker DB152_3 came from the susceptible parent (ICSV 745), while the other two QTLs came from the resistant parent (PB 15520-1).

Colocalization of markers for some of the traits evaluated was detected on site analysis for B. fusca (at Embu and Kabete), and C. partellus (at Kiboko) (Table 12). For example, it was observed that the region flanked by markers (DB172_3, EF184_6, and BF123_9) on chromosomes 3, 6, and 9 conditioned trichome density, deadheart damage, and stem tunneling damage at Kabete. Marker DB169_3 on chromosome 3 conditioned leaf glossiness and deadheart damage at Kabete. Marker CS403_2 on chromosome 2 conditioned seedling vigour and leaf damage at Kabete and Embu, respectively. Marker AD323_5 on chromosome 5 conditioned resistance to deadheart formation at Embu as well as seedling vigour at Embu. Marker CS259_2 on chromosome 2 conditioned resistance to deadheart formation at Kabete as well as bloom waxiness in Kiboko. Marker GH108_7 on chromosome 7 conditioned resistance to seedling vigour at Kiboko as well as stem tunneling in Kabete. Marker AD323_5 on chromosome 5 conditioned resistance to seedling vigour and deadheart damage at Embu. Marker CS389_2 conditioned resistance to leaf damage and deadheart damage at Embu and Kiboko. Marker CS414_2 on chromosome 2 conditioned resistance to exit holes at Kiboko and leaf damage at Embu. Marker CS190_2 on chromosome 2 conditioned resistance to leaf feeding at Embu and exit holes at Kabete.

Joint analysis of B. fusca and C. partellus damage and morphological traits across the three sites is presented on Table 13. Four QTLs due to deadheart damage were detected chromosomes 4, 5, 8, and 10 (flanked by markers BC242-4, m05/023.7, JK399-8, and CPS160-10, respectively) (Table 13). These explained 2.3% of the phenotypic variation, with individual QTLs accounting for 7.5% – 9%. A simultaneous fit with all the four QTLs based on cross validation explained 24.8% of the adjusted genetic variance. The additive gene action ranged from 3.32 to 16.64. One of the QTLs on chromosome 8 (flanked by marker JK399-8) came from the resistant parent (PB 15520-1), while the other three QTLs came from the susceptible parent (ICSV 745).

Six QTLs for leaf feeding damage were detected on markers on chromosomes 1, 3, 6, 7, and 8 (flanked by markers SB255-1, DB160-3, DB169-3, EF222-6, GH70-7, and JK290-8) Table 13). These explained 2.25% of the phenotypic variation, with individual QTLs accounting for 7.6% – 17.4%. A simultaneous fit with all the six QTLs based on cross validation explained 56% of the adjusted genetic variance. The additive gene action ranged from -2.02 to 5.070. Two of the QTLs on chromosomes 3 and 7 (flanked by markers DB169-3 and GH70-7, respectively) came from the resistant parent (PB 15520-1) that contributed to reduce leaf feeding damage, while the other four QTLs came from the susceptible parent (ICSV 745).

Five QTLs for exit hole damage were detected on chromosomes 2, 4, 6, and 9 (flanked by markers CS195-2, BC149-4, EF255-6, EF416-6, and BF138-9) (Table 13). These explained 5.37% of the phenotypic variation, with individual QTLs accounting for 7% – 11.7%. A simultaneous fit with all the five QTLs based on cross validation explained 47.3% of the adjusted genetic variance. The additive gene action ranged from -0.7 to 3.2. Three of the QTLs on chromosomes 2, 4, and 6 (flanked by markers CS195-2, BC149-4, and EF255-6, respectively) came from the resistant parent (PB 15520-1), while the other two QTLs came from the susceptible parent (ICSV 745).

Six QTLs for stem tunneling damage were detected on chromosomes 1, 2, 3, 4, and 9 (flanked by markers SB466_1, CS116-2, CS132-2, DB36_3, BC140_4, and BF140_9) (Table 13). These explained 5.37% of the phenotypic variation, with individual QTLs accounting for 10.1%–14.7%. A simultaneous fit with all the five QTLs based on cross validation explained 65.5% of the adjusted genetic variance. The additive gene action ranged from -3.4 to 4.4. Two of the QTLs on chromosomes 2 and 9 (flanked by markers CS132-2 and BF140_9, respectively) came from the resistant parent (PB 15520-1), while the other four QTLs came from the susceptible parent (ICSV 745). The analyses further revealed that marker CS132-2 was colocalised for leaf toughness and stem tunneling traits on QTL 1 and 2, respectively.

4. Discussion

Heritability is used to predict response to selection and to help plant breeders know if it is more efficient to improve traits of economic importance through selection [31]. Heritability estimates thus explain the level to which genes control expression of the trait of interest. Heritability estimates for deadheart formation, leaf feeding damage, exit holes, stem tunneling, bloom waxiness, and seedling vigour tested for B. fusca at Embu and Kabete were low to high and ranged between 0.21 and 0.65. On the other hand, heritability estimates for the same aforementioned traits evaluated for C. partellus at Kiboko were low and ranged between 0.15 and 0.32. This observation suggested that genetic gain may not be realized in selection for some of the traits evaluated for the two borer species in the three sites. These traits that had low heritabilities are as follows: leaf feeding damage at Kabete and Kiboko; deadheart formation, bloom waxiness, and seedling vigour at Kiboko; and exit holes at all the three sites. Therefore, there is a need to consider other alternatives such as DNA marker technology to fix the desired genes. This observation is supported by other researchers. For example, Oloyede-Kamiyo et al. [31] worked on variability for resistance to the pink stem borer (Sesamia calamistis Hampson) and the sugarcane borer (Eldana saccharina Walker) in two tropical maize populations and observed that narrow-sense heritability was low to moderate and ranged from 1.45% for leaf feeding to 40.6% for stalk breakage.

Path coefficient analyses revealed that leaf feeding, deadheart, exit holes, and stem tunneling had a direct negative effect on grain yield. This observation suggested that either of the damage traits could lead to direct grain yield loss. Singh et al. [6] observed that direct effects of stem tunneling by C. partellus on grain yield loss were greater than leaf feeding and deadhearts in sorghum. However, Starks and Doggett [32] reported that leaf feeding damage by C. partellus in sorghum is a weak pointer of grain yield as production of new fresh leaves compensate for the leaf damage. Sandoya et al. [33] reported that improvement of maize stem tunneling resistance to the Mediterranean corn borer (Sesamia nonagrioides Lef) and the European corn borer (Ostrinia nubilalis Hbn) led to increased grain yield. The significant positive correlations observed between leaf damage and deadheart, and between stem tunneling and exit holes, implied that there is a direct relationship between the two traits, respectively, and either of the traits can be used to predict the other in the sorghum population studied. Generally, tall vigorous plants with bloom waxiness produced longer panicle length and high yield with less stem damage. There was a negative relationship between stem tunneling and bloom waxiness and between seedling vigour and deadheart formation implying that sorghum can be bred with either of the morphological traits for reduced damage with high grain yield to the two stem borers.

Colocalization of the markers for some of the traits evaluated was detected on site analysis for B. fusca (at Embu and Kabete) and C. partellus (at Kiboko). For example, it was observed that the region (flanked by markers DB172_3, EF184_6, and BF123_9) on chromosomes 3, 6, and 9, respectively, conditioned trichome density, deadheart damage, and stem tunneling damage at Kabete for B. fusca. This observation implies that the aforementioned traits can be improved by using the same linked markers. Marker DB169_3 on chromosome 3 conditioned leaf glossiness and deadheart damage at Kabete implying that both traits can be improved by using the same linked marker for B. fusca. Marker CS403_2 on chromosome 2 conditioned seedling vigour and leaf damage for B. fusca at Kabete and C. partellus at Embu, respectively, implying that both traits can be improved by using the same linked marker. Marker AD323_5 on chromosome 5 conditioned resistance to deadheart formation as well as seedling vigour at Embu (for B. fusca) suggesting that both traits can be improved by using the same linked marker. Marker CS259_2 on chromosome 2 conditioned resistance to deadheart formation at Kabete (for B. fusca) as well as bloom waxiness in Kiboko (for C. partellus) implying that the two traits can be improved by using the same linked marker. Marker GH108_7 on chromosome 7 conditioned resistance to seedling vigour at Kiboko (for C. partellus) as well as stem tunneling in Kabete (for B. fusca) implying that both traits can be improved by using the same linked marker. Marker AD323_5 on chromosome 5 conditioned resistance to seedling vigour and deadheart damage at Embu (for B. fusca) implying that both traits can be improved by using the same linked marker. Marker CS389_2 conditioned resistance to leaf damage and deadheart damage at Embu (for B. fusca) and Kiboko (for C. partellus) implying that both traits can be improved by using the same linked marker. Marker CS414_2 on chromosome 2 conditioned resistance to exit holes at Kiboko (for C. partellus) and leaf damage at Embu (for B. fusca) implying that both traits can be improved by using the same linked marker. Marker CS190_2 on chromosome 2 conditioned resistance to leaf feeding at Embu and exit holes at Kabete implying that both traits can be improved by using the same linked marker for B. fusca.

This study identified and mapped genomic regions conditioning dual resistance to B. fusca and C. partellus in sorghum. The study utilized a large number of RILs and high-density evenly spaced genetic markers scored to ensure high power and precision in QTL mapping as described by Bekele et al. [34]. A five-fold cross validation analysis was run in this study as described by Gowda et al. [30] to confirm the frequency of QTL detection that gave an estimation of the precision of QTL localization. Identification of QTLs that influence dual resistance to B. fusca and C. partellus would increase the efficiency of selection and breeding for resistance of these two cereal stem borer species. Majority of the identified loci were highly significant, and they accounted for a substantial amount of phenotypic variation for B. fusca and C. partellus. This implied that fixing these traits (leaf feeding, deadheart formation, exit holes, and stem tunneling damages) would lead to the development of resistant cultivars and thus improve food security since the losses to stem borers will be lowered. A number of reports on linkage and QTL mapping for borer resistance traits have been published in maize [3537]. QTL mapping for C. partellus resistance using the same sorghum mapping population utilized in this study is reported by Vinayan [20]. QTLs governing more than one resistance trait (pleiotropic QTLs) were identified in the present study. For example, it was observed that markers DB172_3, EF184_6, and BF123_9 on chromosomes 3, 6, and 9 conditioned trichome density, deadheart damage, and stem tunneling damage at Kabete implying that all these traits can be improved using the same linked markers. Pleotropic QTLs conditioning resistance to borers have been reported in the maize population [3840].

Chromosome 1 has been associated with multiple resistance to lepidopteran insect pests, mainly, the European corn borer (Ostrinia nubilalis (Hubner)), the southwestern corn borer (Diatrea grandiosella (Dyar)), and the sugarcane borer (Diatraea saccharalis (Fabricius)), and resistance to the maize weevil [41]. Chromosome 4 has been reported to condition resistance to stem tunneling in the tropical maize population against the corn borer [42, 43]. Chromosomes 1, 3, and 8 have been associated with resistance to Mediterranean corn borer Sesamia nonagrioides stem tunneling using maize RILs [44]. Chromosomes 1, 2, 6, and 7 in sorghum have also been reported to regulate resistance to deadheart formation against sorghum shoot fly [45]. Joint analyses across the two species of stem borers and the three sites revealed that 4 QTLs conditioned deadheart damage (chr. 4, 5, 10, and 8) and trichome density (chr. 1, 2, 4, and 6); 6 QTLs conditioned leaf feeding damage (chr. 1, 3, 6, 7, and 8), stem tunneling damages (chr. 1, 2, 3, 4, and 9), and leaf toughness (chr. 2, 4, 6, 9, and 10); 5 QTLs conditioned exit holes (chr. 2, 4, 6, and 9); and 2 QTLs conditioned bloom waxiness (chr. 3, and 9), leaf glossiness (chr. 5, and 7), and seedling vigour (chr. 8, and 9). Genomic regions governing trichome density detected on sorghum chromosomes 2, and 6 in this study have not been reported. QTLs identified for trichome density on chromosomes 1 and 3 in the joint analyses in the current study have also been reported in sorghum against shoot fly (Atherigona soccata Rond) [8, 45, 46]. Foerster et al. [47] reported that chromosome 9 conditioned bloom waxiness in maize. Aruna et al. [45] reported that chromosome 9 controlled seedling vigour in sorghum against shoot fly. Joint analyses further revealed that marker CS132-2 colocalised for leaf toughness and stem tunneling on QTLs 1 and 2, respectively, implying that both traits can be improved with the help of the same linked marker.

5. Conclusion and Recommendation

This study identified and mapped genomic regions that conditioned expression of dual resistance to B. fusca and C. partellus in sorghum. QTLs identified in this study can be used in marker-assisted selection in the development of B. fusca- and C. partellus-resistant sorghum cultivars. Marker CS132-2 colocalised for leaf toughness and stem tunneling on QTLs 1 and 2, respectively; thus, both traits can be improved for B. fusca and C. partellus with the help of the same linked marker. There is a need for further studies to identify gene(s) underlying the mapped QTLs.

Data Availability

Some of the data used to support the findings of this study are included in the article. Additional data are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

Financial support was obtained from Regional Universities Forum for Capacity Building in Agriculture (RUFORUM), Carnegie Corporation of New York, USA and the Harnessing Opportunities for Productivity Enhancement (HOPE) of Sorghum and Millets in Sub-Saharan Africa and South Asia Project funded by the Bill and Melinda Gates Foundation implemented by ICRISAT.

Supplementary Materials

Table 14 Coded markers and their names. (Supplementary Materials)