Search

Abstract As a Brassica crop, Brassica napus typically has single flowers that contain four petals. The double-flower phenotype of rapeseed has been a desirable trait in China because of its potential commercial value in ornamental tourism. However, few double-flowered germplasms have been documented in B. napus, and knowledge of the underlying genes is limited. Here, B. napus D376 was characterized as a double-flowered strain that presented an average of 10.92 ± 1.40 petals and other normal floral organs. F1, F2 and BC1 populations were constructed by crossing D376 with a single-flowered line reciprocally. Genetic analysis revealed that the double-flower trait was a recessive trait controlled by multiple genes. To identify the key genes controlling the double-flower trait, bulk segregant analysis sequencing (BSA-seq) and RNA-seq analyses were conducted on F2 individual bulks with opposite extreme phenotypes. Through BSA-seq, one candidate interval was mapped at the region of chromosome C05: 14.56–16.17 Mb. GO and KEGG enrichment analyses revealed that the DEGs were significantly enriched in carbohydrate metabolic processes, notably starch and sucrose metabolism. Interestingly, five and thirty-six DEGs associated with floral development were significantly up- and down-regulated, respectively, in the double-flowered plants. A combined analysis of BSA-seq and RNA-seq data revealed that five genes were candidates associated with the double flower trait, and BnaC05.ERS2 was the most promising gene. These findings provide novel insights into the breeding of double-flowered varieties and lay a theoretical foundation for unveiling the molecular mechanisms of floral development in B. napus.

The global expansion of rapeseed seed quality has been focused on maintaining glucosinolate (GSL) and erucic acid (EA) contents. However, the influence of seed GSL and EA contents on the germination process under drought stress remains poorly understood. Herein, 114 rapeseed accessions were divided into four groups based on GSL and EA contents to investigate their performance during seed imbibition under drought stress. Our results revealed significant variations in seed germination-related traits, particularly with higher GSL and EA, which exhibited higher germination % (G%) and lower mean germination time (MGT) under drought stress conditions. Moreover, osmoregulation, enzymatic system and hormonal regulation were improved in high GSL and high EA (HGHE) versus low GSL and low EA (LGLE) seeds, indicating the essential protective role of GSL and EA during the germination process in response to drought stress. The transcriptional regulation mechanism for coordinating GSL–EA-related pathways in response to drought stress during seed imbibition was found to involve the differential expression of sugar metabolism-, antioxidant-, and hormone-related genes with higher enrichment in HGHE compared to LGLE seeds. GO enrichment analysis showed higher variations in transcription regulator activity and DNA-binding transcription factors, as well as ATP and microtubule motor activity in GSL–EA-related pathways. Furthermore, KEGG analysis identified cellular processes, environmental information processing, and metabolism categories, with varied gene participation between GSL, EA and GSL–EA-related pathways. For further clarification, QY7 (LGLE) seeds were primed with different concentrations of GSL and EA under drought stress conditions. The results showed that 200 μmol/L of GSL and 400 μmol/L of EA significantly improved G%, MGT, and seedling fresh weight, besides regulating stress and fatty acid responsive genes during the seed germination process under drought stress conditions. Conclusively, exogenous application of GSL and EA is considered a promising method for enhancing the drought tolerance of LGLE seeds. Furthermore, the current investigation could provide a theoretical basis of GSL and EA roles and their underlying mechanisms in stress tolerance during the germination process.

Abstract Background Heterosis is an important biological phenomenon in which the hybrids exceed the parents in many traits. However, the molecular mechanism underlying seedling heterosis remains unclear. Results In the present study, we analyzed the leaf transcriptomes of strong hybrids (AM, HM) and weak hybrids (CM, HW) and their parents (A, C, H, M, and W) at two periods. Phenotypically, hybrids had obvious biomass heterosis at the seedling stage, with statistically significant differences between the strong and weak hybrids. The transcriptomic analysis demonstrated that the number of differentially expressed genes (DEGs) between parents was the highest. Further analysis showed that most DEGs were biased toward parental expression. The biological processes of the two periods were significantly enriched in the plant hormone signal transduction and photosynthetic pathways. In the plant hormone signaling pathway, DEG expression was high in hybrids, with expression differences between strong and weak hybrids. In addition, DEGs related to cell size were identified. Similar changes were observed during photosynthesis. The enhanced leaf area of hybrids generated an increase in photosynthetic products, which was consistent with the phenotype of the biomass. Weighted gene co-expression network analysis of different hybrids and parents revealed that hub genes in vigorous hybrid were mainly enriched in the plant hormone signal transduction and regulation of plant hormones. Conclusion Plant hormone signaling and photosynthesis pathways, as well as differential expression of plant cell size-related genes, jointly regulate the dynamic changes between strong and weak hybrids and the generation of seedling-stage heterosis. This study may elucidate the molecular mechanism underlying early biomass heterosis and help enhance canola yield.

Chloroplasts are essential for plant growth and development, as they play a key role in photosynthesis. The chloroplast biogenesis process is complex and its regulatory mechanism remains elusive. We characterized a spontaneous Brassica napus (rapeseed) mutant, ytg, that showed a delayed greening phenotype in all green organs and retarded growth. We identified BnaA02.YTG1 encoding a chloroplast-localized tetratricopeptide repeat protein widely expressed in rapeseed tissues. We speculated that the ytg phenotype was caused by the deletion of BnaA02.YTG1 based on sequence comparison of 4608 (with normal green leaves, isolated from the elite Chinese rapeseed cultivar ZS11) and ytg combined with transcriptome data and CRISPR/Cas9 gene editing results. The homologous gene (BnaC02.YTG1) restored the phenotype of the mutant. BnaA02.YTG1 interacted with MORF2, MORF8, and OZ1. RNA editing of the ndhD-2, ndhF-290, petL-5, and ndhG-50 plastid transcripts was affected in ytg. These findings suggested that BnaA02.YTG1 participates in RNA editing events. We predicted 29 RNA editing sites in the chloroplast of Brassica napus by comparison with the Arabidopsis chloroplast genome. We conclude that BnaA02.YTG1 affects the posttranscriptional regulation of plastid gene expression and suggest that a tetratricopeptide repeat protein is involved in the chloroplast RNA editing in rapeseed.

Abstract Salinity is a major constraint on crop growth and productivity, limiting sustainable agriculture in arid regions. Understanding the molecular mechanisms of salt-stress adaptation in canola is important to improve salt tolerance and promote its cultivation in saline lands. In this study, roots of control (no salt) and 200 mM NaCl-stressed canola seedlings were collected for RNA-Seq analysis and qRT-PCR validation. A total of 5385, 4268, and 7105 DEGs at the three time points of salt treatment compared to the control were identified, respectively. Several DEGs enriched in plant signal transduction pathways were highly expressed under salt stress, and these genes play an important role in signaling and scavenging of ROS in response to salt stress. Transcript expression in canola roots differed at different stages of salt stress, with the early-stages (2 h) of salt stress mainly related to oxidative stress response and sugar metabolism, while the late-stages (72 h) of salt stress mainly related to transmembrane movement, amino acid metabolism, glycerol metabolism and structural components of the cell wall. Several families of TFs that may be associated with salt tolerance were identified, including ERF, MYB, NAC, WRKY, and bHLH. These results provide a basis for further studies on the regulatory mechanisms of salt stress adaptation in canola.

