Identification of Known and Novel microRNAs and Their Targets in Peach (Prunus persica) Fruit by High-Throughput Sequencing

Characterization of ZAT12 protein from Prunus persica: role in fruit chilling injury tolerance and identification of gene targets

MAIN CONCLUSION

PpZAT12, a transcription factor differentially expressed in peach varieties with distinct susceptibility tochilling injury (CI), is a potential candidate gene for CI tolerance by regulating several identified gene targets. ZAT (zinc finger of Arabidopsis thaliana) proteins play roles in multiple abiotic stress tolerance in Arabidopsis and other species; however, there are few reports on these transcription factors (TFs) in fruit crops. This study aimed to evaluate PpZAT12, a C2H2 TF up-regulated in peach fruit by a heat treatment applied before postharvest cold storage for reducing chilling injury (CI) symptoms. Here, the expression of PpZAT12 in different tissues and fruits subjected to either postharvest heat or cold treatments, was evaluated in peach varieties with differential susceptibility to develop CI. PpZAT12 increased by cold storage in CI-resistant cultivars ('Elegant Lady' and 'Rojo 2'), while it was not modified in a cultivar susceptible to develop CI ('Flordaking'). Besides, we expressed PpZAT12 in Arabidopsis (35S::PpZAT12) and found that these plants show impaired plant growth and development, rendering small plants with senescence delay and aborted seeds. We applied a proteomic approach to decipher the peptides responding to PpZAT12 in Arabidopsis and found 348 differential expressed proteins (DEPs) relative to the wild type. Besides, comparing the DEPs between Arabidopsis plants expressing PpZAT12 or AtZAT12 (35S::AtZAT12) we found common and specific responses to these TFs. Based on the proteomic information obtained here and published data about AtZAT12, we searched ZAT12-targets in peach allowing the identification of a putative ZAT12 regulon in this species. The identified peach ZAT12-protein targets could underlie the differential susceptibility to CI among different peach varieties and can be used as future targets to improve adaptation to refrigeration in fleshy fruits.

The PpPep2-Triggered PTI-like Response in Peach Trees Is Mediated by miRNAs

Plant diseases diminish crop yields and put the world's food supply at risk. Plant elicitor peptides (Peps) are innate danger signals inducing defense responses both naturally and after external application onto plants. Pep-triggered defense networks are compatible with pattern-triggered immunity (PTI). Nevertheless, in complex regulatory pathways, there is crosstalk among different signaling pathways, involving noncoding RNAs in the natural response to pathogen attack. Here, we used Prunus persica, PpPep2 and a miRNA-Seq approach to show for the first time that Peps regulate, in parallel with a set of protein-coding genes, a set of plant miRNAs (~15%). Some PpPep2-regulated miRNAs have been described to participate in the response to pathogens in various plant-pathogen systems. In addition, numerous predicted target mRNAs of PpPep2-regulated miRNAs are themselves regulated by PpPep2 in peach trees. As an example, peach miRNA156 and miRNA390 probably have a role in plant development regulation under stress conditions, while others, such as miRNA482 and miRNA395, would be involved in the regulation of resistance (R) genes and sulfate-mediated protection against oxygen free radicals, respectively. This adds to the established role of Peps in triggering plant defense systems by incorporating the miRNA regulatory network and to the possible use of Peps as sustainable phytosanitary products.

CRISPR/Cas14a-Based Isothermal Amplification for Profiling Plant MicroRNAs

MicroRNAs (miRNAs) play key roles in biological processes in plants, such as stress resistance, yet can hardly be quantified by an enzyme-involved terminal polymerization process due to their 2'-O-methyl modifications at the 3'-terminal. Herein, we proposed a CRISPR/Cas14a-based rolling circle amplification (termed Cas14R) assay, allowing reverse transcription-free and demethylation-free detection of plant miRNAs with single-nucleotide resolution. The employment of target-templated rolling circle amplification circumvents the extension of the unaccessible 2'-O-methyl group at the 3'-terminal. Particularly, the activated Cas14a confers the trans-cleavage activity for identifying target single-stranded DNA sequences without the necessity of the protospacer adjacent motif, generalizing the detection of miRNA sequences and the integration of different isothermal amplification techniques. Ultimately, the Cas14R assay has been applied to profile miR156a to evaluate the ripeness process of banana, indicating its feasibility in analyzing the roles of miRNAs in biological processes of plants.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

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

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.

Genome-Wide Changes of Regulatory Non-Coding RNAs Reveal Pollen Development Initiated at Ecodormancy in Peach

