Apical bulkheads accumulate as adaptive response to impaired bile flow in liver disease.

Hepatocytes form bile canaliculi that dynamically respond to the signalling activity of bile acids and bile flow. Little is known about their responses to intraluminal pressure. During embryonic development, hepatocytes assemble apical bulkheads that increase the canalicular resistance to intraluminal pressure. Here, we investigate whether they also protect bile canaliculi against elevated pressure upon impaired bile flow in adult liver. Apical bulkheads accumulate upon bile flow obstruction in mouse models and patients with primary sclerosing cholangitis (PSC). Their loss under these conditions leads to abnormally dilated canaliculi, resembling liver cell rosettes described in other hepatic diseases. 3D reconstruction reveals that these structures are sections of cysts and tubes formed by hepatocytes. Mathematical modelling establishes that they positively correlate with canalicular pressure and occur in early PSC stages. Using primary hepatocytes and 3D organoids, we demonstrate that excessive canalicular pressure causes the loss of apical bulkheads and formation of rosettes. Our results suggest that apical bulkheads are a protective mechanism of hepatocytes against impaired bile flow, highlighting the role of canalicular pressure in liver diseases.

Genetic variation in TERT modifies the risk of hepatocellular carcinoma in alcohol-related cirrhosis: results from a genome-wide case-control study.

