Response of Sugarcane Leaves to the Infection of Puccinia melanocephala: Transcriptomic Analysis

doi:10.11924/j.issn.1000-6850.casb2020-0622

Abstract

Abstract:

The aim is to analyze the gene expression patterns of sugarcane leaves in response to the infection of P. melanocephala, and to provide data for further research on the molecular mechanism of sugarcane resistance to brown rust. The healthy and infected sugarcane leaves by P. melanocephala were sequenced at transcriptomic level using high-throughput sequencing platform, Illumina Hiseq 2500. The functional annotation and analysis for the differentially expressed genes (DEGs) was carried out with bioinformatics methods. The relative expression levels of five DEGs were further verified by real-time fluorescence quantitative PCR (qRT-PCR). The transcriptome data showed that 4716 DEGs were obtained, including 2979 up-regulated and 1737 down-regulated DEGs. GO annotation results showed that DEGs mainly gathered in the physiological and biochemical processes of cell process, biological process, and catalytic activity, etc. KEGG enrichment analysis revealed that 1919 DEGs were annotated into 114 pathways. The relative genes in the metabolic pathways related to plants’ resistance to biotic stress were significantly up-regulated, and these metabolic pathways included phenylpropane biosynthesis, flavonoids biosynthesis, and glutathione metabolism and so on. In addition, the qRT-PCR results confirmed that the expression profiles of five up-regulated genes were consistent with that of transcriptome sequencing.

Key words: sugarcane, Puccinia melanocephala, compatible interaction, high-throughput sequencing, transcriptomics, differential expression genes

CLC Number:

S432.4

Gao Xiaoning, Wu Zilin, Huang Yonghong, Liu Rui, Qi Yongwen. Response of Sugarcane Leaves to the Infection of Puccinia melanocephala: Transcriptomic Analysis[J]. Chinese Agricultural Science Bulletin, 2021, 37(24): 102-109.

Figures/Tables 7

References 30

