Search

Motivation Meticulous selection of chromatographic peak detection parameters and algorithms is a crucial step in preprocessing liquid chromatography–mass spectrometry (LC-MS) data. However, as mass-to-charge ratio and retention time shifts are larger between batches than within batches, finding apt parameters for all samples of a large-scale multi-batch experiment with the aim of minimizing information loss becomes a challenging task. Preprocessing independent batches individually can curtail said problems but requires a method for aligning and combining them for further downstream analysis. Results We present two methods for aligning and combining individually preprocessed batches in multi-batch LC-MS experiments. Our developed methods were tested on six sets of simulated and six sets of real datasets. Furthermore, by estimating the probabilities of peak insertion, deletion and swap between batches in authentic datasets, we demonstrate that retention order swaps are not rare in untargeted LC-MS data. Availability and implementation kmersAlignment and rtcorrectedAlignment algorithms are made available as an R package with raw data at https://metabocombiner.img.cas.cz Supplementary information Supplementary data are available at Bioinformatics online.

Stress responses are activated by the hypothalamic-pituitary-adrenal axis (HPA axis), culminating in the release of glucocorticoids. During prolonged periods of secretion of glucocorticoids or inappropriate behavioral responses to a stressor, pathologic conditions may occur. Increased glucocorticoid concentration is linked to generalized anxiety, and there are knowledge gaps regarding its regulation. It is known that the HPA axis is under GABAergic control, but the contribution of the individual subunits of the GABA receptor is largely unknown. In this study, we investigated the relationship between the α5 subunit and corticosterone levels in a new mouse model deficient for Gabra5, which is known to be linked to anxiety disorders in humans and phenologs observed in mice. We observed decreased rearing behavior, suggesting lower anxiety in the Gabra5−/− animals; however, such a phenotype was absent in the open field and elevated plus maze tests. In addition to decreased rearing behavior, we also found decreased levels of fecal corticosterone metabolites in Gabra5−/− mice indicating a lowered stress response. Moreover, based on the electrophysiological recordings where we observed a hyperpolarized state of hippocampal neurons, we hypothesize that the constitutive ablation of the Gabra5 gene leads to functional compensation with other channels or GABA receptor subunits in this model.

International audience; Cannabinoid receptor 1 (CB1R), a G protein-coupled receptor, plays a fundamental role in synaptic plasticity. Abnormal activity and deregulation of CB1R signaling result in a broad spectrum of pathological conditions. CB1R signaling is regulated by receptor desensitization including phosphorylation of residues within the intracellular C terminus by G protein-coupled receptor kinases (GRKs) that may lead to endocytosis. Furthermore, CB1R signaling is regulated by the protein Src homology 3-domain growth factor receptor-bound 2-like (SGIP1) that hinders receptor internalization, while enhancing CB1R association with β-arrestin. It has been postulated that phosphorylation of two clusters of serine/threonine residues, 425 SMGDS429 and 460 TMSVSTDTS468 , within the CB1R C-tail controls dynamics of the association between receptor and its interaction partners involved in desensitization. Several molecular determinants of these events are still not well understood. We hypothesized that the dynamics of these interactions are modulated by SGIP1. Using a panel of CB1Rs mutated in the aforementioned serine and threonine residues, together with an array of Bioluminescence energy transfer-based (BRET) sensors, we discovered that GRK3 forms complexes with Gβγ subunits of G proteins that largely independent of GRK3's interaction with CB1R. Furthermore, CB1R interacts only with activated GRK3. Interestingly, phosphorylation of two specific residues on CB1R triggers GRK3 dissociation from the desensitized receptor. SGIP1 increases the association of GRK3 with Gβγ subunits of G proteins, and with CB1R. Altogether, our data suggest that the CB1R signalosome complex is dynamically controlled by sequential phosphorylation of the receptor C-tail and is also modified by SGIP1.

Introduction: Identification of factors involved in early stages of left ventricular systolic function recovery after acute myocardial infarction (AMI) can provide a basis for novel therapies. We aimed to characterize proteomic biomarkers and underlying mechanistic pathways associated with improvement of left ventricular systolic function. Hypothesis: Proteomic biomarkers and protein-protein interaction pathways differ between patients with and without systolic function recovery after AMI. Methods: Patients with ejection fraction (EF) 50% at the follow-up visit 3-6 months after discharge (EF recovery group). Controls were patients with EF Results: In total, 41 cases and 41 controls (mean age 59±10 years, 73% males) were included. At baseline, study groups did not differ in the main clinical variables. Of the 362 biomarkers tested, and after adjustment for age, gender and baseline EF, 44 proteins were differentially expressed, of which 2 (interleukin 6 and stem cell factor) showed association with EF recovery after false discovery rate correction. In the Reactome pathway analysis, Interleukin-6 family signaling was the most significant pathway associated with systolic function recovery. Conclusions: We identified the Interleukin-6 family signaling pathway with IL-6 and leukemia inhibitory factor (LIF) as potential early factors supporting recovery of systolic function after AMI.

