Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid speciesBrassica juncea(L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seededB. junceagenome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported thatB. junceais most likely a single origin in West Asia, 8,000–14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding ofB. juncea.

Brassica juncea, a worldwide cultivated crop plant, produces seeds of different colors. Seed pigmentation is due to the deposition in endothelial cells of proanthocyanidins (PAs), end products from a branch of flavonoid biosynthetic pathway. To elucidate the gene regulatory network of seed pigmentation in B. juncea, transcriptomes in seed coat of a yellow-seeded inbred line and its brown-seeded near- isogenic line were sequenced using the next-generation sequencing platform Illumina/Solexa and de novo assembled. Over 116 million high-quality reads were assembled into 69,605 unigenes, of which about 71.5% (49,758 unigenes) were aligned to Nr protein database with a cut-off E-value of 10⁻⁵. RPKM analysis showed that the brown-seeded testa up-regulated 802 unigenes and down-regulated 502 unigenes as compared to the yellow-seeded one. Biological pathway analysis revealed the involvement of forty six unigenes in flavonoid biosynthesis. The unigenes encoding dihydroflavonol reductase (DFR), leucoantho-cyanidin dioxygenase (LDOX) and anthocyanidin reductase (ANR) for late flavonoid biosynthesis were not expressed at all or at a very low level in the yellow-seeded testa, which implied that these genes for PAs biosynthesis be associated with seed color of B. juncea, as confirmed by qRT-PCR analysis of these genes. To our knowledge, it is the first time to sequence the transcriptome of seed coat in Brassica juncea. The unigene sequences obtained in this study will not only lay the foundations for insight into the molecular mechanisms underlying seed pigmentation in B.juncea, but also provide the basis for further genomics research on this species or its allies.

Plant growth and development require water, but excessive water hinders growth. Sesame (Sesamum indicum L.) is an important oil crop; it is drought-tolerant but sensitive to waterlogging, and its drought tolerance has been extensively studied. However, the waterlogging tolerance of sesame still has relatively few studies. In this study, two kinds of sesame, R (waterlogging-tolerant) and S (waterlogging-intolerant), were used as materials, and they were treated with waterlogging stress for 0, 24, 72, and 120 h. Physiological analysis showed that after waterlogging, sesame plants responded to stress by increasing the contents of ascorbate peroxidase (APX), glutathione (GSH), and some other antioxidants. The results of the multi-omics analysis of sesame under waterlogging stress revealed 15,652 (R) and 12,156 (S) differentially expressed genes (DEGs), 41 (R) and 47 (S) differentially expressed miRNAs (DEMis), and 896 (R) and 1036 (S) differentially accumulated metabolites (DAMs). The combined DEMi-DEG analysis that 24 DEMis regulated 114 DEGs in response to waterlogging stress. In addition, 13 hub genes and three key pathways of plant hormone signal transduction, glutathione metabolism, and glyoxylate and dicarboxylate metabolism were identified by multi-omics analysis under waterlogging stress. The results showed that sesame regulated the content of hormones and antioxidants and promoted energy conversion in the plant through the above pathways to adapt to waterlogging stress. In summary, this study further analyzed the response mechanism of sesame to waterlogging stress and provides helpful information for the breeding of plants for waterlogging tolerance and genetic improvement.

There are three basic Brassica genomes (A, B, and C) and three parallel sets of subgenomes distinguished in the diploid Brassica (i.e.: B. rapa, ArAr; B. nigra, BniBni; B. oleracea, CoCo) and the derived allotetraploid species (i.e.: B. juncea, AjAjBjBj; B. napus, AnAnCnCn; B. carinata, BcBcCcCc). To understand subgenome differentiation in B. juncea in comparison to other A genome-carrying Brassica species (B. rapa and B. napus), we constructed a dense genetic linkage map of B. juncea, and conducted population genetic analysis on diverse lines of the three A-genome carrying Brassica species using a genotyping-by-sequencing approach (DArT-seq). A dense genetic linkage map of B. juncea was constructed using an F2 population derived from Sichuan Yellow/Purple Mustard. The map included 3329 DArT-seq markers on 18 linkage groups and covered 1579 cM with an average density of two markers per cM. Based on this map and the alignment of the marker sequences with the physical genome of Arabidopsis thaliana, we observed strong co-linearity of the ancestral blocks among the different A subgenomes but also considerable block variation. Comparative analyses at the level of genome sequences of B. rapa and B. napus, and marker sequence anchored on the genetic map of B. juncea, revealed a total of 30 potential inversion events across large segments and 20 potential translocation events among the three A subgenomes. Population genetic analysis on 26 accessions of the three A genome-carrying Brassica species showed that the highest genetic distance were estimated when comparing Aj-An than between An-Ar and Aj-Ar subgenome pairs. The development of the dense genetic linkage map of B. juncea with informative DArT-seq marker sequences and availability of the reference sequences of the Ar, and AnCn genomes allowed us to compare the A subgenome structure of B. juncea (Aj) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (Ar, An and Aj) but with apparent subgenomic variation. Population genetic analysis on three A-genome carrying Brassica species support the idea that B. juncea has distinct genomic diversity, and/or evolved from a different A genome progenitor of B. napus.

Strong artificial and natural selection causes the formation of highly conserved haplotypes that harbor agronomically important genes. GWAS combination with haplotype analysis has evolved as an effective method to dissect the genetic architecture of complex traits in crop species. We used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with oleic acid (C18:1) in rapeseed. Six haplotype regions were identified as significantly associated with oleic acid (C18:1) that mapped to chromosomes A02, A07, A08, C01, C02, and C03. Additionally, whole-genome sequencing of 50 rapeseed accessions revealed three genes (BnmtACP2-A02, BnABCI13-A02 and BnECI1-A02) in the A02 chromosome haplotype region and two genes (BnFAD8-C02 and BnSDP1-C02) in the C02 chromosome haplotype region that were closely linked to oleic acid content phenotypic variation. Moreover, the co-expression network analysis uncovered candidate genes from these two different haplotype regions with potential regulatory interrelationships with oleic acid content accumulation. Our results suggest that several candidate genes are closely linked, which provides us with an opportunity to develop functional haplotype markers for the improvement of the oleic acid content in rapeseed.