Comparative Transcriptome Analysis Reveals Key Pathways and Hub Genes in Rapeseed During the Early Stage of Plasmodiophora brassicae Infection

Genome-wide identification of the ICS family genes and its role in resistance to Plasmodiophora brassicae in Brassica napus L

The isochorismate synthase (ICS) proteins are essential regulators of salicylic acid (SA) synthesis, which has been reported to regulate resistance to biotic and abiotic stresses in plants. Clubroot caused by Plasmodiophora brassicae is a common disease that threatens the yield and quality of Oilseed rape (Brassica napus L.). Exogenous application of salicylic acid reduced the incidence of clubroot in oilseed rape. However, the potential importance of the ICS genes family in B. napus and its diploid progenitors has been unclear. Here, we identified 16, 9, and 10 ICS genes in the allotetraploid B. napus, diploid ancestor Brassica rapa and Brassica oleracea, respectively. These ICS genes were classified into three subfamilies (I-III), and member of the same subfamilies showed relatively conserved gene structures, motifs, and protein domains. Furthermore, many hormone-response and stress-related promoter cis-acting elements were observed in the BnaICS genes. Exogenous application of SA delayed the growth of clubroot galls, and the expression of BnaICS genes was significantly different compared to the control groups. Protein-protein interaction analysis identified 58 proteins involved in the regulation of ICS in response to P. brassicae in B. napus. These results provide new clues for understanding the resistance mechanism to P. brassicae.

Plasmodiophora brassicae affects host gene expression by secreting the transcription factor-type effector PbZFE1

The protist pathogen Plasmodiophora brassicae hijacks the metabolism and development of host cruciferous plants and induces clubroot formation, but little is known about its regulatory mechanisms. Previously, the Pnit2int2 sequence, a sequence around the second intron of the nitrilase gene (BrNIT2) involved in auxin biosynthesis in Brassica rapa ssp. pekinensis, was identified as a specific promoter activated during clubroot formation. In this study, we hypothesized that analysis of the transcriptional regulation of Pnit2int2 could reveal how P. brassicae affects the host gene regulatory system during clubroot development. By yeast one-hybrid screening, the pathogen zinc finger protein PbZFE1 was identified to specifically bind to Pnit2int2. Specific binding of PbZFE1 to Pnit2int2 was also confirmed by electrophoretic mobility shift assay. The binding site of PbZFE1 is essential for promoter activity of Pnit2int2 in clubbed roots of transgenic Arabidopsis thaliana (Pnit2int2-2::GUS), indicating that PbZFE1 is secreted from P. brassicae and functions within plant cells. Ectopic expression of PbZEF1 in A. thaliana delayed growth and flowering time, suggesting that PbZFE1 has significant impacts on host development and metabolic systems. Thus, P. brassicae appears to secrete PbZFE1 into host cells as a transcription factor-type effector during pathogenesis.

In-Depth Understanding of the Genetic Control of Stripe Rust Resistance (Puccinia striiformis f. sp. tritici) Induced in Wheat (Triticum aestivum) by Trichoderma asperellum T34

Wheat stripe rust (caused by Puccinia striiformis f. tritici Erikss.) causes severe yield losses worldwide. Due to the continuous appearance of new stripe rust races, resistance has been broken in most of the highly resistant genotypes in Egypt and worldwide. Therefore, looking for new ways to resist such a severe disease is urgently needed. Trichoderma asperellum strain T34 has been known as an effective bioagent against many crop diseases. It exists naturally in Egyptian fields. Therefore, in our study, the effectiveness of strain T34 was tested as a bioagent against wheat stripe rust. For this purpose, 198 spring wheat genotypes were tested for their resistance against two different P. striiformis f. tritici populations collected from the Egyptian fields. The most highly aggressive P. striiformis f. tritici population was used to test the effectiveness of strain T34. Highly significant differences were found between strain T34 and stripe rust, suggesting the effectiveness of strain T34 in stripe rust resistance. A genome-wide association study identified 48 gene models controlling resistance under normal conditions and 46 gene models controlling strain T34-induced resistance. Of these gene models, only one common gene model was found, suggesting the presence of two different genetic systems controlling resistance under each condition. The pathways of the biological processes were investigated under both conditions. This study provided in-depth understanding of genetic control and, hence, will accelerate the future of wheat breeding programs for stripe rust resistance.

