Integrating biological knowledge into variable selection : an empirical Bayes approach with an application in cancer biology

Background: An important question in the analysis of biochemical data is that of identifying subsets of molecular variables that may jointly influence a biological response. Statistical variable selection methods have been widely used for this purpose. In many settings, it may be important to incorporate ancillary biological information concerning the variables of interest. Pathway and network maps are one example of a source of such information. However, although ancillary information is increasingly available, it is not always clear how it should be used nor how it should be weighted in relation to primary data. Results: We put forward an approach in which biological knowledge is incorporated using informative prior distributions over variable subsets, with prior information selected and weighted in an automated, objective manner using an empirical Bayes formulation. We employ continuous, linear models with interaction terms and exploit biochemically-motivated sparsity constraints to permit exact inference. We show an example of priors for pathway- and network-based information and illustrate our proposed method on both synthetic response data and by an application to cancer drug response data. Comparisons are also made to alternative Bayesian and frequentist penalised-likelihood methods for incorporating network-based information. Conclusions: The empirical Bayes method proposed here can aid prior elicitation for Bayesian variable selection studies and help to guard against mis-specification of priors. Empirical Bayes, together with the proposed pathway-based priors, results in an approach with a competitive variable selection performance. In addition, the overall procedure is fast, deterministic, and has very few user-set parameters, yet is capable of capturing interplay between molecular players. The approach presented is general and readily applicable in any setting with multiple sources of biological prior knowledge

Fitness Cost Associated With Enhanced EPSPS Gene Copy Number and Glyphosate Resistance in an Amaranthus tuberculatus Population

The evolution of resistance to pesticides in agricultural systems provides an opportunity to study the fitness costs and benefits of novel adaptive traits. Here, we studied a population of Amaranthus tuberculatus (common waterhemp), which has evolved resistance to glyphosate. The growth and fitness of seed families with contrasting levels of glyphosate resistance was assessed in the absence of glyphosate to determine their ability to compete for resources under intra- and interspecific competition. We identified a positive correlation between the level of glyphosate resistance and gene copy number for the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) glyphosate target, thus identifying gene amplification as the mechanism of resistance within the population. Resistant A. tuberculatus plants were found to have a lower competitive response when compared to the susceptible phenotypes with 2.76 glyphosate resistant plants being required to have an equal competitive effect as a single susceptible plant. A growth trade-off was associated with the gene amplification mechanism under intra-phenotypic competition where 20 extra gene copies were associated with a 26.5 % reduction in dry biomass. Interestingly, this growth trade-off was mitigated when assessed under interspecific competition from maize

A Unique Population of Ventral Tegmental Area Neurons Inhibits the Lateral Habenula to Promote Reward

Lateral habenula (LHb) neurons convey aversive and negative reward conditions through potent indirect inhibition of ventral tegmental area (VTA) dopaminergic neurons. While VTA dopaminergic neurons reciprocally project to the LHb, the electrophysiological properties and the behavioral consequences associated with selective manipulations of this circuit are unknown. Here, we identify a novel inhibitory input to the LHb arising from a unique population of VTA neurons expressing dopaminergic markers. Optogenetic activation of this circuit resulted in no detectable dopamine release in LHb brain slices. Instead, stimulation produced GABA-mediated inhibitory synaptic transmission, which suppressed the firing of postsynaptic LHb neurons in brain slices and increased the spontaneous firing rate of VTA dopaminergic neurons in vivo. Furthermore, in vivo activation of this pathway produced reward-related phenotypes that were dependent on intra-LHb GABAA receptor signaling. These results suggest that non-canonical inhibitory signaling by these hybrid dopaminergic-GABAergic neurons act to suppress LHb output under rewarding conditions

A systems analysis of the chemosensitivity of breast cancer cells to the polyamine analogue PG-11047

