Angela C. M. Luyf

The reach of the genome signature in prokaryotes

BackgroundWith the increased availability of sequenced genomes there have been several initiatives to infer evolutionary relationships by whole genome characteristics. One of these studies suggested good congruence between genome synteny, shared gene content, 16S ribosomal DNA identity, codon usage and the genome signature in prokaryotes. Here we rigorously test the phylogenetic signal of the genome signature, which consists of the genome-specific relative frequencies of dinucleotides, on 334 sequenced prokaryotic genome sequences.ResultsIntrageneric comparisons show that in general the genomic dissimilarity scores are higher than in intraspecific comparisons, in accordance with the suggested phylogenetic signal of the genome signature. Exceptions to this trend, (Bartonella spp., Bordetella spp., Salmonella spp. and Yersinia spp.), which have low average intrageneric genomic dissimilarity scores, suggest that members of these genera might be considered the same species. On the other hand, high genomic dissimilarity values for intraspecific analyses suggest that in some cases (e.g.Prochlorococcus marinus, Pseudomonas fluorescens, Buchnera aphidicola and Rhodopseudomonas palustris) different strains from the same species may actually represent different species. Comparing 16S rDNA identity with genomic dissimilarity values corroborates the previously suggested trend in phylogenetic signal, albeit that the dissimilarity values only provide low resolution.ConclusionThe genome signature has a distinct phylogenetic signal, independent of individual genetic marker genes. A reliable phylogenetic clustering cannot be based on dissimilarity values alone, as bootstrapping is not possible for this parameter. It can however be used to support or refute a given phylogeny and resulting taxonomy.

A Sensitive Assay for Virus Discovery in Respiratory Clinical Samples

In 5–40% of respiratory infections in children, the diagnostics remain negative, suggesting that the patients might be infected with a yet unknown pathogen. Virus discovery cDNA-AFLP (VIDISCA) is a virus discovery method based on recognition of restriction enzyme cleavage sites, ligation of adaptors and subsequent amplification by PCR. However, direct discovery of unknown pathogens in nasopharyngeal swabs is difficult due to the high concentration of ribosomal RNA (rRNA) that acts as competitor. In the current study we optimized VIDISCA by adjusting the reverse transcription enzymes and decreasing rRNA amplification in the reverse transcription, using hexamer oligonucleotides that do not anneal to rRNA. Residual cDNA synthesis on rRNA templates was further reduced with oligonucleotides that anneal to rRNA but can not be extended due to 3′-dideoxy-C6-modification. With these modifications >90% reduction of rRNA amplification was established. Further improvement of the VIDISCA sensitivity was obtained by high throughput sequencing (VIDISCA-454). Eighteen nasopharyngeal swabs were analysed, all containing known respiratory viruses. We could identify the proper virus in the majority of samples tested (11/18). The median load in the VIDISCA-454 positive samples was 7.2 E5 viral genome copies/ml (ranging from 1.4 E3–7.7 E6). Our results show that optimization of VIDISCA and subsequent high-throughput-sequencing enhances sensitivity drastically and provides the opportunity to perform virus discovery directly in patient material.

δρ-Web, an online tool to assess composition similarity of individual nucleic acid sequences

SUMMARY Although whole-genome sequences have been analysed for the presence of anomalous DNA, no dedicated application is currently available to analyse the composition of individual sequence entries, for instance those derived by experimental techniques, such as subtractive hybridization. Since genomic dinucleotide frequency values are conserved between related species, a representative genome sequence can often be found to score for anomalous sequence composition for many of these putative horizontally transferred sequences. We developed the application deltarho-web, which enables the determination of the differences between the dinucleotide composition of an input sequence and that of a selected genome in a size-dependent manner. A feature allowing batch comparisons is included as well. In addition, deltarho-web allows the analysis of the dinucleotide composition of complete genomes. This provides complementary information for the identification of large anomalous gene clusters.

Compositional discordance between prokaryotic plasmids and host chromosomes

