Monitoring the emergence and re-emergence of viral diseases with the goal of containing the spread of viral agents requires both adequate preparedness and quick response. Identifying the causative agent of a new epidemic is one of the most important steps for effective response to disease outbreaks. Traditionally, virus discovery required propagation of the virus in cell culture, a proven technique responsible for the identification of the vast majority of viruses known to date. However, many viruses cannot be easily propagated in cell culture, thus limiting our knowledge of viruses. Viral metagenomic analyses of environmental samples suggest that the field of virology has explored less than 1% of the extant viral diversity. In the last decade, the culture-independent and sequence-independent metagenomic approach has permitted the discovery of many viruses in a wide range of samples. Phylogenetically, some of these viruses are distantly related to previously discovered viruses. In addition, 60-99% of the sequences generated in different viral metagenomic studies are not homologous to known viruses. In this review, we discuss the advances in the area of viral metagenomics during the last decade and their relevance to virus discovery, clinical microbiology and public health. We discuss the potential of metagenomics for characterization of the normal viral population in a healthy community and identification of viruses that could pose a threat to humans through zoonosis. In addition, we propose a new model of the Koch's postulates named the 'Metagenomic Koch's Postulates'. Unlike the original Koch's postulates and the Molecular Koch's postulates as formulated by Falkow, the metagenomic Koch's postulates focus on the identification of metagenomic traits in disease cases. The metagenomic traits that can be traced after healthy individuals have been exposed to the source of the suspected pathogen.

The microbes in the gastrointestinal (GI) tract are of high importance for the health of the host. In this study, Roche 454 pyrosequencing was applied to a pooled set of different 16S rRNA gene amplicons obtained from GI content of common carp (Cyprinus carpio) to make an inventory of the diversity of the microbiota in the GI tract. Compared to other studies, our culture-independent investigation reveals an impressive diversity of the microbial flora of the carp GI tract. The major group of obtained sequences belonged to the phylum Fusobacteria. Bacteroidetes, Planctomycetes and Gammaproteobacteria were other well represented groups of micro-organisms. Verrucomicrobiae, Clostridia and Bacilli (the latter two belonging to the phylum Firmicutes) had fewer representatives among the analyzed sequences. Many of these bacteria might be of high physiological relevance for carp as these groups have been implicated in vitamin production, nitrogen cycling and (cellulose) fermentation.

Vibrio cholerae O1 Amazonia is a pathogen that was isolated from cholera like diarrhea cases in, at least, two countries, Brazil and Ghana. It belongs to a distinct profile by MLSA. The genomic analysis revealed that it contains the Vibrio pathogenicity island-2 and a set of genes related with pathogenesis and fitness, as the type VI secretion system, present in choleragenic V. cholerae strains.

Multiple factors drive the progression from healthy mucosa towards sporadic colorectal carcinomas and accumulating evidence associates intestinal bacteria with disease initiation and progression. Therefore, the aim of this study was to provide a first high-resolution map of colonic dysbiosis that is associated with human colorectal cancer (CRC). To this purpose, the microbiomes colonizing colon tumor tissue and adjacent non-malignant mucosa were compared by deep rRNA sequencing. The results revealed striking differences in microbial colonization patterns between these two sites. Although inter-individual colonization in CRC patients was variable, tumors consistently formed a niche for Coriobacteria and other proposed probiotic bacterial species, while potentially pathogenic Enterobacteria were underrepresented in tumor tissue. As the intestinal microbiota is generally stable during adult life, these findings suggest that CRC-associated physiological and metabolic changes recruit tumor-foraging commensal-like bacteria. These microbes thus have an apparent competitive advantage in the tumor microenvironment and thereby seem to replace pathogenic bacteria that may be implicated in CRC etiology. This first glimpse of the CRC microbiome provides an important step towards full understanding of the dynamic interplay between intestinal microbial ecology and sporadic CRC, which may provide important leads towards novel microbiome-related diagnostic tools and therapeutic interventions.

