IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

Aggregate cultures of foetal rat liver cells: development and maintenance of liver gene expression.

Rotation-mediated aggregate cultures of foetal rat liver cells were prepared and grown in a chemically defined medium. Their capacity for cellular organisation and maturation was studied over a culture period of 3 wk by using both morphologic and biochemical criteria. It was found that within each aggregate, distinct liver cell types were present and attained their normal, differentiated phenotype. Parenchymal cells formed small acini with a central lumen. Within the first 2 wk in culture, albumin and ferritin mRNA levels were maintained, while the alpha-fetoprotein mRNA levels decreased, and tyrosine aminotransferase (TAT) gene expression increased. No significant response to glucocorticoids was observed in early cultures, whereas after 3 wk a marked increase in TAT mRNA levels was elicited by dexamethasone and glucagon (additive stimulatory effects). The results show that foetal rat liver cells cultured in a chemically defined medium are able to rearrange themselves into histotypic structures, and display a developmental pattern of gene expression comparable to that of perinatal rat liver in vivo. This culture system offers therefore a useful model to study the development and function of liver cells.

Liver hyperplasia and paradoxical regulation of glycogen metabolism and glucose-sensitive gene expression in GLUT2-null hepatocytes. Further evidence for the existence of a membrane-based glucose release pathway.

We investigated the impact of GLUT2 gene inactivation on the regulation of hepatic glucose metabolism during the fed to fast transition. In control and GLUT2-null mice, fasting was accompanied by a approximately 10-fold increase in plasma glucagon to insulin ratio, a similar activation of liver glycogen phosphorylase and inhibition of glycogen synthase and the same elevation in phosphoenolpyruvate carboxykinase and glucose-6-phosphatase mRNAs. In GLUT2-null mice, mobilization of glycogen stores was, however, strongly impaired. This was correlated with glucose-6-phosphate (G6P) levels, which remained at the fed values, indicating an important allosteric stimulation of glycogen synthase by G6P. These G6P levels were also accompanied by a paradoxical elevation of the mRNAs for L-pyruvate kinase. Re-expression of GLUT2 in liver corrected the abnormal regulation of glycogen and L-pyruvate kinase gene expression. Interestingly, GLUT2-null livers were hyperplasic, as revealed by a 40% increase in liver mass and 30% increase in liver DNA content. Together, these data indicate that in the absence of GLUT2, the G6P levels cannot decrease during a fasting period. This may be due to neosynthesized glucose entering the cytosol, being unable to diffuse into the extracellular space, and being phosphorylated back to G6P. Because hepatic glucose production is nevertheless quantitatively normal, glucose produced in the endoplasmic reticulum may also be exported out of the cell through an alternative, membrane traffic-based pathway, as previously reported (Guillam, M.-T., Burcelin, R., and Thorens, B. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 12317-12321). Therefore, in fasting, GLUT2 is not required for quantitative normal glucose output but is necessary to equilibrate cytosolic glucose with the extracellular space. In the absence of this equilibration, the control of hepatic glucose metabolism by G6P is dominant over that by plasma hormone concentrations.

Transcriptional profiling reveals divergent roles of PPARalpha and PPARbeta/delta in regulation of gene expression in mouse liver.

Little is known about the role of the transcription factor peroxisome proliferator-activated receptor (PPAR) beta/delta in liver. Here we set out to better elucidate the function of PPARbeta/delta in liver by comparing the effect of PPARalpha and PPARbeta/delta deletion using whole genome transcriptional profiling and analysis of plasma and liver metabolites. In fed state, the number of genes altered by PPARalpha and PPARbeta/delta deletion was similar, whereas in fasted state the effect of PPARalpha deletion was much more pronounced, consistent with the pattern of gene expression of PPARalpha and PPARbeta/delta. Minor overlap was found between PPARalpha- and PPARbeta/delta-dependent gene regulation in liver. Pathways upregulated by PPARbeta/delta deletion were connected to innate immunity and inflammation. Pathways downregulated by PPARbeta/delta deletion included lipoprotein metabolism and various pathways related to glucose utilization, which correlated with elevated plasma glucose and triglycerides and reduced plasma cholesterol in PPARbeta/delta-/- mice. Downregulated genes that may underlie these metabolic alterations included Pklr, Fbp1, Apoa4, Vldlr, Lipg, and Pcsk9, which may represent novel PPARbeta/delta target genes. In contrast to PPARalpha-/- mice, no changes in plasma free fatty acid, plasma beta-hydroxybutyrate, liver triglycerides, and liver glycogen were observed in PPARbeta/delta-/- mice. Our data indicate that PPARbeta/delta governs glucose utilization and lipoprotein metabolism and has an important anti-inflammatory role in liver. Overall, our analysis reveals divergent roles of PPARalpha and PPARbeta/delta in regulation of gene expression in mouse liver.

