Identification of an FXR-modulated liver-intestine hybrid state in iPSC-derived hepatocyte-like cells.

BACKGROUND & AIMS: Pluripotent stem cell (PSC)-derived hepatocyte-like cells (HLC) have enormous potential as a replacement for primary hepatocytes in drug screening, toxicology and cell replacement therapy, but their genome-wide expression patterns differ strongly from primary human hepatocytes (PHH). METHODS: We differentiated human induced pluripotent stem cells (hiPSC) via definitive endoderm to HLC and characterized the cells by single-cell and bulk RNA-seq, with complementary epigenetic analyses. We then compared HLC to PHH and publicly available data on human fetal hepatocytes (FH) ex vivo; we performed bioinformatics-guided interventions to improve HLC differentiation via lentiviral transduction of the nuclear receptor FXR and agonist exposure. RESULTS: Single-cell RNA-seq revealed that transcriptomes of individual HLC display a hybrid state, where hepatocyte-associated genes are expressed in concert with genes that are not expressed in PHH - mostly intestinal genes - within the same cell. Bulk-level overrepresentation analysis, as well as regulon analysis at the single-cell level, identified sets of regulatory factors discriminating HLC, FH, and PHH, hinting at a central role for the nuclear receptor FXR in the functional maturation of HLC. Combined FXR expression plus agonist exposure enhanced the expression of hepatocyte-associated genes and increased the ability of bile canalicular secretion as well as lipid droplet formation, thereby increasing HLCs' similarity to PHH. The undesired non-liver gene expression was reproducibly decreased, although only by a moderate degree. CONCLUSION: In contrast to physiological hepatocyte precursor cells and mature hepatocytes, HLC co-express liver and hybrid genes in the same cell. Targeted modification of the FXR gene regulatory network improves their differentiation by suppressing intestinal traits whilst inducing hepatocyte features. LAY SUMMARY: Generation of human hepatocytes from stem cells represents an active research field but its success is hampered by the fact that the stem cell-derived 'hepatocytes' still show major differences to hepatocytes obtained from a liver. Here, we identified an important reason for the difference, specifically that the stem cell-derived 'hepatocyte' represents a hybrid cell with features of hepatocytes and intestinal cells. We show that a specific protein (FXR) suppresses intestinal and induces liver features, thus bringing the stem cell-derived cells closer to hepatocytes derived from human livers.

LIPG-promoted lipid storage mediates adaptation to oxidative stress in breast cancer.

Endothelial lipase (LIPG) is a cell surface associated lipase that displays phospholipase A1 activity towards phosphatidylcholine present in high-density lipoproteins (HDL). LIPG was recently reported to be expressed in breast cancer and to support proliferation, tumourigenicity and metastasis. Here we show that severe oxidative stress leading to AMPK activation triggers LIPG upregulation, resulting in intracellular lipid droplet accumulation in breast cancer cells, which supports survival. Neutralizing oxidative stress abrogated LIPG upregulation and the concomitant lipid storage. In human breast cancer, high LIPG expression was observed in a limited subset of tumours and was significantly associated with shorter metastasis-free survival in node-negative, untreated patients. Moreover, expression of PLIN2 and TXNRD1 in these tumours indicated a link to lipid storage and oxidative stress. Altogether, our findings reveal a previously unrecognized role for LIPG in enabling oxidative stress-induced lipid droplet accumulation in tumour cells that protects against oxidative stress, and thus supports tumour progression.

Toxicogenomics directory of rat hepatotoxicants in vivo and in cultivated hepatocytes.

