Widespread Contamination of Arabidopsis Embryo and Endosperm Transcriptome Data Sets

Regulation of Parent-of-Origin Allelic Expression in the Endosperm

Genomic imprinting is an epigenetic phenomenon established in the gametes prior to fertilization that causes differential expression of parental alleles, mainly in the endosperm of flowering plants. The overlap between previously identified panels of imprinted genes is limited. To investigate imprinting, we used high-resolution sequencing data acquired with sequence-capture technology. We present a bioinformatics pipeline to assay parent-of-origin allele-specific expression and report more than 300 loci with parental expression bias in Arabidopsis (Arabidopsis thaliana). In most cases, the level of expression from maternal and paternal alleles was not binary, instead supporting a differential dosage hypothesis for the evolution of imprinting in plants. To address imprinting regulation, we systematically employed mutations in regulative epigenetic pathways suggested to be major players in the process. We established the mechanistic mode of imprinting for more than 50 loci regulated by DNA methylation and Polycomb-dependent histone methylation. However, the imprinting patterns of most genes were not affected by these mechanisms. To this end, we also demonstrated that the RNA-directed DNA methylation pathway alone does not substantially influence imprinting patterns, suggesting that more complex epigenetic pathways regulate most of the identified imprinted genes.

Transcriptome landscape of the early Brassica napus seed

Brassica napus L. (canola) is one of the world's most economically important oilseeds. Despite our growing knowledge of Brassica genetics, we still know little about the genes and gene regulatory networks underlying early seed development. In this work, we use laser microdissection coupled with RNA sequencing to profile gene activity of both the maternal and filial subregions of the globular seed. We find subregions of the chalazal end including the chalazal endosperm, chalazal proliferating tissue, and chalazal seed coat, have unique transcriptome profiles associated with hormone biosynthesis and polysaccharide metabolism. We confirm that the chalazal seed coat is uniquely enriched for sucrose biosynthesis and transport, and that the chalazal endosperm may function as an important regulator of the maternal region through brassinosteroid synthesis. The chalazal proliferating tissue, a poorly understood subregion, was specifically enriched in transcripts associated with megasporogenesis and trehalose biosynthesis, suggesting this ephemeral structure plays an important role in both sporophytic development and carbon nutrient balance, respectively. Finally, compartmentalization of transcription factors and their regulatory circuits has uncovered previously unknown roles for the chalazal pole in early seed development.

Family quarrels in seeds and rapid adaptive evolution in Arabidopsis

Evolutionary conflict can drive rapid adaptive evolution, sometimes called an arms race, because each party needs to respond continually to the adaptations of the other. Evidence for such arms races can sometimes be seen in morphology, in behavior, or in the genes underlying sexual interactions of host-pathogen interactions, but is rarely predicted a priori. Kin selection theory predicts that conflicts of interest should usually be reduced but not eliminated among genetic relatives, but there is little evidence as to whether conflict within families can drive rapid adaptation. Here we test multiple predictions about how conflict over the amount of resources an offspring receives from its parent would drive rapid molecular evolution in seed tissues of the flowering plant Arabidopsis As predicted, there is more adaptive evolution in genes expressed in Arabidopsis seeds than in other specialized organs, more in endosperms and maternal tissues than in embryos, and more in the specific subtissues involved in nutrient transfer. In the absence of credible alternative hypotheses, these results suggest that kin selection and conflict are important in plants, that the conflict includes not just the mother and offspring but also the triploid endosperm, and that, despite the conflict-reducing role of kinship, family members can engage in slow but steady tortoise-like arms races.

The embryonic transcriptome of Arabidopsis thaliana

