Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola.

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

A chromosome-scale assembly of Brassica carinata (BBCC) accession HC20 containing resistance to multiple pathogens and an early generation assessment of introgressions into B. juncea (AABB).

Brassica carinata (BBCC) commonly referred to as Ethiopian mustard is a natural allotetraploid containing the genomes of Brassica nigra (BB) and Brassica oleracea (CC). It is an oilseed crop endemic to the northeastern regions of Africa. Although it is under limited cultivation, B. carinata is valuable as it is resistant/highly tolerant to most of the pathogens affecting widely cultivated Brassica species of the U's triangle. We report a chromosome-scale genome assembly of B. carinata accession HC20 using long-read Oxford Nanopore sequencing and Bionano optical maps. The assembly has a scaffold N50 of ~39.8 Mb and covers ~1.11 Gb of the genome. We compared the long-read genome assemblies of the U's triangle species and found extensive gene collinearity between the diploids and allopolyploids with no evidence of major gene losses. Therefore, B. juncea (AABB), B. napus (AACC), and B. carinata can be regarded as strict allopolyploids. We cataloged the nucleotide-binding and leucine-rich repeat immune receptor (NLR) repertoire of B. carinata and, identified 465 NLRs, and compared these with the NLRs in the other Brassica species. We investigated the extent and nature of early-generation genomic interactions between the constituent genomes of B. carinata and B. juncea in interspecific crosses between the two species. Besides the expected recombination between the constituent B genomes, extensive homoeologous exchanges were observed between the A and C genomes. Interspecific crosses, therefore, can be used for transferring disease resistance from B. carinata to B. juncea and broadening the genetic base of the two allotetraploid species.

Genomic asymmetric epigenetic modification of transposable elements is involved in gene expression regulation of allopolyploid Brassica napus.

Polyploids are common and have a wide geographical distribution and environmental adaptability. Allopolyploidy may lead to the activation of transposable elements (TE). However, the mechanism of epigenetic modification of TEs in the establishment and evolution of allopolyploids remains to be explored. We focused on the TEs of model allopolyploid Brassica napus (An An Cn Cn ), exploring the TE characteristics of the genome, epigenetic modifications of TEs during allopolyploidization, and regulation of gene expression by TE methylation. In B. napus, approximately 50% of the genome was composed of TEs. TEs increased with proximity to genes, especially DNA transposons. TE methylation levels were negatively correlated with gene expression, and changes in TE methylation levels were able to regulate the expression of neighboring genes related to responses to light intensity and stress, which promoted powerful adaptation of allopolyploids to new environments. TEs can be synergistically regulated by RNA-directed DNA methylation pathways and histone modifications. The epigenetic modification levels of TEs tended to be similar to those of the diploid parents during the genome evolution of B. napus. The TEs of the An subgenome were more likely to be modified, and the imbalance in TE number and epigenetic modification level in the An and Cn subgenomes may lead to the establishment of subgenome dominance. Our study analyzed the characteristics of TE location, DNA methylation, siRNA, and histone modification in B. napus and highlighted the importance of TE epigenetic modifications during the allopolyploidy process, providing support for revealing the mechanism of allopolyploid formation and evolution.

Nitrate alleviates ammonium toxicity in Brassica napus by coordinating rhizosphere and cell pH and ammonium assimilation.

In natural and agricultural situations, ammonium ( NH 4 + ) is a preferred nitrogen (N) source for plants, but excessive amounts can be hazardous to them, known as NH 4 + toxicity. Nitrate ( NO 3 - ) has long been recognized to reduce NH 4 + toxicity. However, little is known about Brassica napus, a major oil crop that is sensitive to high NH 4 + . Here, we found that NO 3 - can mitigate NH 4 + toxicity by balancing rhizosphere and intracellular pH and accelerating ammonium assimilation in B. napus. NO 3 - increased the uptake of NO 3 - and NH 4 + under high NH 4 + circumstances by triggering the expression of NO 3 - and NH 4 + transporters, while NO 3 - and H+ efflux from the cytoplasm to the apoplast was enhanced by promoting the expression of NO 3 - efflux transporters and genes encoding plasma membrane H+ -ATPase. In addition, NO 3 - increased pH in the cytosol, vacuole, and rhizosphere, and down-regulated genes induced by acid stress. Root glutamine synthetase (GS) activity was elevated by NO 3 - under high NH 4 + conditions to enhance the assimilation of NH 4 + into amino acids, thereby reducing NH 4 + accumulation and translocation to shoot in rapeseed. In addition, root GS activity was highly dependent on the environmental pH. NO 3 - might induce metabolites involved in amino acid biosynthesis and malate metabolism in the tricarboxylic acid cycle, and inhibit phenylpropanoid metabolism to mitigate NH 4 + toxicity. Collectively, our results indicate that NO 3 - balances both rhizosphere and intracellular pH via effective NO 3 - transmembrane cycling, accelerates NH 4 + assimilation, and up-regulates malate metabolism to mitigate NH 4 + toxicity in oilseed rape.

