Polyploidy is a pervasive evolutionary feature of all flowering plants and some animals, leading to genetic and epigenetic changes that affect gene expression and morphology. DNA methylation changes can produce meiotically stable epialleles, which are transmissible through selection and breeding. However, the relationship between DNA methylation and polyploid plant domestication remains elusive.

Cotton fiber development undergoes rapid and dynamic changes in a single cell type, from fiber initiation, elongation, primary and secondary wall biosynthesis, to fiber maturation. Previous studies showed that cotton genes encoding putative MYB transcription factors and phytohormone responsive factors were induced during early stages of ovule and fiber development. Many of these factors are targets of microRNAs (miRNAs) that mediate target gene regulation by mRNA degradation or translational repression.

Seed development in angiosperms requires a 2:1 maternal-to-paternal genome ratio (2m:1p) in the endosperm. When the ratio is disrupted, the seed development is impaired. Rice interploidy crosses result in endosperm failures, but the underlying molecular mechanisms remain unclear. Here, we report that the defective endosperm in rice interploidy crosses was associated with nonadditive expression of small RNAs and protein-coding genes. Interestingly, 24-nt small interfering RNAs were enriched in the 5' and 3' flanking sequences of nonadditively expressed genes in the interploidy crosses and were negatively associated with the expression of imprinted genes. Furthermore, some PRC2 family genes and DNA methylation-related genes including OsMET1b and OsCMT3a were upregulated in the 2×4 cross (pollinating a diploid "mother" with a tetraploid "father") but repressed in the reciprocal cross. These different epigenetic effects could lead to precocious or delayed cellularization during endosperm development. Notably, many endosperm-preferred genes, including starch metabolic and storage protein genes during grain filling, were found to be associated with DNA methylation or H3K27me3, which are repressed in both 2×4 and 4×2 crosses. WUSCHEL homeobox2 (WOX2)-like (WOX2L), an endosperm-preferred gene, was expressed specifically in the rice endosperm, in contrast to WOX2 expression in the Arabidopsis embryo. Disruption of WOX2L in transgenic rice by CRISPR/Cas9-mediated gene editing blocked starch and protein accumulation, resulting in seed abortion. In addition to gene repression, disrupting epigenetic process in the interploidy crosses also induced expression of stress-responsive genes. Thus, maintaining the 2m:1p genome ratio in the endosperm is essential for normal grain development in rice and other cereal crops.

Cotton has been cultivated and used to make fabrics for at least 7000 years. Two allotetraploid species of great commercial importance, Gossypium hirsutum and Gossypium barbadense, were domesticated after polyploidization and are cultivated worldwide. Although the overall genetic diversity between these two cultivated species has been studied with limited accessions, their population structure and genetic variations remain largely unknown.

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.

During evolution successful allopolyploids must overcome 'genome shock' between hybridizing species but the underlying process remains elusive. Here, we report concerted genomic and epigenomic changes in resynthesized and natural Arabidopsis suecica (TTAA) allotetraploids derived from Arabidopsis thaliana (TT) and Arabidopsis arenosa (AA). A. suecica shows conserved gene synteny and content with more gene family gain and loss in the A and T subgenomes than respective progenitors, although A. arenosa-derived subgenome has more structural variation and transposon distributions than A. thaliana-derived subgenome. These balanced genomic variations are accompanied by pervasive convergent and concerted changes in DNA methylation and gene expression among allotetraploids. The A subgenome is hypomethylated rapidly from F1 to resynthesized allotetraploids and convergently to the T-subgenome level in natural A. suecica, despite many other methylated loci being inherited from F1 to all allotetraploids. These changes in DNA methylation, including small RNAs, in allotetraploids may affect gene expression and phenotypic variation, including flowering, silencing of self-incompatibility and upregulation of meiosis- and mitosis-related genes. In conclusion, concerted genomic and epigenomic changes may improve stability and adaptation during polyploid evolution.

