A Pan-plant Protein Complex Map Reveals Deep Conservation and Novel Assemblies.

Plants are foundational for global ecological and economic systems, but most plant proteins remain uncharacterized. Protein interaction networks often suggest protein functions and open new avenues to characterize genes and proteins. We therefore systematically determined protein complexes from 13 plant species of scientific and agricultural importance, greatly expanding the known repertoire of stable protein complexes in plants. By using co-fractionation mass spectrometry, we recovered known complexes, confirmed complexes predicted to occur in plants, and identified previously unknown interactions conserved over 1.1 billion years of green plant evolution. Several novel complexes are involved in vernalization and pathogen defense, traits critical for agriculture. We also observed plant analogs of animal complexes with distinct molecular assemblies, including a megadalton-scale tRNA multi-synthetase complex. The resulting map offers a cross-species view of conserved, stable protein assemblies shared across plant cells and provides a mechanistic, biochemical framework for interpreting plant genetics and mutant phenotypes.

Asymmetric variation in DNA methylation during domestication and de-domestication of rice.

Hundreds of plant species have been domesticated to feed human civilization, while some crops have undergone de-domestication into agricultural weeds, threatening global food security. To understand the genetic and epigenetic basis of crop domestication and de-domestication, we generated DNA methylomes from 95 accessions of wild rice (Oryza rufipogon L.), cultivated rice (Oryza sativa L.) and weedy rice (O. sativa f. spontanea). We detected a significant decrease in DNA methylation over the course of rice domestication but observed an unexpected increase in DNA methylation through de-domestication. Notably, DNA methylation changes occurred in distinct genomic regions for these 2 opposite stages. Variation in DNA methylation altered the expression of nearby and distal genes through affecting chromatin accessibility, histone modifications, transcription factor binding, and the formation of chromatin loops, which may contribute to morphological changes during domestication and de-domestication of rice. These insights into population epigenomics underlying rice domestication and de-domestication provide resources and tools for epigenetic breeding and sustainable agriculture.

Ploidy variation induces butterfly effect on chromatin topology in wheat.

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Imprinting but not cytonuclear interactions determines seed size heterosis in Arabidopsis hybrids.

The parent-of-origin effect on seeds can result from imprinting (unequal expression of paternal and maternal alleles) or combinational effects between cytoplasmic and nuclear genomes, but their relative contributions remain unknown. To discern these confounding factors, we produced cytoplasmic-nuclear substitution (CNS) lines using recurrent backcrossing in Arabidopsis (Arabidopsis thaliana) ecotypes Col-0 and C24. These CNS lines differed only in the nuclear genome (imprinting) or cytoplasm. The CNS reciprocal hybrids with the same cytoplasm displayed ∼20% seed size difference, whereas the seed size was similar between the reciprocal hybrids with fixed imprinting. Transcriptome analyses in the endosperm of CNS hybrids using laser-capture microdissection identified 104 maternally expressed genes (MEGs) and 90 paternally expressed genes (PEGs). These imprinted genes were involved in pectin catabolism and cell wall modification in the endosperm. Homeodomain Glabrous9 (HDG9), an epiallele and one of 11 cross-specific imprinted genes, affected seed size. In the embryo, there were a handful of imprinted genes in the CNS hybrids but only 1 was expressed at higher levels than in the endosperm. AT4G13495 was found to encode a long-noncoding RNA (lncRNA), but no obvious seed phenotype was observed in lncRNA knockout lines. Nuclear RNA Polymerase D1 (NRPD1), encoding the largest subunit of RNA Pol IV, was involved in the biogenesis of small interfering RNAs. Seed size and embryos were larger in the cross using nrpd1 as the maternal parent than in the reciprocal cross, supporting a role of the maternal NRPD1 allele in seed development. Although limited ecotypes were tested, these results suggest that imprinting and the maternal NRPD1-mediated small RNA pathway play roles in seed size heterosis in plant hybrids.

Cytoplasmic genome contributions to domestication and improvement of modern maize.

BACKGROUND: Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated. RESULTS: Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes. CONCLUSIONS: Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean.

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

DNA hypomethylation in tetraploid rice potentiates stress-responsive gene expression for salt tolerance.

