Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea

Pan-genome analysis of 13 Spinacia accessions reveals structural variations associated with sex chromosome evolution and domestication traits in spinach

Structural variations (SVs) are major genetic variants that can be involved in the origin, adaptation and domestication of species. However, the identification and characterization of SVs in Spinacia species are rare due to the lack of a pan-genome. Here, we report eight chromosome-scale assemblies of cultivated spinach and its two wild species. After integration with five existing assemblies, we constructed a comprehensive Spinacia pan-genome and identified 193 661 pan-SVs, which were genotyped in 452 Spinacia accessions. Our pan-SVs enabled genome-wide association study identified signals associated with sex and clarified the evolutionary direction of spinach. Most sex-linked SVs (86%) were biased to occur on the Y chromosome during the evolution of the sex-linked region, resulting in reduced Y-linked gene expression. The frequency of pan-SVs among Spinacia accessions further illustrated the contribution of these SVs to domestication, such as bolting time and seed dormancy. Furthermore, compared with SNPs, pan-SVs act as efficient variants in genomic selection (GS) because of their ability to capture missing heritability information and higher prediction accuracy. Overall, this study provides a valuable resource for spinach genomics and highlights the potential utility of pan-SV in crop improvement and breeding programmes.

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.

Unraveling the mechanism of flower color variation in Brassica napus by integrated metabolome and transcriptome analyses

Brassica napus is one of the most important oil crops in the world. Breeding oilseed rape with colorful flowers can greatly enhance the ornamental value of B. napus and thus improve the economic benefits of planting. As water-soluble flavonoid secondary metabolites, anthocyanins are very important for the synthesis and accumulation of pigments in the petals of plants, giving them a wide range of bright colors. Despite the documentation of over 60 distinct flower shades in B. napus, the intricacies underlying flower color variation remain elusive. Particularly, the mechanisms driving color development across varying flower color backgrounds necessitate further comprehensive investigation. This research undertook a comprehensive exploration through the integration of transcriptome and metabolome analyses to pinpoint pivotal genes and metabolites underpinning an array of flower colors, including beige, beige-red, yellow, orange-red, deep orange-red, white, light-purple, and purple. First, we used a two-way BLAST search to find 275 genes in the reference genome of B. napus Darmor v10 that were involved in making anthocyanins. The subsequent scrutiny of RNA-seq outcomes underscored notable upregulation in the structural genes F3H and UGT, alongside the MYB75, GL3, and TTG1 transcriptional regulators within petals, showing anthocyanin accumulation. By synergizing this data with a weighted gene co-expression network analysis, we identified CHS, F3H, MYB75, MYB12, and MYB111 as the key players driving anthocyanin synthesis in beige-red, orange-red, deep orange-red, light-purple, and purple petals. By integrating transcriptome and weighted gene co-expression network analysis findings with anthocyanin metabolism data, it is hypothesized that the upregulation of MYB75, which, in turn, enhances F3H expression, plays a pivotal role in the development of pigmented oilseed rape flowers. These findings help to understand the transcriptional regulation of anthocyanin biosynthesis in B. napus and provide valuable genetic resources for breeding B. napus varieties with novel flower colors.

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.

Multiomics analysis of a resistant European turnip ECD04 during clubroot infection reveals key hub genes underlying resistance mechanism

The clubroot disease has become a worldwide threat for crucifer crop production, due to its soil-borne nature and difficulty to eradicate completely from contaminated field. In this study we used an elite resistant European fodder turnip ECD04 and investigated its resistance mechanism using transcriptome, sRNA-seq, degradome and gene editing. A total of 1751 DEGs were identified from three time points after infection, among which 7 hub genes including XTH23 for cell wall assembly and two CPK28 genes in PTI pathways. On microRNA, we identified 17 DEMs and predicted 15 miRNA-target pairs (DEM-DEG). We validated two pairs (miR395-APS4 and miR160-ARF) by degradome sequencing. We investigated the miR395-APS4 pair by CRISPR-Cas9 mediated gene editing, the result showed that knocking-out APS4 could lead to elevated clubroot resistance in B. napus. In summary, the data acquired on transcriptional response and microRNA as well as target genes provide future direction especially gene candidates for genetic improvement of clubroot resistance on Brassica species.

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.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

Transcriptomic analyses to summarize gene expression patterns that occur during leaf initiation of Chinese cabbage

