Eivind Coward

SYSTERS, GeneNest, SpliceNest: exploring sequence space from genome to protein

We have integrated the protein families from SYSTERS and the expressed sequence tag (EST) clusters from our database GeneNest with SpliceNest, a new database mapping EST contigs into genomic DNA. The SYSTERS protein sequence cluster set provides an automatically generated classification of all sequences of the SWISS-PROT, TrEMBL and PIR databases into disjoint protein family and superfamily clusters. GeneNest is a database and software package for producing and visualizing gene indices from ESTs and mRNAs. Currently, the database comprises gene indices of human, mouse, Arabidopsis thaliana and zebrafish. SpliceNest is a web-based graphical tool to explore gene structure, including alternative splicing, based on a mapping of the EST consensus sequences from GeneNest to the complete human genome. The integration of SYSTERS, GeneNest and SpliceNest into one framework now permits an overall exploration of the whole sequence space covering protein, mRNA and EST sequences, as well as genomic DNA. The databases are available for querying and browsing at http://cmb.molgen.mpg.de.

Fast sequence clustering using a suffix array algorithm

MOTIVATION Efficient clustering is important for handling the large amount of available EST sequences. Most contemporary methods are based on some kind of all-against-all comparison, resulting in a quadratic time complexity. A different approach is needed to keep up with the rapid growth of EST data. RESULTS A new, fast EST clustering algorithm is presented. Sub-quadratic time complexity is achieved by using an algorithm based on suffix arrays. A prototype implementation has been developed and run on a benchmark data set. The produced clusterings are validated by comparing them to clusterings produced by other methods, and the results are quite promising. AVAILABILITY The source code for the prototype implementation is available under a GPL license from http://www.ii.uib.no/~ketil/bio/.

SpliceNest: visualizing gene structure and alternative splicing based on EST clusters

Abstract Messenger RNA sequences, expressed sequence tags (ESTs), and, in particular, consensus sequences of clustered ESTs provide valuable information about splice variants of genes. We have mapped such (predicted) human mRNA sequences onto human genomic DNA to compute gene structure and splice variants. The results of this computation have been collected in a public database, SpliceNest , with a web-based, interactive graphical user interface. The quality of the presented data profits from previous assembly of ESTs and allows for analysis of gene structure and detection of splice variants in conjunction with the clones and libraries that indicated variants.

Shufflet: shuffling sequences while conserving the k-let counts

Shufflet is a program and a web-application that generates fast random shufflings of sequences (DNA, protein or others), conserving the exact k-let counts for a given k. The sequences are sampled uniformly from all the valid permutations.

Selection in the host structures the microbiota associated with developing cod larvae (Gadus morhua)

Marine fish larvae are immature upon hatching, and share their environment with high numbers of bacteria. The microbial communities associated with developing fish larvae might be structured by other factors than those important in developing terrestrial animals. Here, we analysed the beta (β)-diversity of the microbiota associated with developing cod larvae and compared it with the bacterial communities in water and live feed by applying pyrosequencing of bar coded v4 16S rDNA amplicons. A total of 15 phyla were observed in the cod larval microbiota. Proteobacteria was the most abundant, followed by Firmicutes, Bacteroidetes and Actinobacteria. The composition and diversity of the cod larval microbiota changed considerably with age. The temporal and spatial patterns of β-diversity could not be explained by stochastic processes, and did not coincide with changes in the rearing conditions. Furthermore, the larval microbiota was highly distinct from the water and the live feed microbiota, particularly at early developmental stages. However, the similarity between larval and water microbiota increased with age. This study suggests that strong selection in the host structures the cod larval microbiota. The changes in community structure observed with increasing age can be explained by altered selection pressure due to development of the intestinal system.

A graph based algorithm for generating EST consensus sequences

MOTIVATION EST sequences constitute an abundant, yet error prone resource for computational biology. Expressed sequences are important in gene discovery and identification, and they are also crucial for the discovery and classification of alternative splicing. An important challenge when processing EST sequences is the reconstruction of mRNA by assembling EST clusters into consensus sequences. RESULTS In contrast to the more established assembly tools, we propose an algorithm that constructs a graph over sequence fragments of fixed size, and produces consensus sequences as traversals of this graph. We provide a tool implementing this algorithm, and perform an experiment where the consensus sequences produced by our implementation, as well as by currently available tools, are compared to mRNA. The results show that our proposed algorithm in a majority of the cases produces consensus of higher quality than the established sequence assemblers and at a competitive speed. AVAILABILITY The source code for the implementation is available under a GPL license from http://www.ii.uib.no/~ketil/bioinformatics/ CONTACT ketil@ii.uib.no.

