Pepper (Capsicum spp.) is a widely consumed vegetable with exceptionally large genomes in Solanaceae, yet its genomic evolutionary history remains largely unknown. Here we present 11 high-quality Capsicum genome assemblies, including two gap-free genomes, covering four wild and all five domesticated pepper species. We reconstructed the ancestral karyotype and inferred the evolutionary trajectory of peppers. The expanded and variable genome sizes were attributed to differential transposable element accumulations, which shaped 3D chromatin architecture and introduced mutations associated with traits such as fruit orientation and colour. Using a chromatin accessibility atlas of Capsicum, we highlight the influence of transposable elements on regulatory element evolution. Furthermore, by constructing a haploblock map of 124 pepper core germplasms, we uncover frequent introgressions that facilitate the formation of sweet blocky pepper and the acquisition of important traits such as resistance to pepper mild mottle virus. These findings on the genomic and functional evolution of Capsicum will benefit pepper breeding. This study presents 11 telomere-to-telomere genomes of wild and domesticated pepper, highlights how transposable elements have shaped the evolution of genome structure and regulatory elements and identifies structural variations and introgressions associated with key traits in cultivated pepper.

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.

Cabbage Fusarium wilt (CFW) is a devastating disease caused by the soil-borne fungus Fusarium oxysporum f. sp. conglutinans (Foc). One of the optimal measures for managing CFW is the employment of tolerant/resistant cabbage varieties. However, the interplay between plant genotypes and the pathogen Foc in shaping the rhizosphere microbial community, and the consequent influence of these microbial assemblages on biological resistance, remains inadequately understood. Based on amplicon metabarcoding data, we observed distinct differences in the fungal alpha diversity index (Shannon index) and beta diversity index (unweighted Bray–Curtis dissimilarity) within the rhizosphere of the YR (resistant to Foc) and ZG (susceptible to Foc) cabbage varieties, irrespective of Foc inoculation. Notably, the Shannon diversity shifts in the resistant YR variety were more pronounced following Foc inoculation. Disease-resistant plant variety demonstrate a higher propensity for harboring beneficial microorganisms, such as Pseudomonas, and exhibit superior capabilities in evading harmful microorganisms, in contrast to their disease-susceptible counterparts. Furthermore, the network analysis was performed on rhizosphere-associated microorganisms, including both bacteria and fungi. The networks of association recovered from YR exhibited greater complexity, robustness, and density, regardless of Foc inoculation. Following Foc infection in the YR rhizosphere, there was a notable increase in the dominant bacterium NA13, which is also a hub taxon in the microbial network. Reintroducing NA13 into the soil significantly improved disease resistance in the susceptible ZG variety, by directly inhibiting Foc and triggering defense mechanisms in the roots. The rhizosphere microbial communities of these two cabbage varieties are markedly distinct, with the introduction of the pathogen eliciting significant alterations in their microbial networks which is correlated with susceptibility or resistance to soil-borne pathogens. Furthermore, we identified a rhizobacteria species that significantly boosts disease resistance in susceptible cabbages. Our results indicated that the induction of resistance genes leading to varied responses in microbial communities to pathogens may partly explain the differing susceptibilities of the cabbage varieties tested to CFW.

Constructing inbred lines for self-incompatible species and species with long generation times is challenging, making the use of F1 outcross/segregating populations the main strategy for genetic studies of such species. However, there is a lack of dedicated algorithms/tools for rapid quantitative trait locus (QTL) mapping using the F1 populations. To this end, we have designed and developed an algorithm/tool called OcBSA specifically for QTL mapping of F1 populations. OcBSA transforms the four-haplotype inheritance problem from the two heterozygous diploid parents of the F1 population into the two-haplotype inheritance problem common in current genetic studies by removing the two haplotypes from the heterozygous parent that do not contribute to phenotype segregation in the F1 population. Testing of OcBSA on 1800 simulated F1 populations demonstrated its advantages over other currently available tools in terms of sensitivity and accuracy. In addition, the broad applicability of OcBSA was validated by QTL mapping using seven reported F1 populations of apple, pear, peach, citrus, grape, tea, and rice. We also used OcBSA to map the QTL for flower color in a newly constructed F1 population of potato generated in this study. The OcBSA mapping result was verified by the insertion or deletion markers to be consistent with a previously reported locus harboring the ANTHOCYANIN 2 gene, which regulates potato flower color. Taken together, these results highlight the power and broad utility of OcBSA for QTL mapping using F1 populations and thus a great potential for functional gene mining in outcrossing species. For ease of use, we have developed both Windows and Linux versions of OcBSA, which are freely available at: https://gitee.com/Bioinformaticslab/OcBSA.

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.

