Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing

Daniel Gonzalez-Ibeas
José Blanca
Livia Donaire
Montserrat Saladié
Albert Mascarell-Creus
Ana Cano-Delgado
Jordi Garcia-Mas
Cesar Llave
Miguel A Aranda
Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing
R E S E A R C H A R T I C L E Open Access
Background: Melon (Cucumis melo L.) is a commercially important fruit crop that is cultivated worldwide. The melon research community has recently benefited from the determination of a complete draft genome sequence and the development of associated genomic tools, which have allowed us to focus on small RNAs (sRNAs). These are short, non-coding RNAs 21-24 nucleotides in length with diverse physiological roles. In plants, they regulate gene expression and heterochromatin assembly, and control protection against virus infection. Much remains to be learned about the role of sRNAs in melon. Results: We constructed 10 sRNA libraries from two stages of developing ovaries, fruits and photosynthetic cotyledons infected with viruses, and carried out high-throughput pyrosequencing. We catalogued and analysed the melon sRNAs, resulting in the identification of 26 known miRNA families (many conserved with other species), the prediction of 84 melon-specific miRNA candidates, the identification of trans-acting siRNAs, and the identification of chloroplast, mitochondrion and transposon-derived sRNAs. In silico analysis revealed more than 400 potential targets for the conserved and novel miRNAs. Conclusion: We have discovered and analysed a large number of conserved and melon-specific sRNAs, including miRNAs and their potential target genes. This provides insight into the composition and function of the melon small RNAome, and paves the way towards an understanding of sRNA-mediated processes that regulate melon fruit development and melon-virus interactions.

Background

Melon (Cucumis melo L., family Cucurbitaceae) is an important horticultural species cultivated in temperate, subtropical and tropical regions worldwide, with Spain being the largest producer in Europe and fifth in the world [1]. The melon genome has 12 chromosomes and is thought to contain 450-500 Mb of DNA, which is 3-4 times more than Arabidopsis [2]. Melon is a useful model for the analysis of fruit traits because of the vast morphological, physiological and biochemical diversity within the species, which can be exploited to dissect the biological processes controlling color, flavor and texture and how these properties arise during fruit development [3,4].

Despite the importance of melon, not much was available in the way of genomic sequence information prior to the establishment of a functional genomics consortium in 2004, which developed a range of tools and accumulated more than 33,000 expressed sequence tags (ESTs) and~17,000 tentative consensus sequences (unigenes) [5]. This EST collection has been expanded recently with the addition of 94,000 new ESTs from full-length enriched cDNA and standard cDNA libraries from various melon tissues and cultivars in the framework of the International Cucurbit Genome Initiative [6]. These ESTs as well as other resources are now accessible in a public database [7]. The unigene sequences have also been used to construct an oligonucleotide microarray, which has been applied in the analysis of fruit quality traits, ovary development and pathogen resistance [8]. In addition, a melon sequencing consortium has recently produced a high-quality draft of the melon genome (unpublished data). Although these resources provided significant advances in the analysis of melon gene expression, the small RNA (sRNAs) component of the melon transcriptome has not been studied in detail. These important molecules have been studied in other crop species and have been shown to fulfill a number of critical regulatory roles [9][10][11][12].

sRNAs are short, non-coding RNAs 21-24 nucleotides (nt) in length which are found in protists, fungi, plants and animals [13]. In plants, their roles include maintenance of genome stability, initiation of heterochromatin assembly, post-transcriptional regulation of gene expression and protection against viruses using an RNA-based immune system. The most abundant and best-characterised sRNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs). miRNAs are widely studied because of their regulatory activity, particularly in development, pathogen resistance and stress responses [13]. miRNAs are cleaved from stem-loop precursor molecules that derive from single stranded non-coding transcripts. miRNAs regulate protein-coding genes posttranscriptionally by mediating RNA cleavage or translational repression. Unlike miRNAs, siRNAs are generated from double-stranded RNA precursors and function on cognate RNA or DNA molecules by instigating degradation or promoting RNA-directed DNA methylation, respectively. cis-acting siRNAs (ca-siRNAs) arise from and target endogenous loci such as transposons and DNA repeats to direct cytosine methylation and chromatin modifications [14]. Natural antisense-transcript siRNAs (nat-siRNAs), which derive from pairs of natural-antisense transcripts, guide the cleavage of one of the two parent transcripts, leading to the production of a series of secondary 21-nt siRNAs of unclear function [15,16]. Finally, trans-acting siRNAs (ta-siRNAs) derived from TAS genes, which transcribe long primary noncoding RNAs as precursors for ta-siRNA biogenesis. TAS primary RNAs are cleaved by specific miRNAs and are sequentially processed into 21-nt ta-siRNAs starting from the miRNA-cleaved end, to generate clusters of phased siRNAs [17,18]. In addition to endogenous sRNAs, exogenous siRNAs from virus genomes can be detected in virus-infected plants as a part of the RNAbased immune system [19].

RNA viruses that infect melon are responsible for significant yield losses as well as poor fruit quality [20,21], particularly the widespread Watermelon mosaic virus (WMV, genus Potyvirus, family Potyviridae) [22,23]. Recently, a collection of accessions representing cultivated melon and its wild relatives was screened to identify sources of resistance to mosaic-inducing viruses [24]. TGR-1551 was identified as a resistant accession based on the lower WMV titer compared to susceptible genotypes (e.g. melon cv. Tendral) and the absence or mildness of the mosaic symptoms normally observed in systemically infected leaves [25]. Melon necrotic spot virus (MNSV, genus Carmovirus, family Tombusviridae), although less economically important, may also cause yield losses, and epidemic outbreaks have been reported worldwide [26,23]. In melon, resistance to MNSV is controlled by the single recessive gene nsv, which encodes eukaryotic translation initiation factor 4E (Cm-eIF4E) [27]. This resistance is effective against all MNSV strains (e.g. MNSV-Malfa5) except MNSV-264 [28]. Studies of chimeric viruses have shown that the MNSV 3' untranslated region (3'-UTR) contains the resistance-breaking determinant of MNSV-264, and that it functions as a cap-independent translational enhancer [29,30].

We constructed 10 sRNA libraries from a range of healthy and virus-infected melon tissues, and we sequenced a set of endogenous and exogenous sRNAs using the pyrosequencing-based 454 technology from Roche [31]. To gain insights into the role of sRNAs on key aspects of fruit development, maturation and pathogen defense, samples from two stages of the developing ovary, fruits 15 and 45 days after pollination, and photosynthetic cotyledons from resistant and susceptible melon accessions infected with WMV and MNSV were analysed. In a previous study, we reported the profile of virus-derived sRNAs (viRNAs) from cotyledon samples [32]. Here we report a catalog of endogenous melon sRNAs, including miRNAs from known families and new candidate miRNAs potentially unique to melon, focusing on the number of sequence reads as a reflection of their expression profiles. Potential targets for these miRNAs in the melon transcriptome were identified.

Results

cDNA libraries and sequencing of small RNAs

We used high throughput sequencing data to analyze the composition of the small RNA transcriptome (sRNAome) of melon and compare the results to data in publicly-available RNA and genomic databases. Ten sRNA libraries were constructed from total RNA extracted from fruits, ovaries and healthy and virusinfected melon cotyledons (Table 1). PCR amplification products corresponding to each library were pooled in equal amounts and sequences were obtained by multiplexed high-throughput pyrosequencing (Roche 454). This produced 447,180 raw sequences, each~100 bases in length, 432,743 of which had a complete 3' adaptor in the correct position. Based on these data, we estimated a sequencing error rate of 3.7%. After removing reads where one or the two adaptors could not be identified, 398,450 useful sequences with 3' and 5' adaptors were selected. Only 44 sequences comprising ligated adaptors without an insert were identified. Although we pooled similar amounts of PCR products from each library, different numbers of sequences were obtained according to the 5' adaptor sequence barcode (Table 1). For instance, the fruit and ovary libraries (15d, 45d, c1 and c5) were poorly represented providing a collection of fewer than one third of the number of sequences from the other six libraries. A set of 186,698 non-redundant sRNA sequences was generated for downstream analysis. The representation of sequences with different lengths in the redundant and non-redundant sRNAs datasets is shown in Figure 1. The most abundant sequences were 21, 24, 20 and 22 nts. A few sequences shorter than 20 nt were also retrieved, and these probably represent cloning artifacts and/or degradation products. Sequences > 30 nt in length in our dataset predominantly represented combinations of other melon sRNAs identified in our work. Detailed data are provided in Additional file 1.

Identification of known miRNAs

In order to identify known miRNAs, the melon sRNA data set was used as a BLAST query against the Arabidopsis small RNA database (ASRP) [33] and the microRNA database (miRbase) [34]. We identified 46 melon unique sequences corresponding to 26 miRNA families. Thirty nine sequences were identical to known miRNAs from other plant species, while 7 additional species were sequence variants highly conserved (up to two mismatches allowed). In order to clearly identify each melon sequence, melon miRNAs were named according to the homologous reference miRNA from each database ( Table 2). For each reference miRNA, we found that~3% of the corresponding melon sequences differed at one or two sites with mismatches distributed randomly along the sequence, so these were considered sequencing errors. Only specific sequence variants that represented more than 3% of the total population for each reference miRNA were considered biologically relevant. We identified only two non-conserved miRNAs, corresponding to ath-miR2111a from Arabidopsis and peu-miR2910 from Populus euphratica, respectively ( Table 2). The largest diversity of miRNA species was found in ovary samples and the lowest in fruit samples.   The abundance distribution of different miRNAs in each library was estimated based on sequencing frequencies as shown in Figure 2. We used sequencing data for quantitative profiling of small RNAs, though estimation of abundance based on sequencing frequencies could be misleading due to limited sequencing depth. Many miRNAs differed in abundance according to the source library. Nevertheless, most of the redundancy reflected the accumulation of miR159a, which accounted for more than 14,000 sequences in total.  genotypes Tendral and TGR-1551. Melon miRNA species with similarity to Arabidopsis miR156abcdef, miR160abc and miR168ab, which target mRNAs encoding squamosa promoter binding proteins, auxin response factors (ARFs) and argonaute-like proteins (AGO), respectively, showed different trends in the genotypes tested. For example, miR168ab is more abundant in healthy Tendral tissues compared to infected tissues whereas it is more abundant in WMV-infected TGR-1551 tissues than in healthy tissues. Other known miR-NAs in our sequenced set were generally more abundant in healthy tissues irrespective of the melon variety tested. For example, miRNAs with similarity to Arabidopsis miR159a and miR167d, which target MYB transcription factors and ARFs, respectively, followed this trend in both genotypes albeit with differences in magnitude. Comparison of the two libraries from ovary and fruit samples ( Figure 2C, D) revealed that miRNAs were particularly abundant and diverse in ovaries compared to fruits. Several miRNAs appeared to be temporally regulated during ovary development (e.g. members of the miR160, miR164, miR167, miR169, miR319 and miR390 families) whereas others were equally abundant at both ovary stages (miR156, miR167 and miR171 families). Fruits contained far fewer miRNAs than ovaries, and only miRNAs similar in sequence to Arabidopsis miR159a, miR164ab and miR397a showed significant differences in accumulation (with trends opposite to those seen in ovaries). These findings indicated that miRNAs in melon were expressed in specific tissues and in response to particular physiological conditions. In Arabidopsis, most of these miRNAs target mRNAs encoding transcription factors with roles in development, such as hormone signal transduction and organ identity. Figure 2E, F compares the accumulation of miRNAs in healthy and MNSV-infected tissues. Similar accumulation profiles were observed in both samples for most of the miRNAs identified. Exceptionally, miRNAs similar to Arabidopsis miR396a, miR396b and miR162a, which regulate transcripts encoding GRF transcription factors and DCL proteins, respectively, showed opposite accumulation patterns.

