Unlike animals, plants can pause their life cycle as dormant seeds. In both plants and animals, DNA methylation is involved in the regulation of gene expression and genome integrity. In animals, reprogramming erases and re-establishes DNA methylation during development. However, knowledge of reprogramming or reconfiguration in plants has been limited to pollen and the central cell. To better understand epigenetic reconfiguration in the embryo, which forms the plant body, we compared time-series methylomes of dry and germinating seeds to publicly available seed development methylomes.

During the double fertilization process, pollen tubes deliver two sperm cells to an ovule containing the female gametes. In the pollen tube, the vegetative nucleus and sperm cells move together to the apical region where the vegetative nucleus is thought to play a crucial role in controlling the direction and growth of the pollen tube. Here, we report the generation of pollen tubes in Arabidopsis thaliana whose vegetative nucleus and sperm cells are isolated and sealed by callose plugs in the basal region due to apical transport defects induced by mutations in the WPP domain-interacting tail-anchored proteins (WITs) and sperm cell-specific expression of a dominant mutant of the CALLOSE SYNTHASE 3 protein. Through pollen-tube guidance assays, we show that the physiologically anuclear mutant pollen tubes maintain the ability to grow and enter ovules. Our findings provide insight into the sperm cell delivery mechanism and illustrate the independence of the tip-localized vegetative nucleus from directional growth control of the pollen tube.

A new glutelin gene, designated GluD-1, has been discovered by comparing the seed storage proteins from 48 japonica and indica rice cultivars on SDS-PAGE gels. Evidence that GluD-1 is a member of the glutelin family was provided by Western blots using anti-glutelin antiserum and by mapping the gene to the chromosomal glutelin gene cluster. The limited GluD-1 size polymorphism among the rice varieties is due to amino acid substitutions rather than to post-transcriptional modification. GluD-1 is maximally expressed in the starchy endosperm starting at 5 d after flowering (DAF) and increasing through 30 DAF, a major difference from the other glutelins which are primarily expressed in the subaleurone from 10-16 DAF. Only about 0.2 kb of the GluD-1 promoter was sufficient to confer inner starchy endosperm-specific expression. The 0.2 kb truncated GluD-1 promoter contains a bifactorial endosperm box consisting of a truncated GCN4 motif (TGA(G/C)TCA) and AAAG Prolamin box (P box), and ACGT and AACA motifs as cis-regulatory elements. Gel retardation assays and trans-activation experiments indicated that the truncated GCN4 and P box are specifically recognized by RISBZ1 b-ZIP and RPBF Dof activators in vitro, respectively, and are synergistically transactivated, indicating that combinatorial interactions of these motifs are involved in essential endosperm-specific regulation. Furthermore, deviation from the cognate GCN4 motif alters tissue-specific expression in the inner starchy endosperm to include other endosperm tissues.

Mature seed dormancy is a vital plant trait that prevents germination out of season. In Arabidopsis, the trait can be maternally regulated but the underlying mechanisms sustaining this regulation, its general occurrence and its biological significance among accessions are poorly understood. Upon seed imbibition, the endosperm is essential to repress the germination of dormant seeds. Investigation of genomic imprinting in the mature seed endosperm led us to identify a novel set of imprinted genes that are expressed upon seed imbibition. Remarkably, programs of imprinted gene expression are adapted according to the dormancy status of the seed. We provide direct evidence that imprinted genes play a role in regulating germination processes and that preferential maternal allelic expression can implement maternal inheritance of seed dormancy levels.

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.

RNA extraction has been improved by integration of a variety of materials in the protocol, such as phenol, guanidine thiocyanate, and silica, according to the case-specific demands. However, few methods have been designed for high-throughput RNA preparation for large-scale transcriptome studies. In this study, we established a high-throughput guanidinium thiocyanate and isopropyl alcohol based RNA extraction method (HighGI). HighGI is based on simple and phenol-free homemade buffers and the cost is substantially lower than a column-based commercial kit. We demonstrated that the quality and quantity of RNA extracted with HighGI were comparable to those extracted with a conventional phenol/chloroform-based method and a column-based commercial kit. HighGI retained small RNAs less than 200 bp, which are lost with a commercial column-based kit. We also demonstrated that HighGI is readily applicable to semi-automated RNA extraction. HighGI enables high-throughput RNA extraction for large-scale RNA preparation with high yield and quality.

