Concentrations of individual RNA sequences in polyadenylated nuclear and cytoplasmic RNA populations of Drosophila cells.

Steady state concentrations of individual RNA sequences in poly(A) nuclear and cytoplasmic RNA populations of Drosophila Kc cells were determined using cloned cDNA fragments. These cDNAs represent poly(A) RNA sequences of different abundance in the cytoplasm of Kc cells, but their steady state concentrations in poly(A) hnRNA was always lower. Of ten different sequences analysed, eight showed some four-fold lower concentration in hnRNA mRNA, two were underrepresented in hnRNA relative to the others. The obvious clustering of mRNA/hnRNA ratios is discussed in relation to sequence complexity and turnover rates of these RNA populations.     

Corrections to the human mitochondrial ribosomal RNA sequences

Two new errors and one consensus change were identified in the human mitochondrial Cambridge consensus sequence. The errors are an A to G substitution at nucleotide 750 in the 12S rRNA gene and a single nucleotide deletion at nt 3107 in the 16S rRNA gene. The consensus change is nt 2706 AG in the 16S rRNA gene.     

An analytical approach to the problem of inverse optimization with additive objective functions: an application to human prehension

We consider the problem of what is being optimized in human actions with respect to various aspects of human movements and different motor tasks. From the mathematical point of view this problem consists of finding an unknown objective function given the values at which it reaches its minimum. This problem is called the inverse optimization problem. Until now the main approach to this problems has been the cut-and-try method, which consists of introducing an objective function and checking how it reflects the experimental data. Using this approach, different objective functions have been proposed for the same motor action. In the current paper we focus on inverse optimization problems with additive objective functions and linear constraints. Such problems are typical in human movement science. The problem of muscle (or finger) force sharing is an example. For such problems we obtain sufficient conditions for uniqueness and propose a method for determining the objective functions. To illustrate our method we analyze the problem of force sharing among the fingers in a grasping task. We estimate the objective function from the experimental data and show that it can predict the force-sharing pattern for a vast range of external forces and torques applied to the grasped object. The resulting objective function is quadratic with essentially non-zero linear terms.     

TWARIT: an extremely rapid and efficient approach for phylogenetic classification of metagenomic sequences

Phylogenetic assignment of individual sequence reads to their respective taxa, referred to as 'taxonomic binning', constitutes a key step of metagenomic analysis. Existing binning methods have limitations either with respect to time or accuracy/specificity of binning. Given these limitations, development of a method that can bin vast amounts of metagenomic sequence data in a rapid, efficient and computationally inexpensive manner can profoundly influence metagenomic analysis in computational resource poor settings. We introduce TWARIT, a hybrid binning algorithm, that employs a combination of short-read alignment and composition-based signature sorting approaches to achieve rapid binning rates without compromising on binning accuracy and specificity. TWARIT is validated with simulated and real-world metagenomes and the results demonstrate significantly lower overall binning times compared to that of existing methods. Furthermore, the binning accuracy and specificity of TWARIT are observed to be comparable/superior to them. A web server implementing TWARIT algorithm is available at http://metagenomics.atc.tcs.com/Twarit/     

Predicting RNA secondary structure by the comparative approach: how to select the homologous sequences

Background  

The secondary structure of an RNA must be known before the relationship between its structure and function can be determined. One way to predict the secondary structure of an RNA is to identify covarying residues that maintain the pairings (Watson-Crick, Wobble and non-canonical pairings). This "comparative approach" consists of identifying mutations from homologous sequence alignments. The sequences must covary enough for compensatory mutations to be revealed, but comparison is difficult if they are too different. Thus the choice of homologous sequences is critical. While many possible combinations of homologous sequences may be used for prediction, only a few will give good structure predictions. This can be due to poor quality alignment in stems or to the variability of certain sequences. This problem of sequence selection is currently unsolved.     

Restrictions to RNA virus adaptation: an experimental approach

