Mapping rRNA 2'-O-methylations and identification of C/D snoRNAs in Arabidopsis thaliana plants.

In all eukaryotic cells, the most abundant modification of ribosomal RNA (rRNA) is methylation at the ribose moiety (2'-O-methylation). Ribose methylation at specific rRNA sites is guided by small nucleolar RNAs (snoRNAs) of C/D-box type (C/D snoRNA) and achieved by the methyltransferase Fibrillarin (FIB). Here we used the Illumina-based RiboMethSeq approach for mapping rRNA 2'-O-methylation sites in A. thaliana Col-0 (WT) plants. This analysis detected novel C/D snoRNA-guided rRNA 2'-O-methylation positions and also some orphan sites without a matching C/D snoRNA. Furthermore, immunoprecipitation of Arabidopsis FIB2 identified and demonstrated expression of C/D snoRNAs corresponding to majority of mapped rRNA sites. On the other hand, we show that disruption of Arabidopsis Nucleolin 1 gene (NUC1), encoding a major nucleolar protein, decreases 2'-O-methylation at specific rRNA sites suggesting functional/structural interconnections of 2'-O-methylation with nucleolus organization and plant development. Finally, based on our findings and existent database sets, we introduce a new nomenclature system for C/D snoRNA in Arabidopsis plants.

Searching RNA motifs and their intermolecular contacts with constraint networks.

MOTIVATION: Searching RNA gene occurrences in genomic sequences is a task whose importance has been renewed by the recent discovery of numerous functional RNA, often interacting with other ligands. Even if several programs exist for RNA motif search, none exists that can represent and solve the problem of searching for occurrences of RNA motifs in interaction with other molecules. RESULTS: We present a constraint network formulation of this problem. RNA are represented as structured motifs that can occur on more than one sequence and which are related together by possible hybridization. The implemented tool MilPat is used to search for several sRNA families in genomic sequences. Results show that MilPat allows to efficiently search for interacting motifs in large genomic sequences and offers a simple and extensible framework to solve such problems. New and known sRNA are identified as H/ACA candidates in Methanocaldococcus jannaschii. AVAILABILITY: http://carlit.toulouse.inra.fr/MilPaT/MilPat.pl.

Genome sequence of the plant pathogen Ralstonia solanacearum.

Ralstonia solanacearum is a devastating, soil-borne plant pathogen with a global distribution and an unusually wide host range. It is a model system for the dissection of molecular determinants governing pathogenicity. We present here the complete genome sequence and its analysis of strain GMI1000. The 5.8-megabase (Mb) genome is organized into two replicons: a 3.7-Mb chromosome and a 2.1-Mb megaplasmid. Both replicons have a mosaic structure providing evidence for the acquisition of genes through horizontal gene transfer. Regions containing genetically mobile elements associated with the percentage of G+C bias may have an important function in genome evolution. The genome encodes many proteins potentially associated with a role in pathogenicity. In particular, many putative attachment factors were identified. The complete repertoire of type III secreted effector proteins can be studied. Over 40 candidates were identified. Comparison with other genomes suggests that bacterial plant pathogens and animal pathogens harbour distinct arrays of specialized type III-dependent effectors.

Box C/D RNA guides for the ribose methylation of archaeal tRNAs. The tRNATrp intron guides the formation of two ribose-methylated nucleosides in the mature tRNATrp.

Following a search of the Pyrococcus genomes for homologs of eukaryotic methylation guide small nucleolar RNAs, we have experimentally identified in Pyrococcus abyssi four novel box C/D small RNAs predicted to direct 2'-O-ribose methylations onto the first position of the anticodon in tRNALeu(CAA), tRNALeu(UAA), elongator tRNAMet and tRNATrp, respectively. Remarkably, one of them corresponds to the intron of its presumptive target, pre-tRNATrp. This intron is predicted to direct in cis two distinct ribose methylations within the unspliced tRNA precursor, not only onto the first position of the anticodon in the 5' exon but also onto position 39 (universal tRNA numbering) in the 3' exon. The two intramolecular RNA duplexes expected to direct methylation, which both span an exon-intron junction in pre-tRNATrp, are phylogenetically conserved in euryarchaeotes. We have experimentally confirmed the predicted guide function of the box C/D intron in halophile Haloferax volcanii by mutagenesis analysis, using an in vitro splicing/RNA modification assay in which the two cognate ribose methylations of pre-tRNATrp are faithfully reproduced. Euryarchaeal pre-tRNATrp should provide a unique system to further investigate the molecular mechanisms of RNA-guided ribose methylation and gain new insights into the origin and evolution of the complex family of archaeal and eukaryotic box C/D small RNAs.

Identification of 66 box C/D snoRNAs in Arabidopsis thaliana: extensive gene duplications generated multiple isoforms predicting new ribosomal RNA 2'-O-methylation sites.

Dozens of box C/D small nucleolar RNAs (snoRNAs) have recently been found in eukaryotes (vertebrates, yeast), ancient eukaryotes (trypanosomes) and archae, that specifically target ribosomal RNA sites for 2'-O-ribose methylation. Although early biochemical data revealed that plant rRNAs are among the most highly ribomethylated in eukaryotes, only a handful of methylation guide snoRNAs have been characterized in this kingdom. We report 66 novel box C/D snoRNAs identified by computational screening of Arabidopsis genomic sequences that are expressed in vivo from either single genes, 17 different clusters or three introns. At the structural level, many box C/D snoRNAs have dual antisense elements often matching rRNA regions close to each other on the rRNA secondary structure, which is reminiscent of their archaeal counterparts. Remarkable specimens are found that display two antisense elements having the potential to form an extended snoRNA-rRNA duplex of 23 to 30 nt, in line with the hypothetical function of box C/D snoRNAs in pre-rRNA folding or chaperoning. In contrast to other species, many Arabidopsis snoRNAs are found in multiple isoforms mainly resulting from two different mechanisms: large chromosomal duplications and small tandem duplications producing polycistronic genes. The discovery of numerous different snoRNAs, some of them arising from common ancestors, provide new insights to understand snoRNAs evolution and the birth of new rRNA methylation sites in plants and other organisms.

Archaeal homologs of eukaryotic methylation guide small nucleolar RNAs: lessons from the Pyrococcus genomes.

Ribose methylation is a prevalent type of nucleotide modification in rRNA. Eukaryotic rRNAs display a complex pattern of ribose methylations, amounting to 55 in yeast Saccharomyces cerevisiae and about 100 in vertebrates. Ribose methylations of eukaryotic rRNAs are each guided by a cognate small RNA, belonging to the family of box C/D antisense snoRNAs, through transient formation of a specific base-pairing at the rRNA modification site. In prokaryotes, the pattern of rRNA ribose methylations has been fully characterized in a single species so far, Escherichia coli, which contains only four ribose methylated rRNA nucleotides. However, the hyperthermophile archaeon Sulfolobus solfataricus contains, like eukaryotes, a large number of (yet unmapped) rRNA ribose methylations and homologs of eukaryotic box C/D small nucleolar ribonuclear proteins have been identified in archaeal genomes. We have therefore searched archaeal genomes for potential homologs of eukaryotic methylation guide small nucleolar RNAs, by combining searches for structured motifs with homology searches. We have identified a family of 46 small RNAs, conserved in the genomes of three hyperthermophile Pyrococcus species, which we have experimentally characterized in Pyrococcus abyssi. The Pyrococcus small RNAs, the first reported homologs of methylation guide small nucleolar RNAs in organisms devoid of a nucleus, appear as a paradigm of minimalist box C/D antisense RNAs. They differ from their eukaryotic homologs by their outstanding structural homogeneity, extended consensus box motifs and the quasi-systematic presence of two (instead of one) rRNA antisense elements. Remarkably, for each small RNA the two antisense elements always match rRNA sequences close to each other in rRNA structure, suggesting an important role in rRNA folding. Only a few of the predicted P. abyssi rRNA ribose methylations have been detected so far. Further analysis of these archaeal small RNAs could provide new insights into the origin and functions of methylation guide small nucleolar RNAs and illuminate the still elusive role of rRNA ribose methylations.

A composite genetic map of the parasitoid wasp Trichogramma brassicae based on RAPD markers.

Three linkage maps of the genome of the microhymenopteran Trichogramma brassicae were constructed from the analysis of segregation of random amplified polymorphic DNA markers in three F2 populations. These populations were composed of the haploid male progeny of several virgin F1 females, which resulted from the breeding of four parental lines that were nearly fixed for different random amplified polymorphic DNA markers and that were polymorphic for longevity and fecundity characters. As the order of markers common to the three mapping populations was found to be well conserved, a composite linkage map was constructed. Eighty-four markers were organized into five linkage groups and two pairs. The mean interval between two markers was 17.7 cM, and the map spanned 1330 cM.

ESSA: an integrated and interactive computer tool for analysing RNA secondary structure.

With ESSA, we propose an approach of RNA secondary structure analysis based on extensive viewing within a friendly graphical interface. This computer program is organized around the display of folding models produced by two complementary methods suitable to draw long RNA molecules. Any feature of interest can be managed directly on the display and highlighted by a rich combination of colours and symbols with emphasis given to structural probe accessibilities. ESSA also includes a word searching procedure allowing easy visual identification of structural features even complex and degenerated. Analysis functions make it possible to calculate the thermodynamic stability of any part of a folding using several models and compare homologous aligned RNA both in primary and secondary structure. The predictive capacities of ESSA which brings together the experimental, thermodynamic and comparative methods, are increased by coupling it with a program dedicated to RNA folding prediction based on constraints management and propagation. The potentialities of ESSA are illustrated by the identification of a possible tertiary motif in the LSU rRNA and the visualization of a pseudoknot in S15 mRNA.

CARTHAGENE: constructing and joining maximum likelihood genetic maps.

Genetic mapping is an important step in the study of any organism. An accurate genetic map is extremely valuable for locating genes or more generally either qualitative or quantitative trait loci (QTL). This paper presents a new approach to two important problems in genetic mapping: automatically ordering markers to obtain a multipoint maximum likelihood map and building a multipoint maximum likelihood map using pooled data from several crosses. The approach is embodied in an hybrid algorithm that mixes the statistical optimization algorithm EM with local search techniques which have been developed in the artificial intelligence and operations research communities. An efficient implementation of the EM algorithm provides maximum likelihood recombination fractions, while the local search techniques look for orders that maximize this maximum likelihood. The specificity of the approach lies in the neighborhood structure used in the local search algorithms which has been inspired by an analogy between the marker ordering problem and the famous traveling salesman problem. The approach has been used to build joined maps for the wasp Trichogramma brassicae and on random pooled data sets. In both cases, it compares quite favorably with existing softwares as far as maximum likelihood is considered as a significant criteria.

An interactive framework for RNA secondary structure prediction with a dynamical treatment of constraints.

A novel approach aiding in the prediction of RNA secondary structures is presented. Although phylogenetic methods are the most successful at deriving RNA secondary structures, the are not applicable when the number of sequences or the sequence variability is too low. Methods based on energy minimization are therefore of great interest. However, some of the suboptimal RNA secondary structures computed with classic methods are unsaturated structures, i.e. some structures are included into others. Thus, the incorporation of constraints during the process of folding is not possible, while the incorporation of constraints before the process of folding often introduces a bias into the energy function. This paper describes a new procedure which allows for the incorporation of constraints before and during the process of RNA folding. SAPSSARN is an interactive program which offers a framework, both to specify a secondary structure through a set of folding constraints and to compute all the supoptimal saturated RNA secondary structures which satisfy all the folding constraints. At the start, it relies on the computation of the probabilities of pairing of each base with all others according to McCaskill's algorithm. The constraint satisfaction formulation of the problem deals dynamically with a chosen set of folding constraints and, finally, a search algorithm computes all the suboptimal saturated secondary structures which satisfy those folding constraints. Within such a framework, it is possible to test new ideas about RNA folding and secondary structures, including pseudoknots, can be computed. The program is illustrated with RNA sequences on which we obtained results in agreement with known structures by using a protocol which mimics the hierarchical folding of RNA molecules.

DRAWNA: a program for drawing schematic views of nucleic acids.

A program for drawing automatically exact and schematic views of nucleic acids is described. The program is written in C ANSI and uses the Silicon Graphics GL and Xirisw libraries within the X11/Motif environment. Through menus, the user can choose, specify, and manipulate in real time the three-dimensional views to be displayed. Drawing options include partitioning of structures into differently colored or shaped fragments, representation of backbones as flat or with conic-section ribbons, display of paired or free bases as rods, and display of surfaces as filled or outlined and stereo or depth-cued views.

Automatic display of RNA secondary structures.

A set of programs written in C language with the GL library and under UNIX has been developed for generating compact, pleasant and non-overlapping displays of secondary structures of ribonucleic acids. The first program, rnasearch, implements a new search procedure that dynamically rearranges overlapping portions of the two-dimensional drawing while preserving clear and readable displays of the two-dimensional structure. The algorithm is fast (the execution time for the command rnasearch is 38.6 s for the 16S rRNA of Escherichia coli with 1542 bases), accepts outputs from two-dimensional prediction programs and therefore allows for rapid comparison between the various two-dimensional folds generated. A second program, rnadisplay, allows the graphical display of the computed two-dimensional structures on a graphics workstation. Otherwise, it is possible to obtain a paper output of the two-dimensional structure by using the program print2D which builds a Postscript file. Moreover the two-dimensional drawing can be labelled for representing data coming from chemical modifications and/or enzymatic cleavages. Application to a few secondary structures such as RNaseP, 5S rRNA and 16S rRNA are given.