Bud dormancy is under the regulation of complex mechanisms including genetic and epigenetic factors. To study the function of regulatory non-coding RNAs in winter dormancy release, we analyzed the small RNA and long non-coding RNA (lncRNA) expression from peach (Prunus persica) floral buds in endodormancy, ecodormancy and bud break stages. Small RNAs underwent a major shift in expression primarily between dormancy and flowering with specific pairs of microRNAs and their mRNA target genes undergoing coordinated differential expression. From endodormancy to ecodormancy, ppe-miR6285 was significantly upregulated while its target gene, an ASPARAGINE-RICH PROTEIN involved in the regulation of abscisic acid signaling, was downregulated. At ecodormancy, ppe-miR2275, a homolog of meiosis-specific miR2275 across angiosperms, was significantly upregulated, supporting microsporogenesis in anthers at a late stage of dormancy. The expression of 785 lncRNAs, unlike the overall expression pattern in the small RNAs, demonstrated distinctive expression signatures across all dormancy and flowering stages. We predicted that a subset of lncRNAs were targets of microRNAs and found 18 lncRNA/microRNA target pairs with both differentially expressed across time points. The genome-wide differential expression and network analysis of non-coding RNAs and mRNAs from the same tissues provide new candidate loci for dormancy regulation and suggest complex noncoding RNA interactions control transcriptional regulation across these key developmental time points.

Genetic Transformation in Peach (Prunus persica L.): Challenges and Ways Forward

Almost 30 years have passed since the first publication reporting regeneration of transformed peach plants. Nevertheless, the general applicability of genetic transformation of this species has not yet been established. Many strategies have been tested in order to obtain an efficient peach transformation system. Despite the amount of time and the efforts invested, the lack of success has significantly limited the utility of peach as a model genetic system for trees, despite its relatively short generation time; small, high-quality genome; and well-studied genetic resources. Additionally, the absence of efficient genetic transformation protocols precludes the application of many biotechnological tools in peach breeding programs. In this review, we provide an overview of research on regeneration and genetic transformation in this species and summarize novel strategies and procedures aimed at producing transgenic peaches. Promising future approaches to develop a robust peach transformation system are discussed, focusing on the main bottlenecks to success including the low efficiency of A. tumefaciens-mediated transformation, the low level of correspondence between cells competent for transformation and those that have regenerative competence, and the high rate of chimerism in the few shoots that are produced following transformation.

miRNA-mRNA Integrated Analysis Reveals Roles for miRNAs in a Typical Halophyte, Reaumuria soongorica, during Seed Germination under Salt Stress

MicroRNAs (miRNAs) are endogenous small RNAs that play a crucial role in plant growth, development, and environmental stress responses. Reaumuria soongorica is a typical halophyte that is widely distributed in saline–alkali desert regions. Under salt stress, R. soongorica can complete germination, a critical biological process in the life cycle of seed plants. To identify miRNAs and predict target mRNAs involved in seed germination during salt stress, nine small-RNA libraries were constructed and analyzed from R. soongorica seeds treated with various concentrations of NaCl. We also obtained transcriptome data under the same treatment conditions. Further analysis identified 88 conserved miRNAs representing 25 defined families and discovered 13 novel miRNAs from nine libraries. A co-expression analysis was performed on the same samples to identify putative miRNA–mRNA interactions that were responsive to salt stress. A comparative analysis of expression during germination under 273 (threshold) and 43 mM (optimal) NaCl treatments identified 13 differentially expressed miRNAs and 23 corresponding target mRNAs, while a comparison between 43 mM NaCl and non-salt-stress conditions uncovered one differentially expressed miRNA and one corresponding target mRNA. These results provide basic data for further study of molecular mechanisms involved in the germination of salt-stressed R. soongorica seeds, and also provide a reference for the improvement of salt tolerance during plant germination.

MicroRNA858-mediated regulation of anthocyanin biosynthesis in kiwifruit (Actinidia arguta) based on small RNA sequencing

As important regulators, miRNAs could play pivotal roles in regulation of fruit coloring. Actinidia arguta is a newly emerged fruit tree with extensively application prospects. However, miRNAs involved in A. arguta fruit coloring are unknown. In this study, A. arguta fruit were investigated at three developmental stages by small RNAs high-throughput sequencing. A total of 482 conserved miRNAs corresponding to 526 pre-miRNAs and 581 novel miRNAs corresponding to 619 pre-miRNAs were grouped into 46 miRNA families. Target gene prediction and analysis revealed that miR858, a strongly candidate miRNA, was involved in anthocyanin biosynthesis in which contributes to fruit coloring. The anthocyanin level was determined in three A. arguta cultivars by UPLC-MS/MS (ultra-performance liquid chromatography coupled with tandem mass spectrometry). In addition, qPCR (quantitative real-time PCR), cluster analysis were conducted as well as correlation analysis. All results were combined to propose a model in which describes an association of miRNA and anthocyanin biosynthesis in A. arguta. The data presented herein is the first report on miRNA profile analysis in A. arguta, which can provide valuable information for further research into the regulation of the miRNAs in anthocyanin biosynthesis and fruit coloring.

Characterization of genome-wide microRNAs and their roles in development and biotic stress in pear

