棉花抗黄萎病共表达网络构建与关键基因筛选
安阳工学院, 安阳市生物信息学重点实验室, 河南安阳 455000
| 摘 要： | 大丽轮枝菌引起的棉花黄萎病是一种危害非常严重的土传维管束病害，对棉花产业造成了很大的影响。加权基因共表达网络分析是一种利用基因在不同条件下的表达数据对基因进行聚类的数据挖掘方法，该方法可以将性状和基因进行关联，是一种广泛应用的系统生物学技术。本项目使用具有抗黄萎病特性的瑟伯氏棉、具有高抗黄萎病特性的三裂棉和对黄萎病敏感的雷蒙德棉的根、茎、叶感染大丽轮枝菌后的RNAseq数据，瑟伯氏棉、三裂棉和雷蒙德棉在冷、盐处理12小时候的转录数据，以及雷蒙德棉开花后3天胚珠发育的转录数据构建基因共表达网络并筛选核心基因。分析结果获得1个与抗黄萎病正相关的模块和1个与抗黄萎病负相关的模块。基因功能富集分析结果表明，抗黄萎病性状正相关模块基因主要涉及有机物分解代谢过程、细胞壁组织或生物合成、外部封装结构组织和激素水平的调控等代谢过程。本研究的结果将对棉花黄萎病抗性基因的筛选有重要参考价值。 |
| 关 键 词： | 黄萎病; 棉花; WGCNA; 模块; 聚类 |
| DOI： | 10.57237/j.life.2023.04.002 |
View:
Downloads:
收稿日期: 2023-10-02
接受日期: 2023-11-09
在线出版日期: 2023-11-15
Construction of Cotton Verticillium wilt Resistance Co-expression Network and Identification of Hub Genes
Anyang Key Laboratory of Bioinformatics, Anyang Institute of Technology, Anyang 455000, China
| Abstract: | Cotton verticillium wilt caused by Verticillium dahliae is a serious soil transmission vascular bundle disease, which has a great impact on the yield and quality of cotton. Weighted gene co-expression network analysis (WGCNA) is a data mining method that uses gene expression data under different conditions to cluster genes. This method can associate traits with genes, and it is a widely used system biology technology. The transcriptional data of verticillium dahliae infected root, stem and leaf of 0h, 12h, 48h of Gossypium thurberi with verticillium wilt resistance, Gossypium trilobum with high verticillium wilt resistance and Gossypium raimondii with verticillium wilt sensitivity were used in this project, transcriptional data of Gossypium trilobum, Gossypium trilobum and Gossypium raimondii under cold and salt treatment for 12 hours respectively and transcriptional data of ovule development 3 days after anthesis in Gossypium raimondii were used to construct the gene co expression network and select hub genes. The results showed that one module was positively correlated with verticillium wilt resistance and the other was negatively correlated with verticillium wilt resistance. The results of gene function enrichment analysis showed that the positive correlation module genes of verticillium wilt resistance mainly involved in organic matter catabolism, cell wall tissue or biosynthesis, external encapsulation structure and hormone level regulation. The results of this study have important value to find verticillium wilt resistance genes in cotton. |
| Keywords: | Vrticillium Wilt; Cotton; WGCNA; Modular; Cluster |
参考文献
| [1] | 朱荷琴, 冯自力, 李志芳, 等. 分离自棉花的轮枝菌“种”的鉴定[J]. 中国农业科学, 2013, 46(10): 2032-2040. |
| [2] | 张绪振, 张树琴, 陈吉棣, 等. 我国棉花黄萎病菌“种”的鉴定 [J]. 植物病理学报, 1981(03): 15-20. |
| [3] | 丁晓华. 棉花枯、黄萎病的症状及防治方法 [J]. 河北农业, 2018(07): 30-31. |
| [4] | 伊静. 棉花黄萎病的发生与防治技术 [J]. 现代农业研究, 2018(06): 87-88. |
| [5] | 王坤波, 刘旭. 棉属多倍化研究进展 [J]. 中国农业科技导报, 2013, 15(02): 20-27. |
| [6] | Cai Y, Cai X, Wang Q, et al. Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis [J]. Plant Biotechnol J, 2020, 18(3): 814-828. |
| [7] | Wang K, Wang Z, Li F, et al. The draft genome of a diploid cotton Gossypium raimondii [J]. Nat Genet, 2012, 44(10): 1098-1103. |
| [8] | Hu Y, Chen J, Fang L, et al. Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton [J]. Nat Genet, 2019, 51(4): 739-748. |
| [9] | Li F, Fan G, Lu C, et al. Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution [J]. Nat Biotechnol, 2015, 33(5): 524-530. |
| [10] | Chen Z J, Sreedasyam A, Ando A, et al. Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement [J]. Nat Genet, 2020, 52(5): 525-533. |
| [11] | Li F, Fan G, Wang K, et al. Genome sequence of the cultivated cotton Gossypium arboretum [J]. Nat Genet, 2014, 46(6): 567-572. |
| [12] | 顾本康, 李经仪, 顾萍, 等. 棉属野生种枯萎病黄萎病抗性鉴定初报 [J]. 江苏农业科学, 1993(05): 36-37. |
| [13] | 梁理民, 刘有良, 王增信, 等. 陆地棉×斯特提棉种间杂交创造抗枯黄萎病新种质 [J]. 西北农业学报, 2002(04): 16-18. |
| [14] | 赵凤轩, 戴小枫. 棉花黄萎病菌的侵染过程 [J]. 基因组学与应用生物学, 2009, 28(04): 786-792. |
| [15] | Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis [J]. Stat Appl Genet Mol Biol, 2005, 4: e17. |
| [16] | Bolger A M, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data [J]. Bioinformatics, 2014, 30(15): 2114-2120. |
| [17] | Kim D, Langmead B, Salzberg S L. HISAT: a fast spliced aligner with low memory requirements [J]. Nat Methods, 2015, 12(4): 357-360. |
| [18] | Ghosh S, Chan C K. Analysis of RNA-Seq Data Using TopHat and Cufflinks [J]. Methods Mol Biol, 2016, 1374: 339-361. |
| [19] | Pollier J, Rombauts S, Goossens A. Analysis of RNA-Seq data with TopHat and Cufflinks for genome-wide expression analysis of jasmonate-treated plants and plant cultures [J]. Methods Mol Biol, 2013, 1011: 305-315. |
| [20] | Liao Y, Smyth G K, Shi W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads [J]. Nucleic Acids Res, 2019, 47(8): e47. |
| [21] | Robinson M D, McCarthy D J, Smyth G K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data [J]. Bioinformatics, 2010, 26(1): 139-140. |
| [22] | Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis [J]. BMC Bioinformatics, 2008, 9: 559. |
| [23] | Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks [J]. Genome Res, 2003, 13(11): 2498-2504. |
| [24] | Chin C H, Chen S H, Wu H H, et al. cytoHubba: identifying hub objects and sub-networks from complex interactome [J]. BMC Syst Biol, 2014, 8 Suppl 4(Suppl 4): S11. |
| [25] | Wu T, Hu E, Xu S, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data [J]. Innovation (Camb), 2021, 2(3): 100141. |
| [26] | Yu G, Wang L G, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters [J]. OMICS, 2012, 16(5): 284-287. |
| [27] | Du H, Zhang L, Liu L, et al. Biochemical and molecular characterization of plant MYB transcription factor family [J]. Biochemistry (Mosc), 2009, 74(1): 1-11. |
| [28] | Millard P S, Kragelund B B, Burow M. R2R3 MYB Transcription Factors - Functions outside the DNA-Binding Domain [J]. Trends Plant Sci, 2019, 24(10): 934-946. |
| [29] | Devoto A, Piffanelli P, Nilsson I, et al. Topology, subcellular localization, and sequence diversity of the Mlo family in plants [J]. J Biol Chem, 1999, 274(49): 34993-35004. |
