Abstract：Abstract Soybean cyst nematode (SCN, Heterodera glycines) seriously endangers the yield and quality of crops. SCN as one plant-specific endoparasitic nematode, is the most destructive pathogen of soybean (Glycine max) and other legumes, while tobacco (Nicotiana tabacum) is its non-host. However, in recent years, our research team found a special SCN population (SCNT) in Shandong province, which could infect tobacco, but whose pathogenicity to soybean was very poor. In order to reveal the molecular mechanism of significant difference in pathogenicity of the same host infected by SCNT and the normal SCN population, the Illumina platform HiSeqTM 2500 high-throughput sequencing technique was conducted to detect transcriptional profile of the 2nd Stage Juvenile (second stage juvenile, J2), which were derived from SCN and SCNT, respectively. qRT-PCR was performed to validate the result of transcriptome sequencing. The result demonstrated that there were a total of 1 628 differentially expressed genes (DEGs) between SCN and SCNT, among which, 1 347 DEGs were up-regulated in SCNT with 281 DEGs being down-regulated. GO (Gene Ontology) analysis was carried out to classify functions of DEGs. In the category of biological processes, the relative up-regulated and down-regulated genes were significantly enriched in the translation process, including 139 and 72 DEGs respectively. The two GO terms under cell component representing cytosolic small ribosomal subunit and the cytosolic large ribosomal subunit were notably enriched in the relative up-regulated genes, while other apparent GO terms in the relative down-regulated genes associated with DEGs were ribosome and cytosolic small ribosomal subunit. A notable number of genes (204 up-regulated and 80 down-regulated) were categorized under the GO term of structural constituent of ribosome belonging to molecular function. Furthermore, gene expression patterns were integrated with biochemical pathway from KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis. 639 DEGs might involve in 277 different metabolic pathway branches, among which, a significant number of genes participating ribosome pathway were enriched, including 210 relative up-regulated and 73 relative down-regulated genes. Other pathway of oxidative phosphorylation was notably enriched. 76 DEGs (70 up-regulated and 6 down-regulated) were categorized under phosphorylation. 6 down-regulated genes were localized to cytochrome b gene, cytochrome c oxidase CO I, CO Ⅱ, CO Ⅲ subunit, NADH oxidoreductase ND5 subunit and ATP enzyme subunit, while ND1, ND2, ND3, and subunits involved in NADH oxidoreductase were up regulated. In addition, pathways associated with nematode growth and development and parasitism (such as Nicotinate and nicotinamide metabolism, Glutathione metabolism and Metabolism of xenobiotics by cytochrome P450) were enriched, followed by 13 DEGs encoding esophageal gland cell secretory proteins. There were only two genes encoding secretory proteins (pel2 and Hgg-20) suppressed, while the other 11 genes performed different degrees of down-regulation in SCNT population. The qRT-PCR data were correlated with transcriptome sequencing result, since the correlation coefficient R2=0.98, which indicated transcriptome sequencing data was reliable. This is the first report describing differences on transcriptional levels of different pathogenic SCN and SCNT populations, which will provide a scientific basis for further analysis on molecular mechanisms of pathogenicity differences of SCN.
[1]程子超, 赵洪海, 李建立, 等.山东省寄生烟草的孢囊线虫种类鉴定及种内群体间rDNA-ITS-RFLP分析[J].植物病理学报, 2012, 42(4):387-395[2]郝大程, 马培, 穆军, 等.中药植物虎杖根的高通量转录组测序及转录组特性分析中国科学[J].中国科学, 2012, 42(5):398-412[3]扈丽丽, 卓侃, 林柏荣, 等.植物寄生线虫效应蛋白功能分析方法的研究进展[J].中国生物工程杂志, 2016, 6(02):101-108[4]李海燕, 王芳, 段玉玺, 等.大豆胞囊线虫侵染诱导五寨黑豆早期的转录组分析[J].中国油料作物学报, 2015, 37(02):194-200[5]黎娟华, 赵平娟, 孙海彦, 等.南方根结线虫延长因子2基因cDNA全长克隆和序列分析[J].中国农学通报, 2011, 27(01):223-228[6]刘海勤, 刘泉成, 赵洪海.大豆孢囊线虫大豆和烟草群体寄生性比较分析[J].植物保护, 2016, 42(2):68-73[7]刘维志.2004. 中国检疫性植物线虫[M]. 中国农业科学技术出版社, 北京. pp. 238-239. (Liu W Z. 2004. Quarantine plant nematode in China[M]. China agricultural science and technology press, Bejing. pp. 238-239.)[8]彭德良, 郑经武, 廖金玲, 等.2006. 重要植物线虫致病相关基因研究进展[C]. 中国植物病理学会2006年学术年会论文集. 中国农业科学技术出版社, 1: 221-230. (Peng D L, Zheng J W, Liao J L, et al. 2006. Progress in the research of pathogenic genes of important plant nematodes[C]. Proceedings of the Annual Meeting of Chinese Society for Plant Pathology (2006). China agricultural science and technology press, 1: 221-230.)[9]王慧敏, 赵洪海, 王凤龙.大豆孢囊线虫不同寄主种群的形态变异及对大豆寄主的适应性差异分析[J].植物保护, 2014, 40(6):106-111[10]吴书峰, 王志, 乔治军, 等.大豆抗胞囊线虫4号小种转录组测序及分析[J].农业生物技术学报, 2016, 24(04):519-529[11]杨丹.2017. 几种重要植物线虫高通量测序数据比较分析及靶标基因挖掘[D]. 博士学位论文, 中国农业大学, 导师: 简恒. pp: 37-41. (Yang D. 2017. Comparative analysis of high-throughput sequencing datasets target genes mining for several important plant-nematodes[D]. Thesis for Ph.D., Chinese Agricultural University, Supervisor: Jian H, pp: 37-41.)[12]赵洪海, 程子超, 王凤龙.山东省烟草孢囊线虫的群体动态和世代发生特点[J].植物保护学报, 2013, 40(6):529-532[13]赵洪海, 丁海燕, 王凤龙.大豆和烟草上大豆孢囊线虫田间侵染特征比较分析[J].中国农业科学, 2014, 47(22):4417-4425[14]Andersen J S, Lam Y W, Leung A K, et al.Nucleolar proteome dynamics[J].Nature, 2005, 433(7021):77-83[15]Bakhetia M, Urwin P E, Atkinson H J.Characterisation by RNAi of pioneer genes expressed in the dorsal pharyngeal gland cell of Heterodera glycines and the effects of combinatorial RNAi[J].International journal for parasitology, 2008, 38(13):1589-1597[16]Fraser A G, Kamath R S, Zipperlen P, et al.Functional genomic analysis of C.elegans chromosome I by systematic RNA interference[J].Nature, 2000, 408(6810):325-330[17]Gao B L, Allen R, Davis E L, et al.Developmental expression and biochemical properties of a beta-1, 4-endoglucanase family in the soybean cyst nematode, Heterodera glycines[J].Molecular plant pathology, 2004b, 5(2):93-104[18]Gao B L, Allen R, Maier T, et al.The parasitome of the phytonematode Heterodera glycines[J].Molecular plant-microbe internactions, 2003, 16(8):720-726[19]Gao B L, Allen R, Maier T, et al.Identification of a new beta-1, 4-endoglucanase gene expressed in the esophageal subventral gland cells of Heterodera glycines[J].Journal of nematology, 2002, 34(1):12-15[20]Hamamouch N, Li C Y, Hewezi T, et al.The interaction of the novel 30C02 cyst nematode effector protein with a plant β-1, 3-endoglucanase may suppress host defence to promote parasitism[J].Journal of Experimental Botany, 2012, 63(10):3683-3695[21]Ithal N, Recknor J, Nettleton D, et al.Developmental transcript profiling of cyst nematode feeding cells in soybean roots[J].Molecular Plant-Microbe Interactions, 2007, 20(5):510-525[22]Kapp L D, Lorsch J R.The molecular mechanics of eukaryotic translation[J].Annu Rev Biochem, 2004, 73(1):657-704[23]Mardis E R.Next-generation DNA sequencing methods[J].Annu Rev Genomics Hum Genet, 2008, 9(1):387-402[24]Mason C E, Elemento O.Faster sequencers, larger datasets, new challenges[J].Genomebiology, 2012 , 13(3):314-314[25]Niblack T L, Arelli P R, Noel G R, et al.A revised classification scheme for genetically diverse populations of Heterodera glycines[J].Journal Nematology, 2002, 34(4):279-288[26]Niblack T L, Colgrove A L.Adaptability of Heterodera glycines to new hosts and resistant soybean cultivars[J].Phytopathology, 2006, 96(6):S84-84[27]Noel G R.Heterodera glycines in soybean[J].Nematologia Brasileira, 1993, 17(2):103-121[28]Perkd J.Finding the true $1000 genome[J].BioTechniques, 2013, 54(2):71-74[29]Riggs R D, Schmitt D P.Complete characterization of the race scheme for Heterodera glycines[J].Journal of nematology, 1988, 20(3):392-395[30]Rinke C, Schwientek P, Sczyrba A, et al.Insights into the phylogeny and coding potential of microbial dark matter[J].Nature, 2013, 499(7459):431-437[31]Shi H, zheng J.First report of soybean cyst nematode (Heterodera glycines) on tobacco in Henan, central China[J].Plant Disease, 2013, 97(6):852-852[32]Subbotin SA, Cid Del Prado Vera I, Mundo-Ocampo M, et al.Identification, phylogeny and phylogeography of circumfenestrate cyst nematodes (Nematoda: Heteroderidae) as inferred from analysis of ITS-rDNA.y[J].Nematolog, 2011, 13(1):805-824[33]Tu L C, Yan X, Hood L, et al.Proteomics analysis of the interactome of N-myc downstream regulated gene 1 and its interactions with the androgen response program in prostate cancer cells[J].Mol Cell Proteomics, 2007, 6(4):575-588[34]Tylka G L, Niblack T L.2010. Soybean cyst nematode management guide[M]. Plant Health Initiative, http://www.manitobapulse.ca/wp-content/uploads/SCN_Management.pdf. pp: 1-16.[35]Wang X H, Mitchum M G, Gao B L, et al.A parasitism gene from a plant-parasitic nematode with function similar to CLAVATA3/ESR(CLE) of Arabidopsis thaliana.[J].Molecular plant pathology(2): 187-191., 2005, 6(2):187-191[36]Yan Y T, Smant G, Davis E.Functional screening yields a new beta-1, 4-endoglucanase gene from Heterodera glycines that may be the product of recent gene duplication[J].Molecular plant-microbe internactions, 2001, 14(1):63-71[37]Yu Qing.2011. Soybean cyst nematode (Heterodera glycines Ichinohe) [M]//EI-Shemy H A. Soybean Physiology and Biochemistry. Rijeka, Croatia:InTech Europe, pp:461-474.
