Vincent Moulton

RDP3: a flexible and fast computer program for analyzing recombination

Summary: RDP3 is a new version of the RDP program for characterizing recombination events in DNA-sequence alignments. Among other novelties, this version includes four new recombination analysis methods (3SEQ, VISRD, PHYLRO and LDHAT), new tests for recombination hot-spots, a range of matrix methods for visualizing over-all patterns of recombination within datasets and recombination-aware ancestral sequence reconstruction. Complementary to a high degree of analysis flow automation, RDP3 also has a highly interactive and detailed graphical user interface that enables more focused hands-on cross-checking of results with a wide variety of newly implemented phylogenetic tree construction and matrix-based recombination signal visualization methods. The new RDP3 can accommodate large datasets and is capable of analyzing alignments ranging in size from 1000×10 kilobase sequences to 20×2 megabase sequences within 48 h on a desktop PC. Availability: RDP3 is available for free from its web site http://darwin.uvigo.es/rdp/rdp.html Contact: darrenpatrickmartin@gmail.com Supplementary information: The RDP3 program manual contains detailed descriptions of the various methods it implements and a step-by-step guide describing how best to use these.

Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening

In plants there are several classes of 21-24-nt short RNAs that regulate gene expression. The most conserved class is the microRNAs (miRNAs), although some miRNAs are found only in specific species. We used high-throughput pyrosequencing to identify conserved and nonconserved miRNAs and other short RNAs in tomato fruit and leaf. Several conserved miRNAs showed tissue-specific expression, which, combined with target gene validation results, suggests that miRNAs may play a role in fleshy fruit development. We also identified four new nonconserved miRNAs. One of the validated targets of a novel miRNA is a member of the CTR family involved in fruit ripening. However, 62 predicted targets showing near perfect complementarity to potential new miRNAs did not validate experimentally. This suggests that target prediction of plant short RNAs could have a high false-positive rate and must therefore be validated experimentally. We also found short RNAs from a Solanaceae-specific foldback transposon, which showed a miRNA/miRNA*-like distribution, suggesting that this element may function as a miRNA gene progenitor. The other Solanaceae-specific class of short RNA was derived from an endogenous pararetrovirus sequence inserted into the tomato chromosomes. This study opens a new avenue in the field of fleshy fruit biology by raising the possibility that fruit development and ripening may be under miRNA regulation.

Identification of grapevine microRNAs and their targets using high throughput sequencing and degradome analysis

In plants, microRNAs (miRNAs) comprise one of three classes of small RNAs regulating gene expression at the post-transcriptional level. Many plant miRNAs are conserved, and play a role in development, abiotic stress responses or pathogen responses. However, some miRNAs have only been found in certain species. Here, we use deep-sequencing, computational and molecular methods to identify, profile, and describe conserved and non-conserved miRNAs in four grapevine (Vitis vinifera) tissues. A total of 24 conserved miRNA families were identified in all four tissues, and 26 known but non-conserved miRNAs were also found. In addition to known miRNAs, we also found 21 new grapevine-specific miRNAs together with their star strands. We have also shown that almost all of them originated from single genes. Furthermore, 21 other plausible miRNA candidates have been described. We have found that many known and new miRNAs showed tissue-specific expression. Finally, 112 target mRNAs of known and 44 target mRNAs of new grapevine-specific miRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs.

High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families

BackgroundHigh-throughput sequencing technology is capable to identify novel short RNAs in plant species. We used Solexa sequencing to find new microRNAs in one of the model legume species, barrel medic (Medicago truncatula).Results3,948,871 reads were obtained from two separate short RNA libraries generated from total RNA extracted from M. truncatula leaves, representing 1,563,959 distinct sequences. 2,168,937 reads were mapped to the available M. truncatula genome corresponding to 619,175 distinct sequences. 174,504 reads representing 25 conserved miRNA families showed perfect matches to known miRNAs. We also identified 26 novel miRNA candidates that were potentially generated from 32 loci. Nine of these loci produced eight distinct sequences, for which the miRNA* sequences were also sequenced. These sequences were not described in other plant species and accumulation of these eight novel miRNAs was confirmed by Northern blot analysis. Potential target genes were predicted for most conserved and novel miRNAs.ConclusionDeep sequencing of short RNAs from M. truncatula leaves identified eight new miRNAs indicating that specific miRNAs exist in legume species.

A toolkit for analysing large-scale plant small RNA datasets

UNLABELLED Recent developments in high-throughput sequencing technologies have generated considerable demand for tools to analyse large datasets of small RNA sequences. Here, we describe a suite of web-based tools for processing plant small RNA datasets. Our tools can be used to identify micro RNAs and their targets, compare expression levels in sRNA loci, and find putative trans-acting siRNA loci. AVAILABILITY The tools are freely available for use at http://srna-tools.cmp.uea.ac.uk.

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

