Computational and analytical framework for small RNA profiling by high-throughput sequencing

Bioinformatic analysis of barcoded cDNA libraries for small RNA profiling by next-generation sequencing

The characterization of post-transcriptional gene regulation by small regulatory RNAs of 20-30 nt length, particularly miRNAs and piRNAs, has become a major focus of research in recent years. A prerequisite for the characterization of small RNAs is their identification and quantification across different developmental stages, normal and diseased tissues, as well as model cell lines. Here we present a step-by-step protocol for the bioinformatic analysis of barcoded cDNA libraries for small RNA profiling generated by Illumina sequencing, thereby facilitating miRNA and other small RNA profiling of large sample collections.

Functional analysis of three Arabidopsis ARGONAUTES using slicer-defective mutants

In RNA-directed silencing pathways, ternary complexes result from small RNA-guided ARGONAUTE (AGO) associating with target transcripts. Target transcripts are often silenced through direct cleavage (slicing), destabilization through slicer-independent turnover mechanisms, and translational repression. Here, wild-type and active-site defective forms of several Arabidopsis thaliana AGO proteins involved in posttranscriptional silencing were used to examine several AGO functions, including small RNA binding, interaction with target RNA, slicing or destabilization of target RNA, secondary small interfering RNA formation, and antiviral activity. Complementation analyses in ago mutant plants revealed that the catalytic residues of AGO1, AGO2, and AGO7 are required to restore the defects of Arabidopsis ago1-25, ago2-1, and zip-1 (AGO7-defective) mutants, respectively. AGO2 had slicer activity in transient assays but could not trigger secondary small interfering RNA biogenesis, and catalytically active AGO2 was necessary for local and systemic antiviral activity against Turnip mosaic virus. Slicer-defective AGOs associated with miRNAs and stabilized AGO-miRNA-target RNA ternary complexes in individual target coimmunoprecipitation assays. In genome-wide AGO-miRNA-target RNA coimmunoprecipitation experiments, slicer-defective AGO1-miRNA associated with target RNA more effectively than did wild-type AGO1-miRNA. These data not only reveal functional roles for AGO1, AGO2, and AGO7 slicer activity, but also indicate an approach to capture ternary complexes more efficiently for genome-wide analyses.

MicroRNAs regulate tumor angiogenesis modulated by endothelial progenitor cells

Bone marrow-derived endothelial progenitor cells (EPC) contribute to the angiogenesis-dependent growth of tumors in mice and humans. EPCs regulate the angiogenic switch via paracrine secretion of proangiogenic growth factors and by direct luminal incorporation into sprouting nascent vessels. miRNAs have emerged as key regulators of several cellular processes including angiogenesis; however, whether miRNAs contribute to bone marrow-mediated angiogenesis has remained unknown. Here, we show that genetic ablation of miRNA-processing enzyme Dicer, specifically in the bone marrow, decreased the number of circulating EPCs, resulting in angiogenesis suppression and impaired tumor growth. Furthermore, genome-wide deep sequencing of small RNAs revealed tumor EPC-intrinsic miRNAs including miR-10b and miR-196b, which have been previously identified as key regulators of HOX signaling and adult stem cell differentiation. Notably, we found that both miR-10b and miR-196b are responsive to vascular endothelial growth factor stimulation and show elevated expression in human high-grade breast tumor vasculature. Strikingly, targeting miR-10b and miR-196b led to significant defects in angiogenesis-mediated tumor growth in mice. Targeting these miRNAs may constitute a novel strategy for inhibiting tumor angiogenesis.

Towards an understanding of the epigenetics of schistosomes: a comparative epigenomic study

As in perhaps all eukaryotes, schistosomes use a supplementary information transmitting system, the epigenetic inheritance system, to shape genetic information and to produce different phenotypes. In contrast to other important parasites, the study of epigenetic phenomena in schistosomes is still in its infancy. Nevertheless, we are beginning to grasp what goes on behind the epigenetic scene in this parasite. We have developed techniques of native chromatin immunoprecipitation (N-ChIP) and associated the necessary bioinformatics tools that allow us to run genome-wide comparative chromatin studies on Schistosoma mansoni at different stages of its life cycle, on different strains and on different sexes. We present here an application of such an approach to study the genetic and epigenetic basis for a phenotypic trait, the compatibility of S. mansoni with its invertebrate host Biomphalaria glabrata. We have applied the ChIP procedure to two strains that are either compatible or incompatible with their intermediate host. The precipitated DNA was sequenced and aligned to a reference genome and this information was used to determine regions in which both strands differ in their genomic sequence and/or chromatin structure. This procedure allowed us to identify candidate genes that display either genetic or epigenetic difference between the two strains.

PsRobot: a web-based plant small RNA meta-analysis toolbox

Small RNAs (smRNAs) in plants, mainly microRNAs and small interfering RNAs, play important roles in both transcriptional and post-transcriptional gene regulation. The broad application of high-throughput sequencing technology has made routinely generation of bulk smRNA sequences in laboratories possible, thus has significantly increased the need for batch analysis tools. PsRobot is a web-based easy-to-use tool dedicated to the identification of smRNAs with stem-loop shaped precursors (such as microRNAs and short hairpin RNAs) and their target genes/transcripts. It performs fast analysis to identify smRNAs with stem-loop shaped precursors among batch input data and predicts their targets using a modified Smith–Waterman algorithm. PsRobot integrates the expression data of smRNAs in major plant smRNA biogenesis gene mutants and smRNA-associated protein complexes to give clues to the smRNA generation and functional processes. Besides improved specificity, the reliability of smRNA target prediction results can also be evaluated by mRNA cleavage (degradome) data. The cross species conservation statuses and the multiplicity of smRNA target sites are also provided. PsRobot is freely accessible at http://omicslab.genetics.ac.cn/psRobot/.

PIWI associated siRNAs and piRNAs specifically require the Caenorhabditis elegans HEN1 ortholog henn-1

Small RNAs--including piRNAs, miRNAs, and endogenous siRNAs--bind Argonaute proteins to form RNA silencing complexes that target coding genes, transposons, and aberrant RNAs. To assess the requirements for endogenous siRNA formation and activity in Caenorhabditis elegans, we developed a GFP-based sensor for the endogenous siRNA 22G siR-1, one of a set of abundant siRNAs processed from a precursor RNA mapping to the X chromosome, the X-cluster. Silencing of the sensor is also dependent on the partially complementary, unlinked 26G siR-O7 siRNA. We show that 26G siR-O7 acts in trans to initiate 22G siRNA formation from the X-cluster. The presence of several mispairs between 26G siR-O7 and the X-cluster mRNA, as well as mutagenesis of the siRNA sensor, indicates that siRNA target recognition is permissive to a degree of mispairing. From a candidate reverse genetic screen, we identified several factors required for 22G siR-1 activity, including the chromatin factors mes-4 and gfl-1, the Argonaute ergo-1, and the 3' methyltransferase henn-1. Quantitative RT-PCR of small RNAs in a henn-1 mutant and deep sequencing of methylated small RNAs indicate that siRNAs and piRNAs that associate with PIWI clade Argonautes are methylated by HENN-1, while siRNAs and miRNAs that associate with non-PIWI clade Argonautes are not. Thus, PIWI-class Argonaute proteins are specifically adapted to associate with methylated small RNAs in C. elegans.

