Next Article in Journal
Genome-Wide Profiling of Enterotoxigenic Staphylococcus aureus Strains Used for the Production of Naturally Contaminated Cheeses
Next Article in Special Issue
Characteristics of Microsatellites Mined from Transcriptome Data and the Development of Novel Markers in Paeonia lactiflora
Previous Article in Journal
Transcriptome-Wide Analysis Reveals Key DEGs in Flower Color Regulation of Hosta plantaginea (Lam.) Aschers
Previous Article in Special Issue
Transferability and Polymorphism of SSR Markers Located in Flavonoid Pathway Genes in Fragaria and Rubus Species
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations

1
Directorate of Research, Ministry of Agriculture and Food Security, Ministries Complex, Parliament Road, P.O. Box 33, Juba, South Sudan
2
International Maize and Wheat Improvement Center (CIMMYT), World Agroforestry Centre (ICRAF), United Nations Avenue, Gigiri. P.O. Box 1041–00621, Nairobi, Kenya
3
West Africa Centre for Crop Improvement (WACCI), College of Basic and Applied Sciences, University of Ghana, PMB 30, Legon, Ghana
4
Department of Agricultural Sciences, College of Natural Resources and Environmental Studies, University of Juba, P.O. Box 82, Juba, South Sudan
*
Author to whom correspondence should be addressed.
Genes 2020, 11(1), 32; https://doi.org/10.3390/genes11010032
Submission received: 24 October 2019 / Revised: 17 December 2019 / Accepted: 24 December 2019 / Published: 26 December 2019
(This article belongs to the Special Issue Genetic Diversity Assessment and Marker-Assisted Selection in Crops)

Abstract

:
Maize lethal necrosis (MLN) occurs when maize chlorotic mottle virus (MCMV) and sugarcane mosaic virus (SCMV) co-infect maize plant. Yield loss of up to 100% can be experienced under severe infections. Identification and validation of genomic regions and their flanking markers can facilitate marker assisted breeding for resistance to MLN. To understand the status of previously identified quantitative trait loci (QTL)in diverse genetic background, F3 progenies derived from seven bi-parental populations were genotyped using 500 selected kompetitive allele specific PCR (KASP) SNPs. The F3 progenies were evaluated under artificial MLN inoculation for three seasons. Phenotypic analyses revealed significant variability (P ≤ 0.01) among genotypes for responses to MLN infections, with high heritability estimates (0.62 to 0.82) for MLN disease severity and AUDPC values. Linkage mapping and joint linkage association mapping revealed at least seven major QTL (qMLN3_130 and qMLN3_142, qMLN5_190 and qMLN5_202, qMLN6_85 and qMLN6_157 qMLN8_10 and qMLN9_142) spread across the 7-biparetal populations, for resistance to MLN infections and were consistent with those reported previously. The seven QTL appeared to be stable across genetic backgrounds and across environments. Therefore, these QTL could be useful for marker assisted breeding for resistance to MLN.

1. Introduction

Maize lethal necrosis (MLN) is a major disease in sub-Saharan Africa (SSA) caused by co-infections of maize chlorotic mottle virus (MCMV) and sugarcane mosaic virus (SCMV) [1]. MCMV can able to interact with any member of the Potyviridae family to cause lethal necrosis in maize [2]. Yield loss due to MLN can reach up to 100% under severe infection and MLN favorable environments [1,3]. Breeding for host resistance to MLN is the most effective means of preventing yield losses in farmer’s fields. Application of molecular markers could enhance breeding for resistance to MLN. Although markers are widely used in breeding for crop improvement including maize, the tools are inconsequential unless the linked markers or quantitative trait loci (QTL) are tested for their effectiveness and reproducibility in different genetic backgrounds. QTL validation adds weight to assess the effectiveness of alleles and their linked markers.
QTL mapping approaches are one of the popular genomic tools to dissect the genetic architecture of complex traits [4]. The presence of QTL conferring resistance to several viral diseases in maize has been investigated in numerous linkage and linkage disequilibrium mapping studies [4,5,6,7,8]. QTL mapping or linkage mapping is known for high QTL detection power. Joint linkage association mapping (JLAM) based on different segregating biparental populations is known to provide both the high QTL detection power and high mapping resolution [9,10]. Application of both linkage mapping and JLAM in multiple bi-parental populations is useful to validate earlier findings and to detect possible new sources of resistance for MLN.
Several studies were conducted to validate QTL effects on traits of economic importance in different crops including maize [11,12,13,14,15]. Sukruth et al [16] validated the markers associated with late leaf spot and rust in groundnut using recombinant inbred line (RIL) population and two backcross populations. Zhou et al [17] validated two major QTLs (LEN-3H and LEN-4H) for kernel length in wild barley using a biparental population derived from Fleet (Hordeum vulgare L.) and Awcs276. The authors reported significant association between the two QTL and kernel length. Gowda et al [7] employed four biparental maize populations and adopted linkage mapping and joint linkage mapping options to discover and confirm QTL associated with resistance to MLN. They consistently confirmed in two of the populations that three major QTLs were localized on chromosomes 3, 6, and 9. For effective use of trait linked markers and resources, validation of discovered QTL with large expected impact is crucial, as there are too many QTLs to validate and validation population development and assessment is expensive [18,19,20]. It is a pre-requisite though that for QTL to be effectively used in crop improvement, it should be confirmed that their effects remain consistent across populations and environments [12,16,21].
QTL validation study greatly contributes towards increased resolutions of some of the target QTLs as well as complementarity to previous findings either from genome wide association studies (GWAS) or JLAM or any other approaches. The CIMMYT Global Maize Program has recently identified a number of QTL across maize chromosomes which are associated with resistance to MLN in multiple mapping populations [7,8]. The validation of these QTL using seven different F3 mapping populations would provide better understanding of QTLs associated with resistance to MLN and justify their application for marker assisted breeding towards improvement of maize lines for resistance to MLN. Further this study also helps to find new source of resistance which might not have been reported in earlier studies. This study was aimed to: (i) evaluate seven different F3 populations for their responses to MLN under artificial inoculation; (ii) conduct individual population-based QTL mapping; (iii) apply JLAM with three biometric models and compare with linkage mapping results to identify stable and/or unique QTL; and (iv) assess the potential of genomic prediction for MLN resistance within biparental populations with low marker density.

2. Materials and Methods

2.1. Plant Materials

Seven elite maize parental lines with contrasting response to MLN developed by CIMMYT through pedigree and DH breeding schemes were used in this study. Parents CKDHL120918, CML494, CKLTI0227, and CKDHL312 were tolerant to MLN; CML543 and CKDHL221 were moderately tolerant, and CKDHL0089 was susceptible to MLN. These materials were also known for tolerance to various biotic and abiotic stresses with good agronomic performances. The seven bi-parental populations (F3 pop1–CKDHL120918 × ML494, F3 pop2–CML543 × CML494, F3 pop3–CKDHL120918 × CML543, F3 pop4–CKLTI0227 × CKDHL120918, F3 pop5–CKDHL0089 × CKDHL120918, F3 pop6–CKDHL0221 × CKDHL120312 and F3 pop7–CKDHL0089 × CML494) were used for linkage mapping and JLAM. To develop F3 populations, crossing blocks were established in the nursery in Kiboko, Kenya (37°75′ E; 2°15′ S; 975 m a.s.l.; of 530 mm/year of rain fall and temperature ranges from 14.3 to 35.1 °C) during the 2016/2017 cropping season. Seven bi-parental crosses were made among the seven elite parental lines. Single cross (F1) seeds were grown and F1 plants were selfed. About 300 to 350 F2 plants from each population were randomly selfed and F3 seeds were harvested.

2.2. Phenotypic Evaluation

At least 306 F3 families from each population with their seven parental lines and six commercial checks were evaluated to determine their response to MLN under artificial inoculation in field. Experiments were conducted for three seasons (April 2017, April 2018, and October 2018) in confined MLN facility in Naivasha (36°26′ E; 0°43′ S; 1896 m a.s.l.; 677 mm/year of rain fall and temperature ranges from 12 to 29 °C), hereafter seasons are referred as environments. Trials were planted using alpha lattice design in a 1-row plot of 3.0 m long, with spacing of 0.75 m between rows and 0.25 m between plants. Two seeds were planted per hill and thinned to one plant per hill 3 weeks after germination, making a total of 13 plants per row. Standard agronomic practices were adopted.
MLN inoculum was obtained from separate, sealed greenhouses maintained in Naivasha for each of SCMV and MCMV [7]. Maintenance of virus stocks in the greenhouses was described earlier [22]. In brief, MLN infected leaf tissues were collected from the field and ground in grinding buffer solution at 1:10 dilution ratio (10 mM potassium-phosphate, pH 7.0) [1,22]. The resulting sap extract was centrifuged at 12,000 rpm for 2 min. The sap was decanted and celite powder was added at the rate of 0.02 g/mL. To propagate the viruses in the greenhouses, a susceptible maize hybrid (H614) was grown and infected by rubbing the sap on the leaves at two to four leaf stages. A separate, sealed greenhouse was the maintained for each virus as stock for further use.
Three weeks before inoculation of the experimental materials, enzyme linked immunosorbent assay (ELISA) test was conducted on leaf samples, randomly collected from infected plants in the greenhouses, to determine the presence and purity of the MCMV and SCMV [7]. Separate extracts from the SCMV and MCMV infected plants was prepared and the two extracts were then mixed to form MLN inoculum at the ratio of 4 parts of SCMV to 1 part MCMV (weight/weight). Two inoculations were applied at the 4th and 5th week after planting. In order to keep uniform disease pressure across experiments, a motorized, backpack mist blower (Solo 423 Mist Blower, 12 L capacity) with an open nozzle (2-inches diameter) was used to inoculate the experimental plants at a high inoculum delivery pressure of 10 kg/cm. Drip irrigation was used to provide water and fertilizer. All other agronomic practices relating to maize production were followed according to standard procedures for field practices. Spreader rows of susceptible maize hybrid (H614) were planted along the experiment to enhance disease spread and intensity [16,23].
MLN Disease severity (MLN-DS) score started 10 days after the 2nd inoculations and was recorded four times at 10 days interval using standardized qualitative scale of 1 to 9, where 1 = resistant, clean, no symptoms; 2 = fine or no chlorotic specks, but vigorous plants; 3 = mild chlorotic streaks on emerging leaves; 4 = moderate chlorotic streaks on emerging new leaves; 5 = chlorotic streaks and mottling throughout plants; 6 = intense chlorotic mottling throughout plants, necrosis on leaf margins; 7 = excessive chlorotic mottling, mosaic and leaf necrosis, at times dead heart symptoms; 8 = excessive chlorotic mottling, leaf necrosis, dead heart and premature death of plants; and 9 = susceptible (complete plant necrosis and dead plants) [7,8]. After analyzing MLN-DS for each time score, we chose the third score (40 days post-inoculation) for further analysis because of its higher heritability and full expression of disease symptoms. The area under the disease progress curve (AUDPC) was calculated for each plot using SAS 9.4 (SAS Institute Inc., 2015, Cary, NC, USA) so as to understand the trend of development of MLN severity across the score intervals.