Genome-wide association study reveals the genetic basis for petal-size formation in rapeseed (Brassica napus) and CRISPR-Cas9-mediated mutagenesis of BnFHY3 for petal-size reduction.

Petals in rapeseed (Brassica napus) serve multiple functions, including protection of reproductive organs, nutrient acquisition, and attraction of pollinators. However, they also cluster densely at the top, forming a thick layer that absorbs and reflects a considerable amount of photosynthetically active radiation. Breeding genotypes with large, small, or even petal-less varieties, requires knowledge of primary genes for allelic selection and manipulation. However, our current understanding of petal-size regulation is limited, and the lack of markers and pre-breeding materials hinders targeted petal-size breeding. Here, we conducted a genome-wide association study on petal size using 295 diverse accessions. We identified 20 significant single nucleotide polymorphisms and 236 genes associated with petal-size variation. Through a cross-analysis of genomic and transcriptomic data, we focused on 14 specific genes, from which molecular markers for diverging petal-size features can be developed. Leveraging CRISPR-Cas9 technology, we successfully generated a quadruple mutant of Far-Red Elongated Hypocotyl 3 (q-bnfhy3), which exhibited smaller petals compared to the wild type. Our study provides insights into the genetic basis of petal-size regulation in rapeseed and offers abundant potential molecular markers for breeding. The q-bnfhy3 mutant unveiled a novel role of FHY3 orthologues in regulating petal size in addition to previously reported functions.

Knockout of stigmatic ascorbate peroxidase 1 (APX1) delays pollen rehydration and germination by mediating ROS homeostasis in Brassica napus L.

The successful interaction between pollen and stigma is a critical process for plant sexual reproduction, involving a series of intricate molecular and physiological events. After self-compatible pollination, a significant reduction in reactive oxygen species (ROS) production has been observed in stigmas, which is essential for pollen grain rehydration and subsequent pollen tube growth. Several scavenging enzymes tightly regulate ROS homeostasis. However, the potential role of these ROS-scavenging enzymes in the pollen-stigma interaction in Brassica napus remains unclear. Here, we showed that the activity of ascorbate peroxidase (APX), an enzyme that plays a crucial role in the detoxification of hydrogen peroxide (H2O2), was modulated depending on the compatibility of pollination in B. napus. We then identified stigma-expressed APX1s and generated pentuple mutants of APX1s using CRISPR/Cas9 technology. After compatible pollination, the BnaAPX1 pentuple mutants accumulated higher levels of H2O2 in the stigma, while the overexpression of BnaA09.APX1 resulted in lower levels of H2O2. Furthermore, the knockout of BnaAPX1 delayed the compatible response-mediated pollen rehydration and germination, which was consistent with the effects of a specific APX inhibitor, ρ-Aminophenol, on compatible pollination. In contrast, the overexpression of BnaA09.APX1 accelerated pollen rehydration and germination after both compatible and incompatible pollinations. However, delaying and promoting pollen rehydration and germination did not affect the seed set after compatible and incompatible pollination in APX1 pentuple mutants and overexpression lines, respectively. Our results demonstrate the fundamental role of BnaAPX1 in pollen rehydration and germination by regulating ROS homeostasis during the pollen-stigma interaction in B. napus.

Trehalose-6-phosphate synthase 8 increases photosynthesis and seed yield in Brassica napus.

Trehalose-6-phosphate (T6P) functions as a vital proxy for assessing carbohydrate status in plants. While class II T6P synthases (TPS) do not exhibit TPS activity, they are believed to play pivotal regulatory roles in trehalose metabolism. However, their precise functions in carbon metabolism and crop yield have remained largely unknown. Here, BnaC02.TPS8, a class II TPS gene, is shown to be specifically expressed in mature leaves and the developing pod walls of Brassica napus. Overexpression of BnaC02.TPS8 increased photosynthesis and the accumulation of sugars, starch, and biomass compared to wild type. Metabolomic analysis of BnaC02.TPS8 overexpressing lines and CRISPR/Cas9 mutants indicated that BnaC02.TPS8 enhanced the partitioning of photoassimilate into starch and sucrose, as opposed to glycolytic intermediates and organic acids, which might be associated with TPS activity. Furthermore, the overexpression of BnaC02.TPS8 not only increased seed yield but also enhanced seed oil accumulation and improved the oil fatty acid composition in B. napus under both high nitrogen (N) and low N conditions in the field. These results highlight the role of class II TPS in impacting photosynthesis and seed yield of B. napus, and BnaC02.TPS8 emerges as a promising target for improving B. napus seed yield.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus.

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

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.