Metabonomic profiling of clubroot-susceptible and clubroot-resistant radish and the assessment of disease-resistant metabolites

Plasmodiophora brassicae causes a serious threat to cruciferous plants including radish (Raphanus sativus L.). Knowledge on the pathogenic regularity and molecular mechanism of P. brassicae and radish is limited, especially on the metabolism level. In the present study, clubroot-susceptible and clubroot-resistant cultivars were inoculated with P. brassicae Race 4, root hairs initial infection of resting spores (10⁷ CFU/mL) at 24 h post-inoculation and root galls symptom arising at cortex splitting stage were identified on both cultivars. Root samples of cortex splitting stage of two cultivars were collected and used for untargeted metabonomic analysis. We demonstrated changes in metabolite regulation and pathways during the cortex splitting stage of diseased roots between clubroot-susceptible and clubroot-resistant cultivars using untargeted metabonomic analysis. We identified a larger number of differentially regulated metabolites and heavier metabolite profile changes in the susceptible cultivar than in the resistant counterpart. The metabolites that were differentially regulated in both cultivars were mostly lipids and lipid-like molecules. Significantly regulated metabolites and pathways according to the P value and variable important in projection score were identified. Moreover, four compounds, including ethyl α-D-thioglucopyranoside, imipenem, ginsenoside Rg1, and 6-gingerol, were selected, and their anti-P. brassicae ability and effects on seedling growth were verified on the susceptible cultivar. Except for ethyl α-D-thioglucopyranoside, the remaining could inhibit clubroot development of varing degree. The use of 5 mg/L ginsenoside Rg1 + 5 mg/L 6-gingerol resulted in the lowest disease incidence and disease index among all treatments and enhanced seedling growth. The regulation of pathways or metabolites of carbapenem and ginsenoside was further explored. The results provide a preliminary understanding of the interaction between radish and P. brassicae at the metabolism level, as well as the development of measures for preventing clubroot.

Genome-wide survey of Calcium-Dependent Protein Kinases (CPKs) in five Brassica species and identification of CPKs induced by Plasmodiophora brassicae in B. rapa, B. oleracea, and B. napus

Calcium-dependent protein kinase (CPK) is a class of Ser/Thr protein kinase that exists in plants and some protozoa, possessing Ca²⁺ sensing functions and kinase activity. To better reveal the roles that Brassica CPKs played during plant response to stresses, five Brassica species, namely Brassica rapa (B. rapa), Brassica nigra (B. nigra), Brassica oleracea (B. oleracea), Brassica juncea (B. juncea), and Brassica napus (B. napus) were selected and analyzed. In total, 51 BraCPK, 56 BniCPK, 56 BolCPK, 88 BjuCPK, and 107 BnaCPK genes were identified genome wide and phylogenetics, chromosomal mapping, collinearity, promoter analysis, and biological stress analysis were conducted. The results showed that a typical CPK gene was constituted by a long exon and tandem short exons. They were unevenly distributed on most chromosomes except chromosome A08 in B. napus and B. rapa, and almost all CPK genes were located on regions of high gene density as non-tandem form. The promoter regions of BraCPKs, BolCPKs, and BnaCPKs possessed at least three types of cis-elements, among which the abscisic acid responsive-related accounted for the largest proportion. In the phylogenetic tree, CPKs were clustered into four primary groups, among which group I contained the most CPK genes while group IV contained the fewest. Some clades, like AT5G23580.1(CPK12) and AT2G31500.1 (CPK24) contained much more gene members than others, indicating a possibility that gene expansion occurred during evolution. Furthermore, 4 BraCPKs, 14 BolCPKs, and 31 BnaCPKs involved in the Plasmodiophora brassicae (P. brassicae) defense response in resistant (R) or susceptible (S) materials were derived from online databases, leading to the discovery that some R-specific induced CPKs, such as BnaC02g08720D, BnaA03g03800D, and BolC04g018270.2J.m1 might be ideal candidate genes for P. brassicae resistant research. Overall, these results provide valuable information for research on the function and evolution of CDK genes.

