Polyploidization plays a crucial role in plant evolution and is becoming increasingly important in breeding. Structural variations and epigenomic repatterning have been observed in synthetic polyploidizations. However, the mechanisms underlying the occurrence and their effects on gene expression and phenotype remain unknown. Here, we investigated genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics in leaf organs of progeny from the first eight generations of synthetic tetraploids derived from Chinese cabbage (Brassica rapa L. ssp. pekinensis) and cabbage (Brassica oleracea L. var. capitata). One- or two-copy DelDups, with a mean size of 5.70 Mb (400 kb to 65.85 Mb), occurred from the first generation of selfing and thereafter. The duplication of a fragment in one subgenome consistently coincided with the deletion of its syntenic fragment in the other subgenome, and vice versa, indicating that these DelDups were generated by homoeologous exchanges (HEs). Interestingly, the larger the genomic syntenic region, the higher the frequency of DelDups, further suggesting that the pairing of large homoeologous fragments is crucial for HEs. Moreover, we found that the active transcription of continuously distributed genes in local regions is positively associated with the occurrence of HE breakpoints. In addition, the expression of genes within DelDups exhibited a dosage effect, and plants with extra parental genomic fragments generally displayed phenotypes biased toward the corresponding parent. Genome-wide methylation fluctuated remarkably, which did not clearly affect gene expression on a large scale. Our findings provide insights into the early evolution of polyploid genomes, offering valuable knowledge for polyploidization-based breeding. A comprehensive multi-omics landscape across synthetic Brassica napus tetraploids reveals active transcription of homoeologous exchange-related genes.

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

Polyploidization is important in plant evolution and is becoming increasingly important in crop breeding and material creation. Studies have provided evidence for structural variations and epigenomic repatterning in synthetic polyploidizations, but the relationships between structural and epigenomic variations, as well as their effects on gene expression and phenotypic variations are unknown. Here, we investigated the genome-wide copy number variations of fragments (CVFs) and genomic methylation dynamics, in the leaf organ as a representative, of progenies from eight generations that derived from the synthetic tetraploidization between Brassica rapa and Brassica oleracea. We found that half or complete deletion/duplication of fragments ranging in size from 400 kb to 65.85 Mb, with a mean size of 5.70 Mb, occurred frequently from the first generation of selfing and thereafter. The genes located in these CVFs expressed at levels expected for a positive dosage effect, as indicated by the positive association between expression and the copy number of these genes. Plants containing these CVFs also showed distinct phenotypic variations. The whole genome methylation level experienced significant fluctuations in different generations and eventually decreased in the latter generations. Moreover, the CVFs did not show methylation changes from other individuals of the same generation, and the local regions with methylation alterations did not affect gene expression. Together, the dosage effect of CVFs on gene expression was likely the causal factor underlying phenotypic variations in synthetic tetraploids. Our findings provide new insights into the early evolution of polyploid genomes and guide the use of synthetic polyploidizations in breeding.### Competing Interest StatementThe authors have declared no competing interest.

ABSTRACTPolyploidization is important in plant evolution and is becoming increasingly important in crop breeding and material creation. Studies have provided evidence for structural variations and epigenomic repatterning in synthetic polyploidizations, but the relationships between structural and epigenomic variations, as well as their effects on gene expression and phenotypic variations are unknown. Here, we investigated the genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics, in the leaf organ as a representative, of progenies from eight generations that derived from the synthetic tetraploidization betweenBrassica rapaandBrassica oleracea. We found that half or complete deletion/duplication of fragments ranging in size from 400 kb to 65.85 Mb, with a mean size of 5.70 Mb, occurred frequently from the first generation of selfing and thereafter. The genes located in these DelDups expressed at levels expected for a positive dosage effect, as indicated by the positive association between expression and the copy number of these genes. Plants containing these DelDups also showed distinct phenotypic variations. The whole genome methylation level experienced significant fluctuations in different generations and eventually decreased in the latter generations. Moreover, the DelDups did not show methylation changes from other individuals of the same generation, and the local regions with methylation alterations did not affect gene expression. Our findings provide new insights into the early evolution of polyploid genomes and guide the use of synthetic polyploidizations in breeding.