Arabidopsis embryos possess unique transcriptomes relative to other plant tissues including somatic embryos, and can be partitioned into four transcriptional phases with characteristic biological processes. Cellular differentiation is associated with changes in transcript populations. Accurate quantification of transcriptomes during development can thus provide global insights into differentiation processes including the fundamental specification and differentiation events operating during plant embryogenesis. However, multiple technical challenges have limited the ability to obtain high-quality early embryonic transcriptomes, namely the low amount of RNA obtainable and contamination from surrounding endosperm and seed-coat tissues. We compared the performance of three low-input mRNA sequencing (mRNA-seq) library preparation kits on 0.1 to 5 nanograms (ng) of total RNA isolated from Arabidopsis thaliana (Arabidopsis) embryos and identified a low-cost method with superior performance. This mRNA-seq method was then used to profile the transcriptomes of Arabidopsis embryos across eight developmental stages. By comprehensively comparing embryonic and post-embryonic transcriptomes, we found that embryonic transcriptomes do not resemble any other plant tissue we analyzed. Moreover, transcriptome clustering analyses revealed the presence of four distinct phases of embryogenesis which are enriched in specific biological processes. We also compared zygotic embryo transcriptomes with publicly available somatic embryo transcriptomes. Strikingly, we found little resemblance between zygotic embryos and somatic embryos derived from late-staged zygotic embryos suggesting that somatic and zygotic embryo transcriptomes are distinct from each other. In addition to the biological insights gained from our systematic characterization of the Arabidopsis embryonic transcriptome, we provide a data-rich resource for the community to explore.

Step-by-step protocols for rice gamete isolation

A detailed, step-by-step protocol for isolation of rice gametes for transcriptional profiling, with a general workflow that includes controls for RNA contamination from surrounding cells and tissues is presented. Characterization of the transcriptome and other -omics studies of flowering plant gametes are challenging as a consequence of the small sizes and relative inaccessibility of these cells. Collecting such poorly represented cells is also complicated by potential contamination from surrounding sporophytic, adjacent gametophytic tissues and difficulties in extracting high-quality intact cells. Here we present detailed, step-by-step procedures for collecting intact, unfixed rice (Oryza sativa) egg cells and sperm cells without enzymatic treatments. In addition, we also present a general workflow for assessing sample purity by RT-PCR, using primers specific for marker genes preferentially expressed in surrounding cells and tissues. These protocols should facilitate future studies of genome-scale characterization of gametes in this important model crop.

Isolation of living apical and basal cell lineages of early proembryos for transcriptome analysis

Plant zygotes usually undergo asymmetrical cell division, giving rise to the formation of two daughter cells with distinct developmental cell fate. The small apical cell will develop into the major part of embryo proper, whereas the larger basal cell will divide to form a transient suspensor. Thus, the apical and basal cell lineages are an excellent model to study cell fate determination in relation to zygote polarity. However, the molecular mechanism underlying the differentiation of two distinct cell lineages is not yet understood, possibly due to the technique limitations. Previously, we have established a protocol for isolating apical cell and basal cell for cDNA library construction in tobacco. However, the method for isolating tiny Arabidopsis embryos has long been considered much more difficult. Here, we present a detailed protocol for isolating early Arabidopsis proembryos and separating apical and basal cell lineages of proembryos, which allow us to establish cell lineage-specific transcriptomes of early proembryos.

Endosperm-specific transcriptome analysis by applying the INTACT system

KEY MESSAGE

We report the adaptation of the INTACT method for RNA-sequencing in the endosperm and demonstrate its feasibility for allele-specific expression analysis. Tissue-specific transcriptome analyses provide important insights into the developmental programs of defined cell types. The isolation of nuclei tagged in specific cell types (INTACT) is a versatile method that allows to isolate highly pure nuclei from defined tissue types that can be used for several downstream applications. Here, we describe the adaptation of INTACT from endosperm nuclei for high-throughput RNA-sequencing. By analyzing the ratio of parental reads and tissue-specific gene expression in the endosperm, we could assess the contamination level of our samples. Based on this analysis, we estimate that in most of the samples the contamination level is lower than in previously published datasets. We further show that the nuclear transcriptome and total transcriptome of the endosperm are well correlated. Together, our data show that INTACT of the endosperm is a reliable methodology for endosperm-specific transcriptome analysis that overcomes the limitation of time-consuming manual endosperm dissection that is connected with high levels of maternal tissue contamination. INTACT does not rely on expensive equipment and can be set up in every standard molecular biology laboratory, making it the method of choice for future molecular studies of the endosperm.

