Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor‐bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R‐genes harbour more SNPs than variable genes. More nucleotide binding site‐leucine‐rich repeat (NBS‐LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics‐based improvement and improved disease resistance.     

Precise editing of CLAVATA genes in Brassica napus L. regulates multilocular silique development

Multilocular silique is a desirable agricultural trait with great potential for the development of high‐yield varieties of Brassica. To date, no spontaneous or induced multilocular mutants have been reported in Brassica napus, which likely reflects its allotetraploid nature and the extremely low probability of the simultaneous random mutagenesis of multiple gene copies with functional redundancy. Here, we present evidence for the efficient knockout of rapeseed homologues of CLAVATA3 (CLV3) for a secreted peptide and its related receptors CLV1 and CLV2 in the CLV signalling pathway using the CRISPR/Cas9 system and achieved stable transmission of the mutations across three generations. Each BnCLV gene has two copies located in two subgenomes. The multilocular phenotype can be recovered only in knockout mutations of both copies of each BnCLV gene, illustrating that the simultaneous alteration of multiple gene copies by CRISPR/Cas9 mutagenesis has great potential in generating agronomically important mutations in rapeseed. The mutagenesis efficiency varied widely from 0% to 48.65% in T₀ with different single‐guide RNAs (sgRNAs), indicating that the appropriate selection of the sgRNA is important for effectively generating indels in rapeseed. The double mutation of BnCLV3 produced more leaves and multilocular siliques with a significantly higher number of seeds per silique and a higher seed weight than the wild‐type and single mutant plants, potentially contributing to increased seed production. We also assessed the efficiency of the horizontal transfer of Cas9/gRNA cassettes by pollination. Our findings reveal the potential for plant breeding strategies to improve yield traits in currently cultivated rapeseed varieties.     

The high‐quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi‐winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi‐winter oilseed rape cultivar ‘ZS11’ and its comprehensive genomic comparison with the genomes of the winter‐type cultivar ‘Darmor‐bzh’ as well as two progenitors. The integrated BAC‐to‐BAC and whole‐genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high‐quality genome assembly of B. napus ‘ZS11’. Within a short evolutionary period (~6700 years ago), semi‐winter‐type ‘ZS11’ and the winter‐type ‘Darmor‐bzh’ maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to ‘Darmor‐bzh’, both two subgenomes of ‘ZS11’ are closely related to its progenitors, and the ‘ZS11’ genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi‐winter‐type ‘ZS11’ underwent potential genomic introgressions with B. rapa (A_r). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization‐responsive flowering time in ‘ZS11’ was first experienced HE, and then underwent genomic introgression event with A_r, which potentially has led to genetic differences in controlling vernalization in the semi‐winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi‐winter oilseed rape in Asia.     

Identification and comprehensive analysis of the CLV3/ESR‐related (CLE) gene family in Brassica napus L.

• The CLE (CLAVATA3/ESR) gene family, encoding a group of small secretory peptides, plays important roles in cell‐to‐cell communication, thereby controlling a broad spectrum of development processes. The CLE family has been systematically characterized in some plants, but not in Brassica napus.
• In the present study, 116 BnCLE genes were identified in the B. napus genome, including seven unannotated, six incorrectly predicted and five multi‐CLE domain‐encoding genes. These BnCLE members were separated into seven distinct groups based on phylogenetic analysis, which might facilitate the functional characterization of the peptides.
• Further characterization of CLE pre‐propeptides revealed 31 unique CLE peptides from 45 BnCLE genes, which may give rise to distinct roles of BnCLE and expansion of the gene family. The biological activity of these unique CLE dodecamer peptides was tested further through in vitro peptide assays. Variations in several important residues were identified as key contributors to the functional differentiation of BnCLE and expansion of the gene family in B. napus. Expression profile analysis helped to characterize possible functional redundancy and sub‐functionalization among the BnCLE members.
• This study presents a comprehensive overview of the CLE gene family in B. napus and provides a foundation for future evolutionary and functional studies.

     

Fatty acid distribution and lipid metabolism in developing seeds of laurate-producing rape (Brassica napus L.)

The fatty acid composition and content of membrane and storage lipids of two transgenic laurate-producing rape (Brassica napus L.) lines were monitored during seed development. The two lines, the medium-laurate (ML) line and the high-laurate (HL) line, accumulated 34 mol% and 55 mol% of laurate in their seed triacylglycerols, respectively. The diacylglycerols contained about 17 and 33 mol% of laurate in the ML- and HL-lines, respectively, from the mid-stage of seed development up to seed maturity. The ML-line showed a maximal relative laurate content in phosphatidylcholine (17 mol%) at the mid-stage of seed development whereafter the content decreased to 2.7 mol% with seed maturity. The laurate content in phosphatidylcholine was observed to remain high (26 mol%) in the HL-line from the mid-stage to the end of triacylglycerol deposition. Thereafter, the relative content decreased and reached 6.6 mol% in the mature seeds. There was an enhanced activity of lauroyl-phosphatidylcholine- metabolizing enzymes in the seed membranes from laurate-producing lines compared with control lines, which might explain the decrease seen in laurate content in phosphatidylcholine during seed maturation. A comparison of the laurate distribution in the lipids from developing laurate-producing rape seeds and developing seeds from three species naturally accumulating laurate at similar levels revealed differences in laurate metabolism compared with these species. The results suggest that phospholipids and triacylglycerols are synthesized from the same diacylglycerol pool in rape seeds and that rape lysophosphatidic acid acyltransferase and diacylglycerol acyltransferase do not have the same preference for laurate substrates as the corresponding enzymes in seed tissues naturally accumulating this acyl group. In addition, the mechanisms that specifically remove or exclude laurate from membrane lipids appear less effective in rape seed than in tissues naturally evolved to synthesize laurate-rich oils. Received: 23 December 1996 / Accepted: 16 April 1997     

Induction of telomere‐mediated chromosomal truncation and behavior of truncated chromosomes in Brassica napus

Engineered minichromosomes could be stably inherited and serve as a platform for simultaneously transferring and stably expressing multiple genes. Chromosomal truncation mediated by repeats of telomeric sequences is a promising approach for the generation of minichromosomes. In the present work, direct repetitive sequences of Arabidopsis telomere were used to study telomere‐mediated truncation of chromosomes in Brassica napus. Transgenes containing alien Arabidopsis telomere were successfully obtained, and Southern blotting and fluorescence in situ hybridization (FISH) results show that the transgenes resulted in successful chromosomal truncation in B. napus. In addition, truncated chromosomes were inherited at rates lower than that predicted by Mendelian rules. To determine the potential manipulations and applications of the engineered chromosomes, such as the stacking of multiple transgenes and the Cre/lox and FRT/FLP recombination systems, both amenable to genetic manipulations through site‐specific recombination in somatic cells, were tested for their ability to undergo recombination in B. napus. These results demonstrate that alien Arabidopsis telomere is able to mediate chromosomal truncation in B. napus. This technology would be feasible for chromosomal engineering and for studies on chromosome structure and function in B. napus.     

Development of iFOX‐hunting as a functional genomic tool and demonstration of its use to identify early senescence‐related genes in the polyploid Brassica napus

Functional genomic studies of many polyploid crops, including rapeseed (Brassica napus), are constrained by limited tool sets. Here we report development of a gain‐of‐function platform, termed ‘iFOX (inducible Full‐length cDNA OvereXpressor gene)‐Hunting’, for inducible expression of B. napus seed cDNAs in Arabidopsis. A Gateway‐compatible plant gene expression vector containing a methoxyfenozide‐inducible constitutive promoter for transgene expression was developed. This vector was used for cloning of random cDNAs from developing B. napus seeds and subsequent Agrobacterium‐mediated transformation of Arabidopsis. The inducible promoter of this vector enabled identification of genes upon induction that are otherwise lethal when constitutively overexpressed and to control developmental timing of transgene expression. Evaluation of a subset of the resulting ~6000 Arabidopsis transformants revealed a high percentage of lines with full‐length B. napus transgene insertions. Upon induction, numerous iFOX lines with visible phenotypes were identified, including one that displayed early leaf senescence. Phenotypic analysis of this line (rsl‐1327) after methoxyfenozide induction indicated high degree of leaf chlorosis. The integrated B. napuscDNA was identified as a homolog of an Arabidopsis acyl‐CoA binding protein (ACBP) gene designated BnACBP1‐like. The early senescence phenotype conferred by BnACBP1‐like was confirmed by constitutive expression of this gene in Arabidopsis and B. napus. Use of the inducible promoter in the iFOX line coupled with RNA‐Seq analyses allowed mechanistic clues and a working model for the phenotype associated with BnACBP1‐like expression. Our results demonstrate the utility of iFOX‐Hunting as a tool for gene discovery and functional characterization of Brassica napus genome.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Knockout of two BnaMAX1 homologs by CRISPR/Cas9‐targeted mutagenesis improves plant architecture and increases yield in rapeseed (Brassica napus L.)

Plant height and branch number are essential components of rapeseed plant architecture and are directly correlated with its yield. Presently, improvement of plant architecture is a major challenge in rapeseed breeding. In this study, we first verified that the two rapeseed BnaMAX1 genes had redundant functions resembling those of Arabidopsis MAX1, which regulates plant height and axillary bud outgrowth. Therefore, we designed two sgRNAs to edit these BnaMAX1 homologs using the CRISPR/Cas9 system. The T₀ plants were edited very efficiently (56.30%–67.38%) at the BnaMAX1 target sites resulting in homozygous, heterozygous, bi‐allelic and chimeric mutations. Transmission tests revealed that the mutations were passed on to the T₁ and T₂ progeny. We also obtained transgene‐free lines created by the CRISPR/Cas9 editing, and no mutations were detected in potential off‐target sites. Notably, simultaneous knockout of all four BnaMAX1 alleles resulted in semi‐dwarf and increased branching phenotypes with more siliques, contributing to increased yield per plant relative to wild type. Therefore, these semi‐dwarf and increased branching characteristics have the potential to help construct a rapeseed ideotype. Significantly, the editing resources obtained in our study provide desirable germplasm for further breeding of high yield in rapeseed.     

Melatonin mitigates cadmium and aluminium toxicity through modulation of antioxidant potential in Brassica napus L

• Melatonin has emerged as an essential molecule in plants, due to its role in defence against metal toxicity. Aluminium (Al) and cadmium (Cd) toxicity inhibit rapeseed seedling growth.
• In this study, we applied different doses of melatonin (50 and 100 µm ) to alleviate Al (25 µm ) and Cd (25 µm ) stress in rapeseed seedlings. Results show that Al and Cd caused toxicity in rapeseed seedling, as evidenced by a decrease in height, biomass and antioxidant enzyme activity.
• Melatonin increased the expression of melatonin biosynthesis‐related Brassica napus genes for caffeic acid O‐methyl transferase (BnCOMT) under Al and Cd stress. The genes BnCOMT‐1, BnCOMT‐5 and BnCOMT‐8 showed up‐regulated expression, while BnCOMT‐4 and BnCOMT‐6 were down‐regulated during incubation in water. Melatonin application increased the germination rate, shoot length, root length, fresh and dry weight of seedlings. Melatonin supplementation under Al and Cd stress increased superoxide dismutase, catalase, peroxidase, ascorbate peroxidase, proline, chlorophyll and anthocyanin content, as well as photosynthesis rate. Both Cd and Al treatments significantly increased hydrogen peroxide and malondialdehyde levels in rapeseed seedlings, which were strictly counterbalanced by melatonin. Analysis of Cd and Al in different subcellular compartments showed that melatonin enhanced cell wall and soluble fractions, but reduced the vacuolar and organelle fractions in Al‐ and Cd‐treated seedlings.
• These results suggest that melatonin‐induced improvements in antioxidant potential, biomass, photosynthesis rate and successive Cd and Al sequestration play a pivotal role in plant tolerance to Al and Cd stress. This mechanism may have potential implications in safe food production.

     

Down‐regulation of BnDA1, whose gene locus is associated with the seeds weight,improves the seeds weight and organ size in Brassica napus

Brassica napus L. is an important oil crop worldwide and is the main raw material for biofuel. Seed weight and seed size are the main contributors to seed yield. DA1 (DA means big in Chinese) is an ubiquitin receptor and negatively regulates seed size. Down‐regulation of AtDA1 in Arabidopsis leads to larger seeds and organs by increasing cell proliferation in integuments. In this study, BnDA1 was down‐regulated in B. napus by over expressed of AtDA1^R358K, which is a functional deficiency of DA1 with an arginine‐to‐lysine mutation at the 358th amino acid. The results showed that the biomass and size of the seeds, cotyledons, leaves, flowers and siliques of transgenic plants all increased significantly. In particular, the 1000 seed weight increased 21.23% and the seed yield per plant increased 13.22% in field condition. The transgenic plants had no negative traits related to yield. The candidate gene association analysis demonstrated that the BnDA1 locus was contributed to the seeds weight. Therefore, our study showed that regulation of DA1 in B. napus can increase the seed yield and biomass, and DA1 is a promising target for crop improvement.     

Thiol‐based redox proteins in abscisic acid and methyl jasmonate signaling in Brassica napus guard cells

Reversibly oxidized cysteine sulfhydryl groups serve as redox sensors or targets of redox sensing that are important in various physiological processes. However, little is known about redox‐sensitive proteins in guard cells and how they function in stomatal signaling. In this study, Brassica napus guard‐cell proteins altered by redox in response to abscisic acid (ABA) or methyl jasmonate (MeJA) were identified by complementary proteomics approaches, saturation differential in‐gel electrophoresis and isotope‐coded affinity tagging. In total, 65 and 118 potential redox‐responsive proteins were identified in ABA‐ and MeJA‐treated guard cells, respectively. All the proteins contain at least one cysteine, and over half of them are predicted to form intra‐molecular disulfide bonds. Most of the proteins fall into the functional groups of ‘energy’, ‘stress and defense’ and ‘metabolism’. Based on the peptide sequences identified by mass spectrometry, 30 proteins were common to ABA‐ and MeJA‐treated samples. A total of 44 cysteines were mapped in the identified proteins, and their levels of redox sensitivity were quantified. Two of the proteins, a sucrose non‐fermenting 1‐related protein kinase and an isopropylmalate dehydrogenase, were confirmed to be redox‐regulated and involved in stomatal movement. This study creates an inventory of potential redox switches, and highlights a protein redox regulatory mechanism in ABA and MeJA signal transduction in guard cells.     

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome‐wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue‐specific associated loci were related to total leaf and seed GSL content. Thirty‐six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double‐low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.