Inge Jonassen

The genome sequence of Atlantic cod reveals a unique immune system

Atlantic cod (Gadus morhua) is a large, cold-adapted teleost that sustains long-standing commercial fisheries and incipient aquaculture. Here we present the genome sequence of Atlantic cod, showing evidence for complex thermal adaptations in its haemoglobin gene cluster and an unusual immune architecture compared to other sequenced vertebrates. The genome assembly was obtained exclusively by 454 sequencing of shotgun and paired-end libraries, and automated annotation identified 22,154 genes. The major histocompatibility complex (MHC) II is a conserved feature of the adaptive immune system of jawed vertebrates, but we show that Atlantic cod has lost the genes for MHC II, CD4 and invariant chain (Ii) that are essential for the function of this pathway. Nevertheless, Atlantic cod is not exceptionally susceptible to disease under natural conditions. We find a highly expanded number of MHC I genes and a unique composition of its Toll-like receptor (TLR) families. This indicates how the Atlantic cod immune system has evolved compensatory mechanisms in both adaptive and innate immunity in the absence of MHC II. These observations affect fundamental assumptions about the evolution of the adaptive immune system and its components in vertebrates.

Approaches to the automatic discovery of patterns in biosequences.

This paper surveys approaches to the discovery of patterns in biosequences and places these approaches within a formal framework that systematises the types of patterns and the discovery algorithms. Patterns with expressive power in the class of regular languages are considered, and a classification of pattern languages in this class is developed, covering the patterns that are the most frequently used in molecular bioinformatics. A formulation is given of the problem of the automatic discovery of such patterns from a set of sequences, and an analysis is presented of the ways in which an assessment can be made of the significance of the discovered patterns. It is shown that the problem is related to problems studied in the field of machine learning. The major part of this paper comprises a review of a number of existing methods developed to solve the problem and how these relate to each other, focusing on the algorithms underlying the approaches. A comparison is given of the algorithms, and examples are given of patterns that have been discovered using the different methods.

J-Express: exploring gene expression data using Java

J-Express is a Java application that allows the user to analyze gene expression (microarray) data in a flexible way giving access to multidimensional scaling, clustering, and visualization methods in an integrated manner. Specifically, J-Express includes implementations of hierarchical clustering, k-means, principal component analysis, and self-organizing maps. At present, it does not include methods for comparing two or more experiments for differentially expressed genes. The application is completely portable and requires only that a Java runtime environment 1.2 is installed on the system. Its efficiency allows interactive clustering of thousands of expression profiles on standard personal computers.

New feature subset selection procedures for classification of expression profiles

BackgroundMethods for extracting useful information from the datasets produced by microarray experiments are at present of much interest. Here we present new methods for finding gene sets that are well suited for distinguishing experiment classes, such as healthy versus diseased tissues. Our methods are based on evaluating genes in pairs and evaluating how well a pair in combination distinguishes two experiment classes. We tested the ability of our pair-based methods to select gene sets that generalize the differences between experiment classes and compared the performance relative to two standard methods. To assess the ability to generalize class differences, we studied how well the gene sets we select are suited for learning a classifier.ResultsWe show that the gene sets selected by our methods outperform the standard methods, in some cases by a large margin, in terms of cross-validation prediction accuracy of the learned classifier. We show that on two public datasets, accurate diagnoses can be made using only 15-30 genes. Our results have implications for how to select marker genes and how many gene measurements are needed for diagnostic purposes.ConclusionWhen looking for differential expression between experiment classes, it may not be sufficient to look at each gene in a separate universe. Evaluating combinations of genes reveals interesting information that will not be discovered otherwise. Our results show that class prediction can be improved by taking advantage of this extra information.

The Atlantic salmon genome provides insights into rediploidization

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

Integrated genomic profiling of endometrial carcinoma associates aggressive tumors with indicators of PI3 kinase activation