2.3. Phenotypic Data Analysis

All quantitative genetic parameters were estimated based on the performance of the 2142 F3 lines. To check the quality of the data, first, test for normality of distribution of error terms for MLN-DS and AUDPC was determined using the Shapiro-Wilk test [24]. Secondly, analyses of description statistics (mean, range, skewness and kurtosis) and correlation among phenotypic traits were performed using META-R [25]. We assumed that each environment was a representative of a replication. Therefore, statistical analysis of phenotypic data was conducted using linear mixed model:
Yijk = μ + gi + lj + bkj + εijk,
where Yijk was the disease severity of the ith genotype at the jth environment in the kth incomplete block, μ was an intercept term, gi was the genetic effect of the ith genotype, lj was the effect of the jth environment, bkj was the effect of the kth incomplete block at the jth environment, and εijk was the error term confounding with the genotype-by-environment interaction effect. To determine variance components by the restricted maximum likelihood (REML) method, both block and genotype effects were treated as fixed. Significance of variance component estimates was tested by model comparison with likelihood ratio tests where the halved P-values were used as approximations [26]. Heritability (H2) on an entry-mean basis was estimated as the ratio of genotypic to phenotypic variance. The phenotypic variance comprises genotypic variance and the masking GxE interaction variances divided by the number of environments. Further, the mixed linear model (MLN) established in META-R software (http://hdl.handle.net/11529/10201) was adopted and the best linear unbiased predictor (BLUP) and best linear unbiased estimator (BLUE) for each genotype across environments were generated. The BLUPs were used in linkage mapping and joint linkage association mapping analyses whereas BLUEs were used in genomic prediction studies.

2.4. DNA Extraction, Genotyping, Linkage Map Construction and QTL Analysis

In this study, SNP markers were first screened on parental lines and 500 markers, which showed polymorphism between at least two of the seven parents were selected and used to genotype all the populations. Leaf samples were collected from 306 individuals per population 2 to 3 weeks after emergence, based on CIMMYT laboratory protocols [27]. The leaf samples were sent to LGC Genomics (https://www.biosearchtech.com/services/genotyping-services, Herts, UK) and were genotyped using 500 SNPs. Genotypic data obtained from LGC was subjected to quality check and SNPs were called and filtered using TASSEL version 5.0 software [28].
Polymorphic markers for each population were selected. Segregation of each SNP was verified for deviation from classic Mendelian inheritance using x2-test and SNPs that significantly deviated were discarded [29,30]. Linkage maps were created for the individual populations using the MAP function established in QTL IciMapping v 4.1 software [24,31]. Linkage groups were identified using Group command based on logarithm of odds (LOD) score of 3.0, and recombination rate were converted into centimorgans (cM) using Kosambi mapping function [32]. Ordering of the markers was conducted using the “ordering” instruction with the nnTwoOpt algorithm. Adjustment of the map order was done according to the sum of adjacent recombination frequencies (SARF) and sorted in the “rippling” instruction with a window size of 5 as the amplitude. Instruction generated from the map was used to draw and visualized the map using Map Chart software [33].
The BLUPs obtained from across environments for MLN-DS and AUDPC for each population were subjected to inclusive composite interval mapping (ICIM) analysis so as to determine the QTL linkage. The ICIM is an effective two-step statistical approach that allows separation of co-factor selection from interval mapping process, in order to control the background effects and improve mapping of QTL with additive effects [31]. A LOD threshold of 3.0 with a scanning step of 1 cM were used to declare significant QTL [21,34]. Stepwise regression was adopted to determine the percentages of phenotypic variance explained (R2) by individual QTL and additive effects at LOD peaks. QTL nomenclature [35] was adopted to nominate QTL conferring MLN resistance, where a two or three letter abbreviation of trait, followed by the chromosome number on which the QTL is found and the marker position to distinguish multiple QTLs were employed. The percentages of phenotypic variance explained (% PVE) by individual QTL, and additive and dominant effects at LOD peaks were generated. Sources of favorable alleles were determined depending on signs of the QTL additive effects [7,36]. For each F3 population, estimated additive (a) and dominance (d) effects for each QTL were used to calculate the ratio of dominance level (|d/a|). This ratio was used to classify the nature of QTL as follow [37]: additive (A; 0 ≤ |d/a| ≤ 0.2); partially dominant (PD; 0.2 ≤ |d/a| ≤ 0.8); dominant (D; 0.8 < |d/a| ≤ 1.2) and over dominant (OD; |d/a| > 1.2).
Based on the SNP markers shared by different populations, an integrated map was built by using IciMapping software [31]. In brief, SNPs overlapped across genetic maps were selected as anchor markers and used to integrate corresponding linkage groups on individual linkage maps. The marker order and marker positions were calculated after calculating the order and the relative position (within each genetic map) of the anchored markers, followed by integrating of all the detected markers into one map. Then all QTL identified from the seven populations were projected onto the integrated map based on their confidence interval.

2.5. Joint Linkage Association Mapping (JLAM)

The JLAM method is a combination of linkage mapping and association mapping approaches [38,39]. Here, BLUPs obtained from across seven bi-parental populations and about 420 SNPs with missing values of <5% and minor allele frequency >0.05 were considered for JLAM study. We employed three biometric models to elucidate the QTL-trait relationships [7]. For each model, first, co-factors were selected using stepwise multiple linear regression based on the Schwarz Bayesian Criterion [40], following the Proc GLM SELECT model of SAS 9.2 [41]. Secondly, P-value and F-test were determined based on the full model (with QTL effects) and the reduced model (without QTL effects), followed by QTL scan using R software (version 3.5.3). The approach for Model A involved the trait as a function of the co-factor + marker. Model B incorporated population effect as an additional factor to adjust for the population structure so that trait = population + co-factor + marker; and Model C involved nesting of the co-factor and marker within population such that trait = pop + co-factor (pop) + marker (pop) [10,42]. Significance of association between QTL and MLN resistance was determined at P < 0.05 following Bonferroni-Holm procedure [43]. The adjusted total phenotypic variance explained (R2) values for the detected QTL were generated when the significant QTL were simultaneously fitted in linear model. Further principal component analysis (PCA) of all F3 lines was carried out using TASSEL version 5.2 [28], and the CurlyWhirly version 1.15 software [44] was used to construct PCA biplot for the population structure.

2.6. Genomic Prediction

Due to its advantages over either phenotypic and marker assisted selection, genomic prediction (GP) is becoming a popular approach for breeding, especially for complex traits [45,46]. The technique allows for the selection of superior genotypes without conducting phenotypic evaluation, if a subset of the population has genotypic and phenotypic data. Here, we performed GP analyses using ridge-regression BLUP (RR-BLUP) with five-fold cross-validation. Uniformly distributed, polymorphic, SNPs from each population showing minor allele frequency of less than 5% were used.
Genotypes were sampled from each population and used to constitute both training set and testing set. The sampling of the sets was repeated 100 times [7]. The predictive ability of the GP was computed by dividing the correlation between genomic estimated breeding values (GEBVs) and the observed phenotypic values by the square root of the heritability estimates obtained from the respective populations [46].

3. Results

3.1. Response of Parents and F3 Populations to MLN Infections

The response of F3 lines for MLN-DS and AUDPC showed continuous distribution, ranging from highly resistant or tolerant to completely susceptible in each individual population as well as across populations (Figure 1).
The normality tests indicated that the means were positively and moderately skewed with values ranging from 0.65 to 0.77. Kurtosis values were platykurtotic with maximum of 2.07 for AUDPC. High W-test values were observed for all the traits and ranged from 0.95 to 0.97 (data not shown). Individual means across populations ranged from 4.40 to 5.11 with a mean of 4.70 for MLN-DS, and the AUDPC mean values for individual populations ranged from 125.6 to 134.4 with an across population mean of 132.4. The magnitude of genotypic variance for MLN-DS was lowest (0.17) in F3 pop 2 (CML543 × CML494) and highest (0.69) in F3 pop 4 (CKLTI0227 × CKDHL120918). Genotypic variances were highly significant (P < 0.01) in all seven F3 populations and across populations for both MLN-DS and AUDPC values. Moderate to high broad–sense heritability estimates of 0.62 to 0.79 were found for MLN-DS and AUDPC values, indicative of high-quality phenotypic data for further genetic analysis (Table 1).
Three-dimensional principal component analysis (PCA) of the seven bi-prenatal populations are shown in Figure 2. The results identified five distinct groups including CML543 × CML494, CKDHL120918 × CML543, CKDHL0089 × CKDHL120918, CKDHL0221 × CKDHL120312, and CKDHL0089 × CML494. The distinctiveness of the five populations implies their diversity in their genetic backgrounds. However, populations CKDHL120918 × CML494 and CKLTI0227 × CKDHL120918 were overlapping and CKDHL221 × CKDHO120312 is distinct from other populations. This might indicate their relatedness and distinctness in genetic compositions.

3.2. Molecular Analyses

3.2.1. Linkage Group

Genetic linkage groups consisting of 10 maize chromosomes were constructed for each of the seven F3 populations. Marker density on each map varied among the seven populations. F3 pop 4 carried the largest number of polymorphic markers with 298 SNPs and a total length of 1550.07 cM at an average density of 5.20 cM between the markers. F3 pop 1 had only 112 polymorphic SNPs spanning a total length of 1223.97 cM with a mean spacing of 10.93 cM between markers. Similarly, variation was detected in number of SNPs per linkage group (LG), with marker density per LG ranging from 40 with a mean spacing of 5.71 for LG7 to 153 and an average of 21.86 on LG1. Total linkage map consisted of 389 SNPs with a total length of 2007.85 cM and an average marker interval length of 5.16 cM (Table S1).

3.2.2. QTL Mapping

QTL analyses across environments revealed 60 and 58 significant QTLs for MLN-DS and AUDPC values for the seven populations, respectively (Table 2), with different sizes of QTL effects. Generally, the study revealed that positive alleles for resistance to MLN were donated either by male or female parent in a cross in each population. Some QTL identified in F3 pop 6 and 7 showed large additive effects. QTL detected in the other F3 populations predominantly showed dominant to over-dominant effects. Majority of the QTLs was detected in F3 pop4 (CKLTI0227 × CKDHL120918), whereby 7% of the QTL showed additive effects for MLN-DS and AUDPC values.
In F3 pop 1, for MLN-DS and AUDPC values, together 15 QTL were detected, which explained 30.08% and 34.72% of the total PVE, respectively. Most of the detected QTL were of small effects except qMLN6_85, which had AUDPC values explaining 21.6% of PVE. In F3 pop 2, seven and 10 QTLs were detected which together explained 43.5 and 49.8% of the total phenotypic variance for MLN-DS and AUDPC values, respectively. Among all the QTLs identified, only three QTLs are specific to AUDPC values. In F3 pop3, a total of eight QTLs were found for MLN-DS, explaining 47.9% of the total observed phenotypic variability. Among the QTL with major effects, qMLN3_142 was detected with 28.8% and 11.1% of PVE for MLN-DS and AUDPC values, respectively. There were 19 and 10 QTLs detected which together explained 25.1 and 28.7% of the PVE for MLN-DS and AUDPC values, respectively in F3 pop4. All the detected QTLs were of small effects. For F3 pop5, seven QTL each were found for MLN-DS and AUDPC values which together explained 54.3% and 57.6% of the total phenotypic variance, respectively. QTL qMLN3_142 was found for both MLN-DS and AUDPC values with major effects. F3 pop6 revealed five QTLs explaining 52.1% of the variation for MLN-DS. Whereas for AUDPC values, nine QTLs were found which explained 59.1% of the total phenotypic variance. QTL qMLN3_130 consistently explained >25% of the phenotypic variance. F3 pop7 revealed a total of three and four QTL with 46.7% and 50.8% PVE for MLN-DS and AUDPC values, respectively.
Among several genomic regions identified with stable QTL for MLN resistance, QTL on chromosome 3 was highly consistent. We found two strongly associated SNPs (PHM15449_10 at 125,077,922 bp and PZA00920_1 at 142,821,031 bp) within this region and their allelic effects on MLN resistance in different populations were also prominent. The phenotypic values of the different allele classes of these two major-effect SNPs in five F3 populations and across populations for MLN-DS and AUDPC value were presented in Figure 3.

3.2.3. Consensus Map Construction

Since all populations were genotyped with the same set of SNPs, we integrated the seven maps into a consensus linkage map based on the markers shared by populations (Figure 4). From the total of 118 QTLs (60 for MLN-DS and 58 for AUDPC), similar QTL detected for both traits were treated as single QTL, as a result 77 QTLs were mapped on the consensus map. Among all the QTLs detected from seven populations, 19 QTLs were consistently found in at least two populations; QTL qMLN3_142 on chromosome 3 was consistently identified in six of the seven populations. Two QTLs (qMLN1_265 and qMLN6_157) were consistently found in four populations, seven QTLs (qMLN1_47, qMLN3_130, qMLN4_150, qMLN6_85, qMLN7_158, qMLN8_10, and qMLN10_9) were repeatedly identified in three populations and nine QTLs (qMLN3_146, qMLN4_30, qMLN5_42, qMLN5_160, qMLN5_190, qMLN5_202, qMLN9_109, qMLN9_142, and qMLN10_114) were consistently found in two populations. These clustered QTLs from different populations were placed on the consensus map (Figure 4).

3.2.4. Joint Linkage Association Mapping (JLAM)

The JLAM analyses based on three biometric models together found 27 and 28 main effect QTL for MLN-DS and AUDPC values (Table 3). For MLN-DS, 15, 12 and 18 QTLs were detected, which together explained 34.4%, 27.3% and 29.1% of the total variation in model A, B and C, respectively. Among 27 QTLs, five QTLs were consistently detected in all three models for MLN-DS. For AUDPC values, for model A, B and C, we found 14, 12 and 22 QTL which together explained 33.6%, 29% and 39.8% of the total variance. Seven QTLs were consistently detected in all three models. Most of the QTL detected through JLAM showed minor effects except one QTL qMLN3_148 which was also consistent across models as well as across traits. Genomic prediction (GP) was attempted to predict the genetic values in the seven F3 populations for MLN-DS and AUDPC performances. The mean prediction accuracy ranged from as low as 0.12 in F3 pop1 to 0.75 in F3 pop 7 for MLN-DS (Figure 5). The prediction accuracies were similar for AUDPC values with range of 0.13 in F3 pop1 to 0.75 in F3 pop5.

4. Discussion

4.1. Response of Parents and F3 Populations to MLN Infections

A total of 2142 F3 genotypes and seven parental lines were evaluated for MLN response in an artificially inoculated MLN plots for three seasons (2017–2018). Average mean for the F3 populations was lower than parental average for MLN-DS. Superiority of the populations over parents could be due to dominance and over-dominance effects resulting from combinations of favorable alleles from different parents and therefore, had better heterosis for resistance to MLN. The concept of heterosis has been widely reported in maize [47,48,49]. Significant mid-parent and best-parent heterosis for resistance to MLN across locations were reported in Northern Tanzania [50]. The observed heterosis could also be partially explained by higher vigor in the progenies than the inbred parents. The inbred parents are not as vigorous as the heterozygous progenies and less able to cope with stresses of various kinds including pathogens.
The significant variability observed among the populations for resistance to MLN is possibly due to genetic effects and could be an indication of differences in genetic backgrounds of the genotypes for resistance to MLN (Table 1). Further, the study detected high heritability estimates (0.71 to 0.94) for resistance to MLN in all the seven populations. Previous investigations revealed moderate to high heritability estimates (0.34–0.89) for resistance to MLN [7,22,48,51]. The high heritability estimates detected in this study corroborate earlier reports and strongly suggest that resistance to MLN could be largely influenced by several genes with major effects.

4.2. QTL Analyses

In this study, seven major QTL were identified for resistance to MLN, localized on chromosomes 3 (qMLN3_130 and qMLN3_142), 5 (qMLN5_190 and qMLN5_202), 6 (qMLN6_85 and qMLN6_157), 8 (qMLN8_10) and 9 (qMLN9_142) and spread across the 7-biparetal populations with the percentage of PVE ranging from 10.54% to 44.50% (Table 2). The results indicate that a large portion of the phenotypic variation in MLN resistance can be explained by few QTL with major effects while the remaining portion is due to QTL with minor effects. Similar observations have been made on various quantitative traits [52], including MLN [7,8]. It was observed in the present study that some QTLs showed major effects in some populations and minor effects in other populations. This variability could be attributed to various factors and gene actions including QTL × QTL interactions or QTL × environment interactions which might impact the size of effect of any given QTL in any given biparental population.
To better understand the above interactions and their cause of variability among the genotypes, it is important to consider each component of gene action (additive, dominance and epistasis) separately. Additive effect is realized when each allele is independently contributing to genetic variability. Dominance effect is a variation caused by interaction between alleles at a locus. Epistasis effect is observed when the PVE is due to interactions between alleles at different loci. In this study, it was observed that QTL had varying levels of additive and dominance effects for resistance to MLN (Table 2), suggesting the importance of both modes of gene action in conditioning resistance to the disease. Similarly, some QTL had smaller %PVE values but had larger additive or dominance effects compared to the QTL with larger %PVE. Generally, when epistasis is involved, QTL with moderate and major additive effects are much more affected and they tend to have reduced %PVE. Therefore, epistasis could also be one of the main causes in some QTL having larger additive effects but with reduced %PVE values for resistance to MLN. QTL detected in F3 pop 1, 2, 3, 4, and 5 predominantly showed dominant to over-dominant effects. The dominance effects observed indicates MLN resistance is possibly due to interactions between individual alleles at specific loci. Other factors may include recombination between QTL and markers, leading to change in number of expected resistant alleles in the progeny hence reduced penetrance in QTL effects [53,54]. For QTL to be useful in a breeding program, it is important to first carry out validation studies to confirm their repeatability in different genetic backgrounds and environments. In line with expectations, the major QTL observed across populations in this study such as qMLN3_130, qMLN3_142, qMLN6_85, qMLN5_190, etc., indicated their stability across different genetic backgrounds.
It was observed that favorable alleles for each QTL were contributed by either resistant or susceptible parent, which is indicated here by sign of the additive effect of the QTL. Negative (−) additive effect was considered for resistant parent being the source of favorable allele for resistance to MLN. Similarly, positive (+) additive effect implied that the susceptible parent was the donor of favorable allele for the observed phenotypic variability. This suggests the importance of major QTL in reduction of MLN infections and their stability across different genetic backgrounds [54]. Earlier studies [7,8] also detected contributions of favorable alleles from both resistant and/or susceptible parents.
At least seven major QTLs localized across chromosomes and across populations with significantly larger %PVE were detected with total %PVE ranging from 25.13 to 59.07. Three of these QTLs (qMLN3_130, qMLN6_85 and qMLN5_190) showed effects (%PVE), chromosome positions and confidence intervals (CI) comparable to the previous study [7]. For example, for qMLN3_130, the current CI is 125 to 169 Mb, whereas in the previous study CI was 125–130 Mb; for qMLN6_85 CI is 84–91 Mb and from previous study CI was 85–96 Mb; and for qMLN5_190 CI is 23–191 Mb and previous CI was 190–191 Mb. The larger CI observed in this study, compared to earlier study, could be due to few SNPs used in this study, however the results confirms the stability of the QTL in specific genomic regions. Another major QTL (qMLN6_157) detected on chromosome 6 in the current study was also found in the previous investigations [8]. This observation implies that three QTLs are stable across diverge genetic backgrounds, and hence can be useful in marker assisted breeding for resistance to MLN. Another major effect QTL qMLN3_142 was detected in six out of the seven populations in this study, which also overlapped with QTL for MLN resistance in DH populations from the earlier study [8]. Two major QTLs (qMLN5_202 and qMLN8_10) mapped across populations in this study, have not been reported previously. This implies the novelty of these genomic regions for MLN resistance. Similarly, it could suggest their specificity to the genetic backgrounds used in the current study.

4.3. Joint Linkage Association Mapping (JLAM)

For JLAM, three biometric models were used to detect the maximum number of QTL linked to genes for MLN resistance. The principle of both models A and B resemble the composite interval mapping approach used in bi-parental populations [55] where co-factors were used to adjust for the population structure and background error within the segregating populations, leading to enhanced capacity to detect QTL. In both MLN-DS and AUDPC, we observed that >65% of QTL detected in model A and B overlapped with slightly higher %PVE in Model A (Table 3). This clearly indicates that control of population stratification by using co-factors alone was moderate, so population effect is also important. The slightly higher number of QTL detected by model A compared to model B could imply that model A was able to more effectively detect the variability among populations, and hence was able to reveal more QTL [10]. Model C with the nested marker effects, exploited the linkage disequilibrium (LD) within segregating populations and assumed population-specific QTL effects [10]. Model C was comparable with Model B in terms of the %PVE by these QTLs. This was well supported by all seven populations with relatively bigger population size (n = 306), where the nested model C is able to detect QTL with small to large effect size.
To evaluate the reliability of QTL detected in linkage mapping, we compared the identified 60 and 58 QTLs with 27 and 28 main effect QTLs detected through JLAM for MLN-DS and AUDPC values, respectively. Twenty-three QTLs are common across linkage mapping and JLAM for both the traits (Table 2 and Table 3). Comparison of QTL detected through JLAM and linkage mapping revealed five unique QTLs (qMLN3_199, qMLN4_235, qMLN5_207, qMLN5_209, and qMLN6_13) which were detected only through JLAM. QTLs common across methods helped to pinpoint the specific markers compared to QTLs detected with large CI through linkage mapping. Plausible gene candidates underlie some of the detected QTL regions with significant SNPs found in genes involved in plant defense system, like PZA00447_8 associated with serine/threonine protein phosphatase, PZA01313_2 linked to leucine-rich repeat protein and PZA00726_8 linked to zinc ion binding function which has a role in plant defense responses (Table S2). Among the unique QTL detected in JLAM, QTL on chromosome 6 seems to be overlapping with earlier reported major effect QTL for MLN resistance [7]. Major recessive QTL was reported in this region [56] in F2 populations. Similarly, with different association panels new QTL detected for resistance to SCMV in the same location on chromosome 6 [57]. Unique QTL on chromosomes 3 and 6 found in both model A and B but not in linkage mapping in this study indicates the variation across populations has helped to detect these QTL, which supports the use of multiple populations to find the novel source of resistance. On the contrary, the other three unique QTLs on chromosomes 4 and 5 were detected only in model C and were QTLs with population specific expression.
We employed the same populations to evaluate GP by sampling both training and testing sets. Moderate to high accuracy of GP was observed for MLN-DS and AUDPC within populations. Low accuracy in F3 pop1 is also due to low number of polymorphic SNPs compared to other populations. Application of GP for selection for improvement of MLN resistance has been reported [7]. The previous investigations detected high prediction accuracy within the populations for both MLN-DS and AUDPC values. Resistance for MLN is complex and requires considerable resources and multiple selection cycles to identify resistant lines using phenotypic selection alone. Therefore, GP could provide an alternative practical approach for breeding for resistance to MLN because it accounts both major and minor effect QTL for accuracy, enhanced breeding gain (shortened breeding cycle) and reduced costs.

5. Conclusions

Seven major QTL were identified for resistance to MLN, localized on chromosomes 3 (qMLN3_130 and qMLN3_142), 5 (qMLN5_190 and qMLN5_202), 6 (qMLN6_85 and qMLN6_157), 8 (qMLN8_10) and 9 (qMLN9_142) and spread across the 7-biparetal populations with a percentage of PVE ranging from 10.54% to 44.50%. Two SNPs on chromosome 3 (PHM15449_10 at 125,077,922 bp and PZA00920_1 at 142,821,031) were strongly associated with MLN resistance and their allelic effects in different populations were also prominent. Several QTL on chromosomes 1, 3, 5, and 6 shared similar CI with those reported previously. These QTL could be adopted for marker assisted breeding for improvement of maize against MLN infection. Additional three major QTL were not reported before further research is warranted to confirm their reproducibility and use in breeding for resistance to MLN. Both linkage mapping and JLAM can be incorporated to confirm earlier reported QTLs and discover new sources of resistance. Incorporation of GP can help to capture both large effect and small effect QTL to improve the level of MLN resistance as a result improved genetic gain. Twenty-three QTL were common across linkage mapping and JLAM however, five unique QTL (qMLN3_199, qMLN4_235, qMLN5_207, qMLN5_209, and qMLN6_13) were unique to JLAM with the QTL on chromosome 6 appearing to be overlapping with an earlier reported major effect QTL for SCMV resistance. Both linkage mapping and JLAM can be incorporated to confirm earlier reported QTL and discover new source of resistance. The incorporation of genomic selection can help to capture both large effect and small effect QTL to improve the level of MLN resistance.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/1/32/s1, Table S1: genetic linkage groups constructed for resistance to MLN using 500 SNPs in seven F3 populations; Table S2: Twenty-three common QTLs and specific SNP detected in both linkage mapping and joint linkage association mapping across seven F3 populations and their associated candidate genes.

Author Contributions

Conceptualization and methodology, L.A.O.A.; software and data analyses, L.A.O.A. and M.G.; validation, M.G., P.T., E.D. and B.E.I.; investigation, L.A.O.A.; resources, B.M.P. and M.G.; writing—original draft, L.A.O.A.; writing—review and editing, B.E.I., P.T., E.D., P.W.M.-D., V.O., Y.B., S.L.M., M.B.J., M.O., S.M., B.M.P. and M.G.; supervision, B.E.I., P.T., E.D., P.W.M.-D. and M.G.; project administration, P.T., S.L.M. and M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was part of a PhD program funded by Intra-ACP Mobility Scheme. The research was supported by the Bill and Melinda Gates Foundation, the Howard G. Buffett Foundation, and the United States Agency for International Development (USAID) through Stress Tolerant Maize for Africa (STMA, Grant # OPP1134248), and the CGIAR Research Program MAIZE. MAIZE receives W1 and W2 support from the Governments of Australia, Belgium, Canada, China, France, India, Japan, Korea, Mexico, Netherlands, New Zealand, Norway, Sweden, Switzerland, U.K., U.S., and the World Bank. We thank CIMMYT technicians for data collection at experimental sites.

Acknowledgments

The study was managed by The Global maize Program, CIMMYT and University of Ghana.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mahuku, G.; Lockhart, B.E.; Wanjala, B.; Jones, M.W.; Kimunye, J.N.; Stewart, L.R.; Cassone, B.J.; Subramanian, S.; Nyasani, J.; Kusia, E.; et al. Maize lethal necrosis (MLN), an emerging threat to maize-based food security in sub-Saharan Africa. Phytopathol. 2015, 105, 956–965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wamaitha, M.J.; Nigam, D.; Maina, S.; Stomeo, F.; Wangai, A.; Njuguna, J.N.; Holton, T.A.; Wanjala, B.W.; Wamalwa, M.; Lucas, T.; et al. Metagenomic analysis of viruses associated with maize lethal necrosis in Kenya. Virol. J. 2018, 15, 90. [Google Scholar] [CrossRef] [PubMed]
  3. Wangai, A.W.; Redinbaugh, M.G.; Kinyua, Z.M.; Miano, D.W.; Leley, P.K.; Kasina, M.; Mahuku, G.; Scheets, K.; Jeffers, D. First Report of Maize chlorotic mottle virus and Maize Lethal Necrosis in Kenya. Plant Dis. 2012, 96, 1582. [Google Scholar] [CrossRef] [PubMed]
  4. Redinbaugh, M.; Pratt, R. Virus resistance. In Handbook of Maize: It’s Biology; Bennetzen, J.L., Hake, S.C., Eds.; Springer: New York, NY, USA, 2009; pp. 251–268. [Google Scholar]
  5. Bonamico, N.C.; Di Renzzo, M.; Ibanez, M.; Borghi, M.; Diaz, D.; Salerno, J.; Balzarini, M. QTL analysis of resistance to Mal de Río Cuarto disease in maize using recombinant inbred lines QTL analysis of resistance to Mal de Río Cuarto disease in maize using recombinant inbred lines. J. Agric. Sci. 2012, 150, 5. [Google Scholar] [CrossRef]
  6. Zambrano, J.L.; Jones, M.W.; Brenner, E.; Francis, D.M.; Tomas, A.; Redinbaugh, M.G. Genetic analysis of resistance to six virus diseases in a multiple virus-resistant maize inbred line. Theor. Appl. Genet. 2014, 127, 867–880. [Google Scholar] [CrossRef]
  7. Gowda, M.; Beyene, Y.; Makumbi, D.; Semagn, K.; Olsen, M.S.; Bright, J.M.; Das, B.; Mugo, S.; Suresh, L.M.; Prasanna, B.M. Discovery and validation of genomic regions associated with resistance to maize lethal necrosis in four biparental populations. Mol. Breed. 2018, 38, 66. [Google Scholar] [CrossRef] [Green Version]
  8. Sitonik, C.; Suresh, L.M.; Beyene, Y.; Olsen, M.S.; Makumbi, D.; Oliver, K.; Das, B.; Bright, J.M.; Mugo, S.; Crossa, J.; et al. Genetic architecture of maize chlorotic mottle virus and maize lethal necrosis through GWAS, linkage analysis and genomic prediction in tropical maize germplasm. Theor. Appl. Genet. 2019, 132, 2381–2399. [Google Scholar] [CrossRef] [Green Version]
  9. Yu, J.; Holland, J.B.; McMullen, M.D.; Buckler, E.S. Genetic Design and Statistical Power of Nested Association Mapping in Maize. Genet 2008, 178, 539–551. [Google Scholar] [CrossRef] [Green Version]
  10. Wurschum, T.; Liu, W.; Gowda, M.; Maurer, H.P.; Fischer, S.; Schechert, A.; Reif, J.C. Comparison of biometrical models for joint linkage association mapping. Heredity 2012, 108, 332–340. [Google Scholar] [CrossRef] [Green Version]
  11. Hawkins, L.K.; Warburton, M.L.; Tang, J.D.; Tomashek, J.; Oliveira, D.A.; Ogunola, O.F.; Smith, J.S.; Williams, W.P. Survey of Candidate Genes for Maize Resistance to Infection by Aspergillus flavus and/or Aflatoxin Contamination. Toxins 2018, 10, 61. [Google Scholar] [CrossRef] [Green Version]
  12. Gelli, M.; Konda, A.R.; Liu, K.; Zhang, C.; Clemente, T.E.; Holding, D.R.; Dweikat, I.M. Validation of QTL mapping and transcriptome profiling for identification of candidate genes associated with nitrogen stress tolerance in sorghum. BMC Plant Boil. 2017, 17, 123. [Google Scholar] [CrossRef] [PubMed]
  13. Sudheesh, S.; Rodda, M.S.; Davidson, J.; Javid, M.; Stephens, A.; Slater, A.T.; Cogan, N.O.I.; Forster, J.W.; Kaur, S. SNP-Based Linkage Mapping for Validation of QTLs for Resistance to Ascochyta Blight in Lentil. Front. Plant Sci. 2016, 7, 311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gajjar, K.N.; Mishra, G.P.; Radhakrishnan, T.; Dodia, S.M.; Rathnakumar, A.L.; Kumar, N.; Kumar, K.; Dobaria, J.R.; Kumar, A. Validation of SSR markers linked to the rust and late leaf spot diseases resistance in diverse peanut genotypes. Aust. J. Crop. Sci. 2014, 8, 927–936. [Google Scholar]
  15. Landi, P.; Sanguineti, M.C.; Salvi, S.; Giuliani, S.; Bellotti, M.; Maccaferri, M.; Conti, S.; Tuberosa, R. Validation and characterization of a major QTL affecting leaf ABA concentration in maize. Mol. Breed. 2005, 15, 291–303. [Google Scholar] [CrossRef]
  16. Sukruth, M.; Paratwagh, S.A.; Sujay, V.; Kumari, V.; Gowda, M.V.C.; Nadaf, H.L.; Motagi, B.N.; Lingaraju, S.; Pandey, M.K.; Varshney, R.K.; et al. Validation of markers linked to late leaf spot and rust resistance, and selection of superior genotypes among diverse recombinant inbred lines and backcross lines in peanut (Arachis hypogaea L.). Euphytica 2015, 204, 343–351. [Google Scholar] [CrossRef]
  17. Zhou, H.; Liu, S.; Liu, Y.; Liu, Y.; You, J.; Deng, M.; Ma, J.; Chen, G.; Wei, Y.; Liu, C.; et al. Mapping and validation of major quantitative trait loci for kernel length in wild barley (Hordeum vulgare ssp. spontaneum). BMC Genet. 2016, 17, 130. [Google Scholar] [CrossRef] [Green Version]
  18. De Leon, T.B.; Linscombe, S.; Subudhi, P.K. Identification and validation of QTLs for seedling salinity tolerance in introgression lines of a salt tolerant rice landrace ‘Pokkali’. PLoS ONE 2017, 12, e0175361. [Google Scholar] [CrossRef]
  19. Ishikawa, G.; Nakamura, K.; Ito, H.; Saito, M.; Sato, M.; Jinno, H.; Yoshimura, Y.; Nishimura, T.; Maejima, H.; Uehara, Y.; et al. Association Mapping and Validation of QTLs for Flour Yield in the Soft Winter Wheat Variety Kitahonami. PLoS ONE 2014, 9, e111337. [Google Scholar] [CrossRef]
  20. Bernardo, R. Molecular Markers and Selection for Complex Traits in Plants: Learning from the Last 20 Years. Crop. Sci. 2008, 48, 1649. [Google Scholar] [CrossRef] [Green Version]
  21. Tuberosa, R.; Salvi, S.; Sanguineti, M.C.; Landi, P.; Maccaferri, M.; Conti, S. Mapping QTLs Regulating Morpho-physiological Traits and Yield: Case Studies, Shortcomings and Perspectives in Drought-stressed Maize. Ann. Bot. 2002, 89, 941–963. [Google Scholar] [CrossRef] [Green Version]
  22. Gowda, M.; Das, B.; Makumbi, D.; Babu, R.; Semagn, K.; Mahuku, G.; Olsen, M.S.; Bright, J.M.; Beyene, Y.; Prasanna, B.M. Genome-wide association and genomic prediction of resistance to maize lethal necrosis disease in tropical maize germplasm. Theor. Appl. Genet. 2015, 128, 1957–1968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Tivoli, B.; Baranger, A.; Avila, C.M.; Banniza, S.; Barbetti, M.; Chen, W.; Davidson, J.; Lindeck, K.; Kharrat, M.; Rubiales, D.; et al. Screening techniques and sources of resistance to foliar diseases caused by major necrotrophic fungi in grain legumes. Euphytica 2006, 147, 223–253. [Google Scholar] [CrossRef]
  24. Wang, J.; Li, H.; Zhang, L.; Meng, L. Users’ Manual of QTL IciMapping, 4.1; International Maize and Wheat Improvement Center (CIMMYT): Texcoco, Mexico, 2016. [Google Scholar]
  25. Alvarado, G.; López, M.; Vargas, M.; Pacheco, A.; Rodriguez, F.; Burgueno, J.; Crossa, J. META-R (Multi Environment Trail Analysis with R for Windows) Version 5.0—CIMMYT Research Software Dataverse—CIMMYT Dataverse Network. 2015. Available online: http://hdl.handle.net/11529/10201 (accessed on 5 September 2019).
  26. Stram, D.O.; Lee, J.W. Variance Components Testing in the Longitudinal Mixed Effects Model. Biometrics 1994, 50, 1171–1177. [Google Scholar] [CrossRef] [PubMed]
  27. CIMMYT. Laboratoty Protocols: CIMMYT Applied and Molecular Genetics Laboratoty, 3rd ed.; CIMMYT: Batan, Mexico, 2005. [Google Scholar]
  28. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  29. Xia, X.; Melchinger, A.E.; Kuntze, L.; Lübberstedt, T. Quantitative Trait Loci Mapping of Resistance to Sugarcane Mosaic Virus in Maize. Phytopathology 1999, 89, 660–667. [Google Scholar] [CrossRef] [Green Version]
  30. Horn, F.; Habekuss, A.; Stich, B.; Habekuß, A. Linkage mapping of Barley yellow dwarf virus resistance in connected populations of maize. BMC Plant Boil. 2015, 15, 29. [Google Scholar] [CrossRef] [Green Version]
  31. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop. J. 2015, 3, 269–283. [Google Scholar] [CrossRef] [Green Version]
  32. Kosambi, D. The estimation of map distances from recombination values. Ann. Eugen. 1944, 12, 172–175. [Google Scholar] [CrossRef]
  33. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [Green Version]
  34. Li, H.; Ye, G.; Wang, J. A modified algorithm for the improvement of composite interval mapping. Genetics 2007, 175, 361–374. [Google Scholar] [CrossRef] [Green Version]
  35. Lübberstedt, T.; Melchinger, A.E.; Schön, C.C.; Utz, H.F.; Klein, D. QTL Mapping in Testcrosses of European Flint Lines of Maize: I. Comparison of Different Testers for Forage Yield Traits. Crop. Sci. 1997, 37, 921. [Google Scholar] [CrossRef]
  36. Liu, Y.; Chen, L.; Fu, N.; Lou, Q.; Mei, H.; Xiong, L.; Li, M.; Xu, X.; Mei, X.; Luo, L. Dissection of additive, epistatic effect and QTL × environment interaction of quantitative trait loci for sheath blight resistance in rice. Hereditas 2014, 151, 28–37. [Google Scholar] [CrossRef] [PubMed]
  37. Stuber, C.W.; Edwards, M.D.; Wendel, J.F. Molecular Marker-Facilitated Investigations of Quantitative Trait Loci in Maize. II. Factors Influencing Yield and its Component Traits1. Crop. Sci. 1987, 27, 639. [Google Scholar] [CrossRef] [Green Version]
  38. Carlson, C.H.; Gouker, F.; Crowell, C.R.; Evans, L.; DiFazio, S.P.; Smart, C.D.; Smart, L.B. Joint linkage and association mapping of complex traits in shrub willow (Salix purpurea L.). Ann. Bot. 2019, 124, 701–715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Mammadov, J.; Sun, X.; Gao, Y.; Ochsenfeld, C.; Bakker, E.; Ren, R.; Flora, J.; Wang, X.; Kumpatla, S.; Meyer, D.; et al. Combining powers of linkage and association mapping for precise dissection of QTL controlling resistance to gray leaf spot disease in maize (Zea mays L.). BMC Genom. 2015, 16, 916. [Google Scholar] [CrossRef] [Green Version]
  40. Schwarz, G. Estimating the Dimension of a Model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  41. SAS Institute Inc. SAS® 9.4 Intelligence Platform: Data Administration Guide, 6th ed.; SAS Institute Inc.: Cary, NC, USA, 2016. [Google Scholar]
  42. Reif, J.C.; Zhao, Y.; Würschum, T.; Gowda, M.; Hahn, V. Genomic prediction of sunflower hybrid performance. Plant Breed. 2013, 132, 107–114. [Google Scholar] [CrossRef]
  43. Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 1979, 6, 65–70. [Google Scholar]
  44. Chen, J.; Zavala, C.; Ortega, N.; Petroli, C.; Franco, J.; Burgueño, J.; Costich, D.E.; Hearne, S.J. The Development of Quality Control Genotyping Approaches: A Case Study Using Elite Maize Lines. PLoS ONE 2016, 11, e0157236. [Google Scholar] [CrossRef]
  45. Barría, A.; Christensen, K.A.; Yoshida, G.M.; Correa, K.; Jedlicki, A.; Lhorente, J.P.; Davidson, W.S.; Yanez, J.M. Genomic Predictions and Genome-Wide Association Study of Resistance Against Piscirickettsia salmonis in Coho Salmon (Oncorhynchus kisutch) Using ddRAD Sequencing. G3: Genes|Genomes|Genetics 2018, 8, 1183–1194. [Google Scholar] [CrossRef] [Green Version]
  46. Gao, N.; Teng, J.; Ye, S.; Yuan, X.; Huang, S.; Zhang, H.; Zhang, X.; Li, J.; Zhang, Z. Genomic Prediction of Complex Phenotypes Using Genic Similarity Based Relatedness Matrix. Front. Genet. 2018, 9, 364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Hochholdinger, F.; Hoecker, N. Towards the molecular basis of heterosis. Trends Plant Sci. 2007, 12, 427–432. [Google Scholar] [CrossRef] [PubMed]
  48. Beyene, Y.; Gowda, M.; Suresh, L.M.; Mugo, S.; Olsen, M.; Oikeh, S.O.; Juma, C.; Tarekegne, A.; Prasanna, B.M. Genetic analysis of tropical maize inbred lines for resistance to maize lethal necrosis disease. Euphytica 2017, 213, 224. [Google Scholar] [CrossRef] [Green Version]
  49. Springer, N.M.; Stupar, R.M. Allelic variation and heterosis in maize: How do two halves make more than a whole? Genome Res. 2007, 17, 264–275. [Google Scholar] [CrossRef] [Green Version]
  50. Karume, K.; Nazir, N.; Salmin, I. Assessment of heterosis and combining ability in maize (Zea mays L.) for maize lethal necrosis disease. Afr. J. Plant Breed. 2017, 4, 104–200. [Google Scholar] [CrossRef] [Green Version]
  51. Mbega, E.; Ndakidemi, P.; Mamiro, D.; Mushongi, A.; Kitenge, K.; Ndomba, O. Role of Potyviruses in Synergistic Interaction Leading to Maize Lethal Necrotic Disease on Maize. Int. J. Curr. Microbiol. Appl. Sci. 2016, 5, 85–96. [Google Scholar] [CrossRef]
  52. Mackay, T.F.C. The genetic architecture of quantitative traits: Lessons from Drosophila, Current Opinion in Genetics and Development. Curr. Opin. Genet. Dev. 2004, 14, 253–257. [Google Scholar] [CrossRef]
  53. Forabosco, P.; Falchi, M.; Devoto, M. Statistical tools for linkage analysis and genetic association studies. Expert Rev. Mol. Diagn. 2005, 5, 781–796. [Google Scholar] [CrossRef]
  54. Gallois, J.-L.; Moury, B.; German-Retana, S. Role of the Genetic Background in Resistance to Plant Viruses. Int. J. Mol. Sci. 2018, 19, 2856. [Google Scholar] [CrossRef] [Green Version]
  55. Jansen, R.C.; Stam, P. High Resolution of Quantitative Traits into Multiple Loci via Interval Mapping. Genet. 1994, 136, 1447–1455. [Google Scholar]
  56. Jones, M.W.; Penning, B.W.; Jamann, T.M.; Glaubitz, J.C.; Romay, C.; Buckler, E.S.; Redinbaugh, M.G. Diverse chromosomal locations of quantitative trait loci for tolerance to maize chlorotic mottle virus in five maize populations. Phytopathology 2018, 108, 748–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Gustafson, T.J.; de Leon, N.; Kaeppler, S.M.; Tracy, W.F. Genetic analysis of sugarcane mosaic virus resistance in the Wisconsin diversity panel of maize. Crop. Sci. 2018, 58, 1853–1865. [Google Scholar] [CrossRef]
Figure 1. Phenotypic distribution of MLN disease severity (MLN-DS) and the area under disease progress curve (AUDPC) values recorded in seven individual and across F3 populations. Arrows indicate the performance of parents.
Figure 1. Phenotypic distribution of MLN disease severity (MLN-DS) and the area under disease progress curve (AUDPC) values recorded in seven individual and across F3 populations. Arrows indicate the performance of parents.
Genes 11 00032 g001
Figure 2. PCA (principal component analysis) biplot showing structures of the seven F3 populations.
Figure 2. PCA (principal component analysis) biplot showing structures of the seven F3 populations.
Genes 11 00032 g002
Figure 3. Box plots showing the phenotypic values of the different allele classes of two major-effect SNPs identified in five F3 populations and across populations for MLN-DS and AUDPC value. The SNP names, alleles and the specific F3 population where the effect is witnessed are mentioned above. The black horizontal lines in the middle of the boxes are the median values for the MLN-DS and AUDPC value in the respective allele classes. The vertical size of the boxes represents the inter-quantile range. The upper and lower whiskers represent the minimum and maximum values of data.
Figure 3. Box plots showing the phenotypic values of the different allele classes of two major-effect SNPs identified in five F3 populations and across populations for MLN-DS and AUDPC value. The SNP names, alleles and the specific F3 population where the effect is witnessed are mentioned above. The black horizontal lines in the middle of the boxes are the median values for the MLN-DS and AUDPC value in the respective allele classes. The vertical size of the boxes represents the inter-quantile range. The upper and lower whiskers represent the minimum and maximum values of data.
Genes 11 00032 g003
Figure 4. QTL for MLN resistance in the consensus linkage map of seven bi-parental populations. Different colors represent QTL detected from different populations.
Figure 4. QTL for MLN resistance in the consensus linkage map of seven bi-parental populations. Different colors represent QTL detected from different populations.
Genes 11 00032 g004
Figure 5. Box-whisker-plots for the accuracy of genomic predictions assessed by five-fold cross-validation within each population for MLN-DS (white box) and AUDPC values (grey box).
Figure 5. Box-whisker-plots for the accuracy of genomic predictions assessed by five-fold cross-validation within each population for MLN-DS (white box) and AUDPC values (grey box).
Genes 11 00032 g005
Table 1. Analysis of variance for MLN disease severity (MLN-DS) and area under disease progress curve (AUDPC) values using seven segregating F3 populations evaluated for three seasons under MLN inoculated plots.
Table 1. Analysis of variance for MLN disease severity (MLN-DS) and area under disease progress curve (AUDPC) values using seven segregating F3 populations evaluated for three seasons under MLN inoculated plots.
TraitMean (Range)σ2Gσ2e *H2
CKDHL120918 × CML494 (F3 pop1)
MLN-DS4.52 (3.06–7.12)0.41 **0.440.74
AUDPC131.9 (92.9–185.8)272.20 **230.730.78
CML543 × CML494 (F3 pop2)
MLN-DS5.11 (3.86–6.35)0.17 **0.320.62
AUDPC140.8 (103.8–168.9)141.05 **177.680.70
CKDHL120918 × CML543 (F3 pop3)
MLN-DS4.52 (3.17–6.86)0.46 **0.500.73
AUDPC126.2 (82.1–201.0)282.50 **280.870.75
CKLTI0227 × CKDHL120918 (F3 pop4)
MLN-DS4.60 (2.80–7.65)0.69 **0.570.79
AUDPC129.1 (92.6–194.7)467.78 **299.860.82
CKDHL0089 × CKDHL120918 (F3 pop5)
MLN-DS4.94 (2.85–8.20)0.65 **0.720.73
AUDPC133.9 (79.9–209.3)411.12 **379.390.76
CKDHL0221 × CKDHL120312 (F3 pop6)
MLN-DS4.40 (2.56–6.97)0.50 **0.490.75
AUDPC125.6 (78.0–187.9)374.25 **245.350.82
CKDHL0089 × CML494 (F3 pop7)
MLN-DS4.82 (2.94–7.31)0.44 **0.510.72
AUDPC134.4 (91.0–191.8)300.75 **260.510.78
Across seven populations
MLN-DS4.70 (2.56–8.20)0.40 **0.500.70
AUDPC132.4 (78.0–209.3)265.48 **267.900.75
* and ** indicate significance at P < 0.05 and P < 0.01, respectively MLN-DS = MLN disease severity 42 dpi; AUDPC = area under disease progress curve; σ2G = genotypic variance; σ2e* = error variance confounded with GxE variance; and H2 = broad sense heritability.
Table 2. Quantitative trait loci (QTL) detected by integrated composite intervals mapping analysis for resistance to MLN in seven F3 populations evaluated in MLN inoculated plots over three cropping seasons.
Table 2. Quantitative trait loci (QTL) detected by integrated composite intervals mapping analysis for resistance to MLN in seven F3 populations evaluated in MLN inoculated plots over three cropping seasons.
Trait NameQTL Name aChrPosition (cM)LODPVE (%)AddDomNature of QTLTotal PVE (%)Marker NamePhysical Position (bp)Fav Parent
Left MRight MLeft MRight M
CKDHL120918 × CML494
MLN-DSqMLN1_26511833.992.58−0.46−0.26PD30.08PZB00648_5d8_317,595,139265,199,938CKDHL120918
qMLN2_156234.132.34−0.041.09ODPZA01232_1PHM3055_9155,868,024192,602,324CKDHL120918
qMLN2_1021223.160.690.07−0.23ODPZA00620_3PZA00365_210,429,4051,221,385CML494
qMLN3_14232323.924.270.37−0.21PDPZA00920_1PHM15449_10142,821,031125,077,922CML494
qMLN4_1504904.412.53−0.020.70ODS4_153520131S4_149896839153,520,131149,896,839CKDHL120918
qMLN5_1775884.242.580.030.83ODPZA01410_1S5_177634071172,682,963177,634,071CML494
qMLN6_8567814.113.980.34−0.20PDPZB01009_1PHM8909_1284,664,84091,883,155CML494
qMLN6_8561326.254.560.190.76ODPZB01009_1PHM8909_1284,664,840918,83,155CML494
qMLN7_1587874.712.48−0.071.11ODPHM7898_10S7_157472460161,993,743157,472,460CKDHL120918
qMLN9_1429663.641.02−0.11−0.21ODPZB00221_3PHM229_15142,271,04730,003,189CKDHL120918
qMLN10_91043.822.430.080.84ODPZA00866_2PHM5740_9124,203,56587,73,358CML494
AUDPCqMLN2_1021213.033.231.54−6.16OD34.72PZA00620_3PZA00365_210,429,4051,221,385CML494
qMLN3_14232312.688.925.50−5.57DOPZA00920_1PHM15449_10142,821,031125,077,922CML494
qMLN4_14341402.791.71−2.781.27PDPHM1505_31S4_9850443143,162,7459,850,443CKDHL120918
qMLN6_8567817.9621.649.80−5.38PDPZB01009_1PHM8909_1284,664,84091,883,155CML494
qMLN9_1429684.796.25−4.08−5.16ODPZB00221_3PHM229_15142,271,04730,003,189CKDHL120918
CML543 × CML494
MLN-DSqMLN1_471433.411.10−0.070.04PD43.49csu1138_4S1_46411896119,018,55646,411,896CML494
qMLN3_14631411.794.08−0.720.25PDS3_146966676S3_146363360146,966,676146,363,360CML494
qMLN4_304142.509.94−0.220.52ODPZA02457_1bt2_729,031,20066,290,994CML494
qMLN5_160533.521.22−0.080.04PDS5_42297152PZA01796_142,297,152160,321,846CML494
qMLN6_21645.031.640.09−0.05PDS6_21007530PZA03063_2121,007,53025,335,225CML543
qMLN8_108843.4212.37−0.210.53ODS8_10001165PZA02388_110,001,165169,137CML494
qMLN9_1429331.2112.720.70−0.08ADPZA00832_1PHM7916_4147,131,097132,762,904CML543
AUDPCqMLN1_471426.141.64−2.931.62PD49.84csu1138_4S1_46411896119,018,55646,411,896CML494
qMLN1_26511542.974.17−7.589.26ODd8_3umc128_2265,199,938227,602,208CML494
qMLN3_14631415.654.18−24.017.57PDS3_146966676S3_146363360146,966,676146,363,360CML494
qMLN3_14231073.771.172.730.84PDPZA00920_1PZA02402_1142,821,031169,771,952CML543
qMLN4_304143.887.99−6.9115.27ODPZA02457_1bt2_729,031,20066,290,994CML494
qMLN5_160533.911.01−2.480.94PDS5_42297152PZA01796_142,297,152160,321,846CML494
qMLN6_21628.952.383.84−0.71ADS6_21007530PZA03063_2121,007,53025,335,225CML543
qMLN7_15871194.818.135.12−18.38ODPHM7898_10PHM1912_23161,993,743155,970,264CML543
qMLN8_108853.055.28−6.739.00ODS8_10001165PZA02388_110,001,165169,137CML494
qMLN9_1429334.1810.5421.95−4.51PDPZA00832_1PHM7916_4147,131,097132,762,904CML543
CKDHL120918 × CML543
MLN-DSqMLN1_2651144.122.79−0.04−0.21OD47.88d8_3PHM14475_7265,199,938256,245,118CKDHL120918
qMLN1_2521283.313.08−0.12−0.15ODPZA02269_4PHM4942_12252,721,946226,461,786CKDHL120918
qMLN2_22504.303.45−0.05−0.23ODPHM13440_13PZA00365_22,527,3441,221,385CKDHL120918
qMLN3_14237333.3328.840.46−0.26PDS3_146966676S3_68596995146,966,67668,596,995CML543
qMLN4_14341204.082.870.16−0.04PDbt2_7PHM1505_3166,290,994143,162,745CML543
qMLN5_20251122.894.43−0.14−0.19ODPZB00765_1PHM5484_22202,174,58521,449,633CKDHL120918
qMLN7_158712.623.15−0.39−0.44DOS7_157472460S7_137455469157,472,460137,455,469CKDHL120918
qMLN10_910816.775.10−0.15−0.22ODPHM5740_9PZA01642_18,773,35814,703,451CKDHL120918
AUDPCqMLN1_2651144.852.50−1.32−5.28OD52.80d8_3PHM14475_7265,199,938256,245,118CKDHL120918
qMLN1_2521283.442.40−3.01−3.52DOPZA02269_4PHM4942_12252,721,946226,461,786CKDHL120918
qMLN2_22504.032.46−0.78−5.45ODPHM13440_13PZA00365_22,527,3441,221,385CKDHL120918
qMLN3_154332.935.2310.02−11.50DOS3_154250438S3_150836832154,250,438150,836,832CML543
qMLN3_14237417.6811.097.15−6.24DOS3_146966676S3_68596995146,966,67668,596,995CML543
qMLN3_20731092.732.02−3.50−0.77PDPZA00538_15PZA01931_17206,889,707227,682,081CKDHL120918
qMLN4_14341206.033.304.58−0.86ADbt2_7PHM1505_3166,290,994143,162,745CML543
qMLN5_20251123.163.65−3.50−4.49ODPZB00765_1PHM5484_22202,174,58521,449,633CKDHL120918
qMLN6_15761183.834.7911.14−8.68PDPZA01618_2S6_156386857129,927,781156,386,857CML543
qMLN7_158713.853.30−10.90−11.74DOS7_157472460S7_137455469157,472,460137,455,469CKDHL120918
qMLN9_1099292.811.63−3.06−1.67PDPHM229_15S9_10954923030,003,189109,549,230CKDHL120918
qMLN10_411022.903.469.86−10.22DOPHM4066_11PHM1956_9041,187,56540,187,565CML543
qMLN10_910827.454.61−3.83−5.70ODPHM5740_9PZA01642_18,773,35814,703,451CKDHL120918
CKLTI0227 × CKDHL120918
MLN-DSqMLN1_2821993.211.550.69−0.60DO25.13PHM4752_14PZA03020_8298,874,066282,044,048CKLTI0227
qMLN1_26511093.521.560.68−0.63DOPZA03020_8PZA01254_2282,044,048106,204,446CKLTI0227
qMLN1_4711394.791.380.76−0.82DOS1_22744948S1_8715893022,744,94887,158,930CKLTI0227
qMLN1_2712023.540.910.93−0.94DOS1_15353866S1_269396815,353,8662,693,968CKLTI0227
qMLN3_1423374.441.450.73−0.35PDS3_146966676S3_55444954146,966,67655,444,954CKLTI0227
qMLN3_1303655.681.480.70−0.82DOS3_92864540S3_13616556592,864,540136,165,565CKLTI0227
qMLN3_13031215.321.49−0.68−1.08ODPHM2343_25S3_15425043827,981,649154,250,438CKDHL120918
qMLN4_304364.360.350.72−0.05ADPZA02457_1PZA00726_1029,031,20060,768,063CKLTI0227
qMLN4_1504964.291.590.59−0.74ODS4_155296684S4_9741874155,296,6849,741,874CKLTI0227
qMLN5_16051304.121.050.90−0.83DOS5_154350617S5_198716574154,350,617198,716,574CKLTI0227
qMLN5_4251613.490.25−0.08−1.02ODS5_42297152PZA02164_1642,297,152112,179,855CKDHL120918
qMLN6_1576205.871.460.59−1.12ODS6_156386857PHM4748_16156,386,857158,540,019CKLTI0227
qMLN6_856297.171.680.72−0.34PDPHM4748_16PZB01009_1158,540,01984,664,840CKLTI0227
qMLN6_856838.861.790.81−0.43PDPZB01009_1S6_8982377284,664,84089,823,772CKLTI0227
qMLN6_15761164.731.100.071.37ODPHM2658_129PZA01884_1164,999,578132,316,835CKLTI0227
qMLN7_15871453.561.080.86−0.99DOS7_167230991S7_127970539167,230,991127,970,539CKLTI0227
qMLN8_1081004.341.370.73−0.90ODPZA02746_2PZA02388_1163,067,200169,137CKLTI0227
qMLN9_109993.191.64−0.63−0.57DOPZA00708_3S9_109549230147,381,231109,549,230CKDHL120918
qMLN10_11410365.741.610.69−0.83ODS10_113832226PZA01001_2113,832,226146,538,889CKLTI0227
AUDPCqMLN1_26511636.851.3524.93−8.24PD28.69d8_3PZA02269_4265,199,938252,721,946CKLTI0227
qMLN2_412212.202.3140.41−11.93PDPHM10404_8PZA02450_140,967,99147,575,949CKLTI0227
qMLN3_1423345.673.4214.61−8.98PDS3_146966676S3_55444954146,966,67655,444,954CKLTI0227
qMLN3_146311310.791.98−40.42−5.50ADS3_146250249S3_146026612146,250,249146,026,612CKDHL120918
qMLN3_13031253.493.8112.36−18.09ODPHM2343_25S3_15425043827,981,649154,250,438CKLTI0227
qMLN4_304773.292.4419.01−20.07DOS4_6601124PZA02509_156,601,1243,904,858CKLTI0227
qMLN6_1576328.843.7716.09−5.31PDPHM4748_16PZB01009_1158,540,01984,664,840CKLTI0227
qMLN9_147963.283.17−15.11−9.92PDS9_146012201PZA00708_3146,012,201147,381,231CKDHL120918
qMLN9_1099213.192.6618.29−19.22DOPZA00708_3S9_109549230147,381,231109,549,230CKLTI0227
qMLN10_11410363.543.2016.73−18.91DOS10_113832226PZA01001_2113,832,226146,538,889CKLTI0227
CKDHL0089 × CKDHL120918
MLN-DSqMLN1_18144.427.15−0.110.87OD54.34S1_18838432PZA00175_218,838,4328,510,027CKDHL120918
qMLN2_16921684.256.61−0.061.07ODPZA02727_1PZA00515_10227,921,381169,265,278CKDHL120918
qMLN3_14237344.1424.440.62−0.10ADPZA00920_1S3_133048570142,821,031133,048,570CKDHL0089
qMLN5_19051283.751.56−0.14−0.06PDS5_190675983S5_201226926190,675,983201,226,926CKDHL120918
qMLN6_1576203.218.200.31−0.41ODS6_156386857PHM3466_69156,386,857167,148,384CKDHL0089
qMLN6_8561264.161.780.160.04PDPZA00942_2PHM8909_12102,566,00091,883,155CKDHL0089
qMLN10_910412.066.32−0.28−0.10PDPZA01313_2PHM5740_93,598,2628,773,358CKDHL120918
AUDPCqMLN1_181103.698.65−3.6715.28OD57.62S1_18838432PZA00175_218,838,4328,510,027CKDHL120918
qMLN2_16921663.186.28−1.8121.55ODPZA02727_1PZA00515_10227,921,381169,265,278CKDHL120918
qMLN3_14237350.0227.4616.42−2.21ADPZA00920_1S3_133048570142,821,031133,048,570CKDHL0089
qMLN5_19051275.492.11−4.27−1.58PDPHM7908_25S5_190675983191,075,472190,675,983CKDHL120918
qMLN6_1576212.968.127.34−10.37ODS6_156386857PHM3466_69156,386,857167,148,384CKDHL0089
qMLN6_8561335.031.964.230.36ADPHM8909_12PZA00427_391,883,15579,815,961CKDHL0089
qMLN10_910311.235.20−6.46−2.44PDPZA01313_2PHM5740_93,598,2628,773,358CKDHL120918
CKDHL0221 × CKDHL120312
MLN-DSqMLN1_471715.323.74−0.160.07PD52.13PZA00447_8S1_464118969,024,00546,411,896CKDHL120312
qMLN3_13035644.8344.510.56−0.11ADPZA02402_1PHM15449_10169,771,952125,077,922CKDHL0221
qMLN4_1504372.662.74−0.13−0.09PDPZA01187_1PHM1505_31177,666,738143,162,745CKDHL120312
qMLN4_741284.633.230.150.01ADS4_6544767PHM16788_66,544,76713,581,955CKDHL0221
qMLN8_108454.263.120.150.03ADPHM5235_8PZA00368_194,414,9785,632,308CKDHL0221
AUDPCqMLN1_471714.974.38−4.051.41PD59.07PZA00447_8S1_464118969,024,00546,411,896CKDHL120312
qMLN1_4711002.982.492.21−2.90ODS1_46411896PHM12323_1746,411,89653,357,797CKDHL0221
qMLN3_13035624.6426.009.70−2.76PDPZA02402_1PHM15449_10169,771,952125,077,922CKDHL0221
qMLN3_14236313.1912.716.76−1.85PDPZA00279_2PZA00920_152,804,070142,821,031CKDHL0221
qMLN4_1504393.655.24−4.27−2.01PDPZA01187_1PHM1505_31177,666,738143,162,745CKDHL120312
qMLN4_741283.953.443.710.18ADS4_6544767PHM16788_66,544,76713,581,955CKDHL0221
qMLN5_425342.632.36−0.23−4.17ODPHM16854_3PZA00522_1234,587,02957,933,548CKDHL120312
qMLN8_108455.064.704.16−0.42ADPHM5235_8PZA00368_194,414,9785,632,308CKDHL0221
qMLN10_11410435.164.473.72−2.52PDPZA00814_1PHM1576_2587,194,491124,203,168CKDHL0221
CKDHL0089 × CML494
MLN-DSqMLN5_1905894.405.48−0.170.03AD46.74PZA01427_1PHM7908_2523,135,578191,075,472CML494
qMLN5_20251118.947.97−0.20−0.03ADS5_200938637PHM563_9200,938,637204,993,639CML494
qMLN6_15761073.378.420.090.71ODS6_157568432S6_156386857157,568,432156,386,857CKDHL0089
AUDPCqMLN3_130313.026.640.75−4.18OD50.87PZA01447_1S3_13304857053,549,251133,048,570CKDHL0089
qMLN5_1905895.1415.91−4.780.81ADPZA01427_1PHM7908_2523,135,578191,075,472CML494
qMLN5_20251117.6216.65−4.81−0.76ADS5_200938637PHM563_920,0938,637204,993,639CML494
qMLN6_15761083.2014.952.3416.81ODS6_157568432S6_156386857157,568,432156,386,857CKDHL0089
LOD = logarithm of odds; Add = additive effect; Dom = dominance effect; PVE = phenotypic variance explained; fav parent = parental genotype from where favorable allele for MLN resistance is contributing; a QTL name composed by the trait code followed by the chromosome number in which the QTL was mapped and a physical position of the QTL.
Table 3. Joint linkage association mapping depicting allele substitution effects and total phenotypic variance explained (PVE) using segregating F3 progenies derived from seven bi-parental populations evaluated for three seasons under MLN inoculation.
Table 3. Joint linkage association mapping depicting allele substitution effects and total phenotypic variance explained (PVE) using segregating F3 progenies derived from seven bi-parental populations evaluated for three seasons under MLN inoculation.
MarkerQTL NameChrPosition (Mbp)Model AModel BModel C
α-EffectPVE (%)α-EffectPVE (%)α-EffectPVE (%)
PZA00447_8qMLN1_919.02−0.080.60
PHM5622_21qMLN1_1841183.83−0.120.40
S3_48493677qMLN3_48348.490.327.90−0.020.40
S3_55444954qMLN3_55355.440.140.100.101.00
S3_68596995qMLN3_68368.60−0.060.200.060.40
S3_92694873qMLN3_92392.69−0.280.70−0.021.80
S3_113429913qMLN3_1133113.43−0.381.00
PHM15449_10qMLN3_1253125.080.103.300.110.800.162.20
S3_148291047qMLN3_1483148.29−0.7210.20−0.664.70−0.161.40
S3_151342843qMLN3_1513151.34−0.233.20−0.423.30
PHM2919_23qMLN3_1993199.89−0.120.40−0.120.40
PZA00726_8qMLN4_60460.77−0.040.70
S4_235381719qMLN4_2354235.380.030.90
PHM565_31qMLN5_24524.24−0.030.00−0.470.30
S5_170023563qMLN5_1705170.020.010.10
PHM7908_25qMLN5_1915191.080.040.20
S5_196017729qMLN5_1965196.02−0.110.10−0.030.10−0.010.20
S5_202816906qMLN5_2025202.82−0.100.10
PHM563_9qMLN5_2045204.99−0.080.20
PZA03167_5qMLN5_2075207.600.320.30
S5_209467974qMLN5_2095209.47−0.070.30
S6_13300385qMLN6_13613.300.191.900.201.10
S6_86475982qMLN6_86686.48−0.271.00
S6_89823772qMLN6_90689.82−0.242.70−0.230.90−0.222.50
PHM5235_8qMLN8_94894.410.150.50
PZA01313_2qMLN10_4103.600.111.300.182.900.100.80
PHM5740_9qMLN10_9108.770.090.50
Total PVE (%) 34.40 27.30 29.10
PZA00447_8qMLN1_919.02−1.630.20
PHM5622_21qMLN1_1841183.83−3.540.40
S3_48493677qMLN3_48348.496.636.80
S3_55444954qMLN3_55355.443.531.001.940.60
S3_68596995qMLN3_68368.60−2.450.60−1.552.30
PHM15449_10qMLN3_1253125.083.011.003.653.00
S3_148291047qMLN3_1483148.29−20.1410.20−16.474.80−3.861.50
S3_151342843qMLN3_1513151.34−12.064.90−10.593.70−4.721.80
PHM2919_23qMLN3_1993199.89−4.350.50−4.040.80−1.700.80
PZA00726_8qMLN4_60460.77−1.640.70
S4_155378923qMLN4_1554155.38−8.031.00
S4_235381719qMLN4_2354235.38−10.740.80
S5_170164477qMLN5_1705170.165.530.5010.130.50
PHM7908_25qMLN5_1915191.081.380.20
S5_196017729qMLN5_1965196.02−2.140.20−1.130.40−4.370.20
S5_202816906qMLN5_2025202.82–2.530.10
PHM563_9qMLN5_2045204.99–5.610.30
PZA03167_5qMLN5_2075207.607.570.30
S5_209467974qMLN5_2095209.472.000.40
S6_13300385qMLN6_13613.306.304.305.691.400.541.00
S6_86475982qMLN6_86686.48–6.970.90
S6_89823772qMLN6_90689.82–6.580.90–6.141.00–5.812.00
S8_74144408qMLN8_74874.143.860.30
PHM5235_8qMLN8_94894.414.411.50
S8_102533570qMLN8_1028102.530.830.30
PZA01313_2qMLN10_4103.603.411.204.733.300.660.90
PHM5740_9qMLN10_9108.77–5.510.30
PZA00866_2qMLN10_12410124.201.710.601.440.401.790.90
Total PVE (%) 33.60 29.00 39.80

Share and Cite

MDPI and ACS Style

Awata, L.A.O.; Beyene, Y.; Gowda, M.; L. M., S.; Jumbo, M.B.; Tongoona, P.; Danquah, E.; Ifie, B.E.; Marchelo-Dragga, P.W.; Olsen, M.; et al. Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations. Genes 2020, 11, 32. https://doi.org/10.3390/genes11010032

AMA Style

Awata LAO, Beyene Y, Gowda M, L. M. S, Jumbo MB, Tongoona P, Danquah E, Ifie BE, Marchelo-Dragga PW, Olsen M, et al. Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations. Genes. 2020; 11(1):32. https://doi.org/10.3390/genes11010032

Chicago/Turabian Style

Awata, Luka A. O., Yoseph Beyene, Manje Gowda, Suresh L. M., McDonald B. Jumbo, Pangirayi Tongoona, Eric Danquah, Beatrice E. Ifie, Philip W. Marchelo-Dragga, Michael Olsen, and et al. 2020. "Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations" Genes 11, no. 1: 32. https://doi.org/10.3390/genes11010032

APA Style

Awata, L. A. O., Beyene, Y., Gowda, M., L. M., S., Jumbo, M. B., Tongoona, P., Danquah, E., Ifie, B. E., Marchelo-Dragga, P. W., Olsen, M., Ogugo, V., Mugo, S., & Prasanna, B. M. (2020). Genetic Analysis of QTL for Resistance to Maize Lethal Necrosis in Multiple Mapping Populations. Genes, 11(1), 32. https://doi.org/10.3390/genes11010032

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop