Identification of ta-siRNAs and cis-nat-siRNAs in cassava and their roles in response to cassava bacterial blight

High throughput sRNA sequencing revealed gene regulatory role mediated by pathogen-derived small RNAs during Sri Lankan Cassava Mosaic Virus infection in Cassava

Small RNA (sRNA) mediated gene regulation during Sri Lankan Cassava Mosaic Virus (SLCMV) infection was studied from the Indian Cassava Cultivar H226. Our study generated high throughput sRNA dataset of 23.64 million reads from the control and SLCMV infected H226 leaf libraries. mes-miR9386 was detected as the most prominent miRNA expressed in control and infected leaf. Among the differentially expressed miRNAs, mes-miR156, mes- miR395 and mes-miR535a/b showed significant down regulation in the infected leaf. Genome-wide analysis of the three small RNA profiles revealed critical role of virus-derived small RNAs (vsRNAs) from the infected leaf tissues of H226. The vsRNAs were mapped to the bipartite SLCMV genome and high expression of siRNAs generated from the virus genomic region encoding AV1/AV2 genes in the infected leaf pointed towards the susceptibility of H226 cultivars to SLCMV. Furthermore, the sRNA reads mapped to the antisense strand of the SLCMV ORFs was higher than the sense strand. These vsRNAs were potential to target key host genes involved in virus interaction such as aldehyde dehydrogenase, ADP-ribosylation factor1 and ARF1-like GTP-binding proteins. The sRNAome-assisted analysis also revealed the origin of virus-encoded miRNAs from the SLCMV genome in the infected leaf. These virus-derived miRNAs were predicted to have hair-pin like secondary structures, and have different isoforms. Moreover, our study revealed that the pathogen sRNAs play a critical role in the infection process in H226 plants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03494-2.

The future of cassava in the era of biotechnology in Southern Africa

Cassava (Manihot esculenta) is a major staple food and the world's fourth source of calories. Biotechnological contributions to enhancing this crop, its advances, and present issues must be assessed regularly. Functional genomics, genomic-assisted breeding, molecular tools, and genome editing technologies, among other biotechnological approaches, have helped improve the potential of economically important crops like cassava by addressing some of its significant constraints, such as nutrient deficiency, toxicity, poor starch quality, disease susceptibility, low yield capacity, and postharvest deterioration. However, the development, improvement, and subsequent acceptance of the improved cultivars have been challenging and have required holistic approaches to solving them. This article provides an update of trends and gaps in cassava biotechnology, reviewing the relevant strategies used to improve cassava crops and highlighting the potential risk and acceptability of improved cultivars in Southern Africa.

Recent Advances in Plant Gene Silencing Methods

With the increasing understanding of fundamentals of gene silencing pathways in plants, various tools and techniques for downregulating the expression of a target gene have been developed across multiple plant species. This chapter provides an insight into the molecular mechanisms of gene silencing and highlights the advancements in various gene silencing approaches. The prominent aspects of different gene silencing methods, their advantages and disadvantages have been discussed. A succinct discussion on the newly emerged microRNA-based technologies like microRNA-induced gene silencing (MIGS) and microRNA-mediated virus-induced gene silencing (MIR-VIGS) are also presented. We have also discussed the gene-editing system like CRISPR-Cas. The prominent bottlenecks in gene silencing methods are the off-target effects and lack of universal applicability. However, the tremendous growth in understanding of this field reflects the potentials for improvements in the currently available approaches and the development of new widely applicable methods for easy, fast, and efficient functional characterization of plant genes.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Molecular identification of GAPDHs in cassava highlights the antagonism of MeGAPCs and MeATG8s in plant disease resistance against cassava bacterial blight

KEY MESSAGE