Integrated genetic mapping and transcriptome analysis reveal the BnaA03.IAA7 protein regulates plant architecture and gibberellin signaling in Brassica napus L

A novel mutation in the BnaA03.IAA7 protein reduces plant height and enhances gibberellin signaling in Brassica napus L. Rapeseed (Brassica napus) is an excellent and important source for vegetable oil production, but its production is severely affected by lodging. Lodging hinders mechanization and decreases yield, and an ideal solution is semidwarf breeding. Limited by germplasm resources, semidwarf breeding developed slowly in rapeseed. In the current study, a mutant called sdA03 was isolated from EMS-mutagenized lines of Zhongshuang 11 (ZS11). The inheritance analysis showed that phenotypes of sdA03 were controlled by a single semidominant gene. Genetic mapping, RNA-seq and candidate gene analysis identified BnaA03.IAA7 as a candidate gene, and a function test confirmed that the mutated BnaA03.iaa7 regulates plant architecture in a dose-dependent manner. Yeast two-hybrid and transient expression experiments illustrated the P87L substitution in the GWPPV/I degron motif of BnaA03.iaa7 impaired the interaction between BnaA03.IAA7 and TIR1 proteins, and BnaA03.iaa7 prevented ARF from activating the auxin signaling pathway.The gibberellin (GA) content was higher in sdA03 hypocotyls than in those of ZS11. Further expression analysis showed more active gibberellin signaling in hypocotyl and richer expression of GA synthetic genes in root and cotyledon of sdA03 seedlings. Finally, a marker was developed based on the SNP found in BnaA03.iaa7 and used in molecular breeding. The study enriched our understanding of the architectural regulation of rapeseed and provided germplasm resources for breeding.

Comparative transcriptome analysis reveals different defence responses during the early stage of wounding stress in Chi-Nan germplasm and ordinary Aquilaria sinensis

BACKGROUND

Agarwood is a valuable Chinese medicinal herb and spice that is produced from wounded Aquilaria spp., is widely used in Southeast Asia and is highly traded on the market. The lack of highly responsive Aquilaria lines has seriously restricted agarwood yield and the development of its industry. In this article, a comparative transcriptome analysis was carried out between ordinary A. sinensis and Chi-Nan germplasm, which is a kind of A. sinensis tree with high agarwood-producing capacity in response to wounding stress, to elucidate the molecular mechanism underlying wounding stress in different A. sinensis germplasm resources and to help identify and breed high agarwood-producing strains.

RESULTS

A total of 2427 and 1153 differentially expressed genes (DEGs) were detected in wounded ordinary A. sinensis and Chi-Nan germplasm compared with the control groups, respectively. KEGG enrichment analysis revealed that genes participating in starch metabolism, secondary metabolism and plant hormone signal transduction might play major roles in the early regulation of wound stress. 86 DEGs related to oxygen metabolism, JA pathway and sesquiterpene biosynthesis were identified. The majority of the expression of these genes was differentially induced between two germplasm resources under wounding stress. 13 candidate genes related to defence and sesquiterpene biosynthesis were obtained by WGCNA. Furthermore, the expression pattern of genes were verified by qRT-PCR. The candidate genes expression levels were higher in Chi-Nan germplasm than that in ordinary A. sinensis during early stage of wounding stress, which may play important roles in regulating high agarwood-producing capacity in Chi-Nan germplasm.

CONCLUSIONS

Compared with A. sinensis, Chi-Nan germplasm invoked different biological processes in response to wounding stress. The genes related to defence signals and sesquiterepene biosynthesis pathway were induced to expression differentially between two germplasm resources. A total of 13 candidate genes were identified, which may correlate with high agarwood-producting capacity in Chi-Nan germplasm during the early stage of wounding stress. These genes will contribute to the development of functional molecular markers and the rapid breeding highly of responsive Aquilaria lines.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

