A simple and robust method for connecting small-molecule drugs using gene-expression signatures

Interaction of a drug or chemical with a biological system can result in a gene-expression profile or signature characteristic of the event. Using a suitably robust algorithm these signatures can potentially be used to connect molecules with similar pharmacological or toxicological properties. The Connectivity Map was a novel concept and innovative tool first introduced by Lamb et al to connect small molecules, genes, and diseases using genomic signatures [Lamb et al (2006), Science 313, 1929-1935]. However, the Connectivity Map had some limitations, particularly there was no effective safeguard against false connections if the observed connections were considered on an individual-by-individual basis. Further when several connections to the same small-molecule compound were viewed as a set, the implicit null hypothesis tested was not the most relevant one for the discovery of real connections. Here we propose a simple and robust method for constructing the reference gene-expression profiles and a new connection scoring scheme, which importantly allows the valuation of statistical significance of all the connections observed. We tested the new method with the two example gene-signatures (HDAC inhibitors and Estrogens) used by Lamb et al and also a new gene signature of immunosuppressive drugs. Our testing with this new method shows that it achieves a higher level of specificity and sensitivity than the original method. For example, our method successfully identified raloxifene and tamoxifen as having significant anti-estrogen effects, while Lamb et al's Connectivity Map failed to identify these. With these properties our new method has potential use in drug development for the recognition of pharmacological and toxicological properties in new drug candidates.Comment: 8 pages, 2 figures, and 2 tables; supplementary data supplied as a ZIP fil

How consistent are the transcriptome changes associated with cold acclimation in two species of the Drosophila virilis group?

This work was financially support by a Marie Curie Initial Training Network grant, “Understanding the evolutionary origin of biological diversity” (ITN-2008–213780 SPECIATION), grants from the Academy of Finland to A.H. (project 132619) and M.K. (projects 268214 and 272927), a grant from NERC, UK to M.G.R. (grant NE/J020818/1), and NERC, UK PhD studentship to D.J.P. (NE/I528634/1).For many organisms the ability to cold acclimate with the onset of seasonal cold has major implications for their fitness. In insects, where this ability is widespread, the physiological changes associated with increased cold tolerance have been well studied. Despite this, little work has been done to trace changes in gene expression during cold acclimation that lead to an increase in cold tolerance. We used an RNA-Seq approach to investigate this in two species of the Drosophila virilis group. We found that the majority of genes that are differentially expressed during cold acclimation differ between the two species. Despite this, the biological processes associated with the differentially expressed genes were broadly similar in the two species. These included: metabolism, cell membrane composition, and circadian rhythms, which are largely consistent with previous work on cold acclimation/cold tolerance. In addition, we also found evidence of the involvement of the rhodopsin pathway in cold acclimation, a pathway that has been recently linked to thermotaxis. Interestingly, we found no evidence of differential expression of stress genes implying that long-term cold acclimation and short-term stress response may have a different physiological basis.PostprintPeer reviewe

Siderophore production by Bacillus megaterium : effect of growth-phase and cultural conditions

Siderophore production by Bacillus megaterium was detected, in an iron-deficient culture medium, during the exponential growth phase, prior to the sporulation, in the presence of glucose; these results suggested that the onset of siderophore production did not require glucose depletion and was not related with the sporulation. The siderophore production by B. megaterium was affected by the carbon source used. The growth on glycerol promoted the very high siderophore production (1,182 μmol g−1 dry weight biomass); the opposite effect was observed in the presence of mannose (251 μmol g−1 dry weight biomass). The growth in the presence of fructose, galactose, glucose, lactose, maltose or sucrose, originated similar concentrations of siderophore (546–842 μmol g−1 dry weight biomass). Aeration had a positive effect on the production of siderophore. Incubation of B. megaterium under static conditions delayed and reduced the growth and the production of siderophore, compared with the incubation in stirred conditions.The authors thank Porto University/Totta Bank for their financial support through the project "Microbiological production of chelating agents" (Ref: 180). The authors also thank the Fundacao para a Ciencia e a Tecnologia (FCT) through the Portuguese Government for their financial support of this work through the grants Strategic project-LA23/2013-2014 (IBB) and PEST-C/EQB/LA0006/2011 (REQUIMTE). Manuela D. Machado gratefully acknowledges the postdoctoral (SFRH/BPD/72816/2010) grant from FCT

