Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection

Abstract

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection of Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at 7 days after P. brassicae inoculation, whereas the expression levels of starch degradation-related genes were increased at 14 days after inoculation. Then, genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Previous genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization, as indicated by histological observations and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

References

[1] Wallenhammar AC . Prevalence of Plasmodiophora brassicae in a spring oilseed rape growing area in Central Sweden and factors influencing soil infestation levels. Plant Pathol 1996; 45: 710-9.

[2] Gossen B, Strelkov SE, Manoli VP et al. Spread of Plasmodiophora brassicae on canola in Canada, 2003-2014: old pathogen, new home. Can J Plant Pathol 2015; 37: 403-13.

[3] Botero A, Celsa G, Gossen BD et al. Clubroot disease in Latin America: distribution and management strategies. Plant Pathol 2019; 68: 827-33.

[4] Liu L, Qin L, Zhou Z et al. Refining the life cycle of Plasmodiophora brassicae . Phytopathology. 2020; 110: 1704-12.

[5] Kageyama K, Asano T . Life cycle of Plasmodiophora brassicae . J Plant Growth Regul 2009; 28: 203-11.

[6] McDonald MR, Sharma K, Gossen BD et al. The role of primary and secondary infection in host response to Plasmodiophora brassicae . Phytopathology. 2014; 104: 1078-87.

[7] Buczacki S, Ockendon J . Preliminary observations on variation in susceptibility to clubroot among collections of some wild crucifers. Ann Appl Biol. 1979; 92: 113-8.

[8] Dixon GR . The occurrence and economic impact of Plasmodiophora brassicae and clubroot disease. J Plant Growth Regul 2009; 28: 194-202.

[9] Bi K, He Z, Gao Z et al. Integrated omics study of lipid droplets from Plasmodiophora brassicae . Sci Rep. 2016; 6: 36965.

[10] Niwa R, Nomura Y, Osaki M et al. Suppression of clubroot disease under neutral pH caused by inhibition of spore germination of Plasmodiophora brassicae in the rhizosphere. Plant Pathol 2008; 57: 445-52.

[11] Moxham SE, Buczacki S, Turner R . Problems in interpreting scanning electron micrographs of Plasmodiophora plasmodia . Trans Br Mycol Soc 1981; 76: 153-5.

[12] Schuller A, Ludwig-Müller J . Histological methods to detect the clubroot pathogen Plasmodiophora brassicae during its complex life cycle. Plant Pathol. 2016; 65: 1223-37.

[13] Thalmann M, Santelia D . Starch as a determinant of plant fitness under abiotic stress. New Phytol. 2017; 214: 943-51.

[14] Fladung M, Gieffers W . Resistance reactions of leaves and tubers of rolC transgenic tetraploid potato to bacterial and fungal pathogens. Correlation with sugar, starch and chlorophyll content. Physiol Mol Plant Pathol 1993; 42: 123-32.

[15] Garg H, Li H, Sivasithamparam K et al. The infection processes of Sclerotinia sclerotiorum in cotyledon tissue of a resistant and a susceptible genotype of Brassica napus . Ann Bot 2010; 106: 897-908.

[16] Rafiq L, Shuab R, Kumar V et al. Effect of arbuscular mycorrhizal fungi on growth and development of potato (Solanum tuberosum) plant . Asian J Crop Sci. 2015; 7: 233-43.

[17] Li P, Liu W, Zhang Y et al. Fungal canker pathogens trigger carbon starvation by inhibiting carbon metabolism in poplar stems. Sci Rep. 2019; 9: 10111.

[18] Li X, An M, Xia Z et al. Transcriptome analysis of watermelon (Citrullus lanatus) fruits in response to cucumber green mottle mosaic virus (CGMMV) infection . Sci Rep 2017; 7: 16747.

[19] Brodmann D, Schuller A, Ludwig Muller J et al. Induction of trehalase in Arabidopsis plants infected with the trehalose-producing pathogen Plasmodiophora brassicae . Mol Plant Microbe Interact. 2002; 15: 693-700.

[20] Irani S, Trost B, Waldner M et al. Transcriptome analysis of response to Plasmodiophora brassicae infection in the Arabidopsis shoot and root . BMC Genomics. 2018; 19: 23.

[21] Ikegami H, Watanabe Y, Ito T et al. Ultrastructures of starch grains in Chinese cabbage and turnip root cells infected with Plasmodiophora brassicae (Studies on the Clubroot of Cruciferous Plants 6). Res Bull Fac Agr Gifu Univ. 1984; 49: 37-48.

[22] Tanaka S, Mido H, Ito SI . Colonization by two isolates of Plasmodiophora brassicae with differing pathogenicity on a clubroot-resistant cultivar of Chinese cabbage (Brassica rapa L. subsp. pekinensis) . J Gen Plant Pathol. 2006; 72: 205-9.

[23] Tsai HL, Lue WL, Lu K et al. Starch synthesis in Arabidopsis is achieved by spatial cotranscription of core starch metabolism genes. Plant Physiol. 2009; 151: 1582-95.

[24] Monroe JD, Storm AR . The Arabidopsis β-amylase (BAM) gene family: diversity of form and function. Plant Sci. 2018; 276: 163-70.

