BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics

Phylogenetic inference is fundamental to our understanding of most aspects of the origin and evolution of life, and in recent years, there has been a concentration of interest in statistical approaches such as Bayesian inference and maximum likelihood estimation. Yet, for large data sets and realistic or interesting models of evolution, these approaches remain computationally demanding. High-throughput sequencing can yield data for thousands of taxa, but scaling to such problems using serial computing often necessitates the use of nonstatistical or approximate approaches. The recent emergence of graphics processing units (GPUs) provides an opportunity to leverage their excellent floating-point computational performance to accelerate statistical phylogenetic inference. A specialized library for phylogenetic calculation would allow existing software packages to make more effective use of available computer hardware, including GPUs. Adoption of a common library would also make it easier for other emerging computing architectures, such as field programmable gate arrays, to be used in the future. We present BEAGLE, an application programming interface (API) and library for high-performance statistical phylogenetic inference. The API provides a uniform interface for performing phylogenetic likelihood calculations on a variety of compute hardware platforms. The library includes a set of efficient implementations and can currently exploit hardware including GPUs using NVIDIA CUDA, central processing units (CPUs) with Streaming SIMD Extensions and related processor supplementary instruction sets, and multicore CPUs via OpenMP. To demonstrate the advantages of a common API, we have incorporated the library into several popular phylogenetic software packages. The BEAGLE library is free open source software licensed under the Lesser GPL and available from http://beagle-lib.googlecode.com. An example client program is available as public domain software.     

Revealing Different Roles of the mTOR-Targets S6K1 and S6K2 in Breast Cancer by Expression Profiling and Structural Analysis

Background

The AKT/mTORC1/S6K pathway is frequently overstimulated in breast cancer, constituting a promising therapeutic target. The benefit from mTOR inhibitors varies, likely as a consequence of tumour heterogeneity, and upregulation of several compensatory feed-back mechanisms. The mTORC1 downstream effectors S6K1, S6K2, and 4EBP1 are amplified and overexpressed in breast cancer, associated with a poor outcome and divergent endocrine treatment benefit. S6K1 and S6K2 share high sequence homology, but evidence of partly distinct biological functions is emerging. The aim of this work was to explore possible different roles and treatment target potentials of S6K1 and S6K2 in breast cancer.

Materials and methods

Whole-genome expression profiles were compared for breast tumours expressing high levels of S6K1, S6K2 or 4EBP1, using public datasets, as well as after in vitro siRNA downregulation of S6K1 and/or S6K2 in ZR751 breast cancer cells. In silico homology modelling of the S6K2 kinase domain was used to evaluate its possible structural divergences to S6K1.

Results

Genome expression profiles were highly different in S6K1 and S6K2 high tumours, whereas S6K2 and 4EBP1 profiles showed significant overlaps, both correlated to genes involved in cell cycle progression, among these the master regulator E2F1. S6K2 and 4EBP1 were inversely associated with IGF1 levels, and their prognostic value was shown to be restricted to tumours positive for IGFR and/or HER2. In vitro, S6K1 and S6K2 silencing resulted in upregulation of genes in the mTORC1 and mTORC2 complexes. Isoform-specific silencing also showed distinct patterns, e.g. S6K2 downregulation lead to upregulation of several cell cycle associated genes. Structural analyses of the S6K2 kinase domain showed unique structure patterns, deviating from those of S6K1, facilitating the development of isoform-specific inhibitors. Our data support emerging proposals of distinct biological features of S6K1 and S6K2, suggesting their importance as separate oncogenes and clinical markers, where specific targeting in different breast cancer subtypes could facilitate further individualised therapies.     

Substantial Metabolic Activity of Human Brown Adipose Tissue during Warm Conditions and Cold-Induced Lipolysis of Local Triglycerides

1. Download high-res image (141KB)
2. Download full-size image

     

Decreased oxidative stress during glycolytic inhibition enables maintenance of ATP production and astrocytic survival

