High serum albumin ascites gradient ascites--an atypical presentation of metastatic pancreatic cancer.

Pancreatic adenocarcinoma has less than a 5% 5-year survival rate, and metastatic disease is associated with a median survival of 4.5 months. A typical presentation often includes evidence of biliary obstruction, abdominal pain, jaundice, and weight loss. Significant ascites is not commonly seen at initial presentation and, when present, is typically associated with a low serum albumin ascites gradient (SAAG). We discuss a patient who presented with high-SAAG ascites as her initial presentation, only to be later diagnosed with metastatic pancreatic adenocarcinoma.

Dynamic chromatin states in human ES cells reveal potential regulatory sequences and genes involved in pluripotency.

Pluripotency, the ability of a cell to differentiate and give rise to all embryonic lineages, defines a small number of mammalian cell types such as embryonic stem (ES) cells. While it has been generally held that pluripotency is the product of a transcriptional regulatory network that activates and maintains the expression of key stem cell genes, accumulating evidence is pointing to a critical role for epigenetic processes in establishing and safeguarding the pluripotency of ES cells, as well as maintaining the identity of differentiated cell types. In order to better understand the role of epigenetic mechanisms in pluripotency, we have examined the dynamics of chromatin modifications genome-wide in human ES cells (hESCs) undergoing differentiation into a mesendodermal lineage. We found that chromatin modifications at promoters remain largely invariant during differentiation, except at a small number of promoters where a dynamic switch between acetylation and methylation at H3K27 marks the transition between activation and silencing of gene expression, suggesting a hierarchy in cell fate commitment over most differentially expressed genes. We also mapped over 50 000 potential enhancers, and observed much greater dynamics in chromatin modifications, especially H3K4me1 and H3K27ac, which correlate with expression of their potential target genes. Further analysis of these enhancers revealed potentially key transcriptional regulators of pluripotency and a chromatin signature indicative of a poised state that may confer developmental competence in hESCs. Our results provide new evidence supporting the role of chromatin modifications in defining enhancers and pluripotency.

ChIP-chip designs to interrogate the genome of Xenopus embryos for transcription factor binding and epigenetic regulation.

BACKGROUND: Chromatin immunoprecipitation combined with genome tile path microarrays or deep sequencing can be used to study genome-wide epigenetic profiles and the transcription factor binding repertoire. Although well studied in a variety of cell lines, these genome-wide profiles have so far been little explored in vertebrate embryos. PRINCIPAL FINDINGS: Here we report on two genome tile path ChIP-chip designs for interrogating the Xenopus tropicalis genome. In particular, a whole-genome microarray design was used to identify active promoters by close proximity to histone H3 lysine 4 trimethylation. A second microarray design features these experimentally derived promoter regions in addition to currently annotated 5' ends of genes. These regions truly represent promoters as shown by binding of TBP, a key transcription initiation factor. CONCLUSIONS: A whole-genome and a promoter tile path microarray design was developed. Both designs can be used to study epigenetic phenomena and transcription factor binding in developing Xenopus embryos.

Dynamic histone H3 epigenome marking during the intraerythrocytic cycle of Plasmodium falciparum.

Epigenome profiling has led to the paradigm that promoters of active genes are decorated with H3K4me3 and H3K9ac marks. To explore the epigenome of Plasmodium falciparum asexual stages, we performed MS analysis of histone modifications and found a general preponderance of H3/H4 acetylation and H3K4me3. ChIP-on-chip profiling of H3, H3K4me3, H3K9me3, and H3K9ac from asynchronous parasites revealed an extensively euchromatic epigenome with heterochromatin restricted to variant surface antigen gene families (VSA) and a number of genes hitherto unlinked to VSA. Remarkably, the vast majority of the genome shows an unexpected pattern of enrichment of H3K4me3 and H3K9ac. Analysis of synchronized parasites revealed significant developmental stage specificity of the epigenome. In rings, H3K4me3 and H3K9ac are homogenous across the genes marking active and inactive genes equally, whereas in schizonts, they are enriched at the 5' end of active genes. This study reveals an unforeseen and unique plasticity in the use of the epigenetic marks and implies the presence of distinct epigenetic pathways in gene silencing/activation throughout the erythrocytic cycle.