RBR: library-less repeat detection for ESTs

MOTIVATION Repeat sequences in ESTs are a source of problems, in particular for clustering. ESTs are therefore commonly masked against a library of known repeats. High quality repeat libraries are available for the widely studied organisms, but for most other organisms the lack of such libraries is likely to compromise the quality of EST analysis. RESULTS We present a fast, flexible and library-less method for masking repeats in EST sequences, based on match statistics within the EST collection. The method is not linked to a particular clustering algorithm. Extensive testing on datasets using different clustering methods and a genomic mapping as reference shows that this method gives results that are better than or as good as those obtained using RepeatMasker with a repeat library. AVAILABILITY The implementation of RBR is available under the terms of the GPL from http://www.ii.uib.no/~ketil/bioinformatics CONTACT ketil.malde@bccs.uib.no SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Monitoring multiple myeloma by quantification of recurrent mutations in serum

Circulating tumor DNA is a promising biomarker to monitor tumor load and genome alterations. We explored the presence of circulating tumor DNA in multiple myeloma patients and its relation to disease activity during long-term follow-up. We used digital droplet polymerase chain reaction analysis to monitor recurrent mutations, mainly in mitogen activated protein kinase pathway genes NRAS, KRAS and BRAF. Mutations were identified by next-generation sequencing or polymerase chain reaction analysis of bone marrow plasma cells, and their presence analyzed in 251 archived serum samples obtained from 20 patients during a period of up to 7 years. In 17 of 18 patients, mutations identified in bone marrow during active disease were also found in a time-matched serum sample. The concentration of mutated alleles in serum correlated with the fraction in bone marrow plasma cells (r=0.507, n=34, P<0.002). There was a striking covariation between circulating mutation levels and M protein in ten out of 11 patients with sequential samples. When relapse evaluation by circulating tumor DNA and M protein could be directly compared, the circulating tumor DNA showed relapse earlier in two patients (3 and 9 months), later in one patient (4 months) and in three patients there was no difference. In three patients with transformation to aggressive disease, the concentrations of mutations in serum increased up to 400 times, an increase that was not seen for the M protein. In conclusion, circulating tumor DNA in myeloma is a multi-faceted biomarker reflecting mutated cells, total tumor mass and transformation to a more aggressive disease. Its properties are both similar and complementary to M protein.

BRAF V600E mutation in early-stage multiple myeloma: good response to broad acting drugs and no relation to prognosis.

In this study, we analyzed the prevalence and clone size of BRAF V600E mutation in 209 patients with multiple myeloma and related the results to clinical phenotype, response and survival. Biopsies were screened for BRAF V600E by allele-specific real-time PCR (AS-PCR). Positive results were confirmed by immunohistochemistry, Sanger sequencing and, in three patients from whom we had stored purified myeloma cells, whole-exome sequencing. Eleven patients (5.3%) were BRAF V600E mutation positive by AS-PCR and at least one other method. The fraction of mutated cells varied from 4 to 100%. BRAF V600E-positive patients had no characteristic clinical phenotype except for significantly higher levels of serum creatinine (125 versus 86 μmol/l) Seven of eleven patients responded with at least very good partial response to alkylators, immunomodulatory agents or proteasome inhibitors. Progression-free and overall survival were similar in patients with and without the mutation. By this integrated approach, we found that patients with BRAF V600E mutation responded very well to broad acting drugs and there was no relation to prognosis in early-stage myeloma. In particular, a large mutated cell fraction did not correlate with aggressive disease.

Masking repeats while clustering ESTs

A problem in EST clustering is the presence of repeat sequences. To avoid false matches, repeats have to be masked. This can be a time-consuming process, and it depends on available repeat libraries. We present a fast and effective method that aims to eliminate the problems repeats cause in the process of clustering. Unlike traditional methods, repeats are inferred directly from the EST data, we do not rely on any external library of known repeats. This makes the method especially suitable for analysing the ESTs from organisms without good repeat libraries. We demonstrate that the result is very similar to performing standard repeat masking before clustering.