Genetic Diversity, Population Structure and Inter-Trait Relationships of Combined Heat and Drought Tolerant Early-Maturing Maize Inbred Lines from West and Central Africa

Adequate knowledge and understanding of the genetic diversity and inter-trait relationships among elite maize inbred lines are crucial for determining breeding strategies and predicting hybrid performance. The objectives of this study were to investigate the genetic diversity of 162 early maturing white and yellow tropical maize inbred lines, and to determine the population structure, heterotic groups and inter-trait relationships among the lines. Using 9684 DArT single nucleotide polymorphism (SNP) markers, a gene diversity (GD) of 0.30 was recorded for the inbred lines with polymorphic information content (PIC) ranging from 0.08 to 0.38. The genetic relatedness among the inbred lines evaluated revealed six different groups based on the history of selection, colour of endosperm and pedigree. The genotype-by-trait (GT) biplot analysis identified inbred 1 (TZEI 935) as outstanding in terms of combined heat and drought (HD) tolerance with the base index analysis identifying 15 superior inbreds in the HD environment. A wide range of genetic variability was observed among the inbred lines, indicating that they are an invaluable resource for breeding for HD tolerance in maize breeding programmes, especially in West and Central Africa.


Introduction
Molecular analyses in combination with morphological and agronomic evaluation of germplasm are useful for providing complementary information on the genetic diversity in a breeding programme [1]. Thus, monitoring of genetic diversity in the base material helps breeders to formulate breeding strategies by selecting parental combinations for advance selection [2], determining the level of genetic variability whilst identifying the main clusters selected for specific traits [3]. Information on genetic diversity is also useful in classifying inbred lines into heterotic groups [4,5], and predicting hybrid performance [6]. By definition, a heterotic group encompasses a set of directly related inbred lines which produce high yielding and vigorous hybrids when crossed to inbred lines from different heterotic groups, but not when crossed to inbred lines from the same heterotic group [7].
Additionally, the availability of information on the genetic diversity and population structure facilitates the assessment of the likelihood of loss in genetic diversity during preservation or selection [8], traits are very low. However, this problem could be mitigated through indirect selection of secondary traits that have high genetic correlations with grain yield coupled with high heritability and that are easily measured. Yan and Kang [26] proposed genotype-by-trait (GT) analysis for multiple trait-based assessment of genotypes. Such statistical tools aid in the identification of genotypes with specific desirable traits that could be employed as a genetic resource in a breeding programme or released for commercialization while also presenting a graphical display of the genetic correlations among traits [27,28]. Thus, the objectives of this study were (1) to assess the genetic diversity and population structure of selected heat-and drought-tolerant early maturing white and yellow maize inbreds using DArT SNP markers and (2) examine the interrelationships among traits of the early maturing inbreds in an effort to identify those strongly associated with grain yield in maize under combined HD stress.

Description of Germplasm
A total of 162 early maturing maize inbred lines developed in the IITA-MBP, Ibadan, Nigeria, with tolerance to drought and combined heat and drought stress were used in the present study ( Table 1). The IITA inbred lines were extracted from the broad-based populations (TZE Comp5-Y C 6 , TZE-W P OP STR 104 C 0 , TZE-W P OP STR 107 C 0 , TZE-Y P OP STR 106 C 0 , TZE-W P OP STR 108 C 0 , TZE-W Pop STR C 0, TZE-W Pop × LD C 0 and TZE-Y Pop STR C 0 ) developed from exotic and local germplasm identified as a result of many years of extensive testing for broad adaptation to dry and hot climates of the savannas of WCA. Furthermore, some of the inbred lines were derived from crosses between the broad-based drought-tolerant population, TZE COMP 3 DT C 2 F 2 , and elite Striga resistant inbred lines, TZEI 18, TZEI 4, TZEI 60 and TZEI 98, as well as crosses between the broad-based Striga resistant but drought-susceptible population, TZE COMP 5-W DT C 7 , and the drought tolerant inbred lines TZEI 2, TZEI 18, TZEI 31, TZEI 56 and TZEI 65. In addition to the 34 lines derived from two IITA bi-parental populations [(TZEI 11 × TZEI 8) and (TZEI 7 × TZEI 3)], the genetic materials for this study also included three inbred lines derived from two CIMMYT broad-based populations (M37W/ZM607 and Cuba/Guad C 3 ).

