Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.)

Abstract

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper–pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

References

[1] Lee HA, Yeom SI . Plant NB-LRR proteins: tightly regulated sensors in a complex manner. Brief Funct Genomics. 2015; 14: 233-42.

[2] Wan WL, Frohlich K, Pruitt RN et al. Plant cell surface immune receptor complex signaling. Curr Opin Plant Biol. 2019; 50: 18-28.

[3] Albert I, Hua C, Nurnberger T et al. Surface sensor systems in plant immunity. Plant Physiol. 2020; 182: 1582-96.

[4] Du J, Verzaux E, Chaparro-Garcia A et al. Elicitin recognition confers enhanced resistance to Phytophthora infestans in potato. Nat Plants. 2015; 1: 15034.

[5] Lacombe S, Rougon-Cardoso A, Sherwood E et al. Interfamily transfer of a plant pattern-recognition receptor confers broad-spectrum bacterial resistance. Nat Biotechnol. 2010; 28: 365-9.

[6] Dixon MS, Hatzixanthis K, Jones DA et al. The tomato Cf-5 disease resistance gene and six homologs show pronounced allelic variation in leucine-rich repeat copy number. Plant Cell. 1998; 10: 1915-25.

[7] Dixon MS, Jones DA, Keddie JS et al. The tomato Cf-2 disease resistance locus comprises two functional genes encoding leucine-rich repeat proteins. Cell. 1996; 84: 451-9.

[8] Takken FL, Schipper D, Nijkamp HJ et al. Identification and ds-tagged isolation of a new gene at the Cf-4 locus of tomato involved in disease resistance to Cladosporium fulvum race 5. Plant J. 1998; 14: 401-11.

[9] Thomas CM, Jones DA, Parniske M et al. Characterization of the tomato Cf-4 gene for resistance to Cladosporium fulvum identifies sequences that determine recognitional specificity in Cf-4 and Cf-9. Plant Cell. 1997; 9: 2209-24.

[10] Hegenauer V, Fürst U, Kaiser B et al. Detection of the plant parasite Cuscuta reflexa by a tomato cell surface receptor. Science. 2016; 353: 478-81.

[11] Jehle AK, Lipschis M, Albert M et al. The receptor-like protein ReMAX of Arabidopsis detects the microbe-associated molecular pattern eMax from Xanthomonas. Plant Cell. 2013; 25: 2330-40.

[12] Kawchuk LM, Hachey J, Lynch DR et al. Tomato Ve disease resistance genes encode cell surface-like receptors. Proc Natl Acad Sci U S A. 2001; 98: 6511-5.

[13] Larkan NJ, Lydiate DJ, Parkin IAP et al. The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1. New Phytol. 2013; 197: 595-605.

[14] Ron M, Avni A . The receptor for the fungal elicitor ethylene-inducing xylanase is a member of a resistance-like gene family in tomato. Plant Cell. 2004; 16: 1604-15.

[15] Shen Y, Diener AC . Arabidopsis thaliana resistance to fusarium oxysporum 2 implicates tyrosine-sulfated peptide signaling in susceptibility and resistance to root infection. PLoS Genet. 2013; 9: e1003525.

[16] Steinbrenner AD, Muñoz-Amatriaín M, Chaparro AF et al. A receptor-like protein mediates plant immune responses to herbivore-associated molecular patterns. Proc Natl Acad Sci U S A. 2020; 117: 31510-8.

[17] Vinatzer BA, Patocchi A, Gianfranceschi L et al. Apple contains receptor-like genes homologous to the Cladosporium fulvum resistance gene family of tomato with a cluster of genes cosegregating with Vf apple scab resistance. Mol Plant-Microbe Interact. 2001; 14: 508-15.

[18] Zhang B, Li P, Su T et al. BrRLP48, encoding a receptor-like protein, involved in downy mildew resistance in Brassica rapa. Front Plant Sci. 2018; 9: 1708.

[19] Liebrand TW, van den Burg HA, Joosten MH . Two for all: receptor-associated kinases SOBIR1 and BAK1. Trends Plant Sci. 2014; 19: 123-32.

[20] Jeong S, Trotochaud AE, Clark SE . The Arabidopsis CLAVATA2 gene encodes a receptor-like protein required for the stability of the CLAVATA1 receptor-like kinase. Plant Cell. 1999; 11: 1925-34.

[21] Nadeau JA, Sack FD . Control of stomatal distribution on the Arabidopsis leaf surface. Science. 2002; 296: 1697-700.

[22] Taguchi-Shiobara F, Yuan Z, Hake S et al. The fasciated ear2 gene encodes a leucine-rich repeat receptor-like protein that regulates shoot meristem proliferation in maize. Genes Dev. 2001; 15: 2755-66.

[23] Kang WH, Park B, Lee J et al. Tissue-specific RNA-seq analysis and identification of receptor-like proteins related to plant growth in Capsicum annuum. Plan Theory. 2021; 10: 972.