What Can We Learn from -Omics Approaches to Understand Clubroot Disease?

Clubroot is one of the most economically significant diseases worldwide. As a result, many investigations focus on both curing the disease and in-depth molecular studies. Although the first transcriptome dataset for the clubroot disease describing the clubroot disease was published in 2006, many different pathogen-host plant combinations have only recently been investigated and published. Articles presenting -omics data and the clubroot pathogen Plasmodiophora brassicae as well as different host plants were analyzed to summarize the findings in the richness of these datasets. Although genome data for the protist have only recently become available, many effector candidates have been identified, but their functional characterization is incomplete. A better understanding of the life cycle is clearly required to comprehend its function. While only a few proteome studies and metabolome analyses were performed, the majority of studies used microarrays and RNAseq approaches to study transcriptomes. Metabolites, comprising chemical groups like hormones were generally studied in a more targeted manner. Furthermore, functional approaches based on such datasets have been carried out employing mutants, transgenic lines, or ecotypes/cultivars of either Arabidopsis thaliana or other economically important host plants of the Brassica family. This has led to new discoveries of potential genes involved in disease development or in (partial) resistance or tolerance to P. brassicae. The overall contribution of individual experimental setups to a larger picture will be discussed in this review.

Transcriptomic Analysis Reveals the Regulatory Networks and Hub Genes Controlling the Unsaturated Fatty Acid Contents of Developing Seed in Soybean

Soybean [Glycine max (L.) Merr.] is one of the most important crops, which produces about 25% of the world's edible oil. The nutritional value of soybean oil depends mostly on the relative contents of three unsaturated fatty acids (UFAs), i.e., oleic acid, linoleic acid (LA), and linolenic acid. However, the biosynthetic mechanism of UFAs remains largely unknown, and there are few studies on RNA-seq analysis of developing seeds. To identify the candidate genes and related pathways involved in the regulation of UFA contents during seed development in soybean, two soybean lines with different UFA profiles were selected from 314 cultivars and landraces originated from Southern China, and RNA-seq analysis was performed in soybean seeds at three developmental stages. Using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a series of genes and pathways related to fatty acid metabolism were identified, and 40 days after flowering (DAF) was found to be the crucial period in the formation of UFA profiles. Further, weighted gene co-expression network analysis identified three modules with six genes whose functions were highly associated with the contents of oleic and LA. The detailed functional investigation of the networks and hub genes could further improve the understanding of the underlying molecular mechanism of UFA contents and might provide some ideas for the improvement in fatty acids profiles in soybean.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.

Bioactive compounds and their libraries: An insight into prospective phytotherapeutics approach for oral mucocutaneous cancers

Oral mucocutaneous cancers (OMCs) are cancers that affect both the oral mucosa and perioral cutaneous structures. Common OMCs are squamous cell carcinoma (SCC), basal cell carcinoma (BCC) and malignant melanoma (MM). Anatomical similarities and conventions which categorizes these lesions blur the magnitude of OMCs in diverse populations. The burden of OMC is high in the sub-Saharan Africa and Indian subcontinents, and the cost of management is prohibitive in the resource-limited, developing world. Hence, there is a pressing demand for the use of cost-effective in silico approaches to identify diagnostic tools and treatment targets for diseases with high burdens in these regions. Due to their ubiquitousness and accessibility, the use of therapeutic efficacy of plant bioactive compounds in the management of OMC is both appropriate and plausible. Furthermore, screening known mechanistic disease targets with well annotated plant bioactive compound libraries is poised to improve the routine management of OMCs provided that the requisite access to database resources are available and accessible. Using natural products minimizes the side effects and morbidities associated with conventional therapies. The development of innovative treatments approaches would tremendously benefit the African and Indian populace and reduce the mortalities associated with OMCs in the developing world. Hence, we discuss herein, the potential benefits, opportunities and challenges of using bioactive compound libraries in the management of OMCs.

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

