Validation study of existing gene expression signatures for anti-TNF treatment in patients with rheumatoid arthritis

So far, there are no means of identifying rheumatoid arthritis (RA) patients who will fail to respond to tumour necrosis factor blocking agents (anti-TNF), prior to treatment. We set out to validate eight previously reported gene expression signatures predicting therapy outcome. Genome-wide expression profiling using Affymetrix GeneChip Exon 1.0 ST arrays was performed on RNA isolated from whole blood of 42 RA patients starting treatment with infliximab or adalimumab. Clinical response according to EULAR criteria was determined at week 14 of therapy. Genes that have been reported to be associated with anti-TNF treatment were extracted from our dataset. K-means partition clustering was performed to assess the predictive value of the gene-sets. We performed a hypothesis-driven analysis of the dataset using eight existing gene sets predictive of anti-TNF treatment outcome. The set that performed best reached a sensitivity of 71% and a specificity of 61%, for classifying the patients in the current study. We successfully validated one of eight previously reported predictive expression profile. This replicated expression signature is a good starting point for developing a prediction model for anti-TNF treatment outcome that can be used in a daily clinical setting. Our results confirm that gene expression profiling prior to treatment is a useful tool to predict anti-TNF (non) response.     

IFI35, mir-99a and HCV Genotype to Predict Sustained Virological Response to Pegylated-Interferon Plus Ribavirin in Chronic Hepatitis C

Although, the treatment of chronic hepatitis C (CHC) greatly improved with the use of direct antiviral agents, pegylated-interferon (PEG-IFN) plus ribavirin remains an option for many patients, worldwide. The intra-hepatic level of expression of interferon stimulated genes (ISGs) and the rs12979860 CC genotype located within IFNL3 have been associated with sustained virological response (SVR), in patients with CHC. The aim of the study was to identify micro-RNAs associated with SVR and to build an accurate signature to predict SVR. Pre-treatment liver biopsies from 111 patients, treated with PEG-IFN plus ribavirin, were studied. Fifty-seven patients had SVR, 36 non-response (NR) and 18 relapse (RR). The expression of 851 human miRNAs and 30 selected mRNAs, including ISGs, was assessed by RT-qPCR. In the first group of patients (screen), 20 miRNAs out of the 851 studied were deregulated between NRs and SVRs. From the 4 miRNAs validated (mir-23a, mir-181a*, mir-217 and mir-99a), in the second group of patients (validation), 3 (mir-23a, mir-181a* and mir-99a) were down-regulated in NRs as compared to SVRs. The ISGs, studied, were accumulated in SVRs and IFNL3 rs12979860 CT/TT carriers compared respectively to NRs and CC carriers. Combining, clinical data together with the expression of selected genes and micro-RNAs, we identified a signature (IFI35, mir-99a and HCV genotype) to predict SVR (AUC:0.876) with a positive predictive value of 86.54% with high sensibility (80%) and specificity (80.4%). This signature may help to characterize patients with low chance to respond to PEG-IFN/ribavirin and to elucidate mechanisms of NR.     

The Impact of IL28B Genotype and Liver Fibrosis on the Hepatic Expression of IP10, IFI27, ISG15, and MX1 and Their Association with Treatment Outcomes in Patients with Chronic Hepatitis C

The strong impact of interleukin 28B (IL28B) polymorphisms on sustained virological response (SVR) after peginterferon and ribavirin treatment in patients with chronic hepatitis C (CHC) is well-known. We investigated IL28B variability and hepatic expression of IP10, IFI27, ISG15, and MX1 in CHC patients, the relation of each with their clinical characteristics, and how they associated with responses to combined therapy. Genotyping and gene expression analysis were conducted in a selected cohort of treatment-naïve patients who underwent interferon and ribavirin treatment. Differential expression of IP10, IFI27, ISG15, and MX1 genes was assessed from pretreatment liver biopsies using quantitative PCR. Histopathological evaluation of liver specimens was performed on the basis of the Scheuer’s modified scale. We showed that hepatic IFI27, ISG15, and MX1 expression was lower in the IL28B CC 12979860 and TT rs8099917 groups than in the CT-TT rs12979860 and TG-GG rs8099917 groups (P < 0.001). We found no differences in IP10 expression between the IL28B genotypes (P > 0.05); in contrast, IP10 expression was significantly affected by the progression of fibrosis (P = 0.007). We showed that the rs12979860 CC genotype was associated with successful treatment when compared to the rs12979860 CT-TT genotype (P = 0.004). Additionally, the expression levels of IP10, IFI27 and ISG15, but not MX1, were significantly higher in non-SVR patients than in SVR patients. The effect of variation in IL28B on the results of IFN-based treatment may be associated with changes in IFI27 and ISG15, but not with IP10. Silencing of IP10 is positive and independent from IL28B prediction of SVR, which is strongly associated with liver fibrosis in CHC patients.     

