Comparative mapping between Arabidopsis thaliana and Brassica nigra indicates that Brassica genomes have evolved through extensive genome replication accompanied by chromosome fusions and frequent rearrangements

Genome-Wide Identification of FCS-Like Zinc Finger (FLZ) Family Genes in Three Brassica Plant Species and Functional Characterization of BolFLZs in Chinese Kale Under Abiotic Stresses

FCS-like zinc finger (FLZ) proteins are plant-specific regulatory proteins, which contain a highly conserved FLZ domain, and they play critical roles in plant growth and stress responses. Although the FLZ family has been systematically characterized in certain plants, it remains underexplored in Brassica species, which are vital sources of vegetables, edible oils, and condiments for human consumption and are highly sensitive to various abiotic stresses. Following the whole-genome triplication events (WGT) in Brassica, elucidating how the FLZ genes have expanded, differentiated, and responded to abiotic stresses is valuable for uncovering the genetic basis and functionality of these genes. In this study, we identified a total of 113 FLZ genes from three diploid Brassica species and classified them into four groups on the basis of their amino acid sequences. Additionally, we identified 109 collinear gene pairs across these Brassica species, which are dispersed among different chromosomes, suggesting that whole-genome duplication (WGD) has significantly contributed to the expansion of the FLZ family. Subcellular localization revealed that six representative BolFLZ proteins are located in the nucleus and cytoplasm. Yeast two-hybrid assays revealed that 13 selected BolFLZs interact with BolSnRK1α1 and BolSnRK1α2, confirming the conservation of the SnRK1α-FLZ module in Brassica species. Expression profile analysis revealed differential expression patterns of BolFLZ across various tissues. Notably, the expression levels of seven BolFLZ genes out of the fifteen genes analyzed changed significantly following treatment with various abiotic stressors, indicating that the BolFLZ genes play distinct physiological roles and respond uniquely to abiotic stresses in Brassica species. Together, our results provide a comprehensive overview of the FLZ gene family in Brassica species and insights into their potential applications for enhancing stress tolerance and growth in Chinese kale.

Functional characterisation of BnaA02.TOP1α and BnaC02.TOP1α involved in true leaf biomass accumulation in Brassica napus L

Leaves, as primary photosynthetic organs essential for high crop yield and quality, have attracted significant attention. The functions of DNA topoisomerase 1α (TOP1α) in various biological processes, including leaf development, in Brassica napus remain unknown. Here, four paralogs of BnaTOP1α, namely BnaA01.TOP1α, BnaA02.TOP1α, BnaC01.TOP1α and BnaC02.TOP1α, were identified and cloned in the B. napus inbred line 'K407'. Expression pattern analysis revealed that BnaA02.TOP1α and BnaC02.TOP1α, but not BnaA01.TOP1α and BnaC01.TOP1α, were persistently and highly expressed in B. napus true leaves. Preliminary analysis in Arabidopsis thaliana revealed that BnaA02.TOP1α and BnaC02.TOP1α paralogs, but not BnaA01.TOP1α and BnaC01.TOP1α, performed biological functions. Targeted mutations of four BnaTOP1α paralogs in B. napus using the CRISPR-Cas9 system revealed that BnaA02.TOP1α and BnaC02.TOP1α served as functional paralogs and redundantly promoted true leaf number and size, thereby promoting true leaf biomass accumulation. Moreover, BnaA02.TOP1α modulated the levels of endogenous gibberellins, cytokinins and auxins by indirectly regulating several genes related to their metabolism processes. BnaA02.TOP1α directly activated BnaA03.CCS52A2 and BnaC09.AN3 by facilitating the recruitment of RNA polymerase II and modulating H3K27me3, H3K36me2 and H3K36me3 levels at these loci and indirectly activated the BnaA08.PARL1 expression, thereby positively controlling the true leaf size in B. napus. Additionally, BnaA02.TOP1α indirectly activated the BnaA07.PIN1 expression to positively regulate the true leaf number. These results reveal the important functions of BnaTOP1α and provide insights into the regulatory network controlling true leaf biomass accumulation in B. napus.

BnaC06.WIP2-BnaA09.STM transcriptional regulatory module promotes leaf lobe formation in Brassica napus

The lobed leaves of rapeseed (Brassica napus L.) offer significant advantages in dense planting, leading to increased yield. Although AtWIP2, a C2H2 zinc finger transcription factor, acts as a regulator of leaf development in Arabidopsis thaliana, the function and regulatory mechanisms of BnaWIP2 in B. napus remain unclear. Here, constitutive expression of the BnaC06.WIP2 paralog, predominantly expressed in leaf serrations, produced lobed leaves in both A. thaliana and B. napus. We demonstrated that BnaC06.WIP2 directly repressed the expression of BnaA01.TCP4, BnaA03.TCP4, and BnaC03.TCP4 and indirectly inhibited the expression of BnaA05.BOP1 and BnaC02.AS2 to promote leaf lobe formation. On the other hand, we discovered that BnaC06.WIP2 modulated the levels of endogenous gibberellin, cytokinin, and auxin, and controlled the auxin distribution in B. napus leaves, thus accelerating leaf lobe formation. Meanwhile, we revealed that BnaA09.STM physically interacted with BnaC06.WIP2, and ectopic expression of BnaA09.STM generated smaller and lobed leaves in B. napus. Furthermore, we found that BnaC06.WIP2 and BnaA09.STM synergistically promoted leaf lobe formation through forming transcriptional regulatory module. Collectively, our findings not only facilitate in-depth understanding of the regulatory mechanisms underlying lobed leaf formation, but also are helpful for guiding high-density breeding practices through improving leaf morphology in B. napus.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids

KEY MESSAGE

We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

Comprehensive Evolution and Expression anaLysis of PHOSPHATE 1 Gene Family in Allotetraploid Brassica napus and Its Diploid Ancestors

The members of PHOSPHATE 1 (PHO1) family play important roles in plant phosphate (Pi) transport and adaptation to Pi deficiency. The functions of PHO1 family proteins have been reported in several plant species, with the exception of Brassica species. Here, we identified 23, 23, and 44 putative PHO1 family genes in Brassica rapa, Brassica oleracea, and Brassica napus by whole genome analysis, respectively. The phylogenetic analysis divided PHO1 family proteins into eight groups, which represented the orthologous relationships among PHO1 members. The gene structure and the conserved motif analysis indicated that the most PHO1 family genes had similar gene structures and the PHO1 proteins shared mutual conserved motifs. The chromosome distribution analysis showed that the majority of BnPHO1 family genes distributed analogously at chromosomes with BrPHO1 and BoPHO1 family genes. The data showed that PHO1 family genes were highly conserved during evolution from diploid to tetraploid. Furthermore, the expression analysis showed that PHO1 family genes had different expression patterns in plant tissues, suggesting the diversity of gene functions in Brassica species. Meanwhile, the expression analysis also revealed that some PHO1 family genes were significantly responsive to Pi deficiency, suggesting that PHO1 family genes play critical roles in Pi uptake and homeostasis under low Pi stress. Altogether, the characteristics of PHO1 family genes provide a reliable groundwork for further dissecting their functions in Brassica species.

H2S-mediated balance regulation of stomatal and non-stomatal factors responding to drought stress in Chinese cabbage

