Transcriptome and Coexpression Network Analyses Reveal Hub Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis) During Different Stages of Plasmodiophora brassicae Infection

Transcriptome Analysis of Chinese Cabbage Infected with Plasmodiophora Brassicae in the Primary Stage

Clubroot disease caused by the infection of Plasmodiophora brassicae is widespread in China, and significantly reduces the yield of Chinese cabbage (Brassica rapa L. ssp. pekinensis). However, the resistance mechanism of Chinese cabbage against clubroot disease is still unclear. It is important to exploit the key genes that response to early infection of P. brassicae. In this study, it was found that zoospores were firstly invaded hair roots on the 8th day after inoculating with 1 × 107 spores/mL P. brassicae. Transcriptome analysis found that the early interaction between Chinese cabbage and P. brassicae caused the significant expression change of some defense genes, such as NBS-LRRs and pathogenesis-related genes, etc. The above results were verified by quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Otherwise, peroxidase (POD) salicylic acid (SA) and jasmonic acid (JA) were also found to be important signal molecules in the resistance to clubroot disease in Chinese cabbage. This study provides important clues for understanding the resistance mechanism of Chinese cabbage against clubroot disease.

Comparative Metabolome and Transcriptome Analysis Reveals the Defense Mechanism of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) against Plasmodiophora brassicae Infection

Chinese cabbage (Brassica rapa L. ssp. pekinensis) ranks among the most cultivated and consumed vegetables in China. A major threat to its production is Plasmodiophora brassicae, which causes large root tumors, obstructing nutrient and water absorption and resulting in plant withering. This study used a widely targeted metabolome technique to identify resistance-related metabolites in resistant (DH40R) and susceptible (DH199S) Chinese cabbage varieties after inoculation with P. brassicae. This study analyzed disease-related metabolites during different periods, identifying 257 metabolites linked to resistance, enriched in the phenylpropanoid biosynthesis pathway, and 248 metabolites linked to susceptibility, enriched in the arachidonic acid metabolism pathway. Key metabolites and genes in the phenylpropanoid pathway were upregulated at 5 days post-inoculation (DPI), suggesting their role in disease resistance. In the arachidonic acid pathway, linoleic acid and gamma-linolenic acid were upregulated at 5 and 22 DPI in resistant plants, while arachidonic acid was upregulated at 22 DPI in susceptible plants, leading to the conclusion that arachidonic acid may be a response substance in susceptible plants after inoculation. Many genes enriched in these pathways were differentially expressed in DH40R and DH199S. The research provided insights into the defense mechanisms of Chinese cabbage against P. brassicae through combined metabolome and transcriptome analysis.

Deciphering the role of rhizosphere microbiota in modulating disease resistance in cabbage varieties

BACKGROUND

Cabbage Fusarium wilt (CFW) is a devastating disease caused by the soil-borne fungus Fusarium oxysporum f. sp. conglutinans (Foc). One of the optimal measures for managing CFW is the employment of tolerant/resistant cabbage varieties. However, the interplay between plant genotypes and the pathogen Foc in shaping the rhizosphere microbial community, and the consequent influence of these microbial assemblages on biological resistance, remains inadequately understood.

RESULTS

Based on amplicon metabarcoding data, we observed distinct differences in the fungal alpha diversity index (Shannon index) and beta diversity index (unweighted Bray-Curtis dissimilarity) within the rhizosphere of the YR (resistant to Foc) and ZG (susceptible to Foc) cabbage varieties, irrespective of Foc inoculation. Notably, the Shannon diversity shifts in the resistant YR variety were more pronounced following Foc inoculation. Disease-resistant plant variety demonstrate a higher propensity for harboring beneficial microorganisms, such as Pseudomonas, and exhibit superior capabilities in evading harmful microorganisms, in contrast to their disease-susceptible counterparts. Furthermore, the network analysis was performed on rhizosphere-associated microorganisms, including both bacteria and fungi. The networks of association recovered from YR exhibited greater complexity, robustness, and density, regardless of Foc inoculation. Following Foc infection in the YR rhizosphere, there was a notable increase in the dominant bacterium NA13, which is also a hub taxon in the microbial network. Reintroducing NA13 into the soil significantly improved disease resistance in the susceptible ZG variety, by directly inhibiting Foc and triggering defense mechanisms in the roots.