<p>Abstract</p> <p>Background</p> <p>Polyamines regulate important cellular functions and polyamine dysregulation frequently occurs in cancer. The objective of this study was to use a systems approach to study the relative effects of PG-11047, a polyamine analogue, across breast cancer cells derived from different patients and to identify genetic markers associated with differential cytotoxicity.</p> <p>Methods</p> <p>A panel of 48 breast cell lines that mirror many transcriptional and genomic features present in primary human breast tumours were used to study the antiproliferative activity of PG-11047. Sensitive cell lines were further examined for cell cycle distribution and apoptotic response. Cell line responses, quantified by the GI₅₀(dose required for 50% relative growth inhibition) were correlated with the omic profiles of the cell lines to identify markers that predict response and cellular functions associated with drug sensitivity.</p> <p>Results</p> <p>The concentrations of PG-11047 needed to inhibit growth of members of the panel of breast cell lines varied over a wide range, with basal-like cell lines being inhibited at lower concentrations than the luminal cell lines. Sensitive cell lines showed a significant decrease in S phase fraction at doses that produced little apoptosis. Correlation of the GI₅₀values with the omic profiles of the cell lines identified genomic, transcriptional and proteomic variables associated with response.</p> <p>Conclusions</p> <p>A 13-gene transcriptional marker set was developed as a predictor of response to PG-11047 that warrants clinical evaluation. Analyses of the pathways, networks and genes associated with response to PG-11047 suggest that response may be influenced by interferon signalling and differential inhibition of aspects of motility and epithelial to mesenchymal transition.</p> <p>See the related commentary by Benes and Settleman: <url>http://www.biomedcentral.com/1741-7015/7/78</url></p

Functional Analysis of a Breast Cancer-Associated FGFR2 Single Nucleotide Polymorphism Using Zinc Finger Mediated Genome Editing

The authors thank Barts Charity for funding and Breast Cancer Campaign for funding the Barts Breast Tissue Bank

Suppression of apoptosis inhibitor c-FLIP selectively eliminates breast cancer stem cell activity in response to the anti-cancer agent, TRAIL

