BACKGROUND
durum wheat  is the most widespread crop in the mediterranean basin, since it is the base for the production of pasta  <cit>  and other typical products of this area such as breads and cuscus.
durum wheat is an allotetraploid species  which derives from the cross between t. urartu  and an ancient relative of aegilops speltoides tausch.
.
in spite of the economic importance of durum wheat, there is still a lot to know about the genomics and transcriptomics of this species.
in particular, information on the micrornaome of durum wheat is lacking.
mirnas are small non-coding rna molecules which downregulate gene expression by mrna degradation or translational repression.
several conserved mirnas regulate genes encoding transcription factors with important roles in a variety of processes in plants, including protein degradation, response to biotic and abiotic stresses, signal transduction, and developmental processes .
mirnas are well studied in the reproductive parts of plants, where they regulate flowering time  <cit> , seed  <cit>  and flower  <cit>  development.
the key role of mirnas in developmental processes has been widely demonstrated in model plants and crops including some species of the poaceae family.
in rice, mir <dig> regulates flowering time and temporal expression changes of numerous genes during leaf development  <cit> , whereas mir <dig> and mir <dig> were detected to have an important role in rice seed development  <cit> .
in barley, mir <dig> potentially targets the auxin response factors   <dig> and  <dig>  which have critical roles in plant root development  <cit> .
additionally, barley mir <dig> targets a plant-specific scarecrow-like transcription factor  <dig>  which regulates several plant developmental processes  <cit> .
in maize, it has been proved that mir <dig> downregulates glossy <dig> to promote vegetative phase change  <cit> .
the widespread use of high throughput sequencing has promoted advancements in the genomics of plants, including those species with large, repetitive and polyploid genomes as wheat  <cit> .
the recent availability of the hexaploid bread wheat  draft genome has sped up the discovery of small rnas in this species in the past few years .
thus, so far,  <dig> mature mirnas are annotated for t. aestivum in the mirna repository mirbase .
several studies have uncovered development-associated mirnas in bread wheat mainly in the seeds, thus being potentially involved in the regulation of grain-filling  <cit> , and preferentially expressed in developing or germinating seeds  <cit> .
conversely, a limited number of studies have been carried out on mirnas in tetraploid wheats.
for instance, mirna expression patterns in response to drought stress were assessed in leaf and root of wild emmer wheat , the ancestor of durum wheat, through a plant mirna microarray platform and a homology-based approach  <cit> .
additionally, the same group performed an in silico identification of several mirnas from root transcriptome of t. turgidum ssp.
durum and t. turgidum ssp.
dicoccoides  <cit> .
in spite of this, the discovery of mirnas in durum wheat is still largely unknown, with just one mirna reported in mirbase release  <dig>
the aim of this study was to identify conserved and novel mirnas in durum wheat, and evaluate their possible involvement in plant development.
to this purpose, two small rna libraries from the cultivars ciccio and svevo were generated from a pool of different tissues at the late milk stage  <cit> .
following high throughput sequencing, conserved and novel mirnas and their target genes were identified.
some mirnas were predicted to regulate genes involved in plant development or in important biochemical processes such as photosynthesis.
interestingly, we found three novel mirnas derived from ribosomal rna in durum wheat.
in order to investigate the possible involvement of selected conserved and novel mirnas and target genes in the development of leaves and roots, quantitative pcr  was performed to highlight differences in expression profiles.
this work provides useful information that may be employed in breeding programs aiming at improving durum wheat development.
methods
plant material and growth conditions
seeds of two durum wheat cultivars, ciccio and svevo, were sterilized for 2 min in ethanol and then for 3 min in 2 % sodium hypochlorite.
after 12 h in water aired with an oxygen air pump, seeds were transferred into  <dig> mm petri dishes with a wet filter paper and kept for  <dig> days at  <dig> °c.
once the seed hypocotyl became half centimetre, plantlets were transferred into  <dig> cm diameter net pots containing a mix of 80 % agriperlite and 20 % expanded clay, and placed in a greenhouse at 18–20 °c under  <dig> h photo-period.
plants were grown in a hydroponic system containing  <dig> mm mgso <dig> h2o,  <dig> mm k2so <dig>   <dig> mm kh2po <dig>  2mm ca <dig> h2o,  <dig> μm 6mo7o <dig> h2o, 10μm h3bo <dig>  1μm mnso <dig>  1μm znso <dig>  5μm cuso <dig>  100μm fe-edta and a ph =  <dig>  roots were submerged all the time in the solution and this was refreshed every two days.
at the z <dig> stage of seedling growth  <cit> , roots and leaves were collected from three individual plants from ciccio and svevo separately.
at the late milk developmental stage , roots, leaves, flag leaves, and spikes were collected as above.
all tissues were immediately frozen in liquid nitrogen and stored at −80 °c.
small rna extraction and library preparation
small rnas were isolated from  <dig> mg of each tissue from each durum wheat plant at the stage z <dig>  using the mirpremier microrna isolation kit .
rna quality was verified by agarose gel according to manufacturer’s instructions.
the small rna fractions from the four tissues of the three plants were pooled for library preparation.
truseq small rna sample preparation  protocol was used for the construction of two small rna libraries from ciccio and svevo tissues, respectively.
libraries were sent for sequencing to the institute of applied genomics .
data analysis and identification of conserved and novel mirna
all reads were filtered and trimmed in order to discard low quality sequences and adapters.
unique small rnas were counted eliminating redundancies and each count was normalized and expressed as transcript per million .
in order to identify conserved mirnas, unique reads were subjected to similarity search using all plant mature mirnas from mirbase  as reference database.
to this purpose, blastn was used with parameters optimized for short sequences.
only perfect matching srna were extracted and sequences with at least  <dig> counts in one of the two libraries were considered for further analysis.
potential novel mirnas were identified using a custom pipeline.
the potential precursor sequences of novel durum wheat mirnas were obtained by blastn search against all the ests from the ncbi database.
this similarity search was also extended to t. aestivum unigenes in the dfci gene index , version  <dig> .
all sequences identified were subjected to unafold  <cit>  to predict their secondary structure.
to select conserved and novel mirnas in durum wheat we used the criteria for plant mirna annotation described by meyers et al.
<cit> .
qpcr quantification of micrornas
for the qpcr validation of selected mirnas, svevo was chosen since it is one of the most studied durum wheat varieties, and is well-known for its qualitative properties, e.g.
high protein content and high yellow index  <cit> .
small rnas were extracted using the mirpremier microrna isolation kit  and mirnas were detected using stem-loop rt-pcr method  <cit> .
small rna quality was checked by electrophoresis through a 1 % agarose gel containing 1 % gel red .
quantitative qpcr analyses were performed using the universal probe library  probe assay specific for mature mirna expression, described by varkonyi-gasic and colleagues  <cit> .
each reaction was performed in a final volume of  <dig> μl containing  <dig> μl of rotor-gene probe pcr master mix 2x ,  <dig>  μl of mirna-specific forward primer  and the universal reverse primer  ,  <dig>  μl of upl, probe # <dig> ,  <dig>  μl of nuclease-free water, and  <dig>  µl of the respective retrotranscribed mirna.
primers for gene expression analysis were designed based on sequences obtained from the library analysis .
pcr conditions were: 95 °c for  <dig> min, followed by  <dig> cycles of 95 °c for 5 s, 60 °c for 10 s and 72 °c for 1 s. reactions were performed using a rotor-gene  <dig> machine .
three technical replications of three different plants and a no-template control were included for each experiment.
the housekeeping mirna was chosen based on expression data using qpcr of eight different putative reference mirnas  previously described in wheat  <cit> .
data were analysed using genorme and normfinder software and mir <dig> showed to be the most stable mirna in young and mature tissues of svevo roots and leaves in our experimental conditions.
relative quantification of each single mirna expression was analysed using the comparative ct method as described by livak and schmittgen  <cit> .
differentially expressed mirnas were identified as those showing > 2-fold change between the two developmental stages.
after performing independent student’s t test, fold changes with a p-value <  <dig>  were considered statistically significant.
prediction of ttu-mirna target genes
novel and conserved durum wheat mirnas were used as query in the online tool psrnatarget   <cit> .
the analysis was performed using three different transcript libraries for target search: durum wheat ests from the ncbi database, t. aestivum unigenes from dfci gene index  version  <dig>  and hordeum vulgare unigenes from dfci gene index  version  <dig>  blast2go  <cit>  was used to obtain go annotation of the wheat unigenes based on blastx hits against the ncbi nr database.
go functional classification was obtained executing mapping and annotation steps in blastgo with default parameters.
two graphs representing go terms categories  were determined for all durum wheat mirna targets.
expression analysis of target genes
the expression of mirna putative target genes was analysed using the sybr green dye and the rotor-gene  <dig> machine .
primers for gene expression analysis are found in additional file  <dig>  total rna was isolated using the rneasy mini kit  and reverse transcribed into cdna, using the script cdna synthesis kit  according to the manufacturer’s instructions.
before expression analysis, the same primers were used for sequencing of the target gene fragments including the mirna binding site.
real-time qpcr experiments were carried out using  <dig> ng of cdna and the iq sybr green supermix  following the manufacturer’s protocol.
pcr conditions were: 95 °c for  <dig> min, followed by  <dig> cycles of 95 °c for 10 s and 58 °c for 30 s. the reference gene was chosen by testing three putative reference genes, i.e.
actin, glyceraldehyde-3-phosphate dehydrogenase  and cell division control protein, aaa-superfamily of atpases , through a qpcr expression analysis in leaves and roots of svevo young and mature tissues.
data were analysed using genorme and normfinder software suggesting cdc as the most stable gene in our experimental conditions.
therefore, using cdc as the reference gene, relative quantification of each single gene was performed using the comparative ct method as described by livak and schmittgen  <cit> .
differentially expressed genes were identified as those showing > 2-fold change between the two developmental stages.
after performing independent student’s t test, fold changes with a p-value <  <dig>  were considered statistically significant.
RESULTS
sequencing data analysis and identification of durum wheat conserved and novel mirnas
in order to determine the expression of the micrornaome in durum wheat, two small rna libraries were generated and sequenced from a mix of leaves, flag leaves, roots and spikes of ciccio  and svevo  plants, separately, at the developmental stage of late milk .
conserved and novel mirnas and their target genes were identified using bioinformatics tools.
after removing the low quality reads and adaptors, a total of  <dig> , <dig> and  <dig> , <dig> reads were collected from ciccio and svevo, respectively.
the size distribution of the total clean reads in the two libraries  shows a canonical trend, with the 21– <dig> nucleotide long srnas representing the majority of the total reads.
the  <dig> nucleotide peak was the most abundant one, as observed in most plant species.
reads of the two libraries were searched against ncrnas in rfam  to filter out the snornas, snrnas, trnas, and rrnas before annotation with conserved plant mirnas.
as a result,  <dig> , <dig> and  <dig> , <dig> reads were discarded from the library of ciccio and svevo, respectively.
in order to identify conserved mirnas in durum wheat, all the unique reads longer than  <dig> or shorter than  <dig> nucleotides were mapped to known plant mirnas in mirbase release  <dig>  a total of  <dig> known mirnas were found in the two libraries,  <dig> in ciccio and  <dig> in svevo.
out of the conserved mirnas identified,  <dig> belonged to  <dig> conserved mirna families, whereas  <dig> were unassigned .
all the remaining reads, which were not annotated with mirbase, were used to identify potential new mirnas using bioinformatics tools.
based on the analysis of the predicted secondary structure of durum wheat sequences,  <dig> potential mirnas were considered as novel mirna candidates .
this analysis revealed three non-canonical putative mirnas .
similarity search of their precursor sequences matched to ribosomal rna, which has never been reported as a source of plant mirnas before.
the analysis of the secondary structures and their predicted thermodynamic stability meet all criteria for in silico mirna prediction  <cit> .
moreover, these novel mirnas were detected in independent ngs libraries, both the ones used in this study, and other durum wheat srna proprietary libraries .
two non-canonical mirnas  were mapped to the two opposite arms of the same sequence .
the precursor sequences of the novel mirnas ttu-novel- <dig>  and ttu-novel-42/ttu-novel- <dig>  were mapped to 18s rrna and 28s rrna, respectively.
in order to verify whether other mirnas present in mirbase could derive from rrna, all plant mirna precursors were blasted against the nucleotide database of viridiplantae in ncbi.
at least one mirna, peu-mir <dig> from populus euphratica matched with 18s rna, showing 96 % identity with idesia polycarpa gene for 18s ribosomal rna, as the best hit.fig.
<dig> pre-mirna secondary structure of ttu-novel- <dig>  and ttu-novel- <dig> .
5′-end is marked within a circle
expression analysis of selected micrornas
several of the conserved mirnas found were reported to regulate development in a range of plant species.
in order to highlight a possible involvement of these mirnas in durum wheat cv.
svevo leaf and root development, we selected  <dig> of those for the analysis of relative expression by quantitative pcr  at an early stage  and an advanced late milk developmental stage   <cit> .
additionally, five novel mirnas were selected based on their predicted target or according to their origin from rrna molecules.
this set of novel mirnas was analysed from roots and leaves of svevo plants in the two different stages.
in order to choose a suitable housekeeping mirna, based on the literature  <cit> , eight mirnas were selected and evaluated as putative housekeeping genes in our experiment: ttu-mir <dig>  ttu-mir <dig>  ttu-mir <dig>  ttu-mir <dig>  ttu-mir <dig>  ttu-mir <dig>  ttu-mir <dig>  and ttu-mir <dig>  the qpcr expression profiles of these mirnas were analysed for their stability, and ttu-mir <dig> proved to be the most stable in our experimental conditions.
therefore, ttu-mir <dig> was chosen from this group as the housekeeping mirna for our analyses.
results of mirna expression analysis show that six conserved mirnas  were significantly downregulated in leaves at the mature stage compared to the leaves at the early stage, with fold changes ranging from  <dig>   to  <dig>  .
additionally, ttu-mir167h and ttu-mir393a were found to be significantly upregulated in leaves at the adult stage compared to the early stage, with fold changes of  <dig>  and  <dig> , respectively .
for ttu-mir319f, although a higher expression level was observed in the z <dig> leaves, this value was not significant .
on the other hand, in the roots, five conserved mirnas  were found to be significantly downregulated at the z <dig> stage compared to the z <dig> stage, with fold changes ranging from  <dig>   to  <dig>  .
ttu-mir166k was also strongly downregulated in roots , but with a p-value of  <dig>  .
in contrast, ttu-mir156n, ttu-mir167h and ttu-mir444d were found significantly up regulated in elder roots, with fold changes of  <dig> ,  <dig> , and  <dig> , respectively .fig.
<dig> relative expression analysis of selected mirnas and target genes.
conserved and novel t. turgidum mirnas in leaves  and roots  at the early stage  and at the mature stage .
expression levels are given as fold change of mature stage versus early stage.
data are means ± se.
p values between three early stage replicates and three mature stage replicates for each mirna were calculated using unpaired, 2-tailed student’s t test
the analysis carried out for selected novel mirnas confirmed the expression in durum wheat plants of the two putative rrna-derived mirnas previously identified, ttu-novel- <dig> and ttu-novel- <dig>  in details, quantitative assay revealed that in leaves, two novel mirnas, ttu-novel- <dig>  ttu-novel- <dig> were significantly downregulated at the stage z <dig>  compared to the early stage z <dig>  with fold changes of  <dig>  and  <dig> , respectively .
in roots, ttu-novel- <dig> and ttu-novel- <dig> were strongly downregulated .
ttu-novel- <dig> showed a downregulation with a fold change of  <dig>  and  <dig>  in leaves and roots, respectively, but these values were not significant.
for ttu-novel- <dig> and ttu-novel- <dig> the expression levels observed at the two developmental stages were comparable for both the tissues analysed.
prediction of mirna specific target genes
in order to predict ttu-mirna target genes, durum wheat ests, bread wheat unigenes and barley unigenes were used as distinct datasets, and subsequently results were analysed separately.
sequences carrying a region of high complementarity with the ttu-mirnas were identified in the three databases, which were analysed with different priority starting from t. turgidum followed by t. aestivum and h. vulgare.
a number of  <dig>  target genes were predicted for all durum wheat mirnas identified in this work, for a total of  <dig>  and  <dig>  sequences, matching conserved and novel durum wheat mirnas, respectively .
mirna target genes predicted using wheat dfci gene index were subjected to in silico functional analysis using gene ontology  annotation.
the assignment of go terms highlighted different biological processes and molecular functions for novel and conserved mirna targets.
according to terms description relative to molecular function, the majority of the targets interact with nucleic acids.
in fact, for both novel and conserved durum wheat mirna target genes, nucleotide binding or nucleoside phosphate binding were major classes with approximately 15 % of total sequences.
with regard to the other go terms, there is a certain overlap of function in the comparison between novel and conserved mirna targets, with some significant differences such as cation binding including  <dig>  and  <dig>  % targets for conserved and novel mirnas, respectively .
with regard to the biological processes, the majority of the targets were annotated to terms related to metabolic processes.
go annotations partially overlap for novel and conserved mirna targets, with the exception of the class “response to chemical”, which is missing in novel mirna targets .fig.
<dig> bar charts representing biological processes and molecular function categorization according to go terms associated to durum wheat mirna targets; bars represent the percentage of targets enriched for a specific go term relative to the total number in the category.
a go annotation for conserved mirna targets; b go annotation for novel mirna targets
some of the target genes identified in silico were possibly regulated by mirnas showing a differential expression in leaves and/or roots between the two developmental stages.
a number of these genes were chosen for further analysis.
for instance, our target gene computational approach, identified that ttu-mir169c may regulate the homologue to hap <dig> subunit  of hap complex.
another mirna, ttu-mir156n, was predicted to regulate durum wheat genes encoding the squamosa promoter-binding-like proteins  spl <dig>  spl <dig>  spl <dig>  and spl <dig>  and a protein similar to the spl teosinte glume architecture  <dig>  of zea mays .
additionally, genes encoding for a homologue to the mikc-type mads-box transcription factor wm <dig> , a flowering promoting factor-like  <dig>  and the myb <dig>  transcription factor  were predicted to be the target genes of ttu-mir444d, ttu-mir <dig> and ttu-mir319f, respectively.
two genes encoding for the tira <dig> tir1-like auxin receptor  and a protein similar to the auxin signaling f-box  <dig>  of arabidopsis  were identified as target genes of mir393a.
finally, three genes encoding for a cytochrome p <dig> like-tbp , the oryza sativa early proembryo mrna, similar to rpe <dig> gene  and a photosystem ii p <dig> chlorophyll a apoprotein  were identified to be putative target genes of the ttu-novel- <dig>
expression of predicted target genes
eleven predicted target genes of seven  mirnas differentially expressed between the two developmental stages and possibly involved in plant development, were considered for expression using qpcr analysis.
prior to expression analyses, the target gene fragments used for qpcr experiments were sequenced to confirm the mirna binding site and the sequence of the predicted target genes.
the expression level of each gene was determined in roots and leaves of the late milk stage relative to the young stage.
in general, more interesting results were observed in roots compared to leaves.
the putative target gene squamosa spl <dig>  showed an evident strong downregulation  in roots of the adult plants compared to young roots, while ttu-mir156n was upregulated in the same tissue , evidencing a possible correlation between the expression of the mirna and this target gene.
the other putative target of ttu-mir156n, tga <dig>  showed a non-significant variation in the expression level between the two developmental stages.
on the other hand, both genes were highly upregulated in z <dig> leaves compared to z <dig> leaves, even though ttu-mir156n was not differentially expressed in this tissue in our conditions .fig.
<dig> several mirnas differentially expressed and their target genes  in leaves  and roots  at the mature stage compared to the early stage .
target genes are indicated by abbreviations: tga <dig>  teosinte glume architecture 1; afb <dig>  auxin signaling f-box  <dig> protein; fpf <dig>  flowering promoting factor-like 1; mmw, mikc-type mads-box transcription factor wm30; tira <dig>  tir1-like auxin receptor; cyp <dig>  cytochrome p <dig> like_tbp; ep, early proembryo; psbb, photosystem ii p <dig> chlorophyll a apoprotein; spl <dig>  squamosa promoter-binding-like protein 2; myb <dig>  transcription factor myb <dig>  expression levels are given as fold change of mature stage versus early stage.
data are means ± se.
p values between three early stage replicates and three mature stage replicates for each gene were calculated using unpaired, 2-tailed student’s t test
the relative expression level of myb <dig> , the putative target gene of ttu-mir319f, showed about two-fold increase  in z <dig> roots compared to the z <dig> stage.
this value was inversely correlated to the expression of ttu-mir319f, which showed an approximately two-fold decrease in the late milk stage.
in the leaves, myb <dig> did not vary significantly and ttu-mir319f showed an expression level, whose variation was not significant.
although mikc-type mads-box wm <dig>  tf was highly downregulated in mature leaves, its expression did not appear to be correlated to that of the mirna expected to regulate it, ttu-mir444d.
in the roots, mmw did not show a significant variation in gene expression, while ttu-mir444d was upregulated.
tira <dig>  tir1-like auxin receptor, and the auxin signaling f-box  <dig>  protein are putative targets of ttu-mir393a.
both in leaves and roots, the two genes were not differentially expressed at the late milk stage compared to the early developmental phase, even though a significantly upregulation of ttu-mir393a was observed in the leaves.
the flowering promoting factor-like  <dig>  was predicted to be the putative target of ttu-mir <dig>  however, the two expression analyses did not show any correlation .
the three genes identified as putative targets of the ttu-novel- <dig> were also analysed by qpcr.
in leaves, the photosystem ii p <dig> chlorophyll a apoprotein  was significantly upregulated in adult leaves , and its mirna was significantly downregulated, showing a correlation between the mirna and its putative target.
the other two predicted targets of ttu-novel- <dig>  cytochrome p <dig> like-tbp and early proembryo did not show differential expression in leaves.
in roots, the three genes, as well as the ttu-novel- <dig> were downregulated.
discussion
micrornas have proved to regulate a wide range of physiological processes in plants, including development, transition from one physiological phase to another, flowering time, etc.
<cit> .
in durum wheat, little is known on the mirna presence, abundance, and involvement in plant development.
using massive sequencing of srna libraries of the two different durum wheat cultivars, ciccio and svevo, we obtained over  <dig> million filtered sequences, which allowed us to retrieve  <dig> conserved and  <dig> potential novel mirnas.
in order to understand whether some of the identified mirnas could play a role in durum wheat plant development, especially in leaves and roots, the expression level of selected mirnas was evaluated in a young stage  and in an adult stage .
most of the conserved mirnas identified showed a differential expression between the early and the mature stage, and therefore they may play an important role in durum wheat development.
based on go annotation of target genes of durum wheat mirnas, various biological processes and molecular functions were predicted.
as in most plant species, a large number of durum wheat mirnas were expected to target transcripts annotated as nucleic acid binding proteins, several of which were homologs of known transcription factors controlling gene expression involved in plant development, morphology, and flowering time  <cit> .
in plants, gene regulation accomplished by mirnas can occur via target mrna degradation and/or translational repression.
the first mechanism usually arises when a strong sequence complementarity exists between the mirna and its target, while imperfect pairing causes translational repression without mrna cleavage  <cit> .
due to mirna regulation, a negative correlation is observed between the expression level of a mirna and its corresponding target.
for this reason, qpcr experiments measuring the expression level of mirnas and their targets are a valuable method to validate putative target genes .
in this work, we observed a negative correlation between the expression of some mirnas and their putative target genes, especially in three cases: ttu-mir156n and its target spl <dig>  ttu-mir319f and the gene myb <dig> in the roots, ttu-novel- <dig> and its putative target psbb in the leaves.
squamosa spl genes encode plant-specific tfs involved in the transition from juvenile to adult phase and from vegetative to reproductive phase  <cit> .
several spl genes  are post-transcriptionally regulated by mir <dig>  <cit> .
investigations on mir156/spls relationship published so far, have focused almost exclusively on their roles in plant aerial organs.
however, recently, the important role of mir156/spls in arabidopsis lateral root development has been unravelled  <cit> .
these mir156/spl feedback interactions throughout the development of plant aerial organs and roots may be conserved in other species  <cit> .
therefore, the upregulation of ttu-mir156n in the roots at the mature development stage and the downregulation of its putative target genes  in the same tissue, suggest a role of mir156/spls and mir156/tga <dig> in root development of durum wheat, even though the downregulation of tga <dig>  was not statistically significant.
in arabidopsis leaves, mir <dig> expression levels are high in the juvenile phase and decrease during shoot maturation  <cit> .
it has to be noticed that in our durum wheat adult leaves, both spl <dig> and tga <dig> are strongly and significantly upregulated, however, ttu-mir156n is similar to the control.
the higher expression of the spl genes at the mature stage could be explained because of the role that spl genes play in the control of lateral organs development in association with shoot maturation in the reproductive phase, as it has been demonstrated for spl <dig> in arabidopsis  <cit> .
further analysis should be carried out in order to understand the mechanisms of regulation of spl genes in adult leaves of durum wheat.
our findings suggest that ttu-mir319f might regulate the expression of myb <dig> in the roots of svevo plants.
homology analysis showed that there is a high similarity between wheat myb <dig> and arabidopsis atmyb <dig>  <cit> .
it has been demonstrated that atmyb <dig> expression increases in response to phytohormones, including jasmonic acid and ethylene  <cit> , which inhibit developmental processes such as cell growth.
therefore, the upregulation of ttu-mir <dig> at the immature stage and consequently the downregulation of its target gene myb <dig> may be a consequence of the higher cell growth, cell division and dna synthesis that characterizes an early developmental stage.
for ttu-novel- <dig> mirna, we predicted in silico three putative target genes encoding for the cytochrome p <dig> like_tbp, a protein similar to rpe <dig> , and the photosystem ii p <dig> chlorophyll a apoprotein .
we tested all these potential targets in qpcr and evaluated their expression level compared to that of ttu-novel- <dig> mirna.
the significant upregulation of the psbb gene in the mature leaves suggests that the activity of the photosystem ii p <dig> chlorophyll a apoprotein is higher in the mature leaves of durum wheat, compared to young leaves.
the lower expression of psbb in the leaves of the early stage could be a possible result of a higher target degradation by the ttu-novel- <dig> mirna in the young leaves.
as many processes of the plant, photosynthesis is targeted by a high number of mirnas, as it was demonstrated in rice  <cit> .
in monocots, there is a polarization in cell division and expansion being the oldest cells of the leaves photosynthetically more developed than the younger ones  <cit> .
this phenomenon might explain why in durum wheat mature leaves, a higher expression of the gene was observed.
our data suggest that a possible regulatory mechanism of psbb gene dependent on ttu-novel- <dig> might occur in durum wheat leaves, even though deeper analyses need to be conducted in order to confirm that psbb is the target of this durum wheat novel mirna.
ribosomal rna  is extremely abundant in the eukaryotic cells and constitutes about 80 % of the total cytoplasmatic rna  <cit> .
small rdna-derived rnas  have been identified in mouse and human and they are involved in the regulation of metabolism and other biological processes  <cit> .
high-throughput sequencing of small rna transcriptome reveals the existence of many different rna fragments derived from small rna.
in many cases, these sequences are discarded as degradation products.
nevertheless, experimental evidence suggests that rna fragments derived from small non coding rnas such as nucleolar rna  and transfer rna  are not just the result of rna turnover but rather functional molecules  <cit> .
to our knowledge, this is the first report on mirnas deriving from rrna in plants.
the expression level of the putative target genes of mir444d  and mir393a  was evaluated since a differential expression of these mirnas was observed in roots  and both leaves and roots .
however, the expression analysis of their target genes did not show a correlation with the mirna expression.
mir <dig> is specific to monocots and it targets four mikc-type mads-box genes in rice , homologous to arabidopsis anr <dig> .
in arabidopsis, it has been shown that anr <dig> is a major component in the no3−-signalling pathway that triggers lateral root growth.
in this work, ttu-mir <dig> was predicted to regulate genes encoding the mikc-type mads-box tfs , homologues of the t. aestivum mikc-type mads-box tfs.
ttu-mir444d clearly shows to be highly upregulated in roots of durum wheat in an advanced developmental stage, although no inverse correlation was observed with its putative target wm <dig>  probably, this mirna targets another member of the mikc-type mads-box tfs, or a post-transcriptional regulation might occur.
mir <dig> has been shown to play multiple roles in the rice no3−-signalling pathway involving root development, nitrate accumulation and pi starvation responses  <cit> .
overexpression of mir <dig> resulted in a reduction of rice lateral root elongation, but in a promotion of rice primary and adventitious root growth, in a nitrate-dependent manner.
therefore, it may be that in durum wheat, in the first developmental stages, plants are still more active in developing primary and adventitious roots than lateral roots.
mir <dig> regulates the activity of tir/afb <dig>  auxin receptor  family f-box proteins in arabidopsis leaves.
the perception and signalling of auxin, a hormone implicated in plant development, depends on these receptors  <cit> .
in our work, we detected a tira <dig> and an afb <dig> gene from durum wheat, by blasting our mirna data against t. aestivum and h. vulgare transcripts, respectively.
ttu-mir393a showed to be more expressed in the leaves at an advanced developmental stage compared to the early stage, so that taar genes are expected to be more active when durum wheat leaves are still in an immature stage.
however, in our analysis, no variation was detected in the expression levels of the predicted target genes.
this could be due to the presence of other isoforms of taar genes, which were not considered, or to the fact that notwithstanding perfect complementarity with mrna, in some cases mirna can inhibit translation irrespective of mrna cleavage  <cit> .
in the roots, we detected a very low mir <dig> expression.
this fact is in agreement with what was observed in other species, e.g.
artichoke where this mirna was not detected  <cit> .
in conclusion, some of the mirnas identified in durum wheat are differentially expressed in leaves and roots of the two developmental stages analyzed, indicating a regulation in the plant development.
for three mirnas, a negative correlation with their target was observed.
in other cases, a post-transcriptional regulation could occur.
interestingly, three novel mirnas were predicted to derive from rrna, and this could open new perspectives in the study of plant mirna.
CONCLUSIONS
the analysis of the durum wheat mirnaome from ciccio and svevo tissues at the late milk stage allowed to identify  <dig> conserved and  <dig> potential novel mirnas.
several mirnas identified were differentially expressed in the mature  developmental stage compared to young  tissues.
our results suggest that in roots, the putative genes encoding for the squamosa spl <dig> and tga <dig> proteins are regulated by ttu-mir156n, while myb <dig> transcription factor by ttu-mir319f.
additionally, the photosystem ii p <dig> chlorophyll a apoprotein gene showed an expression level negatively correlated to that of ttu-novel- <dig> in leaves.
these genes might play important roles in root/leaf development in durum wheat.
the prediction of novel mirnas putatively derived from ribosomal rna opens new perspectives on the study of plant mirnas.
abbreviations
arf, auxin response; asf <dig>  auxin signaling f-box 2; cdc, cell division control protein; edta, ethylenediaminetetraacetic acid; gapdh, glyceraldehyde-3-phosphate dehydrogenase; go, gene ontology; hvgi, hordeum vulgare unigenes from dfci gene index; mirna, microrna; psbb, photosystem ii p <dig> chlorophyll a apoprotein; qpcr, quantitative pcr; rrna, ribosomal rna; scl <dig>  scarecrow-like transcription factor 6; spl, squamosa promoter-binding-like proteins; tagi, triticum aestivum unigenes from dfci gene index; tf, transcription factor; tga <dig>  teosinte glume architecture 1; tir/afb <dig>  transport inhibitor response/auxin signalling f-box protein <dig>
additional files
additional file 1: table s <dig>  list of primers used for q-rt-pcr assay for durum wheat mirnas and their putative target genes.
additional file 2: figure s <dig>  length distribution of small rna sequences in two libraries obtained from durum wheat ciccio and svevo cultivars.
additional file 3: table s <dig>  conserved mirnas identified in durum wheat ciccio and svevo cultivars by illumina sequencing.
mirnas were annotated according to mirbase release  <dig>
additional file 4: table s <dig>  novel mirnas from durum wheat identified in durum wheat ciccio and svevo cultivars by illumina sequencing.
mirnas were predicted in silico using similarity search and secondary structure analyses of putative pre-mirnas sequences.
additional file 5: table s <dig>  putative target genes of durum wheat mirnas identified in silico.
the analysis was performed on three datasets: durum wheat ests from the ncbi database, t. aestivum unigenes from dfci gene index  version  <dig>  and hordeum vulgare unigenes from dfci gene index  version  <dig>
diana l. zuluaga and domenico de paola are both first authors.

 durum wheat is the most widespread crop in the mediterranean basin, since it is the base for the production of pasta and other typical products of this area such as breads and cuscus. durum wheat is an allotetraploid species which derives from the cross between t. urartu and an ancient relative of aegilops speltoides tausch.. in spite of the economic importance of durum wheat, there is still a lot to know about the genomics and transcriptomics of this species. in particular, information on the micrornaome of durum wheat is lacking. mirnas are small non-coding rna molecules which downregulate gene expression by mrna degradation or translational repression. several conserved mirnas regulate genes encoding transcription factors with important roles in a variety of processes in plants, including protein degradation, response to biotic and abiotic stresses, signal transduction, and developmental processes. mirnas are well studied in the reproductive parts of plants, where they regulate flowering time, seed and flower development. the key role of mirnas in developmental processes has been widely demonstrated in model plants and crops including some species of the poaceae family. in rice, mir regulates flowering time and temporal expression changes of numerous genes during leaf development, whereas mir and mir were detected to have an important role in rice seed development. in barley, mir potentially targets the auxin response factors and which have critical roles in plant root development. additionally, barley mir targets a plant-specific scarecrow-like transcription factor which regulates several plant developmental processes. in maize, it has been proved that mir downregulates glossy to promote vegetative phase change. the widespread use of high throughput sequencing has promoted advancements in the genomics of plants, including those species with large, repetitive and polyploid genomes as wheat. the recent availability of the hexaploid bread wheat draft genome has sped up the discovery of small rnas in this species in the past few years. thus, so far, mature mirnas are annotated for t. aestivum in the mirna repository mirbase. several studies have uncovered development-associated mirnas in bread wheat mainly in the seeds, thus being potentially involved in the regulation of grain-filling, and preferentially expressed in developing or germinating seeds. conversely, a limited number of studies have been carried out on mirnas in tetraploid wheats. for instance, mirna expression patterns in response to drought stress were assessed in leaf and root of wild emmer wheat, the ancestor of durum wheat, through a plant mirna microarray platform and a homology-based approach. additionally, the same group performed an in silico identification of several mirnas from root transcriptome of t. turgidum ssp. durum and t. turgidum ssp. dicoccoides. in spite of this, the discovery of mirnas in durum wheat is still largely unknown, with just one mirna reported in mirbase release the aim of this study was to identify conserved and novel mirnas in durum wheat, and evaluate their possible involvement in plant development. to this purpose, two small rna libraries from the cultivars ciccio and svevo were generated from a pool of different tissues at the late milk stage. following high throughput sequencing, conserved and novel mirnas and their target genes were identified. some mirnas were predicted to regulate genes involved in plant development or in important biochemical processes such as photosynthesis. interestingly, we found three novel mirnas derived from ribosomal rna in durum wheat. in order to investigate the possible involvement of selected conserved and novel mirnas and target genes in the development of leaves and roots, quantitative pcr was performed to highlight differences in expression profiles. this work provides useful information that may be employed in breeding programs aiming at improving durum wheat development. methods plant material and growth conditions seeds of two durum wheat cultivars, ciccio and svevo, were sterilized for 2 min in ethanol and then for 3 min in 2 % sodium hypochlorite. after 12 h in water aired with an oxygen air pump, seeds were transferred into mm petri dishes with a wet filter paper and kept for days at °c. once the seed hypocotyl became half centimetre, plantlets were transferred into cm diameter net pots containing a mix of 80 % agriperlite and 20 % expanded clay, and placed in a greenhouse at 18–20 °c under h photo-period. plants were grown in a hydroponic system containing mm mgso h2o, mm k2so mm kh2po 2mm ca h2o, μm 6mo7o h2o, 10μm h3bo 1μm mnso 1μm znso 5μm cuso 100μm fe-edta and a ph =  roots were submerged all the time in the solution and this was refreshed every two days. at the z stage of seedling growth, roots and leaves were collected from three individual plants from ciccio and svevo separately. at the late milk developmental stage, roots, leaves, flag leaves, and spikes were collected as above. all tissues were immediately frozen in liquid nitrogen and stored at −80 °c. small rna extraction and library preparation small rnas were isolated from mg of each tissue from each durum wheat plant at the stage z using the mirpremier microrna isolation kit. rna quality was verified by agarose gel according to manufacturer’s instructions. the small rna fractions from the four tissues of the three plants were pooled for library preparation. truseq small rna sample preparation protocol was used for the construction of two small rna libraries from ciccio and svevo tissues, respectively. libraries were sent for sequencing to the institute of applied genomics. data analysis and identification of conserved and novel mirna all reads were filtered and trimmed in order to discard low quality sequences and adapters. unique small rnas were counted eliminating redundancies and each count was normalized and expressed as transcript per million. in order to identify conserved mirnas, unique reads were subjected to similarity search using all plant mature mirnas from mirbase as reference database. to this purpose, blastn was used with parameters optimized for short sequences. only perfect matching srna were extracted and sequences with at least counts in one of the two libraries were considered for further analysis. potential novel mirnas were identified using a custom pipeline. the potential precursor sequences of novel durum wheat mirnas were obtained by blastn search against all the ests from the ncbi database. this similarity search was also extended to t. aestivum unigenes in the dfci gene index, version. all sequences identified were subjected to unafold to predict their secondary structure. to select conserved and novel mirnas in durum wheat we used the criteria for plant mirna annotation described by meyers et al.. qpcr quantification of micrornas for the qpcr validation of selected mirnas, svevo was chosen since it is one of the most studied durum wheat varieties, and is well-known for its qualitative properties, e.g. high protein content and high yellow index. small rnas were extracted using the mirpremier microrna isolation kit and mirnas were detected using stem-loop rt-pcr method. small rna quality was checked by electrophoresis through a 1 % agarose gel containing 1 % gel red. quantitative qpcr analyses were performed using the universal probe library probe assay specific for mature mirna expression, described by varkonyi-gasic and colleagues. each reaction was performed in a final volume of μl containing μl of rotor-gene probe pcr master mix 2x, μl of mirna-specific forward primer and the universal reverse primer, μl of upl, probe #, μl of nuclease-free water, and µl of the respective retrotranscribed mirna. primers for gene expression analysis were designed based on sequences obtained from the library analysis. pcr conditions were: 95 °c for min, followed by cycles of 95 °c for 5 s, 60 °c for 10 s and 72 °c for 1 s. reactions were performed using a rotor-gene machine. three technical replications of three different plants and a no-template control were included for each experiment. the housekeeping mirna was chosen based on expression data using qpcr of eight different putative reference mirnas previously described in wheat. data were analysed using genorme and normfinder software and mir showed to be the most stable mirna in young and mature tissues of svevo roots and leaves in our experimental conditions. relative quantification of each single mirna expression was analysed using the comparative ct method as described by livak and schmittgen. differentially expressed mirnas were identified as those showing > 2-fold change between the two developmental stages. after performing independent student’s t test, fold changes with a p-value <  were considered statistically significant. prediction of ttu-mirna target genes novel and conserved durum wheat mirnas were used as query in the online tool psrnatarget. the analysis was performed using three different transcript libraries for target search: durum wheat ests from the ncbi database, t. aestivum unigenes from dfci gene index version and hordeum vulgare unigenes from dfci gene index version blast2go was used to obtain go annotation of the wheat unigenes based on blastx hits against the ncbi nr database. go functional classification was obtained executing mapping and annotation steps in blastgo with default parameters. two graphs representing go terms categories were determined for all durum wheat mirna targets. expression analysis of target genes the expression of mirna putative target genes was analysed using the sybr green dye and the rotor-gene machine. primers for gene expression analysis are found in additional file total rna was isolated using the rneasy mini kit and reverse transcribed into cdna, using the script cdna synthesis kit according to the manufacturer’s instructions. before expression analysis, the same primers were used for sequencing of the target gene fragments including the mirna binding site. real-time qpcr experiments were carried out using ng of cdna and the iq sybr green supermix following the manufacturer’s protocol. pcr conditions were: 95 °c for min, followed by cycles of 95 °c for 10 s and 58 °c for 30 s. the reference gene was chosen by testing three putative reference genes, i.e. actin, glyceraldehyde-3-phosphate dehydrogenase and cell division control protein, aaa-superfamily of atpases, through a qpcr expression analysis in leaves and roots of svevo young and mature tissues. data were analysed using genorme and normfinder software suggesting cdc as the most stable gene in our experimental conditions. therefore, using cdc as the reference gene, relative quantification of each single gene was performed using the comparative ct method as described by livak and schmittgen. differentially expressed genes were identified as those showing > 2-fold change between the two developmental stages. after performing independent student’s t test, fold changes with a p-value <  were considered statistically significant.  sequencing data analysis and identification of durum wheat conserved and novel mirnas in order to determine the expression of the micrornaome in durum wheat, two small rna libraries were generated and sequenced from a mix of leaves, flag leaves, roots and spikes of ciccio and svevo plants, separately, at the developmental stage of late milk. conserved and novel mirnas and their target genes were identified using bioinformatics tools. after removing the low quality reads and adaptors, a total of, and, reads were collected from ciccio and svevo, respectively. the size distribution of the total clean reads in the two libraries shows a canonical trend, with the 21– nucleotide long srnas representing the majority of the total reads. the nucleotide peak was the most abundant one, as observed in most plant species. reads of the two libraries were searched against ncrnas in rfam to filter out the snornas, snrnas, trnas, and rrnas before annotation with conserved plant mirnas. as a result,, and, reads were discarded from the library of ciccio and svevo, respectively. in order to identify conserved mirnas in durum wheat, all the unique reads longer than or shorter than nucleotides were mapped to known plant mirnas in mirbase release a total of known mirnas were found in the two libraries, in ciccio and in svevo. out of the conserved mirnas identified, belonged to conserved mirna families, whereas were unassigned. all the remaining reads, which were not annotated with mirbase, were used to identify potential new mirnas using bioinformatics tools. based on the analysis of the predicted secondary structure of durum wheat sequences, potential mirnas were considered as novel mirna candidates. this analysis revealed three non-canonical putative mirnas. similarity search of their precursor sequences matched to ribosomal rna, which has never been reported as a source of plant mirnas before. the analysis of the secondary structures and their predicted thermodynamic stability meet all criteria for in silico mirna prediction. moreover, these novel mirnas were detected in independent ngs libraries, both the ones used in this study, and other durum wheat srna proprietary libraries. two non-canonical mirnas were mapped to the two opposite arms of the same sequence. the precursor sequences of the novel mirnas ttu-novel- and ttu-novel-42/ttu-novel- were mapped to 18s rrna and 28s rrna, respectively. in order to verify whether other mirnas present in mirbase could derive from rrna, all plant mirna precursors were blasted against the nucleotide database of viridiplantae in ncbi. at least one mirna, peu-mir from populus euphratica matched with 18s rna, showing 96 % identity with idesia polycarpa gene for 18s ribosomal rna, as the best hit.fig. pre-mirna secondary structure of ttu-novel- and ttu-novel-. 5′-end is marked within a circle expression analysis of selected micrornas several of the conserved mirnas found were reported to regulate development in a range of plant species. in order to highlight a possible involvement of these mirnas in durum wheat cv. svevo leaf and root development, we selected of those for the analysis of relative expression by quantitative pcr at an early stage and an advanced late milk developmental stage. additionally, five novel mirnas were selected based on their predicted target or according to their origin from rrna molecules. this set of novel mirnas was analysed from roots and leaves of svevo plants in the two different stages. in order to choose a suitable housekeeping mirna, based on the literature, eight mirnas were selected and evaluated as putative housekeeping genes in our experiment: ttu-mir ttu-mir ttu-mir ttu-mir ttu-mir ttu-mir ttu-mir and ttu-mir the qpcr expression profiles of these mirnas were analysed for their stability, and ttu-mir proved to be the most stable in our experimental conditions. therefore, ttu-mir was chosen from this group as the housekeeping mirna for our analyses. results of mirna expression analysis show that six conserved mirnas were significantly downregulated in leaves at the mature stage compared to the leaves at the early stage, with fold changes ranging from to. additionally, ttu-mir167h and ttu-mir393a were found to be significantly upregulated in leaves at the adult stage compared to the early stage, with fold changes of and, respectively. for ttu-mir319f, although a higher expression level was observed in the z leaves, this value was not significant. on the other hand, in the roots, five conserved mirnas were found to be significantly downregulated at the z stage compared to the z stage, with fold changes ranging from to. ttu-mir166k was also strongly downregulated in roots, but with a p-value of. in contrast, ttu-mir156n, ttu-mir167h and ttu-mir444d were found significantly up regulated in elder roots, with fold changes of,, and, respectively.fig. relative expression analysis of selected mirnas and target genes. conserved and novel t. turgidum mirnas in leaves and roots at the early stage and at the mature stage. expression levels are given as fold change of mature stage versus early stage. data are means ± se. p values between three early stage replicates and three mature stage replicates for each mirna were calculated using unpaired, 2-tailed student’s t test the analysis carried out for selected novel mirnas confirmed the expression in durum wheat plants of the two putative rrna-derived mirnas previously identified, ttu-novel- and ttu-novel- in details, quantitative assay revealed that in leaves, two novel mirnas, ttu-novel- ttu-novel- were significantly downregulated at the stage z compared to the early stage z with fold changes of and, respectively. in roots, ttu-novel- and ttu-novel- were strongly downregulated. ttu-novel- showed a downregulation with a fold change of and in leaves and roots, respectively, but these values were not significant. for ttu-novel- and ttu-novel- the expression levels observed at the two developmental stages were comparable for both the tissues analysed. prediction of mirna specific target genes in order to predict ttu-mirna target genes, durum wheat ests, bread wheat unigenes and barley unigenes were used as distinct datasets, and subsequently results were analysed separately. sequences carrying a region of high complementarity with the ttu-mirnas were identified in the three databases, which were analysed with different priority starting from t. turgidum followed by t. aestivum and h. vulgare. a number of target genes were predicted for all durum wheat mirnas identified in this work, for a total of and sequences, matching conserved and novel durum wheat mirnas, respectively. mirna target genes predicted using wheat dfci gene index were subjected to in silico functional analysis using gene ontology annotation. the assignment of go terms highlighted different biological processes and molecular functions for novel and conserved mirna targets. according to terms description relative to molecular function, the majority of the targets interact with nucleic acids. in fact, for both novel and conserved durum wheat mirna target genes, nucleotide binding or nucleoside phosphate binding were major classes with approximately 15 % of total sequences. with regard to the other go terms, there is a certain overlap of function in the comparison between novel and conserved mirna targets, with some significant differences such as cation binding including and  % targets for conserved and novel mirnas, respectively. with regard to the biological processes, the majority of the targets were annotated to terms related to metabolic processes. go annotations partially overlap for novel and conserved mirna targets, with the exception of the class “response to chemical”, which is missing in novel mirna targets.fig. bar charts representing biological processes and molecular function categorization according to go terms associated to durum wheat mirna targets; bars represent the percentage of targets enriched for a specific go term relative to the total number in the category. a go annotation for conserved mirna targets; b go annotation for novel mirna targets some of the target genes identified in silico were possibly regulated by mirnas showing a differential expression in leaves and/or roots between the two developmental stages. a number of these genes were chosen for further analysis. for instance, our target gene computational approach, identified that ttu-mir169c may regulate the homologue to hap subunit of hap complex. another mirna, ttu-mir156n, was predicted to regulate durum wheat genes encoding the squamosa promoter-binding-like proteins spl spl spl and spl and a protein similar to the spl teosinte glume architecture of zea mays. additionally, genes encoding for a homologue to the mikc-type mads-box transcription factor wm, a flowering promoting factor-like and the myb transcription factor were predicted to be the target genes of ttu-mir444d, ttu-mir and ttu-mir319f, respectively. two genes encoding for the tira tir1-like auxin receptor and a protein similar to the auxin signaling f-box of arabidopsis were identified as target genes of mir393a. finally, three genes encoding for a cytochrome p like-tbp, the oryza sativa early proembryo mrna, similar to rpe gene and a photosystem ii p chlorophyll a apoprotein were identified to be putative target genes of the ttu-novel- expression of predicted target genes eleven predicted target genes of seven mirnas differentially expressed between the two developmental stages and possibly involved in plant development, were considered for expression using qpcr analysis. prior to expression analyses, the target gene fragments used for qpcr experiments were sequenced to confirm the mirna binding site and the sequence of the predicted target genes. the expression level of each gene was determined in roots and leaves of the late milk stage relative to the young stage. in general, more interesting results were observed in roots compared to leaves. the putative target gene squamosa spl showed an evident strong downregulation in roots of the adult plants compared to young roots, while ttu-mir156n was upregulated in the same tissue, evidencing a possible correlation between the expression of the mirna and this target gene. the other putative target of ttu-mir156n, tga showed a non-significant variation in the expression level between the two developmental stages. on the other hand, both genes were highly upregulated in z leaves compared to z leaves, even though ttu-mir156n was not differentially expressed in this tissue in our conditions.fig. several mirnas differentially expressed and their target genes in leaves and roots at the mature stage compared to the early stage. target genes are indicated by abbreviations: tga teosinte glume architecture 1; afb auxin signaling f-box protein; fpf flowering promoting factor-like 1; mmw, mikc-type mads-box transcription factor wm30; tira tir1-like auxin receptor; cyp cytochrome p like_tbp; ep, early proembryo; psbb, photosystem ii p chlorophyll a apoprotein; spl squamosa promoter-binding-like protein 2; myb transcription factor myb expression levels are given as fold change of mature stage versus early stage. data are means ± se. p values between three early stage replicates and three mature stage replicates for each gene were calculated using unpaired, 2-tailed student’s t test the relative expression level of myb, the putative target gene of ttu-mir319f, showed about two-fold increase in z roots compared to the z stage. this value was inversely correlated to the expression of ttu-mir319f, which showed an approximately two-fold decrease in the late milk stage. in the leaves, myb did not vary significantly and ttu-mir319f showed an expression level, whose variation was not significant. although mikc-type mads-box wm tf was highly downregulated in mature leaves, its expression did not appear to be correlated to that of the mirna expected to regulate it, ttu-mir444d. in the roots, mmw did not show a significant variation in gene expression, while ttu-mir444d was upregulated. tira tir1-like auxin receptor, and the auxin signaling f-box protein are putative targets of ttu-mir393a. both in leaves and roots, the two genes were not differentially expressed at the late milk stage compared to the early developmental phase, even though a significantly upregulation of ttu-mir393a was observed in the leaves. the flowering promoting factor-like was predicted to be the putative target of ttu-mir however, the two expression analyses did not show any correlation. the three genes identified as putative targets of the ttu-novel- were also analysed by qpcr. in leaves, the photosystem ii p chlorophyll a apoprotein was significantly upregulated in adult leaves, and its mirna was significantly downregulated, showing a correlation between the mirna and its putative target. the other two predicted targets of ttu-novel- cytochrome p like-tbp and early proembryo did not show differential expression in leaves. in roots, the three genes, as well as the ttu-novel- were downregulated. discussion micrornas have proved to regulate a wide range of physiological processes in plants, including development, transition from one physiological phase to another, flowering time, etc.. in durum wheat, little is known on the mirna presence, abundance, and involvement in plant development. using massive sequencing of srna libraries of the two different durum wheat cultivars, ciccio and svevo, we obtained over million filtered sequences, which allowed us to retrieve conserved and potential novel mirnas. in order to understand whether some of the identified mirnas could play a role in durum wheat plant development, especially in leaves and roots, the expression level of selected mirnas was evaluated in a young stage and in an adult stage. most of the conserved mirnas identified showed a differential expression between the early and the mature stage, and therefore they may play an important role in durum wheat development. based on go annotation of target genes of durum wheat mirnas, various biological processes and molecular functions were predicted. as in most plant species, a large number of durum wheat mirnas were expected to target transcripts annotated as nucleic acid binding proteins, several of which were homologs of known transcription factors controlling gene expression involved in plant development, morphology, and flowering time. in plants, gene regulation accomplished by mirnas can occur via target mrna degradation and/or translational repression. the first mechanism usually arises when a strong sequence complementarity exists between the mirna and its target, while imperfect pairing causes translational repression without mrna cleavage. due to mirna regulation, a negative correlation is observed between the expression level of a mirna and its corresponding target. for this reason, qpcr experiments measuring the expression level of mirnas and their targets are a valuable method to validate putative target genes. in this work, we observed a negative correlation between the expression of some mirnas and their putative target genes, especially in three cases: ttu-mir156n and its target spl ttu-mir319f and the gene myb in the roots, ttu-novel- and its putative target psbb in the leaves. squamosa spl genes encode plant-specific tfs involved in the transition from juvenile to adult phase and from vegetative to reproductive phase. several spl genes are post-transcriptionally regulated by mir. investigations on mir156/spls relationship published so far, have focused almost exclusively on their roles in plant aerial organs. however, recently, the important role of mir156/spls in arabidopsis lateral root development has been unravelled. these mir156/spl feedback interactions throughout the development of plant aerial organs and roots may be conserved in other species. therefore, the upregulation of ttu-mir156n in the roots at the mature development stage and the downregulation of its putative target genes in the same tissue, suggest a role of mir156/spls and mir156/tga in root development of durum wheat, even though the downregulation of tga was not statistically significant. in arabidopsis leaves, mir expression levels are high in the juvenile phase and decrease during shoot maturation. it has to be noticed that in our durum wheat adult leaves, both spl and tga are strongly and significantly upregulated, however, ttu-mir156n is similar to the control. the higher expression of the spl genes at the mature stage could be explained because of the role that spl genes play in the control of lateral organs development in association with shoot maturation in the reproductive phase, as it has been demonstrated for spl in arabidopsis. further analysis should be carried out in order to understand the mechanisms of regulation of spl genes in adult leaves of durum wheat. our findings suggest that ttu-mir319f might regulate the expression of myb in the roots of svevo plants. homology analysis showed that there is a high similarity between wheat myb and arabidopsis atmyb. it has been demonstrated that atmyb expression increases in response to phytohormones, including jasmonic acid and ethylene, which inhibit developmental processes such as cell growth. therefore, the upregulation of ttu-mir at the immature stage and consequently the downregulation of its target gene myb may be a consequence of the higher cell growth, cell division and dna synthesis that characterizes an early developmental stage. for ttu-novel- mirna, we predicted in silico three putative target genes encoding for the cytochrome p like_tbp, a protein similar to rpe, and the photosystem ii p chlorophyll a apoprotein. we tested all these potential targets in qpcr and evaluated their expression level compared to that of ttu-novel- mirna. the significant upregulation of the psbb gene in the mature leaves suggests that the activity of the photosystem ii p chlorophyll a apoprotein is higher in the mature leaves of durum wheat, compared to young leaves. the lower expression of psbb in the leaves of the early stage could be a possible result of a higher target degradation by the ttu-novel- mirna in the young leaves. as many processes of the plant, photosynthesis is targeted by a high number of mirnas, as it was demonstrated in rice. in monocots, there is a polarization in cell division and expansion being the oldest cells of the leaves photosynthetically more developed than the younger ones. this phenomenon might explain why in durum wheat mature leaves, a higher expression of the gene was observed. our data suggest that a possible regulatory mechanism of psbb gene dependent on ttu-novel- might occur in durum wheat leaves, even though deeper analyses need to be conducted in order to confirm that psbb is the target of this durum wheat novel mirna. ribosomal rna is extremely abundant in the eukaryotic cells and constitutes about 80 % of the total cytoplasmatic rna. small rdna-derived rnas have been identified in mouse and human and they are involved in the regulation of metabolism and other biological processes. high-throughput sequencing of small rna transcriptome reveals the existence of many different rna fragments derived from small rna. in many cases, these sequences are discarded as degradation products. nevertheless, experimental evidence suggests that rna fragments derived from small non coding rnas such as nucleolar rna and transfer rna are not just the result of rna turnover but rather functional molecules. to our knowledge, this is the first report on mirnas deriving from rrna in plants. the expression level of the putative target genes of mir444d and mir393a was evaluated since a differential expression of these mirnas was observed in roots and both leaves and roots. however, the expression analysis of their target genes did not show a correlation with the mirna expression. mir is specific to monocots and it targets four mikc-type mads-box genes in rice, homologous to arabidopsis anr. in arabidopsis, it has been shown that anr is a major component in the no3−-signalling pathway that triggers lateral root growth. in this work, ttu-mir was predicted to regulate genes encoding the mikc-type mads-box tfs, homologues of the t. aestivum mikc-type mads-box tfs. ttu-mir444d clearly shows to be highly upregulated in roots of durum wheat in an advanced developmental stage, although no inverse correlation was observed with its putative target wm probably, this mirna targets another member of the mikc-type mads-box tfs, or a post-transcriptional regulation might occur. mir has been shown to play multiple roles in the rice no3−-signalling pathway involving root development, nitrate accumulation and pi starvation responses. overexpression of mir resulted in a reduction of rice lateral root elongation, but in a promotion of rice primary and adventitious root growth, in a nitrate-dependent manner. therefore, it may be that in durum wheat, in the first developmental stages, plants are still more active in developing primary and adventitious roots than lateral roots. mir regulates the activity of tir/afb auxin receptor family f-box proteins in arabidopsis leaves. the perception and signalling of auxin, a hormone implicated in plant development, depends on these receptors. in our work, we detected a tira and an afb gene from durum wheat, by blasting our mirna data against t. aestivum and h. vulgare transcripts, respectively. ttu-mir393a showed to be more expressed in the leaves at an advanced developmental stage compared to the early stage, so that taar genes are expected to be more active when durum wheat leaves are still in an immature stage. however, in our analysis, no variation was detected in the expression levels of the predicted target genes. this could be due to the presence of other isoforms of taar genes, which were not considered, or to the fact that notwithstanding perfect complementarity with mrna, in some cases mirna can inhibit translation irrespective of mrna cleavage. in the roots, we detected a very low mir expression. this fact is in agreement with what was observed in other species, e.g. artichoke where this mirna was not detected. in conclusion, some of the mirnas identified in durum wheat are differentially expressed in leaves and roots of the two developmental stages analyzed, indicating a regulation in the plant development. for three mirnas, a negative correlation with their target was observed. in other cases, a post-transcriptional regulation could occur. interestingly, three novel mirnas were predicted to derive from rrna, and this could open new perspectives in the study of plant mirna.  the analysis of the durum wheat mirnaome from ciccio and svevo tissues at the late milk stage allowed to identify conserved and potential novel mirnas. several mirnas identified were differentially expressed in the mature developmental stage compared to young tissues. our results suggest that in roots, the putative genes encoding for the squamosa spl and tga proteins are regulated by ttu-mir156n, while myb transcription factor by ttu-mir319f. additionally, the photosystem ii p chlorophyll a apoprotein gene showed an expression level negatively correlated to that of ttu-novel- in leaves. these genes might play important roles in root/leaf development in durum wheat. the prediction of novel mirnas putatively derived from ribosomal rna opens new perspectives on the study of plant mirnas. abbreviations arf, auxin response; asf auxin signaling f-box 2; cdc, cell division control protein; edta, ethylenediaminetetraacetic acid; gapdh, glyceraldehyde-3-phosphate dehydrogenase; go, gene ontology; hvgi, hordeum vulgare unigenes from dfci gene index; mirna, microrna; psbb, photosystem ii p chlorophyll a apoprotein; qpcr, quantitative pcr; rrna, ribosomal rna; scl scarecrow-like transcription factor 6; spl, squamosa promoter-binding-like proteins; tagi, triticum aestivum unigenes from dfci gene index; tf, transcription factor; tga teosinte glume architecture 1; tir/afb transport inhibitor response/auxin signalling f-box protein additional files additional file 1: table s list of primers used for q-rt-pcr assay for durum wheat mirnas and their putative target genes. additional file 2: figure s length distribution of small rna sequences in two libraries obtained from durum wheat ciccio and svevo cultivars. additional file 3: table s conserved mirnas identified in durum wheat ciccio and svevo cultivars by illumina sequencing. mirnas were annotated according to mirbase release additional file 4: table s novel mirnas from durum wheat identified in durum wheat ciccio and svevo cultivars by illumina sequencing. mirnas were predicted in silico using similarity search and secondary structure analyses of putative pre-mirnas sequences. additional file 5: table s putative target genes of durum wheat mirnas identified in silico. the analysis was performed on three datasets: durum wheat ests from the ncbi database, t. aestivum unigenes from dfci gene index version and hordeum vulgare unigenes from dfci gene index version diana l. zuluaga and domenico de paola are both first authors. 
