Microbial stimulation of different Toll-like receptor signalling pathways induces diverse metabolic programmes in human monocytes.

Microbial stimuli such as lipopolysaccharide (LPS) induce robust metabolic rewiring in immune cells known as the Warburg effect. It is unknown whether this increase in glycolysis and decrease in oxidative phosphorylation (OXPHOS) is a general characteristic of monocytes that have encountered a pathogen. Using CD14+ monocytes from healthy donors, we demonstrated that most microbial stimuli increased glycolysis, but that only stimulation of Toll-like receptor (TLR) 4 with LPS led to a decrease in OXPHOS. Instead, activation of other TLRs, such as TLR2 activation by Pam3CysSK4 (P3C), increased oxygen consumption and mitochondrial enzyme activity. Transcriptome and metabolome analysis of monocytes stimulated with P3C versus LPS confirmed the divergent metabolic responses between both stimuli, and revealed significant differences in the tricarboxylic acid cycle, OXPHOS and lipid metabolism pathways following stimulation of monocytes with P3C versus LPS. At a functional level, pharmacological inhibition of complex I of the mitochondrial electron transport chain diminished cytokine production and phagocytosis in P3C- but not LPS-stimulated monocytes. Thus, unlike LPS, complex microbial stimuli and the TLR2 ligand P3C induce a specific pattern of metabolic rewiring that involves upregulation of both glycolysis and OXPHOS, which enables activation of host defence mechanisms such as cytokine production and phagocytosis.

Using bioconductor package BiGGR for metabolic flux estimation based on gene expression changes in brain.

Predicting the distribution of metabolic fluxes in biochemical networks is of major interest in systems biology. Several databases provide metabolic reconstructions for different organisms. Software to analyze flux distributions exists, among others for the proprietary MATLAB environment. Given the large user community for the R computing environment, a simple implementation of flux analysis in R appears desirable and will facilitate easy interaction with computational tools to handle gene expression data. We extended the R software package BiGGR, an implementation of metabolic flux analysis in R. BiGGR makes use of public metabolic reconstruction databases, and contains the BiGG database and the reconstruction of human metabolism Recon2 as Systems Biology Markup Language (SBML) objects. Models can be assembled by querying the databases for pathways, genes or reactions of interest. Fluxes can then be estimated by maximization or minimization of an objective function using linear inverse modeling algorithms. Furthermore, BiGGR provides functionality to quantify the uncertainty in flux estimates by sampling the constrained multidimensional flux space. As a result, ensembles of possible flux configurations are constructed that agree with measured data within precision limits. BiGGR also features automatic visualization of selected parts of metabolic networks using hypergraphs, with hyperedge widths proportional to estimated flux values. BiGGR supports import and export of models encoded in SBML and is therefore interoperable with different modeling and analysis tools. As an application example, we calculated the flux distribution in healthy human brain using a model of central carbon metabolism. We introduce a new algorithm termed Least-squares with equalities and inequalities Flux Balance Analysis (Lsei-FBA) to predict flux changes from gene expression changes, for instance during disease. Our estimates of brain metabolic flux pattern with Lsei-FBA for Alzheimer's disease agree with independent measurements of cerebral metabolism in patients. This second version of BiGGR is available from Bioconductor.

Submembranous recruitment of creatine kinase B supports formation of dynamic actin-based protrusions of macrophages and relies on its C-terminal flexible loop.

Subcellular partitioning of creatine kinase contributes to the formation of patterns in intracellular ATP distribution and the fuelling of cellular processes with a high and sudden energy demand. We have previously shown that brain-type creatine kinase (CK-B) accumulates at the phagocytic cup in macrophages where it is involved in the compartmentalized generation of ATP for actin remodeling. Here, we report that CK-B catalytic activity also helps in the formation of protrusive ruffle structures which are actin-dependent and abundant on the surface of both unstimulated and LPS-activated macrophages. Recruitment of CK-B to these structures occurred transiently and inhibition of the enzyme's catalytic activity with cyclocreatine led to a general smoothening of surface morphology as visualized by scanning electron microscopy. Comparison of the dynamics of distribution of YFP-tagged CK-mutants and isoforms by live imaging revealed that amino acid residues in the C-terminal segment (aa positions 323-330) that forms one of the protein's two mobile loops are involved in partitioning over inner regions of the cytosol and nearby sites where membrane protrusions occur during induction of phagocytic cup formation. Although wt CK-B, muscle-type CK (CK-M), and a catalytically dead CK-B-E232Q mutant with intact loop region were normally recruited from the cytosolic pool, no dynamic transition to the phagocytic cup area was seen for the CK-homologue arginine kinase and a CK-B-D326A mutant protein. Bioinformatics analysis helped us to predict that conformational flexibility of the C-terminal loop, independent of conformational changes induced by substrate binding or catalytic activity, is likely involved in exposing the enzyme for binding at or near the sites of membrane protrusion formation.