Depressed energy metabolism and oxidative stress are common features in many pathological situations in the brain, including stroke. In order to investigate astrocytic responses to such stress, we induced metabolic depression in cultured rat astrocytes. Iodoacetate (IA), an inhibitor of the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used and resulted in a rapid inhibition of GAPDH activity. After 1h of GAPDH inhibition the ATP levels started to decrease and were completely abolished at 4h. In parallel, the activity of reactive oxygen species (ROS) was significantly increased, followed by extensive cell death involving flipping of phosphatidylserine and translocation of apoptosis-inducing factor, but not caspase-3 activation. When IA was combined with azide, a respiratory chain complex IV inhibitor, the ATP levels decreased immediately. Interestingly, with azide present, the ROS activity remained low and the astrocytes remained viable even at very low ATP levels. Addition of exogenous ROS-scavengers prevented the IA-induced ROS activity, the ATP levels were maintained and cell death was prevented. Similar protection could be obtained when astrocytes, prior to addition of IA, were incubated with substances known to activate the nuclear factor erythroid 2-related factor 2 (Nrf2)-regulated endogenous antioxidant system. When IA was washed out, after a relatively moderate ATP depression, massive cell death occurred. This was efficiently prevented by addition of azide or ROS scavengers during the IA treatment or by pre-activation of the Nrf2 system. Our results demonstrate that astrocytes in culture can endure and recover from glycolytic inhibition if the ROS activity remained at a low level and suggest that oxidative stress can be an important component for astrocytic cell death following metabolic stress.     

The RBCC gene RFP2 (Leu5) encodes a novel transmembrane E3 ubiquitin ligase involved in ERAD

RFP2, a gene frequently lost in various malignancies, encodes a protein with RING finger, B-box, and coiled-coil domains that belongs to the RBCC/TRIM family of proteins. Here we demonstrate that Rfp2 is an unstable protein with auto-polyubiquitination activity in vivo and in vitro, implying that Rfp2 acts as a RING E3 ubiquitin ligase. Consequently, Rfp2 ubiquitin ligase activity is dependent on an intact RING domain, as RING deficient mutants fail to drive polyubiquitination in vitro and are stabilized in vivo. Immunopurification and tandem mass spectrometry enabled the identification of several putative Rfp2 interacting proteins localized to the endoplasmic reticulum (ER), including valosin-containing protein (VCP), a protein indispensable for ER-associated degradation (ERAD). Importantly, we also show that Rfp2 regulates the degradation of the known ER proteolytic substrate CD3-delta, but not the N-end rule substrate Ub-R-YFP (yellow fluorescent protein), establishing Rfp2 as a novel E3 ligase involved in ERAD. Finally, we show that Rfp2 contains a C-terminal transmembrane domain indispensable for its localization to the ER and that Rfp2 colocalizes with several ER-resident proteins as analyzed by high-resolution immunostaining. In summary, these data are all consistent with a function for Rfp2 as an ERAD E3 ubiquitin ligase.     

Structural bases of the redox-dependent conformational switch in the serpin PAI-2

Depending on the redox-status, the serpin plasminogen activator inhibitor type 2 (PAI-2) can exist in either a stable monomeric or polymerogenic form. The latter form, which spontaneously forms loop-sheet polymers, has an open beta-sheet A and is stabilized by a disulfide bond between C79 (in the CD-loop) and C161 (at the bottom of PAI-2). Reduction of this bond results in a closing of the beta-sheet A and converts PAI-2 to a stable monomeric form. Here we show that the stable monomeric and polymerogenic forms of PAI-2 are fully interconvertible, depending on redox-status of the environment. Our intramolecular distance measurements indicate that the CD-loop folds mainly on one side of the stable monomeric form of the inhibitor. However, the loop can translocate about 54A to the bottom of PAI-2 so that the C79-C161 disulfide bond can form under oxidizing conditions. We show also that the redox-active C79 can form a disulfide-link to the matrix protein vitronectin, suggesting that vitronectin can stabilize active PAI-2 in extracellular compartments. PAI-2 is therefore a rare example of a redox-sensitive protein for which the activity and polymerization ability are regulated by reversible disulfide bond formation leading to major translocation of a loop and significant conformational changes in the molecule.     