Polycomb Target Genes Are Silenced in Multiple Myeloma

Multiple myeloma (MM) is a genetically heterogeneous disease, which to date remains fatal. Finding a common mechanism for initiation and progression of MM continues to be challenging. By means of integrative genomics, we identified an underexpressed gene signature in MM patient cells compared to normal counterpart plasma cells. This profile was enriched for previously defined H3K27-tri-methylated genes, targets of the Polycomb group (PcG) proteins in human embryonic fibroblasts. Additionally, the silenced gene signature was more pronounced in ISS stage III MM compared to stage I and II. Using chromatin immunoprecipitation (ChIP) assay on purified CD138+ cells from four MM patients and on two MM cell lines, we found enrichment of H3K27me3 at genes selected from the profile. As the data implied that the Polycomb-targeted gene profile would be highly relevant for pharmacological treatment of MM, we used two compounds to chemically revert the H3K27-tri-methylation mediated gene silencing. The S-adenosylhomocysteine hydrolase inhibitor 3-Deazaneplanocin (DZNep) and the histone deacetylase inhibitor LBH589 (Panobinostat), reactivated the expression of genes repressed by H3K27me3, depleted cells from the PRC2 component EZH2 and induced apoptosis in human MM cell lines. In the immunocompetent 5T33MM in vivo model for MM, treatment with LBH589 resulted in gene upregulation, reduced tumor load and increased overall survival. Taken together, our results reveal a common gene signature in MM, mediated by gene silencing via the Polycomb repressor complex. The importance of the underexpressed gene profile in MM tumor initiation and progression should be subjected to further studies.     

Development and validation of a prognostic gene-expression signature for lung adenocarcinoma

Although several prognostic signatures have been developed in lung cancer, their application in clinical practice has been limited because they have not been validated in multiple independent data sets. Moreover, the lack of common genes between the signatures makes it difficult to know what biological process may be reflected or measured by the signature. By using classical data exploration approach with gene expression data from patients with lung adenocarcinoma (n = 186), we uncovered two distinct subgroups of lung adenocarcinoma and identified prognostic 193-gene gene expression signature associated with two subgroups. The signature was validated in 4 independent lung adenocarcinoma cohorts, including 556 patients. In multivariate analysis, the signature was an independent predictor of overall survival (hazard ratio, 2.4; 95% confidence interval, 1.2 to 4.8; p = 0.01). An integrated analysis of the signature revealed that E2F1 plays key roles in regulating genes in the signature. Subset analysis demonstrated that the gene signature could identify high-risk patients in early stage (stage I disease), and patients who would have benefit of adjuvant chemotherapy. Thus, our study provided evidence for molecular basis of clinically relevant two distinct two subtypes of lung adenocarcinoma.     

Identification of CD28 and PTEN as novel prognostic markers for cervical cancer

Cervical cancer (CC) is the most common malignant tumor with poor clinical outcome among women. Identification of novel biomarkers could be beneficial for the clinical diagnosis and treatment of CC. This study aimed to identify prognostic biomarkers for the prediction of prognostic status of CC patients, and explore the effect of the corresponding methylated genes in the occurrence and development of CC. The methylation microarray data of CC was extracted from The Cancer Genome Atlas (TCGA) dataset. The methylation genes associated with the prognostic status were identified based on the information of the relapse-free survival (RFS) of the CC patients. The prognostic gene pairs were further identified. Then, the prognostic signature was identified by the forward search algorithm based on the C-index method. The results were validated by independent dataset. Finally, the functional analysis was performed on the methylation genes. A total of 276 methylation genes and 2508 gene pairs associated with the prognostic status of the CC were identified. A signature composed of eight methylation gene pairs was obtained to predict the prognostic status of cervical patients. A series of genes that played an important role in the occurrence and development of CC were obtained by the functional enrichment analysis. To summary, a prognostic signature consisting of eight methylation gene pairs was obtained. Of note, the CD28 and PTEN gene pair were found to play important roles in the occurrence and development of CC.     