Using a genome-wide analysis of miRNAs in 'Yali' pear (Pyrus bretschneideri) via the next-generation high-throughput sequencing of small RNAs with a bioinformatics analysis, we found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 and their target genes function in viral defense in 'Duli' and 'Hongbaoshi'. pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319 and their targets function in auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Successful fruit production in pear (Pyrus spp.) depends on the use of optimal combinations of rootstocks and scions. Deciphering plant-pathogen defense mechanisms and hormone signaling pathways is an important step towards developing pear rootstocks and varieties with improved qualities. In the current study, we combined next-generation sequencing of small RNAs with a bioinformatics analysis to systematically identify and characterize 298 miRNAs in the pear scion cultivar 'Yali' (Pyrus bretschneideri). We also analyzed miRNAs in three rootstock varieties ('Duli', 'Zhongai 4', and 'Zhongai 5') and one scion cultivar ('Hongbaoshi'). We found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 are induced following infection with the pear virus Apple stem pitting virus (ASPV), and identified their target genes (pbRPS6, pbNAC, pbTLR, and pbRX-CC, respectively), which participate in viral defense pathways in 'Duli' and 'Hongbaoshi'. Furthermore, we identified pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319, and found that the production of these miRNAs was suppressed under low levels of synthetic auxin. The targets of these miRNAs (pbARF, pbAEC, pbSCL, and pbTCP4) respond to auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Our results lay the foundation for breeding improved pear cultivars.

Comparative Analysis of Fruit Ripening-Related miRNAs and Their Targets in Blueberry Using Small RNA and Degradome Sequencing

MicroRNAs (miRNAs) play vital roles in the regulation of fruit development and ripening. Blueberry is an important small berry fruit crop with economical and nutritional value. However, nothing is known about the miRNAs and their targets involved in blueberry fruit ripening. In this study, using high-throughput sequencing of small RNAs, 84 known miRNAs belonging to 28 families and 16 novel miRNAs were identified in white fruit (WF) and blue fruit (BF) libraries, which represent fruit ripening onset and in progress, respectively. Among them, 41 miRNAs were shown to be differentially expressed during fruit maturation, and 16 miRNAs representing 16 families were further chosen to validate the sRNA sequencing data by stem-loop qRT-PCR. Meanwhile, 178 targets were identified for 41 known and 7 novel miRNAs in WF and BF libraries using degradome sequencing, and targets of miR160 were validated using RLM-RACE (RNA Ligase-Mediated (RLM)-Rapid Amplification of cDNA Ends) approach. Moreover, the expression patterns of 6 miRNAs and their targets were examined during fruit development and ripening. Finally, integrative analysis of miRNAs and their targets revealed a complex miRNA-mRNA regulatory network involving a wide variety of biological processes. The findings will facilitate future investigations of the miRNA-mediated mechanisms that regulate fruit development and ripening in blueberry.

Discovery and Annotation of Plant Endogenous Target Mimicry Sequences from Public Transcriptome Libraries: A Case Study of Prunus persica

Novel transcript discovery through RNA sequencing has substantially improved our understanding of the transcriptome dynamics of biological systems. Endogenous target mimicry (eTM) transcripts, a novel class of regulatory molecules, bind to their target microRNAs (miRNAs) by base pairing and block their biological activity. The objective of this study was to provide a computational analysis framework for the prediction of putative eTM sequences in plants, and as an example, to discover previously un-annotated eTMs in Prunus persica (peach) transcriptome. Therefore, two public peach transcriptome libraries downloaded from Sequence Read Archive (SRA) and a previously published set of long non-coding RNAs (lncRNAs) were investigated with multi-step analysis pipeline, and 44 putative eTMs were found. Additionally, an eTM-miRNA-mRNA regulatory network module associated with peach fruit organ development was built via integration of the miRNA target information and predicted eTM-miRNA interactions. My findings suggest that one of the most widely expressed miRNA families among diverse plant species, miR156, might be potentially sponged by seven putative eTMs. Besides, the study indicates eTMs potentially play roles in the regulation of development processes in peach fruit via targeting specific miRNAs. In conclusion, by following the step-by step instructions provided in this study, novel eTMs can be identified and annotated effectively in public plant transcriptome libraries.

Identification and comparative profiling of miRNAs in herbaceous peony (Paeonia lactiflora Pall.) with red/yellow bicoloured flowers

Herbaceous peony (Paeonia lactiflora Pall.) is popular worldwide because of its gorgeous flower colour, and the yellow flower is the rarest. However, its mechanism of yellow formation is still unexplored from the post-translational level. In this study, the anatomy of the petal, cell sap pH and metal elements were investigated in bicoloured flower cultivar 'Jinhui' with red outer-petal and yellow inner-petal, and the yellow formation was influenced by the anatomy of petal, while not by the cell sap pH and metal elements. Subsequently, microRNAs sequencing (miRNA-seq) was used to identify small RNAs (sRNAs). A total of 4,172,810 and 3,565,152 specific unique sRNAs were obtained, 207 and 204 conserved miRNAs and 38 and 42 novel miRNAs were identified from red outer-petal and yellow inner-petal, respectively, which were confirmed by subcloning. Among these miRNAs, 163 conserved and 28 novel miRNAs were differentially expressed in two wheel of petals. And 5 differentially expressed miRNAs and their corresponding target genes related to yellow formation were screened, and their dynamic expression patterns confirmed that the yellow formation might be under the regulation of miR156e-3p-targeted squamosa promoter binding protein-like gene (SPL1). These results improve the understanding of miRNA regulation of the yellow formation in P. lactiflora.