Although utilization of heterosis has largely improved the yield of many crops worldwide, the underlying molecular mechanism of heterosis, particularly for allopolyploids, remains unclear. Here, we compared epigenome and transcriptome data of an elite hybrid and its parental lines in three assessed tissues (seedling, flower bud, and silique) to explore their contribution to heterosis in allopolyploid B. napus. Transcriptome analysis illustrated that a small proportion of non-additive genes in the hybrid compared with its parents, as well as parental expression level dominance, might have a significant effect on heterosis. We identified histone modification (H3K4me3 and H3K27me3) variation between the parents and hybrid, most of which resulted from the differences between parents. H3K4me3 variations were positively correlated with gene expression differences among the hybrid and its parents. Furthermore, H3K4me3 and H3K27me3 were rather stable in hybridization and were mainly inherited additively in the B. napus hybrid. Together, our data revealed that transcriptome reprogramming and histone modification remodeling in the hybrid could serve as valuable resources for better understanding heterosis in allopolyploid crops.

Inheritable albino mutants are excellent models for exploring the mechanism of chloroplast biogenesis and development. However, only a few non-lethal albino mutations have been reported to date in Brassica species. Here, we describe a resynthesized Brassica napus mutant, whose leaf, stem, and silique tissues showed an inheritable albino phenotype under field conditions after the bud stage but green phenotype in the greenhouse during the whole growing season, indicating that the albino phenotype depends on environmental conditions. Compared with the green leaves of the field-grown wild-type (GL) and greenhouse-grown mutant (WGL) plants, white leaves of the field-grown mutant (WL) showed significantly lower chlorophyll contents and structural defects in chloroplasts. Genetic analysis revealed that the albino phenotype of WL is recessive and is controlled by multiple genes. Bulk segregant analysis-sequencing (BSA-Seq) indicated that the candidate regions responsible for the albino phenotype spanned a total physical distance of approximately 49.68 Mb on chromosomes A03, A07, A08, C03, C04, C06, and C07. To gain insights into the molecular mechanisms that control chloroplast development in B. napus, we performed transcriptome (RNA-Seq) analysis of GL, WGL, and WL samples. GO and KEGG enrichment analyses suggested that differentially expressed genes (DEGs) associated with leaf color were significantly enriched in photosynthesis, ribosome biogenesis and chlorophyll metabolism. Further analysis indicated that DEGs involved in chloroplast development and chlorophyll metabolism were likely the main factors responsible for the albino phenotype in B. napus. A total of 59 DEGs were screened in the candidate regions, and four DEGs (BnaC03G0522600NO, BnaC07G0481600NO, BnaC07G0497800NO, and BnaA08G0016300NO) were identified as the most likely candidates responsible for the albino phenotype. Altogether, this study provides clues for elucidating the molecular mechanisms underlying chloroplast development in B. napus.

Abstract Background Cytoplasmic male sterility (CMS), which naturally exists in higher plants, is a useful mechanism for analyzing nuclear and mitochondrial genome functions and identifying the role of mitochondrial genes in the plant growth and development. Polima (pol) CMS is the most universally valued male sterility type in oil-seed rape. Previous studies have described the pol CMS restorer gene Rfp and the sterility-inducing gene orf224 in oil-seed rape, located in mitochondria. However, the mechanism of fertility restoration and infertility remains unknown. Moreover, it is still unknown how the fecundity restorer gene interferes with the sterility gene, provokes the sterility gene to lose its function, and leads to fertility restoration. Result In this study, we used multi-omics joint analysis to discover candidate genes that interact with the sterility gene orf224 and the restorer gene Rfp of pol CMS to provide theoretical support for the occurrence and restoration mechanisms of sterility. Via multi-omics analysis, we screened 24 differential genes encoding proteins related to RNA editing, respiratory electron transport chain, anther development, energy transport, tapetum development, and oxidative phosphorylation. Using a yeast two-hybrid assay, we obtained a total of seven Rfp interaction proteins, with orf224 protein covering five interaction proteins. Conclusions We propose that Rfp and its interacting protein cleave the transcript of atp6/orf224, causing the infertility gene to lose its function and restore fertility. When Rfp is not cleaved, orf224 poisons the tapetum cells and anther development-related proteins, resulting in pol CMS mitochondrial dysfunction and male infertility. The data from the joint analysis of multiple omics provided information on pol CMS’s potential molecular mechanism and will help breed B. napus hybrids.

Most plants are sensitive to salt-alkali stress, and the degree of tolerance to salt-alkali stress varies from different species and varieties. In order to explore the salt-alkali stress adaptability of Brassica napus, we collected the phenotypic data of 505 B. napus accessions at seedling and mature stages under control, low and high salt-alkali soil stress conditions in Inner Mongolia of China. Six resistant and 5 sensitive materials, respectively, have been identified both in Inner Mongolia and Xinjiang Uygur Autonomous Region of China. Genome-wide association studies (GWAS) for 15 absolute values and 10 tolerance coefficients (TCs) of growth and agronomic traits were applied to investigate the genetic basis of salt-alkali tolerance of B. napus. We finally mapped 9 significant QTLs related to salt-alkali stress response and predicted 20 candidate genes related to salt-alkali stress tolerance. Some important candidate genes, including BnABA4, BnBBX14, BnVTI12, BnPYL8, and BnCRR1, were identified by combining sequence variation annotation and expression differences. The identified valuable loci and germplasms could be useful for breeding salt-alkali-tolerant B.napus varieties. This study laid a foundation for understanding molecular mechanism of salt-alkali stress adaptation and provides rich genetic resources for the large-scale production of B. napus on salt-alkali land in the future.

