Please wait a minute...
Journal of Integrative Agriculture  2023, Vol. 22 Issue (11): 3284-3295    DOI: 10.1016/j.jia.2022.09.014
Special Focus: Germplasm and Molecular Breeding in Horticultural Crops Advanced Online Publication | Current Issue | Archive | Adv Search |
Comprehensive analysis of the full-length transcripts and alternative splicing involved in clubroot resistance in Chinese cabbage

SU He-nan1*, YUAN Yu-xiang1*#, YANG Shuang-juan1*, WEI Xiao-chun1, ZHAO Yan-yan1, WANG Zhi-yong1, QIN Liu-yue1, YANG Zhi-yuan2, NIU Liu-jing1, LI Lin1, ZHANG Xiao-wei1# 

1 Institute of Horticulture, Henan Academy of Agricultural Sciences, Zhengzhou 450002, P.R.China 
2 School of Materials Science and Engineering, Central South University, Changsha 410083, P.R.China  
Download:  PDF in ScienceDirect  
Export:  BibTeX | EndNote (RIS)      
摘要  

字花科芸薹属根肿病是由芸薹根肿菌(Plasmodiophora brassicae)侵染引起的世界范围内普遍发生的一种重要土传病害,在全世界多个国家均有分布。以往对大白菜抗根肿病基因的研究主要采用转录组测序技术,不能提供准确的转录本组装和结构信息。本研究采用PacBio RS II SMRT测序技术,获得了大白菜根肿病DH40R接菌02581322天混合根全长转录组。总的来说,从SMRT测序数据中发现了39,376个高质量转录本和26,270开放阅读框。此外,还鉴定出426个注释的长链非编码RNA56个转录因子家族、1883个具有poly(A)位点的基因和1691可变剪切事件。另外,其中有1202个基因在DH40R中至少有一个AS事件。与以往的RNA-seq数据比较显示,6个差异表达的AS基因(1个抗病,5个防御反应)可能参与了根肿病抗性防御。



Abstract  

Chinese cabbage is an economically important Brassica vegetable worldwide, and clubroot, which is caused by the soil-borne protist plant pathogen Plasmodiophora brassicae is regarded as a destructive disease to Brassica crops.  Previous studies on the gene transcripts related to Chinese cabbage resistance to clubroot mainly employed RNA-seq technology, although it cannot provide accurate transcript assembly and structural information.  In this study, PacBio RS II SMRT sequencing was used to generate full-length transcriptomes of mixed roots at 0, 2, 5, 8, 13, and 22 days after Pbrassicae infection in the clubroot-resistant line DH40R.  Overall, 39 376 high-quality isoforms and 26 270 open reading frames (ORFs) were identified from the SMRT sequencing data.  Additionally, 426 annotated long noncoding RNAs (lncRNAs), 56 transcription factor (TF) families, 1 883 genes with poly(A) sites and 1 691 alternative splicing (AS) events were identified.  Furthermore, 1 201 of the genes had at least one AS event in DH40R.  A comparison with RNA-seq data revealed six differentially expressed AS genes (one for disease resistance and five for defensive response) that are potentially involved in Pbrassicae resistance.  The results of this study provide valuable resources for basic research on clubroot resistance in Chinese cabbage.

Keywords:  Chinese cabbage       clubroot       full-length transcriptome       SMRT sequencing       alternative splicing  
Received: 07 May 2022   Accepted: 26 August 2022
Fund: This work was supported by the National Natural Science Foundation of China (31872945 and 31801874), the earmarked fund for China Agricultural Research System (CARS-23-G15), the Funds for Distinguished Young Scientists from Henan Academy of Agricultural Sciences, China (2021JQ03), and the Innovation Team of Henan Academy of Agricultural Sciences, China (2021TD06).
About author:  SU He-nan, E-mail: 18810835083@163.com; #Correspondence YUAN Yu-xiang, E-mail: yuxiangyuan126@126.com; ZHANG Xiao-wei, E-mail: xiaowei5737@163.com * These authors contributed equally to this study.

Cite this article: 

SU He-nan, YUAN Yu-xiang, YANG Shuang-juan, WEI Xiao-chun, ZHAO Yan-yan, WANG Zhi-yong, QIN Liu-yue, YANG Zhi-yuan, NIU Liu-jing, LI Lin, ZHANG Xiao-wei. 2023. Comprehensive analysis of the full-length transcripts and alternative splicing involved in clubroot resistance in Chinese cabbage. Journal of Integrative Agriculture, 22(11): 3284-3295.

