Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

Temperature-dependent jumonji demethylase modulates flowering time by targeting H3K36me2/3 in Brassica rapa

Global warming has a severe impact on the flowering time and yield of crops. Histone modifications have been well-documented for their roles in enabling plant plasticity in ambient temperature. However, the factor modulating histone modifications and their involvement in habitat adaptation have remained elusive. In this study, through genome-wide pattern analysis and quantitative-trait-locus (QTL) mapping, we reveal that BrJMJ18 is a candidate gene for a QTL regulating thermotolerance in thermotolerant B. rapa subsp. chinensis var. parachinensis (or Caixin, abbreviated to Par). BrJMJ18 encodes an H3K36me2/3 Jumonji demethylase that remodels H3K36 methylation across the genome. We demonstrate that the BrJMJ18 allele from Par (BrJMJ18Par) influences flowering time and plant growth in a temperature-dependent manner via characterizing overexpression and CRISPR/Cas9 mutant plants. We further show that overexpression of BrJMJ18Par can modulate the expression of BrFLC3, one of the five BrFLC orthologs. Furthermore, ChIP-seq and transcriptome data reveal that BrJMJ18Par can regulate chlorophyll biosynthesis under high temperatures. We also demonstrate that three amino acid mutations may account for function differences in BrJMJ18 between subspecies. Based on these findings, we propose a working model in which an H3K36me2/3 demethylase, while not affecting agronomic traits under normal conditions, can enhance resilience under heat stress in Brassica rapa.

Genetic characterization and fine mapping of a recessive genic male-sterile gene in flowering Chinese cabbage (Brassica rapa var. parachinensis)

Brassica vegetables exhibit pronounced heterosis; nevertheless, investigations on fertility-related genes are scarce. The present study scrutinized a recessive genic male-sterile line, 7-3A, capable of generating a completely sterile population, holding significant promise for flowering Chinese cabbage breeding. By whole-genome resequencing of sterile and fertile plants, the male-sterile gene was confined to approximately 185 kb on chromosome A07, situated between markers C719 and NP10 in Brassica rapa var. Chiifu-401. Notably, substantial structural variation was identified within this region across diverse Brassica rapa reference genomes. Despite discernible expression level disparities of a homologous gene, Bnams4b, between male sterile and fertile plants, no sequence divergence was detected. Further elucidation is required to pinpoint a novel sterile gene within the candidate interval. This investigation contributes to the advancement of both the molecular-assisted breeding scheme for flowering Chinese cabbage and the comprehension of male sterility mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04005-7.

Expanding the targeting scope of CRISPR/Cas9-mediated genome editing by Cas9 variants in Brassica

CRISPR/Cas9, presently the most widely used genome editing technology, has provided great potential for functional studies and plant breeding. However, the strict requirement for a protospacer adjacent motif (PAM) has hindered the application of the CRISPR/Cas9 system because the number of targetable genomic sites is limited. Recently, the engineered variants Cas9-NG, SpG, and SpRY, which recognize non-canonical PAMs, have been successfully tested in plants (mainly in rice, a monocot). In this study, we evaluated the targeted mutagenesis capabilities of these Cas9 variants in two important Brassica vegetables, Chinese cabbage (Brassica rapa spp. pekinensis) and cabbage (Brassica oleracea var. capitata). Both Cas9-NG and SpG induced efficient mutagenesis at NGN PAMs, while SpG outperformed Cas9-NG at NGC and NGT PAMs. SpRY achieved efficient editing at almost all PAMs (NRN > NYN), albeit with some self-targeting activity at transfer (T)-DNA sequences. And SpRY-induced mutants were detected in cabbage plants in a PAM-less fashion. Moreover, an adenine base editor was developed using SpRY and TadA8e deaminase that induced A-to-G conversions within target sites using non-canonical PAMs. Together, the toolboxes developed here induced successful genome editing in Chinese cabbage and cabbage. Our work further expands the targeting scope of genome editing and paves the way for future basic research and genetic improvement in Brassica.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-024-00155-7.

Development of self-compatible Chinese cabbage lines of Chiifu through marker-assisted selection

The continuously refined genome assembly of the Chinese cabbage accession Chiifu is widely recognized as the reference for Brassica rapa. However, the high self-incompatibility of Chiifu limits its broader utilization. In this study, we report the development of self-compatible Chiifu lines through a meticulous marker-assisted selection (MAS) strategy, involving the substitution of the Chiifu allele of MLPK (M-locus protein kinase) with that from the self-compatible Yellow Sarson (YS). A YS-based marker (SC-MLPK) was employed to screen 841 B. rapa accessions, confirming that all eight accessions with the mlpk/mlpk (mm) genotype exhibited self-compatibility. Additionally, we designed 131 High-Resolution Melting (HRM) markers evenly distributed across the B. rapa genome as genomic background selection (GBS) markers to facilitate the introgression of self-compatibility from YS into Chiifu along with SC-MLPK. Genome background screening revealed that the BC3S1 population had a proportion of the recurrent parent genome (PR) ranging from 93.9% to 98.5%. From this population, we identified self-compatible individuals exhibiting a high number of pollen tubes penetrating stigmas (NPT) (>25) and a maximum compatibility index (CI) value of 7.5. Furthermore, we selected two individuals demonstrating significant similarity to Chiifu in both genetic background and morphological appearance, alongside self-compatibility. These selected individuals were self-pollinated to generate two novel lines designated as SC-Chiifu Lines. The development of these self-compatible Chiifu lines, together with the SC-MLPK marker and the set of HRM markers, represents valuable tools for B. rapa genetics and breeding.

