Jostein Johansen

Assessment of composite motif discovery methods

BackgroundComputational discovery of regulatory elements is an important area of bioinformatics research and more than a hundred motif discovery methods have been published. Traditionally, most of these methods have addressed the problem of single motif discovery – discovering binding motifs for individual transcription factors. In higher organisms, however, transcription factors usually act in combination with nearby bound factors to induce specific regulatory behaviours. Hence, recent focus has shifted from single motifs to the discovery of sets of motifs bound by multiple cooperating transcription factors, so called composite motifs or cis-regulatory modules. Given the large number and diversity of methods available, independent assessment of methods becomes important. Although there have been several benchmark studies of single motif discovery, no similar studies have previously been conducted concerning composite motif discovery.ResultsWe have developed a benchmarking framework for composite motif discovery and used it to evaluate the performance of eight published module discovery tools. Benchmark datasets were constructed based on real genomic sequences containing experimentally verified regulatory modules, and the module discovery programs were asked to predict both the locations of these modules and to specify the single motifs involved. To aid the programs in their search, we provided position weight matrices corresponding to the binding motifs of the transcription factors involved. In addition, selections of decoy matrices were mixed with the genuine matrices on one dataset to test the response of programs to varying levels of noise.ConclusionAlthough some of the methods tested tended to score somewhat better than others overall, there were still large variations between individual datasets and no single method performed consistently better than the rest in all situations. The variation in performance on individual datasets also shows that the new benchmark datasets represents a suitable variety of challenges to most methods for module discovery.

Identification of serum microRNA profiles in colon cancer.

Background:microRNAs (miRNAs) exist in blood in an apparently stable form. We have explored whether serum miRNAs can be used as non-invasive early biomarkers of colon cancer.Methods:Serum samples from 30 patients with colon cancer stage IV and 10 healthy controls were examined for the expression of 375 cancer-relevant miRNAs. Based on the miRNA profile in this study, 34 selected miRNAs were measured in serum from 40 patients with stage I–II colon cancer and from 10 additional controls.Results:Twenty miRNAs were differentially expressed in serum from stage IV patients compared with controls (P<0.01). Unsupervised clustering revealed four subgroups; one corresponding mostly to the control group and the three others to the patient groups. Of the 34 miRNAs measured in the follow-up study of stage I–II patients, 21 showed concordant expression between stage IV and stage I–II patient. Based on the profiles of these 21 miRNAs, a supervised linear regression analysis (Partial Least Squares Regression) was performed. Using this model we correctly assigned stage I–II colon cancer patients based on miRNA profiles of stage IV patients.Conclusion:Serum miRNA expression profiling may be utilised in early detection of colon cancer.

Gene duplications in prokaryotes can be associated with environmental adaptation

BackgroundGene duplication is a normal evolutionary process. If there is no selective advantage in keeping the duplicated gene, it is usually reduced to a pseudogene and disappears from the genome. However, some paralogs are retained. These gene products are likely to be beneficial to the organism, e.g. in adaptation to new environmental conditions. The aim of our analysis is to investigate the properties of paralog-forming genes in prokaryotes, and to analyse the role of these retained paralogs by relating gene properties to life style of the corresponding prokaryotes.ResultsParalogs were identified in a number of prokaryotes, and these paralogs were compared to singletons of persistent orthologs based on functional classification. This showed that the paralogs were associated with for example energy production, cell motility, ion transport, and defence mechanisms. A statistical overrepresentation analysis of gene and protein annotations was based on paralogs of the 200 prokaryotes with the highest fraction of paralog-forming genes. Biclustering of overrepresented gene ontology terms versus species was used to identify clusters of properties associated with clusters of species. The clusters were classified using similarity scores on properties and species to identify interesting clusters, and a subset of clusters were analysed by comparison to literature data. This analysis showed that paralogs often are associated with properties that are important for survival and proliferation of the specific organisms. This includes processes like ion transport, locomotion, chemotaxis and photosynthesis. However, the analysis also showed that the gene ontology terms sometimes were too general, imprecise or even misleading for automatic analysis.ConclusionsProperties described by gene ontology terms identified in the overrepresentation analysis are often consistent with individual prokaryote lifestyles and are likely to give a competitive advantage to the organism. Paralogs and singletons dominate different categories of functional classification, where paralogs in particular seem to be associated with processes involving interaction with the environment.

High coverage sequencing of DNA from microorganisms living in an oil reservoir 2.5 kilometres subsurface

