BACKGROUND
melon  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  <cit> .
the melon genome has  <dig> chromosomes and is thought to contain 450- <dig> mb of dna, which is 3- <dig> times more than arabidopsis  <cit> .
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  <cit> .
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  <dig>  which developed a range of tools and accumulated more than  <dig>  expressed sequence tags  and ~ <dig>  tentative consensus sequences   <cit> .
this est collection has been expanded recently with the addition of  <dig>  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  <cit> .
these ests as well as other resources are now accessible in a public database  <cit> .
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  <cit> .
in addition, a melon sequencing consortium has recently produced a high-quality draft of the melon genome .
although these resources provided significant advances in the analysis of melon gene expression, the small rna  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  <cit> .
srnas are short, non-coding rnas 21- <dig> nucleotides  in length which are found in protists, fungi, plants and animals  <cit> .
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  and small interfering rnas .
mirnas are widely studied because of their regulatory activity, particularly in development, pathogen resistance and stress responses  <cit> .
mirnas are cleaved from stem-loop precursor molecules that derive from single stranded non-coding transcripts.
mirnas regulate protein-coding genes post-transcriptionally 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  arise from and target endogenous loci such as transposons and dna repeats to direct cytosine methylation and chromatin modifications  <cit> .
natural antisense-transcript 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  <cit> .
finally, trans-acting sirnas  derived from tas genes, which transcribe long primary non-coding 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  <cit> .
in addition to endogenous srnas, exogenous sirnas from virus genomes can be detected in virus-infected plants as a part of the rna-based immune system  <cit> .
rna viruses that infect melon are responsible for significant yield losses as well as poor fruit quality  <cit> , particularly the widespread watermelon mosaic virus   <cit> .
recently, a collection of accessions representing cultivated melon and its wild relatives was screened to identify sources of resistance to mosaic-inducing viruses  <cit> .
tgr- <dig> was identified as a resistant accession based on the lower wmv titer compared to susceptible genotypes  and the absence or mildness of the mosaic symptoms normally observed in systemically infected leaves  <cit> .
melon necrotic spot virus , although less economically important, may also cause yield losses, and epidemic outbreaks have been reported worldwide  <cit> .
in melon, resistance to mnsv is controlled by the single recessive gene nsv, which encodes eukaryotic translation initiation factor 4e   <cit> .
this resistance is effective against all mnsv strains  except mnsv- <dig>  <cit> .
studies of chimeric viruses have shown that the mnsv 3' untranslated region  contains the resistance-breaking determinant of mnsv- <dig>  and that it functions as a cap-independent translational enhancer  <cit> .
we constructed  <dig> 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  <dig> technology from roche  <cit> .
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  <dig> and  <dig> 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  from cotyledon samples  <cit> .
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  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 virus-infected melon cotyledons .
pcr amplification products corresponding to each library were pooled in equal amounts and sequences were obtained by multiplexed high-throughput pyrosequencing .
this produced  <dig>  raw sequences, each ~ <dig> bases in length,  <dig>  of which had a complete 3' adaptor in the correct position.
based on these data, we estimated a sequencing error rate of  <dig> %.
after removing reads where one or the two adaptors could not be identified,  <dig>  useful sequences with 3' and 5' adaptors were selected.
only  <dig> 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 .
for instance, the fruit and ovary libraries  were poorly represented providing a collection of fewer than one third of the number of sequences from the other six libraries.
a set of  <dig>  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  <dig>  the most abundant sequences were  <dig>   <dig>   <dig> and  <dig> nts.
a few sequences shorter than  <dig> nt were also retrieved, and these probably represent cloning artifacts and/or degradation products.
sequences >  <dig> nt in length in our dataset predominantly represented combinations of other melon srnas identified in our work.
detailed data are provided in additional file  <dig>
awmv = watermelon mosaic virus; mnsv  = melon necrotic spot virus, alfa <dig> isolate; mnsv  = melon necrotic spot virus, alfa <dig> isolate with 3' utr from  <dig> isolate
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   <cit>  and the microrna database   <cit> .
we identified  <dig> melon unique sequences corresponding to  <dig> mirna families.
thirty nine sequences were identical to known mirnas from other plant species, while  <dig> additional species were sequence variants highly conserved .
in order to clearly identify each melon sequence, melon mirnas were named according to the homologous reference mirna from each database .
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-mir <dig> from populus euphratica, respectively .
the largest diversity of mirna species was found in ovary samples and the lowest in fruit samples.
a sequences with hit in melon genome = 'yes'; sequences with no hit = '--'
the abundance distribution of different mirnas in each library was estimated based on sequencing frequencies as shown in figure  <dig>  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  <dig>  sequences in total.
figure 2a, b compares the accumulation of mirnas in healthy versus wmv-infected melon tissues from genotypes tendral and tgr- <dig>  melon mirna species with similarity to arabidopsis mir156abcdef, mir160abc and mir168ab, which target mrnas encoding squamosa promoter binding proteins, auxin response factors  and argonaute-like proteins , 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- <dig> tissues than in healthy tissues.
other known mirnas 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  revealed that mirnas were particularly abundant and diverse in ovaries compared to fruits.
several mirnas appeared to be temporally regulated during ovary development  whereas others were equally abundant at both ovary stages .
fruits contained far fewer mirnas than ovaries, and only mirnas similar in sequence to arabidopsis mir159a, mir164ab and mir397a showed significant differences in accumulation .
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  while the near-perfect 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  <dig> mirna* sequences complementary to some of the  <dig> mirnas in our dataset , 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 , whereas the other six had a different number of protruding nucleotides and were considered non-typical .
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  <dig> mirna and  <dig> mirna* sequences, as opposed to mir159a for which  <dig>  mirna sequences but no corresponding mirna* sequences were retrieved in the sequenced collections .
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 ,  <dig>  unique melon srnas remained unclassified, from which the most abundant  were 24-nt species.
initial analysis confirmed that  <dig>   of these sequences had a perfect match in the melon genome.
the frequency distribution was highly skewed:  <dig>  sequences had fewer than  <dig> hits , and only  <dig> sequences had more than  <dig> hits.
sequences that were  <dig>   <dig> or  <dig> 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  <dig> 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  <cit> , was calculated for each precursor and the results were sorted accordingly .
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 , up to  <dig> in some cases.
therefore, these were also considered unsuitable mirna candidates.
three other potential mirnas 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,  <dig> loci that fulfilled the structural criteria for annotation of plant mirnas  <cit>  were selected as plausible mirna precursors; we also added to this list  <dig> other loci that had an asymmetric bulge involving  <dig> bases inside the putative mirna duplex.
from them,  <dig>   <dig> and  <dig> corresponded to sequenced srnas of  <dig>   <dig> and  <dig> nt in length, respectively .
six selected sequences are shown as examples in figure  <dig>  by checking the pairing sequence on the stem of the predicted precursors, mirna*s for seven candidate mirnas 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  <cit> .
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  <dig> sequences .
as occurred for conserved mirnas, sequence variants were identified for some novel mirnas  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.
a precursor secondary structures of srnas in bold are represented in figure  <dig>
b miranda score calculated for the identification of the complementary region to the putative mirna.
c mfei index calculated as described  <cit> .
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  <dig> bases was identified in the mirna/mirna* duplex of the precursor, not meeting exactly the structural criteria previously set forth  <cit> .
prediction of potential mirna targets in the melon transcriptome
to identify potential targets of mirnas, we screened melon unigenes in the publicly-available database  <cit> .
two independent searches were performed using miranda  <cit>  and targetfinder  <cit> , 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 ≥ <dig>  and using this value as cutoff we found  <dig> melon unigenes as potential mirna targets, the best of which are summarized in table  <dig> .
the potential mirna targets in table  <dig> 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 mir <dig> regulation in arabidopsis.
many of the melon unigenes identified as potential targets were not annotated, and some had previously been identified as potential mirna precursors  <cit> .
interestingly, the highest miranda scores  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 ~ <dig> nt.
when this region was used as a blast query against the melon srna dataset, a group of  <dig> sequences  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  <cit>  .
both sites  had similar miranda scores, they did not contain mismatches or g:u wobbles involving nucleotides 9- <dig> and were phased  <dig> 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 mir <dig> 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 >  <dig> transcripts with a miranda score >  <dig>  several of these transcripts were annotated as arfs and ubiquitin related gene products .
the remaining unigenes with two predicted mirna target sites listed in table  <dig> were annotated as protein-coding transcripts and no srnas were identified with similarity to the region flanked by the two target sites , 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  <cit> .
twenty-eight unigenes were identified as potential mirna targets using the same criteria described above, most of which were found to be non-annotated .
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 ≥ <dig>  were identified for each of the potential novel melon-specific mirnas 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 .
srnas similar to transfer rna , trans-acting sirna, small nucleolar rna  and transposons were the least abundant, whereas ribosomal rnas  were largely the most abundant non-coding srna species .
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 database , even if up to two mismatches were allowed in blast comparisons.
the only exception were sequences similar to ta-sirnas, for which  <dig> melon srnas with similarity to arabidopsis tas3a|d <dig> and tas3a|d <dig> were identified,  <dig> containing  <dig> mismatch, and  <dig> containing  <dig> mismatches.
all of them mapped very close in the melon genome and in a different region than c15d_05-d02-m13r_c unigene .
to determine if they were authentic melon ta-sirnas, we selected a  <dig> 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  <dig> sequences  mapped in this region and were arranged according to a near 21-nt phase spacing .
many srna sequences also generated hits in the plastid genomes .
when these srnas were mapped onto the melon chloroplast genome  , two clusters of sequences resolved in regions presumably annotated as chloroplast rrna.
these regions lie within two inverted genomic repeats, and srnas were accordingly identified in both the sense and antisense orientations.
some of the chloroplast srnas had previously been cloned in other species  <cit> .
for example, melon srna a33_ <dig>  was the most abundant melon plastid srna , and a matching sequence is present in more than  <dig> 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_  was less abundant , it is present in  <dig> 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 .
for example, in the virus-resistant melon accession tgr- <dig> there was no difference in the number of srnas with hits to melon chloroplast genome between healthy and virus inoculated samples, but in the virus-susceptible accession tendral, more srnas were counted in inoculated samples .
unlike chloroplast srnas, only  <dig>  srnas  matched the melon mitochondrial genome .
these srnas were mapped on the mitochondrial genome sequence and formed three clusters, again corresponding to the sites of rrna genes .
discussion
in this report, we describe the first screen for melon srnas by deep sequencing.
in total,  <dig>  high-quality sequences were generated, representing 90% of the total raw reads.
rna species  <dig>   <dig>   <dig> and  <dig> 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 dcl <dig> are often the most abundant endogenous plant srnas  <cit> , but this may vary according to species.
for example, 24-nt srnas are more abundant in arabidopsis, rice and tomato  <cit> , whereas 21-nt srnas are more abundant in grapevine, wheat and conifers  <cit> .
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 .
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  <cit> .
in melon, only ~ <dig>  srnas matched known transposon sequences, in contrast to ~ <dig>  srna sequences matching ribosomal rna, which may simply reflect the paucity of melon transposon sequence information in databases, as only  <dig> % of the melon genome has been annotated for transposable elements  <cit> .
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 tas <dig> genes.
non-coding transcripts containing two mir <dig> 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' mir <dig> complementary site contains a mismatch and two g:u wobbles involving positions 9- <dig> and, despite it is not cleaved, it binds the silencing complex and is required for full attas <dig> function in vivo  <cit> .
in melon, both 3' and 5' mir <dig> had perfect complementarity at positions 9- <dig>  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 , as proposed for alternatively phased tas <dig> ta-sirnas in arabidopsis  <cit> .
since there are additional tas loci in other plant genomes, it is reasonable that other melon tas loci remain to be discovered.
more than  <dig>  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- <dig> cotyledons  whereas more srna accumulated in healthy tendral  cotyledons than in infected ones.
whether or not this is related with the resistance phenotype is a matter of speculation.
more than  <dig>  melon srnas in our sequenced collections matched known mirnas in other plants, and  <dig> distinct melon srna species could be assigned to  <dig> 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   <cit> .
this probably reflects the generally-accepted high level of expression reported for conserved mirnas.
in addition to known mirnas,  <dig> srna sequences derived from genomic loci with intramolecular folding capacities and not previously described as mirnas in other plant species were predicted as potential melon-specific 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  <cit> .
the candidates listed in table  <dig> include a number of special cases, i.e.
mirnas with miranda scores ≥ <dig> and very strong secondary structures including an internal loop, resembling type iii foldback transposons  <cit> .
although these sequences do not match known melon transposons, they were not considered as mirna candidates because accurate homology-based transposon annotation and prediction occasionally needs to be complemented with ab initio approaches based on structural features  <cit> .
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 over-expression or knockout of mirna genes led to severe developmental defects  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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 argonaute <dig> , the catalytic subunit of the rna-induced silencing complex responsible for slicing of target mrnas  <cit> .
recent work has described the enhanced expression of mir <dig> and ago <dig> mrna in virus-infected plants specifically and independently of other mirnas  <cit> .
the contrasting mirna profiles observed in the tgr- <dig> and tendral varieties suggests that silencing may underly the resistance of tgr- <dig> to wmv, although this is a hypothesis that will require further research.
we have identified more than  <dig> 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  <cit> .
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 non-conserved targets of conserved mirnas can be cleaved at a lower frequency than conserved targets  <cit> .
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  <cit> , 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 mirnas and their target mrnas can be identified  <cit> .
we have used microarrays to monitor gene expression profiles in healthy tgr- <dig> 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 .
the same profile was observed for mir <dig> accumulation, suggesting that mir <dig> may be involved in virus resistance and providing the basis for future experiments.
CONCLUSIONS
we have analysed and catalogued a collection of melon endogenous srna obtained through massive cdna sequencing and have identified known mirnas and ta-sirnas  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  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.

 melon 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. the melon genome has chromosomes and is thought to contain 450- mb of dna, which is 3- times more than arabidopsis. 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. 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 which developed a range of tools and accumulated more than expressed sequence tags and ~ tentative consensus sequences. this est collection has been expanded recently with the addition of 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. these ests as well as other resources are now accessible in a public database. 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. in addition, a melon sequencing consortium has recently produced a high-quality draft of the melon genome. although these resources provided significant advances in the analysis of melon gene expression, the small rna 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. srnas are short, non-coding rnas 21- nucleotides in length which are found in protists, fungi, plants and animals. 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 and small interfering rnas. mirnas are widely studied because of their regulatory activity, particularly in development, pathogen resistance and stress responses. mirnas are cleaved from stem-loop precursor molecules that derive from single stranded non-coding transcripts. mirnas regulate protein-coding genes post-transcriptionally 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 arise from and target endogenous loci such as transposons and dna repeats to direct cytosine methylation and chromatin modifications. natural antisense-transcript 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. finally, trans-acting sirnas derived from tas genes, which transcribe long primary non-coding 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. in addition to endogenous srnas, exogenous sirnas from virus genomes can be detected in virus-infected plants as a part of the rna-based immune system. rna viruses that infect melon are responsible for significant yield losses as well as poor fruit quality, particularly the widespread watermelon mosaic virus. recently, a collection of accessions representing cultivated melon and its wild relatives was screened to identify sources of resistance to mosaic-inducing viruses. tgr- was identified as a resistant accession based on the lower wmv titer compared to susceptible genotypes and the absence or mildness of the mosaic symptoms normally observed in systemically infected leaves. melon necrotic spot virus, although less economically important, may also cause yield losses, and epidemic outbreaks have been reported worldwide. in melon, resistance to mnsv is controlled by the single recessive gene nsv, which encodes eukaryotic translation initiation factor 4e. this resistance is effective against all mnsv strains except mnsv-. studies of chimeric viruses have shown that the mnsv 3' untranslated region contains the resistance-breaking determinant of mnsv- and that it functions as a cap-independent translational enhancer. we constructed 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 technology from roche. 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 and 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 from cotyledon samples. 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.  cdna libraries and sequencing of small rnas we used high throughput sequencing data to analyze the composition of the small rna transcriptome 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 virus-infected melon cotyledons. pcr amplification products corresponding to each library were pooled in equal amounts and sequences were obtained by multiplexed high-throughput pyrosequencing. this produced raw sequences, each ~ bases in length, of which had a complete 3' adaptor in the correct position. based on these data, we estimated a sequencing error rate of %. after removing reads where one or the two adaptors could not be identified, useful sequences with 3' and 5' adaptors were selected. only 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. for instance, the fruit and ovary libraries were poorly represented providing a collection of fewer than one third of the number of sequences from the other six libraries. a set of 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 the most abundant sequences were and nts. a few sequences shorter than nt were also retrieved, and these probably represent cloning artifacts and/or degradation products. sequences > nt in length in our dataset predominantly represented combinations of other melon srnas identified in our work. detailed data are provided in additional file awmv = watermelon mosaic virus; mnsv = melon necrotic spot virus, alfa isolate; mnsv = melon necrotic spot virus, alfa isolate with 3' utr from isolate 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 and the microrna database. we identified melon unique sequences corresponding to mirna families. thirty nine sequences were identical to known mirnas from other plant species, while additional species were sequence variants highly conserved. in order to clearly identify each melon sequence, melon mirnas were named according to the homologous reference mirna from each database. 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-mir from populus euphratica, respectively. the largest diversity of mirna species was found in ovary samples and the lowest in fruit samples. a sequences with hit in melon genome ='yes'; sequences with no hit ='--' the abundance distribution of different mirnas in each library was estimated based on sequencing frequencies as shown in figure 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 sequences in total. figure 2a, b compares the accumulation of mirnas in healthy versus wmv-infected melon tissues from genotypes tendral and tgr- melon mirna species with similarity to arabidopsis mir156abcdef, mir160abc and mir168ab, which target mrnas encoding squamosa promoter binding proteins, auxin response factors and argonaute-like proteins, 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- tissues than in healthy tissues. other known mirnas 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 revealed that mirnas were particularly abundant and diverse in ovaries compared to fruits. several mirnas appeared to be temporally regulated during ovary development whereas others were equally abundant at both ovary stages. fruits contained far fewer mirnas than ovaries, and only mirnas similar in sequence to arabidopsis mir159a, mir164ab and mir397a showed significant differences in accumulation. 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 while the near-perfect 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 mirna* sequences complementary to some of the mirnas in our dataset, 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, whereas the other six had a different number of protruding nucleotides and were considered non-typical. 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 mirna and mirna* sequences, as opposed to mir159a for which mirna sequences but no corresponding mirna* sequences were retrieved in the sequenced collections. 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, unique melon srnas remained unclassified, from which the most abundant were 24-nt species. initial analysis confirmed that of these sequences had a perfect match in the melon genome. the frequency distribution was highly skewed: sequences had fewer than hits, and only sequences had more than hits. sequences that were or 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 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, was calculated for each precursor and the results were sorted accordingly. 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, up to in some cases. therefore, these were also considered unsuitable mirna candidates. three other potential mirnas 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, loci that fulfilled the structural criteria for annotation of plant mirnas were selected as plausible mirna precursors; we also added to this list other loci that had an asymmetric bulge involving bases inside the putative mirna duplex. from them, and corresponded to sequenced srnas of and nt in length, respectively. six selected sequences are shown as examples in figure by checking the pairing sequence on the stem of the predicted precursors, mirna*s for seven candidate mirnas 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. 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 sequences. as occurred for conserved mirnas, sequence variants were identified for some novel mirnas 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. a precursor secondary structures of srnas in bold are represented in figure b miranda score calculated for the identification of the complementary region to the putative mirna. c mfei index calculated as described. 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 bases was identified in the mirna/mirna* duplex of the precursor, not meeting exactly the structural criteria previously set forth. prediction of potential mirna targets in the melon transcriptome to identify potential targets of mirnas, we screened melon unigenes in the publicly-available database. two independent searches were performed using miranda and targetfinder, 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 ≥ and using this value as cutoff we found melon unigenes as potential mirna targets, the best of which are summarized in table. the potential mirna targets in table 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 mir regulation in arabidopsis. many of the melon unigenes identified as potential targets were not annotated, and some had previously been identified as potential mirna precursors. interestingly, the highest miranda scores 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 ~ nt. when this region was used as a blast query against the melon srna dataset, a group of sequences 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. both sites had similar miranda scores, they did not contain mismatches or g:u wobbles involving nucleotides 9- and were phased 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 mir 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 > transcripts with a miranda score > several of these transcripts were annotated as arfs and ubiquitin related gene products. the remaining unigenes with two predicted mirna target sites listed in table were annotated as protein-coding transcripts and no srnas were identified with similarity to the region flanked by the two target sites, 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. twenty-eight unigenes were identified as potential mirna targets using the same criteria described above, most of which were found to be non-annotated. 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 ≥ were identified for each of the potential novel melon-specific mirnas 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. srnas similar to transfer rna, trans-acting sirna, small nucleolar rna and transposons were the least abundant, whereas ribosomal rnas were largely the most abundant non-coding srna species. 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 database, even if up to two mismatches were allowed in blast comparisons. the only exception were sequences similar to ta-sirnas, for which melon srnas with similarity to arabidopsis tas3a|d and tas3a|d were identified, containing mismatch, and containing mismatches. all of them mapped very close in the melon genome and in a different region than c15d_05-d02-m13r_c unigene. to determine if they were authentic melon ta-sirnas, we selected a 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 sequences mapped in this region and were arranged according to a near 21-nt phase spacing. many srna sequences also generated hits in the plastid genomes. when these srnas were mapped onto the melon chloroplast genome, two clusters of sequences resolved in regions presumably annotated as chloroplast rrna. these regions lie within two inverted genomic repeats, and srnas were accordingly identified in both the sense and antisense orientations. some of the chloroplast srnas had previously been cloned in other species. for example, melon srna a33_ was the most abundant melon plastid srna, and a matching sequence is present in more than 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_ was less abundant, it is present in 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. for example, in the virus-resistant melon accession tgr- there was no difference in the number of srnas with hits to melon chloroplast genome between healthy and virus inoculated samples, but in the virus-susceptible accession tendral, more srnas were counted in inoculated samples. unlike chloroplast srnas, only srnas matched the melon mitochondrial genome. these srnas were mapped on the mitochondrial genome sequence and formed three clusters, again corresponding to the sites of rrna genes. discussion in this report, we describe the first screen for melon srnas by deep sequencing. in total, high-quality sequences were generated, representing 90% of the total raw reads. rna species and 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 dcl are often the most abundant endogenous plant srnas, but this may vary according to species. for example, 24-nt srnas are more abundant in arabidopsis, rice and tomato, whereas 21-nt srnas are more abundant in grapevine, wheat and conifers. 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. 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. in melon, only ~ srnas matched known transposon sequences, in contrast to ~ srna sequences matching ribosomal rna, which may simply reflect the paucity of melon transposon sequence information in databases, as only % of the melon genome has been annotated for transposable elements. 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 tas genes. non-coding transcripts containing two mir 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' mir complementary site contains a mismatch and two g:u wobbles involving positions 9- and, despite it is not cleaved, it binds the silencing complex and is required for full attas function in vivo. in melon, both 3' and 5' mir had perfect complementarity at positions 9- 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, as proposed for alternatively phased tas ta-sirnas in arabidopsis. since there are additional tas loci in other plant genomes, it is reasonable that other melon tas loci remain to be discovered. more than 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- cotyledons whereas more srna accumulated in healthy tendral cotyledons than in infected ones. whether or not this is related with the resistance phenotype is a matter of speculation. more than melon srnas in our sequenced collections matched known mirnas in other plants, and distinct melon srna species could be assigned to 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. this probably reflects the generally-accepted high level of expression reported for conserved mirnas. in addition to known mirnas, srna sequences derived from genomic loci with intramolecular folding capacities and not previously described as mirnas in other plant species were predicted as potential melon-specific 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. the candidates listed in table include a number of special cases, i.e. mirnas with miranda scores ≥ and very strong secondary structures including an internal loop, resembling type iii foldback transposons. although these sequences do not match known melon transposons, they were not considered as mirna candidates because accurate homology-based transposon annotation and prediction occasionally needs to be complemented with ab initio approaches based on structural features. 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 over-expression or knockout of mirna genes led to severe developmental defects. 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. 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. 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. 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 argonaute, the catalytic subunit of the rna-induced silencing complex responsible for slicing of target mrnas. recent work has described the enhanced expression of mir and ago mrna in virus-infected plants specifically and independently of other mirnas. the contrasting mirna profiles observed in the tgr- and tendral varieties suggests that silencing may underly the resistance of tgr- to wmv, although this is a hypothesis that will require further research. we have identified more than 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. 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 non-conserved targets of conserved mirnas can be cleaved at a lower frequency than conserved targets. 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, 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 mirnas and their target mrnas can be identified. we have used microarrays to monitor gene expression profiles in healthy tgr- 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. the same profile was observed for mir accumulation, suggesting that mir may be involved in virus resistance and providing the basis for future experiments.  we have analysed and catalogued a collection of melon endogenous srna obtained through massive cdna sequencing and have identified known mirnas and ta-sirnas 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 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. 
in silico analysis revealed more than potential targets for the conserved and novel mirnas. we catalogued and analysed the melon srnas, resulting in the identification of known mirna families, the prediction of melon-specific mirna candidates, the identification of trans-acting sirnas, and the identification of chloroplast, mitochondrion and transposon-derived srnas. we constructed srna libraries from two stages of developing ovaries, fruits and photosynthetic cotyledons infected with viruses, and carried out high-throughput pyrosequencing. in plants, they regulate gene expression and heterochromatin assembly, and control protection against virus infection. 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.
 melon 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. these are short, non-coding rnas 21- 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.  we constructed 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 known mirna families, the prediction of 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 potential targets for the conserved and novel mirnas.  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. 