Increased evidence has shown that hydrogen sulfide (H₂S), a novel gasotransmitter, could enhance drought resistance in plants by inducing stomatal closure, with concurrent enhancement of photosynthetic efficiency, but little is known about the mechanism behind this contradictory phenomenon. This study examined the regulating mechanism of H₂S in response to drought stress from stomatal and non-stomatal factors in Chinese cabbage. The results showed that exogenous H₂S could increase the accumulation of photosynthetic pigments and alleviate the damage caused by drought stress. It also regulated the expression in transcriptional level and the activity of ribulose 1,5-bisphosphate carboxylase/oxygenase (BrRuBisCO) under drought stress. The large subunit of BrRuBisCO was found to be modified by S-sulfhydration, which might be the reason for its increased enzyme activity. The fluxes of Cl^-, K⁺, and H⁺ in the guard cells were detected by non-invasive micro-test techniques while under drought stress. The results indicated that H₂S signaling induced a transmembrane Cl^- and H⁺ efflux and inhibited K⁺ influx, and the Cl^- channel was the main responders for H₂S-regulated stomatal movement. In conclusion, H₂S signal not only activated the ion channel proteins located in the guard cell membrane to induce stomatal closure, but also regulated the transcriptional expression and the activity of RuBisCO, a non-stomatal factor to enhance the photosynthetic efficiency of leaves. There is therefore a beneficial balance between the regulation of H₂S signaling on stomatal factors and non-stomatal factors due to drought stress, which needs to be better understood to apply it practically to increase crop yields.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Overexpression of BnaAGL11, a MADS-Box Transcription Factor, Regulates Leaf Morphogenesis and Senescence in Brassica napus

Previous studies have reported that SEEDSTICK/AGAMOUS-LIKE 11 (AtSTK/AtAGL11), a MADS-box transcription factor, plays important roles in many biological processes in Arabidopsis thaliana. However, the function of BnaAGL11, an AtAGL11 homologous gene from Brassica napus, in leaf development remains unknown. Here, we found that the ectopic expression of any copy of Bna.C09.AGL11, Bna.A03.AGL11, and Bna.A09.AGL11 in A. thaliana led to smaller and curly leaves and promoted leaf senescence. Consistently, the overexpression of Bna.C09.AGL11 in B. napus also caused smaller and curly leaves and accelerated leaf senescence. Furthermore, we demonstrated that Bna.C09.AGL11 controlled leaf morphogenesis by indirectly downregulating the genes of Bna.A01.DWF4 and Bna.C07.PGX3 and promoted leaf senescence by indirectly upregulating the genes of Bna.A04.ABI5, Bna.A05.ABI5, Bna.C04.ABI5-1, and Bna.C01.SEN4 and directly activating the transcription of Bna.C04.ABI5-2 and Bna.C03.SEN4 genes. Our results provide new insights into the underlying regulatory mechanism of BnaAGL11 during leaf development in B. napus.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

Molecular Breeding Strategy and Challenges Towards Improvement of Downy Mildew Resistance in Cauliflower (Brassica oleracea var. botrytis L.)

Cauliflower (Brassica oleracea var. botrytis L.) is one of the important, nutritious and healthy vegetable crops grown and consumed worldwide. But its production is constrained by several destructive fungal diseases and most importantly, downy mildew leading to severe yield and quality losses. For sustainable cauliflower production, developing resistant varieties/hybrids with durable resistance against broad-spectrum of pathogens is the best strategy for a long term and reliable solution. Identification of novel resistant resources, knowledge of the genetics of resistance, mapping and cloning of resistance QTLs and identification of candidate genes would facilitate molecular breeding for disease resistance in cauliflower. Advent of next-generation sequencing technologies (NGS) and publishing of draft genome sequence of cauliflower has opened the flood gate for new possibilities to develop enormous amount of genomic resources leading to mapping and cloning of resistance QTLs. In cauliflower, several molecular breeding approaches such as QTL mapping, marker-assisted backcrossing, gene pyramiding have been carried out to develop new resistant cultivars. Marker-assisted selection (MAS) would be beneficial in improving the precision in the selection of improved cultivars against multiple pathogens. This comprehensive review emphasizes the fascinating recent advances made in the application of molecular breeding approach for resistance against an important pathogen; Downy Mildew (Hyaloperonospora parasitica) affecting cauliflower and Brassica oleracea crops and highlights the QTLs identified imparting resistance against this pathogen. We have also emphasized the critical research areas as future perspectives to bridge the gap between availability of genomic resources and its utility in identifying resistance genes/QTLs to breed downy mildew resistant cultivars. Additionally, we have also discussed the challenges and the way forward to realize the full potential of molecular breeding for downy mildew resistance by integrating marker technology with conventional breeding in the post-genomics era. All this information will undoubtedly provide new insights to the researchers in formulating future breeding strategies in cauliflower to develop durable resistant cultivars against the major pathogens in general and downy mildew in particular.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Expression of a Brassica napus metal transport protein (BnMTP3) in Arabidopsis thaliana confers tolerance to Zn and Mn

The essential micronutrient elements zinc (Zn) and manganese (Mn) are crucial for plant growth and development. As an important oil crop, the yield and quality of rapeseed are affected by Zn and Mn toxicity. The cation diffusion facilitator (CDF) family of proteins play significant roles in maintaining intracellular ionic homeostasis and tolerance in plants. However, research on CDF proteins in rapeseed is lacking. In this study, the function of a Brassica napus cation diffusion facilitator/ metal tolerance protein (CDF/MTP) was investigated. The protein, abbreviated BnMTP3 is homologous to the Arabidopsis thaliana MTP3 (AtMTP3). Heterologous expression of BnMTP3 in yeast enhanced tolerance and intracellular sequestration of Zn and Mn. Expression of BnMTP3 in A. thaliana increased Zn and Mn tolerance and markedly increased Zn accumulation in roots. Quantitative RT-PCR analysis showed that BnMTP3 is primarily expressed in roots, and subcellular localization suggested that BnMTP3 is localized in the trans-Golgi network (TGN) and the prevacuolar compartment (PVC) in Arabidopsis and rape protoplast. After treatment with Zn and Mn, BnMTP3 was observed on the vacuolar membrane in transgenic Arabidopsis lines. These findings suggest that BnMTP3 confers Zn and Mn tolerance by sequestering Zn and/or Mn into the vacuole.

Genome wide association analyses to understand genetic basis of flowering and plant height under three levels of nitrogen application in Brassica juncea (L.) Czern & Coss