Nutrient transfer from mother to embryo is essential for reproduction in viviparous animals. In the viviparous teleost Xenotoca eiseni (family Goodeidae), the intraovarian embryo intakes the maternal component secreted into the ovarian fluid via the trophotaenia. Our previous study reported that the epithelial layer cells of the trophotaenia incorporate a maternal protein via vesicle trafficking. However, the molecules responsible for the absorption were still elusive. Here, we focused on Cubam (Cubilin-Amnionless) as a receptor involved in the absorption, and cathepsin L as a functional protease in the vesicles. Our results indicated that the Cubam receptor is distributed in the apical surface of the trophotaenia epithelium and then is taken into the intracellular vesicles. The trophotaenia possesses acidic organelles in epithelial layer cells and cathepsin L-dependent proteolysis activity. This evidence does not conflict with our hypothesis that receptor-mediated endocytosis and proteolysis play roles in maternal macromolecule absorption via the trophotaenia in viviparous teleosts. Such nutrient absorption involving endocytosis is not a specific trait in viviparous fish. Similar processes have been reported in the larval stage of oviparous fish or the suckling stage of viviparous mammals. Our findings suggest that the viviparous teleost acquired trophotaenia-based viviparity from a modification of the intestinal absorption system common in vertebrates. This is a fundamental study to understand the strategic variation of the reproductive system in vertebrates.

In the evolutionary history of plants, variation in cis-regulatory elements (CREs) resulting in diversification of gene expression has played a central role in driving the evolution of lineage-specific traits. However, it is difficult to predict expression behaviors from CRE patterns to properly harness them, mainly because the biological processes are complex. In this study, we used cistrome datasets and explainable convolutional neural network (CNN) frameworks to predict genome-wide expression patterns in tomato (Solanum lycopersicum) fruit from the DNA sequences in gene regulatory regions. By fixing the effects of trans-acting factors using single cell-type spatiotemporal transcriptome data for the response variables, we developed a prediction model for crucial expression patterns in the initiation of tomato fruit ripening. Feature visualization of the CNNs identified nucleotide residues critical to the objective expression pattern in each gene, and their effects were validated experimentally in ripening tomato fruit. This cis-decoding framework will not only contribute to the understanding of the regulatory networks derived from CREs and transcription factor interactions, but also provides a flexible means of designing alleles for optimized expression.

In angiosperms, endosperm development comprises a series of developmental transitions controlled by genetic and epigenetic mechanisms that are initiated after double fertilization. Polycomb repressive complex 2 (PRC2) is a key component of these mechanisms that mediate histone H3 lysine 27 trimethylation (H3K27me3); the action of PRC2 is well described in Arabidopsis thaliana but remains uncertain in cereals. In this study, we demonstrate that mutation of the rice (Oryza sativa) gene EMBRYONIC FLOWER2a (OsEMF2a), encoding a zinc-finger containing component of PRC2, causes an autonomous endosperm phenotype involving proliferation of the central cell nuclei with separate cytoplasmic domains, even in the absence of fertilization. Detailed cytological and transcriptomic analyses revealed that the autonomous endosperm can produce storage compounds, starch granules, and protein bodies specific to the endosperm. These events have not been reported in Arabidopsis. After fertilization, we observed an abnormally delayed developmental transition in the endosperm. Transcriptome and H3K27me3 ChIP-seq analyses using endosperm from the emf2a mutant identified downstream targets of PRC2. These included >100 transcription factor genes such as type-I MADS-box genes, which are likely required for endosperm development. Our results demonstrate that OsEMF2a-containing PRC2 controls endosperm developmental programs before and after fertilization.

