Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening

Abstract

The perception and signal transduction of the plant hormone abscisic acid (ABA) are crucial for strawberry fruit ripening, but the underlying mechanism of how ABA regulates ripening-related genes has not been well understood. By employing high-throughput sequencing technology, we comprehensively analyzed transcriptomic and miRNA expression profiles simultaneously in ABA- and nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker)-treated strawberry fruits with temporal resolution. The results revealed that ABA regulated many genes in different pathways, including hormone signal transduction and the biosynthesis of secondary metabolites. Transcription factor genes belonging to WRKY and heat shock factor (HSF) families might play key roles in regulating the expression of ABA inducible genes, whereas the KNOTTED1-like homeobox protein and Squamosa Promoter-Binding-like protein 18 might be responsible for ABA-downregulated genes. Additionally, 20 known and six novel differentially expressed miRNAs might be important regulators that assist ABA in regulating target genes that are involved in versatile physiological processes, such as hormone balance regulation, pigments formation and cell wall degradation. Furthermore, degradome analysis showed that one novel miRNA, Fa_novel6, could degrade its target gene HERCULES1, which likely contributed to fruit size determination during strawberry ripening. These results expanded our understanding of how ABA drives the strawberry fruit ripening process as well as the role of miRNAs in this process.

References

[1] Chai, Y. M., Jia, H. F., Li, C. L., Dong, Q. H. & Shen, Y. Y. FaPYR1 is involved in strawberry fruit ripening. J. Exp. Bot. 62, 5079-5089 (2011).

[2] Jia, H. F. et al. Abscisic acid plays an important role in the regulation of strawberry fruit ripening. Plant Physiol. 157, 188-199 (2011).

[3] Hubbard, K. E., Nishimura, N., Hitomi, K., Getzoff, E. D. & Schroeder, J. I. Early abscisic acid signal transduction mechanisms: newly discovered components and newly emerging questions. Gene Dev. 24, 1695-1708 (2010).

[4] Cutler, S. R., Rodriguez, P. L., Finkelstein, R. R. & Abrams, S. R. Abscisic acid: emergence of a core signaling network. Annu. Rev. Plant Biol. 61, 651-679 (2010).

[5] Jia, H. F. et al. Type 2C protein phosphatase ABI1 is a negative regulator of strawberry fruit ripening. J. Exp. Bot. 64, 1677-1687 (2013).

[6] Medina-Puche, L. et al. MYB10 plays a major role in the regulation of flavonoid/phenylpropanoid metabolism during ripening of Fragaria ananassa fruits. J. Exp. Bot. 65, 401-417 (2014).

[7] Vallarino, J. G. et al. Central role of FaGAMYB in the transition of the strawberry receptacle from development to ripening. New Phytol. 208, 482-496 (2015).

[8] Medina-Puche, L. et al. An R2R3-MYB transcription factor regulates eugenol production in ripe strawberry fruit receptacles. Plant Physiol. 168, 598-614 (2015).

[9] Jones-Rhoades, M. W., Bartel, D. P. & Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 57, 19-53 (2006).

[10] Chen, K. & Rajewsky, N. The evolution of gene regulation by transcription factors and microRNAs. Nat. Rev. Genet. 8, 93-103 (2007).

[11] Xia, R., Ye, S. Q., Liu, Z. R., Meyers, B. C. & Liu, Z. C. Novel and recently evolved microRNA clusters regulate expansive F-BOX gene networks through phased small interfering RNAs in wild diploid strawberry. Plant Physiol. 169, 594-610 (2015).

[12] Li, D. D. et al. Developmental and stress regulation on expression of a novel miRNA, Fan-miR73, and its target ABI5 in strawberry. Sci. Rep.-Uk 6, 28358 (2016).

[13] Moxon, S. et al. Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res. 18, 1602-1609 (2008).

[14] Karlova, R. et al. Identification of microRNA targets in tomato fruit development using high-throughput sequencing and degradome analysis. J. Exp. Bot. 64, 1863-1878 (2013).

[15] Symons, G. M. et al. Hormonal changes during non-climacteric ripening in strawberry. J. Exp. Bot. 63, 4741-4750 (2012).

[16] Hu, Y. et al. Identification, isolation, and expression analysis of heat shock transcription factors in the diploid woodland strawberry Fragaria vesca. Front Plant Sci. 6, 00736 (2015).

[17] Ueguchi-Tanaka, M., Nakajima, M., Motoyuki, A. & Matsuoka, M. Gibberellin receptor and its role in gibberellin signaling in plants. Annu. Rev. Plant. Biol 58, 183-198 (2007).