Identification of YUC genes associated with leaf wrinkling trait in Tacai variety of Chinese cabbage

Chinese cabbage (Brassica campestris L. ssp. chinensis (L.) Makino) stands as a widely cultivated leafy vegetable in China, with its leaf morphology significantly influencing both quality and yield. Despite its agricultural importance, the precise mechanisms governing leaf wrinkling development remain elusive. This investigation focuses on 'Wutacai', a representative cultivar of the Tacai variety (Brassica campestris L. ssp. chinensis var. rosularis Tsen et Lee), renowned for its distinct leaf wrinkling characteristics. Within the genome of 'Wutacai', we identified a total of 18 YUCs, designated as BraWTC_YUCs, revealing their conservation within the Brassica genus, and their close homology to YUCs in Arabidopsis. Expression profiling unveiled that BraWTC_YUCs in Chinese Cabbage exhibited organ-specific and leaf position-dependent variation. Additionally, transcriptome sequencing data from the flat leaf cultivar 'Suzhouqing' and the wrinkled leaf cultivar 'Wutacai' revealed differentially expressed genes (DEGs) related to auxin during the early phases of leaf development, particularly the YUC gene. In summary, this study successfully identified the YUC gene family in 'Wutacai' and elucidated its potential function in leaf wrinkling trait, to provide valuable insights into the prospective molecular mechanisms that regulate leaf wrinkling in Chinese cabbage.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Functional Differentiation of the Duplicated Gene BrrCIPK9 in Turnip (Brassica rapa var. rapa)

Gene duplication is a key biological process in the evolutionary history of plants and an important driving force for the diversification of genomic and genetic systems. Interactions between the calcium sensor calcineurin B-like protein (CBL) and its target, CBL-interacting protein kinase (CIPK), play important roles in the plant's response to various environmental stresses. As a food crop with important economic and research value, turnip (Brassica rapa var. rapa) has been well adapted to the environment of the Tibetan Plateau and become a traditional crop in the region. The BrrCIPK9 gene in turnip has not been characterized. In this study, two duplicated genes, BrrCIPK9.1 and BrrCIPK9.2, were screened from the turnip genome. Based on the phylogenetic analysis, BrrCIPK9.1 and BrrCIPK9.2 were found located in different sub-branches on the phylogenetic tree. Real-time fluorescence quantitative PCR analyses revealed their differential expression levels between the leaves and roots and in response to various stress treatments. The differences in their interactions with BrrCBLs were also revealed by yeast two-hybrid analyses. The results indicate that BrrCIPK9.1 and BrrCIPK9.2 have undergone Asparagine-alanine-phenylalanine (NAF) site divergence during turnip evolution, which has resulted in functional differences between them. Furthermore, BrrCIPK9.1 responded to high-pH (pH 8.5) stress, while BrrCIPK9.2 retained its ancestral function (low K+), thus providing further evidence of their functional divergence. These functional divergence genes facilitate turnip's good adaptation to the extreme environment of the Tibetan Plateau. In summary, the results of this study reveal the characteristics of the duplicated BrrCIPK9 genes and provide a basis for further functional studies of BrrCBLs-BrrCIPKs in turnip.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Large-scale gene expression alterations introduced by structural variation drive morphotype diversification in Brassica oleracea

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.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Molecular evolution analysis of MYB5 in Brassicaceae with specific focus on seed coat color of Brassica napus

BACKGROUND

MYB transcription factors are splay a vital role in plant biology, with previous research highlighting the significant impact of the R2R3-MYB-like transcription factor MYB5 on seed mucilage biosynthesis, trichome branching, and seed coat development. However, there is a dearth of studies investigating its role in the regulation of proanthocyanidin (PA) biosynthesis.

RESULTS

In this study, a total of 51 MYB5 homologous genes were identified across 31 species belonging to the Brassicaceae family, with particular emphasis on Brassica napus for subsequent investigation. Through phylogenetic analysis, these genes were categorized into four distinct subclasses. Protein sequence similarity and identity analysis demonstrated a high degree of conservation of MYB5 among species within the Brassicaceae family. Additionally, the examination of selection pressure revealed that MYB5 predominantly underwent purifying selection during its evolutionary history, as indicated by the Ka/Ks values of all MYB5 homologous gene pairs being less than one. Notably, we observed a higher rate of non-synonymous mutations in orthologous genes compared to paralogous genes, and the Ka/Ks value displayed a stronger correlation with Ka. In B. napus, an examination of expression patterns in five tissues revealed that MYB5 exhibited particularly high expression in the black seed coat. The findings from the WGCNA demonstrated a robust correlation between MYB5 and BAN(ANR) associated with PA biosynthesis in the black seed coat, providing further evidence of their close association and co-expression. Furthermore, the results obtained from of the analysis of protein interaction networks offer supplementary support for the proposition that MYB5 possesses the capability to interact with transcriptional regulatory proteins, specifically TT8 and TT2, alongside catalytic enzymes implicated in the synthesis of PAs, thereby making a contribution to the biosynthesis of PAs. These findings imply a plausible and significant correlation between the nuique expression pattern of MYB5 and the pigmentation of rapeseed coats. Nevertheless, additional research endeavors are imperative to authenticate and substantiate these findings.

