Functional identification of lncRNAs in sweet cherry (Prunus avium) pollen tubes via transcriptome analysis using single-molecule long-read sequencing

Abstract

Sweet cherry (Prunus avium) is a popular fruit with high nutritional value and excellent flavor. Although pollen plays an important role in the double fertilization and subsequent fruit production of this species, little is known about its pollen tube transcriptome. In this study, we identified 16,409 transcripts using single-molecule sequencing. After filtering 292 transposable elements, we conducted further analyses including mRNA classification, gene function prediction, alternative splicing (AS) analysis, and long noncoding RNA (lncRNA) identification to gain insight into the pollen transcriptome. The filtered transcripts could be matched with 3,438 coding region sequences from the sweet cherry genome. GO and KEGG analyses revealed complex biological processes during pollen tube elongation. A total of 2043 AS events were predicted, 7 of which were identified in different organs, such as the leaf, pistil and pollen tube. Using BLASTnt and the Coding-Potential Assessment Tool (CPAT), we distinguished a total of 284 lncRNAs, among which 154 qualified as natural antisense transcripts (NATs). As the NATs could be the reverse complements of coding mRNA sequences, they might bind to coding sequences. Antisense transfection assays showed that the NATs could regulate the expression levels of their complementary sequences and even affect the growth conditions of pollen tubes. In summary, this research characterizes the transcripts of P. avium pollen and lays the foundation for elucidating the physiological and biochemical mechanisms underlying sexual reproduction in the male gametes of this species.

References

[1] Williams, J. H. & Reese, J. B. Evolution of development of pollen performance. Curr. Topics Dev. Biol. 131, 299-336 (2019).

[2] Guan, Y., Guo, J., Li, H. & Yang, Z. Signaling in pollen tube growth: crosstalk, feedback, and missing links. Mol. Plant 6, 1053-1064 (2013).

[3] Zhou, Z. et al. Arabidopsis RIC1 severs actin filaments at the apex to regulate pollen tube growth. Plant Cell 27, 1140-1161 (2015).

[4] Loraine, A. E., McCormick, S., Estrada, A., Patel, K. & Qin, P. RNA-seq of Arabidopsis pollen uncovers novel transcription and alternative splicing. Plant Physiol. 162, 1092-1109 (2013).

[5] Vogel, H., Badapanda, C., Knorr, E. & Vilcinskas, A. RNA-sequencing analysis reveals abundant developmental stage-specific and immunity-related genes in the pollen beetle Meligethes aeneus. Insect Mol. Biol. 23, 98-112 (2014).

[6] Huang, L. et al. Systematic identification of long non-coding RNAs during pollen development and fertilization in Brassica rapa. Plant J. 96, 203-222 (2018).

[7] Han, Y. et al. Transcriptomic analysis of differentially expressed genes in flower buds of genetic male sterile and wild type cucumber by RNA sequencing. Physiol. Mol. Biol. Plants 24, 359-367 (2018).

[8] Wang, B. et al. Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat. Commun. 7, 11708 (2016).

[9] Correia, S., Schouten, R., Silva, A. P. & Goncalves, B. Factors affecting quality and health promoting compounds during growth and postharvest life of sweet cherry (Prunus avium L.). Front. Plant Sci. 8, 2166 (2017).

[10] Shirasawa, K. et al. The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding. DNA Res. 24, 499-508 (2017).

[11] Ono, K., Akagi, T., Morimoto, T., Wunsch, A. & Tao, R. Genome re-sequencing of diverse sweet cherry (Prunus avium) individuals reveals a modifier gene mutation conferring pollen-part self-compatibility. Plant Cell Physiol. 59, 1265-1275 (2018).

[12] Black, D. L. Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem. 72, 291-336 (2003).

[13] López-Bigas, N., Audit, B., Ouzounis, C., Parra, G. & Guigó, R. Are splicing mutations the most frequent cause of hereditary disease?. FEBS Lett. 579, 1900-1903 (2005).

[14] Zhang, G. et al. Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. Genome Res. 20, 646-654 (2010).

[15] Marquez, Y., Brown, J. W., Simpson, C., Barta, A. & Kalyna, M. Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res. 22, 1184-1195 (2012).

[16] Li, P. et al. The developmental dynamics of the maize leaf transcriptome. Nat. Genet. 42, 1060-1067 (2010).

[17] Lee, J. H. et al. Regulation of temperature-responsive flowering by MADS-box transcription factor repressors. Science 342, 628-632 (2013).

[18] Lin, F., Zhang, Y. & Jiang, M. Y. Alternative splicing and differential expression of two transcripts of nicotine adenine dinucleotide phosphate oxidase B gene from Zea mays. J. Integr. Plant Biol. 51, 287-298 (2009).

[19] Rinn, J. L. & Chang, H. Y. Genome regulation by long noncoding RNAs. Annu. Rev. Biochem. 81, 145-166 (2012).

[20] Brown, C. J. et al. The human XIST gene: analysis of a 17 kb inactive X-specific RNA that contains conserved repeats and is highly localized within the nucleus. Cell 71, 527-542 (1992).

[21] Liu, J., Wang, H. & Chua, N. H. Long noncoding RNA transcriptome of plants. Plant Biotechnol. J. 13, 319-328 (2015).

[22] Liu, T. T. et al. A global identification and analysis of small nucleolar RNAs and possible intermediate-sized non-coding RNAs in Oryza sativa. Mol. Plant 6, 830-846 (2013).

[23] Wang, Y. et al. Genomic features and regulatory roles of intermediate-sized non-coding RNAs in Arabidopsis. Mol. Plant 7, 514-527 (2014).