Polyploidy is a prominent feature for genome evolution in many animals and all flowering plants. Plant polyploids often show enhanced fitness in diverse and extreme environments, but the molecular basis for this remains elusive. Soil salinity presents challenges for many plants including agricultural crops. Here we report that salt tolerance is enhanced in tetraploid rice through lower sodium uptake and correlates with epigenetic regulation of jasmonic acid (JA)-related genes. Polyploidy induces DNA hypomethylation and potentiates genomic loci coexistent with many stress-responsive genes, which are generally associated with proximal transposable elements (TEs). Under salt stress, the stress-responsive genes including those in the JA pathway are more rapidly induced and expressed at higher levels in tetraploid than in diploid rice, which is concurrent with increased jasmonoyl isoleucine (JA-Ile) content and JA signaling to confer stress tolerance. After stress, elevated expression of stress-responsive genes in tetraploid rice can induce hypermethylation and suppression of the TEs adjacent to stress-responsive genes. These induced responses are reproducible in a recurring round of salt stress and shared between two japonica tetraploid rice lines. The data collectively suggest a feedback relationship between polyploidy-induced hypomethylation in rapid and strong stress response and stress-induced hypermethylation to repress proximal TEs and/or TE-associated stress-responsive genes. This feedback regulation may provide a molecular basis for selection to enhance adaptation of polyploid plants and crops during evolution and domestication.

Temporal Regulation of the Metabolome and Proteome in Photosynthetic and Photorespiratory Pathways Contributes to Maize Heterosis.

Heterosis or hybrid vigor is widespread in plants and animals. Although the molecular basis for heterosis has been extensively studied, metabolic and proteomic contributions to heterosis remain elusive. Here we report an integrative analysis of time-series metabolome and proteome data in maize (Zea mays) hybrids and their inbred parents. Many maize metabolites and proteins are diurnally regulated, and many of these show nonadditive abundance in the hybrids, including key enzymes and metabolites involved in carbon assimilation. Compared with robust trait heterosis, metabolic heterosis is relatively mild. Interestingly, most amino acids display negative mid-parent heterosis (MPH), i.e., having lower values than the average of the parents, while sugars, alcohols, and nucleoside metabolites show positive MPH. From the network perspective, metabolites in the photosynthetic pathway show positive MPH, whereas metabolites in the photorespiratory pathway show negative MPH, which corresponds to nonadditive protein abundance and enzyme activities of key enzymes in the respective pathways in the hybrids. Moreover, diurnally expressed proteins that are upregulated in the hybrids are enriched in photosynthesis-related gene-ontology terms. Hybrids may more effectively remove toxic metabolites generated during photorespiration, and thus maintain higher photosynthetic efficiency. These metabolic and proteomic resources provide unique insight into heterosis and its utilization for high yielding maize and other crop plants.

The Rice Circadian Clock Regulates Tiller Growth and Panicle Development Through Strigolactone Signaling and Sugar Sensing.

Circadian clocks regulate growth and development in plants and animals, but the role of circadian regulation in crop production is poorly understood. Rice (Oryza sativa) grain yield is largely determined by tillering, which is mediated by physiological and genetic factors. Here we report a regulatory loop that involves the circadian clock, sugar, and strigolactone (SL) pathway to regulate rice tiller-bud and panicle development. Rice CIRCADIAN CLOCK ASSOCIATED1 (OsCCA1) positively regulates expression of TEOSINTE BRANCHED1 (OsTB1, also known as FC1), DWARF14 (D14), and IDEAL PLANT ARCHITECTURE1 (IPA1, also known as OsSPL14) to repress tiller-bud outgrowth. Downregulating and overexpressing OsCCA1 increases and reduces tiller numbers, respectively, whereas manipulating PSEUDORESPONSE REGULATOR1 (OsPPR1) expression results in the opposite effects. OsCCA1 also regulates IPA1 expression to mediate panicle and grain development. Genetic analyses using double mutants and overexpression in the mutants show that OsTB1, D14, and IPA1 act downstream of OsCCA1 Sugars repress OsCCA1 expression in roots and tiller buds to promote tiller-bud outgrowth. The circadian clock integrates sugar responses and the SL pathway to regulate tiller and panicle development, providing insights into improving plant architecture and yield in rice and other cereal crops.

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.

LCM and RNA-seq analyses revealed roles of cell cycle and translational regulation and homoeolog expression bias in cotton fiber cell initiation.