OBJECTIVE: Hepatocellular carcinoma (HCC) often develops in patients with alcohol-related cirrhosis at an annual risk of up to 2.5%. Some host genetic risk factors have been identified but do not account for the majority of the variance in occurrence. This study aimed to identify novel susceptibility loci for the development of HCC in people with alcohol related cirrhosis. DESIGN: Patients with alcohol-related cirrhosis and HCC (cases: n=1214) and controls without HCC (n=1866), recruited from Germany, Austria, Switzerland, Italy and the UK, were included in a two-stage genome-wide association study using a case-control design. A validation cohort of 1520 people misusing alcohol but with no evidence of liver disease was included to control for possible association effects with alcohol misuse. Genotyping was performed using the InfiniumGlobal Screening Array (V.24v2, Illumina) and the OmniExpress Array (V.24v1-0a, Illumina). RESULTS: Associations with variants rs738409 in PNPLA3 and rs58542926 in TM6SF2 previously associated with an increased risk of HCC in patients with alcohol-related cirrhosis were confirmed at genome-wide significance. A novel locus rs2242652(A) in TERT (telomerase reverse transcriptase) was also associated with a decreased risk of HCC, in the combined meta-analysis, at genome-wide significance (p=6.41×10-9, OR=0.61 (95% CI 0.52 to 0.70). This protective association remained significant after correction for sex, age, body mass index and type 2 diabetes (p=7.94×10-5, OR=0.63 (95% CI 0.50 to 0.79). Carriage of rs2242652(A) in TERT was associated with an increased leucocyte telomere length (p=2.12×10-44). CONCLUSION: This study identifies rs2242652 in TERT as a novel protective factor for HCC in patients with alcohol-related cirrhosis.

The HSF1-CPT1a Pathway Is Differentially Regulated in NAFLD Progression.

Obesity and obesity-associated diseases represent one of the key health challenges of our time. In this context, aberrant hepatic lipid accumulation is a central pathological aspect of non-alcoholic fatty liver disease (NAFLD). By comparing methylation signatures of liver biopsies before and after bariatric surgery, we recently demonstrated the strong enrichment of differentially methylated heat shock factor 1 (HSF1) binding sites (>400-fold) in the process of liver remodeling, indicating a crucial role of HSF1 in modulating central aspects of NAFLD pathogenesis. Using cellular models of NAFLD, we were able to show that HSF1 is activated during fat accumulation in hepatocytes, mimicking conditions in patients before bariatric surgery. This induction was abolished by starving the cells, mimicking the situation after bariatric surgery. Regarding this connection, carnitine palmitoyltransferase 1 isoform A (CTP1a), a central regulator of lipid beta-oxidation, was identified as a HSF1 target gene by promoter analysis and HSF1 knockdown experiments. Finally, pharmacological activation of HSF1 through celastrol reduced fat accumulation in the cells in a HSF1-dependent manner. In conclusion, we were able to confirm the relevance of HSF1 activity and described a functional HSF1-CPT1a pathway in NAFLD pathogenesis.

Microbiota-dependent activation of the myeloid calcineurin-NFAT pathway inhibits B7H3- and B7H4-dependent anti-tumor immunity in colorectal cancer.

Bacterial sensing by intestinal tumor cells contributes to tumor growth through cell-intrinsic activation of the calcineurin-NFAT axis, but the role of this pathway in other intestinal cells remains unclear. Here, we found that myeloid-specific deletion of calcineurin in mice activated protective CD8+ T cell responses and inhibited colorectal cancer (CRC) growth. Microbial sensing by myeloid cells promoted calcineurin- and NFAT-dependent interleukin 6 (IL-6) release, expression of the co-inhibitory molecules B7H3 and B7H4 by tumor cells, and inhibition of CD8+ T cell-dependent anti-tumor immunity. Accordingly, targeting members of this pathway activated protective CD8+ T cell responses and inhibited primary and metastatic CRC growth. B7H3 and B7H4 were expressed by the majority of human primary CRCs and metastases, which was associated with low numbers of tumor-infiltrating CD8+ T cells and poor survival. Therefore, a microbiota-, calcineurin-, and B7H3/B7H4-dependent pathway controls anti-tumor immunity, revealing additional targets for immune checkpoint inhibition in microsatellite-stable CRC.

Sex-dependent dynamics of metabolism in primary mouse hepatocytes.

The liver is one of the most sexually dimorphic organs. The hepatic metabolic pathways that are subject to sexual dimorphism include xenobiotic, amino acid and lipid metabolism. Non-alcoholic fatty liver disease and hepatocellular carcinoma are among diseases with sex-dependent prevalence, progression and outcome. Although male and female livers differ in their abilities to metabolize foreign compounds, including drugs, sex-dependent treatment and pharmacological dynamics are rarely applied in all relevant cases. Therefore, it is important to consider hepatic sexual dimorphism when developing new treatment strategies and to understand the underlying mechanisms in model systems. We isolated primary hepatocytes from male and female C57BL6/N mice and examined the sex-dependent transcriptome, proteome and extracellular metabolome parameters in the course of culturing them for 96 h. The sex-specific gene expression of the general xenobiotic pathway altered and the female-specific expression of Cyp2b13 and Cyp2b9 was significantly reduced during culture. Sex-dependent differences of several signaling pathways increased, including genes related to serotonin and melatonin degradation. Furthermore, the ratios of male and female gene expression were inversed for other pathways, such as amino acid degradation, beta-oxidation, androgen signaling and hepatic steatosis. Because the primary hepatocytes were cultivated without the influence of known regulators of sexual dimorphism, these results suggest currently unknown modulatory mechanisms of sexual dimorphism in vitro. The large sex-dependent differences in the regulation and dynamics of drug metabolism observed during cultivation can have an immense influence on the evaluation of pharmacodynamic processes when conducting initial preclinical trials to investigate potential new drugs.

SETDB1 is required for intestinal epithelial differentiation and the prevention of intestinal inflammation.

OBJECTIVE: The intestinal epithelium is a rapidly renewing tissue which plays central roles in nutrient uptake, barrier function and the prevention of intestinal inflammation. Control of epithelial differentiation is essential to these processes and is dependent on cell type-specific activity of transcription factors which bind to accessible chromatin. Here, we studied the role of SET Domain Bifurcated Histone Lysine Methyltransferase 1, also known as ESET (SETDB1), a histone H3K9 methyltransferase, in intestinal epithelial homeostasis and IBD. DESIGN: We investigated mice with constitutive and inducible intestinal epithelial deletion of Setdb1, studied the expression of SETDB1 in patients with IBD and mouse models of IBD, and investigated the abundance of SETDB1 variants in healthy individuals and patients with IBD. RESULTS: Deletion of intestinal epithelial Setdb1 in mice was associated with defects in intestinal epithelial differentiation, barrier disruption, inflammation and mortality. Mechanistic studies showed that loss of SETDB1 leads to de-silencing of endogenous retroviruses, DNA damage and intestinal epithelial cell death. Predicted loss-of-function variants in human SETDB1 were considerably less frequently observed than expected, consistent with a critical role of SETDB1 in human biology. While the vast majority of patients with IBD showed unimpaired mucosal SETDB1 expression, comparison of IBD and non-IBD exomes revealed over-representation of individual rare missense variants in SETDB1 in IBD, some of which are predicted to be associated with loss of function and may contribute to the pathogenesis of intestinal inflammation. CONCLUSION: SETDB1 plays an essential role in intestinal epithelial homeostasis. Future work is required to investigate whether rare variants in SETDB1 contribute to the pathogenesis of IBD.

Loss of hepatic Mboat7 leads to liver fibrosis.

OBJECTIVE: The rs641738C>T variant located near the membrane-bound O-acyltransferase domain containing 7 (MBOAT7) locus is associated with fibrosis in liver diseases, including non-alcoholic fatty liver disease (NAFLD), alcohol-related liver disease, hepatitis B and C. We aim to understand the mechanism by which the rs641738C>T variant contributes to pathogenesis of NAFLD. DESIGN: Mice with hepatocyte-specific deletion of MBOAT7 (Mboat7Δhep) were generated and livers were characterised by histology, flow cytometry, qPCR, RNA sequencing and lipidomics. We analysed the association of rs641738C>T genotype with liver inflammation and fibrosis in 846 NAFLD patients and obtained genotype-specific liver lipidomes from 280 human biopsies. RESULTS: Allelic imbalance analysis of heterozygous human liver samples pointed to lower expression of the MBOAT7 transcript on the rs641738C>T haplotype. Mboat7Δhep mice showed spontaneous steatosis characterised by increased hepatic cholesterol ester content after 10 weeks. After 6 weeks on a high fat, methionine-low, choline-deficient diet, mice developed increased hepatic fibrosis as measured by picrosirius staining (p<0.05), hydroxyproline content (p<0.05) and transcriptomics, while the inflammatory cell populations and inflammatory mediators were minimally affected. In a human biopsied NAFLD cohort, MBOAT7 rs641738C>T was associated with fibrosis (p=0.004) independent of the presence of histological inflammation. Liver lipidomes of Mboat7Δhep mice and human rs641738TT carriers with fibrosis showed increased total lysophosphatidylinositol levels. The altered lysophosphatidylinositol and phosphatidylinositol subspecies in MBOAT7Δhep livers and human rs641738TT carriers were similar. CONCLUSION: Mboat7 deficiency in mice and human points to an inflammation-independent pathway of liver fibrosis that may be mediated by lipid signalling and a potentially targetable treatment option in NAFLD.

Shotgun lipidomics-based characterization of the landscape of lipid metabolism in colorectal cancer.

Solid tumors are characterized by global metabolic alterations which contribute to their growth and progression. Altered gene expression profiles and plasma lipid composition suggested a role for metabolic reprogramming in colorectal cancer (CRC) development. However, a conclusive picture of CRC-associated lipidome alterations in the tumor tissue has not emerged. Here, we determined molar abundances of 342 species from 20 lipid classes in matched biopsies of CRC and adjacent normal mucosa. We demonstrate that in contrast to previous reports, CRC shows a largely preserved lipidome composition that resembles that of normal colonic mucosa. Important exceptions include increased levels of lyso-phosphatidylinositols in CRC and reduced abundance of ether phospholipids in advanced stages of CRC. As such, our observations challenge the concept of widespread alterations in lipid metabolism in CRC and rather suggest changes in the cellular lipid profile that are limited to selected lipids involved in signaling and the scavenging of reactive oxygen species.

Heterozygous carriage of the alpha1-antitrypsin Pi*Z variant increases the risk to develop liver cirrhosis.

OBJECTIVE: Homozygous alpha1-antitrypsin (AAT) deficiency increases the risk for developing cirrhosis, whereas the relevance of heterozygous carriage remains unclear. Hence, we evaluated the impact of the two most relevant AAT variants ('Pi*Z' and 'Pi*S'), present in up to 10% of Caucasians, on subjects with non-alcoholic fatty liver disease (NAFLD) or alcohol misuse. DESIGN: We analysed multicentric case-control cohorts consisting of 1184 people with biopsy-proven NAFLD and of 2462 people with chronic alcohol misuse, both cohorts comprising cases with cirrhosis and controls without cirrhosis. Genotyping for the Pi*Z and Pi*S variants was performed. RESULTS: The Pi*Z variant presented in 13.8% of patients with cirrhotic NAFLD but only in 2.4% of counterparts without liver fibrosis (p<0.0001). Accordingly, the Pi*Z variant increased the risk of NAFLD subjects to develop cirrhosis (adjusted OR=7.3 (95% CI 2.2 to 24.8)). Likewise, the Pi*Z variant presented in 6.2% of alcohol misusers with cirrhosis but only in 2.2% of alcohol misusers without significant liver injury (p<0.0001). Correspondingly, alcohol misusers carrying the Pi*Z variant were prone to develop cirrhosis (adjusted OR=5.8 (95% CI 2.9 to 11.7)). In contrast, the Pi*S variant was not associated with NAFLD-related cirrhosis and only borderline with alcohol-related cirrhosis (adjusted OR=1.47 (95% CI 0.99 to 2.19)). CONCLUSION: The Pi*Z variant is the hitherto strongest single nucleotide polymorphism-based risk factor for cirrhosis in NAFLD and alcohol misuse, whereas the Pi*S variant confers only a weak risk in alcohol misusers. As 2%-4% of Caucasians are Pi*Z carriers, this finding should be considered in genetic counselling of affected individuals.

Aberrant DNA methylation of ADAMTS16 in colorectal and other epithelial cancers.

BACKGROUND: ADAMs (a disintegrin and metalloproteinase) have long been associated with tumor progression. Recent findings indicate that members of the closely related ADAMTS (ADAMs with thrombospondin motifs) family are also critically involved in carcinogenesis. Gene silencing through DNA methylation at CpG loci around e.g. transcription start or enhancer sites is a major mechanism in cancer development. Here, we aimed at identifying genes of the ADAM and ADAMTS family showing altered DNA methylation in the development or colorectal cancer (CRC) and other epithelial tumors. METHODS: We investigated potential changes of DNA methylation affecting ADAM and ADAMTS genes in 117 CRC, 40 lung cancer (LC) and 15 oral squamous-cell carcinoma (SCC) samples. Tumor tissue was analyzed in comparison to adjacent non-malignant tissue of the same patients. The methylation status of 1145 CpGs in 51 ADAM and ADAMTS genes was measured with the HumanMethylation450 BeadChip Array. ADAMTS16 protein expression was analyzed in CRC samples by immunohistochemistry. RESULTS: In CRC, we identified 72 CpGs in 18 genes which were significantly affected by hyper- or hypomethylation in the tumor tissue compared to the adjacent non-malignant tissue. While notable/frequent alterations in methylation patterns within ADAM genes were not observed, conspicuous changes were found in ADAMTS16 and ADAMTS2. To figure out whether these differences would be CRC specific, additional LC and SCC tissue samples were analyzed. Overall, 78 differentially methylated CpGs were found in LC and 29 in SCC. Strikingly, 8 CpGs located in the ADAMTS16 gene were commonly differentially methylated in all three cancer entities. Six CpGs in the promoter region were hypermethylated, whereas 2 CpGs in the gene body were hypomethylated indicative of gene silencing. In line with these findings, ADAMTS16 protein was strongly expressed in globlet cells and colonocytes in control tissue but not in CRC samples. Functional in vitro studies using the colorectal carcinoma cell line HT29 revealed that ADAMTS16 expression restrained tumor cell proliferation. CONCLUSIONS: We identified ADAMTS16 as novel gene with cancer-specific promoter hypermethylation in CRC, LC and SCC patients implicating ADAMTS16 as potential biomarker for these tumors. Moreover, our results provide evidence that ADAMTS16 may have tumor suppressor properties.

Evolutionary Distance Predicts Recurrence After Liver Transplantation in Multifocal Hepatocellular Carcinoma.

BACKGROUND: Liver transplantation (LTx) is a potentially curative treatment option for hepatocellular carcinoma (HCC) in cirrhosis. However, patients, where HCC is already a systemic disease, LTx may be individually harmful and has a negative impact on donor organ usage. Thus, there is a need for improved selection criteria beyond nodule morphology to select patients with a favorable outcome for LTx in multifocal HCC. Evolutionary distance measured from genome-wide single-nucleotide polymorphism data between tumor nodules and the cirrhotic liver may be a prognostic marker of survival after LTx for multifocal HCC. METHODS: In a retrospective multicenter study, clinical data and formalin-fixed paraffin-embedded specimens of the liver and 2 tumor nodules were obtained from explants of 30 patients in the discovery and 180 patients in the replication cohort. DNA was extracted from formalin-fixed paraffin-embedded specimens followed by genome wide single-nucleotide polymorphism genotyping. RESULTS: Genotype quality criteria allowed for analysis of 8 patients in the discovery and 17 patients in the replication set. DNA concentrations of a total of 25 patients fulfilled the quality criteria and were included in the analysis. Both, in the discovery (P = 0.04) and in the replication data sets (P = 0.01), evolutionary distance was associated with the risk of recurrence of HCC after transplantation (combined P = 0.0002). In a univariate analysis, evolutionary distance (P = 7.4 × 10) and microvascular invasion (P = 1.31 × 10) were significantly associated with survival in a Cox regression analysis. CONCLUSIONS: Evolutionary distance allows for the determination of a high-risk group of recurrence if preoperative liver biopsy is considered.

Comparison of Gene Expression Patterns Between Mouse Models of Nonalcoholic Fatty Liver Disease and Liver Tissues From Patients.

BACKGROUND & AIMS: Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disorder in industrialized countries. Mouse models of NAFLD have been used in studies of pathogenesis and treatment, and have certain features of the human disease. We performed a systematic transcriptome-wide analysis of liver tissues from patients at different stages of NAFLD progression (ranging from healthy obese individuals to those with steatosis), as well as rodent models of NAFLD, to identify those that most closely resemble human disease progression in terms of gene expression patterns. METHODS: We performed a systematic evaluation of genome-wide messenger RNA expression using liver tissues collected from mice fed a standard chow diet (controls) and 9 mouse models of NAFLD: mice on a high-fat diet (with or without fructose), mice on a Western-type diet, mice on a methionine- and choline-deficient diet, mice on a high-fat diet given streptozotocin, and mice with disruption of Pten in hepatocytes. We compared gene expression patterns with those of liver tissues from 25 patients with nonalcoholic steatohepatitis (NASH), 27 patients with NAFLD, 15 healthy obese individuals, and 39 healthy nonobese individuals (controls). Liver samples were obtained from patients undergoing liver biopsy for suspected NAFLD or NASH, or during liver or bariatric surgeries. Data sets were analyzed using the limma R-package. Overlap of functional profiles was analyzed by gene set enrichment analysis profiles. RESULTS: We found differences between human and mouse transcriptomes to be significantly larger than differences between disease stages or models. Of the 65 genes with significantly altered expression in patients with NASH and 177 genes with significantly altered expression in patients with NAFLD, compared with controls, only 1-18 of these genes also differed significantly in expression between mouse models of NAFLD and control mice. However, expression of genes that regulate pathways associated with the development of NAFLD were altered in some mouse models (such as pathways associated with lipid metabolism). On a pathway level, gene expression patterns in livers of mice on the high-fat diet were associated more closely with human fatty liver disease than other models. CONCLUSIONS: In comparing gene expression profiles between liver tissues from different mouse models of NAFLD and patients with different stages of NAFLD, we found very little overlap. Our data set is available for studies of pathways that contribute to the development of NASH and NAFLD and selection of the most applicable mouse models (http://www.nash-profiler.com).

Metabolomic tissue signature in human non-alcoholic fatty liver disease identifies protective candidate metabolites.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver disorder in industrialized countries, yet its pathophysiology is incompletely understood. Small-molecule metabolite screens may offer new insights into disease mechanisms and reveal new treatment targets. METHODS: Discovery (N = 33) and replication (N = 66) of liver biopsies spanning the range from normal liver histology to non-alcoholic steatohepatitis (NASH) were ascertained ensuring rapid freezing under 30 s in patients. 252 metabolites were assessed using GC/MS. Replicated metabolites were evaluated in a murine high-fat diet model of NAFLD. RESULTS: In a two-stage metabolic screening, hydroquinone (HQ, p(combined) = 3.0 × 10(-4)) and nicotinic acid (NA, p(combined) = 3.9 × 10(-9)) were inversely correlated with histological NAFLD severity. A murine high-fat diet model of NAFLD demonstrated a protective effect of these two substances against NAFLD: Supplementation with 1% HQ reduced only liver steatosis, whereas 0.6% NA reduced both liver fat content and serum transaminase levels and induced a complex regulatory network of genes linked to NALFD pathogenesis in a global expression pathway analysis. Human nutritional intake of NA equivalent was also consistent with a protective effect of NA against NASH progression. CONCLUSION: This first small-molecular screen of human liver tissue identified two replicated protective metabolites. Either the use of NA or targeting its regulatory pathways might be explored to treat or prevent human NAFLD.

Obesity accelerates epigenetic aging of human liver.

Because of the dearth of biomarkers of aging, it has been difficult to test the hypothesis that obesity increases tissue age. Here we use a novel epigenetic biomarker of aging (referred to as an "epigenetic clock") to study the relationship between high body mass index (BMI) and the DNA methylation ages of human blood, liver, muscle, and adipose tissue. A significant correlation between BMI and epigenetic age acceleration could only be observed for liver (r = 0.42, P = 6.8 × 10(-4) in dataset 1 and r = 0.42, P = 1.2 × 10(-4) in dataset 2). On average, epigenetic age increased by 3.3 y for each 10 BMI units. The detected age acceleration in liver is not associated with the Nonalcoholic Fatty Liver Disease Activity Score or any of its component traits after adjustment for BMI. The 279 genes that are underexpressed in older liver samples are highly enriched (1.2 × 10(-9)) with nuclear mitochondrial genes that play a role in oxidative phosphorylation and electron transport. The epigenetic age acceleration, which is not reversible in the short term after rapid weight loss induced by bariatric surgery, may play a role in liver-related comorbidities of obesity, such as insulin resistance and liver cancer.

B Lymphocyte Stimulator (BLyS) is expressed in human adipocytes in vivo and is related to obesity but not to insulin resistance.

Inflammation and metabolism have been shown to be evolutionary linked and increasing evidence exists that pro-inflammatory factors are involved in the pathogenesis of obesity and type 2 diabetes. Until now, most data suggest that within adipose tissue these factors are secreted by cells of the innate immune system, e. g. macrophages. In the present study we demonstrate that B lymphocyte stimulator (BLyS) is increased in human obesity. In contrast to several pro-inflammatory factors, we found the source of BLyS in human adipose tissue to be the adipocytes rather than immune cells. In grade 3 obese human subjects, expression of BLyS in vivo in adipose tissue is significantly increased (p<0.001). Furthermore, BLyS serum levels are elevated in grade 3 human obesity (862.5+222.0 pg/ml vs. 543.7+60.7 pg/ml in lean controls, p<0.001) and are positively correlated to the BMI (r = 0.43, p<0.0002). In the present study, bariatric surgery significantly altered serum BLyS concentrations. In contrast, weight loss due to a very-low-calorie-formula-diet (800 kcal/d) had no such effect. To examine metabolic activity of BLyS, in a translational research approach, insulin sensitivity was measured in human subjects in vivo before and after treatment with the human recombinant anti-BLyS antibody belimumab. Since BLyS is known to promote B-cell proliferation and immunoglobulin secretion, the present data suggest that adipocytes of grade 3 obese human subjects are able to activate the adaptive immune system, suggesting that in metabolic inflammation in humans both, innate and adaptive immunity, are of pathophysiological relevance.

DNA methylation analysis in nonalcoholic fatty liver disease suggests distinct disease-specific and remodeling signatures after bariatric surgery.

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disorder in industrialized countries. Liver samples from morbidly obese patients (n = 45) with all stages of NAFLD and controls (n = 18) were analyzed by array-based DNA methylation and mRNA expression profiling. NAFLD-specific expression and methylation differences were seen for nine genes coding for key enzymes in intermediate metabolism (including PC, ACLY, and PLCG1) and insulin/insulin-like signaling (including IGF1, IGFBP2, and PRKCE) and replicated by bisulfite pyrosequening (independent n = 39). Transcription factor binding sites at NAFLD-specific CpG sites were >1,000-fold enriched for ZNF274, PGC1A, and SREBP2. Intraindividual comparison of liver biopsies before and after bariatric surgery showed NAFLD-associated methylation changes to be partially reversible. Postbariatric and NAFLD-specific methylation signatures were clearly distinct both in gene ontology and transcription factor binding site analyses, with >400-fold enrichment of NRF1, HSF1, and ESRRA sites. Our findings provide an example of treatment-induced epigenetic organ remodeling in humans.

Distinct DNA methylation patterns in cirrhotic liver and hepatocellular carcinoma.

Abberrant DNA methylation is one of the hallmarks of cancerogenesis. Our study aims to delineate differential DNA methylation in cirrhosis and hepatic cancerogenesis. Patterns of methylation of 27,578 individual CpG loci in 12 hepatocellular carcinomas (HCCs), 15 cirrhotic controls and 12 normal liver samples were investigated using an array-based technology. A supervised principal component analysis (PCA) revealed 167 hypomethylated loci and 100 hypermethylated loci in cirrhosis and HCC as compared to normal controls. Thus, these loci show a "cirrhotic" methylation pattern that is maintained in HCC. In pairwise supervised PCAs between normal liver, cirrhosis and HCC, eight loci were significantly changed in all analyses differentiating the three groups (p < 0.0001). Of these, five loci showed highest methylation levels in HCC and lowest in control tissue (LOC55908, CELSR1, CRMP1, GNRH2, ALOX12 and ANGPTL7), whereas two loci showed the opposite direction of change (SPRR3 and TNFSF15). Genes hypermethylated between normal liver to cirrhosis, which maintain this methylation pattern during the development of HCC, are depleted for CpG islands, high CpG content promoters and polycomb repressive complex 2 (PRC2) targets in embryonic stem cells. In contrast, genes selectively hypermethylated in HCC as compared to nonmalignant samples showed an enrichment of CpG islands, high CpG content promoters and PRC2 target genes (p < 0.0001). Cirrhosis and HCC show distinct patterns of differential methylation with regards to promoter structure, PRC2 targets and CpG islands.

Statistical inference of allelic imbalance from transcriptome data.

Next-generation sequencing and the availability of high-density genotyping arrays have facilitated an analysis of somatic and meiotic mutations at unprecedented level, but drawing sensible conclusions about the functional relevance of the detected variants still remains a formidable challenge. In this context, the study of allelic imbalance in intermediate RNA phenotypes may prove a useful means to elucidate the likely effects of DNA variants of unknown significance. We developed a statistical framework for the assessment of allelic imbalance in next-generation transcriptome sequencing (RNA-seq) data that requires neither an expression reference nor the underlying nuclear genotype(s), and that allows for allele miscalls. Using extensive simulation as well as publicly available whole-transcriptome data from European-descent individuals in HapMap, we explored the power of our approach in terms of both genotype inference and allelic imbalance assessment under a wide range of practically relevant scenarios. In so doing, we verified a superior performance of our methodology, particularly at low sequencing coverage, compared to the more simplistic approach of completely ignoring allele miscalls. Because the proposed framework can be used to assess somatic mutations and allelic imbalance in one and the same set of RNA-seq data, it will be particularly useful for the analysis of somatic genetic variation in cancer studies.

Loci from a genome-wide analysis of bilirubin levels are associated with gallstone risk and composition.

BACKGROUND & AIMS: Genome-wide association studies have mapped loci that are associated with serum levels of bilirubin. Bilirubin is a major component of gallstones so we investigated whether these variants predict gallstone bilirubin content and overall risk for gallstones. METHODS: Loci that were identified in a meta-analysis to attain a genome-wide significance level of a P value less than 1.0×10(-7) (UGT1A1, SLCO1B1, LST-3TM12, SLCO1A2) were analyzed in 1018 individuals with known gallstone composition. Gallstone risk was analyzed in 2606 German choleystecomized individuals and 1121 controls and was replicated in 210 cases and 496 controls from South America. RESULTS: By using the presence of bilirubin as a phenotype, variants rs6742078 (UGT1A1; P = .003), rs4149056 (SLCO1B1; P = .003), and rs4149000 (SLCO1A2; P = .015) were associated with gallstone composition. In regression analyses, only UGT1A1 and SLCO1B1 were independently retained in the model. UGT1A1 (rs6742078; P = .018) was associated with overall gallstone risk. In a sex-stratified analysis, only male carriers of rs6742078 had an increased risk for gallstone disease (P = 2.1×10(-7); odds ratio(recessive), 2.34; P(women) = .47). The sex-specific association of rs6742078 was confirmed in samples from South America (P(men) = .046; odds ratio(recessive), 2.19; P(women) = .96). CONCLUSIONS: The UGT1A1 Gilbert syndrome variant rs6742078 is associated with gallstone disease in men; further studies are required regarding the sex-specific physiology of bilirubin and bile acid metabolism. Variants of ABCG8 and UGT1A1 are the 2 major risk factors for overall gallstone disease, they contribute a population attributable risk of 21.2% among men.