Comparative analysis of non-autonomous effects of tasiRNAs and miRNAs in Arabidopsis thaliana

In plants, small interfering RNAs (siRNAs) can trigger a silencing signal that may spread within a tissue to adjacent cells or even systemically to other organs. Movement of the signal is initially limited to a few cells, but in some cases the signal can be amplified and travel over larger distances. How far silencing initiated by other classes of plant small RNAs (sRNAs) than siRNAs can extend has been less clear. Using a system based on the silencing of the CH42 gene, we have tracked the mobility of silencing signals initiated in phloem companion cells by artificial microRNAs (miRNA) and trans-acting siRNA (tasiRNA) that have the same primary sequence. In this system, both the ta-siRNA and the miRNA act at a distance. Non-autonomous effects of the miRNA can be triggered by several different miRNA precursors deployed as backbones. While the tasiRNA also acts non-autonomously, it has a much greater range than the miRNA or hairpin-derived siRNAs directed against CH42, indicating that biogenesis can determine the non-autonomous effects of sRNAs. In agreement with this hypothesis, the silencing signals initiated by different sRNAs differ in their genetic requirements.     

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.     

MicroRNAs in the shoot apical meristem of soybean

Plant microRNAs (miRNAs) play crucial regulatory roles in various developmental processes. In this study, we characterize the miRNA profile of the shoot apical meristem (SAM) of an important legume crop, soybean, by integrating high-throughput sequencing data with miRNA microarray analysis. A total of 8423 non-redundant sRNAs were obtained from two libraries derived from micro-dissected SAM or mature leaf tissue. Sequence analysis allowed the identification of 32 conserved miRNA families as well as 8 putative novel miRNAs. Subsequent miRNA profiling with microarrays verified the expression of the majority of these conserved and novel miRNAs. It is noteworthy that several miRNAs* were expressed at a level similar to or higher than their corresponding mature miRNAs in SAM or mature leaf, suggesting a possible biological function for the star species. In situ hybridization analysis revealed a distinct spatial localization pattern for a conserved miRNA, miR166, and its star speciessuggesting that they serve different roles in regulating leaf development. Furthermore, localization studies showed that a novel soybean miRNA, miR4422a, was nuclear-localized. This study also indicated a novel expression pattern of miR390 in soybean. Our approach identified potential key regulators and provided vital spatial information towards understanding the regulatory circuits in the SAM of soybean during shoot development.     

Global alteration of microRNAs and transposon-derived small RNAs in cotton (Gossypium hirsutum) during Cotton leafroll dwarf polerovirus (CLRDV) infection

Small RNAs (sRNAs) are a class of non-coding RNAs ranging from 20- to 40-nucleotides (nts) that are present in most eukaryotic organisms. In plants, sRNAs are involved in the regulation of development, the maintenance of genome stability and the antiviral response. Viruses, however, can interfere with and exploit the silencing-based regulatory networks, causing the deregulation of sRNAs, including small interfering RNAs (siRNAs) and microRNAs (miRNAs). To understand the impact of viral infection on the plant sRNA pathway, we deep sequenced the sRNAs in cotton leaves infected with Cotton leafroll dwarf virus (CLRDV), which is a member of the economically important virus family Luteoviridae. A total of 60 putative conserved cotton miRNAs were identified, including 19 new miRNA families that had not been previously described in cotton. Some of these miRNAs were clearly misregulated during viral infection, and their possible role in symptom development and disease progression is discussed. Furthermore, we found that the 24-nt heterochromatin-associated siRNAs were quantitatively and qualitatively altered in the infected plant, leading to the reactivation of at least one cotton transposable element. This is the first study to explore the global alterations of sRNAs in virus-infected cotton plants. Our results indicate that some CLRDV-induced symptoms may be correlated with the deregulation of miRNA and/or epigenetic networks.     

Comparative microarray analysis of microRNA expression profiles in primary cutaneous malignant melanoma, cutaneous malignant melanoma metastases, and benign melanocytic nevi

