Search

Abstract Background Rapeseed (Brassica napus L.) is an important oil and industrial crop worldwide. Black rot caused by the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc) is an infectious vascular disease that leads to considerable yield losses in rapeseed. Resistance improvement through genetic breeding is an effective and sustainable approach to control black rot disease in B. napus. However, the molecular mechanisms underlying Brassica-Xcc interactions are not yet fully understood, especially regarding the impact of post-transcriptional gene regulation via alternative splicing (AS). Results In this study, we compared the AS landscapes of a susceptible parental line and two mutagenized B. napus lines with contrasting levels of black rot resistance. Different types of AS events were identified in these B. napus lines at three time points upon Xcc infection, among which intron retention was the most common AS type. A total of 1,932 genes was found to show differential AS patterns between different B. napus lines. Multiple defense-related differential alternative splicing (DAS) hub candidates were pinpointed through an isoform-based co-expression network analysis, including genes involved in pathogen recognition, defense signalling, transcriptional regulation, and oxidation reduction. Conclusion This study provides new insights into the potential effects of post-transcriptional regulation on immune responses in B. napus towards Xcc attack. These findings could be beneficial for the genetic improvement of B. napus to achieve durable black rot resistance in the future.

Oilseed rape (Brassica napus L.) is an allotetraploid (AACC, 2n ​= ​38) crop, valued for its edible oil and protein content. seed yield and nutritional composition of rapeseed are influenced by its yield and oil quality traits. However, the genetic basis of yield-related and oil-quality traits remain ambiguous. A panel of 266 diversified oilseed rape accessions was genotyped using 223 simple sequence repeat (SSR) markers covering all 19 chromosomes to identify significant markers associated with yield and quality traits. Twelve yield-related and six quality traits were investigated in two consecutive years (2014 and 2015), with three replications in two environments (Changshun, CS; and Qinghe, QH). Using the model GLM with population structure and kinship (Q ​+ ​K), a total of 25 significant SSR markers (P ​

Abstract Background Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus. Results Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops. Conclusion Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

To address the global demand for rapeseed while considering farmers’ profit, we face the challenges of making a quantum leap in seed yield and, at the same time, reducing yield loss due to biotic and abiotic stresses. We also face the challenge of efficiently applying new transformative biotechnology tools such as gene editing and breeding by genome design to increase rapeseed productivity and profitability. In this Perspective, we review advances in research on the physiological and genetic bases of both stress factors-affected yield stability and seed yield potential, focusing on source–sink relationships and allocation of photosynthetic assimilates to vegetative growth and seed development. We propose research directions and highlight the role of plant architecture in the relative contributions of the root system, leaves, and pods to seed yield. We call for de novo design of new rapeseed crops. We review trait variation in existing germplasm and biotechnologies available for crop design. Finally, we discuss opportunities to apply fundamental knowledge and key germplasm to rapeseed production and propose an ideotype for de novo design of future rapeseed cultivars.

Rapeseed (Brassica napus L.) is a crucial oil crop cultivated worldwide. First branch height, an essential component of rapeseed plant architecture, has an important effect on yield and mechanized harvesting; however, the underlying genetic mechanism remains unclear. In this study, based on the 60K single nucleotide polymorphism array and a recombinant inbred lines population derived from M083 and 888-5, a total of 19 QTLs were detected in five environments, distributed on linkage groups A02, A09, A10, C06, and C07, which explained phenotypic variation ranging from 4.87 to 29.87%. Furthermore, 26 significant SNPs were discovered on Chr.A02 by genome-wide association study in a diversity panel of 324 re-sequencing accessions. The major QTL of the first branch height trait was co-located on Chr.A02 by integrating linkage mapping and association mapping. Eleven candidate genes were screened via allelic variation analysis, inter-subgenomic synteny analysis, and differential expression of genes in parental shoot apical meristem tissues. Among these genes, BnaA02g13010D, which encodes a TCP transcription factor, was confirmed as the target gene according to gene function annotation, haplotype analysis, and full-length gene sequencing, which revealed that TATA insertion/deletion in the promoter region was closely linked to significantly phenotypic differences BnaA02.TCP1M083 overexpression resulted in decreased branch height and increased branch number in Arabidopsis. These results provide a genetic basis for first branch height and the ideal architecture of B. napus.

Rapeseed (Brassica napus) is an allotetraploid crop that is the main source of edible oils and feed proteins in the world. The ideal plant architecture breeding is a major objective of rapeseed breeding and determining the appropriate plant height is a key element of the ideal plant architecture. Therefore, this study aims to improve the understanding of the genetic controls underlying plant height. The plant heights of 230 rapeseed accessions collected worldwide were investigated in field experiments over two consecutive years in Wuhan, China. Whole-genome resequencing of these accessions yielded a total of 1,707,194 informative single nucleotide polymorphisms (SNPs) that were used for genome-wide association analysis (GWAS). GWAS and haplotype analysis showed that BnaA01g09530D, which encodes BRASSINOSTEROID-INSENSITIVE 2 and belongs to the GLYCOGEN SYNTHASE KINASE 3 (GSK3) family, was significantly associated with plant height in B. napus. Moreover, a total of 31 BnGSK3s with complete domains were identified from B. napus genome and clustered into four groups according to phylogenetic analysis, gene structure, and motif distribution. The expression patterns showed that BnGSK3s exhibited significant differences in 13 developmental tissues in B. napus, suggesting that BnGSK3s may be involved in tissue-specific development. Sixteen BnGSK3 genes were highly expressed the in shoot apical meristem, which may be related to plant height or architecture development. These results are important for providing new haplotypes of plant height in B. napus and for extending valuable genetic information for rapeseed genetic improvement of plant architecture.

Crop genomics and breeding CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) proteins belong to a small peptide family in plants. During plant development, CLE gene family members play a pivotal role in regulating cell-to-cell communication and stem cell maintenance. However, the evolutionary process and functional importance of CLEs are unclear in Brassicaceae. In this study, a total of 70 BnCLEs were identified in Brassica napus (2n = 4x = 38, AnCn): 32 from the An subgenome, 36 from the Cn subgenome, and 2 from the unanchored subgenome. Meanwhile, 29 BrCLE and 32 BoCLE genes were explored in Brassica rapa (2n = 2x = 20, Ar) and Brassica oleracea (2n = 2x = 18, Co). Phylogenetic analysis revealed that 163 CLEs derived from three Brassica species and Arabidopsis thaliana can be divided into seven subfamilies. Homology and synteny analyses indicated whole-genome triplication (WGT) and segmental duplication may be the major contributors to the expansion of CLE family. In addition, RNA-seq and qPCR analysis indicated that 19 and 16 BnCLEs were more highly expressed in immature seeds and roots than in other tissues. Some CLE gene pairs exhibited different expression patterns in the same tissue, which indicated possible functional divergence. Furthermore, genetic variations and regional association mapping analysis indicated that 12 BnCLEs were potential genes for regulating important agronomic traits. This study provided valuable information to understand the molecular evolution and biological function of CLEs in B. napus and its diploid progenitors, which will be helpful for genetic improvement of high-yield breeding in B. napus.

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.