Epigenetic signatures associated with imprinted paternally expressed genes in the Arabidopsis endosperm

BACKGROUND

Imprinted genes are epigenetically modified during gametogenesis and maintain the established epigenetic signatures after fertilization, causing parental-specific gene expression.

RESULTS

In this study, we show that imprinted paternally expressed genes (PEGs) in the Arabidopsis endosperm are marked by an epigenetic signature of Polycomb Repressive Complex2 (PRC2)-mediated H3K27me3 together with heterochromatic H3K9me2 and CHG methylation, which specifically mark the silenced maternal alleles of PEGs. The co-occurrence of H3K27me3 and H3K9me2 on defined loci in the endosperm drastically differs from the strict separation of both pathways in vegetative tissues, revealing tissue-specific employment of repressive epigenetic pathways in plants. Based on the presence of this epigenetic signature on maternal alleles, we are able to predict known PEGs at high accuracy and identify several new PEGs that we confirm using INTACT-based transcriptomes generated in this study.

CONCLUSIONS

The presence of the three repressive epigenetic marks, H3K27me3, H3K9me2, and CHG methylation on the maternal alleles in the endosperm serves as a specific epigenetic signature that allows prediction of genes with parental-specific gene expression. Our study reveals that there are substantially more PEGs than previously identified, indicating that paternal-specific gene expression is of higher functional relevance than currently estimated. The combined activity of PRC2-mediated H3K27me3 together with the heterochromatic H3K9me3 has also been reported to silence the maternal Xist locus in mammalian preimplantation embryos, suggesting convergent employment of both pathways during the evolution of genomic imprinting.

Consistent Reanalysis of Genome-wide Imprinting Studies in Plants Using Generalized Linear Models Increases Concordance across Datasets

Genomic imprinting leads to different expression levels of maternally and paternally derived alleles. Over the last years, major progress has been made in identifying novel imprinted candidate genes in plants, owing to affordable next-generation sequencing technologies. However, reports on sequencing the transcriptome of hybrid F1 seed tissues strongly disagree about how many and which genes are imprinted. This raises questions about the relative impact of biological, environmental, technical, and analytic differences or biases. Here, we adopt a statistical approach, frequently used in RNA-seq data analysis, which properly models count overdispersion and considers replicate information of reciprocal crosses. We show that our statistical pipeline outperforms other methods in identifying imprinted genes in simulated and real data. Accordingly, reanalysis of genome-wide imprinting studies in Arabidopsis and maize shows that, at least for Arabidopsis, an increased agreement across datasets could be observed. For maize, however, consistent reanalysis did not yield a larger overlap between the datasets. This suggests that the discrepancy across publications might be partially due to different analysis pipelines but that technical, biological, and environmental factors underlie much of the discrepancy between datasets. Finally, we show that the set of genes that can be characterized regarding allelic bias by all studies with minimal confidence is small (~8,000/27,416 genes for Arabidopsis and ~12,000/39,469 for maize). In conclusion, we propose to use biologically replicated reciprocal crosses, high sequence coverage, and a generalized linear model approach to identify differentially expressed alleles in developing seeds.

New cues for body axis formation in plant embryos

Plant embryogenesis initiates with the fusion of sperm and egg cell, and completes the generation of the basic outline of the future plant. Here, we summarize the recent findings about the signaling cascade triggering the zygotic transcription, and the intracellular events and regulatory factors involved in the formation of the two major body axes. We highlight the lack of systematic de novo transcriptional activation in the zygote, and emphasize the importance of cytoskeletal reorganization to polarize the zygote and control the first asymmetric division that establishes the apical-basal axis. Finally, the limited knowledge of mechanisms that control the cell divisions separating the inner and outer cell layers is summarized and we propose approaches to enhance our understanding of basic principles of plant embryogenesis.

Maternal small RNAs mediate spatial-temporal regulation of gene expression, imprinting, and seed development in Arabidopsis