SUMMARYLeaf heading is an important and economically valuable horticultural trait in many vegetables. The formation of a leafy head is a specialized leaf morphogenesis characterized by the emergence of the enlarged incurving leaves. However, the transcriptional regulation mechanisms underlying the transition to leaf heading remain unclear. We carried out large‐scale time‐series transcriptome assays covering the major vegetative growth phases of two headingBrassica crops, Chinese cabbage and cabbage, with the non‐heading morphotype Taicai as the control. A regulatory transition stage that initiated the heading process is identified, accompanied by a developmental switch from rosette leaf to heading leaf in Chinese cabbages. This transition did not exist in the non‐heading control. Moreover, we reveal that the heading transition stage is also conserved in the cabbage clade. Chinese cabbage acquired through domestication a leafy head independently from the origins of heading in other cabbages; phylogenetics supports that the ancestor of all cabbages is non‐heading. The launch of the transition stage is closely associated with the ambient temperature. In addition, examination of the biological activities in the transition stage identified the ethylene pathway as particularly active, and we hypothesize that this pathway was targeted for selection for domestication to form the heading trait specifically in Chinese cabbage. In conclusion, our findings on the transcriptome transition that initiated the leaf heading in Chinese cabbage and cabbage provide a new perspective for future studies of leafy head crops.

Motivation: Over the last decade, RNA-Seq whole-genome sequencing has become a widely used method for measuring and understanding transcriptome-level changes in gene expression. Since RNA-Seq is relatively inexpensive, it can be used on multiple genomes to evaluate gene expression across many different conditions, tissues and cell types. Although many tools exist to map and compare RNA-Seq at the genomics level, few web-based tools are dedicated to making data generated for individual genomic analysis accessible and reusable at a gene-level scale for comparative analysis between genes, across different genomes and meta-analyses. Results: To address this challenge, we revamped the comparative gene expression tool qTeller to take advantage of the growing number of public RNA-Seq datasets. qTeller allows users to evaluate gene expression data in a defined genomic interval and also perform two-gene comparisons across multiple user-chosen tissues. Though previously unpublished, qTeller has been cited extensively in the scientific literature, demonstrating its importance to researchers. Our new version of qTeller now supports multiple genomes for intergenomic comparisons, and includes capabilities for both mRNA and protein abundance datasets. Other new features include support for additional data formats, modernized interface and back-end database and an optimized framework for adoption by other organisms' databases.

BACKGROUND:Despite the prevalence and recurrence of polyploidization in the speciation of flowering plants, its impacts on crop intraspecific genome diversification are largely unknown. Brassica rapa is a mesopolyploid species that is domesticated into many subspecies with distinctive morphotypes.RESULTS:Herein, we report the consequences of the whole-genome triplication (WGT) on intraspecific diversification using a pan-genome analysis of 16 de novo assembled and two reported genomes. Among the genes that derive from WGT, 13.42% of polyploidy-derived genes accumulate more transposable elements and non-synonymous mutations than other genes during individual genome evolution. We denote such genes as being "flexible." We construct the Brassica rapa ancestral genome and observe the continuing influence of the dominant subgenome on intraspecific diversification in B. rapa. The gene flexibility is biased to the more fractionated subgenomes (MFs), in contrast to the more intact gene content of the dominant LF (least fractionated) subgenome. Furthermore, polyploidy-derived flexible syntenic genes are implicated in the response to stimulus and the phytohormone auxin; this may reflect adaptation to the environment. Using an integrated graph-based genome, we investigate the structural variation (SV) landscapes in 524 B. rapa genomes. We observe that SVs track morphotype domestication. Four out of 266 candidate genes for Chinese cabbage domestication are speculated to be involved in the leafy head formation.CONCLUSIONS:This pan-genome uncovers the possible contributions of allopolyploidization on intraspecific diversification and the possible and underexplored role of SVs in favorable trait domestication. Collectively, our work serves as a rich resource for genome-based B. rapa improvement.

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