Abdel-Ghany S E, Hamilton M, Jacobi J L, Ngam P, Devitt N, Schilkey F, Ben-Hur A, Reddy A S. 2016. A survey of the sorghum transcriptome using single-molecule long reads. Nature Communications7, 11706.

Chao Y, Yuan J, Li S, Jia S, Han L, Xu L. 2018. Analysis of transcripts and splice isoforms in red clover (Trifolium pratense L.) by single-molecule long-read sequencing. BMC Plant Biology18, 300.

Chen J, Pang W X, Chen B, Zhang C, Piao Z Y. 2016. Transcriptome analysis of Brassica rapa near-isogenic lines carrying clubroot-resistant and -susceptible alleles in response to Plasmodiophora brassicae during early infection. Frontiers in Plant Science6, 1183.

Chen K, Gao J H, Sun S J, Zhang Z J, Yu B, Li J, Xie C G, Li G J, Wang P C, Song C P, Bressan R A, Hua J, Zhu J K, Zhao Y. 2020. BONZAI proteins control global osmotic stress responses in plants. Current Biology30, 4815–4825.

Chen S T, Qiu G L, Yang M L. 2019. SMRT sequencing of full-length transcriptome of seagrasses Zostera japonicaScientific Reports9, 14537.

Chu M, Song T, Falk K C, Zhang X, Liu X, Chang A, Lahlali R, McGregor L, Gossen B D, Peng G, Yu F. 2014. Fine mapping of Rcr1 and analyses of its effect on transcriptome patterns during infection by Plasmodiophora brassicaeBMC Genomics15, 1166.

Dixon G R. 2009. The occurrence and economic impact of Plasmodiophora brassicae and clubroot disease. The Journal of Plant Growth Regulation28, 194–202.

Dong L L, Liu H F, Zhang J C, Yang S J, Kong G Y, Chu J S, Chen N S, Wang D W. 2015. Single-molecule real-time transcript sequencing facilitates common wheat genome annotation and grain transcriptome research. BMC Genomics16, 1039.

Elkon R, Ugalde A P, Agami R. 2013. Alternative cleavage and polyadenylation: extent, regulation and function. Nature Reviews Genetics14, 496–506.

Endt D V, Kijne J W, Memelink J. 2002. Transcription factors controlling plant secondary metabolism: What regulates the regulators? Phytochemistry61, 107–114.

Eulgem T, Somssich I E. 2007. Networks of WRKY transcription factors in defense signaling. Current Opinion in Plant Biology10, 366–371.

Feller A, Machemer K, Braun E L, Grotewold E. 2011. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. The Plant Journal66, 94–116.

Foissac S, Sammeth M. 2007. ASTALAVISTA: dynamic and fexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Research35, W297–W299.

Fu F F, Pen Y S, Wang G B, Kassaby Y A, Cao F L. 2021. Integrative analysis of the metabolome and transcriptome reveals seed germination mechanism in Punica granatum L. Journal of Integrative Agriculture, 20, 132–146.

Galindo-González L, Manolii V, Hwang S F, Strelkov S E. 2020. Response of Brassica napus to Plasmodiophora brassicae involves salicylic acid-mediated immunity: an RNA-seq-based study. Frontiers in Plant Science11, 1025.

Haas B J, Papanicolaou A, Yassour M, Grabherr M, Blood P D, Bowden J, Couger M B, Eccles D, Li B, Lieber M, MacManes M D, Ott M, Orvis J, Pochet N, Strozzi F, Weeks N, Westerman R, William T, Dewey C N, Henschel R, et al. 2013. De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity. Nature Protocols8, 1494–1512.

Hackmann C, Korneli C, Kutyniok M, Köster T, Wiedenlübbert M, Müller C, Staiger D. 2014. Salicylic acid-dependent and -independent impact of an RNA-binding protein on plant immunity. Plant Cell Environment37, 696–706.

Howard B E, Hu Q, Babaoglu A C, Chandra M, Borghi M, Tan X, He L, Winter-Sederoff H, Gassmann W, Veronese P, Heber S. 2013. High-throughput RNA sequencing of pseudomonas-infected Arabidopsis reveals hidden transcriptome complexity and novel splice variants. PLoS ONE8, e74183.