[1]	张学博. 福建省甘蔗锈病的初步调查研究[J]. 植物保护, 1990, 16(s):10-12.
[2]	Rott P, Bailey R A, Comstock J C, et al. A guide to sugarcane diseases[M]. Montepellier: CIRAD-ISSCT, CIRAD Publication Services, 2000.
[3]	李增平, 张淑珍. 海南甘蔗病虫害诊断图谱[M]. 北京: 中国农业出版社, 2014:31-34.
[4]	仓晓燕, 黄应昆, 李文凤, 等. 甘蔗抗褐锈病基因与分子标记研究进展[J]. 中国糖料, 2016, 38(6):62-65.
[5]	Daugrois J H, Grivet L, Roques D. A putative major gene for rust resistance linked with a RFLP marker in sugarcane cultivar R570[J]. Theoretical and Applied Genetics, 1996, 92(8):1059-1064. doi: 10.1007/BF00224049 pmid: 24166636
[6]	Asnaghi C, Dhont A, Glaszmann J, et al. Resistance of sugarcane cultivar R570 to Puccinia melanocephala isolates from different geographic locations[J]. Plant Disease, 2001, 85(3):282-286. doi: 10.1094/PDIS.2001.85.3.282 pmid: 30832043
[7]	Asnaghi C, Roques D, Ruffel S, et al. Targeted mapping of a sugarcane rust resistance gene (Bru1) using bulked segregant analysis and AFLP markers[J]. Theoretical and Applied Genetics, 2004, 108(4):759-764. pmid: 14586507
[8]	Costet L, Toubi L, Le Cunff L. Haplotype structure around Bru1 reveals a narrow genetic basis for brown rust resistance in modern sugarcane cultivars[J]. Theoretical and Applied Genetics, 2012, 125(5):825-836. doi: 10.1007/s00122-012-1875-x pmid: 22572763
[9]	Parco A S, Hale A L, Avellaneda M C, et al. Distribution and frequency of Bru1, a major brown rust resistance gene, in the sugarcane world collection[J]. Plant Breeding, 2017, 136(5):637-651. doi: 10.1111/pbr.2017.136.issue-5 URL
[10]	李文凤, 王晓燕, 黄应昆, 等. 34份甘蔗栽培原种抗褐锈病性鉴定及Bru1基因的分子检测[J]. 分子植物育种, 2015, 13(8):1814-1821.
[11]	李文凤, 王晓燕, 黄应昆, 等. 101份中国甘蔗主要育种亲本褐锈病抗性鉴定及Bru1基因的分子检测[J]. 作物学报, 2016, 42(9):1411-1416.
[12]	李文凤, 单红丽, 张荣跃, 等. 甘蔗优良新品种(系)抗褐锈病基因Bru1的分子检测及自然抗性评价[J]. 植物病理学报, 2017, 47(5):667-674.
[13]	Li W F, Shan H L, Zhang R Y, et al. Identification of field resistance and molecular detection of the brown rust resistance gene Bru1 in new elite sugarcane varieties in China[J]. Crop Protection, 2018, 103:46-50. doi: 10.1016/j.cropro.2017.09.007 URL
[14]	Que Y X, Su Y C, Guo J L, et al. A global view of transcriptome dynamics during Sporisorium scitamineum challenge in sugarcane by RNA-seq[J]. PLoS ONE, 2014, 9(8):e106476. doi: 10.1371/journal.pone.0106476 URL
[15]	Akbar S, Yao W, Yu K, et al. Photosynthetic characterization and expression profiles of sugarcane infected by sugarcane mosaic virus (SCMV) [J]. Photosynjournal Research, 2020, https://doi.org/10.1007/s11120-019-00706-w .
[16]	Chu N, Zhou J R, Fu H Y, et al. Global gene responses of resistant and susceptible sugarcane cultivars to Acidovorax avenae subsp. avenae identified using comparative transcriptome analysis[J]. Mciroorganisms, 2020, 8:10.
[17]	Ntambo M S, Meng J Y, Rott P C, et al. Comparative transcriptome profiling of resistant and susceptible sugarcane cultivars in response to infection by Xanthomonas albilineans[J]. International Journal of Molecular Sciences, 2019, 20:6138. doi: 10.3390/ijms20246138 URL
[18]	Zhang J, Zhang X, Tang H, et al. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L.[J]. Nature Genetics, 2018, 50(11):1565-1573. doi: 10.1038/s41588-018-0237-2 URL
[19]	刘晓妹, 刘文波, 施焕焕, 等. 甘蔗锈病研究进展[J]. 中国植保导刊, 2007, 27(12):12-15.
[20]	Carmona E, Vargas D, Borroto C J, et al. cDNA-AFLP analysis of differential gene expression during the interaction between sugarcane and Puccinia melanocephala[J]. Plant Breed, 2004, 123:499-501. doi: 10.1111/pbr.2004.123.issue-5 URL
[21]	Oloriz M I, Gil V, Rojas L, et al. Sugarcne genes differentially expressed in response to Puccinia melanocephala infection: identification and transcript profiling[J]. Plant Cell Reports, 2012, 31:955-969. doi: 10.1007/s00299-011-1216-6 pmid: 22212461
[22]	骆琰妍, 李魏, 戴良英. 不同激素信号途径在植物抗病中的相互作用研究进展[J]. 中国农学通报, 2013, 29(18):153-157.
[23]	Glazebrook J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens[J]. Annual Review of Phytopathology, 2005, 43:205-227. pmid: 16078883
[24]	田义, 张彩霞, 康国栋, 等. 植物TGA转录因子研究进展[J]. 中国农业科学, 2016, 49(4):632-642.
[25]	郭双元, 任慧文, 张艳琴, 等. 小麦抗病相关基因TaTGA2.2在小麦-条锈菌互作中的表达分析及功能研究[J]. 麦类作物学报, 2020, 40(6):645-655.
[26]	谢政文, 王连军, 陈锦洋, 等. 植物WRKY转录因子及其生物学功能研究进展[J]. 中国农业科技导报, 2016, 18(3):46-54.
[27]	Zheng Z, Qamar S A, Chen Z, et al. Arabidopsis WRKY33 transcription factor is required for resistance to necrotrophic fungal pathogens[J]. Plant Journal, 2006, 48(4):592-605. pmid: 17059405
[28]	Mao G H, Meng X Z, Liu Y D, et al. Phosphorylation of a WRKY transcription factor by two pathogen-responsive MAPKs drives phytoalexin biosynjournal in Arabidopsis[J]. Plant Cell, 2011, 23(4):1639-1653. doi: 10.1105/tpc.111.084996 URL
[29]	余晓丛, 娜仁, 张少英. 钙信号在植物抗病中的作用研究进展[J]. 中国农学通报, 2012, 28(03):12-16.
[30]	孙云南, 许燕, 冉隆珣, 等. 茶树叶片响应茶饼病侵染的转录组分析[J]. 茶叶科学, 2020, 40(1):113-124.