CONCLUSIONS

This study offers valuable insights into the genetic evolution of Brassicaceae plants, thereby serving as a significant reference for the genetic enhancement of Brassicaceae seed coat color.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

Two-Component System Genes in Brassica napus: Identification, Analysis, and Expression Patterns in Response to Abiotic and Biotic Stresses

The two-component system (TCS), consisting of histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs) in eukaryotes, plays pivotal roles in regulating plant growth, development, and responses to environment stimuli. However, the TCS genes were poorly characterized in rapeseed, which is an important tetraploid crop in Brassicaceae. In this work, a total of 182 BnaTCS genes were identified, including 43 HKs, 16 HPs, and 123 RRs, which was more than that in other crops due to segmental duplications during the process of polyploidization. It was significantly different in genetic diversity between the three subfamilies, and some members showed substantial genetic differentiation among the three rapeseed ecotypes. Several hormone- and stress-responsive cis-elements were identified in the putative promoter regions of BnaTCS genes. Furthermore, the expression of BnaTCS genes under abiotic stresses, exogenous phytohormone, and biotic stresses was analyzed, and numerous candidate stress-responsive genes were screened out. Meanwhile, using a natural population with 505 B. napus accessions, we explored the genetic effects of BnaTCS genes on salt tolerance by association mapping analysis and detected some significant association SNPs/genes. The result will help to further understand the functions of TCS genes in the developmental and stress tolerance improvement in B. napus.

Chasing the mechanisms of ecologically adaptive salinity tolerance

Plants adapted to challenging environments offer fascinating models of evolutionary change. Importantly, they also give information to meet our pressing need to develop resilient, low-input crops. With mounting environmental fluctuation-including temperature, rainfall, and soil salinity and degradation-this is more urgent than ever. Happily, solutions are hiding in plain sight: the adaptive mechanisms from natural adapted populations, once understood, can then be leveraged. Much recent insight has come from the study of salinity, a widespread factor limiting productivity, with estimates of 20% of all cultivated lands affected. This is an expanding problem, given increasing climate volatility, rising sea levels, and poor irrigation practices. We therefore highlight recent benchmark studies of ecologically adaptive salt tolerance in plants, assessing macro- and microevolutionary mechanisms, and the recently recognized role of ploidy and the microbiome on salinity adaptation. We synthesize insight specifically on naturally evolved adaptive salt-tolerance mechanisms, as these works move substantially beyond traditional mutant or knockout studies, to show how evolution can nimbly "tweak" plant physiology to optimize function. We then point to future directions to advance this field that intersect evolutionary biology, abiotic-stress tolerance, breeding, and molecular plant physiology.

Deciphering peanut complex genomes paves a way to understand its origin and domestication

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

A complex locus regulates highly lobed-leaf formation in Brassica juncea

KEY MESSAGE

Lineage-specific evolution of RCO was described in Brassicaceae. BjRCO.1 and BjRCO.2 within the complex locus regulated highly lobed-leaf formation in Brassica juncea. RCO regulates the formation of lobed leaves in Brassicaceae species. RCO originated from the duplication of LMI1-type sequences and evolved through gene duplication and loss within the Brassicaceae. However, the evolutionary process and diversification of RCO in different lineages of Brassicaceae remain unclear. Although the RCO locus in B. juncea has been associated with lobed-leaf formation, its complexity has remained largely unknown. This study involved the identification of 55 LMI1-like genes in 16 species of Brassicaceae through syntenic analysis. We classified these LMI1-like genes into two types, namely LMI1-type and RCO-type, based on their phylogenetic relationship. Additionally, we proposed two independent lineage-specific evolution routes for RCO following the divergence of Aethionema. Our findings revealed that the LMI1-like loci responsible for lobed-leaf formation in Brassica species are located on the LF subgenomes. For B. juncea (T84-66V2), we discovered that the complex locus underwent duplication through segments of nucleic acid sequence containing Exostosin-LMI1-RCO (E-R-L), resulting in the tandem presence of two RCO-type and two LMI1-type genes on chromosome A10. As additional evidence, we successfully mapped the complex locus responsible for highly lobed-leaf formation to chromosome A10 using a B. juncea F2 population, which corroborated the results of our evolutionary analysis. Furthermore, through transcriptome analysis, we clarified that BjRCO.1 and BjRCO.2 within the complex locus are functional genes involved in the regulation of highly lobed-leaf formation. The findings of this study offer valuable insights into the regulation of leaf morphology for the breeding of Brassica crops.

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11 Mb, with an N50 scaffold of 46.20 Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40 Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