MeGAPCs were identified as negative regulators of plant disease resistance, and the interaction of MeGAPCs and MeATG8s was highlighted in plant defense response. As an important enzyme of glycolysis metabolic pathway, glyceraldehyde-3-P dehydrogenase (GAPDH) plays important roles in plant development, abiotic stress and immune responses. Cassava (Manihot esculenta) is most important tropical crop and one of the major food crops, however, no information is available about GAPDH gene family in cassava. In this study, 14 MeGAPDHs including 6 cytosol GAPDHs (MeGAPCs) were identified from cassava, and the transcripts of 14 MeGAPDHs in response to Xanthomonas axonopodis pv manihotis (Xam) indicated their possible involvement in immune responses. Further investigation showed that MeGAPCs are negative regulators of disease resistance against Xam. Through transient expression in Nicotiana benthamiana, we found that overexpression of MeGAPCs led to decreased disease resistance against Xam. On the contrary, MeGAPCs-silenced cassava plants through virus-induced gene silencing (VIGS) conferred improved disease resistance. Notably, MeGAPCs physically interacted with autophagy-related protein 8b (MeATG8b) and MeATG8e and inhibited autophagic activity. Moreover, MeATG8b and MeATG8e negatively regulated the activities of NAD-dependent MeGAPDHs, and are involved in MeGAPCs-mediated disease resistance. Taken together, this study highlights the involvement of MeGAPCs in plant disease resistance, through interacting with MeATG8b and MeATG8e.

Plant Responses to Pathogen Attack: Small RNAs in Focus

Small RNAs (sRNA) are a significant group of gene expression regulators for multiple biological processes in eukaryotes. In plants, many sRNA silencing pathways produce extensive array of sRNAs with specialized roles. The evidence on record advocates for the functions of sRNAs during plant microbe interactions. Host sRNAs are reckoned as mandatory elements of plant defense. sRNAs involved in plant defense processes via different pathways include both short interfering RNA (siRNA) and microRNA (miRNA) that actively regulate immunity in response to pathogenic attack via tackling pathogen-associated molecular patterns (PAMPs) and other effectors. In response to pathogen attack, plants protect themselves with the help of sRNA-dependent immune systems. That sRNA-mediated plant defense responses play a role during infections is an established fact. However, the regulations of several sRNAs still need extensive research. In this review, we discussed the topical advancements and findings relevant to pathogen attack and plant defense mediated by sRNAs. We attempted to point out diverse sRNAs as key defenders in plant systems. It is hoped that sRNAs would be exploited as a mainstream player to achieve food security by tackling different plant diseases.

MeWRKY20 and its interacting and activating autophagy-related protein 8 (MeATG8) regulate plant disease resistance in cassava

As a highly conserved mechanism, autophagy is responsible for the transport of cytoplasmic constituents in the vacuoles or lysosomes. Moreover, autophagy is essential for plant development and various stress responses. In this study, 34 MeATGs were systematically identified in cassava, and their transcripts were commonly regulated by Xanthomonas axonopodis pv manihotis (Xam). Through transient expression in Nicotiana benthamiana, the subcellular locations of 4 MeATG8s were revealed. Notably, MeWRKY20 was identified as physical interacting protein of MeATG8a/8f/8h and upstream transcriptional activator of MeATG8a. Through virus-induced gene silencing (VIGS) in cassava, we found that MeATG8-silenced and MeWRKY20-silenced plants resulted in disease sensitive, with less callose depositions and lower autophagic activity. This study may facilitate our understanding of the upstream MeWRKY20 and underlying target as well as interacting proteins of MeATG8s in immune response. Taken together, MeWRKY20 and MeATG8a/8f/8h are essential for disease resistance against bacterial blight by forming various transcriptional modules and interacting complex in cassava.

In silico analysis of the grapefruit sRNAome, transcriptome and gene regulation in response to CTV-CDVd co-infection

BACKGROUND

Small RNA (sRNA) associated gene regulation has been shown to play a significant role during plant-pathogen interaction. In commercial citrus orchards co-infection of Citrus tristeza virus (CTV) and viroids occur naturally.

METHODS