Sphagnum peatlands are important ecosystems in the methane cycle. Methanotrophs in these ecosystems have been shown to reduce methane emissions and provide additional carbon to Sphagnum mosses. However, little is known about the diversity and identity of the methanotrophs present in and on Sphagnum mosses in these peatlands. In this study, we applied a pmoA microarray and high-throughput 454 pyrosequencing to pmoA PCR products obtained from total DNA from Sphagnum mosses from a Dutch peat bog to investigate the presence of methanotrophs and to compare the two different methods. Both techniques showed comparable results and revealed an abundance of Methylomonas and Methylocystis species in the Sphagnum mosses. The advantage of the microarray analysis is that it is fast and cost-effective, especially when many samples have to be screened. Pyrosequencing is superior in providing pmoA sequences of many unknown or uncultivated methanotrophs present in the Sphagnum mosses and, thus, provided much more detailed and quantitative insight into the microbial diversity.

Background The mouse is a valuable model for unravelling the role of hepcidin in iron homeostasis. Here, we aimed to assess mouse hepcidin-1 (Hep-1) and -2 (Hep-2) peptide levels in serum and urine by a novel mass spectrometry (MS)-based approach. Methods We used time-of-flight (TOF) MS to determine Hep-1 and -2 levels and Fourier transform ion cyclotron resonance (FTICR) and tandem-MS for hepcidin identifications. The method was biologically validated by hepcidin assessment in: i) 3 mouse strains (C57Bl/6; DBA/2 and BABL/c) upon stimulation with intravenous iron and LPS, ii) homozygous Hfe knock out, homozygous transferrin receptor 2 (Y245X) mutated mice and double affected mice, and iii) mice treated with a sublethal hepatotoxic dose of paracetamol. Results Hep-1 detection was restricted to serum, while Hep-2 was only found in urine and consisted of several isoforms. Elevations in serum Hep-1 and urine Hep-2 upon intravenous iron or LPS were only moderate and varied considerably between mouse strains. Serum Hep-1 was decreased in all three hemochromatosis models and lowest in the double affected mouse. Serum Hep-1 levels correlated with liver hepcidin-1 gene expression, while acute liver damage by paracetamol depleted Hep-1 from serum. Furthermore, serum Hep-1 appeared to be an excellent indicator of splenic iron accumulation. Conclusion Hep-1 and Hep-2 peptide responses in experimental mouse agree with the known biology of hepcidin mRNA regulators, and their measurement can now be implemented in iron-related experimental mouse models to provide novel insights in post-transcriptional regulation, hepcidin function, and kinetics.

It is generally accepted that hydrogenosomes (hydrogen-producing organelles) evolved from a mitochondrial ancestor. However, until recently, only indirect evidence for this hypothesis was available. Here we present the almost complete genome of the hydrogen-producing mitochondrion of the anaerobic ciliate Nyctotherus ovalis and show that, except for the notable absence of genes encoding electron-transport chain components of Complexes III, IV and V, it has a gene content similar to the mitochondrial genomes of aerobic ciliates. Analysis of the genome of the hydrogen-producing mitochondrion, in combination with that of more than 9,000 gDNA and cDNA sequences, allows a preliminary reconstruction of the organellar metabolism. The sequence data indicate that N. ovalis possesses hydrogen-producing mitochondria that have a truncated, two step (Complex I and II) electron-transport chain that uses fumarate as electron acceptor. In addition, components of an extensive protein network for the metabolism of amino-acids, defense against oxidative stress, mitochondrial protein synthesis, mitochondrial protein import and processing, and transport of metabolites across the mitochondrial membrane were identified. Genes for MPV17 and ACN9, two hypothetical proteins linked to mitochondrial disease in humans, were also found. The inferred metabolism is remarkably similar to the organellar metabolism of the phylogenetically distant anaerobic Stramenopile Blastocystis. Notably, the Blastocystis organelle and that of the related flagellate Proteromonas lacertae also lacks genes encoding components of Complexes III, IV and V. Thus, our data show that the hydrogenosomes of N. ovalis are highly specialized, hydrogen-producing mitochondria.