Timely transition to flowering, maturity and plant height are important for agronomic adaptation and productivity of Indian mustard (B. juncea), which is a major edible oilseed crop of low input ecologies in Indian subcontinent. Breeding manipulation for these traits is difficult because of the involvement of multiple interacting genetic and environmental factors. Here, we report a genetic analysis of these traits using a population comprising 92 diverse genotypes of mustard. These genotypes were evaluated under deficient (N75), normal (N100) or excess (N125) conditions of nitrogen (N) application. Lower N availability induced early flowering and maturity in most genotypes, while high N conditions delayed both. A genotyping-by-sequencing approach helped to identify 406,888 SNP markers and undertake genome wide association studies (GWAS). 282 significant marker-trait associations (MTA's) were identified. We detected strong interactions between GWAS loci and nitrogen levels. Though some trait associated SNPs were detected repeatedly across fertility gradients, majority were identified under deficient or normal levels of N applications. Annotation of the genomic region (s) within ± 50 kb of the peak SNPs facilitated prediction of 30 candidate genes belonging to light perception, circadian, floral meristem identity, flowering regulation, gibberellic acid pathways and plant development. These included over one copy each of AGL24, AP1, FVE, FRI, GID1A and GNC. FLC and CO were predicted on chromosomes A02 and B08 respectively. CDF1, CO, FLC, AGL24, GNC and FAF2 appeared to influence the variation for plant height. Our findings may help in improving phenotypic plasticity of mustard across fertility gradients through marker-assisted breeding strategies.

Hydrogen sulfide promotes flowering in heading Chinese cabbage by S-sulfhydration of BraFLCs

Heading Chinese cabbage (Brassica rapa L. syn. B. campestris L. ssp. chinensis Makino var. pekinensis (Rupr.) J. Cao et Sh. Cao) is a cruciferous Brassica vegetable that has a triplicate genome, owing to an ancient genome duplication event. It is unclear whether the duplicated homologs have conserved or diversified functions. Hydrogen sulfide (H₂S) is a plant gasotransmitter that plays important physiological roles in growth, development, and responses to environmental stresses. The modification of cysteines through S-sulfhydration is an important mechanism of H₂S, which regulates protein functions. H₂S promotes flowering in Arabidopsis and heading Chinese cabbage. Here we investigated the molecular mechanisms of H₂S used to promote flowering in the latter. Four, five, and four BraFLC, BraSOC I, and BraFT homologs were identified in heading Chinese cabbage. Different BraFLC proteins were bound to different CArG boxes in the promoter regions of the BraSOC I and BraFT homologs, producing different binding patterns. Thus, there may be functionally diverse BraFLC homologs in heading Chinese cabbage. Exogenous H₂S at 100 μmol L^-1 significantly promoted flowering by compensating for insufficient vernalization. BraFLC 1 and BraFLC 3 underwent S-sulfhydration by H₂S, after which their abilities to bind most BraSOC I or BraFT promoter probes weakened or even disappeared. These changes in binding ability were consistent with the expression pattern of the BraFT and BraSOC I homologs in seedlings treated with H₂S. These results indicated that H₂S signaling regulates flowering time. In summary, H₂S signaling promoted plant flowering by weakening or eliminating the binding abilities of BraFLCs to downstream promoters through S-sulfhydration.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Identification of Chromosomes and Chromosome Rearrangements in Crop Brassicas and Raphanus sativus: A Cytogenetic Toolkit Using Synthesized Massive Oligonucleotide Libraries

Crop brassicas include three diploid [Brassica rapa (AA; 2n = 2x = 16), B. nigra (BB; 2n = 2x = 18), and B. oleracea (CC; 2n = 2x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa. Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500-1,000 kb) with 1,997-5,420 oligonucleotides localized at each site in B. rapa. For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R (Raphanus sativus) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.

A highly contiguous genome assembly of Brassica nigra (BB) and revised nomenclature for the pseudochromosomes

Background

Brassica nigra (BB), also called black mustard, is grown as a condiment crop in India. B. nigra represents the B genome of U’s triangle and is one of the progenitor species of B. juncea (AABB), an important oilseed crop of the Indian subcontinent. We report the genome assembly of B. nigra variety Sangam.

Results

The genome assembly was carried out using Oxford Nanopore long-read sequencing and optical mapping. A total of 1549 contigs were assembled, which covered ~ 515.4 Mb of the estimated ~ 522 Mb of the genome. The final assembly consisted of 15 scaffolds that were assigned to eight pseudochromosomes using a high-density genetic map of B. nigra. Around 246 Mb of the genome consisted of the repeat elements; LTR/Gypsy types of retrotransposons being the most predominant. The B genome-specific repeats were identified in the centromeric regions of the B. nigra pseudochromosomes. A total of 57,249 protein-coding genes were identified of which 42,444 genes were found to be expressed in the transcriptome analysis. A comparison of the B genomes of B. nigra and B. juncea revealed high gene colinearity and similar gene block arrangements. A comparison of the structure of the A, B, and C genomes of U’s triangle showed the B genome to be divergent from the A and C genomes for gene block arrangements and centromeric regions.

Conclusions

A highly contiguous genome assembly of the B. nigra genome reported here is an improvement over the previous short-read assemblies and has allowed a comparative structural analysis of the A, B, and C genomes of the species belonging to the U’s triangle. Based on the comparison, we propose a new nomenclature for B. nigra pseudochromosomes, taking the B. rapa pseudochromosome nomenclature as the reference.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07271-w.

Improving seed size, seed weight and seedling emergence in Camelina sativa by overexpressing the Atsob3-6 gene variant

Seedling stand establishment is a critical factor affecting crop yield in low-precipitation agricultural regions. This is especially true for small seeded crops, such as Camelina (Camelina sativa) and canola (Brassica napus), that need to be planted shallow. Deeper planting would be desirable so that seeds can access soil moisture and bigger seeds could improve emergence and stand establishment by providing the energy necessary for seedling elongation. AHL (AT-Hook Containing, Nuclear Localized) genes play an important role in seedling growth and development. AHL proteins contain two structural units, the DNA-binding AT-hook motif and the Plant and Prokaryote Conserved (PPC) domain, required for protein-protein interactions. Our previous studies demonstrate that AtAHL29/SOB3 (Suppressor of phytochrome B-4 #3) regulates seedling development in Arabidopsis (Arabidopsis thaliana). Activation-tagged overexpression of AtSOB3 (Atsob3-D) represses the long-hypocotyl phenotype of an Arabidopsis phytochrome B mutant. In contrast, overexpression of the Atsob3-6 variant (Atsob3-6-OX), with a non-functional AT-hook, confers a long-hypocotyl phenotype. In this study, we demonstrate the role of Atsob3-D and Atsob3-6-OX in modulating seed size and hypocotyl length in the brassicas Arabidopsis and Camelina. In Arabidopsis, Atsob3-D reduces seed weight whereas Atsob3-6-OX increases seed weight and size when compared to the wild type. Similarly, Atsob3-6-OX transgenic Camelina seedlings are taller than the wild type, and produce larger and heavier seeds. These larger Atsob3-6-OX Camelina seeds also confer better emergence in deep-soil planting when compared to the wild type. Taken together, Atsob3-6-OX increases seed size, seed weight, seedling hypocotyl length and stand establishment in the oilseed crop Camelina.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Sex Differences in the Recombination Landscape

Sex differences in overall recombination rates are well known, but little theoretical or empirical attention has been given to how and why sexes differ in their recombination landscapes: the patterns of recombination along chromosomes. In the first scientific review of this phenomenon, we find that recombination is biased toward telomeres in males and more uniformly distributed in females in most vertebrates and many other eukaryotes. Notable exceptions to this pattern exist, however. Fine-scale recombination patterns also frequently differ between males and females. The molecular mechanisms responsible for sex differences remain unclear, but chromatin landscapes play a role. Why these sex differences evolve also is unclear. Hypotheses suggest that they may result from sexually antagonistic selection acting on coding genes and their regulatory elements, meiotic drive in females, selection during the haploid phase of the life cycle, selection against aneuploidy, or mechanistic constraints. No single hypothesis, however, can adequately explain the evolution of sex differences in all cases. Sex-specific recombination landscapes have important consequences for population differentiation and sex chromosome evolution.

Evolution of MIR159/319 genes in Brassica campestris and their function in pollen development

MIR159/319 have conserved evolution and diversified function after WGT in Brassica campestris, both of them can lead pollen vitality and germination abnormality, Bra-MIR319c also can function in flower development. MiR159 and miR319 are extensively studied highly conserved microRNAs which play roles in vegetative development, reproduction, and hormone regulation. In this study, the effects of whole-genome triplication (WGT) on the evolution of the MIR159/319 family and the functional diversification of the genes were comprehensively investigated in Brassica campestris. We identified 11 MIR159/319 genes in B. campestris, which produced five mature sequences. After analyzing the precursor sequences and phylogenetic tree, we found that Bra-MIR159/319 have evolutionary conservatism. Furthermore, Bra-MIR159/319 show functional diversification after WGT, as indicated by their expression patterns and the cis-element in their promoter. GUS signal showed that Bra-MIR159a and Bra-MIR319c can be expressed in anther but in different development stages. In B. campestris, overexpressed MIR159a and MIR319c contribute to late anther development and promote pollen abortion. Moreover, Bra-MIR319c can partially assume the function of MIR319a in flower development.

Characterization of differentially expressed genes to Cu stress in Brassica nigra by Arabidopsis genome arrays

Phytoremediation is an efficient and promising cleanup technology to extract or inactivate heavy metals and several organic and inorganic pollutants from soil and water. In this study, different Brassica nigra L. ecotypes, including Diyarbakır, collected from mining areas were exposed to different concentrations of copper and harvested after 72 h of Cu stress for the assessment of phytoremediation capacity. The Diyarbakır ecotype was called as "metallophyte" because of surviving at 500 μM Cu. To better understand Cu stress mechanism, ArabidopsisATH1 genome array was used to compare the gene expression in root and shoot tissues of B. nigra under 25 μM Cu. The response to Cu was much stronger in roots (88 genes showing increased or decreased mRNA levels) than in leaf tissues (24 responding genes). These genes were classified into the metal transport and accumulation-related genes, signal transduction and metabolism-related genes, and transport facilitation genes. Glutathione pathway-related genes (γ-ECS, PC, etc.) mRNAs were identified as differentially expressed in root and shoot tissues. QRT-PCR validation experiments showed that γ-ECS and PC expression was upregulated in the shoot and leaf tissues of the 100 μM Cu-subjected B. nigra-tolerant ecotype. This is the first study showing global expression profiles in response to Cu stress in B. nigra by Arabidopsis genome array. This work presented herein provides a well-illustrated insight into the global gene expression to Cu stress response in plants, and identified genes from microarray data will serve as molecular tools for the phytoremediation applications in the future.

Genome-Wide Analysis of Parent-of-Origin Allelic Expression in Endosperms of Brassicaceae Species, Brassica rapa

Uniparental gene expression, observed in both animals and plants, is termed genomic imprinting. Genomic imprinting is a well-known epigenetic phenomenon regulated through epigenetic modifications such as DNA methylation and histone modifications. Recent genome-wide studies of endosperm transcription have revealed the rapid change of imprinted genes between species, suggesting the flexibility of this phenomenon. Although the functional significance and evolutionary trends of imprinted genes are still obscure, it can be clarified by inter-species comparisons. In this study, we analyzed the pattern of genomic imprinting in Brassica rapa, a species related to Arabidopsis thaliana. Compared with the ancient karyotype of A. thaliana and B. rapa, B. rapa has a triplicated genome. Many imprinted genes, beyond the estimated number previously reported in other species, were observed. Several imprinted genes have been conserved among species in Brassicaceae. We also observed rapid molecular evolution of imprinted genes compared to non-imprinted genes in B. rapa. Especially, imprinted gene overlapping between species showed more rapid molecular evolution and preferential expression in endosperms. It may imply that a small number of imprinted genes have retained functional roles among diverged species and have been the target of natural selection.

Microsynteny analysis to understand evolution and impact of polyploidization on MIR319 family within Brassicaceae

The availability of a large number of whole-genome sequences allows comparative genomic analysis to reveal and understand evolution of regulatory regions and elements. The role played by events such as whole-genome and segmental duplications followed by genome fractionation in shaping genomic landscape and in expansion of gene families is crucial toward developing insights into evolutionary trends and consequences such as sequence and functional diversification. Members of Brassicaceae are known to have experienced several rounds of whole-genome duplication (WGD) that have been termed as paleopolyploidy, mesopolyploidy, and neopolyploidy. Such repeated events led to the creation and expansion of a large number of gene families. MIR319 is reported to be one of the most ancient and conserved plant MIRNA families and plays a role in growth and development including leaf development, seedling development, and embryo patterning. We have previously reported functional diversification of members of miR319 in Brassica oleracea affecting leaf architecture; however, the evolutionary history of the MIR319 gene family across Brassicaceae remains unknown and requires investigation. We therefore identified homologous and homeologous segments of ca. 100 kb, with or without MIR319, performed comparative synteny analysis and genome fractionation studies. We detected variable rates of gene retention across members of Brassicaceae when genomic blocks of MIR319a, MIR319b, and MIR319c were compared either between themselves or against Arabidopsis thaliana genome which was taken as the base genome. The highest levels of shared genes were found between A. thaliana and Capsella rubella in both MIR319b- and MIR319c-containing genomic segments, and with the closest species of A. thaliana, A. lyrata, only in MIR319a-containing segment. Synteny analysis across 12 genomes (with 30 sub-genomes) revealed MIR319c to be the most conserved MIRNA loci (present in 27 genomes/sub-genomes) followed by MIR319a (present in 23 genomes/sub-genomes); MIR319b was found to be frequently lost (present in 20 genomes/sub-genomes) and thus is under least selection pressure for retention. Genome fractionation revealed extensive and differential loss of MIRNA homeologous loci and flanking genes from various sub-genomes of Brassica species that is in accordance with their older history of polyploidy when compared to Camelina sativa, a recent neopolyploid, where the effect of genome fractionation was least. Finally, estimation of phylogenetic relationship using precursor sequences of MIR319 reveals MIR319a and MIR319b form sister clades, with MIR319c forming a separate clade. An intra-species synteny analysis between MIR319a-, MIR319b-, and MIR319c-containing genomic segments suggests segmental duplications at the base of Brassicaceae to be responsible for the origin of MIR319a and MIR319b.

Histone lysine methyltransferases BnaSDG8.A and BnaSDG8.C are involved in the floral transition in Brassica napus

Although increasing experimental evidence demonstrates that histone methylations play important roles in Arabidopsis plant growth and development, little information is available regarding Brassica napus. In this study, we characterized two genes encoding homologues of the Arabidopsis histone 3 lysine 36 (H3K36) methyltransferase SDG8, namely, BnaSDG8.A and BnaSDG8.C. Although no duplication of SDG8 homologous genes had been previously reported to occur during the evolution of any sequenced species, a domain-duplication was uncovered in BnaSDG8.C. This duplication led to the identification of a previously unknown NNH domain in the SDG8 homologues, providing a useful reference for future studies and revealing the finer mechanism of SDG8 function. One NNH domain is present in BnaSDG8.A, while two adjacent NNH domains are present in BnaSDG8.C. Reverse transcriptase-quantitative polymerase chain reaction analysis revealed similar patterns but with varied levels of expression of BnaSDG8.A/C in different plant organs/tissues. To directly investigate their function, BnaSDG8.A/C cDNA was ectopically expressed to complement the Arabidopsis mutant. We observed that the expression of either BnaSDG8.A or BnaSDG8.C could rescue the Arabidopsis sdg8 mutant to the wild-type phenotype. Using RNAi and CRISPR/Cas9-mediated gene editing, we obtained BnaSDG8.A/C knockdown and knockout mutants with the early flowering phenotype as compared with the control. Further analysis of two types of the mutants revealed that BnaSDG8.A/C are required for H3K36 m2/3 deposition and prevent the floral transition of B. napus by directly enhancing the H3K36 m2/3 levels at the BnaFLC chromatin loci. This observation on the floral transition by epigenetic modification in B. napus provides useful information for breeding early-flowering varieties.

Interactions of WRKY15 and WRKY33 transcription factors and their roles in the resistance of oilseed rape to Sclerotinia infection

WRKY transcription factors are known to participate in the defence responses of higher plants. However, little is known about the roles of such proteins, especially regarding their functions in the resistance of oilseed rape (Brassica napus) to Sclerotinia sclerotiorum, a necrotrophic fungal pathogen that causes stem rot. In this study, we identified BnWRKY33 as a S. sclerotiorum-responsive gene that positively regulates resistance to this pathogen by enhancing the expression of genes involved in camalexin synthesis and genes regulated by salicylic acid (SA) and jasmonic acid (JA). We also identified a S. sclerotiorum-responsive region in the promoter of BnWRKY33, which we revealed to be a relatively conserved W-box region in the promoters of homologous genes in different species. Using this S. sclerotiorum-responsive region as bait in a yeast one-hybrid assay, we identified another WRKY transcription factor, BnWRKY15, and observed that both BnWRKY15 and BnWRKY33 could bind to this region. In addition, BnWRKY15 overexpression simultaneously increased the susceptibility of B. napus to S. sclerotiorum and down-regulated BnWRKY33 after different durations of infection. Furthermore, BnWRKY15, which contains a transcriptional repression domain, exhibited reduced transactivation ability and could reduce the transactivation ability of BnWRKY33 in Arabidopsis protoplast assays. Therefore, we suggest that the increased susceptibility of BnWRKY15-overexpressing plants results from reduced BnWRKY33 expression, which is due to the inhibition of BnWRKY33 transcriptional activation by BnWRKY15.

Functional characterization of a heterologously expressed Brassica napus WRKY41-1 transcription factor in regulating anthocyanin biosynthesis in Arabidopsis thaliana

Previous studies have shown that a plant WRKY transcription factor, WRKY41, has multiple functions, and regulates seed dormancy, hormone signaling pathways, and both biotic and abiotic stress responses. However, it is not known about the roles of AtWRKY41 from the model plant, Arabidopsis thaliana, and its ortholog, BnWRKY41, from the closely related and important oil-producing crop, Brassica napus, in the regulation of anthocyanin biosynthesis. Here, we found that the wrky41 mutation in A. thaliana resulted in a significant increase in anthocyanin levels in rosette leaves, indicating that AtWRKY41 acts as repressor of anthocyanin biosynthesis. RNA sequencing and quantitative real-time PCR analysis revealed increased expression of three regulatory genes AtMYB75, AtMYB111, and AtMYBD, and two structural genes, AT1G68440 and AtGSTF12, all of which contribute to anthocyanin biosynthesis, in the sixth rosette leaves of wrky41-2 plants at 20 days after germination. We cloned the full length complementary DNA of BnWRKY41-1 from the C2 subgenome of the B. napus genotype Westar and observed that, when overexpressed in tobacco leaves as a fusion protein with green fluorescent protein, BnWRKY41-1 is localized to the nucleus. We further showed that overexpression of BnWRKY41-1 in the A. thaliana wrky41-2 mutant rescued the higher anthocyanin content phenotype in rosette leaves of the mutant. Moreover, the elevated expression levels in wrky41-2 rosette leaves of several important regulatory and structural genes regulating anthocyanin biosynthesis were not observed in the BnWRKY41-1 overexpressing lines. These results reveal that BnWRKY41-1 has a similar role with AtWRKY41 in regulating anthocyanin biosynthesis when overexpressed in A. thaliana. This gene represents a promising target for genetically manipulating B. napus to increase the amounts of anthocyanins in rosette leaves.

Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus

Background

Deciphering the genetic architecture of a species is a good way to understand its evolutionary history, but also to tailor its profile for breeding elite cultivars with desirable traits. Aligning QTLs from diverse population in one map and utilizing it for comparison, but also as a basis for multiple analyses assure a stronger evidence to understand the genetic system related to a given phenotype.

Results

In this study, 439 genes involved in fatty acid (FA) and triacylglycerol (TAG) biosyntheses were identified in Brassica napus. B. napus genome showed mixed gene loss and insertion compared to B. rapa and B. oleracea, and C genome had more inserted genes. Identified QTLs for oil (OC-QTLs) and fatty acids (FA-QTLs) from nine reported populations were projected on the physical map of the reference genome “Darmor-bzh” to generate a map. Thus, 335 FA-QTLs and OC-QTLs could be highlighted and 82 QTLs were overlapping. Chromosome C3 contained 22 overlapping QTLs with all trait studied except for C18:3. In total, 218 candidate genes which were potentially involved in FA and TAG were identified in 162 QTLs confidence intervals and some of them might affect many traits. Also, 76 among these candidate genes were found inside 57 overlapping QTLs, and candidate genes for oil content were in majority (61/76 genes). Then, sixteen genes were found in overlapping QTLs involving three populations, and the remaining 60 genes were found in overlapping QTLs of two populations. Interaction network and pathway analysis of these candidate genes indicated ten genes that might have strong influence over the other genes that control fatty acids and oil formation.

Conclusion

The present results provided new information for genetic basis of FA and TAG formation in B. napus. A map including QTLs from numerous populations was built, which could serve as reference to study the genome profile of B. napus, and new potential genes emerged which might affect seed oil. New useful tracks were showed for the selection of population or/and selection of interesting genes for breeding improvement purpose.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4176-6) contains supplementary material, which is available to authorized users.

Functional characterization of Brassica napus DNA topoisomerase Iα-1 and its effect on flowering time when expressed in Arabidopsis thaliana

Previous studies have shown that DNA topoisomerase Iα (AtTOP1α) has specific developmental functions during growth and development in Arabidopsis thaliana. However, little is known about the roles of DNA topoisomerases in the closely related and commercially important plant, rapeseed (Brassica napus). Here, the full-length BnTOP1α-1 coding sequence was cloned from the A2 subgenome of the Brassica napus inbred line L111. We determine that all BnTOP1α paralogs showed differing patterns of expression in different organs of L111, and that when expressed in tobacco leaves as a fusion protein with green fluorescent protein, BnTOP1α-1 localized to the nucleus. We further showed that ectopic expression of BnTOP1α-1 in the A. thaliana top1α-7 mutant fully complemented the early flowering phenotype of the mutant. Moreover, altered expression levels in top1α-7 seedlings of several key genes controlling flowering time were restored to wild type levels by ectopic expression of BnTOP1α-1. These results provide valuable insights into the roles of rapeseed DNA topoisomerases in flowering time, and provide a promising target for genetic manipulation of this commercially significant process in rapeseed.

Cytogenetic and Molecular Characterization of B-Genome Introgression Lines of Brassica napus L

Brassica napus introgression lines (ILs), having B-genome segments from B. carinata, were assessed genetically for extent of introgression and phenotypically for siliqua shatter resistance. Introgression lines had 7-9% higher DNA content, were meiotically stable, and had almost normal pollen fertility/seed set. Segment introgressions were confirmed by fluorescent genomic in situ hybridization (fl-GISH), SSR analyses, and SNP studies. Genotyping with 48 B-genome specific SSRs detected substitutions from B3, B4, B6, and B7 chromosomes on 39 of the 69 ILs whereas SNP genotyping detected a total of 23 B-segments (≥3 Mb) from B4, B6, and B7 introgressed into 10 of the 19 (C1, C2, C3, C5, C6, C8, C9, A3, A9, A10) chromosomes in 17 ILs. The size of substitutions varied from 3.0 Mb on chromosome A9 (IL59) to 42.44 Mb on chromosome C2 (IL54), ranging from 7 to 83% of the recipient chromosome. Average siliqua strength in ILs was observed to be higher than that of B. napus parents (2.2-6.0 vs. 1.9-4.0 mJ) while siliqua strength in some of the lines was almost equal to that of the donor parent B. carinata (6.0 vs.7.2 mJ). These ILs, with large chunks of substituted B-genome, can prove to be a useful prebreeding resource for germplasm enhancement in B. napus, especially for siliqua shatter resistance.

Conserved Function of ACYL-ACYL CARRIER PROTEIN DESATURASE 5 on Seed Oil and Oleic Acid Biosynthesis between Arabidopsis thaliana and Brassica napus

Previous studies have shown that several ACYL-ACYL CARRIER PROTEIN DESATURASE (AtAAD) members in Arabidopsis thaliana are responsible for oleic acid (C18:1) biosynthesis. Limited research has been conducted on another member, AtAAD5, and its paralog BnAAD5 in the closely related and commercially important plant, Brassica napus. Here, we found that AtAAD5 was predominantly and exclusively expressed in developing embryos at the whole seed developmental stages. The aad5 mutation caused a significant decrease in the amounts of oil and C18:1, and a considerable increase in the content of stearic acid (C18:0) in mature seeds, suggesting that AtAAD5 functioned as an important facilitator of seed oil biosynthesis. We also cloned the full-length coding sequence of BnAAD5-1 from the A3 subgenome of the B. napus inbred line L111. We showed that ectopic expression of BnAAD5-1 in the A. thaliana aad5-2 mutant fully complemented the phenotypes of the mutant, such as lower oil content and altered contents of C18:0 and C18:1. These results help us to better understand the functions of AAD members in A. thaliana and B. napus and provide a promising target for genetic manipulation of B. napus.

Cloning and Phylogenetic Analysis of Brassica napus L. Caffeic Acid O-Methyltransferase 1 Gene Family and Its Expression Pattern under Drought Stress

For many plants, regulating lignin content and composition to improve lodging resistance is a crucial issue. Caffeic acid O-methyltransferase (COMT) is a lignin monomer-specific enzyme that controls S subunit synthesis in plant vascular cell walls. Here, we identified 12 BnCOMT1 gene homologues, namely BnCOMT1-1 to BnCOMT1-12. Ten of 12 genes were composed of four highly conserved exons and three weakly conserved introns. The length of intron I, in particular, showed enormous diversification. Intron I of homologous BnCOMT1 genes showed high identity with counterpart genes in Brassica rapa and Brassica oleracea, and intron I from positional close genes in the same chromosome were relatively highly conserved. A phylogenetic analysis suggested that COMT genes experience considerable diversification and conservation in Brassicaceae species, and some COMT1 genes are unique in the Brassica genus. Our expression studies indicated that BnCOMT1 genes were differentially expressed in different tissues, with BnCOMT1-4, BnCOMT1-5, BnCOMT1-8, and BnCOMT1-10 exhibiting stem specificity. These four BnCOMT1 genes were expressed at all developmental periods (the bud, early flowering, late flowering and mature stages) and their expression level peaked in the early flowering stage in the stem. Drought stress augmented and accelerated lignin accumulation in high-lignin plants but delayed it in low-lignin plants. The expression levels of BnCOMT1s were generally reduced in water deficit condition. The desynchrony of the accumulation processes of total lignin and BnCOMT1s transcripts in most growth stages indicated that BnCOMT1s could be responsible for the synthesis of a specific subunit of lignin or that they participate in other pathways such as the melatonin biosynthesis pathway.

Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

Plant genomes are characterized by the presence of large miRNA gene families which are few in number. The expansion of miRNA families is thought to be driven by gene and genome duplication. Some members of these miRNA gene families are tandemly arranged and their analysis is of interest because such organisation may indicate origin through tandem duplication and also to investigate whether some such tandem clusters have similar expression patterns, and whether these are regulated through a common set of cis-regulatory elements (eg. promoters and enhancers). As a first step, we undertake a comprehensive study using micro-synteny analyses of tandemly organised miRNA families across the Brassicaceae spanning an evolutionary time scale of ca. 45 million years, among Arabidopsis, Capsella, Brassica and Thellungiella species, to address the following questions: Are most miRNA gene families present as tandem clusters? To what extent are these tandem patterns retained? To what extent can family sizes be ascribed to genome duplication? Our analysis of thirteen tandemly organised miRNA families revealed that synteny is largely conserved among Arabidopsis thaliana, A. lyrata and Capsella rubella, which form a clade spanning approximately between 6.2-9.8 my (Acarkan et al., 2000) [1]. On the other hand, comparison of sequences from these species with Brassica rapa, B. oleracea and Thellungiella halophila, which form a separate clade spanning 31 my (Franzke et al., 2011)[2] reveals many differences. The latter clade reveals several paralogous duplications that probably resulted from whole genome duplication, as well as disrupted synteny. Phylogenetic analyses of precursor sequences generally support the history inferred from synteny analysis. Synteny and phylogenetic analysis of six members of the tandemly organised miR169 family suggest that the Brassicaceae ancestral state consisted of a "dimer as a unit" which may have undergone direct local duplication to retain the transcriptional orientation followed by lineage specific changes. MiR169, to the best of our knowledge, is one of the largest tandemly organised miRNA gene family across plant kingdom and further analysis should reveal the generality of this pattern of evolution. The conserved organisation of miR395A-B-C and miR395 D-E-F as two clusters on same chromosome/scaffold across A. thaliana, B. rapa and salsuginea demonstrates retention of the large chromosomal segment across the two lineages. MiRNA family miR845 was detected only in Arabidopsis species and Thellungiella indicating a complex loss and retention history. MiR447A-B family was only found in A. thaliana indicating that it is a species-specific gene family of recent origin.

Resynthesis of Brassica napus through hybridization between B. juncea and B. carinata

First report for the resynthesis of Brassica napus by recombining A and C genome from B. juncea and B. carinata , respectively. Also documents B genome introgressions in resynthesized B. napus. Resynthesis of Brassica napus (AACC) was achieved by hybridizing Brassica juncea (AABB) with Brassica carinata (BBCC). This was facilitated by spontaneous chromosome doubling in the F1 hybrid (ABBC) to yield octaploid (AABBBBCC), elimination of extra B genome chromosomes in the resulting octaploid and in subsequent selfed generations, aided with directed selection for fertile plants having B. napus morphology. Twenty-five plants with varying degrees of resemblance to natural B. napus were identified from 17 A5 progenies and assayed for cytogenetic stability and genetic diversity. Majority of these plants, except six (2n = 38) were hyperploids (2n = 40-56). The six plants with 2n = 38 were designated as derived B. napus types. These showed an expected meiotic configuration of 19II at metaphase-I, with 19-19 distribution at anaphase-I. Genotyping based on A and C genome specific primers confirmed genetic identity of six derived (2n = 38) B. napus plants with natural types whereas genotyping with B genome specific primers indicated introgression of B genome segments. This was also confirmed by genomic in situ hybridization (GISH). Strong signals of B genome probe were detected, proving hitherto unreported genetic exchanges between B and A/C chromosomes. These introgressions possibly occurred en route five generations of selfing. Derived plants yielded fertile hybrids in crosses with natural B. napus var. GSC 6. The selfed derived plants as evaluated in A6 plant to progeny rows were morphologically similar to natural B. napus, and meiotically stable. Agronomic assessment of these progenies revealed variation for key morpho-physiological traits. Of special interest were the progenies with plants having oil content exceeding 47% as against about 39-41% in existing cultivars.

Co-linearity and divergence of the A subgenome of Brassica juncea compared with other Brassica species carrying different A subgenomes

Background

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, A^rA^r; B. nigra, BⁿⁱBⁿⁱ; B. oleracea, C^oC^o) and the derived allotetraploid species (i.e.: B. juncea, A^jA^jB^jB^j; B. napus, AⁿAⁿCⁿCⁿ; B. carinata, B^cB^cC^cC^c). 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).