DNA Extraction and Genotyping
The 162 early maturing inbred lines were planted at the IITA research field, Ibadan, Nigeria under optimal growing conditions in 2018. At two weeks after planting, leaf samples were collected from each of the inbred lines, dried in a Labconco Freezone 2.5 L System lycophilizer (Marshall Scientific, Hampton, NH, USA) and ground using a Spex™ Sample Prep 2010 Geno/Grinder (Thomas Scientific, Swedesboro, NJ, USA).

Data Filtering and Analysis
Genomic DNA samples were extracted from freeze-dried leaf tissues of each inbred line using the DArT protocol (www.diversityarrys.com/files/DArTDNAisolation.pdf) as previously reported [18]. The quality of the DNA sample was determined using a 1% agarose gel test whilst the quantity was measured using the Nanodrop system (Thermo Scientific, Swedesboro, NJ, USA). Genotyping was performed at the IGSS-BecA-ILRI Platform, Nairobi, Kenya.
Following the methodology previously described [15,29], the DArT-SNP markers with missing data above 5% and 20% heterozygosity and minor allele frequency of below 0.05 were removed. From the initial DArTseq markers obtained from sequencing, a total of 9684 SNP markers distributed across the 10 maize chromosomes were retained for all subsequent analyses. The major allele frequency, polymorphic information content (PIC), number of alleles, heterozygosity and the gene diversity were calculated using Power Marker software V3.2.5 [30]. Analysis of the population structure was done by employing the admixture model using the STRUCTURE 2.3.4 software package [31]. The results from the model were achieved by changing the number of clusters (K) from one to twenty with ten to twenty adjustments for each K and subjecting the data to 10,000 Markov Chain Monte Carlo and 10,000 burn-in as described previously [15,25,29]. The most appropriate K value was identified by running the data obtained into the STRUCTURE HARVESTER using the Evanno method [32,33]. Each inbred was allocated to a group using a cut off threshold of 70%, while inbreds showing less than the designated threshold value were assigned to an additional group described as the mixed group.
Subsequent to verification of the number of sub-populations using the STRUCTURE, the phylogenetic relationship among the inbred lines was accessed using the neighbor-joining (NJ) method under 30,000 bootstraps in the DARwin software [34]. The genetic distance matrix was obtained using the Jaccard similarity test [35] using the formula dij = (a + b)/ (a + (a + c)), where dij is the dissimilarity between units i and j, a is the number of variables where Xi and Xj are present, b is the number of variables where Xi is present and Xj is absent, c is the number of variables where Xi is absent and Xj is present.

