BACKGROUND
micrornas  are a class of small non-protein-coding rnas generated from single-stranded precursors with unique hairpin structures.
they regulate the expression of mrnas by targeting transcripts for cleavage or translational repression  <cit> .
the mirnas were initially isolated in caenorhabditis elegans as developmental timing regulators  <cit> .
since then, they have been found in a broad range of plants, as well as viruses and mammals.
plants mirnas were first identified in arabidopsis by different research groups  <cit> .
the biogenesis of plant mirnas is a complex multi-step enzymatic process  <cit> .
the mirnas are initially transcribed by rna polymerase ii in the cell nucleus as long primary mirnas that are cleaved into the mirna:mirna* duplexes by dicer-like  <dig> .
export of the duplexes into the cell cytoplasm is mediated by hasty.
after methyl groups are added to the 3' ends of the duplexes catalyzed by hen <dig>  one strand of the duplexes is selectively incorporated into the rna-induced silencing complex  to form the mature mirnas, whereas the other strand, designated the mirna*, is typically degraded.
dcl <dig> has also been shown to play a role in the biogenesis of a few mirnas with long hairpin precursors  <cit> .
plant mirnas exhibit high complementarity to their targets and can direct risc-mediated cleavage of target mrnas.
thus, it is widely accepted that plant mirnas induce post-transcriptional gene silencing predominantly through guiding mrna cleavage.
recent studies, however, have showed that translational inhibition is another kind of action mechanism for mirnas in plants  <cit> .
although initial studies have largely demonstrated the role of plant mirnas in development and morphogenesis processes, there are increasing number of reports indicating that plant mirnas also target genes involved in biotic as well as abiotic stress responses  <cit> .
low temperature is one of the most important environmental stimuli that affect plant growth and development.
recently, more and more reports demonstrate the role of plant mirnas in cold stress response.
the first study in this area, performed by sunkar and zhu, revealed the induced expression of mir <dig>  mir397b and mir <dig> in response to cold stress as well as other kinds of stress treatment.
in this study, mir319c appeared to be up-regulated by cold but not by dehydration, nacl, or aba treatment  <cit> .
zhou et al.
identified cold up-regulated mirnas using a computational, transcriptome-based approach  <cit> .
liu et al.
and lu et al.
used microarray analysis to identify cold-responsive mirnas in arabidopsis and populus, respectively  <cit> .
these studies indicate that the expression of several mirnas is affected by the cold treatment.
despite these efforts, our knowledge of the role played by mirnas in plant cold stress response is still limited at the whole-genome level.
the most challenging problem in understanding plant mirnas is to identify more novel mirnas.
three major approaches have been used for mirna discovery in plants: forward genetics, bioinformatic prediction as well as direct cloning and sequencing.
only a few mirnas were identified by forward genetic studies  <cit>  and predicting species-specific mirnas using bioinformatics method is difficult.
thus, direct cloning and sequencing is the most effective method for plant mirna discovery.
only a few hundred mirnas have been identified with this approach, which leads to a premature conclusion that the types of mirnas in plants are limited.
recently, the deep sequencing approach appeared and allowed the identification of numerous small rnas  <cit> .
it not only revealed a lot more species-specific mirnas, but also provided a picture of the genomic landscape of small rnas.
although significant progress has been made in identifying plant mirnas and understanding their action mechanism, the discovery of novel mirnas in plants on a genome-wide scale is still at the preliminary stage.
one limiting factor in mirna discovery is the availability of the whole genome sequence, only with which a comprehensive analysis of the potential hairpin precursor structures of cloned small rnas can be performed to distinguish mirnas from other kinds of small rnas.
thus, most of the studies have been done in arabidopsis, rice and populus, whose whole genome sequences are known.
many evolutionarily or economically important species have not been examined yet.
to further understand the function of plant mirnas, more effort should be directed toward plant species with specific developmental features, which may contain mirnas that are specific for these features  <cit> .
several economically important winter-habit monocots, such as winter wheat, barley, oat, rye, as well as cool-season biofuel and forage grasses, can cold-acclimate and acquire high tolerance to low temperature.
although it has been shown that mirnas are involved in plant cold stress response, little work has been done for monocotyledonous plant, especially for these winter-habit monocots, probably because rice, the model plant for monocots, is a tropical plant and incapable of cold acclimation.
brachypodium distachyon has emerged as a new monocot model plant, especially for temperate cereals and related grasses  <cit> .
brachypodium not only has a closer evolutionary relationship with cool-season temperate cereals and grasses than rice, but also possesses growth and developmental features that are common to these plants.
as a widely distributed winter-habit temperate plant, it has vernalization requirement and is capable of cold acclimation  <cit> .
recently, the draft genome sequence of brachypodium has been released , which makes this plant a good model for performing whole-genome-wide study of mirnas involved in cold response of winter-habit temperate cereals and grasses.
here we sequenced small rna populations in brachypodium with and without cold stress treatment using solexa, the high-throughput sequencing technology.
our studies not only identified species-specific mirnas for the novel monocotyledonous model plant, but also provided useful information for cold-responsive mirnas in temperate cereals and grasses.
RESULTS
deep sequencing of brachypodium small rnas
two small rna libraries, with  and without  cold-treatment, were generated using pooled rna isolated from the aerial parts of brachypodium seedlings.
sequencing of the brachypodium small rna libraries was performed with solexa, a high throughput sequencing technology producing highly accurate, reproducible and quantitative readouts of small rnas  <cit> .
it can be used as a tool for mirna expression profiling  <cit> .
for direct comparison, brachypodium seedlings used for nc and wc library construction were grown under similar conditions except for the cold treatment.
the same amount of rna was used to construct these two libraries and the samples were prepared in a similar manner.
sequencing of these two libraries was performed on the illumina's solexa sequencer and the samples were run side by side.
solexa sequencing of nc and wc libraries generated a total of  <dig> , <dig> and  <dig> , <dig> raw reads, respectively.
analysis of these reads resulted in identification of  <dig> , <dig> and  <dig> , <dig> sequences ranging from  <dig> to  <dig> nucleotides , respectively .
the remaining sequences were either of low quality  or smaller than  <dig> nt, and were excluded from further analysis.
the majority of the obtained small rna sequences from the nc library were 20- <dig> nt in size, which is the typical size range for dicer-derived products .
these small rnas were used for further analysis.
the 20- <dig> nt sequences from the nc and wc libraries were aligned to the draft brachypodium genome using soap  <cit> , and a total of  <dig> , <dig>  and  <dig> , <dig>  sequences were found to match the genome perfectly.
after removing the sequences corresponding to known rrnas, trnas, small nuclear rnas and small nucleolar rnas, a total of  <dig> , <dig>  and  <dig> , <dig>  sequences were obtained .
the wc library contained more non-coding rnas than the nc library, which reduced the population size of non-coding-rna-removed genome-matched sequences in the wc library to half of that in the nc library.
for direct comparison, the nc library was normalized to the wc library according to the population size mentioned above.
in both nc and wc libraries, most of the known mirnas and predicted mirnas were detected at least  <dig> times, enabling the comparison of mirna expression levels between these two libraries.
identification of novel monocot-specific and brachypodium-specific mirnas
one important feature that distinguishes mirnas from other endogenous small rnas is that their surrounding sequences can adopt a hairpin structure  <cit> .
the raw whole genome sequence data of brachypodium distachyon has been released by the us department of energy joint genome institute ;  and gene annotation was accomplished with fgenesh.
in this study, the annotated brachypodium genome sequences were used to predict hairpin structures based on mirna surrounding sequences.
to identify mirnas in brachypodium, the following strategy was adopted.
first, the cloned sequences were aligned with the brachypodium genome sequence using soap  <cit>  to search for perfectly matched sequences.
in the second step, sequences corresponding to known non-coding rnas  were discarded by performing the blastn search  <cit>  against the genbank  and rfam databases  <cit> .
the remaining sequences were used for the fold-back structure prediction and classified as either mirna candidates or other kinds of endogenous small rnas.
our analysis revealed that in the nc library  <dig> , <dig> sequences can be used for fold-back structure prediction and those that fulfilled the hairpin structure criteria described by jones-rhoades et al.
<cit>  were selected as candidate mirnas .
among these candidates, some also matched the plastid or mitochondrial genomes.
these sequences were not considered as mirnas in this study.
after these analyses,  <dig> unique sequences were obtained as mirna candidates .
according to the currently available gene annotation of the brachypodium genome ,  <dig> of these mirna candidates were located in the exons of potentially protein-coding genes, which were not further analyzed.
among the remaining  <dig> brachypodium mirna candidates,  <dig> of them have really low sequencing reads and were not detected in the wc library.
those candidates were excluded from further research according to the updated plant mirna annotation criteria  <cit> .
the remaining  <dig> putative brachypodium mirnas were carefully examined to be sure that they are qualified for the updated plant mirna annotation criteria  <cit> .
the lengths of the  <dig> predicted brachypodium mirnas vary from  <dig> to  <dig> nt, and more than half of them begin with 5' uridine , which is a characteristic feature of mirnas  <cit> .
the minimum free energies of these mirna precursors range from - <dig>  to - <dig>  kcal mol- <dig>  and the average value is about - <dig>  kcal mol- <dig> , which is similar to the free energy values of mirna precursors in rice and arabidopsis  <cit> .
these values are much lower than the folding free energies of trnas  or rrnas   <cit> .
the predicted hairpin structures for the mirna precursors are in the range of 71- <dig> nt , which is similar to what was observed in rice  <cit> .
none of these putative mirnas has been found in arabidopsis or rice before.
identification of conserved mirnas in brachypodium
to identify the conserved mirna homologs in brachypodium, small rna libraries were analyzed for the presence of known mirnas.
blastn search  <cit>  with an e-value cutoff of  <dig> was employed to search for small rnas with predicted hairpin structures against the central mirna registry database .
with this similarity search,  <dig> unique sequences were identified as known mirnas.
they showed higher sequence similarity to their homologs in rice than in arabidopsis and populus .
all of these conserved mirna precursors were identified in the draft brachypodium genome and they can adopt hairpin structures .
based on sequence similarity, these known mirnas could be grouped into  <dig> families .
four conserved mirna families were represented by more than one sequence variant in the small rna libraries.
mir <dig>  mir <dig> and mir <dig> were represented by two variants and mir <dig> was represented by four variants in the library .
some of these mirnas including mir <dig>  mir <dig>  mir <dig> and mi <dig> are deeply conserved, even in lower plants such as physcomitrella patens  <cit> .
mir <dig> and mir <dig>  which are species-specific mirnas, were also found in the nc library.
according to mirbase  and published papers, mir <dig> only appears in oryza sativa and mir <dig> has been identified in eschscholzia californica  <cit> , oryza sativa, physcomitrella patens as well as some pinus species  <cit> .
the mirna* sequences were detected for some bdi-mirnas .
the isolation frequencies for most of these mirna* sequences were about one-tenth of those for mature mirnas, consistent with the conclusion of rajagopalan et al.
<cit> .
for mir <dig>  there were  <dig> reads of the mirna and  <dig> reads of the mirna* in the nc library.
similarly, in the wc library mir <dig> had  <dig> reads of the mirna and  <dig> reads of the mirna*.
the mir529* was reported to be sequenced more often than the mir <dig> in rice and the annotated mir529*, therefore, is suspected to be the true mirna  <cit> .
our data indicate that mir <dig> is the true mirna and the existing annotation should not be changed.
it is also interesting to note that the mir169f was sequenced  <dig> times in the nc library and its potential mirna* was sequenced  <dig> times, which is almost equal to the mature mirna sequencing frequency.
name, the name of conserved brachypodium mirnas and underlines denote mirnas whose mirna* have been detected in solexa sequencing analysis; sequence, mirna sequence cloned in the nc library and; l, the length of mirna; arm/nt, the mirna location in the predicted hairpin structure /the length of the precursor; n, the number of loci in the genome for different mirna families; location, mirna location in the genome; ara, arabidopsis; pop, populus; ric, rice; whe, wheat; for plus symbols in the table: ++++, mirna sequences of brachypodium were identical to those in other plant species; +++, mirna sequences of brachypodium were conserved in other plant species but have variations at  <dig> nucleotide positions; ++, mirna sequences of brachypodium were conserved in other plant species but have variations at  <dig> nucleotide positions; +, mirna sequences of brachypodium were conserved in other plant species but have variations at ≥  <dig> nucleotide positions.
for mirnas whose homologs are not identified in rice, sequence comparison was only performed for the other four plant species.
the sequencing frequencies for mirnas in the library could be used as an index for estimating the relative abundance of mirnas.
solexa sequencing produced a large number of mirna sequences, allowing us to determine the relative abundance of mirnas in brachypodium.
the frequencies of the mirnas varied from  <dig>  to  <dig>  , indicating that the expression of mirnas varied greatly in brachypodium .
the mir <dig> was the most abundant mirna in our two sequencing datasets, accounting for about 20% of all sequence reads matching with the brachypodium genome and about 60% of the sequence reads of the conserved mirnas .
mir <dig> has also been the most abundant mirna in previous high-throughput sequencing studies in rice  <cit> , but not in arabidopsis .
according to the sequencing frequencies, some mirnas  were highly expressed in brachypodium, and were sequenced more than  <dig> times each.
the mir160a, mir164d, mir169f, mir172a, mir172b, mir <dig>  mir <dig>  mir <dig>  mir <dig>  mir395a, mir397a, mir <dig> and mir <dig> were moderately expressed, and were represented by the number of sequences varying between  <dig> and  <dig>  the third category, whose expression level was low , was represented by mir <dig>  mir169a, mir169c, mir169e and mir <dig> .
sequence analysis indicated that the relative abundance of certain members within a mirna family varied greatly in brachypodium.
for instance, the sequencing frequency for mir166a was  <dig> , whereas mir166g's sequencing frequency was only  <dig>  .
these results indicate that in one mirna family different members show clearly different expression levels, probably because their expression is tissue- or developmental-specific.
genomic organization of the mirnas
among the  <dig> conserved bachypodium mirnas,  <dig> were encoded by a single copy in the brachypodium genome, whereas the other  <dig> mirnas had multiple loci , which probably resulted from duplication events that were still active in the brachypodium genome.
for these  <dig> mirnas, most of them only had 2- <dig> loci in the genome, only a few had more than  <dig> loci.
for the predicted brachypodium mirnas, most of them had one single locus in the genome and only a few were represented by more than  <dig> loci  .
according to the currently available gene annotation for the brachypodium genome, most of these predicted mirnas were mapped to the intergenic regions.
interestingly, for almost all the conserved mirna families  with multiple members, their family sizes in brachypodium are much smaller than those in rice and populus, and most of them have sizes similar to those in arabidopsis .
for example, the number of family members for mir <dig>  mir <dig>  mir <dig>  mir <dig>  mir <dig>  mir <dig> and mir <dig> in brachypodium was similar to that in arabidopsis, but much higher in rice and populus .
for mir <dig>  its size in brachypodium was approximately half of that in rice .
name, the name of predicted cold-responsive brachypodium mirnas and underlines denote mirnas whose mirna* have been detected in solexa sequencing analysis; sequence, sequence cloned in small rna libraries; l, the length of mirna; arm/nt, the mirna location in the predicted hairpin structure /the length of the precursor; n, number of loci in the genome; location, mirna location in the genome; amfe, the adjusted minimum free energy  representing the mfe of  <dig> nucleotides, it is calculated by × 100; cd, the effect of cold-treatment on mirna expression, nc/wc, normalized sequencing frequencies in the nc and wc library.
the number of mirna family members in rice, arabidopsis and populus was obtained from the mirbase .
mir <dig> belongs to conserved mirnas that have been found across different plant species.
both mir <dig> and mir395* could be detected in the nc library, although their reads were not high.
according to the detected mirna and mirna* sequences as well as results of blastn search based on sequences similarity,  <dig> precursors with reasonable minimum free energy value were identified for mir <dig> in the draft brachypodium genome .
they were named mir395a~mir395n according to the order of their genomic location.
additionally, two other sequences in the surrounding region showed high similarity to mir <dig> precursors, but they contained no mature mir <dig> sequence and were therefore not considered for further analysis in this study .
for all the mir <dig> precursors, the mature mirnas are all located on the 3' arms.
the precursors of mir395a~c are clustered together and so were the precursors of mir395e~m.
.
it is interesting to note that mir395d in brachypodium is not clustered with other family members  although the members of this family always form clusters in other plant species  <cit> .
no mir395- precursor-like sequence was detected in its surrounding region.
on the basis of size and sequence similarity, two types of mir <dig> precursors exist in brachypodium.
type a includes mir395a~c and mir395f~n, and type b includes mir395d and mir395e.
both kinds of precursors show high sequence similarity to their homologs in rice .
the precursors of mir395a~c and mir395f~n show high sequence similarity, suggesting that they are derived from a series of gene duplication events.
these precursors exhibit different levels of sequence similarity to one another, implying that the duplication events happened at different time points during the evolution history of the mir <dig> family.
the precursors of mir395i and mir395k are nearly identical, implying a recent origin .
the impact of cold stress on conserved and species-specific mirnas in brachypodium
to detect the effect of cold stress on brachypodium mirnas, the expression of mirnas in brachypodium seedlings with  and without cold-treatment  were examined with solexa technology.
these two libraries were compared for their size distribution of small rnas.
in the nc library, the 21- <dig> nt rnas represented the predominant species , with 21-nt and 24-nt rnas being the two most abundant classes , consistent with the distribution patterns of small rnas from other plant species  <cit> .
however, in the wc library, about 25% of the total small rnas were  <dig> nt in size, which was the most abundant class .
the high percentage of the 19-nt small rnas in the wc library may result from the increased amount of rna degradation products caused by the cold treatment.
overall, approximately one-fourth of mirnas showed altered expression  in the cold-stressed brachypodium seedlings, while the other mirnas were expressed almost equally in the nc and wc libraries.
in response to cold treatment,  <dig> known mirnas and  <dig> predicted mirnas showed significant changes  in their expression levels.
the most obvious change was observed for mir <dig>  whose expression level increased about 15-fold in the wc library compared with that the nc library.
the expression of mir169e and mir172b also showed more than 5-fold of increases in the wc library.
the cold-induction of mir172a  was not as significant as that of mir172b.
among the  <dig> predicted mirnas with altered expression under cold stress,  <dig> were up-regulated and  <dig> were down-regulated  .
based on sequence similarity, these  <dig> predicted mirnas were classified into  <dig> mirna families.
the sequences of mir913t, mir914t, mir915t showed high similarity and could be grouped as one family.
it is the same case for mir911t, mir926t and mir927t as well as mir917t, mir922t and mir928t.
the mir918t and mir919t could also be grouped as one family.
for these mirna families, all the members showed decreased expression in response to the cold treatment.
analysis of the isolation frequencies for these  <dig> cold-responsive predicted mirnas indicated that most of them belonged to the moderately expressed mirnas, with isolation frequencies in the range of  <dig> and  <dig> .
the nc library was also searched for mirna* sequences of these cold-responsive predicted mirnas, and only the mirna* for mir912t and mir918t were found.
no mirna* sequence was detected for other predicted mirnas with cold response, probably because their isolation frequencies, about 10% of the frequencies of mature mirnas  <cit> , were too low.
to determine whether the  <dig> cold-related predicted mirnas are conserved among other plant species, the nucleotide databases were searched with blastn  <cit>  to identify their homologs and surrounding sequences.
these cold-related predicted mirnas were also aligned with the genome of the monocotyledonous model plant rice   using patscan  <cit> .
hairpin structures were predicted for these mirna homologs with the rnafold program  using the surrounding sequences.
some of these cold-related predicted mirnas  had homologs in the rice genome, some of them  also had homologs in the wheat and barley genome, but no homolog was found in the arabidopsis or populus genome, suggesting that these are monocot-specific mirnas.
for some of these cold-related predicted mirnas, such as mir901t-906t and mir923t-927t, no homolog was found in other plant species, suggesting that they are brachypodium-specific mirnas .
experimental validation of the cold-responsive expression of brachypodium mirnas
to confirm the expression of the identified mirnas and their response to cold stress, the expression of all the known cold-responsive mirnas and all the predicted mirnas whose expression changed significantly  after cold treatment were analyzed.
three known mirnas and eight predicted mirnas were selected for rna gel blot analysis.
the result showed that the expression of both mir <dig> and mir <dig> was up-regulated under the cold treatment .
the expression of mir169e was difficult to be detected due to its sequence similarity to other members in the mir <dig> family and low expression level.
clear changes were observed for the expression of five analyzed predicted mirnas  after cold treatment .
two of the analyzed predicted mirnas  did not show detectable hybridization signals, probably because their expression level was low .
their cold-responsive expression was then validated by real-time pcr analysis .
target predictions for cold-responsive brachypodium mirnas
it has been shown that plant mirnas exhibit high degree of sequence complementarity to their targets, which allows for effective target prediction  <cit> .
target prediction analysis, therefore, was performed for the cold-related bdi-mirnas.
for brachypodium, only draft genome sequence is now available, which has been annotated according to rice gene models as well as grass transcript assemblies and expressed sequence tags .
the accuracy and reliability of this notation needs further confirmation.
the ncbi  and other database  contain brachypodium ests, but the number is limited.
based on all the available data, the mirna targets were predicted manually using the patscan program  <cit>  with penalty scores ≤  <dig> for mismatched patterns in the mirna/mrna duplexes as described by lu et al.
<cit>   .
number in parentheses show the scores of predicted mirna targets based on rules described in  <cit> .
for conserved mirnas, their target sites should also be conserved across different plant species.
as expected, the targets for conserved cold-related mirnas in our system were similar or functionally related to previously validated plant mirna targets.
for example, bdi-mir <dig> could target laccase and laccase precursor  .
laccase has also been predicted as the target of mir <dig> in arabidopsis and rice  <cit> .
putative targets were also predicted for some newly identified cold-stress-related mirnas  .
the target of this mir901t, a cold-induced predicted mirna, is predicted to be a gene encoding lob domain protein which is proved to be involved in various aspects of plant development  <cit> .
discussion
the mirnas are a group of small non-coding rnas that play an important role in various developmental and stress response processes through negatively controlling of gene expression  <cit> .
they have been identified in diverse plant species, and the large-scale mirna identification in monocots with experimental approaches has mainly been performed in rice.
although mirna identification has also been reported in wheat  <cit> , extensive study was hampered by the unavailability of the whole genome sequence for wheat.
brachypodium, with close evolutionary relationship to wheat, has recently been proposed as a model temperate grass in monocots.
the genome sequence of brachypodium distachyon has been released ; .
to our knowledge, few studies have been done on identifying and analyzing mirnas in this plant species.
in this paper, genome-wide analysis of brachypodium mirnas were performed with high-throughput sequencing technology and their response to cold stress was also analyzed, providing useful information to deepen our understanding of mirnas in monocotyledonous plants.
the mirnas in brachypodium
solexa sequencing of the brachypodium small rna libraries revealed the existence of  <dig> conserved mirnas as well as  <dig> predicted mirnas , indicating that in brachypodium only a small fraction of mirnas are conserved mirnas, the rest of them appear to be species-specific mirnas.
most of the mirna families conserved between arabidopsis, populus and rice were also identified in our dataset   <cit> .
sequence analysis indicated that conserved mirnas in brachypoidium showed higher similarity to their homologs in rice than in arabidopsis and populus .
mir <dig>  previously only identified in rice  <cit> , was also found in brachypodium , indicating that it is not specific to rice.
for predicted mirnas, only these with cold stress response were further analyzed.
some of them had no homologs in other plant species and some appeared to be also conserved in wheat, rice and barley, but not in arabidopsis and populus , indicating that they were either brachypodium-specific or monocot-specific mirnas.
the predicted mirnas exhibited relatively low expression levels , in consistent with the notion that non-conserved mirnas are often expressed at a lower level than conserved mirnas.
these mirnas have not been reported in other plants before, possibly because their expression levels are low and need deeper sequencing to be discovered.
high-throughput sequencing technology and whole-genome-scale data mining enable us to discover probably most of the mirnas in brachypodium, which is an important resource for mirnas analysis.
deep sequencing of different species of brachypodium or different kinds of samples may lead us to identify more novel bdi-mirnas.
in monocotyledonous plants, the brachypodium has close phylogenetic relationships with rice.
both of them belong to the poaceae family, with rice in the ehrhartoideae subfamily and brachypodium in the pooideae subfamily.
the mirnas in rice have been examined extensively.
therefore, deep sequencing analysis of brachypodium provides a good chance to compare mirnas at the whole genome level in these two closely related plant species.
interestingly, although the bdi-mirnas and osa-mirnas show high sequence similarity, their genome organization and family size are quite different.
among the conserved brachypodium mirna families, mir <dig> is distinguished because its members are clustered in several plant genomes.
this family, therefore, has been selected for genome organization analysis.
according to the mirbase , the mir <dig> family has multiple clustered members in both eudicotyledonous and monocotyledonous plants.
it also appears in physcomitrella patens, the model moss species, but only has a single copy.
in rice, this family is encoded by four clusters of  <dig> genome loci, including two 7-gene clusters, one 6-gene cluster and one 4-gene cluster.
in brachypodium, the mir <dig> family is encoded by  <dig> genome loci and two clusters are formed , suggesting that duplication evens in brachypodium are not as active as that in rice.
although up to now, all the mir <dig> gene family members in other eudicotyledonous and monocotyledonous plant genomes are clustered, the brachypodium mir395d was not clustered with other paralogs .
because mir <dig> only has one member in lower plant physcomitrella, mir395d may be the ancient form of this family that is kept unchanged in brachypodium, possibly under some kind of positive selection.
it has been shown that the origin of the plant mirnas is dependent on the occurrence of various duplications, probably followed by chromosomal rearrangements and loss of duplicated genes  <cit> .
thus, an alternative explanation is that the unclustered mir395d may result from the loss of mir <dig> genes happened after duplication events during the evolution of this family.
if this is the case, the loss of mir <dig> genes probably occurred after the divergence of ehrhartoideae  and pooideae , because there is no unclustered mir <dig> gene in rice.
comparative analysis showed that the sizes of conserved multiple-member mirna families in brachypodium were smaller than those in rice and populus, most of them have similar sizes to those in arabidopsis .
for example, seven families  had similar number of members in brachypodium and arabidopsis, but their sizes were much larger in rice and populus .
considering the genome size of brachypodium , rice   <cit> , arabidopsis   <cit>  and populus   <cit> , it seems that the size of plant genomes has little effect on the number of mirna family members.
the genome size of brachypodium is close to that of rice and populus, compared with the size of the arabidopsis genome.
the mirna families in brachypodium, however, had much fewer members than those in rice and populus .
one reasonable explanation is that the duplication events during the evolution of these mirna families in brachypodium are less active or the selection pressure on duplicated mirnas is higher than that in rice and populus.
the great difference in mirna family size and genomic organization between brachypodium and rice indicates that mirna genes change their copy numbers and genome positions actively, probably to provide dosage effects or different regulatory patterns  in target gene regulation.
in addition, these data show that the evolution of mirna families in brachypodium is not as active as that in rice, suggesting that brachypodium could be used as a new platform for mirna studies in monocots.
the cold-responsive mirnas in brachypodium
although mirnas have been shown to play an important role in plant cold stress response, little information is available for monocotyledonous plants in this area  <cit> .
brachypodium is a potential model species for the cool season triticeae crops  <cit> .
to study the cold-responsive mirnas in this plant will provide useful information for improving the cold tolerance of economically important crops.
in this study, the expression levels of mirnas in brachypodium with and without cold treatment were compared and the results indicated that the expression of about one-fourth of mirnas was affected by cold stress.
among them,  <dig> conserved mirnas  and  <dig> predicted mirnas showed significant changes .
these  <dig> known cold-responsive mirnas are conserved among dicotyledonous and monocotyledonous plants.
mir <dig> was shown to be cold responsive in arabidopsis by sequencing of the stress-related small rna library  <cit> .
then, mir <dig> and mir <dig> were found to be responsive to cold stress in arabidopsis both through a computational, transcriptome-based approach and by microarray analysis almost simultaneously  <cit> .
mir <dig> and mir <dig> were also shown to be cold-upregulated in populus  <cit> .
the detection of the cold-induced expression for these mirnas in brachypodium indicates that they are also involved in cold stress response in monocotyledonous plants.
some previously identified cold-responsive conserved mirnas, such as mir <dig> in arabidopsis  <cit> , were not found in our small rna library, nor are they found in rice, suggesting that they probably are species-specific mirnas and do not exist in brachypodium.
some mirnas, including mir <dig>  mir <dig>  mir <dig>  and mir <dig>  respond to cold stress in arabidopsis  <cit> , but in our experiment no obvious change was found for these mirnas after the cold treatment.
one explanation for these discrepancies is that the induction levels of these mirnas are too low to be recognized as significant changes in our experimental system.
it is also possible that these mirnas show cold stress response only in specific tissues or at specific developmental stages in brachypodium.
previous studies put great emphasis on the cold-induced mirnas, whereas the cold-suppressed mirnas have received little research attention.
our study indicated that about one-third of cold-responsive mirnas were up-regulated and two-thirds of them were down-regulated, indicating that both kinds of regulation for mirnas were involved in cold response .
interestingly, in our study all the cold-responsive mirnas conserved in diverse plant species were up-regulated under cold stress condition.
for the  <dig> cold-down-regulated predicted mirnas, their expression was relatively low .
these data suggest that brachypodium responds to cold stress through both up- and down-regulation of mirna expression, in which the up-regulation is conserved and probably more important in cold response.
CONCLUSIONS
we have performed a genome-wide analysis of mirnas in brachypodium, providing a relatively complete view of bdi-mirnas.
taking advantage of the good reproducibility of high-throughput sequencing, cold-responsive conserved as well as predicted brachypodium mirnas were identified, delivering new insights into the role of mirnas in cold response of winter-habit monocotyledonous plants.
comparative studies of mirnas in brachypodium and rice suggest that the composition and location of mirna families are different even in closely related plant species.
our study provides useful information for further analysis in this area.
characterization of mirna targets and exploring mirna regulation mechanism will deepen our understanding of plant mirnas.

 micrornas are a class of small non-protein-coding rnas generated from single-stranded precursors with unique hairpin structures. they regulate the expression of mrnas by targeting transcripts for cleavage or translational repression. the mirnas were initially isolated in caenorhabditis elegans as developmental timing regulators. since then, they have been found in a broad range of plants, as well as viruses and mammals. plants mirnas were first identified in arabidopsis by different research groups. the biogenesis of plant mirnas is a complex multi-step enzymatic process. the mirnas are initially transcribed by rna polymerase ii in the cell nucleus as long primary mirnas that are cleaved into the mirna:mirna* duplexes by dicer-like. export of the duplexes into the cell cytoplasm is mediated by hasty. after methyl groups are added to the 3' ends of the duplexes catalyzed by hen one strand of the duplexes is selectively incorporated into the rna-induced silencing complex to form the mature mirnas, whereas the other strand, designated the mirna*, is typically degraded. dcl has also been shown to play a role in the biogenesis of a few mirnas with long hairpin precursors. plant mirnas exhibit high complementarity to their targets and can direct risc-mediated cleavage of target mrnas. thus, it is widely accepted that plant mirnas induce post-transcriptional gene silencing predominantly through guiding mrna cleavage. recent studies, however, have showed that translational inhibition is another kind of action mechanism for mirnas in plants. although initial studies have largely demonstrated the role of plant mirnas in development and morphogenesis processes, there are increasing number of reports indicating that plant mirnas also target genes involved in biotic as well as abiotic stress responses. low temperature is one of the most important environmental stimuli that affect plant growth and development. recently, more and more reports demonstrate the role of plant mirnas in cold stress response. the first study in this area, performed by sunkar and zhu, revealed the induced expression of mir mir397b and mir in response to cold stress as well as other kinds of stress treatment. in this study, mir319c appeared to be up-regulated by cold but not by dehydration, nacl, or aba treatment. zhou et al. identified cold up-regulated mirnas using a computational, transcriptome-based approach. liu et al. and lu et al. used microarray analysis to identify cold-responsive mirnas in arabidopsis and populus, respectively. these studies indicate that the expression of several mirnas is affected by the cold treatment. despite these efforts, our knowledge of the role played by mirnas in plant cold stress response is still limited at the whole-genome level. the most challenging problem in understanding plant mirnas is to identify more novel mirnas. three major approaches have been used for mirna discovery in plants: forward genetics, bioinformatic prediction as well as direct cloning and sequencing. only a few mirnas were identified by forward genetic studies and predicting species-specific mirnas using bioinformatics method is difficult. thus, direct cloning and sequencing is the most effective method for plant mirna discovery. only a few hundred mirnas have been identified with this approach, which leads to a premature conclusion that the types of mirnas in plants are limited. recently, the deep sequencing approach appeared and allowed the identification of numerous small rnas. it not only revealed a lot more species-specific mirnas, but also provided a picture of the genomic landscape of small rnas. although significant progress has been made in identifying plant mirnas and understanding their action mechanism, the discovery of novel mirnas in plants on a genome-wide scale is still at the preliminary stage. one limiting factor in mirna discovery is the availability of the whole genome sequence, only with which a comprehensive analysis of the potential hairpin precursor structures of cloned small rnas can be performed to distinguish mirnas from other kinds of small rnas. thus, most of the studies have been done in arabidopsis, rice and populus, whose whole genome sequences are known. many evolutionarily or economically important species have not been examined yet. to further understand the function of plant mirnas, more effort should be directed toward plant species with specific developmental features, which may contain mirnas that are specific for these features. several economically important winter-habit monocots, such as winter wheat, barley, oat, rye, as well as cool-season biofuel and forage grasses, can cold-acclimate and acquire high tolerance to low temperature. although it has been shown that mirnas are involved in plant cold stress response, little work has been done for monocotyledonous plant, especially for these winter-habit monocots, probably because rice, the model plant for monocots, is a tropical plant and incapable of cold acclimation. brachypodium distachyon has emerged as a new monocot model plant, especially for temperate cereals and related grasses. brachypodium not only has a closer evolutionary relationship with cool-season temperate cereals and grasses than rice, but also possesses growth and developmental features that are common to these plants. as a widely distributed winter-habit temperate plant, it has vernalization requirement and is capable of cold acclimation. recently, the draft genome sequence of brachypodium has been released, which makes this plant a good model for performing whole-genome-wide study of mirnas involved in cold response of winter-habit temperate cereals and grasses. here we sequenced small rna populations in brachypodium with and without cold stress treatment using solexa, the high-throughput sequencing technology. our studies not only identified species-specific mirnas for the novel monocotyledonous model plant, but also provided useful information for cold-responsive mirnas in temperate cereals and grasses.  deep sequencing of brachypodium small rnas two small rna libraries, with and without cold-treatment, were generated using pooled rna isolated from the aerial parts of brachypodium seedlings. sequencing of the brachypodium small rna libraries was performed with solexa, a high throughput sequencing technology producing highly accurate, reproducible and quantitative readouts of small rnas. it can be used as a tool for mirna expression profiling. for direct comparison, brachypodium seedlings used for nc and wc library construction were grown under similar conditions except for the cold treatment. the same amount of rna was used to construct these two libraries and the samples were prepared in a similar manner. sequencing of these two libraries was performed on the illumina's solexa sequencer and the samples were run side by side. solexa sequencing of nc and wc libraries generated a total of, and, raw reads, respectively. analysis of these reads resulted in identification of, and, sequences ranging from to nucleotides, respectively. the remaining sequences were either of low quality or smaller than nt, and were excluded from further analysis. the majority of the obtained small rna sequences from the nc library were 20- nt in size, which is the typical size range for dicer-derived products. these small rnas were used for further analysis. the 20- nt sequences from the nc and wc libraries were aligned to the draft brachypodium genome using soap, and a total of, and, sequences were found to match the genome perfectly. after removing the sequences corresponding to known rrnas, trnas, small nuclear rnas and small nucleolar rnas, a total of, and, sequences were obtained. the wc library contained more non-coding rnas than the nc library, which reduced the population size of non-coding-rna-removed genome-matched sequences in the wc library to half of that in the nc library. for direct comparison, the nc library was normalized to the wc library according to the population size mentioned above. in both nc and wc libraries, most of the known mirnas and predicted mirnas were detected at least times, enabling the comparison of mirna expression levels between these two libraries. identification of novel monocot-specific and brachypodium-specific mirnas one important feature that distinguishes mirnas from other endogenous small rnas is that their surrounding sequences can adopt a hairpin structure. the raw whole genome sequence data of brachypodium distachyon has been released by the us department of energy joint genome institute; and gene annotation was accomplished with fgenesh. in this study, the annotated brachypodium genome sequences were used to predict hairpin structures based on mirna surrounding sequences. to identify mirnas in brachypodium, the following strategy was adopted. first, the cloned sequences were aligned with the brachypodium genome sequence using soap to search for perfectly matched sequences. in the second step, sequences corresponding to known non-coding rnas were discarded by performing the blastn search against the genbank and rfam databases. the remaining sequences were used for the fold-back structure prediction and classified as either mirna candidates or other kinds of endogenous small rnas. our analysis revealed that in the nc library, sequences can be used for fold-back structure prediction and those that fulfilled the hairpin structure criteria described by jones-rhoades et al. were selected as candidate mirnas. among these candidates, some also matched the plastid or mitochondrial genomes. these sequences were not considered as mirnas in this study. after these analyses, unique sequences were obtained as mirna candidates. according to the currently available gene annotation of the brachypodium genome, of these mirna candidates were located in the exons of potentially protein-coding genes, which were not further analyzed. among the remaining brachypodium mirna candidates, of them have really low sequencing reads and were not detected in the wc library. those candidates were excluded from further research according to the updated plant mirna annotation criteria. the remaining putative brachypodium mirnas were carefully examined to be sure that they are qualified for the updated plant mirna annotation criteria. the lengths of the predicted brachypodium mirnas vary from to nt, and more than half of them begin with 5' uridine, which is a characteristic feature of mirnas. the minimum free energies of these mirna precursors range from - to - kcal mol- and the average value is about - kcal mol-, which is similar to the free energy values of mirna precursors in rice and arabidopsis. these values are much lower than the folding free energies of trnas or rrnas. the predicted hairpin structures for the mirna precursors are in the range of 71- nt, which is similar to what was observed in rice. none of these putative mirnas has been found in arabidopsis or rice before. identification of conserved mirnas in brachypodium to identify the conserved mirna homologs in brachypodium, small rna libraries were analyzed for the presence of known mirnas. blastn search with an e-value cutoff of was employed to search for small rnas with predicted hairpin structures against the central mirna registry database. with this similarity search, unique sequences were identified as known mirnas. they showed higher sequence similarity to their homologs in rice than in arabidopsis and populus. all of these conserved mirna precursors were identified in the draft brachypodium genome and they can adopt hairpin structures. based on sequence similarity, these known mirnas could be grouped into families. four conserved mirna families were represented by more than one sequence variant in the small rna libraries. mir mir and mir were represented by two variants and mir was represented by four variants in the library. some of these mirnas including mir mir mir and mi are deeply conserved, even in lower plants such as physcomitrella patens. mir and mir which are species-specific mirnas, were also found in the nc library. according to mirbase and published papers, mir only appears in oryza sativa and mir has been identified in eschscholzia californica, oryza sativa, physcomitrella patens as well as some pinus species. the mirna* sequences were detected for some bdi-mirnas. the isolation frequencies for most of these mirna* sequences were about one-tenth of those for mature mirnas, consistent with the conclusion of rajagopalan et al.. for mir there were reads of the mirna and reads of the mirna* in the nc library. similarly, in the wc library mir had reads of the mirna and reads of the mirna*. the mir529* was reported to be sequenced more often than the mir in rice and the annotated mir529*, therefore, is suspected to be the true mirna. our data indicate that mir is the true mirna and the existing annotation should not be changed. it is also interesting to note that the mir169f was sequenced times in the nc library and its potential mirna* was sequenced times, which is almost equal to the mature mirna sequencing frequency. name, the name of conserved brachypodium mirnas and underlines denote mirnas whose mirna* have been detected in solexa sequencing analysis; sequence, mirna sequence cloned in the nc library and; l, the length of mirna; arm/nt, the mirna location in the predicted hairpin structure /the length of the precursor; n, the number of loci in the genome for different mirna families; location, mirna location in the genome; ara, arabidopsis; pop, populus; ric, rice; whe, wheat; for plus symbols in the table: ++++, mirna sequences of brachypodium were identical to those in other plant species; +++, mirna sequences of brachypodium were conserved in other plant species but have variations at nucleotide positions; ++, mirna sequences of brachypodium were conserved in other plant species but have variations at nucleotide positions; +, mirna sequences of brachypodium were conserved in other plant species but have variations at ≥ nucleotide positions. for mirnas whose homologs are not identified in rice, sequence comparison was only performed for the other four plant species. the sequencing frequencies for mirnas in the library could be used as an index for estimating the relative abundance of mirnas. solexa sequencing produced a large number of mirna sequences, allowing us to determine the relative abundance of mirnas in brachypodium. the frequencies of the mirnas varied from to, indicating that the expression of mirnas varied greatly in brachypodium. the mir was the most abundant mirna in our two sequencing datasets, accounting for about 20% of all sequence reads matching with the brachypodium genome and about 60% of the sequence reads of the conserved mirnas. mir has also been the most abundant mirna in previous high-throughput sequencing studies in rice, but not in arabidopsis. according to the sequencing frequencies, some mirnas were highly expressed in brachypodium, and were sequenced more than times each. the mir160a, mir164d, mir169f, mir172a, mir172b, mir mir mir mir mir395a, mir397a, mir and mir were moderately expressed, and were represented by the number of sequences varying between and the third category, whose expression level was low, was represented by mir mir169a, mir169c, mir169e and mir. sequence analysis indicated that the relative abundance of certain members within a mirna family varied greatly in brachypodium. for instance, the sequencing frequency for mir166a was, whereas mir166g's sequencing frequency was only. these results indicate that in one mirna family different members show clearly different expression levels, probably because their expression is tissue- or developmental-specific. genomic organization of the mirnas among the conserved bachypodium mirnas, were encoded by a single copy in the brachypodium genome, whereas the other mirnas had multiple loci, which probably resulted from duplication events that were still active in the brachypodium genome. for these mirnas, most of them only had 2- loci in the genome, only a few had more than loci. for the predicted brachypodium mirnas, most of them had one single locus in the genome and only a few were represented by more than loci. according to the currently available gene annotation for the brachypodium genome, most of these predicted mirnas were mapped to the intergenic regions. interestingly, for almost all the conserved mirna families with multiple members, their family sizes in brachypodium are much smaller than those in rice and populus, and most of them have sizes similar to those in arabidopsis. for example, the number of family members for mir mir mir mir mir mir and mir in brachypodium was similar to that in arabidopsis, but much higher in rice and populus. for mir its size in brachypodium was approximately half of that in rice. name, the name of predicted cold-responsive brachypodium mirnas and underlines denote mirnas whose mirna* have been detected in solexa sequencing analysis; sequence, sequence cloned in small rna libraries; l, the length of mirna; arm/nt, the mirna location in the predicted hairpin structure /the length of the precursor; n, number of loci in the genome; location, mirna location in the genome; amfe, the adjusted minimum free energy representing the mfe of nucleotides, it is calculated by × 100; cd, the effect of cold-treatment on mirna expression, nc/wc, normalized sequencing frequencies in the nc and wc library. the number of mirna family members in rice, arabidopsis and populus was obtained from the mirbase. mir belongs to conserved mirnas that have been found across different plant species. both mir and mir395* could be detected in the nc library, although their reads were not high. according to the detected mirna and mirna* sequences as well as results of blastn search based on sequences similarity, precursors with reasonable minimum free energy value were identified for mir in the draft brachypodium genome. they were named mir395a~mir395n according to the order of their genomic location. additionally, two other sequences in the surrounding region showed high similarity to mir precursors, but they contained no mature mir sequence and were therefore not considered for further analysis in this study. for all the mir precursors, the mature mirnas are all located on the 3' arms. the precursors of mir395a~c are clustered together and so were the precursors of mir395e~m.. it is interesting to note that mir395d in brachypodium is not clustered with other family members although the members of this family always form clusters in other plant species. no mir395- precursor-like sequence was detected in its surrounding region. on the basis of size and sequence similarity, two types of mir precursors exist in brachypodium. type a includes mir395a~c and mir395f~n, and type b includes mir395d and mir395e. both kinds of precursors show high sequence similarity to their homologs in rice. the precursors of mir395a~c and mir395f~n show high sequence similarity, suggesting that they are derived from a series of gene duplication events. these precursors exhibit different levels of sequence similarity to one another, implying that the duplication events happened at different time points during the evolution history of the mir family. the precursors of mir395i and mir395k are nearly identical, implying a recent origin. the impact of cold stress on conserved and species-specific mirnas in brachypodium to detect the effect of cold stress on brachypodium mirnas, the expression of mirnas in brachypodium seedlings with and without cold-treatment were examined with solexa technology. these two libraries were compared for their size distribution of small rnas. in the nc library, the 21- nt rnas represented the predominant species, with 21-nt and 24-nt rnas being the two most abundant classes, consistent with the distribution patterns of small rnas from other plant species. however, in the wc library, about 25% of the total small rnas were nt in size, which was the most abundant class. the high percentage of the 19-nt small rnas in the wc library may result from the increased amount of rna degradation products caused by the cold treatment. overall, approximately one-fourth of mirnas showed altered expression in the cold-stressed brachypodium seedlings, while the other mirnas were expressed almost equally in the nc and wc libraries. in response to cold treatment, known mirnas and predicted mirnas showed significant changes in their expression levels. the most obvious change was observed for mir whose expression level increased about 15-fold in the wc library compared with that the nc library. the expression of mir169e and mir172b also showed more than 5-fold of increases in the wc library. the cold-induction of mir172a was not as significant as that of mir172b. among the predicted mirnas with altered expression under cold stress, were up-regulated and were down-regulated. based on sequence similarity, these predicted mirnas were classified into mirna families. the sequences of mir913t, mir914t, mir915t showed high similarity and could be grouped as one family. it is the same case for mir911t, mir926t and mir927t as well as mir917t, mir922t and mir928t. the mir918t and mir919t could also be grouped as one family. for these mirna families, all the members showed decreased expression in response to the cold treatment. analysis of the isolation frequencies for these cold-responsive predicted mirnas indicated that most of them belonged to the moderately expressed mirnas, with isolation frequencies in the range of and. the nc library was also searched for mirna* sequences of these cold-responsive predicted mirnas, and only the mirna* for mir912t and mir918t were found. no mirna* sequence was detected for other predicted mirnas with cold response, probably because their isolation frequencies, about 10% of the frequencies of mature mirnas, were too low. to determine whether the cold-related predicted mirnas are conserved among other plant species, the nucleotide databases were searched with blastn to identify their homologs and surrounding sequences. these cold-related predicted mirnas were also aligned with the genome of the monocotyledonous model plant rice using patscan. hairpin structures were predicted for these mirna homologs with the rnafold program using the surrounding sequences. some of these cold-related predicted mirnas had homologs in the rice genome, some of them also had homologs in the wheat and barley genome, but no homolog was found in the arabidopsis or populus genome, suggesting that these are monocot-specific mirnas. for some of these cold-related predicted mirnas, such as mir901t-906t and mir923t-927t, no homolog was found in other plant species, suggesting that they are brachypodium-specific mirnas. experimental validation of the cold-responsive expression of brachypodium mirnas to confirm the expression of the identified mirnas and their response to cold stress, the expression of all the known cold-responsive mirnas and all the predicted mirnas whose expression changed significantly after cold treatment were analyzed. three known mirnas and eight predicted mirnas were selected for rna gel blot analysis. the result showed that the expression of both mir and mir was up-regulated under the cold treatment. the expression of mir169e was difficult to be detected due to its sequence similarity to other members in the mir family and low expression level. clear changes were observed for the expression of five analyzed predicted mirnas after cold treatment. two of the analyzed predicted mirnas did not show detectable hybridization signals, probably because their expression level was low. their cold-responsive expression was then validated by real-time pcr analysis. target predictions for cold-responsive brachypodium mirnas it has been shown that plant mirnas exhibit high degree of sequence complementarity to their targets, which allows for effective target prediction. target prediction analysis, therefore, was performed for the cold-related bdi-mirnas. for brachypodium, only draft genome sequence is now available, which has been annotated according to rice gene models as well as grass transcript assemblies and expressed sequence tags. the accuracy and reliability of this notation needs further confirmation. the ncbi and other database contain brachypodium ests, but the number is limited. based on all the available data, the mirna targets were predicted manually using the patscan program with penalty scores ≤ for mismatched patterns in the mirna/mrna duplexes as described by lu et al.. number in parentheses show the scores of predicted mirna targets based on rules described in. for conserved mirnas, their target sites should also be conserved across different plant species. as expected, the targets for conserved cold-related mirnas in our system were similar or functionally related to previously validated plant mirna targets. for example, bdi-mir could target laccase and laccase precursor. laccase has also been predicted as the target of mir in arabidopsis and rice. putative targets were also predicted for some newly identified cold-stress-related mirnas. the target of this mir901t, a cold-induced predicted mirna, is predicted to be a gene encoding lob domain protein which is proved to be involved in various aspects of plant development. discussion the mirnas are a group of small non-coding rnas that play an important role in various developmental and stress response processes through negatively controlling of gene expression. they have been identified in diverse plant species, and the large-scale mirna identification in monocots with experimental approaches has mainly been performed in rice. although mirna identification has also been reported in wheat, extensive study was hampered by the unavailability of the whole genome sequence for wheat. brachypodium, with close evolutionary relationship to wheat, has recently been proposed as a model temperate grass in monocots. the genome sequence of brachypodium distachyon has been released;. to our knowledge, few studies have been done on identifying and analyzing mirnas in this plant species. in this paper, genome-wide analysis of brachypodium mirnas were performed with high-throughput sequencing technology and their response to cold stress was also analyzed, providing useful information to deepen our understanding of mirnas in monocotyledonous plants. the mirnas in brachypodium solexa sequencing of the brachypodium small rna libraries revealed the existence of conserved mirnas as well as predicted mirnas, indicating that in brachypodium only a small fraction of mirnas are conserved mirnas, the rest of them appear to be species-specific mirnas. most of the mirna families conserved between arabidopsis, populus and rice were also identified in our dataset. sequence analysis indicated that conserved mirnas in brachypoidium showed higher similarity to their homologs in rice than in arabidopsis and populus. mir previously only identified in rice, was also found in brachypodium, indicating that it is not specific to rice. for predicted mirnas, only these with cold stress response were further analyzed. some of them had no homologs in other plant species and some appeared to be also conserved in wheat, rice and barley, but not in arabidopsis and populus, indicating that they were either brachypodium-specific or monocot-specific mirnas. the predicted mirnas exhibited relatively low expression levels, in consistent with the notion that non-conserved mirnas are often expressed at a lower level than conserved mirnas. these mirnas have not been reported in other plants before, possibly because their expression levels are low and need deeper sequencing to be discovered. high-throughput sequencing technology and whole-genome-scale data mining enable us to discover probably most of the mirnas in brachypodium, which is an important resource for mirnas analysis. deep sequencing of different species of brachypodium or different kinds of samples may lead us to identify more novel bdi-mirnas. in monocotyledonous plants, the brachypodium has close phylogenetic relationships with rice. both of them belong to the poaceae family, with rice in the ehrhartoideae subfamily and brachypodium in the pooideae subfamily. the mirnas in rice have been examined extensively. therefore, deep sequencing analysis of brachypodium provides a good chance to compare mirnas at the whole genome level in these two closely related plant species. interestingly, although the bdi-mirnas and osa-mirnas show high sequence similarity, their genome organization and family size are quite different. among the conserved brachypodium mirna families, mir is distinguished because its members are clustered in several plant genomes. this family, therefore, has been selected for genome organization analysis. according to the mirbase, the mir family has multiple clustered members in both eudicotyledonous and monocotyledonous plants. it also appears in physcomitrella patens, the model moss species, but only has a single copy. in rice, this family is encoded by four clusters of genome loci, including two 7-gene clusters, one 6-gene cluster and one 4-gene cluster. in brachypodium, the mir family is encoded by genome loci and two clusters are formed, suggesting that duplication evens in brachypodium are not as active as that in rice. although up to now, all the mir gene family members in other eudicotyledonous and monocotyledonous plant genomes are clustered, the brachypodium mir395d was not clustered with other paralogs. because mir only has one member in lower plant physcomitrella, mir395d may be the ancient form of this family that is kept unchanged in brachypodium, possibly under some kind of positive selection. it has been shown that the origin of the plant mirnas is dependent on the occurrence of various duplications, probably followed by chromosomal rearrangements and loss of duplicated genes. thus, an alternative explanation is that the unclustered mir395d may result from the loss of mir genes happened after duplication events during the evolution of this family. if this is the case, the loss of mir genes probably occurred after the divergence of ehrhartoideae and pooideae, because there is no unclustered mir gene in rice. comparative analysis showed that the sizes of conserved multiple-member mirna families in brachypodium were smaller than those in rice and populus, most of them have similar sizes to those in arabidopsis. for example, seven families had similar number of members in brachypodium and arabidopsis, but their sizes were much larger in rice and populus. considering the genome size of brachypodium, rice, arabidopsis and populus, it seems that the size of plant genomes has little effect on the number of mirna family members. the genome size of brachypodium is close to that of rice and populus, compared with the size of the arabidopsis genome. the mirna families in brachypodium, however, had much fewer members than those in rice and populus. one reasonable explanation is that the duplication events during the evolution of these mirna families in brachypodium are less active or the selection pressure on duplicated mirnas is higher than that in rice and populus. the great difference in mirna family size and genomic organization between brachypodium and rice indicates that mirna genes change their copy numbers and genome positions actively, probably to provide dosage effects or different regulatory patterns in target gene regulation. in addition, these data show that the evolution of mirna families in brachypodium is not as active as that in rice, suggesting that brachypodium could be used as a new platform for mirna studies in monocots. the cold-responsive mirnas in brachypodium although mirnas have been shown to play an important role in plant cold stress response, little information is available for monocotyledonous plants in this area. brachypodium is a potential model species for the cool season triticeae crops. to study the cold-responsive mirnas in this plant will provide useful information for improving the cold tolerance of economically important crops. in this study, the expression levels of mirnas in brachypodium with and without cold treatment were compared and the results indicated that the expression of about one-fourth of mirnas was affected by cold stress. among them, conserved mirnas and predicted mirnas showed significant changes. these known cold-responsive mirnas are conserved among dicotyledonous and monocotyledonous plants. mir was shown to be cold responsive in arabidopsis by sequencing of the stress-related small rna library. then, mir and mir were found to be responsive to cold stress in arabidopsis both through a computational, transcriptome-based approach and by microarray analysis almost simultaneously. mir and mir were also shown to be cold-upregulated in populus. the detection of the cold-induced expression for these mirnas in brachypodium indicates that they are also involved in cold stress response in monocotyledonous plants. some previously identified cold-responsive conserved mirnas, such as mir in arabidopsis, were not found in our small rna library, nor are they found in rice, suggesting that they probably are species-specific mirnas and do not exist in brachypodium. some mirnas, including mir mir mir and mir respond to cold stress in arabidopsis, but in our experiment no obvious change was found for these mirnas after the cold treatment. one explanation for these discrepancies is that the induction levels of these mirnas are too low to be recognized as significant changes in our experimental system. it is also possible that these mirnas show cold stress response only in specific tissues or at specific developmental stages in brachypodium. previous studies put great emphasis on the cold-induced mirnas, whereas the cold-suppressed mirnas have received little research attention. our study indicated that about one-third of cold-responsive mirnas were up-regulated and two-thirds of them were down-regulated, indicating that both kinds of regulation for mirnas were involved in cold response. interestingly, in our study all the cold-responsive mirnas conserved in diverse plant species were up-regulated under cold stress condition. for the cold-down-regulated predicted mirnas, their expression was relatively low. these data suggest that brachypodium responds to cold stress through both up- and down-regulation of mirna expression, in which the up-regulation is conserved and probably more important in cold response.  we have performed a genome-wide analysis of mirnas in brachypodium, providing a relatively complete view of bdi-mirnas. taking advantage of the good reproducibility of high-throughput sequencing, cold-responsive conserved as well as predicted brachypodium mirnas were identified, delivering new insights into the role of mirnas in cold response of winter-habit monocotyledonous plants. comparative studies of mirnas in brachypodium and rice suggest that the composition and location of mirna families are different even in closely related plant species. our study provides useful information for further analysis in this area. characterization of mirna targets and exploring mirna regulation mechanism will deepen our understanding of plant mirnas. 