NeighborNet: An Agglomerative Method for the Construction of Planar Phylogenetic Networks

We introduce NeighborNet, a network construction and data representation method that combines aspects of the neighbor joining (NJ) and SplitsTree. Like NJ, NeighborNet uses agglomeration: taxa are combined into progressively larger and larger overlapping clusters. Like SPLITSTREE, NeighborNet constructs networks rather than trees, and so can be used to represent multiple phylogenetic hypotheses simultaneously, or to detect complex evolutionary processes like recombination, lateral transfer and hybridization. NeighborNet tends to produce networks that are substantially more resolved than those made with SPLITSTREE. The method is efficient (O(n3) time) and is well suited for the preliminary analyses of complex phylogenetic data. We report results of three case studies: one based on mitochondrial gene order data from early branching eukaryotes, another based on nuclear sequence data from New Zealand alpine buttercups (Ranunculi), and a third on poorly corrected synthetic data.

T -theory: an overview

Abstract T-theory is the name that we adopt for the theory of trees, injective envelopes of metric spaces, and all of the areas that are connected with these topics, which has been developed over the past 10–15 years in Bielefeld. Its motivation was originally—and still is to a large extent—the development of mathematical tools for reconstructing phylogenetic trees.T-theory expanded considerably when its relationships with the theory of affine buildings, valuated matroids, and decompositions of metrics were discovered. In this paper, we give a brief introduction to this theory, which we hope will serve as a useful reference to some of the main results, and also as a guide for further investigations into whatT-theory has to offer.

Reducing ligation bias of small RNAs in libraries for next generation sequencing

BackgroundThe use of nucleic acid-modifying enzymes has driven the rapid advancement in molecular biology. Understanding their function is important for modifying or improving their activity. However, functional analysis usually relies upon low-throughput experiments. Here we present a method for functional analysis of nucleic acid-modifying enzymes using next generation sequencing.FindingsWe demonstrate that sequencing data of libraries generated by RNA ligases can reveal novel secondary structure preferences of these enzymes, which are used in small RNA cloning and library preparation for NGS. Using this knowledge we demonstrate that the cloning bias in small RNA libraries is RNA ligase-dependent. We developed a high definition (HD) protocol that reduces the RNA ligase-dependent cloning bias. The HD protocol doubled read coverage, is quantitative and found previously unidentified microRNAs. In addition, we show that microRNAs in miRBase are those preferred by the adapters of the main sequencing platform.ConclusionsSequencing bias of small RNAs partially influenced which microRNAs have been studied in depth; therefore most previous small RNA profiling experiments should be re-evaluated. New microRNAs are likely to be found, which were selected against by existing adapters. Preference of currently used adapters towards known microRNAs suggests that the annotation of all existing small RNAs, including miRNAs, siRNAs and piRNAs, has been biased.

Structural and functional analysis of viral siRNAs

A large amount of short interfering RNA (vsiRNA) is generated from plant viruses during infection, but the function, structure and biogenesis of these is not understood. We profiled vsiRNAs using two different high-throughput sequencing platforms and also developed a hybridisation based array approach. The profiles obtained through the Solexa platform and by hybridisation were very similar to each other but different from the 454 profile. Both deep sequencing techniques revealed a strong bias in vsiRNAs for the positive strand of the virus and identified regions on the viral genome that produced vsiRNA in much higher abundance than other regions. The hybridisation approach also showed that the position of highly abundant vsiRNAs was the same in different plant species and in the absence of RDR6. We used the Terminator 5′-Phosphate-Dependent Exonuclease to study the 5′ end of vsiRNAs and showed that a perfect control duplex was not digested by the enzyme without denaturation and that the efficiency of the Terminator was strongly affected by the concentration of the substrate. We found that most vsiRNAs have 5′ monophosphates, which was also confirmed by profiling short RNA libraries following either direct ligation of adapters to the 5′ end of short RNAs or after replacing any potential 5′ ends with monophosphates. The Terminator experiments also showed that vsiRNAs were not perfect duplexes. Using a sensor construct we also found that regions from the viral genome that were complementary to non-abundant vsiRNAs were targeted in planta just as efficiently as regions recognised by abundant vsiRNAs. Different high-throughput sequencing techniques have different reproducible sequence bias and generate different profiles of short RNAs. The Terminator exonuclease does not process double stranded RNA, and because short RNAs can quickly re-anneal at high concentration, this assay can be misleading if the substrate is not denatured and not analysed in a dilution series. The sequence profiles and Terminator digests suggest that CymRSV siRNAs are produced from the structured positive strand rather than from perfect double stranded RNA or by RNA dependent RNA polymerase.