Systems Biology Approaches Reveal a Multi-stress Responsive WRKY Transcription Factor and Stress Associated Gene Co-expression Networks in Chickpea

Author(s): Aravind K. Konda* , Parasappa R. Sabale , Khela R. Soren , Shanmugavadivel P. Subramaniam , Pallavi Singh , Santosh Rathod , Sushil K. Chaturvedi , Narendra P. Singh .

Journal Name: Current Bioinformatics

Volume 14 , Issue 7 , 2019

Become EABM
Become Reviewer

Graphical Abstract:


Abstract:

Background: Chickpea is a nutritional rich premier pulse crop but its production encounters setbacks due to various stresses and understanding of molecular mechanisms can be ascribed foremost importance.

Objective: The investigation was carried out to identify the differentially expressed WRKY TFs in chickpea in response to herbicide stress and decipher their interacting partners.

Methods: For this purpose, transcriptome wide identification of WRKY TFs in chickpea was done. Behavior of the differentially expressed TFs was compared between other stress conditions. Orthology based cofunctional gene networks were derived from Arabidopsis. Gene ontology and functional enrichment analysis was performed using Blast2GO and STRING software. Gene Coexpression Network (GCN) was constructed in chickpea using publicly available transcriptome data. Expression pattern of the identified gene network was studied in chickpea-Fusarium interactions.

Results: A unique WRKY TF (Ca_08086) was found to be significantly (q value = 0.02) upregulated not only under herbicide stress but also in other stresses. Co-functional network of 14 genes, namely Ca_08086, Ca_19657, Ca_01317, Ca_20172, Ca_12226, Ca_15326, Ca_04218, Ca_07256, Ca_14620, Ca_12474, Ca_11595, Ca_15291, Ca_11762 and Ca_03543 were identified. GCN revealed 95 hub genes based on the significant probability scores. Functional annotation indicated role in callose deposition and response to chitin. Interestingly, contrasting expression pattern of the 14 network genes was observed in wilt resistant and susceptible chickpea genotypes, infected with Fusarium.

Conclusion: This is the first report of identification of a multi-stress responsive WRKY TF and its associated GCN in chickpea.

Keywords: Chickpea, WRKY transcription factor, Gene Co-expression Networks, multi-stress, stress signaling pathways, Arabidopsis.