Parallel molecular routes to cold adaptation in eight genera of New Zealand stick insects.

The acquisition of physiological strategies to tolerate novel thermal conditions allows organisms to exploit new environments. As a result, thermal tolerance is a key determinant of the global distribution of biodiversity, yet the constraints on its evolution are not well understood. Here we investigate parallel evolution of cold tolerance in New Zealand stick insects, an endemic radiation containing three montane-occurring species. Using a phylogeny constructed from 274 orthologous genes, we show that stick insects have independently colonized montane environments at least twice. We compare supercooling point and survival of internal ice formation among ten species from eight genera, and identify both freeze tolerance and freeze avoidance in separate montane lineages. Freeze tolerance is also verified in both lowland and montane populations of a single, geographically widespread, species. Transcriptome sequencing following cold shock identifies a set of structural cuticular genes that are both differentially regulated and under positive sequence selection in each species. However, while cuticular proteins in general are associated with cold shock across the phylogeny, the specific genes at play differ among species. Thus, while processes related to cuticular structure are consistently associated with adaptation for cold, this may not be the consequence of shared ancestral genetic constraints

Semi-supervised discovery of differential genes

BACKGROUND: Various statistical scores have been proposed for evaluating the significance of genes that may exhibit differential expression between two or more controlled conditions. However, in many clinical studies to detect clinical marker genes for example, the conditions have not necessarily been controlled well, thus condition labels are sometimes hard to obtain due to physical, financial, and time costs. In such a situation, we can consider an unsupervised case where labels are not available or a semi-supervised case where labels are available for a part of the whole sample set, rather than a well-studied supervised case where all samples have their labels. RESULTS: We assume a latent variable model for the expression of active genes and apply the optimal discovery procedure (ODP) proposed by Storey (2005) to the model. Our latent variable model allows gene significance scores to be applied to unsupervised and semi-supervised cases. The ODP framework improves detectability by sharing the estimated parameters of null and alternative models of multiple tests over multiple genes. A theoretical consideration leads to two different interpretations of the latent variable, i.e., it only implicitly affects the alternative model through the model parameters, or it is explicitly included in the alternative model, so that the interpretations correspond to two different implementations of ODP. By comparing the two implementations through experiments with simulation data, we have found that sharing the latent variable estimation is effective for increasing the detectability of truly active genes. We also show that the unsupervised and semi-supervised rating of genes, which takes into account the samples without condition labels, can improve detection of active genes in real gene discovery problems. CONCLUSION: The experimental results indicate that the ODP framework is effective for hypotheses including latent variables and is further improved by sharing the estimations of hidden variables over multiple tests

Antiinflammatory Therapy with Canakinumab for Atherosclerotic Disease