Transient efflux inhibition improves plant regeneration by natural auxins.

Plant genome editing and propagation are important tools in crop breeding and production. Both rely heavily on the development of efficient in vitro plant regeneration systems. Two prominent regeneration systems that are widely employed in crop production are somatic embryogenesis (SE) and de novo shoot regeneration. In many of the protocols for SE or shoot regeneration, explants are treated with the synthetic auxin analog 2,4-dichlorophenoxyacetic acid (2,4-D), since natural auxins, such as indole-3-acetic acid (IAA) or 4-chloroindole-3-acetic acid (4-Cl-IAA), are less effective or even fail to induce regeneration. Based on previous reports that 2,4-D, compared to endogenous auxins, is not effectively exported from plant cells, we investigated whether efflux inhibition of endogenous auxins could convert these auxins into efficient inducers of SE in Arabidopsis immature zygotic embryos (IZEs). We show that natural auxins and synthetic analogs thereof become efficient inducers of SE when their efflux is transiently inhibited by co-application of the auxin transport inhibitor naphthylphthalamic acid (NPA). Moreover, IZEs of auxin efflux mutants pin2 or abcb1 abcb19 show enhanced SE efficiency when treated with IAA or efflux-inhibited IAA, confirming that auxin efflux reduces the efficiency of Arabidopsis SE. Importantly, in contrast to the 2,4-D system, where only 50-60% of the embryos converted to seedlings, all SEs induced by transport-inhibited natural auxins converted to seedlings. Efflux-inhibited IAA, like 2,4-D, also efficiently induced SE from carrot suspension cells, whereas IAA alone could not, and efflux-inhibited 4-Cl-IAA significantly improved de novo shoot regeneration in Brassica napus. Our data provides new insights into the action of 2,4-D as an efficient inducer of plant regeneration but also shows that replacing this synthetic auxin for efflux-inhibited natural auxin significantly improves different types of plant regeneration, leading to a more synchronized and homogenous development of the regenerated plants.

A novel three-layer module BoMYB1R1-BoMYB4b/BoMIEL1-BoDFR1 regulates anthocyanin accumulation in kale.

Anthocyanin is an important pigment responsible for plant coloration and beneficial to human health. Kale (Brassica oleracea var. acephala), a primary cool-season flowers and vegetables, is an ideal material to study anthocyanin biosynthesis and regulation mechanisms due to its anthocyanin-rich leaves. However, the underlying molecular mechanism of anthocyanin accumulation in kale remains poorly understood. Previously, we demonstrated that BoDFR1 is a key gene controlling anthocyanin biosynthesis in kale. Here, we discovered a 369-bp InDel variation in the BoDFR1 promoter between the two kale inbred lines with different pink coloration, which resulted in reduced transcriptional activity of the BoDFR1 gene in the light-pink line. With the 369-bp insertion as a bait, an R2R3-MYB repressor BoMYB4b was identified using the yeast one-hybrid screening. Knockdown of the BoMYB4b gene led to increased BoDFR1 expression and anthocyanin accumulation. An E3 ubiquitin ligase, BoMIEL1, was found to mediate the degradation of BoMYB4b, thereby promoting anthocyanin biosynthesis. Furthermore, the expression level of BoMYB4b was significantly reduced by light signals, which was attributed to the direct repression of the light-signaling factor BoMYB1R1 on the BoMYB4b promoter. Our study revealed that a novel regulatory module comprising BoMYB1R1, BoMIEL1, BoMYB4b, and BoDFR1 finely regulates anthocyanin accumulation in kale. The findings aim to establish a scientific foundation for genetic improvement of leaf color traits in kale, meanwhile, providing a reference for plant coloration studies.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds.

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U.

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