NAMPT-mediated salvage synthesis of NAD+ controls morphofunctional changes of macrophages.

Functional morphodynamic behavior of differentiated macrophages is strongly controlled by actin cytoskeleton rearrangements, a process in which also metabolic cofactors ATP and NAD(H) (i.e. NAD+ and NADH) and NADP(H) (i.e. NADP+ and NADPH) play an essential role. Whereas the link to intracellular ATP availability has been studied extensively, much less is known about the relationship between actin cytoskeleton dynamics and intracellular redox state and NAD+-supply. Here, we focus on the role of nicotinamide phosphoribosyltransferase (NAMPT), found in extracellular form as a cytokine and growth factor, and in intracellular form as one of the key enzymes for the production of NAD+ in macrophages. Inhibition of NAD+ salvage synthesis by the NAMPT-specific drug FK866 caused a decrease in cytosolic NAD+ levels in RAW 264.7 and Maf-DKO macrophages and led to significant downregulation of the glycolytic flux without directly affecting cell viability, proliferation, ATP production capacity or mitochondrial respiratory activity. Concomitant with these differential metabolic changes, the capacity for phagocytic ingestion of particles and also substrate adhesion of macrophages were altered. Depletion of cytoplasmic NAD+ induced cell-morphological changes and impaired early adhesion in phagocytosis of zymosan particles as well as spreading performance. Restoration of NAD+ levels by NAD+, NMN, or NADP+ supplementation reversed the inhibitory effects of FK866. We conclude that direct coupling to local, actin-based, cytoskeletal dynamics is an important aspect of NAD+'s cytosolic role in the regulation of morphofunctional characteristics of macrophages.

Glucose controls morphodynamics of LPS-stimulated macrophages.

Macrophages constantly undergo morphological changes when quiescently surveying the tissue milieu for signs of microbial infection or damage, or after activation when they are phagocytosing cellular debris or foreign material. These morphofunctional alterations require active actin cytoskeleton remodeling and metabolic adaptation. Here we analyzed RAW 264.7 and Maf-DKO macrophages as models to study whether there is a specific association between aspects of carbohydrate metabolism and actin-based processes in LPS-stimulated macrophages. We demonstrate that the capacity to undergo LPS-induced cell shape changes and to phagocytose complement-opsonized zymosan (COZ) particles does not depend on oxidative phosphorylation activity but is fueled by glycolysis. Different macrophage activities like spreading, formation of cell protrusions, as well as phagocytosis of COZ, were thereby strongly reliant on the presence of low levels of extracellular glucose. Since global ATP production was not affected by rewiring of glucose catabolism and inhibition of glycolysis by 2-deoxy-D-glucose and glucose deprivation had differential effects, our observations suggest a non-metabolic role for glucose in actin cytoskeletal remodeling in macrophages, e.g. via posttranslational modification of receptors or signaling molecules, or other effects on the machinery that drives actin cytoskeletal changes. Our findings impute a decisive role for the nutrient state of the tissue microenvironment in macrophage morphodynamics.

The human leukocyte antigen-presented ligandome of B lymphocytes.

Peptides presented by human leukocyte antigen (HLA) molecules on the cell surface play a crucial role in adaptive immunology, mediating the communication between T cells and antigen presenting cells. Knowledge of these peptides is of pivotal importance in fundamental studies of T cell action and in cellular immunotherapy and transplantation. In this paper we present the in-depth identification and relative quantification of 14,500 peptide ligands constituting the HLA ligandome of B cells. This large number of identified ligands provides general insight into the presented peptide repertoire and antigen presentation. Our uniquely large set of HLA ligands allowed us to characterize in detail the peptides constituting the ligandome in terms of relative abundance, peptide length distribution, physicochemical properties, binding affinity to the HLA molecule, and presence of post-translational modifications. The presented B-lymphocyte ligandome is shown to be a rich source of information by the presence of minor histocompatibility antigens, virus-derived epitopes, and post-translationally modified HLA ligands, and it can be a good starting point for solving a wealth of specific immunological questions. These HLA ligands can form the basis for reversed immunology approaches to identify T cell epitopes based not on in silico predictions but on the bona fide eluted HLA ligandome.