Personalized medicine and proteomics: lessons from non-small cell lung cancer

Personalized medicine allows the selection of treatments best suited to an individual patient and disease phenotype. To implement personalized medicine, effective tests predictive of response to treatment or susceptibility to adverse events are needed, and to develop a personalized medicine test, both high quality samples and reliable data are required. We review key features of state-of-the-art proteomic profiling and introduce further analytic developments to build a proteomic toolkit for use in personalized medicine approaches. The combination of novel analytical approaches in proteomic data generation, alignment and comparison permit translation of identified biomarkers into practical assays. We further propose an expanded statistical analysis to understand the sources of variability between individuals in terms of both protein expression and clinical variables and utilize this understanding in a predictive test.     

Escherichia coli TolA tolerates multiple amino-acid substitutions as revealed by screening randomized variants for membrane integrity and phage receptor function

Escherichia coli TolA is a cytoplasmic membrane protein required for outer membrane integrity and the translocation of F-specific filamentous (Ff) bacteriophage DNA. Both phage infection and membrane integrity depend on several TolA interactions, e.g. those of the TolA C-terminal domain (TolAIII). Membrane integrity involves interaction with two host proteins and phage translocation requires direct interaction with the N-terminal domain (N1) of Ff phage protein g3p. Although cocrystallization of TolAIII and N1g3p has identified several contact points, it is still uncertain which residues are selectively involved in the different TolA functions. Thus, four different limited substitution libraries of TolA were created, targeting contacts at positions 415-420. These libraries were introduced into the tolA strain K17DE3tolA/F(+) and several variants, containing complementing, multiple amino-acid substitutions, were identified. However, most randomized variants did not complement the tolA strain K17DE3tolA/F(+). The TolA variants that restored sensitivity to phage infection displayed a considerable sequence variation, while the few variants that restored tolerance to detergent were from the same library. A comparison of the generated residue variation and natural variation, suggests that structural dependence overrides contact residue dependence. Thus, library screening can be efficient in identifying TolA variants with different functionally associated characteristics.     

A Naturally Occurring Variant in Human TLR9, P99L, Is Associated with Loss of CpG Oligonucleotide Responsiveness

The innate immune system employs Toll-like receptors (TLRs) for the detection of invading microorganisms based on distinct molecular patterns. For example, TLR9 is activated by microbial DNA and also by short therapeutic CpG-containing oligonucleotides (CpG-ODN). TLR9 activation leads to the production of interferons and the priming of humoral adaptive immune responses. Unfortunately, the principles of ligand recognition by TLR9 are poorly understood, and genetic variants of TLR9, which may affect its function, have not been characterized systematically on the molecular level. We therefore sought to functionally characterize reported single nucleotide polymorphisms of TLR9 in the HEK293 model system. We discovered that two variants, P99L and M400I, are associated with altered receptor function regarding NF-κB activation and cytokine induction. Our investigations show that for the most functionally impaired variant, P99L, the ability to respond to physiological and therapeutic TLR9 ligands is severely compromised. However, CpG-ODN binding is normal. CpG-ODN recognition by TLR9 thus appears to involve two separate events, CpG-ODN binding and sensing. Our studies highlight Pro-99 as a residue important for the latter process. In genotyping studies, we confirmed that both M400I (rs41308230) and P99L (rs5743844) are relatively rare variants of TLR9. Our data add rs41308230 and rs5743844 to the list of functionally important TLR variants and warrant further research into their relevance for infectious disease susceptibility or responsiveness to CpG-ODN-based therapies.