MicroRNAs expression dynamics reveal post-transcriptional mechanisms regulating seed development in Phaseolus vulgaris L.

Abstract

The knowledge on post-transcriptional regulation mechanisms implicated in seed development (SD) is still limited, particularly in one of the most consumed grain legumes, Phaseolus vulgaris L. We explore for the first time the miRNA expression dynamics in P. vulgaris developing seeds. Seventy-two known and 39 new miRNAs were found expressed in P. vulgaris developing seeds. Most of the miRNAs identified were more abundant at 10 and 40 days after anthesis, suggesting that late embryogenesis/early filling and desiccation were SD stages in which miRNA action is more pronounced. Degradome analysis and target prediction identified targets for 77 expressed miRNAs. While several known miRNAs were predicted to target HD-ZIP, ARF, SPL, and NF-Y transcription factors families, most of the predicted targets for new miRNAs encode for functional proteins. MiRNAs-targets expression profiles evidenced that these miRNAs could tune distinct seed developmental stages. MiRNAs more accumulated at early SD stages were implicated in regulating the end of embryogenesis, postponing the seed maturation program, storage compound synthesis and allocation. MiRNAs more accumulated at late SD stages could be implicated in seed quiescence, desiccation tolerance, and longevity with still uncovered roles in germination. The miRNAs herein described represent novel P. vulgaris resources with potential application in future biotechnological approaches to modulate the expression of genes implicated in legume seed traits with impact in horticultural production systems.

References

[1] Polania, J. A., Chater, C. C. C., Covarrubias, A. A. & Rao, I. M. Phaseolus species responses and tolerance to drought. in The Plant Family Fabaceae (eds Hasanuzzaman, M., Araújo, S. & Gill, S. S.) 319-336 (Springer Singapore, 2020).

[2] Vaz Patto, M. C. & Araújo, S. S. Positioning Portugal into the context of world production and research in grain legumes. Rev. Ciências Agrárias 39, 471-489 (2016).

[3] Kang, Y., Li, M., Sinharoy, S. & Verdier, J. A snapshot of functional genetic studies in Medicago truncatula. Front. Plant Sci. 7, 1175 (2016).

[4] Parreira, J. R. et al. Differential proteomics reveals the hallmarks of seed development in common bean (Phaseolus vulgaris L.). J. Proteom. 143, 188-198 (2016).

[5] Parreira, J. R., Balestrazzi, A., Fevereiro, P. & Araújo, S. Maintaining genome integrity during seed development in Phaseolus vulgaris L.: evidence from a transcriptomic profiling study. Genes (Basel) 9, 463 (2018).

[6] Yu, L. et al. Genome-wide identification and characterization of novel microRNAs in seed development of soybean. Biosci. Biotechnol. Biochem. 83, 233-242 (2019).

[7] DeBoer, K. et al. Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development. BMC Genomics 20, 135 (2019).

[8] Yu, Y., Jia, T. & Chen, X. The ‘how’ and ‘where’ of plant microRNAs. N. Phytol. 216, 1002-1017 (2017).

[9] Arenas-Huertero, C. et al. Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. Plant Mol. Biol. 70, 385-401 (2009).

[10] Valdés-López, O. et al. MicroRNA expression profile in common bean (Phaseolus vulgaris) under nutrient deficiency stresses and manganese toxicity. N. Phytol. 187, 805-818 (2010).

[11] Formey, D. et al. Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing. BMC Genomics 16, 423 (2015).

[12] Pelaez, P. et al. Identification and characterization of microRNAs in Phaseolus vulgaris by high-throughput sequencing. BMC Genomics 13, 83 (2012).

[13] Contreras-Cubas, C. et al. The Phaseolus vulgaris miR159a precursor encodes a second differentially expressed microRNA. Plant Mol. Biol. 80, 103-115 (2012).

[14] Patwa, N., Nithin, C., Bahadur, R. P. & Basak, J. Identification and characterization of differentially expressed Phaseolus vulgaris miRNAs and their targets during mungbean yellow mosaic India virus infection reveals new insight into Phaseolus-MYMIV interaction. Genomics 111, 1333-1342 (2019).

[15] Han, J. et al. Computational identification of miRNAs and their targets in Phaseolus vulgaris. Genet. Mol. Res. 13, 310-322 (2014).