CONCLUSIONS

The rhizosphere microbial communities of these two cabbage varieties are markedly distinct, with the introduction of the pathogen eliciting significant alterations in their microbial networks which is correlated with susceptibility or resistance to soil-borne pathogens. Furthermore, we identified a rhizobacteria species that significantly boosts disease resistance in susceptible cabbages. Our results indicated that the induction of resistance genes leading to varied responses in microbial communities to pathogens may partly explain the differing susceptibilities of the cabbage varieties tested to CFW. Video Abstract.

Comparative Analysis of Transcriptomes Reveals Pathways and Verifies Candidate Genes for Clubroot Resistance in Brassica oleracea

Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.

Comparison of Root Transcriptomes against Clubroot Disease Pathogens in a Resistant Chinese Cabbage Cultivar (Brassica rapa cv. 'Akimeki')

Clubroot, caused by Plasmodiophora brassicae, is one of the diseases that causes major economic losses in cruciferous crops worldwide. Although prevention strategies, including soil pH adjustment and crop rotation, have been used, the disease's long persistence and devastating impact continuously remain in the soil. CR varieties were developed for clubroot-resistant (CR) Chinese cabbage, and 'Akimeki' is one of the clubroot disease-resistant cultivars. However, recent studies have reported susceptibility to several Korean pathotypes in Akimeki and the destruction of the resistance to P. brassicae in many Brassica species against CR varieties, requiring the understanding of more fine-tuned plant signaling by fungal pathogens. In this study, we focused on the early molecular responses of Akimeki during infection with two P. brassicae strains, Seosan (SS) and Hoengseong2 (HS2), using RNA sequencing (RNA-seq). Among a total of 2358 DEGs, 2037 DEGs were differentially expressed following SS and HS2 infection. Gene ontology (GO) showed that 1524 and 513 genes were up-regulated following SS and HS2 inoculations, respectively. Notably, the genes of defense response and jasmonic acid regulations were enriched in the SS inoculation condition, and the genes of water transport and light intensity response were enriched in the HS2 inoculation condition. Moreover, KEGG pathways revealed that the gene expression set were related to pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) mechanisms. The results will provide valuable information for developing CR cultivars in Brassica plants.

Fine mapping and candidate gene analysis of CRA8.1.6, which confers clubroot resistance in turnip (Brassica rapa ssp. rapa)

Clubroot disease poses a significant threat to Brassica crops, necessitating ongoing updates on resistance gene sources. In F2 segregants of the clubroot-resistant inbred line BrT18-6-4-3 and susceptible DH line Y510, the genetic analysis identified a single dominant gene responsible for clubroot resistance. Through bulk segregant sequencing analysis and kompetitive allele-specific polymerase chain reaction assays, CRA8.1.6 was mapped within 110 kb (12,255-12,365 Mb) between markers L-CR11 and L-CR12 on chromosome A08. We identified B raA08g015220.3.5C as the candidate gene of CRA8.1.6. Upon comparison with the sequence of disease-resistant material BrT18-6-4-3, we found 249 single-nucleotide polymorphisms, seven insertions, six deletions, and a long terminal repeat (LTR) retrotransposon (5,310 bp) at 909 bp of the first intron. However, the LTR retrotransposon was absent in the coding sequence of the susceptible DH line Y510. Given the presence of a non-functional LTR insertion in other materials, it showed that the LTR insertion might not be associated with susceptibility. Sequence alignment analysis revealed that the fourth exon of the susceptible line harbored two deletions and an insertion, resulting in a frameshift mutation at 8,551 bp, leading to translation termination at the leucine-rich repeat domain's C-terminal in susceptible material. Sequence alignment of the CDS revealed a 99.4% similarity to Crr1a, which indicate that CRA8.1.6 is likely an allele of the Crr1a gene. Two functional markers, CRA08-InDel and CRA08-KASP1, have been developed for marker-assisted selection in CR turnip cultivars. Our findings could facilitate the development of clubroot-resistance turnip cultivars through marker-assisted selection.

