Visualizing dimensionality reduction of systems biology data

One of the challenges in analyzing high-dimensional expression data is the detection of important biological signals. A common approach is to apply a dimension reduction method, such as principal component analysis. Typically, after application of such a method the data is projected and visualized in the new coordinate system, using scatter plots or profile plots. These methods provide good results if the data have certain properties which become visible in the new coordinate system and which were hard to detect in the original coordinate system. Often however, the application of only one method does not suffice to capture all important signals. Therefore several methods addressing different aspects of the data need to be applied. We have developed a framework for linear and non-linear dimension reduction methods within our visual analytics pipeline SpRay. This includes measures that assist the interpretation of the factorization result. Different visualizations of these measures can be combined with functional annotations that support the interpretation of the results. We show an application to high-resolution time series microarray data in the antibiotic-producing organism Streptomyces coelicolor as well as to microarray data measuring expression of cells with normal karyotype and cells with trisomies of human chromosomes 13 and 21

Concept and optical design of the cross-disperser module for CRIRES

This is the peer reviewed version of the following article: Oliva, Ernesto, A. Tozzi, D. Ferruzzi, L. Origlia, A. Hatzes, R. Follert, T. Loewinger et al. "Concept and optical design of the cross-disperser module for CRIRES+." In SPIE Astronomical Telescopes+ Instrumentation, pp. 91477R-91477R. International Society for Optics and Photonics, 2014, which has been published in final form at 10.1117/12.2054381

An effective all-atom potential for proteins

We describe and test an implicit solvent all-atom potential for simulations of protein folding and aggregation. The potential is developed through studies of structural and thermodynamic properties of 17 peptides with diverse secondary structure. Results obtained using the final form of the potential are presented for all these peptides. The same model, with unchanged parameters, is furthermore applied to a heterodimeric coiled-coil system, a mixed alpha/beta protein and a three-helix-bundle protein, with very good results. The computational efficiency of the potential makes it possible to investigate the free-energy landscape of these 49--67-residue systems with high statistical accuracy, using only modest computational resources by today's standards

B Cells Regulate Neutrophilia during Mycobacterium tuberculosis Infection and BCG Vaccination by Modulating the Interleukin-17 Response

We have previously demonstrated that B cells can shape the immune response to Mycobacterium tuberculosis, including the level of neutrophil infiltration and granulomatous inflammation at the site of infection. The present study examined the mechanisms by which B cells regulate the host neutrophilic response upon exposure to mycobacteria and how neutrophilia may influence vaccine efficacy. To address these questions, a murine aerosol infection tuberculosis (TB) model and an intradermal (ID) ear BCG immunization mouse model, involving both the μMT strain and B cell-depleted C57BL/6 mice, were used. IL (interleukin)-17 neutralization and neutrophil depletion experiments using these systems provide evidence that B cells can regulate neutrophilia by modulating the IL-17 response during M. tuberculosis infection and BCG immunization. Exuberant neutrophilia at the site of immunization in B cell-deficient mice adversely affects dendritic cell (DC) migration to the draining lymph nodes and attenuates the development of the vaccine-induced Th1 response. The results suggest that B cells are required for the development of optimal protective anti-TB immunity upon BCG vaccination by regulating the IL-17/neutrophilic response. Administration of sera derived from M. tuberculosis-infected C57BL/6 wild-type mice reverses the lung neutrophilia phenotype in tuberculous μMT mice. Together, these observations provide insight into the mechanisms by which B cells and humoral immunity modulate vaccine-induced Th1 response and regulate neutrophila during M. tuberculosis infection and BCG immunization. © 2013 Kozakiewicz et al

Cross-species gene expression analysis of species specific differences in the preclinical assessment of pharmaceutical compounds

Animals are frequently used as model systems for determination of safety and efficacy in pharmaceutical research and development. However, significant quantitative and qualitative differences exist between humans and the animal models used in research. This is as a result of genetic variation between human and the laboratory animal. Therefore the development of a system that would allow the assessment of all molecular differences between species after drug exposure would have a significant impact on drug evaluation for toxicity and efficacy. Here we describe a cross-species microarray methodology that identifies and selects orthologous probes after cross-species sequence comparison to develop an orthologous cross-species gene expression analysis tool. The assumptions made by the use of this orthologous gene expression strategy for cross-species extrapolation is that; conserved changes in gene expression equate to conserved pharmacodynamic endpoints. This assumption is supported by the fact that evolution and selection have maintained the structure and function of many biochemical pathways over time, resulting in the conservation of many important processes. We demonstrate this cross-species methodology by investigating species specific differences of the peroxisome proliferatoractivator receptor (PPAR) a response in rat and human

A large-scale study of a poultry trading network in Bangladesh: implications for control and surveillance of avian influenza viruses