Results

A dense genetic linkage map of B. juncea was constructed using an F₂ 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 A^j-Aⁿ than between Aⁿ-A^r and A^j-A^r subgenome pairs.

Conclusions

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 A^r, and AⁿCⁿ genomes allowed us to compare the A subgenome structure of B. juncea (A^j) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (A^r, Aⁿ and A^j) 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2343-1) contains supplementary material, which is available to authorized users.

Molecular Mapping and QTL for Expression Profiles of Flavonoid Genes in Brassica napus

Flavonoids are secondary metabolites that are extensively distributed in the plant kingdom and contribute to seed coat color formation in rapeseed. To decipher the genetic networks underlying flavonoid biosynthesis in rapeseed, we constructed a high-density genetic linkage map with 1089 polymorphic loci (including 464 SSR loci, 97 RAPD loci, 451 SRAP loci, and 75 IBP loci) using recombinant inbred lines (RILs). The map consists of 19 linkage groups and covers 2775 cM of the B. napus genome with an average distance of 2.54 cM between adjacent markers. We then performed expression quantitative trait locus (eQTL) analysis to detect transcript-level variation of 18 flavonoid biosynthesis pathway genes in the seeds of the 94 RILs. In total, 72 eQTLs were detected and found to be distributed among 15 different linkage groups that account for 4.11% to 52.70% of the phenotypic variance atrributed to each eQTL. Using a genetical genomics approach, four eQTL hotspots together harboring 28 eQTLs associated with 18 genes were found on chromosomes A03, A09, and C08 and had high levels of synteny with genome sequences of A. thaliana and Brassica species. Associated with the trans-eQTL hotspots on chromosomes A03, A09, and C08 were 5, 17, and 1 genes encoding transcription factors, suggesting that these genes have essential roles in the flavonoid biosynthesis pathway. Importantly, bZIP25, which is expressed specifically in seeds, MYC1, which controls flavonoid biosynthesis, and the R2R3-type gene MYB51, which is involved in the synthesis of secondary metabolites, were associated with the eQTL hotspots, and these genes might thus be involved in different flavonoid biosynthesis pathways in rapeseed. Hence, further studies of the functions of these genes will provide insight into the regulatory mechanism underlying flavonoid biosynthesis, and lay the foundation for elaborating the molecular mechanism of seed coat color formation in B. napus.