Identification of miRNA/miRNA* duplexes DCL-mediated cleavage of miRNA precursors having the characteristic stem-loop structure gives rise to miRNA duplexes where one of the two strands is the guide miRNA (the functional molecule) while the nearperfect complement sequence is known as the passenger miRNA, or miRNA*. The miRNA* is rapidly degraded but transient species can be cloned and therefore sequenced. We identified 16 miRNA* sequences complementary to some of the 46 miRNAs in our dataset ( Table 2), nine of which had the predicted sequence based on the fold-back structure of their presumptive precursors with internal mismatches and two additional terminal nucleotides forming a 3' tail ( Figure 3A), whereas the other six had a different number of protruding nucleotides and were considered non-typical ( Figure 3B).

The number of sequenced miRNA*s was generally much lower than the number of mature sequences but there were some remarkable exceptions. For example, for miR396a we counted 134 miRNA and 84 miRNA* sequences, as opposed to miR159a for which 14,651 miRNA sequences but no corresponding miRNA* sequences were retrieved in the sequenced collections ( Table 2). The most extreme example was miR157d, for which we recovered the same numbers of miRNA and miRNA* sequences.

Identification of putative melon-specific miRNAs

After identification of known miRNA sequences and other sRNA sequences (see below), 108,454 unique melon sRNAs remained unclassified, from which the most abundant (28.6%) were 24-nt species. Initial analysis confirmed that 36,783 (33.9%) of these sequences had a perfect match in the melon genome. The frequency distribution was highly skewed: 33,621 sequences had fewer than 25 hits (24,488 originated from a single locus), and only 659 sequences had more than 100 hits.

Sequences that were 21, 22 or 24 nt in length with a maximum of six hits in the genome were selected as potential novel miRNAs, and flanking genomic regions were analysed according to three consecutive criteria. First, we used miRanda software to detect sequences complementary to the potential miRNA inside the flanking regions. Second, potential miRNAs with precursors less than 70 nt in length were discarded. Finally, the MFEI index, which is used to distinguish miRNA precursors from other coding and non-coding RNAs and is based on free energy estimates and nucleotide composition [35], was calculated for each precursor and the results were sorted accordingly (the more negative the index, the better the precursor).

Predicted miRNA precursors and their genomic flanking regions that were found to be similar in sequence to previously described transposons were discarded. Other predicted miRNA precursors with intramolecular folding potential showed no similarity to known transposon sequences although their secondary structures were similar to those of known foldback transposons; these were characterised by strong negative MFEI indexes and high miRanda scores, both features consequence of high sequence complementarity in the pairing stem sequence. For some of these precursors, several uncharacterised melon sRNAs mapped on them in both the sense and antisense orientations (e.g. a11_62726 in Figure 4), up to 85 in some cases. Therefore, these were also considered unsuitable miRNA candidates. Three other potential miR-NAs were shown to be the miRNA* sequences of known miRNAs that had not been picked up in our initial screen.

After manually inspecting the remaining secondary structures, 77 loci that fulfilled the structural criteria for annotation of plant miRNAs [36,37] were selected as plausible miRNA precursors; we also added to this list 7 other loci that had an asymmetric bulge involving 3 bases inside the putative miRNA duplex. From them, 43, 20 and 21 corresponded to sequenced sRNAs of 21, 22 and 24 nt in length, respectively (Table 3). Six selected sequences are shown as examples in Figure 4. By checking the pairing sequence on the stem of the predicted precursors, miRNA*s for seven candidate miR-NAs were found in the sequenced set. Therefore these miRNAs were regarded as authentic miRNAs that conformed to the biogenesis and expression criteria for confident miRNA annotation [37]. The remaining sequences, not supported by the complementary passenger strands, were classified as candidate miRNAs. Most of the potential novel miRNA were represented by a small number of sequences, a single sequence in more than half of the cases, but six exceptional candidates were represented by more than 10 sequences (Table 3). As occurred for conserved miRNAs, sequence variants were identified for some novel miRNAs (Table 3) which mapped on the genomic melon sequence with slight variations in length and position relative to the most abundant sequence. In the absence of a reference sequence from any database, these variants were counted. Sequencing errors of differential cleavage of potential miRNA precursors possibly explain these length and positional polymorphisms.  Table 3) were selected based on quality criteria and are shown as examples. Red lines represent regions where miRNA/miRNA* duplexes are located. Identifier a11_72726 represents an example of unsuitable miRNA candidate. Prediction of potential miRNA targets in the melon transcriptome

To identify potential targets of miRNAs, we screened melon unigenes in the publicly-available database [7]. Two independent searches were performed using miRanda [38] and TargetFinder [39], and the results were compared. Each program scores potential targets based on sequence complementarity, with high scores better in miRanda, and low scores better in TargetFinder. Both algorithms identified a common set of presumptive targets albeit with different scores, and the few discrepancies involved targets with low confidence scores.

Targets in Arabidopsis defined by miRanda generally have a score ≥170, and using this value as cutoff we found 150 melon unigenes as potential miRNA targets, the best of which are summarized in Table 4 (a complete list is provided in Additional file 2). The potential miRNA targets in Table 4 generally had similar annotations to their Arabidopsis counterparts, although there are some exceptions. For example, melon unigene cHS_39-F10-M13R_c is a predicted target of miR159a but it is annotated as positive regulator of brassinosteroid signaling rather than a MYB or TCP transcription factor, which is sensitive to miR159 regulation in a Precursor secondary structures of sRNAs in bold are represented in Figure 4. b miRanda score calculated for the identification of the complementary region to the putative miRNA. c MFEI index calculated as described [35]. d Include cases where sequence variants with up to two mismatches were identified and counted e Indicates miRNAs for which a miRNA* sequence was identified. f Indicates miRNAs for which an asymmetric bulge of more than 2 bases was identified in the miRNA/miRNA* duplex of the precursor, not meeting exactly the structural criteria previously set forth [36,37]. Arabidopsis. Many of the melon unigenes identified as potential targets were not annotated, and some had previously been identified as potential miRNA precursors [5]. Interestingly, the highest miRanda scores (> 300) were achieved for transcripts with two separate miRNA targets on the same molecule. For example, unigene c15d_05-D02-M13R_c had two target sites for miR390ab separated by~200 nt. When this region was used as a BLAST query against the melon sRNA dataset, a group of 257 sequences (more than 92% of them being 21 nt long) was identified with nearby clusters of related 21-nt sequences in both the sense and antisense orientations, which is reminiscent of the ta-siRNAs biogenesis mechanism [18] (Figure 5). Both sites (complementary to miR390 family members in unigene c15d_05-D02-M13R_c) had similar miRanda scores, they did not contain mismatches or G:U wobbles involving nucleotides 9-11 and were phased 21 nt one of each other. The number of sRNA copies was different in each cluster and were more abundant in sense orientation compared to antisense orientation. Two registers of phased 21-nt siRNAs were observed. One of them was phased with the miR390 complementary sites but the other one was not. A representative sequence from each cluster was selected and used to search for potential targets in melon unigenes, identifying > 100 transcripts with a miRanda score > 170. Several of these transcripts were annotated as ARFs and ubiquitin related gene products ( Table 5).

The remaining unigenes with two predicted miRNA target sites listed in Table 4 were annotated as protein-coding transcripts and no sRNAs were identified with similarity to the region flanked by the two target sites ( Figure 5), suggesting that they did not account for authentic ta-siRNA-producing loci. Targets were also sought in the reverse-complement sequences of melon unigenes, because a small proportion of the ESTs could be incorrectly oriented as an artifact of the cloning procedure [5]. Twenty-eight unigenes were identified as potential miRNA targets using the same criteria described above, most of which were found to be nonannotated (Table 6). In this new set of data, unigenes with two targets were used again as a BLAST query against the melon sRNA dataset but no hits were obtained, so these unigenes were no longer considered as potential ta-siRNAs.

With some exceptions, several miRNA targets with miRanda scores ≥170 (see Additional File 3) were identified for each of the potential novel melon-specific miR-NAs listed in the previous section.

Characterization of other melon sRNAs

Next, we blasted our sRNA sequences against RNA and genomic databases to search for other sRNA species by sequence similarity (Figure 6). sRNAs similar to transfer RNA (tRNA), trans-acting siRNA, small nucleolar RNA (snoRNA) and transposons were the least abundant, whereas ribosomal RNAs (rRNA) were largely the most abundant non-coding sRNA species ( Figure 6A). Intriguingly, exogenous virus-derived sRNAs were as abundant as other endogenous plant sRNAs, at least in the case of MNSV. Most of the sRNAs identified had complete sequence similarity with the reference RNA from each    Figure 6A), even if up to two mismatches were allowed in BLAST comparisons. The only exception were sequences similar to ta-siRNAs, for which 14 melon sRNAs with similarity to Arabidopsis TAS3a|D7 (+) and TAS3a|D8(+) were identified, 2 containing 1 mismatch, and 12 containing 2 mismatches. All of them mapped very close in the melon genome and in a different region than c15d_05-D02-M13R_c unigene (the other potential source of ta-siRNAs, see above). To determine if they were authentic melon ta-siRNAs, we selected a 600 bp window sequence upstream and downstream from the genomic location determined in the melon genome for each candidate; then, a BLAST query against the melon sRNA dataset was performed, revealing that at least 126 sequences (95 of them being 21-nt in length) mapped in this region and were arranged according to a near 21-nt phase spacing (data not shown). Many sRNA sequences also generated hits in the plastid genomes (30,239 sRNAs corresponding to 4,254 unique plastid sequences). When these sRNAs were mapped onto the melon chloroplast genome (unpublished data) ( Figure 6B), two clusters of sequences resolved in regions presumably annotated as chloroplast rRNA. These regions lie within two inverted genomic repeats, and sRNAs were accordingly Table 6 Best quality miRNA targets identified in reverse-complement sequences of melon unigenes Some of the chloroplast sRNAs had previously been cloned in other species [40]. For example, melon sRNA a33_398374 (sequence AGT TAC TAA TTC ATG ATC TGG C) was the most abundant melon plastid sRNA (18,054 counts), and a matching sequence is present in more than 900 chloroplast genomes. It is located in an intergenic region and may target a methyltrasferase transcript, although there is no direct evidence that it has silencing functions. Melon sRNA a14_392967_ (sequence GGT AGT TCG ATC GTG GAA TTT) was less abundant (166 counts), it is present in 10 different chloroplast genomes, and it may target a transcript encoding an electron carrier protein. Interestingly, different numbers of plastid sequences were obtained from each library ( Figure 6C). For example, in the virus-resistant melon accession TGR-1551 there was no difference in the number of sRNAs with hits to melon chloroplast genome between healthy and virus inoculated samples, but in the virussusceptible accession Tendral, more sRNAs were counted in inoculated samples ( Figure 6C). Unlike chloroplast sRNAs, only 7,854 sRNAs (corresponding to 2,384 unique sequences) matched the melon mitochondrial genome (unpublished data). These sRNAs were mapped on the mitochondrial genome sequence and formed three clusters, again corresponding to the sites of rRNA genes (data not shown).

Discussion

In this report, we describe the first screen for melon sRNAs by deep sequencing. In total, 398,450 high-quality sequences were generated, representing 90% of the total raw reads. RNA species 21, 24, 20 and 22 nt in length dominated the sRNA transcriptome in melon with the 21-nt class being the most abundant in our libraries. Molecules of 24-nt processed by DCL3 are often the most abundant endogenous plant sRNAs [13], but this may vary according to species. For example, 24nt sRNAs are more abundant in Arabidopsis, rice and tomato [41,42,9], whereas 21-nt sRNAs are more abundant in grapevine, wheat and conifers [12,43,44]. It is also possible that the composition of the sRNA population of a given plant species varies according to tissue and physiological conditions, as seems to be the case of melon (see Additional file I). Perhaps the higher proportion of 24-nt sRNAs found in melon ovaries compared to the other tissues reflects the predominance of developmental processes based on epigenetic events in the ovary.

