Irina Mohorianu

The UEA sRNA workbench

Summary: RNA silencing is a complex, highly conserved mechanism mediated by small RNAs (sRNAs), such as microRNAs (miRNAs), that is known to be involved in a diverse set of biological functions including development, pathogen control, genome maintenance and response to environmental change. Advances in next generation sequencing technologies are producing increasingly large numbers of sRNA reads per sample at a fraction of the cost of previous methods. However, many bioinformatics tools do not scale accordingly, are cumbersome, or require extensive support from bioinformatics experts. Therefore, researchers need user-friendly, robust tools, capable of not only processing large sRNA datasets in a reasonable time frame but also presenting the results in an intuitive fashion and visualizing sRNA genomic features. Herein, we present the UEA sRNA workbench, a suite of tools that is a successor to the web-based UEA sRNA Toolkit, but in downloadable format and with several enhanced and additional features. Availability: The program and help pages are available at http://srna-workbench.cmp.uea.ac.uk. Contact: vincent.moulton@cmp.uea.ac.uk

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.

Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions.

MicroRNAs negatively regulate the accumulation of mRNAs therefore when they are expressed in the same cells their expression profiles show an inverse correlation. We previously described one positively correlated miRNA/target pair, but it is not known how widespread this phenomenon is. Here, we investigated the correlation between the expression profiles of differentially expressed miRNAs and their targets during tomato fruit development using deep sequencing, Northern blot and RT-qPCR. We found an equal number of positively and negatively correlated miRNA/target pairs indicating that positive correlation is more frequent than previously thought. We also found that the correlation between microRNA and target expression profiles can vary between mRNAs belonging to the same gene family and even for the same target mRNA at different developmental stages. Since microRNAs always negatively regulate their targets, the high number of positively correlated microRNA/target pairs suggests that mutual exclusion could be as widespread as temporal regulation. The change of correlation during development suggests that the type of regulatory circuit directed by a microRNA can change over time and can be different for individual gene family members. Our results also highlight potential problems for expression profiling-based microRNA target identification/validation.

Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant

BackgroundThe 7B-1 tomato line (Solanum lycopersicum cv. Rutgers) is a photoperiod-sensitive male-sterile mutant, with potential application in hybrid seed production. Small RNAs (sRNAs) in tomato have been mainly characterized in fruit development and ripening, but none have been studied with respect to flower development and regulation of male-sterility. Using sRNA sequencing, we identified miRNAs that are potentially involved in anther development and regulation of male-sterility in 7B-1 mutant.ResultsTwo sRNA libraries from 7B-1 and wild type (WT) anthers were sequenced and thirty two families of known miRNAs and 23 new miRNAs were identified in both libraries. MiR390, miR166, miR159 were up-regulated and miR530, miR167, miR164, miR396, miR168, miR393, miR8006 and two new miRNAs, miR#W and miR#M were down-regulated in 7B-1 anthers. Ta-siRNAs were not differentially expressed and likely not associated with 7B-1 male-sterility. miRNA targets with potential roles in anther development were validated using 5′-RACE. QPCR analysis showed differential expression of miRNA/target pairs of interest in anthers and stem of 7B-1, suggesting that they may regulate different biological processes in these tissues. Expression level of most miRNA/target pairs showed negative correlation, except for few. In situ hybridization showed predominant expression of miR159, GAMYBL1, PMEI and cystatin in tapetum, tetrads and microspores.ConclusionOverall, we identified miRNAs with potential roles in anther development and regulation of male-sterility in 7B-1. A number of new miRNAs were also identified from tomato for the first time. Our data could be used as a benchmark for future studies of the molecular mechanisms of male-sterility in other crops.

Small RNA Profile in Moso Bamboo Root and Leaf Obtained by High Definition Adapters

Moso bamboo (Phyllostachy heterocycla cv. pubescens L.) is an economically important fast-growing tree. In order to gain better understanding of gene expression regulation in this important species we used next generation sequencing to profile small RNAs in leaf and roots of young seedlings. Since standard kits to produce cDNA of small RNAs are biased for certain small RNAs, we used High Definition adapters that reduce ligation bias. We identified and experimentally validated five new microRNAs and a few other small non-coding RNAs that were not microRNAs. The biological implication of microRNA expression levels and targets of microRNAs are discussed.

Global discovery and characterization of small non-coding RNAs in marine microalgae

BackgroundMarine phytoplankton are responsible for 50% of the CO2 that is fixed annually worldwide and contribute massively to other biogeochemical cycles in the oceans. Diatoms and coccolithophores play a significant role as the base of the marine food web and they sequester carbon due to their ability to form blooms and to biomineralise. To discover the presence and regulation of short non-coding RNAs (sRNAs) in these two important phytoplankton groups, we sequenced short RNA transcriptomes of two diatom species (Thalassiosira pseudonana, Fragilariopsis cylindrus) and validated them by Northern blots along with the coccolithophore Emiliania huxleyi.ResultsDespite an exhaustive search, we did not find canonical miRNAs in diatoms. The most prominent classes of sRNAs in diatoms were repeat-associated sRNAs and tRNA-derived sRNAs. The latter were also present in E. huxleyi. tRNA-derived sRNAs in diatoms were induced under important environmental stress conditions (iron and silicate limitation, oxidative stress, alkaline pH), and they were very abundant especially in the polar diatom F. cylindrus (20.7% of all sRNAs) even under optimal growth conditions.ConclusionsThis study provides first experimental evidence for the existence of short non-coding RNAs in marine microalgae. Our data suggest that canonical miRNAs are absent from diatoms. However, the group of tRNA-derived sRNAs seems to be very prominent in diatoms and coccolithophores and maybe used for acclimation to environmental conditions.

Genomic responses to socio-sexual environment in male Drosophila melanogaster exposed to conspecific rivals