Transcriptomics is developing into an invaluable tool in toxicology. The aim of this study was, using a transcriptomics approach, to identify genes that respond similar to many different chemicals (including drugs and industrial compounds) in both rat liver in vivo and in cultivated hepatocytes. For this purpose, we analyzed Affymetrix microarray expression data from 162 compounds that were previously tested in a concentration-dependent manner in rat livers in vivo and in rat hepatocytes cultivated in sandwich culture. These data were obtained from the Japanese Toxicogenomics Project (TGP) and North Rhine-Westphalian (NRW) data sets, which represent 138 and 29 compounds, respectively, and have only 5 compounds in common between them. The in vitro gene expression data from the NRW data set were generated in the present study, while TGP is publicly available. For each of the data sets, the overlap between up- or down-regulated genes in vitro and in vivo was identified, and named in vitro-in vivo consensus genes. Interestingly, the in vivo-in vitro consensus genes overlapped to a remarkable extent between both data sets, and were 21-times (upregulated genes) or 12-times (down-regulated genes) enriched compared to random expectation. Finally, the genes in the TGP and NRW overlap were used to identify the upregulated genes with the highest compound coverage, resulting in a seven-gene set of Cyp1a1, Ugt2b1, Cdkn1a, Mdm2, Aldh1a1, Cyp4a3, and Ehhadh. This seven-gene set was then successfully tested with structural analogues of valproic acid that are not present in the TGP and NRW data sets. In conclusion, the seven-gene set identified in the present study responds similarly in vitro and in vivo to a wide range of different chemicals. Despite these promising results with the seven-gene set, transcriptomics with cultivated rat hepatocytes remains a challenge, because in general many genes are up- or downregulated by in vitro culture per se, respond differently to test compounds in vitro and in vivo, and/or show higher variability in the in vitro system compared to the corresponding in vivo data.

Metabolic profiling of ob/ob mouse fatty liver using HR-MAS 1H-NMR combined with gene expression analysis reveals alterations in betaine metabolism and the transsulfuration pathway.

Metabolic perturbations resulting from excessive hepatic fat accumulation are poorly understood. Thus, in this study, leptin-deficient ob/ob mice, a mouse model of fatty liver disease, were used to investigate metabolic alterations in more detail. Metabolites were quantified in intact liver tissues of ob/ob (n = 8) and control (n = 8) mice using high-resolution magic angle spinning (HR-MAS) 1H-NMR. In addition, after demonstrating that HR-MAS 1H-NMR does not affect RNA integrity, transcriptional changes were measured by quantitative real-time PCR on RNA extracted from the same specimens after HR-MAS 1H-NMR measurements. Importantly, the gene expression changes obtained agreed with those observed by Affymetrix microarray analysis performed on RNA isolated directly from fresh-frozen tissue. In total, 40 metabolites could be assigned in the spectra and subsequently quantified. Quantification of lactate was also possible after applying a lactate-editing pulse sequence that suppresses the lipid signal, which superimposes the lactate methyl resonance at 1.3 ppm. Significant differences were detected for creatinine, glutamate, glycine, glycolate, trimethylamine-N-oxide, dimethylglycine, ADP, AMP, betaine, phenylalanine, and uridine. Furthermore, alterations in one-carbon metabolism, supported by both metabolic and transcriptional changes, were observed. These included reduced demethylation of betaine to dimethylglycine and the reduced expression of genes coding for transsulfuration pathway enzymes, which appears to preserve methionine levels, but may limit glutathione synthesis. Overall, the combined approach is advantageous as it identifies changes not only at the single gene or metabolite level but also deregulated pathways, thus providing critical insight into changes accompanying fatty liver disease. Graphical abstract A Evaluation of RNA integrity before and after HR-MAS 1H-NMR of intact mouse liver tissue. B Metabolite concentrations and gene expression levels assessed in ob/ob (steatotic) and ob/+ (control) mice using HR-MAS 1H-NMR and qRT-PCR, respectively.

The lung-specific proteome defined by integration of transcriptomics and antibody-based profiling.