Deciphering the regulatory networks involved in mild and severe salt stress responses in the roots of wild grapevine Vitis vinifera spp. sylvestris

Transcriptional regulatory networks are pivotal components of plant's response to salt stress. However, plant adaptation strategies varied as a function of stress intensity, which is mainly modulated by climate change. Here, we determined the gene regulatory networks based on transcription factor (TF) TF_gene co-expression, using two transcriptomic data sets generated from the salt-tolerant "Tebaba" roots either treated with 50 mM NaCl (mild stress) or 150 mM NaCl (severe stress). The analysis of these regulatory networks identified specific TFs as key regulatory hubs as evidenced by their multiple interactions with different target genes related to stress response. Indeed, under mild stress, NAC and bHLH TFs were identified as central hubs regulating nitrogen storage process. Moreover, HSF TFs were revealed as a regulatory hub regulating various aspects of cellular metabolism including flavonoid biosynthesis, protein processing, phenylpropanoid metabolism, galactose metabolism, and heat shock proteins. These processes are essentially linked to short-term acclimatization under mild salt stress. This was further consolidated by the protein-protein interaction (PPI) network analysis showing structural and plant growth adjustment. Conversely, under severe salt stress, dramatic metabolic changes were observed leading to novel TF members including MYB family as regulatory hubs controlling isoflavonoid biosynthesis, oxidative stress response, abscisic acid signaling pathway, and proteolysis. The PPI network analysis also revealed deeper stress defense changes aiming to restore plant metabolic homeostasis when facing severe salt stress. Overall, both the gene co-expression and PPI network provided valuable insights on key transcription factor hubs that can be employed as candidates for future genetic crop engineering programs.

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.

bra-miR167a Targets ARF8 and Negatively Regulates Arabidopsis thaliana Immunity against Plasmodiophora brassicae

Clubroot is a soil-borne disease caused by Plasmodiophora brassicae, which can seriously affect the growth and production of cruciferous crops, especially Chinese cabbage crops, worldwide. At present, few studies have been conducted on the molecular mechanism of this disease's resistance response. In this experiment, we analyzed the bioinformation of bra-miR167a, constructed a silencing vector (STTM167a) and an overexpression vector (OE-miR167a), and transformed them to Arabidopsis to confirm the role of miR167a in the clubroot resistance mechanism of Arabidopsis. Afterwards, phenotype analysis and expression level analysis of key genes were conducted on transgenic plants. From the result, we found that the length and number of lateral roots of silence transgenic Arabidopsis STTM167a was higher than that of WT and OE-miR167a. In addition, the STTM167a transgenic Arabidopsis induced up-regulation of disease resistance-related genes (PR1, PR5, MPK3, and MPK6) at 3 days after inoculation. On the other hand, the auxin pathway genes (TIR1, AFB2, and AFB3), which are involved in maintaining the balance of auxin/IAA and auxin response factor (ARF), were down-regulated. These results indicate that bra-miR167a is negative to the development of lateral roots and auxins, but positive to the expression of resistance-related genes. This also means that the STTM167a can improve the resistance of clubroot by promoting lateral root development and the level of auxin, and can induce resistance-related genes by regulating its target genes. We found a positive correlation between miR167a and clubroot disease, which is a new clue for the prevention and treatment of clubroot disease.

Differences in the transcriptional immune response to Albugo candida between white rust resistant and susceptible cultivars in Brassica rapa L