Recent studies have shown that cis-acting siRNAs arising from heterochromatin, transposons and other repeat elements account for the greatest proportion of endogenous sRNA populations in plants [13,[45][46][47][48]. In melon, only~7,000 sRNAs matched known transposon sequences, in contrast to~60,000 sRNA sequences matching ribosomal RNA, which may simply reflect the paucity of melon transposon sequence information in databases, as only 1.5% of the melon genome has been annotated for transposable elements [49]. Transposon sequences in different species show more divergence than rRNA sequences, so the representation of transposon-related sRNAs could increase when a more accurate and complete annotation of the melon genome becomes available. We also identified two sets of ta-siRNAs in our data, which mapped to different loci in the melon genome thus revealing the presence of at least two potential TAS genes. One locus was not represented in the melon unigene database, most likely because of its incomplete coverage. The sequence of the other locus was similar to that of a non-annotated melon transcript, and contained two registers of sRNAs in a 21-nt phase bounded by two target sites to miR390ab, reminiscent to TAS3 genes. Non-coding transcripts containing two miR390 complementary sites that give rise to phased siRNAs have been described in other organisms. In the moss Physcomitrella patens, both 3' and 5' target sites are cleaved. In Arabidopsis, the 5' miR390 complementary site contains a mismatch and two G:U wobbles involving positions 9-11 and, despite it is not cleaved, it binds the silencing complex and is required for full AtTAS3 function in vivo [50]. In melon, both 3' and 5' miR390 had perfect complementarity at positions 9-11, suggesting that both could be cleaved, as opposed to Arabidopsis, to specify a phased register for ta-siRNA biogenesis. Interestingly, an additional siRNA register that is likely independent of miR390-directed cleavage of the putative melon TAS transcript was observed. This alternative register might be determined by the processing activity of TAS transcripts by one of the most abundant melon primary ta-siRNAs during generation of secondary ta-siRNAs ( Figure 5), as proposed for alternatively phased TAS3 ta-siRNAs in Arabidopsis [18,50]. Since there are additional TAS loci in other plant genomes, it is reasonable that other melon TAS loci remain to be discovered.

More than 30,000 of our sRNA sequences matched the plastid genome, suggesting intense sRNA activity in this organelle. Mitochondrion-specific sRNAs were less abundant in comparison. The abundance of plastid sRNAs varied by source, with fewer sequences obtained from the ovary and fruit libraries compared to the cotyledon libraries, perhaps reflecting a relationship between chloroplast sRNA activity and photosynthesis. Interestingly, there was no significant difference in sRNA accumulation when comparing infected and healthy TGR-1551 cotyledons (resistant to WMV) whereas more sRNA accumulated in healthy Tendral (susceptible to WMV and MNSV) cotyledons than in infected ones. Whether or not this is related with the resistance phenotype is a matter of speculation.

More than 28,000 melon sRNAs in our sequenced collections matched known miRNAs in other plants, and 46 distinct melon sRNA species could be assigned to 26 known miRNA families. Although we generated a relatively low number of sequence reads, our data nevertheless were in good harmony with previous studies of miRNA profiling based on exhaustive sequencing of sRNA populations (e.g. in grapevine, 24 million reads, 26 known miRNA families and 26 non-conserved miRNA families; in tomato, 721,874 reads, 30 known miRNA families; and in orange, 13,106,573 reads, 42 highly-conserved miRNA families) [12,11,9]. This probably reflects the generally-accepted high level of expression reported for conserved miRNAs.

In addition to known miRNAs, 84 sRNA sequences derived from genomic loci with intramolecular folding capacities and not previously described as miRNAs in other plant species were predicted as potential melonspecific miRNAs. In most cases, only one sequence was counted from each of these miRNAs, which is consistent with reports suggesting that species-specific miRNAs are usually expressed at low level and in a tissue-specific manner [41]. The candidates listed in Table 2 include a number of special cases, i.e. miRNAs with miRanda scores ≥195 and very strong secondary structures including an internal loop, resembling type III foldback transposons [51,52,9]. Although these sequences do not match known melon transposons, they were not considered as miRNA candidates because accurate homologybased transposon annotation and prediction occasionally needs to be complemented with ab initio approaches based on structural features [53,49]. However, even not considering this particular group, our data indicate that most of the precursors we identified are candidates to encode melon-specific miRNAs.

The accumulation of miRNAs was estimated by census sequencing and this showed that there is more miRNA diversity and that miRNAs are more abundant in ovaries than fruits. Although miRNAs are involved in many processes, 60-70% of known plant miRNAs control the expression of transcription factors that regulate critical developmental processes, such as proper specification of floral organ identity or leaf polarity, and overexpression or knockout of MIRNA genes led to severe developmental defects [48,54,13]. It is likely that the greater abundance of miRNAs in the early ovary stages compared to fruit reflects the more significant developmental activity in ovaries, and confirms that meristems and other developmentally active tissues are good resources for miRNA screening.

The comparison of healthy and virus-infected melon tissues showed that generally miRNAs were less abundant in infected tissues. Viruses interfere with and exploit endogenous RNA-silencing pathways using diverse strategies [55,19]. For example, the potyvirus silencing suppressor HC-Pro has been shown to suppress the miRNA pathway by inhibiting miRNA assembly into AGO1-containing silencing complexes and unwinding of miRNA/miRNA* duplexes, causing accumulation of stable duplexes [56]. Several studies have shown that virus infection can regulate the accumulation of mRNAs targeted by miRNAs without affecting the abundance of the miRNAs themselves, or even by promoting a slight accumulation [57][58][59]. In contrast, we found that miRNA accumulation was generally depressed in infected plants compared to controls subjected to mock inoculations. A notable exception was miR168ab, which was upregulated in the resistant genotype but downregulated in the susceptible one. This miRNA has previously been shown to be involved in controlling the expression of ARGONAUTE1 (AGO1), the catalytic subunit of the RNA-induced silencing complex responsible for slicing of target mRNAs [60]. Recent work has described the enhanced expression of miR168 and AGO1 mRNA in virus-infected plants specifically and independently of other miRNAs [61][62][63]. The contrasting miRNA profiles observed in the TGR-1551 and Tendral varieties suggests that silencing may underly the resistance of TGR-1551 to WMV, although this is a hypothesis that will require further research.

We have identified more than 150 melon unigenes as potential targets for the known and novel miRNA sequences discovered in this investigation. Many animal transcripts are targets for more than one miRNA but this phenomenon is uncommon in plants [64]. Accordingly, most of melon unigenes identified as potential targets featured only a single miRNA site. miRNAs that are conserved across species tend to have conserved targets too, and our data confirm this is the case in melon. However, several unigenes predicted with high confidence as targets for conserved miRNAs had different annotations to the corresponding target genes in Arabidopsis, although these may represent false positives that would fail additional validation. Furthermore, the nonconserved targets of conserved miRNAs can be cleaved at a lower frequency than conserved targets [12]. For these two reasons, the selection of targets for individual validation experiments can be challenging.

An interesting alternative for miRNA target discovery in a genome-scale is the analysis of the small RNA degradome [65], as this avoids the a priori selection of potential targets. High-throughput gene expression profiling techniques such as microarray hybridization can also help to predict miRNA targets because some times a negative correlation between the abundance of miR-NAs and their target mRNAs can be identified [66,67]. We have used microarrays to monitor gene expression profiles in healthy TGR-1551 and Tendral plants and plants infected with WMV. When compared with our miRNA data, we were able to identify two unigenes encoding AGO proteins that were differentially expressed and showed contrasting expression profiles in susceptible and resistant genotypes (Gonzalez-Ibeas and Aranda, unpublished data). The same profile was observed for miR168 accumulation, suggesting that miR168 may be involved in virus resistance and providing the basis for future experiments.

Conclusion

We have analysed and catalogued a collection of melon endogenous sRNA obtained through massive cDNA sequencing and have identified known miRNAs and ta-siRNAs (conserved in other species) as well as potential melon-specific miRNAs with no database matches. We have also identified potential targets for these miRNAs in the melon transcriptome. Census sequencing (i.e. counting the number of sequence reads for each sRNA) was used to profile their expression in different tissues, and in healthy vs. virus-infected cotyledons. By comparing the predicted targets and the differential expression profiles we were able to provide insights into the role of miRNAs in the regulation of fruit development and plant-virus interactions.

Methods

Plant material

Small RNA libraries were prepared using material from three melon accessions: 1) the Tendral cultivar (Semillas Fitó, Barcelona, Spain), which is susceptible to MNSV and WMV, 2) the breeding line T-111 of the cultivar Piel de Sapo (Semillas Fitó, Barcelona, Spain), and 3) the genotype TGR-1551 (germplasm collection of "Estación Experimental La Mayora" (EELM-CSIC), Málaga, Spain), which is resistant to WMV.

Melon plants were grown under greenhouse conditions (~25/20°C, 16-h photoperiod,~70% relative humidity) in 0.5-L pots with substrate (Tendral and TGR-1551) or in soil bags with the capacity for four plants (Piel de Sapo). Fruits of 15 and 45 days after pollination (DAP) were collected and mesocarp tissues were recovered and used for RNA extractions. Virus infected samples were obtained from completely expanded cotyledons rubbed with carborundum (ø = 0.037 mm) and the corresponding viral inoculum. MNSV-infected melon cotyledons exhibiting lesions and marrow leafs systemically infected with WMV were ground in cold inoculation buffer (0.2 M phosphate buffer pH = 8.0, 0.1% (v/v) beta-mercaptoethanol, 0.03 g/ ml activated charcoal) for inoculum preparation. Mockinoculated control cotyledons were rubbed with inoculation buffer and carborundum alone.

Cotyledons were harvested at 1, 3, 5 and 7 days postinoculation (dpi) and pooled for RNA extraction. Fruit samples were prepared as previously described [8].

Ovaries were collected at stages C1 and C5 (Mascarell-Creus et al., unpublished). The C1 stage corresponds to flower emergence from the inflorescence bud, when the outermost perianth organs commence development and no floral whorls are visible. The C5 stage corresponds to anthesis, when the flower is ready to be fertilized and all floral organs are fully formed, including the yellow petals that attract pollinators. Under normal growth conditions, C1 to C5 development takes approximately 5 days.

Small RNA library construction

DNA amplicons were gel-purified using 4% Metaphor Agarose and isolated using the QIAEX II Gel Extraction Kit (Qiagen, Hilden, Germany). The quantity and quality of the DNA amplicons were determined using a ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, Delaware, USA) and an Experion Automated Electrophoresis System (Bio-Rad, Hercules, California, USA). The same quantity of DNA from each library was pooled and sequenced using the 454 Life Science Technology platform (Lifesequencing S.L., Paterna, Valencia, Spain). Sequence data in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE28653 http://www.ncbi.nlm.nih.gov/geo/query/acc. cgi?acc=GSE28653.

Bioinformatics

