Evidence for two domestication lineages supporting a middle-eastern origin for Brassica oleracea crops from diversified kale populations

Metabolomic and Transcriptomic Profiles in Diverse Brassica oleracea Crops Provide Insights into the Genetic Regulation of Glucosinolate Profiles

Glucosinolates (GSLs) are plant secondary metabolites commonly found in the cruciferous vegetables of the Brassicaceae family, offering health benefits to humans and defense against pathogens and pests to plants. In this study, we investigated 23 GSL compounds' relative abundance in four tissues of five different Brassica oleracea morphotypes. Using the five corresponding high-quality B. oleracea genome assemblies, we identified 183 GSL-related genes and analyzed their expression with mRNA-Seq data. GSL abundance and composition varied strongly, among both tissues and morphotypes, accompanied by different gene expression patterns. Interestingly, broccoli exhibited a nonfunctional AOP2 gene due to a conserved 2OG-FeII_Oxy domain loss, explaining the unique accumulation of two health-promoting GSLs. Additionally, transposable element (TE) insertions were found to affect the gene structure of MAM3 genes. Our findings deepen the understanding of GSL variation and genetic regulation in B. oleracea morphotypes, providing valuable insights for breeding with tailored GSL profiles in these crops.

Polyphenols in edible plant leaves: an overview of their occurrence and health properties

Edible plant leaves (EPLs) constitute a major renewable functional plant biomass available all year round, providing an essential source of polyphenols in the global diet. Polyphenols form a large family of antioxidant molecules. They protect against the harmful effects of free radicals, strengthen immunity and stimulate the body's natural defenses thanks to their antibacterial and antiviral functions. This study refers to phenolic compounds from 50 edible plant leaves divided into four categories: green leafy vegetables, underutilized leafy vegetables, leafy spices and leafy drinks. It provides data on the identification, occurrence and pharmacological functions of polyphenols contained in EPLs, and provides a better understanding of trends and gaps in their consumption and study. Certain EPLs, such as moringa (Moringa oleifera Lam.), tea (Camellia sinensis L.) and several leafy spices of the Lamiaceae family, reveal important characteristics and therapeutic potential. The polyphenol composition of EPLs makes them functional plants that offer relevant solutions in the fight against obesity, the management of food insecurity and the prevention of chronic diseases.

Genomic analyses reveal the stepwise domestication and genetic mechanism of curd biogenesis in cauliflower

Cauliflower (Brassica oleracea L. var. botrytis) is a distinctive vegetable that supplies a nutrient-rich edible inflorescence meristem for the human diet. However, the genomic bases of its selective breeding have not been studied extensively. Herein, we present a high-quality reference genome assembly C-8 (V2) and a comprehensive genomic variation map consisting of 971 diverse accessions of cauliflower and its relatives. Genomic selection analysis and deep-mined divergences were used to explore a stepwise domestication process for cauliflower that initially evolved from broccoli (Curd-emergence and Curd-improvement), revealing that three MADS-box genes, CAULIFLOWER1 (CAL1), CAL2 and FRUITFULL (FUL2), could have essential roles during curd formation. Genome-wide association studies identified nine loci significantly associated with morphological and biological characters and demonstrated that a zinc-finger protein (BOB06G135460) positively regulates stem height in cauliflower. This study offers valuable genomic resources for better understanding the genetic bases of curd biogenesis and florescent development in crops.

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.

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.

Cabbage (Brassica oleracea var. capitata) Development in Time: How Differential Parenchyma Tissue Growth Affects Leafy Head Formation