Hu L J, Zhang X W, Yuan Y X, Wang Z Y, Yang S J, Li R N, Nath U K, Zhao Y Y, Tian B M, Shi G, Xie Z Q, Wei F, Wei X C. 2021. Comparative transcriptome identifies gene expression networks regulating developmental pollen abortion in ogura cytoplasmic male sterility in Chinese cabbage (Brassica rapa ssp. pekinensis). Horticulturae7, 157.

Huang F Y, Abbas F, Rothenberg D O, Imran M, Fiaz S, Rehman N, Amanullah S, Younas A, Ding Y, Cai X J, Chen X L, Yu L, Ye X W, Jiang L, Ke Y G, He Y S. 2022. Molecular cloning, characterization and expression analysis of two 12-oxophytodienoate reductases (NtOPR1 and NtOPR2) from Nicotiana tabacumMolecular Biology Reports49, 5379–5387.

Huang G T, Ma S L, Bai L P, Zhang L, Ma H, Jia P, Liu J, Zhong M, Guo Z F. 2012. Signal transduction during cold, salt, and drought stresses in plants. Molecular Biology Reports39, 969–987.

Kong L, Zhang Y, Ye Z Q, Liu X Q, Zhao S Q, Wei L, Gao G. 2007. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Research35, 345–349.

Langfelder P, Horvath S. 2008. WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics9, 559.

Lauvergeat V, Lacomme C, Lacombe E, Lasserre E, Roby D, Grima-Pettenati J. 2001. Two cinnamoyl-CoA reductase (CCR) genes from Arabidopsis thaliana are differentially expressed during development and in response to infection with pathogenic bacteria. Phytochemistry57, 1187–1195.

Liu X, Hao L, Li D, Zhu L, Hu S. 2015. Long non-coding RNAs and their biological roles in plants. Genomics Proteomics Bioinformatics13, 137–147.

Liu X, Zhao C B, Yang L M, Zhuang M, Zhang Y Y, Wang Y, Fang Z Y, Lv H H. 2020. A time-resolved dual transcriptome analysis reveals the molecular regulating network underlying the compatible/incompatible interactions between cabbage (Brassica oleracea) and Fusarium oxysporum f. sp. conglutinans. Plant Soil448, 455–478.

Ma J, Xiang Y, Xiong Y, Lin Z, Xue Y, Mao M, Sun L, Zhou Y, Li X, Huang Z. 2019. SMRT sequencing analysis reveals the full-length transcripts and alternative splicing patterns in Ananas comosus var. bracteatusPeer J7, e7062.

Mandadi K K, Scholthof K B. 2015. Genome-wide analysis of alternative splicing landscapes modulated during plant–virus interactions in Brachypodium distachyonThe Plant Cell27, 71–85.

Meyers B C, Kozik A, Griego A, Kuang H H, Michelmore R W. 2003. Genome-wide analysis of NBS-LRR-encoding genes in ArabidopsisThe Plant Cell15, 809–834.

Ni L J, Wang Z Q, Guo J B, Pei X, Liu L, Li H, Yuan H, Gu C. 2021. Full-length transcriptome sequencing and comparative transcriptome analysis to evaluate drought and salt stress in Iris lactea var. ChinensisGenes12, 434.

Ning Y, Wang Y, Fang Z Y, Zhuang M, Zhang Y Y, Lv H H, Liu Y M, Li Z S, Yang L M. 2019. Comparative transcriptome analysis of cabbage (Brassica oleracea var. capitata) infected by Plasmodiophora brassicae reveals drastic defense response at secondary infection stage. Plant Soil443, 167–183.

Olsen A N, Ernst H A, Leggio L L, Skriver K. 2005. NAC transcription factors: structurally distinct, functionally diverse. Trends in Plant Science10, 79–87.

Pierrugues O, Pierrugues C, Oshiro J, Gouy M, Deveaux Y, Carman G M, Thuriaux P, Kazmaier M. 2001. Lipid phosphate phosphatases in ArabidopsisThe Journal of Biological Chemistry23, 20300–20308.

Song J, Liu H, Zhuang H, Zhao C, Xu Y, Wu S, Qi J, Li J, Hettenhausen C, Wu J. 2017. Transcriptomics and alternative splicing analyses reveal large differences between maize lines B73 and Mo17 in response to aphid Rhopalosiphum padi infestation. Frontiers in Plant Science8, 1738.