A next-generation sequencing-based approach was used to study the sRNA and transcriptional response in grapefruit to the co-infection of CTV and Citrus dwarfing viroid.

RESULTS

The co-infection resulted in a difference in the expression of a number of sRNA species when comparing healthy and infected plants; the majority of these were derived from transcripts processed in a phased manner. Several RNA transcripts were also differentially expressed, including transcripts derived from two genes, predicted to be under the regulation of sRNAs. These genes are involved in plant hormone systems; one in the abscisic acid, and the other in the cytokinin regulatory pathway. Additional analysis of virus- and viroid-derived small-interfering RNAs (siRNAs) showed areas on the pathogen genomes associated with increased siRNA synthesis. Most interestingly, the starting position of the p23 silencing suppressor's sub-genomic RNA generated a siRNA hotspot on the CTV genome.

CONCLUSIONS

This study showed the involvement of various genes, as well as endogenous and exogenous RNA-derived sRNA species in the plant-defence response. The results highlighted the role of sRNA-directed plant hormone regulation during biotic stress, as well as a counter-response of plants to virus suppressors of RNA-silencing.

In Silico Identification and Validation of Potential microRNAs in Kinnow Mandarin (Citrus reticulata Blanco)

MicroRNAs (miRNAs) are a large family of 19-25 nucleotides, regulatory, non-coding RNA molecules that control gene expression by cleaving or inhibiting the translation of target gene transcripts in animals and plants. Despite the important functions of miRNAs related to regulation of plant growth and development processes, metabolism, and abiotic and biotic stresses, little is known about the disease-related miRNA. Here, we present a new pipeline for miRNA analysis using expressed sequence tags (ESTs)-based bioinformatics approach in Kinnow mandarin, a commercially important citrus fruit crop. For this, 56,041 raw EST sequences of Citrus reticulata Blanco were retrieved from EST database in NCBI through step-by-step filtering and processing methods and 130 miRNAs were predicted. Upon blast with Citrus sinensis transcriptome data, these produced potential targets related to disease resistance proteins, pectin lyase-like superfamily proteins, lateral organ boundaries (LOB) domain-containing proteins 11, and protein phosphatase 2C family proteins, protein kinases, dehydrogenases, and methyltransferases. Majority of the predicted miRNAs were of 22, 23, and 24 nucleotides in length. To validate these computationally predicted miRNA, poly(A)-tailed Reverse Transcription-PCR was applied to detect the expression of seven miRNA which showed disease-related potential targets, in citrus greening diseased leaf tissues in comparison to the healthy tissues of Kinnow mandarin. Our study provides information on regulatory roles of these potential miRNAs for the citrus greening disease development, miRNA targets, and would be helpful for future research of miRNA function in citrus.

Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava

Cold and drought stresses seriously affect cassava (Manihot esculenta) plant growth and yield. Recently, long noncoding RNAs (lncRNAs) have emerged as key regulators of diverse cellular processes in mammals and plants. To date, no systematic screening of lncRNAs under abiotic stress and their regulatory roles in cassava has been reported. In this study, we present the first reference catalog of 682 high-confidence lncRNAs based on analysis of strand-specific RNA-seq data from cassava shoot apices and young leaves under cold, drought stress and control conditions. Among them, 16 lncRNAs were identified as putative target mimics of cassava known miRNAs. Additionally, by comparing with small RNA-seq data, we found 42 lncNATs and sense gene pairs can generate nat-siRNAs. We identified 318 lncRNAs responsive to cold and/or drought stress, which were typically co-expressed concordantly or discordantly with their neighboring genes. Trans-regulatory network analysis suggested that many lncRNAs were associated with hormone signal transduction, secondary metabolites biosynthesis, and sucrose metabolism pathway. The study provides an opportunity for future computational and experimental studies to uncover the functions of lncRNAs in cassava.

Two Cassava Basic Leucine Zipper (bZIP) Transcription Factors (MebZIP3 and MebZIP5) Confer Disease Resistance against Cassava Bacterial Blight