Assessment of the Performance of the Inbreds Under Combined Heat and Drought Stress Conditions
The set of 162 inbred lines were planted at Manga (11 • (Table 2). Similarly, at Manga, although some amount of rainfall was recorded in the month of April during 2018 and 2019, the highest day and night temperatures recorded during the growth cycle of the plants were obtained in April ( Table 2). The maize plants were irrigated twice a week using a drip irrigation system that supplied 17 mm of water per irrigation from planting to 2-3 weeks before anthesis at Manga, while at Kadawa, a furrow irrigation system supplied the same amount of water from planting to 2-3 weeks before anthesis when the irrigation water was withdrawn at each test location. A follow-up rescue irrigation was employed 2 weeks after the withdrawal of irrigation water at both test locations and resumed 10-12 days after the "rescue irrigation" until physiological maturity. The rescue irrigation was applied based on the assessment of some physiological stress indicator traits such as severe leaf rolling during the early mornings, leaf senescence, tassel blasting and leaf wilting. Basal fertilizer rate of 60 kg each of N, P 2 O 5 and K 2 O ha −1 of N, P, K (15:15:15) was applied followed by top dressing 2 weeks after planting (WAP) using 60 kg ha −1 of urea. Weed management was carried out using the post emergence herbicide (Roundup) and the pre-emergence herbicides (Caliherb 720 SL and Atrazine) in Ghana, while atrazine and gramozone were applied in Nigeria at 5 l ha -1 followed by manual weeding to ensure weed-free trials.

Data Collection
Agronomic data recorded under combined heat and drought stress at both Kadawa and Manga testing sites were similar to that described by Nelimor et al. [16]. We measured days to anthesis (AD) as the number of days from planting to when 50% of the plants in a plot had extruded pollen, and days to silking (SD) as when 50% of the plants had produced silks. The anthesis-silking interval (ASI) was recorded as the difference between SD and AD. Plant and ear heights were measured at about three weeks after flowering. Leaf death (LD) was determined by estimating the percentage of leaves dead from the base of the plant to the flag leaf using a scale of 1-9 where 1 = 1% to 10% of leaves dead and 9 = 81% to 100% of leaves dead at 70 days after planting (DAP). At about three weeks post-flowering, data were collected on husk cover (HUSKC), plant aspect (PLASP), plant height (PLTH), ear height (EHT), leaf firing and tassel blasting (TB), as described by Badu-Apraku and Fakorede [22]. A few days prior to harvesting, data were recorded on root and stalk lodging and number of plants per plot. At harvesting, records were taken on ear aspect and number of ears per plant (EPP). EPP was obtained by dividing the total number of ears per plot by the number of plants per plot. For the experiments conducted under optimal environments, the moisture content of the grains at harvest was recorded for the shelled grains of five randomly selected ears per plot using a moisture meter. Grain yield (kg ha −1 ) was estimated from field weight of ears per plot, assuming a shelling percentage of 80, adjusted to moisture content of 15%. For the trials conducted under combined HD stress, grain weight on plot basis was obtained by shelling all harvested ears of each plot.

Statistical Analysis
The data were subjected to analysis of variance employing PROCGLM of SAS as well as the RANDOM statement with a TEST option [36] for grain yield and other measured traits. In the statistical model, environments, replicates within environments and incomplete blocks within replicates × environment interaction were considered as random factors, while the inbred lines were regarded as a fixed factor. Block effects on genotype means were adjusted using the lattice design proposed by Cochran and Cox [37]. Mean values of measured characters were separated using the standard error of difference (SED). The multiple trait selection index was adopted for the identification of inbred lines combining tolerance to drought and heat stress environments. The base index was computed using Yield, EPP, ASI, PLASP, EASP and LD [38] as follows: Mean values of the six measured traits were standardized prior to calculation of the MI as the traits were measured in different units. Furthermore, estimated means of the different traits of 15 best and 10 worst inbred lines were adopted for the Genotype × Trait (GT) biplot analysis [39,40] to obtain information on the secondary traits which could be relied on for indirect selection for combined HD stress. The polygon and the vector views of the GT biplot were constructed using the measured characters and the genotype-focused singular value partition ('SVP = 2 ), rendering the polygon view suitable for visualizing the relationship between traits.

Summary Statistics
The summary statistics from the 9684 DArT SNP markers is presented in Table 3. Heterozygosity values were in the range of 0.00 to 0.20 with a mean value of 0.08. The 9684 DArT SNP markers had a major allele frequency ranging from 0.50 to 0.95 with an average value of 0.79. Gene diversity (GD) values varied from 0.09 to 0.50 with a mean of 0.30. For this set of inbred lines, the polymorphic information content (PIC) varied from 0.08 to 0.38 with a mean value of 0.25 (Table 3). The minor allele frequencies were (MAF) 0.05 for minimum and 0.50 for maximum, with an average of 0.21.

Genetic Distance and Population Structure Using DArT SNP Markers
From the pairwise comparisons of the inbred lines, the genetic distance for the lines was in the range of 0.02 to 0.48, with a mean value of 0.42. About 70.32% of the genetic distance among the inbred lines fell within 0.416 and 0.461 ( Figure 3). Lines TZdEI 134 and TZdEI 198 recorded the least genetic distance (Supplementary Table S1), consistent with their genetic background as they shared similar ancestry and history of selection (Supplementary Table S2). Inbred lines TZEI 500 and TZEI 7 recorded the highest (0.48) genetic distance, indicating that the two inbred lines diverged significantly. TZEI 500 was developed from the bi-parental cross involving TZEI 11 × TZEI 8 through repeated selfing up to the 8th generation of inbreeding, while TZEI 7 was developed from repeated self-pollination of WEC STR up to the 8th generation . Moreover, TZEI 500 is characterized by yellow endosperm while TZEI 7 has white endosperm (Supplementary Table S2).  Table S2). About 53.09% (86 of the inbreds) with the probability of association of less than 70% were grouped as mixed individuals (Supplementary Table S2). These six sub-populations were separated based largely on the pedigrees and endosperm colours. Sub-populations 1, 2, 4, 5 and 6 consisted of white endosperm maize inbred lines while sub-population 3 comprised yellow endosperm lines, with about 53% consisting of a mixture of white and yellow endosperm colour (Supplementary Table S2). The results showed that subpopulations 2, 3 and 4 contained only an inbred line, each derived from the different source populations: TZEI 11 × TZEI 8 S 8 , TZE COMP 5-W DT C7 × TZEI 65 S 8 and TZE COMP 3 DT C 2 F 2 × TZEI 98 S 8 , respectively. Sub-population 5 consisted of two source populations with inbred lines derived from TZE W-P OP STR 108 S 8 and TZE -W P OP STR 104 S 8 . Sub-population 1 consisted of four dissimilar source populations, which included TZE COMP 5 W DT C7 × TZEI 2, (TZEI 7 × TZEI 3) S 8 , TZE -Y P OP STR 106 S 5 , and TZE COMP 5 -W DT C7 × TZEI 31 S 8 . Sub-population 6 had the highest genetic diversity and comprised five source populations, namely WEC STR S 8 , TZE-W P OP Co S 8 , TZE -W P OP x LD S 8 , TZE COMP 3 DT C 2 F 2 × TZEI 60 S 8 , TZE-Y P OP STR 106 S 8 and TZE W P OP STR 107 S 8 . Within the six sub-populations, the expected heterozygosity amongst the inbred lines ranged from 0.10 for sub-population 4 to 0.43 for sub-population 6 with a mean of 0.25 (Supplementary Table S3a) Table S3c).