[25] Van Harsselaar JK, Lorenz J, Senning M et al. Genome-wide analysis of starch metabolism genes in potato (Solanum tuberosum L.) . BMC Genomics. 2017; 18: 37.

[26] Wang Y, Li Y, Zhang H et al. A soluble starch synthase I gene, IbSSI, alters the content, composition, granule size and structure of starch in transgenic sweet potato. Sci Rep. 2017; 7: 2315.

[27] Hirose T, Terao T . A comprehensive expression analysis of the starch synthase gene family in rice (Oryza sativa L.) . Planta. 2004; 220: 9-16.

[28] Fan X, Guo M, Li RD et al. Allelic variations in the soluble starch synthase II gene family result in changes of grain quality and starch properties in rice (Oryza sativa L.) . J Agric Sci. 2017; 155: 129.

[29] Malinowski R, Truman W et al. Genius architect or clever thief-how Plasmodiophora brassicae reprograms host development to establish a pathogen-oriented physiological sink. Mol Plant Microbe Interact. 2019; 32: 1259-66.

[30] Schwelm A, Fogelqvist J, Knaust A et al. The Plasmodiophora brassicae genome reveals insights in its life cycle and ancestry of chitin synthases. Sci Rep. 2015; 5: 11153.

[31] Rolfe SA, Strelkov SE, Links MG et al. The compact genome of the plant pathogen Plasmodiophora brassicae is adapted to intracellular interactions with host Brassica spp. BMC Genomics. 2016; 17: 272.

[32] Winkelmüller TM, Entila F, Anver S et al. Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species. Plant Cell. 2021. https://doi.org/10.1093/plcell/koab073.

[33] Graveland R, Dale P, Mithen R . Gall development in hairy root cultures infected with Plasmodiophora brassicae . Mycol Res. 1992; 96: 225-8.

[34] Skryhan K, Gurrieri L, Sparla F et al. Redox regulation of starch metabolism. Front Plant Sci. 2018; 9: 1344.

[35] Krasensky J, Jonak C . Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J Exp Bot. 2012; 63: 1593-608.

[36] Dong S, Beckles DM . Dynamic changes in the starch-sugar interconversion within plant source and sink tissues promote a better abiotic stress response. J Plant Physiol. 2019; 234-235: 80-93.

[37] Trapero C, Alcantara E, Jimenez J et al. Starch hydrolysis and vessel occlusion related to wilt symptoms in olive stems of susceptible cultivars infected by Verticillium dahliae . Front Plant Sci. 2018; 9: 72.

[38] Beilstein MA, Nagalingum NS, Clements MD et al. Dated molecular phylogenies indicate a Miocene origin for Arabidopsis thaliana . Proc Natl Acad Sci USA. 2010; 107: 18724-8.

[39] Lysak MA, Koch MA, Pecinka A et al. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005; 15: 516-25.

[40] Nagaharu U . Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization . Jpn J Bot. 1935; 7: 389-452.

[41] Zhang X, Liu Y, Fang Z et al. Comparative transcriptome analysis between broccoli (Brassica oleracea var. italica) and wild cabbage (Brassica macrocarpa Guss.) in response to Plasmodiophora brassicae during different infection stages . Front Plant Sci. 2016; 7: 1929.

[42] Zhou Q, Galindo-González L, Hwang SF et al. Application of genomics and transcriptomics to accelerate development of clubroot resistant canola. Can J Plant Pathol. 2020.

[43] Ning Y, Wong Y, Yang L et al. Comparative transcriptome analysis of cabbage (Brassica oleracea var. capitata) infected by Plasmodiophora brassicae reveals drastic defense response at secondary infection stage . Plant Soil. 2019; 443: 167-83.

[44] Siemens J, Keller I, Sarx J et al. Transcriptome analysis of Arabidopsis clubroots indicate a key role for cytokinins in disease development. Mol Plant Microbe Interact. 2006; 19: 480-94.

[45] Pang W, Liang Y, Zhan Z et al. Development of a Sinitic clubroot differential set for the pathotype classification of Plasmodiophora brassicae . Front Plant Sci. 2020; 11: 1360.

[46] Canene-Adams K . Preparation of formalin-fixed paraffin-embedded tissue for immunohistochemistry. Methods Enzymol. 2013; 533: 225-33.

[47] Finn RD, Clements J, Eddy SR . HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011; 39: W29-37.

[48] Larkin MA, Blackshields G, Brown NP et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007; 23: 2947-8.

[49] Kumar S, Stecher G, Li M et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018; 35: 1547-9.

[50] Livak KJ, Schmittgen TD . Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001; 25: 402-8.

[51] Chen C, Chen H, He Y et al. TBtools, a toolkit for biologists integrating various biological data handling tools with a user-friendly interface. BioRxiv. 2018; 289660.

[52] Wang Y, Tang H, Debarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49-9.

[53] Wang D, Zhang Y, Zhang Z et al. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010; 8: 77-80.

[54] Yu J, Zhao M, Wang X et al. Bolbase: a comprehensive genomics database for Brassica oleracea . BMC Genomics. 2013; 14: 664.