Arabidopsis seed development involves maternal small interfering RNAs (siRNAs) that induce RNA-directed DNA methylation (RdDM) through the NRPD1-mediated pathway. To investigate their biological functions, we characterized siRNAs in the endosperm and seed coat that were separated by laser-capture microdissection (LCM) in reciprocal genetic crosses with an nrpd1 mutant. We also monitored the spatial-temporal activity of the NRPD1-mediated pathway on seed development using the AGO4:GFP::AGO4 (promoter:GFP::protein) reporter and promoter:GUS sensors of siRNA-mediated silencing. From these approaches, we identified four distinct groups of siRNA loci dependent on or independent of the maternal NRPD1 allele in the endosperm or seed coat. A group of maternally expressed NRPD1-siRNA loci targets endosperm-preferred genes, including those encoding AGAMOUS-LIKE (AGL) transcription factors. Using translational promoter:AGL::GUS constructs as sensors, we demonstrate that spatial and temporal expression patterns of these genes in the endosperm are regulated by the NRPD1-mediated pathway irrespective of complete silencing (AGL91) or incomplete silencing (AGL40) of these target genes. Moreover, altered expression of these siRNA-targeted genes affects seed size. We propose that the corresponding maternal siRNAs could account for parent-of-origin effects on the endosperm in interploidy and hybrid crosses. These analyses reconcile previous studies on siRNAs and imprinted gene expression during seed development.

Genome-Wide Analysis of Parent-of-Origin Allelic Expression in Endosperms of Brassicaceae Species, Brassica rapa

Uniparental gene expression, observed in both animals and plants, is termed genomic imprinting. Genomic imprinting is a well-known epigenetic phenomenon regulated through epigenetic modifications such as DNA methylation and histone modifications. Recent genome-wide studies of endosperm transcription have revealed the rapid change of imprinted genes between species, suggesting the flexibility of this phenomenon. Although the functional significance and evolutionary trends of imprinted genes are still obscure, it can be clarified by inter-species comparisons. In this study, we analyzed the pattern of genomic imprinting in Brassica rapa, a species related to Arabidopsis thaliana. Compared with the ancient karyotype of A. thaliana and B. rapa, B. rapa has a triplicated genome. Many imprinted genes, beyond the estimated number previously reported in other species, were observed. Several imprinted genes have been conserved among species in Brassicaceae. We also observed rapid molecular evolution of imprinted genes compared to non-imprinted genes in B. rapa. Especially, imprinted gene overlapping between species showed more rapid molecular evolution and preferential expression in endosperms. It may imply that a small number of imprinted genes have retained functional roles among diverged species and have been the target of natural selection.

Wild tomato endosperm transcriptomes reveal common roles of genomic imprinting in both nuclear and cellular endosperm

Genomic imprinting is a conspicuous feature of the endosperm, a triploid tissue nurturing the embryo and synchronizing angiosperm seed development. An unknown subset of imprinted genes (IGs) is critical for successful seed development and should have highly conserved functions. Recent genome-wide studies have found limited conservation of IGs among distantly related species, but there is a paucity of data from closely related lineages. Moreover, most studies focused on model plants with nuclear endosperm development, and comparisons with properties of IGs in cellular-type endosperm development are lacking. Using laser-assisted microdissection, we characterized parent-specific expression in the cellular endosperm of three wild tomato lineages (Solanum section Lycopersicon). We identified 1025 candidate IGs and 167 with putative homologs previously identified as imprinted in distantly related taxa with nuclear-type endosperm. Forty-two maternally expressed genes (MEGs) and 17 paternally expressed genes (PEGs) exhibited conserved imprinting status across all three lineages, but differences in power to assess imprinted expression imply that the actual degree of conservation might be higher than that directly estimated (20.7% for PEGs and 10.4% for MEGs). Regardless, the level of shared imprinting status was higher for PEGs than for MEGs, indicating dissimilar evolutionary trajectories. Expression-level data suggest distinct epigenetic modulation of MEGs and PEGs, and gene ontology analyses revealed MEGs and PEGs to be enriched for different functions. Importantly, our data provide evidence that MEGs and PEGs interact in modulating both gene expression and the endosperm cell cycle, and uncovered conserved cellular functions of IGs uniting taxa with cellular- and nuclear-type endosperm.