Alternative splicing (AS) is an important post-transcriptional mechanism promoting the diversity of transcripts and proteins to regulate various life processes in eukaryotes. Sclerotinia stem rot is a major disease of Brassica napus caused by Sclerotinia sclerotiorum, which causes severe yield loss in B. napus production worldwide. Although many transcriptome studies have been carried out on the growth, development, and infection of S. sclerotiorum, the genome-wide AS events of S. sclerotiorum remain poorly understood, particularly at the infection stage. In this study, transcriptome sequencing was performed to systematically explore the genome-scale AS events of S. sclerotiorum at five important infection stages on a susceptible oilseed rape cultivar. A total of 130 genes were predicted to be involved in AS from the S. sclerotiorum genome, among which 98 genes were differentially expressed and may be responsible for AS reprogramming for its successful infection. In addition, 641 differential alternative splicing genes (DASGs) were identified during S. sclerotiorum infection, accounting for 5.76% of all annotated S. sclerotiorum genes, and 71 DASGs were commonly found at all the five infection stages. The most dominant AS type of S. sclerotiorum was found to be retained introns or alternative 3′ splice sites. Furthermore, the resultant AS isoforms of 21 DASGs became pseudogenes, and 60 DASGs encoded different putative proteins with different domains. More importantly, 16 DASGs of S. sclerotiorum were found to have signal peptides and possibly encode putative effectors to facilitate the infection of S. sclerotiorum. Finally, about 69.27% of DASGs were found to be non-differentially expressed genes, indicating that AS serves as another important way to regulate the infection of S. sclerotiorum on plants besides the gene expression level. Taken together, this study provides a genome-wide landscape for the AS of S. sclerotiorum during infection as well as an important resource for further elucidating the pathogenic mechanisms of S. sclerotiorum.