Microorganisms colonize a variety of extreme environments, and based on cultivation studies and analyses of PCR-amplified 16S rDNA sequences, microbial life appears to extend deep into the earth crust. However, none of these studies involved comprehensive characterizations of total DNA. Here we report results of a high-coverage DNA pyrosequencing of an apparently representative and uncontaminated sample from a deep sea oil reservoir located 2.5 km subsurface, attributing a pressure and temperature of 250 bars and 85°C respectively. Bioinformatic analyses of the DNA sequences indicate that the reservoir harbours a rich microbial community dominated by a smaller number of taxa. Comparison of the metagenome with sequences in databases indicated that there may have been contact between the oil reservoir and surface communities late in the sequence of geological events leading to oil reservoir formation. One specific gene, encoding a putative enolase, was synthesized and expressed in Escherichia coli. Enolase activity was confirmed and was found to be much more thermotolerant than for a corresponding E. coli enzyme, consistent with the conditions in the oil reservoir.

A novel POLE mutation associated with cancers of colon, pancreas, ovaries and small intestine.

In some families there is an increased risk for colorectal cancer, caused by heritable, but often unidentified genetic mutations predisposing to the disease. We have identified the likely genetic cause for disease predisposition in a large family with high burden of colorectal adenomas and carcinomas, in addition to extra-colonic cancers. This family had previously been tested for known cancer susceptibility genes, with negative results. Exome sequencing was used to identify a novel mutation, c.1373A>T (p.Tyr458Phe), in the gene for DNA polymerase epsilon catalytic subunit (POLE). This mutation is located in the active site of the exonuclease domain of the enzyme, and affects a residue that has previously been shown to be important for exonuclease activity. The first predisposing mutation identified in POLE (c.1270C>G, p.Leu424Val) was associated with colorectal cancer only, but another mutation with a broader tumour spectrum (c.1089C>A, p.Asn363Lys) has recently been reported. In the family described in the present study, carriers generally have multiple colorectal adenomas and cancer of colon, pancreas, ovaries and small intestine which represents an important broadening of the tumour spectrum of POLE mutation carriers. We also observe a large phenotypic variation among the POLE mutation carriers in this family, most likely explained by modifying variants in other genes. One POLE mutation carrier has a novel variant in EXO1 (c.458C>T, p.Ala153Val), which may contribute to a more severe phenotype. The findings in this study will have important implications for risk assessment and surveillance of POLE mutation carriers.

The microbial communities in two apparently physically separated deep subsurface oil reservoirs show extensive DNA sequence similarities.

It is well established that micro-organisms colonize a variety of extreme environments, including habitats like oil reservoirs deep inside the earth crust. Here, we present the results of a comparative high-coverage DNA sequencing study of metagenomes derived from two different oil reservoirs, both located about 2.5 km subseafloor below the Norwegian Sea. A previously reported bioinformatic analysis of DNA sequence data derived from one of the reservoirs (Well I) indicated that the community is dominated by bacterial species with a smaller fraction of Archaea. Here, we report results of a similar analysis from another reservoir (Well II) located in the same geographical area, however, according to available geological knowledge lacking direct physical contact with Well I. Interestingly, the Well II community is largely dominated by Archaea with a subordinate fraction of Bacteria. Comparison of the two datasets showed that large fractions of the sequences are extremely similar, both with respect to identity (typically above 98%) and gene organization. We therefore conclude that both wells contain essentially the same organisms, but in different relative abundances. Assuming that the communities have been distinct for long timescales because of physical separation, the results also indicate that microbial growth in the reservoirs is extremely slow.

Relationship between operon preference and functional properties of persistent genes in bacterial genomes

BackgroundGenes in bacteria may be organised into operons, leading to strict co-expression of the genes that participate in the same operon. However, comparisons between different bacterial genomes have shown that much of the operon structure is dynamic on an evolutionary time scale. This indicates that there are opposing effects influencing the tendency for operon formation, and these effects may be reflected in properties like evolutionary rate, complex formation, metabolic pathways and gene fusion.ResultsWe have used multi-species protein-protein comparisons to generate a high-quality set of genes that are persistent in bacterial genomes (i.e. they have close to universal distribution). We have analysed these genes with respect to operon participation and important functional properties, including evolutionary rate and protein-protein interactions.ConclusionsGenes for ribosomal proteins show a very slow rate of evolution. This is consistent with a strong tendency for the genes to participate in operons and for their proteins to be involved in essential and well defined complexes. Persistent genes for non-ribosomal proteins can be separated into two classes according to tendency to participate in operons. Those with a strong tendency for operon participation make proteins with fewer interaction partners that seem to participate in relatively static complexes and possibly linear pathways. Genes with a weak tendency for operon participation tend to produce proteins with more interaction partners, but possibly in more dynamic complexes and convergent pathways. Genes that are not regulated through operons are therefore more evolutionary constrained than the corresponding operon-associated genes and will on average evolve more slowly.