in the present work, two small rna libraries from durum wheat ciccio and svevo cultivars were generated from different tissues at the late milk developmental stage. our results suggest that, in durum wheat, these genes may play important roles in root/leaf development and are subjected to mirna regulation. putative target genes were predicted for conserved and novel mirnas, the majority of which interact with nucleic acids, according to go terms relative to molecular function. the allotetraploid durum wheat is a highly economically important species especially in the mediterranean basin. a total of conserved and potential novel mirnas were identified in the two libraries and interestingly, three novel mirnas were found to be derived from ribosomal rna. the prediction of novel mirnas putatively derived from ribosomal rna opens new perspectives on the study of plant mirnas. moreover, target gene expression analysis suggested that in roots, the putative genes encoding for the squamosa spl and tga proteins are regulated by ttu-mir156n, while myb transcription factor by ttu-mir319f. additionally, the photosystem ii p chlorophyll a apoprotein gene showed an expression level negatively correlated to that of ttu-novel- in leaves. micrornasdurum wheatsequencingbioinformatics analysisqpcr expression analysistarget genesrna derived mirnashttp://dx.doi.org/ /501100003407ministero dell’istruzione, dell’università e della ricercapon01_ iscocem “sviluppo tecnologico e innovazione per la sostenibilità e competitività della cerealicoltura meridionale”prin 2010- dd23/10/ n. “identificazione e caratterizzazione di geni utili ad incrementare la produttività e sostenibilità del frumento duro”sonnante gabriella issue-copyright-statement© the author 2016
 the allotetraploid durum wheat is a highly economically important species especially in the mediterranean basin. however, its genomics, transcriptomics and in particular micrornaome are still largely unknown.  in the present work, two small rna libraries from durum wheat ciccio and svevo cultivars were generated from different tissues at the late milk developmental stage. a total of conserved and potential novel mirnas were identified in the two libraries and interestingly, three novel mirnas were found to be derived from ribosomal rna. putative target genes were predicted for conserved and novel mirnas, the majority of which interact with nucleic acids, according to go terms relative to molecular function. quantitative qpcr analysis showed that several mirnas identified were differentially expressed in the mature developmental stage compared to young tissues. moreover, target gene expression analysis suggested that in roots, the putative genes encoding for the squamosa spl and tga proteins are regulated by ttu-mir156n, while myb transcription factor by ttu-mir319f. additionally, the photosystem ii p chlorophyll a apoprotein gene showed an expression level negatively correlated to that of ttu-novel- in leaves.  our results suggest that, in durum wheat, these genes may play important roles in root/leaf development and are subjected to mirna regulation. the prediction of novel mirnas putatively derived from ribosomal rna opens new perspectives on the study of plant mirnas. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords micrornasdurum wheatsequencingbioinformatics analysisqpcr expression analysistarget genesrna derived mirnashttp://dx.doi.org/ /501100003407ministero dell’istruzione, dell’università e della ricercapon01_ iscocem “sviluppo tecnologico e innovazione per la sostenibilità e competitività della cerealicoltura meridionale”prin 2010- dd23/10/ n. “identificazione e caratterizzazione di geni utili ad incrementare la produttività e sostenibilità del frumento duro”sonnante gabriella issue-copyright-statement© the author 2016 