In Chinese cabbage, rosette leaves expose their adaxial side to the light converting light energy into chemical energy, acting as a source for the growth of the leafy head. In the leafy head, the outer heading leaves expose their abaxial side to the light while the inner leaves are shielded from the light and have become a sink organ of the growing Chinese cabbage plant. Interestingly, variation in several ad/abaxial polarity genes is associated with the typical leafy head morphotype. The initiation of leaf primordia and the establishment of leaf ad/abaxial polarity are essential steps in the initiation of marginal meristem activity leading to leaf formation. Understanding the molecular genetic mechanisms of leaf primordia formation, polar differentiation, and leaf expansion is thus relevant to understand leafy head formation. As Brassica's are mesa-hexaploids, many genes have multiple paralogues, complicating analysis of the genetic regulation of leaf development. In this study, we used laser dissection of Chinese cabbage leaf primordia and the shoot apical meristem (SAM) to compare gene expression profiles between both adaxial and abaxial sides and the SAM aiming to capture transcriptome changes underlying leaf primordia development. We highlight genes with roles in hormone pathways and transcription factors. We also assessed gene expression gradients along expanded leaf blades from the same plants to analyze regulatory links between SAM, leaf primordia and the expanding rosette leaf. The catalogue of differentially expressed genes provides insights in gene expression patterns involved in leaf development and form a starting point to unravel leafy head formation.

Studying Salt-Induced Shifts in Gene Expression Patterns of Glucosinolate Transporters and Glucosinolate Accumulation in Two Contrasting Brassica Species

Brassica crops are well known for the accumulation of glucosinolates-secondary metabolites crucial for plants' adaptation to various stresses. Glucosinolates also functioning as defence compounds pose challenges to food quality due to their goitrogenic properties. Their disruption leaves plants susceptible to insect pests and diseases. Hence, a targeted reduction in seed glucosinolate content is of paramount importance to increase food acceptance. GLUCOSINOLATE TRANSPORTERS (GTRs) present a promising avenue for selectively reducing glucosinolate concentrations in seeds while preserving biosynthesis elsewhere. In this study, 54 putative GTR protein sequences found in Brassica were retrieved, employing Arabidopsis GTR1 and GTR2 templates. Comprehensive bioinformatics analyses, encompassing gene structure organization, domain analysis, motif assessments, promoter analysis, and cis-regulatory elements, affirmed the existence of transporter domains and stress-related regulatory elements. Phylogenetic analysis revealed patterns of conservation and divergence across species. Glucosinolates have been shown to increase under stress conditions, indicating a potential role in stress response. To elucidate the role of GTRs in glucosinolate transportation under NaCl stress in two distinct Brassica species, B. juncea and B. napus, plants were subjected to 0, 100, or 200 mM NaCl. Based on the literature, key GTR genes were chosen and their expression across various plant parts was assessed. Both species displayed divergent trends in their biochemical profiles as well as glucosinolate contents under elevated salt stress conditions. Statistical modelling identified significant contributors to glucosinolate variations, guiding the development of targeted breeding strategies for low-glucosinolate varieties. Notably, GTR2A2 exhibited pronounced expressions in stems, contributing approximately 52% to glucosinolate content variance, while GTR2B1/C2 displayed significant expression in flowers. Additionally, GTR2A1 and GTR1A2/B1 demonstrated noteworthy expression in roots. This study enhances our understanding of glucosinolate regulation under stress conditions, offering avenues to improve Brassica crop quality and resilience.

Comparative study of the quality indices, antioxidant substances, and mineral elements in different forms of cabbage

BACKGROUND

As the second largest leafy vegetable, cabbage (Brassica oleracea L. var. capitata) is grown globally, and the characteristics of the different varieties, forms, and colors of cabbage may differ. In this study, five analysis methods-variance analysis, correlation analysis, cluster analysis, principal component analysis, and comprehensive ranking-were used to evaluate the quality indices (soluble protein, soluble sugar, and nitrate), antioxidant content (vitamin C, polyphenols, and flavonoids), and mineral (K, Ca, Mg, Cu, Fe, Mn, and Zn) content of 159 varieties of four forms (green spherical, green oblate, purple spherical, and green cow heart) of cabbage.

RESULTS

The results showed that there are significant differences among different forms and varieties of cabbage. Compared to the other three forms, the purple spherical cabbage had the highest flavonoid, K, Mg, Cu, Mn, and Zn content. A scatter plot of the principal component analysis showed that the purple spherical and green cow heart cabbage varieties were distributed to the same quadrant, indicating that their quality indices and mineral contents were highly consistent, while those of the green spherical and oblate varieties were irregularly distributed. Overall, the green spherical cabbage ranked first, followed by the green cow heart, green oblate, and purple spherical varieties.

CONCLUSIONS

Our results provide a theoretical basis for the cultivation and high-quality breeding of cabbage.

Wheat BREVIS RADIX (BRX) regulates organ size, stomatal density and enhances drought tolerance in Arabidopsis