High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers.

Many genome-wide assays involve the generation of a subset (or representation) of the genome following restriction enzyme digestion. The use of enzymes sensitive to cytosine methylation allows high-throughput analysis of this epigenetic regulatory process. We show that the use of a dual-adapter approach allows us to generate genomic representations that includes fragments of <200 bp in size, previously not possible when using the standard approach of using a single adapter. By expanding the representation to smaller fragments using HpaII or MspI, we increase the representation by these isoschizomers to more than 1.32 million loci in the human genome, representing 98.5% of CpG islands and 91.1% of refSeq promoters. This advance allows the development of a new, high-resolution version of our HpaII-tiny fragment Enrichment by Ligation-mediated PCR (HELP) assay to study cytosine methylation. We also show that the MspI representation generates information about copy-number variation, that the assay can be used on as little as 10 ng of DNA and that massively parallel sequencing can be used as an alternative to microarrays to read the output of the assay, making this a powerful discovery platform for studies of genomic and epigenomic abnormalities.

Histone modifications at human enhancers reflect global cell-type-specific gene expression.

The human body is composed of diverse cell types with distinct functions. Although it is known that lineage specification depends on cell-specific gene expression, which in turn is driven by promoters, enhancers, insulators and other cis-regulatory DNA sequences for each gene, the relative roles of these regulatory elements in this process are not clear. We have previously developed a chromatin-immunoprecipitation-based microarray method (ChIP-chip) to locate promoters, enhancers and insulators in the human genome. Here we use the same approach to identify these elements in multiple cell types and investigate their roles in cell-type-specific gene expression. We observed that the chromatin state at promoters and CTCF-binding at insulators is largely invariant across diverse cell types. In contrast, enhancers are marked with highly cell-type-specific histone modification patterns, strongly correlate to cell-type-specific gene expression programs on a global scale, and are functionally active in a cell-type-specific manner. Our results define over 55,000 potential transcriptional enhancers in the human genome, significantly expanding the current catalogue of human enhancers and highlighting the role of these elements in cell-type-specific gene expression.

Distinct transcriptional outputs associated with mono- and dimethylated histone H3 arginine 2.

Dimethylation of histone H3 Arg2 (H3R2me2) maintains transcriptional silencing by inhibiting Set1 mediated trimethylation of H3K4. Here we demonstrate that Arg2 is also monomethylated (H3R2me1) in yeast but that its functional characteristics are distinct from H3R2me2: (i) H3R2me1 does not inhibit histone H3 Lys4 (H3K4) methylation; (ii) it is present throughout the coding region of genes; and (iii) it correlates with active transcription. Collectively, these results indicate that different H3R2 methylation states have defined roles in gene expression.

Chromosome 12p deletions in TEL-AML1 childhood acute lymphoblastic leukemia are associated with retrotransposon elements and occur postnatally.

TEL-AML1 (ETV6-RUNX1) is the most common translocation in the childhood leukemias, and is a prenatal mutation in most children. This translocation has been detected at a high rate among newborns ( approximately 1%); therefore, the rate-limiting event for leukemia seems to be secondary mutations. One such frequent mutation in this subtype is partial deletion of chromosome 12p, trans from the translocation. Nine del(12p) breakpoints within six leukemia cases were sequenced to explore the etiology of this genetic event, and most involved cryptic sterile translocations. Twelve of 18 del(12p) parent sequences involved in these breakpoints were located in repeat regions (8 of these in long interspersed nuclear elements). This stands in contrast with TEL-AML1, in which only 21 of 110 previously assessed breakpoints (19%) occur in DNA repeats (P=0.0001). An exploratory assessment of archived neonatal blood cards revealed significantly more long interspersed nuclear element CpG methylations in individuals at birth who were later diagnosed with TEL-AML1 leukemia, compared with individuals who did not contract leukemia (P=0.01). Nontemplate nucleotides were also more frequent in del(12p) than in TEL-AML1 junctions (P=0.004), suggesting formation by terminal deoxynucleotidyl transferase. Assessment of six archived neonatal blood cards indicated that no del(12p) rearrangements backtracked to birth, although two of these patients were previously positive for TEL-AML1 using the same assay with comparable sensitivity. These data are compatible with a two-stage natural history: TEL-AML1 occurs prenatally, and del(12p) occurs postnatally in more mature cells with a structure that suggests the involvement of retrotransposon instability.