Perturbations in microRNA (miRNA) expression profiles have been reported for cutaneous malignant melanoma (CMM) predominantly when examined in cell lines. Despite the rapidly growing number of newly discovered human miRNA sequences, the availability of up-to-date miRNA expression profiles for clinical samples of primary cutaneous malignant melanoma (PCMM), cutaneous malignant melanoma metastases (CMMM), and benign melanocytic nevi (BMN) is limited. Specimens excised from the center of tumors (lesional) from patients with PCMM (n=9), CMMM (n=4), or BMN (n=8) were obtained during surgery. An exploratory microarray analysis was performed by miRNA expression profiling based on Agilent platform screening for 1205 human miRNAs. The results from the microarray analysis were validated by TaqMan quantitative real-time polymerase chain reaction. In addition to several miRNAs previously known to be associated with CMM, 19 unidentified miRNA candidates were found to be dysregulated in CMM patient samples. Among the 19 novel miRNA candidates, the genes hsa-miR-22, hsa-miR-130b, hsa-miR-146b-5p, hsa-miR-223, hsa-miR-301a, hsa-miR-484, hsa-miR-663, hsa-miR-720, hsa-miR-1260, hsa-miR-1274a, hsa-miR-1274b, hsa-miR-3663-3p, hsa-miR-4281, and hsa-miR-4286 were upregulated, and the genes hsa-miR-24-1*, hsa-miR-26a, hsa-miR-4291, hsa-miR-4317, and hsa-miR-4324 were downregulated. The results of this study partially confirm previous CMM miRNA profiling studies identifying miRNAs that are dysregulated in CMM. However, we report several novel miRNA candidates in CMM tumors; these miRNA sequences require further validation and functional analysis to evaluate whether they play a role in the pathogenesis of CMM.     

Strategies for profiling microRNA expression

MicroRNAs (miRNAs) are a class of small RNAs ( approximately 22-nt) that play an important role in the control of different cell processes by negative regulation of protein-coding genes. In the last several years, a number of miRNA profiling strategies have been used to document the miRNA expression changes during physiological and pathological processes. Aberrant expression of miRNAs has been linked to developmental defects, cancer, neurological disorders, and heart diseases. Over 540 human miRNAs have been validated to date; however, computer models suggest there may be thousands more. As bench work continue to verify in silico predictions, miRNA profiling will remain a prominent tool for identification of differential expression miRNAs in normal cellular courses and human disorders. This review focuses on current strategies for miRNA expression profiling and discusses their sensitivity and specificity, as well as advantage and disadvantage.     

Acting antisense: plasmid- and chromosome-encoded sRNAs from Gram-positive bacteria

sRNAs that act by base pairing were first discovered in plasmids, phages and transposons, where they control replication, maintenance and transposition. Since 2001, however, computational searches were applied that led to the discovery of a plethora of sRNAs in bacterial chromosomes. Whereas the majority of these chromsome-encoded sRNAs have been investigated in Escherichia coli, Salmonella and other Gram-negative bacteria, only a few well-studied examples are known from Gram-positive bacteria. Here, the author summarizes our current knowledge on plasmid- and chromosome-encoded sRNAs from Gram-positive species, thereby focusing on regulatory mechanisms used by these RNAs and their biological role in complex networks. Furthermore, regulatory factors that control the expression of these RNAs will be discussed and differences between sRNAs from Gram-positive and Gram-negative bacteria highlighted. The main emphasis of this review is on sRNAs that act by base pairing (i.e., by an antisense mechanism). Thereby, both plasmid-encoded and chromosome-encoded sRNAs will be considered.     

Non-coding RNAs in marine Synechococcus and their regulation under environmentally relevant stress conditions

Regulatory small RNAs (sRNAs) have crucial roles in the adaptive responses of bacteria to changes in the environment. Thus far, potential regulatory RNAs have been studied mainly in marine picocyanobacteria in genetically intractable Prochlorococcus, rendering their molecular analysis difficult. Synechococcus sp. WH7803 is a model cyanobacterium, representative of the picocyanobacteria from the mesotrophic areas of the ocean. Similar to the closely related Prochlorococcus it possesses a relatively streamlined genome and a small number of genes, but is genetically tractable. Here, a comparative genome analysis was performed for this and four additional marine Synechococcus to identify the suite of possible sRNAs and other RNA elements. Based on the prediction and on complementary microarray profiling, we have identified several known as well as 32 novel sRNAs. Some sRNAs overlap adjacent coding regions, for instance for the central photosynthetic gene psbA. Several of these novel sRNAs responded specifically to environmentally relevant stress conditions. Among them are six sRNAs changing their accumulation level under cold stress, six responding to high light and two to iron limitation. Target predictions suggested genes encoding components of the light-harvesting apparatus as targets of sRNAs originating from genomic islands and that one of the iron-regulated sRNAs might be a functional homolog of RyhB. These data suggest that marine Synechococcus mount adaptive responses to these different stresses involving regulatory sRNAs.     

Transposable element-associated microRNA hairpins produce 21-nt sRNAs integrated into typical microRNA pathways in rice