Cluster Analysis and Principal Coordinate Analysis (PCoA)
The neighbour joining (NJ) phylogenetic tree constructed assigned the 162 inbred lines into three main clusters followed by several sub-clusters ( Figure 5). In conformity with the biplot output of the STRUCTURE analysis ( Figure 4A,B), the six clusters could be clearly depicted in the phylogenetic tree ( Figure 5). Unique colours were used to represent each sub-population with red, green, blue, yellow, pink and teal denoting 1-6, respectively, and the black colour represented the mixed population when 70% probability threshold was considered.  Figure 5). The second cluster included 29 inbred lines, of which 27 clearly represented sub-population 2 (green bars). All the inbred lines in this cluster were derived from the yellow endosperm bi-parental cross (TZEI 11 × TZEI 8) S 8 (Supplementary Table S2). The third cluster consisted of 16 inbred lines with 14 being similar to those generated from the population structure analysis (sub-population 3; blue bars). All the inbred lines from this cluster were derived from one white kernel source population (TZE COMP 5-W DT C 7 X TZEI 65) S 8 . Cluster 4 had 23 inbred lines out of which 5 were derived from the white endosperm population, (TZE COMP 3 DT C 2 F 2 × TZEI 98) S 8 , and were similar to those generated from the population structure (sub-population 4; orange bars). The fifth cluster obtained from the phylogenetic analysis comprised a total of 18 inbred lines, with 10 out of these belonging to the white endosperm maize population, TZE-W Pop STR 104 S 5 (sub-population 5; pink colour). The results of the grouping of the 10 inbreds lines into sub-population 5 was consistent with the results obtained from the population structure analysis. The remaining eight inbred lines of the fifth cluster were derived from the source population, TZE-W Pop STR 108, and were also characterized by white endosperm kernels. The sixth cluster consisted of 12 inbred lines, with seven of them similar to those generated from the population structure analysis (indicated with teal colour lines) and represented by sub-population 6. Inbred lines placed in this cluster were derived from six different source populations with each inbred line extracted from a different source population, with the exception of the inbred lines TZdEI 47 and TZdEI 52 derived from the same source population, TZE-Y Pop STR 106 (Supplementary Table S2).
To understand and confirm the dynamics of the inbred population structure, we also adopted the principal coordinate analysis (PCoA). Based on the pairwise genetic distance matrix among the 162 inbred lines, a clear distinction among the six groups of inbred lines could be visualized, and this was in concordance with the results of the population STRUCTURE biplot ( Figure 6).