Background: Experimental and clinical data suggest that reducing inflammation without affecting lipid levels may reduce the risk of cardiovascular disease. Yet, the inflammatory hypothesis of atherothrombosis has remained unproved. Methods: We conducted a randomized, double-blind trial of canakinumab, a therapeutic monoclonal antibody targeting interleukin-1β, involving 10,061 patients with previous myocardial infarction and a high-sensitivity C-reactive protein level of 2 mg or more per liter. The trial compared three doses of canakinumab (50 mg, 150 mg, and 300 mg, administered subcutaneously every 3 months) with placebo. The primary efficacy end point was nonfatal myocardial infarction, nonfatal stroke, or cardiovascular death. RESULTS: At 48 months, the median reduction from baseline in the high-sensitivity C-reactive protein level was 26 percentage points greater in the group that received the 50-mg dose of canakinumab, 37 percentage points greater in the 150-mg group, and 41 percentage points greater in the 300-mg group than in the placebo group. Canakinumab did not reduce lipid levels from baseline. At a median follow-up of 3.7 years, the incidence rate for the primary end point was 4.50 events per 100 person-years in the placebo group, 4.11 events per 100 person-years in the 50-mg group, 3.86 events per 100 person-years in the 150-mg group, and 3.90 events per 100 person-years in the 300-mg group. The hazard ratios as compared with placebo were as follows: in the 50-mg group, 0.93 (95% confidence interval [CI], 0.80 to 1.07; P = 0.30); in the 150-mg group, 0.85 (95% CI, 0.74 to 0.98; P = 0.021); and in the 300-mg group, 0.86 (95% CI, 0.75 to 0.99; P = 0.031). The 150-mg dose, but not the other doses, met the prespecified multiplicity-adjusted threshold for statistical significance for the primary end point and the secondary end point that additionally included hospitalization for unstable angina that led to urgent revascularization (hazard ratio vs. placebo, 0.83; 95% CI, 0.73 to 0.95; P = 0.005). Canakinumab was associated with a higher incidence of fatal infection than was placebo. There was no significant difference in all-cause mortality (hazard ratio for all canakinumab doses vs. placebo, 0.94; 95% CI, 0.83 to 1.06; P = 0.31). Conclusions: Antiinflammatory therapy targeting the interleukin-1β innate immunity pathway with canakinumab at a dose of 150 mg every 3 months led to a significantly lower rate of recurrent cardiovascular events than placebo, independent of lipid-level lowering. (Funded by Novartis; CANTOS ClinicalTrials.gov number, NCT01327846.

In Vivo Assessment of Cold Adaptation in Insect Larvae by Magnetic Resonance Imaging and Magnetic Resonance Spectroscopy

Background Temperatures below the freezing point of water and the ensuing ice crystal formation pose serious challenges to cell structure and function. Consequently, species living in seasonally cold environments have evolved a multitude of strategies to reorganize their cellular architecture and metabolism, and the underlying mechanisms are crucial to our understanding of life. In multicellular organisms, and poikilotherm animals in particular, our knowledge about these processes is almost exclusively due to invasive studies, thereby limiting the range of conclusions that can be drawn about intact living systems. Methodology Given that non-destructive techniques like 1H Magnetic Resonance (MR) imaging and spectroscopy have proven useful for in vivo investigations of a wide range of biological systems, we aimed at evaluating their potential to observe cold adaptations in living insect larvae. Specifically, we chose two cold-hardy insect species that frequently serve as cryobiological model systems–the freeze-avoiding gall moth Epiblema scudderiana and the freeze-tolerant gall fly Eurosta solidaginis. Results In vivo MR images were acquired from autumn-collected larvae at temperatures between 0°C and about -70°C and at spatial resolutions down to 27 µm. These images revealed three-dimensional (3D) larval anatomy at a level of detail currently not in reach of other in vivo techniques. Furthermore, they allowed visualization of the 3D distribution of the remaining liquid water and of the endogenous cryoprotectants at subzero temperatures, and temperature-weighted images of these distributions could be derived. Finally, individual fat body cells and their nuclei could be identified in intact frozen Eurosta larvae. Conclusions These findings suggest that high resolution MR techniques provide for interesting methodological options in comparative cryobiological investigations, especially in vivo

Enzymatic capacities of metabolic fuel use in cuttlefish (Sepia officinalis) and responses to food deprivation: insight into the metabolic organization and starvation survival strategy of cephalopods

Food limitation is a common challenge for animals. Cephalopods are sensitive to starvation because of high metabolic rates and growth rates related to their "live fast, die young" life history. We investigated how enzymatic capacities of key metabolic pathways are modulated during starvation in the common cuttlefish (Sepia officinalis) to gain insight into the metabolic organization of cephalopods and their strategies for coping with food limitation. In particular, lipids have traditionally been considered unimportant fuels in cephalopods, yet, puzzlingly, many species (including cuttlefish) mobilize the lipid stores in their digestive gland during starvation. Using a comprehensive multi-tissue assay of enzymatic capacities for energy metabolism, we show that, during long-term starvation (12 days), glycolytic capacity for glucose use is decreased in cuttlefish tissues, while capacities for use of lipid-based fuels (fatty acids and ketone bodies) and amino acid fuels are retained or increased. Specifically, the capacity to use the ketone body acetoacetate as fuel is widespread across tissues and gill has a previously unrecognized capacity for fatty acid catabolism, albeit at low rates. The capacity for de novo glucose synthesis (gluconeogenesis), important for glucose homeostasis, likely is restricted to the digestive gland, contrary to previous reports of widespread gluconeogenesis among cephalopod tissues. Short-term starvation (3-5 days) had few effects on enzymatic capacities. Similar to vertebrates, lipid-based fuels, putatively mobilized from fat stores in the digestive gland, appear to be important energy sources for cephalopods, especially during starvation when glycolytic capacity is decreased perhaps to conserve available glucose

On the Wegener granulomatosis associated region on chromosome 6p21.3

BACKGROUND: Wegener granulomatosis (WG) belongs to the heterogeneous group of systemic vasculitides. The multifactorial pathophysiology of WG is supposedly caused by yet unknown environmental influence(s) on the basis of genetic predisposition. The presence of anti-neutrophil cytoplasmic antibodies (ANCA) in the plasma of patients and genetic involvement of the human leukocyte antigen system reflect an autoimmune background of the disease. Strong associations were revealed with WG by markers located in the major histocompatibility complex class II (MHC II) region in the vicinity of human leukocyte antigen (HLA)-DPB1 and the retinoid X receptor B (RXRB) loci. In order to define the involvement of the 6p21.3 region in WG in more detail this previous population-based association study was expanded here to the respective 3.6 megabase encompassing this region on chromosome 6. The RXRB gene was analysed as well as a splice-site variation of the butyrophilin-like (BTNL2) gene which is also located within the respective region. The latter polymorphism has been evaluated here as it appears as a HLA independent susceptibility factor in another granulomatous disorder, sarcoidosis. METHODS: 150–180 German WG patients and a corresponding cohort of healthy controls (n = 100–261) were used in a two-step study. A panel of 94 microsatellites was designed for the initial step using a DNA pooling approach. Markers with significantly differing allele frequencies between patient and control pools were individually genotyped. The RXRB gene was analysed for single strand conformation polymorphisms (SSCP) and restriction fragment length polymorphisms (RFLP). The splice-site polymorphism in the BTNL2 gene was also investigated by RFLP analysis. RESULTS: A previously investigated microsatellite (#1.0.3.7, Santa Cruz genome browser (UCSC) May 2004 Freeze localisation: chr6:31257596-34999883), which was used as a positive control, remained associated throughout the whole two-step approach. Yet, no additional evidence for association of other microsatellite markers was found in the entire investigated region. Analysis of the RXRB gene located in the WG associated region revealed associations of two variations (rs10548957 p(allelic )= 0.02 and rs6531 p(allelic )= 5.20 × 10(-5), OR = 1.88). Several alleles of markers located between HLA-DPB1, SNP rs6531 and microsatellite 1.0.3.7 showed linkage disequilibrium with r(2 )values exceeding 0.10. Significant differences were not demonstrable for the sarcoidosis associated splice-site variation (rs2076530 p(allelic )= 0.80) in our WG cohort. CONCLUSION: Since a microsatellite flanking the RXRB gene and two intragenic polymorphisms are associated significantly with WG on chromosome 6p21.3, further investigations should be focussed on extensive fine-mapping in this region by densely mapping with additional markers such as SNPs. This strategy may reveal even deeper insights into the genetic contributions of the respective region for the pathogenesis of WG