Overlapping euchromatin/heterochromatin- associated marks are enriched in imprinted gene regions and predict allele-specific modification.

Most genome-level analysis treats the two parental alleles equivalently, yet diploid genomes contain two parental genomes that are often epigenetically distinct. While single nucleotide polymorphisms (SNPs) can be used to distinguish these genomes, it would be useful to develop a generalized strategy for identifying candidate genes or regions showing allele-specific differences, independent of SNPs. We have explored this problem by looking for overlapping marks in the genome related to both euchromatin (histone H3 dimethyl lysine-4 [H3K4Me2]) and heterochromatin (DNA methylation [DNAm]). "Double hits" were defined by the intersection of H3K4Me2 and DNAm. For the top 5% of marks, defined by a sliding window, imprinted gene regions were enriched for double hits 5.4-fold. When the location information of CTCF binding sites were integrated, the "triple hits" were enriched 76-fold for known imprinted genes in the regions studied. The double hits in imprinted genes were found to occur usually at the site of alternative or antisense transcripts. In addition, four of four imprinted genes tested showing double hits also showed allele-specific methylation. We suggest that overlapping euchromatin/heterochromatin marks are common and are enriched for epigenetically distinct parental chromosome regions. Furthermore, we developed a novel approach to identifying allele-specific marks that is SNP independent, by fractionating using H3K4Me2 antibodies followed by DNA methylation analysis.

The genomic distribution and function of histone variant HTZ-1 during C. elegans embryogenesis.

In all eukaryotes, histone variants are incorporated into a subset of nucleosomes to create functionally specialized regions of chromatin. One such variant, H2A.Z, replaces histone H2A and is required for development and viability in all animals tested to date. However, the function of H2A.Z in development remains unclear. Here, we use ChIP-chip, genetic mutation, RNAi, and immunofluorescence microscopy to interrogate the function of H2A.Z (HTZ-1) during embryogenesis in Caenorhabditis elegans, a key model of metazoan development. We find that HTZ-1 is expressed in every cell of the developing embryo and is essential for normal development. The sites of HTZ-1 incorporation during embryogenesis reveal a genome wrought by developmental processes. HTZ-1 is incorporated upstream of 23% of C. elegans genes. While these genes tend to be required for development and occupied by RNA polymerase II, HTZ-1 incorporation does not specify a stereotypic transcription program. The data also provide evidence for unexpectedly widespread independent regulation of genes within operons during development; in 37% of operons, HTZ-1 is incorporated upstream of internally encoded genes. Fewer sites of HTZ-1 incorporation occur on the X chromosome relative to autosomes, which our data suggest is due to a paucity of developmentally important genes on X, rather than a direct function for HTZ-1 in dosage compensation. Our experiments indicate that HTZ-1 functions in establishing or maintaining an essential chromatin state at promoters regulated dynamically during C. elegans embryogenesis.

Transcriptional enhancement by GATA1-occupied DNA segments is strongly associated with evolutionary constraint on the binding site motif.

Tissue development and function are exquisitely dependent on proper regulation of gene expression, but it remains controversial whether the genomic signals controlling this process are subject to strong selective constraint. While some studies show that highly constrained noncoding regions act to enhance transcription, other studies show that DNA segments with biochemical signatures of regulatory regions, such as occupancy by a transcription factor, are seemingly unconstrained across mammalian evolution. To test the possible correlation of selective constraint with enhancer activity, we used chromatin immunoprecipitation as an approach unbiased by either evolutionary constraint or prior knowledge of regulatory activity to identify DNA segments within a 66-Mb region of mouse chromosome 7 that are occupied by the erythroid transcription factor GATA1. DNA segments bound by GATA1 were identified by hybridization to high-density tiling arrays, validated by quantitative PCR, and tested for gene regulatory activity in erythroid cells. Whereas almost all of the occupied segments contain canonical WGATAR binding site motifs for GATA1, in only 45% of the cases is the motif deeply preserved (found at the orthologous position in placental mammals or more distant species). However, GATA1-bound segments with high enhancer activity tend to be the ones with an evolutionarily preserved WGATAR motif, and this relationship was confirmed by a loss-of-function assay. Thus, GATA1 binding sites that regulate gene expression during erythroid maturation are under strong selective constraint, while nonconstrained binding may have only a limited or indirect role in regulation.