基因	正向引物(5’-3’)	反向引物 (5’-3’)
PR1	GACTGCCAGCTGGTTCACTC	AACTGCTTCTCGGACACCCA
PR10	GTTGGCAGCATCAGGCAGTT	TGTCCAGGTCCAGGAAATCG
CDPK	ATCTGCGTGGGTTAGTTCAC	ATCTGACTCAGGACCAGACC
NPR1	GCTTCTTAACTTTGTAGACAAG	CTAGCAATCCTTTGAATACATTTG
CTR1	GAGGATACTCAATCTGGACTAA	TGCTCCATTAGAATCTTACTG
GAPDH	CACGGCCACTGGAAGCA	TCCTCAGGGTTCCTGATGCC

基因	正向引物(5’-3’)	反向引物 (5’-3’)
PR1	GACTGCCAGCTGGTTCACTC	AACTGCTTCTCGGACACCCA
PR10	GTTGGCAGCATCAGGCAGTT	TGTCCAGGTCCAGGAAATCG
CDPK	ATCTGCGTGGGTTAGTTCAC	ATCTGACTCAGGACCAGACC
NPR1	GCTTCTTAACTTTGTAGACAAG	CTAGCAATCCTTTGAATACATTTG
CTR1	GAGGATACTCAATCTGGACTAA	TGCTCCATTAGAATCTTACTG
GAPDH	CACGGCCACTGGAAGCA	TCCTCAGGGTTCCTGATGCC

GO ID	GO 条目	功能分类	上调差异表达基因数目	注释的基因数目	KS值
0005524	ATP binding ATP结合	分子功能	474	10796	1.80E-26
0004364	glutathione transferase activity谷胱甘肽转移酶活性	分子功能	71	465	2.10E-17
0006749	glutathione metabolic process谷胱甘肽代谢过程	生物过程	80	496	1.20E-16
0004386	helicase activity解旋酶活性	分子功能	9	862	2.80E-12
0009407	toxin catabolic process毒素分解过程	生物过程	57	190	3.80E-12
0005737	cytoplasm细胞质	细胞组分	919	34128	6.00E-12
0030244	cellulose biosynthetic process纤维素生物合成过程	生物过程	9	193	1.70E-11
0004854	xanthine dehydrogenase activity黄嘌呤脱氢酶活性	分子功能	13	43	4.90E-11
0009115	xanthine catabolic process黄嘌呤分解过程	生物过程	13	43	5.70E-11
0004459	L-lactate dehydrogenase activity L-乳酸脱氢酶活性	分子功能	16	20	8.50E-11
0016760	cellulose synthase (UDP-forming) activity纤维素合成酶活性	分子功能	9	174	3.20E-10
0004803	transposase activity转座酶活性	分子功能	1	35	8.00E-10
0016210	naringenin-chalcone synthase activity查耳酮合成酶活性	分子功能	21	38	1.90E-9
0006357	regulation of transcription from RNA polymerase II promoter RNA聚合酶II启动子转录调控	生物过程	11	817	4.60E-9
0055114	oxidation-reduction process氧化还原过程	生物过程	511	7139	6.60E-9