Brassica napus L. (canola, oil seed rape) is one of the world’s most important oil seed crops. In the last four decades, the discovery of cytoplasmic male-sterility (CMS) systems and the restoration of fertility (Rf) genes in B. napus has improved the crop traits by heterosis. The homologs of Rf genes, known as the restoration of fertility-like (RFL) genes, have also gained importance because of their similarities with Rf genes. Such as a high non-synonymous/synonymous codon replacement ratio (dN/dS), autonomous gene duplications, and a possible engrossment in fertility restoration. B. napus contains 53 RFL genes on chromosomes A9 and C8. Our research aims to study the function of BnaRFL11 in fertility restoration using the CRISPR/Cas9 genome editing technique. A total of 88/108 (81.48%) T0 lines, and for T1, 110/145 (75%) lines carried T-DNA insertions. Stable mutations were detected in the T0 and T1 generations, with an average allelic mutation transmission rate of 81%. We used CRISPR-P software to detect off-target 50 plants sequenced from the T0 generation that showed no off-target mutation, signifying that if the designed sgRNA is specific for the target, the off-target effects are negligible. We also concluded that the mutagenic competence of the designed sgRNAs mediated by U6-26 and U6-29 ranged widely from 31% to 96%. The phenotypic analysis of bnarfl11 revealed defects in the floral structure, leaf size, branch number, and seed production. We discovered a significant difference between the sterile line and fertile line flower development after using a stereomicroscope and scanning electron microscope. The pollen visibility test showed that the pollen grain had utterly degenerated. The cytological observations of homozygous mutant plants showed an anther abortion stage similar to nap-CMS, with a Orf222, Orf139, Ap3, and nad5c gene upregulation. The bnarfl11 shows vegetative defects, including fewer branches and a reduced leaf size, suggesting that PPR-encoding genes are essential for the plants’ vegetative and reproductive growth. Our results demonstrated that BnaRFL11 has a possible role in fertility restoration. The current study’s findings suggest that CRISPR/Cas9 mutations may divulge the functions of genes in polyploid species and provide agronomically desirable traits through a targeted mutation.

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. The planting area and output of rapeseed are affected by the flowering time, which is a critical agronomic feature. COL9 controls growth and development in many different plant species as a member of the zinc finger transcription factor family. However, BnaCOL9 in rapeseed has not been documented. The aim of this study was to apply CRISPR/Cas9 technology to create an early-flowering germplasm resource to provide useful material for improving the early-maturing breeding of rapeseed. We identified four COL9 homologs in rapeseed that were distributed on chromosomes A05, C05, A03, and C03. We successfully created quadruple BnaCOL9 mutations in rapeseed using the CRISPR/Cas9 platform. The quadruple mutants of BnaCOL9 flowered earlier than the wild-type. On the other hand, the flowering time of the BnaCOL9 overexpression lines was delayed. An analysis of the expression patterns revealed that these genes were substantially expressed in the leaves and flowers. A subcellular localization experiment demonstrated that BnaCOL9 was in the nucleus. Furthermore, we discovered that two key flowering-related genes, BnaCO and BnaFT, were highly elevated in the BnaCOL9 mutants, but dramatically downregulated in the BnaCOL9 overexpression lines. Our findings demonstrate that BnaCOL9 is a significant flowering inhibitor in rapeseed and may be employed as a crucial gene for early-maturing breeding.

AbstractRapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for Bnapus

Despite the economic importance of the allotetraploid crop Brassica napus, our knowledge of its phylogenomic relationships, genetic structure, and diversification is limited. Here, the authors show its complex origin and diversification by analyzing transcriptome and organellar sequencing data.

Abstract Background Helitron is a rolling-circle DNA transposon; it plays an important role in plant evolution. However, Helitron distribution and contribution to evolution at the family level have not been previously investigated. Results We developed the software easy-to-annotate Helitron (EAHelitron), a Unix-like command line, and used it to identify Helitrons in a wide range of 53 plant genomes (including 13 Brassicaceae species). We determined Helitron density (abundance/Mb) and visualized and examined Helitron distribution patterns. We identified more than 104,653 Helitrons, including many new Helitrons not predicted by other software. Whole genome Helitron density is independent from genome size and shows stability at the species level. Using linear discriminant analysis, de novo genomes (next-generation sequencing) were successfully classified into Arabidopsis thaliana groups. For most Brassicaceae species, Helitron density negatively correlated with gene density, and Helitron distribution patterns were similar to those of A. thaliana. They preferentially inserted into sequence around the centromere and intergenic region. We also associated 13 Helitron polymorphism loci with flowering-time phenotypes in 18 A. thaliana ecotypes. Conclusion EAHelitron is a fast and efficient tool to identify new Helitrons. Whole genome Helitron density can be an informative character for plant classification. Helitron insertion polymorphism could be used in association analysis.

Abstract Background Alloplasmic lines, in which the nuclear genome is combined with wild cytoplasm, are often characterized by cytoplasmic male sterility (CMS), regardless of whether it was derived from sexual or somatic hybridization with wild relatives. In this study, we sequenced and analyzed the mitochondrial genomes of five such alloplasmic lines in Brassica juncea. Results The assembled and annotated mitochondrial genomes of the five alloplasmic lines were found to have virtually identical gene contents. They preserved most of the ancestral mitochondrial segments, and the same candidate male sterility gene (orf108) was found harbored in mitotype-specific sequences. We also detected promiscuous sequences of chloroplast origin that were conserved among plants of the Brassicaceae, and found the RNA editing profiles to vary across the five mitochondrial genomes. Conclusions On the basis of our characterization of the genetic nature of five alloplasmic mitochondrial genomes, we speculated that the putative candidate male sterility gene orf108 may not be responsible for the CMS observed in Brassica oxyrrhina and Diplotaxis catholica. Furthermore, we propose the potential coincidence of CMS in alloplasmic lines. Our findings lay the foundation for further elucidation of male sterility gene.

As a desirable agricultural trait, multi-inflorescence (MI) fulfills the requirement of mechanized harvesting and yield increase in rapeseed (Brassica napus L.). However, the genetic mechanism underlying the multi-inflorescence trait remain poorly understood. We previously identified a difference of one pair of dominant genes between the two mapping parental materials. In this study, phenotype and expression analysis indicated that the imbalance of the CLAVATA (CLV)-WUSCHEL (WUS) feedback loop may contribute to the abnormal development of the shoot apical meristem (SAM). BnaMI was fine-mapped to a 55 kb genomic region combining with genotype and phenotype of 5768 BCF1 individuals using a traditional mapping approach. Through comparative and expression analyses, combined with the annotation in Arabidopsis, five genes in this interval were identified as candidate genes. The present findings may provide assistance in functional analysis of the mechanism associated with multi-inflorescence and yield increase in rapeseed.

Increasing rapeseed yield has always been a primary goal of rapeseed research and breeding. However, flowering time is a prerequisite for stable rapeseed yield and determines its adaptability to ecological regions. MIKC-type MADS-box (MICK) genes are a class of transcription factors that are involved in various physiological and developmental processes in plants. To understand their role in floral transition-related pathways, a genome-wide screening was conducted with Brassica napus (B. napus), which revealed 172 members. Using previous data from a genome-wide association analysis of flowering traits, BnaSVP and BnaSEP1 were identified as candidate flowering genes. Therefore, we used the CRISPR/Cas9 system to verify the function of BnaSVP and BnaSEP1 in B. napus. T0 plants were edited efficiently at the BnaSVP and BnaSEP1 target sites to generate homozygous and heterozygous mutants with most mutations stably inherited by the next generation. Notably, the mutant only showed the early flowering phenotype when all homologous copies of BnaSVP were edited, indicating functional redundancy between homologous copies. However, no changes in flowering were observed in the BnaSEP1 mutant. Quantitative analysis of the pathway-related genes in the BnaSVP mutant revealed the upregulation of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and FLOWERING LOCUS T (FT) genes, which promoted early flowering in the mutant. In summary, our study created early flowering mutants, which provided valuable resources for early maturing breeding, and provided a new method for improving polyploid crops.

