Pathotype Characterization of Plasmodiophora brassicae by European Clubroot Differential and Williams Sets in China.

Clubroot disease caused by the soil-borne Plasmodiophora brassicae is devastating to Brassicaceae crops and spreading rapidly in China in recent years, resulting in great yield losses annually. Virulence of P. brassicae populations specializes and is in dynamic change in the fields. Information on the pathotypes and their distributions is crucial to control the clubroot disease. Presently, the pathotypes of P. brassicae prevalent in China, however, are not well determined. In this study, we used 16 Brassica hosts, including the European Clubroot Differential (ECD) and Williams sets, to designate the pathotypes of 33 P. brassicae populations from 13 provinces. The 33 P. brassicae populations could be divided into 26 pathotypes by the ECD set or seven pathotypes by the Williams set, revealing ECD16/15/31 and ECD16/31/31 or P4 and P2 as the predominant pathotypes. We found that the Brassica rapa differentials ECD01 to ECD04 showed stable and high levels of resistance to most pathotypes of P. brassicae in China, thereby providing valuable resources for clubroot-resistance breeding of Brassicaceae crops. The ECD set exhibited much higher discernibility and further divided the isolates that belonged to the P4 pathotype into 10 ECD pathotypes. Isolates of ECD16/23/31 and ECD16/15/31 were strongly virulent on Huashuang 5R, the first and widely used clubroot-resistant cultivar of oilseed rape in China. As we learn, 26 pathotypes are the most diverse populations of P. brassicae characterized until now in China. Our study provides new insights into virulence specialization of P. brassicae and their geographical distributions, contributing to exploitation of clubroot-resistant resources and the field layout of the present resistant Brassica crops in China.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L.

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.

The plant trans-Golgi network component ECHIDNA regulates defense, cell death, and endoplasmic reticulum stress.

The trans-Golgi network (TGN) acts as a central platform for sorting and secreting various cargoes to the cell surface, thus being essential for the full execution of plant immunity. However, the fine-tuned regulation of TGN components in plant defense and stress response has been not fully elucidated. Our study revealed that despite largely compromising penetration resistance, the loss-of-function mutation of the TGN component protein ECHIDNA (ECH) induced enhanced postinvasion resistance to powdery mildew in Arabidopsis thaliana. Genetic and transcriptome analyses and hormone profiling demonstrated that ECH loss resulted in salicylic acid (SA) hyperaccumulation via the ISOCHORISMATE SYNTHASE 1 biosynthesis pathway, thereby constitutively activating SA-dependent innate immunity that was largely responsible for the enhanced postinvasion resistance. Furthermore, the ech mutant displayed accelerated SA-independent spontaneous cell death and constitutive POWDERY MILDEW RESISTANCE 4-mediated callose depositions. In addition, ECH loss led to a chronically prolonged endoplasmic reticulum stress in the ech mutant. These results provide insights into understanding the role of TGN components in the regulation of plant immunity and stress responses.

Genome-wide association study reveals a GLYCOGEN SYNTHASE KINASE 3 gene regulating plant height in Brassica 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.

Genome-wide identification and comparative analysis of CLE family in rapeseed and its diploid progenitors.

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.

Differential alternative splicing genes and isoform co-expression networks of Brassica napus under multiple abiotic stresses.

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 reprogramming in fungal pathogen Sclerotinia sclerotiorum at different infection stages on Brassica napus.

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.

Comparative transcriptome analysis of compatible and incompatible Brassica napus-Xanthomonas campestris interactions.

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.

Genome-Wide Association Analysis Combined With Quantitative Trait Loci Mapping and Dynamic Transcriptome Unveil the Genetic Control of Seed Oil Content in Brassica napus L.

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.

Genome-Wide Identification and Analysis of Ariadne Gene Family Reveal Its Genetic Effects on Agronomic Traits of Brassica napus.

E3 ligases promote protein ubiquitination and degradation, which regulate every aspect of eukaryotic life. The Ariadne (ARI) proteins of RBR (ring between ring fingers) protein subfamily has been discovered as a group of potential E3 ubiquitin ligases. Only a few available research studies show their role in plant adaptations processes against the external environment. Presently, the functions of ARI proteins are largely unknown in plants. Therefore, in this study, we performed genome-wide analysis to identify the ARI gene family and explore their potential importance in B. napus. A total of 39 ARI genes were identified in the B. napus genome and were classified into three subfamilies (A, B and C) based on phylogenetic analysis. The protein-protein interaction networks and enrichment analysis indicated that BnARI genes could be involved in endoreduplication, DNA repair, proteasome assembly, ubiquitination, protein kinase activity and stress adaptation. The transcriptome data analysis in various tissues provided us an indication of some BnARI genes' functional importance in tissue development. We also identified potential BnARI genes that were significantly responsive towards the abiotic stresses. Furthermore, eight BnARI genes were identified as candidate genes for multiple agronomic traits through association mapping analysis in B. napus; among them, BnaA02g12100D, which is the ortholog of AtARI8, was significantly associated with ten agronomic traits. This study provided useful information on BnARI genes, which could aid targeted functional research and genetic improvement for breeding in B. napus.

Genome-Wide Characterization of Serine/Arginine-Rich Gene Family and Its Genetic Effects on Agronomic Traits of Brassica 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.

Genome-Wide Identification and Characterization of SET Domain Family Genes in Brassica napus L.

SET domain group encoding proteins function as histone lysine methyltransferases. These proteins are involved in various biological processes, including plant development and adaption to the environment by modifying the chromatin structures. So far, the SET domain genes (SDGs) have not been systematically investigated in Brassica napus (B. napus). In the current study, through genome-wide analysis, a total of 122 SDGs were identified in the B. napus genome. These BnSDGs were subdivided into seven (I-VII) classes based on phylogeny analysis, domain configurations, and motif distribution. Segmental duplication was involved in the evolution of this family, and the duplicated genes were under strong purifying selection. The promoter sequence of BnSDGs consisted of various growth, hormones, and stress-related cis-acting elements along with transcription factor binding sites (TFBSs) for 20 TF families in 59 of the 122 BnSDGs. The gene ontology (GO) analysis revealed that BnSDGs were closely associated with histone and non-histone methylation and metal binding capacity localized mostly in the nucleus. The in silico expression analysis at four developmental stages in leaf, stem root, floral organ, silique, and seed tissues showed a broad range of tissue and stage-specific expression pattern. The expression analysis under four abiotic stresses (dehydration, cold, ABA, and salinity) also provided evidence for the importance of BnSDGs in stress environments. Based on expression analysis, we performed reverse transcription-quantitative PCR for 15 target BnSDGs in eight tissues (young leaf, mature leaf, root, stem, carpel, stamen, sepal, and petals). Our results were in accordance with the in silico expression data, suggesting the importance of these genes in plant development. In conclusion, this study lays a foundation for future functional studies on SDGs in B. napus.

Integrating GWAS, linkage mapping and gene expression analyses reveal the genetic control of first branch height in Brassica napus L.

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.TCP1 M083 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.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica 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.

Mapping of a major QTL controlling plant height using a high-density genetic map and QTL-seq methods based on whole-genome resequencing in Brassica 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.

Characterization and Fine Mapping of a Yellow-Virescent Gene Regulating Chlorophyll Biosynthesis and Early Stage Chloroplast Development in Brassica napus.

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.