In late puberty, estrogen decelerates bone growth by stimulating growth plate maturation. Here, we studied the mechanism of estrogen action using two pubertal growth plate specimens of one girl at Tanner stage B2 and Tanner stage B3. Histological analysis showed that progression of puberty coincided with characteristic morphological changes; a decrease in total growth plate height (p=0.002), height of the individual zones (p<0.001) and an increase in intercolumnar space (p<0.001). Microarray analysis of the specimens identified 394 genes (72% upregulated, 28% downregulated) that changed with the progression of puberty. Overall changes in gene expression were small (average 1.38-fold upregulated and 1.36-fold downregulated genes). The 394 genes mapped to 13 significantly changing pathways (p<0.05) associated with growth plate maturation (e.g., extracellular matrix, cell cycle and cell death). We next scanned the upstream promoter regions of the 394 genes for the presence of evolutionarily conserved binding sites for transcription factors implicated in growth plate maturation such as Estrogen Receptor, Androgen Receptor, Elk1, Stat5b, CREB and RUNX2. High quality motif sites for RUNX2 (87 genes), Elk1 (43 genes) and Stat5b (31 genes), but not estrogen receptor, were evolutionarily conserved, indicating their functional relevance across primates. Moreover, we show that some of these sites are direct target genes of these transcription factors as shown by ChIP assays.

Heterozygous mutations in p63 are associated with split hand/foot malformations (SHFM), orofacial clefting and ectodermal abnormalities. Elucidation of the p63 gene network that includes target genes and regulatory elements may reveal new genes for other malformation disorders. We performed genome-wide DNA-binding profiling by chromatin immunoprecipitation (ChIP) followed by deep sequencing (ChIP-seq) in primary human keratinocytes, and identified potential target genes and regulatory elements controlled by p63. We show that p63 binds to an enhancer element in the SHFM1 locus on chromosome 7q and that this element controls expression of DLX6 and possibly DLX5, both of which are important for limb development. A unique microdeletion including this enhancer element but not the DLX5/DLX6 genes was identified in a patient with SHFM. Our study strongly indicates disruption of a non-coding cis-regulatory element located more than 250 kb from the DLX5/DLX6 genes as a novel disease mechanism in SHFM1. These data provide a proof-of-concept that the catalogue of p63 binding sites identified in this study may be of relevance to the studies of SHFM and other congenital malformations that resemble the p63-associated phenotypes.

This award from the Dutch Scientific Organization (NWO) enables me to do 3 years of independent research. I will interpret the functionality of metagenomes at the level of individual micro-organisms. The award has been highlighted by several news sources, including Gezondheidskrant.nl, Medicalfacts.nl.

Only three biological pathways are known to produce oxygen: photosynthesis, chlorate respiration and the detoxification of reactive oxygen species. Here we present evidence for a fourth pathway, possibly of considerable geochemical and evolutionary importance. The pathway was discovered after metagenomic sequencing of an enrichment culture that couples anaerobic oxidation of methane with the reduction of nitrite to dinitrogen. The complete genome of the dominant bacterium, named 'Candidatus Methylomirabilis oxyfera', was assembled. This apparently anaerobic, denitrifying bacterium encoded, transcribed and expressed the well-established aerobic pathway for methane oxidation, whereas it lacked known genes for dinitrogen production. Subsequent isotopic labelling indicated that 'M. oxyfera' bypassed the denitrification intermediate nitrous oxide by the conversion of two nitric oxide molecules to dinitrogen and oxygen, which was used to oxidize methane. These results extend our understanding of hydrocarbon degradation under anoxic conditions and explain the biochemical mechanism of a poorly understood freshwater methane sink. Because nitrogen oxides were already present on early Earth, our finding opens up the possibility that oxygen was available to microbial metabolism before the evolution of oxygenic photosynthesis.

Background There are thousands of very diverse ciliates species from which only a handful mitochondrial genomes have been studied so far. These genomes are rather similar because the ciliates analysed (Tetrahymena spp. and Paramecium aurelia) are closely related. Here we study the mitochondrial genomes of the hypotrichous ciliates Euplotes minuta and Euplotes crassus. These ciliates are only distantly related to Tetrahymena spp. and Paramecium aurelia, but more closely related to Nyctotherus ovalis, which possesses a hydrogenosomal (mitochondrial) genome. Results The linear mitochondrial genomes of the hypotrichous ciliates Euplotes minuta and Euplotes crassus were sequenced and compared with the mitochondrial genomes of several Tetrahymena species, Paramecium aurelia and the partially sequenced mitochondrial genome of the anaerobic ciliate Nyctotherus ovalis. This study reports new features such as long 5'gene extensions of several mitochondrial genes, extremely long cox1 and cox2 open reading frames and a large repeat in the middle of the linear mitochondrial genome. The repeat separates the open reading frames into two blocks, each having a single direction of transcription, from the repeat towards the ends of the chromosome. Although the Euplotes mitochondrial gene content is almost identical to Paramecium and Tetrahymena, the order of the genes is completely different. In contrast, the 33273 bp (excluding the repeat region) piece of the mitochondrial genome that has been sequenced in both Euplotes species exhibits no difference in gene order. Unexpectedly, many of the mitochondrial genes of E. minuta encoding ribosomal proteins possess N-terminal extensions that are similar to mitochondrial targeting signals. Conclusions The mitochondrial genomes of the hypotrichous ciliates Euplotes minuta and Euplotes crassus are rather different from the previously studied genomes. Many genes are extended in size compared to mitochondrial genes from other sources.