Drought is a major threat to plant growth and crop productivity. Reduced level of the gibberellin would result in increased drought tolerance, but the underlying mechanism is still unclear. In Brassica napus, there are four BnaRGA genes that code for DELLA proteins, negative regulators of GA signaling. Among them, expression of BnaA6.RGA was greatly induced by drought and abscisic acid (ABA). Previously, we created the gain-of-function mutant of BnaA6.RGA, bnaa6.rga-D, and the loss-of-function quadruple mutant, bnarga by CRISPR/Cas9, respectively. Here we show that bnaa6.rga-D displayed enhanced drought tolerance, and its stomatal closure was hypersensitive to ABA treatment. By contrast, bnarga displayed reduced drought tolerance and was less sensitive to ABA treatment, but there is no difference in drought tolerance between single BnaRGA mutant and WT, suggesting a functional redundancy between the BnaRGA genes in this process. Furthermore, we found that BnaRGAs were able to interact physically with BnaA10.ABF2, an essential transcription factor in ABA signaling. The BnaA10.ABF2-BnaA6.RGA protein complex greatly increased the expression level of the drought responsive gene BnaC9.RAB18. Taken together, this work highlighted the fundamental roles of DELLA proteins in drought tolerance in B. napus, and provide desirable germplasm for further breeding of drought tolerance in rapeseed.

Brassica juncea is an important dietary vegetable cultivated and consumed in China for its edible stalks and leaves. The purple leaf mustard, which is rich in anthocyanins, is eye-catching and delivers valuable nutrition. However, the molecular mechanism involved in anthocyanin biosynthesis has not been well studied in B. juncea. Here, histological and transcriptome analyses were used to characterize the purple leaf color and gene expression profiles. Free-hand section analysis showed that the anthocyanin was mainly accumulated in the adaxial epidermal leaf cells. The anthocyanin content in the purple leaves was significantly higher than that in the green leaves. To investigate the critical genes and pathways involved in anthocyanin biosynthesis and accumulation, the transcriptome analysis was used to identify the differentially expressed genes (DEGs) between the purple and green leaves from the backcrossed BC3 segregation population in B. juncea. A total of 2,286 different expressed genes were identified between the purple and green leaves. Among them, 1,593 DEGs were up-regulated and 693 DEGs were down-regulated. There were 213 differently expressed transcription factors among them. The MYB and bHLH transcription factors, which may regulate anthocyanin biosynthesis, were up-regulated in the purple leaves. Interestingly, most of the genes involved in plant–pathogen interaction pathway were also up-regulated in the purple leaves. The late biosynthetic genes involved in anthocyanin biosynthesis were highly up-regulated in the purple leaves of B. juncea. The up regulation of BjTT8 and BjMYC2 and anthocyanin biosynthetic genes (BjC4H, BjDFR, and BjANS) may activate the purple leaf formation in B. juncea. This study may help to understand the transcriptional regulation of anthocyanin biosynthesis in B. juncea.

Plant architecture involves important agronomic traits affecting crop yield, resistance to lodging, and fitness for mechanical harvesting in Brassica napus. Breeding high-yield varieties with plant architecture suitable for mechanical harvesting is the main goal of rapeseed breeders. Here, we report an accession of B. napus (4942C-5), which has a dwarf and compact plant architecture in contrast to cultivated varieties. A BC8 population was constructed by crossing a normal plant architecture line, 8008, with the recurrent parent 4942C-5. To investigate the molecular mechanisms underlying plant architecture, we performed phytohormone profiling, bulk segregant analysis sequencing (BSA-Seq), and RNA sequencing (RNA-Seq) in BC8 plants with contrasting plant architecture. Genetic analysis indicated the plant architecture traits of 4942C-5 were recessive traits controlled by multiple genes. The content of auxin (IAA), gibberellin (GA), and abscisic acid (ABA) differed significantly between plants with contrasting plant architecture in the BC8 population. Based on BSA-Seq analysis, we identified five candidate intervals on chromosome A01, namely those of 0 to 6.33 Mb, 6.45 to 6.48 Mb, 6.51 to 6.53 Mb, 6.77 to 6.79 Mb, and 7 to 7.01 Mb regions. The RNA-Seq analysis revealed a total of 4378 differentially expressed genes (DEGs), of which 2801 were up-regulated and 1577 were down-regulated. There, further analysis showed that genes involved in plant hormone biosynthesis and signal transduction, cell structure, and the phenylpropanoid pathway might play a pivotal role in the morphogenesis of plant architecture. Association analysis of BSA-Seq and RNA-Seq suggested that seven DEGs involved in plant hormone signal transduction and a WUSCHEL-related homeobox (WOX) gene (BnaA01g01910D) might be candidate genes responsible for the dwarf and compact phenotype in 4942C-5. These findings provide a foundation for elucidating the mechanisms underlying rapeseed plant architecture and should contribute to breed new varieties suitable for mechanization.

Brassica napus (oilseed rape) is one of the most important oil crops worldwide, but its growth is seriously threatened by Sclerotinia sclerotiorum. The mechanism of oilseed rape response to this pathogen has rarely been studied. Here, it was identified that BnaA03.MKK5 whose expression was induced by S. sclerotiorum infection was involved in plant immunity. BnaA03.MKK5 overexpression lines exhibited decreased disease symptoms compared to wild-type plants, accompanied by the increased expression of camalexin-biosynthesis-related genes, including BnPAD3 and BnCYP71A13. In addition, two copies of BnMPK3 (BnA06.MPK3 and BnC03.MPK3) were induced by Sclerotinia incubation, and BnaA03.MKK5 interacted with BnaA06.MPK3/BnaC03.MPK3 in yeast. These interactions were confirmed using in vivo co-immunoprecipitation assays. In vitro phosphorylation assays showed that BnaA06.MPK3 and BnaC03.MPK3 were the direct phosphorylation substrates of BnaA03.MKK5. The transgenic oilseed rape plants including BnaA06.MPK3 and BnaC03.MPK3 overexpression lines and BnMPK3 gene editing lines mediated by CRISPR/Cas9 were generated; the results of the genetic transformation of BnaA06.MPK3/BnaC03.MPK3 indicate that BnMPK3 also has a positive role in Sclerotinia resistance. This study provides information about the potential mechanism of B. napus defense against S. Sclerotiorum mediated by a detailed BnaA03.MKK5-BnaA06.MPK3/BnaC03.MPK3 module.

Soil salinity is one of the most significant abiotic stresses affecting crop yield around the world. To explore the molecular mechanism of salt tolerance in rapeseed (Brassica napus L.), the transcriptome analysis and metabolomics analysis were used to dissect the differentially expressed genes and metabolites in two rapeseed varieties with significant differences in salt tolerance; one is an elite rapeseed cultivar, Huayouza 62. A total of 103 key differentially expressed metabolites (DEMs) and 53 key differentials expressed genes (DEGs) that might be related to salt stress were identified through metabolomics and transcriptomics analysis. GO and KEGG analysis revealed that the DEGs were mainly involved in ion transport, reactive oxygen scavenging, osmotic regulation substance synthesis, and macromolecular protein synthesis. The DEMs were involved in TCA cycle, proline metabolism, inositol metabolism, carbohydrate metabolic processes, and oxidation-reduction processes. In addition, overexpression of BnLTP3, which was one of the key DEGs, could increase tolerance to salt stress in Arabidopsis plants. This study reveals that the regulation mechanism of salt tolerance in rapeseed at the transcriptome and metabolism level and provides abundant data for further in-depth identification of essential salt tolerance genes.