Pan-genome analysis sheds light on structural variation-based dissection of agronomic traits in melon crops

Sweetness and appearance of fresh fruits are key palatable and preference attributes for consumers and are often controlled by multiple genes. However, fine-mapping the key loci or genes of interest by single genome-based genetic analysis is challenging. Herein, we present the chromosome-level genome assembly of 1 landrace melon accession (Cucumis melo ssp. agrestis) with wild morphologic features and thus construct a melon pan-genome atlas via integrating sequenced melon genome datasets. Our comparative genomic analysis reveals a total of 3.4 million genetic variations, of which the presence/absence variations (PAVs) are mainly involved in regulating the function of genes for sucrose metabolism during melon domestication and improvement. We further resolved several loci that are accountable for sucrose contents, flesh color, rind stripe, and suture using a structural variation (SV)-based genome-wide association study. Furthermore, via bulked segregation analysis (BSA)-seq and map-based cloning, we uncovered that a single gene, (CmPIRL6), determines the edible or inedible characteristics of melon fruit exocarp. These findings provide important melon pan-genome information and provide a powerful toolkit for future pan-genome-informed cultivar breeding of melon.

Role of BraRGL1 in regulation of Brassica rapa bolting and flowering

Gibberellin (GA) plays a major role in controlling Brassica rapa stalk development. As an essential negative regulator of GA signal transduction, DELLA proteins may exert significant effects on stalk development. However, the regulatory mechanisms underlying this regulation remain unclear. In this study, we report highly efficient and inheritable mutagenesis using the CRISPR/Cas9 gene editing system in BraPDS (phytoene desaturase) and BraRGL1 (key DELLA protein) genes. We observed a loss-of-function mutation in BraRGL1 due to two amino acids in GRAS domain. The flower bud differentiation and bolting time of BraRGL1 mutants were significantly advanced. The expression of GA-regulatory protein (BraGASA6), flowering related genes (BraSOC1, BraLFY), expansion protein (BraEXPA11) and xyloglucan endotransferase (BraXTH3) genes was also significantly upregulated in these mutants. BraRGL1-overexpressing plants displayed the contrasting phenotypes. BraRGL1 mutants were more sensitive to GA signaling. BraRGL1 interacted with BraSOC1, and the interaction intensity decreased after GA3 treatment. In addition, BraRGL1 inhibited the transcription-activation ability of BraSOC1 for BraXTH3 and BraLFY genes, but the presence of GA3 enhanced the activation ability of BraSOC1, suggesting that the BraRGL1-BraSOC1 module regulates bolting and flowering of B. rapa through GA signal transduction. Thus, we hypothesized that BraRGL1 is degraded, and BraSOC1 is released in the presence of GA3, which promotes the expression of BraXTH3 and BraLFY, thereby inducing stalk development in B. rapa. Further, the BraRGL1-M mutant promoted the flower bud differentiation without affecting the stalk quality. Thus, BraRGL1 can serve as a valuable target for the molecular breeding of early maturing varieties.

A BrLINE1-RUP insertion in BrCER2 alters cuticular wax biosynthesis in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Glossiness is an important quality-related trait of Chinese cabbage, which is a leafy vegetable crop in the family Brassicaceae. The glossy trait is caused by abnormal cuticular wax accumulation. In this study, on the basis of a bulked segregant analysis coupled with next-generation sequencing (BSA-seq) and fine-mapping, the most likely candidate gene responsible for the glossy phenotype of Chinese cabbage was identified. It was subsequently named Brcer2 because it is homologous to AtCER2 (At4g24510). A bioinformatics analysis indicated a long interspersed nuclear element 1 (LINE-1) transposable element (named BrLINE1-RUP) was inserted into the first exon of Brcer2 in HN19-G via an insertion-mediated deletion mechanism, which introduced a premature termination codon. Gene expression analysis showed that the InDel mutation of BrCER2 reduced the transcriptional expression levels of Brcer2 in HN19-G. An analysis of cuticular waxes suggested that a loss-of-function mutation to BrCER2 in Chinese cabbage leads to a severe decrease in the abundance of very-long-chain-fatty-acids (> C28), resulting in the production of a cauline leaf, inflorescence stem, flower, and pistil with a glossy phenotype. These findings imply the insertion of the LINE-1 transposable element BrLINE1-RUP into BrCER2 can modulate the waxy traits of Chinese cabbage plants.

Role of Brassica orphan gene BrLFM on leafy head formation in Chinese cabbage (Brassica rapa)