BREVIS RADIX (BRX) is a small plant-specific and evolutionary conserved gene family with divergent yet partially redundant biological functions including root and shoot growth, stomatal development and tiller angle in plants. We characterized a BRX family gene from wheat (Triticum aestivum) by gain-of-function in Arabidopsis. Overexpression of TaBRXL2A resulted in longer primary roots with increased root meristem size and higher root growth under control and exogenous hormone treatments as compared to wild type (Col-0) plants. Overexpression lines also exhibited significant differences with the wild type such as increased rosette size, higher leaf number and leaf size. At reproductive stage, overexpression lines exhibited wider siliques and higher grain weight per plant. Under drought stress, overexpression lines exhibited enhanced drought tolerance in terms of higher chlorophyll retention and lower oxidative stress, thereby leading to significant recovery from drought stress. The analysis suggests that the inherent lower stomatal density in the leaves of overexpression lines and higher stomatal closure in response to ABA might contribute to lower water loss from the overexpression lines. Furthermore, TaBRXL2A protein showed membrane localization, presence of conserved residues at N-terminal for palmitoylation, and phosphosites in the linker region which are prescribed for its potential role in protophloem differentiation and stomatal lineage. Thus, we identified a TaBRX family gene which is involved in developmental pathways essential for plant growth, and also enhances drought tolerance in Arabidopsis.

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.

Cuticular proteins in codling moth (Cydia pomonella) respond to insecticide and temperature stress

The insect cuticle consists of chitin and cuticular proteins (CPs), which stabilize the body shape and provide an effective physical barrier against the external environment. They are also potential target sites for developing environmentally friendly insect management through the utilization of physiology-based methods. The codling moth, Cydia pomonella, is a pest afflicting fruit orchards worldwide. This study used a comparative genomic approach, whole-genome resequencing, and transcriptome data to understand the role that CPs played in the environmental adaptation of the codling moth. A total of 182 putative CPs were identified in the codling moth genome, which were classified into 12 CP families. 119 CPR genes, including 54 RR-1, 60 RR-2, and 5 RR-3 genes were identified and accounted for 65.4% of the total CPs. Eight and seven gene clusters are formed in RR1 and RR2 subfamily and the ancestor-descendant relationship was explained. Five CPAP genes were highly expressed during the egg stage and exposed to high temperature, which indicated their potential role in aiding codling moth eggs in acclimating to varying external heat conditions. Moreover, six CPs belonging to the CPR and CPLCP families were identified in association with insecticide resistance by population resequencing. Their expression levels increased after exposure to insecticides, suggesting they might be involved in codling moth resistance to the insecticides azinphos-methyl or deltamethrin. Our results provide insight into the evolution of codling moth CPs and their association with high temperature adaptation and insecticide resistance, and provide an additional information required for further analysis of CPs in environmental adaptation.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

The common carp (Cyprinus carpio) might best represent the domesticated allopolyploid animals. Although subgenome divergence which is well-known to be a key to allopolyploid domestication has been comprehensively characterized in common carps, the link between genetic architecture underlying agronomic traits and subgenome divergence is unknown in the selective breeding of common carps globally.

RESULTS

We utilized a comprehensive SNP dataset in 13 representative common carp strains worldwide to detect genome-wide genetic variations associated with scale reduction, vibrant skin color, and high growth rate in common carp domestication. We identified numerous novel candidate genes underlie the three agronomically most desirable traits in domesticated common carps, providing potential molecular targets for future genetic improvement in the selective breeding of common carps. We found that independently selective breeding of the same agronomic trait (e.g., fast growing) in common carp domestication could result from completely different genetic variations, indicating the potential advantage of allopolyploid in domestication. We observed that candidate genes associated with scale reduction, vibrant skin color, and/or high growth rate are repeatedly enriched in the immune system, suggesting that domestication of common carps was often accompanied by the disease resistance improvement.

CONCLUSIONS

In common carp domestication, asymmetric subgenome selection is prevalent, while parallel subgenome selection occurs in selective breeding of common carps. This observation is not due to asymmetric gene retention/loss between subgenomes but might be better explained by reduced pleiotropy through transposable element-mediated expression divergence between ohnologs. Our results demonstrate that domestication benefits from polyploidy not only in plants but also in animals.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

Physiological Control and Genetic Basis of Leaf Curvature and Heading in Brassica rapa L

BACKGROUND

Heading is an important agronomic feature for Chinese cabbage, cabbage, and lettuce. The heading leaves function as nutrition storage organs, which contribute to the high quality and economic worth of leafy heads. Leaf development is crucial during the heading stage, most genes previously predicted to be involved in the heading process are based on Arabidopsis leaf development studies.

AIM OF REVIEW