The Solanaceae family contains many agriculturally important crops, including tomato, potato, pepper, and tobacco, as well as other species with potential for agricultural development, such as the orphan crops groundcherry, wolfberry, and pepino. Research progress varies greatly among these species, with model crops like tomato being far ahead. This disparity limits the broader agricultural application of other Solanaceae species. In this study, we constructed an interspecies pan-genome for the Solanaceae family and identified various gene retention patterns. Our findings reveal that the activity of specific transposable elements is closely associated with gene fractionation and transposition. The pan-genome was further resolved at the level of T subgenomes, which were generated by Solanaceae-specific paleo-hexaploidization (T event). We demonstrate substantial gene fractionation (loss) and divergence events following ancient duplications. For example, all class A and E flower model genes in Solanaceae originated from two tandemly duplicated genes, which expanded through the γ and T events before fractionating into 10 genes in tomato, each acquiring distinct functions critical for fruit development. Based on these results, we developed the Solanaceae Pan-Genome Database (SolPGD, http://www.bioinformaticslab.cn/SolPGD), which integrates datasets from both inter- and intra-species pan-genomes of Solanaceae. These findings and resources will facilitate future studies of solanaceous species, including orphan crops.

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.

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.

Chinese cabbage (Brassica rapa ssp. pekinensis) is an economically important vegetable crop throughout the world, especially in Asia. High-quality genome sequences are available for Chinese cabbage, but gene functional studies remain challenging. To promote functional genomic studies of Chinese cabbage, we generated an ethyl methane sulfonate (EMS) mutant population of ∼8000 M2 plants using the double haploid inbred line A03 as the parent. The genome of A03 was sequenced and used as a reference for high-throughput functional characterization of gene mutations at the whole-genome level. A total of 300 M2 to M5 EMS mutants were phenotypically screened and then sequenced, revealing 750 629 SNPs and 46 272 InDel mutations that cover 98.27% of all predicted genes in the A03 genome. A forward-genetics approach was successfully used to identify two genes with chloroplast-related functions that are responsible for the yellow leaf mutant trait. A reverse-genetics approach was also used to identify associations between mutations in five genes of the glucosinolate biosynthetic pathway and variations in glucosinolate content of the mutant plants. In addition, we built the Chinese cabbage EMS mutation database (CCEMD, www.bioinformaticslab.cn/EMSmutation/home) to increase the usability of this mutant population resource. In summary, we performed large-scale screening of a heading Chinese cabbage EMS mutant collection at the phenotypic and genotypic levels, which will facilitate gene mining of Chinese cabbage and might also be useful for the study of other Brassica crops.

Whole genome duplication (WGD) and tandem duplication (TD) are important modes of gene amplification and functional innovation, and they are common in plant genome evolution. We analyzed the genomes of three Solanaceae species (Solanum lycopersicum, Capsicum annuum, and Petunia inflata), which share a common distant ancestor with Vitis vinifera, Theobroma cacao, and Coffea canephora but have undergone an extra whole genome triplication (WGT) event. The analysis was used to investigate the phenomenon of tandem gene evolution with (S. lycopersicum) or without WGT (V. vinifera). Among the tandem gene arrays in these genomes, we found that V. vinifera, which has not experienced the WGT event, retained relatively more and larger tandem duplicated gene (TDG) clusters than the Solanaceae species that experienced the WGT event. Larger TDG clusters tend to be derived from older TD events, so this indicates that continuous TDGs (absolute dosage) accumulated during long-term evolution. In addition, WGD and TD show a significant bias in the functional categories of the genes retained. WGD tends to retain dose-sensitive genes related to biological processes, including DNA-binding and transcription factor activity, while TD tends to retain genes involved in stress resistance. WGD and TD also provide more possibilities for gene functional innovation through gene fusion and fission. The TDG cluster containing the tomato fusarium wilt resistance gene I3 contains 15 genes, and one of these genes, Solyc07g055560, has undergone a fusion event after the duplication events. These data provide evidence that helps explain the new functionalization of TDGs in adapting to environmental changes.

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.

DNA methylation plays an important role in plant growth and development, and in regulating the activity of transposable elements (TEs). Research on DNA methylation-related (DMR) genes has been reported in Arabidopsis, but little research on DMR genes has been reported in Brassica rapa and Brassica oleracea, the genomes of which exhibit significant differences in TE content. In this study, we identified 78 and 77 DMR genes in Brassica rapa and Brassica oleracea, respectively. Detailed analysis revealed that the numbers of DMR genes in different DMR pathways varied in B. rapa and B. oleracea. The evolutionary selection pressure of DMR genes in B. rapa and B. oleracea was compared, and the DMR genes showed differential evolution between these two species. The nucleotide diversity (π) and selective sweep (Tajima's D) revealed footprints of selection in the B. rapa and B. oleracea populations. Transcriptome analysis showed that most DMR genes exhibited similar expression characteristics in B. rapa and B. oleracea. This study dissects the evolutionary differences and genetic variations of the DMR genes in B. rapa and B. oleracea, and will provide valuable resources for future research on the divergent evolution of DNA methylation between B. rapa and B. oleracea.