A Small RNA Pathway Mediates Allelic Dosage in Endosperm

Balance between maternal and paternal genomes within the triploid endosperm is necessary for normal seed development. The majority of endosperm genes are expressed in a 2:1 maternal:paternal ratio, reflecting genomic DNA content. Here, we find that the 2:1 transcriptional ratio is, unexpectedly, actively regulated. In A. thaliana and A. lyrata, endosperm 24-nt small RNAs are reduced in transposable elements and enriched in genes compared with the embryo. We find an inverse relationship between the parent of origin of sRNAs and mRNAs, with genes more likely to be associated with maternally than paternally biased sRNAs. Disruption of the Pol IV sRNA pathway causes a shift toward maternal allele mRNA expression for many genes. Furthermore, paternal inheritance of an RNA Pol IV mutation is sufficient to rescue seed abortion caused by excess paternal genome dosage. Thus, RNA Pol IV mediates the transcriptional balance between maternally and paternally inherited genomes in endosperm.

FACT complex is required for DNA demethylation at heterochromatin during reproduction in Arabidopsis

The DEMETER (DME) DNA glycosylase catalyzes genome-wide DNA demethylation and is required for endosperm genomic imprinting and embryo viability. Targets of DME-mediated DNA demethylation reside in small, euchromatic, AT-rich transposons and at the boundaries of large transposons, but how DME interacts with these diverse chromatin states is unknown. The STRUCTURE SPECIFIC RECOGNITION PROTEIN 1 (SSRP1) subunit of the chromatin remodeler FACT (facilitates chromatin transactions), was previously shown to be involved in the DME-dependent regulation of genomic imprinting in Arabidopsis endosperm. Therefore, to investigate the interaction between DME and chromatin, we focused on the activity of the two FACT subunits, SSRP1 and SUPPRESSOR of TY16 (SPT16), during reproduction in Arabidopsis We found that FACT colocalizes with nuclear DME in vivo, and that DME has two classes of target sites, the first being euchromatic and accessible to DME, but the second, representing over half of DME targets, requiring the action of FACT for DME-mediated DNA demethylation genome-wide. Our results show that the FACT-dependent DME targets are GC-rich heterochromatin domains with high nucleosome occupancy enriched with H3K9me2 and H3K27me1. Further, we demonstrate that heterochromatin-associated linker histone H1 specifically mediates the requirement for FACT at a subset of DME-target loci. Overall, our results demonstrate that FACT is required for DME targeting by facilitating its access to heterochromatin.

Sequential gene activation and gene imprinting during early embryo development in maize

Gene imprinting is a widely observed epigenetic phenomenon in maize endosperm; however, whether it also occurs in the maize embryo remains controversial. Here, we used high-throughput RNA sequencing on laser capture microdissected and manually dissected maize embryos from reciprocal crosses between inbred lines B73 and Mo17 at six time points (3-13 days after pollination, DAP) to analyze allelic gene expression patterns. Co-expression analysis revealed sequential gene activation during maize embryo development. Gene imprinting was observed in maize embryos, and a greater number of imprinted genes were identified at early embryo stages. Sixty-four strongly imprinted genes were identified (at the threshold of 9:1) on manually dissected embryos 5-13 DAP (more imprinted genes at 5 DAP). Forty-one strongly imprinted genes were identified from laser capture microdissected embryos at 3 and 5 DAP (more imprinted genes at 3 DAP). Furthermore, of the 56 genes that were completely imprinted (at the threshold of 99:1), 36 were not previously identified as imprinted genes in endosperm or embryos. In situ hybridization demonstrated that most of the imprinted genes were expressed abundantly in maize embryonic tissue. Our results shed lights on early maize embryo development and provide evidence to support that gene imprinting occurs in maize embryos.

Parental-genome dosage effects on the transcriptome of F1 hybrid triploid embryos of Arabidopsis thaliana