Analysis of the mechanisms mediating tumor-specific changes in gene expression in human liver tumors.

There is widespread interest in efficient characterization of differences between tumor and normal samples. Here, we show an effective methodology for genome-scale characterization of tumors. Using matched normal and tumor samples from liver cancer patients, as well as non-cancer-related normal liver tissue, we first determined changes in gene expression as monitored on RNA expression arrays. We identified several hundred mRNAs that were consistently changed in the tumor samples. To characterize the mechanisms responsible for creation of the tumor-specific transcriptome, we performed chromatin immunoprecipitation on microarray experiments to assay binding of RNA polymerase II, H3me3K27, and H3me3K9 and DNA methylation in 25,000 promoter regions. These experiments identified changes in active and silenced regions of the genome in the tumor cells. Finally, we used a "virtual comparative genomic hybridization" method to identify copy number alterations in the tumor samples. Through comparison of RNA polymerase II binding, chromatin structure, DNA methylation, and copy number changes, we suggest that the major contributor to creation of the liver tumor transcriptome was changes in gene copy number.

Using ChIP-chip technology to reveal common principles of transcriptional repression in normal and cancer cells.

We compared 12 different cell populations, including embryonic stem cells before and during differentiation into embryoid bodies as well as various types of normal and tumor cells to determine if pluripotent versus differentiated cell types use different mechanisms to establish their transcriptome. We first identified genes that were not expressed in the 12 different cell populations and then determined which of them were regulated by histone methylation, DNA methylation, at the step of productive elongation, or by the inability to establish a preinitiation complex. For these experiments, we performed chromatin immunoprecipitation using antibodies to H3me3K27, H3me3K9, 5-methyl-cytosine, and POLR2A. We found that (1) the percentage of low expressed genes bound by POLR2A, H3me3K27, H3me3K9, or 5-methyl-cytosine is similar in all 12 cell types, regardless of differentiation or neoplastic state; (2) a gene is generally repressed by only one mechanism; and (3) distinct classes of genes are repressed by certain mechanisms. We further characterized two transitioning cell populations, 3T3 cells progressing from G0/G1 into S phase and mES cells differentiating into embryoid bodies. We found that the transient regulation through the cell cycle was achieved predominantly by changes in the recruitment of the general transcriptional machinery or by post-POLR2A recruitment mechanisms. In contrast, changes in chromatin silencing were critical for the permanent changes in gene expression in cells undergoing differentiation.

An integrative genomic and epigenomic approach for the study of transcriptional regulation.

The molecular heterogeneity of acute leukemias and other tumors constitutes a major obstacle towards understanding disease pathogenesis and developing new targeted-therapies. Aberrant gene regulation is a hallmark of cancer and plays a central role in determining tumor phenotype. We predicted that integration of different genome-wide epigenetic regulatory marks along with gene expression levels would provide greater power in capturing biological differences between leukemia subtypes. Gene expression, cytosine methylation and histone H3 lysine 9 (H3K9) acetylation were measured using high-density oligonucleotide microarrays in primary human acute myeloid leukemia (AML) and acute lymphocytic leukemia (ALL) specimens. We found that DNA methylation and H3K9 acetylation distinguished these leukemias of distinct cell lineage, as expected, but that an integrative analysis combining the information from each platform revealed hundreds of additional differentially expressed genes that were missed by gene expression arrays alone. This integrated analysis also enhanced the detection and statistical significance of biological pathways dysregulated in AML and ALL. Integrative epigenomic studies are thus feasible using clinical samples and provide superior detection of aberrant transcriptional programming than single-platform microarray studies.

SNP-specific array-based allele-specific expression analysis.