Like many important crops, peanut is a polyploid that underwent polyploidization, evolution, and domestication. The wild allotetraploid peanut species Arachis monticola (A. monticola) is an important and unique link from the wild diploid species to cultivated tetraploid species in the Arachis lineage. However, little is known about A. monticola and its role in the evolution and domestication of this important crop. A fully annotated sequence of ≈2.6 Gb A. monticola genome and comparative genomics of the Arachis species is reported. Genomic reconstruction of 17 wild diploids from AA, BB, EE, KK, and CC groups and 30 tetraploids demonstrates a monophyletic origin of A and B subgenomes in allotetraploid peanuts. The wild and cultivated tetraploids undergo asymmetric subgenome evolution, including homoeologous exchanges, homoeolog expression bias, and structural variation (SV), leading to subgenome functional divergence during peanut domestication. Significantly, SV-associated homoeologs tend to show expression bias and correlation with pod size increase from diploids to wild and cultivated tetraploids. Moreover, genomic analysis of disease resistance genes shows the unique alleles present in the wild peanut can be introduced into breeding programs to improve some resistance traits in the cultivated peanuts. These genomic resources are valuable for studying polyploid genome evolution, domestication, and improvement of peanut production and resistance.

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.

Both ploidy and hybridity affect cell size and growth vigor in plants and animals, but the relative effects of genome dosage and hybridization on biomass, fitness, and gene expression changes have not been systematically examined. Here we performed the first comparative analysis of seed, cell, and flower sizes, starch and chlorophyll content, biomass, and gene expression changes in diploid, triploid, and tetraploid hybrids and their respective parents in three Arabidopsis thaliana ecotypes: Columbia, C24, and Landsberg erecta (Ler). Ploidy affects many morphological and fitness traits, including stomatal size, flower size, and seed weight, whereas hybridization between the ecotypes leads to altered expression of central circadian clock genes and increased starch and chlorophyll content, biomass, and seed weight. However, varying ploidy levels has subtle effects on biomass, circadian clock gene expression, and chlorophyll and starch content. Interestingly, biomass, starch content, and seed weight are significantly different between the reciprocal hybrids at all ploidy levels tested, with the lowest and highest levels found in the reciprocal triploid hybrids, suggesting parent-of-origin effects on biomass, starch content, and seed weight. These findings provide new insights into molecular events of polyploidy and heterosis, as well as complex agronomic traits that are important to biomass and seed production in hybrid and polyploid crops.

Each cotton fiber is a single-celled seed trichome or hair, and over 20,000 fibers may develop semi-synchronously on each seed. The molecular basis for seed hair development is unknown but is likely to share many similarities with leaf trichome development in Arabidopsis. Leaf trichome initiation in Arabidopsis thaliana is activated by GLABROUS1 (GL1) that is negatively regulated by TRIPTYCHON (TRY). Using laser capture microdissection and microarray analysis, we found that many putative MYB transcription factor and structural protein genes were differentially expressed in fiber and non-fiber tissues. Gossypium hirsutum MYB2 (GhMYB2), a putative GL1 homolog, and its downstream gene, GhRDL1, were highly expressed during fiber cell initiation. GhRDL1, a fiber-related gene with unknown function, was predominately localized around cell walls in stems, sepals, seed coats, and pollen grains. GFP:GhRDL1 and GhMYB2:YFP were co-localized in the nuclei of ectopic trichomes in siliques. Overexpressing GhRDL1 or GhMYB2 in A. thaliana Columbia-0 (Col-0) activated fiber-like hair production in 4-6% of seeds and had on obvious effects on trichome development in leaves or siliques. Co-overexpressing GhRDL1 and GhMYB2 in A. thaliana Col-0 plants increased hair formation in ∼8% of seeds. Overexpressing both GhRDL1 and GhMYB2 in A. thaliana Col-0 try mutant plants produced seed hair in ∼10% of seeds as well as dense trichomes inside and outside siliques, suggesting synergistic effects of GhRDL1 and GhMYB2 with try on development of trichomes inside and outside of siliques and seed hair in A. thaliana. These data suggest that a different combination of factors is required for the full development of trichomes (hairs) in leaves, siliques, and seeds. A. thaliana can be developed as a model a system for discovering additional genes that control seed hair development in general and cotton fiber in particular.