Genomic imprinting in the seed endosperm could be due to unequal parental-genome contribution effects in triploid endosperm tissue that trigger parent-of-origin specific activation and/or silencing of loci prone to genomic imprinting. To determine whether genomic imprinting is triggered by unequal parental-genome contribution effects, we generated a whole-genome transcriptome dataset of F1 hybrid triploid embryos (as mimics of F1 hybrid triploid endosperm). For the vast majority of genes, the parental contributions to their expression levels in the F1 triploid hybrid embryos follow a biallelic and linear expression pattern. While allele-specific expression (ASE) bias was detected, such effects were predominantly parent-of-origin independent. We demonstrate that genomic imprinting is largely absent from F1 triploid embryos, strongly suggesting that neither triploidy nor unequal parental-genome contribution are key triggers of genomic imprinting in plants. However, extensive parental-genome dosage effects on gene expression were observed between the reciprocal F1 hybrid embryos, particularly for genes involved in defence response and nutrient reservoir activity, potentially leading to the seed size differences between reciprocal triploids. We further determined that unequal parental-genome contribution in F1 triploids can lead to overexpression effects that are parent-of-origin dependent, and which are not observed in diploid or tetraploid embryos in which the parental-genome dosage is balanced. Overall, our study demonstrates that neither triploidy nor unequal parental-genome contribution is sufficient to trigger imprinting in plant tissues, suggesting that genomic imprinting is an intrinsic and unique feature of the triploid seed endosperm.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Arabidopsis LORELEI, a Maternally Expressed Imprinted Gene, Promotes Early Seed Development

In flowering plants, the female gametophyte controls pollen tube reception immediately before fertilization and regulates seed development immediately after fertilization, although the controlling mechanisms remain poorly understood. Previously, we showed that LORELEI (LRE), which encodes a putative glycosylphosphatidylinositol-anchored membrane protein, is critical for pollen tube reception by the female gametophyte before fertilization and the initiation of seed development after fertilization. Here, we show that LRE is expressed in the synergid, egg, and central cells of the female gametophyte and in the zygote and proliferating endosperm of the Arabidopsis (Arabidopsis thaliana) seed. Interestingly, LRE expression in the developing seeds was primarily from the matrigenic LRE allele, indicating that LRE expression is imprinted. However, LRE was biallelically expressed in 8-d-old seedlings, indicating that the patrigenic allele does not remain silenced throughout the sporophytic generation. Regulation of imprinted LRE expression is likely novel, as LRE was not expressed in pollen or pollen tubes of mutants defective for MET1, DDM1, RNA-dependent DNA methylation, or MSI-dependent histone methylation. Additionally, the patrigenic LRE allele inherited from these mutants was not expressed in seeds. Surprisingly, and contrary to the predictions of the parental conflict hypothesis, LRE promotes growth in seeds, as loss of the matrigenic but not the patrigenic LRE allele caused delayed initiation of seed development. Our results showed that LRE is a rare imprinted gene that functions immediately after double fertilization and supported the model that a passage through the female gametophyte establishes monoalleleic expression of LRE in seeds and controls early seed development.

Zygotic Genome Activation Occurs Shortly after Fertilization in Maize

The formation of a zygote via the fusion of an egg and sperm cell and its subsequent asymmetric division herald the start of the plant's life cycle. Zygotic genome activation (ZGA) is thought to occur gradually, with the initial steps of zygote and embryo development being primarily maternally controlled, and subsequent steps being governed by the zygotic genome. Here, using maize (Zea mays) as a model plant system, we determined the timing of zygote development and generated RNA-seq transcriptome profiles of gametes, zygotes, and apical and basal daughter cells. ZGA occurs shortly after fertilization and involves ∼10% of the genome being activated in a highly dynamic pattern. In particular, genes encoding transcriptional regulators of various families are activated shortly after fertilization. Further analyses suggested that chromatin assembly is strongly modified after fertilization, that the egg cell is primed to activate the translational machinery, and that hormones likely play a minor role in the initial steps of early embryo development in maize. Our findings provide important insights into gamete and zygote activity in plants, and our RNA-seq transcriptome profiles represent a comprehensive, unique RNA-seq data set that can be used by the research community.