PAREsnip: a tool for rapid genome-wide discovery of small RNA/target interactions evidenced through degradome sequencing

Small RNAs (sRNAs) are a class of short (20–25 nt) non-coding RNAs that play important regulatory roles in gene expression. An essential first step in understanding their function is to confidently identify sRNA targets. In plants, several classes of sRNAs such as microRNAs (miRNAs) and trans-acting small interfering RNAs have been shown to bind with near-perfect complementarity to their messenger RNA (mRNA) targets, generally leading to cleavage of the mRNA. Recently, a high-throughput technique known as Parallel Analysis of RNA Ends (PARE) has made it possible to sequence mRNA cleavage products on a large-scale. Computational methods now exist to use these data to find targets of conserved and newly identified miRNAs. Due to speed limitations such methods rely on the user knowing which sRNA sequences are likely to target a transcript. By limiting the search to a tiny subset of sRNAs it is likely that many other sRNA/mRNA interactions will be missed. Here, we describe a new software tool called PAREsnip that allows users to search for potential targets of all sRNAs obtained from high-throughput sequencing experiments. By searching for targets of a complete ‘sRNAome’ we can facilitate large-scale identification of sRNA targets, allowing us to discover regulatory interaction networks.

Maternal siRNAs as regulators of parental genome imbalance and gene expression in endosperm of Arabidopsis seeds

Seed size is important to crop domestication and natural selection and is affected by the balance of maternal and paternal genomes in endosperm. Endosperm, like placenta in mammals, provides reserves to the developing embryo. Interploidy crosses disrupt the genome balance in endosperm and alter seed size. Specifically, paternal-excess crosses (2 × 4) delay endosperm cellularization (EC) and produce larger seeds, whereas maternal-excess crosses (4 × 2) promote precocious EC and produce smaller seeds. The mechanisms for responding to the parental genome dosage imbalance and for gene expression changes in endosperm are unknown. In plants, RNA polymerase IV (PolIV or p4) encoded by NRPD1a is required for biogenesis of a major class of 24-nt small interfering RNAs (also known as p4-siRNAs), which are predominately expressed in developing endosperm. Here we show that p4-siRNA accumulation depends on the maternal genome dosage, and maternal p4-siRNAs target transposable elements (TEs) and TE-associated genes (TAGs) in seeds. The p4-siRNAs correlate negatively with expression levels of AGAMOUS-LIKE (AGL) genes in endosperm of interploidy crosses. Moreover, disruption of maternal NRPD1a expression is associated with p4-siRNA reduction and AGL up-regulation in endosperm of reciprocal crosses. This is unique genetic evidence for maternal siRNAs in response to parental genome imbalance and in control of transposons and gene expression during endosperm development.

Library preparation and data analysis packages for rapid genome sequencing

High-throughput sequencing (HTS) has quickly become a valuable tool for comparative genetics and genomics and is now regularly carried out in laboratories that are not connected to large sequencing centers. Here we describe an updated version of our protocol for constructing single- and paired-end Illumina sequencing libraries, beginning with purified genomic DNA. The present protocol can also be used for "multiplexing," i.e. the analysis of several samples in a single flowcell lane by generating "barcoded" or "indexed" Illumina sequencing libraries in a way that is independent from Illumina-supported methods. To analyze sequencing results, we suggest several independent approaches but end users should be aware that this is a quickly evolving field and that currently many alignment (or "mapping") and counting algorithms are being developed and tested.

Deep sequencing of small RNAs in plants: applied bioinformatics

Small RNAs, including microRNA and short-interfering RNAs, play important roles in plants. In recent years, developments in sequencing technology have enabled the large-scale discovery of sRNAs in various cells, tissues and developmental stages and in response to various stresses. This review describes the bioinformatics challenges to analysing these large datasets of short-RNA sequences and some of the solutions to those challenges.

High-throughput sequencing of RNA silencing-associated small RNAs in olive (Olea europaea L.)

Small RNAs (sRNAs) of 20 to 25 nucleotides (nt) in length maintain genome integrity and control gene expression in a multitude of developmental and physiological processes. Despite RNA silencing has been primarily studied in model plants, the advent of high-throughput sequencing technologies has enabled profiling of the sRNA component of more than 40 plant species. Here, we used deep sequencing and molecular methods to report the first inventory of sRNAs in olive (Olea europaea L.). sRNA libraries prepared from juvenile and adult shoots revealed that the 24-nt class dominates the sRNA transcriptome and atypically accumulates to levels never seen in other plant species, suggesting an active role of heterochromatin silencing in the maintenance and integrity of its large genome. A total of 18 known miRNA families were identified in the libraries. Also, 5 other sRNAs derived from potential hairpin-like precursors remain as plausible miRNA candidates. RNA blots confirmed miRNA expression and suggested tissue- and/or developmental-specific expression patterns. Target mRNAs of conserved miRNAs were computationally predicted among the olive cDNA collection and experimentally validated through endonucleolytic cleavage assays. Finally, we use expression data to uncover genetic components of the miR156, miR172 and miR390/TAS3-derived trans-acting small interfering RNA (tasiRNA) regulatory nodes, suggesting that these interactive networks controlling developmental transitions are fully operational in olive.

The ERI-6/7 helicase acts at the first stage of an siRNA amplification pathway that targets recent gene duplications

