Diversity and conservation of plant small secreted proteins associated with arbuscular mycorrhizal symbiosis

Abstract

Arbuscular mycorrhizal symbiosis (AMS) is widespread mutualistic association between plants and fungi, which plays an essential role in nutrient exchange, enhancement in plant stress resistance, development of host, and ecosystem sustainability. Previous studies have shown that plant small secreted proteins (SSPs) are involved in beneficial symbiotic interactions. However, the role of SSPs in the evolution of AMS has not been well studied yet. In this study, we performed computational analysis of SSPs in 60 plant species and identified three AMS-specific ortholog groups containing SSPs only from at least 30% of the AMS species in this study and three AMS-preferential ortholog groups containing SSPs from both AMS and non-AMS species, with AMS species containing significantly more SSPs than non-AMS species. We found that independent lineages of monocot and eudicot plants contained genes in the AMS-specific ortholog groups and had significant expansion in the AMS-preferential ortholog groups. Also, two AMS-preferential ortholog groups showed convergent changes, between monocot and eudicot species, in gene expression in response to arbuscular mycorrhizal fungus Rhizophagus irregularis. Furthermore, conserved cis-elements were identified in the promoter regions of the genes showing convergent gene expression. We found that the SSPs, and their closely related homologs, in each of three AMS-preferential ortholog groups, had some local variations in the protein structural alignment. We also identified genes co-expressed with the Populus trichocarpa SSP genes in the AMS-preferential ortholog groups. This first plant kingdom-wide analysis on SSP provides insights on plant-AMS convergent evolution with specific SSP gene expression and local diversification of protein structures.

References

[1] Lease KA, Walker JC . The Arabidopsis unannotated secreted peptide database, a resource for plant peptidomics. Plant Physiol. 2006; 142: 831-8.

[2] Tabata R, Sawa S . Maturation processes and structures of small secreted peptides in plants. Front Plant Sci. 2014; 5: 311.

[3] Tavormina P, De Coninck B, Nikonorova N et al. The plant peptidome: an expanding repertoire of structural features and biological functions. Plant Cell. 2015; 27: 2095-118.

[4] Chen YL, Fan KT, Hung SC, Chen YR . The role of peptides cleaved from protein precursors in eliciting plant stress reactions. New Phytol. 2020; 225: 2267-82.

[5] Fukuda H, Ohashi-Ito K . Vascular tissue development in plants. Curr Top Dev Biol. 2019; 131: 141-60.

[6] Araya T, Miyamoto M, Wibowo J et al. CLE-CLAVATA1 peptide-receptor signaling module regulates the expansion of plant root systems in a nitrogen-dependent manner. Proc Natl Acad Sci U S A. 2014; 111: 2029-34.

[7] Takahashi F, Suzuki T, Osakabe Y et al. A small peptide modulates stomatal control via abscisic acid in long-distance signalling. Nature. 2018; 556: 235-8.

[8] Mohd-Radzman NA, Binos S, Truong TT et al. Novel MtCEP1 peptides produced in vivo differentially regulate root development in Medicago truncatula . J Exp Bot. 2015; 66: 5289-300.

[9] Bonfante P, Genre A . Mechanisms underlying beneficial plant-fungus interactions in mycorrhizal symbiosis. Nat Commun. 2010; 1: 1-11.

[10] de Bang TC, Lundquist PK, Dai X et al. Genome-wide identification of Medicago peptides involved in macronutrient responses and nodulation. Plant Physiol. 2017; 175: 1669-89.

[11] Plett JM, Yin H, Mewalal R et al. Populus trichocarpa encodes small, effector-like secreted proteins that are highly induced during mutualistic symbiosis . Sci Rep. 2017; 7: 1-13.

[12] Peeters MK, Menschaert G. The hunt for sORFs: a multidisciplinary strategy. Exp Cell Res. 2020;391:111923.

[13] Hanada K, Akiyama K, Sakurai T et al. sORF finder: a program package to identify small open reading frames with high coding potential. Bioinformatics. 2010; 26: 399-400.

[14] Zhou P, Silverstein KA, Gao L et al. Detecting small plant peptides using SPADA (small peptide alignment discovery application). BMC Bioinformatics. 2013; 14: 1-16.

[15] Guo J-C, Fang S-S, Wu Y et al. CNIT: a fast and accurate web tool for identifying protein-coding and long non-coding transcripts based on intrinsic sequence composition. Nucleic Acids Res. 2019; 47: W516-22.

[16] Zhu M, Gribskov M . MiPepid: MicroPeptide identification tool using machine learning. BMC Bioinformatics. 2019; 20: 1-11.