Diversification and evolution of the SDG gene family in Brassica rapa after the whole genome triplication

Histone lysine methylation, controlled by the SET Domain Group (SDG) gene family, is part of the histone code that regulates chromatin function and epigenetic control of gene expression. Analyzing the SDG gene family in Brassica rapa for their gene structure, domain architecture, subcellular localization, rate of molecular evolution and gene expression pattern revealed common occurrences of subfunctionalization and neofunctionalization in BrSDGs. In comparison with Arabidopsis thaliana, the BrSDG gene family was found to be more divergent than AtSDGs, which might partly explain the rich variety of morphotypes in B. rapa. In addition, a new evolutionary pattern of the four main groups of SDGs was presented, in which the Trx group and the SUVR subgroup evolved faster than the E(z), Ash groups and the SUVH subgroup. These differences in evolutionary rate among the four main groups of SDGs are perhaps due to the complexity and variability of the regions that bind with biomacromolecules, which guide SDGs to their target loci.

Speeding up chromosome evolution in Phaseolus: multiple rearrangements associated with a one-step descending dysploidy

The genus Phaseolus L. has been subject of extensive cytogenetic studies due to its global economic importance. It is considered karyotypically stable, with most of its ca. 75 species having 2n = 22 chromosomes, and only three species (Phaseolus leptostachyus, Phaseolus macvaughii, and Phaseolus micranthus), which form the Leptostachyus clade, having 2n = 20. To test whether a simple chromosomal fusion was the cause of this descending dysploidy, mitotic chromosomes of P. leptostachyus (2n = 20) were comparatively mapped by fluorescent in situ hybridization (FISH) using bacterial artificial chromosomes (BACs) and ribosomal DNA (rDNA) probes. Our results corroborated the conservation of the 5S and 45S rDNA sites on ancestral chromosomes 10 and 6, respectively. The reduction from x = 11 to x = 10 was the result of the insertion of chromosome 10 into the centromeric region of chromosome 11, supporting a nested chromosome fusion (NCF) as the main cause of this dysploidy. Additionally, the terminal region of the long arm of chromosome 6 was translocated to this larger chromosome. Surprisingly, the NCF was accompanied by several additional translocations and inversions previously unknown for the genus, suggesting that the dysploidy may have been associated to a burst of genome reorganization in this otherwise stable, diploid plant genus.

