Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane
- Almir S Zanca†1,
- Renato Vicentini†1,
- Fausto A Ortiz-Morea2,
- Luiz EV Del Bem1,
- Marcio J da Silva1,
- Michel Vincentz1, 3 and
- Fabio TS Nogueira2, 4Email author
© Zanca et al; licensee BioMed Central Ltd. 2010
Received: 27 April 2010
Accepted: 24 November 2010
Published: 24 November 2010
MicroRNAs (miRNAs) are small regulatory RNAs, some of which are conserved in diverse plant genomes. Therefore, computational identification and further experimental validation of miRNAs from non-model organisms is both feasible and instrumental for addressing miRNA-based gene regulation and evolution. Sugarcane (Saccharum spp.) is an important biofuel crop with publicly available expressed sequence tag and genomic survey sequence databases, but little is known about miRNAs and their targets in this highly polyploid species.
In this study, we have computationally identified 19 distinct sugarcane miRNA precursors, of which several are highly similar with their sorghum homologs at both nucleotide and secondary structure levels. The accumulation pattern of mature miRNAs varies in organs/tissues from the commercial sugarcane hybrid as well as in its corresponding founder species S. officinarum and S. spontaneum. Using sugarcane MIR827 as a query, we found a novel MIR827 precursor in the sorghum genome. Based on our computational tool, a total of 46 potential targets were identified for the 19 sugarcane miRNAs. Several targets for highly conserved miRNAs are transcription factors that play important roles in plant development. Conversely, target genes of lineage-specific miRNAs seem to play roles in diverse physiological processes, such as SsCBP1. SsCBP1 was experimentally confirmed to be a target for the monocot-specific miR528. Our findings support the notion that the regulation of SsCBP1 by miR528 is shared at least within graminaceous monocots, and this miRNA-based post-transcriptional regulation evolved exclusively within the monocots lineage after the divergence from eudicots.
Using publicly available nucleotide databases, 19 sugarcane miRNA precursors and one new sorghum miRNA precursor were identified and classified into 14 families. Comparative analyses between sugarcane and sorghum suggest that these two species retain homologous miRNAs and targets in their genomes. Such conservation may help to clarify specific aspects of miRNA regulation and evolution in the polyploid sugarcane. Finally, our dataset provides a framework for future studies on sugarcane RNAi-dependent regulatory mechanisms.
MicroRNAs (miRNAs) are small regulatory RNAs (19-21 nt) that play crucial roles in diverse aspects of plant development [1–3], biotic and abiotic stress responses [4, 5], signal transduction and protein degradation [6, 7]. MiRNAs are generated by stepwise processing of RNA polymerase II (Pol II)-dependent primary miRNA transcripts (pri-miRNAs). The pri-miRNAs typically form an imperfect fold-back structure, which is processed into a stem-loop precursor (pre-miRNA) and further excised as an RNA duplex by the DICER-LIKE1 (DCL1) enzyme. Partial or complete base-pairing between the miRNA and its target RNA allows the miRNA-associated RNA-induced silencing complexes (RISCs) to promote translational inhibition, accelerated exonucleolytic mRNA decay, and/or mRNA cleavage through slicing within miRNA-mRNA base-pairing (for review, see ). The majority of the target genes of highly conserved miRNAs are transcription factors that play important roles in development . Conversely, lineage-specific miRNAs seem to regulate the expression of a broader type of genes, including those involved in cellular metabolism, stress response, and post-translational modifications [7, 10].
The identification of miRNAs and their targets in a large number of plant species is an important step to understand the function and evolution of miRNAs and miRNA-dependent gene regulation. Over 1,300 miRNAs from eudicotyledoneous and 832 miRNAs from monocotyledonous plants have been deposited in the latest release of miRBase (release 14.0 September 2009). Although deep sequencing methods have substantially contributed to the identification of conserved and lineage-specific miRNAs in model species , these approaches are time-consuming and relatively expensive. In this regard, public EST databases and genomic survey sequences (GSSs) have become attractive alternatives to identify non-coding sequences through computational approaches in non-model plants.
The fact that most known miRNAs are evolutionarily conserved raises the possibility of identifying new miRNA homologs in other species using computer-based strategies , and such in silico approaches have been reviewed and classified not only as homology-based but also as structure similarity-based searches [7, 13]. Therefore, recent computational methods provide an accurate, fast, inexpensive, and consequently convenient way to retrieve miRNA precursor sequences from publicly available sequence databases. Finally, target mRNAs of conserved miRNAs can be searched using web-based  or in-house algorithms and analyzed across plant species.
Most identified miRNAs and their targets have been predicted in plants for which whole genome information is available such as Arabidopsis thaliana and rice. Currently, there is no experimental and only scarce computational information about miRNAs and their targets in sugarcane (Saccharum spp.). Sugarcane is an economically important biofuel crop. Recently, it has become a target for improvement of sustainable biomaterial production due to its high biomass productivity and built-in containment features . Modern sugarcane cultivars are highly polyploid, aneuploid hybrids between S. officinarum L. (octoploid, with 2n = 80 chromosomes) and S. spontaneum L. (ploidy level of 5-16, with 2n = 40-128 chromosomes). Modern sugarcane cultivars typically have 2n = 100-130 chromosomes, of which approximately 15-20% are derived from S. spontaneum and 5% are recombinants derived from both species. Therefore, the genome of modern sugarcane cultivars has at least 10 copies of most homo(eo)logous loci, contributing to the high complexity of its genome .
In this study, we used conserved miRNAs to systematically search public EST and GSS databases for sugarcane pri-miRNAs or miRNA precursors. A total of 19 distinct sugarcane pri-miRNAs were identified by our computational protocol, of which nine are monocot-specific. The expression profiles of selected sugarcane miRNAs were monitored by pulsed stem-loop RT-PCR  in organs/tissues of a modern cultivar as well as in S. officinarum and S. spontaneum. To identify target genes of the identified miRNAs, we developed a BLAST-based computational tool to search the NCBI EST and BAC sequences of sugarcane, rice, and sorghum. By using this method, we predicted several target messages, of which one novel target was experimentally tested and confirmed. Finally, we integrated sugarcane miRNA primary precursor and target information into a web-based database (http://sysbiol.cbmeg.unicamp.br/SCmiRNA), which is publicly available. The identification of miRNAs and their targets is important not only to help us learn more about the roles of miRNAs in sugarcane development and physiology, but also to provide a framework for further studies on RNAi-based regulation mechanisms in this highly polyploid species.
Results and Discussion
Identification of miRNA primary transcripts in Saccharum spp
Distinct sugarcane pri-miRNAs identified in this study
miRNA Mature Sequence
Conserved in rice
The sugarcane miRNA families identified in this study include the six families already deposited in the miRbase (Table 1), indicating the robustness of our approach. Nonetheless, careful inspection of the sugarcane miRNA precursor sequences deposited in the miRbase v.14 and comparison with our analysis revealed some divergences between these databases. For example, the SsMIR156b/c (Table 1) was previously annotated as a single stem-loop MIR156 precursor (miRbase v.14). However, our analyses revealed that this precursor belongs to a cluster representing a two-tandem microRNA precursor, which is highly similar to its sorghum homolog (90% nt identity) and to the maize Corngrass1 microRNA (84% nt identity) . Moreover, genomic DNA PCR amplification from sugarcane hybrid RB 83-5486 using specific primers and subsequent sequencing indicate that SsMIR156b/c locus encodes tandem MIR156 genes (data not shown). Comparison among the miRbase-derived precursor sequences and with those identified in this study suggests that the 16 previously annotated sugarcane miRNA precursors represent only eight different precursors (Table 1). For instance, we identified only two distinct precursors of miR408, SsMIR408a and SsMIR408b (Table 1), instead of five (miRbase v.14). Closer inspection suggests that SsMIR408a and SsMIR408b are likely different alleles of the same locus. This observation is supported by the fact that MIR408 genes have been found only as one copy in all plant genomes evaluated to date (miRbase v.14). The discrepancies between our data and previous annotation in the miRbase may be due to the use of SoGI Release 2.2 (July, 2008) that contains substantially more Tentative Consensus (TCs) than the earlier releases, which likely reflect differences in EST clustering or assembling.
In agreement with previous results [7, 10], most sugarcane miRNA sequences have uracil as their first nucleotide (13 out of 19 mature sugarcane miRNAs; Table 1). Moreover, sugarcane miRNA precursors displayed high minimal free energy index (MFEI) values (average 1.02 ± 0.22), which is a criterion used for distinguishing miRNA precursors from other types of RNAs. MFEI is a parameter that considers not only the minimal free energy (MFE) value of a particular sequence but also its length and G+C content. MFEI values were calculated as described by .
MiRNAs are located either in the 5'-arm or 3'-arm of the stem-loop hairpin pre-miRNA sequences (Table 1). All new identified miRNA precursors could fold into stem-loop structures (see additional file 1), following the rules and parameters reported by . One exception was the EST TC87836, which displays high similarity (e-value 0.0 and 89.5% nt identity) with one of the MIR319 precursors present in the draft of the Sorghum bicolor genome  (see additional file 2). It could not form a suitable stem-loop structure and thus it was not validated by our in-house MIRcheck-based script. This might be due to the fact that the miR319* is located at the 5'end of the sequence, which is not present in the TC87836 sugarcane EST. Nevertheless, based on its extensive homology with sorghum MIR319 precursor (see additional file 2), we annotated this TC as a potential SsMIR319 precursor (Table 1). That most ESTs do not contain their entire 5'-end sequence information undermines EST databases as sources for miRNA precursor searching. Based on the example given in this study, it may be interesting to develop rules and parameters to assign EST sequences as miRNA precursors based only upon extensive nucleotide identity with precursors from highly closely related species.
Not only SsMIR319, but several sugarcane pre-miRNAs show high sequence similarity with their sorghum homologs (values between 86% and 94% nt identity). Sorghum and sugarcane are each other's closest relatives among cultivated crops. They belong to the Andropogoneae tribe and diverged from a common ancestor around 8-9 Myr ago . Based on genomic sequence comparisons [15, 20], it is likely that sugarcane and sorghum did not have sufficient time to diverge, which reflects the high degree of identity observed between their miRNA precursors.
Among the monocot-specific miRNA precursors, we have identified three potential precursors of microRNA444 (Table 1). Interestingly, SsMIR444b and SsMIR444c contain tandem and overlapping mature miRNA sequences (additional file 1), similar to MIR444 precursors identified in rice and sorghum [20, 22]. At least in rice, such precursors are able to generate natural antisense miRNAs, or nat-miRNAs. The production of nat-miRNAs depends upon sense/antisense transcription and alternative splicing of the precursors prior to DCL1 cleavage . These nat-miRNAs seem to be restricted to monocot graminae, indicating this new pathway is less than 50 million years old.
Given the limited number of sugarcane EST and GSS sequences available as well as existent sequencing errors, the frequency of candidate miRNA precursors identified in this study is comparable to others using such databases [29, 30]. It is noteworthy that all miRNAs reported in this study have been identified using previously known miRNAs from several plant species. Therefore, we did not uncover miRNAs that are specific to sugarcane. Further investigations that employ small RNA libraries combined with computational approaches are needed to identify sugarcane-specific miRNAs.
Expression patterns of sugarcane miRNAs
We also compared the expression profiles of these miRNAs between S. officinarum and S. spontaneum to evaluate whether both species produce detectable mature miRNA molecules. All miRNAs are detected in the evaluated organs/tissues from these two closely related species (Figure 3B). Although most miRNAs seem to accumulate similarly in both species, some presented variations in abundance when comparing the same organs/tissues at similar developmental stages of S. officinarum and S. spontaneum. For example, miR444 is slightly more abundant in vegetative apex of S. officinarum. In contrast, miR408 accumulates at higher levels in leaf blade tissues of S. spontaneum (Figure 3B). Similar data was observed for miRNAs accumulating in some organs/tissues of stable Arabidopsis allopolyploids . The relatively low variation in miRNA accumulation between these species is likely a reflection of their level of ploidy. Highly polyploid species might have developed a genetic buffering against extensive miRNA expression variation in particular organs/tissues or developmental stages to maintain target gene expression stability across generations of ploidy . Our data also present the possibility that both ancient species contributed similarly to the miRNA-based regulatory pathways present in modern sugarcane hybrids. It will be interesting to test whether all target loci in hybrid modern cultivars are down-regulated by miRNAs from one ancient progenitor or from both.
Although we did not detect any tandem repeats, the CpG/CpNpG islands found in the three MIR1432 loci overlap broadly with the possibly conserved regions upstream of pre-miR1432. These regions also included common predicted TFB sites for the investigated species, such as an auxin response element (AuxRE) (Figure 4). Taken together, these findings suggest the promoter regions of the sugarcane MIR1432 locus share conserved elements with its sorghum and maize homologs. Such elements might be biologically important for the final organ/tissue localization of miR1432 mature species and, consequently, for target down-regulation. It has recently been reported an evolutionary sequencing comparison for the MIR319a locus in Arabidopsis and related Brassicaceae. Reporter experiments have demonstrated that regions under stronger evolutionary constraints contain important information for MIR319a transcription . As more sugarcane genomic sequences become available, it will be interesting to verify whether most, if not all homologous miRNAs between sorghum and sugarcane also share conserved elements in their promoters.
Potential targets of sugarcane miRNAs
Target genes for 14 sugarcane miRNA families
TC96571(1), TC76232(1), CF571975(1), CA188863(1), TC98485(1), TC95509(1), TC87521(1)
Squamosa promoter binding protein (SBP)
Myb protein-like protein
CA232593(4), CA292372(4), CA201181(4), CA070734(4), TC109486(4), TC94174(4), TC82328(4), TC76219(4), TC74847(4)
Auxin response factor
TC85238(4), CA099906(4), TC85149(4), TC74893(4)
ARGONAUTE1 (AGO1)-like protein
TC106305(4), TC103899(4), TC100837(4)
CCAAT-box transcription factor complex
TCP family transcription factor
CA133327(2), TC105472(2), TC82603(2)
Basic blue copper protein-like protein
SPX domain-containing protein
protein kinase-like protein
MIKC MADS-box transcription factor
CA208272 (2), CA282606 (2)
CA105902 (2), TC90577 (2)
TC90826(1), TC99902(1), CA213231(1)
Cu2+-binding domain-containing protein
To gain more insight into the evolutionary history of the CBPs, we performed a phylogenetic analysis using SsCBP1 sequence as a query to search for homologous proteins within genomic and EST databanks of a set of green plants, including angiosperms, basal land plants, and green algae (Viridiplantae 1.0; see Methods). Our analysis revealed that SsCBP1 belongs to a Possible Group of Orthologous (PoGO A; for a definition and criteria for PoGO, see ) that integrates only angiosperm sequences (Figure 5C). The simplest explanation is that these genes share a common origin within the last common ancestor of angiosperms. Interestingly, the miR528-target recognition site is only present within monocot genes from PoGO A (data not shown). All eudicots orthologous to SsCBP1 from Arabidopsis, poplar, grape, and soybean genomes completely lack the miR528-target recognition site, suggesting that miR528 is indeed a monocot-specific microRNA. Taken together, our findings support the notion that the regulation of SsCBP1 by miR528 is shared at least within graminaceous monocots, and this miRNA-based post-transcriptional regulation evolved exclusively within the monocots lineage after the divergence from eudicots. Further studies on plant CBPs are needed to define their physiological role(s) and the possible evolutionary advantages given by the miR528-based post-transcriptional regulation of monocot SsCBP1 orthologs.
Our findings indicate that several sugarcane miRNA precursors share high homology with their sorghum's possible orthologous beyond miRNA mature sequence. In the case of pre-miR1432, which was obtained from genomic sequences, we found precursor-surrounding regions conserved among sugarcane, sorghum, and maize homologs. This finding indicates these genes may share common genetic and epigenetic regulatory programs. However, further work that includes additional homologous sequences from other closely related species is required to confirm such conservation. Our data also indicate that sugarcane miRNAs are expressed in commercial hybrids as well as in the ancient progenitors S. officinarum and S. spontaneum. Our approach leads to the prediction of several conserved and non-conserved sugarcane miRNA targets in the available EST and genomic databases. The data is available in the public website (http://sysbiol.cbmeg.unicamp.br/SCmiRNA) that will be continuously updated to incorporate future miRBase updates. Our findings will be a useful resource toward tracing the evolution of small RNA-based regulation in sugarcane and related species. Most importantly, this study will serve as a foundation for future research into the functional roles of miRNAs and their target genes in this important biofuel crop.
Plant material and RNA extraction
Leaf blade and sheath tissues were collected from three-week-old sugarcane seedlings (hybrid RB 83-5486) grown in greenhouse conditions. Mature six-month-old plants of the same hybrid were used to obtain lateral buds and leaf roll (apical meristem plus leaf primordia) tissues. We also collected tissues from Saccharum officinarum (accession Muntok, Java) and S. spontaneum (accession SES205A). Five-month-old plantlets cultivated in vitro were transferred to greenhouse conditions. After one month, leaf blade and vegetative apex tissues (pool of four plantlets) were harvested from both species. Tissues were also collected from whole three-week-old seedlings of sugarcane hybrid (RB83-5486) and Sorghum bicolor (BTx623), four-week-old seedlings of Zea mays, one-month-old plantlets of Oryza sativa (ssp. japonica cv Nipponbare), and from one-month-old plantlets of Arabidopsis thaliana (Columbia). Total RNA was extracted using Trizol reagent according to manufacturer's instructions.
Stem-loop reverse transcriptase (RT)-PCR
Stem-loop RT and PCR primers for sugarcane miR408, miR156, miR444, miR528, miR1128, and miR1432 were designed according to  (see Additional file 3). Total RNA was treated with DNAse I (Promega) to eliminate any residue of genomic DNA. Six-hundred nanograms of DNAse-treated RNA were used to generate the first strand cDNA . Oligo(dT) primer was added to the reaction for further normalization with the endogenous control gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH; ). The reaction mixture was placed at GeneAmp9700 thermocycler (Applied Biosystems) and incubated at 16°C for 30 minutes, followed by 60 cycles of pulsed reverse transcription at 30°C for 30 seconds, 42°C for 30 seconds, and 50°C for one second.
cDNA dilutions were used for PCR reactions as following: 1.0 μL of cDNA, 1.5 mM Magnesium Sulfate, 0.25 mM each dNTP, 10 pmol each primer, and 1 U of Taq DNA Polymerase (Fermentas). The reactions were placed in the thermocycler with the following conditions: 94°C for two minutes and appropriate cycle numbers of 94°C for 20 seconds, 60°C for 30 seconds, and 72°C for 45 seconds. All reactions were repeated at least three times.
Prior use of SsGAPDH (accession TC77224) as a control to evaluate miRNA accumulation in sugarcane ancient wild species, the efficiency of its primers was tested in genomic DNA from leaves of S. officinarum (accession Muntok, Java) and S. spontaneum (accession SES205A; see additional file 4). Thirty nanograms of genomic DNA were used as a template for PCR reactions. The reactions were placed in the thermocycler with the following conditions: 94°C for three minutes and 32 cycles of 94°C for 20 seconds, 58°C for 30 seconds, and 72°C for 45 seconds. The reactions were repeated twice.
Analysis of 5'RACE
Five micrograms of total RNA from sugarcane plantlets (hybrid RB 83-5486) were ligated to a RNA adapter, in a reaction mixture containing 0.5 U/μL of T4 RNA Ligase, 4 U/μL RNAse inhibitor, and 1 mM ATP. The subsequent steps were performed according to the manufacturer's guide of the GeneRacer kit (Invitrogen). The first PCR was done using the following SsCBP1 specific primer: 5'-GAAAGCCCTCTCCGCCAGC. The PCR reaction was subsequently used as a template for a semi-NESTED PCR with an internal SsCBP1 specific primer (5'-GCGCCGTCGCCGCACCC). After amplification, 5'RACE products were gel-purified and cloned, and at least 13 independent clones were randomly chosen and sequenced.
Sugarcane miRNA precursor identification
Sugarcane ESTs and GSSs were retrieved from The Gene Index Program (116,588 unique sequences; Release 2.2, July 2008) and NCBI, respectively. The sequences were used as drivers for a BLASTX search (e-value e-10) against the NCBI protein sequence database (September 2008). All potential no hit sequences were recorded as a distinct dataset. Recorded miRNAs from plants were obtained from the miRBase (over 2,300 miRNA sequences; Release 14.0)  and used as drivers for BLASTN search of sugarcane miRNA precursors in the aforementioned dataset, similarly as described by Zhang et al.. We allowed 0-3 nt mismatches or gaps between drivers and database sequences. The BLASTN parameters were adjusted to expected values of 1000 and number of descriptions and alignments of 1000. The default word-match size between the query and the database sequences was seven with a low complexity filtering ability. We also employed BLAST searches to remove sugarcane sequences similar to tRNAs, ncRNAs (http://biobases.ibch.poznan.pl/ncRNA), snoRNAs (http://bioinf.scri.sari.ac.uk/cgi-bin/plant_snorna/home) or other RNAs found in the Rfam database .
Wherever available, precursor sequences of approximately 620 nt were extracted (300 nt upstream of and 300 nt downstream from the BLAST hits) and used for hairpin structure predictions using MFOLD3.2 algorithm . Number of structures, free energy, miRNA-like helicity, number of arms per structure, size of helices within arms, and size/symmetry of internal loops within arms were analyzed by our in-house MIRcheck-based script , following manual inspection. RNA sequences were considered miRNA precursor candidates only if they fitted the following criteria: (1) the RNA sequence could fold into an appropriate stem-loop hairpin secondary structure; (2) mature miRNA site was located in one arm of the hairpin structure; (3) the mature miRNA sequence was located in the same arm of the hairpin as its homolog in other plant species; (4) mature miRNA had six or fewer, and one or more, mismatches with the miRNA* sequence in the opposite arm; (5) no break in miRNA* sequences; (6) predicted secondary structures had MFEI values higher than 0.65 [12, 19], negative MFEs, and 30-70% G + C contents; (7) two base pairs of maximum consecutive mismatches between miRNA and miRNA*; (8) a minimum of two bases pairing after the alignment between the predicted miRNA sequence and its opposite miRNA* sequence within the secondary structure; (9) and a final stem loop with a minimum of 60 nt.
Predicting sugarcane miRNA targets
Sugarcane miRNA mature sequences were used to BLAST search for possible gene targets present in the SoGI database and available BAC sequences. To minimize the number of false positives, 21-nt miRNA sequences were initially divided into three blocks of eight (block 1), three (block 2), and 10 bp (block 3). The maximum mismatches permitted in each block for the mRNA:miRNA duplex were two, zero, and three, respectively. To more thoroughly assess the mRNA::miRNA potential pairing, we additionally developed a more sensitive computational approach to identify target candidates. Each miRNA complementary site was scored, with perfect matches given a score of zero. Points were added for each G:U bulge (0.5), non-G:U mismatch (one), and bulged nucleotide in the miRNA or target strand (1.5). Only SoGI/BAC sequences that scored ≤3.5 points were further considered as potential miRNA targets. Closely related sequences were blasted against each other and analyzed. Sequences sharing ≥95% of identity at nucleotide level were considered as one gene target.
SsCBP1 comparative sequence analysis
Comparative analysis of sugarcane SsCBP1's possible orthologous in green plants was done by constructing a phylogenetic tree containing highly similar plant sequences. A BLASTX search was performed using SsCBP1 as query against a green plant protein dataset of 365,187 protein sequences obtained from several completed genomes (Arabidopsis thaliana, version 7.0 - http://www.arabidopsis.org;Populus trichocarpa, version 1.1 - http://genome.jgi-psf.org/Poptr1_1/Poptr1_1.home.html; Glycine max, version 0.1 - http://www.phytozome.net/soybean.php;Oryza sativa, version 5.0 - http://www.tigr.org/tdb/e2k1/osa1/pseudomolecules/info.shtml;Sorghum bicolor, version 1.4 - http://genome.jgi-psf.org/Sorbi1/Sorbi1.home.html; Selaginella moellendorffii, version 1.0 - http://genome.jgi-psf.org/Selmo1/Selmo1.home.html; Physcomitrella patens patens, version 1.1 - http://genome.jgi-psf.org/Phypa1_1/Phypa1_1.home.html; Volvox carteri, version 1.0 - http://genome.jgi-psf.org/Volca1/Volca1.home.html; Chlamydomonas reinhardtii, version 3.0 - http://genome.jgi-psf.org/chlre3/chlre3.home.html; Ostreococcus lucimarinus, version 2.0 - http://genome.jgi-psf.org/Ostta4/Ostta4.home.html; Ostreococcus tauri, version 2.0 - http://genome.jgi-psf.org/Ostta4/Ostta4.home.html; Micromonas pusilla CCMP1545, version 2.0 - http://genome.jgi-psf.org/MicpuC2/MicpuC2.home.html; Micromonas strain RCC299, version 2.0 - http://genome.jgi-psf.org/MicpuN2/MicpuN2.home.html). The conserved domains found among protein sequences were aligned using ClustalW  to produce ungapped alignments. The phylogenetic relationship of these aligned sequences was then constructed using the Neighbor-Joining method. Phylogenetic analysis was conducted in MEGA4 software . This process allowed identifying the most probable orthologous sequences of the SsCBP1. EST sequences from barley, wheat, and sugarcane were obtained from "TIGR Plant Transcript Assemblies Database" , and cDNA sequences from maize were obtained from MAGI (http://magi.plantgenomics.iastate.edu/). The accession numbers of genes shown in Figure 5 are as follows: SsCBP1 - TC67256; SbCBP1 - Sb02g036870; ZmCBP1 - MAGIv4.0 54669; HvCBP1 - BI947163; TaCBP1 - CK217219; OsCBP1 - Os07g38290; AtCBP1 - At5g26330; VvCBP1 - Sim4.aln-TCVV023209; PtCBP1 - 821987; GmCBP1 - Gm0010 × 00014; GmCBP2 - Gm0133 × 00019; PtCBP2 - 415490; PtCBP3 - 195948; PtCBP4 - 410618; PtCBP5 - 561943; PtCBP6 - 173259; ZmCBP2 - MAGIv4.0 158060; SbCBP2 - Sb01g004320; SmCBP1 - 27471; AtCBP2 - At2g26720; AtCBP3 - At2g31050.
The authors are indebted to Dr. Silvana Creste Souza (Instituto Agronomico de Campinas) for providing plant materials of S. officinarum and S spontaneum. This work was supported by the State of Sao Paulo Research Foundation, FAPESP (grant no. 07/58289-5) and partially by the National Council for Scientific and Technological Development, CNPq (grant no. 474635/2008-2).
- Nogueira FTS, Madi S, Chitwood DH, Juarez MT, Timmermans MCP: Two small regulatory RNAs establish opposing fates of a developmental axis. Genes Development. 2007, 21: 750-755.PubMedPubMed CentralView ArticleGoogle Scholar
- Chitwood DH, Nogueira FTS, Howell MD, Montgomery TA, Carrington JC, Timmermans MCP: Pattern formation via small RNA mobility. Genes & Development. 2009, 23: 549-554.View ArticleGoogle Scholar
- Rubio-Somoza I, Cuperus JT, Weigel D, Carrington JC: Regulation and functional specialization of small RNA-target nodes during plant development. Current Opinion in Plant Biology. 2009, 12 (5): 622-627.PubMedView ArticleGoogle Scholar
- Shukla LI, Chinnusamy V, Sunkar R: The role of microRNAs and other endogenous small RNAs in plant stress responses. Biochimica et Biophysica Acta - Gene Regulatory Mechanisms. 2008, 1779 (11): 1874-9399.View ArticleGoogle Scholar
- Ruiz-Ferrer V, Voinnet O: Roles of Plant Small RNAs in Biotic Stress Responses. Annual Review of Plant Biology. 2009, 60: 485-510.PubMedView ArticleGoogle Scholar
- Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 2005, 17: 1376-1386.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang B, Pan X, Wang Q, Cobb GP, Anderson TA: Computational identification of microRNAs and their targets. Computational Biology and Chemistry. 2006, 30 (6): 395-407.PubMedView ArticleGoogle Scholar
- Voinnet O: Origin, Biogenesis, and Activity of Plant MicroRNAs. Cell. 2009, 136 (4): 669-687.PubMedView ArticleGoogle Scholar
- Rhoades MW, Reinhardt BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520.PubMedView ArticleGoogle Scholar
- Zhao CZ, Xia H, Frazier TP, Yao YY, Bi YP, Li AQ, Li MJ, Li CS, Zhang BH, Wang XJ: Deep sequencing identifies novel and conserved microRNAs in peanuts (Arachis hypogaea L.). BMC Plant Biology. 2010, 5 (10): 3-View ArticleGoogle Scholar
- Lu C, Meyers BC, Green PJ: Construction of small RNA cDNA libraries for deep sequencing. Methods. 2007, 43 (2): 110-117.PubMedView ArticleGoogle Scholar
- Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species. BMC Plant Biology. 2008, 8: 37-PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang BH, Pan XP, Wang QL, Cobb GP, Anderson TA: Identification and characterization of new plant microRNAs using EST analysis. Cell Research. 2005, 15: 336-360.PubMedView ArticleGoogle Scholar
- Birch RG: Metabolic Engineering in Sugarcane: Assisting the Transition to a Bio-based Economy. Applications of Plant Metabolic Engineering. Springer Netherlands; 2007,249-281. 2007View ArticleGoogle Scholar
- Jannoo N, Grivet L, Chantret N, Garsmeur O, Glaszmann JC, Arruda P, D'Hont A: Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome. Plant Journal. 2007, 50: 574-585.PubMedView ArticleGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP: A highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007, 3: 12-PubMedPubMed CentralView ArticleGoogle Scholar
- Lam E, Shine J, Da Silva J, Lawton M, Bonos S, Calvino M, Carrer H, Silva-Filho MC, Glynn N, Helsel Z, Ma J, Richard E, Souza MG, Ming R: Improving sugarcane for biofuel: engineering for an even better feedstock. GCB Bioenergy. 2009, 1: 251-255.View ArticleGoogle Scholar
- Chuck G, Cigan AM, Saeteurn K, Hake S: The heterochronic maize mutant Corngrass1 results from overexpression of a tandem microRNA. Nature Genetics. 2007, 39: 544-549.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cellular and Molecular Life Sciences. 2006, 63: 246-254.PubMedView ArticleGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, Mehboob-ur-Rahman , Ware D, Westhoff P, Mayer KFX, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457 (7229): 551-556.PubMedView ArticleGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. Journal Molecular Biology. 2000, 302: 205-217.View ArticleGoogle Scholar
- Lu C, Jeong DH, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, Liu B, Cao X, Meyers BC, Green PJ: Genome-wide analysis for discovery of rice microRNAs revels natural antisense microRNAs (nat-miRNAs). Proceedings of the National Academy of Sciences. 2008, 105 (12): 4951-4956.View ArticleGoogle Scholar
- Piriyapongsa J, Jordan IK: A family of human microRNA genes from miniature inverted-repeat transposable elements. PLoS One. 2007, 2 (2): e203-PubMedPubMed CentralView ArticleGoogle Scholar
- Piriyapongsa J, Jordan IK: Dual coding of siRNAs and miRNAs by plant transposable elements. RNA. 2008, 14 (5): 814-21.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang L, Chia JM, Kumari S, Stein JC, Liu Z, Narechania A, Maher CA, Guill K, McMullen MD, Ware D: A Genome-Wide Characterization of MicroRNA Genes in Maize. PLoS Genetics. 2009, 5 (11): e1000716-PubMedPubMed CentralView ArticleGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Sung GH, Spatafora JW, Carrington JC: Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nature Genetics. 2004, 36: 1282-1290.PubMedView ArticleGoogle Scholar
- Xie F, Frazier TP, Zhang B: Identification and characterization of microRNAs and their targets in the bioenergy plant switchgrass (Panicum virgatum). Planta. 2010, 232 (2): 417-434.PubMedView ArticleGoogle Scholar
- Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu JK, Sun Q: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biology. 2007, 8 (6): R96-PubMedPubMed CentralView ArticleGoogle Scholar
- Unver T, Budak H: Conserved microRNAs and their targets in model grass species Brachypodium distachyon. Planta. 2009, 230: 659-669.PubMedView ArticleGoogle Scholar
- Song C, Fang J, Li X, Liu H, Thomas Chao C: Identification and characterization of 27 conserved microRNAs in citrus. Planta. 2009, 230: 671-685.PubMedPubMed CentralView ArticleGoogle Scholar
- Ha M, Lu J, Tian L, Ramachandran V, Kasschau KD, Chapman EJ, Carrington JC, Chen X, Wang XJ, Chen ZJ: Small RNAs serve as a genetic buffer against genomic shock in Arabidopsis interspecific hybrids and allopolyploids. Proceedings of the National Academy of Sciences. 2009, 106: 17835-17840.View ArticleGoogle Scholar
- Nogueira FTS, Chitwood DH, Madi S, Ohtsu K, Schnable PS, Scanlon MJ, Timmermans MC: Regulation of small RNA accumulation in the maize shoot apex. PLoS Genetics. 2009, 5: e1000320-PubMedPubMed CentralView ArticleGoogle Scholar
- Ovcharenko I, Boffelli D, Loots GG: eShadow: a tool for comparing closely related sequences. Genome Research. 2004, 14: 1191-1198.PubMedPubMed CentralView ArticleGoogle Scholar
- Chang WC, Lee TY, Huang HD, Huang HY, Pan RL: PlantPAN: Plant promoter analysis navigator, for identifying combinatorial cis-regulatory elements with distance constraint in plant gene groups. BMC Genomics. 2008, 26: 561-View ArticleGoogle Scholar
- Vaucheret H, Fagard M: Transcriptional gene silencing in plants: targets, inducers and regulators. Trends in Genetics. 2001, 17 (1): 29-35.PubMedView ArticleGoogle Scholar
- Warthmann N, Das S, Lanz C, Weigel D: Comparative Analysis of the MIR319a MicroRNA Locus in Arabidopsis and Related Brassicaceae. Molecular Biology and Evolution. 2008, 25 (5): 892-902.PubMedView ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Molecular Cell. 2004, 14: 787-799.PubMedView ArticleGoogle Scholar
- Llave C, Kasschau KD, Rector MA, Carrington JC: Endogenous and silencing-associated small-RNAs in plants. Plant Cell. 2002, 14: 1605-1619.PubMedPubMed CentralView ArticleGoogle Scholar
- Sunkar R, Girke T, Jain PK, Zhu JK: Cloning and characterization of microRNAs from rice. Plant Cell. 2005, 17: 1397-1411.PubMedPubMed CentralView ArticleGoogle Scholar
- Adai A, Johnson C, Mlotshwa S, Archer-Evans S, Manocha V, Vance V, Sundaresan V: Computational prediction of miRNAs in Arabidopsis thaliana. Genome Research. 2005, 15: 78-91.PubMedPubMed CentralView ArticleGoogle Scholar
- Axtell MJ, Bowman JL: Evolution of plant microRNAs and their targets. Trends Plant Sci. 2008, 13: 343-349.PubMedView ArticleGoogle Scholar
- Ding D, Zhang L, Wang H, Liu Z, Zhang Z, Zheng Y: Differential expression of miRNAs in response to salt stress in maize roots. Annals of Botany (London). 2009, 103: 29-38.View ArticleGoogle Scholar
- Vincentz M, Bandeira-Kobarg C, Gauer L, Schlögl P, Leite A: Evolutionary pattern of angiosperm bZIP factors homologous to the maize Opaque2 regulatory protein. Journal of Molecular Evolution. 2003, 56 (1): 105-116.PubMedView ArticleGoogle Scholar
- Iskandar HM, Simpson RS, Casu RE, Bonnett GD, MacLean DJ, Manners JM: Comparison of Reference Genes for Quantitative Real-Time Polymerase Chain Reaction Analysis of Gene Expression in Sugarcane. Plant Molecular Biology Reporter. 2004, 22: 325-337.View ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Research. 2006, 34: D140-D144.PubMedPubMed CentralView ArticleGoogle Scholar
- Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR, Bateman A: Rfam: Updates to the RNA families database. Nucleic Acids Research. 2009, 37: D136-140.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Research. 2003, 31 (13): 3406-3415.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Research. 1994, 22: 4673-4680.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution. 2007, 24 (8): 1596-1599.PubMedView ArticleGoogle Scholar
- Childs KL, Hamilton JP, Zhu W, Ly E, Cheung F, Wu H, Rabinowicz PD, Town CD, Buell CR, Chan AP: The TIGR Plant Transcript Assemblies database. Nucleic Acids Research. 2007, 35: D846-D851.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.