Intracellular NAD(H) levels control motility and invasion of glioma cells.

Oncogenic transformation involves reprogramming of cell metabolism, whereby steady-state levels of intracellular NAD(+) and NADH can undergo dramatic changes while ATP concentration is generally well maintained. Altered expression of nicotinamide phosphoribosyltransferase (NAMPT), the rate-limiting enzyme of NAD(+)-salvage, accompanies the changes in NAD(H) during tumorigenesis. Here, we show by genetic and pharmacological inhibition of NAMPT in glioma cells that fluctuation in intracellular [NAD(H)] differentially affects cell growth and morphodynamics, with motility/invasion capacity showing the highest sensitivity to [NAD(H)] decrease. Extracellular supplementation of NAD(+) or re-expression of NAMPT abolished the effects. The effects of NAD(H) decrease on cell motility appeared parallel coupled with diminished pyruvate-lactate conversion by lactate dehydrogenase (LDH) and with changes in intracellular and extracellular pH. The addition of lactic acid rescued and knockdown of LDH-A replicated the effects of [NAD(H)] on motility. Combined, our observations demonstrate that [NAD(H)] is an important metabolic component of cancer cell motility. Nutrient or drug-mediated modulation of NAD(H) levels may therefore represent a new option for blocking the invasive behavior of tumors.

ProRepeat: an integrated repository for studying amino acid tandem repeats in proteins.