Root system architecture plays a crucial role in nutrient and water absorption during rice production. Genetic improvement of the rice root system requires elucidating its genetic control. Genome-wide association studies (GWASs) have identified genomic regions responsible for rice root phenotypes. However, candidate gene prioritization around the peak region often suffers from low statistical power and resolution. Transcriptomics enables other statistical mappings, such as transcriptome-wide association study (TWAS) and expression GWAS (eGWAS), which improve candidate gene identification by leveraging the natural variation of the expression profiles. To explore the genes responsible for root phenotypes, we conducted GWAS, TWAS, and eGWAS for 12 root phenotypes in 57 rice accessions using 427,751 single nucleotide polymorphisms (SNPs) and the expression profiles of 16,901 genes expressed in the roots. The GWAS identified three significant peaks, of which the most significant peak responsible for seven root phenotypes (crown root length, crown root surface area, number of crown root tips, lateral root length, lateral root surface area, lateral root volume, and number of lateral root tips) was detected at 6,199,732 bp on chromosome 8. In the most significant GWAS peak region, OsENT1 was prioritized as the most plausible candidate gene because its expression profile was strongly negatively correlated with the seven root phenotypes. In addition to OsENT1, OsEXPA31, OsSPL14, OsDEP1, and OsDEC1 were identified as candidate genes responsible for root phenotypes using TWAS. Furthermore, a cis-eGWAS peak SNP was detected for OsDjA6, which showed the eighth strongest association with lateral root volume in the TWAS. The cis-eGWAS peak SNP for OsDjA6 was in strong linkage disequilibrium (LD) with a GWAS peak SNP on the same chromosome for lateral root volume and in perfect LD with another SNP variant in a putative cis-element at the 518 bp upstream of the gene. These candidate genes provide new insights into the molecular breeding of root system architecture.

The repression of transcription from transposable elements (TEs) by DNA methylation is necessary to maintain genome integrity and prevent harmful mutations. However, under certain circumstances, TEs may escape from the host defense system and reactivate their transcription. In Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), DNA demethylases target the sequences derived from TEs in the central cell, the progenitor cell for the endosperm in the female gametophyte. Genome-wide DNA demethylation is also observed in the endosperm after fertilization. In the present study, we used a custom microarray to survey the transcripts generated from TEs during rice endosperm development and at selected time points in the embryo as a control. The expression patterns of TE transcripts are dynamically up- and downregulated during endosperm development, especially those of miniature inverted-repeat TEs (MITEs). Some TE transcripts were directionally controlled, whereas the other DNA transposons and retrotransposons were not. We also discovered the NUCLEAR FACTOR Y binding motif, CCAAT, in the region near the 5' terminal inverted repeat of Youren, one of the transcribed MITEs in the endosperm. Our results uncover dynamic changes in TE activity during endosperm development in rice.

Salinity critically limits rice metabolism, growth, and productivity worldwide. Improvement of the salt resistance of locally grown high-yielding cultivars is a slow process. The objective of this study was to develop a new salt-tolerant rice germplasm using speed-breeding. Here, we precisely introgressed the hst1 gene, transferring salinity tolerance from "Kaijin" into high-yielding "Yukinko-mai" (WT) rice through single nucleotide polymorphism (SNP) marker-assisted selection. Using a biotron speed-breeding technique, we developed a BC3F3 population, named "YNU31-2-4", in six generations and 17 months. High-resolution genotyping by whole-genome sequencing revealed that the BC3F2 genome had 93.5% similarity to the WT and fixed only 2.7% of donor parent alleles. Functional annotation of BC3F2 variants along with field assessment data indicated that "YNU31-2-4" plants carrying the hst1 gene had similar agronomic traits to the WT under normal growth condition. "YNU31-2-4" seedlings subjected to salt stress (125 mM NaCl) had a significantly higher survival rate and increased shoot and root biomasses than the WT. At the tissue level, quantitative and electron probe microanalyzer studies indicated that "YNU31-2-4" seedlings avoided Na+ accumulation in shoots under salt stress. The "YNU31-2-4" plants showed an improved phenotype with significantly higher net CO2 assimilation and lower yield decline than WT under salt stress at the reproductive stage. "YNU31-2-4" is a potential candidate for a new rice cultivar that is highly tolerant to salt stress at the seedling and reproductive stages, and which might maintain yields under a changing global climate.