The sRNA sequences were parsed from FASTA-formatted files containing 447,180 reads from two independent 454 sequencing runs and assigned to specific libraries by identifying the sRNA/adaptor boundaries and barcode analysis. Sequences were analysed with standard Python scripts [68] and the BioPython library [69]. Only sequences with the 3' and 5' adaptors in the correct position were considered. Known sRNAs were identified by searching public databases using BLAST version 2.2.19 [70] and allowing up to two mismatches. The following databases and sequences were searched: Transfer RNA Database (version 2009) [71], Plant Small Nucleolar RNA Database (v1.2) [72], SILVA (ribosomal RNA database, v100) [73], The Arabidopsis Small RNA Project (ASRP) Database [33], Rfam Database 10.0 [74], miRBase (release 16) [34], The Plant Repeat Database [75], Cucumis melo chloroplast genome (unpublished data), Cucumis melo mitochondrial genome (unpublished data), MNSV genome (GenBank accession

Figure 1
1Representation of sequences with different lengths in the melon sRNA data set. (a) All sequences. (b) Non-redundant sequences. The number of sequences is expressed as a percentage of the total number of sequences.

Figure 2A, B compares the accumulation of miRNAs in healthy versus WMV-infected melon tissues from

Figure 2
2Relative accumulation of conserved miRNAs in melon samples used for sRNA library construction. Total reads for each miRNA in each library were normalised relative to the total number of reads from the library, and expressed per 10,000 reads.(a) Cotyledons from melon cv. Tendral inoculated with WMV-M116 compared to mock inoculated cotyledons of the same cultivar. (b) Cotyledons from the melon accession TGR-1551 inoculated with WMV-M116 compared to mock inoculated cotyledons of the same accession. (c) Stage C1 and C5 ovaries from melon cv. Piel de Sapo. (d) Fruit from melon cv. Piel de Sapo 15 days after pollination (15d) compared to fruit from the same cultivar 45 days after pollination (45d). (e) Cotyledons from melon cv. Tendral inoculated with MNSV-alfa5 compared to mock-inoculated cotyledons of the same cultivar. (f) Cotyledons from melon cv. Tendral inoculated with MNSV (chimeric virus) compared to mock inoculated cotyledons of the same cultivar.

Figure 3
3Duplexes of mature miRNA and passenger (miRNA*) sequences identified in the melon sRNA collections. (a) Typical duplex structure. (b) Non-typical duplex structure (number of protruding nucleotides ≠ 2).

Figure 4
4Secondary structures of putative novel melon-specific miRNA precursors. Some of the sRNAs identified as potentially novel and melon-specific miRNAs (listed in

Figure 5
5Identification of potential ta-siRNA transcripts. Four unigenes were found to have two miRNA target sites each. For each unigene, the region delimited by the two miRNA suites was used as a BLAST query against the sRNA data set and only c15d_05-D02-M13R_c generated hits. Unigenes are represented by black horizontal lines. The miRNA site boundaries are represented by vertical arrows. Potential ta-siRNAs are represented in the inset by colored short horizontal lines mapped onto the unigene.

both the sense and antisense orientations.

Figure 6
6Types and frequencies of known melon sRNAs. (a) Number of sRNAs with similarity to known sequences in public databases: tRNA = transfer RNA; snoRNA = small nucleolar RNA; LSUrRNA = large subunit of ribosomal RNA; SSUrRNA = small subunit of ribosomal RNA; taRNA = trans-acting-si-RNAs; transposon = transposon sequences; chloroplast = melon chloroplast genome; mitochondrion = melon mitochondrial genome; MNSV = Melon necrotic spot virus genome; WMV = Watermelon mosaic virus genome; miRNA = microRNA. (b) Mapping of small RNAs onto the melon chloroplast genome. X-axis represents genome nucleotide positions. Y-axis represents the number of sRNAs mapped at each position. The chloroplast genome is represented by a horizontal black line (156,018 nt in length) at y = 0, and the 5' ends of sRNAs are represented by points (blue: sense sRNAs, pink: antisense sRNAs). The frequencies of antisense sRNAs are shown as negative numbers. (c) Number of melon sRNAs with similarity to the melon chloroplast genome.

Table 1
1Description of small RNA libraries from different melon tissuesLibrary
Cultivar/accession
Tissue
Physiological condition
Virus a
Reads
Unique sequences

Wtm
cv. Tendral
Cotyledon
Mock-inoculated
-
33123
15624

Wt
cv. Tendral
Cotyledon
Virus-infected
WMV-M116
35860
12840

Cwm
accession TGR-1551
Cotyledon
Mock-inoculated
-
41039
21122

Cw
accession TGR-1551
Cotyledon
Virus-infected
WMV-M116
36330
24100

15d
cv. Piel de Sapo
Fruit
Healthy, 15 days after pollination
-
21662
14620

45d
cv. Piel de Sapo
Fruit
Healthy, 45 days after pollination
-
9942
8167

c1
cv. Piel de Sapo
Ovary
Healthy
-
18764
15269

c5
cv. Piel de Sapo
Ovary
Healthy
-
14529
12608
Ta5
cv. Tendral
Cotyledon
Virus-infected
MNSV-alfa5
43170
22869

3'T
cv. Tendral
Cotyledon
Virus-infected
MNSV (chimeric)
56425
56425

a

WMV = Watermelon mosaic virus; MNSV (alfa5) = Melon necrotic spot virus, alfa5 isolate; MNSV (chimeric) = Melon necrotic spot virus, alfa5 isolate with 3' UTR
from 264 isolate

Table 2
2Known plant miRNAs identified in melonAnnotation
Melon sRNA sequence (5'-3') Similarity
Number of miRNA
sequences

miRNA*
sequences

Hit in melon
genome a

miR156|a, b, c, d, e, f
UGACAGAAGAGAGUGAGCAC
100%
469
60
YES

miR157|a, b, c
UUGACAGAAGAUAGAGAGCAC
100%
269
0
YES
miR157|d
UGACAGAAGAUAGAGAGCAC
100%
19
21
YES

miR158|a
UCCCAAAUGUAGACAAAGCA
100%
1
0
-

miR159|a
UUUGGAUUGAAGGGAGCUCUA
100%
14651
0
YES

miR159|b
UUUGGAUUGAAGGGAGCUCUU
100%
18
0
-

miR159|c
UUUGGAUUGAAGGGAGCUCCU
100%
1
0
-

miR160|a, b, c
UGCCUGGCUCCCUGUAUGCCA
100%
537
0
YES

miR161|a.1
UUGAAAGUGACUACAUCGGGG
100%
6
0
-

miR161|a.2
UCAAUGCAUUGAAAGUGACUA
100%
1
0
-

miR162|a, b
UCGAUAAACCUCUGCAUCCAG
100%
825
0
YES

miR164|a, b
UGGAGAAGCAGGGCACGUGCA
100%
172
1
YES

miR165|a, b
UCGGACCAGGCUUCAUCCCCC
100%
4
0
-

miR166|a, b, c, d, e, f,
g

UCGGACCAGGCUUCAUUCCCC
100%
65
27
YES

miR167|a, b
UGAAGCUGCCAGCAUGAUCUA
100%
136
0
YES

miR167|d
UGAAGCUGCCAGCAUGAUCUGG
100%
16
1
-

miR168|a, b
UCGCUUGGUGCAGGUCGGGAA
100%
967
0
YES

miR169|a
CAGCCAAGGAUGACUUGCCGA
100%
3
1
-

miR169|b, c
CAGCCAAGGAUGACUUGCCGG
100%
76
1
YES
miR169|h, i, j, k, l, m, n UAGCCAAGGAUGACUUGCCUG
100%
83
1
YES

miR170|a
UGAUUGAGCCGUGUCAAUAUC
100%
3
0
-

miR171|a
UGAUUGAGCCGCGCCAAUAUC
100%
85
5
YES

miR171|b, c
UUGAGCCGUGCCAAUAUCACG
100%
64
0
YES

miR172|a
AGAAUCUUGAUGAUGCUGCAU
100%
85
58
YES

miR172|c, d
AGAAUCUUGAUGAUGCUGCAG
100%
4
0
YES

miR172|e
GGAAUCUUGAUGAUGCUGCAU
100%
3
0
YES

miR319|a, b
UUGGACUGAAGGGAGCUCCC
100%
2
3
YES
miR390|a, b
AAGCUCAGGAGGGAUAGCGCC
100%
32
6
YES

miR391|a
UUCGCAGGAGAGAUAGCGCCA
100%
1
0
-

miR393|a, b
UCCAAAGGGAUCGCAUUGAUC
100%
18
0
YES

miR394|a, b
UUGGCAUUCUGUCCACCUCC
100%
4
0
YES

miR396|a
UUCCACAGCUUUCUUGAACUG
100%
134
84
YES

miR396|b
UUCCACAGCUUUCUUGAACUU
100%
82
16
YES

miR397|a
UCAUUGAGUGCAGCGUUGAUG
100%
26
0
YES

miR408|a
AUGCACUGCCUCUUCCCUGGC
100%
14
1
YES
ath-miR2111a
UAAUCUGCAUCCUGAGGUUUA
100%
1
0
YES

peu-miR2910
UAGUUGGUGGAGCGAUUUGUC
100%
8
0
YES

osa-miR167d
UGAAGCUGCCAGCAUGAUCUG
100%
3401
1
YES

tae-miR395b
UGAAGUGUUUGGGGGAACUC
100%
1
0
YES

bna-miR397a
CAUUGAGUGCAGCGUUGAUGU
95%
77
0
YES

miR156|h
UUGACAGAAGAGAGUGAGCAC
95%
91
0
YES

miR156|g
ACAGAAGAGAGUGAGCACA
90%
5
0
YES

miR169|d, e, f, g
UGAGCCAAGGAUGACUUGCCU
95%
130
0
YES
miR169|d, e, f, g
UGAGCCAAAGAUGACUUGCCU
90%
112
0
YES

miR399|a
UGCCAAAAGAGACUUGCCCUG
95%
3
0
YES

miR403|a
CUAGAUUCACGCACAAGCUCG
90%
1
0
-

Table 3
3Potential novel melon specific miRNAsMelon sRNA
name a

nt
sRNA sequence
Number of
sequences

Hits in
genome

miRanda
score b

Potential precursor
lenght (nt)

MFEI c

a34_130677_
21
AUAGAUAUUGAUAUGCUUUUA
1
4
163
94
-1.0573

a24_2602_
21
UGCUACAUGGUUUAUCAGUGA
2
5
115
72
-1.2524
a24_177791_
21
UCGCAGAAGAGAUGGCGCCGA
7
1
143
91
-0.8587

a23_118111_
21
CAUUGAUAGACACUAAUAGAA
1
5
167
90
-1.475

a33_14294_
21
AUAGACUUCUAUUGGUGUCUA
1
1
154
74
-1.5105

a32_31625_
21
GUUCCCACGGUAAUGAUAAUA
2
1
127
72
-1.1045

a32_1324_
21
AGGUGUCAUCUUGCUGCGAUA
1
1
179
99
-1.2

a22_190223_
21
UGUUAUGCAUGGCGUCGGGAG
1
5
187
157
-1.2915

a14_98657_
21
AUAGCGAAGUAUAUCAGUGAU
1
1
154
127
-1.2423

a14_51701_
21
UGAGCCGUGCCAAUAUCGACG
1
1
159
97
-1.1
a14_180374_
21
UAAAUAUUUAGAAAGUCAAUC
1
4
159
79
-1.855

a21_80766_ f
21
UAAUAUUCAUUUUCACUUCUU
1
1
116
232
-1.2494

a21_78863_
21
UGCAUCCUGAGGUUUAGGGAG
3
1
159
194
-0.9776

a21_244426_
21
AGUAACCACUAAGCUAAUGGC
1
1
187
597
-1.0096

a23_1441_ f
21
UAUAGCAAAGUCUAUCGAUGG
5 d
1
107
77
-1.1435

a13_84447_ e
21
GGUCAUUCUAGCAGCUUCAAU
13 d
1
139
201
-0.96

a23_370_ e
21
UGGUGUGCAUGUGAUGGAAUA
13 d
1
163
113
-0.9061

a21_134553_ f 21
GUUAUACGAGUUGGGUUGGGU
1
1
155
114
-0.9571
a11_95657_ e 21
UUGUGUCAUUGACAUUGUGGU
1
2
195
191
-1.1708

a14_9427_
21
UCGUCCUGAGAAUACAUGUCA
35 d
1
159
97
-0.9409

a14_668_
21
UGAGUUAUCGGUGAAUUCAAG
5 d
3
159
516
-0.9675

a13_252112_
21
ACUGCUGCUUGUACUAUUGAA
1
1
191
255
-1.9670

a11_33177_
21
UUUAGUUUAGCCUAUUGCUUU
1
3
187
139
-1.1035

a11_191362_ e 21
UUCUAUUGUCUUCAUUUGUGA
1
1
191
119
-1.2718

a34_224062_
21
UGAAAUGACUUGUCAAGUGCU
1
1
151
97
-0.975

a13_228150_
21
CUUGUACUUGAUUUUGUUGCC
1
1
191
115
-1.4469
a12_32299_ e
21
AAUUUGUUGGUCAAAUGAUUG
2
1
195
107
-1.7552

a12_272161_
21
UUGUAUGGUGGAAAGAUGGAA
1
1
162
96
-1.4167

a11_33986_
21
GCUGACUUGCUGAUUGAGUUA
2
3
179
189
-1.4852

a13_33760_
21
UGAAUUAUCUGCUUAAGUUUU
1
1
187
95
-1.2889

a11_389198_
21
ACACGCAGAAGAGACGAUUGA
1
1
191
120
-1.5575

a13_357842_
21
UGGAGCAAUAUUGAUGCAUAU
1
1
195
220
-0.9548

a11_364692_
21
UUGGGUCUAUUUAAUGGGAGC
1
1
155
107
-1.2308

a13_281334_
21
ACUUUCUGUCAAUAUAAUCAG
1
1
175
115
-1.3943
a12_71107_
21
UAUCAUAGUUGGUGGUUCAGG
3
1
143
115
-1.2167

a13_120551_
21
UCAACGAUAGACAUUGAUAGA
1
1
171
107
-1.2875

a12_123886_
21
UUAUCAUUGAUAGACUAGUAU
1
2
155
174
-1.0612

a11_227522_
21
CAAGCCCAUGACAAAGCAAGC
1
1
187
225
-1.0929

a14_133932_ f 21
UCAACACGAUCGUCUAGCAUG
1
2
173
113
-1.2295

a11_203340_
21
UUUGAGUGUCCUACUCACCUC
1
1
191
411
-1.0833

a11_11135_
21
UAGUGCCGCGCUGCGUGCGUC
85
1
147
102
-0.98

a11_31022_
21
UUUCGCUUUUCCUCUUUCGUG
1
1
191
454
-1.1523
a12_144938_
21
UCGUGGAUAUUGCUCUUUUCU
2
1
171
504
-1.3303

a33_181157_
22 GAUAGAUACUAAUAUGCUUCUA
1
2
188
87
-1.3227

a33_37151_
22 AGAUUAAUUUAUUGGGCGUUAU
1
2
144
94
-1.3313

a12_72169_
22 UUGAGCUAUGCUCAGGUUGACA
30
1
176
174
-1.2933

a24_96796_ e
22 UGAGCUAUGCUCGCUUUGGCAA
21 d
1
175
169
-1.2855

a11_378153_
22 GAGUUCCUAAGUUUUGAUGAAU
1
1
144
351
-1.2221

a11_378297_
22 UUUUGGAUUCUAUCGAUGAAAG
1
1
155
123
-1.4857

a23_244052_ e 22 GGGCAGCCCCACGUUGGGCAUG
5 d
1
175
353
-0.9263
a11_85662_
22 AAAUAUAUCGGUGUCUAUCAAU
1
2
132
85
-1.2208

Table 3
3Potential novel melon specific miRNAs (Continued)a21_388555_
22 GAUAGACGCUGAUAGAUAGACA
3
1
124
76
-0.8963

a11_84237_
22
CGGCCAAAAAUGACUUGCCCGG
2
1
150
105
-0.8902

a23_124460_
22 AGGUGAGUUCUUUUUAUAGGCU
1
1
184
179
-1.6306

a23_163065_
22 AUUUGAUUAGCCAAAUUUAAAC
2
2
148
130
-1.0514
a23_71826_
22 CAAUAGUCAGAUGUAAACGAUC
1
1
180
228
-1.4282

a22_52587_
22 AAAAUUAUUGGGUGAAUUAGUU
2
1
179
162
-0.9013

a13_234225_
22 UGAAUUUUGUUAUGUUUUGUAA
1
2
168
80
-2.0417

a13_286453_
22
AGUCUAUCACCGAUAGAAGCCU
1
6
182
430
-1.4372

a21_339397_
22 CGCGAGGUUCUUUGUUUGUCUU
1
1
192
413
-1.3341

a14_283014_
22 UACCUAGUGAUGCCAUUGUCAA
1
1
184
533
-1.2829

a21_170878_
22 CUAAGGUUGCCCAGAGAUGUUC
1
1
163
271
-1.2341

a21_63125_ f
22 GAAUAAUUAUCAAGUGUGUAGC
1
1
165
210
-1.7804
a11_146182_
24 UUAAAAUGUUGCUAUAUAAUUAAU
1
6
202
390
-1.6207

a21_169735_
24 UAUACGGGCCGUAAAUAGUUUGAU
2
6
132
369
-0.8885

a23_3672_
24 UUGCUCAUUGCUAACUGCAAAGAG
1
3
198
183
-1.7594

a14_260703_ f 24 UUAAAAGUAUGAGACGAAAAUGAA
1
1
190
111
-1.4204

a14_218837_
24 UUAAAAAGAACUACACGAACGUGC
1
1
194
342
-1.1314

a14_148530_ f 24 AAAUAGCGUCGGGGAAAGGUGUCU
3
1
152
73
-1.0167

a21_127379_
24 UGGGACAAAAGAAAACUGUGGGUC
2
1
162
586
-1.32

a11_2899_
24 AUGAUGCUUUGGUGCUAAGGAGGU
1
1
198
470
-1.6094
a11_248538_
24 AUUUUUGGCAUUUUACACGGUGAG
2
2
192
218
-1.5206

a13_59670_
24 AGUGGAGUGGGCUAUUUUAGUCCA
6
1
166
570
-1.1249

a32_72333_
24 AACUAUUUUAUUGGAACAUGUUGA
3
1
170
96
-1.0391

a33_22103_
24 AGUAUGAUCUCGGGCUAAGGUUGC
1
4
202
246
-1.4138

a24_224684_
24 AACAAAACGAAUGAUCAAAAUGGU
3
1
134
80
-0.8871

a13_225824_
24 ACCAAAUGGAUCUAUUCUUAUAAU
1
1
198
311
-1.3393

a24_82972_
24 AACGAUCGGGUUGACUACGUAAAU
3
1
190
146
-0.9354

a21_98074_
24 AAUUUUUAGUGGUCCCGAAAUGCA
3
1
166
112
-0.9033
a11_350684_
24 UGAGUGUAUCAUCGAGAUAGUGCG
1
1
190
307
-1.1642

a21_216237_
24 AAAUUUCAGGGUCUAAAUUGAUGC
2
1
138
447
-1.2096

a33_49599_
24 CAGCGUGAUUGAUGGGGCAUUUUU
3
1
118
287
-1.4678

a33_58155_
24 AUGGUCGAUCUCAACCGAGAUUGA
1
1
194
298
-1.3188

a23_103110_
24 CGUGAAGAUUGUGGAUAUUGGAGA
1
1
190
414
-1.5507

Table 4
4Best quality miRNA targets identified in melon unigenesmiRNA
annotation

Unigene
Score
(miRanda)

Score
(TargetFinder)

Unigene annotation

miR390|a, b
c15d_05-D02-
M13R_c

362
2.5
non-annotated unigene

miR390|a, b
c15d_05-D02-
M13R_c

362
4
non-annotated unigene

miR390|a, b
c15d_21-G08-
M13R_c

362
2.5
non-annotated unigene

miR390|a, b
c15d_21-G08-
M13R_c

362
4
non-annotated unigene

miR391|a
cCL286Contig1
325
-
histone H1, putative

miR164|a, b
cA_04-D07-
M13R_c

319
-
non-annotated unigene

miR390|a, b
cCL384Contig1
316
-
ATSK11, SK 11ATSK11; protein kinase/protein serine/threonine kinase

miR391|a
cPSI_29-G09-
M13R_c

313
-
UVR8UVR8 (UVB-RESISTANCE 8); chromatin binding/guanyl-nucleotide
exchange factor

miR167|d
cCL1653Contig1
200
-
non-annotated unigene
miR167|d
cCL2516Contig1
200
-
non-annotated unigene

miR167|a, b
cCL1653Contig1
195
1
non-annotated unigene

miR167|a, b
cCL2516Contig1
195
0
non-annotated unigene

miR168|a, b
c46d_19-A03-
M13R_c

195
0
non-annotated unigene

miR171|a
cHS_39-C12-
M13R_c

195
0
scarecrow-like transcription factor 6 (SCL6)

miR397|a
c15d_32-E08-
M13R_c

195
0
potential miR397a precursor

miR160|a, b, c cCL5073Contig1
191
0.5
ARF17ARF17 (AUXIN RESPONSE FACTOR 17); transcription factor

miR164|a, b
cPSI_18-H09-
M13R_c

190
3
ANAC100, ATNAC5ANAC100 (ARABIDOPSIS NAC DOMAIN CONTAINING
PROTEIN 100)

miR393|a, b
cCL3757Contig1
190
2
AFB2AFB2 (AUXIN SIGNALING F-BOX 2); auxin binding/ubiquitin-protein
ligasechr3

miR393|a, b
cCL4853Contig1
190
1
TIR1TIR1 (TRANSPORT INHIBITOR RESPONSE 1); auxin binding/protein binding/
ubiquitin-protein ligase

miR408|a
cCL975Contig1
190
2.5
ARPNARPN (PLANTACYANIN); copper ion binding/electron carrierchr2

miR408|a
cHS_18-D07-
M13R_c

190
2.5
ARPNARPN (PLANTACYANIN); copper ion binding/electron carrierchr2

miR157|a, b, c
c46d_26-C05-
M13R_c

187
3
SPL4SPL4 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 4); DNA binding/
transcription factor

miR157|a, b, c
cCI_30-A09-
M13R_c

187
2
SPL9SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9); transcription
factor