Brassica orphan gene BrFLM, identified by two allelic mutants, was involved in leafy head formation in Chinese cabbage. Leafy head formation is a unique agronomic trait of Chinese cabbage that determines its yield and quality. In our previous study, an EMS mutagenesis Chinese cabbage mutant library was constructed using the heading Chinese cabbage double haploid (DH) line FT as the wild-type. Here, we screened two extremely similar leafy head deficiency mutants lfm-1 and lfm-2 with geotropic growth leaves from the library to investigate the gene(s) related to leafy head formation. Reciprocal crossing results showed that these two mutants were allelic. We utilized lfm-1 to identify the mutant gene(s). Genetic analysis showed that the mutated trait was controlled by a single nuclear gene Brlfm. Mutmap analysis showed that Brlfm was located on chromosome A05, and BraA05g012440.3C or BraA05g021450.3C were the candidate gene. Kompetitive allele-specific PCR analysis eliminated BraA05g012440.3C from the candidates. Sanger sequencing identified an SNP from G to A at the 271st nucleotide on BraA05g021450.3C. The sequencing of lfm-2 detected another non-synonymous SNP (G to A) located at the 266st nucleotide on BraA05g021450.3C, which verified its function on leafy head formation. We blasted BraA05g021450.3C on database and found that it belongs to a Brassica orphan gene encoding an unknown 13.74 kDa protein, named BrLFM. Subcellular localization showed that BrLFM was located in the nucleus. These findings reveal that BrLFM is involved in leafy head formation in Chinese cabbage.

The NLRomes of Zea mays NAM founder lines and Zea luxurians display presence-absence variation, integrated domain diversity, and mobility

Plant pathogens cause significant crop loss worldwide, and new resistance genes deployed to combat diseases can be overcome quickly. Understanding the existing resistance gene diversity within the germplasm of major crops, such as maize, is crucial for the development of new disease-resistant varieties. We analysed the nucleotide-binding leucine-rich repeat receptors (NLRs) of 26 recently sequenced diverse founder lines from the maize nested association mapping (NAM) population and compared them to the R gene complement present in a wild relative of maize, Zea luxurians. We found that NLRs in both species contain a large diversity of atypical integrated domains, including many domains that have not previously been found in the NLRs of other species. Additionally, the single Z. luxurians genome was found to have greater integrated atypical domain diversity than all 26 NAM founder lines combined, indicating that this species may represent a rich source of novel resistance genes. NLRs were also found to have very high sequence diversity and presence-absence variation among the NAM founder lines, with a large NLR cluster on Chr10 representing a diversity hotspot. Additionally, NLRs were shown to be mobile within maize genomes, with several putative interchromosomal translocations identified.

A High-Continuity Genome Assembly of Chinese Flowering Cabbage (Brassica rapa var. parachinensis) Provides New Insights into Brassica Genome Structure Evolution

Chinese flowering cabbage (Brassica rapa var. parachinensis) is a popular and widely cultivated leaf vegetable crop in Asia. Here, we performed a high quality de novo assembly of the 384 Mb genome of 10 chromosomes of a typical cultivar of Chinese flowering cabbage with an integrated approach using PacBio, Illumina, and Hi-C technology. We modeled 47,598 protein-coding genes in this analysis and annotated 52% (205.9/384) of its genome as repetitive sequences including 17% in DNA transposons and 22% in long terminal retrotransposons (LTRs). Phylogenetic analysis reveals the genome of the Chinese flowering cabbage has a closer evolutionary relationship with the AA diploid progenitor of the allotetraploid species, Brassica juncea. Comparative genomic analysis of Brassica species with different subgenome types (A, B and C) reveals that the pericentromeric regions on chromosome 5 and 6 of the AA genome have been significantly expanded compared to the orthologous genomic regions in the BB and CC genomes, largely driven by LTR-retrotransposon amplification. Furthermore, we identified a large number of structural variations (SVs) within the B. rapa lines that could impact coding genes, suggesting the functional significance of SVs on Brassica genome evolution. Overall, our high-quality genome assembly of the Chinese flowering cabbage provides a valuable genetic resource for deciphering the genome evolution of Brassica species and it can potentially serve as the reference genome guiding the molecular breeding practice of B. rapa crops.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.

Resistance strategies for defense against Albugo candida causing white rust disease

Albugo candida, the causal organism of white rust, is an oomycete obligate pathogen infecting crops of Brassicaceae family occurred on aerial part, including vegetable and oilseed crops at all growth stages. The disease expression is characterized by local infection appearing on the abaxial region developing white or creamy yellow blister (sori) on leaves and systemic infections cause hypertrophy and hyperplasia leading to stag-head of reproductive organ. To overcome this problem, several disease management strategies like fungicide treatments were used in the field and disease-resistant varieties have also been developed using conventional and molecular breeding. Due to high variability among A. candida isolates, there is no single approach available to understand the diverse spectrum of disease symptoms. In absence of resistance sources against pathogen, repetitive cultivation of genetically-similar varieties locally tends to attract oomycete pathogen causing heavy yield losses. In the present review, a deep insight into the underlying role of the non-host resistance (NHR) defence mechanism available in plants, and the strategies to exploit available gene pools from plant species that are non-host to A. candida could serve as novel sources of resistance. This work summaries the current knowledge pertaining to the resistance sources available in non-host germ plasm, the understanding of defence mechanisms and the advance strategies covers molecular, biochemical and nature-based solutions in protecting Brassica crops from white rust disease.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