Gene co-expression network analysis reveals key pathways and hub genes in Chinese cabbage (Brassica rapa L.) during vernalization

BACKGROUND

Vernalization is a type of low temperature stress used to promote rapid bolting and flowering in plants. Although rapid bolting and flowering promote the reproduction of Chinese cabbages (Brassica rapa L. ssp. pekinensis), this process causes their commercial value to decline. Clarifying the mechanisms of vernalization is essential for its further application. We performed RNA sequencing of gradient-vernalization in order to explore the reasons for the different bolting process of two Chinese cabbage accessions during vernalization.

RESULTS

There was considerable variation in gene expression between different-bolting Chinese cabbage accessions during vernalization. Comparative transcriptome analysis and weighted gene co-expression network analysis (WGCNA) were performed for different-bolting Chinese cabbage during different vernalization periods. The biological function analysis and hub gene annotation of highly relevant modules revealed that shoot system morphogenesis and polysaccharide and sugar metabolism caused early-bolting 'XBJ' to bolt and flower faster; chitin, ABA and ethylene-activated signaling pathways were enriched in late-bolting 'JWW'; and leaf senescence and carbohydrate metabolism enrichment were found in the two Chinese cabbage-related modules, indicating that these pathways may be related to bolting and flowering. The high connectivity of hub genes regulated vernalization, including MTHFR2, CPRD49, AAP8, endoglucanase 10, BXLs, GATLs, and WRKYs. Additionally, five genes related to flower development, BBX32 (binds to the FT promoter), SUS1 (increases FT expression), TSF (the closest homologue of FT), PAO and NAC029 (plays a role in leaf senescence), were expressed in the two Chinese cabbage accessions.

CONCLUSION

The present work provides a comprehensive overview of vernalization-related gene networks in two different-bolting Chinese cabbages during vernalization. In addition, the candidate pathways and hub genes related to vernalization identified here will serve as a reference for breeders in the regulation of Chinese cabbage production.

Clubroot in Brassica: recent advances in genomics, breeding, and disease management

Clubroot disease, caused by Plasmodiophora brassicae, affects Brassica oilseed and vegetable production worldwide. This review is focused on various aspects of clubroot disease and its management, including understanding the pathogen and resistance in the host plants. Advances in genetics, molecular biology techniques, and omics research have helped to identify several major loci, QTL, and genes from the Brassica genomes involved in the control of clubroot resistance. Transcriptomic studies have helped to extend our understanding of the mechanism of infection by the pathogen and the molecular basis of resistance/susceptibility in the host plants. A comprehensive understanding of the clubroot disease and host resistance would allow developing a better strategy by integrating the genetic resistance with cultural practices to manage this disease from a long-term perspective.

Root Transcriptome and Metabolome Profiling Reveal Key Phytohormone-Related Genes and Pathways Involved Clubroot Resistance in Brassica rapa L

Plasmodiophora brassicae, an obligate biotrophic pathogen-causing clubroot disease, can seriously affect Brassica crops worldwide, especially Chinese cabbage. Understanding the transcriptome and metabolome profiling changes during the infection of P. brassicae will provide key insights in understanding the defense mechanism in Brassica crops. In this study, we estimated the phytohormones using targeted metabolome assays and transcriptomic changes using RNA sequencing (RNA-seq) in the roots of resistant (BrT24) and susceptible (Y510-9) plants at 0, 3, 9, and 20 days after inoculation (DAI) with P. brassicae. Differentially expressed genes (DEGs) in resistant vs. susceptible lines across different time points were identified. The weighted gene co-expression network analysis of the DEGs revealed six pathways including "Plant-pathogen interaction" and "Plant hormone signal transduction" and 15 hub genes including pathogenic type III effector avirulence factor gene (RIN4) and auxin-responsive protein (IAA16) to be involved in plants immune response. Inhibition of Indoleacetic acid, cytokinin, jasmonate acid, and salicylic acid contents and changes in related gene expression in R-line may play important roles in regulation of clubroot resistance (CR). Based on the combined metabolome profiling and hormone-related transcriptomic responses, we propose a general model of hormone-mediated defense mechanism. This study definitely enhances our current understanding and paves the way for improving CR in Brassica rapa.