Staiger D, Brown J W. 2013. Alternative splicing at the intersection of biological timing, development, and stress responses. The Plant Cell25, 3640–3656.

Streitner C, Köster T, Simpson C G, Shaw P, Danisman S, Brown J W, Staiger D. 2012. An hnRNP-like RNA-binding protein affects alternative splicing by in vivo interaction with transcripts in Arabidopsis thalianaNucleic Acids Research40, 11240–11255.

Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, Liu Y, Chen R, Zhao Y. 2013. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Research41, e166.

Tang F, Yang S, Gao M, Zhu H. 2013. Alternative splicing is required for RCT1-mediated disease resistance in Medicago truncatulaPlant Molecular Biology82, 367–374.

Tsuda K, Somssich I E. 2015. Transcriptional networks in plant immunity. New Phytologist206, 932–947.

Wang J, Lin J, Kan J, Wang H, Li X, Yang Q, Li H, Chang Y. 2018. Genome-wide identification and functional prediction of novel drought-responsive lncRNAs in Pyrus betulifoliaGenes9, 311.

Wang L, Park H J, Dasari S, Wang S, Kocher J P, Li W. 2013. CPAT: Coding-Potential Assessment Tool using an alignment-free logistic regression model. Nucleic Acids Research41, e74.

Wang Y, Ji J L, Tong L, Fang Z Y, Yang L M, Zhuang M, Zhang Y Y, Lv H H. 2021. Global survey of the full-length cabbage transcriptome (Brassica oleracea var. capitata L.) reveals key alternative splicing events involved in growth and disease response. International Journal of Molecular Sciences22, 10443.

Wu Q, Wang J, Mao S, Xu H, Wu Q, Liang M, Yuan Y, Liu M, Huang K. 2019. Comparative transcriptome analyses of genes involved in sulforaphane metabolism at different treatment in Chinese kale using full-length transcriptome sequencing. BMC Genomics20, 377.

Wu Y W, Xu J, Han X M, Qiao G, Yang K, Wen Z, Wen X P. 2020. Comparative transcriptome analysis combining SMRT- and Illumina-based RNA-Seq identifies potential candidate genes involved in betalain biosynthesis in pitaya fruit. International Journal of Molecular Sciences21, 3288.

Yang B, Jiang Y Q, Rahman M H, Deyholos M K, Kav N N. 2009. Identification and expression analysis of WRKY transcription factor genes in canola (Brassica napus L.) in response to fungal pathogens and hormone treatments. BMC Plant Biology9, 68.

Yang S, Tang F, Zhu H. 2014. Alternative splicing in plant immunity. International Journal of Molecular Sciences15, 10424–10445.

Yang Z, Yang Z, Yang C, Wang Z, Chen D, Xie Y, Wu Y. 2019. Identification and genetic analysis of alternative splicing of long non-coding RNAs in tomato initial flowering stage. Genomics112, 897–907.

Yuan Y X, Qin L Y, Su H N, Yang S J, Wei X C, Wang Z Y, Zhao Y Y, Li L, Liu H L, Tian B M, Zhang X W. 2021. Transcriptome and coexpression network analyses reveal hub genes in Chinese cabbage (Brassica rapa L. ssp. pekinensis) during different stages of Plasmodiophora brassicae infection. Frontiers in Plant Science12, 650252.

Zhang X C, Gassmann W. 2007. Alternative splicing and mRNA levels of the disease resistance gene RPS4 are induced during defense responses. Plant Physiology145, 1577–1587.

Zhang X L, Liu Y M, Fang Z Y, Li Z S, Yang L M, Zhuang M, Zhang Y Y, Lv H H. 2016. Comparative transcriptome analysis between Broccoli (Brassica oleracea varitalica) and wild cabbage (Brassica macrocarpa Guss.) in response to Plasmodiophora brassicae during different infection stages. Frontiers in Plant Science7, 1929.

Zhao J, Liu T, Liu W C, Zhang D P, Dong D, Wu H L, Zhang T T, Liu D W. 2021. Transcriptomic insights into growth promotion effect of Trichoderma afroharzianum TM2–4 microbial agent on tomato plants. Journal of Integrative Agriculture20, 1266–1276.