Albugo candida causing white rust disease decreases the yield of Brassica rapa vegetables greatly. Resistant and susceptible cultivars in B. rapa vegetables have different immune responses against A. candida inoculation, however, the mechanism of how host plants respond to A. candida is still unknown. Using RNA-sequencing, we identified differentially expressed genes (DEGs) between A. candida inoculated [48 and 72 h after inoculation (HAI)] and non-inoculated samples in resistant and susceptible cultivars of komatsuna (B. rapa var. perviridis). Functional DEGs differed between the resistant and susceptible cultivars in A. candida inoculated samples. Salicylic acid (SA) responsive genes tended to be changed in their expression levels by A. candida inoculation in both resistant and susceptible cultivars, but different genes were identified in the two cultivars. SA-dependent systemic acquired resistance (SAR) involving genes were upregulated following A. candida inoculation in the resistant cultivar. Particular genes categorized as SAR that changed expression levels overlapped between A. candida and Fusarium oxysporum f. sp. conglutinans inoculated samples in resistant cultivar, suggesting a role for SAR in defense response to both pathogens particularly in the effector-triggered immunity downstream pathway. These findings will be useful for understanding white rust resistance mechanisms in B. rapa.

The root meristem growth factor BrRGF6 positively regulates Chinese cabbage to infection of clubroot disease caused by Plasmodiophora Brassicae

Chinese cabbage has a high annual demand in China. However, clubroot disease caused by the infection of Plasmodiophora brassicae seriously affects its yield. Transcriptome analysis identified a root meristem growth factor 6 (BrRGF6) as significantly up-regulated in Chinese cabbage roots infected with Plasmodiophora brassicae. Quantitative reverse-transcription polymerase chain reaction and in situ hybridization analysis showed higher BrRGF6 expression in susceptible materials than in resistant materials. After Plasmodiophora brassicae infection, BrRGF6 expression was significantly up-regulated, especially in susceptible materials. Gene function analysis showed that the roots of Arabidopsis mutant rgf6 grew faster than the wild-type, and delayed the infection progress of Plasmodiophora brassicae. A Protein, nuclear transcription factor Y subunit C (BrNF-YC), was screened from yeast two-hybrid library of Chinese cabbage induced by Plasmodiophora brassicae, and verified to interact with BrRGF6 by yeast two-hybrid co-transfer. Yeast one-hybrid and β-Glucuronidase activity analysis showed that BrNF-YC could directly bind to and strongly activate the promoter of BrRGF6. Transgenic verification showed that BrRGF6 or BrNF-YC silenced Chinese cabbage significantly decreased the expression of BrRGF6, accelerated root development, and reduced incidence of clubroot disease. However, after overexpression of BrRGF6 or BrNF-YC, the phenotype showed a reverse trend. Therefore, BrRGF6 silencing accelerated root growth and enhanced resistance to clubroot disease, which was regulated by BrNF-YC.

Integrative Transcriptome, miRNAs, Degradome, and Phytohormone Analysis of Brassica rapa L. in Response to Plasmodiophora brassicae

Clubroot is an infectious root disease caused by Plasmodiophora brassicae in Brassica crops, which can cause immeasurable losses. We analyzed integrative transcriptome, small RNAs, degradome, and phytohormone comprehensively to explore the infection mechanism of P. brassicae. In this study, root samples of Brassica rapa resistant line material BrT24 (R-line) and susceptible line material Y510-9 (S-line) were collected at four different time points for cytological, transcriptome, miRNA, and degradome analyses. We found the critical period of disease resistance and infection were at 0-3 DAI (days after inoculation) and 9-20 DAI, respectively. Based on our finding, we further analyzed the data of 9 DAI vs. 20 DAI of S-line and predicted the key genes ARF8, NAC1, NAC4, TCP10, SPL14, REV, and AtHB, which were related to clubroot disease development and regulating disease resistance mechanisms. These genes are mainly related to auxin, cytokinin, jasmonic acid, and ethylene cycles. We proposed a regulatory model of plant hormones under the mRNA-miRNA regulation in the critical period of P. brassicae infection by using the present data of the integrative transcriptome, small RNAs, degradome, and phytohormone with our previously published results. Our integrative analysis provided new insights into the regulation relationship of miRNAs and plant hormones during the process of disease infection with P. brassicae.