Abstract Background Accumulation of storage compounds during seed development plays an important role in the life cycle of oilseed plants; these compounds provide carbon and energy resources to support the establishment of seedlings. Results In this study, we show that BnCIPK9 has a broad expression pattern in Brassica napus L. tissues and that wounding stress strongly induces its expression. The overexpression of BnCIPK9 during seed development reduced oil synthesis in transgenic B. napus compared to that observed in wild-type (WT) plants. Functional analysis revealed that seed oil content (OC) of complementation lines was similar to that of WT plants, whereas OC in Arabidopsis thaliana (L.) Heynh. Atcipk9 knockout mutants (cipk9) was higher than that of WT plants. Seedling of cipk9 mutants failed to establish roots on a sugar-free medium, but root establishment could be rescued by supplementation of sucrose or glucose. The phenotype of complementation transgenic lines was similar to that of WT plants when grown on sugar-free medium. Mutants, cipk9, cbl2, and cbl3 presented similar phenotypes, suggesting that CIPK9, CBL2, and CBL3 might work together and play similar roles in root establishment under sugar-free condition. Conclusion This study showed that BnCIPK9 and AtCIPK9 encode a protein kinase that is involved in sugar-related response and plays important roles in the regulation of energy reserves. Our results suggest that AtCIPK9 negatively regulates lipid accumulation and has a significant effect on early seedling establishment in A. thaliana. The functional characterization of CIPK9 provides insights into the regulation of OC, and might be used for improving OC in B. napus. We believe that our study makes a significant contribution to the literature because it provides information on how CIPKs coordinate stress regulation and energy signaling.

Most flowering plants have evolved a self-incompatibility (SI) system to maintain genetic diversity by preventing self-pollination. The Brassica species possesses sporophytic self-incompatibility (SSI), which is controlled by the pollen- and stigma-determinant factors SP11/SCR and SRK. However, the mysterious molecular mechanism of SI remains largely unknown. Here, a new class II S haplotype, named BrS-325, was identified in a pak choi line ‘325’, which was responsible for the completely self-compatible phenotype. To obtain the entire S locus sequences, a complete pak choi genome was gained through Nanopore sequencing and de novo assembly, which provided a good reference genome for breeding and molecular research in B. rapa. S locus comparative analysis showed that the closest relatives to BrS-325 was BrS-60, and high sequence polymorphism existed in the S locus. Meanwhile, two duplicated SRKs (BrSRK-325a and BrSRK-325b) were distributed in the BrS-325 locus with opposite transcription directions. BrSRK-325b and BrSCR-325 were expressed normally at the transcriptional level. The multiple sequence alignment of SCRs and SRKs in class II S haplotypes showed that a number of amino acid variations were present in the contact regions (CR II and CR III) of BrSCR-325 and the hypervariable regions (HV I and HV II) of BrSRK-325s, which may influence the binding and interaction between the ligand and the receptor. Thus, these results suggested that amino acid variations in contact sites may lead to the SI destruction of a new class II S haplotype BrS-325 in B. rapa. The complete SC phenotype of ‘325’ showed the potential for practical breeding application value in B. rapa.

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

Self-incompatibility (SI) is a pollen-stigma recognition system controlled by a single and highly polymorphic genetic locus known as the S-locus. The S-locus exists in all Brassica napus (B. napus, AACC), but natural B. napus accessions are self-compatible. About 100 and 50 S haplotypes exist in Brassica rapa (AA) and Brassica oleracea (CC), respectively. However, S haplotypes have not been detected in B. napus populations. In this study, we detected the S haplotype distribution in B. napus and ascertained the function of a common S haplotype BnS-6 through genetic transformation. BnS-1/BnS-6 and BnS-7/BnS-6 were the main S haplotypes in 523 B. napus cultivars and inbred lines. The expression of SRK in different S haplotypes was normal (the expression of SCR in the A subgenome affected the SI phenotype) while the expression of BnSCR-6 in the C subgenome had no correlation with the SI phenotype in B. napus. The BnSCR-6 protein in BnSCR-6 overexpressed lines was functional, but the self-compatibility of overexpressed lines did not change. The low expression of BnSCR-6 could be a reason for the inactivation of BnS-6 in the SI response of B. napus. This study lays a foundation for research on the self-compatibility mechanism and the SI-related breeding in B. napus.

Despite heterosis contributing to genetic improvements in crops, root growth heterosis in rapeseed plants is poorly understood at the molecular level. The current study was performed to discover key differentially expressed genes (DEGs) related to heterosis in two hybrids with contrasting root growth performance (FO; high hybrid and FV; low hybrid) based on analysis of the root heterosis effect. Based on comparative transcriptomic analysis, we believe that the overdominance at the gene expression level plays a critical role in hybrid roots’ early biomass heterosis. Our findings imply that a considerable increase in up-regulation of gene expression underpins heterosis. In the FO hybrid, high expression of DEGs overdominant in the starch/sucrose and galactose metabolic pathways revealed a link between hybrid vigor and root growth. DEGs linked to auxin, cytokinin, brassinosteroids, ethylene, and abscisic acid were also specified, showing that these hormones may enhance mechanisms of root growth and the development in the FO hybrid. Moreover, transcription factors such as MYB, ERF, bHLH, NAC, bZIP, and WRKY are thought to control downstream genes involved in root growth. Overall, this is the first study to provide a better understanding related to the regulation of the molecular mechanism of heterosis, which assists in rapeseed growth and yield improvement.

Timely flowering is important for seed formation and maximization of rapeseed (Brassica napus) yield. Here, we performed flowering-time quantitative trait loci (QTL) mapping using a double haploid (DH) population grown in three environments to study the genetic architecture. Brassica 60 K Illumina Infinium™ single nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers were used for genotyping of the DH population, and a high-density genetic linkage map was constructed. QTL analysis of flowering time from the three environments revealed five consensus QTLs, including two major QTLs. A major QTL located on chromosome A03 was detected specifically in the semi-winter rapeseed growing region, and the one on chromosome C08 was detected in all environments. Ribonucleic acid sequencing (RNA-seq) was performed on the parents’ leaves at seven time-points in a day to determine differentially expressed genes (DEGs). The biological processes and pathways with significant enrichment of DEGs were obtained. The DEGs in the QTL intervals were analyzed, and four flowering time-related candidate genes were found. These results lay a foundation for the genetic regulation of rapeseed flowering time and create a rapeseed gene expression library for seven time-points in a day.