Real-time recording of circadian liver gene expression in freely moving mice reveals the phase-setting behavior of hepatocyte clocks.

The mammalian circadian timing system consists of a master pacemaker in the suprachiasmatic nucleus (SCN) in the hypothalamus, which is thought to set the phase of slave oscillators in virtually all body cells. However, due to the lack of appropriate in vivo recording technologies, it has been difficult to study how the SCN synchronizes oscillators in peripheral tissues. Here we describe the real-time recording of bioluminescence emitted by hepatocytes expressing circadian luciferase reporter genes in freely moving mice. The technology employs a device dubbed RT-Biolumicorder, which consists of a cylindrical cage with reflecting conical walls that channel photons toward a photomultiplier tube. The monitoring of circadian liver gene expression revealed that hepatocyte oscillators of SCN-lesioned mice synchronized more rapidly to feeding cycles than hepatocyte clocks of intact mice. Hence, the SCN uses signaling pathways that counteract those of feeding rhythms when their phase is in conflict with its own phase.

Requirement of glucose metabolism for regulation of glucose transporter type 2 (GLUT2) gene expression in liver.

Previous studies have shown that glucose increases the glucose transporter (GLUT2) mRNA expression in the liver in vivo and in vitro. Here we report an analysis of the effects of glucose metabolism on GLUT2 gene expression. GLUT2 mRNA accumulation by glucose was not due to stabilization of its transcript but rather was a direct effect on gene transcription. A proximal fragment of the 5' regulatory region of the mouse GLUT2 gene linked to a reporter gene was transiently transfected into liver GLUT2-expressing cells. Glucose stimulated reporter gene expression in these cells, suggesting that glucose-responsive elements were included within the proximal region of the promoter. A dose-dependent effect of glucose on GLUT2 expression was observed over 10 mM glucose irrespective of the hexokinase isozyme (glucokinase K(m) 16 mM; hexokinase I K(m) 0.01 mM) present in the cell type used. This suggests that the correlation between extracellular glucose and GLUT2 mRNA concentrations is simply a reflection of an activation of glucose metabolism. The mediators and the mechanism responsible for this response remain to be determined. In conclusion, glucose metabolism is required for the proper induction of the GLUT2 gene in the liver and this effect is transcriptionally regulated.

A role for altered TLR gene expression in association with increased expression of CD200R in the induction of mucosal tissue CD4(+) Treg in aged mice following gavage with a liver extract along with intramuscular monophosphoryl lipid A (MPLA) injection.

Previous studies showed a fetal sheep liver extract (FSLE), in association with LPS, injected into aged (>20 months) mice reversed the altered polarization (increased IL-4 and IL-10 with decreased IL-2 and IFN-gamma) in cytokine production seen from ConA stimulated lymphoid cells of those mice. Aged mice show a >60% decline in numbers and suppressive function of both CD4(+)CD25(+)Foxp3(+)Treg and so-called Tr3 (CD4(+)TGFbeta(+)). Their number/function is restored to levels seen in control (8-week-old) mice by FSLE. We have reported at length on the ability of a novel pair of immunoregulatory molecules, members of the TREM family, namely CD200:CD200R, to control development of dendritic cells (DCs) which themselves regulate production of Foxp3(+) Treg. The latter express a distinct subset of TLRs which control their function. We report that a feature of the altered Treg expression following combined treatment with FSLE and monophosphoryl lipid A, MPLA (a bioactive component of lipid A of LPS) is the altered gene expression both of distinct subsets of TLRs and of CD200Rs. We speculate that this may represent one of the mechanisms by which FSLE and MPLA alter immunity in aged mice.

IFN Stimulated Gene Expression in the Liver is a Better Predictor of Treatment Response in Chronic Hepatitis C than the IL28B (IFN lambda 3) Genotype