Some basic properties of RNA viruses are their high mutation rate, their enormous population sizes and their short generation time. These properties allow RNA virus populations to quickly explore fitness landscapes. A great adaptability has been amply demonstrated in experimental, as well as in natural, populations of RNA viruses. However, at least from a theoretical point of view, a limit to the extent of viral adaptation may exist as a consequence of adaptive trade-offs arising during evolution in changing environmental conditions. Here, I review previously published results searching for such fitness trade-offs. The following scenario has been explored: the cost of host-range expansion, the cost of resistance to antiviral drugs, and the adaptation to different population densities. Despite the environmental conditions tested, results show a common pattern: whenever a virus adapt to a simple environmental situation it pays a cost in terms of adaptation to alternative situations. However, in those cases where the virus has been simultaneously adapted to different environmental conditions, this cost disappears or, at least, is greatly reduced. Finally, and as another factor imposing a limit to their speed of adaptation, I review results showing that clonal interference also plays an important role during viral evolution. This revised version was published online in August 2006 with corrections to the Cover Date.     

Reverse taxonomy: an approach towards determining the diversity of meiobenthic organisms based on ribosomal RNA signature sequences

Organisms living in or on the sediment layer of water bodies constitute the benthos fauna, which is known to harbour a large number of species of diverse taxonomic groups. The benthos plays a significant role in the nutrient cycle and it is, therefore, of high ecological relevance. Here, we have explored a DNA-taxonomic approach to access the meiobenthic organismic diversity, by focusing on obtaining signature sequences from a part of the large ribosomal subunit rRNA (28S), the D3-D5 region. To obtain a broad representation of taxa, benthos samples were taken from 12 lakes in Germany, representing different ecological conditions. In a first approach, we have extracted whole DNA from these samples, amplified the respective fragment by PCR, cloned the fragments and sequenced individual clones. However, we found a relatively large number of recombinant clones that must be considered PCR artefacts. In a second approach we have, therefore, directly sequenced PCR fragments that were obtained from DNA extracts of randomly picked individual organisms. In total, we have obtained 264 new unique sequences, which can be readily placed into taxon groups, based on phylogenetic comparison with currently available database sequences. The group with the highest taxon abundance were nematodes and protozoa, followed by chironomids. However, we find also that we have by far not exhausted the diversity of organisms in the samples. Still, our data provide a framework within which a meiobenthos DNA signature sequence database can be constructed, that will allow to develop the necessary techniques for studying taxon diversity in the context of ecological analysis. Since many taxa in our analysis are initially only identified via their signature sequences, but not yet their morphology, we propose to call this approach 'reverse taxonomy'.     

Entropy of protein sequences: an integral approach

Several classifications of protein spatial structures and their structural elements are known. This makes revealing of the relation between these structural elements and sequence fragments rather topical. The most important move in this direction would be the determination of positional sensitivity levels and ranges between the residues in protein sequences. In this work the Shannon-Weaver informational entropy was used as a disorder criterion for solving this problem. This entropy was computed as function of the distance between the amino acid residues in different sets of unhomological protein sequences. Similarity of this function for different sets of protein sequences was shown. Analysis of informational entropy allows detecting a long-range positional correlation (> or =30) between the amino acid residues and oscillations with periods of 3.6 and 2.9. These oscillation periods correspond to periodicity of alpha- and 3(10)-helices.     

Dynalign: an algorithm for finding the secondary structure common to two RNA sequences

With the rapid increase in the size of the genome sequence database, computational analysis of RNA will become increasingly important in revealing structure-function relationships and potential drug targets. RNA secondary structure prediction for a single sequence is 73 % accurate on average for a large database of known secondary structures. This level of accuracy provides a good starting point for determining a secondary structure either by comparative sequence analysis or by the interpretation of experimental studies. Dynalign is a new computer algorithm that improves the accuracy of structure prediction by combining free energy minimization and comparative sequence analysis to find a low free energy structure common to two sequences without requiring any sequence identity. It uses a dynamic programming construct suggested by Sankoff. Dynalign, however, restricts the maximum distance, M, allowed between aligned nucleotides in the two sequences. This makes the calculation tractable because the complexity is simplified to O(M(3)N(3)), where N is the length of the shorter sequence.The accuracy of Dynalign was tested with sets of 13 tRNAs, seven 5 S rRNAs, and two R2 3' UTR sequences. On average, Dynalign predicted 86.1 % of known base-pairs in the tRNAs, as compared to 59.7 % for free energy minimization alone. For the 5 S rRNAs, the average accuracy improves from 47.8 % to 86.4 %. The secondary structure of the R2 3' UTR from Drosophila takahashii is poorly predicted by standard free energy minimization. With Dynalign, however, the structure predicted in tandem with the sequence from Drosophila melanogaster nearly matches the structure determined by comparative sequence analysis.     