A GDSL motif-containing lipase modulates Sclerotinia sclerotiorum resistance in Brassica napus.

Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum (Lib.) De Bary is a devastating disease infecting hundreds of plant species. It also restricts the yield, quality, and safe production of rapeseed (Brassica napus) worldwide. However, the lack of resistance sources and genes to S. sclerotiorum has greatly restricted rapeseed SSR-resistance breeding. In this study, a previously identified GDSL motif-containing lipase gene, Brassica napus GDSL LIPASE-LIKE 1 (BnaC07.GLIP1), encoding a protein localized to the intercellular space, was characterized as functioning in plant immunity to S. sclerotiorum. The BnaC07.GLIP1 promoter is S. sclerotiorum-inducible and the expression of BnaC07.GLIP1 is substantially enhanced after S. sclerotiorum infection. Arabidopsis (Arabidopsis thaliana) heterologously expressing and rapeseed lines overexpressing BnaC07.GLIP1 showed enhanced resistance to S. sclerotiorum, whereas RNAi suppression and CRISPR/Cas9 knockout B. napus lines were hyper-susceptible to S. sclerotiorum. Moreover, BnaC07.GLIP1 affected the lipid composition and induced the production of phospholipid molecules, such as phosphatidylethanolamine, phosphatidylcholine and phosphatidic acid, which were correlated with decreased levels of reactive oxygen species (ROS) and enhanced expression of defense-related genes. A B. napus bZIP44 transcription factor specifically binds the CGTCA motif of the BnaC07.GLIP1 promoter to positively regulate its expression. BnbZIP44 responded to S. sclerotiorum infection, and its heterologous expression inhibited ROS accumulation, thereby enhancing S. sclerotiorum resistance in Arabidopsis. Thus, BnaC07.GLIP1 functions downstream of BnbZIP44 and is involved in S. sclerotiorum resistance by modulating the production of phospholipid molecules and ROS homeostasis in B. napus, providing insights into the potential roles and functional mechanisms of BnaC07.GLIP1 in plant immunity and for improving rapeseed SSR disease-resistance breeding.

Transcriptional regulation of transcription factor genes WRI1 and LAFL during Brassica napus seed development.

The WRINKLED1 (WRI1) and LAFL [LEAFY COTYLEDON1 (LEC1), ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2] transcription factors play essential roles in governing seed development and oil biosynthesis. To gain a comprehensive understanding of the transcriptional regulation of WRI1 and LAFL, we conducted genome-wide association studies for the expression profiles of WRI1 and LAFL in developing seeds at 20 and 40 days after flowering (DAF) using 302 rapeseed (Brassica napus) accessions. We identified a total of 237 expression quantitative trait nucleotides (eQTNs) and 51 expression QTN-by-environment interactions (eQEIs) associated with WRI1 and LAFL. Around these eQTNs and eQEIs, we pinpointed 41 and 8 candidate genes with known transcriptional regulations or protein interactions with their expression traits, respectively. Based on RNA-seq and ATAC-seq data, we employed the XGBoost and Basenji models which predicted 15 candidate genes potentially regulating the expression of WRI1 and LAFL. We further validated the predictions via tissue expression profile, haplotype analysis, and expression correlation analysis, and verified the transcriptional activation activity of BnaC03.MYB56 (R2R3-MYB transcription factor 56) on the expression of BnaA09.LEC1 by dual-luciferase reporter and yeast one-hybrid assays. BnaA10.AGL15 (AGAMOUS-LIKE 15), BnaC04.VAL1 (VIVIPAROUS1/ABSCISIC ACID INSENSITIVE3-LIKE 1), BnaC03.MYB56, and BnaA10.MYB56 were co-expressed with WRI1 and LAFL at 20 DAF in M35, a key module for seed development and oil biosynthesis. We further validated the positive regulation of MYB56 on seed oil accumulation using Arabidopsis (Arabidopsis thaliana) mutants. This study not only delivers a framework for future eQEI identification but also offers insights into the developmental regulation of seed oil accumulation.

The roles of cell wall polysaccharides in response to waterlogging stress in Brassica napus L. root.

BACKGROUND: Brassica napus L. (B. napus) is susceptible to waterlogging stress during different cultivation periods. Therefore, it is crucial to enhance the resistance to waterlogging stress to achieve a high and stable yield of B. napus. RESULTS: Here we observed significant differences in the responses of two B. napus varieties in root under waterlogging stress. The sensitive variety (23651) exhibited a more pronounced and rapid reduction in cell wall thickness and root integrity compared with the tolerant variety (Santana) under waterlogging stress. By module clustering analysis based on transcriptome data, we identified that cell wall polysaccharide metabolism responded to waterlogging stress in root. It was found that pectin content was significantly reduced in the sensitive variety compared with the tolerant variety. Furthermore, transcriptome analysis revealed that the expression of two homologous genes encoding polygalacturonase-inhibiting protein 2 (PGIP2), involved in polysaccharide metabolic pathways, was highly upregulated in root of the tolerant variety under waterlogging stress. BnaPGIP2s probably confer waterlogging resistance by inhibiting the activity of polygalacturonases (PGs), which in turn reduces the degradation of the pectin backbone polygalacturonic acid. CONCLUSIONS: Our findings demonstrate that cell wall polysaccharides in root plays a vital role in response to the waterlogging stress and provide a theoretical foundation for breeding waterlogging resistance in B. napus varieties.