[24] Kung, J. T., Colognori, D. & Lee, J. T. Long noncoding RNAs: past, present, and future. Genetics 193, 651-669 (2013).

[25] Khalil, A. M. et al. Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc. Natl Acad. Sci. USA 106, 11667-11672 (2009).

[26] Bardou, F. et al. Long noncoding RNA modulates alternative splicing regulators in Arabidopsis. Dev. Cell 30, 166-176 (2014).

[27] Borsani, O., Zhu, J., Verslues, P. E., Sunkar, R. & Zhu, J. K. Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell 123, 1279-1291 (2005).

[28] Kim, D. H. & Sung, S. Vernalization-triggered intragenic chromatin loop formation by long noncoding RNAs. Dev. Cell 40, 302-312 (2017).

[29] Wang, Y. et al. Arabidopsis noncoding RNA mediates control of photomorphogenesis by red light. Proc. Natl Acad. Sci. USA 111, 10364 (2014).

[30] Cui, J., Luan, Y., Jiang, N., Bao, H. & Meng, J. Comparative transcriptome analysis between resistant and susceptible tomato allows the identification of lncRNA16397 conferring resistance to Phytophthora infestans by co-expressing glutaredoxin. Plant J. 89, 577-589 (2017).

[31] Qin, T., Zhao, H., Cui, P., Albesher, N. & Xiong, L. A nucleus-localized long non-coding RNA enhances drought and salt stress tolerance. Plant Physiol. 175, 1321-1336 (2017).

[32] Ding, J. et al. A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice. Proc. Natl Acad. Sci. USA 109, 2654-2659 (2012).

[33] Hackl, T., Hedrich, R., Schultz, J. & Förster, F. proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 30, 3004-3011 (2014).

[34] Jaime, H. et al. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, 286-293 (2016).

[35] Jia, Y. et al. WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update. Nucleic Acids Res. 46, 71-75 (2018).

[36] Wang, M. et al. A global survey of alternative splicing in allopolyploid cotton: landscape, complexity and regulation. New Phytol. 217, 163-178 (2018).

[37] Rogers, M. F., Thomas, J., Reddy, A. S. & Ben-Hur, A. SpliceGrapher: detecting patterns of alternative splicing from RNA-Seq data in the context of gene models and EST data. Genome Biol. 13, R4 (2012).

[38] Wang, L. et al. CPAT: coding-potential assessment tool using an alignment-free logistic regression model. Nucleic Acids Res. 41, e74 (2013).

[39] Zhao, X. et al. Global identification of Arabidopsis lncRNAs reveals the regulation of MAF4 by a natural antisense RNA. Nat. Commun. 9, 5056 (2018).

[40] Daccord, N. et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat. Genet. 49, 1099-1106 (2017).

[41] Wang, H. et al. Alternative splicing during Arabidopsis flower development results in constitutive and stage-regulated isoforms. Front. Genet. 5, 25 (2014).

[42] Thatcher, S. R. et al. Genome-wide analysis of alternative splicing during development and drought stress in maize. Plant Physiol. 170, 586-599 (2016).

[43] Mangano, S., Juárez, S. P. D. & Estevez, J. M. ROS regulation of polar growth in plant cells. Plant Physiol. 171, 1593-1605 (2016).

[44] Oxenoid, K. et al. Architecture of the mitochondrial calcium uniporter. Nature 533, 269-273 (2016).

[45] Fu, Y. The cytoskeleton in the pollen tube. Curr. Opin. Plant Biol. 28, 111-119 (2015).

[46] Lin, Y. et al. Comparative transcriptome profiling analysis of red- and white-fleshed strawberry (Fragaria × ananassa) provides new insight into the regulation of anthocyanins pathway. Plant Cell Physiol. 59, 1844-1859 (2018).

[47] Lv, Y. et al. Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize (Zea mays L.). BMC Genomics 17, 350 (2016).

[48] Wang, J. et al. Genome-wide identification and functional prediction of novel drought-responsive lncRNAs in Pyrus betulifolia. Genes 9, 311 (2018).

[49] Chekanova, J. A. Long non-coding RNAs and their functions in plants. Curr. Opin. Plant Biol. 27, 207-216 (2015).

[50] Lagarde, J. et al. High-throughput annotation of full-length long noncoding RNAs with capture long-read sequencing. Nat. Genet. 49, 1731-1740 (2017).

[51] González-Buendía, E., Saldaña-Meyer, R., Meier, K. & Recillas-Targa, F. Transcriptome-wide identification of in vivo interactions between RNAs and RNA-binding proteins by RIP and PAR-CLIP assays. Methods Mol. Biol. 1288, 413-428 (2015).

[52] Wu, T. et al. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

[53] Larkin, M. et al. Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947-2948 (2007).

[54] Jin, W. et al. The R2R3 MYB transcription factor PavMYB10.1 involves in anthocyanin biosynthesis and determines fruit skin colour in sweet cherry (Prunus avium L.). Plant Biotechnol. J. 14, 2120-2133 (2016).

[55] Stephen, F. et al. Basic local alignment search tool. J. Mol. Biol. 215, 403-410 (1990).

[56] Altschul, S. et al. Issues in searching molecular sequence databases. Nat. Genet. 6, 119-129 (1994).

[57] Olson, S. EMBOSS opens up sequence analysis. Brief. Bioinform. 3, 87-91 (2002).

[58] Li, W. et al. Apple S-RNase triggers inhibition of tRNA aminoacylation by interacting with a soluble inorganic pyrophosphatase in growing self-pollen tubes in vitro. New Phytol. 218, 579-593 (2018).