Molecular breeding in Brassica for salt tolerance: importance of microsatellite (SSR) markers for molecular breeding in Brassica

Salinity is one of the important abiotic factors for any crop management in irrigated as well as rainfed areas, which leads to poor harvests. This yield reduction in salt affected soils can be overcome by improving salt tolerance in crops or by soil reclamation. Salty soils can be reclaimed by leaching the salt or by cultivation of salt tolerance crops. Salt tolerance is a quantitative trait controlled by several genes. Poor knowledge about mechanism of its inheritance makes slow progress in its introgression into target crops. Brassica is known to be a good reclamation crop. Inter and intra specific variation within Brassica species shows potential of molecular breeding to raise salinity tolerant genotypes. Among the various molecular markers, SSR markers are getting high attention, since they are randomly sparsed, highly variable and show co-dominant inheritance. Furthermore, as sequencing techniques are improving and softwares to find SSR markers are being developed, SSR markers technology is also evolving rapidly. Comparative SSR marker studies targeting Arabidopsis thaliana and Brassica species which lie in the same family will further aid in studying the salt tolerance related QTLs and subsequent identification of the "candidate genes" and finding out the origin of important QTLs. Although, there are a few reports on molecular breeding for improving salt tolerance using molecular markers in Brassica species, usage of SSR markers has a big potential to improve salt tolerance in Brassica crops. In order to obtain best harvests, role of SSR marker driven breeding approaches play important role and it has been discussed in this review especially for the introgression of salt tolerance traits in crops.