miR157|a, b, c cCL2877Contig1
187
3
SPL3SPL3 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 3); DNA binding/
transcription factor

miR164|a, b
cCI_64-A04-
M13R_c

187
2
NAC1, ANAC022NAC1; transcription factorchr1

miR170|a
cHS_39-C12-
M13R_c

187
1.5
scarecrow-like transcription factor 6 (SCL6)

miR159|a
cHS_39-F10-
M13R_c

183
3
brassinosteroid signaling positive regulator-related

miR161|a.2
cA_16-D06-
M13R_c

183
3
pentatricopeptide (PPR) repeat-containing protein

miR169|a
cA_37-E12-
M13R_c

183
4
NF-YA9NF-YA9 (NUCLEAR FACTOR Y, SUBUNIT A9); specific transcriptional
repressor/transcription factor

miR156|a, b, c,
d, e, f

c46d_26-C05-
M13R_c

182
2
SPL4SPL4 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 4); DNA binding/
transcription factor

miR156|a, b, c,
d, e, f

cCI_30-A09-
M13R_c

182
1
SPL9SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9); transcription
factor
miR156|a, b, c,
d, e, f

cCL2877Contig1
182
2
SPL3SPL3 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 3); DNA binding/
transcription factor

Table 4
4Best quality miRNA targets identified in melon unigenes (Continued)miR157|d
c46d_26-C05-
M13R_c

182
3
SPL4SPL4 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 4); DNA binding/
transcription factor

miR157|d
cCI_30-A09-
M13R_c

182
2
SPL9SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9); transcription
factor