Separating male and female sex organs is one of the main strategies used to maintain genetic diversity within a species. However, the genetic determinants and their regulatory mechanisms have been identified in only a few species. In dioecious persimmons, the homeodomain transcription factor, MeGI, which is the target of a Y chromosome-encoded small-RNA, OGI, can determine floral sexuality. The basic features of this system are conserved in the monoecious hexaploid Oriental persimmon, in which an additional epigenetic regulation of MeGI determines floral sexuality. The downstream regulatory pathways of MeGI remain uncharacterized. In this study, we examined transcriptomic data for male and female flowers from monoecious persimmon cultivars to unveil the gene networks orchestrated by MeGI. A network visualization and cistrome assessment suggested that class-1 KNOTTED-like homeobox (KNOX)/ovate family protein (OFP)/growth regulating factors (GRFs) and short vegetative phase (SVP) genes mediate the differences in gynoecium and androecium development between male and female flowers, respectively. The expression of these genes is directly controlled by MeGI. The gene networks also suggested that some cytokinin, auxin, and gibberellin signaling genes function cooperatively in the KNOX/OFP/GRF pathway during gynoecium differentiation. Meanwhile, SVP may repress PI expression in developing androecia. Overall, our results suggest that MeGI evolved the ability to promote gynoecium development and suppress androecium development by regulating KNOX/OFP/GRF and SVP expression levels, respectively. These insights may help to clarify the molecular mechanism underlying the production of unisexual flowers, while also elucidating the physiological background enabling a single-factor system to establish dioecy in plants.

IR64 is a rice variety with high-yield that has been widely cultivated around the world. IR64 has been replaced by modern varieties in most growing areas. Given that modern varieties are mostly progenies or relatives of IR64, genetic analysis of IR64 is valuable for rice functional genomics. However, chromosome-level genome sequences of IR64 have not been available previously. Here, we sequenced the IR64 genome using synthetic long reads obtained by linked-read sequencing and ultra-long reads obtained by nanopore sequencing. We integrated these data and generated the de novo assembly of the IR64 genome of 367 Mb, equivalent to 99% of the estimated size. Continuity of the IR64 genome assembly was improved compared with that of a publicly available IR64 genome assembly generated by short reads only. We annotated 41,458 protein-coding genes, including 657 IR64-specific genes, that are missing in other high-quality rice genome assemblies IRGSP-1.0 of japonica cultivar Nipponbare or R498 of indica cultivar Shuhui498. The IR64 genome assembly will serve as a genome resource for rice functional genomics as well as genomics-driven and/or molecular breeding.

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

To induce transcriptional gene silencing (TGS) of endogenous genes of rice (Oryza sativa L.), we expressed double-strand RNA of each promoter region and thus induced RNA-directed DNA methylation (RdDM). We targeted constitutively expressed genes encoding calnexin (CNX), protein disulphide isomerase (PDIL1-1) and luminal binding protein (BiP1); an endoplasmic reticulum stress-inducible gene (OsbZIP50); and genes with seed-specific expression encoding α-globulin (Glb-1) and glutelin-B4 (GluB4). TGS of four genes was obtained with high efficiency (CNX, 66.7% of regenerated plants; OsBiP1, 67.4%; OsbZIP50, 63.4%; GluB4, 66.1%), whereas the efficiency was lower for PDIL1-1 (33.3%) and Glb-1 TGS lines (10.5%). The heredity of TGS, methylation levels of promoter regions and specificity of silencing of the target gene were investigated in some of the TGS lines. In progeny of CNX and OsbZIP50 TGS lines, suppression of the target genes was preserved (except in the endosperm) even after the removal of trigger genes (T-DNA) by segregation. TGS of CNX was reverted by demethylation treatment, and a significant difference in CG and CHG methylation levels in the -1 to -250 bp region of the CNX promoter was detected between the TGS and revertant lines, suggesting that TGS is closely related to the methylation levels of promoter. TGS exhibited specific suppression towards the target gene compared with post-transcriptional gene silencing when GluB4 gene from glutelin multigene family was targeted. Based on these results, future perspectives and problems to be solved in the application of RdDM to new plant breeding techniques in rice are discussed.