Black rot caused by the vascular pathogenic bacterium Xanthomonas campestris pv. campestris (Xcc) is widespread in Brassicaceae plants and an infectious disease that causes large yield losses in oil seed rape (Brassica napus L.). Improvement of resistance through breeding is a crucial strategy to prevent black rot disease in B. napus, but presently hampered by insufficient understanding of Xcc-Brassica interactions. This study compares two EMS-mutagenized B. napus lines that show contrasting resistance levels to their susceptible progenitor. Patterns of differential gene expression between these B. napus lines were evaluated at three time points post inoculation by comparative RNA-seq analysis. In line with the observed disease phenotypes, the susceptible line ZS9mXccS-1 displayed a steady amount of differentially expressed genes (DEGs) at different time points of infection, whereas the resistant line ZS9mXccR-1 displayed a gradual increase in DEGs throughout the course of infection. Weighted gene co-expression network analysis (WGCNA) pinpointed multiple defense-related hub genes with potential central roles in immunity, including the cell surface receptor genes CRK11 and BIR1, and the associated downstream regulatory genes WRKY11 and PBL30. KEGG analysis of DEGs belonging to two distinct co-expression modules revealed enriched pathways associated with defense, including Ca2+-signaling, receptor-mediated immunity, and phytohormone balance. Taken together, our comparative transcriptome analysis provides new avenues to unravel the mechanisms underlying black rot resistance in B. napus.

Rapeseed, an allotetraploid oil crop, provides vegetable oil for human consumption. The growing demand for oilseeds has necessitated the development of rapeseed varieties with improved quality. Therefore, a clear understanding of the genetic basis underlying the seed oil content (SOC) is required. In this study, a natural population comprising 204 diverse accessions and recombinant inbred lines (RILs) derived from Brassica napus and Sinapis alba via distant hybridization were collected for genome-wide association analysis (GWAS) and quantitative trait loci (QTL) mapping of the SOC trait, respectively. The variable coefficient of the RIL and natural populations ranged from 7.43 to 10.43% and 8.40 to 10.91%. Then, a high-density linkage map was constructed based on whole genome re-sequencing (WGS); the map harbored 2,799 bin markers and covered a total distance of 1,835.21 cM, with an average marker interval of 0.66 cM. The QTLs for SOC on chromosome A07 were stably detected in both single and multiple environments. Finally, a novel locus qA07.SOC was identified as the major QTL for SOC based on the GWAS and RIL populations. In addition, the RNA-seq results showed that photosynthesis, lipid biosynthesis proteins, fatty acid metabolism, and unsaturated fatty acid biosynthesis were significantly different between the developed seeds of the two parents of the RIL population. By comparing the variation information and expression levels of the syntenic genes within qA07.SOC and its syntenic genomic regions, as well as through haplotype analysis via GWAS, BnaA07.STR18, BnaA07.NRT1, and BnaA07g12880D were predicted as candidate genes in the qA07.SOC interval. These stable QTLs containing candidate genes and haplotypes can potentially provide a reliable basis for marker-assisted selection in B. napus breeding for SOC.