BACKGROUND: Therapy of chronic hepatitis C (CHC) with pegIFNα/ribavirin achieves a sustained virologic response (SVR) in ∼55%. Pre-activation of the endogenous interferon system in the liver is associated with non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. METHODS: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. RESULTS: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. CONCLUSIONS: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

Differential involvement of peroxisome-proliferator-activated receptors alpha and delta in fibrate and fatty-acid-mediated inductions of the gene encoding liver fatty-acid-binding protein in the liver and the small intestine.

Liver fatty-acid-binding protein (L-FABP) is a cytoplasmic polypeptide that binds with strong affinity especially to long-chain fatty acids (LCFAs). It is highly expressed in both the liver and small intestine, where it is thought to have an essential role in the control of the cellular fatty acid (FA) flux. Because expression of the gene encoding L-FABP is increased by both fibrate hypolipidaemic drugs and LCFAs, it seems to be under the control of transcription factors, termed peroxisome-proliferator-activated receptors (PPARs), activated by fibrate or FAs. However, the precise molecular mechanism by which these regulations take place remain to be fully substantiated. Using transfection assays, we found that the different PPAR subtypes (alpha, gamma and delta) are able to mediate the up-regulation by FAs of the gene encoding L-FABP in vitro. Through analysis of LCFA- and fibrate-mediated effects on L-FABP mRNA levels in wild-type and PPARalpha-null mice, we have found that PPARalpha in the intestine does not constitute a dominant regulator of L-FABP gene expression, in contrast with what is known in the liver. Only the PPARdelta/alpha agonist GW2433 is able to up-regulate the gene encoding L-FABP in the intestine of PPARalpha-null mice. These findings demonstrate that PPARdelta can act as a fibrate/FA-activated receptor in tissues in which it is highly expressed and that L-FABP is a PPARdelta target gene in the small intestine. We propose that PPARdelta contributes to metabolic adaptation of the small intestine to changes in the lipid content of the diet.

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Infrequent promoter methylation of the MGMT gene in liver metastases from uveal melanoma.

Uveal melanoma is associated with a high mortality rate once metastases occur, with over >90% of metastatic patients dying within less than 1 year from metastases to the liver. The intraarterial hepatic (iah) administration of the alkylating agent fotemustine holds some promise with response rates of 36% and median survival of 15 months. Here, we investigated whether the DNA-repair-protein MGMT may be involved in the variability of response to fotemustine and temozolomide in uveal melanoma. Epigenetic inactivation of MGMT has been demonstrated to be a predictive marker for benefit from alkylating agent therapy in glioblastoma. We found a methylated MGMT promoter in 6% of liver metastases from 34 uveal melanoma patients. The mean MGMT activity measured in liver metastases with negligible liver tissue content was significantly lower than in liver tissue (146 versus 523 fmol/mg protein, p = 0.002). Expression of the MGMT protein was detectable in 50% of 88 metastases by immunohistochemistry on a tissue microarray. Expression was heterogeneous, and in accordance with MGMT activity data, usually lower than in the surrounding liver. Differential MGMT activity/expression between metastasis and liver tissue and more efficient depletion of MGMT with higher doses of alkylating agent therapy using iah delivery may provide the pharmacologic window for the higher response rate. However, these results do not support MGMT methylation status or protein expression as predictive markers for treatment outcome to iah chemotherapy with alkylating agents.

An autoregulatory loop controlling CYP1A1 gene expression: role of H(2)O(2) and NFI.

Cytochrome P450 1A1 (CYP1A1), like many monooxygenases, can produce reactive oxygen species during its catalytic cycle. Apart from the well-characterized xenobiotic-elicited induction, the regulatory mechanisms involved in the control of the steady-state activity of CYP1A1 have not been elucidated. We show here that reactive oxygen species generated from the activity of CYP1A1 limit the levels of induced CYP1A1 mRNAs. The mechanism involves the repression of the CYP1A1 gene promoter activity in a negative-feedback autoregulatory loop. Indeed, increasing the CYP1A1 activity by transfecting CYP1A1 expression vectors into hepatoma cells elicited an oxidative stress and led to the repression of a reporter gene driven by the CYP1A1 gene promoter. This negative autoregulation is abolished by ellipticine (an inhibitor of CYP1A1) and by catalase (which catalyzes H(2)O(2) catabolism), thus implying that H(2)O(2) is an intermediate. Down-regulation is also abolished by the mutation of the proximal nuclear factor I (NFI) site in the promoter. The transactivating domain of NFI/CTF was found to act in synergy with the arylhydrocarbon receptor pathway during the induction of CYP1A1 by 2,3,7,8-tetrachloro-p-dibenzodioxin. Using an NFI/CTF-Gal4 fusion, we show that NFI/CTF transactivating function is decreased by a high activity of CYP1A1. This regulation is also abolished by catalase or ellipticine. Consistently, the transactivating function of NFI/CTF is repressed in cells treated with H(2)O(2), a novel finding indicating that the transactivating domain of a transcription factor can be targeted by oxidative stress. In conclusion, an autoregulatory loop leads to the fine tuning of the CYP1A1 gene expression through the down-regulation of NFI activity by CYP1A1-based H(2)O(2) production. This mechanism allows a limitation of the potentially toxic CYP1A1 activity within the cell.

