Dynamic Patterns of Gene Expression and Regulatory Variation in the Maize Seed Coat

Dynamic Patterns of Gene Expression and Regulatory Variation in the Maize Seed Coat | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Dynamic Patterns of Gene Expression and Regulatory Variation in the Maize Seed Coat Juan Li, Liangfa Wang, Jiong Wan, Kuntai Dang, Yuan Lin, Shujun Meng, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-1768077/v2 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 07 Feb, 2023 Read the published version in BMC Plant Biology → Version 2 posted 7 You are reading this latest preprint version Show more versions Abstract Background Seed size is an important factor contributing to maize yield, but its molecular mechanism remains unclear. The seed coat, which serves as one of the three components of the maize grain, determines seed size to a certain extent. The seed coat also shares the maternal genotype and is an ideal material for studying heterosis. Results In this study, the self-pollinated seeds of the maize hybrid Yudan888 and its parental lines were continuously collected from 0 day after pollination (DAP) to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP and 8 DAP are the best time points to study maize seed coat heterosis. Cytological observations indicated that maize seed coat heterosis might be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results showed that the nonadditive genes changed significantly between 3 DAP and 8 DAP. However, the number of genes expressed additively was not significantly different. Our findings suggest that seed coat heterosis in hybrid is the result of nonadditive expression caused by dynamic changes in genes at different time points during seed expansion and seed coat development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment indicated that genes related to DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation contributed significantly to hybrid seed coat heterosis. Conclusion Maize seed coat phenotyping allowed us to infer that 3 DAP and 8 DAP are important time points in the study of seed coat heterosis. Our findings provide evidence for genes involved in DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation in hybrid with high or low parental expression as major contributors to hybrid seed coat heterosis. Maize Seed coat Heterosis RNA-seq. Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Background Heterosis has been widely used in plant breeding[1]. Due to heterosis, hybrid F 1 generations are superior to both parents in terms of viability, plant height, yield, stress and disease resistance [2]. It is hypothesized that the difference in gene expression between a hybrid and its parents is one of the main causes of heterosis [3, 4]. Genetic models to explain the increased yield of hybrids consider the interaction of alleles at many loci generating altered expression levels and patterns [5, 6]. Whole-genome transcriptomes have shown that in developing rice leaves and panicles, differentially expressed genes (DEGs) between hybrids and their parents were overrepresented during energy metabolism and transport [7]. Changes in epigenetic states such as DNA methylation, small RNA production and histone modification have been found in hybrid genomes [8, 9]. Genetic mapping has been used to identify the detailed genetic loci contributing to hybrid performances and investigate their mechanisms, in a number of plant species including Arabidopsis, rice, maize, sorghum, and tomato [10-14]. Rice GAI homologues ectopically expressed in the F 1 generation of Arabidopsis thaliana caused significant plant height changes [10]. The major quantitative trait gene Ghd7 , which encodes a CCT domain protein in rice, was identified as a key factor for improving the yield and adaptability of the excellent hybrid Shanyou63 and other indica rice varieties [11]. The orthologue of the flowering locus SFT for tomato [12] and Hd3a for rice [13], has been shown to have single-gene overdominance in yield traits. However, there have been many studies on heterosis, and the genetic mechanism and molecular basis of heterosis remain unclear. Seed size is an important component of grain yield and a key trait for crop domestication [15]. Seed size is coordinately determined by the development of the triploid endosperm, diploid embryo, and seed coat [7, 16]. The seed coat which is developed from the integument and shares the maternal genotype protects the embryo and endosperm and grows coordinately with them[17]. It plays an important role in the seed filling process [18]. Therefore, the seed coat determines the yield of cereal crops to a certain extent. The development of the seed coat and endosperm takes precedence over that of the embryo. In the early stage of seed development, from double fertilization to 6 days after pollination (DAP), the growth of the seed coat is accompanied by enlargement of the endosperm. Within 72 hours after fertilization, the endosperm converts to the multinuclear stage [19]. At 3-6 DAP, the endosperm of maize is in the cellularization stage, during which the cytoskeleton recognizes the nuclei and forms the microtubule system between the nuclei. All nuclei are separated and cytoplasmic division is carried out at the same time to form normal cells [20]. Eight to 12 DAP is the endosperm cell differentiation stage, during which nutrients and storage materials accumulate alongside cell proliferation and enlargement [21]. The morphological changes of the seed coat occur only when the endosperm begins to develop, and the development of the embryo does not promote the growth and differentiation of the seed coat [22, 23]. Thus, the seed coat coincides with the enlargement of the seed in the early stage of development. In addition, the seed coat regulates the accumulation of phloem assimilates in seeds during seed development [24]. At the beginning of the filling stage at 12 DAP, the seed coat thickness is continuously compressed until the seed coat thickness decreases to only a few layers of cells. At this stage, the seed coat cell wall thickens, providing physical protection for future seed dormancy. The function of the seed coat during development may be to use its own starch for endosperm expansion and to provide more space for endosperm development. An impressive array of genomic and regulatory variations in kernel traits has been observed between maize inbred lines. However, how these variations affect gene expression in hybrids remains unknown. Since the integument does not participate in the fertilization process, the seed coat differentiates primarily from the integument is determined by the maternal plant and shares the maternal genotype[25]. Therefore, the genotype of the seed coat of the F 1 hybrid and its parents is still the same as that of the maternal plant after fertilization. Because seed size is a heterotic trait, the seed coat of F2 seeds obtained by self-pollination of F 1 plants shows heterozygosity, which is a good model for studying heterosis. However, the genetic basis of heterosis in the seed coat has remained largely uncharacterized. To this end, this study attempted to better define gene expression variations between seed coats of maize inbred lines, how these variations affect gene expression levels in F 1 hybrids, and the relative contributions of cis- and trans-variantsto seed coat heterosis. In this study, the self-pollinated kernels of the hybrid Yudan888 and its parental lines were continuously collected from 0 DAP to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP is an important time point for the expansion of the seed coat area in the hybrid, and 8 DAP is the key time point when the seed coat area of the hybrid significantly extends in the parents.Three DAP and 8 DAP are the best time points to study maize seed coat heterosis. According cytological observations at 8 DAP, we speculate that the heterosis of the seed coat may be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results at 3 DAP and 8 DAP indicated that the number of genes expressed nonadditively changed significantly in the hybrid. However, the number of genes expressed additively was not significantly different. Our findings suggest that seed coat heterosis is the result of nonadditive expression caused by dynamic changes in genes at different time points. These dynamic change suggests that a dose effect of gene expression results in heterosis. GO enrichment and KEGG analysis indicated that at 3 DAP, the early stage of seed coat development, genes related to DNA replication and cell cycle regulation contributed significantly to heterosis. At 8 DAP, the hybrid weakened the cell cycle by enhancing metabolite accumulation combined with increased kernel weight and stabilizing heterosis. We expect that this transcriptome comparison will provide a new step towards understanding the causative mechanisms of altered gene expression in hybrids and the molecular mechanisms of maize seed coat heterosis. Results The seed coat area of hybrid and its parents was significantly different at 8 DAP. The seed coat area and kernel weight from 0 to 15 DAP in the hybrid and parental lines were measured continuously (Fig.1A). At 0 DAP, the seed coat area of the hybrid was lower than that of both the mid-parental value (MPV) and the maternal line. At 6 DAP, the hybrid values began to exceed the MPV but were still lower than those of the maternal line. With the expansion of the seed coat, the hybrid value continued to be higher than the MPV at 7 DAP and slightly exceeded the value of the maternal line. At 8 DAP, there was a significant difference (P value = 0.034) between the hybrid and the MPV regarding seed coat area (Fig. 1B). From 9 to 15 DAP, the seed coat area of the hybrid was increasingly different from the MPV and paternal line values (Table S1). In conclusion, 8 DAP is the inflection point and the entry point of heterosis in the seed coat of maize. The kernel weight of the hybrid was significantly different from that of the MPV after 10 DAP (Fig. 1D). These results indicated that the heterosis of the seed coat area occurred prior to changes in seed weight. In the seed coat area, the relative growth rate of the hybrid and paternal line peaked at 3 DAP, and the maternal line grew most rapidly at 2 DAP (Fig. 1C). In terms of the relative increasing rate of kernel weight, the value of the hybrid was the highest at 3 DAP, while it was highest at 2 DAP for the maternal line and at 4 DAP for the paternal line (Fig. 1E). In general, seed coat area and kernel weight showed a large increase at 2-4 DAP. More importantly, 3 DAP was found to be a critical time point for hybrid seed coat development. Relative cell numbers in the seed coat To understand whether the increased seed coat area in the hybrid was due to changes in cell number or cell size, we counted the number of cells in the seed coat per unit area (0.1 mm 2 ) (Fig. 2). The middle longitudinal section of seed embryos collected at 8 DAP was obtained (Fig. 2A-C). The number of cells per unit area for the paternal line (319 cells per 0.1 mm 2 ) was higher than that for the hybrid (149.7 cells per 0.1 mm 2 ), the maternal line (138 cells per 0.1 mm 2 ) and the MPV (228.5 cells per 0.1 mm 2 ) at 8 DAP (Fig. 2D). Although there was no difference in the number of cells per unit area between the hybrid and maternal line at 8 DAP, there was a significant difference in seed coat surface area (Fig. 2D, E). Thus, compared with the maternal line, the cell number of the hybrid was determined by its larger seed coat area. The differences in cell number and seed coat area between the hybrid and MPV were significant (Fig. 2D), suggesting that the difference in seed coat size between the hybrid and parents might be the result of the coordination between cell number and cell size. Analysis of expressed gene in seed coat To investigate the expression of genes before and after the time points at which the seed coat of the hybrid line and its parents showed significant differences, we performed RNA-seq on 18 samples at 6 time points (2-4 DAP and 8-10 DAP) from the hybrid and parental lines, with 3 biological replicates for each sample. In total, 54 libraries and 2.4 billion high-quality reads were generated. The reads were then mapped to the maize B73 reference genome (RefGen_V4 [26]) by Hisat [27]. Approximately (87%) of reads were uniquely mapped (Table S2) and used to calculate the normalized gene expression level as fragments per kb of transcript per million mapped reads (FPKM). There was a high correlation among the three biological replicates at each time point (average R 2 =0.93) (Table S3). All the above results indicated that the sampling quality of this survey was high, and that the RNA-seq sequencing data were accurate and reliable. Thus, we took the average FPKM value of the three replicates as the expression level for the sample at each time point. The Principal Component Analysis (PCA) revealed that the 18 samples at six time points from three genotypes were assigned to two stages, and that 2-4 DAP were clearly distinguished from 8-10 DAP (Fig. S1). Moreover, the RNA-seq data of 2-4 DAP were overlapping (Fig. S1). There were 39005 annotated genes in the B73 AGP_v4 genome[26]. The transcriptomes of both the hybrid and parents exhibited very similar distributions in the number of expressed genes at each time point. In detail, the total numbers of genes expressed by the hybrid, maternal and paternal lines at six time points were 23516, 23596 and 23247, respectively. A total of 17691, 18328 and 17604 genes overlapped at these 6 time points, accounting for 75.23%, 77.67% and 75.73% of the total expressed genes (Fig. 3A-C), respectively. Taken together, the results of overlapping expressed genes in the hybrid and its parental lines showed that the growth and development of the seed coat was mainly due to gene expression level rather than the change in gene expression. Differentially expressed gene identification It was suggested that the hybrid showed significant differences from its parents at 8 DAP, while the relative growth rate of the hybrid seed coat area and kernel weight peaked at 3 DAP. Therefore, we focused on RNA-seq data at 3 DAP and 8 DAP for subsequent analysis. In DEG analysis, the gene expression levels in the hybrid were compared one-by-one with its parental lines, with a difference in the fold change (FC) ≥ 2 and (FC) ≤ -2 denoted as upregulated expression and downregulated expression, respectively. Some DEGs were highly or weakly expressed in either the hybrid or parental lines. At 3 DAP, a total of 3915 DEGs varied between the hybrid and maternal line (2134 (54.50%) upregulated and 1781 (45.50%) downregulated), and 2647 DEGs varied between the hybrid and paternal line (1248 (47.15%) upregulated and 1399 (52.85%) downregulated). Compared with that at 3 DAP, the proportion of DEGs upregulated at 8 DAP increased sharply. A total of 2182 DEGs varied between the hybrid and maternal line (1499 (68.70%) upregulated and 683 (31.30%) downregulated), and 2724 DEGs varied between the hybrid and paternal line (1847 (67.80%) upregulated and 877 (32.20%) downregulated) (Fig. 4A). The differences in the upregulation and/or downregulation of gene expression between the hybrid and its parents varied across comparisons at 3 DAP and 8 DAP. Some DEGs had genotype-specific expression and were thus identified as genotype-specific unigenes (Fig. 4B). Venn diagram analysis indicated that 1156 genes were DEGs (highlighted as black bold font) between the hybrid and its parents at both time points (Fig. 4B). Out of these 1156 genes, 834 DEGs (sum of bold font numerals) were commonly observed between M vs. H and P vs. H at 3 DAP. However, 430 DEGs (sum of bold font numerals) were commonly observed between M vs. H and P vs. H at 8 DAP (Fig. 4B). In the hybrids at both time points, 108 DEGs were found to be commonly differentially expressed (Fig. 4B). Expression patterns of DEGs To explore and categorize the trends for changes in expression, DEGs in 3 DAP and 8 DAP samples were divided into 13 possible hybrid and parental expression patterns (Table 1). Gene expression patterns were divided into additive expression and nonadditive expression. We found no significant changes in the number of additively expressed DEGs between the two time points. The expression pattern of most genes in the hybrid at the two time points was additive (85.18% at 3 DAP and 86.45% at 8 DAP). However, the effects of nonadditive DEGs were quite different, especially the dominantly expressed DEGs (classes 5, 7, 10 and 12). That is, the gene expression in the hybrid was similar to that of one parent, which accounted for more than 13% of all DEGs at both time points. At 3 DAP, the number of dominantly expressed genes with higher paternal expression accounted for the highest proportion in the nonadditive classes, while at 8 DAP, the number of dominantly expressed genes with higher maternal expression was the most abundant. The number of overdominantly expressed genes with transitional downregulation (0.35% at 3 DAP and 0.14% at 8 DAP) was relatively higher than that with transitional upregulation (0.10% at 3 DAP and 0.08% at 8 DAP). We summarized 13 possible hybrid and parental expression patterns into mid-parent, high-parent and low-parent categories to analyse the expression changes at both time points. By comparing the number of DEGs in three expression patterns at both time points, it was found that the number of low-parent DEGs changed the most, followed by the number of high-parent DEGs, while the number of mid-parents DEGs changed the least (Table S4). According to the analysis of the fold changes at both time points, it was suggested that the change in the number of low-parent and high-parent DEGs may play a role in hybrid seed coat heterosis. Functional classification by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis The DEGs with high and low parental expression at 3 DAP and 8 DAP were annotated by GO. The enrichment of GO terms with low-parental expression at 3 DAP showed that a large number of DEG categories related to stress were significantly enriched in biological processes, followed by cellular components, and categories related to the plasma membrane (Fig. 5A). The highly expressed DEG GO enrichment results at 3 DAP suggested that categories related to chromosome assembly and processing were enriched in cellular components, and categories related to the cell cycle and DNA replication were enriched in biological processes (Fig. 5B). The GO enrichment results for low-parental expression at 8 DAP demonstrated that the categories related to cytoskeleton and microtubule were enriched in cellular component, and the categories enriched in biological process were mainly related to cell cycle, nuclear division and microtubule (Fig. 5C). The GO enrichment results for highly expressed DEGs at 8 DAP indicated that categories related to metabolism were enriched (Fig. 5D). The KEGG pathway enrichment analysis was carried out to further understand the biological function of genes and their interactions. We performed KEGG enrichment analysis on DEGs with high and low-parental expression at 3 DAP and 8 DAP, respectively (Fig. 6). Among the 3 DAP low-parental expressed DEGs, metabolic pathways, biosynthesis of secondary metabolites, protein processing in the endoplasmic reticulum and glutathione metabolism were significantly enriched (Fig. 6A). For the 3 DAP high-parental expressed DEGs, the biosynthesis of secondary metals, DNA replication and photosynthesis were significantly enriched (Fig. 6B). Metabolic pathways, the biosynthesis of secondary metabolites and alpha linolenic acid metabolism were significantly enriched in 8 DAP low-parental expressed DEGs (Fig. 6C). In contrast to DEGs with high parental expression at 3 DAP, DNA replication was also enriched in the low-parental expression DEGs at 8 DAP. For the high-parental expression DEGs at 8 DAP, circadian rhythms, the biosynthesis of secondary metals, glycolysis/gluconeogenesis, thiamine metabolism, metabolic pathways, galactose metabolism and other pathways were significantly enriched (Fig. 6D). Weighted gene co-expression network analysis (WGCNA) To identify the specific genes highly related to seed coat area, we constructed co-expression networks by using the transcriptome data from comparisons of the hybrid and its parents at 6 time points, and associated the co-expression module with seed coat area traits. The results showed that these genes should be divided into 41 modules (Fig. S2). By observing the correlation between modules and samples, it was found that the absolute value of the correlation coefficient between the MEyellow (0.86) and MEsienna3 (0.82) modules was the highest (Fig. S2). In other words, these two modules have the highest correlation with seed coat area. In these two modules, the genes with the top 100 weight values were screened to construct the regulation network diagram (Fig. 7, Fig. S3). First, the module MEyellow revealed that Zm00001d019306 was a hub gene (Fig. 7). However, in module MEsienna3, the regulatory network was not as simple and clear as that in MEyellow. We chose genes with a degree value of more than 10 as hub genes. Therefore, the hub genes in module MEsienna3 were Zm00001d006330 , Zm00001d043289 , Zm00001d007181 , Zm00001d043290 , Zm00001d033510 , Zm00001d024732 , Zm00001d015025 , Zm00001d006331 and Zm00001d047796 (Fig. S3, Table S5). Validation of DEGs by quantitative real-time PCR The expression patterns of nine DEGs were further confirmed by quantitative real-time PCR (qRT‒PCR). We compared the transcriptional profiles of each gene with the seed coat RNA samples from the hybrid and its parents at different time points. The data confirmed that the expression patterns of all 9 DEGs were consistent with the expression levels obtained by RNA-seq analysis (Fig. 8, Fig. S4) Discussion Heterosis has been known for centuries [28]. Maize seeds are characterized as having heterosis, and seed size is determined by the maize seed coat area to a certain extent. The maize seed coat develops from the ovary wall and carries a maternal genotype. Due to its genetic stability and fixed heterozygosity, the seed coat of F2 seeds obtained by self-pollination of F 1 plants carries a heterozygous genotype, which makes it a good model for studying heterosis. Therefore, exploring the seed coat area of seeds obtained by self-pollination of maize F 1 plants and their self-crossed parents can deepen our understanding of maize seed coat heterosis, and in turn seed size and yield heterosis. The initiation of seed coat development depends on the fusion of the central cell with one of the sperm cells [29]. In our study, the relative growth rate of the seed coat area of the hybrid and paternal line peaked at 3 DAP, while the maternal line grew most rapidly with the peak at 2 DAP. This phenomenon may result from the fact that the synthesis of auxin production in the endosperm after fertilization removes the PRC2 block on seed coat development, to promote the coordinated development of the seed coat and endosperm [30]. There was a significant difference between the hybrid and mid-parental value for seed coat area at 8 DAP, suggesting that heterosis of the seed coat area might be the result of the coordination between cell number and cell size (Fig. 2). These results were consistent with the GO enrichment results from the transcriptome comparison between the hybrid and parents in the early stage of seed coat growth and development, suggesting that heterosis of the seed coat area resulted from the enhanced cell cycle and cell division activity in the early stage. Previous studies have reported that plant growth heterosis basically results from changes in cell number, similar to most plant characteristics [31]. However, the elongation of seed coat cells also affects seed size in Arabidopsis [32]. Our results showed that the maize seed coat developmental program of the hybrid was not dramatically altered, instead, cell proliferation and cell elongation changed. Transcriptome analysis is a well-established method to identify differential expression levels and regulatory mechanisms of genes from different samples at the transcriptional level [33]. Previous studies in plants have identified heterotic genes associated with many traits between hybrids and parents and have proved that gene dose effects play a role in heterosis [34-36]. RNA-seq was performed on the seed coats of the 2-4 DAP and 8-10 DAP hybrid and parental lines. The expression results indicated that the heterosis of seed coat extension was mainly due to the dosage effect of gene expression rather than none-or-all expression. The differential expression pattern of DEGs between the hybrid and its parents plays a significant role in heterosis [8, 37]. In hybrids, where 2 different alleles of a gene are combined, the allelic allelic expression may deviate from that of either parent or mid-parent predictions [38]. Concerning the relative gene expression level within a hybrid–parent triad, there are usually two scenarios about the relative gene expression level between the triples of hybrid parents. In the first case, the gene expression of the hybrid showed a cumulative pattern, which was contributed by each allele of parents. This pattern was additive and mainly regulated by cis action. In another case, the expression deviates from the mid-parental level, and other trans acting factors may lead to the change in expression for the corresponding allele in the hybrid [38]. In the present study, the majority (85.18% at 3 DAP and 86.45% at 8 DAP) of genes in the hybrid were additively expressed, suggesting that cis-acting elements played a major role in the regulation of gene expression. The nonadditive genes expressed throughout our entire dataset accounted for less than 15% at each time point but accounted for 54.35–56.24% of the DEGs. This result indicates that the difference in gene expression between the hybrid and parental lines was mainly due to nonadditive effects. Analysis of the expression of nonadditive genes indicated that the allelic expression pattern of the hybrid may not be a simple combination of alleles from both parents, but may be regulated by other genes or trans-action factors. By comparing the expression patterns of hybrids in 3 DAP and 8 DAP, the number and type of additive expressed genes changed slightly (19714 or 85.18% at 3 DAP to 19460 or 86.45% at 8 DAP DEGs), but the number and type of highly-expressed and weakly-expressed genes at 3 DAP and 8 DAP varied greatly (Table S4, Fig. 5). In conclusion, it was suggested that the expression of nonadditive genes in hybrids may be the main cause of heterosis in hybrid seed coats. Three DAP and 8 DAP are the key stages of maize endosperm cellularization and differentiation, respectively [20, 21]. In this study, the results of GO enrichment analysis and KEGG pathway analysis of different expression patterns for the seed coat in hybrid and parental lines at 3 DAP and 8 DAP matched the endosperm development process. In other words, at 3 DAP, compared with the parents, the hybrid increased its seed coat size mainly by enhancing cell cycle activity, weakening the response to stress and codeveloping with the endosperm (Fig. 5, A-B). However, at 8 DAP, the hybrid seed coat showed a reduced rate of cell division, instead by strengthening the metabolic process. At this stage, the seed coat changed mainly in coordination with the endosperm to enhance the synthesis and accumulation of nutrients in the seed, to enrich the grain and increase the kernel weight (Fig.5, C-D). This result was supported by the significant difference in kernel weight between the hybrid and MPV at 10 DAP (Fig. 1, D). The KEGG analysis results of genes with high parental expression at 8 DAP showed that, pathways related to circadian rhythm and metabolite accumulation were significantly enriched. Changing the circadian rhythm can improve the viability and biomass of hybrids [39]. In addition, evidence of changes in metabolic profiles has been recorded in hybrids [40, 41]. Moreover, seed development is accompanied by metabolic activities for the synthesis and accumulation of stored products, including protein and carbohydrates [42]. Therefore, in addition to genetic control, the final size of seeds may also be affected by metabolic activities [43]. Zm00001d018415 and Zm00001d051995 , two genes encoding proliferating cell nuclear antigen（PCNA）protein, showed high parental expression patterns in the GO and KEGG analyses at 3 DAP. PCNA is a key protein in the mitotic DNA replication mechanism of all eukaryotes [44]. It acts as a DNA polymerase δ cofactors and is thus necessary for the synthesis of new DNA strands by forming a homotrimeric complex [28, 45] that binds to various cyclin-CDK complexes [46-48]. In addition, as a sliding clamp, PCNA also plays a role in regulating DNA metabolism, processing DNA damage and controlling the cell cycle [49]. Therefore, the high parental expression of the PCNA gene and its related proteins at 3 DAP may play a role in enhancing DNA replication, repair and metabolism in the early stage of seed coat development, to increase the cell cycle rate and cell number in the hybrid. Moreover, the excessive transcript abundance of these genes in hybrids in the early stage of seed coat development may eventually help to maintain heterosis through dominant gene effects to a great extent. Zm00001d039131 , encodes ADP-glucose pyrophosphorylase, which catalyses the key steps of the starch synthesis pathway. This gene was enriched in GO terms and KEGG pathways related to metabolism and the biosynthesis of secondary metabolites at 8 DAP when high parental expression was present. This gene has been reported to be associated with seed development and filling and kernel weight [50-52]. Eight to 12 DAP is the stage of endosperm cell differentiation, during which nutrients and storage materials accumulate alongside cell proliferation and enlargement [21]. Seed development is a coordinated process among the seed coat, embryo and endosperm. ADP-glucose pyrophosphorylase was highly parentally expressed at 8 DAP in the hybrid, while the cell cycle-related genes showed a low-parental expression level. This phenomenon may be the result of the rapid development and nutrient filling of the endosperm through enhanced accumulation of metabolites in the hybrid. In other words, during this period, the gene expression pattern of the hybrid changed, that is, cell cycle activity was weakened and metabolic processes were enhanced, to stabilize heterosis by cooperating with the development and filling processes in the endosperm. WGCNA can be used to analyse the relationship between a gene set and sample phenotype, draw the regulatory network between genes in a gene set and identify key regulatory genes [53]. Zm00001d019306 , the hub gene of WGCNA in the highest correlation coefficient modules with the highest correlation coefficients, was significantly positively correlated with the seed coat phenotype and encoded, an indole-3-acetyl-leu hydrolase. This enzyme can catalyse the hydrolysis of indole-3-acetyl-L-leucine and release free indole-3-acetic acid (IAA) [54]. Seed development is regulated by phytohormones [43]. Auxin, as a central plant hormone, regulates many biological processes, including the maintenance of meristems, cell division and cell expansion [55-57]. In Arabidopsis, indole-3-acetyl-leu hydrolase regulates the auxin response by regulating auxin homeostasis in the endoplasmic reticulum [58]. Auxin is the key regulator of the plant defence response and the main control factor of plant growth and development. The change in gene activity related to the IAA pathway has potential importance in producing a heterosis phenotype [59, 60]. Previous studies showed that in Arabidopsis, the utilization of the IAA pathway caused differences in the control of auxin on cell proliferation by auxin, resulting in different levels of heterosis [61]. Interestingly, these views are consistent with our results. On the one hand, at 3 DAP, genes related to the defence response in hybrids were significantly downregulated. On the other hand, heterosis of the maize seed coat area is determined by the number of cells. In summary, indole-3-acetyl-leu hydrolase may play a significant role in of the IAA pathway and seed coat heterosis. It may regulate the auxin response by regulating auxin homeostasis in different seed coat development processes, to effectively maintain the balance of hybrid immunity and growth and development regarding resource allocation at different developmental time points of the seed coat. Finally, the expression levels of 9 genes identified by RNA-seq analysis were verified by qRT‒PCR, and the expression patterns of all tested genes were consistent with the results of RNA-seq, indicating the reliability of RNA-seq data in Yudan888 and its parents (Fig. 8). In conclusion, the genes identified in our study may provide some insights into broader aspects of seed coat heterosis in maize. However, further validation is needed to confirm the association between gene expression patterns and target agronomic traits in maize. Materials And Methods Field experiment and phenotyping In summer 2020, maize hybrid Yudan888, a widely cultivated maize variety in the Huang-Huai-Hai area of China, together with its parental lines 15S717 and T4691, was planted on the farm of Henan Agricultural University (Yuanyang, China, 113°16’ E, 35°41’ N). The planting density was 75,000 plants per hectare, with 3 replicates each containing 600 seedlings. A randomized block design was adopted in the field experiment. All the plants of the hybrid and parental lines were artificially self-crossed. Seed coat phenotyping began at 0 DAP and continued daily until 15 DAP, phenotyping was completed between 9 am and 11 am each day. The kernel weight was measured immediately after sampling by an analytical balance (ME204E, METTLER TOLEDO) with an accuracy of 0.0001 g. The whole seed coat was peeled off with tweezers and photographed, on a black background. Finally, the seed coat area was calculated by ImageJ software. At least six individual plants per DAP with no less than six kernels per plant were sampled for phenotyping. Univariate ANOVA was performed using IBM SPSS statistics 24 to test the significant differences in the measured the area of the seed coat between the hybrid and parental lines. Seeds from the centre of three different ears were fixed in FAA buffer (50% ethanol: formaldehyde: acetic acid = 90:5:5 volume) and then vacuumed at 4 °C for further imaging. Kernels from each sample were removed from the FAA solution and rinsed thoroughly in deionized water. Kernels were dissected at the longest part with a scalpel. Each dissected kernel was embedded in paraffin. Sections from kernels ranged from 8 to 10 μm in thickness. The maximum longitudinal cross section of intact kernels was dyed with 1% toluidine blue after dewaxing. Images of the sections were captured with a microscope (Axio Scope A1, Carl Zeiss) fitted with a digital camera (Axiocam 503 colour, Carl Zeiss). ImageJ was used to measure the number of seed coat cells per unit area. Sample collection and RNA extraction RNA-seq samples were taken during the period when hybrid and parental lines showed significant differences in the seed coat areas. The seed coat samples for RNA-seq were collected by manual dissection with a scalpel. Three biological replicates were set for each sample at each time point. Each replicate sample was collected from at least three individual ears. The samples collected were frozen immediately in liquid nitrogen and stored at -80 °C before RNA extraction. Total RNA was extracted from seed coat samples of the hybrid and its parental lines at 2-4 DAP and 8-10 DAP using TRIzol reagent (Invitrogen, United States). Three biological replicates were used for at each DAP. The quality and concentration of each RNA sample were determined using both gel electrophoresis and a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE). Only RNAs that met the criterion of an OD260/280 ratio of 1.9-2.1 were stored in a −80 °C freezer for further use. RNA-seq and data analysis A total of 1 μg of RNA was used as the input material for the construction of each RNA library, with the NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA). Briefly, mRNA was purified from total RNA by using poly-T oligo-attached magnetic beads. First strand cDNA was synthesized using random hexamer primers and M-MuLV Reverse Transcriptase (with RNase H). After second-strand cDNA synthesis, terminal repair, poly(A) tail addtion sequencing oligonucleotide adaptor ligation, the fragments were purified and subsequently amplified by PCR. The PCR products were assessed with the Agilent Bioanalyzer 2100 system. Finally, the libraries containing 150 bp cDNA inserts in size were generated and sequenced by a commercial service company (Guoke Biotechnology, Beijing, China) on an Illumia Noveseq platform. All data analyses were based on clean data with high quality [62]. Feature Counts v1.5.0-p3 was used to count the read numbers mapped to each gene (ftp://ftp.ensemblgenomes.org/pub/release-40/plants/fasta/zea_mays/dna/). Multiple testing with the Benjamini‒Hochberg approach for controlling the false discovery rate (FDR) was taken into account by using an adjusted P value. A corrected P value of 0.05 and an absolute fold change of two times were set as the thresholds for significantly differential expression. The normalized gene expression values with FPKM were used for PCA. PCA was performed in R 4.2.1. The GO enrichment analyses were obtained using agriGO v2.0 (http://systemsbiology.cau.edu.cn/agriGOv2/index.php) [63]. The KEGG annotations were obtained on the KEGG website (https://www.kegg.jp/kegg/) [64]. The GO and KEGG enrichment analyses were performed using TBtools [65]. Classification of differential gene expression patterns Based on previous studies on gene expression patterns [66-70], the FPKM values of genes were compared between the hybrid and its parents. The expression patterns were divided into additive and nonadditive (dominance or overdominance). A relative expression level of a hybrid gene that was located between parental values was classified as additive (regardless of whether there were DEGs between the parents or not). A relative expression level of DEGs similar to that of the dominant parents was classified as dominant. A relative expression level of DEGs that was higher than that of the parents (excessive upregulation) or lower than that of the parents (excessive downregulation) was classified as overdominance. Weighted gene co-expression network analysis (WGCNA) The assignments of gene co-expression modules using the WGCNA protocol were based on FPKM data (Langfelder and Horvath, 2008). To determine the module related to seed coat surface area, we associated the genes in each module with the seed coat area trait from 2-4 DAP and 8-10 DAP. The genes with an average FPKM＞1 in 54 samples were analyzed. The soft threshold power β for network construction was set to 11. The dynamic tree cutting algorithm with a minimum module size of 50 genes was used for hierarchical cluster cutting; a value of 0.15 was used to merge similar modules. If the p value of the module trait association was less than 0.05, it was defined as a meaningful module. In R 4.0.1, the WGCNA package was used to calculate the correlation weight of each gene and all other genes in the module. The top 100 gene sets in the module were visualized using Cytoscape V3.9.1, and the gene with the highest degree was considered to be a hub gene. qRT ‒ PCR Nine differentially expressed genes were selected for qRT‒PCR verification, and the first-strand cDNA was synthesized by reverse transcription using total RNA from the seed coat of the hybrid, the maternal line and the paternal line as templates. The primer 3 online toolbox was used to design gene specific primers for qRT‒PCR. The primers used for quantitative PCR are listed in Supplemental Table S6. The reverse transcription and qRT‒PCR reaction system and reaction procedures were in accordance with the instructions for the PrimeScriptTM RT Reagent Kit with gDNA Eraser (TaKaRa) and TB GreenTM Premix Ex TaqTM II (TaKaRa) Kit. Each experiment had 3 biological replicates and 3 technical replicates. The relative expression levels of these genes were analysed by the 2 - ΔΔ Ct method with Zm-actin-1 serving as the internal reference [71]. Abbreviations DAP: days after pollination MPV: mid-parental value PCA: Principal Component Analysis DEGs: differentially expressed genes GO: Gene Ontology KEGG: Kyoto Encyclopedia of Genes and Genomes WGCNA: Weighted gene co-expression network analysis qRT‒PCR: Real-time quantitative PCR Declarations Ethics approval and consent to participate All experimental studies on plants complied with relevant institutional, national, and international guidelines and legislation. Consent for publication Not applicable. Availability of data and materials The sequence datasets in fastq format of the current study are available in the NCBI Sequence Read Archive (SRA) database under Bioproject PRJNA828527 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA828527). Competing interests The authors declare that they have no competing interests. Funding This work was supported by a grant from the National Natural Foundation of China (31971961, 31871641, 32060460), the Science and Technique Foundation of Henan Province of China (202102110036) and the Basic Research Plan of Guizhou Province of China (zk [2022] ‒ 236). Author Contributions DD, ZHC and JHT designed the research. JL and LFW performed the experiments. JL, LFW and JW analysed the data. JL wrote the manuscript, and DD, ZHC and JHT revised the manuscript. All authors approved the final manuscript. Acknowledgements The authors would like to thank Huili Yang, Yanan Lin and Guoqiang Xu for assisting with this article. Authors' information 1 College of Agronomy, Hunan Agricultural University, Changsha 410128, China 2 National Key Laboratory of Wheat and Maize Crop Science; Collaborative Innovation Center of Henan Grain Crops; College of Agronomy; Henan Agricultural University, Zhengzhou 450002, China 3 Institute of Crop Germplasm Resources, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China 4 Hebi Academy of Agricultural Sciences, Hebi 458030, China 5 The Shennong laboratory, Zhengzhou 450002, China 6 Institute of Upland Food Crops, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China Juan Li 1,2,3# , Liangfa Wang 1,2,4# , Jiong Wan 2 , Kuntai Dang 2 , Yuan Lin 4 , Shujun Meng 2 , Xiaoqian Qiu 2 , Qiyue Wang 4 , Jiawen Zhao 2 , Liqin Mu 2 , Hongbing Luo 1 , Dong Ding 2* , Zehui Chen 1,6* , and Jihua Tang 2,5* References Hochholdinger F, Hoecker N: Towards the molecular basis of heterosis . Trends Plant Sci 2007, 12 (9):427-432. Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA: Heterosis . Plant Cell 2010, 22 (7):2105-2112. Hubner N, Wallace CA, Zimdahl H, Petretto E, Schulz H, Maciver F, Mueller M, Hummel O, Monti J, Zidek V et al : Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease . Nat Genet 2005, 37 (3):243-253. Song RT, Messing J: Gene expression of a gene family in maize based on noncollinear haplotypes . Proc Natl Acad Sci U S A 2003, 100 (15):9055-9060. Lippman ZB, Zamir D: Heterosis: revisiting the magic . Trends Genet 2007, 23 (2):60-66. Charlesworth D, Willis JH: The genetics of inbreeding depression . Nat Rev Genet 2009, 10 (11):783-796. 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-7701. He GM, Zhu XP, Elling AA, Chen LB, Wang XF, Guo L, Liang MZ, He H, Zhang HY, Chen FF et al : Global Epigenetic and Transcriptional Trends among Two Rice Subspecies and Their Reciprocal Hybrids . Plant Cell 2010, 22 (1):17-33. Groszmann M, Greaves IK, Albertyn ZI, Scofield GN, Peacock WJ, Dennis ES: Changes in 24-nt siRNA levels in Arabidopsis hybrids suggest an epigenetic contribution to hybrid vigor . Proc Natl Acad Sci U S A 2011, 108 (6):2617-2622. Su N, Sullivan JA, Deng XW: Modulation of F1 hybrid stature without altering parent plants through trans-activated expression of a mutated rice GAI homologue . Plant Biotechnol J 2005, 3 (2):157-164. Xue WY, Xing YZ, Weng XY, Zhao Y, Tang WJ, Wang L, Zhou HJ, Yu SB, Xu CG, Li XH et al : Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice . Nat Genet 2008, 40 (6):761-767. Krieger U, Lippman ZB, Zamir D: The flowering gene SINGLE FLOWER TRUSS drives heterosis for yield in tomato . Nat Genet 2010, 42 (5):459-463. 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-633. Li X, Li XR, Fridman E, Tesso TT, Yu JM: Dissecting repulsion linkage in the dwarfing gene Dw3 region for sorghum plant height provides insights into heterosis . Proc Natl Acad Sci U S A 2015, 112 (38):11823-11828. Zhang X, Hirsch CN, Sekhon RS, de Leon N, Kaeppler SM: Evidence for maternal control of seed size in maize from phenotypic and transcriptional analysis . J Exp Bot 2016, 67 (6):1907-1917. Sundaresan V: Control of seed size in plants . Proc Natl Acad Sci U S A 2005, 102 (50):17887-17888. Adamski NM, Anastasiou E, Eriksson S, O'Neill CM, Lenhard M: Local maternal control of seed size by KLUH/CYP78A5-dependent growth signaling . Proc Natl Acad Sci U S A 2009, 106 (47):20115-20120. Verdier J, Dessaint F, Schneider C, Abirached-Darmency M: A combined histology and transcriptome analysis unravels novel questions on Medicago truncatula seed coat . J Exp Bot 2013, 64 (2):459-470. Bennetzen JL, Hake SC: Handbook of maize: its biology . Springer New York 2009, 10.1007/978-0-387-79418-1 . Sabelli PA, Larkins BA: The contribution of cell cycle regulation to endosperm development . Sex Plant Reprod 2009, 22 (4):207-219. Ottaviano E, Petroni D, Pe ME: Gametophytic Expression of Genes-Controlling Endosperm Development in Maize . Theor Appl Genet 1988, 75 (2):252-258. Hehenberger E, Kradolfer D, Kohler C: Endosperm cellularization defines an important developmental transition for embryo development . Development 2012, 139 (11):2031-2039. Nowack MK, Shirzadi R, Dissmeyer N, Dolf A, Endl E, Grini PE, Schnittger A: Bypassing genomic imprinting allows seed development . Nature 2007, 447 (7142):312-315. Thorne JH: Phloem Unloading of C-Assimilates and N-Assimilates in Developing Seeds . Annu Rev Plant Physiol Plant Mol Biol 1985, 36 :317-343. Western TL, Skinner DJ, Haughn GW: Differentiation of mucilage secretory cells of the Arabidopsis seed coat . Plant Physiol 2000, 122 (2):345-356. Jiao Y, Peluso P, Shi J, Liang T, Stitzer MC, Wang B, Campbell MS, Stein JC, Wei X, Chin CS et al : Improved maize reference genome with single-molecule technologies . Nature 2017, 546 (7659):524-527. Kim D, Langmead B, Salzberg SL: HISAT: a fast spliced aligner with low memory requirements . Nat Methods 2015, 12 (4):357-360. Darwin C: The effects of cross and self fertilisation in the vegetable kingdom . New York,: D. Appleton and company; 1902. Roszak P, Kohler C: Polycomb group proteins are required to couple seed coat initiation to fertilization . Proc Natl Acad Sci U S A 2011, 108 (51):20826-20831. Figueiredo DD, Batista RA, Roszak PJ, Kohler C: Auxin production couples endosperm development to fertilization . Nat Plants 2015, 1 (12). East EM: Heterosis . Genetics 1936, 21 (4):375-397. Garcia D, Fitz Gerald JN, Berger F: Maternal control of integument cell elongation and zygotic control of endosperm growth are coordinated to determine seed size in Arabidopsis . Plant Cell 2005, 17 (1):52-60. Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics . Nat Rev Genet 2009, 10 (1):57-63. Yao H, Dogra Gray A, Auger DL, Birchler JA: Genomic dosage effects on heterosis in triploid maize . Proc Natl Acad Sci U S A 2013, 110 (7):2665-2669. Birchler JA, Riddle NC, Auger DL, Veitia RA: Dosage balance in gene regulation: biological implications . Trends Genet 2005, 21 (4):219-226. Birchler JA, Johnson AF, Veitia RA: Kinetics genetics: Incorporating the concept of genomic balance into an understanding of quantitative traits . Plant Sci (Amsterdam, Neth) 2016, 245 :128-134. Zhang HY, He H, Chen LB, Li L, Liang MZ, Wang XF, Liu XG, He GM, Chen RS, Ma LG et al : A genome-wide transcription analysis reveals a close correlation of promoter INDEL polymorphism and heterotic gene expression in rice hybrids . Mol Plant 2008, 1 (5):720-731. Birchler JA, Auger DL, Riddle NC: In search of the molecular basis of heterosis . Plant Cell 2003, 15 (10):2236-2239. Ni Z, Kim ED, Ha M, Lackey E, Liu J, Zhang Y, Sun Q, Chen ZJ: Altered circadian rhythms regulate growth vigour in hybrids and allopolyploids . Nature 2009, 457 (7227):327-331. Gartner T, Steinfath M, Andorf S, Lisec J, Meyer RC, Altmann T, Willmitzer L, Selbig J: Improved heterosis prediction by combining information on DNA- and metabolic markers . PLoS ONE 2009, 4 (4):e5220. Fievet JB, Dillmann C, de Vienne D: Systemic properties of metabolic networks lead to an epistasis-based model for heterosis . Theor Appl Genet 2010, 120 (2):463-473. Sreenivasulu N, Wobus U: Seed-Development Programs: A Systems Biology-Based Comparison Between Dicots and Monocots . Annual Review of Plant Biology, Vol 64 2013, 64 :189-+. Borisjuk L, Rolletschek H, Radchuk R, Weschke W, Wobus U, Weber H: Seed development and differentiation: a role for metabolic regulation . Plant Biol (Stuttg) 2004, 6 (4):375-386. Tan CK, Castillo C, So AG, Downey KM: An auxiliary protein for DNA polymerase-delta from fetal calf thymus . J Biol Chem 1986, 261 (26):12310-12316. de la Paz Sanchez M, Torres A, Boniotti MB, Gutierrez C, Vazquez-Ramo JM: PCNA protein associates to Cdk-A type protein kinases in germinating maize . Plant Mol Biol 2002, 50 (2):167-175. Celis JE, Madsen P, Celis A, Nielsen HV, Gesser B: Cyclin (Pcna, Auxiliary Protein of DNA Polymerase-Delta) Is a Central Component of the Pathway(S) Leading To DNA-Replication and Cell-Division . FEBS Lett 1987, 220 (1):1-7. Prelich G, Tan CK, Kostura M, Mathews MB, So AG, Downey KM, Stillman B: Functional Identity of Proliferating Cell Nuclear Antigen and a DNA Polymerase-Delta Auxiliary Protein . Nature 1987, 326 (6112):517-520. Sanchez MD, Gurusinghe SH, Bradford KJ, Vazquez-Ramos JM: Differential response of PCNA and Cdk-A proteins and associated kinase activities to benzyladenine and abscisic acid during maize seed germination . J Exp Bot 2005, 56 (412):515-523. Strzalka W, Ziemienowicz A: Proliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation . Ann Bot-London 2011, 107 (7):1127-1140. Kawagoe Y, Kubo A, Satoh H, Takaiwa F, Nakamura Y: Roles of isoamylase and ADP-glucose pyrophosphorylase in starch granule synthesis in rice endosperm . Plant J 2005, 42 (2):164-174. Na G, Aryal N, Fatihi A, Kang J, Lu C: Seed-specific suppression of ADP-glucose pyrophosphorylase in Camelina sativa increases seed size and weight . Biotechnol Biofuels 2018, 11 :330. Prathap V, Tyagi A: Correlation between expression and activity of ADP glucose pyrophosphorylase and starch synthase and their role in starch accumulation during grain filling under drought stress in rice . Plant Physiol Biochem 2020, 157 :239-243. Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis . BMC Bioinformatics 2008, 9 :559. Bartel B, Leclere S, Magidin M, Zolman BK: Inputs to the Active Indole-3-Acetic Acid Pool: De Novo Synthesis, Conjugate Hydrolysis, and Indole-3-Butyric Acid b-Oxidation . J Plant Growth Regul 2001, 20 (3):198-216. Chapman EJ, Estelle M: Mechanism of auxin-regulated gene expression in plants . Annu Rev Genet 2009, 43 :265-285. Vanneste S, Friml J: Auxin: a trigger for change in plant development . Cell 2009, 136 (6):1005-1016. Weijers D, Wagner D: Transcriptional Responses to the Auxin Hormone . Annu Rev Plant Biol 2016, 67 :539-574. Rampey, Rebekah A: The roles of Arabidopsis IAA-amino acid conjugate hydrolases in auxin homeostasis . Kazan K, Manners JM: Linking development to defense: auxin in plant-pathogen interactions . Trends Plant Sci 2009, 14 (7):373-382. Busov VB, Brunner AM, Strauss SH: Genes for control of plant stature and form . New Phytol 2008, 177 (3):589-607. Groszmann M, Gonzalez-Bayon R, Lyons RL, Greaves IK, Kazan K, Peacock WJ, Dennis ES: Hormone-regulated defense and stress response networks contribute to heterosis in Arabidopsis F1 hybrids . Proc Natl Acad Sci U S A 2015, 112 (46):E6397-E6406. Bolger AM, Lohse M, Usadel B: Trimmomatic: a flexible trimmer for Illumina sequence data . Bioinformatics 2014, 30 (15):2114-2120. Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, Xu W, Su Z: agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update . Nucleic Acids Res 2017, 45 (W1):W122-W129. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M: KEGG: integrating viruses and cellular organisms . Nucleic Acids Res 2021, 49 (D1):D545-D551. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R: TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data . Mol Plant 2020, 13 (8):1194-1202. Yoo MJ, Szadkowski E, Wendel JF: Homoeolog expression bias and expression level dominance in allopolyploid cotton . Heredity (Edinb) 2013, 110 (2):171-180. Zhou P, Hirsch CN, Briggs SP, Springer NM: Dynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development . Mol Plant 2019, 12 (3):410-425. Ding HP, Qin C, Luo XR, Li LJ, Chen Z, Liu HJ, Gao J, Lin HJ, Shen YO, Zhao MJ et al : Heterosis in Early Maize Ear Inflorescence Development: A Genome-Wide Transcription Analysis for Two Maize Inbred Lines and Their Hybrid . Int J Mol Sci 2014, 15 (8):13892-13915. Li AL, Liu DC, Wu J, Zhao XB, Hao M, Geng SF, Yan J, Jiang XX, Zhang LQ, Wu JY et al : mRNA and Small RNA Transcriptomes Reveal Insights into Dynamic Homoeolog Regulation of Allopolyploid Heterosis in Nascent Hexaploid Wheat . Plant Cell 2014, 26 (5):1878-1900. Howlader J, Robin AHK, Natarajan S, Biswas MK, Sumi KR, Song CY, Park JI, Nou IS: Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum) . Sci Rep-Uk 2020, 10 (1). Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method . Methods 2001, 25 (4):402-408. Howlader J, Robin AHK, Natarajan S, Biswas MK, Sumi KR, Song CY, Park JI, Nou IS: Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum) . Sci Rep 2020, 10 (1):9082. Tables Table 1 is available in the Supplementary Files section. Additional Declarations No competing interests reported. Supplementary Files Fig.S2.pdf TableS1.xlsx TableS2.xlsx TableS4.xlsx TableS5.xlsx TableS6.xlsx Table1.docx Cite Share Download PDF Status: Published Journal Publication published 07 Feb, 2023 Read the published version in BMC Plant Biology → Version 2 posted Editorial decision: Accepted 19 Jan, 2023 Reviews received at journal 17 Jan, 2023 Reviewers agreed at journal 17 Jan, 2023 Reviewers agreed at journal 16 Jan, 2023 Reviewers invited by journal 13 Jan, 2023 Submission checks completed at journal 12 Jan, 2023 First submitted to journal 10 Jan, 2023 You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-1768077","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[{"code":1,"date":"2022-06-27 14:03:38","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":175483414,"identity":"a8965e36-4e44-437b-8e8f-908be58b3d2f","order_by":0,"name":"Juan Li","email":"","orcid":"","institution":"Hunan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Juan","middleName":"","lastName":"Li","suffix":""},{"id":175483415,"identity":"f3e2cb11-1c89-47ba-89e9-c9b3fdfbf21e","order_by":1,"name":"Liangfa Wang","email":"","orcid":"","institution":"Hunan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Liangfa","middleName":"","lastName":"Wang","suffix":""},{"id":175483416,"identity":"624afe43-6095-4a3f-a427-192c72beb3ac","order_by":2,"name":"Jiong Wan","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jiong","middleName":"","lastName":"Wan","suffix":""},{"id":175483417,"identity":"4e4d3869-4a3b-490b-80c5-916170c4e741","order_by":3,"name":"Kuntai Dang","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Kuntai","middleName":"","lastName":"Dang","suffix":""},{"id":175483418,"identity":"794285d1-d060-45c5-8994-2e4fbe6a83d1","order_by":4,"name":"Yuan Lin","email":"","orcid":"","institution":"Hebi Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yuan","middleName":"","lastName":"Lin","suffix":""},{"id":175483419,"identity":"9ca256dd-126f-4809-9f78-36d31d4a698d","order_by":5,"name":"Shujun Meng","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Shujun","middleName":"","lastName":"Meng","suffix":""},{"id":175483420,"identity":"8c3b7092-aa05-443d-92a8-3780a9f07991","order_by":6,"name":"Xiaoqian Qiu","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xiaoqian","middleName":"","lastName":"Qiu","suffix":""},{"id":175483421,"identity":"7f0330e4-d5fb-49ff-8a45-0ad77570c73e","order_by":7,"name":"Qiyue Wang","email":"","orcid":"","institution":"Hebi Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Qiyue","middleName":"","lastName":"Wang","suffix":""},{"id":175483422,"identity":"c9d0b1a5-05ff-46bc-bc25-ec8d97cb2573","order_by":8,"name":"Jiawen Zhao","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jiawen","middleName":"","lastName":"Zhao","suffix":""},{"id":175483423,"identity":"2de2e284-63cf-4fe1-81cd-9cdbbf0e21f0","order_by":9,"name":"Liqin Mu","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Liqin","middleName":"","lastName":"Mu","suffix":""},{"id":175483424,"identity":"c4a3af3d-0f11-421a-a0a2-47fc0ce90334","order_by":10,"name":"Hongbing Luo","email":"","orcid":"","institution":"Hunan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Hongbing","middleName":"","lastName":"Luo","suffix":""},{"id":175483425,"identity":"1be26532-9fa9-4311-a669-0fcd76022530","order_by":11,"name":"Dong Ding","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABCElEQVRIiWNgGAWjYDCCw0DMA8Rs7I2NDz42gIQYGw8QpYWf53Cz4cwGBgmglgb8Wg5AtUjOSG8T5gVrgQjiBHzHeQ+/eFNxx27DgcQ2ZtsdNnW67YeBttTYROPSInmYL81yzplnyRsOHGx7nHsmTcLsTCJQy7G03AYcWgwO85gZ87YdTjY42NhunNt2WMLsAFALY8NhIrQcZmyTtgRpOf+QoBbjx0AtdpJtQC2MIC03CNgiCbSFcc6Zwwn8PIzNhr1taZLbbgBtScDjF77zZ4w/vKk4bM8m//zhg59tNvxm59MfPvhQY4NTCxCwgaIiEVVBAm7lIMD8AUjY41czCkbBKBgFIxoAAL6zalL1M//lAAAAAElFTkSuQmCC","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Dong","middleName":"","lastName":"Ding","suffix":""},{"id":175483426,"identity":"26804ef8-2aa1-4847-9cdd-e330af8b90d5","order_by":12,"name":"Zehui Chen","email":"","orcid":"","institution":"Institute of Upland Food Crops, Guizhou Academy of Agricultural Sciences","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Zehui","middleName":"","lastName":"Chen","suffix":""},{"id":175483427,"identity":"f8a50c37-440e-4b29-92e5-04536b9711e2","order_by":13,"name":"Jihua Tang","email":"","orcid":"","institution":"Henan Agricultural University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jihua","middleName":"","lastName":"Tang","suffix":""}],"badges":[],"createdAt":"2022-06-17 09:44:26","currentVersionCode":2,"declarations":"","doi":"10.21203/rs.3.rs-1768077/v2","doiUrl":"https://doi.org/10.21203/rs.3.rs-1768077/v2","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12870-023-04078-1","type":"published","date":"2023-02-07T18:45:45+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":32873896,"identity":"64220be1-a9f1-4e74-ba73-60efb0142f1b","added_by":"auto","created_at":"2023-02-13 22:56:03","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":546196,"visible":true,"origin":"","legend":"\u003cp\u003eComparisons of phenotypic variations in the hybrid and its parental lines during seed coat development.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eA \u003c/strong\u003eImage of seeds from the hybrid (H) and its parental lines (M for maternal line and P for paternal line) at 0-15 DAP. Scale bar, 2 millimetres. \u003cstrong\u003eB\u003c/strong\u003e The variation in seed coat area at 0-15 DAP of the hybrid and the mid-parent value. \u003cstrong\u003eC\u003c/strong\u003e The relative growth rate of seed coat area in the hybrid and its parental lines on 0-15 DAP. \u003cstrong\u003eD\u003c/strong\u003e The variation in kernel weight during 0-15 DAP of the hybrid and the mid-parent value. \u003cstrong\u003eE\u003c/strong\u003e The relative growth rate of kernel weight in the hybrid and its parental lines at 0-15 DAP.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/a310cc5b27ea629f121a0169.png"},{"id":32875712,"identity":"176fa3e3-83e9-4405-9e72-a3f3bd7d8dde","added_by":"auto","created_at":"2023-02-13 23:20:03","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":991185,"visible":true,"origin":"","legend":"\u003cp\u003eSeed coat cell number per unit area at 8 DAP in the hybrid and its parental lines. \u003cstrong\u003eA\u003c/strong\u003e Micrographs of the seed coat cell number per unit area in the maternal line. \u003cstrong\u003eB\u003c/strong\u003e Micrograph of seed coat cell numbers per unit area in the hybrid. \u003cstrong\u003eC\u003c/strong\u003e Micrograph of seed coat cell numbers per unit area in the paternal line. \u003cstrong\u003eD\u003c/strong\u003e Comparison of the number of seed coat cells per unit area between the hybrid (H) and maternal line (M), paternal line (P) and mid-parent value (MPV). The cell numbers were counted in a 0.1 mm\u003csup\u003e2\u003c/sup\u003e area, with 20 um scale bars. \u003cstrong\u003eE\u003c/strong\u003e Comparison of seed coat area between the hybrid (H) and maternal line (M).\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/211dd4da621f895e4bb19ae0.png"},{"id":32875092,"identity":"b2b53d3b-bfe5-4d18-b5e3-61b11c186d6f","added_by":"auto","created_at":"2023-02-13 23:12:03","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":155098,"visible":true,"origin":"","legend":"\u003cp\u003eVenn diagram analyses of expressed genes at 6 time points in the hybrid and its parental lines. \u003cstrong\u003eA\u003c/strong\u003e hybrid Yudan888. \u003cstrong\u003eB\u003c/strong\u003ematernal line. \u003cstrong\u003eC\u003c/strong\u003e paternal line\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/dd12de140cac4c7c77ff09a5.png"},{"id":32873911,"identity":"a1a62618-f564-45d3-9360-4b798463b3b5","added_by":"auto","created_at":"2023-02-13 22:56:03","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":174974,"visible":true,"origin":"","legend":"\u003cp\u003eIdentification of differentially expressed genes (DEGs) between the hybrid and parents on 3 DAP and 8 DAP. \u003cstrong\u003eA\u003c/strong\u003e Statistics of up- or downregulated genes between the hybrid and parents are shown by coloured columns at 3 DAP and 8 DAP. \u003cstrong\u003eB\u003c/strong\u003e Total numbers of DEGs in the maternal line (M) vs. hybrid (H) and paternal line (P) vs. hybrid(H) on 3 DAP, maternal line (M) vs. hybrid (H) and paternal line (P) vs. hybrid (H) on 8 DAP by Venn diagram analysis. Bold numerals indicate commonly expressed DEGs between the hybrid and parents.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/fd671b9793afa3051a919e53.png"},{"id":32876775,"identity":"347d086f-be04-4424-aa24-8a71979f70bd","added_by":"auto","created_at":"2023-02-13 23:28:03","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":296380,"visible":true,"origin":"","legend":"\u003cp\u003eGO enrichment analysis of genes with low- parental expressionand high -parental expression at 3 DAP and 8 DAP. \u003cstrong\u003eA\u003c/strong\u003e The top 20 GO terms of the shared genes with low- parental expression at 3 DAP. \u003cstrong\u003eB\u003c/strong\u003e The top 20 GO terms of the shared genes with high- parental expression of 3 DAP. \u003cstrong\u003eC\u003c/strong\u003eThe top 20 GO terms of the shared genes with low- parental expression at 8 DAP. \u003cstrong\u003eD\u003c/strong\u003e The top 20 GO terms of the shared genes with high- parental expression at 8 DAP.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/f38021de34729d9a03cd9693.png"},{"id":32874676,"identity":"acfd8316-462d-4f62-b62d-6b943664dce1","added_by":"auto","created_at":"2023-02-13 23:04:03","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":293433,"visible":true,"origin":"","legend":"\u003cp\u003eKEGG analysis of genes with low- parental expression and high- parental expression at 3 DAP and 8 DAP. \u003cstrong\u003eA\u003c/strong\u003eThe top 20 KEGG pathways of the shared genes with low- parental expression at 3 DAP. \u003cstrong\u003eB\u003c/strong\u003e The top 20 KEGG pathways of the shared genes with high- parental expression at 3 DAP. \u003cstrong\u003eC\u003c/strong\u003e The top 20 KEGG pathways of the shared genes with low- parental expression at 8DAP. \u003cstrong\u003eD\u003c/strong\u003e The top 20 KEGG pathways of the shared genes with high-parental expression at 8 DAP.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/b1e15e68618db184e3e8d8a7.png"},{"id":32875099,"identity":"2ba7f959-3db4-417e-bc38-fcd5ecdab199","added_by":"auto","created_at":"2023-02-13 23:12:03","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":325582,"visible":true,"origin":"","legend":"\u003cp\u003eThe gene network of \u003cem\u003eZm00001d019306\u003c/em\u003e.\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/bd09a9f22bcf6ce183ed9175.png"},{"id":32876776,"identity":"d25d9b15-0999-4080-89d6-d39f4b187142","added_by":"auto","created_at":"2023-02-13 23:28:03","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":386848,"visible":true,"origin":"","legend":"\u003cp\u003eThe relative expression levels of 9 differentially expressed genes at 8 DAP as assessed by qRT‒PCR. The qRT‒PCR relative expression levels are shown in the column diagram in the figures, and the FPKM value obtained by RNA-seq is represented by the lines in the figure.\u003c/p\u003e","description":"","filename":"8.png","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/92b2d1b78df559efe98f7937.png"},{"id":44719315,"identity":"5244057e-11db-4f74-91d7-09bf93955c22","added_by":"auto","created_at":"2023-10-16 18:54:26","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4629548,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/473c55ac-258e-4ec8-ad6b-08627fb8467a.pdf"},{"id":32873897,"identity":"7e2955d1-e5d8-4bd3-ab48-a4550d9a088b","added_by":"auto","created_at":"2023-02-13 22:56:03","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":10594,"visible":true,"origin":"","legend":"","description":"","filename":"Fig.S2.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/60b652d1bfb7c6d22d6102ac.pdf"},{"id":32874670,"identity":"7bdd43f3-9092-4dee-b1e1-c683ae2cc585","added_by":"auto","created_at":"2023-02-13 23:04:03","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":11247,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/2149505536b1f8f416f42c7f.xlsx"},{"id":32875095,"identity":"63d67831-50be-4d24-8a14-2172d11204cd","added_by":"auto","created_at":"2023-02-13 23:12:03","extension":"xlsx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":15897,"visible":true,"origin":"","legend":"","description":"","filename":"TableS2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/55311f21036b95acf9aa4308.xlsx"},{"id":32875101,"identity":"4c779307-1625-4b82-b514-dd837520e60a","added_by":"auto","created_at":"2023-02-13 23:12:03","extension":"xlsx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":9342,"visible":true,"origin":"","legend":"","description":"","filename":"TableS4.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/36cf351db67557101fba6e8f.xlsx"},{"id":32875713,"identity":"f4ced53d-2e30-4f98-84b7-c7dd7d793b10","added_by":"auto","created_at":"2023-02-13 23:20:03","extension":"xlsx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":10488,"visible":true,"origin":"","legend":"","description":"","filename":"TableS5.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/4500079aa16f7b88ea088334.xlsx"},{"id":32875718,"identity":"a13e0df1-e4ac-4fd4-a8c8-2a6dd1f167ca","added_by":"auto","created_at":"2023-02-13 23:20:03","extension":"xlsx","order_by":8,"title":"","display":"","copyAsset":false,"role":"supplement","size":10915,"visible":true,"origin":"","legend":"","description":"","filename":"TableS6.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/f35b9f3402a6a0a506f0a517.xlsx"},{"id":32873902,"identity":"2d8300fd-eaed-433f-9ffd-a500020372c1","added_by":"auto","created_at":"2023-02-13 22:56:03","extension":"docx","order_by":9,"title":"","display":"","copyAsset":false,"role":"supplement","size":47357,"visible":true,"origin":"","legend":"","description":"","filename":"Table1.docx","url":"https://assets-eu.researchsquare.com/files/rs-1768077/v2/a94ceb0c9e644ba8ce8b1344.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Dynamic Patterns of Gene Expression and Regulatory Variation in the Maize Seed Coat","fulltext":[{"header":"Background","content":"\u003cp\u003eHeterosis has been widely used in plant breeding[1].\u0026nbsp;Due to heterosis, hybrid F\u003csub\u003e1\u003c/sub\u003e generations are superior to both parents in terms of viability, plant height, yield, stress and disease resistance\u003cem\u003e\u0026nbsp;\u003c/em\u003e[2]. It is hypothesized that the difference in gene expression between a hybrid and its parents is one of the main causes of heterosis\u0026nbsp;[3, 4]. Genetic models to explain the increased yield of hybrids consider the interaction of alleles at many loci generating altered expression levels and patterns\u0026nbsp;[5, 6].\u0026nbsp;Whole-genome transcriptomes have shown that in developing rice leaves and panicles, differentially expressed genes (DEGs) between hybrids and their parents were overrepresented during energy metabolism and transport\u0026nbsp;[7].\u0026nbsp;Changes in epigenetic states such as DNA methylation, small RNA production and histone modification have been found in hybrid genomes\u0026nbsp;[8, 9]. Genetic mapping has been used to identify the detailed genetic loci contributing to hybrid performances and investigate their mechanisms, in a number of plant species including Arabidopsis, rice, maize, sorghum, and tomato\u003cem\u003e\u0026nbsp;\u003c/em\u003e[10-14]. Rice GAI homologues ectopically expressed in the F\u003csub\u003e1\u003c/sub\u003e generation of \u003cem\u003eArabidopsis thaliana\u003c/em\u003e caused significant plant height changes\u0026nbsp;[10].\u0026nbsp;The major quantitative trait gene \u003cem\u003eGhd7\u003c/em\u003e, which encodes a CCT domain protein in rice, was identified as a key factor for improving the yield and adaptability of the excellent hybrid Shanyou63 and other indica rice varieties\u0026nbsp;[11]. The orthologue of \u0026nbsp;the flowering locus \u003cem\u003eSFT\u003c/em\u003e for tomato\u0026nbsp;[12]\u0026nbsp;and \u003cem\u003eHd3a\u003c/em\u003e for rice\u0026nbsp;[13], has been shown to have single-gene overdominance\u0026nbsp;in yield traits.\u0026nbsp;However,\u0026nbsp;there have been many studies on heterosis, and the genetic mechanism and molecular basis of heterosis remain unclear.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSeed size is an important component of grain yield and a key trait for crop domestication\u0026nbsp;[15].\u0026nbsp;Seed size is\u0026nbsp;coordinately\u0026nbsp;determined by the development of the triploid endosperm, diploid embryo, and seed coat\u0026nbsp;[7, 16]. The seed coat\u0026nbsp;which\u0026nbsp;is developed from the integument and shares the maternal genotype\u0026nbsp;protects the embryo and endosperm and grows\u0026nbsp;coordinately\u0026nbsp;with them[17]. It plays an important role in the seed filling process\u0026nbsp;[18]. Therefore, the seed coat determines the yield of cereal crops to a certain extent.\u0026nbsp;The development of the seed coat and endosperm takes precedence over that of the embryo. In the early stage of seed development, from double fertilization to 6 days after pollination (DAP), the growth of the seed coat is accompanied by enlargement of the endosperm. Within 72 hours after fertilization, the endosperm converts to the multinuclear stage\u0026nbsp;[19]. At 3-6 DAP, the endosperm of maize is in the cellularization stage, during which the cytoskeleton recognizes the nuclei and forms the microtubule system between the nuclei. All nuclei are separated and cytoplasmic division is carried out at the same time to form normal cells\u0026nbsp;[20]. Eight to 12 DAP is the endosperm cell differentiation stage, during which nutrients and storage materials accumulate alongside cell proliferation and enlargement\u0026nbsp;[21].\u0026nbsp;The morphological changes of the seed coat occur only when the endosperm begins to develop, and the development of the embryo does not promote the growth and differentiation of the seed coat\u0026nbsp;[22, 23].\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eThus, the seed coat coincides with the enlargement of the seed in the early stage of development.\u0026nbsp;In addition, the seed coat regulates the accumulation of phloem assimilates in seeds during seed development\u0026nbsp;[24]. At the beginning of the filling stage at 12 DAP, the seed coat thickness is continuously compressed until the seed coat thickness decreases to only a few layers of cells. At this stage, the seed coat cell wall thickens, providing physical protection for future seed dormancy. The function of the seed coat during development may be to use its own starch for endosperm expansion and to provide more space for endosperm development.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAn impressive array of genomic and regulatory variations in kernel traits has been observed between maize inbred lines. However, how these variations affect gene expression in hybrids remains unknown.\u0026nbsp;Since the integument does not participate in the fertilization process, the seed coat differentiates primarily from the integument is determined by the maternal plant and shares the maternal genotype[25]. \u0026nbsp;Therefore, the genotype of the seed coat of the F\u003csub\u003e1\u003c/sub\u003e hybrid and its parents is still the same as that of the maternal plant after fertilization.\u0026nbsp;Because\u0026nbsp;seed size is\u0026nbsp;a\u0026nbsp;heterotic\u0026nbsp;trait, the seed coat of F2 seeds obtained by self-pollination of F\u003csub\u003e1\u003c/sub\u003e plants shows heterozygosity, which is a good model for studying heterosis. However, the genetic basis of heterosis in the seed coat has remained largely uncharacterized.\u0026nbsp;To this end, this study attempted to better define gene expression variations between seed coats of maize inbred lines, how these variations affect gene expression levels in F\u003csub\u003e1\u003c/sub\u003e hybrids, and the relative contributions of cis- and trans-variantsto seed coat heterosis.\u0026nbsp;In this study, the self-pollinated kernels of \u0026nbsp;the hybrid Yudan888 and its parental lines were continuously collected from 0 DAP to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP is an important time point for the expansion of the seed coat area in the hybrid, and 8 DAP is the key time point when the seed coat area of the hybrid significantly extends in the parents.Three DAP and 8 DAP are the best time points to study maize seed coat heterosis. According cytological observations at 8 DAP, we speculate that the heterosis of the seed coat may be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results at 3 DAP and 8 DAP indicated that the number of genes expressed nonadditively changed significantly in the hybrid. However, the number of genes expressed additively was not significantly different.\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eOur findings suggest that seed coat heterosis is the result of nonadditive expression caused by dynamic changes in genes at different time points. These dynamic change suggests that a dose effect of gene expression results in heterosis. GO enrichment and KEGG analysis indicated that at 3 DAP, the early stage of seed coat development, genes related to DNA replication and cell cycle regulation contributed significantly to heterosis. At 8 DAP, the hybrid weakened the cell cycle by enhancing metabolite accumulation combined with increased kernel weight and stabilizing heterosis. We expect that this transcriptome comparison will provide a new step towards understanding the causative mechanisms of altered gene expression in hybrids and the molecular mechanisms of maize seed coat heterosis.\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eThe seed coat area of hybrid and its parents was significantly different at 8 DAP.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe seed coat area and kernel weight from 0 to 15 DAP in the hybrid and parental lines were measured continuously (Fig.1A). At 0 DAP, the seed coat area of the hybrid was lower than that of both the mid-parental value (MPV)\u0026nbsp;and the maternal line. At 6 DAP,\u0026nbsp;the hybrid values began to exceed the MPV but were still lower than those of the maternal line. With the expansion of the seed coat, the hybrid value continued to be higher than the MPV at 7 DAP and slightly exceeded the value of the maternal line. At 8 DAP, there was a significant difference\u0026nbsp;(P value\u0026nbsp;= 0.034)\u0026nbsp;between the hybrid and\u0026nbsp;the MPV regarding seed coat area (Fig. 1B).\u0026nbsp;From 9 to 15 DAP, the seed coat area of the hybrid was increasingly different from the MPV and paternal line values (Table S1). In conclusion, 8 DAP is the inflection point and the entry point of heterosis in the seed coat of maize. The kernel weight of the hybrid was significantly different from that of the MPV after 10 DAP (Fig. 1D). These results indicated that the heterosis of the seed coat area occurred prior to changes in seed weight.\u003c/p\u003e\n\u003cp\u003eIn the seed coat area, the relative growth rate of the hybrid and paternal line peaked at 3 DAP, and the maternal line grew most rapidly at 2 DAP\u0026nbsp;(Fig. 1C).\u0026nbsp;In terms of the relative increasing rate of kernel weight,\u0026nbsp;the value of the hybrid was the highest at 3 DAP, while it was highest at 2 DAP for the maternal line and at 4 DAP for the paternal line\u0026nbsp;(Fig. 1E).\u0026nbsp;In general,\u0026nbsp;seed coat area and kernel weight showed a large increase at 2-4 DAP. More importantly, 3 DAP was found to be a critical time point for hybrid seed coat development.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRelative cell numbers in the seed coat\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo understand whether the increased seed coat area in the hybrid was due to changes in cell number or cell size, we counted the number of cells in the seed coat per unit area (0.1 mm\u003csup\u003e2\u003c/sup\u003e) (Fig. 2). The middle longitudinal section of seed embryos collected at 8 DAP was obtained (Fig. 2A-C). The number of cells per unit area for the paternal line (319 cells per\u0026nbsp;0.1 mm\u003csup\u003e2\u003c/sup\u003e) was higher than that for the hybrid (149.7 cells per\u0026nbsp;0.1 mm\u003csup\u003e2\u003c/sup\u003e), the maternal line (138 cells per\u0026nbsp;0.1 mm\u003csup\u003e2\u003c/sup\u003e) and the MPV (228.5 cells per\u0026nbsp;0.1 mm\u003csup\u003e2\u003c/sup\u003e) at 8 DAP (Fig. 2D). Although there was no difference in the number of cells per unit area between the hybrid and maternal line at 8 DAP, there was a significant difference in seed coat surface area (Fig. 2D, E). Thus, compared with the maternal line, the cell number of the hybrid was determined by its larger seed coat area. The differences in cell number and seed coat area between the hybrid and MPV were significant (Fig. 2D), suggesting that the difference in seed coat size between the hybrid and parents might be the result of the coordination between cell number and cell size. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAnalysis of expressed gene in seed coat\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo investigate the expression of genes before and after the time points at which the seed coat of the hybrid line and its parents showed significant differences, we performed RNA-seq on 18 samples at 6 time points (2-4 DAP and 8-10 DAP) from the hybrid and parental lines, with 3 biological replicates for each sample. In total, 54 libraries and 2.4 billion high-quality reads were generated. The reads were then mapped to the maize B73 reference genome (RefGen_V4\u0026nbsp;[26]) by Hisat\u0026nbsp;[27].\u0026nbsp;Approximately (87%) of reads were uniquely mapped (Table S2) and used to calculate the normalized gene expression level as fragments per kb of transcript per million mapped reads (FPKM). There was a high correlation among the three biological replicates at each time point (average R\u003csup\u003e2\u003c/sup\u003e=0.93) (Table S3).\u0026nbsp;All the above results indicated that the sampling quality of this survey was high, and that the RNA-seq sequencing data were accurate and reliable.\u0026nbsp;Thus,\u0026nbsp;we took the average FPKM value of the three replicates as the expression level for the sample at each time point. The Principal Component Analysis (PCA) revealed that the 18\u0026nbsp;samples at six time points from three genotypes were assigned to two stages, and that 2-4 DAP were clearly distinguished from 8-10 DAP (Fig. S1). Moreover, the RNA-seq data of 2-4 DAP were overlapping (Fig. S1).\u003c/p\u003e\n\u003cp\u003eThere were 39005 annotated genes in the B73 AGP_v4 genome[26].\u0026nbsp;The transcriptomes of both the hybrid and parents exhibited very similar distributions in the number of expressed genes at each time point.\u0026nbsp;In detail,\u0026nbsp;the total numbers of genes expressed by the hybrid, maternal and paternal lines at six time points were 23516, 23596 and 23247, respectively. A total of 17691, 18328 and 17604 genes overlapped at these 6 time points, accounting for 75.23%, 77.67% and 75.73% of the total expressed genes (Fig. 3A-C), respectively. Taken together, the results of overlapping expressed genes in the hybrid and its parental lines showed that the growth and development of the seed coat was mainly due to\u0026nbsp;gene expression level rather than the change in gene expression. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDifferentially expressed gene identification\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIt was suggested that the hybrid showed significant differences from its parents at 8 DAP, while the relative growth rate of the hybrid seed coat area and kernel weight peaked at 3 DAP. Therefore, we focused on RNA-seq data at 3 DAP and 8 DAP for subsequent analysis. In DEG analysis, the gene expression levels in the hybrid were compared one-by-one with its parental lines, with a difference in the fold change (FC) \u0026ge; 2 and (FC) \u0026le; -2 denoted as upregulated expression and downregulated expression, respectively. Some DEGs were highly or weakly expressed in either the hybrid or parental lines. At 3 DAP, a total of 3915 DEGs varied between the hybrid and maternal line (2134 (54.50%) upregulated and 1781 (45.50%) downregulated), and 2647 DEGs varied between the hybrid and paternal line (1248 (47.15%) upregulated and 1399 (52.85%) downregulated). Compared with that at 3 DAP, the proportion of DEGs upregulated at 8 DAP increased sharply. A total of 2182 DEGs varied between the hybrid and maternal line (1499 (68.70%) upregulated and 683 (31.30%) downregulated), and 2724 DEGs varied between the hybrid and paternal line (1847 (67.80%) upregulated and 877 (32.20%) downregulated) (Fig. 4A). The differences in the upregulation and/or downregulation of gene expression between the hybrid and its parents varied across comparisons at 3 DAP and 8 DAP. Some DEGs had genotype-specific expression and were thus identified as genotype-specific unigenes (Fig. 4B). Venn diagram analysis indicated that 1156 genes were DEGs (highlighted as black bold font) between the hybrid and its parents at both time points (Fig. 4B). Out of these 1156 genes, 834 DEGs (sum of bold font numerals) were commonly observed between M vs. H and P vs. H at 3 DAP. However, 430 DEGs (sum of bold font numerals) were commonly observed between M vs. H and P vs. H at 8 DAP (Fig. 4B). In the hybrids at both time points, 108 DEGs were found to be commonly differentially expressed (Fig. 4B).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eExpression patterns of DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo explore and categorize the trends for changes in expression, DEGs in 3 DAP and 8 DAP samples were divided into 13 possible hybrid and parental expression patterns (Table 1). Gene expression patterns were divided into additive expression and nonadditive expression.\u0026nbsp;We found no significant changes in the number of additively expressed DEGs between the two time points.\u0026nbsp;The expression pattern of most genes in the hybrid at the two time points was additive (85.18% at 3 DAP and 86.45% at 8 DAP).\u0026nbsp;However, the effects of nonadditive DEGs were quite different, especially the dominantly expressed DEGs (classes 5, 7, 10 and 12). That is, the gene expression in the hybrid was similar to that of one parent, which accounted for more than 13% of all DEGs at both time points.\u0026nbsp;At 3 DAP, the number of dominantly expressed genes with higher paternal expression accounted for the highest proportion in the nonadditive classes, while at 8 DAP, the number of dominantly expressed genes with higher maternal expression was the most abundant.\u0026nbsp;The number of overdominantly expressed genes with transitional downregulation (0.35% at 3 DAP and 0.14% at 8 DAP) was relatively higher than that with transitional upregulation (0.10% at 3 DAP and 0.08% at 8 DAP).\u003c/p\u003e\n\u003cp\u003eWe summarized 13 possible hybrid and parental expression patterns into mid-parent, high-parent and low-parent categories to analyse the expression changes at both time points. By comparing the number of DEGs in three expression patterns at both time points, it was found that the number of low-parent DEGs changed the most, followed by the number of high-parent DEGs, while the number of mid-parents DEGs changed the least (Table S4). According to the analysis of the fold changes at both time points, it was suggested that the change in the number of low-parent and high-parent DEGs may play a role in hybrid seed coat heterosis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunctional classification by Gene Ontology (GO) and\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eKyoto Encyclopedia of Genes and Genomes (KEGG)\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003epathway enrichment analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe DEGs with high and low parental expression at 3 DAP and 8 DAP were annotated by GO. The enrichment of GO terms with low-parental expression at 3 DAP showed that a large number of DEG categories related to stress were significantly enriched in biological processes, followed by cellular components, and categories related to the plasma membrane (Fig. 5A). The highly expressed DEG GO enrichment results at 3 DAP suggested that categories related to chromosome assembly and processing were enriched in cellular components, and categories related to the cell cycle and DNA replication were enriched in biological processes (Fig. 5B). The GO enrichment results for low-parental expression at 8 DAP demonstrated that the categories related to cytoskeleton and microtubule were enriched in cellular component, and the categories enriched in biological process were mainly related to cell cycle, nuclear division and microtubule (Fig. 5C). The GO enrichment results for highly expressed DEGs at 8 DAP indicated that categories related to metabolism were enriched (Fig. 5D). \u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe KEGG pathway enrichment analysis was carried out to further understand the biological function of genes and their interactions. We performed KEGG enrichment analysis on DEGs with high and low-parental expression at 3 DAP and 8 DAP, respectively (Fig. 6). Among the 3 DAP low-parental expressed DEGs, metabolic pathways, biosynthesis of secondary metabolites, protein processing in the endoplasmic reticulum and glutathione metabolism were significantly enriched (Fig. 6A). For the 3 DAP high-parental expressed DEGs, the biosynthesis of secondary metals, DNA replication and photosynthesis were significantly enriched (Fig. 6B). Metabolic pathways, the biosynthesis of secondary metabolites and alpha linolenic acid metabolism were significantly enriched in 8 DAP low-parental expressed DEGs (Fig. 6C). In contrast to DEGs with high parental expression at 3 DAP, DNA replication was also enriched in the low-parental expression DEGs at 8 DAP. For the high-parental expression DEGs at 8 DAP, circadian rhythms, the biosynthesis of secondary metals, glycolysis/gluconeogenesis, thiamine metabolism, metabolic pathways, galactose metabolism and other pathways were significantly enriched (Fig. 6D).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWeighted gene co-expression network analysis (WGCNA)\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo identify the specific genes highly related to seed coat area, we constructed co-expression networks by using the transcriptome data from comparisons of the hybrid and its parents at 6 time points, and associated the co-expression module with seed coat area traits. The results showed that these genes should be divided into 41 modules (Fig. S2). By observing the correlation between modules and samples, it was found that the absolute value of the correlation coefficient between the MEyellow (0.86) and MEsienna3 (0.82) modules was the highest (Fig. S2). In other words, these two modules have the highest correlation with seed coat area. In these two modules, the genes with the top 100 weight values were screened to construct the regulation network diagram (Fig. 7, Fig. S3). First, the module MEyellow revealed that \u003cem\u003eZm00001d019306\u003c/em\u003e was a hub gene (Fig. 7). However, in module MEsienna3, the regulatory network was not as simple and clear as that in MEyellow. We chose genes with a degree value of more than 10 as hub genes. Therefore, the hub genes in module MEsienna3 were \u003cem\u003eZm00001d006330\u003c/em\u003e, \u003cem\u003eZm00001d043289\u003c/em\u003e, \u003cem\u003eZm00001d007181\u003c/em\u003e, \u003cem\u003eZm00001d043290\u003c/em\u003e, \u003cem\u003eZm00001d033510\u003c/em\u003e, \u003cem\u003eZm00001d024732\u003c/em\u003e, \u003cem\u003eZm00001d015025\u003c/em\u003e, \u003cem\u003eZm00001d006331\u003c/em\u003e and \u003cem\u003eZm00001d047796\u0026nbsp;\u003c/em\u003e(Fig. S3, Table S5).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eValidation of DEGs by quantitative real-time PCR\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe expression patterns of nine DEGs were further confirmed by quantitative real-time PCR (qRT‒PCR). We compared the transcriptional profiles of each gene with the seed coat RNA samples from the hybrid and its parents at different time points. The data confirmed that the expression patterns of all 9 DEGs were consistent with the expression levels obtained by RNA-seq analysis (Fig. 8, Fig. S4)\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eHeterosis has been known for centuries\u0026nbsp;[28].\u0026nbsp;Maize seeds are characterized as having heterosis, and seed size is determined by the maize seed coat area to a certain extent. The maize seed coat develops from the ovary wall and carries a maternal genotype. Due to its genetic stability and fixed heterozygosity, the seed coat of F2 seeds obtained by self-pollination of F\u003csub\u003e1\u003c/sub\u003e plants carries a heterozygous genotype, which makes it a good model for studying heterosis. Therefore, exploring the seed coat area of seeds obtained by self-pollination of maize F\u003csub\u003e1\u003c/sub\u003e plants and their self-crossed parents can deepen our understanding of maize seed coat heterosis, and in turn seed size and yield heterosis.\u003c/p\u003e\n\u003cp\u003eThe initiation of seed coat development depends on the fusion of the central cell with one of the sperm cells\u0026nbsp;[29].\u0026nbsp;In our study, the relative growth rate of the seed coat area of the hybrid and paternal line peaked at 3 DAP,\u0026nbsp;while the maternal line grew most rapidly with the peak at 2 DAP. This phenomenon may result from the fact that the synthesis of auxin production in the endosperm after fertilization removes the PRC2 block on seed coat development, to promote the coordinated development of the seed coat and endosperm\u0026nbsp;[30]. There was a significant difference between the hybrid and mid-parental value for seed coat area at 8 DAP,\u0026nbsp;suggesting\u0026nbsp;that heterosis of the seed coat area might\u0026nbsp;be\u0026nbsp;the result of the coordination between cell number and cell size (Fig. 2). These results were consistent with the GO enrichment results from the transcriptome comparison between the hybrid and parents in the early stage of seed coat growth and development, suggesting that heterosis of the seed coat area resulted from the enhanced cell cycle and cell division activity in the early stage. Previous studies have reported that \u0026nbsp;plant growth heterosis basically results from \u0026nbsp;changes in cell number, similar to most plant characteristics\u0026nbsp;[31]. However, the elongation of seed coat cells also affects seed size in Arabidopsis\u0026nbsp;[32]. Our results showed that the maize seed coat developmental program of the hybrid was not dramatically altered, instead,\u0026nbsp;cell proliferation and cell elongation changed.\u003c/p\u003e\n\u003cp\u003eTranscriptome\u0026nbsp;analysis is a well-established method to identify differential expression levels and regulatory mechanisms of genes from different samples at the transcriptional level\u0026nbsp;[33]. Previous studies in plants have identified heterotic genes associated with many traits between hybrids and parents and have proved that gene dose effects play a role in heterosis\u0026nbsp;[34-36]. RNA-seq was performed on the seed coats of the 2-4 DAP and 8-10 DAP hybrid and parental lines. The expression results indicated that the heterosis of seed coat extension was mainly due to the dosage effect of gene expression rather than none-or-all expression.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe differential expression pattern of DEGs between the hybrid and its parents plays a significant role in heterosis\u0026nbsp;[8, 37]. In hybrids, where 2 different alleles of a gene are combined, the allelic allelic expression may deviate from that of either parent or mid-parent predictions\u0026nbsp;[38]. Concerning the relative gene expression level within a hybrid\u0026ndash;parent triad, there are usually two scenarios about the relative gene expression level between the triples of hybrid parents. In the first case, the gene expression of the hybrid showed a cumulative pattern, which was contributed by each allele of parents. This pattern was additive and mainly regulated by cis action. In another case, the expression deviates from the mid-parental level, and other trans acting factors may lead to the change in expression for the corresponding allele in the hybrid\u0026nbsp;[38]. In the present study, the majority (85.18% at 3 DAP and 86.45% at 8 DAP) of genes in the hybrid were additively expressed, suggesting that cis-acting elements played a major role in the regulation of gene expression. The nonadditive genes expressed throughout our entire dataset accounted for less than 15% at each time point but accounted for 54.35\u0026ndash;56.24% of the DEGs. This result indicates that the difference in gene expression between the hybrid and parental lines was mainly due to nonadditive effects. Analysis of the expression of nonadditive genes indicated that the allelic expression pattern of the hybrid may not be a simple combination of alleles from both parents, but may be regulated by other genes or trans-action factors. By comparing the expression patterns of hybrids in 3 DAP and 8 DAP, the number and type of additive expressed genes changed slightly (19714 or 85.18% at 3 DAP to 19460 or 86.45% at 8 DAP DEGs), but the number and type of highly-expressed and weakly-expressed genes at 3 DAP and 8 DAP varied greatly (Table S4,\u0026nbsp;Fig. 5). In conclusion, it was suggested that the expression of nonadditive genes in hybrids may be the main cause of heterosis in hybrid seed coats.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThree DAP and 8 DAP are the key stages of maize endosperm cellularization and differentiation, respectively\u0026nbsp;[20, 21]. In this study, the results of GO enrichment analysis and KEGG pathway analysis of different expression patterns for the seed coat in hybrid and parental lines at 3 DAP and 8 DAP matched the endosperm development process. In other words, at 3 DAP, compared with the parents, the hybrid increased its seed coat size mainly by enhancing cell cycle activity, weakening the response to stress and codeveloping with the endosperm (Fig. 5, A-B). However, at 8 DAP, the hybrid seed coat showed a reduced rate of cell division, instead by strengthening the metabolic process. At this stage, the seed coat changed mainly in coordination with the endosperm to enhance the synthesis and accumulation of nutrients in the seed, to enrich the grain and increase the kernel weight (Fig.5, C-D). This result was supported by the significant difference in kernel weight between the hybrid and MPV at 10 DAP (Fig. 1, D). The KEGG analysis results of genes with high parental expression at 8 DAP showed that, pathways related to circadian rhythm and metabolite accumulation were significantly enriched. Changing the circadian rhythm can improve the viability and biomass of hybrids\u0026nbsp;[39]. In addition, evidence of changes in metabolic profiles has been recorded in hybrids\u0026nbsp;[40, 41]. Moreover, seed development is accompanied by metabolic activities for the synthesis and accumulation of stored products, including protein and carbohydrates\u0026nbsp;[42]. Therefore, in addition to genetic control, the final size of seeds may also be affected by metabolic activities\u0026nbsp;[43].\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eZm00001d018415\u003c/em\u003e and \u003cem\u003eZm00001d051995\u003c/em\u003e,\u0026nbsp;two genes encoding proliferating cell nuclear antigen（PCNA）protein, showed high parental expression patterns in the GO and KEGG analyses at 3 DAP. PCNA is a key protein in the mitotic DNA replication mechanism of all eukaryotes\u0026nbsp;[44]. It acts as a DNA polymerase \u0026delta; cofactors and is thus necessary for the synthesis of new DNA strands by forming a homotrimeric complex\u0026nbsp;[28, 45]\u0026nbsp;that binds to various cyclin-CDK complexes\u0026nbsp;[46-48]. In addition, as a sliding clamp, PCNA also plays a role in regulating DNA metabolism, processing DNA damage and controlling the cell cycle\u0026nbsp;[49].\u0026nbsp;Therefore, the high parental expression of the PCNA gene and its related proteins at 3 DAP may play a role in enhancing DNA replication, repair and metabolism in the early stage of seed coat development, to increase the cell cycle rate and cell number in the hybrid. Moreover, the excessive transcript abundance of these genes in hybrids in the early stage of seed coat development may eventually help to maintain heterosis through dominant gene effects to a great extent.\u0026nbsp;\u003cem\u003eZm00001d039131\u003c/em\u003e, encodes ADP-glucose pyrophosphorylase, which catalyses the key steps of the starch synthesis pathway. This gene was enriched in GO terms and KEGG pathways related to metabolism and the biosynthesis of secondary metabolites at 8 DAP when high parental expression was present.\u0026nbsp;This gene has been reported to be associated with seed development and filling and kernel weight\u0026nbsp;[50-52]. Eight to 12 DAP is the stage of endosperm cell differentiation, during which nutrients and storage materials accumulate alongside cell proliferation and enlargement\u0026nbsp;[21].\u0026nbsp;Seed development is a coordinated process among the seed coat, embryo and endosperm. ADP-glucose pyrophosphorylase was highly parentally expressed at 8 DAP in the hybrid, while the cell cycle-related genes showed a low-parental expression level. This phenomenon may be the result of the rapid development and nutrient filling of the endosperm through enhanced accumulation of metabolites in the hybrid. In other words, during this period, the gene expression pattern of the hybrid changed, that is, cell cycle activity was weakened and metabolic processes were enhanced, to stabilize heterosis by cooperating with the development and filling processes in the endosperm.\u003c/p\u003e\n\u003cp\u003eWGCNA can be used to analyse the relationship between a gene set and sample phenotype, draw the regulatory network between genes in a gene set and identify key regulatory genes\u0026nbsp;[53]. \u003cem\u003eZm00001d019306\u003c/em\u003e, the hub gene of WGCNA in the highest correlation coefficient modules\u0026nbsp;with the highest correlation coefficients, was significantly positively correlated with the seed coat phenotype and encoded, an indole-3-acetyl-leu hydrolase. This enzyme can catalyse the hydrolysis of indole-3-acetyl-L-leucine and release free indole-3-acetic acid (IAA)\u0026nbsp;[54]. Seed development is regulated by phytohormones\u0026nbsp;[43]. Auxin, as a central plant hormone, regulates many biological processes, including the maintenance of meristems, cell division and cell expansion\u0026nbsp;[55-57].\u0026nbsp;In Arabidopsis, indole-3-acetyl-leu hydrolase regulates the auxin response by regulating auxin homeostasis in the endoplasmic reticulum\u0026nbsp;[58]. Auxin is the key regulator of the plant defence response and the main control factor of plant growth and development. The change in gene activity related to the IAA pathway has potential importance in producing a heterosis phenotype\u0026nbsp;[59, 60]. Previous studies showed that in Arabidopsis, the utilization of the IAA pathway caused differences in the control of auxin on cell proliferation by auxin, resulting in different levels of heterosis\u0026nbsp;[61]. Interestingly, these views are consistent with our results. On the one hand, at 3 DAP, genes related to the defence response in hybrids were significantly downregulated. On the other hand, heterosis of the maize seed coat area is determined by the number of cells. In summary, indole-3-acetyl-leu hydrolase may play a significant role in of the IAA pathway and seed coat heterosis. It may regulate the auxin response by regulating auxin homeostasis in different seed coat development processes, to effectively maintain the balance of hybrid immunity and growth and development regarding resource allocation at different developmental time points of the seed coat.\u003c/p\u003e\n\u003cp\u003eFinally, the expression levels of 9 genes identified by RNA-seq analysis were verified by qRT‒PCR, and the expression patterns of all tested genes were consistent with the results of RNA-seq, indicating the reliability of RNA-seq data in Yudan888 and its parents (Fig. 8). In conclusion, the genes identified in our study may provide some insights into broader aspects of seed coat heterosis in maize. However, further validation is needed to confirm the association between gene expression patterns and target agronomic traits in maize.\u0026nbsp;\u003c/p\u003e"},{"header":"Materials And Methods","content":"\u003cp\u003e\u003cstrong\u003eField experiment and\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003ephenotyping\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn summer 2020, maize hybrid Yudan888, a widely cultivated maize variety in the Huang-Huai-Hai area of China, together with its parental lines 15S717 and T4691, was planted on the farm of\u0026nbsp;Henan Agricultural University\u0026nbsp;(Yuanyang, China,\u0026nbsp;113\u0026deg;16\u0026rsquo; E, 35\u0026deg;41\u0026rsquo;\u0026nbsp;N).\u0026nbsp;The planting density was 75,000 plants per hectare, with 3 replicates each containing 600 seedlings.\u0026nbsp;A\u0026nbsp;randomized block design was adopted in the field experiment. All the plants of the hybrid and parental lines were artificially\u0026nbsp;self-crossed.\u0026nbsp;\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSeed coat phenotyping began at 0 DAP and continued daily until 15 DAP, phenotyping was completed between 9 am and 11 am each day.\u0026nbsp;The kernel weight was measured immediately after sampling by an analytical balance (ME204E,\u0026nbsp;METTLER TOLEDO) with an accuracy of 0.0001\u0026nbsp;g. The whole seed coat was peeled off with tweezers and photographed, on a black background. Finally, the seed coat area was calculated by\u0026nbsp;ImageJ\u0026nbsp;software.\u0026nbsp;At least six individual plants per DAP with no less than six kernels per plant were sampled for phenotyping.\u0026nbsp;Univariate ANOVA was performed using IBM SPSS statistics 24 to test the significant differences in the measured the area of the seed coat between the hybrid and parental lines.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSeeds from the centre of three different ears were fixed in FAA buffer (50% ethanol: formaldehyde: acetic acid = 90:5:5 volume) and then vacuumed at 4 \u0026deg;C for further imaging.\u0026nbsp;Kernels from each sample were removed from the FAA solution and rinsed thoroughly in deionized water.\u0026nbsp;Kernels were dissected at the longest part with a scalpel. Each dissected kernel was embedded in paraffin. Sections from kernels ranged from 8 to 10 \u0026mu;m in thickness. The maximum longitudinal cross section of intact kernels was dyed with 1% toluidine blue after dewaxing. Images of the sections were captured with a\u0026nbsp;microscope\u0026nbsp;(Axio Scope A1,\u0026nbsp;Carl Zeiss) fitted with a digital camera (Axiocam 503 colour,\u0026nbsp;Carl Zeiss).\u0026nbsp;ImageJ was used to measure the number of seed coat cells per unit area.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSample collection and\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;RNA\u0026nbsp;\u003c/strong\u003e\u003cstrong\u003eextraction\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eRNA-seq samples were taken during the period when hybrid and parental lines showed significant differences in the seed coat areas.\u0026nbsp;The\u0026nbsp;seed coat samples for RNA-seq were collected by manual dissection with a scalpel.\u0026nbsp;Three biological replicates were set for each sample at each time point. Each replicate sample was collected from at least three individual ears. The\u0026nbsp;samples collected were frozen immediately in liquid nitrogen and\u0026nbsp;stored at -80 \u0026deg;C before RNA extraction.\u0026nbsp;Total RNA was extracted from seed coat samples of the hybrid and its parental lines at 2-4 DAP and 8-10 DAP using TRIzol reagent (Invitrogen, United States). Three biological replicates were used for at each DAP. The\u0026nbsp;quality\u0026nbsp;and concentration of each RNA sample were determined using\u0026nbsp;both gel electrophoresis and a NanoDrop 2000 spectrophotometer\u0026nbsp;(Thermo Fisher Scientific, Wilmington, DE). Only RNAs that met the criterion of an OD260/280 ratio of 1.9-2.1 were stored in a \u0026minus;80 \u0026deg;C freezer for further\u0026nbsp;use.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRNA-seq and data analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA total of 1 \u0026mu;g\u0026nbsp;of RNA was used as the input material for the construction of each RNA library, with the NEBNext\u0026reg; UltraTM RNA Library Prep Kit for Illumina\u0026reg; (NEB, USA). Briefly, mRNA was purified from total RNA by using poly-T oligo-attached magnetic beads. First strand cDNA was synthesized using random hexamer primers and M-MuLV Reverse Transcriptase (with RNase H). After second-strand cDNA synthesis, terminal repair, poly(A) tail addtion sequencing oligonucleotide adaptor ligation, the fragments were purified and subsequently amplified by PCR. The PCR products were assessed with the Agilent Bioanalyzer 2100 system. Finally, the libraries containing 150 bp cDNA inserts in size were generated and sequenced by a commercial service company (Guoke Biotechnology, Beijing, China) on an Illumia Noveseq platform.\u003c/p\u003e\n\u003cp\u003eAll data analyses were based on clean data with high quality [62]. Feature Counts v1.5.0-p3 was used to count the read numbers mapped to each gene (ftp://ftp.ensemblgenomes.org/pub/release-40/plants/fasta/zea_mays/dna/). Multiple testing with the Benjamini‒Hochberg approach for controlling the false discovery rate (FDR) was taken into account by using an adjusted P value. A corrected P value of 0.05 and an absolute fold change of two times were set as the thresholds for significantly differential expression. The normalized gene expression values with FPKM were used for PCA. PCA was performed in R 4.2.1. The GO enrichment analyses were obtained using agriGO v2.0 (http://systemsbiology.cau.edu.cn/agriGOv2/index.php) [63]. The KEGG annotations were obtained on the KEGG website (https://www.kegg.jp/kegg/) [64]. The GO and KEGG enrichment analyses were performed using TBtools [65].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eClassification of differential gene expression patterns\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBased on previous studies on gene expression patterns\u0026nbsp;[66-70], the FPKM values of genes were compared between the hybrid and its parents. The expression patterns were divided into additive and nonadditive (dominance or overdominance). A relative expression level of a hybrid gene that was located between parental values was classified as additive (regardless of whether there were DEGs between the parents or not). A relative expression level of DEGs similar to that of the dominant parents was classified as dominant. A relative expression level of DEGs that was higher than that of the parents (excessive upregulation) or lower than that of the parents (excessive downregulation) was classified as overdominance.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWeighted gene co-expression network analysis\u003c/strong\u003e\u003cstrong\u003e\u0026nbsp;(WGCNA)\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe assignments of gene co-expression modules using the WGCNA protocol were based on FPKM data (Langfelder and Horvath, 2008). To determine the module related to seed coat surface area, we associated the genes in each module with the seed coat area trait from 2-4 DAP and 8-10 DAP. The genes with an average FPKM＞1 in 54 samples were\u0026nbsp;analyzed. The soft threshold power \u0026beta; for network construction was set to 11. The dynamic tree cutting algorithm with a minimum module size of 50 genes was used for hierarchical cluster cutting; a value of 0.15 was used to merge similar modules. If the p value of the module trait association was less than 0.05, it was defined as a meaningful module. In R 4.0.1, the WGCNA package was used to calculate the correlation weight of each gene and all other genes in the module. The top 100 gene sets in the module were visualized using Cytoscape V3.9.1, and the gene with the highest degree was considered to be a hub gene.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eqRT\u003c/strong\u003e‒\u003cstrong\u003ePCR\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNine differentially expressed genes were selected for qRT‒PCR verification, and the first-strand cDNA was synthesized by reverse transcription using total RNA from the seed coat of the hybrid, the maternal line and the paternal line as templates. The primer 3 online toolbox was used to design gene specific primers for qRT‒PCR. The primers used for quantitative PCR are listed in Supplemental Table S6. The reverse transcription and qRT‒PCR reaction system and reaction procedures were in accordance with the instructions for the PrimeScriptTM RT Reagent Kit with gDNA Eraser (TaKaRa) and TB GreenTM Premix Ex TaqTM II (TaKaRa) Kit. Each experiment had 3 biological replicates and 3 technical replicates. The relative expression levels of these genes were analysed by the 2\u003csup\u003e-\u003c/sup\u003e\u003csup\u003e\u0026Delta;\u0026Delta;\u003c/sup\u003e\u003csup\u003eCt\u0026nbsp;\u003c/sup\u003emethod with \u003cem\u003eZm-actin-1\u0026nbsp;\u003c/em\u003eserving as the internal reference [71].\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003e\u003cstrong\u003eDAP:\u0026nbsp;\u003c/strong\u003edays after pollination\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMPV:\u0026nbsp;\u003c/strong\u003emid-parental value\u003c/p\u003e\n\u003cp\u003ePCA: Principal Component Analysis\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDEGs:\u003c/strong\u003e differentially expressed genes\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGO:\u003c/strong\u003e Gene Ontology\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eKEGG:\u0026nbsp;\u003c/strong\u003eKyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWGCNA:\u0026nbsp;\u003c/strong\u003eWeighted gene co-expression network analysis\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eqRT‒PCR:\u0026nbsp;\u003c/strong\u003eReal-time quantitative PCR\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll experimental studies on plants complied with relevant institutional, national, and international guidelines and legislation.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe sequence datasets in fastq format of the current study are available in the NCBI Sequence Read Archive (SRA) database under Bioproject PRJNA828527 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA828527).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by a grant from the\u0026nbsp;National Natural Foundation of China (31971961, 31871641, 32060460), the Science and Technique Foundation of Henan Province of China (202102110036) and the Basic Research Plan of Guizhou Province of China (zk [2022] ‒ 236).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDD, ZHC and JHT designed the research. JL and LFW performed the experiments. JL, LFW and JW analysed the data. JL wrote the manuscript, and DD, ZHC and JHT revised the manuscript. All authors approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors would like to thank Huili Yang, Yanan Lin and Guoqiang Xu for assisting with this article.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; information\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e1\u003c/sup\u003eCollege of Agronomy, Hunan Agricultural University, Changsha 410128, China\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e2\u003c/sup\u003eNational Key Laboratory of Wheat and Maize Crop Science; Collaborative Innovation Center of Henan Grain Crops; College of Agronomy; Henan Agricultural University, Zhengzhou 450002, China\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e3\u003c/sup\u003e Institute of Crop Germplasm Resources, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e4\u003c/sup\u003e Hebi Academy of Agricultural Sciences, Hebi 458030, China\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e5\u003c/sup\u003eThe Shennong laboratory, Zhengzhou 450002, China\u003c/p\u003e\n\u003cp\u003e\u003csup\u003e6\u003c/sup\u003eInstitute of Upland Food Crops, Guizhou Academy of Agricultural Sciences, Guiyang 550006, China\u003c/p\u003e\n\u003cp\u003eJuan Li \u003csup\u003e1,2,3#\u003c/sup\u003e, Liangfa Wang\u003csup\u003e\u0026nbsp;1,2,4#\u003c/sup\u003e, Jiong Wan\u003csup\u003e\u0026nbsp;2\u003c/sup\u003e, Kuntai Dang\u003csup\u003e2\u003c/sup\u003e, Yuan Lin\u003csup\u003e\u0026nbsp;4\u003c/sup\u003e, Shujun Meng\u003csup\u003e\u0026nbsp;2\u003c/sup\u003e, Xiaoqian Qiu\u003csup\u003e\u0026nbsp;2\u003c/sup\u003e, Qiyue Wang\u003csup\u003e\u0026nbsp;4\u003c/sup\u003e, Jiawen Zhao\u003csup\u003e\u0026nbsp;2\u003c/sup\u003e, Liqin Mu\u003csup\u003e\u0026nbsp;2\u003c/sup\u003e, Hongbing Luo\u003csup\u003e\u0026nbsp;1\u003c/sup\u003e, Dong Ding\u003csup\u003e\u0026nbsp;2*\u003c/sup\u003e, Zehui Chen\u003csup\u003e\u0026nbsp;1,6*\u003c/sup\u003e, and Jihua Tang\u003csup\u003e\u0026nbsp;2,5*\u003c/sup\u003e\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eHochholdinger F, Hoecker N: \u003cstrong\u003eTowards the molecular basis of heterosis\u003c/strong\u003e. \u003cem\u003eTrends Plant Sci \u003c/em\u003e2007, \u003cstrong\u003e12\u003c/strong\u003e(9):427-432.\u003c/li\u003e\n\u003cli\u003eBirchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA: \u003cstrong\u003eHeterosis\u003c/strong\u003e. \u003cem\u003ePlant Cell \u003c/em\u003e2010, \u003cstrong\u003e22\u003c/strong\u003e(7):2105-2112.\u003c/li\u003e\n\u003cli\u003eHubner N, Wallace CA, Zimdahl H, Petretto E, Schulz H, Maciver F, Mueller M, Hummel O, Monti J, Zidek V\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eIntegrated transcriptional profiling and linkage analysis for identification of genes underlying disease\u003c/strong\u003e. \u003cem\u003eNat Genet \u003c/em\u003e2005, \u003cstrong\u003e37\u003c/strong\u003e(3):243-253.\u003c/li\u003e\n\u003cli\u003eSong RT, Messing J: \u003cstrong\u003eGene expression of a gene family in maize based on noncollinear haplotypes\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2003, \u003cstrong\u003e100\u003c/strong\u003e(15):9055-9060.\u003c/li\u003e\n\u003cli\u003eLippman ZB, Zamir D: \u003cstrong\u003eHeterosis: revisiting the magic\u003c/strong\u003e. \u003cem\u003eTrends Genet \u003c/em\u003e2007, \u003cstrong\u003e23\u003c/strong\u003e(2):60-66.\u003c/li\u003e\n\u003cli\u003eCharlesworth D, Willis JH: \u003cstrong\u003eThe genetics of inbreeding depression\u003c/strong\u003e. \u003cem\u003eNat Rev Genet \u003c/em\u003e2009, \u003cstrong\u003e10\u003c/strong\u003e(11):783-796.\u003c/li\u003e\n\u003cli\u003eWei G, Tao Y, Liu G, Chen C, Luo R, Xia H, Gan Q, Zeng H, Lu Z, Han Y\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eA transcriptomic analysis of superhybrid rice LYP9 and its parents\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2009, \u003cstrong\u003e106\u003c/strong\u003e(19):7695-7701.\u003c/li\u003e\n\u003cli\u003eHe GM, Zhu XP, Elling AA, Chen LB, Wang XF, Guo L, Liang MZ, He H, Zhang HY, Chen FF\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eGlobal Epigenetic and Transcriptional Trends among Two Rice Subspecies and Their Reciprocal Hybrids\u003c/strong\u003e. \u003cem\u003ePlant Cell \u003c/em\u003e2010, \u003cstrong\u003e22\u003c/strong\u003e(1):17-33.\u003c/li\u003e\n\u003cli\u003eGroszmann M, Greaves IK, Albertyn ZI, Scofield GN, Peacock WJ, Dennis ES: \u003cstrong\u003eChanges in 24-nt siRNA levels in Arabidopsis hybrids suggest an epigenetic contribution to hybrid vigor\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2011, \u003cstrong\u003e108\u003c/strong\u003e(6):2617-2622.\u003c/li\u003e\n\u003cli\u003eSu N, Sullivan JA, Deng XW: \u003cstrong\u003eModulation of F1 hybrid stature without altering parent plants through trans-activated expression of a mutated rice GAI homologue\u003c/strong\u003e. \u003cem\u003ePlant Biotechnol J \u003c/em\u003e2005, \u003cstrong\u003e3\u003c/strong\u003e(2):157-164.\u003c/li\u003e\n\u003cli\u003eXue WY, Xing YZ, Weng XY, Zhao Y, Tang WJ, Wang L, Zhou HJ, Yu SB, Xu CG, Li XH\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eNatural variation in Ghd7 is an important regulator of heading date and yield potential in rice\u003c/strong\u003e. \u003cem\u003eNat Genet \u003c/em\u003e2008, \u003cstrong\u003e40\u003c/strong\u003e(6):761-767.\u003c/li\u003e\n\u003cli\u003eKrieger U, Lippman ZB, Zamir D: \u003cstrong\u003eThe flowering gene SINGLE FLOWER TRUSS drives heterosis for yield in tomato\u003c/strong\u003e. \u003cem\u003eNat Genet \u003c/em\u003e2010, \u003cstrong\u003e42\u003c/strong\u003e(5):459-463.\u003c/li\u003e\n\u003cli\u003eHuang X, Yang S, Gong J, Zhao Q, Feng Q, Zhan Q, Zhao Y, Li W, Cheng B, Xia J\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eGenomic architecture of heterosis for yield traits in rice\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2016, \u003cstrong\u003e537\u003c/strong\u003e(7622):629-633.\u003c/li\u003e\n\u003cli\u003eLi X, Li XR, Fridman E, Tesso TT, Yu JM: \u003cstrong\u003eDissecting repulsion linkage in the dwarfing gene Dw3 region for sorghum plant height provides insights into heterosis\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2015, \u003cstrong\u003e112\u003c/strong\u003e(38):11823-11828.\u003c/li\u003e\n\u003cli\u003eZhang X, Hirsch CN, Sekhon RS, de Leon N, Kaeppler SM: \u003cstrong\u003eEvidence for maternal control of seed size in maize from phenotypic and transcriptional analysis\u003c/strong\u003e. \u003cem\u003eJ Exp Bot \u003c/em\u003e2016, \u003cstrong\u003e67\u003c/strong\u003e(6):1907-1917.\u003c/li\u003e\n\u003cli\u003eSundaresan V: \u003cstrong\u003eControl of seed size in plants\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2005, \u003cstrong\u003e102\u003c/strong\u003e(50):17887-17888.\u003c/li\u003e\n\u003cli\u003eAdamski NM, Anastasiou E, Eriksson S, O\u0026apos;Neill CM, Lenhard M: \u003cstrong\u003eLocal maternal control of seed size by KLUH/CYP78A5-dependent growth signaling\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2009, \u003cstrong\u003e106\u003c/strong\u003e(47):20115-20120.\u003c/li\u003e\n\u003cli\u003eVerdier J, Dessaint F, Schneider C, Abirached-Darmency M: \u003cstrong\u003eA combined histology and transcriptome analysis unravels novel questions on Medicago truncatula seed coat\u003c/strong\u003e. \u003cem\u003eJ Exp Bot \u003c/em\u003e2013, \u003cstrong\u003e64\u003c/strong\u003e(2):459-470.\u003c/li\u003e\n\u003cli\u003eBennetzen JL, Hake SC: \u003cstrong\u003eHandbook of maize: its biology\u003c/strong\u003e. \u003cem\u003eSpringer New York \u003c/em\u003e2009, \u003cstrong\u003e10.1007/978-0-387-79418-1\u003c/strong\u003e.\u003c/li\u003e\n\u003cli\u003eSabelli PA, Larkins BA: \u003cstrong\u003eThe contribution of cell cycle regulation to endosperm development\u003c/strong\u003e. \u003cem\u003eSex Plant Reprod \u003c/em\u003e2009, \u003cstrong\u003e22\u003c/strong\u003e(4):207-219.\u003c/li\u003e\n\u003cli\u003eOttaviano E, Petroni D, Pe ME: \u003cstrong\u003eGametophytic Expression of Genes-Controlling Endosperm Development in Maize\u003c/strong\u003e. \u003cem\u003eTheor Appl Genet \u003c/em\u003e1988, \u003cstrong\u003e75\u003c/strong\u003e(2):252-258.\u003c/li\u003e\n\u003cli\u003eHehenberger E, Kradolfer D, Kohler C: \u003cstrong\u003eEndosperm cellularization defines an important developmental transition for embryo development\u003c/strong\u003e. \u003cem\u003eDevelopment \u003c/em\u003e2012, \u003cstrong\u003e139\u003c/strong\u003e(11):2031-2039.\u003c/li\u003e\n\u003cli\u003eNowack MK, Shirzadi R, Dissmeyer N, Dolf A, Endl E, Grini PE, Schnittger A: \u003cstrong\u003eBypassing genomic imprinting allows seed development\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2007, \u003cstrong\u003e447\u003c/strong\u003e(7142):312-315.\u003c/li\u003e\n\u003cli\u003eThorne JH: \u003cstrong\u003ePhloem Unloading of C-Assimilates and N-Assimilates in Developing Seeds\u003c/strong\u003e. \u003cem\u003eAnnu Rev Plant Physiol Plant Mol Biol \u003c/em\u003e1985, \u003cstrong\u003e36\u003c/strong\u003e:317-343.\u003c/li\u003e\n\u003cli\u003eWestern TL, Skinner DJ, Haughn GW: \u003cstrong\u003eDifferentiation of mucilage secretory cells of the Arabidopsis seed coat\u003c/strong\u003e. \u003cem\u003ePlant Physiol \u003c/em\u003e2000, \u003cstrong\u003e122\u003c/strong\u003e(2):345-356.\u003c/li\u003e\n\u003cli\u003eJiao Y, Peluso P, Shi J, Liang T, Stitzer MC, Wang B, Campbell MS, Stein JC, Wei X, Chin CS\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eImproved maize reference genome with single-molecule technologies\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2017, \u003cstrong\u003e546\u003c/strong\u003e(7659):524-527.\u003c/li\u003e\n\u003cli\u003eKim D, Langmead B, Salzberg SL: \u003cstrong\u003eHISAT: a fast spliced aligner with low memory requirements\u003c/strong\u003e. \u003cem\u003eNat Methods \u003c/em\u003e2015, \u003cstrong\u003e12\u003c/strong\u003e(4):357-360.\u003c/li\u003e\n\u003cli\u003eDarwin C: \u003cstrong\u003eThe effects of cross and self fertilisation in the vegetable kingdom\u003c/strong\u003e. New York,: D. Appleton and company; 1902.\u003c/li\u003e\n\u003cli\u003eRoszak P, Kohler C: \u003cstrong\u003ePolycomb group proteins are required to couple seed coat initiation to fertilization\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2011, \u003cstrong\u003e108\u003c/strong\u003e(51):20826-20831.\u003c/li\u003e\n\u003cli\u003eFigueiredo DD, Batista RA, Roszak PJ, Kohler C: \u003cstrong\u003eAuxin production couples endosperm development to fertilization\u003c/strong\u003e. \u003cem\u003eNat Plants \u003c/em\u003e2015, \u003cstrong\u003e1\u003c/strong\u003e(12).\u003c/li\u003e\n\u003cli\u003eEast EM: \u003cstrong\u003eHeterosis\u003c/strong\u003e. \u003cem\u003eGenetics \u003c/em\u003e1936, \u003cstrong\u003e21\u003c/strong\u003e(4):375-397.\u003c/li\u003e\n\u003cli\u003eGarcia D, Fitz Gerald JN, Berger F: \u003cstrong\u003eMaternal control of integument cell elongation and zygotic control of endosperm growth are coordinated to determine seed size in Arabidopsis\u003c/strong\u003e. \u003cem\u003ePlant Cell \u003c/em\u003e2005, \u003cstrong\u003e17\u003c/strong\u003e(1):52-60.\u003c/li\u003e\n\u003cli\u003eWang Z, Gerstein M, Snyder M: \u003cstrong\u003eRNA-Seq: a revolutionary tool for transcriptomics\u003c/strong\u003e. \u003cem\u003eNat Rev Genet \u003c/em\u003e2009, \u003cstrong\u003e10\u003c/strong\u003e(1):57-63.\u003c/li\u003e\n\u003cli\u003eYao H, Dogra Gray A, Auger DL, Birchler JA: \u003cstrong\u003eGenomic dosage effects on heterosis in triploid maize\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2013, \u003cstrong\u003e110\u003c/strong\u003e(7):2665-2669.\u003c/li\u003e\n\u003cli\u003eBirchler JA, Riddle NC, Auger DL, Veitia RA: \u003cstrong\u003eDosage balance in gene regulation: biological implications\u003c/strong\u003e. \u003cem\u003eTrends Genet \u003c/em\u003e2005, \u003cstrong\u003e21\u003c/strong\u003e(4):219-226.\u003c/li\u003e\n\u003cli\u003eBirchler JA, Johnson AF, Veitia RA: \u003cstrong\u003eKinetics genetics: Incorporating the concept of genomic balance into an understanding of quantitative traits\u003c/strong\u003e. \u003cem\u003ePlant Sci (Amsterdam, Neth) \u003c/em\u003e2016, \u003cstrong\u003e245\u003c/strong\u003e:128-134.\u003c/li\u003e\n\u003cli\u003eZhang HY, He H, Chen LB, Li L, Liang MZ, Wang XF, Liu XG, He GM, Chen RS, Ma LG\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eA genome-wide transcription analysis reveals a close correlation of promoter INDEL polymorphism and heterotic gene expression in rice hybrids\u003c/strong\u003e. \u003cem\u003eMol Plant \u003c/em\u003e2008, \u003cstrong\u003e1\u003c/strong\u003e(5):720-731.\u003c/li\u003e\n\u003cli\u003eBirchler JA, Auger DL, Riddle NC: \u003cstrong\u003eIn search of the molecular basis of heterosis\u003c/strong\u003e. \u003cem\u003ePlant Cell \u003c/em\u003e2003, \u003cstrong\u003e15\u003c/strong\u003e(10):2236-2239.\u003c/li\u003e\n\u003cli\u003eNi Z, Kim ED, Ha M, Lackey E, Liu J, Zhang Y, Sun Q, Chen ZJ: \u003cstrong\u003eAltered circadian rhythms regulate growth vigour in hybrids and allopolyploids\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2009, \u003cstrong\u003e457\u003c/strong\u003e(7227):327-331.\u003c/li\u003e\n\u003cli\u003eGartner T, Steinfath M, Andorf S, Lisec J, Meyer RC, Altmann T, Willmitzer L, Selbig J: \u003cstrong\u003eImproved heterosis prediction by combining information on DNA- and metabolic markers\u003c/strong\u003e. \u003cem\u003ePLoS ONE \u003c/em\u003e2009, \u003cstrong\u003e4\u003c/strong\u003e(4):e5220.\u003c/li\u003e\n\u003cli\u003eFievet JB, Dillmann C, de Vienne D: \u003cstrong\u003eSystemic properties of metabolic networks lead to an epistasis-based model for heterosis\u003c/strong\u003e. \u003cem\u003eTheor Appl Genet \u003c/em\u003e2010, \u003cstrong\u003e120\u003c/strong\u003e(2):463-473.\u003c/li\u003e\n\u003cli\u003eSreenivasulu N, Wobus U: \u003cstrong\u003eSeed-Development Programs: A Systems Biology-Based Comparison Between Dicots and Monocots\u003c/strong\u003e. \u003cem\u003eAnnual Review of Plant Biology, Vol 64 \u003c/em\u003e2013, \u003cstrong\u003e64\u003c/strong\u003e:189-+.\u003c/li\u003e\n\u003cli\u003eBorisjuk L, Rolletschek H, Radchuk R, Weschke W, Wobus U, Weber H: \u003cstrong\u003eSeed development and differentiation: a role for metabolic regulation\u003c/strong\u003e. \u003cem\u003ePlant Biol (Stuttg) \u003c/em\u003e2004, \u003cstrong\u003e6\u003c/strong\u003e(4):375-386.\u003c/li\u003e\n\u003cli\u003eTan CK, Castillo C, So AG, Downey KM: \u003cstrong\u003eAn auxiliary protein for DNA polymerase-delta from fetal calf thymus\u003c/strong\u003e. \u003cem\u003eJ Biol Chem \u003c/em\u003e1986, \u003cstrong\u003e261\u003c/strong\u003e(26):12310-12316.\u003c/li\u003e\n\u003cli\u003ede la Paz Sanchez M, Torres A, Boniotti MB, Gutierrez C, Vazquez-Ramo JM: \u003cstrong\u003ePCNA protein associates to Cdk-A type protein kinases in germinating maize\u003c/strong\u003e. \u003cem\u003ePlant Mol Biol \u003c/em\u003e2002, \u003cstrong\u003e50\u003c/strong\u003e(2):167-175.\u003c/li\u003e\n\u003cli\u003eCelis JE, Madsen P, Celis A, Nielsen HV, Gesser B: \u003cstrong\u003eCyclin (Pcna, Auxiliary Protein of DNA Polymerase-Delta) Is a Central Component of the Pathway(S) Leading To DNA-Replication and Cell-Division\u003c/strong\u003e. \u003cem\u003eFEBS Lett \u003c/em\u003e1987, \u003cstrong\u003e220\u003c/strong\u003e(1):1-7.\u003c/li\u003e\n\u003cli\u003ePrelich G, Tan CK, Kostura M, Mathews MB, So AG, Downey KM, Stillman B: \u003cstrong\u003eFunctional Identity of Proliferating Cell Nuclear Antigen and a DNA Polymerase-Delta Auxiliary Protein\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e1987, \u003cstrong\u003e326\u003c/strong\u003e(6112):517-520.\u003c/li\u003e\n\u003cli\u003eSanchez MD, Gurusinghe SH, Bradford KJ, Vazquez-Ramos JM: \u003cstrong\u003eDifferential response of PCNA and Cdk-A proteins and associated kinase activities to benzyladenine and abscisic acid during maize seed germination\u003c/strong\u003e. \u003cem\u003eJ Exp Bot \u003c/em\u003e2005, \u003cstrong\u003e56\u003c/strong\u003e(412):515-523.\u003c/li\u003e\n\u003cli\u003eStrzalka W, Ziemienowicz A: \u003cstrong\u003eProliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation\u003c/strong\u003e. \u003cem\u003eAnn Bot-London \u003c/em\u003e2011, \u003cstrong\u003e107\u003c/strong\u003e(7):1127-1140.\u003c/li\u003e\n\u003cli\u003eKawagoe Y, Kubo A, Satoh H, Takaiwa F, Nakamura Y: \u003cstrong\u003eRoles of isoamylase and ADP-glucose pyrophosphorylase in starch granule synthesis in rice endosperm\u003c/strong\u003e. \u003cem\u003ePlant J \u003c/em\u003e2005, \u003cstrong\u003e42\u003c/strong\u003e(2):164-174.\u003c/li\u003e\n\u003cli\u003eNa G, Aryal N, Fatihi A, Kang J, Lu C: \u003cstrong\u003eSeed-specific suppression of ADP-glucose pyrophosphorylase in Camelina sativa increases seed size and weight\u003c/strong\u003e. \u003cem\u003eBiotechnol Biofuels \u003c/em\u003e2018, \u003cstrong\u003e11\u003c/strong\u003e:330.\u003c/li\u003e\n\u003cli\u003ePrathap V, Tyagi A: \u003cstrong\u003eCorrelation between expression and activity of ADP glucose pyrophosphorylase and starch synthase and their role in starch accumulation during grain filling under drought stress in rice\u003c/strong\u003e. \u003cem\u003ePlant Physiol Biochem \u003c/em\u003e2020, \u003cstrong\u003e157\u003c/strong\u003e:239-243.\u003c/li\u003e\n\u003cli\u003eLangfelder P, Horvath S: \u003cstrong\u003eWGCNA: an R package for weighted correlation network analysis\u003c/strong\u003e. \u003cem\u003eBMC Bioinformatics \u003c/em\u003e2008, \u003cstrong\u003e9\u003c/strong\u003e:559.\u003c/li\u003e\n\u003cli\u003eBartel B, Leclere S, Magidin M, Zolman BK: \u003cstrong\u003eInputs to the Active Indole-3-Acetic Acid Pool: De Novo Synthesis, Conjugate Hydrolysis, and Indole-3-Butyric Acid b-Oxidation\u003c/strong\u003e. \u003cem\u003eJ Plant Growth Regul \u003c/em\u003e2001, \u003cstrong\u003e20\u003c/strong\u003e(3):198-216.\u003c/li\u003e\n\u003cli\u003eChapman EJ, Estelle M: \u003cstrong\u003eMechanism of auxin-regulated gene expression in plants\u003c/strong\u003e. \u003cem\u003eAnnu Rev Genet \u003c/em\u003e2009, \u003cstrong\u003e43\u003c/strong\u003e:265-285.\u003c/li\u003e\n\u003cli\u003eVanneste S, Friml J: \u003cstrong\u003eAuxin: a trigger for change in plant development\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2009, \u003cstrong\u003e136\u003c/strong\u003e(6):1005-1016.\u003c/li\u003e\n\u003cli\u003eWeijers D, Wagner D: \u003cstrong\u003eTranscriptional Responses to the Auxin Hormone\u003c/strong\u003e. \u003cem\u003eAnnu Rev Plant Biol \u003c/em\u003e2016, \u003cstrong\u003e67\u003c/strong\u003e:539-574.\u003c/li\u003e\n\u003cli\u003eRampey, Rebekah A: \u003cstrong\u003eThe roles of Arabidopsis IAA-amino acid conjugate hydrolases in auxin homeostasis\u003c/strong\u003e.\u003c/li\u003e\n\u003cli\u003eKazan K, Manners JM: \u003cstrong\u003eLinking development to defense: auxin in plant-pathogen interactions\u003c/strong\u003e. \u003cem\u003eTrends Plant Sci \u003c/em\u003e2009, \u003cstrong\u003e14\u003c/strong\u003e(7):373-382.\u003c/li\u003e\n\u003cli\u003eBusov VB, Brunner AM, Strauss SH: \u003cstrong\u003eGenes for control of plant stature and form\u003c/strong\u003e. \u003cem\u003eNew Phytol \u003c/em\u003e2008, \u003cstrong\u003e177\u003c/strong\u003e(3):589-607.\u003c/li\u003e\n\u003cli\u003eGroszmann M, Gonzalez-Bayon R, Lyons RL, Greaves IK, Kazan K, Peacock WJ, Dennis ES: \u003cstrong\u003eHormone-regulated defense and stress response networks contribute to heterosis in Arabidopsis F1 hybrids\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2015, \u003cstrong\u003e112\u003c/strong\u003e(46):E6397-E6406.\u003c/li\u003e\n\u003cli\u003eBolger AM, Lohse M, Usadel B: \u003cstrong\u003eTrimmomatic: a flexible trimmer for Illumina sequence data\u003c/strong\u003e. \u003cem\u003eBioinformatics \u003c/em\u003e2014, \u003cstrong\u003e30\u003c/strong\u003e(15):2114-2120.\u003c/li\u003e\n\u003cli\u003eTian T, Liu Y, Yan H, You Q, Yi X, Du Z, Xu W, Su Z: \u003cstrong\u003eagriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res \u003c/em\u003e2017, \u003cstrong\u003e45\u003c/strong\u003e(W1):W122-W129.\u003c/li\u003e\n\u003cli\u003eKanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M: \u003cstrong\u003eKEGG: integrating viruses and cellular organisms\u003c/strong\u003e. \u003cem\u003eNucleic Acids Res \u003c/em\u003e2021, \u003cstrong\u003e49\u003c/strong\u003e(D1):D545-D551.\u003c/li\u003e\n\u003cli\u003eChen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R: \u003cstrong\u003eTBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data\u003c/strong\u003e. \u003cem\u003eMol Plant \u003c/em\u003e2020, \u003cstrong\u003e13\u003c/strong\u003e(8):1194-1202.\u003c/li\u003e\n\u003cli\u003eYoo MJ, Szadkowski E, Wendel JF: \u003cstrong\u003eHomoeolog expression bias and expression level dominance in allopolyploid cotton\u003c/strong\u003e. \u003cem\u003eHeredity (Edinb) \u003c/em\u003e2013, \u003cstrong\u003e110\u003c/strong\u003e(2):171-180.\u003c/li\u003e\n\u003cli\u003eZhou P, Hirsch CN, Briggs SP, Springer NM: \u003cstrong\u003eDynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development\u003c/strong\u003e. \u003cem\u003eMol Plant \u003c/em\u003e2019, \u003cstrong\u003e12\u003c/strong\u003e(3):410-425.\u003c/li\u003e\n\u003cli\u003eDing HP, Qin C, Luo XR, Li LJ, Chen Z, Liu HJ, Gao J, Lin HJ, Shen YO, Zhao MJ\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eHeterosis in Early Maize Ear Inflorescence Development: A Genome-Wide Transcription Analysis for Two Maize Inbred Lines and Their Hybrid\u003c/strong\u003e. \u003cem\u003eInt J Mol Sci \u003c/em\u003e2014, \u003cstrong\u003e15\u003c/strong\u003e(8):13892-13915.\u003c/li\u003e\n\u003cli\u003eLi AL, Liu DC, Wu J, Zhao XB, Hao M, Geng SF, Yan J, Jiang XX, Zhang LQ, Wu JY\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003emRNA and Small RNA Transcriptomes Reveal Insights into Dynamic Homoeolog Regulation of Allopolyploid Heterosis in Nascent Hexaploid Wheat\u003c/strong\u003e. \u003cem\u003ePlant Cell \u003c/em\u003e2014, \u003cstrong\u003e26\u003c/strong\u003e(5):1878-1900.\u003c/li\u003e\n\u003cli\u003eHowlader J, Robin AHK, Natarajan S, Biswas MK, Sumi KR, Song CY, Park JI, Nou IS: \u003cstrong\u003eTranscriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)\u003c/strong\u003e. \u003cem\u003eSci Rep-Uk \u003c/em\u003e2020, \u003cstrong\u003e10\u003c/strong\u003e(1).\u003c/li\u003e\n\u003cli\u003eLivak KJ, Schmittgen TD: \u003cstrong\u003eAnalysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method\u003c/strong\u003e. \u003cem\u003eMethods \u003c/em\u003e2001, \u003cstrong\u003e25\u003c/strong\u003e(4):402-408.\u003c/li\u003e\n\u003cli\u003eHowlader J, Robin AHK, Natarajan S, Biswas MK, Sumi KR, Song CY, Park JI, Nou IS: \u003cstrong\u003eTranscriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)\u003c/strong\u003e. \u003cem\u003eSci Rep \u003c/em\u003e2020, \u003cstrong\u003e10\u003c/strong\u003e(1):9082.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003eTable 1 is available in the Supplementary Files section.\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-plant-biology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pbio","sideBox":"Learn more about [BMC Plant Biology](http://bmcplantbiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/pbio/default.aspx","title":"BMC Plant Biology","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Maize, Seed coat, Heterosis, RNA-seq. ","lastPublishedDoi":"10.21203/rs.3.rs-1768077/v2","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-1768077/v2","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSeed size is an important factor contributing to maize yield, but its molecular mechanism remains unclear. The seed coat, which serves as one of the three components of the maize grain, determines seed size to a certain extent. The seed coat also shares the maternal genotype and is an ideal material for studying heterosis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResults\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn this study, the self-pollinated seeds of the maize hybrid Yudan888 and its parental lines were continuously collected from 0 day after pollination (DAP) to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP and 8 DAP are the best time points to study maize seed coat heterosis. Cytological observations indicated that maize seed coat heterosis might be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results showed that the nonadditive genes changed significantly between 3 DAP and 8 DAP. However, the number of genes expressed additively was not significantly different. Our findings suggest that seed coat heterosis in hybrid is the result of nonadditive expression caused by dynamic changes in genes at different time points during seed expansion and seed coat development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment indicated that genes related to DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation contributed significantly to hybrid \u0026nbsp;seed coat heterosis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConclusion\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eMaize seed coat phenotyping allowed us to infer that 3 DAP and 8 DAP are important time points in the study of seed coat heterosis. Our findings provide evidence for genes involved in DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation in hybrid with high or low parental expression as major contributors to hybrid seed coat heterosis.\u003c/p\u003e","manuscriptTitle":"Dynamic Patterns of Gene Expression and Regulatory Variation in the Maize Seed Coat","msid":"","msnumber":"","nonDraftVersions":[{"code":2,"date":"2023-02-13 22:55:57","doi":"10.21203/rs.3.rs-1768077/v2","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Accepted","date":"2023-01-19T07:56:55+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2023-01-17T08:39:36+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"44372362-d5b2-4d20-92fc-b2181f92e3f2","date":"2023-01-17T08:36:54+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"e0dfe215-e3cf-4446-a74a-e50c176ec9c9","date":"2023-01-16T23:17:50+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-01-13T15:11:22+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2023-01-12T12:16:03+00:00","index":"","fulltext":""},{"type":"submitted","content":"BMC Plant Biology","date":"2023-01-11T01:16:52+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-plant-biology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pbio","sideBox":"Learn more about [BMC Plant Biology](http://bmcplantbiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/pbio/default.aspx","title":"BMC Plant Biology","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"b0dc1f60-b5ba-4468-be3f-42755ce7a45a","owner":[],"postedDate":"February 13th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2023-10-16T18:53:54+00:00","versionOfRecord":{"articleIdentity":"rs-1768077","link":"https://doi.org/10.1186/s12870-023-04078-1","journal":{"identity":"bmc-plant-biology","isVorOnly":false,"title":"BMC Plant Biology"},"publishedOn":"2023-02-07 18:45:45","publishedOnDateReadable":"February 7th, 2023"},"versionCreatedAt":"2023-02-13 22:55:57","video":"","vorDoi":"10.1186/s12870-023-04078-1","vorDoiUrl":"https://doi.org/10.1186/s12870-023-04078-1","workflowStages":[]},"version":"v2","identity":"rs-1768077","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-1768077","identity":"rs-1768077","version":["v2"]},"buildId":"GqpaHPwrfC8PjnIFayRh5","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}