Small RNAs are abundant in plant reproductive tissues, especially 24-nucleotide (nt) small interfering RNAs (siRNAs). Most 24-nt siRNAs are dependent on RNA Pol IV and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and establish DNA methylation at thousands of genomic loci in a process called RNA-directed DNA methylation (RdDM). In Brassica rapa , RdDM is required in the maternal sporophyte for successful seed development. Here, we demonstrate that a small number of siRNA loci account for over 90% of siRNA expression during B. rapa seed development. These loci exhibit unique characteristics with regard to their copy number and association with genomic features, but they resemble canonical 24-nt siRNA loci in their dependence on RNA Pol IV/RDR2 and role in RdDM. These loci are expressed in ovules before fertilization and in the seed coat, embryo, and endosperm following fertilization. We observed a similar pattern of 24-nt siRNA expression in diverse angiosperms despite rapid sequence evolution at siren loci. In the endosperm, siren siRNAs show a marked maternal bias, and siren expression in maternal sporophytic tissues is required for siren siRNA accumulation. Together, these results demonstrate that seed development occurs under the influence of abundant maternal siRNAs that might be transported to, and function in, filial tissues.

In the version of this article originally published, author Joshua R. Puzey was incorrectly listed as having affiliation 7 (School of Plant Sciences, University of Arizona, Tucson, AZ, USA); affiliation 6 (Department of Biology, College of William and Mary, Williamsburg, VA, USA) is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