A Pan-Transcriptome Analysis Indicates Efficient Downregulation of the FIB Genes Plays a Critical Role in the Response of Alfalfa to Cold Stress

Alfalfa (Medicago sativa L.) is a perennial forage legume that is widely distributed throughout the world, and cold stress is an important environmental factor limiting the growth and production of alfalfa in cold regions. However, little is known of the molecular mechanisms regarding cold tolerance in alfalfa. Here, we conducted physiological metabolism assays and pan-transcriptome sequencing on eight cultivars of alfalfa under cold stress conditions. The results of the RNA-seq analysis showed that the genes are "oxidoreductase activity" and "transcription regulator activity", suggesting that genes with such functions are more likely to play important roles in the response to cold stress by alfalfa. In addition, to identify specific gene modules and hub genes in response to alfalfa cold stress, we applied weighted gene co-expression network (WGCNA) analyses to the RNA-seq data. Our results indicate that the modules of genes that focus on the ATPase complex, ribosome biogenesis, are more likely to be involved in the alfalfa response to cold stress. It is important to note that we identified two fibronectin (FIB) genes as hub genes in alfalfa in response to cold stress and that they negatively regulate alfalfa response to chilling stress, and it is possible that dormant alfalfa is more effective at down-regulating FIB expression and therefore more resistant to cold stress.

MYB Transcription Factors Becoming Mainstream in Plant Roots

The function of the root system is crucial for plant survival, such as anchoring plants, absorbing nutrients and water from the soil, and adapting to stress. MYB transcription factors constitute one of the largest transcription factor families in plant genomes with structural and functional diversifications. Members of this superfamily in plant development and cell differentiation, specialized metabolism, and biotic and abiotic stress processes are widely recognized, but their roles in plant roots are still not well characterized. Recent advances in functional studies remind us that MYB genes may have potentially key roles in roots. In this review, the current knowledge about the functions of MYB genes in roots was summarized, including promoting cell differentiation, regulating cell division through cell cycle, response to biotic and abiotic stresses (e.g., drought, salt stress, nutrient stress, light, gravity, and fungi), and mediate phytohormone signals. MYB genes from the same subfamily tend to regulate similar biological processes in roots in redundant but precise ways. Given their increasing known functions and wide expression profiles in roots, MYB genes are proposed as key components of the gene regulatory networks associated with distinct biological processes in roots. Further functional studies of MYB genes will provide an important basis for root regulatory mechanisms, enabling a more inclusive green revolution and sustainable agriculture to face the constant changes in climate and environmental conditions.

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.

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.

Genetic Diversity and DNA Fingerprinting in Broccoli Carrying Multiple Clubroot Resistance Genes Based on SSR Markers

To identify cultivars quickly and accurately, DNA fingerprinting of 10 broccoli varieties was performed by using simple sequence repeat (SSR) marker technology. Highly informative and polymorphic SSR markers were screened using broccoli and rapeseed. Out of the 93 SSR marker pairs, 21 pairs were selected and found to have good polymorphism. Each marker pair generated 1 to 10 polymorphic bands with an average of 4.29. The average polymorphism information content (PIC) was 0.41 with a range from 0.16 to 0.95. Six selected marker pairs established the fingerprinting of the 10 accessions and their unique fingerprints. Cluster analysis of 10 accessions showed that the genetic similarity coefficient was between 0.57 and 0.91. They can be divided into 3 groups at the genetic similarity coefficient (GSC) of 0.73. The above results indicated that DNA fingerprinting could provide a scientific basis for the identification of broccoli polymerized multiple clubroot resistance genes. Research shows that SSR marker-based DNA fingerprinting further ensures plant seed purity.