[18] Perez, A. C. et al. The non-JAZ TIFY protein TIFY8 from Arabidopsis thaliana is a transcriptional repressor. PLoS ONE 9, e84891 (2014).

[19] Kong, Y. Z. et al. Molecular analysis of a family of Arabidopsis genes related to Galacturonosyltransferases. Plant Physiol. 155, 1791-1805 (2011).

[20] Ge, A. J. et al. Deep sequencing discovery of novel and conserved microRNAs in strawberry (Fragariaxananassa). Physiol. Plant. 148, 387-396 (2013).

[21] Wu, J. X., Zheng, S. S., Feng, G. Z. & Yi, H. L. Comparative analysis of miRNAs and their target transcripts between a spontaneous late-ripening sweet orange mutant and its wild-type using small RNA and degradome sequencing. Front. Plant Sci. 7, 1416 (2016).

[22] Guo, H. Q. et al. Three related receptor-like kinases are required for optimal cell elongation in Arabidopsis thaliana. Proc. Natl Acad. Sci. USA 106, 7648-7653 (2009).

[23] Loake, G. & Grant, M. Salicylic acid in plant defence-the players and protagonists. Curr. Opin. Plant Biol. 10, 466-472 (2007).

[24] Jiang, C. J. et al. Abscisic acid interacts antagonistically with salicylic acid signaling pathway in rice-magnaporthe grisea interaction. Mol. Plant Microbe. 23, 791-798 (2010).

[25] Larkindale, J. & Knight, M. R. Protection against heat stress-induced oxidative damage in arabidopsis involves calcium, abscisic acid, ethylene, and salicylic acid. Plant Physiol. 128, 682-695 (2002).

[26] Seo, P. J. & Park, C. M. MYB96-mediated abscisic acid signals induce pathogen resistance response by promoting salicylic acid biosynthesis in Arabidopsis. New Phytol. 186, 471-483 (2010).

[27] Ji, K. et al. Non-climacteric ripening in strawberry fruit is linked to ABA, FaNCED2 and FaCYP707A1. Funct. Plant Biol. 39, 351-357 (2012).

[28] Merchante, C. et al. Ethylene is involved in strawberry fruit ripening in an organ-specific manner. J. Exp. Bot. 64, 4421-4439 (2013).

[29] Kang, C. Y. et al. Genome-scale transcriptomic insights into early-stage fruit development in woodland strawberry Fragaria vesca. Plant Cell 25, 1960-1978 (2013).

[30] Sanchez-Sevilla, J. F. et al. Gene expression atlas of fruit ripening and transcriptome assembly from RNA-seq data in octoploid strawberry (Fragaria x ananassa). Sci. Rep.-Uk 7, 13737 (2017).

[31] McAtee, P., Karim, S., Schaffer, R. & David, K. A dynamic interplay between phytohormones is required for fruit development, maturation, and ripening. Front. Plant Sci. 4, 79 (2013).

[32] Vernoux, T. et al. The auxin signalling network translates dynamic input into robust patterning at the shoot apex. Mol. Syst. Biol. 7, 508 (2011).

[33] Piya, S., Shrestha, S. K., Binder, B., Stewart, C. N. & Hewezi, T. Protein−protein interaction and gene co-expression maps of ARFs and Aux/IAAs in Arabidopsis. Front. Plant Sci. 5, 744 (2014).

[34] Li, P., Chen, X., Sun, F. & Dong, H. Tobacco TTG2 and ARF8 function concomitantly to control flower colouring by regulating anthocyanin synthesis genes. Plant Biol. 19, 525-532 (2017).

[35] Trainotti, L., Pavanello, A. & Casadoro, G. Different ethylene receptors show an increased expression during the ripening of strawberries: does such an increment imply a role for ethylene in the ripening of these non-climacteric fruits?. J. Exp. Bot. 56, 2037-2046 (2005).

[36] Iannetta, P. P. M. et al. Ethylene and carbon dioxide production by developing strawberries show a correlative pattern that is indicative of ripening climacteric fruit. Physiol. Plant. 127, 247-259 (2006).

[37] Aharoni, A. & O’Connell, A. P. Gene expression analysis of strawberry achene and receptacle maturation using DNA microarrays. J. Exp. Bot. 53, 2073-2087 (2002).

[38] Albert, N. W. et al. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell 26, 962-980 (2014).

[39] Kui, L. W. et al. An R2R3 MYB transcription factor associated with regulation of the anthocyanin biosynthetic pathway in Rosaceae. BMC Plant Biol. 10, 50 (2010).