Comparative Transcriptome Analysis of Rutabaga (Brassica napus) Cultivars Indicates Activation of Salicylic Acid and Ethylene-Mediated Defenses in Response to Plasmodiophora brassicae

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important soilborne disease of Brassica napus L. and other crucifers. To improve understanding of the mechanisms of resistance and pathogenesis in the clubroot pathosystem, the rutabaga (B. napus subsp. rapifera Metzg) cultivars ‘Wilhelmsburger’ (resistant) and ‘Laurentian’ (susceptible) were inoculated with P. brassicae pathotype 3A and their transcriptomes were analyzed at 7, 14, and 21 days after inoculation (dai) by RNA sequencing (RNA-seq). Thousands of transcripts with significant changes in expression were identified in each host at each time-point in inoculated vs. non-inoculated plants. Molecular responses at 7 and 14 dai supported clear differences in the clubroot response mechanisms of the two genotypes. Both the resistant and the susceptible cultivars activated receptor-like protein (RLP) genes, resistance (R) genes, and genes involved in salicylic acid (SA) signaling as clubroot defense mechanisms. In addition, genes related to calcium signaling and genes encoding leucine-rich repeat (LRR) receptor kinases, the respiratory burst oxidase homolog (RBOH) protein, and transcription factors such as WRKYs, ethylene responsive factors, and basic leucine zippers (bZIPs), appeared to be upregulated in ‘Wilhelmsburger’ to restrict P. brassicae development. Some of these genes are essential components of molecular defenses, including ethylene (ET) signaling and the oxidative burst. Our study highlights the importance of activation of genes associated with SA- and ET-mediated responses in the resistant cultivar. A set of candidate genes showing contrasting patterns of expression between the resistant and susceptible cultivars was identified and includes potential targets for further study and validation through approaches such as gene editing.

Stem lodging resistance in hulless barley: Transcriptome and metabolome analysis of lignin biosynthesis pathways in contrasting genotypes

Hulless barley is an important economic and food crop for local population in the Qinghai-Tibet plateau. However, due to extreme weather conditions, its production suffers from stem lodging stress, inflicting significant yield losses. Herein, we selected five lodging resistant and five non-resistant genotypes to investigate changes in concentration of lignin related metabolites and expression levels of related genes in node samples. The lodging resistant genotypes displayed high content of lignin intermediate metabolites. 57% of the expressed genes were differentially expressed (DEG) between the two groups. 31 DEGs participate in the lignin pathways and we found that 65% of these DEGs were strongly up-regulated in the lodging resistant group, indicating a mechanism towards high lignin synthesis within said group. The candidate structural genes as well as the co-expressed TFs identified in this study represent important molecular tools for functional characterization and exploitation in molecular breeding programmes.

Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes

The contribution of surrounding plant microbiota to disease development has led to the 'pathobiome' concept, which represents the interaction between the pathogen, the host plant and the associated biotic microbial community, resulting or not in plant disease. The aim herein is to understand how the soil microbial environment may influence the functions of a pathogen and its pathogenesis, and the molecular response of the plant to the infection, with a dual-RNAseq transcriptomics approach. We address this question using Brassica napus and Plasmodiophora brassicae, the pathogen responsible for clubroot. A time-course experiment was conducted to study interactions between P. brassicae, two B. napus genotypes and three soils harbouring high, medium or low microbiota diversities and levels of richness. The soil microbial diversity levels had an impact on disease development (symptom levels and pathogen quantity). The P. brassicae and B. napus transcriptional patterns were modulated by these microbial diversities, these modulations being dependent on the host genotype plant and the kinetic time. The functional analysis of gene expressions allowed the identification of pathogen and plant host functions potentially involved in the change of plant disease level, such as pathogenicity-related genes (NUDIX effector) in P. brassicae and plant defence-related genes (glucosinolate metabolism) in B. napus.