Structural variation (SV)-based pan-genome and GWAS reveal the impacts of SVs on the speciation and diversification of allotetraploid cottons

Structural variations (SVs) have long been described as being involved in the origin, adaption, and domestication of species. However, the underlying genetic and genomic mechanisms are poorly understood. Here, we report a high-quality genome assembly of Gossypium barbadense acc. Tanguis, a landrace that is closely related to formation of extra-long-staple (ELS) cultivated cotton. An SV-based pan-genome (Pan-SV) was then constructed using a total of 182 593 non-redundant SVs, including 2236 inversions, 97 398 insertions, and 82 959 deletions from 11 assembled genomes of allopolyploid cotton. The utility of this Pan-SV was then demonstrated through population structure analysis and genome-wide association studies (GWASs). Using segregation mapping populations produced through crossing ELS cotton and the landrace along with an SV-based GWAS, certain SVs responsible for speciation, domestication, and improvement in tetraploid cottons were identified. Importantly, some of the SVs presently identified as associated with the yield and fiber quality improvement had not been identified in previous SNP-based GWAS. In particular, a 9-bp insertion or deletion was found to associate with elimination of the interspecific reproductive isolation between Gossypium hirsutum and G. barbadense. Collectively, this study provides new insights into genome-wide, gene-scale SVs linked to important agronomic traits in a major crop species and highlights the importance of SVs during the speciation, domestication, and improvement of cultivated crop species.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

Genome assembly of the Brassicaceae diploid Orychophragmus violaceus reveals complex whole-genome duplication and evolution of dihydroxy fatty acid metabolism

Orychophragmus violaceus is a Brassicaceae species widely cultivated in China, particularly as a winter cover crop in northern China because of its low-temperature tolerance and low water demand. Recently, O. violaceus has also been cultivated as a potential industrial oilseed crop because of its abundant 24-carbon dihydroxy fatty acids (diOH-FAs), which contribute to superior high-temperature lubricant properties. In this study, we performed de novo assembly of the O. violaceus genome. Whole-genome synteny analysis of the genomes of its relatives demonstrated that O. violaceus is a diploid that has undergone an extra whole-genome duplication (WGD) after the Brassicaceae-specific α-WGD event, with a basic chromosome number of x = 12. Formation of diOH-FAs is hypothesized to have occurred after the WGD event. Based on the genome and the transcriptome data from multiple stages of seed development, we predicted that OvDGAT1-1 and OvDGAT1-2 are candidate genes for the regulation of diOH-FA storage in O. violaceus seeds. These results may greatly facilitate the development of heat-tolerant and eco-friendly plant-based lubricants using O. violaceus seed oil and improve our understanding of the genomic evolution of Brassicaceae.

Brassica rapa orphan gene BR1 delays flowering time in Arabidopsis

Orphan genes are essential to the emergence of species-specific traits and the process of evolution, lacking sequence similarity to any other identified genes. As they lack recognizable domains or functional motifs, however, efforts to characterize these orphan genes are often difficult. Flowering is a key trait in Brassica rapa, as premature bolting can have a pronounced adverse impact on plant quality and yield. Bolting resistance-related orphan genes, however, have yet to be characterized. In this study, an orphan gene designated BOLTING RESISTANCE 1 (BR1) was identified and found through gene structural variation analyses to be more highly conserved in Chinese cabbage than in other available accessions. The expression of BR1 was increased in bolting resistant Chinese cabbage and decreased in bolting non-resistant type, and the expression of some mark genes were consist with bolting resistance phenotype. BR1 is primarily expressed in leaves at the vegetative growth stage, and the highest BR1 expression levels during the flowering stage were observed in the flower buds and silique as compared to other tissue types. The overexpression of BR1 in Arabidopsis was associated with enhanced bolting resistance under long day (LD) conditions, with these transgenic plants exhibiting significant decreases in stem height, rosette radius, and chlorophyll content. Transcriptomic sequencing of WT and BR1OE plants showed the association of BR1 with other bolting resistance genes. Transcriptomic sequencing and qPCR revealed that six flowering integrator genes and one chlorophyll biosynthesis-related gene were downregulated following BR1 overexpression. Six key genes in photoperiodic flowering pathway exhibited downward expression trends in BR1OE plants, while the expression of floral repressor AtFLC gene was upregulated. The transcripts of these key genes were consistent with observed phenotypes in BR1OE plants, and the results indicated that BR1 may function through vernalization and photoperiodic pathway. Instead, the protein encoded by BR1 gene was subsequently found to localize to the nucleus. Taken together, we first propose that orphan gene BR1 functions as a novel regulator of flowering time, and these results suggested that BR1 may represent a promising candidate gene to support the selective breeding of Chinese cabbage cultivars with enhanced bolting resistance.

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Whole-Genome Comparison Reveals Structural Variations behind Heading Leaf Trait in Brassica oleracea