Till date, there is no published review article that demonstrated a complete layout of all the identified regulators of leaf curvature and heading. In this review, we have summarized all the identified physiological and genetic regulators that are directly or indirectly involved in leaf curvature and heading in Brassica crops. By integrating all identified regulators that provide a coherent logic of leaf incurvature and heading, we proposed a molecular mechanism in Brassica crops with graphical illustrations. This review adds value to future breeding of distinct heading kinds of cabbage and Chinese cabbage by providing unique insights into leaf development.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Leaf curvature and heading are established by synergistic interactions among genes, transcription factors, microRNAs, phytohormones, and environmental stimuli that regulate primary and secondary morphogenesis. Various genes have been identified using transformation and genome editing that are responsible for the formation of leaf curvature and heading in Brassica crops. A range of leaf morphologies have been observed in Brassica, which are established because of the mutated determinants that are responsible for cell division and leaf polarity.

Genetic dissection of morphological variation between cauliflower and a rapid cycling Brassica oleracea line

To improve resolution to small genomic regions and sensitivity to small-effect loci in the identification of genetic factors conferring the enlarged inflorescence and other traits of cauliflower while also expediting further genetic dissection, 104 near-isogenic introgression lines (NIILs) covering 78.56% of the cauliflower genome, were selected from an advanced backcross population using cauliflower [Brassica oleracea var. botrytis L., mutant for Orange gene (ORG)] as the donor parent and a rapid cycling line (TO1434) as recurrent parent. Subsets of the advanced backcross population and NIILs were planted in the field for 8 seasons, finding 141 marker-trait associations for 15 leaf-, stem-, and flower-traits. Exemplifying the usefulness of these lines, we delineated the previously known flower color gene to a 4.5 MB interval on C3; a gene for small plant size to a 3.4 MB region on C8; and a gene for large plant size and flowering time to a 6.1 MB region on C9. This approach unmasked closely linked QTL alleles with opposing effects (on chr. 8) and revealed both alleles with expected phenotypic effects and effects opposite the parental phenotypes. Selected B. oleracea NIILs with short generation time add new value to widely used research and teaching materials.

Comparative pangenome analyses provide insights into the evolution of Brassica rapa resistance gene analogues (RGAs)

Brassica rapa is grown worldwide as economically important vegetable and oilseed crop. However, its production is challenged by yield-limiting pathogens. The sustainable control of these pathogens mainly relies on the deployment of genetic resistance primarily driven by resistance gene analogues (RGAs). While several studies have identified RGAs in B. rapa, these were mainly based on a single genome reference and do not represent the full range of RGA diversity in B. rapa. In this study, we utilized the B. rapa pangenome, constructed from 71 lines encompassing 12 morphotypes, to describe a comprehensive repertoire of RGAs in B. rapa. We show that 309 RGAs were affected by presence-absence variation (PAV) and 223 RGAs were missing from the reference genome. The transmembrane leucine-rich repeat (TM-LRR) RGA class had more core gene types than variable genes, while the opposite was observed for nucleotide-binding site leucine-rich repeats (NLRs). Comparative analysis with the B. napus pangenome revealed significant RGA conservation (93%) between the two species. We identified 138 candidate RGAs located within known B. rapa disease resistance QTL, of which the majority were under negative selection. Using blackleg gene homologues, we demonstrated how these genes in B. napus were derived from B. rapa. This further clarifies the genetic relationship of these loci, which may be useful in narrowing-down candidate blackleg resistance genes. This study provides a novel genomic resource towards the identification of candidate genes for breeding disease resistance in B. rapa and its relatives.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

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.

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.

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.

Differences in the transcriptional immune response to Albugo candida between white rust resistant and susceptible cultivars in Brassica rapa L

Albugo candida causing white rust disease decreases the yield of Brassica rapa vegetables greatly. Resistant and susceptible cultivars in B. rapa vegetables have different immune responses against A. candida inoculation, however, the mechanism of how host plants respond to A. candida is still unknown. Using RNA-sequencing, we identified differentially expressed genes (DEGs) between A. candida inoculated [48 and 72 h after inoculation (HAI)] and non-inoculated samples in resistant and susceptible cultivars of komatsuna (B. rapa var. perviridis). Functional DEGs differed between the resistant and susceptible cultivars in A. candida inoculated samples. Salicylic acid (SA) responsive genes tended to be changed in their expression levels by A. candida inoculation in both resistant and susceptible cultivars, but different genes were identified in the two cultivars. SA-dependent systemic acquired resistance (SAR) involving genes were upregulated following A. candida inoculation in the resistant cultivar. Particular genes categorized as SAR that changed expression levels overlapped between A. candida and Fusarium oxysporum f. sp. conglutinans inoculated samples in resistant cultivar, suggesting a role for SAR in defense response to both pathogens particularly in the effector-triggered immunity downstream pathway. These findings will be useful for understanding white rust resistance mechanisms in B. rapa.

Comparative transcriptome analysis provides insights into molecular pathway and genes associated with head-type formation and phenotypic divergence in Chinese cabbage