RNAPack: an integrated NMR approach to RNA structure determination

Over the last decade, a vast number of useful nuclear magnetic resonance (NMR) experiments have been developed and successfully employed to determine the structure and dynamics of RNA oligonucleotides. Despite this progress, high-resolution RNA structure determination by NMR spectroscopy still remains a lengthy process and requires programming and extensive calibrations to perform NMR experiments successfully. To accelerate RNA structure determination by NMR spectroscopy, we have designed and programmed a package of RNA NMR experiments, called RNAPack. The user-friendly package contains a set of semiautomated single, double, and triple resonance NMR experiments, which are fully optimized for high-resolution RNA solution structure determination on Varian NMR spectrometers. RNAPack provides an autocalibration feature that allows rapid calibration of all NMR experiments in a single step and thereby speeds up the NMR data collection and eliminates user errors. In our laboratory, we have successfully employed this technology to solve RNA solution structures of domains of the internal ribosome entry site of the genomic hepatitis C viral RNA in less than 3 months. RNAPack therefore makes NMR spectroscopy an attractive and rapid structural tool and allows integration of atomic resolution structural information into biochemical studies of large RNA systems.     

Evolutionary conservation of the human 7 S RNA sequences

The cloning and characterization of the cytoplasmic 7 S RNAs of HeLa cells has provided pure probes to study the organization of the corresponding genomic DNA sequences. Such analysis has shown that the 7 S L and K RNAs are derived from families of middle repetitive DNA (Ullu & Melli, 1982; Ullu et al., 1982). In this work we analyze the evolutionary conservation of these sequences in the RNA and DNA of distantly related species. Hybridization of the 7 S recombinants to the RNA of rodents, birds, amphibians and echinoderms suggests high conservation of these sequences throughout evolution. Southern blot analysis of genomic DNAs from the same species shows the presence of families of repeated sequences homologous to the 7 S recombinants and Alu DNAs in the genomes of the same species. We were unable to hybridize the 7 S probes to the RNAs of Drosophila melanogaster or Dictyostelium discoideum, although sequence(s) homologous to the 7 S L probe were found in the genome of D. discoideum and to both 7 S L and K probes in the genome of D. melanogaster.     

A universal protocol to generate consensus level genome sequences for foot-and-mouth disease virus and other positive-sense polyadenylated RNA viruses using the Illumina MiSeq

Background

Next-Generation Sequencing (NGS) is revolutionizing molecular epidemiology by providing new approaches to undertake whole genome sequencing (WGS) in diagnostic settings for a variety of human and veterinary pathogens. Previous sequencing protocols have been subject to biases such as those encountered during PCR amplification and cell culture, or are restricted by the need for large quantities of starting material. We describe here a simple and robust methodology for the generation of whole genome sequences on the Illumina MiSeq. This protocol is specific for foot-and-mouth disease virus (FMDV) or other polyadenylated RNA viruses and circumvents both the use of PCR and the requirement for large amounts of initial template.

Results

The protocol was successfully validated using five FMDV positive clinical samples from the 2001 epidemic in the United Kingdom, as well as a panel of representative viruses from all seven serotypes. In addition, this protocol was successfully used to recover 94% of an FMDV genome that had previously been identified as cell culture negative. Genome sequences from three other non-FMDV polyadenylated RNA viruses (EMCV, ERAV, VESV) were also obtained with minor protocol amendments. We calculated that a minimum coverage depth of 22 reads was required to produce an accurate consensus sequence for FMDV O. This was achieved in 5 FMDV/O/UKG isolates and the type O FMDV from the serotype panel with the exception of the 5′ genomic termini and area immediately flanking the poly(C) region.

Conclusions

We have developed a universal WGS method for FMDV and other polyadenylated RNA viruses. This method works successfully from a limited quantity of starting material and eliminates the requirement for genome-specific PCR amplification. This protocol has the potential to generate consensus-level sequences within a routine high-throughput diagnostic environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-828) contains supplementary material, which is available to authorized users.     