Brassica oleracea displays remarkable morphological variations. It intrigued researchers to study the underlying cause of the enormous diversification of this organism. However, genomic variations in complex heading traits are less known in B. oleracea. Herein, we performed a comparative population genomics analysis to explore structural variations (SVs) responsible for heading trait formation in B. oleracea. Synteny analysis showed that chromosomes C1 and C2 of B. oleracea (CC) shared strong collinearity with A01 and A02 of B. rapa (AA), respectively. Two historical events, whole genome triplication (WGT) of Brassica species and differentiation time between AA and CC genomes, were observed clearly by phylogenetic and Ks analysis. By comparing heading and non-heading populations of B. oleracea genomes, we found extensive SVs during the diversification of the B. oleracea genome. We identified 1205 SVs that have an impact on 545 genes and might be associated with the heading trait of cabbage. Overlapping the genes affected by SVs and the differentially expressed genes identified by RNA-seq analysis, we identified six vital candidate genes that may be related to heading trait formation in cabbage. Further, qRT-PCR experiments also verified that six genes were differentially expressed between heading leaves and non-heading leaves, respectively. Collectively, we used available genomes to conduct a comparison population genome analysis and identify candidate genes for the heading trait of cabbage, which provides insight into the underlying reason for heading trait formation in B. oleracea.

Stigma receptors control intraspecies and interspecies barriers in Brassicaceae

Flowering plants have evolved numerous intraspecific and interspecific prezygotic reproductive barriers to prevent production of unfavourable offspring¹. Within a species, self-incompatibility (SI) is a widely utilized mechanism that rejects self-pollen^2,3 to avoid inbreeding depression. Interspecific barriers restrain breeding between species and often follow the SI × self-compatible (SC) rule, that is, interspecific pollen is unilaterally incompatible (UI) on SI pistils but unilaterally compatible (UC) on SC pistils^1,4-6. The molecular mechanisms underlying SI, UI, SC and UC and their interconnections in the Brassicaceae remain unclear. Here we demonstrate that the SI pollen determinant S-locus cysteine-rich protein/S-locus protein 11 (SCR/SP11)^2,3 or a signal from UI pollen binds to the SI female determinant S-locus receptor kinase (SRK)^2,3, recruits FERONIA (FER)^7-9 and activates FER-mediated reactive oxygen species production in SI stigmas^10,11 to reject incompatible pollen. For compatible responses, diverged pollen coat protein B-class^12-14 from SC and UC pollen differentially trigger nitric oxide, nitrosate FER to suppress reactive oxygen species in SC stigmas to facilitate pollen growth in an intraspecies-preferential manner, maintaining species integrity. Our results show that SRK and FER integrate mechanisms underlying intraspecific and interspecific barriers and offer paths to achieve distant breeding in Brassicaceae crops.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Identification of NPF Family Genes in Brassica rapa Reveal Their Potential Functions in Pollen Development and Response to Low Nitrate Stress

Nitrate Transporter 1/Peptide Transporter Family (NPF) genes encode membrane transporters involved in the transport of diverse substrates. However, little is known about the diversity and functions of NPFs in Brassica rapa. In this study, 85 NPFs were identified in B. rapa (BrNPFs) which comprised eight subfamilies. Gene structure and conserved motif analysis suggested that BrNFPs were conserved throughout the genus. Stress and hormone-responsive cis-acting elements and transcription factor binding sites were identified in BrNPF promoters. Syntenic analysis suggested that tandem duplication contributed to the expansion of BrNPFs in B. rapa. Transcriptomic profiling analysis indicated that BrNPF2.6, BrNPF2.15, BrNPF7.6, and BrNPF8.9 were expressed in fertile floral buds, suggesting important roles in pollen development. Thirty-nine BrNPFs were responsive to low nitrate availability in shoots or roots. BrNPF2.10, BrNPF2.19, BrNPF2.3, BrNPF5.12, BrNPF5.16, BrNPF5.8, and BrNPF6.3 were only up-regulated in roots under low nitrate conditions, indicating that they play positive roles in nitrate absorption. Furthermore, many genes were identified in contrasting genotypes that responded to vernalization and clubroot disease. Our results increase understanding of BrNPFs as candidate genes for genetic improvement studies of B. rapa to promote low nitrate availability tolerance and for generating sterile male lines based on gene editing methods.

Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium

Phenotypic diversity and evolutionary innovation ultimately trace to variation in genomic sequence and rewiring of regulatory networks. Here, we constructed a pan-genome of the Gossypium genus using ten representative diploid genomes. We document the genomic evolutionary history and the impact of lineage-specific transposon amplification on differential genome composition. The pan-3D genome reveals evolutionary connections between transposon-driven genome size variation and both higher-order chromatin structure reorganization and the rewiring of chromatin interactome. We linked changes in chromatin structures to phenotypic differences in cotton fiber and identified regulatory variations that decode the genetic basis of fiber length, the latter enabled by sequencing 1,005 transcriptomes during fiber development. We showcase how pan-genomic, pan-3D genomic and genetic regulatory data serve as a resource for delineating the evolutionary basis of spinnable cotton fiber. Our work provides insights into the evolution of genome organization and regulation and will inform cotton improvement by enabling regulome-based approaches.