Endogenous small interfering RNAs (siRNAs) are a class of naturally occuring regulatory RNAs found in fungi, plants, and animals. Some endogenous siRNAs are required to silence transposons or function in chromosome segregation; however, the specific roles of most endogenous siRNAs are unclear. The helicase gene eri-6/7 was identified in the nematode Caenorhabditis elegans by the enhanced response to exogenous double-stranded RNAs (dsRNAs) of the null mutant. eri-6/7 encodes a helicase homologous to small RNA factors Armitage in Drosophila, SDE3 in Arabidopsis, and Mov10 in humans. Here we show that eri-6/7 mutations cause the loss of 26-nucleotide (nt) endogenous siRNAs derived from genes and pseudogenes in oocytes and embryos, as well as deficiencies in somatic 22-nucleotide secondary siRNAs corresponding to the same loci. About 80 genes are eri-6/7 targets that generate the embryonic endogenous siRNAs that silence the corresponding mRNAs. These 80 genes share extensive nucleotide sequence homology and are poorly conserved, suggesting a role for these endogenous siRNAs in silencing of and thereby directing the fate of recently acquired, duplicated genes. Unlike most endogenous siRNAs in C. elegans, eri-6/7–dependent siRNAs require Dicer. We identify that the eri-6/7–dependent siRNAs have a passenger strand that is ∼19 nt and is inset by ∼3–4 nts from both ends of the 26 nt guide siRNA, suggesting non-canonical Dicer processing. Mutations in the Argonaute ERGO-1, which associates with eri-6/7–dependent 26 nt siRNAs, cause passenger strand stabilization, indicating that ERGO-1 is required to separate the siRNA duplex, presumably through endonucleolytic cleavage of the passenger strand. Thus, like several other siRNA–associated Argonautes with a conserved RNaseH motif, ERGO-1 appears to be required for siRNA maturation.

Endogenous small interfering RNAs (siRNAs) are a class of small RNAs present in fungi, plants, and animals. Small RNAs, including microRNAs, are known to regulate the expression levels of genes, silence invading elements such as transposons, and act in cell division. However, the function of many endogenous siRNAs is unknown. We have found that the ERI-6/7 helicase is required for a subset of endogenous siRNAs present in the nematode Caenorhabditis elegans. The ERI-6/7 helicase acts in a pathway together with the Argonaute protein ERGO-1 to produce two types of siRNAs: a primary class of 26 nucleotides in length present in oocytes and embryos, and a class of 22 nucleotide siRNAs present in later stages of development. These siRNAs correspond to only about one hundred genes. Interestingly, we found that these genes fall into groups of genes that contain nearly identical DNA sequences. The sequences of these genes are not conserved in other organisms, not even in related nematodes. These results point to a potential function of these endogenous siRNAs: silencing of recently acquired, duplicated genes. Our work demonstrates a new role of small RNAs, different from known functions in transposon silencing and regulation of gene expression.

GENE-counter: a computational pipeline for the analysis of RNA-Seq data for gene expression differences

GENE-counter is a complete Perl-based computational pipeline for analyzing RNA-Sequencing (RNA-Seq) data for differential gene expression. In addition to its use in studying transcriptomes of eukaryotic model organisms, GENE-counter is applicable for prokaryotes and non-model organisms without an available genome reference sequence. For alignments, GENE-counter is configured for CASHX, Bowtie, and BWA, but an end user can use any Sequence Alignment/Map (SAM)-compliant program of preference. To analyze data for differential gene expression, GENE-counter can be run with any one of three statistics packages that are based on variations of the negative binomial distribution. The default method is a new and simple statistical test we developed based on an over-parameterized version of the negative binomial distribution. GENE-counter also includes three different methods for assessing differentially expressed features for enriched gene ontology (GO) terms. Results are transparent and data are systematically stored in a MySQL relational database to facilitate additional analyses as well as quality assessment. We used next generation sequencing to generate a small-scale RNA-Seq dataset derived from the heavily studied defense response of Arabidopsis thaliana and used GENE-counter to process the data. Collectively, the support from analysis of microarrays as well as the observed and substantial overlap in results from each of the three statistics packages demonstrates that GENE-counter is well suited for handling the unique characteristics of small sample sizes and high variability in gene counts.

Genome-wide evidence for local DNA methylation spreading from small RNA-targeted sequences in Arabidopsis

Transposable elements (TEs) and their relics play major roles in genome evolution. However, mobilization of TEs is usually deleterious and strongly repressed. In plants and mammals, this repression is typically associated with DNA methylation, but the relationship between this epigenetic mark and TE sequences has not been investigated systematically. Here, we present an improved annotation of TE sequences and use it to analyze genome-wide DNA methylation maps obtained at single-nucleotide resolution in Arabidopsis. We show that although the majority of TE sequences are methylated, ∼26% are not. Moreover, a significant fraction of TE sequences densely methylated at CG, CHG and CHH sites (where H = A, T or C) have no or few matching small interfering RNA (siRNAs) and are therefore unlikely to be targeted by the RNA-directed DNA methylation (RdDM) machinery. We provide evidence that these TE sequences acquire DNA methylation through spreading from adjacent siRNA-targeted regions. Further, we show that although both methylated and unmethylated TE sequences located in euchromatin tend to be more abundant closer to genes, this trend is least pronounced for methylated, siRNA-targeted TE sequences located 5' to genes. Based on these and other findings, we propose that spreading of DNA methylation through promoter regions explains at least in part the negative impact of siRNA-targeted TE sequences on neighboring gene expression.

Genome sequencing and comparative analysis of the carrot bacterial blight pathogen, Xanthomonas hortorum pv. carotae M081, for insights into pathogenicity and applications in molecular diagnostics

Xanthomonas hortorum pv. carotae (Xhc) is an economically important pathogen of carrots. Its ability to epiphytically colonize foliar surfaces and infect seeds can result in bacterial blight of carrots when grown in warm and humid regions. We used high-throughput sequencing to determine the genome sequence of isolate M081 of Xhc. The short reads were de novo assembled and the resulting contigs were ordered using a syntenic reference genome sequence from X. campestris pv. campestris ATCC 33913. The improved, high-quality draft genome sequence of Xhc M081 is the first for its species. Despite its distance from other sequenced xanthomonads, Xhc M081 still shared a large inventory of orthologous genes, including many clusters of virulence genes common to other foliar pathogenic species of Xanthomonas. We also mined the genome sequence and identified at least 21 candidate type III effector genes. Two were members of the avrBs2 and xopQ families that demonstrably elicit effector-triggered immunity. We showed that expression in planta of these two type III effectors from Xhc M081 was sufficient to elicit resistance gene-mediated hypersensitive responses in heterologous plants, indicating a possibility for resistance gene-mediated control of Xhc. Finally, we identified regions unique to the Xhc M081 genome sequence, and demonstrated their potential in the design of molecular diagnostics for this pathogen.

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.

MicroRNA expression profiling reveals miRNA families regulating specific biological pathways in mouse frontal cortex and hippocampus