Polyploidy is an evolutionary innovation, providing extra sets of genetic material for phenotypic variation and adaptation. It is predicted that changes of gene expression by genetic and epigenetic mechanisms are responsible for novel variation in nascent and established polyploids (Liu and Wendel, 2002; Osborn et al., 2003; Pikaard, 2001). Studying gene expression changes in allopolyploids is more complicated than in autopolyploids, because allopolyploids contain more than two sets of genomes originating from divergent, but related, species. Here we describe two methods that are applicable to the genome-wide analysis of gene expression differences resulting from genome duplication in autopolyploids or interactions between homoeologous genomes in allopolyploids. First, we describe an amplified fragment length polymorphism (AFLP)--complementary DNA (cDNA) display method that allows the discrimination of homoeologous loci based on restriction polymorphisms between the progenitors. Second, we describe microarray analyses that can be used to compare gene expression differences between the allopolyploids and respective progenitors using appropriate experimental design and statistical analysis. We demonstrate the utility of these two complementary methods and discuss the pros and cons of using the methods to analyze gene expression changes in autopolyploids and allopolyploids. Furthermore, we describe these methods in general terms to be of wider applicability for comparative gene expression in a variety of evolutionary, genetic, biological, and physiological contexts.

Polyploidy provides new genetic material that facilitates evolutionary novelty, species adaptation, and crop domestication. Polyploidy often leads to an increase in cell or organism size, which may affect transcript abundance or transcriptome size, but the relationship between polyploidy and transcriptome changes remains poorly understood. Plant cells often undergo endoreduplication, confounding the polyploid effect.

Epigenetic regulation of gene expression is critical for controlling embryonic properties during the embryo-to-seedling phase transition. Here we report that a histone deacetylase19 (HDA19)-associated regulator, scarecrow-like15 (SCL15), is essential for repressing the seed maturation programme in vegetative tissues. SCL15 is expressed in and GFP-tagged SCL15 predominantly localizes to, the vascular bundles particularly in the phloem companion cells and neighbouring specialized cells. Mutation of SCL15 leads to a global shift in gene expression in seedlings to a profile resembling late embryogenesis in seeds. In scl15 seedlings, many genes involved in seed maturation are markedly derepressed with concomitant accumulation of seed 12S globulin; this is correlated with elevated levels of histone acetylation at a subset of seed-specific loci. SCL15 physically interacts with HDA19 and direct targets of HDA19-SCL15 association are identified. These studies reveal that SCL15 acts as an HDA19-associated regulator to repress embryonic traits in seedlings.

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.

Inbreeding depression is widespread across plant and animal kingdoms and may arise from the exposure of deleterious alleles and/or loss of overdominant alleles resulting from increased homozygosity, but these genetic models cannot fully explain the phenomenon. Here, we report epigenetic links to inbreeding depression in maize. Teosinte branched1/cycloidea/proliferating cell factor (TCP) transcription factors control plant development. During successive inbreeding among inbred lines, thousands of genomic regions across TCP-binding sites (TBS) are hypermethylated through the H3K9me2-mediated pathway. These hypermethylated regions are accompanied by decreased chromatin accessibility, increased levels of the repressive histone marks H3K27me2 and H3K27me3, and reduced binding affinity of maize TCP-proteins to TBS. Consequently, hundreds of TCP-target genes involved in mitochondrion, chloroplast, and ribosome functions are down-regulated, leading to reduced growth vigor. Conversely, random mating can reverse corresponding hypermethylation sites and TCP-target gene expression, restoring growth vigor. These results support a unique role of reversible epigenetic modifications in inbreeding depression.

