Global transcriptional profiling between inbred parents and hybrids provides comprehensive insights into ear-length heterosis of maize (Zea mays)
BMC Plant Biology volume 21, Article number: 118 (2021)
Maize (Zea mays) ear length, which is an important yield component, exhibits strong heterosis. Understanding the potential molecular mechanisms of ear-length heterosis is critical for efficient yield-related breeding.
Here, a joint netted pattern, including six parent-hybrid triplets, was designed on the basis of two maize lines harboring long (T121 line) and short (T126 line) ears. Global transcriptional profiling of young ears (containing meristem) was performed. Multiple comparative analyses revealed that 874 differentially expressed genes are mainly responsible for the ear-length variation between T121 and T126 lines. Among them, four key genes, Zm00001d049958, Zm00001d027359, Zm00001d048502 and Zm00001d052138, were identified as being related to meristem development, which corroborated their roles in the superior additive genetic effects on ear length in T121 line. Non-additive expression patterns were used to identify candidate genes related to ear-length heterosis. A non-additively expressed gene (Zm00001d050649) was associated with the timing of meristematic phase transition and was determined to be the homolog of tomato SELF PRUNING, which assists SINGLE FLOWER TRUSS in driving yield-related heterosis, indicating that Zm00001d050649 is a potential contributor to drive heterotic effect on ear length.
Our results suggest that inbred parents provide genetic and heterotic effects on the ear lengths of their corresponding F1 hybrids through two independent pathways. These findings provide comprehensive insights into the transcriptional regulation of ear length and improve the understanding of ear-length heterosis in maize.
Heterosis is a phenomenon in which F1 hybrid progenies exhibit superior performances compared with those of their parents [1,2,3]. It is used for hybrid crop breeding, which has greatly increased the productivity of many crops worldwide [4, 5]. The successful exploitation of heterosis has also led researchers to determine its essential features. From the formation of two hypotheses (dominance and over-dominance) [6, 7] to the identification of genetic components [8,9,10], as well as comprehensive analyses of genomes, transcriptomes and metabolomes [11,12,13,14,15,16], tremendous efforts have been made to elucidate the mechanisms responsible for heterosis. However, the underlying molecular mechanisms remain poorly understood.
Transcriptional regulation plays roles in various aspects of plant growth and development. Variation in transcriptional regulation promotes phenotypic diversity in all species  and, thus, is a potential source of heterosis that could explain the differences between F1 hybrids and their parental lines. Many transcriptome analyses between hybrids and inbred lines have been carried out in both maize (Zea mays) and rice (Oryza sativa), and a great number of differentially expressed genes (DEGs) were found in F1 hybrids compared with their parents [4, 18,19,20]. Thus, heterosis was assumed to result from the global variation in gene expression between hybrids and inbred lines. However, several kinds of gene expression modes were observed in F1 hybrids: mid-parent (MP) (additivity), high and low parent (high and low parent dominance, respectively), above the high parent (over-dominance) and below the low parent (under-dominance) . These data revealed that in hybrids, some genes exhibit non-additive expression patterns (not the expected MP level), which suggested a potential association with heterosis [22,23,24]. These expression differences may be caused by allele-specific expression (ASE), which refers to the characteristic of preferentially expressing one parental allele in the hybrid owing to variations in regulatory sequences from the parental genome [25, 26]. Trans- and cis-regulation frequently exist in different parental lines, and they might be responsible for inducing ASE in hybrids [17, 27]. Consequently, analyzing transcriptional regulation is a valuable strategy for untangling the molecular basis of heterosis.
The measurement of heterosis involves a specific trait. Moreover, heterotic level is highly variable depending on the species, the cross parents and the trait(s) of interest . Maize is an important food crop worldwide, and it exhibits superior heterosis for a wide range of traits. In addition, its inbred lines have been classified into several “heterotic groups” on the basis of their heterotic level [28, 29]. Generally, crosses of parents within heterotic groups produce less heterosis than crosses of parents in different groups. This suggested that the inbred lines in each group may have specific exclusive properties that contribute to heterosis. Thus, between-group crosses are more likely to produce greater heterosis. Maize ear length is a representative trait with a superior heterotic level, and it contributes greatly to grain yield . A thorough knowledge of the transcriptional regulation of ear-length heterosis will aid in understanding the molecular basis of heterosis.
Maize ear length is predetermined, to some extent, by the activity of the ear primordium. As the ear primordium (meristem) differentiates, the visible young ear gradually elongates, revealing heterosis, and it is positively correlated with the final ear length . In this study, two specific maize lines, T121 and T126, with long and short ears, respectively, were identified. When crossed with other lines, T121 line, compared with T126 line, produced a series of hybrids with longer ears. However, equal ear-length heterosis was observed in their corresponding hybrids. With these two specific maize lines, we performed comprehensive transcriptional profiling of young ears using a joint netted pattern to determine the underlying cause of long ear in T121 line and to gain insights into ear-length heterosis.
The performance of ear length in different inbred-hybrid triplets
Ear-length heterosis in maize is a very striking phenomenon resulting from a cross of two distinct inbred lines. To explore ear-length heterosis, we selected two specific inbred lines (T121 and T126) with long and short ears, respectively. Additionally, two other inbred lines (PH4CV and PH6WC) were used to form a joint netted pattern (Fig. S1) that included six parent-hybrid triplets to adequately analyze ear-length heterosis. During maize ear differentiation, the young ear gradually elongates and becomes visible. Moreover, the elongation capability of the growth cone determines the final ear length to some extent. Here, we compared the morphologies of young ears of hybrids and their inbred parents at the 13-leaf stage when young ears were initially apparent. The young ears of T121 line were longer than those of T126 line. Moreover, the F1 hybrids generated by T121 crosses (T121 × PH4CV and T121 × PH6WC) had longer young ears than those generated by corresponding T126 crosses (T126 × PH4CV and T126 × PH6WC) (Table 1, Fig. 1a). In addition, the lengths of young ears from each inbred parent were less than those of their F1 hybrids (Table 1, Fig. 1a), indicating that ear-length heterosis had already emerged.
At the maturation stage, we measured the final ear lengths of all lines (inbred and hybrid). T121 line, a long-ear inbred line, had an ear length that reached 19.52 ± 1.18 cm and was much longer than that of the T126 (13.52 ± 0.83 cm) line (Table 1, Fig. 1b). The six F1 hybrids, T121 × PH4CV, T121 × PH6WC, T126 × PH4CV, T126 × PH6WC, T121 × T126 and PH6WC × PH4CV, exhibited MP heterosis for ear length (Table 1). Interestingly, the hybrids produced by T121 line (long ear), T121 × PH4CV and T121 × PH6WC, had longer ears than those produced by the short ear line T126 (Fig. 1b). However, for the MP heterosis, there was no significant difference between the other corresponding hybrids, such as T121 × PH4CV vs T126 × PH4CV (Table 1). These data indicate that T121 line makes a superior contribution to ear length than T126 line, but this is not a result of ear-length heterosis.
Transcriptome profiles of maize young ears among four inbred parents and six F1 hybrids
To understand the comprehensive transcriptional regulation of maize ear-length heterosis, young ears of four inbred parents and six F1 hybrids were used to perform RNA-sequencing (RNA-seq) analysis at the 13-leaf stage. In total, 332,543,189 raw reads were generated, ranging from 13.87 million to 21.76 million per library (Table 2). After filtering, 320,828,384 clean reads, accounting for 96.48% of the total, were maintained (Table 2). Based on the B73 maize reference genome (Version 4), the average unique mapping rate was 85.47%, with a range from 79.03 to 87.93% (Table 2). Moreover, two biological replicates were in close agreement (Fig. S2). Finally, 25,199 unique genes were identified in all lines (Table S1). The RNA-seq data is available for further analyses of transcriptional regulation.
Global transcriptome changes from inbred parents to their hybrids
Variation in gene expression is closely associated with phenotypic diversity. Thus, a series of transcriptional changes should occur from two inbred parents to one hybrid. For the T121–T126–T121 × T126 triplet, 64.97% of the genes in T121 × T126 hybrid kept their expression levels within the parental range, whereas the expression levels of the remaining genes (35.03%) were out of this rang (Table 3). This data indicated that the hybrids had the sufficient potential to surpass the two parents. Using a differential expression analysis, 5027 DEGs were identified between T121 and T126, 2547 DEGs were identified between T121 × T126 and T121, and 2431 DEGs were identified between T121 × T126 and T126 (Fig. 2a; Table S2). Thus, the number of DEGs between a hybrid and one parent (T121 or T126) was less than that between the two parents. Moreover, similar scenarios, including the ranges of the gene expression levels and the numbers of DEGs, were found in other parent-hybrid triplets (Tables 3, S2; Fig. 2a). Thus, some transcriptional regulatory mechanisms appeared to be universal and common in the production of hybrids from inbred parents.
In F1 hybrids, there are two gene expression patterns: additive and non-additive. We performed a non-additive expression analysis for the six parent-hybrid triplets. In total, 1375, 1349, 638, 566, 712 and 1257 non-additive genes were identified for T121–PH4CV–T121 × PH4CV, T121–PH6WC–T121 × PH6WC, T121–T126–T121 × T126, T126–PH4CV–T126 × PH4CV, T126–PH6WC–T126 × PH6WC and PH4CV–PH6WC–PH4CV × PH6WC, respectively (Fig. 2b; Table S3). A small proportion (< 6%) of non-additive patterns appeared in all F1 hybrids (Fig. 2b), indicating that most genes displayed an additive pattern in F1 hybrids. The additive pattern represents the expected MP level, whereas the non-additive pattern significantly deviates from the MP level. Thus, the non-additive genes in each triplet may have contributed to ear-length heterosis.
In hybrids, ASE frequently exists, increasing the plasticity of gene expression governed by diverse alleles from two parents, and this may be the reason that non-additive patterns appear in F1 hybrids. Thus, we analyzed genes having ASE in all the parent-hybrid triplets and then compared them with non-additively expressed genes identified in the same triplet. Quite a number of genes harboring ASE were detected in hybrids (Table S4). However, few of them were non-additively expressed (Fig. 2c–h). For instance, in the T121 × PH4CV hybrid, 1702 genes having ASE were identified, but only 71 exhibited non-additive expression patterns (Fig. 2c). These results indicated that in F1 hybrids, ASE might have a limited contribution to the production of non-additive expression-related variation.
The major genes responsible for the ear-length variation between T121 and T126 lines
T121 line produces longer ears than T126 line at the 13-leaf stage. A transcriptional level analysis of young ears revealed a large number of DEGs (5027) between T121 and T126 lines. However, it was difficult to determine the major genes responsible for the ear-length variation. Nevertheless, compared with T126 line, T121 had a longer ear and might pass this advantage to its F1 hybrid. When T121 and T126 were hybridized with the other parents (PH4CV and PH6WC), the former produced F1 hybrids with longer ears compared with the latter. Consequently, we performed a differential expression analysis between the corresponding F1 hybrids, T121 × PH4CV vs T126 × PH4CV and T121 × PH6WC vs T126 × PH6WC (Table S2). In total, 890 DEGs were found to overlap between the two groups (Fig. 3a). We compared these overlapped genes with the DEGs identified between lines T121 and T126. As expected, they shared many common genes (874) (Fig. 3b), which suggested that these genes take part in the regulation of ear elongation and are mainly responsible for the ear-length variation between T121 and T126 lines.
Gene ontology (GO) enrichment analysis was performed to identify some major terms related to ear length, as well as the key genes implicated in ear-length heterosis. A total of 1672 GO terms were enriched for these genes in biological process (Table S5). Furthermore, the top 10 GO terms were investigated, and they revealed several terms related to development, such as GO:0048582 (regulation of post-embryonic development) and GO:0048831 (regulation of shoot system development) (Fig. 3c; Table S5). Among these terms, four genes (Fig. 3d), Zm00001d027359 (FUSCA homolog, FUS6), Zm00001d048502 (COP9 signalosome complex subunit 1, CNS1), Zm00001d052138 (E3 ubiquitin-protein ligase, COP1) and Zm00001d049958 (WD40 repeat domain family protein, CYP71), were found to also belong to GO:0048507 (meristem development), and they may make major contributions to ear-length variation.
Non-additively expressed genes contributing to ear-length heterosis
Non-additively expressed genes may be potential sources of heterosis . To identify promising potential genes that contribute to ear-length heterosis derived from T121 (or T126) line, we made multiple comparisons of non-additively expressed genes in these parent-hybrid triplets. For T121 line, 47 non-additively expressed genes overlapped among hybrids produced by T121 × PH4CV, T121 × PH6WC and T121 × T126 (Fig. 4a). Whereas, for 126 line, 50 common non-additively expressed genes were identified among hybrids produced by T126 × PH4CV, T126 × PH6WC and T121 × T126 (Fig. 4b). These genes should be involved in ear-length heterosis, because the ear lengths of these F1 hybrids all surpassed the MP values. Moreover, 19 genes were shared (Fig. 4c, d), and these genes displayed non-additive expression patterns in all hybrids, suggesting that they had a potential to contribute to ear-length heterosis. GO enrichment analysis revealed that the top 10 GO terms for T121 and T126 lines were highly similar (Fig. 4e, f; Tables S6, 7), suggesting that there are some common components of the mechanism underlying ear-length heterosis. Among the common GO terms, GO:0048506 (regulation of timing of meristematic phase transition) and GO:0048510 (regulation of timing of transition from vegetative to reproductive phase) were associated with meristem, and a shared gene, Zm00001d050649 (ZCN2), may be responsible for the ear-length heterosis.
Validation of candidate gene expression by quantitative real-time PCR
The application of RNA-seq technology has greatly enhanced the global understanding of transcriptional regulatory networks. To verify the accuracy of RNA-seq analysis, we performed a quantitative real-time PCR (qRT-PCR) analysis of five candidate genes, including four DEGs having additive genetic effects on ear length, Zm00001d027359, Zm00001d048502, Zm00001d052138 and Zm00001d049958, and one non-additively expressed gene having heterotic effects on ear length, Zm00001d050649. Primers were designed to specifically amplify each of the five genes (Table S8). These primers were used to conduct qRT-PCR on three biological replications of RNA from re-prepared samples. All the assayed genes showed expression patterns similar to those determined by RNA-seq (Fig. 5), verifying the reliability of our RNA-seq analysis.
Two independent pathways of genetic and heterotic effects on ear length in maize
Three major hypotheses, dominance [6, 32], over-dominance [7, 33] and epistasis [34, 35], have served as the foundation for exploring the genetic and molecular causes of heterosis. The core notions are complementation within alleles, interactions within alleles and interactions between alleles, respectively . All of them highlight the potential contributions of alleles (preferably considered as genetic loci) from two parental inbred lines to their F1 hybrid. Indeed, conventional genetic loci refer to quantitative trait loci (QTL) for per se traits that control the additive genetic effect for trait performance, while heterotic effects are determined by special genetic loci, defined as QTL for trait heterosis, based on MP heterosis . Using recombinant inbred lines and immortalized F2 populations, some researchers have identified large numbers of QTL for per se traits and QTL for trait heterosis, respectively. Interestingly, very few of the two kinds overlaped, indicating that two independent pathways are responsible for their respective contributions to trait performance [10, 37]. However, quite a few overlapped loci were revealed in another study . More information is needed to determine the relationships between genomics of a per se trait and those involved in heterosis of the same trait.
Gene expression is a complex process involving a series of transcriptional regulations that affect an individual’s phenotype . Transcriptional regulation plays a role in explaining the molecular mechanisms of heterosis that benefit F1 hybrid individuals with a superior trait performance compared with the parental inbred individuals [19, 38]. Transcriptome profiles of two inbred lines and their hybrids have been determined to investigate variations in global gene expression. Consistent with many previous studies [39,40,41,42], there were more DEGs between the two parents than between each parent and the hybrid in all triplets, and only a few genes (< 6%) displayed non-additive patterns (Fig. 2b). Thus, the prevailing additive pattern appears to limit the difference between one parent and the hybrids and makes a limited contribution to heterosis in hybrids. Additionally, T121 line produced longer ears than T126 line, which indicated that T121 line had a superior additive genetic effect on ear length. However, their corresponding hybrids, such as T121 × PH4CV vs T126 × PH4CV, produced nearly equal levels of heterosis for ear length (Table 1). Thus, the additive genetic effect (parental variation) and heterotic effect appeared to be uncorrelated. Similarly, the genetic distance between two parents is a limited predictor of heterosis in their F1 hybrids .
Indeed, the ear lengths of T121 × T126 hybrids easily surpassed that of T121 line (extremely long ear). This was attributed to ear-length heterosis. Thus, the DEGs between T121 × T126 and T121 should be involved in ear-length heterosis, which increased the ear length compared with that of T121 line. However, some of these genes were also identified as being responsible for the ear-length variation between T121 and T126 lines, but their variation trends were opposite (up−/down-regulation) (Fig. S3). This indicated that these shared genes may not play roles in the ear-length heterosis. In this context, two inbred parents would provide the genetic and heterotic effects on the ear lengths of their F1 hybrids through two independent pathways. The superior performance of the hybrid over that of the better inbred parent benefits from the altered regulation of specific genes having non-additive expression patterns.
A possible pathway resulting in the non-additive expression pattern from inbred parents to their hybrids
From two inbred parents to one hybrid, it is tempting to infer that some new transcriptional regulations (non-additive expression patterns) contribute to heterosis. Determining the pathway responsible for the non-additive expression pattern may help to elucidate the mechanisms of heterosis.
Most transcriptional variation may be caused by sequence variation in regulatory regions of genes (cis-regulation) or by functional variations in a regulator (trans-regulation) . In hybrids, all genes consist of a pair of alleles derived from two parental inbred lines, respectively. Usually, the biased expression of alleles (ASE) takes place in some genes. If cis-regulation is present, then the allelic expression is expected to be additive in the hybrids, resulting in an additive expression pattern. In this study, we found that few genes having ASE present the non-additive expression pattern in all parent-hybrid triplets (Fig. 2c–h). Thus, this indicated that ASE is mainly caused by cis-regulation and makes a limited contribution to the non-additive expression pattern. Indeed, numerous studies [44,45,46,47] have revealed that cis-regulation plays major roles in the regulation of allelic expression in the hybrids of many species, and this may result in the production of large numbers of additively expressed genes.
Trans-regulation also occurs widely among regulatory networks, and it regulates the expression of many genes . A deficiency in trans-acting factors in one parent leads to the differential expression of their target genes between parents. However, because F1 hybrids have the same genetic background, the sharing of trans-acting factors facilitates the balanced expression of allelic genes, resulting in non-additive expression patterns. Here, we compared the two blocks: non-additively expressed genes in F1 hybrids and DEGs between two parents. Several common genes were discovered in each triplet (Fig. S4). For example, in the T121–PH4CV–T121 × PH4CV triplet, 395 genes overlapped between the two blocks (Fig. S4a). This revealed that trans-regulation may cause the variation in gene expression between two parents and produce the non-additive expression pattern in F1 hybrids. However, such genes would be responsible for the variation in ear length between the two parents (additive genetic effect) but would not contribute to ear-length heterosis.
Indeed, numerous isolated non-additively expressed genes were not differentially expressed between two inbred parents and presented non-ASE patterns in corresponding F1 hybrids (Figs. 2b, S4). One possible scenario is that cis- and trans-interactions occur. If the cis-regulated alleles also harbor variations in functions (trans-regulation), then non-additive expression patterns would be produced. An excellent example has been reported in tomato hybrids, in which yield was improved by fine-tuning the expression of a transcription factor (MADS-box) and its trans-effects on the target alleles . Alternatively, two or more trans-acting factors may combine to activate or suppress the expression of the target alleles in the F1 hybrids, leading to a non-additive expression pattern. For example, two maize transcription factors, B and Pl, interact to up-regulate the expression of genes A1, A2 and Bz1, which control anthocyanin production . An inbred line with a nonfunctional b or pl allele displays a green phenotype owing to the low, or absent, expression of genes A1, A2 and Bz1, but a hybrid with B/b Pl/pl alleles has high expression levels for genes A1, A2 and Bz1 and a red phenotype . Thus, once two parents are crossed, some specific interactions, rather than either cis- or trans-regulation, play major roles in generating the non-additive expression patterns found in hybrids, contributing to heterosis.
Key genes having superior additive genetic effects on ear length in T121 line
One inbred line often transfers its excellent characteristic to its progeny, including F1 hybrids, and this can be attributed to its superior additive genetic effect. In this study, T121 line produced an extremely long ear, and its ear length was far greater than that of T126 line. Likewise, its F1 hybrids exhibited longer ears compared with those of T126 line, indicating that T121 line had a superior additive genetic effect on ear length.
Maize ear development arises from axillary meristem, which can greatly influence ear size . Here, we identified four candidate genes related to meristem development, Zm00001d027359 (FUS6), Zm00001d048502 (CNS1), Zm00001d052138 (COP1) and Zm00001d049958 (CYP71). Interestingly, FUS6/CNS1 (FUS6 also called CNS1) and COP1 were originally found to act together in the photomorphogenesis of shoot apical meristem in Arabidopsis [51,52,53]. The COP9 signalosome (CSN) is a multifunctional protein complex composed of eight subunits (CSN1–8) in higher eukaryotes, such as plants, which regulates the activity levels of the cullin-RING ligase families of E3 ubiquitin ligase complexes . COP1 acts as an E3 ubiquitin ligase and translocates into the nucleus in a CSN-dependent manner, where it suppresses photomorphogenesis by promoting the degradation of positive of photomorphogenic regulators in the darkness [54, 55]. Upon light exposure, the nuclear COP1 is rapidly depleted, thus alleviating its suppression of photomorphogenic development . CSN and COP1 are involved in a range of plant growth and developmental processes [53, 55,56,57]. Therefore, the lower expression levels of Zm00001d027359 (FUS6/CNS1), Zm00001d048502 (CNS1) and Zm00001d052138 (COP1) in T121, T121 × PH4CV and T121 × PH6WC individuals (Figs. 3d, 5) accelerate ear growth and the rapid growth of plants after releasing photomorphogenetic suppression.
Additionally, Arabidopsis CYP71 is a unique immunophilin with a WD40 domain, and it interacts with histone H3 to regulate gene expression patterns that determine plant organogenesis [58, 59]. The CYP71 gene is preferentially expressed in meristem and other actively dividing tissues, and a loss of CYP71 function causes the arrest of apical meristem development [58, 59]. Similarly, higher expression levels of Zm00001d049958 (CYP71) in T121 line and its progeny (Figs. 3d and 5) was conducive to ear growth. In this context, a possible scenario is that the expression of Zm00001d049958 (CYP71) in maize axillary meristem suppresses the expression of Zm00001d027359 (FUS6/CNS1), Zm00001d048502 (CNS1) and Zm00001d052138 (COP1) genes and then promotes ear elongation. Overall, these four genes may be responsible for the additive genetic effects on ear length in T121 line.
The potential contributors to heterosis for the ear length
The performance of F1 hybrids mainly benefits from two aspects: the genetic and heterotic effects derived from two parents. Heterosis is specific to different traits and may be attributed to specific loci for a particular trait . Independent of the loci, drastic transcriptional variations in key genes must take place in hybrids. Genes harboring non-additive expression patterns have been studied owing to their huge contributing potential to heterosis [31, 60, 61].
Maize ear length is an important agronomic trait that often exhibits super heterosis . In this study, we analyzed global transcriptomes of young maize ears from six parent-hybrid triplets derived from four inbred lines. The Zm00001d050649 (ZCN2) gene displayed a non-additive expression pattern in all the triplets and belonged to GO:0048506 (regulation of timing of meristematic phase transition), implying its contribution to ear-length heterosis. Maize ZCN2 is a member of the TERMINAL FLOWER1 (TFL1)-like gene family, which is highly conserved in plants and is thought to function in the maintenance of meristem indeterminacy . In Arabidopsis, TFL1 and FLOWERING LOCUS T (FT) are two antagonistic integrators of the floral transition pathways that function in repressing and promoting flowering, respectively [63,64,65]. The tomato SINGLE FLOWER TRUSS gene, an ortholog of Arabidopsis FT, drives the heterosis for yield in an over-dominant pattern . These heterotic effects depend on the genetic background having a mutation in SELF PRUNING (SP), an ortholog of Arabidopsis TFL1. If plants carry a functional SP gene, then heterosis is eliminated [66, 67]. This suggested that the sp gene is a required contributor that drives heterosis for yield in tomato. Like the tomato sp gene, lack of expression of the homologous maize ZCN2 gene in hybrids (Fig. 5) contributes to ear-length heterosis.
In this study, multiple comparative analyses of the transcriptional profiles of six parent-hybrid triplets revealed that the genetic and heterotic effects on ear length in maize contribute to the performance of F1 hybrids through two independent pathways. Four key genes, Zm00001d049958 (CYP71), Zm00001d027359 (FUS6/CNS1), Zm00001d048502 (CNS1) and Zm00001d052138 (COP1), were identified as being responsible for the superior additive genetic effects on ear length in T121 line. Cis- and trans-regulatory interactions mainly caused the emergence of non-additive expression patterns in F1 hybrids, providing the potential to drive ear-length heterosis. The lack of expression of a non-additively expressed gene, Zm00001d050649 (ZCN2) was identified as potentially contributing to ear-length heterosis just as its homologous tomato SP gene contributes to yield heterosis. This will lead to investigations of the mechanism behind the silencing of Zm00001d050649 (ZCN2) in F1 hybrids, which will help further elucidate the mechanisms of heterosis. The present work provides insights into the transcriptional regulation of the maize ear-length characteristic from two parents to one hybrid. These findings improve our understanding of ear-length heterosis in maize.
Materials and methods
Experimental design and plant growth
Four maize inbred lines were selected in this study. T121 and T126 are inbred lines having long and short ears, respectively, which were derived from our breeding lines. The other two inbred lines, PH4CV and PH6WC, are the two parents of the excellent hybrid ‘Xianyu 335’. Then, inbreds were crosses to each other following a half-diallel (without reciprocals) design that resulted in a joint netted pattern comprising six F1 hybrids (Fig. S1).
The four parental inbred lines and six F1 hybrids were planted in a specially designed plot. The plot consisted of 4-m-long rows separated by 1-m spaces between each row, with 15 individuals planted per row. The five-row interval planting was performed for each line with two replicates. In addition, inbred and hybrid lines were separated. The 5th and 10th leaves of uniform individuals were labelled in the field. These experiments were carried out at the Scientific & Educational Park of Henan Agricultural University, Yuanyang, China.
Phenotypic heterosis for ear length
In maize, an axillary meristem forms at each stalk node beginning at the base of the stalk and continuing toward the top except for the upper six to eight nodes of the plant. The maize axillary meristem initiates ear development and only the upper one or two ear shoots ultimately become the harvestable (final) ears. The uppermost (final) ear is normally located at the 12th to 14th stalk node, corresponding to the 12th to 14th leaf . In practice, at the 13-leaf stage, young ears approximately 2–5 mm in length were initially visible and easily segregated. Young ears of 10 individuals per line were fixed in FAA composed of 5% formaldehyde (40% v/v), 5% acetate and 90% alcohol (75% v/v). Their morphologies were observed using a stereomicroscope. Then, their lengths were determined and used as phenotypic values. Phenotypic heterosis was evaluated using MP heterosis, which was calculated by the following equation: (F1 − MP)/MP × 100%, where F1 represents the phenotypic value of the F1 hybrid and MP represents the average of the two inbred parents. Moreover, at the maturation stage, mature ears of 10 individuals per line were also harvested and their ear lengths measured.
Sample preparation, RNA isolation and RNA sequencing
At the 13-leaf stage, young ears were isolated and prepared for RNA-seq. A total of 20 individuals per line from each replicate were mixed for the RNA isolation. Total RNA was extracted using TRIzol reagent. RNA quantity and purity were determined using a Nanodrop 2000 and capillary electrophoresis. All samples with an RNA integrity number greater than 7 were considered of good quality. A total amount of 1 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA). In total, 20 RNA samples (10 varieties × 2 replicates) were supplied for deep sequencing using Illumina NovaSeq (150-bp paired-end) at BerryGenomics (Beijing, China), and the raw data included approximately 13–21 million reads per sample.
Analysis of RNA-Seq data: mapping and quantifying
When the sequencing was completed, Cutadapt 1.10 and in-house Perl scripts were used for quality control . Raw reads were filtered to remove adapters and low-quality bases, as well as reads less than 50 bp in length. Then, HISAT 2.0  was used to map clean reads to the B73 maize reference genome (Version 4), and StringTie 1.3  was used to assemble mapped reads. Finally, using StringTie 1.3 together with Ballgown , fragments per kilobase of exon per million (FPKM) mapped sequence reads values were calculated for each sample to estimate the level of gene expression. In addition, the correlation coefficient of the two biological replicates was calculated to evaluate the repeatability of the experiment. The averages of two replicate samples were regarded as the gene expression levels in each line.
The criteria (statistical significance) of p-value < 0.05 and abs (log2 (fold-change) > 1 were used to identify DEGs between two lines with the DESeq package (http://www.bioconductor.org/packages/). Between each two lines (inbred parents and F1 hybrids), differential expression analysis was performed. Moreover, the non-additive expression pattern was analyzed in each parent-hybrid triplet. For each triplet, the average expression levels of two parents were calculated as the MP value. Based on the above criteria, non-additively expressed genes were defined as having differential expression levels between those of the F1 hybrid and the MP value.
GO enrichment analysis was conducted to determine the essential functions of the DEGs (https://www.omicshare.com/tools). The top 10 GO terms were investigated to determine the major candidate genes. The threshold p-value < 0.05 were used for the analysis.
A specific filter was required for mapping reads. Using a customized Perl script, desired reads that were perfectly mapped to one parental sequence and had single nucleotide polymorphisms mapped to the other were retained. Then, the refiltered reads were assembled according to the previously reported criteria . In each triplet, the refiltered reads from F1 hybrid were divided into two sets: set 1, reads aligned against one parent, and set 2, reads aligned against the other parent, to distinguish parent-specific reads in the single nucleotide polymorphism calling step. The normalization of these read numbers was performed using the function estimateSize Factors from the DESeq package . For each gene, ASE was called if the reads of each set deviated significantly from 1:1 by simple random sampling, which was validated by 1000 permutations at a false discovery rate < 0.05.
Quantitative real-time PCR analysis
The same samples were re-prepared for a quantitative real-time PCR (qRT-PCR) analysis in an attempt to validate the expression patterns of key genes. The qRT-PCR was performed in a Bio-Rad CFX96 Real-Time PCR System with SYBR Green PCR Master Mix (Takara Bio). Three technical replicates were included in each plate for qRT-PCR. The Zm00001d013873 (ACTIN-2) gene was used as an internal standard to normalize gene expression, and the relative gene expression levels were measured using the 2–ΔΔCt method . Primers were designed online (https://www.ncbi.nlm.nih.gov/) for the candidate genes, and the primer information was provided in Table S8.
Availability of data and materials
The datasets generated and analysed during the current study are available in the NCBI Sequence Read Archive (SRA) database under Bioproject PRJNA682653 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA682653).
Differentially expressed genes
Quantitative real-time pcr
Quantitative trait loci
Flowering locus t
Fragments per kilobase of exon per million
Chen ZJ. Genomic and epigenetic insights into the molecular bases of heterosis. Nat Rev Genet. 2013;14(7):471–82.
Schnable PS, Springer NM. Progress toward understanding heterosis in crop plants. Annu Rev Plant Biol. 2013;64(1):71–88.
Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA. Heterosis. Plant Cell. 2010;22(7):2105.
Shahzad K, Zhang X, Guo L, Qi T, Bao L, Zhang M, Zhang B, Wang H, Tang H, Qiao X, et al. Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton. BMC Plant Biol. 2020;20(1):239.
Yu Y, Zhu M, Cui Y, Liu Y, Li Z, Jiang N, Xu Z, Xu Q, Sui G. Genome sequence and QTL analyses using backcross recombinant inbred lines (BILs) and BILF1 lines uncover multiple heterosis-related loci. Int J Mol Sci. 2020;21(3):780.
Bruce AB. The mendelian theory of heredity and the augmentation of vigor. Science. 1910;32(827):627–8.
East EM. Heterosis. Genetics. 1936;21(4):375.
Zhou G, Chen Y, Yao W, Zhang C, Xie W, Hua J, Xing Y, Xiao J, Zhang Q. Genetic composition of yield heterosis in an elite rice hybrid. Proc Natl Acad Sci U S A. 2012;109(39):15847–52.
Guo T, Yang N, Tong H, Pan Q, Yang X, Tang J, Wang J, Li J, Yan J. Genetic basis of grain yield heterosis in an "immortalized F2" maize population. Theor Appl Genet. 2014;127(10):2149–58.
Tang J, Yan J, Ma X, Teng W, Wu W, Dai J, Dhillon BS, Melchinger AE, Li J. Dissection of the genetic basis of heterosis in an elite maize hybrid by QTL mapping in an immortalized F2 population. Theor Appl Genet. 2010;120(2):333–40.
Huang X, Yang S, Gong J, Zhao Q, Feng Q, Zhan Q, Zhao Y, Li W, Cheng B, Xia J, et al. Genomic architecture of heterosis for yield traits in rice. Nature. 2016;537(7622):629–33.
Huang X, Yang S, Gong J, Zhao Y, Feng Q, Gong H, Li W, Zhan Q, Cheng B, Xia J, et al. Genomic analysis of hybrid rice varieties reveals numerous superior alleles that contribute to heterosis. Nat Commun. 2015;6:6258.
Yang M, Wang X, Ren D, Huang H, Xu M, He G, Deng XW. Genomic architecture of biomass heterosis in Arabidopsis. Proc Natl Acad Sci U S A. 2017;114(30):8101.
Li D, Huang Z, Song S, Xin Y, Mao D, Lv Q, Zhou M, Tian D, Tang M, Wu Q, et al. Integrated analysis of phenome, genome, and transcriptome of hybrid rice uncovered multiple heterosis-related loci for yield increase. Proc Natl Acad Sci U S A. 2016;113(41):E6026–e6035.
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.
Yi G, Shin H, Park HR, Park JE, Ahn JH, Lim S, Lee JG, Lee EJ, Huh JH. Revealing biomass heterosis in the allodiploid x Brassicoraphanus, a hybrid between Brassica rapa and Raphanus sativus, through integrated transcriptome and metabolites analysis. BMC Plant Biol. 2020;20(1):252.
Botet R, Keurentjes JJB. The role of transcriptional regulation in hybrid vigor. Front Plant Sci. 2020;11:410.
Stupar RM, Gardiner JM, Oldre AG, Haun WJ, Chandler VL, Springer NM. Gene expression analyses in maize inbreds and hybrids with varying levels of heterosis. BMC Plant Biol. 2008;8:33.
Zhiguo E, Huang S, Zhang Y, Ge L, Wang L. Genome-wide transcriptome profiles of rice hybrids and their parents. Int J Mol Sci. 2014;15(11):20833–45.
Hu X, Wang H, Diao X, Liu Z, Li K, Wu Y, Liang Q, Wang H, Huang C. Transcriptome profiling and comparison of maize ear heterosis during the spikelet and floret differentiation stages. BMC Genomics. 2016;17(1):959.
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.
Song GS, Zhai HL, Peng YG, Zhang L, Wei G, Chen XY, Xiao YG, Wang L, Chen YJ, Wu B, et al. Comparative transcriptional profiling and preliminary study on heterosis mechanism of super-hybrid rice. Mol Plant. 2010;3(6):1012–25.
Wei G, Tao Y, Liu G, Chen C, Luo R, Xia H, Gan Q, Zeng H, Lu Z, Han Y, et al. A transcriptomic analysis of superhybrid rice LYP9 and its parents. Proc Natl Acad Sci U S A. 2009;106(19):7695–701.
Bao J, Lee S, Chen C, Zhang X, Zhang Y, Liu S, Clark T, Wang J, Cao M, Yang H, et al. Serial analysis of gene expression study of a hybrid rice strain (LYP9) and its parental cultivars. Plant Physiol. 2005;138(3):1216–31.
Shao L, Xing F, Xu C, Zhang Q, Che J, Wang X, Song J, Li X, Xiao J, Chen L-L, et al. Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis. Proc Natl Acad Sci U S A. 2019;116(12):5653.
Gaur U, Li K, Mei S, Liu G. Research progress in allele-specific expression and its regulatory mechanisms. J Appl Genet. 2013;54(3):271–83.
Metzger BPH, Duveau F, Yuan DC, Tryban S, Yang B, Wittkopp PJ. Contrasting frequencies and effects of cis- and trans-regulatory mutations affecting gene expression. Mol Biol Evol. 2016;33(5):1131–46.
Fan X, Zhang Y, Yao W, Chen H, Tan J, Xu C, Han X, Luo L, Kang M. Classifying maize inbred lines into heterotic groups using a factorial mating design. Agron J. 2009;101(1):106–12.
Suwarno WB, Pixley KV, Palacios-Rojas N, Kaeppler SM, Babu R. Formation of heterotic groups and understanding genetic effects in a provitamin a biofortified maize breeding program. Crop Sci. 2014;54(1):14–24.
Wang X, Cao H, Zhang D, Li B, He Y, Li J, Wang S. Relationship between differential gene expression and heterosis during ear development in maize (Zea mays L.). J Genet Genomics. 2007;34(2):160–70.
Auger DL, Gray AD, Ream TS, Kato A, Coe EH, Birchler JA. Nonadditive gene expression in diploid and triploid hybrids of maize. Genetics. 2005;169(1):389–97.
Jones DF. Dominance of linked factors as a means of accounting for heterosis. Proc Natl Acad Sci U S A. 1917;3(4):310–2.
Shull GH. The composition of a field of maize. J Hered. 1908;os-4(1):296–301.
Minvielle F. Dominance is not necessary for heterosis: a two-locus model. Genet Res. 1987;49:245–7.
Schnell FW, Cockerham CC. Multiplicative vs. arbitrary gene action in heterosis. Genetics. 1992;131(2):461–9.
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.
Chen J, Zhou H, Xie W, Xia D, Gao G, Zhang Q, Wang G, Lian X, Xiao J, He Y. Genome-wide association analyses reveal the genetic basis of combining ability in rice. Plant Biotechnol J. 2019;17(11):2211–22.
Stupar RM, Springer NM. Cis-transcriptional variation in maize inbred lines B73 and Mo17 leads to additive expression patterns in the F1 hybrid. Genetics. 2006;173(4):2199.
Hoecker N, Keller B, Muthreich N, Chollet D, Descombes P, Piepho HP, Hochholdinger F. Comparison of maize (Zea mays L.) F1-hybrid and parental inbred line primary root transcriptomes suggests organ-specific patterns of nonadditive gene expression and conserved expression trends. Genetics. 2008;179(3):1275.
Jahnke S, Sarholz B, Thiemann A, Kühr V, Gutiérrez-Marcos JF, Geiger HH, Piepho HP, Scholten S. Heterosis in early seed development: a comparative study of F1 embryo and endosperm tissues 6 days after fertilization. Theor Appl Genet. 2009;120(2):389.
Li A, Fang MD, Song WQ, Chen CB, Qi LW, Wang CG. Gene expression profiles of two intraspecific Larix lines and their reciprocal hybrids. Mol Biol Rep. 2012;39(4):3773–84.
Meyer RC, Witucka-Wall H, Becher M, Blacha A, Boudichevskaia A, Dörmann P, Fiehn O, Friedel S, von Korff M, Lisec J, et al. Heterosis manifestation during early Arabidopsis seedling development is characterized by intermediate gene expression and enhanced metabolic activity in the hybrids. Plant J. 2012;71(4):669–83.
Springer NM, Stupar RM. Allele-specific expression patterns reveal biases and embryo-specific parent-of-origin effects in hybrid maize. Plant Cell. 2007;19(8):2391–402.
Doss S, Schadt EE, Drake TA, Lusis AJ. Cis-acting expression quantitative trait loci in mice. Genome Res. 2005;15(5):681–91.
Gibson G, Riley-Berger R, Harshman L, Kopp A, Vacha S, Nuzhdin S, Wayne M. Extensive sex-specific nonadditivity of gene expression in Drosophila melanogaster. Genetics. 2004;167(4):1791–9.
Springer NM, Stupar RM. Allelic variation and heterosis in maize: how do two halves make more than a whole? Genome Res. 2007;17(3):264–75.
Wittkopp PJ, Haerum BK, Clark AG. Evolutionary changes in cis and trans gene regulation. Nature. 2004;430(6995):85–8.
Soyk S, Lemmon ZH, Oved M, Fisher J, Liberatore KL, Park SJ, Goren A, Jiang K, Ramos A, van der Knaap E, et al. Bypassing negative epistasis on yield in tomato imposed by a domestication gene. Cell. 2017;169(6):1142–1155.e1112.
Dooner HK, Robbins TP, Jorgensen RA. Genetic and developmental control of anthocyanin biosynthesis. Annu Rev Genet. 1991;25:173–99.
Lejeune P, Bernier G. Effect of environment on the early steps of ear initiation in maize (Zea mays L.). Plant Cell Environ. 2006;19:217–24.
Cope GA, Deshaies RJ. COP9 signalosome: a multifunctional regulator of SCF and other cullin-based ubiquitin ligases. Cell. 2003;114(6):663–71.
Staub JM, Wei N, Deng XW. Evidence for FUS6 as a component of the nuclear-localized COP9 complex in Arabidopsis. Plant Cell. 1996;8(11):2047–56.
Kim JY, Song JT, Seo HS. COP1 regulates plant growth and development in response to light at the post-translational level. J Exp Bot. 2017;68(17):4737–48.
Osterlund MT, Ang LH, Deng XW. The role of COP1 in repression of Arabidopsis photomorphogenic development. Trends Cell Biol. 1999;9(3):113–8.
Qin N, Xu D, Li J, Deng XW. COP9 signalosome: discovery, conservation, activity, and function. J Integr Plant Biol. 2020;62(1):90–103.
Serino G, Deng XW. The COP9 signalosome: regulating plant development through the control of proteolysis. Annu Rev Plant Biol. 2003;54:165–82.
Tuller T, Diament A, Yahalom A, Zemach A, Atar S, Chamovitz DA. The COP9 signalosome influences the epigenetic landscape of Arabidopsis thaliana. Bioinformatics. 2019;35(16):2718–23.
Li H, Luan S. The cyclophilin AtCYP71 interacts with CAF-1 and LHP1 and functions in multiple chromatin remodeling processes. Mol Plant. 2011;4(4):748–58.
Li H, He Z, Lu GH, Lee SC, Alonso J, Ecker J, Luan S. A WD40 domain cyclophilin interacts with histone H3 and functions in gene repression and organogenesis in Arabidopsis. Plant Cell. 2007;19:2403–16.
Baldauf JA, Marcon C, Paschold A, Hochholdinger F. Nonsyntenic genes drive tissue-specific dynamics of differential, nonadditive, and allelic expression patterns in maize hybrids. Plant Physiol. 2016;171(2):1144–55.
Hoecker N, Lamkemeyer T, Sarholz B, Paschold A, Fladerer C, Madlung J, Wurster K, Stahl M, Piepho HP, Nordheim A, et al. Analysis of nonadditive protein accumulation in young primary roots of a maize (Zea mays L.) F1-hybrid compared to its parental inbred lines. Proteomics. 2008;8(18):3882–94.
Danilevskaya ON, Meng X, Ananiev EV. Concerted modification of flowering time and inflorescence architecture by ectopic expression of TFL1-like genes in maize. Plant Physiol. 2010;153(1):238–51.
Kardailsky I, Shukla VK, Ahn JH, Dagenais N, Christensen SK, Nguyen JT, Chory J, Harrison MJ, Weigel D. Activation tagging of the floral inducer FT. Science. 1999;286(5446):1962.
Kobayashi Y, Kaya H, Goto K, Iwabuchi M, Araki T. A pair of related genes with antagonistic roles in mediating flowering signals. Science. 1999;286(5446):1960–2.
Hanzawa Y, Money T, Bradley D. A single amino acid converts a repressor to an activator of flowering. Proc Natl Acad Sci U S A. 2005;102(21):7748.
Krieger U, Lippman ZB, Zamir D. The flowering gene SINGLE FLOWER TRUSS drives heterosis for yield in tomato. Nat Genet. 2010;42(5):459–63.
Shalit A, Rozman A, Goldshmidt A, Alvarez JP, Bowman JL, Eshed Y, Lifschitz E. The flowering hormone florigen functions as a general systemic regulator of growth and termination. Proc Natl Acad Sci U S A. 2009;106(20):8392.
Nielsen RL. Ear initiation & size determination in corn. Corny News Network. 2003; https://www.agry.purdue.edu/ext/corn/news/articles.03/earsize-0609.html. Accessed 9 June 2003.
Martin M. Cutasapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):200.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290–5.
Frazee AC, Pertea G, Jaffe AE, Langmead B, Salzberg SL, Leek JT. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat Biotechnol. 2015;33(3):243–6.
Liu H, Yang X, Liao X, Zuo T, Qin C, Cao S, Dong L, Zhou H, Zhang Y, Liu S, et al. Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development. Genomics. 2015;106(1):52–60.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
We are deeply indebted to the entire research team for their support during this research work. We also thank Lesley Benyon, PhD, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
This research work was supported by funds received from NSFC Henanunion fund project (U1604231), Zhengzhou major science and technology innovation project (188PCXZX803), NSFC Henanunion fund project (U1604113), Basic and Frontier Technology Research Program of Henan Province (162300410165), Henan Province Natural Fund (202300410215).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Schematic drawing of a joint netted pattern including four maize inbred parents and six F1 hybrids. Thereinto, two specific maize lines T121 and T126 displayed long ear and short ear, respectively
. Correlation analysis for each RNA-seq replicate
. Parental variation and greater than better-parental variation in gene expression for triplet T121–T126–T121 × T126. a Venn diagram comparison between the DEGs in T121 × T126 vs T121 and the key genes responsible for the ear length variation in T121 vs T126. b The expression profile of random 10 of those shared genes in a
. Venn diagram comparison between non-additively expressed genes in the F1 hybrids and DEGs between the two parents in the six triplets. a-f represent triplet T121–PH4CV, T121–PH6WC, T121–T126, T126–PH4CV, T126–PH6WC and PH4CV–PH6WC, respectively. Non- represents non-additive expression
. The mapped genes and their expression levels
. The differential expression analysis between different varieties
. The non-additive expression analysis for the six parent-hybrid triplets
. The ASE analysis for the six parent-hybrid triplets
. The GO enrichment analysis of the major genes responsible for the ear length variation between T121 and T126 lines
. The GO enrichment analysis of the non-additive expression genes for T121 line
. The GO enrichment analysis of the non-additive expression genes for T126 line
. The sequence of primers used for qRT-PCR
About this article
Cite this article
Zhang, X., Ma, C., Wang, X. et al. Global transcriptional profiling between inbred parents and hybrids provides comprehensive insights into ear-length heterosis of maize (Zea mays). BMC Plant Biol 21, 118 (2021). https://doi.org/10.1186/s12870-021-02890-1
- Ear length
- Transcriptional regulation
- Non-additive expression patterns