Investigating the viability of photo-identification as an objective tool to study endangered sea turtle populations

We assessed the potential of using natural facial markings to identify individuals in an endangered breeding population of loggerhead sea turtles (Caretta caretta). We divided individual turtles into ten groups based on facial (post-ocular) scale patterns to facilitate rapid comparison of new images in a large photographic catalogue of known turtles (exceeding 400 unique individuals). The matching process was validated by using turtles marked with external flipper tags. An experienced observer achieved a mean 99% success in identifying individuals using photo-id. The reliability and wider utility of the technique was assessed through testing the ability of naïve and trained observers to (1) consistently allocate known (i.e. flipper tagged) individuals into the correct groups (2) correctly match known individuals within one group. In all trials the mean success rate in photographic sorting and matching ranged from 68-100%. A 20 minute training session was found to significantly improve observer ability, i.e the photo-id skills were rapidly acquired by inexperienced workers. Photo-id has the benefit of being suitable for male turtles, which do not come ashore to allow conventional tagging, and so are rarely identified. Photo-id may facilitate the assessment of the numbers of male and female turtles at breeding areas and allow adult sex ratios to be measured.     

Immunotyping of human immunodeficiency virus type 1 (HIV): an approach to immunologic classification of HIV

Because immunologic classification of human immunodeficiency virus type 1 (HIV) might be more relevant than genotypic classification for designing polyvalent vaccines, studies were undertaken to determine whether immunologically defined groups of HIV ("immunotypes") could be identified. For these experiments, the V3 region of the 120-kDa envelope glycoprotein (gp120) was chosen for study. Although antibodies (Abs) to V3 may not play a major protective role in preventing HIV infection, identification of a limited number of immunologically defined structures in this extremely variable region would set a precedent supporting the hypothesis that, despite its diversity, the HIV family, like the V3 region, might be divisible into immunotypes. Consequently, the immunochemical reactivities of 1,176 combinations of human anti-V3 monoclonal Abs (MAbs) and V3 peptides, derived from viruses of several clades, were studied. Extensive cross-clade reactivity was observed. The patterns of reactivities of 21 MAbs with 50 peptides from clades A through H were then analyzed by a multivariate statistical technique. To test the validity of the mathematical approach, a cluster analysis of the 21 MAbs was performed. Five groups were identified, and these MAb clusters corresponded to classifications of these same MAbs based on the epitopes which they recognize. The concordance between the MAb clusters identified by mathematical analysis and by their specificities supports the validity of the mathematical approach. Therefore, the same mathematical technique was used to identify clusters within the 50 peptides. Seven groups of peptides, each containing peptides from more than one clade, were defined. Inspection of the amino acid sequences of the peptides in each of the mathematically defined peptide clusters revealed unique "signature sequences" that suggest structural motifs characteristic of each V3-based immunotype. The results suggest that cluster analysis of immunologic data can define immunotypes of HIV. These immunotypes are distinct from genotypic classifications. The methods described pave the way for identification of immunotypes defined by immunochemical and neutralization data generated with anti-HIV Env MAbs and intact, viable HIV virions.     

Organization of sequences related to U6 RNA in the human genome.

Small nuclear RNAs were isolated from human placenta and fractionated into individual molecular species. They were then iodinated with 125I and used as probes to screen the human genome. Of 2 x 10(4) recombinant phage clones screened, 22 clones hybridized with U6 RNA, suggesting that there were about 200 copies of this sequence family per haploid genome. Southern blots of these cloned DNAs digested with several restriction enzymes gave the following results: 1, each clone had only one fragment that carried the U6 sequence, 2, the lengths of these fragments varied from clone to clone. These observations indicate that U6 sequences exist as dispersed middle repetitive DNA, and that the sequences surrounding these loci vary. Two of the loci and their flanking regions were subcloned into plasmid and sequenced. Both of the loci showed microheterogeneity of mainly A/G and T/C, but had closely related sequences to U6 RNAs of rat or mouse. The divergence of the flanking regions begins immediately outside the loci. The implication on the microheterogeneity of the U6-related sequences is discussed.