Chloroplasts, which perform photosynthesis, are one of the most important organelles in green plants and algae. Chloroplasts maintain an independent genome that includes important genes encoding their photosynthetic machinery and various housekeeping functions. Owing to its non-recombinant nature, low mutation rates, and uniparental inheritance, the chloroplast genome (plastome) can give insights into plant evolution and ecology and in the development of biotechnological and breeding applications. However, efficient methods to obtain high-quality chloroplast DNA (cpDNA) are currently not available, impeding powerful sequencing and further functional genomics research. To investigate effects on rice chloroplast genome quality, we compared cpDNA extraction by three extraction protocols: liquid nitrogen coupled with sucrose density gradient centrifugation, high-salt buffer, and Percoll gradient centrifugation. The liquid nitrogen-sucrose gradient method gave a high yield of high-quality cpDNA with reliable purity. The cpDNA isolated by this technique was evaluated, resequenced, and assembled de novo to build a robust framework for genomic and genetic studies. Comparison of this high-purity cpDNA with total DNAs revealed the read coverage of the sequenced regions; next-generation sequencing data showed that the high-quality cpDNA eliminated noise derived from contamination by nuclear and mitochondrial DNA, which frequently occurs in total DNA. The assembly process produced highly accurate, long contigs. We summarize the extent to which this improved method of isolating cpDNA from rice can provide practical progress in overcoming challenges related to chloroplast genomes and in further exploring the development of new sequencing technologies.

DNA cytosine methylation is an epigenetic mark associated with silencing of transposable elements (TEs) and heterochromatin formation. In plants, it occurs in three sequence contexts: CG, CHG, and CHH (where H is A, T, or C). The latter does not allow direct inheritance of methylation during DNA replication due to lack of symmetry, and methylation must therefore be re-established every cell generation. Genome-wide association studies (GWAS) have previously shown that CMT2 and NRPE1 are major determinants of genome-wide patterns of TE CHH methylation. Here we instead focus on CHH methylation of individual TEs and TE-families, allowing us to identify the pathways involved in CHH methylation simply from natural variation and confirm the associations by comparing them with mutant phenotypes. Methylation at TEs targeted by the RNA-directed DNA methylation (RdDM) pathway is unaffected by CMT2 variation, but is strongly affected by variation at NRPE1, which is largely responsible for the longitudinal cline in this phenotype. In contrast, CMT2-targeted TEs are affected by both loci, which jointly explain 7.3% of the phenotypic variation (13.2% of total genetic effects). There is no longitudinal pattern for this phenotype, however, because the geographic patterns appear to compensate for each other in a pattern suggestive of stabilizing selection.

A central question in genomic imprinting is how a specific sequence is recognized as the target for epigenetic marking. In both mammals and plants, imprinted genes are often associated with tandem repeats and transposon-related sequences, but the role of these elements in epigenetic gene silencing remains elusive. FWA is an imprinted gene in Arabidopsis thaliana expressed specifically in the female gametophyte and endosperm. Tissue-specific and imprinted expression of FWA depends on DNA methylation in the FWA promoter, which is comprised of two direct repeats containing a sequence related to a SINE retroelement. Methylation of this element causes epigenetic silencing, but it is not known whether the methylation is targeted to the SINE-related sequence itself or the direct repeat structure is also necessary. Here we show that the repeat structure in the FWA promoter is highly diverse in species within the genus Arabidopsis. Four independent tandem repeat formation events were found in three closely related species. Another related species, A. halleri, did not have a tandem repeat in the FWA promoter. Unexpectedly, even in this species, FWA expression was imprinted and the FWA promoter was methylated. In addition, our expression analysis of FWA gene in vegetative tissues revealed high frequency of intra-specific variation in the expression level. In conclusion, we show that the tandem repeat structure is dispensable for the epigenetic silencing of the FWA gene. Rather, SINE-related sequence is sufficient for imprinting, vegetative silencing, and targeting of DNA methylation. Frequent independent tandem repeat formation events in the FWA promoter led us to propose that they may be a consequence, rather than cause, of the epigenetic control. The possible significance of epigenetic variation in reproductive strategies during evolution is also discussed.