Genome-wide identification and expression analysis of BrAGC genes in Brassica rapa reveal their potential roles in sexual reproduction and abiotic stress tolerance

AGC protein kinases play important roles in regulating plant growth, immunity, and cell death. However, the function of AGC in Brassica rapa has not yet been clarified. In this study, 62 BrAGC genes were identified, and these genes were distributed on 10 chromosomes and divided into six subfamilies. Analysis of gene structure and conserved motifs showed that the activation segment of BrAGC genes was highly conserved, and genes of the same subfamily showed higher sequence and structural similarity. Collinearity analysis revealed that BrAGCs were more closely related to AtAGCs than to OsAGCs. Expression profile analysis revealed that BrAGCs were preferentially expressed in flowers and BrAGC26, BrAGC33, and BrAGC04 were preferentially expressed in the stigma; the expression of these genes was significantly upregulated after self-incompatibility pollination, and the expression of BrAGC13 and BrAGC32 was significantly upregulated after cross-pollination. In addition, several typical cis-elements involved in the stress response were identified in BrAGC promoters. The expression levels of BrAGC37 and BrAGC44 significantly varied under different types of abiotic stress. Collectively, we identified that BrAGC26, BrAGC33, and BrAGC44 have the greatest potential in regulating pollen-pistil interaction and abiotic stress tolerance, respectively. Our findings will aid future functional investigations of BrAGCs in B. rapa.

The Global Assessment of Oilseed Brassica Crop Species Yield, Yield Stability and the Underlying Genetics

The global demand for oilseeds is increasing along with the human population. The family of Brassicaceae crops are no exception, typically harvested as a valuable source of oil, rich in beneficial molecules important for human health. The global capacity for improving Brassica yield has steadily risen over the last 50 years, with the major crop Brassica napus (rapeseed, canola) production increasing to ~72 Gt in 2020. In contrast, the production of Brassica mustard crops has fluctuated, rarely improving in farming efficiency. The drastic increase in global yield of B. napus is largely due to the demand for a stable source of cooking oil. Furthermore, with the adoption of highly efficient farming techniques, yield enhancement programs, breeding programs, the integration of high-throughput phenotyping technology and establishing the underlying genetics, B. napus yields have increased by >450 fold since 1978. Yield stability has been improved with new management strategies targeting diseases and pests, as well as by understanding the complex interaction of environment, phenotype and genotype. This review assesses the global yield and yield stability of agriculturally important oilseed Brassica species and discusses how contemporary farming and genetic techniques have driven improvements.

Genome-wide identification and expression analysis of WRKY family genes under soft rot in Chinese cabbage

Complex transcriptional networks regulate plant defense against pathogen attack, and plant transcription factors act as key regulators of the plant immune responses. The differences between transcription factor expression and regulation in Chinese cabbage soft rot (Pectobacterium carotovorum; Pc) have not been revealed. In this study, a total of 148 putative Chinese cabbage WRKY genes (BrWRKYs) were identified from the Chinese cabbage genome (v3.0). These genes were divided into seven subgroups (groups I, IIa–e, and III) based on phylogenomic analysis, with distinct motif compositions in each subgroup. Time-series RNA-seq was carried out to elucidate the dynamic expression patterns of the BrWRKYs on the resistant mutant (sr) and the susceptible wild-type (inbred WT) challenged by Pc. Transcriptional analysis showed that 48 WRKY transcription genes at 0–24 hpi were significantly upregulated in sr under soft rot stress. At the 12-h post-inoculation critical time point, we identified three specifically upregulated genes and two downregulated genes in the resistant mutant, which may provide potential applications for genetic improvement against soft rot. The findings improved our understanding of the WRKY-mediated soft rot stress response regulation in Chinese cabbage. The study thus lays a foundation for the genetic improvement of soft rot resistance.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Comprehensive In Silico Characterization and Expression Pro-Filing of DA1/DAR Family Genes in Brassica rapa

The DA1/DAR family genes have been shown to play important roles in regulating organ size and plant biomass in the model plant Arabidopsis and several crops. However, this family has not been characterized in Brassica rapa (B. rapa). In this study, we identified 17 DA1&DAR genes from B. rapa. Phylogenetic analysis indicated that these genes are classified into four groups. Structural and motif analysis of BrDA1&DARs discovered that the genes within the same group have similar exon-intron structures and share an equal number of conserved motifs except for BrDAR6.3 from group IV, which contains two conserved motifs. Cis-regulatory elements identified four phytohormones (Salicylic acid, Abscisic acid, Gibberellin, and auxin) and three major abiotic (Light, Low temperature, and drought) responsive elements. Further, six br-miRNAs named br-miR164a, br-miR164b, br-miR164c, br-miR164d, br-miRN360, and br-miRN366 were found which target BrDAR6.1, BrDA1.4, and BrDA1.5. BrDA1&DAR genes were highly expressed in stem, root, silique, flower, leaf, and callus tissues. Moreover, qRT-PCR analyses indicated that some of these genes were responsive to abiotic stresses or phytohormone treatments. Our findings provide a foundation for further genetic and physiological studies of BrDA1&DARs in B. rapa.