This study aims to categorize the morphological changes during cabbage (B. oleracea ssp. capitata) development, seedling, rosette, folding, and heading, and to elucidate the cellular mechanisms of the leaf curvature, essential for the formation of the leafy head. We followed the growth of two cabbage cultivars with distinct head shapes (round and pointed) and one non-heading collard cultivar; we phenotyped the size and volume of the whole plant as well as the size, shape, and curvature of the leaves during growth. By integrating these phenotypic data, we determined the four vegetative stages for both cabbages. The histological phenotypes of microtome sections from five distinct leaf positions of the rosette, folding, and heading leaves at two timepoints during leaf growth were quantified and revealed variations in cellular parameters among leaf types, between leaf positions, and between the adaxial and abaxial sides. We identified two synergistic cellular mechanisms contributing to the curvature of heading leaves: differential growth across the leaf blade, with increased growth at the leaf's center relative to the margins; and the increased expansion of the spongy parenchyma layer compared to the palisade parenchyma layer, resulting in the direction of the curvature, which is inwards. These two processes together contribute to the typical leafy heads of cabbages.

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.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

Cauliflowers or how the perseverance of a plant to make flowers produces an amazing fractal structure

Biological organisms have an immense diversity of forms. Some of them exhibit conspicuous and fascinating fractal structures that present self-similar patterns at all scales. How such structures are produced by biological processes is intriguing. In a recent publication, we used a multi-scale modelling approach to understand how gene activity can produce macroscopic cauliflower curds. Our work provides a plausible explanation for the appearance of fractal-like structures in plants, linking gene activity with development.

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.

Construction of an SNP fingerprinting database and population genetic analysis of 329 cauliflower cultivars

Cauliflower is one of the most important vegetable crops grown worldwide. However, the lack of genetic diversity information and efficient molecular markers hinders efforts to improve cauliflower. This study aims to construct DNA fingerprints for 329 cauliflower cultivars based on SNP markers and the KASP system. After rigorous filtering, a total of 1662 candidate SNPs were obtained from nearly 17.9 million SNP loci. The mean values of PIC, MAF, heterozygosity and gene diversity of these SNPs were 0.389, 0.419, 0.075, and 0.506, respectively. We developed a program for in silico simulations on 153 core germplasm samples to generate ideal SNP marker sets from the candidates. Finally, 41 highly polymorphic KASP markers were selected and applied to identify 329 cauliflower cultivars, mainly collected from the public market. Furthermore, based on the KASP genotyping data, we performed phylogenetic analysis and population structure analysis of the 329 cultivars. As a result, these cultivars could be classified into three major clusters, and the classification patterns were significantly related to their curd solidity and geographical origin. Finally, fingerprints of the 329 cultivars and 2D barcodes with the genetic information of each sample were generated. The fingerprinting database developed in this study provides a practical tool for identifying the authenticity and purity of cauliflower seeds and valuable genetic information about the current cauliflower cultivars.

Genetic dissection of morphological variation in rosette leaves and leafy heads in cabbage (Brassica oleracea var. capitata)

Correlations between morphological traits of cabbage rosette leaves and heads were found. Genome-wide association studies of these traits identified 50 robust quantitative trait loci in multiple years. Half of these loci affect both organs. Cabbage (Brassica oleracea var. capitata) is an economically important vegetable crop cultivated worldwide. Cabbage plants go through four vegetative stages: seedling, rosette, folding and heading. Rosette leaves are the largest leaves of cabbage plants and provide most of the energy needed to produce the leafy head. To understand the relationship and the genetic basis of leaf development and leafy head formation, 308 cabbage accessions were scored for rosette leaf and head traits in three-year field trials. Significant correlations were found between morphological traits of rosette leaves and heads, namely leaf area with the head area, height and width, and leaf width with the head area and head height, when heads were harvested at a fixed number of days after sowing. Fifty robust quantitative trait loci (QTLs) for rosette leaf and head traits distributed over all nine chromosomes were identified with genome-wide association studies. All these 50 loci were identified in multiple years and generally affect multiple traits. Twenty-five of the QTL were associated with both rosette leaf and leafy head traits. We discuss thirteen candidate genes identified in these QTL that are expressed in heading leaves, with an annotation related to auxin and other phytohormones, leaf development, and leaf polarity that likely play a role in leafy head development or rosette leaf expansion.