BACKGROUND: Cotton fibers provide a powerful model for studying cell differentiation and elongation. Each cotton fiber is a singular and elongated cell derived from epidermal-layer cells of a cotton seed. Efforts to understand this dramatic developmental shift have been impeded by the difficulty of separation between fiber and epidermal cells. RESULTS: Here we employed laser-capture microdissection (LCM) to separate these cell types. RNA-seq analysis revealed transitional differences between fiber and epidermal-layer cells at 0 or 2 days post anthesis. Specifically, down-regulation of putative cell cycle genes was coupled with upregulation of ribosome biosynthesis and translation-related genes, which may suggest their respective roles in fiber cell initiation. Indeed, the amount of fibers in cultured ovules was increased by cell cycle progression inhibitor, Roscovitine, and decreased by ribosome biosynthesis inhibitor, Rbin-1. Moreover, subfunctionalization of homoeologs was pervasive in fiber and epidermal cells, with expression bias towards 10% more D than A homoeologs of cell cycle related genes and 40-50% more D than A homoeologs of ribosomal protein subunit genes. Key cell cycle regulators were predicted to be epialleles in allotetraploid cotton. MYB-transcription factor genes displayed expression divergence between fibers and ovules. Notably, many phytohormone-related genes were upregulated in ovules and down-regulated in fibers, suggesting spatial-temporal effects on fiber cell development. CONCLUSIONS: Fiber cell initiation is accompanied by cell cycle arrest coupled with active ribosome biosynthesis, spatial-temporal regulation of phytohormones and MYB transcription factors, and homoeolog expression bias of cell cycle and ribosome biosynthesis genes. These valuable genomic resources and molecular insights will help develop breeding and biotechnological tools to improve cotton fiber production.

B-BOX DOMAIN PROTEIN28 Negatively Regulates Photomorphogenesis by Repressing the Activity of Transcription Factor HY5 and Undergoes COP1-Mediated Degradation.

Plants have evolved a delicate molecular system to fine-tune their growth and development in response to dynamically changing light environments. In this study, we found that BBX28, a B-box domain protein, negatively regulates photomorphogenic development in a dose-dependent manner in Arabidopsis thaliana BBX28 interferes with the binding of transcription factor HY5 to the promoters of its target genes through physical interactions, thereby repressing its activity and negatively affecting HY5-regulated gene expression. In darkness, BBX28 associates with CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1) and undergoes COP1-mediated degradation via the 26S proteasome system. Collectively, these results demonstrate that BBX28 acts as a key factor in the COP1-HY5 regulatory hub by maintaining proper HY5 activity to ensure normal photomorphogenic development in plants.

COP1 SUPPRESSOR 4 promotes seedling photomorphogenesis by repressing CCA1 and PIF4 expression in Arabidopsis.

CONSTITUTIVELY PHOTOMORPHOGENIC 1 (COP1) and DE-ETIOLATED 1 (DET1) are founding components of two central repressor complexes of photomorphogenesis that trigger the degradation of a larger number of photomorphogenic-promoting factors in darkness. Here, we identify COP1 SUPPRESSOR 4 (CSU4) as a genetic suppressor of the cop1-6 mutation. Mutations in CSU4 largely rescued the constitutively photomorphogenic phenotype of cop1-6 and det1-1 in darkness. Loss of CSU4 function resulted in significantly longer hypocotyl in the light. Further biochemical studies revealed that CSU4 physically interacts with CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1) and negatively regulates its transcriptional repression activity toward its targets. CSU4 represses the expression of CCA1 in the early morning and of PHYTOCHROME INTERACTING FACTOR 4 (PIF4) in the early evening. Our study suggests that CSU4 acts as a negative regulator of CCA1 via physically associating with CCA1, which in turn, likely serves to repress expression of CCA1 and PIF4 to promote photomorphogenesis.

Diurnal down-regulation of ethylene biosynthesis mediates biomass heterosis.

Heterosis is widely applied in agriculture; however, the underlying molecular mechanisms for superior performance are not well understood. Ethylene biosynthesis and signaling genes are shown to be down-regulated in Arabidopsis interspecific hybrids. Ethylene is a plant hormone that promotes fruit ripening and maturation but inhibits hypocotyl elongation. Here we report that application of exogenous ethylene could eliminate biomass vigor in Arabidopsis thaliana F1 hybrids, suggesting a negative role of ethylene in heterosis. Ethylene biosynthesis is mediated by the rate-limiting enzyme, 1-aminocyclopropane-1-carboxylate synthase (ACS). Down-regulation of ACS genes led to the decrease of ethylene production, which was associated with the high-vigor F1 hybrids, but not with the low-vigor ones. At the mechanistic level, expression of ACS genes was down-regulated diurnally and indirectly by Circadian Clock Associated 1 (CCA1) during the day and directly by Phyotochrome-Interacting Factor 5 (PIF5) at night. Consistent with the negative role of ethylene in plant growth, biomass vigor was higher in the acs mutants than in wild-type plants, while increasing endogenous ethylene production in the hybridizing parents reduced growth vigor in the hybrids. Thus, integrating circadian rhythms and light signaling into ethylene production is another regulatory module of complex biological networks, leading to biomass heterosis in plants.

