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

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

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

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

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

Effects of elevated seawater pCO2 on gene expression patterns in the gills of the green crab, Carcinus maenas

Background: The green crab Carcinus maenas is known for its high acclimation potential to varying environmental abiotic conditions. A high ability for ion and acid-base regulation is mainly based on an efficient regulation apparatus located in gill epithelia. However, at present it is neither known which ion transport proteins play a key role in the acid-base compensation response nor how gill epithelia respond to elevated seawater pCO2 as predicted for the future. In order to promote our understanding of the responses of green crab acid-base regulatory epithelia to high pCO2, Baltic Sea green crabs were exposed to a pCO2 of 400 Pa. Gills were screened for differentially expressed gene transcripts using a 4,462-feature microarray and quantitative real-time PCR. Results: Crabs responded mainly through fine scale adjustment of gene expression to elevated pCO2. However, 2% of all investigated transcripts were significantly regulated 1.3 to 2.2-fold upon one-week exposure to CO2 stress. Most of the genes known to code for proteins involved in osmo- and acid-base regulation, as well as cellular stress response, were were not impacted by elevated pCO2. However, after one week of exposure, significant changes were detected in a calcium-activated chloride channel, a hyperpolarization activated nucleotide-gated potassium channel, a tetraspanin, and an integrin. Furthermore, a putative syntaxin-binding protein, a protein of the transmembrane 9 superfamily, and a Cl-/HCO3 - exchanger of the SLC 4 family were differentially regulated. These genes were also affected in a previously published hypoosmotic acclimation response study. Conclusions: The moderate, but specific response of C. maenas gill gene expression indicates that (1) seawater acidification does not act as a strong stressor on the cellular level in gill epithelia; (2) the response to hypercapnia is to some degree comparable to a hypoosmotic acclimation response; (3) the specialization of each of the posterior gill arches might go beyond what has been demonstrated up to date; and (4) a re-configuration of gill epithelia might occur in response to hypercapnia

Oxygen Reperfusion Damage in an Insect

The deleterious effects of anoxia followed by reperfusion with oxygen in higher animals including mammals are well known. A convenient and genetically well characterized small-animal model that exhibits reproducible, quantifiable oxygen reperfusion damage is currently lacking. Here we describe the dynamics of whole-organism metabolic recovery from anoxia in an insect, Drosophila melanogaster, and report that damage caused by oxygen reperfusion can be quantified in a novel but straightforward way. We monitored CO2 emission (an index of mitochondrial activity) and water vapor output (an index of neuromuscular control of the spiracles, which are valves between the outside air and the insect's tracheal system) during entry into, and recovery from, rapid-onset anoxia exposure with durations ranging from 7.5 to 120 minutes. Anoxia caused a brief peak of CO2 output followed by knock-out. Mitochondrial respiration ceased and the spiracle constrictor muscles relaxed, but then re-contracted, presumably powered by anaerobic processes. Reperfusion to sustained normoxia caused a bimodal re-activation of mitochondrial respiration, and in the case of the spiracle constrictor muscles, slow inactivation followed by re-activation. After long anoxia durations, both the bimodality of mitochondrial reactivation and the recovery of spiracular control were impaired. Repeated reperfusion followed by episodes of anoxia depressed mitochondrial respiratory flux rates and damaged the integrity of the spiracular control system in a dose-dependent fashion. This is the first time that physiological evidence of oxygen reperfusion damage has been described in an insect or any invertebrate. We suggest that some of the traditional approaches of insect respiratory biology, such as quantifying respiratory water loss, may facilitate using D. melanogaster as a convenient, well-characterized experimental model for studying the underlying biology and mechanisms of ischemia and reperfusion damage and its possible mitigation