Genotype-by-Trait (GT) Interaction and the Multiple Trait Selection Index
The GT biplot analysis is one of the most reliable statistical tools available for selecting inbred lines for hybrid development using secondary traits. This analysis can identify secondary traits which are directly or indirectly related to traits of major breeding interest such as grain yield and ears per plant. The inbred lines were evaluated under combined heat and drought stress in this study to select desirable seed parents that could be used in productive hybrid combinations. Moreover, the evaluation aimed to examine the inter-relationships among traits measured for the inbreds as a way of predicting outstanding hybrid performance. The present study identified 15 inbred lines with outstanding performance as well as those with poorest performance (worst 10 lines) under combined heat and drought environments using the GT biplot (Figure 7). The inbred lines were selected based on the premise that the best selected lines should not record 50% or more percentage count for leaf firing and tassel blasting under combined HD stress. Approximately 52.56% of the total variation among the measured traits of the inbred lines was accounted for by the principal component axes PCI and PC2. In the biplot view, high values of grain yield and EPP for inbreds were considered as desirable while high values for the secondary traits such as LF, TB, EASP, PLASP, LD, ASI, HUSKC, AD, SD, and SL of inbreds were considered undesirable. The inbred(s) at each vertex of the polygon (vertex inbred) had the highest values for the traits within its sector. The results revealed that the sector with inbreds 21 and 8 as the vertex inbreds contained the traits LD, TB, LF, ASI, HUSKC, SL, EASP and PLASP, indicating that these inbreds had high values for these traits and that they were the low-yielding inbreds based on the high values of these traits. Inbred lines 25 and 20 were positioned at the vertex that contained the RL traits, signifying that these inbreds recorded high values for root lodging. In contrast, inbred 1 (the best performing inbred line) was the vertex inbred in the segment containing EPP and grain yield, suggesting the superiority of the inbred line in terms of grain yield and EPP. Inbred line 7 was the vertex inbred that recorded no trait within its sector, implying that it was not outstanding in any of the measured traits. Fifteen inbred lines were identified as outstanding in performance based on the measured traits. TZEI 935 was identified as the best performing inbred line with a yield of 3461 kg/ha. Contrarily, 10 inbred lines were identified as the worst in performance based on the measured traits, with TZdEI 109 as the worst performing inbred line with a yield of 424 kg/ha (Table 4).

