Phenotypic variation and heritability in kernel lysine content
Four DH populations were developed bu using eight inbred lines, which had the lysine content with a range of 0.18–0.44%. Totally, the four DH populations included 190–316 lines, respectively (Table 1). The lysine content exhibited continuously and approximately normal distribution in each DH population with a range of 0.12–0.50% (Fig. 1, Table 1). ANOVA variance analysis exhibited that genotype variance was greater than environmental variance in nearly all populations. It was demonstrated that lysine content variations were mainly controlled by genetic factors and the alleles responsible for increasing the phenotype reside in both parents. High broad-sense heritability estimates were calculated for lysine content in the four DH populations, with a range of 93–96% (Table 1), which indicated that the majority of lysine content variations are controlled by genetic factors and suitable for further QTL mapping.
Table 1
Phenotypic performance, variance, and broad-sense heritability of lysine content in the four DH populations
Trait a | Populations |
AF109 | AF116 | AF129 | AF170 |
Parents | | | | | |
means ± SD (%) | KB319003 | 0.21 ± 0.03 | KB519007 | 0.34 ± 0.04 | KB320005 | 0.25 ± 0.04 | JinQingWL2 | 0.18 ± 0.01 |
AJ317001 | 0.38 ± 0.02 | KB717001 | 0.43 ± 0.03 | KB120001 | 0.37 ± 0.03 | KB319004 | 0.44 ± 0.02 |
p value b | <0.0001**** | <0.0002*** | <0.0001**** | <0.0001**** |
DHs | | | | | | | | |
Size | 248 | 190 | 316 | 265 |
means ± SD (%) | 0.32 ± 0.07 | 0.32 ± 0.07 | 0.33 ± 0.06 | 0.27 ± 0.08 |
Range (%) | 0.17–0.49 | 0.18–0.49 | 0.18–0.49 | 0.12–0.50 |
σg2 c | 0.018 | 0.015 | 0.013 | 0.016 |
σe2 d | 0.012 | 0.014 | 0.025 | 0.043 |
σε2 e | 0.082 | 0.073 | 0.134 | 0.068 |
h2 (%) f | 95.00% | 93.00% | 93.00% | 96.00% |
a lysine content; b P value based on a t-test evaluating two parental lines; c genetic variance; d environmental variance; e residual variance ; f broad-sense heritability (h2); *** p < 0.001, **** p < 0.0001.
Genotyping and genetic linkage map
All DH lines in four populations were genotyped using GenoBaits Maize 2K marker panel containing 10,378 SNP markers, and further refined by eliminating SNPs with MAF < 0.1 or missing rate > 0.6. This resulted in a total of 8,377 SNPs, with their precise physical positions based on the B73 reference sequence Version 4, were polymorphic between their respective parents for AF109, AF116, AF129 and AF170, respectively. In each DH population, the missing rate in most lines were less than 2%. The average genetic distance between every two adjacent markers was 0.98, 0.72, 0.66, and 0.78 cM in each DH population, respectively.
Identification of QTLs for lysine content in four DH populations
The QTL mapping for lysine content in the four DH populations and the related genetic features were summarized in Table 2. In total, 12 QTLs were identified with a LOD threshold of above 3.00 in the four populations (Fig. 2). These QTLs were located on chromosomes 1, 2, 3, 4, 5, 6, 9 and 10. The average of QTL physical intervals was 15.91 Mb in a range of 0.85–32.17 Mb. The average of the total PVE explained by all identified QTLs in a population was 17.68 and ranged from 8.70 (AF116) to 26.06% (AF129). This was less than broad-sense heritability (Table 1 and Table 2), suggesting that only part of QTLs have been detected in bi-parent populations.
In AF109, three QTLs (qLYS-1-1, qLYS-1-2 and qSC-1-3) were detected on chromosome 3, 4 and 9, respectively. Phenotypic variation explained by these QTLs was 4.42–4.98% and explained by the additive effect of all QTLs was 11.94%. The parental AJ317001 allele at qLYS-1-2 on chromosome 4 and qLYS-1-3 on chromosome 9 increased lysine content. The QTL qLYS-1-1 was located on chromosome 3 and explained 4.42% of phenotypic variance. The parental KB319003 allele at this locus decreased lysine content.
In AF116, a total of two QTLs (qLYS-2-1 and qLYS-2-2) were identified and accounted for 8.70% of the total phenotypic variance. The QTL qLYS-2-1 was located on chromosome 3 and contributed to 4.46% of the explained phenotypic variance. The QTL qLYS-2-2 on chromosome 9 explained 5.60% of phenotypic variance. Parental KB717001 allele at qLYS-2-2 increased lysine content, while parental KB519007 allele at qLYS-2-1 decreased lysine content.
In AF129, a total of three QTLs (qLYS-3-1, qLYS-3-2 and qLYS-3-3) were detected on chromosome 1, 2 and 6, respectively and explained by the additive effect of all QTLs was 26.06%. Phenotypic variation explained by each individual QTL was 5.29–12.66%. Highly effective QTL qLYS-3-3 were detected on chromosome 6 which explained a phenotypic variance of more than 10%, suggesting that qLYS-3-3 was the major QTL controlling lysine content in AF129. The parental KB320005 allele at qLYS-3-1 on chromosome 1 with a phenotypic variance of 5.57% and QTL qLYS-3-3 on chromosome 6 increased lysine content. The QTL qLYS-3-2 was located on chromosome 2 and explained 5.29% of phenotypic variance. The parental KB120001 allele at this locus decreased lysine content.
In AF170, four QTLs (qLYS-4-1, qLYS-4-2, qLYS-4-3 and qLYS-4-4) distributed on chromosome 1, 5, 9 and 10, respectively. Phenotypic variation explained by these QTLs was 5.23–8.65% and explained by the additive effect of all QTLs was 24.03%. The parental JinQingWL2 allele at qLYS-4-1, qLYS-4-2 and qLYS-4-3 increased the lysine content with the phenotypic variance of 6.92–8.65%, whereas parental KB319004 at qLYS-4-4 decreased lysine content with a phenotypic variance of 5.23%.
Table 2
Individual QTL for lysine content in the four DH populations
Populations | QTL | Chr.a | P-Peak (Mb)_V4b | P-Range (Mb)_V4c | LOD | PVE%d | Add.e | Parentf+ | PVE(%) -ALLg |
AF109 | qLYS-1-1 | 3 | 69.67 | 59.07–97.50 | 3.994 | 4.42 | -0.041 | KB319003 | 11.94 |
qLYS-1-2 | 4 | 44.01 | 37.64–65.21 | 3.867 | 4.67 | 0.017 | AJ317001 |
qLYS-1-3 | 9 | 130.80 | 129.70-137.80 | 4.034 | 4.98 | 0.021 | AJ317001 |
AF116 | qLYS-2-1 | 3 | 161.43 | 161.43–169.00 | 3.892 | 4.46 | -0.022 | KB519007 | 8.70 |
qLYS-2-2 | 9 | 132.53 | 130.80-149.39 | 4.661 | 5.60 | 0.031 | KB717001 |
AF129 | qLYS-3-1 | 1 | 267.92 | 267.37-271.96 | 4.226 | 5.57 | 0.083 | KB320005 | 26.06 |
qLYS-3-2 | 2 | 214.01 | 210.90-219.60 | 4.552 | 5.29 | -0.014 | KB120001 |
qLYS-3-3 | 6 | 131.71 | 117.49-141.72 | 6.731 | 12.66 | 0.016 | KB320005 |
AF170 | qLYS-4-1 | 1 | 150.61 | 143.21-159.07 | 3.529 | 7.57 | 0.031 | JinQingWL2 | 24.03 |
qLYS-4-2 | 5 | 12.70 | 12.70-13.55 | 4.111 | 8.65 | 0.020 | JinQingWL2 |
qLYS-4-3 | 9 | 111.07 | 93.07-125.24 | 4.990 | 6.92 | 0.028 | JinQingWL2 |
qLYS-4-4 | 10 | 14.80 | 13.61–17.92 | 3.333 | 5.23 | -0.020 | KB319004 |
a Chromosome; b Physical position of QTL based on the B73 reference sequence (V4); c Physical position range of QTL based on the B73 reference sequence (V4); d Percentage of the phenotypic variation explained by the additive effect of QTL; e Additive effect of QTL; f which parental allele increased lysine content; g Percentage of the phenotypic variation explained by the additive effect of all QTL.
Genetic overlap of QTLs in the four DH populations
To evaluate genetic overlaps among different mapping populations, the 1.5-LOD support interval of QTLs in the four DH populations and other populations for lysine content previously reported were compared (Fig. 3). QTLs with overlapping support intervals were considered as common QTLs. Among the four DH population, a 7.00 Mb overlap was observed between qLYS-1-3 and qLYS-2-2 on chromosome 9. Furthermore, there were only a few of overlaps detected after comparing the results with all types of other populations reported. To the end, there were 8.70 Mb and 7.57 Mb overlap with the QTLs in F2:3 progeny derived from high free amino acids (FAA) parents Oh545o2 and Oh51Ao2 [15].