The combined action of multiple cell types is essential for the physiological function of the lung, and increased awareness of the molecular constituents characterizing each cell type is likely to advance the understanding of lung biology and disease. In the current study, we used genome-wide RNA sequencing of normal lung parenchyma and 26 additional tissue types, combined with antibody-based protein profiling, to localize the expression to specific cell types. Altogether, 221 genes were found to be elevated in the lung compared with their expression in other analyzed tissues. Among the gene products were several well-known markers, but also several proteins previously not described in the context of the lung. To link the lung-specific molecular repertoire to human disease, survival associations of pneumocyte-specific genes were assessed by using transcriptomics data from 7 non-small-cell lung cancer (NSCLC) cohorts. Transcript levels of 10 genes (SFTPB, SFTPC, SFTPD, SLC34A2, LAMP3, CACNA2D2, AGER, EMP2, NKX2-1, and NAPSA) were significantly associated with survival in the adenocarcinoma subgroup, thus qualifying as promising biomarker candidates. In summary, based on an integrated omics approach, we identified genes with elevated expression in lung and localized corresponding protein expression to different cell types. As biomarker candidates, these proteins may represent intriguing starting points for further exploration in health and disease.

Identification of sample annotation errors in gene expression datasets.

The comprehensive transcriptomic analysis of clinically annotated human tissue has found widespread use in oncology, cell biology, immunology, and toxicology. In cancer research, microarray-based gene expression profiling has successfully been applied to subclassify disease entities, predict therapy response, and identify cellular mechanisms. Public accessibility of raw data, together with corresponding information on clinicopathological parameters, offers the opportunity to reuse previously analyzed data and to gain statistical power by combining multiple datasets. However, results and conclusions obviously depend on the reliability of the available information. Here, we propose gene expression-based methods for identifying sample misannotations in public transcriptomic datasets. Sample mix-up can be detected by a classifier that differentiates between samples from male and female patients. Correlation analysis identifies multiple measurements of material from the same sample. The analysis of 45 datasets (including 4913 patients) revealed that erroneous sample annotation, affecting 40 % of the analyzed datasets, may be a more widespread phenomenon than previously thought. Removal of erroneously labelled samples may influence the results of the statistical evaluation in some datasets. Our methods may help to identify individual datasets that contain numerous discrepancies and could be routinely included into the statistical analysis of clinical gene expression data.

Glycerophospholipid profile in oncogene-induced senescence.

Alterations in lipid metabolism and in the lipid composition of cellular membranes are linked to the pathology of numerous diseases including cancer. However, the influence of oncogene expression on cellular lipid profile is currently unknown. In this work we analyzed changes in lipid profiles that are induced in the course of ERBB2-expression mediated premature senescence. As a model system we used MCF-7 breast cancer cells with doxycycline-inducible expression of NeuT, an oncogenic ERBB2 variant. Affymetrix gene array data showed NeuT-induced alterations in the transcription of many enzymes involved in lipid metabolism, several of which (ACSL3, CHPT1, PLD1, LIPG, MGLL, LDL and NPC1) could be confirmed by quantitative realtime PCR. A study of the glycerophospholipid and lyso-glycerophospholipid profiles, obtained by high performance liquid chromatography coupled to Fourier-transform ion cyclotron resonance-mass spectrometry revealed senescence-associated changes in numerous lipid species, including mitochondrial lipids. The most prominent changes were found in PG(34:1), PG(36:1) (increased) and LPE(18:1), PG(40:7) and PI(36:1) (decreased). Statistical analysis revealed a general trend towards shortened phospholipid acyl chains in senescence and a significant trend to more saturated acyl chains in the class of phosphatidylglycerol. Additionally, the cellular cholesterol content was elevated and accumulated in vacuoles in senescent cells. These changes were accompanied by increased membrane fluidity. In mitochondria, loss of membrane potential along with altered intracellular distribution was observed. In conclusion, we present a comprehensive overview of altered cholesterol and glycerophospholipid patterns in senescence, showing that predominantly mitochondrial lipids are affected and lipid species less susceptible to peroxidation are increased.

Proteome analysis of monocytes implicates altered mitochondrial biology in adults reporting adverse childhood experiences.