Direct access to millions of mutations by whole genome sequencing of an oilseed rape mutant population.

Induced mutations are an essential source of genetic variation in plant breeding. Ethyl methanesulfonate (EMS) mutagenesis has been frequently applied, and mutants have been detected by phenotypic or genotypic screening of large populations. In the present study, a rapeseed M2 population was derived from M1 parent cultivar 'Express' treated with EMS. Whole genomes were sequenced from fourfold (4×) pools of 1988 M2 plants representing 497 M2 families. Detected mutations were not evenly distributed and displayed distinct patterns across the 19 chromosomes with lower mutation rates towards the ends. Mutation frequencies ranged from 32/Mb to 48/Mb. On average, 284 442 single nucleotide polymorphisms (SNPs) per M2 DNA pool were found resulting from EMS mutagenesis. 55% of the SNPs were C → T and G → A transitions, characteristic for EMS induced ('canonical') mutations, whereas the remaining SNPs were 'non-canonical' transitions (15%) or transversions (30%). Additionally, we detected 88 725 high confidence insertions and deletions per pool. On average, each M2 plant carried 39 120 canonical mutations, corresponding to a frequency of one mutation per 23.6 kb. Approximately 82% of such mutations were located either 5 kb upstream or downstream (56%) of gene coding regions or within intergenic regions (26%). The remaining 18% were located within regions coding for genes. All mutations detected by whole genome sequencing could be verified by comparison with known mutations. Furthermore, all sequences are accessible via the online tool 'EMSBrassica' (http://www.emsbrassica.plantbreeding.uni-kiel.de), which enables direct identification of mutations in any target sequence. The sequence resource described here will further add value for functional gene studies in rapeseed breeding.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae.

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5 million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

Fine mapping of meiotic crossovers in Brassica oleracea reveals patterns and variations depending on direction and combination of crosses.

Meiotic recombination is crucial for assuring proper segregation of parental chromosomes and generation of novel allelic combinations. As this process is tightly regulated, identifying factors influencing rate, and distribution of meiotic crossovers (COs) is of major importance, notably for plant breeding programs. However, high-resolution recombination maps are sparse in most crops including the Brassica genus and knowledge about intraspecific variation and sex differences is lacking. Here, we report fine-scale resolution recombination landscapes for 10 female and 10 male crosses in Brassica oleracea, by analyzing progenies of five large four-way-cross populations from two reciprocally crossed F1s per population. Parents are highly diverse inbred lines representing major crops, including broccoli, cauliflower, cabbage, kohlrabi, and kale. We produced approximately 4.56T Illumina data from 1248 progenies and identified 15 353 CO across the 10 reciprocal crosses, 51.13% of which being mapped to <10 kb. We revealed fairly similar Mb-scale recombination landscapes among all cross combinations and between the sexes, and provided evidence that these landscapes are largely independent of sequence divergence. We evidenced strong influence of gene density and large structural variations on CO formation in B. oleracea. Moreover, we found extensive variations in CO number depending on the direction and combination of the initial parents crossed with, for the first time, a striking interdependency between these factors. These data improve our current knowledge on meiotic recombination and are important for Brassica breeders.

Assassination tango: an NLR/NLR-ID immune receptors pair of rapeseed co-operates inside the nucleus to activate cell death.

Plant immunity largely relies on intracellular nucleotide-binding domain leucine-rich repeat (NLR) immune receptors. Some plant NLRs carry integrated domains (IDs) that mimic authentic pathogen effector targets. We report here the identification of a genetically linked NLR-ID/NLR pair: BnRPR1 and BnRPR2 in Brassica napus. The NLR-ID carries two ID fusions and the mode of action of the pair conforms to the proposed "integrated sensor/decoy" model. The two NLRs interact and the heterocomplex localizes in the plant-cell nucleus and nucleolus. However, the BnRPRs pair does not operate through a negative regulation as it was previously reported for other NLR-IDs. Cell death is induced only upon co-expression of the two proteins and is dependent on the helper genes, EDS1 and NRG1. The nuclear localization of both proteins seems to be essential for cell death activation, while the IDs of BnRPR1 are dispensable for this purpose. In summary, we describe a new pair of NLR-IDs with interesting features in relation to its regulation and the cell death activation.