miR157|d
cCL2877Contig1
182
3
SPL3SPL3 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 3); DNA binding/
transcription factor

miR319|a, b
c15d_24-H05-
M13R_c

182
4
potential miR319|a, b precursor

miR167|
d_melon

cCL1653Contig1
-
0
non-annotated unigene

miR167|
d_melon

cCL2516Contig1
-
0.5
non-annotated unigene

miR172|e
c15d_13-C08-
M13R_c

-
2
non-annotated unigene

miR172|e
cA_04-D07-
M13R_c

-
2
non-annotated unigene

miR167|a, b
cCL2288Contig1
-
2.5
unknown protein

miR156|a, b, c,
d, e, f

cCL5542Contig1
173
2.5
kelch repeat-containing F-box family protein

miR157|a, b, c cCL2547Contig1
175
2.5
unknown protein

miR164|a, b
cCL2655Contig1
175
2.5
non-annotated unigene

Table 5
5miRNA targets identified in melon transcripts for potential ta-siRNAs derived from unigene c15d_05-D02-M13R_cMelon sRNA
name

sRNA sequence
Targeted
unigene

Unigene
sense

Score
(miRanda)

Unigene annotation

a11_156739_
AGTTTGCTTCTTGGGCTCTTC
cA_05-B09-
M13R_c

Forward
175
IAA16; transcription factor

a11_156739_
AGTTTGCTTCTTGGGCTCTTC
cPS_07-G03-
M13R_c

Forward
175
IAA16; transcription factor

a14_55988_ AGAGCCCAAGAAGCAAACTGG cCL678Contig1
Forward
172
auxin efflux carrier family protein

a24_92242_
AGAGCCCAAGAAGCAAACTG cCL678Contig1
Forward
172
auxin efflux carrier family protein

a33_151240_
CAGTTTGCTTCTTGGGCTCTT
c15d_39-H01-
M13R_c

Forward
171
ARF6 (AUXIN RESPONSE FACTOR 6); transcription
factor

a14_362833_
CGATGGTGATGGGATTTTTGA cCL1479Contig1
Reverse
171
IAA9 (INDOLE-3-ACETIC ACID INDUCIBLE 9);
transcription factor

a14_362833_
CGATGGTGATGGGATTTTTGA cCL4756Contig1
Reverse
175
ATAUX2-11 (AUXIN INDUCIBLE 2-11); DNA binding/
transcription factor

a14_362833_
CGATGGTGATGGGATTTTTGA
cP5.72_c
Reverse
175
IAA7 (INDOLE-3-ACETIC ACID 7); transcription factor

a33_203464_
CATTTTTTACGATGGTGATGG cCL3310Contig1
Forward
179
ATUBP3 (ARABIDOPSIS THALIANA UBIQUITIN-
SPECIFIC PROTEASE 3)

a14_362833_
CGATGGTGATGGGATTTTTGA cCL3310Contig1
Forward
175
ATUBP3 (ARABIDOPSIS THALIANA UBIQUITIN-
SPECIFIC PROTEASE 3)

a14_20904_
TACGATGGTGATGGGATTTTT cCL4210Contig1
Forward
174
ubiquitin-associated (UBA)/TS-N domain-containing
protein

a11_156739_
AGTTTGCTTCTTGGGCTCTTC cCL1290Contig1
Forward
171
binding/ubiquitin-protein ligase

AcknowledgementsThis work was supported by grants AGL2009-07552/AGR, BIO2006-13107 (Ministerio de Ciencia e Innovación, Spain) and MELONOMICS (Fundación Genoma España, Spain). We thank Mari Carmen Montesinos and Blanca Gosalvez for their technical assistance and Richard M. Twyman email@example.com for editorial assistance.AY122286.1), WMV genome (GenBank accession AY437609.1). In the case of miRNAs, melon sequences were named with the reference miRNA from each database in order to distinguish miRNA species of each family. miRNA targets were identified using miRanda v3.0[38]and TargetFinder Perl script 1.5[39]. Putative novel melon-specific miRNA genes were identified by using the candidate miRNA as a BLAST query against the melon genome (unpublished data). For each hit, 600 bp of sequence upstream and downstream of the alignment was used to search for a near-perfect reverse complement (miRNA*) sequence with the miRanda algorithm. Regions lacking a corresponding miRNA* sequence were discarded. Minimum genomic regions (> 70 nt) containing miRNA and miRNA* sequences were selected as potential precursors. Those corresponding to protein-coding genes were identified by BLAST searches against the Arabidopsis protein database (TAIR) and were discarded, whereas non-coding potential precursors were manually inspected and used to predict the RNA secondary structure with Mfold[76]and for calculation of the MFEI index[35]. Precursors that met structural miRNA criteria were selected for further evaluation[36,37].Additional materialAdditional file 1: Length distribution of the small RNA data set for each library. Length distribution of melon sRNAs for each library (listed inTable 1). Sequence numbers are shown as a percentage of the total number of sequences obtained from every library. Data are given for total (with redundancy) and unique (no redundancy) sequences.Additional file 2: Known miRNA targets identified in melon unigenes. Complete set of all known miRNA targets identified in melon unigenes.Additional file 3: Novel melon-specific miRNA targets identified in melon unigenes. Complete set of all novel melon-specific miRNA targets identified in melon unigenes.Authors' contributions DGI prepared RNA from infected and mock-inoculated samples, constructed all the sRNA libraries, carried out the trimming and analysis of the sRNA sequences and wrote the manuscript. JB provided bioinformatics analysis support. LD and CL provided guidance for the preparation of the sRNA libraries and additional technical support. ACD, AMC, MS and JGM prepared RNA from fruit and ovary samples. MAA supervised DGI including writing of the manuscript, and conceived this study together with DGI, JGM and CL. All authors read and approved the final manuscript.
. Apdo. correos. Departamento de Biología del Estrés y Patología Vegetal, Centro de Edafología y Biología Aplicada del Segura (CEBAS) -CSIC164Departamento de Biología del Estrés y Patología Vegetal, Centro de Edafología y Biología Aplicada del Segura (CEBAS) -CSIC, Apdo. correos 164,

Espinardo, Madrid, Spain. 4 IRTA, Center for Research in Agricultural Genomics CSIC-IRTA-UAB. Murcia), Spain. 2 Departamento de Biotecnología, Instituto de Conservación y Mejora de la Agrodiversidad Valenciana (COMAV) -UPV, Camino de Vera s/n, 46022 Valencia, Spain. 3 Departamento de Biología MedioambientalCampus UAB, Edifici CRAG, Bellaterra (Cerdanyola del Vallès; Barcelona), Spain; Campus UAB, Edifici CRAG, Bellaterra (Cerdanyola del Vallès; Barcelona), Spain8193Centro de Investigaciones Biológicas (CIB) -CSIC. References 1. FAOSTAT.Espinardo (Murcia), Spain. 2 Departamento de Biotecnología, Instituto de Conservación y Mejora de la Agrodiversidad Valenciana (COMAV) -UPV, Camino de Vera s/n, 46022 Valencia, Spain. 3 Departamento de Biología Medioambiental, Centro de Investigaciones Biológicas (CIB) -CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain. 4 IRTA, Center for Research in Agricultural Genomics CSIC-IRTA-UAB, Campus UAB, Edifici CRAG, Bellaterra (Cerdanyola del Vallès), 08193 (Barcelona), Spain. 5 Molecular Genetics Department, Center for Research in Agricultural Genomics (CRAG) CSIC-IRTA-UAB, Campus UAB, Edifici CRAG, Bellaterra (Cerdanyola del Vallès), 08193 (Barcelona), Spain. References 1. FAOSTAT. [http://faostat.fao.org/site/339/default.aspx].

The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. E Huala, A W Dickerman, M Garcia-Hernandez, D Weems, L Reiser, F Lafond, D Hanley, D Kiphart, M Zhuang, W Huang, L A Mueller, D Bhattacharyya, D Bhaya, B W Sobral, W Beavis, D W Meinke, C D Town, C Somerville, S Y Rhee, Nucleic Acids Res. 29Huala E, Dickerman AW, Garcia-Hernandez M, Weems D, Reiser L, LaFond F, Hanley D, Kiphart D, Zhuang M, Huang W, Mueller LA, Bhattacharyya D, Bhaya D, Sobral BW, Beavis W, Meinke DW, Town CD, Somerville C, Rhee SY: The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic Acids Res 2001, 29:102-105.

Characterization of six varieties of Cucumis melo L. based on morphological and physiological characters, including shelf-life of fruit. L Liu, F Kakihara, M Kato, Euphytica. 135Liu L, Kakihara F, Kato M: Characterization of six varieties of Cucumis melo L. based on morphological and physiological characters, including shelf-life of fruit. Euphytica 2004, 135:305-313.

Temporal Sequence of Cell Wall Disassembly in Rapidly Ripening Melon Fruit. J K Rose, K A Hadfield, J M Labavitch, A B Bennett, Plant Physiol. 117Rose JK, Hadfield KA, Labavitch JM, Bennett AB: Temporal Sequence of Cell Wall Disassembly in Rapidly Ripening Melon Fruit. Plant Physiol 1998, 117:345-361.

MELOGEN: an EST database for melon functional genomics. D Gonzalez-Ibeas, J Blanca, C Roig, M Gonzalez-To, B Pico, V Truniger, P Gomez, W Deleu, A Cano-Delgado, P Arus, F Nuez, J Garcia-Mas, P Puigdomenech, M Aranda, BMC Genomics. 8306Gonzalez-Ibeas D, Blanca J, Roig C, Gonzalez-To M, Pico B, Truniger V, Gomez P, Deleu W, Cano-Delgado A, Arus P, Nuez F, Garcia-Mas J, Puigdomenech P, Aranda M: MELOGEN: an EST database for melon functional genomics. BMC Genomics 2007, 8:306.

Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. C Clepet, T Joobeur, Y Zheng, D Jublot, M Huang, V Truniger, A Boualem, M E Hernandez-Gonzalez, R Dolcet-Sanjuan, V Portnoy, A Mascarell-Creus, A I Cano-Delgado, N Katzir, A Bendahmane, J J Giovannoni, M A Aranda, J Garcia-Mas, Z Fei, BMC Genomics. 12252Clepet C, Joobeur T, Zheng Y, Jublot D, Huang M, Truniger V, Boualem A, Hernandez-Gonzalez ME, Dolcet-Sanjuan R, Portnoy V, Mascarell-Creus A, Cano-Delgado AI, Katzir N, Bendahmane A, Giovannoni JJ, Aranda MA, Garcia-Mas J, Fei Z: Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. BMC Genomics 2011, 12:252.

. International Cucurbit Genome Initiative, ICuGI. International Cucurbit Genome Initiative, ICuGI. [http://www.icugi.org].

An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L. A Mascarell-Creus, J Cañizares, J Vilarrasa-Blasi, S Mora-García, J Blanca, D Gonzalez-Ibeas, M Saladié, C Roig, W Deleu, B Picó-Silvent, N López-Bigas, M A Aranda, J Garcia-Mas, F Nuez, P Puigdomènech, A I Caño-Delgado, BMC Genomics. 10467Mascarell-Creus A, Cañizares J, Vilarrasa-Blasi J, Mora-García S, Blanca J, Gonzalez-Ibeas D, Saladié M, Roig C, Deleu W, Picó-Silvent B, López-Bigas N, Aranda MA, Garcia-Mas J, Nuez F, Puigdomènech P, Caño-Delgado AI: An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L.). BMC Genomics 2009, 10:467.

Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. S Moxon, R Jing, G Szittya, F Schwach, Rusholme Pilcher, R L Moulton, V Dalmay, T , Genome Res. 18Moxon S, Jing R, Szittya G, Schwach F, Rusholme Pilcher RL, Moulton V, Dalmay T: Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res 2008, 18:1602-1609.

Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. C Arenas-Huertero, B Pérez, F Rabanal, D Blanco-Melo, C De La Rosa, G Estrada-Navarrete, F Sanchez, A A Covarrubias, J L Reyes, Plant Mol Biol. 70Arenas-Huertero C, Pérez B, Rabanal F, Blanco-Melo D, De la Rosa C, Estrada-Navarrete G, Sanchez F, Covarrubias AA, Reyes JL: Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. Plant Mol Biol 2009, 70:385-401.

Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). C Song, C Wang, C Zhang, N Korir, H Yu, Z Ma, J Fang, BMC Genomics. 11431Song C, Wang C, Zhang C, Korir N, Yu H, Ma Z, Fang J: Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics 2010, 11:431.