Chlorophyll biosynthesis and chloroplast development are crucial to photosynthesis and plant growth, but their regulatory mechanism remains elusive in many crop species. We isolated a Brassica napus yellow-virescent leaf (yvl) mutant, which exhibited yellow-younger-leaf and virescent-older-leaf with decreased chlorophyll accumulation and delayed chloroplast development. We mapped yvl locus to a 70-kb interval between molecular markers yvl-O10 and InDel-O6 on chromosome A03 in BC2F2 population using whole genome re-sequencing and bulked segregant analysis. The mutant had a ‘C’ to ‘T’ substitution in the coding sequence of BnaA03.CHLH, which encodes putative H subunit of Mg-protoporphyrin IX chelatase (CHLH). The mutation resulted in an imperfect protein structure and reduced activity of CHLH. It also hampered the plastid encoded RNA polymerase which transcribes regulatory genes of photosystem II and I. Consequently, the chlorophyll a/b and carotenoid contents were reduced and the chloroplast ultrastructure was degraded in yvl mutant. These results explain that a single nucleotide mutation in BnaA03.CHLH impairs PEP activity to disrupt chloroplast development and chlorophyll biosynthesis in B. napus.

Serine/arginine-rich (SR) proteins are indispensable factors for RNA splicing, and they play important roles in development and abiotic stress responses. However, little information on SR genes in Brassica napus is available. In this study, 59 SR genes were identified and classified into seven subfamilies: SR, SCL, RS2Z, RSZ, RS, SR45, and SC. In each subfamily, the genes showed relatively conserved structures and motifs, but displayed distinct expression patterns in different tissues and under abiotic stress, which might be caused by the varied cis-acting regulatory elements among them. Transcriptome datasets from Pacbio/Illumina platforms showed that alternative splicing of SR genes was widespread in B. napus and the majority of paralogous gene pairs displayed different splicing patterns. Protein-protein interaction analysis indicated that SR proteins were involved in the regulation of the whole lifecycle of mRNA, from synthesis to decay. Moreover, the association mapping analysis suggested that 12 SR genes were candidate genes for regulating specific agronomic traits, which indicated that SR genes could affect the development and hence influence the important agronomic traits of B. napus. In summary, this study provided elaborate information on SR genes in B. napus, which will aid further functional studies and genetic improvement of agronomic traits in B. napus.

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Plant height is a crucial element related to plant architecture that influences the seed yield of oilseed rape (Brassica napus L.). In this study, we isolated a natural B. napus mutant, namely a semi-dwarf mutant (sdw-e), which exhibits a 30% reduction in plant height compared with Zhongshuang 11-HP (ZS11-HP). Quantitative trait locus sequencing (QTL-seq) was conducted using two extreme DNA bulks in F2 populations in Wuchang-2017 derived from ZS11-HP ×sdw-e to identify QTLs associated with plant height. The result suggested that two QTL intervals were located on chromosome A10. The F2 population consisting of 200 individuals in Yangluo-2018 derived from ZS11-HP × sdw-e was used to construct a high-density linkage map using whole-genome resequencing. The high-density linkage map harbored 4323 bin markers and covered a total distance of 2026.52 cM with an average marker interval of 0.47 cM. The major QTL for plant height named qPHA10 was identified on linkage group A10 by interval mapping and composite interval mapping methods. The major QTL qPHA10 was highly consistent with the QTL-seq results. And then, we integrated the variation sites and expression levels of genes in the major QTL interval to predict the candidate genes. Thus, the identified QTL and candidate genes could be used in marker-assisted selection for B. napus breeding in the future. [ABSTRACT FROM AUTHOR]