Abstract As a desirable agricultural trait, multilocular trait of rapeseed (Brassica rapa; Brassica napus; Brassica juncea), always represents higher yield per plant compared with bilocular plants. We previously isolated a trilocular gene locus, Bjmc1, and identified a set of molecular markers linked to the trilocular gene. With a map-based cloning, we identified that the BjMc1 was located in B genome of Brassica juncea, and it was a CLAVATA1 (CLV1) gene homologue. The insertion of a copia-LTR retrotransposable element 1 (RTE1) into the coding region of BjMc1 interrupted its transcription in rapeseed, leading to the trilocular phenotype. Phylogenetic analysis showed that Mc1 genes were conserved and widespread in land plants. Two amino acid sites had undergone positive selection in the ancestor of Mc1 genes, and then purifying selection was the dominant force after the divergence of dicots and monocots from their common ancestor in the evolutionary process, indicating that Mc1 genes are conserved in modern land plants. Our results provided new insights in molecular regulatory mechanism of multilocularity in rapeseed, and better understanding of molecular mechanism in crop yield improvement.

Monoacylglycerol lipase (MAGL) hydrolyzes monoacylglycerol, producing free fatty acid and glycerol. Although this enzyme has been shown to play important roles in mammal, its potential function in plants remains poorly understood. In a survey of the MAGL genes in Brassica napus, we found tapetal expression of BnaC.MAGL8.a, a homolog of AtMAGL8, results in male sterility in Arabidopsis thaliana. Retarded tapetal PCD and defective pollen wall were observed in the transgenic plants. The tapetal cells became vacuolated at stage 9, and then degenerated at stage 11. Most microspores degenerated with the tapetal cells, and only few pollen grains with an irregular-shaped exine layer were produced in the transgenic plants. Transcriptome analysis identified 398 differentially expressed genes. Most of them are involved in pollen development and stress response. ABORTED MICROSPORES and its downstream pollen wall biosynthesis genes were down-regulated, but genes related with reactive oxygen species homeostasis and jasmonates signaling were up-regulated in the transgenic plants. These results suggest that expression of BnaC.MAGL8.a in tapetum invokes stress response and impairs pollen development. The apparent phenotypic similarity between atgpat1 mutant and BnA9::BnaC.MAGL8.a transgenic plants lead us to propose a role for monoacylglycerol (MAG) in pollen development in Arabidopsis. Our study provides insights on not only the biological function of plant MAGL genes but also the role of MAG in pollen development.

Photosystem II (PSII) is an important component of the chloroplast. The PSII repair cycle is crucial for the relief of photoinhibition and may be advantageous when improving stress resistance and photosynthetic efficiency. Lethal genes are widely used in the efficiency detection and method improvement of gene editing. In the present study, we identified the naturally occurring lethal mutant 7-521Y with etiolated cotyledons in Brassica napus, controlled by double-recessive genes (named cyd1 and cyd2). By combining whole-genome resequencing and map-based cloning, CYD1 was fine-mapped to a 29 kb genomic region using 15,167 etiolated individuals. Through cosegregation analysis and functional verification of the transgene, BnaC06.FtsH1 was determined to be the target gene; it encodes an filamentation temperature sensitive protein H 1 (FtsH1) hydrolase that degrades damaged PSII D1 in Arabidopsis thaliana. The expression of BnaC06.FtsH1 was high in the cotyledons, leaves, and flowers of B. napus, and localized in the chloroplasts. In addition, the expression of EngA (upstream regulation gene of FtsH) increased and D1 decreased in 7-521Y. Double mutants of FtsH1 and FtsH5 were lethal in A. thaliana. Through phylogenetic analysis, the loss of FtsH5 was identified in Brassica, and the remaining FtsH1 was required for PSII repair cycle. CYD2 may be a homologous gene of FtsH1 on chromosome A07 of B. napus. Our study provides new insights into lethal mutants, the findings may help improve the efficiency of the PSII repair cycle and biomass accumulation in oilseed rape.

Cadmium is a potentially toxic heavy metal to human health. Rapeseed (Brassica napus L.), a vegetable and oilseed crop, might also be a Cd hyperaccumulator, but there is little information on this trait in rapeseed. We evaluated Cd accumulation in different oilseed accessions and employed a genome-wide association study to identify quantitative trait loci (QTLs) related to Cd accumulation. A total of 419 B. napus accessions and inbred lines were genotyped with a 60K Illumina Infinium SNP array of Brassica. Wide genotypic variations in Cd concentration and translocation were found. Twenty-five QTLs integrated with 98 single-nucleotide polymorphisms (SNPs) located at 15 chromosomes were associated with Cd accumulation traits. These QTLs explained 3.49–7.57% of the phenotypic variation observed. Thirty-two candidate genes were identified in these genomic regions, and they were 0.33–497.97 kb away from the SNPs. We found orthologs of Arabidopsis thaliana located near the significant SNPs on the B. napus genome, including NRAMP6 (natural resistance-associated macrophage protein 6), IRT1 (iron-regulated transporter 1), CAD1 (cadmium-sensitive 1), and PCS2 (phytochelatin synthase 2). Of them, four candidate genes were verified by qRT-PCR, the expression levels of which were significantly higher after exposure to Cd than in the controls. Our results might facilitate the study of the genetic basis of Cd accumulation and the cloning of candidate Cd accumulation genes, which could be used to help reduce Cd levels in edible plant parts and/or create more efficient hyperaccumulators.

Small peptides secreted to the extracellular matrix control many aspects of the plant’s physiological activities which were identified in Arabidopsis thaliana, called ATSPs. Here, we isolated and characterized the small peptide gene Bna.SP6 from Brassica napus. The BnaC.SP6 promoter was cloned and identified. Promoter deletion analysis suggested that the -447 to -375 and -210 to -135 regions are crucial for the silique septum and pollen expression of BnaC.SP6, respectively. Furthermore, the minimal promoter region of p158 (-210 to -52) was sufficient for driving gene expression specifically in pollen and highly conserved in Brassica species. In addition, BnaA.bZIP1 was predominantly expressed in anthers where BnaC.SP6 was also expressed, and was localized to the nuclei. BnaA.bZIP1 possessed transcriptional activation activity in yeast and protoplast system. It could specifically bind to the C-box in p158 in vitro, and negatively regulate p158 activity in vivo. BnaA.bZIP1 functions as a transcriptional repressor of BnaC.SP6 in pollen activity. These results provide novel insight into the transcriptional regulation of BnaC.SP6 in pollen activity and the pollen/anther-specific promoter regions of BnaC.SP6 may have their potential agricultural application for new male sterility line generation.