Socio-sexual environments have profound effects on fitness. Local sex ratios can alter the threat of sexual competition, to which males respond via plasticity in reproductive behaviors and ejaculate composition. In Drosophila melanogaster, males detect the presence of conspecific, same-sex mating rivals prior to mating using multiple, redundant sensory cues. Males that respond to rivals gain significant fitness benefits by altering mating duration and ejaculate composition. Here we investigated the underlying genome-wide changes involved. We used RNA-seq to analyze male transcriptomic responses 2, 26, and 50 h after exposure to rivals, a time period that was previously identified as encompassing the major facets of male responses to rivals. The results showed a strong early activation of multiple sensory genes in the head-thorax (HT), prior to the expression of any phenotypic differences. This gene expression response was reduced by 26 h, at the time of maximum phenotypic change, and shut off by 50 h. In the abdomen (A), fewer genes changed in expression and gene expression responses appeared to increase over time. The results also suggested that different sets of functionally equivalent genes might be activated in different replicates. This could represent a mechanism by which robustness is conferred upon highly plastic traits. Overall, our study reveals that mRNA-seq can identify subtle genomic signatures characteristic of flexible behavioral phenotypes.

The Pestivirus N Terminal Protease N pro Redistributes to Mitochondria and Peroxisomes Suggesting New Sites for Regulation of IRF3 by N pro

The N-terminal protease of pestiviruses, Npro is a unique viral protein, both because it is a distinct autoprotease that cleaves itself from the following polyprotein chain, and also because it binds and inactivates IRF3, a central regulator of interferon production. An important question remains the role of Npro in the inhibition of apoptosis. In this study, apoptotic signals induced by staurosporine, interferon, double stranded RNA, sodium arsenate and hydrogen peroxide were inhibited by expression of wild type Npro, but not by mutant protein Npro C112R, which we show is less efficient at promoting degradation of IRF3, and led to the conclusion that Npro inhibits the stress-induced intrinsic mitochondrial pathway through inhibition of IRF3-dependent Bax activation. Both expression of Npro and infection with Bovine Viral Diarrhea Virus (BVDV) prevented Bax redistribution and mitochondrial fragmentation. Given the role played by signaling platforms during IRF3 activation, we have studied the subcellular distribution of Npro and we show that, in common with many other viral proteins, Npro targets mitochondria to inhibit apoptosis in response to cell stress. Npro itself not only relocated to mitochondria but in addition, both Npro and IRF3 associated with peroxisomes, with over 85% of Npro puncta co-distributing with PMP70, a marker for peroxisomes. In addition, peroxisomes containing Npro and IRF3 associated with ubiquitin. IRF3 was degraded, whereas Npro accumulated in response to cell stress. These results implicate mitochondria and peroxisomes as new sites for IRF3 regulation by Npro, and highlight the role of these organelles in the anti-viral pathway.

Small RNA Analysis in Sindbis Virus Infected Human HEK293 Cells

Introduction In contrast to the defence mechanism of RNA interference (RNAi) in plants and invertebrates, its role in the innate response to virus infection of mammals is a matter of debate. Since RNAi has a well-established role in controlling infection of the alphavirus Sindbis virus (SINV) in insects, we have used this virus to investigate the role of RNAi in SINV infection of human cells. Results SINV AR339 and TR339-GFP were adapted to grow in HEK293 cells. Deep sequencing of small RNAs (sRNAs) early in SINV infection (4 and 6 hpi) showed low abundance (0.8%) of viral sRNAs (vsRNAs), with no size, sequence or location specific patterns characteristic of Dicer products nor did they possess any discernible pattern to ascribe to a specific RNAi biogenesis pathway. This was supported by multiple variants for each sequence, and lack of hot spots along the viral genome sequence. The abundance of the best defined vsRNAs was below the limit of Northern blot detection. The adaptation of the virus to HEK293 cells showed little sequence changes compared to the reference; however, a SNP in E1 gene with a preference from G to C was found. Deep sequencing results showed little variation of expression of cellular microRNAs (miRNAs) at 4 and 6 hpi compared to uninfected cells. Twelve miRNAs exhibiting some minor differential expression by sequencing, showed no difference in expression by Northern blot analysis. Conclusions We show that, unlike SINV infection of invertebrates, generation of Dicer-dependent svRNAs and change in expression of cellular miRNAs were not detected as part of the Human response to SINV.

FiRePat—Finding Regulatory Patterns between sRNAs and Genes

Small RNAs are regulatory RNA fragments which, through RNA silencing, can regulate the expression of genes. Because sRNAs are negative regulators it is generally assumed that expression profiles of sRNAs and their targets are negatively correlated. Recently, examples of positive correlation between the expression of sRNAs and their targets have been discovered. At the moment, it is not known how many sRNA‐target pairs are positively and negatively correlated, and it is also not clear in what situations (e.g., under which treatments) any of these correlations can be observed. To determine this, one of the first steps is to develop tools to carry out a genome wide characterization of covariation of expression levels of sRNAs and genes. We present FiRePat—Finding Regulatory Patterns—an unsupervised data mining tool applicable to large datasets, typically produced by high throughput sequencing of sRNAs and mRNAs or microarray experiments, that detects sRNA‐gene pairs with correlated expression levels. The method consists of three steps: first, we select differentially expressed sRNAs and genes; second, we compute the correlation between sRNA and gene series for all possible sRNA–gene pairs; and third, we cluster the sRNA or gene expression series, simultaneously inducing clusters in the other series. Potential uses of FiRePat are presented using publicly available sRNA and mRNA datasets for both plants and animals. The standard output of FiRePat, a list of correlated pairs formed with sRNAs and mRNAs, can be used to investigate the cause and consequences of the respective expression patterns.