Genome-wide identification of drought-responsive microRNAs in two sets of Malus from interspecific hybrid progenies

Abstract

Drought stress can negatively impact apple fruit quality and yield. Apple microRNAs (miRNAs) participate in apple tree and fruit development, as well as in biotic stress tolerance; however, it is largely unknown whether these molecules are involved in the drought response. To identify drought-responsive miRNAs in Malus, we first examined the drought stress tolerance of ten F₁ progenies of R3 (M. × domestica) × M. sieversii. We performed Illumina sequencing on pooled total RNA from both drought-tolerant and drought-sensitive plants. The sequencing results identified a total of 206 known miRNAs and 253 candidate novel miRNAs from drought-tolerant plants and drought-sensitive plants under control or drought conditions. We identified 67 miRNAs that were differentially expressed in drought-tolerant plants compared with drought-sensitive plants under drought conditions. Under drought stress, 61 and 35 miRNAs were differentially expressed in drought-tolerant and drought-sensitive plants, respectively. We determined the expression levels of seven out of eight miRNAs by stem-loop qPCR analysis. We also predicted the target genes of all differentially expressed miRNAs and identified the expression of some genes. Gene Ontology analyses indicated that the target genes were mainly involved in stimulus response and cellular and metabolic processes. Finally, we confirmed roles of two miRNAs in apple response to mannitol. Our results reveal candidate miRNAs and their associated mRNAs that could be targeted for improving drought tolerance in Malus species, thus providing a foundation for understanding the molecular networks involved in the response of apple trees to drought stress.

References

[1] Gu, S. & Kay, M. A. How do miRNAs mediate translational repression?. Silence 1, 11 (2010).

[2] Zamore, P. D. & Haley, B. Ribo-gnome: the big world of small RNAs. Science 309, 1519-1524 (2005).

[3] Bonnet, E., Van de Peer, Y. & Rouzé, P. The small RNA world of plants. New Phytol. 171, 451-468 (2006).

[4] Zhu, J. K. Reconstituting plant miRNA biogenesis. Proc. Natl Acad. Sci. USA 105, 9851-9852 (2008).

[5] Carrington, J. C. & Ambros, V. Role of microRNAs in plant and animal development. Science 301, 336 (2003).

[6] Sunkar, R., Li, Y. F. & Jagadeeswaran, G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 17, 196-203 (2012).

[7] Sunkar, R. & Zhu, J.-K. Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 16, 2001-2019 (2004).

[8] Liu, H. H., Tian, X., Li, Y. J., Wu, C. A. & Zheng, C. C. Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA 14, 836-843 (2008).

[9] Li, W.-X. et al. The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and posttranscriptionally to promote drought resistance. Plant Cell 20, 2238-2251 (2008).

[10] Zhou, L. et al. Genome-wide identification and analysis of drought-responsive microRNAs in Oryza sativa. J. Exp. Bot. 61, 4157-4168 (2010).

[11] Zhao, B. et al. Identification of drought-induced microRNAs in rice. Biochem. Biophys. Res. Commun. 354, 585-590 (2007).

[12] Cheah, B. H., Nadarajah, K., Divate, M. D. & Wickneswari, R. Identification of four functionally important microRNA families with contrasting differential expression profiles between drought-tolerant and susceptible rice leaf at vegetative stage. BMC Genom. 16, 692 (2015).

[13] Shuai, P., Liang, D., Zhang, Z., Yin, W. & Xia, X. Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis. BMC Genom. 14, 233 (2013).

[14] Lu, S., Sun, Y. H. & Chiang, V. L. Stress‐responsive microRNAs in Populus. Plant J. 55, 131-151 (2008).

[15] Candar-Cakir, B., Arican, E. & Zhang, B. Small RNA and degradome deep sequencing reveals drought-and tissue-specific microRNAs and their important roles in drought-sensitive and drought-tolerant tomato genotypes. Plant Biotechnol. J. 14, 1727-1746 (2016).

[16] Xie, F., Wang, Q., Sun, R. & Zhang, B. Deep sequencing reveals important roles of microRNAs in response to drought and salinity stress in cotton. J. Exp. Bot. 66, 789-804 (2015).

[17] Ma, C., Burd, S. & Lers, A. miR408 is involved in abiotic stress responses in Arabidopsis. Plant J. 84, 169-187 (2015).

[18] Trindade, I., Capitão, C., Dalmay, T., Fevereiro, M. P. & Santos, D. Md miR398 and miR408 are up-regulated in response to water deficit in Medicago truncatula. Planta 231, 705-716 (2010).

[19] Bohnert, H. J., Nelson, D. E. & Jensen, R. G. Adaptations to environmental stresses. Plant Cell 7, 1099-1111 (1995).

[20] Han, Y., Zhang, X., Wang, Y. & Ming, F. The suppression of WRKY44 by GIGANTEA-miR172 pathway is involved in drought response of Arabidopsis thaliana. PLoS ONE 8, e73541 (2013).

[21] Chen, X. et al. SQUAMOSA promoter‐binding protein‐like transcription factors: star players for plant growth and development. J. Integr. Plant Biol. 52, 946-951 (2010).