BackgroundMost plasmids depend on the host replication machinery and possess partitioning genes. These properties confine plasmids to a limited range of hosts, yielding a close and presumably stable relationship between plasmid and host. Hence, it is anticipated that due to amelioration the dinucleotide composition of plasmids is similar to that of the genome of their hosts. However, plasmids are also thought to play a major role in horizontal gene transfer and thus are frequently exchanged between hosts, suggesting dinucleotide composition dissimilarity between plasmid and host genome. We compared the dinucleotide composition of a large collection of plasmids with that of their host genomes to shed more light on this enigma.ResultsThe dinucleotide frequency, coined the genome signature, facilitates the identification of putative horizontally transferred DNA in complete genome sequences, since it was found to be typical for a certain genome, and similar between related species. By comparison of the genome signature of 230 plasmid sequences with that of the genome of each respective host, we found that in general the genome signature of plasmids is dissimilar from that of their host genome.ConclusionOur results show that the genome signature of plasmids does not resemble that of their host genome. This indicates either absence of amelioration or a less stable relationship between plasmids and their host. We propose an indiscriminate lifestyle for plasmids preserving the genome signature discordance between these episomes and host chromosomes.

Performance of VIDISCA-454 in Feces-Suspensions and Serum

Virus discovery combining sequence unbiased amplification with next generation sequencing is now state-of-the-art. We have previously determined that the performance of the unbiased amplification technique which is operational at our institute, VIDISCA-454, is efficient when respiratory samples are used as input. The performance of the assay is, however, not known for other clinical materials like blood or stool samples. Here, we investigated the sensitivity of VIDISCA-454 with feces-suspensions and serum samples that are positive and that have been quantified for norovirus and human immunodeficiency virus type 1, respectively. The performance of VIDISCA-454 in serum samples was equal to its performance in respiratory material, with an estimated lower threshold of 1,000 viral genome copies. The estimated threshold in feces-suspension is around 200,000 viral genome copies. The decreased sensitivity in feces suspension is mainly due to sequences that share no recognizable identity with known sequences. Most likely these sequences originate from bacteria and phages which are not completely sequenced.

Defining functional classes of Barth syndrome mutation in humans

The X-linked disease Barth syndrome (BTHS) is caused by mutations in TAZ; TAZ is the main determinant of the final acyl chain composition of the mitochondrial-specific phospholipid, cardiolipin. To date, a detailed characterization of endogenous TAZ has only been performed in yeast. Further, why a given BTHS-associated missense mutation impairs TAZ function has only been determined in a yeast model of this human disease. Presently, the detailed characterization of yeast tafazzin harboring individual BTHS mutations at evolutionarily conserved residues has identified seven distinct loss-of-function mechanisms caused by patient-associated missense alleles. However, whether the biochemical consequences associated with individual mutations also occur in the context of human TAZ in a validated mammalian model has not been demonstrated. Here, utilizing newly established monoclonal antibodies capable of detecting endogenous TAZ, we demonstrate that mammalian TAZ, like its yeast counterpart, is localized to the mitochondrion where it adopts an extremely protease-resistant fold, associates non-integrally with intermembrane space-facing membranes and assembles in a range of complexes. Even though multiple isoforms are expressed at the mRNA level, only a single polypeptide that co-migrates with the human isoform lacking exon 5 is expressed in human skin fibroblasts, HEK293 cells, and murine heart and liver mitochondria. Finally, using a new genome-edited mammalian BTHS cell culture model, we demonstrate that the loss-of-function mechanisms for two BTHS alleles that represent two of the seven functional classes of BTHS mutation as originally defined in yeast, are the same when modeled in human TAZ.

Betulinic acid induces a novel cell death pathway that depends on cardiolipin modification