The experience of adversity in childhood has been associated with poor health outcomes in adulthood. In search of the biological mechanisms underlying these effects, research so far focused on alterations of DNA methylation or shifts in transcriptomic profiles. The level of protein, however, has been largely neglected. We utilized mass spectrometry to investigate the proteome of CD14+ monocytes in healthy adults reporting childhood adversity and a control group before and after psychosocial stress exposure. Particular proteins involved in (i) immune processes, such as neutrophil-related proteins, (ii) protein metabolism, or (iii) proteins related to mitochondrial biology, such as those involved in energy production processes, were upregulated in participants reporting exposure to adversity in childhood. This functional triad was further corroborated by protein interaction- and co-expression analyses, was independent of stress exposure, i.e. observed at both pre- and post-stress time points, and became evident especially in females. In line with the mitochondrial allostatic load model, our findings provide evidence for the long-term effects of childhood adversity on mitochondrial biology.

Inhibiting the glycerophosphodiesterase EDI3 in ER-HER2+ breast cancer cells resistant to HER2-targeted therapy reduces viability and tumour growth.

BACKGROUND: Intrinsic or acquired resistance to HER2-targeted therapy is often a problem when small molecule tyrosine kinase inhibitors or antibodies are used to treat patients with HER2 positive breast cancer. Therefore, the identification of new targets and therapies for this patient group is warranted. Activated choline metabolism, characterized by elevated levels of choline-containing compounds, has been previously reported in breast cancer. The glycerophosphodiesterase EDI3 (GPCPD1), which hydrolyses glycerophosphocholine to choline and glycerol-3-phosphate, directly influences choline and phospholipid metabolism, and has been linked to cancer-relevant phenotypes in vitro. While the importance of choline metabolism has been addressed in breast cancer, the role of EDI3 in this cancer type has not been explored. METHODS: EDI3 mRNA and protein expression in human breast cancer tissue were investigated using publicly-available Affymetrix gene expression microarray datasets (n = 540) and with immunohistochemistry on a tissue microarray (n = 265), respectively. A panel of breast cancer cell lines of different molecular subtypes were used to investigate expression and activity of EDI3 in vitro. To determine whether EDI3 expression is regulated by HER2 signalling, the effect of pharmacological inhibition and siRNA silencing of HER2, as well as the influence of inhibiting key components of signalling cascades downstream of HER2 were studied. Finally, the influence of silencing and pharmacologically inhibiting EDI3 on viability was investigated in vitro and on tumour growth in vivo. RESULTS: In the present study, we show that EDI3 expression is highest in ER-HER2 + human breast tumours, and both expression and activity were also highest in ER-HER2 + breast cancer cell lines. Silencing HER2 using siRNA, as well as inhibiting HER2 signalling with lapatinib decreased EDI3 expression. Pathways downstream of PI3K/Akt/mTOR and GSK3ß, and transcription factors, including HIF1α, CREB and STAT3 were identified as relevant in regulating EDI3 expression. Silencing EDI3 preferentially decreased cell viability in the ER-HER2 + cells. Furthermore, silencing or pharmacologically inhibiting EDI3 using dipyridamole in ER-HER2 + cells resistant to HER2-targeted therapy decreased cell viability in vitro and tumour growth in vivo. CONCLUSIONS: Our results indicate that EDI3 may be a potential novel therapeutic target in patients with HER2-targeted therapy-resistant ER-HER2 + breast cancer that should be further explored.

Expression of aurora kinase A is associated with metastasis-free survival in node-negative breast cancer patients.