GO ID	GO 条目	功能分类	上调差异表达基因数目	注释的基因数目	KS值
0005524	ATP binding ATP结合	分子功能	474	10796	1.80E-26
0004364	glutathione transferase activity谷胱甘肽转移酶活性	分子功能	71	465	2.10E-17
0006749	glutathione metabolic process谷胱甘肽代谢过程	生物过程	80	496	1.20E-16
0004386	helicase activity解旋酶活性	分子功能	9	862	2.80E-12
0009407	toxin catabolic process毒素分解过程	生物过程	57	190	3.80E-12
0005737	cytoplasm细胞质	细胞组分	919	34128	6.00E-12
0030244	cellulose biosynthetic process纤维素生物合成过程	生物过程	9	193	1.70E-11
0004854	xanthine dehydrogenase activity黄嘌呤脱氢酶活性	分子功能	13	43	4.90E-11
0009115	xanthine catabolic process黄嘌呤分解过程	生物过程	13	43	5.70E-11
0004459	L-lactate dehydrogenase activity L-乳酸脱氢酶活性	分子功能	16	20	8.50E-11
0016760	cellulose synthase (UDP-forming) activity纤维素合成酶活性	分子功能	9	174	3.20E-10
0004803	transposase activity转座酶活性	分子功能	1	35	8.00E-10
0016210	naringenin-chalcone synthase activity查耳酮合成酶活性	分子功能	21	38	1.90E-9
0006357	regulation of transcription from RNA polymerase II promoter RNA聚合酶II启动子转录调控	生物过程	11	817	4.60E-9
0055114	oxidation-reduction process氧化还原过程	生物过程	511	7139	6.60E-9

KO号	代谢途径	差异基因数量	Q值
上调基因
ko00940	Phenylpropanoid biosynthesis苯丙烷生物合成	80	0.00E+00
ko00941	Flavonoid biosynthesis类黄酮生物合成	51	0.00E+00
ko00480	Glutathione metabolism谷胱甘肽代谢	88	0.00E+00
ko00960	Tropane, piperidine and pyridine alkaloid biosynthesis 莨菪碱、哌啶和吡啶类生物碱的生物合成	30	0.00E+00
ko00908	Zeatin biosynthesis玉米素生物合成	22	2.74E-11
ko00592	Alpha-linolenic acid metabolism α-亚油酸代谢	29	4.28E-09
ko00910	Nitrogen metabolism氮代谢	22	2.08E-06
ko00360	Phenylalanine metabolism苯丙氨酸代谢	26	2.43E-06
ko00640	Propanoate metabolism氟吡甲禾灵代谢	19	4.68E-05
ko00944	Flavone and flavonol biosynthesis黄酮和黄酮醇生物合成	10	5.00E-05
下调基因
ko00950	Isoquinoline alkaloid biosynthesis异喹啉生物碱生物合成	16	1.67E-08
ko00073	Cutin, suberine and wax biosynthesis角质素、亚黄素和蜡的生物合成	16	6.05E-08
ko00909	Sesquiterpenoid and triterpenoid biosynthesis倍半萜和三萜生物合成	12	9.77E-08
ko00410	Beta-alanine metabolism β-丙氨酸代谢	16	4.81E-06
ko00904	Diterpenoid biosynthesis双萜生物合成	11	9.74E-06
ko00900	Terpenoid backbone biosynthesis萜类生物合成	18	4.16E-05
ko00350	Tyrosine metabolism酪氨酸代谢	17	4.27E-05
ko00730	Thiamine metabolism硫胺素代谢	8	2.97E-04
ko00591	Linoleic acid metabolism亚油酸代谢	9	5.03E-04
ko00564	Glycerophospholipid metabolism甘油磷脂代谢	25	6.18E-04