Introduction It is postulated that breast cancer stem cells (bCSCs) mediate disease recurrence and drive formation of distant metastases - the principal cause of mortality in breast cancer patients. Therapeutic targeting of bCSCs however, is hampered by their heterogeneity and resistance to existing therapeutics. In order to identify strategies to selectively remove bCSCs from breast cancers, irrespective of their clinical subtype, we sought an apoptosis mechanism that would target bCSCs yet would not kill normal cells. Suppression of the apoptosis inhibitor cellular FLICE-Like Inhibitory Protein (c-FLIP) partially sensitizes breast cancer cells to the anti-cancer agent Tumour Necrosis Factor-Related Apoptosis Inducing Ligand (TRAIL). Here we demonstrate in breast cancer cell lines that bCSCs are exquisitely sensitive to the de-repression of this pro-apoptotic pathway, resulting in a dramatic reduction in experimental metastases and the loss of bCSC self-renewal. Methods Suppression c-FLIP was performed by siRNA (FLIPi) in four breast cancer cell lines and by conditional gene-knockout in murine mammary glands. Sensitivity of these cells to TRAIL was determined by complementary cell apoptosis assays, including a novel heterotypic cell assay, while tumour-initiating potential of cancer stem cell subpopulations was determined by mammosphere cultures, aldefluor assay and in vivo transplantation. Results Genetic suppression of c-FLIP resulted in the partial sensitization of TRAIL-resistant cancer lines to the pro-apoptotic effects of TRAIL, irrespective of their cellular phenotype, yet normal mammary epithelial cells remained refractory to killing. While 10%-30% of the cancer cell populations remained viable after TRAIL/FLIPi treatment, subsequent mammosphere and aldefluor assays demonstrated that this pro-apoptotic stimulus selectively targeted the functional bCSC pool, eliminating stem cell renewal. This culminated in an 80% reduction in primary tumours and a 98% reduction in metastases following transplantation. The recurrence of residual tumour initiating capacity was consistent with the observation that post-treated adherent cultures re-acquired bCSC-like properties in vitro. Importantly however this recurrent bCSC activity was attenuated following repeated TRAIL/FLIPi treatment. Conclusions We describe an apoptotic mechanism that selectively and repeatedly removes bCSC activity from breast cancer cell lines and suggest that a combined TRAIL/FLIPi therapy could prevent metastatic disease progression in a broad range of breast cancer subtypes. [PROVISIONAL

Representativeness of microsatellite distributions in genomes, as revealed by 454 GS-FLX Titanium pyrosequencing

<p>Abstract</p> <p>Background</p> <p>Microsatellites are markers of choice in population genetics and genomics, as they provide useful insight into patterns and processes as diverse as genome evolutionary dynamics and demographic processes. The acquisition of microsatellites through multiplex-enriched libraries and 454 GS-FLX Titanium pyrosequencing is a promising new tool for the isolation of new markers in unknown genomes. This approach can also be used to evaluate the extent to which microsatellite-enriched libraries are representative of the genome from which they were isolated. In this study, we deciphered potential discrepancies in microsatellite content recovery for two reference genomes (<it>Apis mellifera </it>and <it>Danio rerio</it>), selected on the basis of their extreme heterogeneity in terms of the proportions and distributions of microsatellites on chromosomes.</p> <p>Results</p> <p>The <it>A. mellifera </it>genome, in particular, was found to be highly heterogeneous, due to extremely high rates of recombination, with hotspots, but the only bias consistently introduced into pyrosequenced multiplex-enriched libraries concerned sequence length, with the overrepresentation of sequences 160 to 320 bp in length. Other deviations from expected proportions or distributions of motifs on chromosomes were observed, but the significance and intensity of these deviations was mostly limited. Furthermore, no consistent adverse competition between multiplexed probes was observed during the motif enrichment phase.</p> <p>Conclusions</p> <p>This approach therefore appears to be a promising strategy for improving the development of microsatellites, as it introduces no major bias in terms of the proportions and distribution of microsatellites.</p

PD 0332991, a selective cyclin D kinase 4/6 inhibitor, preferentially inhibits proliferation of luminal estrogen receptor-positive human breast cancer cell lines in vitro

Abstract Introduction Alterations in cell cycle regulators have been implicated in human malignancies including breast cancer. PD 0332991 is an orally active, highly selective inhibitor of the cyclin D kinases (CDK)4 and CDK6 with ability to block retinoblastoma (Rb) phosphorylation in the low nanomolar range. To identify predictors of response, we determined the in vitro sensitivity to PD 0332991 across a panel of molecularly characterized human breast cancer cell lines. Methods Forty-seven human breast cancer and immortalized cell lines representing the known molecular subgroups of breast cancer were treated with PD 0332991 to determine IC50 values. These data were analyzed against baseline gene expression data to identify genes associated with PD 0332991 response. Results Cell lines representing luminal estrogen receptor-positive (ER+) subtype (including those that are HER2 amplified) were most sensitive to growth inhibition by PD 0332991 while nonluminal/basal subtypes were most resistant. Analysis of variance identified 450 differentially expressed genes between sensitive and resistant cells. pRb and cyclin D1 were elevated and CDKN2A (p16) was decreased in the most sensitive lines. Cell cycle analysis showed G0/G1 arrest in sensitive cell lines and Western blot analysis demonstrated that Rb phosphorylation is blocked in sensitive lines but not resistant lines. PD 0332991 was synergistic with tamoxifen and trastuzumab in ER+ and HER2-amplified cell lines, respectively. PD 0332991 enhanced sensitivity to tamoxifen in cell lines with conditioned resistance to ER blockade. Conclusions These studies suggest a role for CDK4/6 inhibition in some breast cancers and identify criteria for patient selection in clinical studies of PD 0332991

The Effect of Selenium Supplementation in the Prevention of DNA Damage in White Blood Cells of Hemodialyzed Patients: A Pilot Study

Patients with chronic kidney disease (CKD) have an increased incidence of cancer. It is well known that long periods of hemodialysis (HD) treatment are linked to DNA damage due to oxidative stress. In this study, we examined the effect of selenium (Se) supplementation to CKD patients on HD on the prevention of oxidative DNA damage in white blood cells. Blood samples were drawn from 42 CKD patients on HD (at the beginning of the study and after 1 and 3 months) and from 30 healthy controls. Twenty-two patients were supplemented with 200 μg Se (as Se-rich yeast) per day and 20 with placebo (baker's yeast) for 3 months. Se concentration in plasma and DNA damage in white blood cells expressed as the tail moment, including single-strand breaks (SSB) and oxidative bases lesion in DNA, using formamidopyrimidine glycosylase (FPG), were measured. Se concentration in patients was significantly lower than in healthy subjects (P < 0.0001) and increased significantly after 3 months of Se supplementation (P < 0.0001). Tail moment (SSB) in patients before the study was three times higher than in healthy subjects (P < 0.01). After 3 months of Se supplementation, it decreased significantly (P < 0.01) and was about 16% lower than in healthy subjects. The oxidative bases lesion in DNA (tail moment, FPG) of HD patients at the beginning of the study was significantly higher (P < 0.01) compared with controls, and 3 months after Se supplementation it was 2.6 times lower than in controls (P < 0.01). No changes in tail moment was observed in the placebo group. In conclusion, our study shows that in CKD patients on HD, DNA damage in white blood cells is higher than in healthy controls, and Se supplementation prevents the damage of DNA