[40] Hawkins, C., Caruana, J., Schiksnis, E. & Liu, Z. C. Genome-scale DNA variant analysis and functional validation of a SNP underlying yellow fruit color in wild strawberry. Sci. Rep.-UK 6, 29017 (2016).

[41] Aharoni, A. et al. The strawberry FaMYB1 transcription factor suppresses anthocyanin and flavonol accumulation in transgenic tobacco. Plant J. 28, 319-332 (2001).

[42] Kadomura-Ishikawa, Y., Miyawaki, K., Takahashi, A. & Noji, S. RNAi-mediated silencing and overexpression of the FaMYB1 gene and its effect on anthocyanin accumulation in strawberry fruit. Biol. Plant. 59, 677-685 (2015).

[43] Zhang, Y. C. et al. Transcript quantification by RNA-Seq reveals differentially expressed genes in the red and yellow fruits of Fragaria vesca. PLoS ONE 10, e0144356 (2015).

[44] Schaart, J. G. et al. Identification and characterization of MYB-bHLH-WD40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry (Fragaria x ananassa) fruits. New Phytol. 197, 454-467 (2013).

[45] Zhou, H. Y. et al. Genome-wide analysis of the expression of WRKY family genes in different developmental stages of wild strawberry (Fragaria vesca) fruit. PLoS ONE 11, e0154312 (2016).

[46] Wei, W. et al. The WRKY transcription factors in the diploid woodland strawberry Fragaria vesca: identification and expression analysis under biotic and abiotic stresses. Plant Physiol. Bioch. 105, 129-144 (2016).

[47] Encinas-Villarejo, S. et al. Evidence for a positive regulatory role of strawberry (Fragariaxananassa) Fa WRKY1 and Arabidopsis At WRKY75 proteins in resistance. J. Exp. Bot. 60, 3043-3065 (2009).

[48] Jia, X. Y. et al. Differential and dynamic regulation of miR398 in response to ABA and salt stress in Populus tremula and Arabidopsis thaliana. Plant Mol. Biol. 71, 51-59 (2009).

[49] Dugas, D. V. & Bartel, B. Sucrose induction of Arabidopsis miR398 represses two Cu/Zn superoxide dismutases. Plant Mol. Biol. 67, 403-417 (2008).

[50] Han, Y. et al. SUCROSE NONFERMENTING1-RELATED PROTEIN KINASE2.6, an ortholog of OPEN STOMATA1, is a negative regulator of strawberry fruit development and ripening. Plant Physiol. 167, 915-930 (2015).

[51] Sunkar, R., Kapoor, A. & Zhu, J. K. Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell 18, 2415- 2415 (2006).

[52] Daminato, M., Guzzo, F. & Casadoro, G. A SHATTERPROOF-like gene controls ripening in non-climacteric strawberries, and auxin and abscisic acid antagonistically affect its expression. J. Exp. Bot. 64, 3775-3786 (2013).

[53] Baima, S. et al. Negative feedback regulation of Auxin signaling by ATHB8/ACL5-BUD2 transcription module. Mol. Plant 7, 1006-1025 (2014).

[54] Cui, L. G., Shan, J. X., Shi, M., Gao, J. P. & Lin, H. X. The miR156-SPL9-DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants (vol 80, pg 1108, 2014). Plant J. 82, 901- 901 (2015).

[55] Feder, A. et al. A Kelch domain-containing F-Box coding gene negatively regulates flavonoid accumulation in muskmelon. Plant Physiol. 169, 1714-1726 (2015).

[56] Steinwand, B. J. & Kieber, J. J. The role of receptor-like kinases in regulating cell wall function. Plant Physiol. 153, 479-484 (2010).

[57] Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644-U130 (2011).

[58] Shulaev, V. et al. The genome of woodland strawberry (Fragaria vesca). Nat. Genet. 43, 109-116 (2011).

[59] Hirakawa, H. et al. Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species. DNA Res. 21, 169-181 (2014).

[60] Mortazavi, A., Williams, B. A., Mccue, K., Schaeffer, L. & Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621-628 (2008).

[61] German, M. A. et al. Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat. Biotechnol. 26, 941-946 (2008).

[62] Addo-Quaye, C. et al. Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome. RNA 15, 2112-2121 (2009).

[63] Yu, L. J. et al. Comparative transcriptome analysis of transporters, phytohormone and lipid metabolism pathways in response to arsenic stress in rice (Oryza sativa). New Phytol. 195, 97-112 (2012).