The most widely cultivated cotton (Gossypium hirsutum L., AD-genome) is derived from tetraploidization between A- and D-genome species. G. arboreum L. (A-genome) and G. raimondii Ulbr. (D-genome) are two of closely-related extant progenitors. Gene expression studies in allotetraploid cotton are complicated by the homoeologous loci of A- and D-genome origins. To develop genomic resources for gene expression and cotton breeding, we sequenced and assembled expressed sequence tags (ESTs) derived from G. arboreum and G. raimondii.

Genomic interactions in allopolyploids create expression variation of homoeologous alleles through protein-protein and protein-DNA interactions. However, the molecular basis for this is largely unknown. Here we investigated the protein-protein and protein-DNA interactions among homoeologous transcription factors in the circadian-clock feedback loop, consisting of CCA1 HIKING EXPEDITION (CHE), CIRCADIAN CLOCK ASSOCIATED1 (CCA1), and TIMING OF CAB EXPRESSION1 (TOC1), plus the interaction with a chromatin factor, HISTONE DEACETYLASE1 (HD1). In the allotetraploids formed between A. thaliana (At) and Arabidopsis arenosa (Aa), AtCCA1 is expressed at lower levels than AaCCA1, which could alter clock output traits. The reduced AtCCA1 expressions in the allotetraploids are consistent with the biochemical data that AaCHE showed preferential binding to the AtCCA1 promoter, in which AaCHE interacts with a higher affinity to AtHD1 than AtCHE. AaCHE also showed a higher affinity to TOC1 than AtCHE, consistent with the effect of TOC1 on repressing CCA1. Thus, stronger AaCHE-TOC1 and AaCHE-AtHD1 interactions reduce AtCC1 allelic expression. Our current data suggest a biochemical basis for protein interactions in trans with a preference to the cis-acting elements in heterologous combinations to reduce AtCCA1 expression, while altered CCA1 expression has been shown to affect metabolic and biomass heterosis in interspecific hybrids or allotetraploids.

Interspecific hybrids often increase the levels of heterozygosity and hybrid vigor, but some interspecific hybrid seeds are aborted shortly after fertilization. The mechanism behind this postzygotic seed abortion is poorly understood. Here, we report genome-wide analysis of allelic expression changes in developing siliques and seeds in three F1 interspecific crosses between Arabidopsis thaliana (Col, Ler, or C24) and Arabidopsis arenosa. The majority of maternally expressed genes (MEGs) were shared among all three F1 interspecific crosses, whereas ∼90% of 272 paternally expressed genes (PEGs) were found only in one or two F1 crosses, suggesting a role for disrupted paternal gene expression in seed abortion that varies in different crosses. Consistent with this notion, 12 PEGs in the infertile interspecific hybrids matched MEGs in fertile intraspecific hybrids. This disruption of PEGs in the interspecific hybrids was consistent with the upregulation of the genes in the paternal-excess interploidy cross (2X6) between a diploid mother and a hexaploid father, leading to the seed abortion. Moreover, a subset of PEGs in the interspecific crosses were also upregulated in the intraspecific hybrid met1XWT or meaXWT, in which the mutant of MET1 (DNA METHYLTRANSFERASE1) or MEDEA, a Polycomb Repressive Complex2 gene, was used as the maternal parent. These data suggest that maternal epigenetic factors and paternal gene expression play important roles in the postzygotic seed abortion in interspecific hybrids or neo-allopolyploids.

The morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber.

Circadian rhythms modulate growth and development in all organisms through interlocking transcriptional-translational feedback loops. The transcriptional loop involves chromatin modifications of central circadian oscillators in mammals and plants. However, the molecular basis for rhythmic epigenetic modifications and circadian regulation is poorly understood.