We have developed an optimized array-based approach for customizable allele-specific gene expression (ASE) analysis. The central features of the approach are the ability to select SNPs at will for detection, and the absence of need to PCR amplify the target. A surprisingly long probe length (39-49 nt) was needed for allelic discrimination. Reconstitution experiments demonstrate linearity of ASE over a broad range. Using this approach, we have discovered at least two novel imprinted genes, NLRP2, which encodes a member of the inflammasome, and OSBPL1A, which encodes a presumed oxysterol-binding protein, were both preferentially expressed from the maternal allele. In contrast, ERAP2, which encodes an aminopeptidase, did not show preferential parent-of-origin expression, but rather, cis-acting nonimprinted differential allelic control. The approach is scalable to the whole genome and can be used for discovery of functional epigenetic modifications in patient samples.

Systematic evaluation of variability in ChIP-chip experiments using predefined DNA targets.

The most widely used method for detecting genome-wide protein-DNA interactions is chromatin immunoprecipitation on tiling microarrays, commonly known as ChIP-chip. Here, we conducted the first objective analysis of tiling array platforms, amplification procedures, and signal detection algorithms in a simulated ChIP-chip experiment. Mixtures of human genomic DNA and "spike-ins" comprised of nearly 100 human sequences at various concentrations were hybridized to four tiling array platforms by eight independent groups. Blind to the number of spike-ins, their locations, and the range of concentrations, each group made predictions of the spike-in locations. We found that microarray platform choice is not the primary determinant of overall performance. In fact, variation in performance between labs, protocols, and algorithms within the same array platform was greater than the variation in performance between array platforms. However, each array platform had unique performance characteristics that varied with tiling resolution and the number of replicates, which have implications for cost versus detection power. Long oligonucleotide arrays were slightly more sensitive at detecting very low enrichment. On all platforms, simple sequence repeats and genome redundancy tended to result in false positives. LM-PCR and WGA, the most popular sample amplification techniques, reproduced relative enrichment levels with high fidelity. Performance among signal detection algorithms was heavily dependent on array platform. The spike-in DNA samples and the data presented here provide a stable benchmark against which future ChIP platforms, protocol improvements, and analysis methods can be evaluated.

Genome-wide mapping and analysis of active promoters in mouse embryonic stem cells and adult organs.

By integrating genome-wide maps of RNA polymerase II (Polr2a) binding with gene expression data and H3ac and H3K4me3 profiles, we characterized promoters with enriched activity in mouse embryonic stem cells (mES) as well as adult brain, heart, kidney, and liver. We identified approximately 24,000 promoters across these samples, including 16,976 annotated mRNA 5' ends and 5153 additional sites validating cap-analysis of gene expression (CAGE) 5' end data. We showed that promoters with CpG islands are typically non-tissue specific, with the majority associated with Polr2a and the active chromatin modifications in nearly all the tissues examined. By contrast, the promoters without CpG islands are generally associated with Polr2a and the active chromatin marks in a tissue-dependent way. We defined 4396 tissue-specific promoters by adapting a quantitative index of tissue-specificity based on Polr2a occupancy. While there is a general correspondence between Polr2a occupancy and active chromatin modifications at the tissue-specific promoters, a subset of them appear to be persistently marked by active chromatin modifications in the absence of detectable Polr2a binding, highlighting the complexity of the functional relationship between chromatin modification and gene expression. Our results provide a resource for exploring promoter Polr2a binding and epigenetic states across pluripotent and differentiated cell types in mammals.

A comprehensive ChIP-chip analysis of E2F1, E2F4, and E2F6 in normal and tumor cells reveals interchangeable roles of E2F family members.