both kinds of mirnas were involved in cold stress response, but all the conserved mirnas were up-regulated, implying an important role for cold-induced mirnas. although mirnas are proved to be involved in plant cold stress responses, little is known for winter-habit monocots. the different size and genome organization of mirna families in brachypodium and rice suggest that the frequency of duplication events or the selection pressure on duplicated mirnas are different between these two closely related plant species. brachypodium distachyon, with close evolutionary relationship to cool-season cereals, has recently emerged as a novel model plant. there are few reports of brachypodium mirnas. our results suggest that brachypodium mirnas are composed of a set of conserved mirnas and a large proportion of non-conserved mirnas with low expression levels. the genome organization of mir family in brachypodium was quite different from that in rice. among these mirnas, some were cold-induced and some were cold-suppressed, but all the conserved mirnas were up-regulated under cold stress condition. the expression of conserved mirnas and predicted mirnas showed significant changes in response to cold stress. for multiple-member conserved mirna families, their sizes in brachypodium were much smaller than those in rice and populus. micrornas are endogenous small rnas having large-scale regulatory effects on plant development and stress responses.
 micrornas are endogenous small rnas having large-scale regulatory effects on plant development and stress responses. extensive studies of mirnas have only been performed in a few model plants. although mirnas are proved to be involved in plant cold stress responses, little is known for winter-habit monocots. brachypodium distachyon, with close evolutionary relationship to cool-season cereals, has recently emerged as a novel model plant. there are few reports of brachypodium mirnas.  high-throughput sequencing and whole-genome-wide data mining led to the identification of conserved mirnas, as well as predicted mirnas in brachypodium. for multiple-member conserved mirna families, their sizes in brachypodium were much smaller than those in rice and populus. the genome organization of mir family in brachypodium was quite different from that in rice. the expression of conserved mirnas and predicted mirnas showed significant changes in response to cold stress. among these mirnas, some were cold-induced and some were cold-suppressed, but all the conserved mirnas were up-regulated under cold stress condition.  our results suggest that brachypodium mirnas are composed of a set of conserved mirnas and a large proportion of non-conserved mirnas with low expression levels. both kinds of mirnas were involved in cold stress response, but all the conserved mirnas were up-regulated, implying an important role for cold-induced mirnas. the different size and genome organization of mirna families in brachypodium and rice suggest that the frequency of duplication events or the selection pressure on duplicated mirnas are different between these two closely related plant species. 