Search algorithms as a framework for the optimization of drug combinations

Combination therapies are often needed for effective clinical outcomes in the management of complex diseases, but presently they are generally based on empirical clinical experience. Here we suggest a novel application of search algorithms, originally developed for digital communication, modified to optimize combinations of therapeutic interventions. In biological experiments measuring the restoration of the decline with age in heart function and exercise capacity in Drosophila melanogaster, we found that search algorithms correctly identified optimal combinations of four drugs with only one third of the tests performed in a fully factorial search. In experiments identifying combinations of three doses of up to six drugs for selective killing of human cancer cells, search algorithms resulted in a highly significant enrichment of selective combinations compared with random searches. In simulations using a network model of cell death, we found that the search algorithms identified the optimal combinations of 6-9 interventions in 80-90% of tests, compared with 15-30% for an equivalent random search. These findings suggest that modified search algorithms from information theory have the potential to enhance the discovery of novel therapeutic drug combinations. This report also helps to frame a biomedical problem that will benefit from an interdisciplinary effort and suggests a general strategy for its solution.Comment: 36 pages, 10 figures, revised versio

A genome-wide association study identifies protein quantitative trait loci (pQTLs)

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts - cis effects, and elsewhere in the genome - trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8×10 -57), CCL4L1 (p = 3.9×10-21), IL18 (p = 6.8×10-13), LPA (p = 4.4×10-10), GGT1 (p = 1.5×10-7), SHBG (p = 3.1×10-7), CRP (p = 6.4×10-6) and IL1RN (p = 7.3×10-6) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8×10-40), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways. © 2008 Melzer et al

Population-Based Resequencing of Experimentally Evolved Populations Reveals the Genetic Basis of Body Size Variation in Drosophila melanogaster

Body size is a classic quantitative trait with evolutionarily significant variation within many species. Locating the alleles responsible for this variation would help understand the maintenance of variation in body size in particular, as well as quantitative traits in general. However, successful genome-wide association of genotype and phenotype may require very large sample sizes if alleles have low population frequencies or modest effects. As a complementary approach, we propose that population-based resequencing of experimentally evolved populations allows for considerable power to map functional variation. Here, we use this technique to investigate the genetic basis of natural variation in body size in Drosophila melanogaster. Significant differentiation of hundreds of loci in replicate selection populations supports the hypothesis that the genetic basis of body size variation is very polygenic in D. melanogaster. Significantly differentiated variants are limited to single genes at some loci, allowing precise hypotheses to be formed regarding causal polymorphisms, while other significant regions are large and contain many genes. By using significantly associated polymorphisms as a priori candidates in follow-up studies, these data are expected to provide considerable power to determine the genetic basis of natural variation in body size

Testing the Mid-Holocene Relative Sea-Level Highstand Hypothesis in North Wales, United Kingdom

Accurate Holocene relative sea-level curves are vital for modelling future sea-level changes, particularly in regions where relative sea-level changes are dominated by isostatically induced vertical land movements. In North Wales, various glacial isostatic adjustment (GIA) models predict a mid-Holocene relative sea-level highstand between 4 and 6 ka, which is unsubstantiated by any geological sea-level data but affects the ability of geophysical models to model accurately past and future sea levels. Here, we use a newly developed foraminifera-based sea-level transfer function to produce a 3300-year-long late-Holocene relative sea-level reconstruction from a salt marsh in the Malltraeth estuary on the south Anglesey coast in North Wales. This is the longest continuous late-Holocene relative sea-level reconstruction in Northwest Europe. We combine this record with two new late-Holocene sea-level index points (SLIPs) obtained from a freshwater marsh at Rhoscolyn, Anglesey, and with previously published regional SLIPs, to produce a relative sea-level record for North Wales that spans from ca. 13,000 BP to the present. This record leaves no room for a mid-Holocene relative sea-level highstand in the region. We conclude that GIA models that include a mid-Holocene sea-level highstand for North Wales need revision before they are used in the modelling of past and future relative sea-level changes around the British Isles

A Cross-Study Transcriptional Analysis of Parkinson's Disease

The study of Parkinson's disease (PD), like other complex neurodegenerative disorders, is limited by access to brain tissue from patients with a confirmed diagnosis. Alternatively the study of peripheral tissues may offer some insight into the molecular basis of disease susceptibility and progression, but this approach still relies on brain tissue to benchmark relevant molecular changes against. Several studies have reported whole-genome expression profiling in post-mortem brain but reported concordance between these analyses is lacking. Here we apply a standardised pathway analysis to seven independent case-control studies, and demonstrate increased concordance between data sets. Moreover data convergence increased when the analysis was limited to the five substantia nigra (SN) data sets; this highlighted the down regulation of dopamine receptor signaling and insulin-like growth factor 1 (IGF1) signaling pathways. We also show that case-control comparisons of affected post mortem brain tissue are more likely to reflect terminal cytoarchitectural differences rather than primary pathogenic mechanisms. The implementation of a correction factor for dopaminergic neuronal loss predictably resulted in the loss of significance of the dopamine signaling pathway while axon guidance pathways increased in significance. Interestingly the IGF1 signaling pathway was also over-represented when data from non-SN areas, unaffected or only terminally affected in PD, were considered. Our findings suggest that there is greater concordance in PD whole-genome expression profiling when standardised pathway membership rather than ranked gene list is used for comparison