BACKGROUND: Inhibitors targeting the cell cycle-regulated aurora kinase A (AURKA) are currently being developed. Here, we examine the prognostic impact of AURKA in node-negative breast cancer patients without adjuvant systemic therapy (n = 766). METHODS: AURKA was analyzed using microarray-based gene-expression data from three independent cohorts of node-negative breast cancer patients. In multivariate Cox analyses, the prognostic impact of age, histological grade, tumor size, estrogen receptor (ER), and HER2 were considered. RESULTS: Patients with higher AURKA expression had a shorter metastasis-free survival (MFS) in the Mainz (HR 1.93; 95% CI 1.34 - 2.78; P < 0.001), Rotterdam (HR 1.95; 95% CI 1.45- 2.63; P<0.001) and Transbig (HR 1.52; 95% CI 1.14-2.04; P=0.005) cohorts. AURKA was also associated with MFS in the molecular subtype ER+/HER2- carcinomas (HR 2.10; 95% CI 1.70-2.59; P<0.001), but not in ER-/HER2- nor in HER2+ carcinomas. In the multivariate Cox regression adjusted to age, grade and tumor size, AURKA showed independent prognostic significance in the ER+/HER2- subtype (HR 1.73; 95% CI 1.24-2.42; P=0.001). Prognosis of patients in the highest quartile of AURKA expression was particularly poor. In addition, AURKA correlated with the proliferation metagene (R=0.880; P<0.001), showed a positive association with grade (P<0.001), tumor size (P<0.001) and HER2 (P<0.001), and was inversely associated with ER status (P<0.001). CONCLUSIONS: AURKA is associated with worse prognosis in estrogen receptor positive breast carcinomas. Patients with the highest AURKA expression (>75% percentile) have a particularly bad prognosis and may profit from therapy with AURKA inhibitors.

Optimal strategies for sequential validation of significant features from high-dimensional genomic data.

High-dimensional genomic studies play a key role in identifying critical features that are significantly associated with a phenotypic outcome. The two most important examples are the detection of (1) differentially expressed genes from genome-wide gene expression studies and (2) single-nucleotide polymorphisms (SNPs) from genome-wide association studies. Such experiments are often associated with high noise levels, and the validity of statistical conclusions suffers from low sample size compared to large number of features. The corresponding multiple testing problem calls for the identification of optimal strategies for controlling the numbers of false discoveries and false nondiscoveries. In addition, a frequent validation problem is that features identified as important in one study are often less so in another study. Adjustment for multiple testing in both studies separately increases the risk of missing the crucial features even further. These problems can be addressed by sequential validation strategies, where only significant features identified in one study enter as candidates in the next study. The quality associated with different studies, for example, in terms of noise levels, may vary considerably. By performing simulation studies it is possible to demonstrate that the optimal order for this stepwise procedure is to sort experimental studies according to their quality in descending order. The impact of the method for multiple testing adjustment (Bonferroni-Holm, FDR) was also analyzed. Finally, the sequential validation strategy was applied to three large breast cancer studies with gene expression measurements, confirming the crucial impact of the order of the validation steps in a real-world application.

Comparison of scores for bimodality of gene expression distributions and genome-wide evaluation of the prognostic relevance of high-scoring genes.

BACKGROUND: A major goal of the analysis of high-dimensional RNA expression data from tumor tissue is to identify prognostic signatures for discriminating patient subgroups. For this purpose genome-wide identification of bimodally expressed genes from gene array data is relevant because distinguishability of high and low expression groups is easier compared to genes with unimodal expression distributions.Recently, several methods for the identification of genes with bimodal distributions have been introduced. A straightforward approach is to cluster the expression values and score the distance between the two distributions. Other scores directly measure properties of the distribution. The kurtosis, e.g., measures divergence from a normal distribution. An alternative is the outlier-sum statistic that identifies genes with extremely high or low expression values in a subset of the samples. RESULTS: We compare and discuss scores for bimodality for expression data. For the genome-wide identification of bimodal genes we apply all scores to expression data from 194 patients with node-negative breast cancer. Further, we present the first comprehensive genome-wide evaluation of the prognostic relevance of bimodal genes. We first rank genes according to bimodality scores and define two patient subgroups based on expression values. Then we assess the prognostic significance of the top ranking bimodal genes by comparing the survival functions of the two patient subgroups. We also evaluate the global association between the bimodal shape of expression distributions and survival times with an enrichment type analysis.Various cluster-based methods lead to a significant overrepresentation of prognostic genes. A striking result is obtained with the outlier-sum statistic (p < 10-12). Many genes with heavy tails generate subgroups of patients with different prognosis. CONCLUSIONS: Genes with high bimodality scores are promising candidates for defining prognostic patient subgroups from expression data. We discuss advantages and disadvantages of the different scores for prognostic purposes. The outlier-sum statistic may be particularly valuable for the identification of genes to be included in prognostic signatures. Among the genes identified as bimodal in the breast cancer data set several have not yet previously been recognized to be prognostic and bimodally expressed in breast cancer.