Basic domain-leucine zipper (bZIP) transcription factor, one type of conserved gene family, plays an important role in plant development and stress responses. Although 77 MebZIPs have been genome-wide identified in cassava, their in vivo roles remain unknown. In this study, we analyzed the expression pattern and the function of two MebZIPs (MebZIP3 and MebZIP5) in response to pathogen infection. Gene expression analysis indicated that MebZIP3 and MebZIP5 were commonly regulated by flg22, Xanthomonas axonopodis pv. manihotis (Xam), salicylic acid (SA), and hydrogen peroxide (H₂O₂). Subcellular localization analysis showed that MebZIP3 and MebZIP5 are specifically located in cell nucleus. Through overexpression in tobacco, we found that MebZIP3 and MebZIP5 conferred improved disease resistance against cassava bacterial blight, with more callose depositions. On the contrary, MebZIP3- and MebZIP5-silenced plants by virus-induced gene silencing (VIGS) showed disease sensitive phenotype, lower transcript levels of defense-related genes and less callose depositions. Taken together, this study highlights the positive role of MebZIP3 and MebZIP5 in disease resistance against cassava bacterial blight for further utilization in genetic improvement of cassava disease resistance.

High-resolution identification and abundance profiling of cassava (Manihot esculenta Crantz) microRNAs

BACKGROUND

Small RNAs (sRNAs) are endogenous sRNAs that play regulatory roles in plant growth, development, and biotic and abiotic stress responses. In plants, one subset of sRNAs, microRNAs (miRNAs) exhibit tissue-differential expression and regulate gene expression mainly through direct cleavage of mRNA or indirectly via production of secondary phased siRNAs (phasiRNAs) that silence cognate target transcripts in trans.

RESULTS

Here, we have identified cassava (Manihot esculenta Crantz) miRNAs using high resolution sequencing of sRNA libraries from leaf, stem, callus, male and female flower tissues. To analyze the data, we built a cassava genome database and, via sequence analysis and secondary structure prediction, 38 miRNAs not previously reported in cassava were identified. These new cassava miRNAs included two miRNAs not previously been reported in any plant species. The miRNAs exhibited tissue-differential accumulation as confirmed by quantitative RT-PCR and Northern blot analysis, largely reflecting levels observed in sequencing data. Some of the miRNAs identified were predicted to trigger production of secondary phased siRNAs (phasiRNAs) from 80 PHAS loci.

CONCLUSIONS

Cassava is a woody perennial shrub, grown principally for its starch-rich storage roots, which are rich in calories. In this study, new miRNAs were identified and their expression was validated using qRT-PCR of RNA from five different tissues. The data obtained expand the list of annotated miRNAs and provide additional new resources for cassava improvement research.

Virus tolerance and recovery from viral induced-symptoms in plants are associated with transcriptome reprograming

Plant recovery from viral infection is characterized by initial severe systemic symptoms which progressively decrease, leading to reduced symptoms or symptomless leaves at the apices. A key feature to plant recovery from invading nucleic acids such as viruses is the degree of the host's initial basal immunity response. We review current links between RNA silencing, recovery and tolerance, and present a model in which, in addition to regulation of resistance (R) and other defence-related genes by RNA silencing, viral infections incite perturbations of the host physiological state that trigger reprogramming of host responses to by-pass severe symptom development, leading to partial or complete recovery. Recovery, in particular in perennial hosts, may trigger tolerance or virus accommodation. We discuss evidence suggesting that plant viruses can avoid total clearance but persistently replicate at low levels, thereby modulating the host transcriptome response which minimizes fitness cost and triggers recovery from viral-symptoms. In some cases a susceptible host may fail to recover from initial viral systemic symptoms, yet, accommodates the persistent virus throughout the life span, a phenomenon herein referred to as non-recovery accommodation, which differs from tolerance in that there is no distinct recovery phase, and differs from susceptibility in that the host is not killed. Recent advances in plant recovery from virus-induced symptoms involving host transcriptome reprogramming are discussed.