MicroRNAs (miRNAs) are small regulatory molecules that cause post-transcriptional gene silencing. Although some miRNAs are known to have region-specific expression patterns in the adult brain, the functional consequences of the region-specificity to the gene regulatory networks of the brain nuclei are not clear. Therefore, we studied miRNA expression patterns by miRNA-Seq and microarrays in two brain regions, frontal cortex (FCx) and hippocampus (HP), which have separate biological functions. We identified 354 miRNAs from FCx and 408 from HP using miRNA-Seq, and 245 from FCx and 238 from HP with microarrays. Several miRNA families and clusters were differentially expressed between FCx and HP, including the miR-8 family, miR-182|miR-96|miR-183 cluster, and miR-212|miR-312 cluster overexpressed in FCx and miR-34 family overexpressed in HP. To visualize the clusters, we developed support for viewing genomic alignments of miRNA-Seq reads in the Chipster genome browser. We carried out pathway analysis of the predicted target genes of differentially expressed miRNA families and clusters to assess their putative biological functions. Interestingly, several miRNAs from the same family/cluster were predicted to regulate specific biological pathways. We have developed a miRNA-Seq approach with a bioinformatic analysis workflow that is suitable for studying miRNA expression patterns from specific brain nuclei. FCx and HP were shown to have distinct miRNA expression patterns which were reflected in the predicted gene regulatory pathways. This methodology can be applied for the identification of brain region-specific and phenotype-specific miRNA-mRNA-regulatory networks from the adult and developing rodent brain.

Wheat hybridization and polyploidization results in deregulation of small RNAs

Speciation via interspecific or intergeneric hybridization and polyploidization triggers genomic responses involving genetic and epigenetic alterations. Such modifications may be induced by small RNAs, which affect key cellular processes, including gene expression, chromatin structure, cytosine methylation and transposable element (TE) activity. To date, the role of small RNAs in the context of wide hybridization and polyploidization has received little attention. In this work, we performed high-throughput sequencing of small RNAs of parental, intergeneric hybrid, and allopolyploid plants that mimic the genomic changes occurring during bread wheat speciation. We found that the percentage of small RNAs corresponding to miRNAs increased with ploidy level, while the percentage of siRNAs corresponding to TEs decreased. The abundance of most miRNA species was similar to midparent values in the hybrid, with some deviations, as seen in overrepresentation of miR168, in the allopolyploid. In contrast, the number of siRNAs corresponding to TEs strongly decreased upon allopolyploidization, but not upon hybridization. The reduction in corresponding siRNAs, together with decreased CpG methylation, as shown here for the Veju element, represent hallmarks of TE activation. TE-siRNA downregulation in the allopolyploid may contribute to genome destabilization at the initial stages of speciation. This phenomenon is reminiscent of hybrid dysgenesis in Drosophila.

HC-Pro silencing suppressor significantly alters the gene expression profile in tobacco leaves and flowers

BACKGROUND

RNA silencing is used in plants as a major defence mechanism against invasive nucleic acids, such as viruses. Accordingly, plant viruses have evolved to produce counter defensive RNA-silencing suppressors (RSSs). These factors interfere in various ways with the RNA silencing machinery in cells, and thereby disturb the microRNA (miRNA) mediated endogene regulation and induce developmental and morphological changes in plants. In this study we have explored these effects using previously characterized transgenic tobacco plants which constitutively express (under CaMV 35S promoter) the helper component-proteinase (HC-Pro) derived from a potyviral genome. The transcript levels of leaves and flowers of these plants were analysed using microarray techniques (Tobacco 4 × 44 k, Agilent).

RESULTS

Over expression of HC-Pro RSS induced clear phenotypic changes both in growth rate and in leaf and flower morphology of the tobacco plants. The expression of 748 and 332 genes was significantly changed in the leaves and flowers, respectively, in the HC-Pro expressing transgenic plants. Interestingly, these transcriptome alterations in the HC-Pro expressing tobacco plants were similar as those previously detected in plants infected with ssRNA-viruses. Particularly, many defense-related and hormone-responsive genes (e.g. ethylene responsive transcription factor 1, ERF1) were differentially regulated in these plants. Also the expression of several stress-related genes, and genes related to cell wall modifications, protein processing, transcriptional regulation and photosynthesis were strongly altered. Moreover, genes regulating circadian cycle and flowering time were significantly altered, which may have induced a late flowering phenotype in HC-Pro expressing plants. The results also suggest that photosynthetic oxygen evolution, sugar metabolism and energy levels were significantly changed in these transgenic plants. Transcript levels of S-adenosyl-L-methionine (SAM) were also decreased in these plants, apparently leading to decreased transmethylation capacity. The proteome analysis using 2D-PAGE indicated significantly altered proteome profile, which may have been both due to altered transcript levels, decreased translation, and increased proteosomal/protease activity.

CONCLUSION

Expression of the HC-Pro RSS mimics transcriptional changes previously shown to occur in plants infected with intact viruses (e.g. Tobacco etch virus, TEV). The results indicate that the HC-Pro RSS contributes a significant part of virus-plant interactions by changing the levels of multiple cellular RNAs and proteins.

Heterochromatin is required for normal distribution of Neurospora crassa CenH3

Centromeres serve as platforms for the assembly of kinetochores and are essential for nuclear division. Here we identified Neurospora crassa centromeric DNA by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) of DNA associated with tagged versions of the centromere foundation proteins CenH3 (CENP-A) and CEN-C (CENP-C) and the kinetochore protein CEN-T (CENP-T). On each chromosome we found an ∼150- to 300-kbp region of enrichment for all three proteins. These regions correspond to intervals predicted to be centromeric DNA by genetic mapping and DNA sequence analyses. By ChIP-seq we found extensive colocalization of CenH3, CEN-C, CEN-T, and histone H3K9 trimethylation (H3K9me3). In contrast, H3K4me2, which has been found at the cores of plant, fission yeast, Drosophila, and mammalian centromeres, was not enriched in Neurospora centromeric DNA. DNA methylation was most pronounced at the periphery of centromeric DNA. Mutation of dim-5, which encodes an H3K9 methyltransferase responsible for nearly all H3K9me3, resulted in altered distribution of CenH3-green fluorescent protein (GFP). Similarly, CenH3-GFP distribution was altered in the absence of HP1, the chromodomain protein that binds to H3K9me3. We conclude that eukaryotes with regional centromeres make use of different strategies for maintenance of CenH3 at centromeres, and we suggest a model in which centromere proteins nucleate at the core but require DIM-5 and HP1 for spreading.

Investigation of the microRNAs in safflower seed, leaf, and petal by high-throughput sequencing

Safflower (Carthamus tinctorius L.) is one of the most important crop plants that has been utilized for production of oleosins. miRNAs (microRNAs) are a class of small and non-coding RNAs that negatively regulate gene expression at post-transcriptional level thus playing a role in plant growth, development, and stress response. In this study, high-throughput Illumina sequencing technology has been used to comprehensively investigate the small RNA transcriptomes of safflower seed, flower, and leaf. It is found that there are at least 236 known miRNAs expressed in safflower, of which 100 miRNAs with relatively high expression abundance exhibited evolutionary conservation across multiple plants. Comparison of their expression abundance among different tissues shows that a total of 116, 133, and 128 miRNAs are significantly differentially expressed with higher abundance or lower abundance between safflower seed/leaf, seed/petal, and leaf/petal. The majority of the most significant differences in miRNA abundance between tissues are tissue-specific miRNAs. In addition, 13 putative novel miRNAs have been identified in safflower. The small RNA transcriptomes obtained in this study provide a basis for further investigation of the physiological roles of identified miRNAs in safflower.