Maintenance of genome stability is essential for every living cell as genetic information is repeatedly challenged during DNA replication in each cell division event. Errors, defects, delays, and mistakes that arise during mitosis or meiosis lead to an activation of DNA repair processes and in case of their failure, programmed cell death, i.e. apoptosis, could be initiated. Fam208a is a protein whose importance in heterochromatin maintenance has been described recently. In this work, we describe the crucial role of Fam208a in sustaining genome stability during cellular division. The targeted depletion of Fam208a in mice using CRISPR/Cas9 led to embryonic lethality before E12.5. We also used the siRNA approach to downregulate Fam208a in zygotes to avoid the influence of maternal RNA in the early stages of development. This early downregulation increased arresting of the embryonal development at the two-cell stage and the occurrence of multipolar spindles formation. To investigate this further, we used the yeast two-hybrid (Y2H) system and identified new putative interaction partners Gpsm2, Svil, and Itgb3bp. Their co-expression with Fam208a was assessed by RT-qPCR profiling and in situ hybridization [1] in multiple murine tissues. Based on these results we proposed that Fam208a functions within the HUSH complex by interaction with Mphosph8 as these proteins are not only able to physically interact but also co-localise. We are bringing new evidence that Fam208a is a multi-interacting protein affecting genome stability on the cell division level at the earliest stages of development and by interaction with methylation complex in adult tissues. In addition to its epigenetic functions, Fam208a appears to have an important role in the zygotic division, possibly via interaction with newly identified putative partners Gpsm2, Svil, and Itgb3bp.

Type 2 diabetes mellitus represents a major health problem with increasing prevalence worldwide. Limited efficacy of current therapies has prompted a search for novel therapeutic options. Here we show that treatment of pre-diabetic mice with mitochondrially targeted tamoxifen, a potential anti-cancer agent with senolytic activity, improves glucose tolerance and reduces body weight with most pronounced reduction of visceral adipose tissue due to reduced food intake, suppressed adipogenesis and elimination of senescent cells. Glucose-lowering effect of mitochondrially targeted tamoxifen is linked to improvement of type 2 diabetes mellitus-related hormones profile and is accompanied by reduced lipid accumulation in liver. Lower senescent cell burden in various tissues, as well as its inhibitory effect on pre-adipocyte differentiation, results in lower level of circulating inflammatory mediators that typically enhance metabolic dysfunction. Targeting senescence with mitochodrially targeted tamoxifen thus represents an approach to the treatment of type 2 diabetes mellitus and its related comorbidities, promising a complex impact on senescence-related pathologies in aging population of patients with type 2 diabetes mellitus with potential translation into the clinic.

MicroRNA (miRNA) profiling represents a promising source of cancer-related biomarkers. miRNA signatures are specific for each cancer type and subgroups of patients with diverse treatment sensitivity. Yet this miRNA potential has not been satisfactorily explored in rectal cancer (RC). The aim of the study was to identify the specific miRNA signature with clinical and therapeutic relevance for RC. Expressions of 2555 miRNA were examined in 20 pairs of rectal tumors and matched non-malignant tissues by 3D-Gene Toray microarray. Candidate miRNAs were validated in an independent cohort of 100 paired rectal tissues and in whole plasma and exosomes of 100 RC patients. To study the association of miRNA profile with therapeutic outcomes, plasma samples were taken repeatedly over a time period of 1 year reflecting thus patients' treatment responses. Finally, the most prominent miRNAs were investigated in vitro for their involvement in cell growth. We identified RC-specific miRNA signature that distinguishes responders from non-responders to adjuvant chemotherapy. A predominant part of identified miRNAs was represented by the members of miR-17/92 cluster. Upregulation of miRNA-17, -18a, -18b, -19a, -19b, -20a, -20b and -106a in tumor was associated with higher risk of tumor relapse and their overexpression in RC cell lines stimulated cellular proliferation. Examination of these miRNAs in plasma exosomes showed that their levels differed between RC patients and healthy controls and correlated with patient's treatment response. miRNAs from miR-17/92 cluster represent a non-invasive biomarker to predict posttreatment prognosis in RC patients.