Biology has become an information-dense, data-intensive enterprise that requires the ability to access, integrate, and compute on large amounts of data. To enable better use of data, stakeholders (scientists, curators, librarians, and journal publishers) have established standards to ensure that data are findable, accessible, interoperable, and reusable (FAIR; Wilkinson et al., 2016; https://www.force11.org/group/fairgroup/fairprinciples). In a nutshell, FAIR means:

All natural plant species are evolved from ancient polyploids. Polyloidization plays an important role in plant genome evolution, species divergence and crop domestication. We review how the pattern of polyploidy within the plant phylogenetic tree has engendered hypotheses involving mass extinctions, lag-times following polyploidy, and epochs of asexuality. Polyploidization has happened repeatedly in plant evolution and, we conclude, is important for crop domestication. Once duplicated, the effect of purifying selection on any one duplicated gene is relaxed, permitting duplicate gene and regulatory element loss (fractionation). We review the general topic of fractionation, and how some gene categories are retained more than others. Several explanations, including neofunctionalization, subfunctionalization and gene product dosage balance, have been shown to influence gene content over time. For allopolyploids, genetic differences between parental lines immediately manifest as subgenome dominance in the wide-hybrid, and persist and propagate for tens of millions of years. While epigenetic modifications are certainly involved in genome dominance, it has been difficult to determine which came first, the chromatin marks being measured or gene expression. Data support the conclusion that genome dominance and heterosis are antagonistic and mechanically entangled; both happen immediately in the synthetic wide-cross hybrid. Also operating in this hybrid are mechanisms of ‘paralogue interference’. We present a foundation model to explain gene expression and vigour in a wide hybrid/new allotetraploid. This Review concludes that some mechanisms operate immediately at the wide-hybrid, and other mechanisms begin their operations later. Direct interaction of new paralogous genes, as measured using high-resolution chromatin conformation capture, should inform future research and single cell transcriptome sequencing should help achieve specificity while studying gene sub- and neo-functionalization.

Some organisms deploy small RNAs from accessory cells to maintain genome integrity in the zygote, a mechanism that has been proposed but not demonstrated in plants. Here we show that maternal mutations in the Pol IV-dependent small RNA pathway cause abortion of developing seeds in Brassica rapa . Surprisingly, small RNA production is required in maternal somatic tissues, but not in maternal gametes or the developing zygote. We propose that parental influence over zygotic genomes is a common strategy in eukaryotes and that outbreeding species such as B. rapa are key to understanding the role of small RNAs during reproduction.

Van de Peer and colleagues ([Fawcett et al., 2009][1]; [Lohaus and Van de Peer, 2016][2]) have written thought-provoking reviews and commentaries on the nonrandom accumulation of ancient polyploidies at the inferred Cretaceous/Tertiary (K/Pg) boundary in the inferred phylogenetic tree of sequenced

While prokaryotic pan-genomes have been shown to contain many more genes than any individual organism, the prevalence and functional significance of differentially present genes in eukaryotes remains poorly understood. Whole-genome de novo assembly and annotation of 54 lines of the grass Brachypodium distachyon yield a pan-genome containing nearly twice the number of genes found in any individual genome. Genes present in all lines are enriched for essential biological functions, while genes present in only some lines are enriched for conditionally beneficial functions (e.g., defense and development), display faster evolutionary rates, lie closer to transposable elements and are less likely to be syntenic with orthologous genes in other grasses. Our data suggest that differentially present genes contribute substantially to phenotypic variation within a eukaryote species, these genes have a major influence in population genetics, and transposable elements play a key role in pan-genome evolution.

Subgenome dominance is an important phenomenon observed in allopolyploids after whole genome duplication, in which one subgenome retains more genes as well as contributes more to the higher expressing gene copy of paralogous genes. To dissect the mechanism of subgenome dominance, we systematically investigated the relationships of gene expression, transposable element (TE) distribution and small RNA targeting, relating to the multicopy paralogous genes generated from whole genome triplication in Brassica rapa. The subgenome dominance was found to be regulated by a relatively stable factor established previously, then inherited by and shared among B. rapa varieties. In addition, we found a biased distribution of TEs between flanking regions of paralogous genes. Furthermore, the 24-nt small RNAs target TEs and are negatively correlated to the dominant expression of individual paralogous gene pairs. The biased distribution of TEs among subgenomes and the targeting of 24-nt small RNAs together produce the dominant expression phenomenon at a subgenome scale. Based on these findings, we propose a bucket hypothesis to illustrate subgenome dominance and hybrid vigor. Our findings and hypothesis are valuable for the evolutionary study of polyploids, and may shed light on studies of hybrid vigor, which is common to most species.

Oropetium thomaeum is a resurrection plant that can survive extreme water stress through desiccation to complete dryness, providing a model for drought tolerance; here, whole-genome sequencing and assembly of the Oropetium genome using single-molecule real-time sequencing is reported. Oropetium thomaeum is a resurrection plant that can endure extreme water stress through desiccation to complete dryness whilst retaining the ability to revive when water is available, providing a model for drought tolerance. These authors report whole-genome sequencing and assembly of the O. thomaeum genome, using only single-molecule real-time (SMRT) long-read sequencing. Understanding the genomic mechanisms of extreme desiccation tolerance in resurrection plants such as Oropetium may provide targets for engineering drought and stress tolerance in crop plants. Plant genomes, and eukaryotic genomes in general, are typically repetitive, polyploid and heterozygous, which complicates genome assembly1. The short read lengths of early Sanger and current next-generation sequencing platforms hinder assembly through complex repeat regions, and many draft and reference genomes are fragmented, lacking skewed GC and repetitive intergenic sequences, which are gaining importance due to projects like the Encyclopedia of DNA Elements (ENCODE)2. Here we report the whole-genome sequencing and assembly of the desiccation-tolerant grass Oropetium thomaeum. Using only single-molecule real-time sequencing, which generates long (>16 kilobases) reads with random errors, we assembled 99% (244 megabases) of the Oropetium genome into 625 contigs with an N50 length of 2.4 megabases. Oropetium is an example of a ‘near-complete’ draft genome which includes gapless coverage over gene space as well as intergenic sequences such as centromeres, telomeres, transposable elements and rRNA clusters that are typically unassembled in draft genomes. Oropetium has 28,466 protein-coding genes and 43% repeat sequences, yet with 30% more compact euchromatic regions it is the smallest known grass genome. The Oropetium genome demonstrates the utility of single-molecule real-time sequencing for assembling high-quality plant and other eukaryotic genomes, and serves as a valuable resource for the plant comparative genomics community.

The plant gene model remains largely an extrapolation from animals, with the cis functional unit, the gene, cast as a dynamic looping structure. Molecular genetics with model plants continues to make advances; highlighted here are quantitative-occupancy results from the Arabidopsis thaliana (Arabidopsis) Phytochrome-Interacting bHLH transcription Factors (PIF) quartet. Compared to this complex snapshot, results from chromatin occupancy and other Encyclopedia of DNA Elements (ENCODE)-like approaches increase our transcription factor-motif cognate library, but regulation cannot by itself be inferred from binding. Complementary published Arabidopsis conserved noncoding sequence lists are compared, evaluated, merged, and released. Comparative genomic approaches have identified a cis modifier of a gene's expression-hypothetically, a transposon-based 'rheostat'-that works in all cells, times and places.