Dynamic and reversible DNA methylation changes induced by genome separation and merger of polyploid wheat.

BACKGROUND: Wheat is a powerful genetic model for studying polyploid evolution and crop domestication. Hexaploid bread wheat was formed by two rounds of interspecific hybridization and polyploidization, processes which are often accompanied by genetic and epigenetic changes, including DNA methylation. However, the extent and effect of such changes during wheat evolution, particularly from tetraploid-to-hexaploid wheat, are currently elusive. RESULTS: Here we report genome-wide DNA methylation landscapes in extracted tetraploid wheat (ETW, AABB), natural hexaploid wheat (NHW, AABBDD), resynthesized hexaploid wheat (RHW, AABBDD), natural tetraploid wheat (NTW, AABB), and diploid (DD). In the endosperm, levels of DNA methylation, especially in CHG (H=A, T, or C) context, were dramatically decreased in the ETW relative to natural hexaploid wheat; hypo-differentially methylated regions (DMRs) (850,832) were 24-fold more than hyper-DMRs (35,111). Interestingly, those demethylated regions in ETW were remethylated in the resynthesized hexaploid wheat after the addition of the D genome. In ETW, hypo-DMRs correlated with gene expression, and TEs were demethylated and activated, which could be silenced in the hexaploid wheat. In NHW, groups of TEs were dispersed in genic regions of three subgenomes, which may regulate the expression of TE-associated genes. Further, hypo-DMRs in ETW were associated with reduced H3K9me2 levels and increased expression of histone variant genes, suggesting concerted epigenetic changes after separation from the hexaploid. CONCLUSION: Genome merger and separation provoke dynamic and reversible changes in chromatin and DNA methylation. These changes correlate with altered gene expression and TE activity, which may provide insights into polyploid genome and wheat evolution.

Asymmetrical changes of gene expression, small RNAs and chromatin in two resynthesized wheat allotetraploids.

Polyploidy occurs in some animals and all flowering plants, including important crops such as wheat. The consequences of polyploidy in crops remain elusive, partly because their progenitors are unknown. Using two resynthesized wheat allotetraploids Sl Sl AA and AADD with known diploid progenitors, we analyzed mRNA and small RNA transcriptomes in the endosperm, compared transcriptomes between endosperm and root in AADD, and examined chromatin changes in the allotetraploids. In the endosperm, there were more non-additively expressed genes in Sl Sl AA than in AADD. In AADD, non-additively expressed genes were developmentally regulated, and the majority (62-70%) were repressed. The repressed genes in AADD included a group of histone methyltransferase gene homologs, which correlated with reduced histone H3K9me2 levels and activation of various transposable elements in AADD. In Sl Sl AA, there was a tendency for expression dominance of Sl over A homoeologs, but the histone methyltransferase gene homologs were additively expressed, correlating with insignificant changes in histone H3K9me2 levels. Moreover, more 24-nucleotide small inferring RNAs (siRNAs) in the A subgenome were disrupted in AADD than in Sl Sl AA, which were associated with expression changes of siRNA-associated genes. Our results indicate that asymmetrical changes in siRNAs, chromatin modifications, transposons and gene expression coincide with unstable AADD genomes and stable Sl Sl AA genomes, which could help explain the evolutionary trajectories of wheat allotetraploids formed by different progenitors.

Genomic and epigenetic insights into the molecular bases of heterosis.

Heterosis, also known as hybrid vigour, is widespread in plants and animals, but the molecular bases for this phenomenon remain elusive. Recent studies in hybrids and allopolyploids using transcriptomic, proteomic, metabolomic, epigenomic and systems biology approaches have provided new insights. Emerging genomic and epigenetic perspectives suggest that heterosis arises from allelic interactions between parental genomes, leading to altered programming of genes that promote the growth, stress tolerance and fitness of hybrids. For example, epigenetic modifications of key regulatory genes in hybrids and allopolyploids can alter complex regulatory networks of physiology and metabolism, thus modulating biomass and leading to heterosis. The conceptual advances could help to improve plant and animal productivity through the manipulation of heterosis.