Resistance to systemic therapy is a major problem in metastatic breast cancer (MBC) that can be explained by initial tumor heterogeneity as well as by evolutionary changes during therapy and tumor progression. Circulating tumor cells (CTCs) detected in a liquid biopsy can be sampled and characterized repeatedly during therapy in order to monitor treatment response and disease progression. Our aim was to investigate how CTC derived gene expression of treatment predictive markers (ESR1/HER2) and other cancer associated markers changed in patient blood samples during six months of first-line systemic treatment for MBC. CTCs from 36 patients were enriched using CellSearch (Janssen Diagnostics) and AdnaTest (QIAGEN) before gene expression analysis was performed with a customized gene panel (TATAA Biocenter). Our results show that antibodies against HER2 and EGFR were valuable to isolate CTCs unidentified by CellSearch and possibly lacking EpCAM expression. Evaluation of patients with clinically different breast cancer subgroups demonstrated that gene expression of treatment predictive markers changed over time. This change was especially prominent for HER2 expression. In conclusion, we found that changed gene expression during first-line systemic therapy for MBC could be a possible explanation for treatment resistance. Characterization of CTCs at several time-points during therapy could be informative for treatment selection.

Zinc finger 644 (Zfp644 in mouse, ZNF644 in human) gene is a transcription factor whose mutation S672G is considered a potential genetic factor of inherited high myopia. ZNF644 interacts with G9a/GLP complex, which functions as a H3K9 methyltransferase to silence transcription. In this study, we generated mouse models to unravel the mechanisms leading to symptoms associated with high myopia. Employing TALEN technology, two mice mutants were generated, either with the disease-carrying mutation (Zfp644S673G) or with a truncated form of Zfp644 (Zfp644Δ8). Eye morphology and visual functions were analysed in both mutants, revealing a significant difference in a vitreous chamber depth and lens diameter, however the physiological function of retina was preserved as found under the high-myopia conditions. Our findings prove that ZNF644/Zfp644 is involved in the development of high-myopia, indicating that mutations such as, Zfp644S673G and Zfp644Δ8 are causative for changes connected with the disease. The developed models represent a valuable tool to investigate the molecular basis of myopia pathogenesis and its potential treatment. Electronic supplementary material The online version of this article (10.1186/s13578-019-0280-4) contains supplementary material, which is available to authorized users.

Efficient unbiased data analysis is a major challenge for laboratories handling large cytometry datasets. We present EmbedSOM, a non-linear embedding algorithm based on FlowSOM that improves the analyses by providing high-performance visualization of complex single cell distributions within cellular populations and their transition states. The algorithm is designed for linear scaling and speed suitable for interactive analyses of millions of cells without downsampling. At the same time, the visualization quality is competitive with current state-of-art algorithms. We demonstrate the properties of EmbedSOM on workflows that improve two essential types of analyses: The native ability of EmbedSOM to align population positions in embedding is used for comparative analysis of multi-sample data, and the connection to FlowSOM is exploited for simplifying the supervised hierarchical dissection of cell populations. Additionally, we discuss the visualization of the trajectories between cellular states facilitated by the local linearity of the embedding.

Rhabdomyosarcomas (RMS) are a heterogeneous group of mesodermal tumors, the most common sub-types are embryonal (eRMS) and alveolar (aRMS) rhabdomyosarcoma. Immunohistochemical analysis revealed c-Myb expression in both eRMS and aRMS. c-Myb has been reported to be often associated with malignant human cancers. We therefore investigated the c-Myb role in RMS using cellular models of RMS. Specific suppression of c-Myb by a lentiviral vector expressing doxycycline (Dox)-inducible c-Myb shRNA inhibited proliferation, colony formation, and migration of the eRMS cell line (RD), but not of the aRMS cell line (RH30). Upon c-Myb knockdown in eRMS cells, cells accumulated in G0/G1 phase, the invasive behaviour of cells was repressed, and elevated levels of myosin heavy chain, marker of muscle differentiation, was detected. Next, we used an RD-based xenograft model to investigate the role of c-Myb in eRMS tumorigenesis in vivo. We found that Dox administration did not result in efficient suppression of c-Myb in growing tumors. However, when c-Myb-deficient RD cells were implanted into SCID mice, we observed inefficient tumor grafting and attenuation of tumor growth during the initial stages of tumor expansion. The presented study suggests that c-Myb could be a therapeutic target in embryonal rhabdomyosarcoma assuming that its expression is ablated.