Background: The heading type of Chinese cabbage is a significant commercial trait with high economic value. At present, research on the phenotypic divergence and formation mechanism of heading type is limited. Results: Through comparative-transcriptome analysis, the formation and phenotypic divergence mechanism of the leafy head of diploid overlapping type cabbage, diploid outward-curling type cabbage, tetraploid overlapping type cabbage, and tetraploid outward-curling type cabbage were systematically and comprehensively investigated, and the phenotype-specific genes of four varieties were revealed. These phenotype-specific differentially expressed genes (DEGs) were considered crucial for cabbage heading type through WGCNA. Some transcription factors have been predicted as significant genes for phenotypic divergence, including the members of the bHLH, AP2/ERF-ERF, WRKY, MYB, NAC, and C2CH2 families. Phytohormone-related genes, including abscisic acid/auxin hormone, may play an important role in the phenotypic divergence of head type in cabbage. Conclusion: Comparative-transcriptome analysis supports a role for phytohormone-related genes and some transcription factors in head-type formation and divergence for four cultivars. These findings increase our understanding of the molecular basis for pattern formation and divergence of the leafy heads of Chinese cabbage and will contribute to developing more desirable leafy head patterns.

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.

Chinese cabbage: an emerging model for functional genomics in leafy vegetable crops

Leafy vegetable crops (LVCs) are consumed worldwide and offer essential nutrients for humans. Unlike model plant species, systematic characterisation of gene function is lacking, although whole-genome sequences (WGSs) are available for various LVCs. Several recent studies in Chinese cabbage have reported high-density mutant populations linking genotype to phenotype, providing blueprints for functional LVC genomics and beyond.

Genome wide analysis of BREVIS RADIX gene family from wheat (Triticum aestivum): A conserved gene family differentially regulated by hormones and abiotic stresses

BREVIS RADIX is a plant specific gene family with unique protein-protein interaction domain. It regulates developmental processes viz. root elongation and tiller angle which are pertinent for crop improvement. In the present study, five BRX family genes were identified in wheat genome and clustered into five sub-groups. Phylogenetic and synteny analyses revealed evolutionary conservation among BRX proteins from monocot species. Expression analyses showed abundance of TaBRXL1 transcripts in vegetative and reproductive tissues except flag leaf. TaBRXL2, TaBRXL3 and TaBRXL4 showed differential, tissue specific and lower level expression as compared to TaBRXL1. TaBRXL5-A expressed exclusively in stamens. TaBRXL1 was upregulated under biotic stresses while TaBRXL2 expression was enhanced under abiotic stresses. TaBRXL2 and TaBRXL3 were upregulated by ABA and IAA in roots. In shoot, TaBRXL2 was upregulated by ABA while TaBRXL3 and TaBRXL4 were upregulated by IAA. Expression levels, tissue specificity and response time under different conditions suggest distinct as well as overlapping functions of TaBRX genes. This was also evident from global co-expression network of these genes. Further, TaBRX proteins exhibited homotypic and heterotypic interactions which corroborated with the role of BRX domain in protein-protein interaction. This study provides leads for functional characterization of TaBRX genes.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Auxin participates in regulating the leaf curl development of Wucai (Brassica campestris L.)

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) belongs to the Brassica genus of the Cruciferae family, and its leaf curl is a typical feature that distinguishes Wucai from other nonheading cabbage subspecies. Our previous research found that plant hormones were involved in the development of the leaf curl in Wucai. However, the molecular mechanisms and the hormones regulating the formation of leaf curl in Wucai have not yet been reported. This study aimed to understand the molecular functions related to hormone metabolism during the formation of leaf curl in Wucai. A total of 386 differentially expressed genes (DEGs) were identified by transcriptome sequencing of two different morphological parts of the same leaf of Wucai germplasm W7-2, and 50 DEGs were found to be related to plant hormones, which were mainly involved in the auxin signal transduction pathway. Then, we measured the content of endogenous hormones in two different forms of the same leaf of Wucai germplasm W7-2. A total of 17 hormones with differential content were identified, including auxin, cytokinins, jasmonic acids, salicylic acids, and abscisic acid. And we found that treatment with auxin transport inhibitor N-1-naphthylphthalamic acid can affect the leaf curl phenotype of Wucai and pak choi (Brassica rapa L. subsp. Chinensis). These results indicated that plant hormones, especially auxin, are involved in developing the leaf curl of Wucai. Our findings provide a potentially valuable reference for future research on the development of leaf curls.

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.

Identification and in vitro enzymatic activity analysis of the AOP2 gene family associated with glucosinolate biosynthesis in Tumorous stem mustard (Brassica juncea var. tumida)

