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.

SUMMARYDeep learning offers new approaches to investigate the mechanisms underlying complex biological phenomena, such as subgenome dominance. Subgenome dominance refers to the dominant expression and/or biased fractionation of genes in one subgenome of allopolyploids, which has shaped the evolution of a large group of plants. However, the underlying cause of subgenome dominance remains elusive. Here, we adopt deep learning to construct two convolutional neural network (CNN) models, binary expression model (BEM) and homoeolog contrast model (HCM), to investigate the mechanism underlying subgenome dominance using DNA sequence and methylation sites. We apply these CNN models to analyze three representative polyploidization systems, Brassica, Gossypium, and Cucurbitaceae, each with available ancient and neo/synthetic polyploidized genomes. The BEM shows that DNA sequence of the promoter region can accurately predict whether a gene is expressed or not. More importantly, the HCM shows that the DNA sequence of the promoter region predicts dominant expression status between homoeologous gene pairs retained from ancient polyploidizations, thus predicting subgenome dominance associated with these events. However, HCM fails to predict gene expression dominance between new homoeologous gene pairs arising from the neo/synthetic polyploidizations. These results are consistent across the three plant polyploidization systems, indicating broad applicability of our models. Furthermore, the two models based on methylation sites produce similar results. These results show that subgenome dominance is associated with long‐term sequence differentiation between the promoters of homoeologs, suggesting that subgenome expression dominance precedes and is the driving force or even the determining factor for sequence divergence between subgenomes following polyploidization.

Brassica oleracea , globally cultivated for its vegetable crops, consists of very diverse morphotypes, characterized by specialized enlarged organs as harvested products. This makes B. oleracea an ideal model for studying rapid evolution and domestication. We constructed a B. oleracea pan-genome from 27 high-quality genomes representing all morphotypes and their wild relatives. We identified structural variations (SVs) among these genomes and characterized these in 704 B. oleracea accessions using graph-based genome tools. We show that SVs exert bidirectional effects on the expression of numerous genes, either suppressing through DNA methylation or promoting probably by harboring transcription factor-binding elements. The following examples illustrate the role of SVs modulating gene expression: SVs promoting BoPNY and suppressing BoCKX3 in cauliflower/broccoli, suppressing BoKAN1 and BoACS4 in cabbage and promoting BoMYBtf in ornamental kale. These results provide solid evidence for the role of SVs as dosage regulators of gene expression, driving B. oleracea domestication and diversification.

Brassica oleracea has been developed into many important crops, including cabbage, kale, cauliflower, broccoli and so on. The genome and gene annotation of cabbage (cultivar JZS), a representative morphotype of B. oleracea, has been widely used as a common reference in biological research. Although its genome assembly has been updated twice, the current gene annotation still lacks information on untranslated regions (UTRs) and alternative splicing (AS). Here, we constructed a high-quality gene annotation (JZSv3) using a full-length transcriptome acquired by nanopore sequencing, yielding a total of 59 452 genes and 75 684 transcripts. Additionally, we re-analyzed the previously reported transcriptome data related to the development of different tissues and cold response using JZSv3 as a reference, and found that 3843 out of 11 908 differentially expressed genes (DEGs) underwent AS during the development of different tissues and 309 out of 903 cold-related genes underwent AS in response to cold stress. Meanwhile, we also identified many AS genes, including BolLHCB5 and BolHSP70, that displayed distinct expression patterns within variant transcripts of the same gene, highlighting the importance of JZSv3 as a pivotal reference for AS analysis. Overall, JZSv3 provides a valuable resource for exploring gene function, especially for obtaining a deeper understanding of AS regulation mechanisms.

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.

Dear Editor, Potato is a vital food security crop and is ranked as the world's third most important food crop after rice and wheat.In 2011,the first genome assembly of a doubled monoploid potato DM1-3 516 R44(DM)was released(Potato Genome Sequencing Consortium,2011),which has been widely used as one of the most popular reference genomes in the last decade and served as a valuable resource in plant genomics and potato genetics community(Leisner et al.,2018;Yang et al.,2020;Zheng et al.,2020).

SUMMARY Actinidia (‘ Mihoutao ’ in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2 n = 174) with diploid A. chinensis (2 n = 58) as the putative parent. A CCS‐based assembly anchored to a high‐resolution linkage map provided a chromosome‐resolved genome for hexaploid A. deliciosa yielded a 3.91‐Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis ; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele‐defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

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.

In order to trace the evolution of chromosomes of Solanaceae species，a systematic comparative genomic analysis was carried out on seven Solanaceae species. Using the genome of model plant tomato as a reference template，we comprehensively identified the collinear relationship among seven species of Solanaceae，and established the genome fragment unit system of Solanaceae. Based on this system，we further reconstructed the ancestral genome of Solanaceae with 12 chromosomes. The results showed that the ancestral genome gradually formed the current genome structure of Solanaceae mainly through translocation and inversion of genome fragments，and these rearrangement events mostly occurred in the regions of centromere and other repeated sequences were enriched. This study provided an important theoretical reference for further understanding the genome evolution and species differentiation of Solanaceae and other plants.

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

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.