The G Protein-Coupled Bile Acid Receptor TGR5 (Gpbar1) Modulates Endothelin-1 Signaling in Liver.

TGR5 (Gpbar1) is a G protein-coupled receptor responsive to bile acids (BAs), which is expressed in different non-parenchymal cells of the liver, including biliary epithelial cells, liver-resident macrophages, sinusoidal endothelial cells (LSECs), and activated hepatic stellate cells (HSCs). Mice with targeted deletion of TGR5 are more susceptible towards cholestatic liver injury induced by cholic acid-feeding and bile duct ligation, resulting in a reduced proliferative response and increased liver injury. Conjugated lithocholic acid (LCA) represents the most potent TGR5 BA ligand and LCA-feeding has been used as a model to rapidly induce severe cholestatic liver injury in mice. Thus, TGR5 knockout (KO) mice and wildtype (WT) littermates were fed a diet supplemented with 1% LCA for 84 h. Liver injury and gene expression changes induced by the LCA diet revealed an enrichment of pathways associated with inflammation, proliferation, and matrix remodeling. Knockout of TGR5 in mice caused upregulation of endothelin-1 (ET-1) expression in the livers. Analysis of TGR5-dependent ET-1 signaling in isolated LSECs and HSCs demonstrated that TGR5 activation reduces ET-1 expression and secretion from LSECs and triggers internalization of the ET-1 receptor in HSCs, dampening ET-1 responsiveness. Thus, we identified two independent mechanisms by which TGR5 inhibits ET-1 signaling and modulates portal pressure.

Assessment of stem cell differentiation based on genome-wide expression profiles.

In recent years, protocols have been established to differentiate stem and precursor cells into more mature cell types. However, progress in this field has been hampered by difficulties to assess the differentiation status of stem cell-derived cells in an unbiased manner. Here, we present an analysis pipeline based on published data and methods to quantify the degree of differentiation and to identify transcriptional control factors explaining differences from the intended target cells or tissues. The pipeline requires RNA-Seq or gene array data of the stem cell starting population, derived 'mature' cells and primary target cells or tissue. It consists of a principal component analysis to represent global expression changes and to identify possible problems of the dataset that require special attention, such as: batch effects; clustering techniques to identify gene groups with similar features; over-representation analysis to characterize biological motifs and transcriptional control factors of the identified gene clusters; and metagenes as well as gene regulatory networks for quantitative cell-type assessment and identification of influential transcription factors. Possibilities and limitations of the analysis pipeline are illustrated using the example of human embryonic stem cell and human induced pluripotent cells to generate 'hepatocyte-like cells'. The pipeline quantifies the degree of incomplete differentiation as well as remaining stemness and identifies unwanted features, such as colon- and fibroblast-associated gene clusters that are absent in real hepatocytes but typically induced by currently available differentiation protocols. Finally, transcription factors responsible for incomplete and unwanted differentiation are identified. The proposed method is widely applicable and allows an unbiased and quantitative assessment of stem cell-derived cells.This article is part of the theme issue 'Designer human tissue: coming to a lab near you'.

Epsin Family Member 3 and Ribosome-Related Genes Are Associated with Late Metastasis in Estrogen Receptor-Positive Breast Cancer and Long-Term Survival in Non-Small Cell Lung Cancer Using a Genome-Wide Identification and Validation Strategy.