Temporal Shift of Circadian-Mediated Gene Expression and Carbon Fixation Contributes to Biomass Heterosis in Maize Hybrids.

Heterosis has been widely used in agriculture, but the molecular mechanism for this remains largely elusive. In Arabidopsis hybrids and allopolyploids, increased photosynthetic and metabolic activities are linked to altered expression of circadian clock regulators, including CIRCADIAN CLOCK ASSOCIATED1 (CCA1). It is unknown whether a similar mechanism mediates heterosis in maize hybrids. Here we report that higher levels of carbon fixation and starch accumulation in the maize hybrids are associated with altered temporal gene expression. Two maize CCA1 homologs, ZmCCA1a and ZmCCA1b, are diurnally up-regulated in the hybrids. Expressing ZmCCA1 complements the cca1 mutant phenotype in Arabidopsis, and overexpressing ZmCCA1b disrupts circadian rhythms and biomass heterosis. Furthermore, overexpressing ZmCCA1b in maize reduced chlorophyll content and plant height. Reduced height stems from reduced node elongation but not total node number in both greenhouse and field conditions. Phenotypes are less severe in the field than in the greenhouse, suggesting that enhanced light and/or metabolic activities in the field can compensate for altered circadian regulation in growth vigor. Chromatin immunoprecipitation-sequencing (ChIP-seq) analysis reveals a temporal shift of ZmCCA1-binding targets to the early morning in the hybrids, suggesting that activation of morning-phased genes in the hybrids promotes photosynthesis and growth vigor. This temporal shift of ZmCCA1-binding targets correlated with nonadditive and additive gene expression in early and late stages of seedling development. These results could guide breeding better hybrid crops to meet the growing demand in food and bioenergy.

Dynamic Roles for Small RNAs and DNA Methylation during Ovule and Fiber Development in Allotetraploid Cotton.

DNA methylation is essential for plant and animal development. In plants, methylation occurs at CG, CHG, and CHH (H = A, C or T) sites via distinct pathways. Cotton is an allotetraploid consisting of two progenitor genomes. Each cotton fiber is a rapidly-elongating cell derived from the ovule epidermis, but the molecular basis for this developmental transition is unknown. Here we analyzed methylome, transcriptome, and small RNAome and revealed distinct changes in CHH methylation during ovule and fiber development. In ovules, CHH hypermethylation in promoters correlated positively with siRNAs, inducing RNA-dependent DNA methylation (RdDM), and up-regulation of ovule-preferred genes. In fibers, the ovule-derived cells generated additional heterochromatic CHH hypermethylation independent of RdDM, which repressed transposable elements (TEs) and nearby genes including fiber-related genes. Furthermore, CHG and CHH methylation in genic regions contributed to homoeolog expression bias in ovules and fibers. Inhibiting DNA methylation using 5-aza-2'-deoxycytidine in cultured ovules has reduced fiber cell number and length, suggesting a potential role for DNA methylation in fiber development. Thus, RdDM-dependent methylation in promoters and RdDM-independent methylation in TEs and nearby genes could act as a double-lock feedback mechanism to mediate gene and TE expression, potentiating the transition from epidermal to fiber cells during ovule and seed development.

Both maternally and paternally imprinted genes regulate seed development in rice.

Genetic imprinting refers to the unequal expression of paternal and maternal alleles of a gene in sexually reproducing organisms, including mammals and flowering plants. Although many imprinted genes have been identified in plants, the functions of these imprinted genes have remained largely uninvestigated. We report genome-wide analysis of gene expression, DNA methylation and small RNAs in the rice endosperm and functional tests of five imprinted genes during seed development using Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated gene9 (CRISPR/Cas9) gene editing technology. In the rice endosperm, we identified 162 maternally expressed genes (MEGs) and 95 paternally expressed genes (PEGs), which were associated with miniature inverted-repeat transposable elements, imprinted differentially methylated loci and some 21-22 small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs). Remarkably, one-third of MEGs and nearly one-half of PEGs were associated with grain yield quantitative trait loci. Most MEGs and some PEGs were expressed specifically in the endosperm. Disruption of two MEGs increased the amount of small starch granules and reduced grain and embryo size, whereas mutation of three PEGs reduced starch content and seed fertility. Our data indicate that both MEGs and PEGs in rice regulate nutrient metabolism and endosperm development, which optimize seed development and offspring fitness to facilitate parental-offspring coadaptation. These imprinted genes and mechanisms could be used to improve the grain yield of rice and other cereal crops.