Original Research ARTICLE
Large-Scale Analysis of Combining Ability and Heterosis for Development of Hybrid Maize Breeding Strategies Using Diverse Germplasm Resources
- 1College of Agriculture, Northeast Agricultural University, Harbin, China
- 2Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing, China
- 3Qiqihar Branch of Heilongjiang Academy of Agricultural Sciences, Qiqihar, China
- 4CIMMYT-China Specialty Maize Research Center, Shanghai Academy of Agricultural Sciences, Shanghai, China
- 5CIMMYT-China Tropical Maize Research Center, Foshan University, Foshan, China
- 6International Maize and Wheat Improvement Center (CIMMYT), Texcoco, Mexico
Understanding combining ability and heterosis among diverse maize germplasm resources is important for breeding hybrid maize (Zea mays L.). Using 28 temperate and 23 tropical maize inbreds that represent different ecotypes and worldwide diversity of maize germplasm, we first developed a large-scale multiple-hybrid population (MHP) with 724 hybrids, which could be divided into three subsets, 325 temperate diallel hybrids and 136 tropical diallel hybrids generated in Griffing IV, and 263 temperate by tropical hybrids generated in NCD II. All the parental lines and hybrids were evaluated for 11 traits in replicated tests across two locations and three years. Several widely used inbreds showed strong general combining ability (GCA), and their derived hybrids showed strong specific combining ability (SCA). Heterosis is a quantifiable, trait-dependent and environment-specific phenotype, and the response of parental lines and their hybrids to environments resulted in various levels of heterosis. For all the tested traits except plant height and hundred grain weight (HGW), NCD II (temperate × tropical) hybrids showed higher average heterosis than the temperate and tropical diallel hybrids, with higher hybrid performance for ear length, ear diameter, and HGW. Tropical maize germplasm can be used to improve the yield potential for temperate lines. Grain number per row and grain number per ear were two most important traits that determined yield heterosis, which can be used as direct selection criteria for yield heterosis. The hybrids from heterotic groups, Reid × SPT, Reid × LRC, SPT × PA, and Lancaster × LRC, contributed highly significant positive SCA effects and strong heterosis to yield-related traits, and the heterotic patterns identified in this study were potentially useful for commercial maize breeding. Heterosis was more significantly and positively correlated with SCA than GCA, indicating that SCA can be used in heterosis prediction to develop potential hybrids in commercial maize breeding. The results of the present study not only contribute to developing breeding strategies, but also improve targeted breeding efficiency by using both temperate and tropical maize to broaden genetic basis. Large sets of parental lines with available genotypic information can be shared and used in worldwide hybrid breeding programs through an open-source breeding strategy. Potential applications of the reported results in developing hybrid maize breeding strategies were also discussed.
Maize (Zea mays L.) is one of the most important crops for staple food, livestock feed, edible oil, and biofuel (Mackay, 2009). Its cradle is in America’s tropical and subtropical areas, subsequently generating a variety of ecotypes and diverse germplasm through a process of evolution and domestication. Compared to temperate maize, those developed in tropical and subtropical zones usually have more diverse genetic variation with resistance to diseases and pests, flourishing roots, toughness stalk, lodging resistance, drought tolerance, and higher levels of heterosis in their hybrids with temperate inbreds (Vasal et al., 1992). The phenomenon of heterosis or hybrid vigor was perceived by Darwin (1876) and described as hybrid F1 offspring exhibiting phenotypic superiority than both parents (East, 1908; Shull, 1908). In plants, heterosis could be attributed to the interaction among multiple loci, depending on hybrids and traits (Schnable and Springer, 2013), as shown in the magnitude and ratio of heterosis for biomass (Li et al., 2001), flowering related traits (Krieger et al., 2010), yield (Luo et al., 2001), and resistance to abiotic and biotic stresses (Miller et al., 2015).
Utilization of heterosis is of great importance for agricultural production and one of the most successful examples in crops is from maize (Duvick, 2001). Breeding practice indicates that the performance of parents per se is not consistent with the hybrid performance. Excellent hybrid varieties are not necessarily derived from elite parents. Therefore, breeders should judge a parental line by its potential to produce superior hybrids, not only by its performance per se (Riedelsheimer et al., 2012). To identify the parental lines with great potential in making hybrids, combining ability has been estimated and used to select desirable parents and thus their hybrids. Two types of genetic parameters, GCA and SCA, have been used, which may be primarily caused by additive and non-additive gene actions, respectively (Sprague and Tatum, 1942). GCA for an inbred line is measured as the average performance for all the hybrids produced with that inbred line as the common parent, and SCA for a specific cross or hybrid is measured by the deviation of the hybrid performance from what can be predicted by the parental GCA (Sprague and Tatum, 1942). In generally, GCA evaluation is performed at early generations or breeding stages in order to save time and money in hybrid breeding, particularly in hybrid maize breeding (Sprague, 1946). The relative contribution of GCA and SCA effects to hybrid performance depends on traits and hybrids, and in some cases, for example, aflatoxin and grain yield (GY) (Meseka et al., 2018), GCA effects are more important than SCA effects. Evaluating GCA is inevitably a cumbersome and time-consuming task, becoming one of the major constraints in hybrid breeding programs.
To breed ideal hybrids with high GY (Yong et al., 2019), good quality (Machida et al., 2010) and strong resistance to biotic (Sibiya et al., 2013) and abiotic stresses (Makumbi et al., 2018), heterosis and combining ability have been analyzed for available germplasm with limited numbers of parental lines. Heterosis and combining ability are usually estimated by populations derived from special genetic designs, such as diallel (Griffing, 1956) and the NCD II (North Carolina design II) (Comstock and Robinson, 1948), which are two most powerful genetic designs for combining ability analysis and have been applied extensively. Using an NCD II with 6 × 18 parents, the combining ability analysis indicated that at least one parent with higher GCA is required for producing a hybrid cross with high SCA for nitrogen use efficiency (Cui et al., 2014). Heterotic grouping among 378 hybrids derived from diallel crosses of 28 early inbreds was evaluated for their tolerance to Striga hermonthica, indicating that grouping based on SCA and GCA was the most effective in classifying early maturing maize inbreds for tropical maize breeding programs (Akinwale et al., 2014). Therefore, the better understanding of the genetic basis of heterosis and combining ability we can get, the more effective maize improvement programs and hybrid performance prediction can be achieved (Dhillon and Singh, 1977). However, only limited numbers of parents and their hybrids have been used so far, with a few of exceptions that over 100 hybrids were used (Akinwale et al., 2014). The population for a hybrid crop with a large number of hybrids generated from mating designs can be simply called as multiple-hybrid population (MHP) (Wang et al., 2017). Using a specific mating design or a combination of multiple designs, large-scale MHPs can be produced for more effective analyses of combining ability and heterosis. Diverse germplasm from different ecotypes, including temperate and tropical maize, should have been used for identifying genetic variation for both basic research and commercial breeding. So far, very few analyses of combining ability and heterosis have been performed using the hybrids between different ecotypes in maize (Fan et al., 2016). Considering great genetic diversity existing in tropical maize germplasm that could contribute to further genetic improvement, more studies are required by using between-ecotype hybrids. Therefore, large-scale analysis of heterosis and combining ability using diverse germplasm resources will improve our understanding of hybrid performance significantly, thus contributing to increased genetic gain in maize hybrid breeding.
Here we report a large-scale analysis of combining ability and heterosis using an MHP consisting of 724 hybrids derived from 28 temperate and 23 tropical maize inbred lines. Our objective was to measure heterosis in both intra- and between-ecotype hybrids, estimate GCA and SCA effects and compare combining ability and heterosis across ecotypes and environments. Our results will facilitate our future maize breeding through improved combining ability analysis, defined heterotic patterns and broaden genetic basis using different maize ecotypes. The large-scale phenotypic data, combining with high-density genotypes, which can be shared and used in worldwide hybrid breeding programs through open-source breeding strategy, will provide a great opportunity for whole genome prediction of heterosis and hybrid performance.
Materials and Methods
A maize MHP was developed by using diallel and NCD II mating designs consisting of 28 temperate and 23 tropical inbred lines, representing a broad selection of breeding germplasm from temperate and tropical regions (Wang et al., 2017). The 724 hybrids were developed, which were divided into three subsets, 325 temperate hybrids derived in Griffing IV involving 26 parental lines, 136 tropical diallel hybrids involving 17 parental lines, and 263 NCD II hybrids generated between 13 temperate and 21 tropical parental lines. Among the 28 temperate inbred lines, 15 were from China and 13 were from United States, which included six heterotic groups, Reid, SPT, LRC, Lancaster, PA, and PB (Wang et al., 2017). Among the 15 Chinese temperate lines, six are common testers, including Ye478, HZ4, Dan340, Mo17, Tie7922, and Qi319, which have been widely used in Chinese maize breeding programs. The rest nine Chinese temperate inbreds have been also playing a very important role in hybrid breeding across maize zones in China. The 23 tropical lines have been widely used as parents across worldwide breeding programs in China and CIMMYT, three of which, Jiao51, Chuan29 Female, and 18-599, are from China.
Field Experiments and Data Collection
The temperate diallel hybrids, NCD II hybrids and their parental inbred lines were phenotyped in 3 years (2013–2015) and two locations, Shunyi, Beijing (116.6°E, 40.2°N) and Xinxiang, Henan (113.8°E, 35.1°N) in randomized block design with two replications. The tropical diallel hybrids were phenotyped in Jinghong, Yunnan (100.8°E, 22.0°N) in 2014 and Sanya, Hainan (109.2°E, 18.4°N) in 2015. In each replication, each entry was planted with two 4-m rows with 25 cm between plants and 60 cm between rows. Thinning was done at the fifth leaf stage to maintain a density of 66,600 plants/ha. Traditional agronomical practices for local maize production were adopted in each trial to manage the experimental plots.
Eleven traits of agronomic importance for hybrid performance were investigated. Two flowering related traits, days to silk (DTS) and days to anthesis (DTA) were scored. At harvesting stage, plant height (PH) and ear height (EH) were measured and seven yield-related traits, including ear length (EL), ear diameter (ED), row number (RN), grain number per row (GNPR), grain number per ear (GNPE), hundred grain weight (HGW) and grain weight per plant (GWPP), were measured after the harvested ears were air-dried. PH and EH were measured as the average height from ground to the top of the tassel and from ground to the node of the ear, respectively, each with five consecutive plants after excluding the edge ones. DTS and DTA were measured as the number of days from sowing to 50% silk and 50% anthesis, respectively. EL, ED, and RN were measured as the length from the ear bottom to tip, the diameter in the ear middle, and RN per ear, respectively. HGW was estimated with three samples of 100 kernels randomly selected from the total kernels and measured to give the average.
where gi and gj are the GCA effects for i-th and j-th lines, respectively; sij is the SCA effect for ij-th hybrid; yij is the trait value of ij-th hybrid; is the average of the hybrids among i-th line crossed with a series of parents; is the overall mean.
The genetic variances of GCA and SCA effects were obtained in a joint linear mixed model analysis of MHP over all tested environments by following Riedelsheimer et al. (2012):
where yijkl is the phenotypic observation for the i-th environment, μ is the overall mean, Li is the i-th fixed environment effect, bj(i) is the effect of j-th block within the i-th environment, gcak and gcal are the random GCA effects of the k-th female and the l-th male, scakl is the random SCA effects of the k-th and the l-th parents, L × gcaik and L × gcail are the random GCA by environment interaction effects, L × scaikl is the random SCA by location interaction effect, and eijkl is the random error.
where is the variance of GCA effects derived from the mean square of GCA and is the variance of SCA effects derived from the mean square of SCA. Since the total genetic variance among F1 hybrids is equal to twice the GCA component plus the SCA component, the closer this ratio is to unity, the greater the proportion of a specific hybrid’s performance can be predicted based on GCA alone (Baker, 1978).
(2) Heterosis was estimated based on two criteria, mid-parent heterosis (MPH) and high-parent heterosis (HPH), using the following formulae:
where F1 is the mean performance of F1 hybrids, MP is the parental mean, and HP is higher parent values for all tested traits.
(3) Analysis of variance (ANOVA) for phenotypic performance was performed based on combined data and linear mixed model by following Hallauer et al. (2010):
where yijkl is the i-th phenotypic observation for the j-th year, k-th environment, and l-th block, μ is the overall mean, vi is the effect of i-th cross, yj is the effect of j-th year, sk is the effect of k-th environment, (vy)ij is the interaction effect between i-th cross and j-th year, (vs)ik is the interaction effect between i-th cross and k-th environment, (ys)jk is the interaction effect between j-th year and k-th environment, (vys)ijk is the interaction effect among i-th cross, j-th year, and k-th environment, rl(jk) is the effect of l-th block within j-th year and k-th environment, and eijkl is the error. The genetic effect vi is considered as fixed effect while all other effects as random.
Broad sense heritability (H2) was calculated for each trait by following Knapp et al. (1985):
where is genotypic variance of the hybrids, is genotype × environment interaction variance, is error variance, L is the number of environments and R is the number of replications per location.
The ANOVA was performed to evaluate the effects of genotype (G), environment (E), and the interaction (G × E) using the PROC MIXED procedure of SAS®.1 Best linear unbiased predictions (BLUPs) were used to estimate phenotypic traits across multiple environments based on a linear model (Brown et al., 2011). Calculation of GCA, SCA, MPH and HPH were conducted based on the BLUP value for each trait. The correlation coefficients were assessed by the “cor.test” function in R, and the significance of the correlation coefficient was tested with t-test.
Significant Effects of Genotypes and Environments on Phenotypic Traits
Phenotypic performance across environments and years with tested traits for inbreds and hybrids is shown in Table 1. Significant genotype effects were found for all measured traits (P < 0.01), indicating significant genetic variation among parental lines and hybrids (Supplementary Table 1). For the temperate diallel and NCD II hybrids, environment and year had significant effects on all traits except DTS and DTA, and genotype by environment effect was significant for all traits except RN. Significant genotype by year interaction (P < 0.01) were also found for all traits but PH and HGW. Therefore, the phenotypic performance of inbreds and hybrids was significantly affected by genotype, environment and genotype by environment interaction (Supplementary Table 1).
Diverse Combining Ability Effects Contributed by Ecotypes and Heterotic Groups
Highly significant GCA and SCA effects were found for all tested traits (Table 2). Significant GCA × E and SCA × E interaction effects were revealed for all the traits but RN. GCA/(GCA + SCA) ratios indicate that the tested traits were predominantly controlled by additive gene effects. Broad-sense heritability estimated was high (0.77–0.93) for all traits except GWPP (0.59), suggesting that phenotypic variation observed in flowering time and yield-related traits was highly inheritable (Table 2). GCA effects were highly variable across the traits (Supplementary Table 2). In the temperate diallel, 13 temperate inbreds showed negative GCA effects for PH and EH, suggesting that these inbreds had genetic potential for reducing plant and EHs. Two inbreds, P11 (HZ1) and P22 (Qi319), showed positive GCA effects for GNPE, HGW, and GWPP. P26 (Zheng58) exhibited positive GCA effects for HGW and GWPP but negative effects for PH and EH, indicating its value in increasing GY and deceasing plant sizes. The line P2 (AS6103) showed negative effects on PH, EH and flowering related traits, with stress tolerance at early flowering stages. P19 (PH4CV) and P20 (PH6WC) had positive effects on yield traits. For 28 temperate inbred lines, American inbred lines had desirable GCA effects on HGW and GWPP, while Chinese ones had desirable GCA effects on other yield-related traits with shorter plant and EHs. For tropical parental lines, CIMMYT inbred lines had much better GCA and performance than Chinese lines. P38 (TR0415), P41 (622016-ZCN-2), P43 (CML312SR), and P48 (DTMA227B) contributed to shorter plant and EHs. P38 (TR0415), P43 (CML312SR), and P45 (CML330) contributed to early flowering. P39 (TR0423), P40 (18-599), P46 (CML504), and P50 (TR0582) contributed to higher yield-related traits. In the NCD II, temperate inbred lines had desirable GCA effects on EL, ED, HGW, DTT, and DTS, while tropical ones had more favorable effects on EH, RN, GNPR, GNPE, and GWPP.
Table 2. Analysis of variance for combining ability and heritability of tested traits in a maize multiple-hybrid population.
Several widely used inbreds showed higher desirable GCA effects, and SCA analysis demonstrates that the hybrids with elite materials as parents have been identified as combinations with higher desirable SCA effects (Supplementary Table 3). As shown in Figure 1 for the hybrids from different heterotic groups, the hybrids from LRC × PA showed the lowest SCA effects on PH, EH, GNPR, and GNPE but the highest SCA effects on DTS and DTA. PA × PB hybrids had the lowest SCA effects on yield-related traits, including ED, HGW, and GWPP. Reid × LRC hybrids had the highest SCA effects on PH, EH, and RN. Reid × SPT and Lan × PB hybrids had the lowest SCA effects on flowering related traits than the hybrids derived from other heterotic groups. There was no significant difference found in SCA effects for the hybrids between tropical and temperate groups. From the aforementioned analysis, Reid × SPT, Reid × PB, Reid × LRC, and SPT × PA had favorable SCA effects for yield-related traits, while Reid × SPT and Lan × PB had lower SCA effects for flowering related traits.
Figure 1. Specific combining ability (SCA) for hybrids between different heterotic groups in temperate diallel and NCD II designs. SCA is color-coded according to the color keys plotted on the bottom right. The blue and red boxes indicate positive and negative SCA effects, respectively. Different lowercase and uppercase letters indicate significant differences by Tukey test at 0.05 probability in temperate diallel and NCD II hybrids, respectively. Tropical: tropical inbred lines; Six maize heterotic groups are represented by Reid, Lan (Lancaster), SPT, LRC, PA, and PB. PH, plant height; EH, ear height; DTS, days to silk; DTA, days to anthesis; EL, ear length; ED, ear diameter; RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: temperate diallel (N = 325); NCD II (N = 263).
Prevalent Heterosis Observed for Difference Traits
Significant MPH and HPH were observed for all the tested traits (Tables 3, 4). Two flowering related traits, DTS and DTA, exhibited negative heterosis, while others exhibited positive heterosis. For all the tested traits except HGW, a higher level of heterosis was observed in NCD II hybrids than in temperate and tropical diallel hybrids, indicating a higher level of heterosis in between-ecotype hybrids. NCD II hybrids showed a higher level of negative heterosis for flowering related traits, indicating that compared to tropical by tropical hybrids, NCD II (temperate × tropical) hybrids, with shorten growth period and improved stress tolerance, have advantages allowing them to be planted in temperate conditions. Overall, NCD II hybrids showed high levels of heterosis for EH and yield-related traits, indicating that temperate by tropical hybrids showed obviously stronger heterosis than intra-ecotype hybrids. Temperate diallel hybrids showed the lowest level of heterosis on flowering related traits, illustrating that they flowered earlier than other hybrids. The result revealed that tropical hybrids showed stronger heterosis on vegetative traits, while between-ecotype hybrids, i.e., temperate by tropical hybrids, showed a higher level of yield-related heterosis.
Highly Correlated Heterosis for Related Traits
The correlation of heterosis among traits varied greatly (Figure 2). MPH showed highly significant correlation (0.87) between two flowering traits, DTS and DTA. Conversely, the correlations of the flowering traits with other traits were lower than 0.30. Yield-related traits exhibited highly significant MPH correlation with each other except HGW. A very similar trend was observed for HPH (Figure 2). Highly significant correlation was observed between MPH and HPH for all the tested traits, and the correlation for yield-related traits was higher than 0.77 (Table 5). Thus, observed heterosis largely depends on the genotype and could be correlated with related traits. Heterosis per se (MPH and HPH) can be treated as a phenotypic trait and used for genome selection and heterosis prediction.
Figure 2. The correlation of heterosis among 11 tested traits in a maize multiple-hybrid population. Correlation coefficients for mid-parent heterosis (MPH) are listed in the upper-right, while those for high-parent heterosis (HPH) are listed in the lower-left. The correlation level is color-coded according to the color keys plotted on the right. The blue and red boxes indicate positive and negative correlation coefficients, respectively. The symbols * and ** indicate significance at 0.05 and 0.01 probability levels, respectively. PH, plant height; EH, ear height; DTS, days to silk; DTA, days to anthesis; EL, ear length; ED, ear diameter; RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: multiple-hybrid population (N = 724).
Heterosis Varied Significantly Across Environments
In temperate diallel hybrids and NCD II hybrids, significantly different levels of MPH were observed between two environments, Shunyi and Xinxiang, for all tested traits except ED, GNPR, and HGW (Figure 3A). Significantly different levels of HPH were also observed for all tested traits except GNPR and HGW (Figure 3B). Parents, hybrids and their combined effects contributed to the varied heterosis levels across environments. The response of maize parents and their hybrids to environmental factors may result in different levels of heterosis. We used the coefficient of variation (CV) to evaluate the stability of parents and their hybrids across the two environments. For 8 of the 11 traits, parents exhibited significantly higher CVs than hybrids (Supplementary Figure 1), indicating that for these traits the unstable heterosis across environments was driven more significantly by unstable parental lines. The rest three traits, including GNPR, HGW, and GWPP, did not exhibit significantly different CVs between hybrids and parents. These results indicate that a great variability existed in hybrids and parental lines across environments for all tested traits, and hybrids were more stable than parental lines. From the above results, environmental variables affect heterosis and hybrid performance greatly when maize hybrids and inbred lines respond differently to environmental stimuli.
Figure 3. Heterosis performance in two environments, Shunyi and Xinxiang. (A) Mid-parent heterosis (MPH); (B) High-parent heterosis (HPH). In each figure red and green colors indicate Shunyi and Xinxiang, respectively. The symbols * and ** indicate significance at 0.05 and 0.01 probability levels, respectively. PH, plant height; EH, ear height; DTS, days to silk; DTA, days to anthesis; EL, ear length; ED, ear diameter; RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: temperate diallel (N = 325); NCD II (N = 263).
Heterotic Groups Contributed to Trait-Specific Advantages in Heterosis
We evaluated the average levels of heterosis across heterotic groups for the tested traits (Figures 4, 5). The Reid × LRC hybrids showed higher MPH for PH, but higher HPH for RN. The Lancaster × LRC hybrids had higher MPH for EH, and GNPE. The SPT × LRC hybrids had higher MPH for RN. Therefore, LRC group can be used to increase heterosis for yield-related traits. Both Reid × PA and Reid × PB hybrids had lower MPH and HPH for yield-related traits. SPT × LRC and Reid × SPT hybrids had lower MPH for DTS and DTA. Nonsignificant difference was found in MPH for DTA, EL, HGW, GNPR, and GWPP, or in HPH for EL, ED, HGW, DTA, GNPE, GNPR, and GWPP, among heterotic groups. Tropical × Lancaster, tropical × SPT, tropical × LRC hybrids showed higher MPH for GNPR, GNPE, and GWPP, but higher HPH for HGW, GNPE, and GWPP. Therefore, the hybrids from specific heterotic groups tended to exhibit specific trait advantages. For example, LRC group showed favorable heterosis for yield-related traits, while LRC, SPT and Lancaster groups had lower heterosis for flowering related traits.
Figure 4. Heterosis performance in different heterotic groups for mid-parent heterosis (MPH) in temperate diallel and NCD II designs. Different lowercase and uppercase letters indicate significant differences by Tukey test at 0.05 probability in temperate diallel and NCD II hybrids, respectively. Tro, tropical inbred lines; Six maize heterotic groups are represented by Reid, Lan (Lancaster), SPT, LRC, PA, and PB. PH, plant height; EH, ear height; DTS, days to silk; DTA, days to anthesis; EL, ear length; ED, ear diameter; RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: temperate diallel (N = 325); NCD II (N = 263).
Figure 5. Heterosis performance in different heterotic groups for high-parent heterosis (HPH) in temperate diallel and NCD II designs. Different lowercase and uppercase letters indicate significant differences by Tukey test at 0.05 probability in temperate diallel and NCD II hybrids, respectively. Tro, tropical inbred lines; Six maize heterotic groups are represented by Reid, Lan (Lancaster), SPT, LRC, PA, and PB. PH, plant height; EH, ear height; DTS, days to silk; DTA, days to anthesis; EL, ear length; ED, ear diameter; RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: temperate diallel (N = 325); NCD II (N = 263).
GNPR and GNPE as Two Important Traits for Yield Heterosis
Yield was evaluated based on GWPP. The tested hybrids showed significant MPH and HPH for both yield and six yield-related traits (Tables 3, 4). GWPP in temperate diallel hybrids showed 141.30% (MPH) and 116.21% (HPH) on average over parental lines. Moreover, the hybrids showed significantly higher yield-component traits than both parents. The means and ranges of seven yield-related traits in three subset populations were shown in Supplementary Figure 4. The result indicates that apparent heterosis for yield-related traits was observed in the three types of populations. The average performance of hybrids in the temperate diallel was higher than that in the NCD II and tropical diallel hybrids for ED, RN, GNPR, GNPE, while the performance in NCD II was higher than that in other two populations for EL, ED, and HGW. Temperate × tropical hybrids exhibited better yield performance and heterosis than within-ecotype hybrids for most of yield-related traits. Maize yield is determined by GNPR, GNPE, RN, and HGW. Improvement of GNPE related traits is an effective way for breeding yield heterosis. As GNPE is the multiplication of GNPR and RN, in-depth analysis of the traits related to GNPE has important theoretical and practical impacts on hybrid breeding. However, GNPR and RN are two traits that interact with each other, and GNPE is determined by both GNPR and RN. MPH showed highly significant correlation between GNPR and GNPE (0.95), between GNPE and RN (0.75) and between RN and GNPR (0.55) (Figure 2). GNPR and GNPE were two most important traits which determine yield heterosis. Heterosis for GNPR and GNPE in the temperate diallel was 58.80 and 77.18%, respectively. Heterosis for RN and HGW (13.39 and 12.14%, respectively), was relatively low, even with negative effects (Supplementary Figures 2, 3). When compared across years and locations, heterosis for GY and yield-related traits were varied significantly. Correlation analysis of heterosis between the two environments, Shunyi and Xinxiang, indicates that heterosis for yield-related traits was significantly affected by environments. The correlation of heterosis between Shunyi and Xinxiang was more significant for GNPR and GNPE than for RN and HGW (Supplementary Figures 2, 3). As heterosis for RN and HGW was relatively low (Tables 3, 4), no steadily MPH or HPH across the environments was observed. When these observations are considered together, yield heterosis can be largely explained by two major yield-related components, GNPE and GNPR, which in most cases have to offset the negative effects of RN and HGW. By examining the top 5% hybrids ranked by MPH and HPH of GWPP and four representative commercial hybrids from China, Xianyu335 (P19 × P20), Zhengdan958 (P21 × P26), Yedan13 (P10 × P16), Ludan981 (P06 × P22) (Figure 6), significant yield heterosis (MHP and HPH) was observed in all the selected hybrids, and GNPE was a major contributor to yield heterosis. GNPR was also an important determinant for yield heterosis in most hybrids. In contrast, no consistent heterosis (MPH and HPH) across environments or hybrids was observed for RN and HGW. Thus, the results suggest that both MPH and HPH for GNPE and GNPR were highly stable or consistent across these high-yielding and commercial hybrids, whereas those for RN and HGW were case-dependent.
Figure 6. The yield-related heterosis for GWPP ranked in the top 5% hybrids and four representative commercial maize hybrids in a maize multiple-hybrid population. (A) mid-parent heterosis (MPH); (B) high-parent heterosis (HPH). The rank is in accordance with the order for heterosis of GWPP from high to low. RN, row number; GNPR, grain number per row; GNPE, grain number per ear; HGW, hundred grain weight; GWPP, grain weight per plant. Sample sizes: multiple-hybrid population (N = 724).
Significant Correlation Observed Between Heterosis and Combining Ability
Simple linear correlation coefficients were used to reveal the relationship among MPH, HPH, GCA, SCA, and F1 hybrid performance (Table 5). MPH was highly significantly correlated with SCA, and also positively correlated with GCA for all tested traits expect PH, EH, and ED. HPH was positively correlated with the sum of parental GCAs for all the tested traits except EH, DTA, ED, and GWPP. Hybrid performance showed stronger correlation with the sum of parental GCAs than with hybrid SCA. The strong correlation between heterosis (both MPH and HPH) and SCA suggests that SCA could be used to predict hybrid performance and heterosis. In contrast, no significant correlation was found between heterosis and GCA for most tested traits. The correlation of hybrid performance with SCA was higher than that with heterosis (MPH and HPH) for PH, EH, DTA, DTS, ED, and GWPP. We found that both MPH and HPH were significantly correlated with SCA, while their correlation with the sum of parental GCAs was not consistent across the tested traits. Meanwhile, we also found that both SCA and the sum of parental GCAs were highly significantly correlated with hybrid performance, and SCA was highly significantly correlated with heterosis (MPH and HPH).
Potential Application of Multiple-Hybrid Populations and Open-Source Breeding Programs
Combining ability and heterosis among maize lines and correlation between combining ability, heterosis and hybrid performance provide important insights for developing breeding strategies, defining heterotic groups, and predicting hybrid performance. Complex trait dissection and crop improvement for combining ability and heterosis need to use diverse germplasm resources and large populations. Although required phenotypic variability exists in diverse maize germplasm, most researches have been using relatively small sets of inbred lines (Fan et al., 2014; Badu-Apraku et al., 2015), largely due to the fact that the number of hybrids that can be produced increases exponentially with the increase of parental lines. The MHP used in this study, consisting of 724 hybrids, was developed with Griffing IV diallel and NCD II designs using temperate and tropical elite maize inbreds as parental lines, which is suitable for combining ability and heterosis analysis and can be used for breeding different ecotypes by taking the advantages of different maize germplasm resources. Both diallel and NCD II designs can provide detailed genetic information, including dominance-recessiveness relationships and genetic interactions. In this study, an MHP could be used to reveal useful information about combining ability, heterosis, hybrid performance and genotype × environment interaction.
On the other hand, the parental lines, which have been genotyped using 55K SNP markers and resequencing (Wang et al., 2017), can be shared with international collaborators. By inferring hybrid genotypes from their parental lines, various sets of MHPs can be developed for a specific target environment or research purpose, by sharing the 51 parental lines, from which any set of hybrids, up to 1275, can be generated. If the number of parental lines increases to 200, which is manageable for many breeding programs, up to 19,900 potential hybrids can be generated to meet the requirement of worldwide breeding programs. The genotyped parental lines can be used worldwide as proposed for open-source breeding programs (Xu et al., 2017, 2020).
Utilization of Combining Ability and Heterosis in Hybrid Maize Breeding
A full understanding of genetic basis of heterosis and combining ability remains elusive (Birchler, 2015), which, however, does not affect the vital role heterosis and combining ability play in maize breeding. Combining heterosis in different traits such as yield-related traits and stress tolerance could improve gain yield (Fujimoto et al., 2018). As shown in the present study, heterosis varied across environments as maize hybrids and inbred lines responded differently to environmental stimuli. Although heterosis was greatly affected by environmental variables, it is a quantifiable, trait-specific phenotype. In the temperate maize with six heterotic groups, the hybrids from different heterotic groups tended to exhibit trait-specific advantages. To be commercially advantageous, a hybrid should outperform its parents with respect to agronomic traits, especially the traits related to GY. Normally, GY heterosis is an important indicator of yield potential. In the present study, we demonstrated that heterosis was contributed mainly by two outperformed yield components, GNPE and GNPR. Moreover, heterosis was compromised in few cases by the negative effects of other component traits, RN and HGW. Therefore, GNPE and GNPR can be used as direct selection criteria for yield heterosis. Based on our analysis, nine elite inbred lines, HZ1, PH4CV, PH6WC, Qi319, Zheng58, TR0423, 18-599, CML504, and TR0582, had highly significant positive GCA effects for yield-related traits, which should have contributed to the improved hybrid yield. Six inbred lines, AS6103, PH4CV, PH6WC, TR0415, CML312SR, and CML330, manifested negative GCA effects for DTS and DTA, responsible for early flowering. The hybrids between heterotic groups exhibited various levels of dominance across traits, which were inconsistent among identified heterotic loci. In the elite maize hybrid Yuyu22 (Zong3 × Yu87-1), 13 heterotic loci were identified, including three for GY, seven for EL, one for RN and two for HGW (Tang et al., 2010). Several QTL were identified for seedling weight (SW), number of kernels per plant (NK), and GY in a cross between two elite inbred lines, B73 and H99 (Frascaroli et al., 2007). Using the elite maize hybrid Zhengdan958 (Zheng58 × Chang7-2), 38 heterotic loci for ear-related traits were identified, suggesting that the combination of heterotic loci in tested hybrids was genotype-dependent (Li et al., 2017). In another report, 156 QTL, 28 pairs of epistatic loci, and 10 QTL × environment interaction regions were identified, and the inheritance of yield-related traits and their MPH in Reid (PA) × Tem-tropic I (PB) hybrids with improved heterotic pattern is trait-dependent (Yi et al., 2019). The hybrids from different heterotic groups, Reid × SPT, Reid × LRC, SPT × PA, and Lancaster × LRC, showed highly significant positive SCA effect and heterosis. Considering combining ability (GCA and SCA), performance and heterosis together will help identify the hybrid combinations with comparative advantages for maize breeding. P05 × P21 (Lancaster × SPT), P07 × P21 (Reid × SPT), P32 × P23 (Tropical × SPT), P37 × P21 (Tropical × SPT), P39 × P21 (Tropical × SPT), P45 × P27 (Tropical × LRC), and P51 × P26 (Tropical × PA) contributed a highly significant positive SCA effect to GWPP and relatively high yield, MPH and HPH (Supplementary Table 3). As a result, these hybrids could be used in breeding for high yielding and wide-adaptability. Consequently, elite inbred lines with improved combining ability and associated heterotic patterns could be explored for efficient hybrid breeding.
Use of Diverse Germplasm From Different Ecotypes in Maize Breeding
Introduction of exotic and diverse germplasms into breeding programs is of great importance in broadening genetic variation and thus improving breeding efficiency. Tropical maize germplasm, which come from the location of maize origin, host rich genetic diversity that can be used for temperate maize breeding. With the development of molecular markers and high-efficient genotyping technologies, researches have been introgressing favorable alleles from tropical into temperate maize, as shown in GY, grain moisture content and lodging resistance (Lewis and Goodman, 2003) and resistance to maize lethal necrosis disease (Gowda et al., 2015). Introgression of temperate maize germplasm into tropical lines did not disrupt the existing heterotic groups, because the introgressed lines remained genetically inclined towards the original heterotic groups from which they were derived (Musundire et al., 2019). In the present study, we observed higher average heterosis in NCD II (temperate × tropical) hybrids than that for within-ecotype hybrids for all yield-related traits except HGW, with higher hybrid performance for EL, ED, and HGW. Introgression of favorable genes and alleles from tropical maize germplasm can be explored to broaden the genetic basis of temperate maize, improve biotic and abiotic stress tolerance, optimize heterotic patterns, and develop improved temperate-tropical hybrids (Teixeira et al., 2015). Thus, large-scale analysis of combining ability and heterosis can facilitate improving targeted breeding efficiency using different ecotypes by broadening the genetic base of commercial hybrids.
Prediction of Heterosis and Hybrid Performance Using Combining Ability
Prediction of hybrid performance and heterosis is one of the breeders’ dreams in hybrid breeding. Although there are still some gaps in our understanding of the mechanism of heterosis, great progress has been made in predicting hybrid performance (Andorf et al., 2019). The prediction can be done through integrated information including the performance of parental lines per se, combining ability and heterosis along with the phenotypic data collected from field evaluation (Josue and Brewbaker, 2018). In addition, different types of omics data can be explored to improve the prediction in maize (Schrag et al., 2018). There was no association found between the performance of parental lines and heterosis in F1 hybrids, because hybrid performance depends on the nature of genetic variation (Lee et al., 2007). Hence, it can be deduced that performance of parents per se may not be necessarily a reliable indicator of heterosis prediction. Heterosis can be treated as a single trait in genome prediction of heterosis. We observed that SCA was significantly correlated with MPH and HPH for all the tested traits, supporting that non-additive effects are the main effect for heterosis in maize (Falconer and Mackay, 1996). In contrast, the sum of parental GCAs was either negatively or not correlated with MPH or HPH. However, the correlation of hybrid performance with the sum of parental GCAs was higher than that with heterosis. We also found that hybrid performance had stronger correlation with the sum of parental GCAs than that with SCA for all the tested traits. The GCA/(GCA + SCA) ratio indicates that all the tested traits were affected by both additive and non-additive effects, but the additive effect was more important, as reported in several maize studies for kernel rows per ear (Dehghanpour and Ehdaie, 2013), PH (Pswarayi and Vivek, 2008) and GY (Sibiya et al., 2013). Therefore, the GCA/(GCA + SCA) ratio, which is close to one, suggests that GCA is much more important than SCA in hybrid maize breeding. The inbred lines with high GCA can be used as parents to develop superior hybrids in maize breeding. Strong correlation between GCA effects and hybrid performance per se for most traits suggests that parental GCA may be a good indicator for predicting hybrid performance. It has been reported that the accuracy of predicted SCA effects was considerably lower than that of GCA effects (Technow, 2019). As discussed above, hybrid performance was dependent on the effects of GCA and SCA, while heterosis was dependent on SCA effects. The significant and positive correlation among SCA, MPH and HPH indicates that SCA can be used to predict MPH and HPH, which is also supported by several previous studies (Solomon et al., 2012; Ndhlela et al., 2015; Tian et al., 2017). The genetic information generated from a subset of the hybrids developed in the present study can be used to predict the whole set of hybrids with specific heterotic patterns. In genomic selection, the tested parental lines and their hybrids can be used as a training population to predict the rest potential hybrids.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
YX conceived and designed the study. KY, HW, XL, and CX analyzed the data. KY, HW, XL, CX, ZL, XX, and JL performed the experiments. KY, ZW, and YX wrote the manuscript.
National Key Research and Development Program of China (2016YFD0101803); The Agricultural Science and Technology Innovation Program (ASTIP) of CAAS; Fundamental Research Funds for Central Non-Profit of Institute of Crop Sciences, CAAS (1610092016124); and The Bill and Melinda Gates Foundation and the CGIAR Research Program MAIZE.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00660/full#supplementary-material
Akinwale, R. O., Badu-Apraku, B., Fakorede, M. A. B., and Vroh-Bi, I. (2014). Heterotic grouping of tropical early-maturing maize inbred lines based on combining ability in Striga-infested and Striga-free environments and the use of SSR markers for genotyping. Field Crops Res. 156, 48–62. doi: 10.1016/j.fcr.2013.10.015
Andorf, C., Beavis, W. D., Hufford, M., Smith, S., Suza, W. P., Wang, K., et al. (2019). Technological advances in maize breeding: past, present and future. Theor. Appl. Genet. 132, 817–849. doi: 10.1007/s00122-019-03306-3
Badu-Apraku, B., Annor, B., Oyekunle, M., Akinwale, R. O., Fakorede, M. A. B., Talabi, A. O., et al. (2015). Grouping of early maturing quality protein maize inbreds based on SNP markers and combining ability under multiple environments. Field Crops Res. 183, 169–183. doi: 10.1016/j.fcr.2015.07.015
Brown, P. J., Upadyayula, N., Mahone, G. S., Tian, F., Bradbury, P. J., Myles, S., et al. (2011). Distinct genetic architectures for male and female inflorescence traits of maize. PLoS Genet. 7:e1002383. doi: 10.1371/journal.pgen.1002383
Comstock, R. E., and Robinson, H. F. (1948). The components of genetic variance in populations of biparental progenies and their uses in estimating the average degree of dominance. Biometrics 4, 254–266. doi: 10.2307/3001412
Cui, C., Gao, J., Yu, X., Su, Z., Wang, Z., Sun, J., et al. (2014). Combining ability of traits related to nitrogen use efficiency in eighteen maize inbred lines. Acta Agron. Sin. 40, 838–849. doi: 10.3724/sp.j.1006.2014.00838
Dehghanpour, Z., and Ehdaie, B. (2013). Stability of general and specific combining ability effects for grain yield in elite Iranian maize inbred lines. J. Crop Improv. 27, 137–152. doi: 10.1080/15427528.2012.745822
Fan, X., Yin, X., Zhang, Y., Bi, Y., Liu, L., Chen, H., et al. (2016). Combining ability estimation for grain yield of maize exotic germplasm using testers from three heterotic groups. Crop Sci. 56, 2527–2535. doi: 10.2135/cropsci2016.01.0032
Fan, X., Zhang, Y., Yao, W., Bi, Y., Liu, L., Chen, H., et al. (2014). Reciprocal diallel crosses impact combining ability, variance estimation, and heterotic group classification. Crop Sci. 54, 89–97. doi: 10.2135/cropsci2013.06.0393
Frascaroli, E., Canè, M. A., Landi, P., Pea, G., Gianfranceschi, L., Villa, M., et al. (2007). Classical genetic and quantitative trait loci analyses of heterosis in a maize hybrid between two elite inbred lines. Genetics 176, 625–644. doi: 10.1534/genetics.106.064493
Fujimoto, R., Uezono, K., Ishikura, S., Osabe, K., Peacock, W. J., and Dennis, E. S. (2018). Recent research on the mechanism of heterosis is important for crop and vegetable breeding systems. Breed. Sci. 68, 145–158. doi: 10.1270/jsbbs.17155
Gowda, M., Das, B., Makumbi, D., Babu, R., Semagn, K., Mahuku, G., et al. (2015). Genome-wide association and genomic prediction of resistance to maize lethal necrosis disease in tropical maize germplasm. Theor. Appl. Genet. 128, 1957–1968. doi: 10.1007/s00122-015-2559-0
Lee, E. A., Ash, M. J., and Good, B. (2007). Re-examining the relationship between degree of relatedness, genetic effects, and heterosis in maize. Crop Sci. 47, 629–635. doi: 10.2135/cropsci2006.04.0275
Lewis, R. S., and Goodman, M. M. (2003). Incorporation of tropical maize germplasm into inbred lines derived from temperate× temperate-adapted tropical line crosses: agronomic and molecular assessment. Theor. Appl. Genet. 107, 798–805. doi: 10.1007/s00122-003-1341-x
Li, H., Yang, Q., Gao, L., Zhang, M., Ni, Z., and Zhang, Y. (2017). Identification of heterosis-associated stable QTLs for ear-weight-related traits in an elite maize hybrid Zhengdan 958 by design III. Front. Plant Sci. 8:561. doi: 10.3389/fpls.2017.00561
Li, Z., Luo, L., Mei, H., Wang, D., Shu, Q., Tabien, R., et al. (2001). Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield. Genetics 158, 1737–1753. doi: 10.3410/f.1002133.8203
Luo, L., Li, Z., Mei, H., Shu, Q., Tabien, R., Zhong, D., et al. (2001). Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. II. Grain yield components. Genetics 158, 1755–1771.
Machida, L., Derera, J., Tongoona, P., and MacRobert, J. (2010). Combining ability and reciprocal cross effects of elite quality protein maize inbred lines in subtropical environments. Crop Sci. 50, 1708–1717. doi: 10.2135/cropsci2009.09.0538
Makumbi, D., Assanga, S., Diallo, A., Magorokosho, C., Asea, G., Worku, M., et al. (2018). Genetic analysis of tropical midaltitude-adapted maize populations under stress and nonstress conditions. Crop Sci. 58, 1492–1507. doi: 10.2135/cropsci2017.09.0531
Meseka, S., Williams, W. P., Warburton, M. L., Brown, R. L., Augusto, J., Ortega-Beltran, A., et al. (2018). Heterotic affinity and combining ability of exotic maize inbred lines for resistance to aflatoxin accumulation. Euphytica 214:184. doi: 10.1007/s10681-018-2254-8
Miller, M., Song, Q., Shi, X., Juenger, T. E., and Chen, Z. J. (2015). Natural variation in timing of stress-responsive gene expression predicts heterosis in intraspecific hybrids of Arabidopsis. Nat. Commun. 6:7453. doi: 10.1038/ncomms8453
Musundire, L., Derera, J., Dari, S., Tongoona, P., and Cairns, J. E. (2019). Molecular characterisation of maize introgressed inbred lines bred in different environments. Euphytica 215:46. doi: 10.1007/s10681-019-2367-8
Ndhlela, T., Herselman, L., Semagn, K., Magorokosho, C., Mutimaamba, C., and Labuschagne, M. T. (2015). Relationships between heterosis, genetic distances and specific combining ability among CIMMYT and Zimbabwe developed maize inbred lines under stress and optimal conditions. Euphytica 204, 635–647. doi: 10.1007/s10681-015-1353-z
Pswarayi, A., and Vivek, B. (2008). Combining ability amongst CIMMYT’s early maturing maize (Zea mays L.) germplasm under stress and non-stress conditions and identification of testers. Euphytica 162, 353–362. doi: 10.1007/s10681-007-9525-0
Riedelsheimer, C., Czedik-Eysenberg, A., Grieder, C., Lisec, J., Technow, F., Sulpice, R., et al. (2012). Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat. Genet. 44, 217–220. doi: 10.1038/ng.1033
Schrag, T. A., Westhues, M., Schipprack, W., Seifert, F., Thiemann, A., Scholten, S., et al. (2018). Beyond genomic prediction: combining different types of omics data can improve prediction of hybrid performance in maize. Genetics 208, 1373–1385. doi: 10.1534/genetics.117.300374
Sibiya, J., Tongoona, P., and Derera, J. (2013). Combining ability and GGE biplot analyses for resistance to northern leaf blight in tropical and subtropical elite maize inbred lines. Euphytica 191, 245–257. doi: 10.1007/s10681-012-0806-x
Solomon, K. F., Zeppa, A., and Mulugeta, S. D. (2012). Combining ability, genetic diversity and heterosis in relation to F1 performance of tropically adapted shrunken (sh2) sweet corn lines. Plant Breed. 131, 430–436. doi: 10.1111/j.1439-0523.2012.01965.x
Tang, J., Yan, J., Ma, X., Teng, W., Wu, W., Dai, J., et al. (2010). Dissection of the genetic basis of heterosis in an elite maize hybrid by QTL mapping in an immortalized F2 population. Theor. Appl. Genet. 120, 333–340. doi: 10.1007/s00122-009-1213-0
Teixeira, J. E. C., Weldekidan, T., De Leon, N., Flint-Garcia, S., Holland, J. B., Lauter, N., et al. (2015). Hallauer’s Tusón: a decade of selection for tropical-to-temperate phenological adaptation in maize. Heredity 114, 229–240. doi: 10.1038/hdy.2014.90
Vasal, S. K., Srinivasan, G., Crossa, J., and Beck, D. L. (1992). Heterosis and combining ability of CIMMYT’s subtropical and temperate early-maturity maize germplasm. Crop Sci. 32, 884–890. doi: 10.2135/cropsci1992.0011183X003200040010x
Wang, H., Xu, C., Liu, X., Guo, Z., Xu, X., Wang, S., et al. (2017). Development of a multiple-hybrid population for genome-wide association studies: theoretical consideration and genetic mapping of flowering traits in maize. Sci. Rep. 7:40239. doi: 10.1038/srep40239
Xu, Y., Liu, X., Fu, J., Wang, H., Wang, J., Huang, C., et al. (2020). Enhancing genetic gain through genomic selection: from livestock to plants. Plant Commun. 1:100005. doi: 10.1016/j.xplc.2019.100005
Yi, Q., Liu, Y. H., Hou, X. B., Zhang, X. G., Li, H., Zhang, J. J., et al. (2019). Genetic dissection of yield-related traits and mid-parent heterosis for those traits in maize (Zea mays L.). BMC Plant Biol. 19:392. doi: 10.1186/s12870-019-2009-2
Keywords: maize, multiple-hybrid population, heterosis, heterotic groups, combining ability
Citation: Yu K, Wang H, Liu X, Xu C, Li Z, Xu X, Liu J, Wang Z and Xu Y (2020) Large-Scale Analysis of Combining Ability and Heterosis for Development of Hybrid Maize Breeding Strategies Using Diverse Germplasm Resources. Front. Plant Sci. 11:660. doi: 10.3389/fpls.2020.00660
Received: 19 December 2019; Accepted: 28 April 2020;
Published: 01 June 2020.
Edited by:Hanwei Mei, Shanghai Agrobiological Gene Center, China
Reviewed by:Haiming Xu, Zhejiang University, China
Yongzhong Xing, Huazhong Agricultural University, China
Copyright © 2020 Yu, Wang, Liu, Xu, Li, Xu, Liu, Wang and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.