[17] Zhang Y, Jia C, Fullwood MJ, Kwoh CK . DeepCPP: a deep neural network based on nucleotide bias information and minimum distribution similarity feature selection for RNA coding potential prediction. Brief Bioinform. 2020; 22: 2073-84.

[18] Armenteros JJA, Tsirigos KD, Sønderby CK et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019; 37: 420-3.

[19] Möller S, Croning MD, Apweiler R . Evaluation of methods for the prediction of membrane spanning regions. Bioinformatics. 2001; 17: 646-53.

[20] Nielsen H, Petsalaki EI, Zhao L, Stühler K . Predicting eukaryotic protein secretion without signals. Biochim Biophys Acta, Proteins Proteomics. 2019; 1867: 140174.

[21] Sperschneider J, Dodds PN, Singh KB, Taylor JM . ApoplastP: prediction of effectors and plant proteins in the apoplast using machine learning. New Phytol. 2018; 217: 1764-78.

[22] Savojardo C, Martelli PL, Fariselli P et al. BUSCA: an integrative web server to predict subcellular localization of proteins. Nucleic Acids Res. 2018; 46: W459-66.

[23] Sahu SS, Loaiza CD, Kaundal R . Plant-mSubP: a computational framework for the prediction of single-and multi-target protein subcellular localization using integrated machine-learning approaches. AoB Plants. 2020; 12: plz068.

[24] MacLean AM, Bravo A, Harrison MJ . Plant signaling and metabolic pathways enabling arbuscular mycorrhizal symbiosis. Plant Cell. 2017; 29: 2319-35.

[25] Hu XL, Lu H, Hassan MM et al. Advances and perspectives in discovery and functional analysis of small secreted proteins in plants. Hortic Res. 2021; 8: 130.

[26] Lanfranco L, Fiorilli V, Gutjahr C . Partner communication and role of nutrients in the arbuscular mycorrhizal symbiosis. New Phytol. 2018; 220: 1031-46.

[27] Yang G, Liu N, Lu W et al. The interaction between arbuscular mycorrhizal fungi and soil phosphorus availability influences plant community productivity and ecosystem stability. J Ecol. 2014; 102: 1072-82.

[28] Bona E, Cantamessa S, Massa N et al. Arbuscular mycorrhizal fungi and plant growth-promoting pseudomonads improve yield, quality and nutritional value of tomato: a field study. Mycorrhiza. 2017; 27: 1-11.

[29] Rillig MC, Aguilar-Trigueros CA, Camenzind T et al. Why farmers should manage the arbuscular mycorrhizal symbiosis. New Phytol. 2019; 222: 1171-5.

[30] Delaux PM . Comparative phylogenomics of symbiotic associations. New Phytol. 2017; 213: 89-94.

[31] Handa Y, Nishide H, Takeda N et al. RNA-seq transcriptional profiling of an arbuscular mycorrhiza provides insights into regulated and coordinated gene expression in Lotus japonicus and Rhizophagus irregularis . Plant Cell Physiol. 2015; 56: 1490-511.

[32] Mewalal R, Yin H, Hu R et al. Identification of Populus small RNAs responsive to mutualistic interactions with mycorrhizal fungi, Laccaria bicolor and Rhizophagus irregularis . Front Microbiol. 2019; 10: 515.

[33] Liu B, Hu J, Zhang J . Evolutionary divergence of duplicated Hsf genes in Populus . Cells. 2019; 8: 438.

[34] Jiang B, Shi Y, Peng Y et al. Cold-induced CBF-PIF3 interaction enhances freezing tolerance by stabilizing the phyB thermosensor in Arabidopsis. Mol Plant. 2020; 13: 894-906.

[35] Shikata M, Matsuda Y, Ando K et al. Characterization of Arabidopsis ZIM, a member of a novel plant-specific GATA factor gene family. J Exp Bot. 2004; 55: 631-9.

[36] Lei J, Jayaprakasha GK, Singh J et al. CIRCADIAN CLOCK-ASSOCIATED1 controls resistance to aphids by altering indole glucosinolate production. Plant Physiol. 2019; 181: 1344-59.

[37] Carter ME, Helm M, Chapman AVE et al. Convergent evolution of effector protease recognition by Arabidopsis and barley. Mol Plant-Microbe Interact. 2019; 32: 550-65.

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

[39] de Vries S, Stukenbrock EH, Rose LE . Rapid evolution in plant-microbe interactions - an evolutionary genomics perspective. New Phytol. 2020; 226: 1256-62.

[40] Labbé J, Muchero W, Czarnecki O et al. Mediation of plant-mycorrhizal interaction by a lectin receptor-like kinase. Nat Plants. 2019; 5: 676-80.

[41] Qiao Z, Yates TB, Shrestha HK et al. Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase. Plant Biotechnol J. 2021; 19: 2454-68.