Cortical glial cells contain both ionotropic and metabotropic glutamate receptors. Despite several efforts, a comprehensive analysis of the entire family of glutamate receptors and their subunits present in glial cells is still missing. Here, we provide an overall picture of the gene expression of ionotropic (AMPA, kainate, NMDA) and the main metabotropic glutamate receptors in cortical glial cells isolated from GFAP/EGFP mice before and after focal cerebral ischemia. Employing single-cell RT-qPCR, we detected the expression of genes encoding subunits of glutamate receptors in GFAP/EGFP-positive (GFAP/EGFP(+)) glial cells in the cortex of young adult mice. Most of the analyzed cells expressed mRNA for glutamate receptor subunits, the expression of which, in most cases, even increased after ischemic injury. Data analyses disclosed several classes of GFAP/EGFP(+) glial cells with respect to glutamate receptors and revealed in what manner their expression correlates with the expression of glial markers prior to and after ischemia. Furthermore, we also examined the protein expression and functional significance of NMDA receptors in glial cells. Immunohistochemical analyses of all seven NMDA receptor subunits provided direct evidence that the GluN3A subunit is present in GFAP/EGFP(+) glial cells and that its expression is increased after ischemia. In situ and in vitro Ca(2+) imaging revealed that Ca(2+) elevations evoked by the application of NMDA were diminished in GFAP/EGFP(+) glial cells following ischemia. Our results provide a comprehensive description of glutamate receptors in cortical GFAP/EGFP(+) glial cells and may serve as a basis for further research on glial cell physiology and pathophysiology.

We have examined the imprecision in the estimation of PCR efficiency by means of standard curves based on strategic experimental design with large number of technical replicates. In particular, how robust this estimation is in terms of a commonly varying factors: the instrument used, the number of technical replicates performed and the effect of the volume transferred throughout the dilution series. We used six different qPCR instruments, we performed 1–16 qPCR replicates per concentration and we tested 2–10μl volume of analyte transferred, respectively. We find that the estimated PCR efficiency varies significantly across different instruments. Using a Monte Carlo approach, we find the uncertainty in the PCR efficiency estimation may be as large as 42.5% (95% CI) if standard curve with only one qPCR replicate is used in 16 different plates. Based on our investigation we propose recommendations for the precise estimation of PCR efficiency: (1) one robust standard curve with at least 3–4 qPCR replicates at each concentration shall be generated, (2) the efficiency is instrument dependent, but reproducibly stable on one platform, and (3) using a larger volume when constructing serial dilution series reduces sampling error and enables calibration across a wider dynamic range.