[16] Valdés-López, O. et al. Essential role of MYB transcription factor: PvPHR1 and microRNA: PvmiR399 in phosphorus-deficiency signalling in common bean roots. Plant. Cell Environ. 31, 1834-1843 (2008).

[17] Zhang, B. H., Pan, X. P., Cox, S. B., Cobb, G. P. & Anderson, T. A. Evidence that miRNAs are different from other RNAs. Cell. Mol. Life Sci. 63, 246-254 (2006).

[18] Verdier, J. & Thompson, R. D. Transcriptional regulation of storage protein synthesis during dicotyledon seed filling. Plant Cell Physiol. 49, 1263-1271 (2008).

[19] Smith, Z. R. & Long, J. A. Control of Arabidopsis apical-basal embryo polarity by antagonistic transcription factors. Nature 464, 423-426 (2010).

[20] Manan, S. et al. Soybean LEC2 regulates subsets of genes involved in controlling the biosynthesis and catabolism of seed storage substances and seed development. Front. Plant Sci. 8, 1-16 (2017).

[21] Wang, X. G. et al. Deciphering transcriptome and complex alternative splicing transcripts in mammary gland tissues from cows naturally infected with Staphylococcus aureus mastitis. PLoS ONE 11, e0159719 (2016).

[22] Frontini, M., Imbriano, C., Manni, I. & Mantovani, R. Cell-cycle regulation of NF-YC nuclear localization. Cell Cycle 3, 205-210 (2004).

[23] Jo, L., Pelletier, J. M. & Harada, J. J. Central role of the LEAFY COTYLEDON1 transcription factor in seed development. J. Integr. Plant Biol. 61, 564-580 (2019).

[24] Mu, J., Tan, H., Hong, S., Liang, Y. & Zuo, J. Arabidopsis transcription factor genes NF-YA1, 5, 6, and 9 play redundant roles in male gametogenesis, embryogenesis, and seed development. Mol. Plant 6, 188-201 (2013).

[25] Nodine, M. D. & Bartel, D. P. MicroRNAs prevent precocious gene expression and enable pattern formation during plant embryogenesis. Genes Dev. 24, 2678-2692 (2010).

[26] Wang, F. & Perry, S. E. Identification of direct targets of FUSCA3, a key regulator of Arabidopsis seed development. Plant Physiol. 161, 1251-1264 (2013).

[27] Xu, M. et al. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLOS Genet. 12, e1006263 (2016).

[28] Zheng, Z. et al. miR159 represses a constitutive pathogen defense response in tobacco. Plant Physiol. 182, 2182-2198 (2020).

[29] Tsuzuki, M. et al. Profiling and characterization of small RNAs in the liverwort, Marchantia polymorpha, belonging to the first diverged land plants. Plant Cell Physiol. 57, 359-372 (2016).

[30] Zhou, Y., Qu, H., Dibley, K. E., Offler, C. E. & Patrick, J. W. A suite of sucrose transporters expressed in coats of developing legume seeds includes novel pH-independent facilitators. Plant J. 49, 750-764 (2007).

[31] Coelho, C. M. M., Benedito, V. A., Figueira, A., Vitorello, V. A. & Azevedo, R. A. Variation in the enzyme activity and gene expression of myo-inositol-3-phosphate synthase and phytate accumulation during seed development in common bean (Phaseolus vulgaris L.). Acta Physiol. Plant. 29, 265-271 (2007).

[32] Lin, S.-I. et al. Regulatory network of microRNA399 and PHO2 by systemic signaling. Plant Physiol. 147, 732-746 (2008).

[33] Liang, G., Yang, F. & Yu, D. MicroRNA395 mediates regulation of sulfate accumulation and allocation in Arabidopsis thaliana. Plant J. 62, 1046-1057 (2010).

[34] Awazuhara, M. et al. The function of SULTR2;1 sulfate transporter during seed development in Arabidopsis thaliana. Physiol. Plant. 125, 95-105 (2005).

[35] Marsolais, F. et al. Proteomic analysis of common bean seed with storage protein deficiency reveals up-regulation of sulfur-rich proteins and starch and raffinose metabolic enzymes, and down-regulation of the secretory pathway. J. Proteom. 73, 1587-1600 (2010).

[36] Tian, R. et al. Direct and indirect targets of the Arabidopsis seed transcription factor ABSCISIC ACID INSENSITIVE3. Plant J. 103, 1679-1694 (2020).