The major enzyme encoded by the glucosinolate biosynthetic gene AOP2 is involved in catalyzing the conversion of glucoiberin (GIB) into sinigrin (SIN) in Brassicaceae crops. The AOP2 proteins have previously been identified in several Brassicaceae species, but not in Tumorous stem mustard. As per this research, the five identified members of the AOP2 family from the whole genome of Brassica juncea named BjuAOP2.1-BjuAOP2.5 were found to be evenly distributed on five chromosomes. The subcellular localization results implied that BjuAOP2 proteins were mainly concentrated in the cytoplasm. Phylogenetic analysis of the AOP2 proteins from the sequenced Brassicaceae species in BRAD showed that BjuAOP2 genes were more closely linked to Brassica carinata and Brassica rapa than Arabidopsis. In comparison with other Brassicaceae plants, the BjuAOP2 members were conserved in terms of gene structures, protein sequences, and motifs. The light response and hormone response elements were included in the BjuAOP2 genes' cis-regulatory elements. The expression pattern of BjuAOP2 genes was influenced by the different stages of development and the type of tissue being examined. The BjuAOP2 proteins were used to perform the heterologous expression experiment. The results showed that all the five BjuAOP2 proteins can catalyze the conversion of GIB to SIN with different catalytic activity. These results provide the basis for further investigation of the functional study of BjuAOP2 in Tumorous stem mustard glucosinolate biosynthesis.

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.

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.

Their C-termini divide Brassica rapa FT-like proteins into FD-interacting and FD-independent proteins that have different effects on the floral transition

Members of the FLOWERING LOCUS T (FT)-like clade of phosphatidylethanolamine-binding proteins (PEBPs) induce flowering by associating with the basic leucine zipper (bZIP) transcription factor FD and forming regulatory complexes in angiosperm species. However, the molecular mechanism of the FT-FD heterocomplex in Chinese cabbage (Brassica rapa ssp. pekinensis) is unknown. In this study, we identified 12 BrPEBP genes and focused our functional analysis on four BrFT-like genes by overexpressing them individually in an FT loss-of-function mutant in Arabidopsis thaliana. We determined that BrFT1 and BrFT2 promote flowering by upregulating the expression of floral meristem identity genes, whereas BrTSF and BrBFT, although close in sequence to their Arabidopsis counterparts, had no clear effect on flowering in either long- or short-day photoperiods. We also simultaneously genetically inactivated BrFT1 and BrFT2 in Chinese cabbage using CRISPR/Cas9-mediated genome editing, which revealed that BrFT1 and BrFT2 may play key roles in inflorescence organogenesis as well as in the transition to flowering. We show that BrFT-like proteins, except for BrTSF, are functionally divided into FD interactors and non-interactors based on the presence of three specific amino acids in their C termini, as evidenced by the observed interconversion when these amino acids are mutated. Overall, this study reveals that although BrFT-like homologs are conserved, they may have evolved to exert functionally diverse functions in flowering via their potential to be associated with FD or independently from FD in Brassica rapa.

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.

Comprehensive Analyses of the Histone Deacetylases Tuin (HDT) Gene Family in Brassicaceae Reveals Their Roles in Stress Response

Histone deacetylases tuin (HDT) is a plant-specific protein subfamily of histone deacetylation enzymes (HDAC) which has a variety of functions in plant development, hormone signaling and stress response. Although the HDT family's genes have been studied in many plant species, they have not been characterized in Brassicaceae. In this study, 14, 8 and 10 HDT genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the HDTs were divided into four groups: HDT1(HD2A), HDT2(HD2B), HDT3(HD2C) and HDT4(HD2D). There was an expansion of HDT2 orthologous genes in Brassicaceae. Most of the HDT genes were intron-rich and conserved in gene structure, and they coded for proteins with a nucleoplasmin-like (NPL) domain. Expression analysis showed that B. napus, B. rapa, and B. oleracea HDT genes were expressed in different organs at different developmental stages, while different HDT subgroups were specifically expressed in specific organs and tissues. Interestingly, most of the Bna/Br/BoHDT2 members were expressed in flowers, buds and siliques, suggesting they have an important role in the development of reproductive organs in Brassicaceae. Expression of BnaHDT was induced by various hormones, such as ABA and ethylene treatment, and some subgroups of genes were responsive to heat treatment. The expression of most HDT members was strongly induced by cold stress and freezing stress after non-cold acclimation, while it was slightly induced after cold acclimation. In this study, the HDT gene family of Brassicaceae was analyzed for the first time, which helps in understanding the function of BnaHDT in regulating plant responses to abiotic stresses, especially freezing stresses.

GADMA2: more efficient and flexible demographic inference from genetic data

BACKGROUND