[24] Wang G, Ellendorff U, Kemp B et al. A genome-wide functional investigation into the roles of receptor-like proteins in Arabidopsis. Plant Physiol. 2008; 147: 503-17.

[25] Petre B, Hacquard S, Duplessis S et al. Genome analysis of poplar LRR-RLP gene clusters reveals RISP, a defense-related gene coding a candidate endogenous peptide elicitor. Front Plant Sci. 2014; 5: 111.

[26] Kang WH, Yeom SI . Genome-wide identification, classification, and expression analysis of the receptor-like protein family in tomato. Plant Pathol J. 2018; 34: 435-44.

[27] Fritz-Laylin LK, Krishnamurthy N, Tor M et al. Phylogenomic analysis of the receptor-like proteins of rice and Arabidopsis. Plant Physiol. 2005; 138: 611-23.

[28] Yang H, Bayer PE, Tirnaz S et al. Genome-wide identification and evolution of receptor-like kinases (RLKs) and receptor like proteins (RLPs) in Brassica juncea. Biology. 2020; 10: 17.

[29] Restrepo-Montoya D, Brueggeman R, McClean PE et al. Computational identification of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes. BMC Genomics. 2020; 21: 459.

[30] Li Y, Pearl SA, Jackson SA . Gene networks in plant biology: approaches in reconstruction and analysis. Trends Plant Sci. 2015; 20: 664-75.

[31] Sun BM, Zhou X, Chen C et al. Coexpression network analysis reveals an MYB transcriptional activator involved in capsaicinoid biosynthesis in hot peppers. Hortic Res. 2020; 7: 162.

[32] Zheng J, He C, Qin Y et al. Co-expression analysis aids in the identification of genes in the cuticular wax pathway in maize. Plant J. 2019; 97: 530-42.

[33] Ruprecht C, Vaid N, Proost S et al. Beyond genomics: studying evolution with gene coexpression networks. Trends Plant Sci. 2017; 22: 298-307.

[34] Amrine KC, Blanco-Ulate B, Cantu D . Discovery of core biotic stress responsive genes in Arabidopsis by weighted gene co-expression network analysis. PLoS One. 2015; 10: e0118731.

[35] Dossa K, Mmadi MA, Zhou R et al. Depicting the core transcriptome modulating multiple abiotic stresses responses in sesame (Sesamum indicum L.). Int J Mol Sci. 2019; 20: 3930.

[36] Zhu M, Xie H, Wei X et al. WGCNA analysis of salt-responsive core transcriptome identifies novel hub genes in rice. Genes. 2019; 10: 719.

[37] Ruprecht C, Mendrinna A, Tohge T et al. FamNet: a framework to identify multiplied modules driving pathway expansion in plants. Plant Physiol. 2016; 170: 1878-94.

[38] Ichihashi Y, Aguilar-Martínez JA, Farhi M et al. Evolutionary developmental transcriptomics reveals a gene network module regulating interspecific diversity in plant leaf shape. Proc Natl Acad Sci U S A. 2014; 111: E2616-21.

[39] Chang YM, Lin HH, Liu WY et al. Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes. Proc Natl Acad Sci U S A. 2019; 116: 3091-9.

[40] Gerstein MB, Rozowsky J, Yan KK et al. Comparative analysis of the transcriptome across distant species. Nature. 2014; 512: 445-8.

[41] Mutwil M, Klie S, Tohge T et al. PlaNet: combined sequence and expression comparisons across plant networks derived from seven species. Plant Cell. 2011; 23: 895-910.

[42] Kang WH, Sim YM, Koo N et al. Transcriptome profiling of abiotic responses to heat, cold, salt, and osmotic stress of Capsicum annuum L. Sci Data. 2020; 7: 17.

[43] Lee J, Nam JY, Jang H et al. Comprehensive transcriptome resource for response to phytohormone-induced signaling in Capsicum annuum L. BMC Res Notes. 2020; 13: 440.

[44] Kim S, Park J, Yeom SI et al. New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication. Genome Biol. 2017; 18: 210.

[45] Kim MS, Kim S, Jeon J et al. Global gene expression profiling for fruit organs and pathogen infections in the pepper, Capsicum annuum L. Sci Data. 2018; 5: 180103.

[46] Seo E, Kim S, Yeom SI et al. Genome-wide comparative analyses reveal the dynamic evolution of nucleotide-binding leucine-rich repeat gene family among Solanaceae plants. Front Plant Sci. 2016; 7: 1205.

[47] Kang WH, Kim S, Lee HA et al. Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper. Sci Rep. 2016; 6: 33332.

[48] Kim N, Kang WH, Lee J et al. Development of clustered resistance gene analogs-based markers of resistance to Phytophthora capsici in chili pepper. Biomed Res Int. 2019; 2019: 1093186.

[49] Yeom SI, Baek HK, Oh SK et al. Use of a secretion trap screen in pepper following Phytophthora capsici infection reveals novel functions of secreted plant proteins in modulating cell death. Mol Plant-Microbe Interact. 2011; 24: 671-84.