Since its first report in 2007, avian influenza (AI) has been endemic in Bangladesh. While live poultry marketing is widespread throughout the country and known to influence AI dissemination and persistence, trading patterns have not been described. The aim of this study is to assess poultry trading practices and features of the poultry trading networks which could promote AI spread, and their potential implications for disease control and surveillance. Data on poultry trading practices was collected from 849 poultry traders during a cross-sectional survey in 138 live bird markets (LBMs) across 17 different districts of Bangladesh. The quantity and origins of traded poultry were assessed for each poultry type in surveyed LBMs. The network of contacts between farms and LBMs resulting from commercial movements of live poultry was constructed to assess its connectivity and to identify the key premises influencing it

A finite strain fibre-reinforced viscoelasto-viscoplastic model of plant cell wall growth

A finite strain fibre-reinforced viscoelasto-viscoplastic model implemented in a finite element (FE) analysis is presented to study the expansive growth of plant cell walls. Three components of the deformation of growing cell wall, i.e. elasticity, viscoelasticity and viscoplasticity-like growth, are modelled within a consistent framework aiming to present an integrative growth model. The two aspects of growth—turgor-driven creep and new material deposition—and the interplay between them are considered by presenting a yield function, flow rule and hardening law. A fibre-reinforcement formulation is used to account for the role of cellulose microfibrils in the anisotropic growth. Mechanisms in in vivo growth are taken into account to represent the corresponding biologycontrolled behaviour of a cell wall. A viscoelastic formulation is proposed to capture the viscoelastic response in the cell wall. The proposed constitutive model provides a unique framework for modelling both the in vivo growth of cell wall dominated by viscoplasticity-like behaviour and in vitro deformation dominated by elastic or viscoelastic responses. A numerical scheme is devised, and FE case studies are reported and compared with experimental data

Quantum dynamics in strong fluctuating fields

A large number of multifaceted quantum transport processes in molecular systems and physical nanosystems can be treated in terms of quantum relaxation processes which couple to one or several fluctuating environments. A thermal equilibrium environment can conveniently be modelled by a thermal bath of harmonic oscillators. An archetype situation provides a two-state dissipative quantum dynamics, commonly known under the label of a spin-boson dynamics. An interesting and nontrivial physical situation emerges, however, when the quantum dynamics evolves far away from thermal equilibrium. This occurs, for example, when a charge transferring medium possesses nonequilibrium degrees of freedom, or when a strong time-dependent control field is applied externally. Accordingly, certain parameters of underlying quantum subsystem acquire stochastic character. Herein, we review the general theoretical framework which is based on the method of projector operators, yielding the quantum master equations for systems that are exposed to strong external fields. This allows one to investigate on a common basis the influence of nonequilibrium fluctuations and periodic electrical fields on quantum transport processes. Most importantly, such strong fluctuating fields induce a whole variety of nonlinear and nonequilibrium phenomena. A characteristic feature of such dynamics is the absence of thermal (quantum) detailed balance.Comment: review article, Advances in Physics (2005), in pres

Next-generation mitogenomics: A comparison of approaches applied to caecilian amphibian phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a 'traditional' Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina's HiSeq and MiSeq, Roche's 454 GS FLX, and Life Technologies' Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case

Weak temperature dependence of P (+) H A (-) recombination in mutant Rhodobacter sphaeroides reaction centers

International audienceIn contrast with findings on the wild-type Rhodobacter sphaeroides reaction center, biexponential P (+) H A (-) → PH A charge recombination is shown to be weakly dependent on temperature between 78 and 298 K in three variants with single amino acids exchanged in the vicinity of primary electron acceptors. These mutated reaction centers have diverse overall kinetics of charge recombination, spanning an average lifetime from ~2 to ~20 ns. Despite these differences a protein relaxation model applied previously to wild-type reaction centers was successfully used to relate the observed kinetics to the temporal evolution of the free energy level of the state P (+) H A (-) relative to P (+) B A (-) . We conclude that the observed variety in the kinetics of charge recombination, together with their weak temperature dependence, is caused by a combination of factors that are each affected to a different extent by the point mutations in a particular mutant complex. These are as follows: (1) the initial free energy gap between the states P (+) B A (-) and P (+) H A (-) , (2) the intrinsic rate of P (+) B A (-) → PB A charge recombination, and (3) the rate of protein relaxation in response to the appearance of the charge separated states. In the case of a mutant which displays rapid P (+) H A (-) recombination (ELL), most of this recombination occurs in an unrelaxed protein in which P (+) B A (-) and P (+) H A (-) are almost isoenergetic. In contrast, in a mutant in which P (+) H A (-) recombination is relatively slow (GML), most of the recombination occurs in a relaxed protein in which P (+) H A (-) is much lower in energy than P (+) H A (-) . The weak temperature dependence in the ELL reaction center and a YLH mutant was modeled in two ways: (1) by assuming that the initial P (+) B A (-) and P (+) H A (-) states in an unrelaxed protein are isoenergetic, whereas the final free energy gap between these states following the protein relaxation is large (~250 meV or more), independent of temperature and (2) by assuming that the initial and final free energy gaps between P (+) B A (-) and P (+) H A (-) are moderate and temperature dependent. In the case of the GML mutant, it was concluded that the free energy gap between P (+) B A (-) and P (+) H A (-) is large at all times