Inference of complex demographic histories is a source of information about events that happened in the past of studied populations. Existing methods for demographic inference typically require input from the researcher in the form of a parameterized model. With an increased variety of methods and tools, each with its own interface, the model specification becomes tedious and error-prone. Moreover, optimization algorithms used to find model parameters sometimes turn out to be inefficient, for instance, by being not properly tuned or highly dependent on a user-provided initialization. The open-source software GADMA addresses these problems, providing automatic demographic inference. It proposes a common interface for several likelihood engines and provides global parameters optimization based on a genetic algorithm.

RESULTS

Here, we introduce the new GADMA2 software and provide a detailed description of the added and expanded features. It has a renovated core code base, new likelihood engines, an updated optimization algorithm, and a flexible setup for automatic model construction. We provide a full overview of GADMA2 enhancements, compare the performance of supported likelihood engines on simulated data, and demonstrate an example of GADMA2 usage on 2 empirical datasets.

CONCLUSIONS

We demonstrate the better performance of a genetic algorithm in GADMA2 by comparing it to the initial version and other existing optimization approaches. Our experiments on simulated data indicate that GADMA2's likelihood engines are able to provide accurate estimations of demographic parameters even for misspecified models. We improve model parameters for 2 empirical datasets of inbred species.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce

The mechanisms underlying leafy heads in vegetables are poorly understood. Here, we cloned a quantitative trait locus (QTL) controlling leafy heads in lettuce (Lactuca sativa). The QTL encodes a transcription factor, SAWTOOTH 1 (LsSAW1), which has a BEL1-like homeodomain and is a homolog of Arabidopsis thaliana. A 1-bp deletion in Lssaw1 contributes to the development of leafy heads. Laser-capture microdissection and RNA-sequencing showed that LsSAW1 regulates leaf dorsiventrality and loss-of-function of Lssaw1 downregulates the expression of many adaxial genes but upregulates abaxial genes. LsSAW1 binds to the promoter region of the adaxial gene ASYMMETRIC LEAVES 1 (LsAS1) to upregulate its expression. Overexpression of LsAS1 compromised the effects of Lssaw1 on heading. LsSAW1 also binds to the promoter region of the abaxial gene YABBY 1 (LsYAB1), but downregulates its expression. Overexpression of LsYAB1 led to bending leaves in LsSAW1 genotypes. LsSAW1 directly interacts with KNOTTED 1 (LsKN1), which is necessary for leafy heads in lettuce. RNA-seq data showed that LsSAW1 and LsKN1 exert antagonistic effects on the expression of thousands of genes. LsSAW1 compromises the ability of LsKN1 to repress LsAS1. Our results suggest that downregulation or loss-of-function of adaxial genes and upregulation of abaxial genes allow for the development of leafy heads.

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.

An identical-by-descent segment harbors a 12-bp insertion determining fruit softening during domestication and speciation in Pyrus

BACKGROUND

Although the wild relatives of pear originated in southwest China, this fruit crop was independently domesticated and improved in Asia and Europe, and there are major phenotypic differences (e.g., maturity and fruit firmness) between Asian and European pears.  RESULTS: In this study, we examined the genomes of 113 diverse pear accessions using an identity-by-descent (IBD) approach to investigate how historical gene flow has shaped fruit firmness traits in Asian and European pears. We found a 3-Mbp IBD-enriched region (IBD-ER) that has undergone "convergent domestication" in both the Asian and European pear lineages, and a genome-wide association study (GWAS) of fruit firmness phenotypes strongly implicated the TRANSLOCON AT THE INNER CHLOROPLAST ENVELOPE55 (TIC55) locus within this 3-Mbp IBD-ER. Furthermore, we identified a tandem duplication that includes a 12-bp insertion located in the first exon of TIC55 that is uniquely present in Asian pears, and expression analysis showed that the pear TIC55 gene is highly expressed in Asian pear, while it is weakly or not expressed in European pear; this could contribute to the differences in fruit firmness between Asian and European pear fruits.

CONCLUSIONS

Our findings provide insights into how pear fruit softening has been impacted during domestication, and we identified candidate genes associated with fruit softening that can contribute to the breeding and improvement of pear and other fruit crops.

Genome-wide analysis of the JAZ subfamily of transcription factors and functional verification of BnC08.JAZ1-1 in Brassica napus

BACKGROUND

JAZ subfamily plays crucial roles in growth and development, stress, and hormone responses in various plant species. Despite its importance, the structural and functional analyses of the JAZ subfamily in Brassica napus are still limited.

RESULTS