[50] Kim SB, Kang WH, Huy HN et al. Divergent evolution of multiple virus-resistance genes from a progenitor in capsicum spp. New Phytol. 2017; 213: 886-99.

[51] Frederickson Matika DE, Loake GJ . Redox regulation in plant immune function. Antioxid Redox Signal. 2014; 21: 1373-88.

[52] Mittler R . Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002; 7: 405-10.

[53] Wu L, Wang S, Wu J et al. Phosphoproteomic analysis of the resistant and susceptible genotypes of maize infected with sugarcane mosaic virus. Amino Acids. 2015; 47: 483-96.

[54] Saijo Y, Loo EP, Yasuda S . Pattern recognition receptors and signaling in plant-microbe interactions. Plant J. 2018; 93: 592-613.

[55] Kazan K, Gardiner DM . Targeting pathogen sterols: defence and counterdefence? PLoS Pathog. 2017; 13: e1006297.

[56] van Loon LC, Geraats BP, Linthorst HJ . Ethylene as a modulator of disease resistance in plants. Trends Plant Sci. 2006; 11: 184-91.

[57] Erb M, Kliebenstein DJ . Plant secondary metabolites as defenses, regulators, and primary metabolites: the blurred functional trichotomy. Plant Physiol. 2020; 184: 39-52.

[58] Yi SY, Lee DJ, Yeom SI et al. A novel pepper (Capsicum annuum) receptor-like kinase functions as a negative regulator of plant cell death via accumulation of superoxide anions. New Phytol. 2010; 185: 701-15.

[59] Kwon JS, Nam JY, Yeom SI et al. Leaf-to-whole plant spread bioassay for pepper and Ralstonia solanacearum interaction determines inheritance of resistance to bacterial wilt for further breeding. Int J Mol Sci. 2021; 22: 2279.

[60] Andolfo G, Sanseverino W, Rombauts S et al. Overview of tomato (Solanum lycopersicum) candidate pathogen recognition genes reveals important solanum R locus dynamics. New Phytol. 2013; 197: 223-37.

[61] Kruijt M, DE Kock MJD, de Wit. Receptor-like proteins involved in plant disease resistance. Mol Plant Pathol. 2005; 6: 85-97.

[62] Buscaill P, Rivas S . Transcriptional control of plant defence responses. Curr Opin Plant Biol. 2014; 20: 35-46.

[63] Eulgem T. Regulation of the Arabidopsis defense transcriptome. Trends Plant Sci. 2005; 10: 71-8.

[64] Tsuda K, Somssich IE . Transcriptional networks in plant immunity. New Phytol. 2015; 206: 932-47.

[65] Javed T, Shabbir R, Ali A et al. Transcription factors in plant stress responses: challenges and potential for sugarcane improvement. Plan Theory. 2020; 9: 491.

[66] Yang Y, Zhang Y, Ding P et al. The ankyrin-repeat transmembrane protein BDA1 functions downstream of the receptor-like protein SNC2 to regulate plant immunity. Plant Physiol. 2012; 159: 1857-65.

[67] Zhang Y, Yang Y, Fang B et al. Arabidopsis snc2-1D activates receptor-like protein-mediated immunity transduced through WRKY70. Plant Cell. 2010; 22: 3153-63.

[68] Kothari SL, Joshi A, Kachhwaha S et al. Chilli peppers-a review on tissue culture and transgenesis. Biotechnol Adv. 2010; 28: 35-48.

[69] Woo JY, Jeong KJ, Kim YJ et al. CaLecRK-S.5, a pepper L-type lectin receptor kinase gene, confers broad-spectrum resistance by activating priming. J Exp Bot. 2016; 67: 5725-41.

[70] Finn RD, Bateman A, Clements J et al. Pfam: the protein families database. Nucleic Acids Res. 2014; 42: D222-30.

[71] Schultz J, Milpetz F, Bork P et al. SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A. 1998; 95: 5857-64.

[72] Li L, Stoeckert CJ Jr, Roos DS . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003; 13: 2178-89.

[73] Voorrips RE . MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002; 93: 77-8.

[74] Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet journal. 2011; 17: 10.

[75] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[76] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[77] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 550.

[78] Gu Z, Eils R, Schlesner M . Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016; 32: 2847-9.

[79] Ma C, Xin M, Feldmann KA et al. Machine learning-based differential network analysis: a study of stress-responsive transcriptomes in Arabidopsis. Plant Cell. 2014; 26: 520-37.

[80] Shannon P, Markiel A, Ozier O et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13: 2498-504.

[81] Young MD, Wakefield MJ, Smyth GK et al. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010; 11: R14.

[82] Gotz S, García-Gómez JM, Terol J et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008; 36: 3420-35.

[83] Yeom SI, Seo E, Oh SK et al. A common plant cell-wall protein HyPRP1 has dual roles as a positive regulator of cell death and a negative regulator of basal defense against pathogens. Plant J. 2012; 69: 755-68.