Pepper ( Capsicum spp.) is one of the earliest domesticated crops, providing a unique pungent sensation when eaten. Through the construction of the first pepper variome, we describe the main groups that emerged during domestication and breeding of C. annuum , their relations and temporal succession, and the molecular events underlying the main transitions. The initial differentiation in fruit shape and pungency, increase in fruit weight, and transition from erect to pendent fruits, and the recent appearance of blocky, large, sweet fruits (bell peppers), were accompanied by strong selection/fixation of key alleles and introgressions in two large genomic regions. Furthermore, we describe the identification of Up , a key domestication gene controlling erect vs pendent fruit orientation, encoding a BIG GRAIN protein involved in auxin transport, and Flip1 associated with capsaicinoid content, encoding a protein involved in phospholipid flipping. The function of Up was confirmed by virus-induced gene silencing. These findings constitute a cornerstone for understanding the domestication and differentiation of a key horticultural crop.

Rubus corchorifolius (Shanmei or mountain berry, 2n = 14) is widely distributed in China, and its fruits possess high nutritional and medicinal values. Here, we reported a high-quality chromosome-scale genome assembly of Shanmei, with contig size of 215.69 Mb and 26,696 genes. Genome comparison among Rosaceae species showed that Shanmei and Fupenzi (Rubus chingii Hu) were most closely related, followed by blackberry (Rubus occidentalis), and that environmental adaptation-related genes were expanded in the Shanmei genome. Further resequencing of 101 samples of Shanmei collected from four regions in the provinces of Yunnan, Hunan, Jiangxi, and Sichuan in China revealed that among these samples, the Hunan population of Shanmei possessed the highest diversity and represented the more ancestral population. Moreover, the Yunnan population underwent strong selection based on the nucleotide diversity, linkage disequilibrium, and historical effective population size analyses. Furthermore, genes from candidate genomic regions that showed strong divergence were significantly enriched in the flavonoid biosynthesis and plant hormone signal transduction pathways, indicating the genetic basis of adaptation of Shanmei to the local environment. The high-quality assembled genome and the variome dataset of Shanmei provide valuable resources for breeding applications and for elucidating the genome evolution and ecological adaptation of Rubus species.

Miscanthus , a member of the Saccharinae subtribe that includes sorghum and sugarcane, has been widely studied as a feedstock for cellulosic biofuel production. Here, we report the sequencing and assembly of the Miscanthus floridulus genome by the integration of PacBio sequencing and Hi-C mapping, resulting in a chromosome-scale, high-quality reference genome of the genus Miscanthus . Comparisons among Saccharinae genomes suggest that Sorghum split first from the common ancestor of Saccharum and Miscanthus , which subsequently diverged from each other, with two successive whole-genome duplication events occurring independently in the Saccharum genus and one whole-genome duplication occurring in the Miscanthus genus. Fusion of two chromosomes occurred during rediploidization in M. floridulus and no significant subgenome dominance was observed. A survey of cellulose synthases (CesA) in M. floridulus revealed quite high expression of most CesA genes in growing stems, which is in agreement with the high cellulose content of this species. Resequencing and comparisons of 75 Miscanthus accessions suggest that M. lutarioriparius is genetically close to M. sacchariflorus and that M. floridulus is more distantly related to other species and is more genetically diverse. This study provides a valuable genomic resource for molecular breeding and improvement of Miscanthus and Saccharinae crops.

Miscanthus is recognized as a promising lignocellulosic crop for the production of bioethanol and bioproducts worldwide. To facilitate the identification of agronomical important traits and establish genetics knowledge, two genetic maps were developed from a controlled interspecific cross between M. floridulus and M. sacchariflorus. A total of 650 SSR markers were mapped in M. floridulus, spanning 19 linkage groups and 2053.31 cM with an average interval of 3.25 cM. The map of M. sacchariflorus comprised 495 SSR markers in 19 linkage groups covering 1684.86 cM with an average interval of 3.54 cM. The estimation on genome length indicated that the genome coverage of parental genetic maps were 93.87% and 89.91%, respectively. Eighty-eight bi-parental common markers were allowed to connect the two maps, and six pairs of syntenic linkage groups were recognized. Furthermore, quantitative trait loci (QTL) mapping of three agronomic traits, namely, plant height (PH), heading time (HT), and flowering time (FT), demonstrated that a total of 66 QTLs were identified in four consecutive years using interval mapping and multiple-QTL model. The LOD value of these QTLs ranged from 2.51 to 10.60, and the phenotypic variation explained varied from 9.50 to 37.10%. QTL cluster in syntenic groups MF19/MS7 contained six stable QTLs associated with PH, HT, and FT. In conclusion, we report for the first time the genetic mapping of biomass traits in M. floridulus and M. sacchariflorus. These results will be a valuable genetic resource, facilitating the discovery of essential genes and breeding of Miscanthus.