Identification of grapevine microRNAs and their targets using highthroughput sequencing and degradome analysis. V Pantaleo, G Szittya, S Moxon, L Miozzi, V Moulton, T Dalmay, J Burgyan, Plant J. 62Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, Burgyan J: Identification of grapevine microRNAs and their targets using high- throughput sequencing and degradome analysis. Plant J 2010, 62:960-976.

Post-transcriptional small RNA pathways in plants: mechanisms and regulations. H Vaucheret, Genes & Development. 20Vaucheret H: Post-transcriptional small RNA pathways in plants: mechanisms and regulations. Genes & Development 2006, 20:759-771.

Heterochromatic siRNAs and DDM1 Independently Silence Aberrant 5S rDNA Transcripts in Arabidopsis. T Blevins, O Pontes, C S Pikaard, F Meins, PLoS ONE. 4Blevins T, Pontes O, Pikaard CS, Meins F: Heterochromatic siRNAs and DDM1 Independently Silence Aberrant 5S rDNA Transcripts in Arabidopsis. PLoS ONE 2009, 4.

Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. O Borsani, J Zhu, P E Verslues, R Sunkar, J Zhu, Cell. 123Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu J: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell 2005, 123:1279-1291.

The diversity of RNA silencing pathways in plants. P Brodersen, O Voinnet, Trends Genet. 22Brodersen P, Voinnet O: The diversity of RNA silencing pathways in plants. Trends Genet 2006, 22:268-280.

Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. F Vazquez, H Vaucheret, R Rajagopalan, C Lepers, V Gasciolli, A C Mallory, J Hilbert, D P Bartel, P Crété, Mol Cell. 16Vazquez F, Vaucheret H, Rajagopalan R, Lepers C, Gasciolli V, Mallory AC, Hilbert J, Bartel DP, Crété P: Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell 2004, 16:69-79.

. Gonzalez-Ibeas, BMC Genomics. 12Gonzalez-Ibeas et al. BMC Genomics 2011, 12:393 http://www.biomedcentral.com/1471-2164/12/393

microRNA-directed phasing during trans-acting siRNA biogenesis in plants. E Allen, Z Xie, A M Gustafson, J C Carrington, Cell. 121Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 2005, 121:207-221.

Virus-derived small interfering RNAs at the core of plant-virus interactions. C Llave, Trends Plant Sci. 15Llave C: Virus-derived small interfering RNAs at the core of plant-virus interactions. Trends Plant Sci 2010, 15:701-707.

Safeguarding production-losses in major crops and the role of crop protection. E Oerke, H Dehne, Crop Protection. 23Oerke E, Dehne H: Safeguarding production-losses in major crops and the role of crop protection. Crop Protection 2004, 23:275-285.

Emerging pathogens: the epidemiology and evolution of species jumps. M E Woolhouse, D T Haydon, R Antia, Trends in Ecology & Evolution. 20Woolhouse ME, Haydon DT, Antia R: Emerging pathogens: the epidemiology and evolution of species jumps. Trends in Ecology & Evolution 2005, 20:238-244.

Cucurbit aphid-borne yellows virus Is Prevalent in Field-Grown Cucurbit Crops of Southeastern Spain. M A Kassem, R N Sempere, M Juárez, M A Aranda, V Truniger, Plant Disease. 91Kassem MA, Sempere RN, Juárez M, Aranda MA, Truniger V: Cucurbit aphid-borne yellows virus Is Prevalent in Field-Grown Cucurbit Crops of Southeastern Spain. Plant Disease 2007, 91:232-238.

Sources of resistance to viruses in the Potyviridae. R Provvidenti, R O Hampton, Arch Virol Suppl. 5Provvidenti R, Hampton RO: Sources of resistance to viruses in the Potyviridae. Arch Virol Suppl 1992, 5:189-211.

Potential Sources of Resistance for Melon to Nonpersistently Aphid-borne Viruses. J A Díaz, C Mallor, C Soria, R Camero, E Garzo, A Fereres, J M Alvarez, M L Gómez-Guillamón, M Luis-Arteaga, E Moriones, Plant Disease. 87Díaz JA, Mallor C, Soria C, Camero R, Garzo E, Fereres A, Alvarez JM, Gómez- Guillamón ML, Luis-Arteaga M, Moriones E: Potential Sources of Resistance for Melon to Nonpersistently Aphid-borne Viruses. Plant Disease 2003, 87:960-964.

Inheritance of Resistance to Watermelon mosaic virus in Cucumis melo that Impairs Virus Accumulation, Symptom Expression, and Aphid Transmission. J A Diaz-Pendon, R Fernandez-Munoz, M L Gomez-Guillamon, E Moriones, Phytopathology. 95Diaz-Pendon JA, Fernandez-Munoz R, Gomez-Guillamon ML, Moriones E: Inheritance of Resistance to Watermelon mosaic virus in Cucumis melo that Impairs Virus Accumulation, Symptom Expression, and Aphid Transmission. Phytopathology 2005, 95:840-846.

Descr Pl Viruses No. 302 CMI/AAB. T Hibi, I Furuki, Hibi T, Furuki I: Descr Pl Viruses No. 302 CMI/AAB; 1985.

An eIF4E allele confers resistance to an uncapped and non-polyadenylated RNA virus in melon. C Nieto, M Morales, G Orjeda, C Clepet, A Monfort, B Sturbois, P Puigdomènech, M Pitrat, M Caboche, C Dogimont, J Garcia-Mas, M A Aranda, A Bendahmane, Plant J. 48Nieto C, Morales M, Orjeda G, Clepet C, Monfort A, Sturbois B, Puigdomènech P, Pitrat M, Caboche M, Dogimont C, Garcia-Mas J, Aranda MA, Bendahmane A: An eIF4E allele confers resistance to an uncapped and non-polyadenylated RNA virus in melon. Plant J 2006, 48:452-462.

Spanish Melon necrotic spot virus Isolate Overcomes the Resistance Conferred by the Recessive nsv Gene of Melon. J A Díaz, C Nieto, E Moriones, M A Aranda, Plant Disease. 86Díaz JA, Nieto C, Moriones E, Aranda MA: Spanish Melon necrotic spot virus Isolate Overcomes the Resistance Conferred by the Recessive nsv Gene of Melon. Plant Disease 2002, 86:694-694.

Molecular characterization of a Melon necrotic spot virus strain that overcomes the resistance in melon and nonhost plants. J A Díaz, C Nieto, E Moriones, V Truniger, M A Aranda, Mol Plant Microbe Interact. 17Díaz JA, Nieto C, Moriones E, Truniger V, Aranda MA: Molecular characterization of a Melon necrotic spot virus strain that overcomes the resistance in melon and nonhost plants. Mol Plant Microbe Interact 2004, 17:668-675.

Mechanism of plant eIF4E-mediated resistance against a Carmovirus (Tombusviridae): capindependent translation of a viral RNA controlled in cis by an (a) virulence determinant. V Truniger, C Nieto, D González-Ibeas, M Aranda, Plant J. 56Truniger V, Nieto C, González-Ibeas D, Aranda M: Mechanism of plant eIF4E-mediated resistance against a Carmovirus (Tombusviridae): cap- independent translation of a viral RNA controlled in cis by an (a) virulence determinant. Plant J 2008, 56:716-727.

M Margulies, M Egholm, W E Altman, S Attiya, J S Bader, L A Bemben, J Berka, M S Braverman, Y Chen, Z Chen, S B Dewell, L Du, J M Fierro, X V Gomes, B C Goodwin, W He, S Helgesen, C H Ho, G P Irzyk, S C Jando, M L Alenquer, T P Jarvie, K B Jirage, J Kim, J R Knight, J R Lanza, J H Leamon, S M Lefkowitz, M Lei, J Li, Genome Sequencing in Open Microfabricated High Density Picoliter Reactors. 437Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen Y, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Goodwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim J, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, et al: Genome Sequencing in Open Microfabricated High Density Picoliter Reactors. Nature 2005, 437:376-380.

Deepsequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes. L Donaire, Y Wang, D Gonzalez-Ibeas, K F Mayer, M A Aranda, C Llave, Virology. 392Donaire L, Wang Y, Gonzalez-Ibeas D, Mayer KF, Aranda MA, Llave C: Deep- sequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes. Virology 2009, 392:203-214.

ASRP: the Arabidopsis Small RNA Project Database. A M Gustafson, E Allen, S Givan, D Smith, J C Carrington, K D Kasschau, Nucleic Acids Res. 33Gustafson AM, Allen E, Givan S, Smith D, Carrington JC, Kasschau KD: ASRP: the Arabidopsis Small RNA Project Database. Nucleic Acids Res 2005, 33: D637-640.

miRBase: tools for microRNA genomics. S Griffiths-Jones, H K Saini, S Van Dongen, A J Enright, Nucleic Acids Research. 36Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Research 2007, 36:D154-D158.

Evidence that miRNAs are different from other RNAs. B H Zhang, X P Pan, S B Cox, G P Cobb, T A Anderson, Cell Mol Life Sci. 63Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci 2006, 63:246-254.

A uniform system for microRNA annotation. Ambros V Bartel, B Bartel, D P Burge, C B Carrington, J C Chen, X Dreyfuss, G Eddy, S R Griffiths-Jones, S Marshall, M Matzke, M Ruvkun, G Tuschl, T , RNA. 9Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T: A uniform system for microRNA annotation. RNA 2003, 9:277-279.

B C Meyers, M J Axtell, B Bartel, D P Bartel, D Baulcombe, J L Bowman, X Cao, J C Carrington, X Chen, P J Green, S Griffiths-Jones, S E Jacobsen, A C Mallory, R A Martienssen, R S Poethig, Y Qi, H Vaucheret, O Voinnet, Y Watanabe, D Weigel, J Zhu, Criteria for Annotation of Plant MicroRNAs. 20Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu J: Criteria for Annotation of Plant MicroRNAs. Plant Cell 2008, 20:3186-3190.

MicroRNA targets in Drosophila. A J Enright, B John, U Gaul, T Tuschl, C Sander, D S Marks, Genome Biol. 51Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol 2003, 5:R1.

High-Throughput Sequencing of Arabidopsis microRNAs: Evidence for Frequent Birth and Death of MIRNA Genes. N Fahlgren, M D Howell, K D Kasschau, E J Chapman, C M Sullivan, J S Cumbie, S A Givan, T F Law, S R Grant, J L Dangl, J C Carrington, PLoS ONE. 2219Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-Throughput Sequencing of Arabidopsis microRNAs: Evidence for Frequent Birth and Death of MIRNA Genes. PLoS ONE 2007, 2:e219.

Identification of small non-coding RNAs from mitochondria and chloroplasts. B Lung, A Zemann, M J Madej, M Schuelke, S Techritz, S Ruf, R Bock, A Hüttenhofer, Nucleic Acids Res. 34Lung B, Zemann A, Madej MJ, Schuelke M, Techritz S, Ruf S, Bock R, Hüttenhofer A: Identification of small non-coding RNAs from mitochondria and chloroplasts. Nucleic Acids Res 2006, 34:3842-3852.

