A radish root sRNA library sequenced on the Illumina/Solexa platform produced 15,983,864 raw reads, of which 15,779,290 (98.72%) were clean reads spanning 4,615,663 unique sequences. After filtering rRNAs, tRNAs, snRNAs, and snoRNAs, 42,308 unique reads (0.92%) matched known miRNAs while 96.13% remained unannotated and available for novel-miRNA discovery .
Using miRBase 18.0 alignment (maximum two mismatches), the authors identified 545 conserved miRNA families. Abundance varied over four orders of magnitude: miR408 (413,125 reads), miR157 (204,590), miR166 (187,287), miR168 (149,489), miR158 (119,374), and miR156 (108,477) dominated, whereas miR1437 had only one read. Secondary structures and negative folding free energies (-18.4 to -86.7 kcal/mol, mean -44.7) supported 15 novel miRNAs with miRNA* strands plus 64 potentially novel miRNAs .
qRT-PCR of four conserved (miR156, miR166, miR172, miR396) and six novel (rsa-miRn1, n9, n13, n19, n25, n33) miRNAs across six tissues broadly agreed with sequencing and showed tissue-preferential expression, e.g., rsa-miRn1/n19 favoring root and stem; rsa-miRn9 favoring flower and leaf . Target prediction against limited radish unigenes yielded 196 targets for 42 novel miRNAs (mean 4.67 targets/miRNA), enriching transcription factors (AP2, SBP, ARF, zinc finger) and stress-related proteins; 33 novel miRNAs had no predicted targets .
Strengths: first genome-wide radish miRNA inventory; adherence to Meyers 2008 annotation criteria (miRNA* evidence); explicit use of Mireap/Mfold; qRT-PCR confirmation. Key limitations acknowledged by the authors: the radish genome was not yet sequenced, so no genomic loci could be verified for conserved families and novel predictions relied on sparse EST/GSS; only a single tissue (root) was sequenced while expression was profiled across six tissues; targets were purely computational with no degradome or 5'RACE validation. The unannotated fraction (96.13%) also suggests many putative "novel" miRNAs could be siRNAs, degradation fragments, or repeats β a known confound in small-RNA discovery without genomic anchoring . Subsequent small-RNA tooling work also emphasizes that adapter-ligation bias and non-miRNA small-RNA species can inflate novel-miRNA calls in such catalogs .
Confidence note: The headline counts are well-supported by the reported data and internal consistency checks (Table 2 sums, qRT-PCR correlation), but the "novel" miRNA count carries more uncertainty than the conserved family count because loci could not be independently confirmed. The target-gene findings are the weakest layer β predicted only, never experimentally validated.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.