Discussion
Molecular markers have proven to be among the most reliable tools for identifying the genetic components of inbred lines [15,29]. In the present study, 9684 DArT-based SNP markers were used to assess the genetic diversity of 162 early maturing yellow and white endosperm tropical maize inbreds possessing varying degrees of tolerance to drought and/or combined HD stress. Moreover, the population structure, heterotic groups and inter-trait relationships were also examined in the studied materials.
Generally, high PIC values indicate the effectiveness of markers for linkage analysis when estimating the inheritance between parental lines and the derived hybrids [41], whereas the GD or expected heterozygosity (He) reveals gene diversity for haploid markers and gives information on the proportion of the genetic distance and average heterozygosity within a given population [41,42]. In the present study, an average PIC of 0.25 was obtained, which is comparable to that reported earlier by Dao et al. [27] and Zhang et al. [28], as well as the 0.24 of Yu et al. [43] and 0.29 of Wu et al. [44], but higher compared with the results of Adu-Boakyewaa et al. [15] and Silva et al. [45] who reported average PIC values of 0.19 and 0.17, respectively. The high mean PIC value observed implied that majority of the selected SNPs were informative and polymorphic enough to bring out the differences among the 162 inbred lines studied. The average GD of 0.30 obtained in this study is also comparable to 0.32 previously reported in maize [46,47], but higher than the value (0.22) reported by other researchers [15,44,48,49]. These results indicated a broad genetic variability among the inbred lines and the relatively lower number of pairwise individuals with low genetic distance.
Of the 9684 SNPs deployed across the 162 lines, approximately 65% had MAF of between 0.05 and 0.25 (Supplementary Table S2), which is lower compared to the range reported by Semagn et al. [17] involving 450 CIMMYT lines studied using 1065 SNPs. Nevertheless, it is comparable to the values obtained by Adewale et al. [19] in early maturing white tropical maize. Consistent with the results of Adu-Boakyewaa et al. [15] and Adewale et al. [29], the average heterozygosity of 0.08 obtained in the present study is close to zero as expected for inbred lines, thus making them a useful genetic resource for association mapping and genetic studies in maize in which uniformity of the genetic materials is invaluable.
A higher F ST value for the sub-populations was reported in the present study. Together with the relatively moderate allele frequency divergence among sub-populations, the results indicated that the alleles of the inbred lines were fixed and could be categorized into groups with distinct characteristics, and this finding is in conformity with previous results [15]. The six clusters observed in this study indicated that individuals from the different groups were therefore expected to harbour different favourable alleles that could be exploited for breeding for combined HD-tolerant hybrids and synthetics. The results also demonstrated that the DArT-based SNP-derived markers were informative in providing genome profiles which were very useful for the identification of unique characteristics among the inbred lines [17,50,51].
Population structure analysis provides a guide for assigning a set of inbred lines into heterotic groups based on genotypic information of individual ancestry of the lines [52]. In this study, six sub-populations (K = 6) were established among the 162 inbred lines. The inbred lines were placed in distinct groups based on their relatedness in terms of pedigree records and selection history, such as endosperm colour, tolerance or susceptibility to drought and combined HD. This finding was affirmed by the high F ST values ranging from 0.28 to 0.81 for the six sub-populations and that inbred lines from the different groups could harbour unique favourable alleles for development of productive hybrids for population improvement programmes in the tropics. The high F ST values obtained for the sub-populations in the present study were comparable to the results of other researchers. For example, Kashiani et al. [53], Aci et al. [54] and Adu-Boakyewaa et al. [15] indicated F ST values of as high as 0.94, 0.33 and 0.83, respectively, for the sub-populations obtained in their genetic diversity studies for different sets of tropical maize lines. In the present study, cluster 1 contained all the inbred lines with white endosperm and were generally tolerant to drought, cluster 2 had inbred lines with yellow endosperm as well as tolerance to drought, while cluster 3 had highly drought-tolerant white endosperm lines. However, all the inbred lines within these three sub-groups showed moderate tolerance to combined HD stress. This observation is not surprising because the inbred lines have long been selected under drought conditions compared to the results of the more recent evaluations under combined HD conditions.
The multivariate clustering methods showed high consistency regarding the number of groups and individuals assigned to each group, indicating the existence of six genetically distinct groups. This result corroborated the findings of Adu-Boakyewaa et al. [15] and Obeng-Bio et al. [51] who reported a high consistency among the model-based structure analysis, PCoA and the Neighbour-Joining phylogeny using the Roger's and Euclidian genetic distance methods, respectively. However, this result disagreed with the reports of Dao et al. [27] and Semagn et al. [17] who found a high consistency between the model-based structure and the PCoA but very low concordance with the Neighbour-Joining phylogeny generated using the Roger's genetic distance method. The differences in the results could be due to the variation in the inbred lines used in the different studies, the method used in deriving the genetic distances between the inbreds, as well as the different clustering algorithms used. In concordance with earlier researchers, it was evident that the clustering of tropical maize populations is largely consistent with the pedigree information [15,21,27,48,55].
Generally, it is advocated that source genotypes with favourable alleles for combined heat and drought tolerance should possess good agronomic characteristics necessary for combined HD-tolerant hybrid development without placing too much emphasis on the per se grain yield. The GT analysis identified ASI, EASP, PLASP, LF, TB and EPP as important secondary traits that could be utilized to accelerate genetic gains from selection for enhanced grain yield under combined HD conditions. Nelimor et al. [16] suggested that lines that scored 50% or more for LF and TB under combined HD conditions should be rejected, implying the importance of the information on inter-relationships among secondary traits. From the GT biplot view, the best 15 performing lines were placed in five different heterotic groups. As a guide, crosses among the lines from different heterotic groups should be given priority in terms of hybrid combinations towards maximization of heterosis. For example, the best performing inbred line (TZEI 935) placed in the mixed group could be crossed to inbred lines TdZEI 311 and TZEI 1228 which belonged to groups six and three, respectively, to increase the chances of developing productive hybrids. The significant and positive association observed between PLASP and EASP implied that there was a close relationship among these traits and that either of them could be eliminated from the base index without sacrificing the precision in selecting for combined HD-tolerant genotypes. Badu-Apraku et al. [38] also found a positive correlation between these traits and advocated that either PLASP or EASP could be discarded in the base index in selecting for drought-and low-N-tolerant genotypes. Similarly, there was high positive correlation between AD, SD, and ASI, but none of them were found to be redundant as the ASI is dependent on the SD and AD [38]. These three traits were very important in selecting stress-tolerant maize genotypes [56,57], as a relatively constant value (lower ASI) is desirable in selection of early maturing maize [38].