Proteomics and Co-expression Network Analysis Reveal the Importance of Hub Proteins and Metabolic Pathways in Nicotine Synthesis and Accumulation in Tobacco (Nicotiana tabacum L.)

Nicotine is a unique alkaloid present in tobacco that is widely used in cigarettes and in the agricultural, chemical, and pharmaceutical industries. However, the research on nicotine is mostly limited to its synthesis pathways, and only a few studies have explored the effects of other metabolic pathways on nicotine precursors. Regulating the nicotine content in tobacco can greatly promoting the application of nicotine in other fields. In this study, we performed global data-independent acquisition proteomics analysis of four tobacco varieties. Of the four varieties, one had high nicotine content and three had a low nicotine content. A total of 31,259 distinct peptides and 6,018 proteins across two samples were identified. A total of 45 differentially expressed proteins (DEPs) co-existed in the three comparison groups and were mainly involved in the transport and metallic processes of the substances. Most DEPs were enriched in the biosynthesis of secondary metals, glutathione metabolism, carbon metabolism, and glycolysis/gluconeogenesis. In addition, the weighted gene co-expression network analysis identified an expression module closely related to the nicotine content (Brown, r = 0.74, P = 0.006). Gene Ontology annotation and Kyoto Encyclopaedia of Genes and Genomes enrichment analysis showed that the module proteins were mainly involved in the synthesis and metabolism of nicotine precursors such as arginine, ornithine aspartate, proline, and glutathione. The increased levels of these precursors lead to the synthesis and accumulation of nicotine in plants. More importantly, these proteins regulate nicotine synthesis by affecting the formation of putrescine, which is the core intermediate product in nicotine anabolism. Our results provide a reference for tobacco variety selection with a suitable nicotine content and regulation of the nicotine content. Additionally, the results highlight the importance of other precursor metabolism in nicotine synthesis.

Histopathology of the Plasmodiophora brassicae-Chinese Cabbage Interaction in Hosts Carrying Different Sources of Resistance

Clubroot is a serious soil-borne disease of crucifers caused by the obligate parasite Plasmodiophora brassicae. The genetic basis and histopathology of clubroot resistance in two Chinese cabbage (Brassica rapa ssp. pekinensis) inbred lines Bap055 and Bap246, challenged with pathotype 4 of P. brassicae, was evaluated. The Chinese cabbage cultivar "Juxin" served as a susceptible check. The resistance in Bap055 was found to be controlled by the CRa gene, while resistance in Bap246 fit a model of control by unknown recessive gene. Infection of the roots by P. brassicae was examined by inverted microscopy. Despite their resistance, primary and secondary infection were observed to occur in Bap055 and Bap246. Primary infection was detected at 2 days post-inoculation (DPI) in "Juxin," at 4 DPI in Bap055, and at 6 DPI in Bap246. Infection occurred most quickly on "Juxin," with 60% of the root hairs infected at 10 DPI, followed by Bap055 (31% of the root hairs infected at 12 DPI) and Bap246 (20% of the root hairs infected at 14 DPI). Secondary infection of "Juxin" was first observed at 8 DPI, while in Bap055 and Bap246, secondary infection was first observed at 10 DPI. At 14 DPI, the percentage of cortical infection in "Juxin," Bap055 and Bap246 was 93.3, 20.0, and 11.1%, respectively. Although cortical infection was more widespread in Bap055 than in Bap246, secondary infection in both of these hosts was restricted relative to the susceptible check, and the vascular system remained intact. A large number of binucleate secondary plasmodia were observed in "Juxin" and the vascular system was disrupted at 16 DPI; in Bap055 and Bap246, only a few secondary plasmodia were visible, with no binucleate secondary plasmodia. The defense mechanisms and expression of resistance appears to differ between Chinese cabbage cultivars carrying different sources of resistance.

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.