Experimental design, preprocessing, normalization and differential expression analysis of small RNA sequencing experiments

Prior to the advent of new, deep sequencing methods, small RNA (sRNA) discovery was dependent on Sanger sequencing, which was time-consuming and limited knowledge to only the most abundant sRNA. The innovation of large-scale, next-generation sequencing has exponentially increased knowledge of the biology, diversity and abundance of sRNA populations. In this review, we discuss issues involved in the design of sRNA sequencing experiments, including choosing a sequencing platform, inherent biases that affect sRNA measurements and replication. We outline the steps involved in preprocessing sRNA sequencing data and review both the principles behind and the current options for normalization. Finally, we discuss differential expression analysis in the absence and presence of biological replicates. While our focus is on sRNA sequencing experiments, many of the principles discussed are applicable to the sequencing of other RNA populations.

Computational dissection of Arabidopsis smRNAome leads to discovery of novel microRNAs and short interfering RNAs associated with transcription start sites

The profiling of small RNAs by high-throughput sequencing (smRNA-Seq) has revealed the complexity of the RNA world. Here, we describe a computational scheme for dissecting the plant smRNAome by integrating smRNA-Seq datasets in Arabidopsis thaliana. Our analytical approach first defines ab initio the genomic loci that produce smRNAs as basic units, then utilizes principal component analysis (PCA) to predict novel miRNAs. Secondary structure prediction of candidates' putative precursors discovered a group of long hairpin double-stranded RNAs (lh-dsRNAs) formed by inverted duplications of decayed coding genes. These gene remnants produce miRNA-like small RNAs which are predominantly 21- and 22-nt long, dependent of DCL1 but independent of RDR2 and DCL2/3/4, and associated with AGO1. Additionally, we found two classes of transcription start site associated (TSSa) RNAs located at sense (+) and antisense (-) approximately 100-200 bp downstream of TSSs, but are differentially incorporated into AGO1 and AGO4, respectively.

mut-16 and other mutator class genes modulate 22G and 26G siRNA pathways in Caenorhabditis elegans

Argonaute-associated siRNAs and Piwi-associated piRNAs have overlapping roles in silencing mobile genetic elements in animals. In Caenorhabditis elegans, mutator (mut) class genes mediate siRNA-guided repression of transposons as well as exogenous RNAi, but their roles in endogenous RNA silencing pathways are not well-understood. To characterize the endogenous small RNAs dependent on mut class genes, small RNA populations from a null allele of mut-16 as well as a regulatory mut-16(mg461) allele that disables only somatic RNAi were subjected to deep sequencing. Additionally, each of the mut class genes was tested for a requirement in 26G siRNA pathways. The results indicate that mut-16 is an essential factor in multiple endogenous germline and somatic siRNA pathways involving several distinct Argonautes and RNA-dependent RNA polymerases. The results also reveal essential roles for mut-2 and mut-7 in the ERGO-1 class 26G siRNA pathway and less critical roles for mut-8, mut-14, and mut-15. We show that transposons are hypersusceptible to mut-16-dependent silencing and identify a requirement for the siRNA machinery in piRNA biogenesis from Tc1 transposons. We also show that the soma-specific mut-16(mg461) mutant allele is present in multiple C. elegans laboratory strains.

Large-scale sequencing of plant small RNAs

Deep sequencing technologies have become very powerful tools in the identification and quantification of small RNAs involved in gene regulation. Small interfering RNA (siRNA) and miRNA are two classes of DCL-dependent small RNAs known to affect phenotype, developmental regulation, and various traits in plants. These small RNAs function by selectively repressing gene expression mainly by guiding cleavage, resulting in degradation of target transcripts. In this chapter, we describe a method for preparation of 5(')-phosphate-dependent small RNA libraries, a hallmark of RNase III-like DCL products, for high-throughput sequencing, and recommendations for small RNA analysis. This method is useful for determining small RNA involvement in critical pathways in plants, identifying and quantifying novel small RNAs, and examining small RNA global expression patterns.

Construction of small RNA cDNA libraries for high-throughput sequencing

Small RNAs (smRNAs) play an essential role in virtually every aspect of growth and development, by regulating gene expression at the post-transcriptional and/or transcriptional level. New high-throughput sequencing technology allows for a comprehensive coverage of smRNAs in any given biological sample, and has been widely used for profiling smRNA populations in various developmental stages, tissue and cell types, or normal and disease states. In this article, we describe the method used in our laboratory to construct smRNA cDNA libraries for high-throughput sequencing.

Digital gene expression signatures for maize development

Genome-wide expression signatures detect specific perturbations in developmental programs and contribute to functional resolution of key regulatory networks. In maize (Zea mays) inflorescences, mutations in the RAMOSA (RA) genes affect the determinacy of axillary meristems and thus alter branching patterns, an important agronomic trait. In this work, we developed and tested a framework for analysis of tag-based, digital gene expression profiles using Illumina's high-throughput sequencing technology and the newly assembled B73 maize reference genome. We also used a mutation in the RA3 gene to identify putative expression signatures specific to stem cell fate in axillary meristem determinacy. The RA3 gene encodes a trehalose-6-phosphate phosphatase and may act at the interface between developmental and metabolic processes. Deep sequencing of digital gene expression libraries, representing three biological replicate ear samples from wild-type and ra3 plants, generated 27 million 20- to 21-nucleotide reads with frequencies spanning 4 orders of magnitude. Unique sequence tags were anchored to 3'-ends of individual transcripts by DpnII and NlaIII digests, which were multiplexed during sequencing. We mapped 86% of nonredundant signature tags to the maize genome, which associated with 37,117 gene models and unannotated regions of expression. In total, 66% of genes were detected by at least nine reads in immature maize ears. We used comparative genomics to leverage existing information from Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) in functional analyses of differentially expressed maize genes. Results from this study provide a basis for the analysis of short-read expression data in maize and resolved specific expression signatures that will help define mechanisms of action for the RA3 gene.

An improved, high-quality draft genome sequence of the Germination-Arrest Factor-producing Pseudomonas fluorescens WH6

Background

Pseudomonas fluorescens is a genetically and physiologically diverse species of bacteria present in many habitats and in association with plants. This species of bacteria produces a large array of secondary metabolites with potential as natural products. P. fluorescens isolate WH6 produces Germination-Arrest Factor (GAF), a predicted small peptide or amino acid analog with herbicidal activity that specifically inhibits germination of seeds of graminaceous species.

Results