Comparing to the existence of 12 JAZ genes (AtJAZ1-AtJAZ12) in Arabidopsis, there are 28, 31, and 56 JAZ orthologues in the reference genome of B. rapa, B. oleracea, and B. napus, respectively, in accordance with the proven triplication events during the evolution of Brassicaceae. The phylogenetic analysis showed that 127 JAZ proteins from A. thaliana, B. rapa, B. oleracea, and B. napus could fall into five groups. The structure analysis of all 127 JAZs showed that these proteins have the common motifs of TIFY and Jas, indicating their conservation in Brassicaceae species. In addition, the cis-element analysis showed that the main motif types are related to phytohormones, biotic and abiotic stresses. The qRT-PCR of the representative 11 JAZ genes in B. napus demonstrated that different groups of BnJAZ individuals have distinct patterns of expression under normal conditions or treatments with distinctive abiotic stresses and phytohormones. Especially, the expression of BnJAZ52 (BnC08.JAZ1-1) was significantly repressed by abscisic acid (ABA), gibberellin (GA), indoleacetic acid (IAA), polyethylene glycol (PEG), and NaCl treatments, while induced by methyl jasmonate (MeJA), cold and waterlogging. Expression pattern analysis showed that BnC08.JAZ1-1 was mainly expressed in the vascular bundle and young flower including petal, pistil, stamen, and developing ovule, but not in the stem, leaf, and mature silique and seed. Subcellular localization showed that the protein was localized in the nucleus, in line with its orthologues in Arabidopsis. Overexpression of BnC08.JAZ1-1 in Arabidopsis resulted in enhanced seed weight, likely through regulating the expression of the downstream response genes involved in the ubiquitin-proteasome pathway and phospholipid metabolism pathway.

CONCLUSIONS

The systematic identification, phylogenetic, syntenic, and expression analyses of BnJAZs subfamily improve our understanding of their roles in responses to stress and phytohormone in B. napus. In addition, the preliminary functional validation of BnC08.JAZ1-1 in Arabidopsis demonstrated that this subfamily might also play a role in regulating seed weight.

Mapping and identification of a new potential dominant resistance gene to turnip mosaic virus in Brassica rapa

By constructing an F₂ population, a new potential dominant resistance gene to TuMV in Brassica rapa was mapped and identified. Brassica rapa is the most widely grown vegetable crop in China, and turnip mosaic virus (TuMV) is a great threat to its production. Hence, it is a very important work to excavate more and novel resistance genes in B. rapa. In this study, the resistant line B80124 and the susceptible line B80450 were used to construct the F₂ populations, and through genetic analysis, the resistance to TuMV was found to be controlled by a dominant gene. Bulked segregant analysis sequence (BSA-seq) was used for the primary mapping, and an intersection (22.25-25.03 Mb) was obtained. After fine mapping using single nucleotide polymorphisms (SNP) markers, the candidate region was narrowed to 330 kb between the SNP markers A06S11 and A06S14, including eight genes relating to disease resistance. Using the transcriptome analysis and sequence identification, BraA06g035130.3C was screened as the final candidate gene, and it contained two deletion mutations, leading to frameshift in the susceptible line B80450. In addition, the phylogenetic analysis, hydrophilia and hydrophobicity analysis, subcellular location prediction analysis, amino acid bias analysis, and 3D modeling structures of BraA06g035130.3C were conducted to predict its functions. This study was conducive to the identification of a new TuMV resistance gene in B. rapa, which is of important scientific significance and application value for the improvement of TuMV resistance traits and molecular design breeding for Brassica crops.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Stomatal morphological variation contributes to global ecological adaptation and diversification of Brassica napus

Stomatal density and guard cell length of 274 global core germplasms of rapeseed reveal that the stomatal morphological variation contributes to global ecological adaptation and diversification of Brassica napus. Stomata are microscopic structures of plants for the regulation of CO₂ assimilation and transpiration. Stomatal morphology has changed substantially in the adaptation to the external environment during land plant evolution. Brassica napus is a major crop to produce oil, livestock feed and biofuel in the world. However, there are few studies on the regulatory genes controlling stomatal development and their interaction with environmental factors as well as the genetic mechanism of adaptive variation in B. napus. Here, we characterized stomatal density (SD) and guard cell length (GL) of 274 global core germplasms at seedling stage. It was found that among the significant phenotypic variation, European germplasms are mostly winter rapeseed with high stomatal density and small guard cell length. However, the germplasms from Asia (especially China) are semi-winter rapeseed, which is characterized by low stomatal density and large guard cell length. Through selective sweep analysis and homology comparison, we identified several candidate genes related to stomatal density and guard cell length, including Epidermal Patterning Factor2 (EPF2; BnaA09g23140D), Epidermal Patterning Factor Like4 (EPFL4; BnaC01g22890D) and Suppressor of LLP1 (SOL1 BnaC01g22810D). Haplotype and phylogenetic analysis showed that natural variation in EPF2, EPFL4 and SOL1 is closely associated with the winter, spring, and semi-winter rapeseed ecotypes. In summary, this study demonstrated for the first time the relation between stomatal phenotypic variation and ecological adaptation in rapeseed, which is useful for future molecular breeding of rapeseed in the context of evolution and domestication of key stomatal traits and global climate change.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.