Background Nyctotherus ovalis is a single-celled eukaryote that has hydrogen-producing mitochondria and lives in the hindgut of cockroaches. Like all members of the ciliate taxon, it has two types of nuclei, a micronucleus and a macronucleus. N. ovalis generates its macronuclear chromosomes by forming polytene chromosomes that subsequently develop into macronuclear chromosomes by DNA elimination and rearrangement. Results We examined the structure of these gene-sized macronuclear chromosomes in N. ovalis. We determined the telomeres, subtelomeric regions, UTRs, coding regions and introns by sequencing a large set of macronuclear DNA sequences (4,242) and cDNAs (5,484) and comparing them with each other. The telomeres consist of repeats CCC(AAAACCCC)n, similar to those in spirotrichous ciliates such as Euplotes, Sterkiella (Oxytricha) and Stylonychia. Per sequenced chromosome we found evidence for either a single protein-coding gene, a single tRNA, or the complete ribosomal RNAs cluster. Hence the chromosomes appear to encode single transcripts. In the short subtelomeric regions we identified a few over-represented motifs that could be involved in gene regulation, but there is no consensus polyadenylation site. The introns are short (21-29 nucleotides), and a significant fraction (1/3) of the tiny introns is conserved in the distantly related ciliate Paramecium tetraurelia. As has been observed in P. tetraurelia, the N. ovalis introns tend to contain in-frame stop codons or have a length that is not dividable by three. This pattern causes premature termination of mRNA translation in the event of intron retention, and potentially degradation of unspliced mRNAs by the nonsense-mediated mRNA decay pathway. Conclusions The combination of short leaders, tiny introns and single genes leads to very minimal macronuclear chromosomes. The smallest we identified contained only 150 nucleotides.

The gap between the amount of genome information released by genome sequencing projects and our knowledge about the proteins' functions is rapidly increasing. To fill this gap, various 'genomic-context' methods have been proposed that exploit sequenced genomes to predict the functions of the encoded proteins. One class of methods, phylogenetic profiling, predicts protein function by correlating the phylogenetic distribution of genes with that of other genes or phenotypic characteristics. The functions of a number of proteins, including ones of medical relevance, have thus been predicted and subsequently confirmed experimentally. Additionally, various approaches to measure the similarity of phylogenetic profiles and to account for the phylogenetic bias in the data have been proposed. We review the successful applications of phylogenetic profiling and analyse the performance of various profile similarity measures with a set of one microsporidial and 25 fungal genomes. In the fungi, phylogenetic profiling yields high-confidence predictions for the highest and only the highest scoring gene pairs illustrating both the power and the limitations of the approach. Both practical examples and theoretical considerations suggest that in order to get a reliable and specific picture of a protein's function, results from phylogenetic profiling have to be combined with other sources of evidence.

Background The common marmoset monkey (Callithrix jacchus), a small non-endangered New World primate native to eastern Brazil, is becoming increasingly used as a non-human primate model in biomedical research, drug development and safety assessment. In contrast to the growing interest for the marmoset as an animal model, the molecular tools for genetic analysis are extremely limited. Results Here we report the development of the first marmoset-specific oligonucleotide microarray (EUMAMA) containing probe sets targeting 1541 different marmoset transcripts expressed in hippocampus. These 1541 transcripts represent a wide variety of different functional gene classes. Hybridisation of the marmoset microarray with labelled RNA from hippocampus, cortex and a panel of 7 different peripheral tissues resulted in high detection rates of 85% in the neuronal tissues and on average 70% in the non-neuronal tissues. The expression profiles of the 2 neuronal tissues, hippocampus and cortex, were highly similar, as indicated by a correlation coefficient of 0.96. Several transcripts with a tissue-specific pattern of expression were identified. Besides the marmoset microarray we have generated 3215 ESTs derived from marmoset hippocampus, which have been annotated and submitted to GenBank [GenBank: EF214838 - EF215447, EH380242 - EH382846]. Conclusions We have generated the first marmoset-specific DNA microarray and demonstrated its use to characterise large-scale gene expression profiles of hippocampus but also of other neuronal and non-neuronal tissues. In addition, we have generated a large collection of ESTs of marmoset origin, which are now available in the public domain. These new tools will facilitate molecular genetic research into this non-human primate animal model.