We used a hybrid next-generation sequencing approach to develop a high-quality draft genome sequence for P. fluorescens WH6. We employed automated, manual, and experimental methods to further improve the draft genome sequence. From this assembly of 6.27 megabases, we predicted 5876 genes, of which 3115 were core to P. fluorescens and 1567 were unique to WH6. Comparative genomic studies of WH6 revealed high similarity in synteny and orthology of genes with P. fluorescens SBW25. A phylogenomic study also placed WH6 in the same lineage as SBW25. In a previous non-saturating mutagenesis screen we identified two genes necessary for GAF activity in WH6. Mapping of their flanking sequences revealed genes that encode a candidate anti-sigma factor and an aminotransferase. Finally, we discovered several candidate virulence and host-association mechanisms, one of which appears to be a complete type III secretion system.

Conclusions

The improved high-quality draft genome sequence of WH6 contributes towards resolving the P. fluorescens species, providing additional impetus for establishing two separate lineages in P. fluorescens. Despite the high levels of orthology and synteny to SBW25, WH6 still had a substantial number of unique genes and represents another source for the discovery of genes with implications in affecting plant growth and health. Two genes are demonstrably necessary for GAF and further characterization of their proteins is important for developing natural products as control measure against grassy weeds. Finally, WH6 is the first isolate of P. fluorescens reported to encode a complete T3SS. This gives us the opportunity to explore the role of what has traditionally been thought of as a virulence mechanism for non-pathogenic interactions with plants.

Global effects of the small RNA biogenesis machinery on the Arabidopsis thaliana transcriptome

In Arabidopsis thaliana, four different dicer-like (DCL) proteins have distinct but partially overlapping functions in the biogenesis of microRNAs (miRNAs) and siRNAs from longer, noncoding precursor RNAs. To analyze the impact of different components of the small RNA biogenesis machinery on the transcriptome, we subjected dcl and other mutants impaired in small RNA biogenesis to whole-genome tiling array analysis. We compared both protein-coding genes and noncoding transcripts, including most pri-miRNAs, in two tissues and several stress conditions. Our analysis revealed a surprising number of common targets in dcl1 and dcl2 dcl3 dcl4 triple mutants. Furthermore, our results suggest that the DCL1 is not only involved in miRNA action but also contributes to silencing of a subset of transposons, apparently through an effect on DNA methylation.

NOVOMIR: De Novo Prediction of MicroRNA-Coding Regions in a Single Plant-Genome

MicroRNAs (miRNA) are small regulatory, noncoding RNA molecules that are transcribed as primary miRNAs (pri-miRNA) from eukaryotic genomes. At least in plants, their regulatory activity is mediated through base-pairing with protein-coding messenger RNAs (mRNA) followed by mRNA degradation or translation repression. We describe NOVOMIR, a program for the identification of miRNA genes in plant genomes. It uses a series of filter steps and a statistical model to discriminate a pre-miRNA from other RNAs and does rely neither on prior knowledge of a miRNA target nor on comparative genomics. The sensitivity and specificity of NOVOMIR for detection of premiRNAs from Arabidopsis thaliana is ~0.83 and ~0.99, respectively. Plant pre-miRNAs are more heterogeneous with respect to size and structure than animal pre-miRNAs. Despite these difficulties, NOVOMIR is well suited to perform searches for pre-miRNAs on a genomic scale. NOVOMIR is written in Perl and relies on two additional, free programs for prediction of RNA secondary structure (RNALFOLD, RNASHAPES).

Multiple distinct small RNAs originate from the same microRNA precursors

Background

MicroRNAs (miRNAs), which originate from precursor transcripts with stem-loop structures, are essential gene expression regulators in eukaryotes.

Results

We report 19 miRNA precursors in Arabidopsis that can yield multiple distinct miRNA-like RNAs in addition to miRNAs and miRNA*s. These miRNA precursor-derived miRNA-like RNAs are often arranged in phase and form duplexes with an approximately two-nucleotide 3'-end overhang. Their production depends on the same biogenesis pathway as their sibling miRNAs and does not require RNA-dependent RNA polymerases or RNA polymerase IV. These miRNA-like RNAs are methylated, and many of them are associated with Argonaute proteins. Some of the miRNA-like RNAs are differentially expressed in response to bacterial challenges, and some are more abundant than the cognate miRNAs. Computational and expression analyses demonstrate that some of these miRNA-like RNAs are potentially functional and they target protein-coding genes for silencing. The function of some of these miRNA-like RNAs was further supported by their target cleavage products from the published small RNA degradome data. Our systematic examination of public small-RNA deep sequencing data from four additional plant species (Oryza sativa, Physcomitrella patens, Medicago truncatula and Populus trichocarpa) and four animals (Homo sapiens, Mus musculus, Caenorhabditis elegans and Drosophila) shows that such miRNA-like RNAs exist broadly in eukaryotes.

Conclusions

We demonstrate that multiple miRNAs could derive from miRNA precursors by sequential processing of Dicer or Dicer-like proteins. Our results suggest that the pool of miRNAs is larger than was previously recognized, and miRNA-mediated gene regulation may be broader and more complex than previously thought.

Transcription factors in light and circadian clock signaling networks revealed by genomewide mapping of direct targets for neurospora white collar complex

Light signaling pathways and circadian clocks are inextricably linked and have profound effects on behavior in most organisms. Here, we used chromatin immunoprecipitation (ChIP) sequencing to uncover direct targets of the Neurospora crassa circadian regulator White Collar Complex (WCC). The WCC is a blue-light receptor and the key transcription factor of the circadian oscillator. It controls a transcriptional network that regulates ∼20% of all genes, generating daily rhythms and responses to light. We found that in response to light, WCC binds to hundreds of genomic regions, including the promoters of previously identified clock- and light-regulated genes. We show that WCC directly controls the expression of 24 transcription factor genes, including the clock-controlled adv-1 gene, which controls a circadian output pathway required for daily rhythms in development. Our findings provide links between the key circadian activator and effectors in downstream regulatory pathways.

Whole-genome in-silico subtractive hybridization (WISH)--using massive sequencing for the identification of unique and repetitive sex-specific sequences: the example of Schistosoma mansoni

Background

Emerging methods of massive sequencing that allow for rapid re-sequencing of entire genomes at comparably low cost are changing the way biological questions are addressed in many domains. Here we propose a novel method to compare two genomes (genome-to-genome comparison). We used this method to identify sex-specific sequences of the human blood fluke Schistosoma mansoni.

Results