Interferon-Induced Gene Expression Is a Stronger Predictor of Treatment Response Than IL28B Genotype in Patients With Hepatitis C.

BACKGROUND & AIMS: The host immune response during the chronic phase of hepatitis C virus infection varies among individuals; some patients have a no interferon (IFN) response in the liver, whereas others have full activation IFN-stimulated genes (ISGs). Preactivation of this endogenous IFN system is associated with nonresponse to pegylated IFN-α (pegIFN-α) and ribavirin. Genome-wide association studies have associated allelic variants near the IL28B (IFNλ3) gene with treatment response. We investigated whether IL28B genotype determines the constitutive expression of ISGs in the liver and compared the abilities of ISG levels and IL28B genotype to predict treatment outcome. METHODS: We genotyped 109 patients with chronic hepatitis C for IL28B allelic variants and quantified the hepatic expression of ISGs and of IL28B. Decision tree ensembles, in the form of a random forest classifier, were used to calculate the relative predictive power of these different variables in a multivariate analysis. RESULTS: The minor IL28B allele was significantly associated with increased expression of ISG. However, stratification of the patients according to treatment response revealed increased ISG expression in nonresponders, irrespective of IL28B genotype. Multivariate analysis of ISG expression, IL28B genotype, and several other factors associated with response to therapy identified ISG expression as the best predictor of treatment response. CONCLUSIONS: IL28B genotype and hepatic expression of ISGs are independent predictors of response to treatment with pegIFN-α and ribavirin in patients with chronic hepatitis C. The most accurate prediction of response was obtained with a 4-gene classifier comprising IFI27, ISG15, RSAD2, and HTATIP2.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

A 338-bp proximal fragment of the glucose transporter type 2 (GLUT2) promoter drives reporter gene expression in the pancreatic islets of transgenic mice.

The high Km glucose transporter GLUT2 is a membrane protein expressed in tissues involved in maintaining glucose homeostasis, and in cells where glucose-sensing is necessary. In many experimental models of diabetes, GLUT2 gene expression is decreased in pancreatic beta-cells, which could lead to a loss of glucose-induced insulin secretion. In order to identify factors involved in pancreatic beta-cell specific expression of GLUT2, we have recently cloned the murine GLUT2 promoter and identified cis-elements within the 338-bp of the proximal promoter capable of binding islet-specific trans-acting factors. Furthermore, in transient transfection studies, this 338-bp fragment could efficiently drive the expression of the chloramphenicol acetyl transferase (CAT) gene in cell lines derived from the endocrine pancreas, but displayed no promoter activity in non-pancreatic cells. In this report, we tested the cell-specific expression of a CAT reporter gene driven by a short (338 bp) and a larger (1311 bp) fragment of the GLUT2 promoter in transgenic mice. We generated ten transgenic lines that integrated one of the constructs. CAT mRNA expression in transgenic tissues was assessed using the RNAse protection assay and the quantitative reverse transcribed polymerase chain reaction (RT-PCR). Overall CAT mRNA expression for both constructs was low compared to endogenous GLUT2 mRNA levels but the reporter transcript could be detected in all animals in the pancreatic islets and the liver, and in a few transgenic lines in the kidney and the small intestine. The CAT protein was also present in Langerhans islets and in the liver for both constructs by immunocytochemistry. These findings suggest that the proximal 338 bp of the murine GLUT2 promoter contain cis-elements required for the islet-specific expression of GLUT2.