Drought stress is one of the major abiotic factors affecting Brassica napus (B. napus) productivity. In order to identify genes of potential importance to drought stress and obtain a deeper understanding of the molecular mechanisms regarding the responses of B. napus to dehydration stress, we performed large-scale transcriptome sequencing of B. napus plants under dehydration stress using the Illumina sequencing technology. In this work, a relatively drought tolerant B. napus line, Q2, identified in our previous study, was used. Four cDNA libraries constructed from mRNAs of control and dehydration-treated root and leaf were sequenced by Illumina technology. A total of 6018 and 5377 differentially expressed genes (DEGs) were identified in root and leaf. In addition, 1745 genes exhibited a coordinated expression profile between the two tissues under drought stress, 1289 (approximately 74%) of which showed an inverse relationship, demonstrating different regulation patterns between the root and leaf. The gene ontology (GO) enrichment test indicated that up-regulated genes in root were mostly involved in “stimulus” “stress” biological process, and activated genes in leaf mainly functioned in “cell” “cell part” components. Furthermore, a comparative network related to plant hormone signal transduction and AREB/ABF, AP2/EREBP, NAC, WRKY and MYC/MYB transcription factors (TFs) provided a view of different stress tolerance mechanisms between root and leaf. Some of the DEGs identified may be candidates for future research aimed at detecting drought-responsive genes and will be useful for understanding the molecular mechanisms of drought tolerance in root and leaf of B. napus.

Yellow seeds are a favorable trait for Brassica crops breeding due to better quality than their black-seeded counterparts. Here, we compared the Brassica napus seed coat transcriptomes between yellow- and brown-seeded near-isogenic lines (Y-NIL and B-NIL) that were developed from the resynthesized yellow-seeded line No. 2127-17. A total of 4,974 differentially expressed genes (DEG) were identified during seed development, involving 3,128 up-regulated and 1,835 down-regulated genes in yellow seed coats. Phenylpropanoid and flavonoid biosynthesis pathways were enriched in down-regulated genes, whereas the top two pathways for up-regulated genes were plant–pathogen interaction and plant hormone signal transduction. Twelve biosynthetic genes and three regulatory genes involved in the flavonoid pathway exhibited similar expression patterns in seed coats during seed development, of which the down-regulation mainly contributed to the reduction of proanthocyanidins (PAs) in yellow seed coats, indicating that these genes associated with PA biosynthesis may be regulated by an unreported common regulator, possibly corresponding to the candidate for the dominant black-seeded gene D in the NILs. Three transcription factor (TF) genes, including one bHLH gene and two MYB-related genes that are located within the previous seed coat color quantitative trait locus (QTL) region on chromosome A09, also showed similar developmental expression patterns to the key PA biosynthetic genes and they might thus potentially involved participate in flavonoid biosynthesis regulation. Our study identified novel potential TFs involved in PAs accumulation and will provide pivotal information for identifying the candidate genes for seed coat color in B. napus.

Brassica species exhibit both compatible and incompatible pollen-stigma interactions, however, the underlying molecular mechanisms remain largely unknown. Here, RNA-seq technology was applied in a comprehensive time-course experiment (2, 5, 10, 20, and 30 min) to explore gene expression during compatible/incompatible pollen-stigma interactions in stigma. Moderate changes of gene expression were observed both in compatible pollination (PC) and incompatible pollination (PI) within 10 min, whereas drastic changes showed up by 30 min, especially in PI. Stage specific DEGs [Differentially Expressed Gene(s)] were identified, and signaling pathways such as stress response, defense response, cell wall modification and others were found to be over-represented. In addition, enriched genes in all samples were analyzed as well, 293 most highly expressed genes were identified and annotated. Gene Ontology and metabolic pathway analysis revealed 10 most highly expressed genes and 37 activated metabolic pathways. According to the data, downstream components were activated in signaling pathways of both compatible and incompatible responses, and incompatible response had more complicated signal transduction networks. This study provides more detailed molecular information at different time points after compatible and incompatible pollination, deepening our knowledge about pollen-stigma interactions.

Soil salinity is a serious threat to agriculture sustainability worldwide. Salt tolerance at the seedling stage is crucial for plant establishment and high yield in saline soils; however, little information is available on rapeseed (Brassica napus L.) salt tolerance. We evaluated salt tolerance in different rapeseed accessions and conducted a genome-wide association study (GWAS) to identify salt tolerance-related quantitative trait loci (QTL). A natural population comprising 368 B. napus cultivars and inbred lines was genotyped with a Brassica 60K Illumina Infinium SNP array. The results revealed that 75 single-nucleotide polymorphisms (SNPs) distributed across 14 chromosomes were associated with four salt tolerance-related traits. These SNPs integrated into 25 QTLs that explained 4.21–9.23% of the phenotypic variation in the cultivars. Additionally, 38 possible candidate genes were identified in genomic regions associated with salt tolerance indices. These genes fell into several functional groups that are associated with plant salt tolerance, including transcription factors, aquaporins, transporters, and enzymes. Thus, salt tolerance in rapeseed involves complex molecular mechanisms. Our results provide valuable information for studying the genetic control of salt tolerance in B. napus seedlings and may facilitate marker-based breeding for rapeseed salt tolerance.

Self-incompatibility (SI) is a widespread mechanism in angiosperms that prevents inbreeding by rejecting self-pollen. However, the regulation of the SI response in Brassica napus is not well understood. Here, we report that the M-locus protein kinase (MLPK) BnaMLPKs, the functional homolog of BrMLPKs in Brassica rapa, controls SI in B. napus. We identified four paralogue MLPK genes in B. napus, including BnaA3.MLPK, BnaC3.MLPK, BnaA4.MLPK, and BnaC4.MLPK. Two transcripts of BnaA3.MLPK, BnaA3.MLPKf1 and BnaA3.MLPKf2, were generated by alternative splicing. Tissue expression pattern analysis demonstrated that BnaA3.MLPK, especially BnaA3.MLPKf2, is highly expressed in reproductive organs, particularly in stigmas. We subsequently created RNA-silencing lines and CRISPR/Cas9-induced quadruple mutants of BnaMLPKs in B. napus SI line S-70. Phenotypic analysis revealed that SI response is partially suppressed in RNA-silencing lines and is completely blocked in quadruple mutants. These results indicate the importance of BnaMLPKs in regulating the SI response of B. napus. We found that the expression of SI positive regulators S-locus receptor kinase (SRK) and Arm-Repeat Containing 1 (ARC1) are suppressed in bnmlpk mutant, whereas the self-compatibility (SC) element Glyoxalase I (GLO1) maintained a high expression level. Overall, our findings reveal a new regulatory mechanism of MLPK in the SI of B. napus.

The most important Brassica species, B. rapa, is naturally self-incompatible. Self-compatible mutants would be useful for dissecting the molecular mechanism of self-incompatibility (SI), a process that promotes outcrossing by recognizing and refusing self-pollens. The S haplotype in a new self-compatible B. rapa cultivar, Dahuangyoucai, was characterized for the first time in this study. Sequence analysis of the S-locus genes, SLG (S-locus glycoprotein), SRK (S-locus receptor kinase) and SCR (S-locus cysteine-rich protein) revealed that Dahuangyoucai contained S haplotype highly similar to S-f2, a non-functional class I S haplotype identified in another self-compatible B. rapa cultivar, Yellow Sarson. Mutations of MLPK (M-locus protein kinase) and non-transcription of the male determinant, SCR, were observed in this cultivar, which is similar to the situation reported in Yellow Sarson. With respect to the female determinant, SRK, no transcript was detected in Yellow Sarson but two fragments were detected in Dahuangyoucai. One fragment was highly similar to SRK-f2, but the other fragment was different from the signal factors previously identified in the SI reaction. The results suggest that Dahuangyoucai and Yellow Sarson have the same origin and a similar mechanism of self-compatibility, but diverge after mutations in SRK, SCR and MLPK. Further studying the self-compatibility of Dahuangyoucai might identify novel factors involved in the SI signalling cascade and provide new insights into the mechanisms of SI in Brassicaceae.