Genomic DNA was extracted from male and female (heterogametic) S. mansoni adults and sequenced with a Genome Analyzer (Illumina). Sequences are available at the NCBI sequence read archive http://www.ncbi.nlm.nih.gov/Traces/sra/ under study accession number SRA012151.6. Sequencing reads were aligned to the genome, and a pseudogenome composed of known repeats. Straightforward comparative bioinformatics analysis was performed to compare male and female schistosome genomes and identify female-specific sequences. We found that the S. mansoni female W chromosome contains only few specific unique sequences (950 Kb i.e. about 0.2% of the genome). The majority of W-specific sequences are repeats (10.5 Mb i.e. about 2.5% of the genome). Arbitrarily selected W-specific sequences were confirmed by PCR. Primers designed for unique and repetitive sequences allowed to reliably identify the sex of both larval and adult stages of the parasite.

Conclusion

Our genome-to-genome comparison method that we call "whole-genome in-silico subtractive hybridization" (WISH) allows for rapid identification of sequences that are specific for a certain genotype (e.g. the heterogametic sex). It can in principle be used for the detection of any sequence differences between isolates (e.g. strains, pathovars) or even closely related species.

Unique functionality of 22-nt miRNAs in triggering RDR6-dependent siRNA biogenesis from target transcripts in Arabidopsis

RNA interference pathways can involve amplification of secondary siRNAs by RNA-dependent RNA polymerases. In plants, RDR6-dependent secondary siRNAs arise from transcripts targeted by some microRNAs (miRNAs). Here, Arabidopsis thaliana secondary siRNAs from mRNA as well as trans-acting siRNAs are shown to be triggered through initial targeting by a 22-nucleotide (nt) miRNA that associates with AGO1. In contrast to canonical 21-nt miRNAs, 22-nt miRNAs primarily arise from foldback precursors containing asymmetric bulges. Using artificial miRNA constructs, conversion of asymmetric foldbacks to symmetric foldbacks resulted in the production of 21-nt forms of miR173, miR472 and miR828. Both 21- and 22-nt forms associated with AGO1 and guided accurate slicer activity, but only 22-nt forms were competent to trigger RDR6-dependent siRNA production from target RNA. These data suggest that AGO1 functions differentially with 21- and 22-nt miRNAs to engage the RDR6-associated amplification apparatus.

mirTools: microRNA profiling and discovery based on high-throughput sequencing

miRNAs are small, non-coding RNA that negatively regulate gene expression at post-transcriptional level, which play crucial roles in various physiological and pathological processes, such as development and tumorigenesis. Although deep sequencing technologies have been applied to investigate various small RNA transcriptomes, their computational methods are far away from maturation as compared to microarray-based approaches. In this study, a comprehensive web server mirTools was developed to allow researchers to comprehensively characterize small RNA transcriptome. With the aid of mirTools, users can: (i) filter low-quality reads and 3/5′ adapters from raw sequenced data; (ii) align large-scale short reads to the reference genome and explore their length distribution; (iii) classify small RNA candidates into known categories, such as known miRNAs, non-coding RNA, genomic repeats and coding sequences; (iv) provide detailed annotation information for known miRNAs, such as miRNA/miRNA*, absolute/relative reads count and the most abundant tag; (v) predict novel miRNAs that have not been characterized before; and (vi) identify differentially expressed miRNAs between samples based on two different counting strategies: total read tag counts and the most abundant tag counts. We believe that the integration of multiple computational approaches in mirTools will greatly facilitate current microRNA researches in multiple ways. mirTools can be accessed at http://centre.bioinformatics.zj.cn/mirtools/ and http://59.79.168.90/mirtools.

Introduction into the analysis of high-throughput-sequencing based epigenome data

Sequencing-based approaches now allow high-resolution, genome-scale investigation of cellular epigenetic landscapes. For example, mapping of open chromatin regions, post-translational histone modifications and DNA methylation across a whole genome is now feasible, and new non-coding regulatory RNAs can be sensitively identified via RNA sequencing. The resulting large-scale data sets promise to contribute towards a more precise and complete understanding of gene regulation and to yield insights into the interplay between genomes and the environment. In this article, I review some of the conceptual issues and currently available software tools for the analysis of sequencing-based whole-genome epigenetics data.

Methodological framework for functional characterization of plant microRNAs

Since the beginning of this century, microRNAs (miRNAs), which are tiny RNA molecules, have become one of the major research topics on gene expression regulation in both animals and plants. The major task of miRNA study is to elucidate how the miRNAs are expressed in vivo, how they exert regulatory effects on their targets, and how they can be qualitatively or quantitatively cloned. For these purposes, the methodology of miRNA study has been developed and significantly improved in recent years. The focus here is on a number of powerful methods for plant miRNA research including bioinformatics tools and experimental approaches being used for upstream or downstream analysis of miRNAs or miRNA cloning. Some discrepancies exist in the miRNA research methodology between plants and animals, for example, 5' modified RACE (Rapid Amplification of cDNA Ends) can be used for cleavage target validation only in plants. However, numerous common methods are shared by these two miRNA research areas. Thus, this review will enhance our understanding of miRNA research methodology in organisms.

Small RNAs in angiosperms: sequence characteristics, distribution and generation