Cancer is associated with strong changes in lipid metabolism. For instance, normal cells take up fatty acids (FAs) from the circulation, while tumour cells generate their own and become dependent on de novo FA synthesis, which could provide a vulnerability to target tumour cells. Betulinic acid (BetA) is a natural compound that selectively kills tumour cells through an ill-defined mechanism that is independent of BAX and BAK, but depends on mitochondrial permeability transition-pore opening. Here we unravel this pathway and show that BetA inhibits the activity of steroyl-CoA-desaturase (SCD-1). This enzyme is overexpressed in tumour cells and critically important for cells that utilize de novo FA synthesis as it converts newly synthesized saturated FAs to unsaturated FAs. Intriguingly, we find that inhibition of SCD-1 by BetA or, alternatively, with a specific SCD-1 inhibitor directly and rapidly impacts on the saturation level of cardiolipin (CL), a mitochondrial lipid that has important structural and metabolic functions and at the same time regulates mitochondria-dependent cell death. As a result of the enhanced CL saturation mitochondria of cancer cells, but not normal cells that do not depend on de novo FA synthesis, undergo ultrastructural changes, release cytochrome c and quickly induce cell death. Importantly, addition of unsaturated FAs circumvented the need for SCD-1 activity and thereby prevented BetA-induced CL saturation and subsequent cytotoxicity, supporting the importance of this novel pathway in the cytotoxicity induced by BetA.

Visualizing metabolic activity on a genome-wide scale

MOTIVATION To enhance the exploration of gene expression data in a metabolic context, one requires an application that allows the integration of this data and which represents this data in a (genome-wide) metabolic map. The layout of this metabolic map must be highly flexible to enable discoveries of biological phenomena. Moreover, it must allow the simultaneous representation of additional information about genes and enzymes. Since the layout and properties of existing maps did not fulfill our requirements, we developed a new way of representing gene expression data in metabolic charts. RESULTS ViMAc generates user-specified (genome-wide) metabolic maps to explore gene expression data. To enhance the interpretation of these maps information such as sub-cellular localization is included. ViMAc can be used to analyse human or yeast expression data obtained with DNA microarrays or SAGE. We introduce our metabolic map method and demonstrate how it can be applied to explore DNA microarray data for yeast. AVAILABILITY ViMAc is freely available for academic institutions on request from the authors.

A sensitive mass spectrometry platform identifies metabolic changes of life history traits in C. elegans

Abnormal nutrient metabolism is a hallmark of aging, and the underlying genetic and nutritional framework is rapidly being uncovered, particularly using C. elegans as a model. However, the direct metabolic consequences of perturbations in life history of C. elegans remain to be clarified. Based on recent advances in the metabolomics field, we optimized and validated a sensitive mass spectrometry (MS) platform for identification of major metabolite classes in worms and applied it to study age and diet related changes. Using this platform that allowed detection of over 600 metabolites in a sample of 2500 worms, we observed marked changes in fatty acids, amino acids and phospholipids during worm life history, which were independent from the germ-line. Worms underwent a striking shift in lipid metabolism after early adulthood that was at least partly controlled by the metabolic regulator AAK-2/AMPK. Most amino acids peaked during development, except aspartic acid and glycine, which accumulated in aged worms. Dietary intervention also influenced worm metabolite profiles and the regulation was highly specific depending on the metabolite class. Altogether, these MS-based methods are powerful tools to perform worm metabolomics for aging and metabolism-oriented studies.

Initial steps towards a production platform for DNA sequence analysis on the grid

BackgroundBioinformatics is confronted with a new data explosion due to the availability of high throughput DNA sequencers. Data storage and analysis becomes a problem on local servers, and therefore it is needed to switch to other IT infrastructures. Grid and workflow technology can help to handle the data more efficiently, as well as facilitate collaborations. However, interfaces to grids are often unfriendly to novice users.ResultsIn this study we reused a platform that was developed in the VL-e project for the analysis of medical images. Data transfer, workflow execution and job monitoring are operated from one graphical interface. We developed workflows for two sequence alignment tools (BLAST and BLAT) as a proof of concept. The analysis time was significantly reduced. All workflows and executables are available for the members of the Dutch Life Science Grid and the VL-e Medical virtual organizations All components are open source and can be transported to other grid infrastructures.ConclusionsThe availability of in-house expertise and tools facilitates the usage of grid resources by new users. Our first results indicate that this is a practical, powerful and scalable solution to address the capacity and collaboration issues raised by the deployment of next generation sequencers. We currently adopt this methodology on a daily basis for DNA sequencing and other applications. More information and source code is available via http://www.bioinformaticslaboratory.nl/