An erythroid differentiation signature predicts response to lenalidomide in myelodysplastic syndrome

Background

Lenalidomide is an effective new agent for the treatment of patients with myelodysplastic syndrome (MDS), an acquired hematopoietic disorder characterized by ineffective blood cell production and a predisposition to the development of leukemia. Patients with an interstitial deletion of Chromosome 5q have a high rate of response to lenalidomide, but most MDS patients lack this deletion. Approximately 25% of patients without 5q deletions also benefit from lenalidomide therapy, but response in these patients cannot be predicted by any currently available diagnostic assays. The aim of this study was to develop a method to predict lenalidomide response in order to avoid unnecessary toxicity in patients unlikely to benefit from treatment.

Methods and Findings

Using gene expression profiling, we identified a molecular signature that predicts lenalidomide response. The signature was defined in a set of 16 pretreatment bone marrow aspirates from MDS patients without 5q deletions, and validated in an independent set of 26 samples. The response signature consisted of a cohesive set of erythroid-specific genes with decreased expression in responders, suggesting that a defect in erythroid differentiation underlies lenalidomide response. Consistent with this observation, treatment with lenalidomide promoted erythroid differentiation of primary hematopoietic progenitor cells grown in vitro.

Conclusions

These studies indicate that lenalidomide-responsive patients have a defect in erythroid differentiation, and suggest a strategy for a clinical test to predict patients most likely to respond to the drug. The experiments further suggest that the efficacy of lenalidomide, whose mechanism of action in MDS is unknown, may be due to its ability to induce erythroid differentiation.     

Involvement of MAP3K8 and miR-17-5p in Poor Virologic Response to Interferon-Based Combination Therapy for Chronic Hepatitis C

Despite advances in chronic hepatitis C treatment, a proportion of patients respond poorly to treatment. This study aimed to explore hepatic mRNA and microRNA signatures involved in hepatitis C treatment resistance. Global hepatic mRNA and microRNA expression profiles were compared using microarray data between treatment responses. Quantitative real-time polymerase chain reaction validated the gene signatures from 130 patients who were infected with hepatitis C virus genotype 1b and treated with pegylated interferon-alpha and ribavirin combination therapy. The correlation between mRNA and microRNA was evaluated using in silico analysis and in vitro siRNA and microRNA inhibition/overexpression experiments. Multivariate regression analysis identified that the independent variables IL28B SNP rs8099917, hsa-miR-122-5p, hsa-miR-17-5p, and MAP3K8 were significantly associated with a poor virologic response. MAP3K8 and miR-17-5p expression were inversely correlated with treatment response. Furthermore, miR-17-5p repressed HCV production by targeting MAP3K8. Collectively, the data suggest that several molecules and the inverse correlation between mRNA and microRNA contributed to a host genetic refractory hepatitis C treatment response.     

Local synthesis of interferon-alpha in lupus nephritis is associated with type I interferons signature and LMP7 induction in renal tubular epithelial cells