Viroids, infectious long non-coding RNAs with autonomous replication

Transcriptome deep-sequencing studies performed during the last years confirmed that the vast majority of the RNAs transcribed in higher organisms correspond to several types of non-coding RNAs including long non-coding RNAs (lncRNAs). The study of lncRNAs and the identification of their functions, is still an emerging field in plants but the characterization of some of them indicate that they play an important role in crucial regulatory processes like flowering regulation, and responses to abiotic stress and plant hormones. A second group of lncRNAs present in plants is formed by viroids, exogenous infectious subviral plant pathogens well known since many years. Viroids are composed of circular RNA genomes without protein-coding capacity and subvert enzymatic activities of their hosts to complete its own biological cycle. Different aspects of viroid biology and viroid-host interactions have been elucidated in the last years and some of them are the main topic of this review together with the analysis of the state-of-the-art about the growing field of endogenous lncRNAs in plants.

Manipulating microRNAs for improved biomass and biofuels from plant feedstocks

Petroleum-based fuels are nonrenewable and unsustainable. Renewable sources of energy, such as lignocellulosic biofuels and plant metabolite-based drop-in fuels, can offset fossil fuel use and reverse environmental degradation through carbon sequestration. Despite these benefits, the lignocellulosic biofuels industry still faces many challenges, including the availability of economically viable crop plants. Cell wall recalcitrance is a major economic barrier for lignocellulosic biofuels production from biomass crops. Sustainability and biomass yield are two additional, yet interrelated, foci for biomass crop improvement. Many scientists are searching for solutions to these problems within biomass crop genomes. MicroRNAs (miRNAs) are involved in almost all biological and metabolic process in plants including plant development, cell wall biosynthesis and plant stress responses. Because of the broad functions of their targets (e.g. auxin response factors), the alteration of plant miRNA expression often results in pleiotropic effects. A specific miRNA usually regulates a biologically relevant bioenergy trait. For example, relatively low miR156 overexpression leads to a transgenic feedstock with enhanced biomass and decreased recalcitrance. miRNAs have been overexpressed in dedicated bioenergy feedstocks such as poplar and switchgrass yielding promising results for lignin reduction, increased plant biomass, the timing of flowering and response to harsh environments. In this review, we present the status of miRNA-related research in several major biofuel crops and relevant model plants. We critically assess published research and suggest next steps for miRNA manipulation in feedstocks for increased biomass and sustainability for biofuels and bioproducts.

Small RNA sequencing identifies miRNA roles in ovule and fibre development

MicroRNAs (miRNAs) have been found to be differentially expressed during cotton fibre development. However, which specific miRNAs and how they are involved in fibre development is unclear. Here, using deep sequencing, 65 conserved miRNA families were identified and 32 families were differentially expressed between leaf and ovule. At least 40 miRNAs were either leaf or ovule specific, whereas 62 miRNAs were shared in both leaf and ovule. qRT-PCR confirmed these miRNAs were differentially expressed during fibre early development. A total of 820 genes were potentially targeted by the identified miRNAs, whose functions are involved in a series of biological processes including fibre development, metabolism and signal transduction. Many predicted miRNA-target pairs were subsequently validated by degradome sequencing analysis. GO and KEGG analyses showed that the identified miRNAs and their targets were classified to 1027 GO terms including 568 biological processes, 324 molecular functions and 135 cellular components and were enriched to 78 KEGG pathways. At least seven unique miRNAs participate in trichome regulatory interaction network. Eleven trans-acting siRNA (tasiRNA) candidate genes were also identified in cotton. One has never been found in other plant species and two of them were derived from MYB and ARF, both of which play important roles in cotton fibre development. Sixteen genes were predicted to be tasiRNA targets, including sucrose synthase and MYB2. Together, this study discovered new miRNAs in cotton and offered evidences that miRNAs play important roles in cotton ovule/fibre development. The identification of tasiRNA genes and their targets broadens our understanding of the complicated regulatory mechanism of miRNAs in cotton.

