Genome-wide analysis of long noncoding RNAs affecting floral bud dormancy in pears in response to cold stress

Abstract The versatile role of long noncoding RNAs (lncRNAs) in plant growth and development has been established, but a systematic identification and analysis of lncRNAs in the pear has not been reported. Bud dormancy is a crucial and complicated protective mechanism for plants in winter. The roles of lncRNAs in the dormancy process remain largely unclear. In this study, we induced pear floral buds to enter into different dormant statuses by simulating four different chilling accumulation conditions. Then, a time series of RNA-seq analysis was performed and we identified 7594 lncRNAs in Pyrus pyrifolia (Burm. F.) Nakai that have not been identified. The sequence and expression of the lncRNAs were confirmed by PCR analysis. In total, 6253 lncRNAs were predicted to target protein-coding genes including 692 cis-regulated pairs (596 lncRNAs) and 13,158 trans-regulated pairs (6181 lncRNAs). Gene Ontology analysis revealed that most of lncRNAs' target genes were involved in catalytic activity, metabolic processes and cellular processes. In the trend analysis, 124 long-term cold response lncRNAs and 80 short-term cold response lncRNAs were predicted. Regarding the lncRNA-miRNA regulatory networks, 59 lncRNAs were identified as potential precursors for miRNA members of 20 families, 586 lncRNAs were targets of 261 pear miRNAs and 53 lncRNAs were endogenous target mimics for 26 miRNAs. In addition, three cold response lncRNAs, two miRNAs and their target genes were selected for expression confirmed. The trend of their expression was consistent with the predicted relationships among them and suggested possible roles of lncRNAs in ABA metabolic pathway. Our findings not only suggest the potential roles of lncRNAs in regulating the dormancy of pear floral buds but also provide new insights into the lncRNA-miRNA-mRNA regulatory network in plants.