[1]
Camacho DM, Collins JJ. Systems biology strikes gold. Cell 2009; 137(1): 24-6.
[2]
Lee I, Ambaru B, Thakkar P, Marcotte EM, Rhee SY. Rational association of genes with traits using a genome-scale gene network for Arabidopsis thaliana. Nat Biotechnol 2010; 28(2): 149-56.
[3]
Jain M, Misra G, Patel RK, et al. A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.). Plant J 2013; 74(5): 715-29.
[4]
Varshney RK, Song C, Saxena RK, et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat Biotechnol 2013; 31(3): 240-6.
[5]
Grigoriev A. A relationship between gene expression and protein interactions on the proteome scale: analysis of the bacteriophage T7 and the yeast Saccharomyces cerevisiae. Nucleic Acids Res 2001; 29(17): 3513-9.
[6]
Carter SL, Brechbühler CM, Griffin M, Bond AT. Gene co-expression network topology provides a framework for molecular characterization of cellular state. Bioinformatics 2004; 20(14): 2242-50.
[7]
Lee T, Yang S, Kim E, et al. AraNet v2: an improved database of co-functional gene networks for the study of Arabidopsis tha-liana and 27 other nonmodel plant species. Nucleic Acids Res 2015; 43: D996-D1002.
[8]
Agarwal P, Reddy MP, Chikara J. WRKY: its structure, evolutionary relationship, DNA-binding selectivity, role in stress tolerance and development of plants. Mol Biol Rep 2011; 38(6): 3883-96.
[9]
Chen L, Song Y, Li S, Zhang L, Zou C, Yu D. The role of WRKY transcription factors in plant abiotic stresses. Biochim Biophys Acta 2012; 1819(2): 120-8.
[10]
Eulgem T, Somssich IE. Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol 2007; 10(4): 366-71.
[11]
Kumar K, Srivastava V, Purayannur S, Kaladhar VC, Cheruvu PJ, Verma PK. WRKY domain-encoding genes of a crop legume chickpea (Cicer arietinum): comparative analysis with Medicago truncatula WRKY family and characterization of group-III gene(s). DNA Res 2016; 23(3): 225-39.
[12]
Boratyn GM, Schäffer AA, Agarwala R, Altschul SF, Lipman DJ, Madden TL. Domain enhanced lookup time accelerated BLAST. Biol Direct 2012; 7(1): 12.
[13]
Iquebal MA, Soren KR, Gangwar P, et al. Discovery of putative herbicide re-sistance genes and its regulatory network in chickpea using transcriptome sequencing. Front Plant Sci 2017; 8: 958.
[14]
Li J, Dai X, Liu T, Zhao PX. LegumeIP: an integrative database for comparative genomics and transcriptomics of model legumes. Nucleic Acids Res 2012; 40(Database issue): D1221-9.
[15]
Das S, Meher PK, Pradhan UK, Paul AK. Inferring gene regu-latory networks using Kendall’s tau correlation coefficient and identification of salinity stress responsive genes in rice. Curr Sci 2017; 112(6): 1257.
[16]
Das S, Meher PK, Rai A, Bhar LM, Mandal BN. Statistical Approaches for Gene Selection, Hub Gene Identification and Module Interaction in Gene Co-Expression Network Analysis: An Application to Aluminum Stress in Soybean (Glycine max L.). PLoS One 2017; 12(1)e0169605
[17]
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005; 21(18): 3674-6.
[18]
Szklarczyk D, Franceschini A, Wyder S, et al. STRING v10: protein–protein interaction net-works, integrated over the tree of life. Nucleic Acids Res 2015; 43(Database issue): D447-52.
[19]
Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res 2011; 39(Suppl. 2).W155-9
[20]
Kohli D, Joshi G, Deokar AA, et al. Identification and characterization of Wilt and salt stress-responsive microRNAs in chickpea through high-throughput sequencing. PLoS One 2014; 9(10)e108851
[21]
Srivastava S, Zheng Y, Kudapa H, et al. High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea. Plant Sci 2015; 235: 46-57.
[22]
Garg R, Bhattacharjee A, Jain M. Genome-scale transcriptomic insights into molecular aspects of abiotic stress responses in chickpea. Plant Mol Biol Report 2015; 33(3): 388-400.
[23]
Gupta S, Bhar A, Chatterjee M, Ghosh A, Das S. Transcriptomic dissection reveals wide spread differential expression in chickpea during early time points of Fusarium oxysporum f. sp. ciceri Race 1 attack. PLoS One 2017; 12(5)e0178164
[24]
Romeis T, Ludwig AA, Martin R, Jones JD. Calcium-dependent protein kinases play an essential role in a plant defence response. EMBO J 2001; 20(20): 5556-67.
[25]
AbuQamar S, Chen X, Dhawan R, et al. Expression profiling and mutant analysis reveals complex regulatory networks involved in Arabidopsis response to Botrytis infection. Plant J 2006; 48(1): 28-44.
[26]
Cai Y, Zhuang X, Gao C, Wang X, Jiang L. The Arabidopsis endosomal sorting complex required for transport III regu-lates internal vesicle formation of the prevacuolar compart-ment and is required for plant development. Plant Physiol 2014; 165(3): 1328-43.
[27]
Katsiarimpa A, Kalinowska K, Anzenberger F, et al. The deubiquitinating enzyme AMSH1 and the ESCRT-III subunit VPS2.1 are required for autophagic degradation in Arabidopsis. Plant Cell 2013; 25(6): 2236-52.
[28]
Qiang X, Zechmann B, Reitz MU, Kogel KH, Schäfer P. The mutualistic fungus Piriformospora indica colonizes Arabidopsis roots by inducing an endoplasmic reticulum stress-triggered caspase-dependent cell death. Plant Cell 2012; 24(2): 794-809.
[29]
Park CJ, Peng Y, Chen X, et al. Rice XB15, a protein phosphatase 2C, negatively regulates cell death and XA21-mediated innate immunity. PLoS Biol 2008; 6(9)e231
[30]
Ross A, Yamada K, Hiruma K, et al. The Arabidopsis PEPR pathway couples local and systemic plant immunity. EMBO J 2014; 33(1): 62-75.
[31]
Cao J, Jiang M, Li P, Chu Z. Genome-wide identification and evolutionary analyses of the PP2C gene family with their expression profiling in response to multiple stresses in Brachypodium distachyon. BMC Genomics 2016; 17(1): 175.
[32]
Kunz BA, Dando PK, Grice DM, Mohr PG, Schenk PM, Cahill DM. UV-induced DNA damage promotes resistance to the biotrophic pathogen Hyaloperonospora parasitica in Arabidopsis. Plant Physiol 2008; 148(2): 1021-31.
[33]
Zeng Q, Sritubtim S, Ellis BE. AtMKK6 and AtMPK13 are required for lateral root formation in Arabidopsis. Plant Signal Behav 2011; 6(10): 1436-9.
[34]
Zhou C, Yin Y, Dam P, Xu Y. Identification of novel proteins involved in plant cell-wall synthesis based on protein-protein interaction data. J Proteome Res 2010; 9(10): 5025-37.
[35]
Li N, Gügel IL, Giavalisco P, et al. FAX1, a novel membrane protein mediating plastid fatty acid export. PLoS Biol 2015; 13(2)e1002053
[36]
Vanholme R, Cesarino I, Rataj K, et al. Caffeoyl shikimate esterase (CSE) is an enzyme in the lignin biosynthetic pathway in Arabidopsis. Science 2013; 341(6150): 1103-6.
[37]
Bhuiyan NH, Selvaraj G, Wei Y, King J. Gene expression profiling and silencing reveal that monolignol biosynthesis plays a critical role in penetration defence in wheat against powdery mildew invasion. J Exp Bot 2009; 60(2): 509-21.
[38]
Kumar Y, Zhang L, Panigrahi P, et al. Fusarium oxysporum mediates systems metabolic repro-gramming of chickpea roots as revealed by a combination of proteomics and metabolomics. Plant Biotechnol J 2016; 14(7): 1589-603.
[39]
Funnell DL, Pedersen JF. Reaction of sorghum lines genetical-ly modified for reduced lignin content to infection by Fusari-um and Alternaria spp. Plant Dis 2006; 90(3): 331-8.


Rights & PermissionsPrintExport Cite as

Article Details

VOLUME: 14
ISSUE: 7
Year: 2019
Page: [591 - 601]
Pages: 11
DOI: 10.2174/1574893614666190204152500
Price: $58

Article Metrics

PDF: 17
HTML: 1