ProRepeat (http://prorepeat.bioinformatics.nl/) is an integrated curated repository and analysis platform for in-depth research on the biological characteristics of amino acid tandem repeats. ProRepeat collects repeats from all proteins included in the UniProt knowledgebase, together with 85 completely sequenced eukaryotic proteomes contained within the RefSeq collection. It contains non-redundant perfect tandem repeats, approximate tandem repeats and simple, low-complexity sequences, covering the majority of the amino acid tandem repeat patterns found in proteins. The ProRepeat web interface allows querying the repeat database using repeat characteristics like repeat unit and length, number of repetitions of the repeat unit and position of the repeat in the protein. Users can also search for repeats by the characteristics of repeat containing proteins, such as entry ID, protein description, sequence length, gene name and taxon. ProRepeat offers powerful analysis tools for finding biological interesting properties of repeats, such as the strong position bias of leucine repeats in the N-terminus of eukaryotic protein sequences, the differences of repeat abundance among proteomes, the functional classification of repeat containing proteins and GC content constrains of repeats' corresponding codons.

MADMAX - Management and analysis database for multiple ~omics experiments.

The rapid increase of ~omics datasets generated by microarray, mass spectrometry and next generation sequencing technologies requires an integrated platform that can combine results from different ~omics datasets to provide novel insights in the understanding of biological systems. MADMAX is designed to provide a solution for storage and analysis of complex ~omics datasets. In addition, analysis results (such as lists of genes) will be merged to reveal candidate genes supported by all datasets. The system constitutes an ISA-Tab compliant LIMS part which is independent of different analysis pipelines. A pilot study of different type of ~omics data in Brassica rapa demonstrates the possible use of MADMAX. The web-based user interface provides easy access to data and analysis tools on top of the database.

HSPVdb--the Human Short Peptide Variation Database for improved mass spectrometry-based detection of polymorphic HLA-ligands.

T cell epitopes derived from polymorphic proteins or from proteins encoded by alternative reading frames (ARFs) play an important role in (tumor) immunology. Identification of these peptides is successfully performed with mass spectrometry. In a mass spectrometry-based approach, the recorded tandem mass spectra are matched against hypothetical spectra generated from known protein sequence databases. Commonly used protein databases contain a minimal level of redundancy, and thus, are not suitable data sources for searching polymorphic T cell epitopes, either in normal or ARFs. At the same time, however, these databases contain much non-polymorphic sequence information, thereby complicating the matching of recorded and theoretical spectra, and increasing the potential for finding false positives. Therefore, we created a database with peptides from ARFs and peptide variation arising from single nucleotide polymorphisms (SNPs). It is based on the human mRNA sequences from the well-annotated reference sequence (RefSeq) database and associated variation information derived from the Single Nucleotide Polymorphism Database (dbSNP). In this process, we removed all non-polymorphic information. Investigation of the frequency of SNPs in the dbSNP revealed that many SNPs are non-polymorphic "SNPs". Therefore, we removed those from our dedicated database, and this resulted in a comprehensive high quality database, which we coined the Human Short Peptide Variation Database (HSPVdb). The value of our HSPVdb is shown by identification of the majority of published polymorphic SNP- and/or ARF-derived epitopes from a mass spectrometry-based proteomics workflow, and by a large variety of polymorphic peptides identified as potential T cell epitopes in the HLA-ligandome presented by the Epstein-Barr virus cells.

ATP changes the fluorescence lifetime of cyan fluorescent protein via an interaction with His148.

Recently, we described that ATP induces changes in YFP/CFP fluorescence intensities of Fluorescence Resonance Energy Transfer (FRET) sensors based on CFP-YFP. To get insight into this phenomenon, we employed fluorescence lifetime spectroscopy to analyze the influence of ATP on these fluorescent proteins in more detail. Using different donor and acceptor pairs we found that ATP only affected the CFP-YFP based versions. Subsequent analysis of purified monomers of the used proteins showed that ATP has a direct effect on the fluorescence lifetime properties of CFP. Since the fluorescence lifetime analysis of CFP is rather complicated by the existence of different lifetimes, we tested a variant of CFP, i.e. Cerulean, as a monomer and in our FRET constructs. Surprisingly, this CFP variant shows no ATP concentration dependent changes in the fluorescence lifetime. The most important difference between CFP and Cerulean is a histidine residue at position 148. Indeed, changing this histidine in CFP into an aspartic acid results in identical fluorescence properties as observed for the Cerulean fluorescent based FRET sensor. We therefore conclude that the changes in fluorescence lifetime of CFP are affected specifically by possible electrostatic interactions of the negative charge of ATP with the positively charged histidine at position 148. Clearly, further physicochemical characterization is needed to explain the sensitivity of CFP fluorescence properties to changes in environmental (i.e. ATP concentrations) conditions.

Evidence for RNA recombination between distinct isolates of Pepino mosaic virus.

Genetic recombination plays an important role in the evolution of virus genomes. In this study we analyzed publicly available genomic sequences of Pepino mosaic virus (PepMV) for recombination events using several bioinformatics tools. The genome-wide analyses not only confirm the presence of previously found recombination events in PepMV but also provide the first evidence for double recombinant origin of the US2 isolate.

Multi-netclust: an efficient tool for finding connected clusters in multi-parametric networks.

UNLABELLED: Multi-netclust is a simple tool that allows users to extract connected clusters of data represented by different networks given in the form of matrices. The tool uses user-defined threshold values to combine the matrices, and uses a straightforward, memory-efficient graph algorithm to find clusters that are connected in all or in either of the networks. The tool is written in C/C++ and is available either as a form-based or as a command-line-based program running on Linux platforms. The algorithm is fast, processing a network of > 10(6) nodes and 10(8) edges takes only a few minutes on an ordinary computer. AVAILABILITY: http://www.bioinformatics.nl/netclust/.

Gene expression in chicken reveals correlation with structural genomic features and conserved patterns of transcription in the terrestrial vertebrates.

BACKGROUND: The chicken is an important agricultural and avian-model species. A survey of gene expression in a range of different tissues will provide a benchmark for understanding expression levels under normal physiological conditions in birds. With expression data for birds being very scant, this benchmark is of particular interest for comparative expression analysis among various terrestrial vertebrates. METHODOLOGY/PRINCIPAL FINDINGS: We carried out a gene expression survey in eight major chicken tissues using whole genome microarrays. A global picture of gene expression is presented for the eight tissues, and tissue specific as well as common gene expression were identified. A Gene Ontology (GO) term enrichment analysis showed that tissue-specific genes are enriched with GO terms reflecting the physiological functions of the specific tissue, and housekeeping genes are enriched with GO terms related to essential biological functions. Comparisons of structural genomic features between tissue-specific genes and housekeeping genes show that housekeeping genes are more compact. Specifically, coding sequence and particularly introns are shorter than genes that display more variation in expression between tissues, and in addition intergenic space was also shorter. Meanwhile, housekeeping genes are more likely to co-localize with other abundantly or highly expressed genes on the same chromosomal regions. Furthermore, comparisons of gene expression in a panel of five common tissues between birds, mammals and amphibians showed that the expression patterns across tissues are highly similar for orthologous genes compared to random gene pairs within each pair-wise comparison, indicating a high degree of functional conservation in gene expression among terrestrial vertebrates. CONCLUSIONS: The housekeeping genes identified in this study have shorter gene length, shorter coding sequence length, shorter introns, and shorter intergenic regions, there seems to be selection pressure on economy in genes with a wide tissue distribution, i.e. these genes are more compact. A comparative analysis showed that the expression patterns of orthologous genes are conserved in the terrestrial vertebrates during evolution.

A pipeline for high throughput detection and mapping of SNPs from EST databases.

Single nucleotide polymorphisms (SNPs) represent the most abundant type of genetic variation that can be used as molecular markers. The SNPs that are hidden in sequence databases can be unlocked using bioinformatic tools. For efficient application of these SNPs, the sequence set should be error-free as much as possible, targeting single loci and suitable for the SNP scoring platform of choice. We have developed a pipeline to effectively mine SNPs from public EST databases with or without quality information using QualitySNP software, select reliable SNP and prepare the loci for analysis on the Illumina GoldenGate genotyping platform. The applicability of the pipeline was demonstrated using publicly available potato EST data, genotyping individuals from two diploid mapping populations and subsequently mapping the SNP markers (putative genes) in both populations. Over 7000 reliable SNPs were identified that met the criteria for genotyping on the GoldenGate platform. Of the 384 SNPs on the SNP array approximately 12% dropped out. For the two potato mapping populations 165 and 185 SNPs segregating SNP loci could be mapped on the respective genetic maps, illustrating the effectiveness of our pipeline for SNP selection and validation. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-009-9377-5) contains supplementary material, which is available to authorized users.