Background The horizontal transfer of expressed genes from Bacteria into Ciliates which live in close contact with each other in the rumen (the foregut of ruminants) was studied using ciliate Expressed Sequence Tags (ESTs). About 4000 ESTs were sequenced from the two main types of rumen Cilates: Entodiniomorphs (Entodinium simplex, Entodinium caudatum, Eudiplodinium maggii, Metadinium medium, Diploplastron affine, Polyplastron multivesiculatum and Epidinium ecaudatum) and Vestibuliferida, previously called Holotrichs (Isotricha prostoma, Isotricha intestinalis and Dasytricha ruminantium). Results A comparison of the sequences with the completely sequenced genomes of Eukaryotes and Prokaryotes, followed by large scale construction and analysis of phylogenies, identified 148 ciliate genes that specifically cluster with genes from the Bacteria. Of these genes, 34 cluster with genes from the Firmicutes, a phylum of Bacteria that is well represented in the rumen. The phylogenetic clustering with bacterial genes, coupled with the absence of close relatives of these genes in the Ciliate Tetrahymena thermophila, indicates that they have recently been acquired via Horizontal Gene Transfer (HGT). Conclusions Among the HGT candidates, we found an over representation (>75%) of genes involved in metabolism, specifically in the catabolism of complex carbohydrates (>45%), a rich food source in the rumen. We propose that the acquisition of these genes has facilitated the Ciliates' colonization of the rumen and provides evidence for the role of HGT in the adaptation to new niches.

Background The massive scale of microarray derived gene expression data allows for a global view of cellular function. Thus far, comparative studies of gene expression between species have been based on the level of expression of the gene across corresponding tissues, or on the co-expression of the gene with another gene. Results To compare gene expression between distant species on a global scale, we introduce the "expression context". The expression context of a gene is based on the co-expression with all other genes that have unambiguous counterparts in both genomes. Employing this new measure, we show 1) that the expression context is largely conserved between orthologs, and 2) that sequence identity shows little correlation with expression context conservation after gene duplication and speciation. Conclusions This means that the degree of sequence identity has a limited predictive quality for differential expression context conservation between orthologs, and thus presumably also for other facets of gene function.

Genome trees are a means to capture the overwhelming amount of phylogenetic information that is present in genomes. Different formalisms have been introduced to reconstruct genome trees on the basis of various aspects of the genome. On the basis of these aspects, we separate genome trees into five classes: (a) alignment-free trees based on statistic properties of the genome, (b) gene content trees based on the presence and absence of genes, (c) trees based on chromosomal gene order, (d) trees based on average sequence similarity, and (e) phylogenomics-based genome trees. Despite their recent development, genome tree methods have already had some impact on the phylogenetic classification of bacterial species. However, their main impact so far has been on our understanding of the nature of genome evolution and the role of horizontal gene transfer therein. An ideal genome tree method should be capable of using all gene families, including those containing paralogs, in a phylogenomics framework capitalizing on existing methods in conventional phylogenetic reconstruction. We expect such sophisticated methods to help us resolve the branching order between the main bacterial phyla.

In this literature thesis, I attempt to give an update on the state of the art in computerised network reconstruction techniques, and explicitly relate this to actual biological gene networks. I will go into the model formalisms used to describe genetic networks, and explain their specific advantages and disadvantages. Also, a separate chapter will be dedicated to several experimental results obtained in the research of genetic networks, and finally, a short discussion and some hypothesising is added.

[1] S. Gupta, N. Ferguson & R. Anderson: "Chaos, Persistance, and Evolution of Strain Structure in Antigenically Diverse Infectious Agents", Science, 280: 912-915, 1998.