A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. R Rajagopalan, H Vaucheret, J Trejo, D P Bartel, Genes Dev. 20Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 2006, 20:3407-3425.

Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. R D Morin, G Aksay, E Dolgosheina, H A Ebhardt, V Magrini, E R Mardis, S C Sahinalp, P J Unrau, Genome Res. 18Morin RD, Aksay G, Dolgosheina E, Ebhardt HA, Magrini V, Mardis ER, Sahinalp SC, Unrau PJ: Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. Genome Res 2008, 18:571-584.

Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Y Yao, G Guo, Z Ni, R Sunkar, J Du, J Zhu, Q Sun, Genome Biol. 896Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu J, Sun Q: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biol 2007, 8:R96.

Conifers have a unique small RNA silencing signature. E V Dolgosheina, R D Morin, G Aksay, S C Sahinalp, V Magrini, E R Mardis, J Mattsson, P J Unrau, RNA. 14Dolgosheina EV, Morin RD, Aksay G, Sahinalp SC, Magrini V, Mardis ER, Mattsson J, Unrau PJ: Conifers have a unique small RNA silencing signature. RNA 2008, 14:1508-1515.

Elucidation of the Small RNA Component of the Transcriptome. C Lu, S S Tej, S Luo, C D Haudenschild, B C Meyers, P J Green, Science. 309Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the Small RNA Component of the Transcriptome. Science 2005, 309:1567-15691.

Novel and Stress-Regulated MicroRNAs and Other Small RNAs from Arabidopsis. R Sunkar, J Zhu, Plant Cell. 16Sunkar R, Zhu J: Novel and Stress-Regulated MicroRNAs and Other Small RNAs from Arabidopsis. Plant Cell 2004, 16:2001-2019.

Act in microRNA Metabolism in Arabidopsis thaliana. W Park, J Li, R Song, J Messing, X Chen, CARPEL FACTORY, a Dicer Homolog, and HEN1, a Novel Protein. 12Park W, Li J, Song R, Messing J, Chen X: CARPEL FACTORY, a Dicer Homolog, and HEN1, a Novel Protein, Act in microRNA Metabolism in Arabidopsis thaliana. Current Biology 2002, 12:1484-1495.

Genetic and Functional Diversification of Small RNA Pathways in Plants. Z Xie, L K Johansen, A M Gustafson, K D Kasschau, A D Lellis, D Zilberman, S E Jacobsen, J C Carrington, PLoS Biol. 2104Xie Z, Johansen LK, Gustafson AM, Kasschau KD, Lellis AD, Zilberman D, Jacobsen SE, Carrington JC: Genetic and Functional Diversification of Small RNA Pathways in Plants. PLoS Biol 2004, 2:e104.

Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy. V Gonzalez, A Benjak, E Henaff, G Mir, J Casacuberta, J Garcia-Mas, P Puigdomenech, BMC Plant Biology. 10246Gonzalez V, Benjak A, Henaff E, Mir G, Casacuberta J, Garcia-Mas J, Puigdomenech P: Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy. BMC Plant Biology 2010, 10:246.

A two-hit trigger for siRNA biogenesis in plants. M J Axtell, C Jan, R Rajagopalan, D P Bartel, Cell. 127Axtell MJ, Jan C, Rajagopalan R, Bartel DP: A two-hit trigger for siRNA biogenesis in plants. Cell 2006, 127:565-577.

Foldback transposable elements in plants. D Rebatchouk, J O Narita, Plant Mol Biol. 34Rebatchouk D, Narita JO: Foldback transposable elements in plants. Plant Mol Biol 1997, 34:831-835.

Sequence and Analysis of the Tomato JOINTLESS Locus. L Mao, D Begum, S A Goff, R A Wing, Plant Physiol. 126Mao L, Begum D, Goff SA, Wing RA: Sequence and Analysis of the Tomato JOINTLESS Locus. Plant Physiol 2001, 126:1331-1340.

The DAWGPAWS pipeline for the annotation of genes and transposable elements in plant genomes. J Estill, J Bennetzen, Plant Methods. 58Estill J, Bennetzen J: The DAWGPAWS pipeline for the annotation of genes and transposable elements in plant genomes. Plant Methods 2009, 5:8.

MicroRNAS and their regulatory roles in plants. M W Jones-Rhoades, D P Bartel, B Bartel, Annu Rev Plant Biol. 57Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol 2006, 57:19-53.

The complex interplay between plant viruses and host RNA-silencing pathways. P Dunoyer, O Voinnet, Current Opinion in Plant Biology. 8Dunoyer P, Voinnet O: The complex interplay between plant viruses and host RNA-silencing pathways. Current Opinion in Plant Biology 2005, 8:415-423.

Viral RNA silencing suppressors inhibit the microRNA pathway at an intermediate step. E J Chapman, A I Prokhnevsky, K Gopinath, V V Dolja, J C Carrington, Genes & Development. 18Chapman EJ, Prokhnevsky AI, Gopinath K, Dolja VV, Carrington JC: Viral RNA silencing suppressors inhibit the microRNA pathway at an intermediate step. Genes & Development 2004, 18:1179-1186.

Genome-Wide Profiling and Analysis of Arabidopsis siRNAs. K D Kasschau, N Fahlgren, E J Chapman, C M Sullivan, J S Cumbie, S A Givan, J C Carrington, PLoS Biol. 557Kasschau KD, Fahlgren N, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Carrington JC: Genome-Wide Profiling and Analysis of Arabidopsis siRNAs. PLoS Biol 2007, 5:e57.

A viral suppressor of RNA silencing differentially regulates the accumulation of short interfering RNAs and micro-RNAs in tobacco. A C Mallory, B J Reinhart, D Bartel, V B Vance, L H Bowman, Proceedings of the National Academy of Sciences of the United States of America. the National Academy of Sciences of the United States of America99Mallory AC, Reinhart BJ, Bartel D, Vance VB, Bowman LH: A viral suppressor of RNA silencing differentially regulates the accumulation of short interfering RNAs and micro-RNAs in tobacco. Proceedings of the National Academy of Sciences of the United States of America 2002, 99:15228-15233.

Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. A A Bazzini, H E Hopp, R N Beachy, S Asurmendi, Proceedings of the National Academy of Sciences. 104Bazzini AA, Hopp HE, Beachy RN, Asurmendi S: Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. Proceedings of the National Academy of Sciences 2007, 104:12157-12162.

AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. H Vaucheret, A C Mallory, D P Bartel, Mol Cell. 22Vaucheret H, Mallory AC, Bartel DP: AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell 2006, 22:129-136.

Defective Interfering RNA Hinders the Activity of a Tombusvirus-Encoded Posttranscriptional Gene Silencing Suppressor. Z Havelda, C Hornyik, A Válóczi, J Burgyán, J Virol. 79Havelda Z, Hornyik C, Válóczi A, Burgyán J: Defective Interfering RNA Hinders the Activity of a Tombusvirus-Encoded Posttranscriptional Gene Silencing Suppressor. J Virol 2005, 79:450-457.

Cucumber mosaic virus-encoded 2b suppressor inhibits Arabidopsis Argonaute1 cleavage activity to counter plant defense. X Zhang, Y Yuan, Y Pei, S Lin, T Tuschl, D J Patel, N Chua, Genes Dev. 20Zhang X, Yuan Y, Pei Y, Lin S, Tuschl T, Patel DJ, Chua N: Cucumber mosaic virus-encoded 2b suppressor inhibits Arabidopsis Argonaute1 cleavage activity to counter plant defense. Genes Dev 2006, 20:3255-3268.

Plant virus-mediated induction of miR168 is associated with repression of ARGONAUTE1 accumulation. E Várallyay, A Válóczi, A Agyi, J Burgyán, Z Havelda, EMBO J. 29Várallyay E, Válóczi A, Agyi A, Burgyán J, Havelda Z: Plant virus-mediated induction of miR168 is associated with repression of ARGONAUTE1 accumulation. EMBO J 2010, 29:3507-3519.

. Gonzalez-Ibeas, BMC Genomics. 12Gonzalez-Ibeas et al. BMC Genomics 2011, 12:393 http://www.biomedcentral.com/1471-2164/12/393

Prediction of Plant MicroRNA Targets. M W Rhoades, B J Reinhart, L P Lim, C B Burge, B Bartel, D P Bartel, Cell. 110Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of Plant MicroRNA Targets. Cell 2002, 110:513-520.

Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. M A German, M Pillay, D Jeong, A Hetawal, S Luo, P Janardhanan, V Kannan, L A Rymarquis, K Nobuta, R German, De Paoli, E Lu, C Schroth, G Meyers, B C Green, P J , Nat Biotech. 26German MA, Pillay M, Jeong D, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotech 2008, 26:941-946.

Computational challenges in miRNA target predictions: to be or not to be a true target?. C Barbato, I Arisi, M E Frizzo, R Brandi, Da Sacco, L Masotti, A , J Biomed Biotechnol. 803069Barbato C, Arisi I, Frizzo ME, Brandi R, Da Sacco L, Masotti A: Computational challenges in miRNA target predictions: to be or not to be a true target? J Biomed Biotechnol 2009, 2009:803069.

Inferring MicroRNA Activities by Combining Gene Expression with MicroRNA Target Prediction. C Cheng, L M Li, PLoS ONE. 31989Cheng C, Li LM: Inferring MicroRNA Activities by Combining Gene Expression with MicroRNA Target Prediction. PLoS ONE 2008, 3:e1989.

Python Programming Language -Official Website. Python Programming Language -Official Website. [http://www.python. org/].

. Biopython -Biopython, Biopython -Biopython. [http://biopython.org/wiki/Main_Page].

Basic local alignment search tool. S F Altschul, W Gish, W Miller, E W Myers, D J Lipman, J Mol Biol. 215Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215:403-410.

tRNAdb 2009: compilation of tRNA sequences and tRNA genes. F Jühling, M Mörl, R K Hartmann, M Sprinzl, P F Stadler, J Pütz, Nucleic Acids Res. 37Jühling F, Mörl M, Hartmann RK, Sprinzl M, Stadler PF, Pütz J: tRNAdb 2009: compilation of tRNA sequences and tRNA genes. Nucleic Acids Res 2009, 37:D159-162.

Plant snoRNA database. Jws Brown, M Echeverria, L Qu, T M Lowe, J Bachellerie, A Hüttenhofer, J P Kastenmayer, P J Green, P Shaw, D F Marshall, Nucleic Acids Res. 31Brown JWS, Echeverria M, Qu L, Lowe TM, Bachellerie J, Hüttenhofer A, Kastenmayer JP, Green PJ, Shaw P, Marshall DF: Plant snoRNA database. Nucleic Acids Res 2003, 31:432-435.

SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. E Pruesse, C Quast, K Knittel, B M Fuchs, W Ludwig, J Peplies, F O Glöckner, Nucleic Acids Research. 35Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Research 2007, 35:7188-7196.

Rfam: updates to the RNA families database. P P Gardner, J Daub, J G Tate, E P Nawrocki, D L Kolbe, S Lindgreen, A C Wilkinson, R D Finn, S Griffiths-Jones, S R Eddy, A Bateman, Nucleic Acids Res. 37Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR, Bateman A: Rfam: updates to the RNA families database. Nucleic Acids Res 2009, 37: D136-140.

The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants. S Ouyang, C R Buell, Nucleic Acids Res. 32Ouyang S, Buell CR: The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants. Nucleic Acids Res 2004, 32:D360-D363.

Mfold web server for nucleic acid folding and hybridization prediction. M Zuker, 10.1186/1471-2164-12-393Nucleic Acids Research. 31Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Research 2003, 31:3406-3415. doi:10.1186/1471-2164-12-393

Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing. Gonzalez-Ibeas, BMC Genomics. 12393Cite this article asCite this article as: Gonzalez-Ibeas et al.: Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing. BMC Genomics 2011 12:393.