Endogenous small-noncoding RNAs and their roles in chilling response and stress acclimation in Cassava

BACKGROUND

Small noncoding RNA (sncRNA), including microRNAs (miRNAs) and endogenous small-interfering RNAs (endo-siRNAs) are key gene regulators in eukaryotes, playing critical roles in plant development and stress tolerance. Trans-acting siRNAs (ta-siRNAs), which are secondary siRNAs triggered by miRNAs, and siRNAs from natural antisense transcripts (nat-siRNAs) are two well-studied classes of endo-siRNAs.

RESULTS

In order to understand sncRNAs' roles in plant chilling response and stress acclimation, we performed a comprehensive study of miRNAs and endo-siRNAs in Cassava (Manihot esculenta), a major source of food for the world populations in tropical regions. Combining Next-Generation sequencing and computational and experimental analyses, we profiled and characterized sncRNA species and mRNA genes from the plants that experienced severe and moderate chilling stresses, that underwent further severe chilling stress after chilling acclimation at moderate stress, and that grew under the normal condition. We also included castor bean (Ricinus communis) in our study to understand conservation of sncRNAs. In addition to known miRNAs, we identified 32 (22 and 10) novel miRNAs as well as 47 (26 and 21) putative secondary siRNA-yielding and 8 (7 and 1) nat-siRNA-yielding candidate loci in Cassava and castor bean, respectively. Among the expressed sncRNAs, 114 miRNAs, 12 ta-siRNAs and 2 nat-siRNAs showed significant expression changes under chilling stresses.

CONCLUSION

Systematic and computational analysis of microRNAome and experimental validation collectively showed that miRNAs, ta-siRNAs, and possibly nat-siRNAs play important roles in chilling response and chilling acclimation in Cassava by regulating stress-related pathways, e.g. Auxin signal transduction. The conservation of these sncRNA might shed lights on the role of sncRNA-mediated pathways affected by chilling stress and stress acclimation in Euphorbiaceous plants.

De novo transcriptome assembly from inflorescence of Orchis italica: analysis of coding and non-coding transcripts

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.

Epigenetic control of mobile DNA as an interface between experience and genome change

Mobile DNA in the genome is subject to RNA-targeted epigenetic control. This control regulates the activity of transposons, retrotransposons and genomic proviruses. Many different life history experiences alter the activities of mobile DNA and the expression of genetic loci regulated by nearby insertions. The same experiences induce alterations in epigenetic formatting and lead to trans-generational modifications of genome expression and stability. These observations lead to the hypothesis that epigenetic formatting directed by non-coding RNA provides a molecular interface between life history events and genome alteration.

Extending the sRNAome of apple by next-generation sequencing

The global importance of apple as a fruit crop necessitates investigations into molecular aspects of the processes that influence fruit quality and yield, including plant development, fruit ripening and disease resistance. In order to study and understand biological processes it is essential to recognise the range of molecules, which influence these processes. Small non-coding RNAs are regulatory agents involved in diverse plant activities, ranging from development to stress response. The occurrence of these molecules in apple leaves was studied by means of next-generation sequencing. 85 novel microRNA (miRNA) gene loci were predicted and characterized along with known miRNA loci. Both cis- and trans-natural antisense transcript pairs were identified. Although the trans-overlapping regions were enriched in small RNA (sRNA) production, cis-overlaps did not seem to agree. More than 150 phased regions were also identified, and for a small subset of these, potential miRNAs that could initiate phasing, were revealed. Repeat-associated siRNAs, which are generated from repetitive genomic regions such as transposons, were also analysed. For this group almost all available repeat sequences, associated with the apple genome and present in Repbase, were found to produce siRNAs. Results from this study extend our current knowledge on apple sRNAs and their precursors significantly. A rich molecular resource has been created and is available to the research community to serve as a baseline for future studies.