The radish genome and comprehensive gene expression profile of tuberous root formation and development

Understanding the processes that regulate plant sink formation and development at the molecular level will contribute to the areas of crop breeding, food production and plant evolutionary studies. We report the annotation and analysis of the draft genome sequence of the radish Raphanus sativus var. hortensis (long and thick root radish) and transcriptome analysis during root development. Based on the hybrid assembly approach of next-generation sequencing, a total of 383 Mb (N50 scaffold: 138.17 kb) of sequences of the radish genome was constructed containing 54,357 genes. Syntenic and phylogenetic analyses indicated that divergence between Raphanus and Brassica coincide with the time of whole genome triplication (WGT), suggesting that WGT triggered diversification of Brassiceae crop plants. Further transcriptome analysis showed that the gene functions and pathways related to carbohydrate metabolism were prominently activated in thickening roots, particularly in cell proliferating tissues. Notably, the expression levels of sucrose synthase 1 (SUS1) were correlated with root thickening rates. We also identified the genes involved in pungency synthesis and their transcription factors.

Neofunctionalization of duplicated Tic40 genes caused a gain-of-function variation related to male fertility in Brassica oleracea lineages

Gene duplication followed by functional divergence in the event of polyploidization is a major contributor to evolutionary novelties. The Brassica genus evolved from a common ancestor after whole-genome triplication. Here, we studied the evolutionary and functional features of Brassica spp. homologs to Tic40 (for translocon at the inner membrane of chloroplasts with 40 kDa). Four Tic40 loci were identified in allotetraploid Brassica napus and two loci in each of three basic diploid Brassica spp. Although these Tic40 homologs share high sequence identities and similar expression patterns, they exhibit altered functional features. Complementation assays conducted on Arabidopsis thaliana tic40 and the B. napus male-sterile line 7365A suggested that all Brassica spp. Tic40 homologs retain an ancestral function similar to that of AtTic40, whereas BolC9.Tic40 in Brassica oleracea and its ortholog in B. napus, BnaC9.Tic40, in addition, evolved a novel function that can rescue the fertility of 7365A. A homologous chromosomal rearrangement placed bnac9.tic40 originating from the A genome (BraA10.Tic40) as an allele of BnaC9.Tic40 in the C genome, resulting in phenotypic variation for male sterility in the B. napus near-isogenic two-type line 7365AB. Assessment of the complementation activity of chimeric B. napus Tic40 domain-swapping constructs in 7365A suggested that amino acid replacements in the carboxyl terminus of BnaC9.Tic40 cause this functional divergence. The distribution of these amino acid replacements in 59 diverse Brassica spp. accessions demonstrated that the neofunctionalization of Tic40 is restricted to B. oleracea and its derivatives and thus occurred after the divergence of the Brassica spp. A, B, and C genomes.

