Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U’s triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U’s triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U’s triangle.

Background: Regulatory network of cytoplasmic male sterility (CMS) occurrence is still largely unknown in plants, although numerous researches have been attempted to isolate genes involved in CMS. Here, we employed high-throughput sequencing and degradome analysis to identify microRNAs and their targets using high-throughput sequencing in CMS and its maintainer fertile (MF) lines of Brassica juncea. Results: We identified 197 known and 78 new candidate microRNAs during reproductive development of B. juncea. A total of 47 differentially expressed microRNAs between CMS and its MF lines were discovered, according to their sequencing reads number. Different expression levels of selected microRNAs were confirmed by using real-time quantitative PCR between CMS and MF lines. Furthermore, we observed that the transcriptional patterns of these microRNAs could be mimicked by artificially inhibiting mitochondrial F₁F₀-ATPase activity and its function in MF line by using treatment with oligomycin. Targeted genes of the microRNAs were identified by high-throughput sequencing and degradome approaches, including auxin response factor, NAC (No Apical Meristem) domain transcription factor, GRAS family transcription factor, MYB transcription factor, squamosa promoter binding protein, AP2-type transcription factor, homeobox/homeobox-leucine zipper family and TCP family transcription factors, which were observed to be differentially expressed between CMS and MF. Conclusion: Taken together, from these findings we suggested microRNA might participate in the regulatory network of CMS by tuning fork in gene expressions in CMS B. juncea. The differential expression of miRNAs observed between CMS and MF lines suggested that biogenesis of miRNAs could be influenced in the CMS.

The elucidation of heat tolerance mechanisms is required to combat the challenges of global warming. This study aimed to determine the antioxidant enzyme responses to heat stress, at the enzymatic activity and gene expression levels, and to investigate the antioxidative alterations associated with heat tolerance in the stems and roots of squashes using three genotypes differing in heat tolerance. Plants of heat-tolerant “C. moschata”, thermolabile “C. maxima” and moderately heat-tolerant interspecific inbred line “Maxchata” genotypes were exposed to moderate (37 °C) and severe (42 °C) heat shocks. “C. moschata” exhibited comparatively little oxidative damage, with the lowest hydrogen peroxide (H₂O₂), superoxide (O₂⁻) and malondialdehyde (MDA) contents in the roots compared to stems, followed by “Maxchata”. The enzyme activities of superoxide dismutase (SOD), ascorbate peroxidase (APX), catalase (CAT) and peroxidase (POD) were found to be increased with heat stress in tolerant genotypes. The significant inductions of FeSOD, MnSOD, APX2, CAT1 and CAT3 isoforms in tolerant genotypes suggested their participation in heat tolerance. The differential isoform patterns of SOD, APX and CAT between stems and roots also indicated their tissue specificity. Furthermore, despite the sequence similarity of the studied antioxidant genes among “C. maxima” and “Maxchata”, most of these genes were highly induced under heat stress in “Maxchata”, which contributed to its heat tolerance. This phenomenon also indicated the involvement of other unknown genetic and/or epigenetic factors in controlling the expression of these antioxidant genes in squashes, which demands further exploration.

Salicylic acid (SA) is a ubiquitous plant hormone with a long history in human civilization^1,2. Because of the central role of SA in orchestrating plant pathogen defence, understanding SA biosynthesis is fundamental to plant immunity research and crop improvement. Isochorismate-derived SA biosynthesis has been well defined in Arabidopsis. However, increasing evidence suggests a crucial function for phenylalanine-derived SA biosynthesis in many other plant species¹. Here we reveal the phenylalanine-derived SA biosynthetic pathway in rice by identifying three dedicated enzymes - peroxisomal benzoyl-CoA:benzyl alcohol benzoyltransferase (BEBT), the endoplasmic reticulum-associated cytochrome P450 enzyme benzylbenzoate hydroxylase (BBH), and cytosolic benzylsalicylate esterase (BSE) that sequentially convert benzoyl-CoA to benzylbenzoate, benzylsalicylate and SA. The pathogen-induced gene expression pattern and SA biosynthetic functions of this triple-enzyme module are conserved in diverse plants. This work fills a major knowledge gap in the biosynthesis of a key plant defence hormone, establishing a foundation for new strategies to create disease-resistant crops.

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.