An evaluation of the performance of three semantic background knowledge sources in comparative anatomy.

In this paper we evaluate the performance and usefulness of three semantic background knowledge sources for predicting synonymous anatomical terms across species boundaries. The reference sources under evaluation are UMLS, FMA-OBO and WordNet, which are applied to the anatomical ontologies of mouse and zebrafish. Our results show that the use of specialized knowledge sources leads to highly accurate predictions, verified through complete manual curation, which can be further improved by combining multiple of said sources. We argue that these three references complement each other in terms of granularity and specificity. From our results we conclude that these references can be used to create reliable ontology mappings with minimal human supervision.

A comparison of database systems for XML-type data.

BACKGROUND: In the field of bioinformatics interchangeable data formats based on XML are widely used. XML-type data is also at the core of most web services. With the increasing amount of data stored in XML comes the need for storing and accessing the data. In this paper we analyse the suitability of different database systems for storing and querying large datasets in general and Medline in particular. RESULTS: All reviewed database systems perform well when tested with small to medium sized datasets, however when the full Medline dataset is queried a large variation in query times is observed. CONCLUSIONS: There is not one system that is vastly superior to the others in this comparison and, depending on the database size and the query requirements, different systems are most suitable. The best all-round solution is the Oracle 11~g database system using the new binary storage option. Alias-i's Lingpipe is a more lightweight, customizable and sufficiently fast solution. It does however require more initial configuration steps. For data with a changing XML structure Sedna and BaseX as native XML database systems or MySQL with an XML-type column are suitable.

A tail-anchored myotonic dystrophy protein kinase isoform induces perinuclear clustering of mitochondria, autophagy, and apoptosis.

BACKGROUND: Studies on the myotonic dystrophy protein kinase (DMPK) gene and gene products have thus far mainly concentrated on the fate of length mutation in the (CTG)n repeat at the DNA level and consequences of repeat expansion at the RNA level in DM1 patients and disease models. Surprisingly little is known about the function of DMPK protein products. METHODOLOGY/PRINCIPAL FINDINGS: We demonstrate here that transient expression of one major protein product of the human gene, the hDMPK A isoform with a long tail anchor, results in mitochondrial fragmentation and clustering in the perinuclear region. Clustering occurred in a variety of cell types and was enhanced by an intact tubulin cytoskeleton. In addition to morphomechanical changes, hDMPK A expression induces physiological changes like loss of mitochondrial membrane potential, increased autophagy activity, and leakage of cytochrome c from the mitochondrial intermembrane space accompanied by apoptosis. Truncation analysis using YFP-hDMPK A fusion constructs revealed that the protein's tail domain was necessary and sufficient to evoke mitochondrial clustering behavior. CONCLUSION/SIGNIFICANCE: Our data suggest that the expression level of the DMPK A isoform needs to be tightly controlled in cells where the hDMPK gene is expressed. We speculate that aberrant splice isoform expression might be a codetermining factor in manifestation of specific DM1 features in patients.