[22] Wang, J.-W., Czech, B. & Weigel, D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana. Cell 138, 738-749 (2009).

[23] Ding, Y., Tao, Y. & Zhu, C. Emerging roles of microRNAs in the mediation of drought stress response in plants. J. Exp. Bot. 64, 3077-3086 (2013).

[24] Li, C. et al. Silencing the SpMPK1, SpMPK2, and SpMPK3 genes in tomato reduces abscisic acid-mediated drought tolerance. Int. J. Mol. Sci. 14, 21983-21996 (2013).

[25] Naor, A., Naschitz, S., Peres, M. & Gal, Y. Responses of apple fruit size to tree water status and crop load. Tree Physiol. 28, 1255-1261 (2008).

[26] Xing, L. et al. Shoot bending promotes flower bud formation by miRNA-mediated regulation in apple (Malus domestica Borkh.). Plant Biotechnol. J. 14, 749-770 (2016).

[27] Guo, X. et al. Small RNA-sequencing links physiological changes and RdDM process to vegetative-to-floral transition in apple. Front. Plant Sci. 8, 873 (2017).

[28] Jia‐Long, Y. et al. A microRNA allele that emerged prior to apple domestication may underlie fruit size evolution. Plant J. 84, 417-427 (2015).

[29] Ma, C. et al. Cloning and characterization of miRNAs and their targets, including a novel miRNA-targeted NBS-LRR protein class gene in apple (Golden Delicious). Mol. Plant 7, 218-230 (2014).

[30] Kaja, E. et al. Identification of apple miRNAs and their potential role in fire blight resistance. Tree Genet. Genomes 11, 812 (2014).

[31] Feng, H. et al. microRNAs and their targets in apple (Malus domestica cv. “Fuji”) involved in response to infection of pathogen Valsa mali. Front. Plant Sci. 8, 2081 (2017).

[32] Liu, B. et al. Influence of rootstock on antioxidant system in leaves and roots of young apple trees in response to drought stress. Plant Growth Regul. 67, 247-256 (2012).

[33] Abe, K., Kotoda, N., Kato, H. & Soejima, J.-i Genetic studies on resistance to Valsa canker in apple: genetic variance and breeding values estimated from intra- and inter-specific hybrid progeny populations. Tree Genet. Genomes 7, 363-372 (2011).

[34] Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[35] Li, X. et al. Characterization and comparative profiling of the small RNA transcriptomes in two phases of flowering in Cymbidium ensifolium. BMC Genom. 16, 622 (2015).

[36] Friedländer, M. R., Mackowiak, S. D., Li, N., Chen, W. & Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37-52 (2012).

[37] Moxon, S. et al. A toolkit for analysing large-scale plant small RNA datasets. Bioinformatics 24, 2252-3 (2008).

[38] Wen, M., Shen, Y., Shi, S. & Tang, T. miREvo: an integrative microRNA evolutionary analysis platform for next-generation sequencing experiments. BMC Bioinformatics 13, 140 (2012).

[39] Zhou, L. et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PLoS ONE 5, e15224 (2011).

[40] Dai, X. & Zhao, P. X. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 39, W155-W159 (2011).

[41] Tian, T. et al. agriGOv2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 45, W122-W129 https://doi.org/10.1093/nar/gkx382 (2017).

[42] Bustin, S. A. Why the need for qPCR publication guidelines?-The case for MIQE. Methods 50, 217-226 (2010).

[43] Varkonyi-Gasic, E., Wu, R., Wood, M., Walton, E. F. & Hellens, R. P. Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 3, 12 (2007).

[44] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402-408 (2001).

[45] Wang, H. et al. Comprehensive genomic analysis of the TYROSINE AMINOTRANSFERASE (TAT) genes in apple (Malus domestica) allows the identification of MdTAT2 conferring tolerance to drought and osmotic stresses in plants. Plant Physiol. Biochem. 133, 81-91 (2018).

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

[47] Guan, Q., Lu, X., Zeng, H., Zhang, Y. & Zhu, J. Heat stress induction of miR398 triggers a regulatory loop that is critical for thermotolerance in Arabidopsis. Plant J. 74, 840-851 (2013).

[48] Jones-Rhoades, M. W. & Bartel, D. P. Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol. Cell 14, 787-799 (2004).

[49] Li, B., Qin, Y., Duan, H., Yin, W. & Xia, X. Genome-wide characterization of new and drought stress responsive microRNAs in Populus euphratica. J. Exp. Bot. 62, 3765-3779 (2011).

[50] Yao, Y. et al. Non-coding small RNAs responsive to abiotic stress in wheat (Triticum aestivum L.). Funct. Integr. Genom. 10, 187-190 (2010).

[51] Zhang, L. et al. A genome-wide characterization of microRNA genes in maize. PLoS Genet. 5, e1000716 (2009).

[52] Kulcheski, F. R. et al. Identification of novel soybean microRNAs involved in abiotic and biotic stresses. BMC Genom. 12, 307 (2011).