Photosynthesis competent autotrophy should be established during the postgerminative stage of plant growth. Among the multiple factors, light plays a decisive role in the switch from heterotrophic to autotrophic growth. Under dark condition, the rapeseed hypocotyl extends quickly with an apical hook, and the cotyledon is yellow and folded, and maintains high level of the isocitrate lyase (ICL). By contrast, in the light, the hypocotyl extends slowly, the cotyledon unfolds and turns green, the ICL content changes in parallel with the cotyledon greening. To reveal metabolic adaptations during the establishment of postgerminative autotrophy in rapeseed, we conducted comparative proteomic and metabolomic analyses of the cotyledons of seedlings grown under light versus dark conditions. Under both conditions, the increase of the protease, fatty acids β-oxidation and glyoxylate-cycle related proteins was accompanied with rapid degradation of the stored proteins and lipids with an accumulation of the amino acids, while light condition partially retarded these conversions. Light significantly induced the expression of chlorophyll-binding and photorespiration related proteins, resulting in an increase of the reducing-sugar. However, the levels of some chlide conversion, Calvin-cycle and photorespiration related proteins also accumulated in dark grown cotyledons, implying that the transition from heterotrophy to autotrophy is programmed in the seed rather than induced by light. Various anti-stress systems, e.g., redox related proteins, salicylic acid, proline and chaperones, were employed to release the oxidative stress, which was mainly derived from lipid oxidation or photorespiration, under both conditions. This study provides a comprehensive understanding of the differential molecular responses of rapeseed cotyledons to light and dark conditions, which will facilitate further study on the complex mechanism underlying the transition from heterotrophy to autotrophy.

Plant height is a key morphological trait of rapeseed. In this study, we measured plant height of a rapeseed population across six environments. This population contains 476 inbred lines representing the major Chinese rapeseed genepool and 44 lines from other countries. The 60K Brassica Infinium® SNP array was utilized to genotype the association panel. A genome-wide association study (GWAS) was performed via three methods, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Consequently, 68 loci were identified as significantly associated with plant height (P < 5.22 × 10-5), and more than 70% of the loci (48) overlapped the confidence intervals of reported QTLs from nine mapping populations. Moreover, 24 GWAS loci were detected with selective sweep signals, which reflected the signatures of historical semi-dwarf breeding. In the linkage disequilibrium (LD) decay range up- and downstream of 65 loci (r2 > 0.1), we found plausible candidates orthologous to the documented Arabidopsis genes involved in height regulation. One significant association found by GWAS colocalized with the established height locus BnRGA in rapeseed. Our results provide insights into the genetic basis of plant height in rapeseed and may facilitate marker-based breeding.

DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during plant development. However, few studies have determined the DNA methylation profiles of male-sterile rapeseed. Here, we conducted a global comparison of DNA methylation patterns between the rapeseed genic male sterile line 7365A and its near-isogenic fertile line 7365B by whole-genome bisulfite sequencing (WGBS). Profiling of the genome-wide DNA methylation showed that the methylation level in floral buds was lower than that in leaves and roots. Besides, a total of 410 differentially methylated region-associated genes (DMGs) were identified in 7365A relative to 7365B. Traditional bisulfite sequencing polymerase chain reaction (PCR) was performed to validate the WGBS data. Eleven DMGs were found to be involved in anther and pollen development, which were analyzed by quantitative PCR. In particular, Bnams4 was hypo-methylated in 7365A, and its expression was up-regulated, which might affect other DMGs and thus control the male sterility. This study provided genome-wide DNA methylation profiles of floral buds and important clues for revealing the molecular mechanism of genic male sterility in rapeseed.

Self-incompatibility (SI) in plants genetically prevents self-fertilization to promote outcrossing and genetic diversity. Its hybrids in Brassica have been widely cultivated due to the propagation of SI lines by spraying a salt solution. We demonstrated that suppression of Brassica napus SI from edible salt solution treatment was ascribed to sodium chloride and independent of S haplotypes, but it did not obviously change the expression of SI-related genes. Using the isobaric tags for relative and absolute quantitation (iTRAQ) technique, we identified 885 differentially accumulated proteins (DAPs) in Brassica napus stigmas of un-pollinated (UP), pollinated with compatible pollen (PC), pollinated with incompatible pollen (PI), and pollinated with incompatible pollen after edible salt solution treatment (NA). Of the 307 DAPs in NA/UP, 134 were unique and 94 were shared only with PC/UP. In PC and NA, some salt stress protein species, such as glyoxalase I, were induced, and these protein species were likely to participate in the self-compatibility (SC) pathway. Most of the identified protein species were related to metabolic pathways, biosynthesis of secondary metabolites, ribosome, and so on. A systematic analysis implied that salt treatment-overcoming SI in B. napus was likely conferred by at least five different physiological mechanisms: (i) the use of Ca2+ as signal molecule; (ii) loosening of the cell wall to allow pollen tube penetration; (iii) synthesis of compatibility factor protein species for pollen tube growth; (iv) depolymerization of microtubule networks to facilitate pollen tube movement; and (v) inhibition of protein degradation pathways to restrain the SI response.

Yellow seed is a desirable quality trait of the Brassica oilseed species. Previously, several seed coat color genes have been mapped in the Brassica species, but the molecular mechanism is still unknown. In the present investigation, map-based cloning method was used to identify a seed coat color gene, located on A9 in B. rapa. Blast analysis with the Arabidopsis genome showed that there were 22 Arabidopsis genes in this region including at4g09820 to at4g10620. Functional complementation test exhibited a phenotype reversion in the Arabidopsis thaliana tt8-1 mutant and yellow-seeded plant. These results suggested that the candidate gene was a homolog of TRANSPARENT TESTA8 (TT8) locus. BrTT8 regulated the accumulation of proanthocyanidins (PAs) in the seed coat. Sequence analysis of two alleles revealed a large insertion of a new class of transposable elements, Helitron in yellow sarson. In addition, no mRNA expression of BrTT8 was detected in the yellow-seeded line. It indicated that the natural transposon might have caused the loss in function of BrTT8. BrTT8 encodes a basic/helix-loop-helix (bHLH) protein that shares a high degree of similarity with other bHLH proteins in the Brassica. Further expression analysis also revealed that BrTT8 was involved in controlling the late biosynthetic genes (LBGs) of the flavonoid pathway. Our present findings provided with further studies could assist in understanding the molecular mechanism involved in seed coat color formation in Brassica species, which is an important oil yielding quality trait.