Comparative proteomic analysis reveals that the Heterosis of two maize hybrids is related to enhancement of stress response and photosynthesis respectively
BMC Plant Biology volume 21, Article number: 34 (2021)
Heterosis refers to superior traits exhibiting in a hybrid when compared with both parents. Generally, the hybridization between parents can change the expression pattern of some proteins such as non-additive proteins (NAPs) which might lead to heterosis. ‘Zhongdan808’ (ZD808) and ‘Zhongdan909’ (ZD909) are excellent maize hybrids in China, however, the heterosis mechanism of them are not clear. Proteomics has been wildly used in many filed, and comparative proteomic analysis of hybrid and its parents is helpful for understanding the mechanism of heterosis in the two maize hybrids.
Over 2000 protein groups were quantitatively identified from second seedling leaves of two hybrids and their parents by label-free quantification. Statistical analysis of total identified proteins, differentially accumulated proteins (DAPs) and NAPs of the two hybrids revealed that both of them were more similar to their female parents. In addition, most of DAPs were up-regulated and most of NAPs were high parent abundance or above-high parent abundance in ZD808, while in ZD909, most of DAPs were down-regulated and most of NAPs were low parent abundance or below-low parent abundance. Pathway enrichment analysis showed that more of stress response-related NAPs in ZD808 were high parent abundance or above-high parent abundance, and most of PS related NAPs in ZD909 were high parent abundance or above-high parent abundance. Finally, four stress response-related proteins and eight proteins related to PS were verified by PRM, ten of them had significant differences between hybrid and midparent value.
Even though every one of the two hybrids were more similar to its female parent at proteome level, the biological basis of heterosis is different in the two maize hybrids. In comparison with their parents, the excellent agronomic traits of hybrid ZD808 is mainly correlated with the high expression levels of some proteins related to stress responses and metabolic functions, while traits of ZD909 is mainly correlated with high expressed proteins related to photosynthesis. Our proteomics results support previous physiological and morphological research and have provided useful information in understanding the reason of valuable agronomic traits.
Nearly 200 NAPs were identified in maize hybrids ZD808 and ZD909 respectively.
Most heterosis related proteins and pathways were different between two hybrids.
Pathways related to stress response were significantly enhanced in hybrid ZD808.
Pathways related to photosynthesis were significantly enhanced in hybrid ZD909.
Heterosis refers to the phenomenon that of superior traits that exhibiting in hybrid progeny of species or varieties compared with both parents. This can include important agronomic traits, such as biomass, developmental speed, fertility, yield, and abiotic/biotic stress tolerance. Thus, F1 hybrids displaying heterosis are widely used in the production of many crops and vegetables [1, 2]. However, to obtain a hybrid with heterosis, a large number of field tests are needed, which are labor intensive, expensive, time consuming and long-term processes, and the results are unpredictable. To improve the efficiency of crossbreeding, it is necessary to understand the molecular mechanisms of heterosis.
Some methods such as physiological, genetic, molecular, transcriptomic and proteomic approaches, have been used to study heterosis for many years [1, 3, 4]. Among of them, physiological studies revealed that heterosis is affected, to some extent, by the magnitude and rate of vegetative growth, the flowering time, the yield, and the resistance levels to biotic and abiotic environmental stresses . A genetic study showed that one showing primarily nonadditive gene action explained most of the heterotic variation of grain weight in rice . Recently, much attention has been paid to heterosis associated genes [1, 7]. As we known, various genes associated with heterosis have been located or been deeply studied, such as GA-related genes , ErbB-3 binding Protein 1 , calcium-dependent protein kinase [10, 11], xyloglucan endotransglucosylase/hydrolase , cell number regulator gene . With the development of life sciences-related technology, some new methods have also been used to study heterosis, including transcriptomic methods, especially proteomics [1, 4]. Transcriptomics technologies also have been used to compare parental lines with their F1 hybrids to identify genes potentially involved in heterosis, and fond that the mode of gene action in F1 hybrids is classified as having additive and non-additive expression . In proteomics analysis of crop hybrids, some heterosis-related differentially accumulated proteins (DAPs) were identified and non-additive proteins (NAPs) were believed as the important DAPs [15, 16]. In analysis of heterosis-related proteins, proteins that deviated from the midparent value are termed non-additive proteins (NAPs), while others are termed additive proteins (APs). The concept of NAPs is widely used in heterosis research and NAPs are generally classed into four types, ‘above-high parent abundance’ (++), ‘high parent abundance’ (+), ‘below-low parent abundance’ (− −) and ‘low parent abundance’ (−) [16, 17]. Above-high or below-low parent abundance indicates that the NAP content in hybrid is significantly greater or lower than any of parents, high or low parent abundance indicates that the NAP content in hybrid is significantly greater or lower than one of parental line but not another parental line.
Maize is one of the most important cereal crops worldwide and maize hybrids with heterosis are wildly planted in the world . Quantitative trait locus analyses at earlier stage concluded that overdominance is the main cause of heterosis , and at recently, multiparent populations have been used to identify quantitative trait locus associated with the genetic architecture of hybrid performances [20, 21]. Transcriptome analyses also have been performed using seedlings, roots, embryos, and immature ears in maize [22, 23]. Calcium-dependent protein kinases are calcium-binding serine/threonine protein kinases, play essential roles in cellular calcium signaling processes. A putative CDPK gene, ZmCPK16 was reported to had different expressions between hybrids and their parents in maize, indicated it plays a part in heterosis . Xyloglucan endotransglucosylase/hydrolase play an essential role in the construction and restructuring of xyloglucan cross-links, a family of xyloglucan modifying enzymes , a xyloglucan endotransglucosylase/hydrolase (GRMZM2G091118) down-regulated at V3 stem and shoot apical meristem, although it up-regulated at 6 days after sowing primary root, and had consistent expression at others tissues in different heterosis groups . The cell number regulator super family negative regulators across diverse species. Cell number regulator 2 expression was found to be negatively correlated with tissue growth activity and hybrid seedling vigor in maize . Several heterosis-related genes in maize have been analyzed, but it is not easy to understand the causes or the results of heterosis.
In addition to using nucleic acid research methods, proteomic approaches have successfully identified some maize heterosis-related DAPs in different tissues and/or stages. In earlier research, by using two-dimensional gel electrophoresis (2DE) and ESI-MS/MS, it was found that NAP accumulation in primary roots of maize hybrids might be associated with the manifestation of heterosis . In 2010, 2DE and ESI-MS/MS were also used to study NAP accumulation patterns in maize hybrids during embryo development . An analysis of the heterosis of ZD909 determined that half of the DAPs in immature maize ears were related to dominance expression patterns, and high parental dominance proteins mainly participated in carbon metabolism and nitrogen assimilation processes . In 2018, a proteomics analysis of DAPs in the primary roots of the popcorn hybrid UENF/UEM01 showed that the heterosis was associated with an up-regulation of proteins related to synthesis and energy metabolism at an early stage of plant development . A proteomics study of auxin homeostasis’ influence on the heterosis of the eighth internodal length in maize found that proteins involved in phenylpropanoid and benzoxazinoid metabolic pathways were correlated with the extent of the heterosis . Recently, 2DE and MALDI-TOF-MS were also used to analyze the seed germination trait of hybrid DHM 117 and its parental inbred lines. The DAPs between the hybrid and its parental inbred lines were related to metabolism and energy, followed by storage proteins, stress response, transcription and translation . These results indicated that hybridization between parents can change the expression types of some proteins in maize hybrids, which might lead to heterosis. Understanding the proteomic characteristics of maize hybrids and their parents, and determining the protein types associated with the advantages of hybrids, may provide important guidance for efficient breeding.
‘Zhongdan808’ (ZD808) and ‘Zhongdan909’ (ZD909) are excellent maize hybrids with high yield, good quality, multi-resistance, and wide adaptability. The hybrid ZD808 was bred using CLll as the female parent and NG5 as the male parent, and it is mainly planted in southwestern China. Another hybrid, ZD909, was bred using Z58 as the female parent and HD568 as the male parent, and it is mainly planted in the Huang-Huai-Hai area of China. The two hybrids have been grown on over 6.67 million hectares of cropland in the past in China. They both exhibit strong heterotic phenotypes compared with their parents for agronomic traits, such as more vigorous seedlings, larger ears, bigger grain size, and increased stress resistance. Under water stress conditions, ZD808 had a higher relative water content, lower malondialdehyde content and better osmotic adjustment ability compared with other varieties . However, the mean leaf area index and duration, as well as the mean net assimilation rate, were higher in ZD909 compared with other varieties . In the present study, we used LC/MS/MS to analyze DAPs and NAPs in seeding leaves of hybrids ZD808 and ZD909. This may provide useful information for understanding of the molecular mechanisms involved in the heterosis of hybrid maize.
By using a comparative proteomics approach (Fig. 1), proteomes in seeding leaves of maize hybrids ZD808 and ZD909, and their parents were analyzed and 12 heterosis-related proteins were verified.
LC-MS analysis and statistical analysis of quantitative identified proteins
The total ion chromatograms of those samples showed that the nano liquid phase peaked within 10 min, and the distribution of the chromatographic peaks was uniform and smooth (Fig. S1, S2, S3, S4, S5 and S6: Total ion chromatograms of all the samples). Thus, the mass spectrometry system appeared stable, the parameters were suitable, and the repeatability of the sample preparation was good. Mass spectrometry raw files from three biological repeats of six varieties were used for the quantitative identification through a label-free MS-based qualitative proteomics technique. A total of 2222 and 2486 proteins were quantitatively identified from ZD808 and its parents (Table S1) and from ZD909 and its parents (Table S2), respectively. According to the Pearson’s correlation analysis of the LFQ quantitative data, the correlations between the three independent biological replicates of most varieties were high (Fig. 2a, b). In the ZD808 group, except for the first biological replicate of ZD808, the other two biological replicates had high similarity levels with the female parent and low similarity levels with the male parent. The difference between the ZD808 male and female parents was greater than that between the hybrid and the parents (Fig. 2a). A similar situation was found in the ZD909 sample group (Fig. 2b). Thus, the protein expression characteristics of the two hybrids came more from their female parents than their male parents. In addition, principal component analysis was conducted in all of samples, the PC1 and PC2 represented 94.22 and 2.66% of the total variance, respectively, thus making 96.88% for the total variance accounted for in the principal component analysis model (Fig. S7).
To identify the DAPs among samples, a 2-fold change was chosen as an empirical cutoff. In total, there were 833 up-regulated proteins and 459 down-regulated proteins in hybrid ZD808 when compared with its female parent CL11 or male parent NG5, and there were 537 up-regulated proteins and 1167 down-regulated proteins in hybrid ZD909 when compared with its female parent Z58 or male parent HD568 (Fig. 2c). Among the DAPs, there were 325 common up-regulated proteins and 54 common down-regulated proteins in ZD808 when compared with both parents, and there were 309 common up-regulated proteins and 143 common down-regulated proteins in ZD909 when compared with both parents (Fig. 2d, e). NAPs may be highly associated with heterosis. There were 188 and 263 NAPs identified in the seedlings of hybrids ZD808 and ZD909, respectively. Among them, 65 proteins were classified as ‘++’, 80 as ‘+’, 23 as ‘− −’, and 20 as ‘−’ in ZD808 (Table S3). While in ZD909, 22 proteins were classified as ‘++’, 52 were classified as ‘+’, 81 as ‘− −’, and 108 as ‘−’ (Table S4). Among all the NAPs identified from ZD808 and ZD909, only four were identified at the same time in both hybrids (Fig. 2f). The NAPs in two hybrids were clustered, it can be clearly seen that the distribution of NAPs between two hybrids is difference, and the difference between male and female parents is larger than the difference between hybrid and one of the parents in ZD808, the same to ZD909 (Fig. 2g, h).
A comparison of the DAPs in ZD808 and ZD909 showed that the number of common up-regulated DAPs was greater than common down-regulated DAPs. Among the NAPs, the number of above-high and high parental abundance NAPs was greater than below-low and low parental abundance NAPs in ZD808, but to the opposite was true in ZD909. A comparison of the DAPs and NAPs in ZD808 and ZD909 showed that the number of common up-regulated proteins in ZD808 when compared with both parents was higher than the number of NAPs in ZD808 classified as ‘++’ and ‘+’. This was the same for the common down-regulated proteins in both ZD808 and ZD909.
Pathway enrichment analysis of DAPs and NAPs
The pathways of DAPs and NAPs identified from ZD808 and ZD909 were enriched using MapMan software and divided into 34 classes (Fig. 3). DAPs with at least 2-fold change in the hybrids when compared to any parent, were highly enriched in protein, RNA, amino acid metabolism, signaling, lipid metabolism, secondary metabolism, Photosynthesis (PS), stress, redox, hormone metabolism, and development. up-regulated DAPs with at least 2-fold change in the hybrids when compared to both parents were commonly enriched in many pathways, except development and tetrapyrrole synthesis pathways, which were only enriched in ZD808. down-regulated DAPs with at least 2-fold change in the hybrids when compared to both parents were associated with protein pathways, and in ZD909 were especially enriched in secondary metabolism. NAPs of ZD808 classified as ‘++’ and ‘+’ were associated with protein, RNA, miscellaneous, secondary metabolism, amino acid metabolism, redox, signaling, and tricarboxylic acid (TCA)/org transformation pathways, while NAPs of ZD909 classified as ‘++’ and ‘+’ were mainly associated with the PS pathway. NAPs of ZD808 classified as ‘− −‘and ‘−’ were associated with protein, PS, and miscellaneous pathways. NAPs of ZD909 classified as ‘− −‘and ‘−’ were associated with protein, RNA, amino acid metabolism, transport, secondary metabolism, tetrapyrrole synthesis, miscellaneous, and redox pathways.
Stress related pathways enhanced in ZD808
A more detailed analysis showed that a total of 76 NAPs from ZD808 and ZD909 could be enriched in stress response-related pathways (Fig. 4), including secondary metabolism, proteolysis, TCA, abiotic stress, redox state, signaling and so on. Among them, 44 NAPs were identified from ZD808 and its parents, and 32 NAPs were up-regulated in ZD808. At the same time, 36 NAPs were identified from ZD909 and its parents, and only 11 NAPs were up-regulated in ZD909.
PS related pathways enhanced in ZD909
A total of 27 NAPs from ZD808 and ZD909 could be enriched in PS-related pathways (Fig. 5), including light reactions, Calvin cycle and photorespiration. Most of PS-related NAPs were identified from ZD909 and its parents, and many of them were up-regulated, except 2 NAPs were down-regulated in ZD909, they were mainly involved in light reactions. Ten of them were identified from ZD808 and its parents, and only 3 NAPs were up-regulated in ZD808, them were mainly enriched in light reactions and Calvin cycle. The result showed that half of the NAPs of ZD808 involved in PS were classified as ‘− −’ or ‘−’, and unlike in ZD808, most of the NAPs related to PS in ZD909 were classified as ‘+’ or ‘++’.
The PRM assay was also used to verify a subset of the NAPs obtained in the label-free analysis. Because the signature peptides for target proteins must exhibit uniqueness, only proteins that possessed unique peptide sequences were selected for the PRM analysis. Some proteins related to the PS and stress pathways were selected for PRM verification (Fig. 6). The Skyline analysis results of unique peptides of 12 proteins verified by PRM are shown in supplementary files (Fig.S8, S9, Table S7), and the detail data results including peak area and retention time of each peptide fragment and the peptide information are shown in supplementary files (Table S5, S6). The verification results showed that the intensity of most PS-related proteins in hybrid ZD808 were lower than in its parents, while the opposite was true in ZD909.
Heterosis is ubiquitous and has attracted much scientific attention. The mechanism of heterosis is very complex. Many methods, like genomics and transcriptomics, have been used in heterosis studies , and now proteomics is also being used. ZD808 is the main cultivated variety in southwestern China, which has a complicated eco-environment with barren soil and frequent drought and disease occurrences, while ZD909 is the main cultivated variety in the Huang-Hui-Hai area of China. ZD808 can adapt to adverse environmental conditions and exhibits excellent resistance levels to stresses, such as drought , while ZD909 has a strong photosynthetic capacity . So that maize hybrids ZD808 and ZD909 have excellent yield traits compared with those of their parents [26, 29] (Table 1). Here, proteomic techniques were used to analyze the two hybrids and their parental lines to provide new clues regarding the heterosis mechanisms. However, further exploration is needed in the aspects of the breadth, proteomics technology and gene function verification.
Heterosis is affected, to some extent, by many factors [3, 5], and NAPs may be a major source of heterosis. By using LFQ, 2222 and 2486 proteins were quantified from seedling leaves of ZD808 group and ZD909 group, respectively, means that there is a difference in the analysis results between groups. The Pearson’s correlation analysis revealed that the correlations of the quantitative proteins of three independent biological replicates were relatively high, and the two hybrids were more similar to their female parents. In our experiment, there has some difference exist in the proteome analysis results of ZD808–1 and the other two biological repeat samples. However, all the data of ZD808 biological repeats were more similar to the female parent but not to the male parent. The results showed that the correlation pattern of ZD808–1 and the other two biological repeats were similar (Fig. 2a, b). Principal component analysis also showed that correlation pattern of two hybrids were more similar to their female parents, and the correlations of ZD909 group was higher than ZD808 group (Fig. S7). Venn diagrams independently illustrated the relationships of DAPs and NAPs between two hybrids. Some DAPs of ZD808 and ZD909 were simultaneously identified as up-regulated proteins and/or down-regulated proteins when compared with any parent, but there were no significant differences between the numbers of simultaneously identified DAPs in most intersections of the two data sets (Fig. 2c). Although many DAPs in ZD808 were up-regulated and DAPs in ZD909 were down-regulated, the comparative analysis of common DAPs indicated that there were more common up-regulated DAPs than common down-regulated DAPs in both ZD808 and ZD909 (Fig. 2d, e). However, the NAP analysis showed that their patterns were significantly different in ZD808 and ZD909 (Fig. 2f). The Cluster analysis showed a difference distribution of NAPs between two hybrids, but similarity of two hybrids and their female parents is higher than hybrids and their male parents, it indicated female parent had more influence on hybrids NAPs.
All the identified proteins were classified into 33 pathways as well as a group in which proteins were not annotated (Fig. 3). The pathway enrichment analysis showed that most of the up-regulated DAPs in ZD808 and down-regulated DAPs in ZD909, when compared with any parent, were enriched in many pathways. In ZD808, DAPs up-regulated by at least 2-fold change in the hybrids when compared with both parents were enriched in almost all of the pathways, especially stress-related pathways such as amino acid metabolism, signaling, lipid metabolism, secondary metabolism, stress, redox, and hormone metabolism. In ZD909, DAPs up-regulated by at least 2-fold change in the hybrids when compared with both parents were mainly enriched in cell, stress, and mitochondrial electron transport/ATP synthesis and transport. Thus, the functions of the DAPs in ZD808 and ZD909 appeared to be different. However, in ZD808, most of the ‘++’ and ‘+’ NAPs were only enriched in stress-related pathways, while ‘− −’ and ‘−’ NAPs were only obviously enriched in PS pathways. However, in ZD909, although ‘− −’ and ‘−’ NAPs were a large majority of all the NAPs, the ‘++’ and ‘+’ NAPs were highly enriched in PS pathway, while ‘− −’ and ‘−’ NAPs were enrichment in RNA, amino acid metabolism, secondary metabolism, and redox. Thus, stress-related pathways were enhanced in ZD808, and PS-related pathways were enhanced in ZD909.
Because non-additively expressed proteins were predicted to play important roles in heterosis [3, 16, 31], all the identified NAPs were further analyzed using MapMan (Table S3, S4). Most of the NAPs identified from the two hybrids could be divided into four classes of pathways, stress responses, PS, amino acid metabolism and protein. We believe that if the majority of the NAPs in a pathway were up- or down-regulated, then the pathway was enhanced or impeded, respectively. Stress responses play crucial roles when plants are subjected to abiotic or biotic stresses. The MapMan analysis showed that some up-regulated NAPs of ZD808 were enriched in pathways involved in stress responses (Fig. 4). Among these proteins, 70-kDa heat shock proteins (Hsp70s) are a widely expressed family with similar structures that exist in almost all organisms. Hsp70s are important components of the cell’s machinery for protein folding, and they help to protect cells from cold, heat, and drought stresses [32,33,34]. Other proteins shown in Fig. 4 are also related to basal metabolism and stress responses. The basal metabolism and stress responses were enhanced in hybrid ZD808 compared with its parents. Some studies reported the relationship between heterosis and stress tolerance, but a great difference exist on different species and different stress . Some evidence showed that there is negative relationship between expression of defense response genes and growth heterosis , there are also some evidences of a positive or no relationship between heterosis and stress tolerance [36, 37]. In this study, stress responses were enhanced at three-leaf stage, and yield increased at maturity stage, it suggested that there had a positive relationship between heterosis and anti-stress in ZD808. Increased vigor or yield is usually the focus of heterosis research. PS is a key process in promoting plant growth and development. The main PS-related pathways enriched with NAPs of ZD808 and ZD909 are shown in Fig. 5. There is a great difference in the PS activity between the two hybrids at the three-leaf stage. The enhanced of PS-related pathways at three-leaf stage could be the bases on heterosis of yield at maturity stage.
Amino acid metabolism also plays important roles in the growth processes of plants. A total of 14 NAPs in ZD808 and ZD909 were enriched in amino acid metabolism-related pathways. Among them, most proteins identified as ‘++’ or ‘+’ in ZD808, except D-3-phosphoglycerate dehydrogenase. However, in ZD909, all of these NAPs were classified as ‘− −’ or ‘−’. It indicated that amino acid metabolism was enhanced in hybrid ZD808 and impeded in ZD909.
Moreover, protein pathways were greatly affected by hybridization, including amino acid activation, protein degradation (serine protease and ubiquitin), protein folding, protein synthesis (elongation, initiation and ribosomal protein), protein targeting (in the chloroplast and nucleus), assembly and cofactor ligation, and posttranslational modification. In ZD808, most of the NAPs related to protein synthesis, folding, amino acid activation, and targeting were classified as ‘+’ or ‘++’. However, in ZD909, most of the NAPs in protein-related pathways were classified as ‘− −’ or ‘−’. This indicated that the effects of protein synthesis and other protein-related pathways in the two hybrids were very different.
In the PRM verification of some key proteins, the abundance levels of 12 NAPs related to PS and stress were analyzed (Fig. 6, Table S7). The PRM results in parents and hybrids were similar to the label-free quantitative proteomics results. Most of them had significant differences between hybrid and midparent value both in LFQ data and PRM verification.
Four stress response-related proteins, chaperone protein ClpB3, D-3-phosphoglycerate dehydrogenase (PGDH), glutathione S-transferase 6 (GST6), and membrane steroid-binding protein 1 (MSBP1), were verified by PRM. The direct interaction of ClpB3 with Hsp70 may eventually result in the refolding, and hence the reactivation, of deoxyxylulose 5-phosphate synthase, which promotes the repair of the methylerythritol 4-phosphate pathway [38,39,40]. This pathway synthesizes the metabolic precursors of isoprenoids, such as carotenoids, and the prenyl chains of chlorophylls, tocopherols, or plastoquinone. PRM verification classified ClpB3 as ‘++’ in ZD808, while it was classified as ‘+’ in label-free quantitative proteomics, there had significant difference between ZD808 and midparent value in PRM verification. The high abundance level of ClpB3 in hybrid ZD808 indicates that its stress-response capability is greater than that of its parents. PGDH, which originates from glycolysis and the Calvin cycle, participates in the biosynthesis of L-serine . Arabidopsis thaliana PGDH1 is also essential for plant adaptation to high CO2, as well as for embryogenesis and pollen development . Recently, a study showed that PGDH activity increased in sugar beet after exposure to salt stress . Even though PGDH was down-regulated in ZD808, the response to salt stress-related PGDH in ZD808 may be greater than the responses of its parents. GST is a multifunctional enzyme that participates in reducing the formation of reactive oxygen species and cell toxicity [44, 45]. GST6 is up-regulated under salt-stress conditions . GST6 was classified as a ‘−’ NAP in ZD808, and a significant difference exist in ZD808 and its midparent value in PRM verification. Thus, the response to salt stress-related GST6 may be high in ZD808. MSBP1 inhibits cell elongation  and negatively regulates brassinosteroid signaling . MSBP1 expression is also closely related to salinity tolerance . Like many stress-related NAPs in ZD909, this protein was down-regulated, which was verified by PRM.
Eight proteins related to PS, fructose-bisphosphate aldolase 7 (FBA7), ferredoxin-NADP reductase (FNR), photosystem II subunit PsbS1 (PsbS1), post-illumination chlorophyll fluorescence increase (PIFI), chlorophyll a-b binding proteins (CAB1 and CAB2), photosystem I P700 chlorophyll a apoprotein A2 (psaB), photosynthetic NDH subunit of lumenal location 1 (PNSL1), were also subjected to PRM verification (Fig. 6). Except PIFI and PNSL1, the other six proteins related to PS were verified by PRM. These proteins had important roles in PS. PIFI, CAB1 and 2, psaB, and PNSL1 were identified in ZD909 as ‘+’ NAPs, and most of them were classified as ‘+‘or ‘++’ in the PRM verification. PIFI is essential for NDH-mediated nonphotochemical reduction of the plastoquinone pool in chlororespiratory electron transport . CABs are a class of important proteins in thylakoid membranes of higher plants, and take part in both photosystems I and II. In addition, CABs may be involved in abiotic stress responses . It was also found that CABs were highly expressed in the F1 hybrid of Chinese cabbage at two days after sowing, means that early developmental events in the germinating seedling of the hybrid may be important for later developmental vigor and yield advantage . PsaB a primary electron donor of photosystem I. PNSL1 is a lumen subcomplex L subunit of chloroplast NDH. Unlike in ZD909, most proteins related to PS identified by label-free quantitative proteomics were down-regulated, and FBA7, FNR, and PsbS1 classifications were verified by the PRM. FBAs participate in carbon fixation, glycolysis, and Calvin cycle pathways. FNR catalyzes reversible electron transport between ferredoxin and NADP. PsbS1 is a subunit of photosystem II. These roles indicated that PS was decreased in ZD808 compared with its parents. Based on the label-free quantitative proteomics and PRM analyses of NAPs, the heterosis of maize hybrids ZD808 and ZD909 may be mainly related to PS and/or stress resistance.
In addition to the above-mentioned proteins which related to stress responses and photosynthesis, a protein expressed one of reported genes associated with heterosis, ERBB-3 binding protein 1 (EBP1), was up-regulated in ZD909 (Table S2). EBP1 regulates organ size through cell growth and proliferation . And more, some DAPs/NAPs in ZD808 and ZD909 are not yet functionally annotated (Table S1, S2). And even though analysis of DAPs and NAPs in the germinating seedling of hybrids ZD808 and ZD909 has given us some useful information in understanding the reason of valuable agronomic traits, proteomics research throughout the growth stage is still lacking. All these results indicate that the mechanism of maize heterosis needs further study.
In summary, nearly 200 NAPs were identified from seedling leaves of maize hybrids ZD808 and ZD909, and most of them were enriched in stress responses, PS, amino acid metabolism and protein pathways. Pathways related to stress response, proteins, and amino acid metabolism were significantly enhanced in hybrid ZD808 but impeded in ZD909. However, pathways related to PS were significantly enhanced in hybrid ZD909. The results show that the biological basis of heterosis in the two maize hybrids is different. In comparison with their parents, the excellent agronomic traits of hybrid ZD808 was correlated with the high expression levels of some proteins related to stress responses and metabolic functions, while those of ZD909 was correlated with photosynthesis. Our proteomics results support previous physiological and morphological research and have provided useful information in understanding the reason of valuable agronomic traits.
Plant materials and growth conditions
The seeds of maize hybrid ZD808 (NG5 × CL11), along with its female parent CL11 (from 78,599) and male parent NG5 (from 95,236 × 95,167, Iowa State University), and those of the hybrid ZD909 (Z58 × HD568), along with its female parent Z58 and male parent HD568 (from DKC79 × Chang7–2) were acquired from professor Changling Huang’s lab (Center for crop genetics and breeding, Institute of crop science, Chinese Academy of Agricultural Sciences) in this study and some yield traits data are shown in Table 1. The seeds of each variety were washed separately with distilled water, soaked in 2% NaClO for 20 min, and washed again with distilled water. Then, ~ 200 seeds were germinated and cultured in sand for 15 days with a 16 h light at 28 °C, and 8 h dark at 24 °C photocycle and 70% relative humidity. Then, the seedlings with uniform growth were transplanted into new sand supplemented with 1/2 nutrient solution. Finally, the second leaves were collected at three-leaf stage, frozen in liquid nitrogen, and then stored at − 80 °C for the proteome analysis.
Protein extraction and digestion
Proteins were extracted according to Zhu et al. . Maize leaves of each sample (100 mg, three biological repeats and three technical repetitions) were used for protein extraction. The extracted proteins were digested by 50 μL trypsin and incubated overnight at 37 °C, centrifuged at 13,000×g for 30 min and stored at − 20 °C until further use.
Mass spectrometric analysis
Polypeptide fragments were separated and analyzed using an Easy nLC 1000 coupled to a Q-Exactive Plus (Thermo Fisher Scientific). Parameter setting of mass spectrometer were the same to Zhu et al. . The experiments with each sample were repeated at least three times, and then 54 mass spectrometry raw files with high reliability were obtained (Supplementary data: 54 mass spectrometry raw files, Partner repository with the dataset identifier PXD015819.).
Protein identification and LFQ
The label-free quantitative (LFQ) of proteins was carried out using MaxQuant software (version 1.6.0, https://www.coxdocs.org) and based on the Uniprot database (Maize, August 28, 2018). The mass tolerance of the precursor was set as 20 ppm, 0.02 Da fragment ion tolerance; up to two missed cleavages, carbamidomethyl cysteine as a fixed modification, N-terminal acetylation and methionine oxidation of proteins were defined as variable modifications. The cutoff global false discovery rate (FDR) for peptide and protein identification was set at 0.01.
The LFQ data were preprocessed using a standardized method before use in screening related proteins. A Pearson’s correlation analysis was used to evaluate the quality of the original data (SPSS 22.0, https://www.ibm.com/analytics). A Venn analysis (https://bioinfogp.cnb.csic.es/tools/venny/) was used to group the DAPs and NAPs between hybrids and their parents. Based on previous literature, the selection criterion for the DAPs was fold change > 2, and the criterion for the NAPs was fold change > 1.5 and p < 0.05. Principal component analysis was visualized using R package ggplot2_3.3.0 and ggfortify_0.4.9. Heatmap was visualized using R package pheatmap_1.0.12, and scaled row using Z-score. Finally, pathway enrichment analyses of DAPs and NAPs were performed using MapMan (https://MapMan.gabipd.org/).
Parallel reaction monitoring (PRM)
A PRM analysis were performed on an Easy nLC 1000 coupled to a Q-Exactive Plus (Thermo Fisher Scientific). Leaves were collected at three-leaf stage. Proteins were extracted and digested, and then used for PRM analysis. The mass spectrometry parameters included an isolation width of ± 1 Thomson/window of 2 Thomson for target precursor ions, and precursor ions were fragmented through higher energy collisional dissociation with normalized collision energies of 27 eV. In all the experiments, a full mass spectrum at 70,000 resolution relatives to m/z 200 was followed by up to 30 PRM scans at a 17,500 resolution. The other separation conditions and mass spectrometry parameters were as described above. Finally, raw data obtained were analyzed using Proteome Discoverer 2.1 (Thermo Electron, Germany) and Skyline 4.1 software (downloaded from ftp://ftp.ensembl.org/pub/mnt2/).
Availability of data and materials
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange. org) via the iProX partner repository with the dataset identifier PXD015819 (http://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD015819).
differentially accumulated protein
automatic gain control
parallel reaction monitoring
Liquid chromatography-tandem mass spectrometry
two-dimensional gel electrophoresis
Matrix-Assisted Laser Desorption/ Ionization Time of Flight Mass Spectrometry
mitogen-activated protein kinase
chaperone protein ClpB3
D- 3-phosphoglycerate dehydrogenase
glutathione S-transferase 6
membrane steroid-binding protein 1
Photosystem II subunit PsbS1
post-illumination chlorophyll fluorescence increase
chlorophyll a-b binding protein1
Photosystem I P700 chlorophyll a apoprotein A2
photosynthetic NDH subunit of lumenal location 1
Fujimoto R, Uezono K, Ishikura S, Osabe K, Peacock WJ, Dennis ES. Recent research on the mechanism of heterosis is important for crop and vegetable breeding systems. Breed Sci. 2018;68(2):145–58.
Kempe K, Rubtsova M, Gils M. Split-gene system for hybrid wheat seed production. Proc Natl Acad Sci U S A. 2014;111(25):9097–102.
Schnable PS, Springer NM. Progress toward understanding heterosis in crop plants. Annu Rev Plant Biol. 2013;64:71–88.
Marianyela Petrizzelli DdV, Christine Dillmann∗. Decoupling the Variances of Heterosis and Inbreeding Effects Is Evidenced in Yeast's Life-History and Proteomic Traits. Genetics. 2018;211(2).
Lippman ZB, Zamir D. Heterosis: revisiting the magic. Trends Genet. 2007;23(2):60–6.
Luo LJ, Li ZK, Mei HW, Shu QY, Tabien R, Zhong DB, et al. Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. II. Grain yield components. Genetics. 2001;158(4):1755–71.
Liu J, Li M, Zhang Q, Wei X, Huang X. Exploring the molecular basis of heterosis for plant breeding. J Integr Plant Biol. 2020;62(3):287–98.
Ma Q, Hedden P, Zhang Q. Heterosis in rice seedlings: its relationship to gibberellin content and expression of gibberellin metabolism and signaling genes. Plant Physiol. 2011;156(4):1905–20.
Wang T, Sui Z, Liu X, Li Y, Li H, Xing J, et al. Ectopic expression of a maize hybrid up-regulated gene, ErbB-3 binding protein 1 (ZmEBP1), increases organ size by promoting cell proliferation in Arabidopsis. Plant science : an international journal of experimental plant biology. 2016;243:23–34.
Liu D, Li S, Wang L, Li Q, Cui Y, Dai X, et al. Cloning and expression analysis of SiCDPK4, a gene related to Heterosis in foxtail millet [(Setaria italica (L.) P. Beauv.)]. J Plant Growth Regul. 2018;38(2):513–22.
Hushuai N. Cloning and functional analysis of genes related to maize heterosis Changchun: Jilin University; 2015.
Nguyen-Phan TC, Fry SC. Functional and chemical characterization of XAF: a heat-stable plant polymer that activates xyloglucan endotransglucosylase/hydrolase (XTH). Ann Bot. 2019;124(1):131–47.
Guo M, Rupe MA, Dieter JA, Zou JJ, Spielbauer D, Duncan KE, et al. Cell number Regulator1 affects plant and organ size in maize: implications for crop yield enhancement and Heterosis. Plant Cell. 2010;22(4):1057–73.
Birchler JA, Auger DL, Riddle NC. In search of the molecular basis of heterosis. Plant Cell. 2003;15(10):2236–9.
Marcon C, Schutzenmeister A, Schutz W, Madlung J, Piepho HP, Hochholdinger F. Nonadditive protein accumulation patterns in maize (Zea mays L.) hybrids during embryo development. J Proteome Res. 2010;9(12):6511–22.
Rockenbach MF, Correa CCG, Heringer AS, Freitas ILJ, Santa-Catarina C, do Amaral-Junior AT, et al. Differentially abundant proteins associated with heterosis in the primary roots of popcorn. PloS one. 2018;13(5):e0197114.
Vale EM, Reis RS, Santa-Catarina R, Pereira MG, Santa-Catarina C, Silveira V. Comparative proteomic analysis of the heterosis phenomenon in papaya roots. Sci Hortic-Amsterdam. 2016;209:178–86.
Klopfenstein TJ, Erickson GE, Berger LL. Maize is a critically important source of food, feed, energy and forage in the USA. Field Crop Res. 2013;153:5–11.
Charles W. Stuber SEL, David W. Wolff: Tim Helentjarisn, Eric S. Lander. Identification of Genetic Factors Contributing to Heterosis in a Hybrid From Two Elite Maize Inbred Lines Using Molecular Markers. Genetics. 1992;132(3).
Giraud H, Bauland C, Falque M, Madur D, Combes V, Jamin P, et al. Reciprocal genetics: identifying QTL for general and specific combining abilities in hybrids between multiparental populations from two maize (Zea mays L.) Heterotic groups. Genetics. 2017;207(3):1167–80.
Li H, Yang Q, Fan N, Zhang M, Zhai H, Ni Z, et al. Quantitative trait locus analysis of heterosis for plant height and ear height in an elite maize hybrid zhengdan 958 by design III. BMC Genet. 2017;18(1):36.
Jahnke S, Sarholz B, Thiemann A, Kuhr V, Gutierrez-Marcos JF, Geiger HH, et al. Heterosis in early seed development: a comparative study of F1 embryo and endosperm tissues 6 days after fertilization. TAG Theoretical and applied genetics Theoretische und angewandte Genetik. 2010;120(2):389–400.
Hu X, Wang H, Diao X, Liu Z, Li K, Wu Y, et al. Transcriptome profiling and comparison of maize ear heterosis during the spikelet and floret differentiation stages. BMC Genomics. 2016;17(1):959.
Lin Feng GM, Zhou Ling, Zhao Han. Genome-Wide Identification of Glyco-hydro-16 Family in Maize and Differentiation Analysis. Scientia Agricultura Sinica. 2016;49(11).
Hoecker N, Lamkemeyer T, Sarholz B, Paschold A, Fladerer C, Madlung J, et al. Analysis of nonadditive protein accumulation in young primary roots of a maize (Zea mays L.) F(1)-hybrid compared to its parental inbred lines. Proteomics. 2008;8(18):3882–94.
Hu X, Wang H, Li K, Wu Y, Liu Z, Huang C. Genome-wide proteomic profiling reveals the role of dominance protein expression in heterosis in immature maize ears. Sci Rep. 2017;7(1):16130.
Chen Y, Zhou Q, Tian R, Ma Z, Zhao X, Tang J, et al. Proteomic analysis reveals that auxin homeostasis influences the eighth internode length heterosis in maize (Zea mays). Sci Rep. 2018;8(1):7159.
Meena RK, Pullaiahgari D, Gudipalli P. Proteomic analysis of heterotic seed germination in maize using F1 hybrid DHM 117 and its parental inbreds. Turkish journal of biology = Turk biyoloji dergisi. 2018;42(4):345–63.
Bai XLQH, Liu M, Zhang ZY. Study on the relationship between drought resistance and physiological index of maize. Journal o fMaize Sciences. 2007;015(005):79–83.
Hou H-P, Ding Z-S, Ma W, Li C-F, Zhao M. Yield performance characteristics and regulation effects of plant density and sub-soiling tillage system for high yield population of summer maize. Acta Agron Sin. 2013;39(6):1069.
Swanson-Wagner RA, Jia Y, DeCook R, Borsuk LA, Nettleton D, Schnable PS. All possible modes of gene action are observed in a global comparison of gene expression in a maize F1 hybrid and its inbred parents. Proc Natl Acad Sci U S A. 2006;103(18):6805–10.
Wang X, Yan B, Shi M, Zhou W, Zekria D, Wang H, et al. Overexpression of a Brassica campestris HSP70 in tobacco confers enhanced tolerance to heat stress. Protoplasma. 2016;253(3):637–45.
Fei J, Wang YS, Zhou Q, Gu JD. Cloning and expression analysis of HSP70 gene from mangrove plant Kandelia obovata under cold stress. Ecotoxicology. 2015;24(7–8):1677–85.
Tang T, Yu A, Li P, Yang H, Liu G, Liu L. Sequence analysis of the Hsp70 family in moss and evaluation of their functions in abiotic stress responses. Sci Rep. 2016;6:33650.
Miller M, Song QX, Shi XL, Juenger TE, Chen ZJ. Natural variation in timing of stress-responsive gene expression predicts heterosis in intraspecific hybrids of Arabidopsis. Nat Commun. 2015;6.
Yang L, Li BS, Zheng XY, Li JG, Yang M, Dong XN, et al. Salicylic acid biosynthesis is enhanced and contributes to increased biotrophic pathogen resistance in Arabidopsis hybrids (vol 6, 7309, 2015). Nat Commun. 2015;6.
Rohde P, Hincha DK, Heyer AG. Heterosis in the freezing tolerance of crosses between two Arabidopsis thaliana accessions (Columbia-0 and C24) that show differences in non-acclimated and acclimated freezing tolerance. Plant J. 2004;38(5):790–9.
Pulido P, Llamas E, Llorente B, Ventura S, Wright LP, Rodriguez-Concepcion M. Specific Hsp100 chaperones determine the fate of the first enzyme of the Plastidial Isoprenoid pathway for either refolding or degradation by the stromal Clp protease in Arabidopsis. PLoS Genet. 2016;12(1):e1005824.
Llamas E, Pulido P, Rodriguez-Concepcion M. Interference with plastome gene expression and Clp protease activity in Arabidopsis triggers a chloroplast unfolded protein response to restore protein homeostasis. PLoS Genet. 2017;13(9):e1007022.
Perello C, Llamas E, Burlat V, Ortiz-Alcaide M, Phillips MA, Pulido P, et al. Differential Subplastidial localization and turnover of enzymes involved in Isoprenoid biosynthesis in chloroplasts. PLoS One. 2016;11(2):e0150539.
Grant GA. D-3-Phosphoglycerate dehydrogenase. Front Mol Biosci. 2018;5:110.
Okamura E, Hirai MY. Novel regulatory mechanism of serine biosynthesis associated with 3-phosphoglycerate dehydrogenase in Arabidopsis thaliana. Sci Rep. 2017;7(1):3533.
Kito K, Tsutsumi K, Rai V, Theerawitaya C, Cha-Um S, Yamada-Kato N, et al. Isolation and functional characterization of 3-phosphoglycerate dehydrogenase involved in salt responses in sugar beet. Protoplasma. 2017;254(6):2305–13.
Goto S, Kawakatsu M, Izumi S, Urata Y, Kageyama K, Ihara Y, et al. Glutathione S-transferase pi localizes in mitochondria and protects against oxidative stress. Free Radic Biol Med. 2009;46(10):1392–403.
Chasseaud LF. The Role of Glutathione and Glutathione S-Transferases in the Metabolism of Chemical Carcinogens and Other Electrophilic Agents. 1979;29:175–274.
Shinde H, Dudhate A, Tsugama D, Gupta SK, Liu S, Takano T. Pearl millet stress-responsive NAC transcription factor PgNAC21 enhances salinity stress tolerance in Arabidopsis. Plant physiology and biochemistry : PPB. 2019;135:546–53.
Yang XH, Xu ZH, Xue HW. Arabidopsis membrane steroid binding protein 1 is involved in inhibition of cell elongation. Plant Cell. 2005;17(1):116–31.
Song L, Shi QM, Yang XH, Xu ZH, Xue HW. Membrane steroid-binding protein 1 (MSBP1) negatively regulates brassinosteroid signaling by enhancing the endocytosis of BAK1. Cell Res. 2009;19(7):864–76.
Witzel K, Matros A, Moller ALB, Ramireddy E, Finnie C, Peukert M, et al. Plasma membrane proteome analysis identifies a role of barley membrane steroid binding protein in root architecture response to salinity. Plant Cell Environ. 2018;41(6):1311–30.
Wang D, Portis AR Jr. A novel nucleus-encoded chloroplast protein, PIFI, is involved in NAD(P)H dehydrogenase complex-mediated chlororespiratory electron transport in Arabidopsis. Plant Physiol. 2007;144(4):1742–52.
Silva J KY, Sukweenadhi J, Rahimi S, Kwon WS, Yang DC. Molecular characterization of 5 chlorophyll a-b-binding protein genes from Panax ginseng Meyer and their expression analysis during abiotic stresses. Photosynthetica. 2016;54(3).
Saeki N, Kawanabe T, Ying H, Shimizu M, Kojima M, Abe H, et al. Molecular and cellular characteristics of hybrid vigour in a commercial hybrid of Chinese cabbage. BMC Plant Biol. 2016;16.
Zhu L, Wang D, Sun J, Mu Y, Pu W, Ma B, et al. Phenotypic and proteomic characteristics of sorghum (Sorghum bicolor) albino lethal mutant sbe6-a1. Plant physiology and biochemistry : PPB. 2019;139:400–10.
We thank International Science Editing (https://www.internationalscienceediting.com) for editing this manuscript.
This work was jointly supported by the “National Key Research and Development Plan” (2016YFD0101005), the National Research Program of China (No. 2009ZX08012-011B), the “Beijing Municipal Science and Technology Project” (D171100007717002), the “Fundamental Research Funds for Central Non-profit Scientific Institution” (Y2017PT25) and the “Agricultural Science and Technology Innovation Program” (Development and Application of Molecular Markers in Crops). The funder was not involved in the experimental design of the study, data collection, analysis and interpretation, and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
We declare that all authors have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Total ion chromatograms of ZD808.
Total ion chromatograms of ZD808 female parent.
Total ion chromatograms of ZD808 male parent.
Total ion chromatograms of ZD909.
Total ion chromatograms of ZD909 female parent.
Total ion chromatograms of ZD909 male parent.
Principal component analysis of six samples.
Skyline analysis results of 6 PRM verified proteins unique peptides in ZD808.
Skyline analysis results of 6 PRM verified proteins unique peptides in ZD909.
Quality detection of mass spectrometry of ZD808 and its parents.
Quality detection of mass spectrometry of ZD909 and its parents.
Details of non-additive proteins of ZD808.
Details of non-additive proteins of ZD909.
Details of PRM verification data.
Peptide information of PRM verification.
Details of 12 heterosis-related proteins.
About this article
Cite this article
Wang, D., Mu, Y., Hu, X. et al. Comparative proteomic analysis reveals that the Heterosis of two maize hybrids is related to enhancement of stress response and photosynthesis respectively. BMC Plant Biol 21, 34 (2021). https://doi.org/10.1186/s12870-020-02806-5