Novel Antimicrobial Peptides EeCentrocins 1, 2 and EeStrongylocin 2 from the Edible Sea Urchin Echinus esculentus Have 6-Br-Trp Post-Translational Modifications

The global problem of microbial resistance to antibiotics has resulted in an urgent need to develop new antimicrobial agents. Natural antimicrobial peptides are considered promising candidates for drug development. Echinoderms, which rely on innate immunity factors in the defence against harmful microorganisms, are sources of novel antimicrobial peptides. This study aimed to isolate and characterise antimicrobial peptides from the Edible sea urchin Echinus esculentus. Using bioassay-guided purification and cDNA cloning, three antimicrobial peptides were characterised from the haemocytes of the sea urchin; two heterodimeric peptides and a cysteine-rich peptide. The peptides were named EeCentrocin 1 and 2 and EeStrongylocin 2, respectively, due to their apparent homology to the published centrocins and strongylocins isolated from the green sea urchin Strongylocentrotus droebachiensis. The two centrocin-like peptides EeCentrocin 1 and 2 are intramolecularly connected via a disulphide bond to form a heterodimeric structure, containing a cationic heavy chain of 30 and 32 amino acids and a light chain of 13 amino acids. Additionally, the light chain of EeCentrocin 2 seems to be N-terminally blocked by a pyroglutamic acid residue. The heavy chains of EeCentrocins 1 and 2 were synthesised and shown to be responsible for the antimicrobial activity of the natural peptides. EeStrongylocin 2 contains 6 cysteines engaged in 3 disulphide bonds. A fourth peptide (Ee4635) was also discovered but not fully characterised. Using mass spectrometric and NMR analyses, EeCentrocins 1 and 2, EeStrongylocin 2 and Ee4635 were all shown to contain post-translationally brominated Trp residues in the 6 position of the indole ring.

Comparative genomics to delineate pathogenic potential in Non-O157 Shiga toxin-producing Escherichia coli (STEC) from patients with and without haemolytic uremic syndrome (HUS) in Norway

Shiga toxin-producing Escherichia coli (STEC) cause infections in humans ranging from asymptomatic carriage to bloody diarrhoea and haemolytic uremic syndrome (HUS). Here we present whole genome comparison of Norwegian non-O157 STEC strains with the aim to distinguish between strains with the potential to cause HUS and less virulent strains. Whole genome sequencing and comparisons were performed across 95 non-O157 STEC strains. Twenty-three of these were classified as HUS-associated, including strains from patients with HUS (n = 19) and persons with an epidemiological link to a HUS-case (n = 4). Genomic comparison revealed considerable heterogeneity in gene content across the 95 STEC strains. A clear difference in gene profile was observed between strains with and without the Locus of Enterocyte Effacement (LEE) pathogenicity island. Phylogenetic analysis of the core genome showed high degree of diversity among the STEC strains, but all HUS-associated STEC strains were distributed in two distinct clusters within phylogroup B1. However, non-HUS strains were also found in these clusters. A number of accessory genes were found to be significantly overrepresented among HUS-associated STEC, but none of them were unique to this group of strains, suggesting that different sets of genes may contribute to the pathogenic potential in different phylogenetic STEC lineages. In this study we were not able to clearly distinguish between HUS-associated and non-HUS non-O157 STEC by extensive genome comparisons. Our results indicate that STECs from different phylogenetic backgrounds have independently acquired virulence genes that determine pathogenic potential, and that the content of such genes is overlapping between HUS-associated and non-HUS strains.

Use of multigene-panel identifies pathogenic variants in several CRC-predisposing genes in patients previously tested for Lynch Syndrome

Many families with a high burden of colorectal cancer fulfil the clinical criteria for Lynch Syndrome. However, in about half of these families, no germline mutation in the mismatch repair genes known to be associated with this disease can be identified. The aim of this study was to find the genetic cause for the increased colorectal cancer risk in these unsolved cases.