Genome-Wide Fitness and Expression Profiling Implicate Mga2 in Adaptation to Hydrogen Peroxide

Caloric restriction extends lifespan, an effect once thought to involve attenuation of reactive oxygen species (ROS) generated by aerobic metabolism. However, recent evidence suggests that caloric restriction may in fact raise ROS levels, which in turn provides protection from acute doses of oxidant through a process called adaptation. To shed light on the molecular mechanisms of adaptation, we designed a series of genome-wide deletion fitness and mRNA expression screens to identify genes involved in adaptation to hydrogen peroxide. Combined with known transcriptional interactions, the integrated data implicate Yap1 and Skn7 as central transcription factors of both the adaptive and acute oxidative responses. They also identify the transcription factors Mga2 and Rox1 as active exclusively in the adaptive response and show that Mga2 is essential for adaptation. These findings are striking because Mga2 and Rox1 have been thought to control the response to hypoxic, not oxidative, conditions. Expression profiling of mga2Δ and rox1Δ knockouts shows that these factors most strongly regulate targets in ergosterol, fatty-acid, and zinc metabolic pathways. Direct quantitation of ergosterol reveals that its basal concentration indeed depends on Mga2, but that Mga2 is not required for the decrease in ergosterol observed during adaptation

Functional Copy-Number Alterations in Cancer

Understanding the molecular basis of cancer requires characterization of its genetic defects. DNA microarray technologies can provide detailed raw data about chromosomal aberrations in tumor samples. Computational analysis is needed (1) to deduce from raw array data actual amplification or deletion events for chromosomal fragments and (2) to distinguish causal chromosomal alterations from functionally neutral ones. We present a comprehensive computational approach, RAE, designed to robustly map chromosomal alterations in tumor samples and assess their functional importance in cancer. To demonstrate the methodology, we experimentally profile copy number changes in a clinically aggressive subtype of soft-tissue sarcoma, pleomorphic liposarcoma, and computationally derive a portrait of candidate oncogenic alterations and their target genes. Many affected genes are known to be involved in sarcomagenesis; others are novel, including mediators of adipocyte differentiation, and may include valuable therapeutic targets. Taken together, we present a statistically robust methodology applicable to high-resolution genomic data to assess the extent and function of copy-number alterations in cancer

Using lithium as a neuroprotective agent in patients with cancer

Neurocognitive impairment is being increasingly recognized as an important issue in patients with cancer who develop cognitive difficulties either as part of direct or indirect involvement of the nervous system or as a consequence of either chemotherapy-related or radiotherapy-related complications. Brain radiotherapy in particular can lead to significant cognitive defects. Neurocognitive decline adversely affects quality of life, meaningful employment, and even simple daily activities. Neuroprotection may be a viable and realistic goal in preventing neurocognitive sequelae in these patients, especially in the setting of cranial irradiation. Lithium is an agent that has been in use for psychiatric disorders for decades, but recently there has been emerging evidence that it can have a neuroprotective effect.This review discusses neurocognitive impairment in patients with cancer and the potential for investigating the use of lithium as a neuroprotectant in such patients.<br /

Using lithium as a neuroprotective agent in patients with cancer

Neurocognitive impairment is being increasingly recognized as an important issue in patients with cancer who develop cognitive difficulties either as part of direct or indirect involvement of the nervous system or as a consequence of either chemotherapy-related or radiotherapy-related complications. Brain radiotherapy in particular can lead to significant cognitive defects. Neurocognitive decline adversely affects quality of life, meaningful employment, and even simple daily activities. Neuroprotection may be a viable and realistic goal in preventing neurocognitive sequelae in these patients, especially in the setting of cranial irradiation. Lithium is an agent that has been in use for psychiatric disorders for decades, but recently there has been emerging evidence that it can have a neuroprotective effect.This review discusses neurocognitive impairment in patients with cancer and the potential for investigating the use of lithium as a neuroprotectant in such patients.<br /

Transcriptional Profiling of the Dose Response: A More Powerful Approach for Characterizing Drug Activities

The dose response curve is the gold standard for measuring the effect of a drug treatment, but is rarely used in genomic scale transcriptional profiling due to perceived obstacles of cost and analysis. One barrier to examining transcriptional dose responses is that existing methods for microarray data analysis can identify patterns, but provide no quantitative pharmacological information. We developed analytical methods that identify transcripts responsive to dose, calculate classical pharmacological parameters such as the EC50, and enable an in-depth analysis of coordinated dose-dependent treatment effects. The approach was applied to a transcriptional profiling study that evaluated four kinase inhibitors (imatinib, nilotinib, dasatinib and PD0325901) across a six-logarithm dose range, using 12 arrays per compound. The transcript responses proved a powerful means to characterize and compare the compounds: the distribution of EC50 values for the transcriptome was linked to specific targets, dose-dependent effects on cellular processes were identified using automated pathway analysis, and a connection was seen between EC50s in standard cellular assays and transcriptional EC50s. Our approach greatly enriches the information that can be obtained from standard transcriptional profiling technology. Moreover, these methods are automated, robust to non-optimized assays, and could be applied to other sources of quantitative data