High-throughput sequencing (HTS) has opened up a new era for small RNA (sRNA) exploration. Using HTS data for a global survey of sRNAs in 26 angiosperms, elevated GC contents were detected in the monocots, whereas the 5(')-terminal compositions were quite uniform among the angiosperms. Chromosome-wide distribution patterns of sRNAs were investigated by using scrolling-window analysis. We performed de novo natural antisense transcript (NAT) prediction, and found that the overlapping regions of trans-NATs, but not cis-NATs, were hotspots for sRNA generation. One cis-NAT generates phased natural antisense short interfering RNAs (nat-siRNAs) specifically from flowers in Arabidopsis, while one in rice produces phased nat-siRNAs from grains, suggesting their organ-specific regulatory roles.

MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana

MicroRNAs (miRNAs) are short regulatory RNAs processed from partially self-complementary foldbacks within longer MIRNA primary transcripts. Several MIRNA families are conserved deeply through land plants, but many are present only in closely related species or are species specific. The finding of numerous evolutionarily young MIRNA, many with low expression and few if any targets, supports a rapid birth-death model for MIRNA evolution. A systematic analysis of MIRNA genes and families in the close relatives, Arabidopsis thaliana and Arabidopsis lyrata, was conducted using both whole-genome comparisons and high-throughput sequencing of small RNAs. Orthologs of 143 A. thaliana MIRNA genes were identified in A. lyrata, with nine having significant sequence or processing changes that likely alter function. In addition, at least 13% of MIRNA genes in each species are unique, despite their relatively recent speciation (approximately 10 million years ago). Alignment of MIRNA foldbacks to the Arabidopsis genomes revealed evidence for recent origins of 32 families by inverted or direct duplication of mostly protein-coding gene sequences, but less than half of these yield miRNA that are predicted to target transcripts from the originating gene family. miRNA nucleotide divergence between A. lyrata and A. thaliana orthologs was higher for young MIRNA genes, consistent with reduced purifying selection compared with deeply conserved MIRNA genes. Additionally, target sites of younger miRNA were lost more frequently than for deeply conserved families. In summary, our systematic analyses emphasize the dynamic nature of the MIRNA complement of plant genomes.

SeqBuster, a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells

High-throughput sequencing technologies enable direct approaches to catalog and analyze snapshots of the total small RNA content of living cells. Characterization of high-throughput sequencing data requires bioinformatic tools offering a wide perspective of the small RNA transcriptome. Here we present SeqBuster, a highly versatile and reliable web-based toolkit to process and analyze large-scale small RNA datasets. The high flexibility of this tool is illustrated by the multiple choices offered in the pre-analysis for mapping purposes and in the different analysis modules for data manipulation. To overcome the storage capacity limitations of the web-based tool, SeqBuster offers a stand-alone version that permits the annotation against any custom database. SeqBuster integrates multiple analyses modules in a unique platform and constitutes the first bioinformatic tool offering a deep characterization of miRNA variants (isomiRs). The application of SeqBuster to small-RNA datasets of human embryonic stem cells revealed that most miRNAs present different types of isomiRs, some of them being associated to stem cell differentiation. The exhaustive description of the isomiRs provided by SeqBuster could help to identify miRNA-variants that are relevant in physiological and pathological processes. SeqBuster is available at http://estivill_lab.crg.es/seqbuster.

Exploring plant transcriptomes using ultra high-throughput sequencing

Ultra high-throughput sequencing (UHTS) technologies offer the potential to interrogate transcriptomes in detail that has traditionally been restricted to single gene surveys. For instance, it is now possible to globally define transcription start sites, polyadenylation signals, alternative splice sites and generate quantitative data on gene transcript accumulation in single tissues or cell types. These technologies are thus paving the way for whole genome transcriptomics and will undoubtedly lead to novel insights into plant development and biotic and abiotic stress responses. However, several challenges exist to making this technology broadly accessible to the plant research community. These include the current need for a computationally intensive analysis of data sets, a lack of standardized alignment and formatting procedures and a relatively small number of analytical software packages to interpret UHTS outputs. In this review we summarize recent findings from UHTS and discuss potential opportunities and challenges for broad adoption of these technologies in the plant science community.

Arabidopsis RNA-dependent RNA polymerases and dicer-like proteins in antiviral defense and small interfering RNA biogenesis during Turnip Mosaic Virus infection

Plants respond to virus infections by activation of RNA-based silencing, which limits infection at both the single-cell and system levels. Viruses encode RNA silencing suppressor proteins that interfere with this response. Wild-type Arabidopsis thaliana is immune to silencing suppressor (HC-Pro)-deficient Turnip mosaic virus, but immunity was lost in the absence of DICER-LIKE proteins DCL4 and DCL2. Systematic analysis of susceptibility and small RNA formation in Arabidopsis mutants lacking combinations of RNA-dependent RNA polymerase (RDR) and DCL proteins revealed that the vast majority of virus-derived small interfering RNAs (siRNAs) were dependent on DCL4 and RDR1, although full antiviral defense also required DCL2 and RDR6. Among the DCLs, DCL4 was sufficient for antiviral silencing in inoculated leaves, but DCL2 and DCL4 were both involved in silencing in systemic tissues (inflorescences). Basal levels of antiviral RNA silencing and siRNA biogenesis were detected in mutants lacking RDR1, RDR2, and RDR6, indicating an alternate route to form double-stranded RNA that does not depend on the three previously characterized RDR proteins.

Loqs and R2D2 act sequentially in the siRNA pathway in Drosophila

In Drosophila melanogaster, the small interfering RNA (siRNA) pathway is triggered by exogenous double-stranded RNA (dsRNA) or upon viral infection. This pathway requires Dicer-2 (Dcr-2) in association with a dsRNA-binding protein (dsRBP) called R2D2. A potentially distinct siRNA pathway, which requires Dcr-2 in association with a different dsRBP, called Loquacious (Loqs), is activated by endogenous dsRNA derived from transposons, structured loci and overlapping transcripts. Here we show that different sources of dsRNA enter a common siRNA pathway that requires R2D2 and Loqs. R2D2 and loqs mutants show impaired silencing triggered by injection of exogenous dsRNA or by artificial and natural expression of endogenous dsRNA. In addition, we show that these dsRBPs function sequentially and nonredundantly in collaboration with Dcr-2. Loqs is primarily required for dsRNA processing, whereas R2D2 is essential for the subsequent loading of siRNAs into effector Ago-RISC complexes.

Identification of MIR390a precursor processing-defective mutants in Arabidopsis by direct genome sequencing

Transacting siRNA (tasiRNA) biogenesis in Arabidopsis is initiated by microRNA (miRNA) -guided cleavage of primary transcripts. In the case of TAS3 tasiRNA formation, ARGONAUTE7 (AGO7)-miR390 complexes interact with primary transcripts at two sites, resulting in recruitment of RNA-DEPENDENT RNA POLYMERASE6 for dsRNA biosynthesis. An extensive screen for Arabidopsis mutants with specific defects in TAS3 tasiRNA biogenesis or function was done. This yielded numerous ago7 mutants, one dcl4 mutant, and two mutants that accumulated low levels of miR390. A direct genome sequencing-based approach to both map and rapidly identify one of the latter mutant alleles was developed. This revealed a G-to-A point mutation (mir390a-1) that was calculated to stabilize a relatively nonpaired region near the base of the MIR390a foldback, resulting in misprocessing of the miR390/miR390* duplex and subsequent reduced TAS3 tasiRNA levels. Directed substitutions, as well as analysis of variation at paralogous miR390-generating loci (MIR390a and MIR390b), indicated that base pair properties and nucleotide identity within a region 4-6 bases below the miR390/miR390* duplex region contributed to the efficiency and accuracy of precursor processing.

Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing

Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5' end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(-) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs.

Deciphering the diversity of small RNAs in plants: the long and short of it

RNA silencing is a complex and highly conserved regulatory mechanism that is now known to be involved in such diverse processes as development, pathogen control, genome maintenance and response to environmental changes. Since its recent discovery, RNA silencing has become a fast moving key area of research in plant and animal molecular biology. Research in this field has greatly profited from recent developments in novel sequencing technologies that allow massive parallel sequencing of small RNA (sRNA) molecules, the key players of all RNA silencing phenomena. As researchers are beginning to decipher the complexity of RNA silencing, novel methodologies have to be developed to make sense of the large amounts of data that are currently being generated. In this review we present an overview of RNA silencing pathways in plants and the current challenges in analysing sRNA data, with a special focus on computational approaches.