Conclusions
To ensure the selection of appropriate parental lines for crossing and devising strategies for significant gains from selection, adequate knowledge and understanding of the genetic diversity of available inbred lines in a breeding programme is essential. The population structure analysis, PCoA and phylogenetic clustering methods confirmed the existence of six different groups for the panel of 162 maize inbred lines studied. The lines were grouped largely based on the pedigree information, endosperm colour, selection history and similarity of ancestry. There were high genetic distances between the paired inbred lines, demonstrating the distinctiveness of the inbred lines and the presence of high genetic variability which could be exploited in tropical maize breeding programmes for hybrid development. The low heterozygosity as well as relatively low divergence among the sub-populations indicated that the present set of inbred lines had great potential for genetic studies and could serve as future sources of promising inbred lines by contributing new and favourable alleles. The GT biplot analysis revealed that EPP, PLASP, EASP, LF and TB were the most reliable traits for selecting genotypes for tolerance to combined HD stress. Inbred 1 (TZEI 935) was identified as the ideal under combined HD stress environments. Using the base index, the 15 top-yielding inbred lines were identified for use in the maize breeding programmes of SSA for the development of productive hybrids and synthetics with combined HD tolerance. Among the inbred lines studied, inbred TZEI 935 was identified as one of the ideal lines under combined HD stress environments and should be further explored. Finally, the results of the present study will serve as an important guide for appropriate decision-making for future development of productive hybrids, particularly for SSA.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/9/1324/s1, Table S1. Genetic distance matrix of the 162 yellow and white inbred lines using Euclidean method in the present study. Table S2. Summary of the 9684 SNP markers used in the current study. Table S3