[42] Liu D, Chen M, Mendoza B et al. CRISPR/Cas9-mediated targeted mutagenesis for functional genomics research of crassulacean acid metabolism plants. J Exp Bot. 2019; 70: 6621-9.

[43] Liu D, Hu R, Palla KJ et al. Advances and perspectives on the use of CRISPR/Cas9 systems in plant genomics research. Curr Opin Plant Biol. 2016; 30: 70-7.

[44] Yang X, Medford JI, Markel K et al. Plant biosystems design research roadmap 1.0. BioDes Res. 2020; 2020: 8051764.

[45] Hassan MM, Zhang Y, Yuan G et al. Construct design for CRISPR/Cas-based genome editing in plants. Trends Plant Sci. 2021; 26: 1133-52.

[46] Kinoshita T, Fujita M . Biosynthesis of GPI-anchored proteins: special emphasis on GPI lipid remodeling. J Lipid Res. 2016; 57: 6-24.

[47] Zhou K . Glycosylphosphatidylinositol-anchored proteins in Arabidopsis and one of their common roles in signaling transduction. Front Plant Sci. 2019; 10: 1022.

[48] Wang Q, Liu J, Zhu H . Genetic and molecular mechanisms underlying symbiotic specificity in legume-rhizobium interactions. Front Plant Sci. 2018; 9: 313.

[49] Yang X, Kalluri UC, DiFazio SP et al. Poplar genomics: state of the science. Crit Rev Plant Sci. 2009; 28: 285-308.

[50] Dharmawardhana P, Brunner A, Strauss S . Poplar as a tree model for horticulture and beyond: a case study of genome-scale changes in gene expression during bud entry and release from dormancy. Acta Hortic. 2009; 859: 43-7.

[51] Javot H, Pumplin N, Harrison MJ . Phosphate in the arbuscular mycorrhizal symbiosis: transport properties and regulatory roles. Plant Cell Environ. 2007; 30: 310-22.

[52] Pellegrin C, Daguerre Y, Ruytinx J et al. Laccaria bicolor MiSSP8 is a small-secreted protein decisive for the establishment of the ectomycorrhizal symbiosis . Environ Microbiol. 2019; 21: 3765-79.

[53] Paniagua C, Bilkova A, Jackson P et al. Dirigent proteins in plants: modulating cell wall metabolism during abiotic and biotic stress exposure. J Exp Bot. 2017; 68: 3287-301.

[54] Wang B, Qiu Y-L . Phylogenetic distribution and evolution of mycorrhizas in land plants. Mycorrhiza. 2006; 16: 299-363.

[55] Brundrett MC . Mycorrhizal associations and other means of nutrition of vascular plants: understanding the global diversity of host plants by resolving conflicting information and developing reliable means of diagnosis. Plant Soil. 2009; 320: 37-77.

[56] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 1-14.

[57] Katoh K, Rozewicki J, Yamada KD . MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019; 20: 1160-6.

[58] Minh BQ, Schmidt HA, Chernomor O et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020; 37: 1530-4.

[59] Kalyaanamoorthy S, Minh BQ, Wong TK et al. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017; 14: 587-9.

[60] Mirarab S, Reaz R, Bayzid M et al. ASTRAL: genome-scale coalescent-based species tree estimation. Bioinformatics. 2014; 30: i541-8.

[61] Käll L, Krogh A, Sonnhammer EL . Advantages of combined transmembrane topology and signal peptide prediction-the Phobius web server. Nucleic Acids Res. 2007; 35: W429-32.

[62] Armenteros JJA, Salvatore M, Emanuelsson O et al. Detecting sequence signals in targeting peptides using deep learning. Life Sci Alliance. 2019; 2: e201900429.

[63] Nugent T, Jones DT . Detecting pore-lining regions in transmembrane protein sequences. BMC Bioinformatics. 2012; 13: 1-9.

[64] Almagro Armenteros JJ, Sønderby CK, Sønderby SK et al. DeepLoc: prediction of protein subcellular localization using deep learning. Bioinformatics. 2017; 33: 3387-95.

[65] Uhlén M, Fagerberg L, Hallström BM et al. Tissue-based map of the human proteome. Science. 2015; 347: 1260419.

[66] Dobin A, Davis CA, Schlesinger F et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29: 15-21.

[67] Leng N, Dawson JA, Thomson JA et al. EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics. 2013; 29: 1035-43.

[68] Chow C-N, Lee TY, Hung YC et al. PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants. Nucleic Acids Res. 2019; 47: D1155-63.

[69] Kelley LA, Mezulis S, Yates CM et al. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015; 10: 845-58.

[70] 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.

[71] Thimm O, Bläsing O, Gibon Y et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004; 37: 914-39.