Characterization of the small RNA component of leaves and fruits from four different cucurbit species
- Guru Jagadeeswaran†1,
- Padma Nimmakayala†2,
- Yun Zheng†3,
- Kanchana Gowdu1,
- Umesh K Reddy3 and
- Ramanjulu Sunkar1Email author
© Jagadeeswaran et al.; licensee BioMed Central Ltd. 2012
Received: 1 November 2011
Accepted: 29 June 2012
Published: 23 July 2012
MicroRNAs (miRNAs) are a class of non-coding small RNAs involved in post-transcriptional regulation of gene expression critical for plant growth and development, stress responses and other diverse biological processes in plants. The Cucurbitaceae or cucurbit family represents some of economically important species, particularly those with edible and medicinal fruits. Genomic tools for the molecular analysis of members of this family are just emerging. Partial draft genome sequence became available recently for cucumber and watermelon facilitating investigation of the small RNA component of the transcriptomes in cucurbits.
We generated four small RNA libraries from bottle gourd (Lagenaria siceraria), Cucurbita moschata, Cucurbita pepo, and, watermelon (Citrullus lanatus var. lanatus) in order to identify conserved and novel lineage specific miRNAs in these cucurbits. Deep sequencing of small RNA libraries from these species resulted in 1,597,263, 532,948, 601,388, and 493,384 unique sRNA reads from bottle gourd, moschata, pepo and watermelon, respectively. Sequence analysis of these four libraries resulted in identification of 21 miRNA families that are highly conserved and 8 miRNA families that are moderately conserved in diverse dicots. We also identified 4 putative novel miRNAs in these plant species. Furthermore, the tasiRNAs were identified and their biogenesis was determined in these cucurbits. Small RNA blot analysis or q-PCR analyses of leaf and fruit tissues of these cucurbits showed differential expression of several conserved miRNAs. Interestingly, the abundance of several miRNAs in leaves and fruits of closely related C. moschata and C. pepo was also distinctly different. Target genes for the most conserved miRNAs are also predicted.
High-throughput sequencing of small RNA libraries from four cucurbit species has provided a glimpse of small RNA component in their transcriptomes. The analysis also showed considerable variation within four cucurbit species with regards to expression of individual miRNAs.
Although transcriptional gene regulation is the most important mode of gene regulation, miRNA-dependent post-transcriptional gene regulation is also absolutely necessary for completion of life cycle in plants as the disruption of either the miRNA biogenesis or function is lethal. Thus, identification of miRNAs and their targets from various plant species of agricultural importance is not only vital as part of basic biology but might also have biotechnological applications. MicroRNA genes are transcribed by RNA Polymerase II and the resulting pre-miRNA transcript can adopt a hairpin-like structure. These structures serve as substrates for Dicer-like 1 (DCL1), an endoribonuclease that precisely severs the hairpin-like structure and releases 21-nt long miRNA:miRNA* duplex with 2-nt overhangs at the 3′ends. Although DCL1 is a major determinant for processing of miRNAs, several other proteins such as Hyponastic Leaves1 (HYL1), Serrate (SE), Dawdle (DDL) and Cap-binding protein (CBP) are also important for processing of miRNAs from its hairpin-like precursor ([1–6]. The 21 nt miRNA:miRNA* duplex is methylated at the 3′ends by the Hua Enhancer 1 (HEN1). The duplex is then exported to the cytosol, where it can be loaded into RNA-Induced Silencing Complex (RISC). The miRNA in the RISC serves as guide molecule to identify its target mRNA, which will be either degraded or prevented from being translated [7–10].
Approximately 21 miRNA families are known as highly conserved among the angiosperms (dicots and monocots) and 8 and 6 of these miRNA families conserved even in gymnosperms and bryophytes, respectively [11, 12]. Because of the sequence conservation, identification of conserved miRNAs is not a difficult task provided sufficient genomic/EST resources are available for a plant species in question [13, 14]. Even in the absence of genomic/EST resources, conserved miRNAs can be identified by expression analysis using probes designed to detect conserved miRNAs. Besides conserved miRNAs, a large number of miRNA families have been identified in diverse plant species, some of which are conserved in closely related species or species-specific, which are referred to as ‘young miRNAs’ . Only a minor portion of these young miRNAs appears to be functional, whereas the majority of them appear to be non-functional and eventually dissipate, and only few of them will be integrating into gene regulatory networks [12, 15]. Although functionally such non-conserved miRNAs seem less significant, their identification is important to trace the birth and death of miRNAs in a specific plant lineage .
MicroRNAs have been reported from several model or crop plant species [16–23], but little is known about the miRNA component in the family of Cucurbitaceae. Cucurbitaceae is known to have 90 genera and 700 species, out of which there are several domesticated species for food that includes Citrullus lanatus (watermelon), Cucumis sativus (cucumber), Cucumis melo (melon), Lagenaria siceraria (bottle gourd) and Cucurbita (squash & pumpkin) [24–27]. Little is known about the small RNA component of the transcriptome in a Cucurbitaceae member. Only recently, miRNAs in cucumber (Cucumis sativus) and melon (Cucumis melo) have been analyzed using relatively small number of reads [28, 29]. Here we report high-throughput sequencing of small RNA libraries from the four different cucurbit species. The analysis identified 21 highly conserved miRNA families, and 8 miRNA families that are only moderately conserved in various dicots. Analysis of small RNA libraries from four cucurbit species facilitated identification of 4 putative novel miRNAs based on their recovery in at least two different cucurbit genera. Cucurbit small RNAs with similarity to Arabidopsis TAS3 siRNAs were identified and their biogenesis was also determined. Furthermore, we predicted target genes for the most conserved miRNAs with the available watermelon, pepo and melon genomes.
Results and discussion
Sequence analysis and annotation
Small RNA libraries in different cucurbit species
With regards to unique reads, the distribution of size classes was uniform across the species, where 24-nt size class and 23-nt size class ranked first and second respectively in all 4 libraries (Figure 1). Small RNAs of 20-21-nt size class represents largely miRNAs. Recently, Martinez et al. reported small RNAs from cucumber (Cucumis sativus belonging to Cucurbitaceae family) by analyzing 38,787 unique small RNAs , which is far less compared to the current study in which several millions of small RNAs from 4 different cucurbits were analyzed. Size distribution of small RNAs in their library was typical of a plant small RNA library with two peaks one each at 21-nt and at 24-nt. Martinez et al. constructed small RNA libraries from the leaves and phloem exudates , whereas we used pooled RNA from leaves, stems and three different fruit tissues. Mix of diverse tissues used in this study, particularly the inclusion of fruit tissues might have contributed to the variation in size distribution including the observed peak at 19-nt size class and deserves further study.
Identification of conserved miRNAs in different cucurbit species
Based on their functional analysis, conserved miRNAs appear to be involved in almost all aspects of plant growth and development , as well as biotic and abiotic stress responses [36–38]. Our analysis of small RNAs has identified all 21 conserved miRNA families in four different cucurbits as expected (Additional file 1). Additionally, miR894 and miR2111 were recovered in all 4 cucurbits, whereas miR158, miR824, miR827, miR858, miR2916 and miR2950 that are moderately conserved in some of the dicots could be identified only in watermelon (Additional file 1).
Many conserved miRNA families have multiple loci and the members often varied by one or two nucleotides. Most miRNA families were represented by multiple members; miR156 (5 members), miR159 (7 members), miR164 (6 members), miR165/miR166 (13 members), miR167 (7 members), miR169 (9 members), miR170/171 (11 members), miR172 (9 members), miR396 (7 members) and miR319 (6 members) (Additional file 1). However, 10 conserved miRNA families, miR160 (2 members, miR162 (1 member), miR168 (2 members), miR393 (2 members), miR390 (2 members), miR394 (1 member), miR395 (2 members), miR397 (1 member), miR398 (2 members) and miR408 (2 members) were represented by one or two members only. Overall, we recovered 83, 82, 82 and 81 miRNAs belonging to 22, 21, 22 and 29 miRNA families in bottle gourd, moschata, pepo and watermelon respectively. The highest number of miRNA families in watermelon is because of identification of less conserved miR824, miR827, miR858, miR1515, miR2916 and miR2950.
The number of miRNA variants within the miRNA families was similar in the most cases. However, we also observed preferential expression of specific members within the miRNA families. For instance, miR168 represented by two members (miR168a and 168b) that differ by one nucleotide and only one of these two loci preferentially expressed in all 4 cucurbit species (Additional file 1). This was previously shown both in Arabidopsis and rice [44, 45]. In addition, we also found variation within the cucurbits in this study. For example, only one variant of miR396 was expressed at high level in watermelon (Additional file 1).
Few interesting observations emerged from the analysis of small RNA libraries are worth mentioning here. One such observation was the retrieval of miR158 homolog from watermelon but not from the other 3 cucurbits (Additional file 1). miR158 has been previously considered to be specific to Brassicaceae family . In corroboration with our finding, miR158 was also identified in (Cucumis melo) small RNA libraries . Another interesting observation was identification of miR894 homologs in these cucurbits. miR894 was reported from Physcomitrella, and in a few dicots such as sorghum  tomato  and cotton . Although miR894 homologs were found in all four different cucurbits, watermelon has the highest frequency (Additional file 1). Identification of miR894 in all 4 cucurbits and its absence in several other dicots raises several interesting evolutionary questions including whether homologs of miR894 were selectively lost in certain lineages.
Distinct miRNA profiles in four cucurbits
Normalized miRNA abundance indicated that several miRNA families (miR156, miR159 and miR160, miR164, miR167 and miR172) are expressed at high levels in watermelon compared to 3 other cucurbits. The differential abundance was strikingly high in case of miR156, where nearly 50 to 100 fold greater expression was noticed in watermelon as compared to 3 other cucurbits (Figure 2a). Similarly, expression of miR167 and miR172 varied considerably among the different cucurbits. Expression level of miR167 in watermelon was nearly 12-fold and 18-fold greater than in bottle gourd and Cucurbita genus, respectively (Figure 2a; Additional file 1). A similar trend was observed for miR172 (Figure 2a). In general, the miRNA abundance of most conserved miRNAs was the highest in watermelon and the least in bottle gourd (Additional file 1). However, a few miRNAs (miR162, miR168, miR397 and miR408) showed higher levels in bottle gourd compared to the 3 other cucurbit species studied here (Additional file 1). Of these, expression of miR168 in bottle gourd was nearly 10-fold higher when compared to moschata and pepo and 5-fold higher compared to watermelon (Figure 2a).
Overall, miRNA profiles are remarkably similar between moschata and pepo (Figure 2a; Additional file 1). Identical miRNA profiles between moschata and pepo is not surprising given the fact that these two belong to the same genus Cucurbita. Technically, this also confirms the uniform pooling of RNA samples from different tissues for library construction. The read count of miR166 family accounts for nearly 75% of the total reads, thus, is the most highly expressed in bottle gourd, moschata and pepo, whereas miR167 is the most highly expressed in watermelon. The second most abundant miRNA family was miR159 in moschata and pepo, miR164 in bottle gourd and miR156 in watermelon (Figure 2, Additional file 1). Low expression levels of miR162 and miR393 has been reported in diverse plant species [50, 51]. Overall, our data demonstrated that considerable variabilty exists within the four species of Cucurbitaceae with regard to abundance of miRNA families.
Dynamic regulation of conserved miRNAs in cucurbit tissues
Differential miRNA abundance in leaf and fruit tissues of moschata and pepo
C. moschata and C. pepo are two closely related species in Cucurbitacea family that show a high level of macro-synteny in their genomes . In order to examine miRNA expression profiles in leaves and fruit tissues of these two species, we carried out blot analysis for several conserved miRNAs. Expression pattern of nine miRNAs analyzed in these two tissues had a similar profile in both species although some differences were noticed (Figure 4). For instance, miR166, miR167 and miR168 showed almost similar abundance in leaves and fruits of moschata, but their levels differed between leaves and fruits of pepo. miR166 and miR168 abundance was higher in leaves whereas miR167 abundance was higher in fruits of pepo. Expression level of miR159 was particularly high in leaves of both moschata and pepo. By contrast, miR156, miR171 and miR172 expression levels were relatively higher in fruits of both species (Figure 4). miR156 is generally considered to be an abundantly expressed miRNA in leaves. However, it has recently been noted that expression of miR156 increases during fruit ripening in tomato as observed here . Abundance of miR396 was similar in leaves and fruits of pepo but slightly higher in fruits of moschata. Of the different miRNAs studied, miR169 had the least expression in both tissues and only a faint signal was detected, which is again consistent with the normalized read count (Figure 4 and Additional file 1). Expression level of miRNAs determined from q-PCR and small RNA blot analysis did not correlate in certain cases and such differences between different methodologies have been noticed earlier when measuring the miRNA abundance [57–59]. The observed differences could also be due to the differential sensitivity of each of these two methods’ i.e., q-PCR measures the abundance of individual miRNA, whereas small RNA blot analysis measures the abundance of multiple miRNAs belonging to a specific miRNA family.
Putative novel miRNAs in cucurbits
Normalized abundance (TPM) of putative novel miRNAs and their predicted targets in cucurbits
Number of normalized reads
PU022623 (ubiquitin-protein ligase) ;
WMU58406 (acetolactate synthase)
PU061865 (unknown protein) ;
WMU41511 (unknown protein)
PU134802 (unknown protein)
PU119447 (unknown protein) ;
PU062878 (unknown protein) ;
PU005456 (zinc finger-like superfamily protein)
PU036654 (glycine rich protein);
PU014282 (unknown protein)
PU057169 (glyceraldehyde-3-phosphate dehydrogenase)
PU043518 (unknown protein)
PU084326 (cytochrome c oxidase)
WMU42945 (cytochrome c oxidase)
Identification of TASi locus and tasiRNAs in cucurbits
Endogenous siRNAs comprise trans-acting siRNAs (tasiRNAs), natural anti-sense transcript-derived siRNAs (NATsiRNAs) and repeat-associated siRNAs (rasiRNAs) in plants . Accurate annotation of siRNAs requires the knowledge of complete genomes. TasiRNAs are the only siRNAs that are conserved among diverse plant species and therefore these can be identified easily. Transcript from TAS3 locus harbors a tandem repeat of 21-nt almost identical sequences and possesses dual miR390 complimentary sequences. miR390 guided cleavage of TAS3 transcript at the 3′ end leads to production of dsRNA which is further processed by DCL4 to release tasiRNAs .
miRNA targets in cucurbit transcriptome
Potential targets for conserved miRNAs in cucurbits
Targets showing highest homology
Squamosa promoter-binding protein-like (SPL) proteins
MYB-like binding factors
Auxin response factor (ARFs)
NAC domain protein
Scarecrow-like (SCL1) and GRAS family transcription factors
PU002536,PU074309, PU067878,PU002536, PU115341
TCP/DNA binding proteins
Auxin receptors and BLH transcription factors
Putative blue copper binding protein
Plastocyanin-like domain-containing protein
Comparative miRNA profiling in four cucurbits led to the identification of conserved (highly conserved as well as somewhat less conserved) miRNAs in cucurbits. Additionally, 4 putative new miRNAs are identified in these four cucurbits. Expression analysis showed differential regulation of several conserved miRNAs between leaves and fruit tissues of cucurbits. Our analysis also demonstrates considerable variability within four cucurbits with regards to expression of individual miRNAs. Even more strikingly, several miRNAs expression patterns differed between C. moschata and C. pepo, the two closely related cucurbits analyzed in this study. The predicted targets for conserved miRNAs suggested the involvement of miRNAs in regulating growth and development as well as other important physiological processes in cucurbits.
Small RNA library construction and sequence analysis
Total RNA was isolated from different tissues (leaf, stem and flesh, rind and placenta of the fruits) using TRIzol reagent and small RNA libraries were generated from four cucurbit species: bottle gourd (Lagenaria siceraria (accession Grif 1617 collection from India)), Cucurbita moschata (accession Grif 14244 Early Butternut) Cucurbita pepo (accession NSL98075 Table King), and watermelon (Citrullus lanatus var. lanatus) (PI 438676 Charleston Grey) by pooling equimolar amounts of total RNA from the aforesaid tissues. Construction of small RNA libraries from size fractionated RNA was carried out as described previously . In brief, small RNA fractions of 18–28 nt were isolated from 15% denaturing polyacrylamide gels and sequentially ligated to 5′ and 3′ RNA adapters. Small RNAs ligated with adapters were converted to DNA by RT-PCR following Solexa protocol. The final PCR product was gel purified and sequenced by Genome Analyser II (Illumina).
Analysis of sequence reads from four cucurbit libraries was performed as described previously [13, 16]. Briefly, adaptor sequences of the raw reads were removed and the small RNAs in between the adaptors were extracted. Unique small RNAs were obtained after eliminating redundant sequences. Reads mapped to known non-coding RNAs (rRNAs, tRNAs, snRNAs, snoRNAs) and repeats were removed from unique RNAs by aligning to databases Rfam  and Repbase , and the Plant Repeats database . Remaining small RNAs were searched in the miRBase database  (release 16, obtained from http://microrna.sanger.ac.uk/sequences/ftp.shtml) to identify conserved microRNAs. Small RNAs mapping to known miRNAs of other plant species were designated conserved miRNAs in cucurbits. The genome sequence of cucumber was obtained from the Cucumber Genome DataBase . To identify novel miRNAs, unique small RNAs with more than 10 genomic hits in cucumber genome were removed because they might possibly come from repeat-rich loci. The potential candidate miRNAs were identified by folding the flanking genome sequence of unique small RNAs using the RNAfold program . To identify TAS genes and tasiRNAs in cucurbits, the genomes of watermelon and pepo genomes were searched using Hitsensor scores for possible miR390 binding sites  separated by 200 to 300-nt.
Bioinformatic prediction of miRNA targets
To predict potential targets for cucurbit miRNAs, partially annotated coding sequences (unigenes) in watermelon, pepo and melon genomes were used for searching sequences complementary to the miRNAs (Cucurbit Genomics Database: http://www.icugi.org) . A scoring matrix that allows a maximum of 3.5 mismatches (G:U accounts for half mismatch) between the miRNA and its target mRNA was used in the analysis . In order to get putative annotation for the target genes, sequences identified as targets were BLAST searched and the genes showing highest homology in other plants were assigned as putative annotation in cucurbits.
Expression profile of miRNAs in cucurbits was verified by performing stem loop RT-PCR as described previously . One microgram total RNA in leaf and fruit tissue was used to perform RT reaction with miRNA specific RT primers and 1:10 diluted cDNA was used in q-PCR experiments. RT stem-loop primers and all other primers used are given in Additional file 3. Real-Time PCR was carried out using MaximaTM SYBR Green q-PCR Master Mix in an ABI7500 Real-Time PCR System. The relative expression was obtained using delta-Ct method and actin as reference gene. The data shown was mean of 3 replicates. Student t-test was used to determine the significant differences between different tissues as well as in the same tissue of different cucurbits.
Small RNA blot analysis
Total RNA (100 μg) was resolved on a 15% polyacrylamide gel containing 7 M urea in TBE buffer. Size-fractionated small RNAs were then transferred to Hybond-N + (Amersham) membranes, UV cross-linked and baked for 1 h at 80°C. After pre-hybridization, blots were hybridized at 38°C with a 32P-labeled DNA oligo probe complementary to miRNA sequence. After washing, blots were exposed to phosphor screen and scanned using a phosphorimager. Blots were used for reprobing after thorough stripping.
The authors are grateful to Dr. Jarret, Plant Genetic Resources Conservation Unit, USDA-ARS, Griffin, GA, 30223 for providing the seeds of germplasm accessions. We thank the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences (IVF-CAAS) for sharing the genome sequence of cucumber. This research was supported by the USDA-NIFA Research grant (2007-38814-18472) to UKR and RS and a start-up grant of Fudan University and Grant No. 10ZR1403000 of STCSM to YZ. This work also has been supported, in part, by the NSF EPSCoR award EPS0814361 to RS.
- Dong Z, Han MH, Fedoroff N: The RNA-binding proteins HYL1 and SE promote accurate in vitro processing of pri-miRNA by DCL1. Proc Natl Acad Sci U S A. 2008, 105 (29): 9970-9975. 10.1073/pnas.0803356105.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim S, Yang JY, Xu J, Jang IC, Prigge MJ, Chua NH: Two cap-binding proteins CBP20 and CBP80 are involved in processing primary MicroRNAs. Plant Cell Physiol. 2008, 49 (11): 1634-1644. 10.1093/pcp/pcn146.PubMed CentralView ArticlePubMedGoogle Scholar
- Lobbes D, Rallapalli G, Schmidt DD, Martin C, Clarke J: SERRATE: a new player on the plant microRNA scene. EMBO Rep. 2006, 7 (10): 1052-1058. 10.1038/sj.embor.7400806.PubMed CentralView ArticlePubMedGoogle Scholar
- Song L, Axtell MJ, Fedoroff NV: RNA secondary structural determinants of miRNA precursor processing in Arabidopsis. Curr Biol. 2010, 20 (1): 37-41. 10.1016/j.cub.2009.10.076.View ArticlePubMedGoogle Scholar
- Yang L, Liu Z, Lu F, Dong A, Huang H: SERRATE is a novel nuclear regulator in primary microRNA processing in Arabidopsis. Plant J. 2006, 47 (6): 841-850. 10.1111/j.1365-313X.2006.02835.x.View ArticlePubMedGoogle Scholar
- Yu B, Bi L, Zheng B, Ji L, Chevalier D, Agarwal M, Ramachandran V, Li W, Lagrange T, Walker JC, Chen X: The FHA domain proteins DAWDLE in Arabidopsis and SNIP1 in humans act in small RNA biogenesis. Proc Natl Acad Sci U S A. 2008, 105 (29): 10073-10078. 10.1073/pnas.0804218105.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.View ArticlePubMedGoogle Scholar
- Ramachandran V, Chen X: Small RNA metabolism in Arabidopsis. Trends Plant Sci. 2008, 13 (7): 368-374. 10.1016/j.tplants.2008.03.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.View ArticlePubMedGoogle Scholar
- Axtell MJ, Bartel DP: Antiquity of microRNAs and their targets in land plants. Plant Cell. 2005, 17 (6): 1658-1673. 10.1105/tpc.105.032185.PubMed CentralView ArticlePubMedGoogle Scholar
- Cuperus JT, Fahlgren N, Carrington JC: Evolution and functional diversification of MIRNA genes. Plant Cell. 2011, 23 (2): 431-442. 10.1105/tpc.110.082784.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species. BMC Plant Biol. 2008, 8: 37-10.1186/1471-2229-8-37.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang B, Pan X, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes. Plant J. 2006, 46 (2): 243-259. 10.1111/j.1365-313X.2006.02697.x.View ArticlePubMedGoogle Scholar
- Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One. 2007, 2 (2): e219-10.1371/journal.pone.0000219.PubMed CentralView ArticlePubMedGoogle Scholar
- Jagadeeswaran G, Zheng Y, Li YF, Shukla LI, Matts J, Hoyt P, Macmil SL, Wiley GB, Roe BA, Zhang W, Sunkar R: Cloning and characterization of small RNAs from Medicago truncatula reveals four novel legume-specific microRNA families. New Phytol. 2009, 184 (1): 85-98. 10.1111/j.1469-8137.2009.02915.x.View ArticlePubMedGoogle Scholar
- Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome. Science. 2005, 309 (5740): 1567-1569. 10.1126/science.1114112.View ArticlePubMedGoogle Scholar
- Subramanian S, Fu Y, Sunkar R, Barbazuk WB, Zhu JK, Yu O: Novel and nodulation-regulated microRNAs in soybean roots. BMC Genomics. 2008, 9: 160-10.1186/1471-2164-9-160.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunkar R, Girke T, Jain PK, Zhu JK: Cloning and characterization of microRNAs from rice. Plant Cell. 2005, 17 (5): 1397-1411. 10.1105/tpc.105.031682.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC plant biology. 2008, 8: 25-10.1186/1471-2229-8-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16 (8): 2001-2019. 10.1105/tpc.104.022830.PubMed CentralView ArticlePubMedGoogle Scholar
- Szittya G, Moxon S, Santos DM, Jing R, Fevereiro MP, Moulton V, Dalmay T: High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families. BMC Genomics. 2008, 9: 593-10.1186/1471-2164-9-593.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465. 10.1101/gr.075572.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeffrey C: Further notes onCucurbitaceae: III. Some African taxa. Kew Bull. 1975, 30: 475-493. 10.2307/4103075.View ArticleGoogle Scholar
- Jeffrey C: A review of the Cucurbitaceae. Bot J Linn Soc. 1980, 81: 233-247. 10.1111/j.1095-8339.1980.tb01676.x.View ArticleGoogle Scholar
- Jeffrey C: Appendix: an outline classification of the Cucurbitaceae. Biology and Utilization of the Cucurbitaceae. Edited by: Bates DM, Robinson RW, Jeffrey C. 1990, Cornell Univ. Press, Ithaca, 449-463.Google Scholar
- Schaefer H, Heibl C, Renner SS: Gourds afloat: a dated phylogeny reveals an Asian origin of the gourd family (Cucurbitaceae) and numerous oversea dispersal events. Proc Biol Sci. 2009, 276 (1658): 843-851. 10.1098/rspb.2008.1447.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzalez-Ibeas D, Blanca J, Donaire L, Saladie M, Mascarell-Creus A, Cano-Delgado A, Garcia-Mas J, Llave C, Aranda MA: Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing. BMC Genomics. 2011, 12: 393-10.1186/1471-2164-12-393.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez G, Forment J, Llave C, Pallas V, Gomez G: High-throughput sequencing, characterization and detection of new and conserved cucumber miRNAs. PLoS One. 2011, 6 (5): e19523-10.1371/journal.pone.0019523.PubMed CentralView ArticlePubMedGoogle Scholar
- Matts J, Jagadeeswaran G, Roe BA, Sunkar R: Identification of microRNAs and their targets in switchgrass, a model biofuel plant species. J Plant Physiol. 2010, 167 (11): 896-904. 10.1016/j.jplph.2010.02.001.View ArticlePubMedGoogle Scholar
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.PubMed CentralView ArticlePubMedGoogle Scholar
- Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, Burgyan J: Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis. Plant J. 2010, 62 (6): 960-976.PubMedGoogle Scholar
- Hsieh LC, Lin SI, Shih AC, Chen JW, Lin WY, Tseng CY, Li WH, Chiou TJ: Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol. 2009, 151 (4): 2120-2132. 10.1104/pp.109.147280.PubMed CentralView ArticlePubMedGoogle Scholar
- Cole C, Sobala A, Lu C, Thatcher SR, Bowman A, Brown JW, Green PJ, Barton GJ, Hutvagner G: Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs. RNA. 2009, 15 (12): 2147-2160. 10.1261/rna.1738409.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee YS, Shibata Y, Malhotra A, Dutta A: A novel class of small RNAs: tRNA-derived RNA fragments (tRFs). Genes Dev. 2009, 23 (22): 2639-2649. 10.1101/gad.1837609.PubMed CentralView ArticlePubMedGoogle Scholar
- Katiyar-Agarwal S, Jin H: Role of small RNAs in host-microbe interactions. Annu Rev Phytopathol. 2010, 48: 225-246. 10.1146/annurev-phyto-073009-114457.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunkar R: MicroRNAs with macro-effects on plant stress responses. Semin Cell Dev Biol. 2010, 21 (8): 805-811. 10.1016/j.semcdb.2010.04.001.View ArticlePubMedGoogle Scholar
- Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12 (7): 301-309. 10.1016/j.tplants.2007.05.001.View ArticlePubMedGoogle Scholar
- Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su CL: Regulation of phosphate homeostasis by MicroRNA in Arabidopsis. Plant Cell. 2006, 18 (2): 412-421. 10.1105/tpc.105.038943.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK: A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol. 2005, 15 (22): 2038-2043. 10.1016/j.cub.2005.10.016.View ArticlePubMedGoogle Scholar
- Kawashima CG, Yoshimoto N, Maruyama-Nakashita A, Tsuchiya YN, Saito K, Takahashi H, Dalmay T: Sulphur starvation induces the expression of microRNA-395 and one of its target genes but in different cell types. Plant J. 2009, 57 (2): 313-321. 10.1111/j.1365-313X.2008.03690.x.View ArticlePubMedGoogle Scholar
- Abdel-Ghany SE, Pilon M: MicroRNA-mediated systemic down-regulation of copper protein expression in response to low copper availability in Arabidopsis. J Biol Chem. 2008, 283 (23): 15932-15945. 10.1074/jbc.M801406200.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamasaki H, Abdel-Ghany SE, Cohu CM, Kobayashi Y, Shikanai T, Pilon M: Regulation of copper homeostasis by micro-RNA in Arabidopsis. J Biol Chem. 2007, 282 (22): 16369-16378. 10.1074/jbc.M700138200.View ArticlePubMedGoogle Scholar
- Vaucheret H: AGO1 homeostasis involves differential production of 21-nt and 22-nt miR168 species by MIR168a and MIR168b. PLoS One. 2009, 4 (7): e6442-10.1371/journal.pone.0006442.PubMed CentralView ArticlePubMedGoogle Scholar
- Vaucheret H, Mallory AC, Bartel DP: AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell. 2006, 22 (1): 129-136. 10.1016/j.molcel.2006.03.011.PubMed CentralView ArticlePubMedGoogle Scholar
- Fattash I, Voss B, Reski R, Hess WR, Frank W: Evidence for the rapid expansion of microRNA-mediated regulation in early land plant evolution. BMC Plant Biol. 2007, 7: 13-10.1186/1471-2229-7-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang L, Zheng Y, Jagadeeswaran G, Li Y, Gowdu K, Sunkar R: Identification and temporal expression analysis of conserved and novel microRNAs in Sorghum. Genomics. 2011, 98 (6): 460-468. 10.1016/j.ygeno.2011.08.005.View ArticlePubMedGoogle Scholar
- Moxon S, Jing R, Szittya G, Schwach F, Rusholme Pilcher RL, Moulton V, Dalmay T: Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res. 2008, 18 (10): 1602-1609. 10.1101/gr.080127.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Kwak PB, Wang QQ, Chen XS, Qiu CX, Yang ZM: Enrichment of a set of microRNAs during the cotton fiber development. BMC Genomics. 2009, 10: 457-10.1186/1471-2164-10-457.PubMed CentralView ArticlePubMedGoogle Scholar
- Hirsch J, Lefort V, Vankersschaver M, Boualem A, Lucas A, Thermes C, d’Aubenton-Carafa Y, Crespi M: Characterization of 43 non-protein-coding mRNA genes in Arabidopsis, including the MIR162a-derived transcripts. Plant Physiol. 2006, 140 (4): 1192-1204. 10.1104/pp.105.073817.PubMed CentralView ArticlePubMedGoogle Scholar
- Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JD: A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006, 312 (5772): 436-439. 10.1126/science.1126088.View ArticlePubMedGoogle Scholar
- Xue LJ, Zhang JJ, Xue HW: Characterization and expression profiles of miRNAs in rice seeds. Nucleic Acids Res. 2009, 37 (3): 916-930. 10.1093/nar/gkn998.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin RC, Liu PP, Goloviznina NA, Nonogaki H: microRNA, seeds, and Darwin?: diverse function of miRNA in seed biology and plant responses to stress. J Exp Bot. 2010, 61 (9): 2229-2234. 10.1093/jxb/erq063.View ArticlePubMedGoogle Scholar
- Nonogaki H: MicroRNA gene regulation cascades during early stages of plant development. Plant Cell Physiol. 2010, 51 (11): 1840-1846. 10.1093/pcp/pcq154.View ArticlePubMedGoogle Scholar
- Gong L, Pachner M, Kalai K, Lelley T: SSR-based genetic linkage map of Cucurbita moschata and its synteny with Cucurbita pepo. Genome. 2008, 51 (11): 878-887. 10.1139/G08-072.View ArticlePubMedGoogle Scholar
- Mohorianu I, Schwach F, Jing R, Lopez-Gomollon S, Moxon S, Szittya G, Sorefan K, Moulton V, Dalmay T: Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J. 2011, 67 (2): 232-246. 10.1111/j.1365-313X.2011.04586.x.View ArticlePubMedGoogle Scholar
- Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C: Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010, 16 (5): 991-1006. 10.1261/rna.1947110.PubMed CentralView ArticlePubMedGoogle Scholar
- Jayaprakash AD, Jabado O, Brown BD, Sachidanandam R: Identification and remediation of biases in the activity of RNA ligases in small-RNA deep sequencing. Nucleic Acids Res. 2011, 39 (21): e141-10.1093/nar/gkr693.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang W, Gao S, Zhou X, Chellappan P, Chen Z, Zhou X, Zhang X, Fromuth N, Coutino G, Coffey M, Jin H: Bacteria-responsive microRNAs regulate plant innate immunity by modulating plant hormone networks. Plant Molecular Biology. 2011, 75 (1–2): 93-105.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.PubMed CentralView ArticlePubMedGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121 (2): 207-221. 10.1016/j.cell.2005.04.004.View ArticlePubMedGoogle Scholar
- Cucurbit Genomics Database.http://www.icugi.org,
- Plant Repeats Database.http://plantrepeats.plantbiology.msu.edu/downloads.html,
- miRBase (release 16).http://microrna.sanger.ac.uk/sequences/ftp.shtml,
- Hofacker IL: Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31 (13): 3429-3431. 10.1093/nar/gkg599.PubMed CentralView ArticlePubMedGoogle Scholar
- Zheng Y, Zhang W: Animal microRNA target prediction by incorporating diverse sequence-specific determinants. J Bioinform and Comput Biol. 2010, 8: 763-788. 10.1142/S0219720010004896.View ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.View ArticlePubMedGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP: Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007, 3: 12-10.1186/1746-4811-3-12.PubMed CentralView ArticlePubMedGoogle Scholar