Although 75% of endometrial cancers are treated at an early stage, 15% to 20% of these recur. We performed an integrated analysis of genome-wide expression and copy-number data for primary endometrial carcinomas with extensive clinical and histopathological data to detect features predictive of recurrent disease. Unsupervised analysis of the expression data distinguished 2 major clusters with strikingly different phenotypes, including significant differences in disease-free survival. To identify possible mechanisms for these differences, we performed a global genomic survey of amplifications, deletions, and loss of heterozygosity, which identified 11 significantly amplified and 13 significantly deleted regions. Amplifications of 3q26.32 harboring the oncogene PIK3CA were associated with poor prognosis and segregated with the aggressive transcriptional cluster. Moreover, samples with PIK3CA amplification carried signatures associated with in vitro activation of PI3 kinase (PI3K), a signature that was shared by aggressive tumors without PIK3CA amplification. Tumors with loss of PTEN expression or PIK3CA overexpression that did not have PIK3CA amplification also shared the PI3K activation signature, high protein expression of the PI3K pathway member STMN1, and an aggressive phenotype in test and validation datasets. However, mutations of PTEN or PIK3CA were not associated with the same expression profile or aggressive phenotype. STMN1 expression had independent prognostic value. The results affirm the utility of systematic characterization of the cancer genome in clinically annotated specimens and suggest the particular importance of the PI3K pathway in patients who have aggressive endometrial cancer.

Sequencing the genome of the Atlantic salmon (Salmo salar)

The International Collaboration to Sequence the Atlantic Salmon Genome (ICSASG) will produce a genome sequence that identifies and physically maps all genes in the Atlantic salmon genome and acts as a reference sequence for other salmonids.

Structure Comparison and Structure Patterns

This article investigates aspects of pairwise and multiple structure comparison, and the problem of automatically discover common patterns in a set of structures. Descriptions and representation of structures and patterns are described, as well as scoring and algorithms for comparison and discovery. A framework and nomenclature is developed for classifying different methods, and many of these are reviewed and placed into this framework.

Angiogenesis-independent tumor growth mediated by stem-like cancer cells

In this work, highly infiltrative brain tumors with a stem-like phenotype were established by xenotransplantation of human brain tumors in immunodeficient nude rats. These tumors coopted the host vasculature and presented as an aggressive disease without signs of angiogenesis. The malignant cells expressed neural stem cell markers, showed a migratory behavior similar to normal human neural stem cells, and gave rise to tumors in vivo after regrafting. Serial passages in animals gradually transformed the tumors into an angiogenesis-dependent phenotype. This process was characterized by a reduction in stem cells markers. Gene expression profiling combined with high throughput immunoblotting analyses of the angiogenic and nonangiogenic tumors identified distinct signaling networks in the two phenotypes. Furthermore, proinvasive genes were up-regulated and angiogenesis signaling genes were down-regulated in the stem-like tumors. In contrast, proinvasive genes were down-regulated in the angiogenesis-dependent tumors derived from the stem-like tumors. The described angiogenesis-independent tumor growth and the uncoupling of invasion and angiogenesis, represented by the stem-like cancer cells and the cells derived from them, respectively, point at two completely independent mechanisms that drive tumor progression. This article underlines the need for developing therapies that specifically target the stem-like cell pools in tumors.

Characteristics of 454 pyrosequencing data—enabling realistic simulation with flowsim

Motivation: The commercial launch of 454 pyrosequencing in 2005 was a milestone in genome sequencing in terms of performance and cost. Throughout the three available releases, average read lengths have increased to ∼500 base pairs and are thus approaching read lengths obtained from traditional Sanger sequencing. Study design of sequencing projects would benefit from being able to simulate experiments. Results: We explore 454 raw data to investigate its characteristics and derive empirical distributions for the flow values generated by pyrosequencing. Based on our findings, we implement Flowsim, a simulator that generates realistic pyrosequencing data files of arbitrary size from a given set of input DNA sequences. We finally use our simulator to examine the impact of sequence lengths on the results of concrete whole-genome assemblies, and we suggest its use in planning of sequencing projects, benchmarking of assembly methods and other fields. Availability: Flowsim is freely available under the General Public License from http://blog.malde.org/index.php/flowsim/ Contact: susanne.balzer@imr.no; ketil.malde@imr.no