BACKGROUND: In breast cancer, gene signatures that predict the risk of metastasis after surgical tumor resection are mainly indicative of early events. The purpose of this study was to identify genes linked to metastatic recurrence more than three years after surgery. METHODS: Affymetrix HG U133A and Plus 2.0 array datasets with information on metastasis-free, disease-free or overall survival were accessed via public repositories. Time restricted Cox regression models were used to identify genes associated with metastasis during or after the first three years post-surgery (early- and late-type genes). A sequential validation study design, with two non-adjuvantly treated discovery cohorts (n = 409) and one validation cohort (n = 169) was applied and identified genes were further evaluated in tamoxifen-treated breast cancer patients (n = 923), as well as in patients with non-small cell lung (n = 1779), colon (n = 893) and ovarian (n = 922) cancer. RESULTS: Ten late- and 243 early-type genes were identified in adjuvantly untreated breast cancer. Adjustment to clinicopathological factors and an established proliferation-related signature markedly reduced the number of early-type genes to 16, whereas nine late-type genes still remained significant. These nine genes were associated with metastasis-free survival (MFS) also in a non-time restricted model, but not in the early period alone, stressing that their prognostic impact was primarily based on MFS more than three years after surgery. Four of the ten late-type genes, the ribosome-related factors EIF4B, RPL5, RPL3, and the tumor angiogenesis modifier EPN3 were significantly associated with MFS in the late period also in a meta-analysis of tamoxifen-treated breast cancer cohorts. In contrast, only one late-type gene (EPN3) showed consistent survival associations in more than one cohort in the other cancer types, being associated with worse outcome in two non-small cell lung cancer cohorts. No late-type gene was validated in ovarian and colon cancer. CONCLUSIONS: Ribosome-related genes were associated with decreased risk of late metastasis in both adjuvantly untreated and tamoxifen-treated breast cancer patients. In contrast, high expression of epsin (EPN3) was associated with increased risk of late metastasis. This is of clinical relevance considering the well-understood role of epsins in tumor angiogenesis and the ongoing development of epsin antagonizing therapies.

Interferon-inducible guanylate binding protein (GBP2) is associated with better prognosis in breast cancer and indicates an efficient T cell response.

BACKGROUND: Recently, interferon-inducible guanylate binding protein (GBP2) has been discussed as a possible control factor in tumor development, which is controlled by p53, and inhibits NF-Kappa B and Rac protein as well as expression of matrix metalloproteinase 9. However, the potential role that GBP2 plays in tumor development and prognosis has not yet been studied. METHODS: We analyzed whether GBP2 mRNA levels are associated with metastasis-free interval in 766 patients with node negative breast carcinomas who did not receive systemic chemotherapy. Furthermore, response to anthracycline-based chemotherapy was studied in 768 breast cancer patients. RESULTS: High expression of GBP2 in breast carcinomas was associated with better prognosis in the univariate (P < 0.001, hazard ratio 0.763, 95 % CI 0.650-0.896) as well as in the multivariate Cox analysis (P = 0.008, hazard ratio 0.731, 95 % CI 0.580-0.920) adjusted to the established clinical factors age, pT stage, grading, hormone and ERBB2 receptor status. The association was particularly strong in subgroups with high proliferation and positive estrogen receptor status but did not reach significance in carcinomas with low expression of proliferation associated genes. Besides its prognostic capacity, GBP2 also predicted pathologically complete response to anthracycline-based chemotherapy (P = 0.0037, odds ratio 1.39, 95 % CI 1.11-1.74). Interestingly, GBP2 correlated with a recently established T cell signature, indicating tumor infiltration with T cells (R = 0.607, P < 0.001). CONCLUSION: GBP2 is associated with better prognosis in fast proliferating tumors and probably represents a marker of an efficient T cell response.

Loss of circadian clock gene expression is associated with tumor progression in breast cancer.