Zheng Y, Wang Y, Ding B, Fei Z. 2017. Comprehensive transcriptome analyses reveal that potato spindle tuber viroid triggers genome-wide changes in alternative splicing, inducible trans-acting activity of phased secondary small interfering RNAs, and immune responses. Journal of Virology91, e00247.

Zhu C, Li X, Zheng J. 2018. Transcriptome profiling using Illumina- and SMRT-based RNA-seq of hot pepper for in-depth understanding of genes involved in CMV infection. Gene666, 123–133.

Zhu Z, Chen H, Xie K, Liu C, Li L, Liu L, Han X, Jiao C, Wan Z, Sha A. 2020. Characterization of drought-responsive transcriptome during seed germination in Adzuki bean (L.) by PacBio SMRT and Illumina sequencing. Frontiers in Genetics11, 996.

[1] MA Xiao-wen, MA Qiu-xiang, MA Mu-qing, CHEN Yan-hang, GU Jin-bao, LI Yang, HU Qing, LUO Qing-wen, WEN Ming-fu, ZHANG Peng, LI Cong, WANG Zhen-yu.

Cassava MeRS40 is required for the regulation of plant salt tolerance [J]. >Journal of Integrative Agriculture, 2023, 22(5): 1396-1411.

[2] ZHANG Rui-xing, ZHANG Ni-nan, WANG Ya-xiu, Khan ABID, MA Shuai, BAI Xue, ZENG Qi, PAN Qi-ming, LI Bao-hua, ZHANG Lu-gang. Blue light induces leaf color change by modulating carotenoid metabolites in orange-head Chinese cabbage (Brassica rapa L. ssp. pekinensis)[J]. >Journal of Integrative Agriculture, 2023, 22(11): 3296-3311.
[3] YANG Xiao-xiang, HUANG Xiao-qin, WU Wen-xian, XIANG Yun-jia, DU Lei, ZHANG Lei, LIU Yong. Effects of different rotation patterns on the occurrence of clubroot disease and diversity of rhizosphere microbes[J]. >Journal of Integrative Agriculture, 2020, 19(9): 2265-2273.
[4] XU Yu-chao, HOU Xi-lin, XU Wei-wei, SHEN Lu-lu, Lü Shan-wu, ZHANG Shi-lin, HU Chun-mei. Isolation and characterization of an ERF-B3 gene associated with flower abnormalities in non-heading Chinese cabbage[J]. >Journal of Integrative Agriculture, 2016, 15(3): 528-536.
[5] JU Xiang-hong, XU Han-jin, YONG Yan-hong, AN Li-long, XU Ying-mei, JIAO Pei-rong , LIAO Ming. Heat Stress Upregulates the Expression of TLR4 and Its Alternative Splicing Variant in Bama Miniature Pigs[J]. >Journal of Integrative Agriculture, 2014, 13(11): 2479-2487.
[6] ZHANG Peng, DENG Heng, XIAO Fang-ming , LIU Yong-sheng. Alterations of Alternative Splicing Patterns of Ser/Arg-Rich (SR) Genes in Response to Hormones and Stresses Treatments in Different Ecotypes of Rice (Oryza sativa)[J]. >Journal of Integrative Agriculture, 2013, 12(5): 737-748.
[7] WANG Li, GE Ting-ting, PENG Hai-tao, WANG Cheng, LIU Tong-kun, HOU Xi-lin. Molecular Cloning, Expression Analysis and Localization of Exo70A1 Related to Self Incompatibility in Non-Heading Chinese Cabbage (Brassica campestris ssp. chinensis)[J]. >Journal of Integrative Agriculture, 2013, 12(12): 2149-2156.
[8] LIU Yi-peng, LIN Ke-jian, LIU Yang, GUI Fu-rong , WANG Gui-rong. Nicotinic Acetylcholine Receptor Gene Family of the Pea Aphid, Acyrthosiphon pisum[J]. >Journal of Integrative Agriculture, 2013, 12(11): 2083-2091.
[9] LIU Yan, YIN You-ping, WANG Zhong-kang , LUO Yuan-li. Expression of Nitrilases in Brassica juncea var. tumida Tsen in Root Galls Caused by Plasmodiophora brassicae[J]. >Journal of Integrative Agriculture, 2012, 12(1): 100-108.
No Suggested Reading articles found!