Using ChIP-chip assays (employing ENCODE arrays and core promoter arrays), we examined the binding patterns of three members of the E2F family in five cell types. We determined that most E2F1, E2F4, and E2F6 binding sites are located within 2 kb of a transcription start site, in both normal and tumor cells. In fact, the majority of promoters that are active (as defined by TAF1 or POLR2A binding) in GM06990 B lymphocytes and Ntera2 carcinoma cells were also bound by an E2F. This very close relationship between E2F binding sites and binding sites for general transcription factors in both normal and tumor cells suggests that a chromatin-bound E2F may be a signpost for active transcription initiation complexes. In general, we found that several E2Fs bind to a given promoter and that there is only modest cell type specificity of the E2F family. Thus, it is difficult to assess the role of any particular E2F in transcriptional regulation, due to extreme redundancy of target promoters. However, Ntera2 carcinoma cells were exceptional in that a large set of promoters were bound by E2F6, but not by E2F1 or E2F4. It has been proposed that E2F6 contributes to gene silencing by recruiting enzymes involved in methylating histone H3. To test this hypothesis, we created Ntera2 cell lines harboring shRNAs to E2F6. We found that reduction of E2F6 only induced minimal alteration of the transcriptome of Ntera2 transcriptome. Our results support the concept of functional redundancy in the E2F family and suggest that E2F6 is not critical for histone methylation.

Arginine methylation at histone H3R2 controls deposition of H3K4 trimethylation.

Modifications on histones control important biological processes through their effects on chromatin structure. Methylation at lysine 4 on histone H3 (H3K4) is found at the 5' end of active genes and contributes to transcriptional activation by recruiting chromatin-remodelling enzymes. An adjacent arginine residue (H3R2) is also known to be asymmetrically dimethylated (H3R2me2a) in mammalian cells, but its location within genes and its function in transcription are unknown. Here we show that H3R2 is also methylated in budding yeast (Saccharomyces cerevisiae), and by using an antibody specific for H3R2me2a in a chromatin immunoprecipitation-on-chip analysis we determine the distribution of this modification on the entire yeast genome. We find that H3R2me2a is enriched throughout all heterochromatic loci and inactive euchromatic genes and is present at the 3' end of moderately transcribed genes. In all cases the pattern of H3R2 methylation is mutually exclusive with the trimethyl form of H3K4 (H3K4me3). We show that methylation at H3R2 abrogates the trimethylation of H3K4 by the Set1 methyltransferase. The specific effect on H3K4me3 results from the occlusion of Spp1, a Set1 methyltransferase subunit necessary for trimethylation. Thus, the inability of Spp1 to recognize H3 methylated at R2 prevents Set1 from trimethylating H3K4. These results provide the first mechanistic insight into the function of arginine methylation on chromatin.

Genome-wide analysis of KAP1 binding suggests autoregulation of KRAB-ZNFs.

We performed a genome-scale chromatin immunoprecipitation (ChIP)-chip comparison of two modifications (trimethylation of lysine 9 [H3me3K9] and trimethylation of lysine 27 [H3me3K27]) of histone H3 in Ntera2 testicular carcinoma cells and in three different anatomical sources of primary human fibroblasts. We found that in each of the cell types the two modifications were differentially enriched at the promoters of the two largest classes of transcription factors. Specifically, zinc finger (ZNF) genes were bound by H3me3K9 and homeobox genes were bound by H3me3K27. We have previously shown that the Polycomb repressive complex 2 is responsible for mediating trimethylation of lysine 27 of histone H3 in human cancer cells. In contrast, there is little overlap between H3me3K9 targets and components of the Polycomb repressive complex 2, suggesting that a different histone methyltransferase is responsible for the H3me3K9 modification. Previous studies have shown that SETDB1 can trimethylate H3 on lysine 9, using in vitro or artificial tethering assays. SETDB1 is thought to be recruited to chromatin by complexes containing the KAP1 corepressor. To determine if a KAP1-containing complex mediates trimethylation of the identified H3me3K9 targets, we performed ChIP-chip assays and identified KAP1 target genes using human 5-kb promoter arrays. We found that a large number of genes of ZNF transcription factors were bound by both KAP1 and H3me3K9 in normal and cancer cells. To expand our studies of KAP1, we next performed a complete genomic analysis of KAP1 binding using a 38-array tiling set, identifying ~7,000 KAP1 binding sites. The identified KAP1 targets were highly enriched for C2H2 ZNFs, especially those containing Krüppel-associated box (KRAB) domains. Interestingly, although most KAP1 binding sites were within core promoter regions, the binding sites near ZNF genes were greatly enriched within transcribed regions of the target genes. Because KAP1 is recruited to the DNA via interaction with KRAB-ZNF proteins, we suggest that expression of KRAB-ZNF genes may be controlled via an auto-regulatory mechanism involving KAP1.