[53] Eldem, V. et al. Genome-wide identification of miRNAs responsive to drought in peach (Prunus persica) by high-throughput deep sequencing. PLoS ONE 7, e50298 (2012).

[54] Hackenberg, M., Gustafson, P., Langridge, P. & Shi, B. J. Differential expression of microRNAs and other small RNAs in barley between water and drought conditions. Plant Biotechnol. J. 13, 2-13 (2015).

[55] Yates, S. A. et al. De novo assembly of red clover transcriptome based on RNA-Seq data provides insight into drought response, gene discovery and marker identification. BMC Genom. 15, 453 (2014).

[56] Bai, Y., Dougherty, L., Cheng, L., Zhong, G.-Y. & Xu, K. Uncovering co-expression gene network modules regulating fruit acidity in diverse apples. BMC Genom. 16, 612 (2015).

[57] Wang, N. et al. Comparative transcriptomes analysis of red- and white-fleshed apples in an F1 population of Malus sieversii f. niedzwetzkyana crossed with M. domestica ‘Fuji’. PLoS ONE 10, e0133468 (2015).

[58] Xia, R., Zhu, H., An, Y.-q, Beers, E. P. & Liu, Z. Apple miRNAs and tasiRNAs with novel regulatory networks. Genome Biol. 13, R47 (2012).

[59] Wang, T., Chen, L., Zhao, M., Tian, Q. & Zhang, W.-H. Identification of drought-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. BMC Genom. 12, 367 (2011).

[60] Fang, Y., Xie, K. & Xiong, L. Conserved miR164-targeted NAC genes negatively regulate drought resistance in rice. J. Exp. Bot. 65, 2119-2135 (2014).

[61] Jung, Y. et al. BrRZFP1 a Brassica rapa C3HC4‐type RING zinc finger protein involved in cold, salt and dehydration stress. Plant Biol. 15, 274-283 (2013).

[62] Abogadallah, G. M., Nada, R. M., Malinowski, R. & Quick, P. Overexpression of HARDY, an AP2/ERF gene from Arabidopsis, improves drought and salt tolerance by reducing transpiration and sodium uptake in transgenic Trifolium alexandrinum L.. Planta 233, 1265-1276 (2011).

[63] Licausi, F., Ohme‐Takagi, M. & Perata, P. APETALA2/ethylene responsive factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 199, 639-649 (2013).

[64] Quan, R. et al. Overexpression of an ERF transcription factor TSRF1 improves rice drought tolerance. Plant Biotechnol. J. 8, 476-488 (2010).

[65] Bakhshi, B. et al. The contrasting microRNA content of a drought tolerant and a drought susceptible wheat cultivar. J. Plant Physiol. 216, 35-43 (2017).

[66] Hameed, A., Goher, M. & Iqbal, N. Drought induced programmed cell death and associated changes in antioxidants, proteases, and lipid peroxidation in wheat leaves. Biol. Plant 57, 370-374 (2013).

[67] Borrell, A. K. et al. Drought adaptation of stay-green sorghum is associated with canopy development, leaf anatomy, root growth, and water uptake. J. Exp. Bot. 65, 6251-6263 (2014).

[68] Ding, S., Zhang, B. & Qin, F. Arabidopsis RZFP34/CHYR1, a ubiquitin E3 ligase, regulates stomatal movement and drought tolerance via SnRK2.6-mediated phosphorylation. Plant Cell 27, 3228-3244 (2015).

[69] Golldack, D., Lüking, I. & Yang, O. Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep. 30, 1383-1391 (2011).

[70] Li, S., Fu, Q., Chen, L., Huang, W. & Yu, D. Arabidopsis thaliana WRKY25, WRKY26, and WRKY33 coordinate induction of plant thermotolerance. Planta 233, 1237-1252 (2011).

[71] Ouyang, S. Q. et al. Receptor‐like kinase OsSIK1 improves drought and salt stress tolerance in rice (Oryza sativa) plants. Plant J. 62, 316-329 (2010).

[72] Xie, K., Wu, C. & Xiong, L. Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice. Plant Physiol. 142, 280-293 (2006).

[73] Song, C. et al. miRNA and degradome sequencing reveal miRNA and their target genes that may mediate shoot growth in spur type mutant “Yanfu 6”. Front. Plant Sci. 8, 441 (2017).

[74] Jung, J.-H., Lee, H.-J., Ryu, J. Y. & Park, C.-M. SPL3/4/5 integrate developmental aging and photoperiodic signals into the FT-FD module in Arabidopsis flowering. Mol. Plant 9, 1647-1659 (2016).

[75] Padmanabhan, M. S. et al. Novel positive regulatory role for the SPL6 transcription factor in the N TIR-NB-LRR receptor-mediated plant innate immunity. PLoS Pathog. 9, e1003235 (2013).

[76] Xu, P. et al. Virus infection improves drought tolerance. New Phytol. 180, 911-921 (2008).

[77] Arshad, M., Feyissa, B. A., Amyot, L., Aung, B. & Hannoufa, A. MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13. Plant Sci. 258, 122-136 (2017).