[37] Changan, S. S., Ali, K., Kumar, V., Garg, N. K. & Tyagi, A. Abscisic acid biosynthesis under water stress: anomalous behavior of the 9-cis-epoxycarotenoid dioxygenase1 (NCED1) gene in rice. Biol. Plant. 62, 663-670 (2018).

[38] Jo, L. et al. Combinatorial interactions of the LEC1 transcription factor specify diverse developmental programs during soybean seed development. Proc. Natl Acad. Sci. USA 117, 1223-1232 (2020).

[39] Battaglia, M., Olvera-Carrillo, Y., Garciarrubio, A., Campos, F. & Covarrubias, A. A. The enigmatic LEA proteins and other hydrophilins. Plant Physiol. 148, 6-24 (2008).

[40] Hundertmark, M., Buitink, J., Leprince, O . & Hincha, D. K. The reduction of seed-specific dehydrins reduces seed longevity in Arabidopsis thaliana. Seed Sci. Res. 21, 165-173 (2011).

[41] Martin, R. C., Liu, P.-P., Goloviznina, N. A. & Nonogaki, H. microRNA, seeds, and Darwin?: diverse function of miRNA in seed biology and plant responses to stress. J. Exp. Bot. 61, 2229-2234 (2010).

[42] Xu, M. Y. et al. Stress-induced early flowering is mediated by miR169 in Arabidopsis thaliana. J. Exp. Bot. 65, 89-101 (2014).

[43] Stephenson, T. J., McIntyre, C. L., Collet, C. & Xue, G.-P. Genome-wide identification and expression analysis of the NF-Y family of transcription factors in Triticum aestivum. Plant Mol. Biol. 65, 77-92 (2007).

[44] Yin, H. et al. miR156/SPL9 regulates reactive oxygen species accumulation and immune response in Arabidopsis thaliana. Phytopathology 109, 632-642 (2019).

[45] Rocak, S. & Linder, P. DEAD-box proteins: the driving forces behind RNA metabolism. Nat. Rev. Mol. Cell Biol. 5, 232-241 (2004).

[46] Lorković, Z. J., Herrmann, R. G. & Oelmüller, R. PRH75, a new nucleus-localized member of the DEAD-box protein family from higher plants. Mol. Cell. Biol. 17, 2257-2265 (1997).

[47] Nayak, N. R. et al. An Arabidopsis ATP-dependent, DEAD-box RNA helicase loses activity upon IsoAsp formation but is restored by PROTEIN ISOASPARTYL METHYLTRANSFERASE. Plant Cell 25, 2573-2586 (2013).

[48] Bai, B. et al. Seed-stored mRNAs that are specifically associated to monosomes are translationally regulated during germination. Plant Physiol. 182, 378-392 (2020).

[49] Wei, Z., Liu, X., Feng, T. & Chang, Y. Novel and conserved microRNAs in dalian purple urchin (Strongylocentrotus nudus) identified by next generation sequencing. Int. J. Biol. Sci. 7, 180-192 (2011).

[50] Li, M. et al. MicroRNAome of porcine pre- and postnatal development. PLoS ONE 5, e11541 (2010).

[51] Griffiths-Jones, S. The microRNA Registry. Nucleic Acids Res. 32, D109-D111 (2004).

[52] Meyers, B. C. et al. Criteria for annotation of plant microRNAs. Plant Cell 20, 3186-3190 (2008).

[53] Zhong, M. et al. The effect of cadmium on the microRNAome, degradome and transcriptome of rice seedlings. Plant Growth Regul. 90, 15-27 (2020).

[54] Zhu, H. et al. MicroRNA528, a hub regulator modulating ROS homeostasis via targeting of a diverse set of genes encoding copper‐containing proteins in monocots. N. Phytol. 225, 385-399 (2020).

[55] Martin, M . Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10 (2011).

[56] Addo-Quaye, C., Miller, W. & Axtell, M. J. CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics 25, 130-131 (2009).

[57] Lohse, M. et al. Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant. Cell Environ. 37, 1250-1258 (2014).

[58] Shannon, P. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498-2504 (2003).

[59] Androvic, P., Valihrach, L., Elling, J., Sjoback, R. & Kubista, M. Two-tailed RT-qPCR: a novel method for highly accurate miRNA quantification. Nucleic Acids Res. 45, e144- e144 (2017).

[60] Ruijter, J. M. et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 37, e45- e45 (2009).

[61] Pfaffl, M. W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29, 45e- 45e (2001).