Several studies suggest a link between circadian rhythm disturbances and tumorigenesis. However, the association between circadian clock genes and prognosis in breast cancer has not been systematically studied. Therefore, we examined the expression of 17 clock components in tumors from 766 node-negative breast cancer patients that were untreated in both neoadjuvant and adjuvant settings. In addition, their association with metastasis-free survival (MFS) and correlation to clinicopathological parameters were investigated. Aiming to estimate functionality of the clockwork, we studied clock gene expression relationships by correlation analysis. Higher expression of several clock genes (e.g., CLOCK, PER1, PER2, PER3, CRY2, NPAS2 and RORC) was found to be associated with longer MFS in univariate Cox regression analyses (HR<1 and FDR-adjusted P < 0.05). Stratification according to molecular subtype revealed prognostic relevance for PER1, PER3, CRY2 and NFIL3 in the ER+/HER2- subgroup, CLOCK and NPAS2 in the ER-/HER2- subtype, and ARNTL2 in HER2+ breast cancer. In the multivariate Cox model, only PER3 (HR = 0.66; P = 0.016) and RORC (HR = 0.42; P = 0.003) were found to be associated with survival outcome independent of established clinicopathological parameters. Pairwise correlations between functionally-related clock genes (e.g., PER2-PER3 and CRY2-PER3) were stronger in ER+, HER2- and low-grade carcinomas; whereas, weaker correlation coefficients were observed in ER- and HER2+ tumors, high-grade tumors and tumors that progressed to metastatic disease. In conclusion, loss of clock genes is associated with worse prognosis in breast cancer. Coordinated co-expression of clock genes, indicative of a functional circadian clock, is maintained in ER+, HER2-, low grade and non-metastasizing tumors but is compromised in more aggressive carcinomas.

The prognostic relevance of tumour-infiltrating plasma cells and immunoglobulin kappa C indicates an important role of the humoral immune response in non-small cell lung cancer.

A prognostic impact of immunoglobulin kappa C (IGKC) expression has been described in cancer. We analysed the influence of B-cell and plasma cell markers, as well as IGKC expression, in non-small lung cancer (NSCLC) using immunohistochemistry on a tissue microarray. IGKC protein expression was independently associated with longer survival, with particular impact in the adenocarcinoma subgroup. Moreover, a correlation was seen with CD138+ cells, but not with CD20. CD138 expression revealed a comparable association with survival. In conclusion, IGKC expression in stroma-infiltrating plasma cells is a prognostic marker in NSCLC, supporting emerging treatment concepts that exploit the humoral immune response.

Biomarker discovery in non-small cell lung cancer: integrating gene expression profiling, meta-analysis, and tissue microarray validation.

PURPOSE: Global gene expression profiling has been widely used in lung cancer research to identify clinically relevant molecular subtypes as well as to predict prognosis and therapy response. So far, the value of these multigene signatures in clinical practice is unclear, and the biologic importance of individual genes is difficult to assess, as the published signatures virtually do not overlap. EXPERIMENTAL DESIGN: Here, we describe a novel single institute cohort, including 196 non-small lung cancers (NSCLC) with clinical information and long-term follow-up. Gene expression array data were used as a training set to screen for single genes with prognostic impact. The top 450 probe sets identified using a univariate Cox regression model (significance level P < 0.01) were tested in a meta-analysis including five publicly available independent lung cancer cohorts (n = 860). RESULTS: The meta-analysis revealed 14 genes that were significantly associated with survival (P < 0.001) with a false discovery rate <1%. The prognostic impact of one of these genes, the cell adhesion molecule 1 (CADM1), was confirmed by use of immunohistochemistry on tissue microarrays from 2 independent NSCLC cohorts, altogether including 617 NSCLC samples. Low CADM1 protein expression was significantly associated with shorter survival, with particular influence in the adenocarcinoma patient subgroup. CONCLUSIONS: Using a novel NSCLC cohort together with a meta-analysis validation approach, we have identified a set of single genes with independent prognostic impact. One of these genes, CADM1, was further established as an immunohistochemical marker with a potential application in clinical diagnostics.