IntroductionType I interferons are pivotal in the activation of autoimmune response in systemic lupus erythematous. However, the pathogenic role of interferon-alpha in patients affected by lupus nephritis remains uncertain. The aim of our study was to investigate the presence of a specific interferon signature in lupus nephritis and the effects of interferon-alpha at renal level.MethodsWe performed immunohistochemical analysis for MXA-protein and in situ hybridization to detect interferon-alpha signature and production in human lupus nephritis. Through microarray studies, we analyzed the gene expression profile of renal tubular epithelial cells, stimulated with interferon-alpha. We validated microarray results through real-time polymerase chain reaction, flow cytometry on renal tubular epithelial cells, and through immunohistochemical analysis and confocal microscopy on renal biopsies.ResultsType I interferons signature was characterized by MXA-specific staining in renal tubular epithelial cells; in addition, in situ hybridization showed that renal tubular epithelial cells were the major producers of interferon-alpha, indicating a potential autocrine effect. Whole-genome expression profile showed interferon-alpha induced up-regulation of genes involved in innate immunity, protein ubiquitination and switching to immunoproteasome. In accordance with the in vitro data, class IV lupus nephritis showed up-regulation of the immunoproteasome subunit LMP7 in tubular epithelial cells associated with type I interferon signature.ConclusionsOur data indicate that type I interferons might have a pathogenic role in lupus nephritis characterized by an autocrine effect of interferon-alpha on renal tubular epithelial cells. Therefore we hypothesize that inhibition of type I interferons might represent a therapeutic target to prevent tubulo-interstitial damage in patients with lupus nephritis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0588-3) contains supplementary material, which is available to authorized users.     

A Translatable Predictor of Human Radiation Exposure

Terrorism using radiological dirty bombs or improvised nuclear devices is recognized as a major threat to both public health and national security. In the event of a radiological or nuclear disaster, rapid and accurate biodosimetry of thousands of potentially affected individuals will be essential for effective medical management to occur. Currently, health care providers lack an accurate, high-throughput biodosimetric assay which is suitable for the triage of large numbers of radiation injury victims. Here, we describe the development of a biodosimetric assay based on the analysis of irradiated mice, ex vivo-irradiated human peripheral blood (PB) and humans treated with total body irradiation (TBI). Interestingly, a gene expression profile developed via analysis of murine PB radiation response alone was inaccurate in predicting human radiation injury. In contrast, generation of a gene expression profile which incorporated data from ex vivo irradiated human PB and human TBI patients yielded an 18-gene radiation classifier which was highly accurate at predicting human radiation status and discriminating medically relevant radiation dose levels in human samples. Although the patient population was relatively small, the accuracy of this classifier in discriminating radiation dose levels in human TBI patients was not substantially confounded by gender, diagnosis or prior exposure to chemotherapy. We have further incorporated genes from this human radiation signature into a rapid and high-throughput chemical ligation-dependent probe amplification assay (CLPA) which was able to discriminate radiation dose levels in a pilot study of ex vivo irradiated human blood and samples from human TBI patients. Our results illustrate the potential for translation of a human genetic signature for the diagnosis of human radiation exposure and suggest the basis for further testing of CLPA as a candidate biodosimetric assay.     

Met Kinetic Signature Derived from the Response to HGF/SF in a Cellular Model Predicts Breast Cancer Patient Survival

To determine the signaling pathways leading from Met activation to metastasis and poor prognosis, we measured the kinetic gene alterations in breast cancer cell lines in response to HGF/SF. Using a network inference tool we analyzed the putative protein-protein interaction pathways leading from Met to these genes and studied their specificity to Met and prognostic potential. We identified a Met kinetic signature consisting of 131 genes. The signature correlates with Met activation and with response to anti-Met therapy (p<0.005) in in-vitro models. It also identifies breast cancer patients who are at high risk to develop an aggressive disease in six large published breast cancer patient cohorts (p<0.01, N>1000). Moreover, we have identified novel putative Met pathways, which correlate with Met activity and patient prognosis. This signature may facilitate personalized therapy by identifying patients who will respond to anti-Met therapy. Moreover, this novel approach may be applied for other tyrosine kinases and other malignancies.     

Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features