A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.

Doubling down on genomes: polyploidy and crop plants

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.

Intergenomic single nucleotide polymorphisms as a tool for bacterial artificial chromosome contig building of homoeologous Brassica napus regions

BACKGROUND

Homoeologous sequences pose a particular challenge if bacterial artificial chromosome (BAC) contigs shall be established for specific regions of an allopolyploid genome. Single nucleotide polymorphisms (SNPs) differentiating between homoeologous genomes (intergenomic SNPs) may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.

RESULTS

Sequence alignments between Brassica rapa (AA) and Brassica oleracea (CC) sequences mapping to corresponding regions on chromosomes A1 and C1, respectively were used to identify single nucleotide polymorphisms between the A and C genomes. A large fraction of these polymorphisms was also present in Brassica napus (AACC), an allopolyploid species that originated from hybridisation of A and C genome species. Intergenomic SNPs mapping throughout homoeologous chromosome segments spanning approximately one Mbp each were included in Illumina's GoldenGate® Genotyping Assay and used to screen multidimensional pools of a Brassica napus bacterial artificial chromosome library with tenfold genome coverage. Based on the results of 50 SNP assays, a BAC contig for the Brassica napus A subgenome was established that spanned the entire region of interest. The C subgenome region was represented in three BAC contigs.

CONCLUSIONS

This proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic SNPs suitable for multiplex polymerase chain reaction (PCR)-based screening of multidimensional BAC pools of a polyploid organism. Owing to their high abundance and ease of identification, intergenomic SNPs represent a versatile tool to establish BAC contigs for homoeologous regions of a polyploid genome.

Constructing a dense genetic linkage map and mapping QTL for the traits of flower development in Brassica carinata

An integrated dense genetic linkage map was constructed in a B. carinata population and used for comparative genome analysis and QTL identification for flowering time. An integrated dense linkage map of Brassica carinata (BBCC) was constructed in a doubled haploid population based on DArT-Seq(TM) markers. A total of 4,031 markers corresponding to 1,366 unique loci were mapped including 639 bins, covering a genetic distance of 2,048 cM. We identified 136 blocks and islands conserved in Brassicaceae, which showed a feature of hexaploidisation representing the suggested ancestral crucifer karyotype. The B and C genome of B. carinata shared 85 % of commonly conserved blocks with the B genome of B. nigra/B. juncea and 80 % of commonly conserved blocks with the C genome of B. napus, and shown frequent structural rearrangements such as insertions and inversions. Up to 24 quantitative trait loci (QTL) for flowering and budding time were identified in the DH population. Of these QTL, one consistent QTL (qFT.B4-2) for flowering time was identified in all of the environments in the J block of the B4 linkage group, where a group of genes for flowering time were aligned in A. thaliana. Another major QTL for flowering time under a winter-cropped environment was detected in the E block of C6, where the BnFT-C6 gene was previously localised in B. napus. This high-density map would be useful not only to reveal the genetic variation in the species with QTL analysis and genome sequencing, but also for other applications such as marker-assisted selection and genomic selection, for the African mustard improvement.