The interplay between dietary habits and individual genetic make-up is assumed to influence risk of cancer, via modulation of DNA integrity. Our aim was to characterize internal and external factors that underlie inter-individual variability in DNA damage and repair and to identify dietary habits beneficial for maintaining DNA integrity. Habitual diet was estimated in 340 healthy individuals using a food frequency questionnaire and biomarkers of antioxidant status were quantified in fasting blood samples. Markers of DNA integrity were represented by DNA strand breaks, oxidized purines, oxidized pyrimidines and a sum of all three as total DNA damage. DNA repair was characterized by genetic variants and functional activities of base and nucleotide excision repair pathways. Sex, fruit-based food consumption and XPG genotype were factors significantly associated with the level of DNA damage. DNA damage was higher in women (p=0.035). Fruit consumption was negatively associated with the number of all measured DNA lesions, and this effect was mediated mostly by β-cryptoxanthin and β-tocopherol (p

Maintenance of genome stability is essential for every living cell as genetic information is repeatedly challenged during DNA replication in each cell division event. Errors, defects, delays, and mistakes that arise during mitosis or meiosis lead to an activation of DNA repair processes and in case of their failure, programmed cell death, i.e. apoptosis, could be initiated. Fam208a is a protein whose importance in heterochromatin maintenance has been described recently. In this work, we describe the crucial role of Fam208a in sustaining the genome stability during the cellular division. The targeted depletion of Fam208a in mice using CRISPR/Cas9 leads to embryonic lethality before E12.5. We also used the siRNA approach to downregulate Fam208a in zygotes to avoid the influence of maternal RNA in the early stages of development. This early downregulation increased arresting of the embryonal development at the two-cell stage and occurrence of multipolar spindles formation. To investigate this further, we used the yeast two-hybrid (Y2H) system and identified new putative interaction partners Gpsm2, Amn1, Eml1, Svil, and Itgb3bp. Their co-expression with Fam208a was assessed by qRT-PCR profiling and in situ hybridisation [1] in multiple murine tissues. Based on these results we proposed that Fam208a functions within the HUSH complex by interaction with Mphosph8 as these proteins are not only able to physically interact but also co-localise. We are bringing new evidence that Fam208a is multi-interacting protein affecting genome stability on the level of cell division at the earliest stages of development and also by interaction with methylation complex in adult tissues. In addition to its epigenetic functions, Fam208a appears to have an additional role in zygotic division, possibly via interaction with newly identified putative partners Gpsm2, Amn1, Eml1, Svil, and Itgb3bp.

Successful molecular analyses of human solid tissues require intact biological material with well-preserved nucleic acids, proteins, and other cell structures. Pre-analytical handling, comprising of the collection of material at the operating theatre, is among the first critical steps that influence sample quality. The aim of this study was to compare the experimental outcomes obtained from samples collected and stored by the conventional means of snap freezing and by PAXgene Tissue System (Qiagen). These approaches were evaluated by measuring rRNA and mRNA integrity of the samples (RNA Quality Indicator and Differential Amplification Method) and by gene expression profiling. The collection procedures of the biological material were implemented in two hospitals during colon cancer surgery in order to identify the impact of the collection method on the experimental outcome. Our study shows that the pre-analytical sample handling has a significant effect on the quality of RNA and on the variability of qPCR data. PAXgene collection mode proved to be more easily implemented in the operating room and moreover the quality of RNA obtained from human colon tissues by this method is superior to the one obtained by snap freezing.

BACKGROUND Transcriptome analysis of circulating tumor cells (CTCs) holds great promise to unravel the biology of cancer cell dissemination and identify expressed genes and signaling pathways relevant to therapeutic interventions. METHODS CTCs were enriched based on their EpCAM expression (CellSearch®) or by size and deformability (ParsortixTM), identified by EpCAM and/or pan-keratin–specific antibodies, and isolated for single cell multiplex RNA profiling. RESULTS Distinct breast and prostate CTC expression signatures could be discriminated from RNA profiles of leukocytes. Some CTCs positive for epithelial transcripts (EpCAM and KRT19) also coexpressed leukocyte/mesenchymal associated markers (PTPRC and VIM). Additional subsets of CTCs within individual patients were characterized by divergent expression of genes involved in epithelial–mesenchymal transition (e.g., CDH2, MMPs, VIM, or ZEB1 and 2), DNA repair (RAD51), resistance to cancer therapy (e.g., AR, AR-V7, ERBB2, EGFR), cancer stemness (e.g., CD24 and CD44), activated signaling pathways involved in tumor progression (e.g., PIK3CA and MTOR) or cross talks between tumors and immune cells (e.g., CCL4, CXCL2, CXCL9, IL15, IL1B, or IL8). CONCLUSIONS Multimarker RNA profiling of single CTCs reveals distinct CTC subsets and provides important insights into gene regulatory networks relevant for cancer progression and therapy.

Background With the introduction of the first high-throughput qPCR instrument on the market it became possible to perform thousands of reactions in a single run compared to the previous hundreds. In the high-throughput reaction, only limited volumes of highly concentrated cDNA or DNA samples can be added. This necessity can be solved by pre-amplification, which became a part of the high-throughput experimental workflow. Here, we focused our attention on the limits of the specific target pre-amplification reaction and propose the optimal, general setup for gene expression experiment using BioMark instrument (Fluidigm). Results For evaluating different pre-amplification factors following conditions were combined: four human blood samples from healthy donors and five transcripts having high to low expression levels; each cDNA sample was pre-amplified at four cycles (15, 18, 21, and 24) and five concentrations (equivalent to 0.078 ng, 0.32 ng, 1.25 ng, 5 ng, and 20 ng of total RNA). Factors identified as critical for a success of cDNA pre-amplification were cycle of pre-amplification, total RNA concentration, and type of gene. The selected pre-amplification reactions were further tested for optimal Cq distribution in a BioMark Array. The following concentrations combined with pre-amplification cycles were optimal for good quality samples: 20 ng of total RNA with 15 cycles of pre-amplification, 20x and 40x diluted; and 5 ng and 20 ng of total RNA with 18 cycles of pre-amplification, both 20x and 40x diluted. Conclusions We set up upper limits for the bulk gene expression experiment using gene expression Dynamic Array and provided an easy-to-obtain tool for measuring of pre-amplification success. We also showed that variability of the pre-amplification, introduced into the experimental workflow of reverse transcription-qPCR, is lower than variability caused by the reverse transcription step. Electronic supplementary material The online version of this article (doi:10.1186/s12867-015-0033-9) contains supplementary material, which is available to authorized users.

Variations in the TP53 gene have been suggested to play a role in many cancers, including breast. We previously observed an association between TP53 haplotypes based on four polymorphisms (rs17878362, rs1042522, rs12947788, and rs17884306) and the risk of colorectal and pancreatic cancer. Based on these results, in the present study, we have investigated the same polymorphisms and their haplotypes in 705 breast cancer cases and 611 healthy controls in relation to the disease risk, histopathological features of the tumor and clinical outcomes. In comparison to the most common haplotype A1-G-C-G, all the other identified haplotypes were globally associated with a significantly decreased breast cancer risk (P = 0.006). In particular, the A2-G-C-G haplotype was associated with a marked decreased risk of breast cancer when compared with the common haplotype (P = 0.0001). Moreover, rs1042522 in patients carrying the GC genotype and receiving only the anthracycline-based chemotherapy was associated with both overall and disease-free survival (recessive model for overall survival HR = 0.30 95% CI 0.11–0.80, P = 0.02 and for disease-free survival HR = 0.42 95% CI 0.21–0.84, P = 0.01). Present results suggest common genetic features in the susceptibility to breast and gastrointestinal cancers in respect to TP53 variations. In fact, similar haplotype distributions were observed for breast, colorectal, and pancreatic patients in associations with cancer risk. Rs1042522 polymorphism (even after applying the Dunn-Bonferroni correction for multiple testing) appears to be an independent prognostic marker in breast cancer patients.

Background: Circulating tumor cells (CTC) are discussed to be an ideal surrogate marker to monitor disease progression in metastatic breast cancer (MBC). Besides CTC characterization for targeted therapies, it would also be desirable to know where these cells derive from or to which organ site they are going to. Here we investigated whether it is possible to predict the origin of metastatic lesion based on the expression of 46 genes in CTC of MBC patients (pts). Materials and Methods: 2×5 ml blood of 45 MBC pts and 20 healthy controls was collected at the time of disease progression (T0) and at two consecutive clinical staging (T1 and T2) for the detection of CTC applying immunomagnetic enrichment using the AdnaTest EMT-2/Stem Cell Select (QIAGEN Hannover GmbH, Germany). Pts were grouped a) into overall responders (OR) and overall non-responders (ONR), thus responding or not responding at T1 and T2 and b) according to sites of metastasis. PCR assays, targeting 46 transcripts and reference markers were used for a workflow based on pre-amplification and high throughput profiling (each samples in duplicates) with the full set of markers including also ValidPrime to correct for genomic background and InterPlate Calibrator to even out variations between runs. Data were analyzed using GenEx (MultiD, Sweden) and SAS. qPCR as well as technical reads were normalized using several normalization strategies. Results: The multidrug resistant protein gene MRP1 was significantly differently expressed if OR and ONR groups were compared. In the following order of significance, VEGFR1, Keratin (KRT) 19, EGFR, MET1, ALDH, progesterone receptor (PR), UPA, Cathepsin D, KIT1 and Ki67 were differentially expressed in CTC of pts who had already developed liver metastasis as compared to pts without liver metastasis. Interestingly, a small group of pts, developing liver metastases in the course of disease, showed the estrogen receptor (ER), PR, HER2, mammaglobin, KRT19 on a significantly lower level as compared to the other pts. Similarly, once CTC were ER and PR positive, the probability of bone metastasis development decreased. Conclusion: Our preliminary results indicate that the development of metastatic lesions is associated with site-specific CTC. Thus, besides using CTC as a monitoring tool to guide therapy, they might also indicate the site of metastasis which will allow a more precise decision concerning treatment strategy. Citation Format: Maren Bredemeier, Philippos Edimiris, Pawel Mach, Mikael Kubista, Robert Sjoback, Marie Jindrichova, Eva Rohlova, Vendula Novosadova, Katarina Kolostova, Siegfried Hauch, Bahriye Aktas, Mitra Tewes, Rainer Kimmig, Sabine Kasimir-Bauer. Gene expression signatures in circulating tumor cells are prognostic for metastatic lesions in breast cancer patients and correlate with response to therapy. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 502.