With the goal of studying epigenetic alterations in fibrolamellar hepatocellular carcinoma (FLC) and establish an associated DNA methylation signature, we analyzed LINE-1 methylation in a cohort of FLC and performed next-generation sequencing of DNA methylation in a training set of pure-FLCs and non-cirrhotic hepatocellular carcinomas (nc-HCC). DNA methylation was correlated with gene expression. Furthermore, we established and validated an epigenetic signature differentiating pure-FLC from other HCCs. LINE-1 methylation correlated with shorter recurrence-free survival and overall survival in resected pure-FLC patients. Unsupervised clustering using CG sites located in islands distinguished pure-FLC from nc-HCC. Major DNA methylation changes occurred outside promoters, mainly in gene bodies and intergenic regions located in the vicinity of liver developmental genes (i.e., SMARCA4 and RXRA). Partially methylated domains were more prone to DNA methylation changes. Furthermore, we identified several putative tumor suppressor genes (e.g., DLEU7) and oncogenes (e.g., DUSP4). While ∼70% of identified gene promoters gaining methylation were marked by bivalent histone marks (H3K4me3/H3K27me3) in embryonic stem cells, ∼70% of those losing methylation were marked by H3K4me3. Finally, we established a pure FLC DNA methylation signature and validated it in an independent dataset. Our analysis reveals a distinct epigenetic signature of pure FLC as compared to nc-HCC, with DNA methylation changes occurring in the vicinity of liver developmental genes. These data suggest new options for targeting FLC based on cancer epigenome aberrations.     

Up-regulation of the novel proinflammatory adipokines lipocalin-2, chitinase-3 like-1 and osteopontin as well as angiogenic-related factors in visceral adipose tissue of patients with colon cancer

Background

Obesity is widely recognised as an important risk factor for colorectal cancer (CC).

Aim

The study aimed to evaluate the effect of CC on circulating concentrations and gene expression levels of inflammatory and angiogenesis-related factors in human visceral adipose tissue (VAT).

Methods

VAT biopsies were obtained from 18 healthy individuals and 11 patients with CC. Real-time polymerase chain reactions were performed to quantify gene expression levels and zymographic analyses were used to determine the activity of matrix metalloproteinases (MMPs).

Results

Patients with CC exhibited increased mRNA expression levels of lipocalin-2 (P=.014), osteopontin (P=.027), tumor necrosis factor-α (TNF-α) (P=.016) and chitinase-3 like-1 (P=.006) compared to control subjects in VAT. Gene expression levels of hypoxia-inducible factor-1 α, vascular endothelial growth factor and MMP-2 were significantly higher (P<.05) in VAT of patients with CC. The expression of insulin-like growth factor I, insulin growth factor binding protein 3 and MMP-9 followed the same trend, although no significant differences were reached. The enzymatic activity of MMP-9 was increased (P<.001) in patients with CC. Furthermore, individuals with CC showed increased (P<.05) circulating concentrations of the inflammatory markers interleukin-6, tumour necrosis factor α and hepatocyte growth factor, whereas levels of the anti-inflammatory adipokine adiponectin were decreased (P<.01).

Conclusion

These findings represent the first observation that mRNA levels of the novel inflammatory factors lipocalin-2, chitinase-3 like-1 and osteopontin are increased in human VAT of subjects with CC. This observation together with the up-regulation of angiogenic factors suggests that adipokines secreted by VAT may be involved in the development of colon cancer.     

Expression Profiling of Rectal Tumors Defines Response to Neoadjuvant Treatment Related Genes

To date, no effective method exists that predicts the response to preoperative chemoradiation (CRT) in locally advanced rectal cancer (LARC). Nevertheless, identification of patients who have a higher likelihood of responding to preoperative CRT could be crucial in decreasing treatment morbidity and avoiding expensive and time-consuming treatments. The aim of this study was to identify signatures or molecular markers related to response to pre-operative CRT in LARC. We analyzed the gene expression profiles of 26 pre-treatment biopsies of LARC (10 responders and 16 non-responders) without metastasis using Human WG CodeLink microarray platform. Two hundred and fifty seven genes were differentially over-expressed in the responder patient subgroup. Ingenuity Pathway Analysis revealed a significant ratio of differentially expressed genes related to cancer, cellular growth and proliferation pathways, and c-Myc network. We demonstrated that high Gng4, c-Myc, Pola1, and Rrm1 mRNA expression levels was a significant prognostic factor for response to treatment in LARC patients (p<0.05). Using this gene set, we were able to establish a new model for predicting the response to CRT in rectal cancer with a sensitivity of 60% and 100% specificity. Our results reflect the value of gene expression profiling to gain insight about the molecular pathways involved in the response to treatment of LARC patients. These findings could be clinically relevant and support the use of mRNA levels when aiming to identify patients who respond to CRT therapy.