microRNAs (miRNAs) are a class of small RNAs (sRNAs) of ～21 nucleotides (nt) in length processed from foldback hairpins by DICER-LIKE1 (DCL1) or DCL4. They regulate the expression of target mRNAs by base pairing through RNA-induced silencing complex (RISC). In the RISC, ARGONAUTE1 (AGO1) is the key protein that cleaves miRNA targets at position ten of a miRNA:target duplex. The authenticity of many annotated rice miRNA hairpins is under debate because of their homology to repeat sequences. Some of them, like miR1884b, have been removed from the current release of miRBase based on incomplete information. In this study, we investigated the association of transposable element (TE)-derived miRNAs with typical miRNA pathways (DCL1/4- and AGO1-dependent) using publicly available deep sequencing datasets. Seven miRNA hairpins with 13 unique sRNAs were specifically enriched in AGO1 immunoprecipitation samples and relatively reduced in DCL1/4 knockdown genotypes. Interestingly, these species are ～21-nt long, instead of 24-nt as annotated in miRBase and the literature. Their expression profiles meet current criteria for functional annotation of miRNAs. In addition, diagnostic cleavage tags were found in degradome datasets for predicted target mRNAs. Most of these miRNA hairpins share significant homology with miniature inverted-repeat transposable elements, one type of abundant DNA transposons in rice. Finally, the root-specific production of a 24-nt miRNA-like sRNA was confirmed by RNA blot for a novel EST that maps to the 3′-UTR of a candidate pseudogene showing extensive sequence homology to miR1884b hairpin. Our data are consistent with the hypothesis that TEs can serve as a driving force for the evolution of some MIRNAs, where co-opting of DICER-LIKE1/4 processing and integration into AGO1 could exapt transcribed TE-associated hairpins into typical miRNA pathways.     

Genome-wide profiling of novel and conserved <Emphasis Type="Italic">Populus</Emphasis> microRNAs involved in pathogen stress response by deep sequencing

MicroRNAs (miRNAs) are small RNAs, generally of 20–23 nt, that down-regulate target gene expression during development, differentiation, growth, and metabolism. In Populus, extensive studies of miRNAs involved in cold, heat, dehydration, salinity, and mechanical stresses have been performed; however, there are few reports profiling the miRNA expression patterns during pathogen stress. We obtained almost 38 million raw reads through Solexa sequencing of two libraries from Populus inoculated and uninoculated with canker disease pathogen. Sequence analyses identified 74 conserved miRNA sequences belonging to 37 miRNA families from 154 loci in the Populus genome and 27 novel miRNA sequences from 35 loci, including their complementary miRNA* strands. Intriguingly, the miRNA* of three conserved miRNAs were more abundant than their corresponding miRNAs. The overall expression levels of conserved miRNAs increased when subjected to pathogen stress, and expression levels of 33 miRNA sequences markedly changed. The expression trends determined by sequencing and by qRT-PCR were similar. Finally, nine target genes for three conserved miRNAs and 63 target genes for novel miRNAs were predicted using computational analysis, and their functions were annotated. Deep sequencing provides an opportunity to identify pathogen-regulated miRNAs in trees, which will help in understanding the regulatory mechanisms of plant defense responses during pathogen infection.     

Normalization strategies for microRNA profiling experiments: a ‘normal’ way to a hidden layer of complexity?

MicroRNA (miRNA) profiling is a first important step in elucidating miRNA functions. Real time quantitative PCR (RT-qPCR) and microarray hybridization approaches as well as ultra high throughput sequencing of miRNAs (small RNA-seq) are popular and widely used profiling methods. All of these profiling approaches face significant introduction of bias. Normalization, often an underestimated aspect of data processing, can minimize systematic technical or experimental variation and thus has significant impact on the detection of differentially expressed miRNAs. At present, there is no consensus normalization method for any of the three miRNA profiling approach. Several normalization techniques are currently in use, of which some are similar to mRNA profiling normalization methods, while others are specifically modified or developed for miRNA data. The characteristic nature of miRNA molecules, their composition and the resulting data distribution of profiling experiments challenges the selection of adequate normalization techniques. Based on miRNA profiling studies and comparative studies on normalization methods and their performances, this review provides a critical overview of commonly used and newly developed normalization methods for miRNA RT-qPCR, miRNA hybridization microarray, and small RNA-seq datasets. Emphasis is laid on the complexity, the importance and the potential for further optimization of normalization techniques for miRNA profiling datasets.