Comparison of sample preparation methods for ChIP-chip assays

Nanoparticles and trained immunity: Glimpse into the future

Emerging evidences show that innate immune cells can display changes in their functional programs after infection or vaccination, which lead to immunomodulation (increased or reduced responsiveness) upon secondary activation to the same stimuli or even to a different one. Innate cells acquire features of immunological memory, nowadays using the new term of "trained immunity" or "innate immune memory", which is different from the specific memory immune response elicited by B and T lymphocytes. The review focused on the concept of trained immunity, mostly on myeloid cells. Special attention is dedicated to the pathogen recognition along the evolution (bacteria, plants, invertebrate and vertebrate animals), and to techniques used to study epigenetic reprogramming and metabolic rewiring. Nanomaterials can be recognized by immune cells offering a very promising way to learn about trained immunity. Nanomaterials could be modified in order to immunomodulate the responses ad hoc. Many therapeutic possibilities are opened, and they should be explored.

Methods and platforms for analysis of nucleic acids from single-cell based on microfluidics

Single-cell nucleic acid analysis aims at discovering the genetic differences between individual cells which is well known as the cellular heterogeneity. This technology facilitates cancer diagnosis, stem cell research, immune system analysis, and other life science applications. The conventional platforms for single-cell nucleic acid analysis more rely on manual operation or bulky devices. Recently, the emerging microfluidic technology has provided a perfect platform for single-cell nucleic acid analysis with the characteristic of accurate and automatic single-cell manipulation. In this review, we briefly summarized the procedure of single-cell nucleic acid analysis including single-cell isolation, single-cell lysis, nucleic acid amplification, and genetic analysis. And then, three representative microfluidic platforms for single-cell nucleic acid analysis are concluded as valve-, microwell-, and droplet-based platforms. Furthermore, we described the state-of-the-art integrated single-cell nucleic acid analysis systems based on the three platforms. Finally, the future development and challenges of microfluidics-based single-cell nucleic acid analysis are discussed as well.

E2F4 regulates transcriptional activation in mouse embryonic stem cells independently of the RB family

E2F transcription factors are central regulators of cell division and cell fate decisions. E2F4 often represents the predominant E2F activity in cells. E2F4 is a transcriptional repressor implicated in cell cycle arrest and whose repressive activity depends on its interaction with members of the RB family. Here we show that E2F4 is important for the proliferation and the survival of mouse embryonic stem cells. In these cells, E2F4 acts in part as a transcriptional activator that promotes the expression of cell cycle genes. This role for E2F4 is independent of the RB family. Furthermore, E2F4 functionally interacts with chromatin regulators associated with gene activation and we observed decreased histone acetylation at the promoters of cell cycle genes and E2F targets upon loss of E2F4 in RB family-mutant cells. Taken together, our findings uncover a non-canonical role for E2F4 that provide insights into the biology of rapidly dividing cells.

E2F transcription factors are regulators of cell division and cell fate decisions. Here the authors show that E2F4 is important for proliferation and survival of mouse ESCs, independent of the RB family, and that E2F4 interacts with chromatin regulators associated with gene activation.

CRISPR-mediated HDAC2 disruption identifies two distinct classes of target genes in human cells

The transcriptional functions of the class I histone deacetylases (HDACs) HDAC1 and HDAC2 are mainly viewed as both repressive and redundant based on murine knockout studies, but they may have additional independent roles and their physiological functions in human cells are not as clearly defined. To address the individual epigenomic functions of HDAC2, here we utilized CRISPR-Cas9 to disrupt HDAC2 in human cells. We find that while HDAC2 null cells exhibited signs of cross-regulation between HDAC1 and HDAC2, specific epigenomic phenotypes were still apparent using RNA-seq and ChIP assays. We identified specific targets of HDAC2 repression, and defined a novel class of genes that are actively expressed in a partially HDAC2-dependent manner. While HDAC2 was required for the recruitment of HDAC1 to repressed HDAC2-gene targets, HDAC2 was dispensable for HDAC1 binding to HDAC2-activated targets, supporting the notion of distinct classes of targets. Both active and repressed classes of gene targets demonstrated enhanced histone acetylation and methylation in HDAC2-null cells. Binding of the HDAC1/2-associated SIN3A corepressor was altered at most HDAC2-targets, but without a clear pattern. Overall, our study defines two classes of HDAC2 targets in human cells, with a dependence of HDAC1 on HDAC2 at one class of targets, and distinguishes unique functions for HDAC2.

LEC1 sequentially regulates the transcription of genes involved in diverse developmental processes during seed development

LEAFY COTYLEDON1 (LEC1), an atypical subunit of the nuclear transcription factor Y (NF-Y) CCAAT-binding transcription factor, is a central regulator that controls many aspects of seed development including the maturation phase during which seeds accumulate storage macromolecules and embryos acquire the ability to withstand desiccation. To define the gene networks and developmental processes controlled by LEC1, genes regulated directly by and downstream of LEC1 were identified. We compared the mRNA profiles of wild-type and lec1-null mutant seeds at several stages of development to define genes that are down-regulated or up-regulated by the lec1 mutation. We used ChIP and differential gene-expression analyses in Arabidopsis seedlings overexpressing LEC1 and in developing Arabidopsis and soybean seeds to identify globally the target genes that are transcriptionally regulated by LEC1 in planta Collectively, our results show that LEC1 controls distinct gene sets at different developmental stages, including those that mediate the temporal transition between photosynthesis and chloroplast biogenesis early in seed development and seed maturation late in development. Analyses of enriched DNA sequence motifs that may act as cis-regulatory elements in the promoters of LEC1 target genes suggest that LEC1 may interact with other transcription factors to regulate distinct gene sets at different stages of seed development. Moreover, our results demonstrate strong conservation in the developmental processes and gene networks regulated by LEC1 in two dicotyledonous plants that diverged ∼92 Mya.

Transcriptional regulatory logic of the diurnal cycle in the mouse liver

Many organisms exhibit temporal rhythms in gene expression that propel diurnal cycles in physiology. In the liver of mammals, these rhythms are controlled by transcription-translation feedback loops of the core circadian clock and by feeding-fasting cycles. To better understand the regulatory interplay between the circadian clock and feeding rhythms, we mapped DNase I hypersensitive sites (DHSs) in the mouse liver during a diurnal cycle. The intensity of DNase I cleavages cycled at a substantial fraction of all DHSs, suggesting that DHSs harbor regulatory elements that control rhythmic transcription. Using chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq), we found that hypersensitivity cycled in phase with RNA polymerase II (Pol II) loading and H3K27ac histone marks. We then combined the DHSs with temporal Pol II profiles in wild-type (WT) and Bmal1-/- livers to computationally identify transcription factors through which the core clock and feeding-fasting cycles control diurnal rhythms in transcription. While a similar number of mRNAs accumulated rhythmically in Bmal1-/- compared to WT livers, the amplitudes in Bmal1-/- were generally lower. The residual rhythms in Bmal1-/- reflected transcriptional regulators mediating feeding-fasting responses as well as responses to rhythmic systemic signals. Finally, the analysis of DNase I cuts at nucleotide resolution showed dynamically changing footprints consistent with dynamic binding of CLOCK:BMAL1 complexes. Structural modeling suggested that these footprints are driven by a transient heterotetramer binding configuration at peak activity. Together, our temporal DNase I mappings allowed us to decipher the global regulation of diurnal transcription rhythms in the mouse liver.

Genome-wide DNA promoter methylation and transcriptome analysis in human adipose tissue unravels novel candidate genes for obesity

Objective/methods

DNA methylation plays an important role in obesity and related metabolic complications. We examined genome-wide DNA promoter methylation along with mRNA profiles in paired samples of human subcutaneous adipose tissue (SAT) and omental visceral adipose tissue (OVAT) from non-obese vs. obese individuals.

Results

We identified negatively correlated methylation and expression of several obesity-associated genes in our discovery dataset and in silico replicated ETV6 in two independent cohorts. Further, we identified six adipose tissue depot-specific genes (HAND2, HOXC6, PPARG, SORBS2, CD36, and CLDN1). The effects were further supported in additional independent cohorts. Our top hits might play a role in adipogenesis and differentiation, obesity, lipid metabolism, and adipose tissue expandability. Finally, we show that in vitro methylation of SORBS2 directly represses gene expression.

Conclusions

Taken together, our data show distinct tissue specific epigenetic alterations which associate with obesity.

•

Obesity-associated differences in DNA promoter methylation and transcriptome in human adipose tissue (ETV6).

•

Depot-specific analyses revealed novel/known genes (HAND2, HOXC6, PPARG, SORBS2, CD36, CLDN1).

•

EWAS revealed SSPN and CCDC125 associated to BMI in SAT or OVAT, respectively.

•

Differentially methylated genes overlap in part with GWAS hits for obesity and fat distribution.

Down-regulatory mechanism of mammea E/BB from Mammea siamensis seed extract on Wilms' Tumor 1 expression in K562 cells

Background

Wilms’ tumor 1 (WT1) is a biological marker for predicting leukemia progression. In this study, mammea E/BB, an active compound from Saraphi (Mammea siamensis) seed extract was examined for its effect on down-regulatory mechanism of WT1 gene expression, WT1 protein and mRNA stability, and cell proliferation in K562 cell line.

Methods

M. siamensis seeds were obtained from the region of Chiang Mai (North of Thailand). Mammea E/BB was extracted from seeds of M. siamensis. WT1 protein expression and stability were evaluated by Western blot analysis. WT1 mRNA stability was assessed by qRT-PCR. WT1-DNA binding and WT1 promoter activity were assayed by ChIP assay and luciferase-reporter assay, respectively. Cell cycle arrest was studied by flow cytometry.

Results

Treatment with mammea E/BB led to down-regulation of WT1 expression. The suppression of WT1 expression did not involve protein and mRNA degradation. Rather, WT1 protein was down-regulated through disruption of transcriptional auto-regulation of the WT1 gene. Mammea E/BB inhibited WT1-DNA binding at the WT1 promoter and decreased luciferase activity. It also disrupted c-Fos/AP-1 binding to the WT1 promoter via ERK1/2 signaling pathway and induced S phase cell cycle arrest in K562 cells.

Conclusion

Mammea E/BB had pleotropic effects on kinase signaling pathways, resulting in inhibition of leukemia cell proliferation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12906-016-1107-z) contains supplementary material, which is available to authorized users.

Reciprocal changes in DNA methylation and hydroxymethylation and a broad repressive epigenetic switch characterize FMR1 transcriptional silencing in fragile X syndrome

Background

Fragile X syndrome (FXS) is the most common form of inherited intellectual disability, resulting from the loss of function of the fragile X mental retardation 1 (FMR1) gene. The molecular pathways associated with FMR1 epigenetic silencing are still elusive, and their characterization may enhance the discovery of novel therapeutic targets as well as the development of novel clinical biomarkers for disease status.

Results

We have deployed customized epigenomic profiling assays to comprehensively map the FMR1 locus chromatin landscape in peripheral mononuclear blood cells (PBMCs) from eight FXS patients and in fibroblast cell lines derived from three FXS patient. Deoxyribonucleic acid (DNA) methylation (5-methylcytosine (5mC)) and hydroxymethylation (5-hydroxymethylcytosine (5hmC)) profiling using methylated DNA immunoprecipitation (MeDIP) combined with a custom FMR1 microarray identifies novel regions of DNA (hydroxy)methylation changes within the FMR1 gene body as well as in proximal flanking regions. At the region surrounding the FMR1 transcriptional start sites, increased levels of 5mC were associated to reciprocal changes in 5hmC, representing a novel molecular feature of FXS disease. Locus-specific validation of FMR1 5mC and 5hmC changes highlighted inter-individual differences that may account for the expected DNA methylation mosaicism observed at the FMR1 locus in FXS patients. Chromatin immunoprecipitation (ChIP) profiling of FMR1 histone modifications, together with 5mC/5hmC and gene expression analyses, support a functional relationship between 5hmC levels and FMR1 transcriptional activation and reveal cell-type specific differences in FMR1 epigenetic regulation. Furthermore, whilst 5mC FMR1 levels positively correlated with FXS disease severity (clinical scores of aberrant behavior), our data reveal for the first time an inverse correlation between 5hmC FMR1 levels and FXS disease severity.

Conclusions

We identify novel, cell-type specific, regions of FMR1 epigenetic changes in FXS patient cells, providing new insights into the molecular mechanisms of FXS. We propose that the combined measurement of 5mC and 5hmC at selected regions of the FMR1 locus may significantly enhance FXS clinical diagnostics and patient stratification.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0181-x) contains supplementary material, which is available to authorized users.

CTCF regulates NELF, DSIF and P-TEFb recruitment during transcription

CTCF is a versatile transcription factor with well-established roles in chromatin organization and insulator function. Recent findings also implicate CTCF in the control of elongation by RNA polymerase (RNAP) II. Here we show that CTCF knockdown abrogates RNAP II pausing at the early elongation checkpoint of c-myc by affecting recruitment of DRB-sensitivity-inducing factor (DSIF). CTCF knockdown also causes a termination defect on the U2 snRNA genes (U2), by affecting recruitment of negative elongation factor (NELF). In addition, CTCF is required for recruitment of positive elongation factor b (P-TEFb), which phosphorylates NELF, DSIF, and Ser2 of the RNAP II CTD to activate elongation of transcription of c-myc and recognition of the snRNA gene-specific 3' box RNA processing signal. These findings implicate CTCF in a complex network of protein:protein/protein:DNA interactions and assign a key role to CTCF in controlling RNAP II transcription through the elongation checkpoint of the protein-coding c-myc and the termination site of the non-coding U2, by regulating the recruitment and/or activity of key players in these processes.

Age-Related Increase of EED Expression in Early Hematopoietic Progenitor Cells is Associated with Global Increase of the Histone Modification H3K27me3

Human hematopoietic stem and progenitor cells (HSPCs) from umbilical cord blood exhibit higher differentiation potential and repopulation capacity compared to adult HSPCs. The molecular basis for these functional differences is currently unknown. Upon screening for epigenetic effector genes being differentially expressed in neonatal and adult HSPC subpopulations, the Polycomb Repressive Complex 2 (PRC2) member EED was identified. Even though EED is expressed at comparable amounts in neonatal and adult multipotent HSPCs, early adult lineage committed progenitors of the lymphomyeloid (LM) and erythromyeloid lineages expressed higher EED amounts than neonatal HPCs. We demonstrate that EED overexpression directly leads to higher H3K27me3 levels, a repressive histone modification that is mediated by the PRC2 complex. Quantitative analysis of H3K27me3 levels by FPLC-based ELISA revealed elevated levels in primary blood cells from adults. Besides quantitative changes, gene ontology analysis of the genome-wide H3K27me3 distribution revealed qualitative changes in adult HSPCs with elevated levels in genes associated with nonhematopoietic development pathways. In contrast, H3K4me3 which labels active chromatin was enriched on hematopoietic genes. In vitro differentiation of EED-transfected neonatal HSPCs revealed aberrant expression of the myelopoietic marker CD14, suggesting that EED affects the lymphoid versus myeloid decision processes within the lymphomyeloid lineage. This is in line with LM progenitors having the most pronounced differences in EED expression. Highlighting the dynamic roles of epigenetic modifications in human hematopoiesis, the present data demonstrate shifts in the PRC2-associated histone modification H3K27me3 from birth to adulthood.

ChIP on chip and ChIP-Seq assays: genome-wide analysis of transcription factor binding and histone modifications

Deregulation of transcriptional activity of many genes has been causatively linked to human diseases including cancer. Altered patterns of gene expression in normal and cancer cells are the result of inappropriate expression of transcription factors and chromatin modifying proteins. Chromatin immunoprecipitation assay is a well-established tool for investigating the interactions between regulatory proteins and DNA at distinct stages of gene activation. ChIP coupled with DNA microarrays, known as ChIP on chip, or sequencing of DNA associated with the factors (ChIP-Seq) allow us to determine the entire spectrum of in vivo DNA binding sites for a given protein. This has been of immense value because ChIP on chip assays and ChIP-Seq experiments can provide a snapshot of the transcriptional regulatory mechanisms on a genome-wide scale. This chapter outlines the general strategies used to carry out ChIP-chip assays to study the differential recruitment of regulatory molecules based on the studies conducted in our lab as well as other published protocols; these can be easily modified to a ChIP-Seq analysis.

Profiling gene promoter occupancy of Sox2 in two phenotypically distinct breast cancer cell subsets using chromatin immunoprecipitation and genome-wide promoter microarrays

Introduction

Aberrant expression of the embryonic stem cell marker Sox2 has been reported in breast cancer (BC). We previously identified two phenotypically distinct BC cell subsets separated based on their differential response to a Sox2 transcription activity reporter, namely the reporter-unresponsive (RU) and the more tumorigenic reporter-responsive (RR) cells. We hypothesized that Sox2, as a transcription factor, contributes to their phenotypic differences by mediating differential gene expression in these two cell subsets.

Methods

We used chromatin immunoprecipitation and a human genome-wide promoter microarray (ChIP-chip) to determine the promoter occupancies of Sox2 in the MCF7 RU and RR breast cancer cell populations. We validated our findings with conventional chromatin immunoprecipitation, quantitative reverse transcription polymerase chain reaction (qPCR), and western blotting using cell lines, and also performed qPCR using patient RU and RR samples.

Results

We found a largely mutually exclusive profile of gene promoters bound by Sox2 between RU and RR cells derived from MCF7 (1830 and 456 genes, respectively, with only 62 overlapping genes). Sox2 was bound to stem cell- and cancer-associated genes in RR cells. Using quantitative RT-PCR, we confirmed that 15 such genes, including PROM1 (CD133), BMI1, GPR49 (LGR5), and MUC15, were expressed significantly higher in RR cells. Using siRNA knockdown or enforced expression of Sox2, we found that Sox2 directly contributes to the higher expression of these genes in RR cells. Mucin-15, a novel Sox2 downstream target in BC, contributes to the mammosphere formation of BC cells. Parallel findings were observed in the RU and RR cells derived from patient samples.

Conclusions

In conclusion, our data supports the model that the Sox2 induces differential gene expression in the two distinct cell subsets in BC, and contributes to their phenotypic differences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-014-0470-2) contains supplementary material, which is available to authorized users.

Global loss of DNA methylation uncovers intronic enhancers in genes showing expression changes

BACKGROUND

Gene expression is epigenetically regulated by a combination of histone modifications and methylation of CpG dinucleotides in promoters. In normal cells, CpG-rich promoters are typically unmethylated, marked with histone modifications such as H3K4me3, and are highly active. During neoplastic transformation, CpG dinucleotides of CG-rich promoters become aberrantly methylated, corresponding with the removal of active histone modifications and transcriptional silencing. Outside of promoter regions, distal enhancers play a major role in the cell type-specific regulation of gene expression. Enhancers, which function by bringing activating complexes to promoters through chromosomal looping, are also modulated by a combination of DNA methylation and histone modifications.

RESULTS

Here we use HCT116 colorectal cancer cells with and without mutations in DNA methyltransferases, the latter of which results in a 95% reduction in global DNA methylation levels. These cells are used to study the relationship between DNA methylation, histone modifications, and gene expression. We find that the loss of DNA methylation is not sufficient to reactivate most of the silenced promoters. In contrast, the removal of DNA methylation results in the activation of a large number of enhancer regions as determined by the acquisition of active histone marks.

CONCLUSIONS

Although the transcriptome is largely unaffected by the loss of DNA methylation, we identify two distinct mechanisms resulting in the upregulation of distinct sets of genes. One is a direct result of DNA methylation loss at a set of promoter regions and the other is due to the presence of new intragenic enhancers.

Herpes Simplex Virus 1 (HSV-1) ICP22 protein directly interacts with cyclin-dependent kinase (CDK)9 to inhibit RNA polymerase II transcription elongation

The Herpes Simplex Virus 1 (HSV-1)-encoded ICP22 protein plays an important role in viral infection and affects expression of host cell genes. ICP22 is known to reduce the global level of serine (Ser)2 phosphorylation of the Tyr1Ser2Pro3Thr4Ser5Pro6Ser7 heptapeptide repeats comprising the carboxy-terminal domain (CTD) of the large subunit of RNA polymerase (pol) II. Accordingly, ICP22 is thought to associate with and inhibit the activity of the positive-transcription elongation factor b (P-TEFb) pol II CTD Ser2 kinase. We show here that ICP22 causes loss of CTD Ser2 phosphorylation from pol II engaged in transcription of protein-coding genes following ectopic expression in HeLa cells and that recombinant ICP22 interacts with the CDK9 subunit of recombinant P-TEFb. ICP22 also interacts with pol II in vitro. Residues 193 to 256 of ICP22 are sufficient for interaction with CDK9 and inhibition of pol II CTD Ser2 phosphorylation but do not interact with pol II. These results indicate that discrete regions of ICP22 interact with either CDK9 or pol II and that ICP22 interacts directly with CDK9 to inhibit expression of host cell genes.

Organ size control is dominant over Rb family inactivation to restrict proliferation in vivo

In mammals, a cell's decision to divide is thought to be under the control of the Rb/E2F pathway. We previously found that inactivation of the Rb family of cell cycle inhibitors (Rb, p107, and p130) in quiescent liver progenitors leads to uncontrolled division and cancer initiation. Here, we show that, in contrast, deletion of the entire Rb gene family in mature hepatocytes is not sufficient for their long-term proliferation. The cell cycle block in Rb family mutant hepatocytes is independent of the Arf/p53/p21 checkpoint but can be abrogated upon decreasing liver size. At the molecular level, we identify YAP, a transcriptional regulator involved in organ size control, as a factor required for the sustained expression of cell cycle genes in hepatocytes. These experiments identify a higher level of regulation of the cell cycle in vivo in which signals regulating organ size are dominant regulators of the core cell cycle machinery.

Global analysis of ZNF217 chromatin occupancy in the breast cancer cell genome reveals an association with ERalpha

BACKGROUND

The ZNF217 gene, encoding a C2H2 zinc finger protein, is located at 20q13 and found amplified and overexpressed in greater than 20% of breast tumors. Current studies indicate ZNF217 drives tumorigenesis, yet the regulatory mechanisms of ZNF217 are largely unknown. Because ZNF217 associates with chromatin modifying enzymes, we postulate that ZNF217 functions to regulate specific gene signaling networks. Here, we present a large-scale functional genomic analysis of ZNF217, which provides insights into the regulatory role of ZNF217 in MCF7 breast cancer cells.

RESULTS

ChIP-seq analysis reveals that the majority of ZNF217 binding sites are located at distal regulatory regions associated with the chromatin marks H3K27ac and H3K4me1. Analysis of ChIP-seq transcription factor binding sites shows clustering of ZNF217 with FOXA1, GATA3 and ERalpha binding sites, supported by the enrichment of corresponding motifs for the ERalpha-associated cis-regulatory sequences. ERalpha expression highly correlates with ZNF217 in lysates from breast tumors (n = 15), and ERalpha co-precipitates ZNF217 and its binding partner CtBP2 from nuclear extracts. Transcriptome profiling following ZNF217 depletion identifies differentially expressed genes co-bound by ZNF217 and ERalpha; gene ontology suggests a role for ZNF217-ERalpha in expression programs associated with ER+ breast cancer studies found in the Molecular Signature Database. Data-mining of expression data from breast cancer patients correlates ZNF217 with reduced overall survival.

CONCLUSIONS

Our genome-wide ZNF217 data suggests a functional role for ZNF217 at ERalpha target genes. Future studies will investigate whether ZNF217 expression contributes to aberrant ERalpha regulatory events in ER+ breast cancer and hormone resistance.

Convergent and divergent cellular responses by ErbB4 isoforms in mammary epithelial cells

Associations of ErbB4 (ERBB4/HER4), the fourth member of the EGFR family, with cancer are variable, possibly as a result of structural diversity of this receptor. There are multiple structural isoforms of ERBB4 arising by alternative mRNA splicing, and a subset undergo proteolysis that releases membrane-anchored and soluble isoforms that associate with transcription factors and coregulators to modulate transcription. To compare the differential and common signaling activities of full-length (FL) and soluble intracellular isoforms of ERBB4, four JM-a isoforms (FL and soluble intracellular domain (ICD) CYT-1 and CYT-2) were expressed in isogenic MCF10A cells and their biologic activities were analyzed. Both FL and ICD CYT-2 promoted cell proliferation and invasion, and CYT-1 suppressed cell growth. Transcriptional profiling revealed several new and underexplored ERBB4-regulated transcripts, including: proteases/protease inhibitors (MMP3 and SERPINE2), the YAP/Hippo pathway (CTGF, CYR61, and SPARC), the mevalonate/cholesterol pathway (HMGCR, HMGCS1, LDLR, and DHCR7), and cytokines (IL8, CCL20, and CXCL1). Many of these transcripts were subsequently validated in a luminal breast cancer cell line that normally expresses ERBB4. Furthermore, ChIP-seq experiments identified ADAP1, APOE, SPARC, STMN1, and MXD1 as novel molecular targets of ERBB4. These findings clarify the diverse biologic activities of ERBB4 isoforms, and reveal new and divergent functions.

IMPLICATIONS

ErbB4 as a regulator of Hippo and mevalonate pathways provides new insight into milk production and anabolic processes in normal mammary epithelia and cancer.

Promoter methylation profiles between human lung adenocarcinoma multidrug resistant A549/cisplatin (A549/DDP) cells and its progenitor A549 cells

Although aberrant DNA methylation has been implicated in the pathophysiology of lung cancer, the role of methylation in multidrug resistance (MDR) of lung cancer has remained unclear. To investigate whether certain distinct DNA methylation pattern is associated with acquired MDR of lung adenocarcinoma, methylated-DNA immunoprecipitation-chromatin immunoprecipitation (MeDIP-ChIP) was utilised to compare the genome-wide promoter methylation of the human lung adenocarcinoma MDR A549/cisplatin (A549/DDP) cells with its progenitor A549 cells. The comparison identified 3617 genes with differentially methylated promoter, of which 1581 were hypermethylated and 2036 were hypomethylated. Then, bisulphite sequencing polymerase chain reaction (PCR) (BSP) and quantitative reverse transcription (RT)-PCR (Q-PCR) were used to validate the promoter methylation of five candidate genes and to determine whether the expression of genes was associated with the promoter methylation. BSP confirmed that the promoter methylation incidence of the hypermethylated genes, G protein-coupled receptor 56 isoform 3 (GPR56), metallothionein 1G (MT1G), and RAS association domain family gene 1 (RASSF1), was significantly higher in A549/DDP cells compared with A549 cells (p<0.001, p=0.0099, and p=0.0165), whereas no significant difference was found in that of the other two genes, CCNL2 and BAD (p=0.0594 and p=0.5546). Additionally, Q-PCR showed that the mRNA expression of the three hypermethylated genes was significantly lower in A549/DDP cells compared with A549 cells (all p<0.001). In conclusion, this study reported for the first time that a distinct promoter methylation pattern is associated with MDR of lung adenocarcinoma A549/DDP cells and suggested that GPR56, MT1G, and RASSF1 might be the potential methylation markers associated with acquired MDR of lung adenocarcinoma.

Spatial separation of Xist RNA and polycomb proteins revealed by superresolution microscopy

In female mammals, one of the two X chromosomes is transcriptionally silenced to equalize X-linked gene dosage relative to XY males, a process termed X chromosome inactivation. Mechanistically, this is thought to occur via directed recruitment of chromatin modifying factors by the master regulator, X-inactive specific transcript (Xist) RNA, which localizes in cis along the entire length of the chromosome. A well-studied example is the recruitment of polycomb repressive complex 2 (PRC2), for which there is evidence of a direct interaction involving the PRC2 proteins Enhancer of zeste 2 (Ezh2) and Supressor of zeste 12 (Suz12) and the A-repeat region located at the 5' end of Xist RNA. In this study, we have analyzed Xist-mediated recruitment of PRC2 using two approaches, microarray-based epigenomic mapping and superresolution 3D structured illumination microscopy. Making use of an ES cell line carrying an inducible Xist transgene located on mouse chromosome 17, we show that 24 h after synchronous induction of Xist expression, acquired PRC2 binding sites map predominantly to gene-rich regions, notably within gene bodies. Paradoxically, these new sites of PRC2 deposition do not correlate with Xist-mediated gene silencing. The 3D structured illumination microscopy was performed to assess the relative localization of PRC2 proteins and Xist RNA. Unexpectedly, we observed significant spatial separation and absence of colocalization both in the inducible Xist transgene ES cell line and in normal XX somatic cells. Our observations argue against direct interaction between Xist RNA and PRC2 proteins and, as such, prompt a reappraisal of the mechanism for PRC2 recruitment in X chromosome inactivation.

In vivo mapping of arabidopsis scaffold/matrix attachment regions reveals link to nucleosome-disfavoring poly(dA:dT) tracts

Scaffold or matrix attachment regions (S/MARs) are found in all eukaryotes. The pattern of distribution and genomic context of S/MARs is thought to be important for processes such as chromatin organization and modulation of gene expression. Despite the importance of such processes, much is unknown about the large-scale distribution and sequence content of S/MARs in vivo. Here, we report the use of tiling microarrays to map 1358 S/MARs on Arabidopsis thaliana chromosome 4 (chr4). S/MARs occur throughout chr4, spaced much more closely than in the large plant and animal genomes that have been studied to date. Arabidopsis S/MARs can be divided into five clusters based on their association with other genomic features, suggesting a diversity of functions. While some Arabidopsis S/MARs may define structural domains, most occur near the transcription start sites of genes. Genes associated with these S/MARs have an increased probability of expression, which is particularly pronounced in the case of transcription factor genes. Analysis of sequence motifs and 6-mer enrichment patterns show that S/MARs are preferentially enriched in poly(dA:dT) tracts, sequences that resist nucleosome formation, and the majority of S/MARs contain at least one nucleosome-depleted region. This global view of S/MARs provides a framework to begin evaluating genome-scale models for S/MAR function.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

The RNA polymerase II C-terminal domain-interacting domain of yeast Nrd1 contributes to the choice of termination pathway and couples to RNA processing by the nuclear exosome

The RNA polymerase II (RNApII) C-terminal domain (CTD)-interacting domain (CID) proteins are involved in two distinct RNApII termination pathways and recognize different phosphorylated forms of CTD. To investigate the role of differential CTD-CID interactions in the choice of termination pathway, we altered the CTD-binding specificity of Nrd1 by domain swapping. Nrd1 with the CID from Rtt103 (Nrd1(CID(Rtt103))) causes read-through transcription at many genes, but can also trigger termination where multiple Nrd1/Nab3-binding sites and the Ser(P)-2 CTD co-exist. Therefore, CTD-CID interactions target specific termination complexes to help choose an RNApII termination pathway. Interactions of Nrd1 with both CTD and nascent transcripts contribute to efficient termination by the Nrd1 complex. Surprisingly, replacing the Nrd1 CID with that from Rtt103 reduces binding to Rrp6/Trf4, and RNA transcripts terminated by Nrd1(CID(Rtt103)) are predominantly processed by core exosome. Thus, the Nrd1 CID couples Ser(P)-5 CTD not only to termination, but also to RNA processing by the nuclear exosome.

Coordinated effects of sequence variation on DNA binding, chromatin structure, and transcription

DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.

Gene duplication and neofunctionalization: POLR3G and POLR3GL

RNA polymerase III (Pol III) occurs in two versions, one containing the POLR3G subunit and the other the closely related POLR3GL subunit. It is not clear whether these two Pol III forms have the same function, in particular whether they recognize the same target genes. We show that the POLR3G and POLR3GL genes arose from a DNA-based gene duplication, probably in a common ancestor of vertebrates. POLR3G- as well as POLR3GL-containing Pol III are present in cultured cell lines and in normal mouse liver, although the relative amounts of the two forms vary, with the POLR3G-containing Pol III relatively more abundant in dividing cells. Genome-wide chromatin immunoprecipitations followed by high-throughput sequencing (ChIP-seq) reveal that both forms of Pol III occupy the same target genes, in very constant proportions within one cell line, suggesting that the two forms of Pol III have a similar function with regard to specificity for target genes. In contrast, the POLR3G promoter—not the POLR3GL promoter—binds the transcription factor MYC, as do all other promoters of genes encoding Pol III subunits. Thus, the POLR3G/POLR3GL duplication did not lead to neo-functionalization of the gene product (at least with regard to target gene specificity) but rather to neo-functionalization of the transcription units, which acquired different mechanisms of regulation, thus likely affording greater regulation potential to the cell.

EED and KDM6B coordinate the first mammalian cell lineage commitment to ensure embryo implantation

The first mammalian cell lineage commitment is the formation of the trophectoderm (TE) and the inner cell mass (ICM) lineages during preimplantation development. Proper development of the TE and ICM lineages is dependent upon establishment of specific transcriptional programs. However, the epigenetic mechanisms that functionally contribute to establish TE- and ICM-specific transcriptional programs are poorly understood. Here, we show that proper development of the TE and ICM lineages is coordinated via combinatorial regulation of embryonic ectoderm development (EED) and lysine-specific demethylase 6B (KDM6B). During blastocyst formation, the relative levels of EED and KDM6B expression determine altered polycomb repressor 2 (PRC2) complex recruitment and incorporation of the repressive histone H3 lysine 27 trimethylation (H3K27Me3) mark at the chromatin domains of TE-specific master regulators CDX2 and GATA3, leading to their activation in the TE lineage and repression in the ICM lineage. Furthermore, ectopic gain of EED along with depletion of KDM6B in preimplantation mouse embryos abrogates CDX2 and GATA3 expression in the nascent TE lineage. The loss of CDX2 and GATA3 in the nascent TE lineage results in improper TE development, leading to failure in embryo implantation to the uterus. Our study delineates a novel epigenetic mechanism that orchestrates proper development of the first mammalian cell lineages.

Female bias in Rhox6 and 9 regulation by the histone demethylase KDM6A

The Rhox cluster on the mouse X chromosome contains reproduction-related homeobox genes expressed in a sexually dimorphic manner. We report that two members of the Rhox cluster, Rhox6 and 9, are regulated by de-methylation of histone H3 at lysine 27 by KDM6A, a histone demethylase with female-biased expression. Consistent with other homeobox genes, Rhox6 and 9 are in bivalent domains prior to embryonic stem cell differentiation and thus poised for activation. In female mouse ES cells, KDM6A is specifically recruited to Rhox6 and 9 for gene activation, a process inhibited by Kdm6a knockdown in a dose-dependent manner. In contrast, KDM6A occupancy at Rhox6 and 9 is low in male ES cells and knockdown has no effect on expression. In mouse ovary where Rhox6 and 9 remain highly expressed, KDM6A occupancy strongly correlates with expression. Our study implicates Kdm6a, a gene that escapes X inactivation, in the regulation of genes important in reproduction, suggesting that KDM6A may play a role in the etiology of developmental and reproduction-related effects of X chromosome anomalies.

Monitoring homology search during DNA double-strand break repair in vivo

Homologous recombination (HR) is crucial for genetic exchange and accurate repair of DNA double-strand breaks and is pivotal for genome integrity. HR uses homologous sequences for repair, but how homology search, the exploration of the genome for homologous DNA sequences, is conducted in the nucleus remains poorly understood. Here, we use time-resolved chromatin immunoprecipitations of repair proteins to monitor homology search in vivo. We found that homology search proceeds by a probing mechanism, which commences around the break and samples preferentially on the broken chromosome. However, elements thought to instruct chromosome loops mediate homology search shortcuts, and centromeres, which cluster within the nucleus, may facilitate homology search on other chromosomes. Our study thus reveals crucial parameters for homology search in vivo and emphasizes the importance of linear distance, chromosome architecture, and proximity for recombination efficiency.

Distinct signal transduction pathways downstream of the (P)RR revealed by microarray and ChIP-chip analyses

The (pro)renin receptor ((P)RR) signaling is involved in different pathophysiologies ranging from cardiorenal end-organ damage via diabetic retinopathy to tumorigenesis. We have previously shown that the transcription factor promyelocytic leukemia zinc finger (PLZF) is an adaptor protein of the (P)RR. Furthermore, recent publications suggest that major functions of the (P)RR are mediated ligand-independently by its transmembrane and intracellular part, which acts as an accessory protein of V-ATPases. The transcriptome and recruitmentome downstream of the V-ATPase function and PLZF in the context of the (P)RR are currently unknown. Therefore, we performed a set of microarray and chromatin-immunoprecipitation (ChIP)-chip experiments using siRNA against the (P)RR, stable overexpression of PLZF, the PLZF translocation inhibitor genistein and the specific V-ATPase inhibitor bafilomycin to dissect transcriptional pathways downstream of the (P)RR. We were able to identify distinct and overlapping genetic signatures as well as novel real-time PCR-validated target genes of the different molecular functions of the (P)RR. Moreover, bioinformatic analyses of our data confirm the role of (P)RŔs signal transduction pathways in cardiovascular disease and tumorigenesis.

Chromatin immunoprecipitation assays for Myc and N-Myc

Myc and N-Myc have widespread impacts on the chromatin state within cells, both in a gene-specific and genome-wide manner. Our laboratory uses functional genomic methods including chromatin immunoprecipitation (ChIP), ChIP-chip, and, more recently, ChIP-seq to analyze the binding and genomic location of Myc. In this chapter, we describe an effective ChIP protocol using specific validated antibodies to Myc and N-Myc. We discuss the application of this protocol to several types of stem and cancer cells, with a focus on aspects of sample preparation prior to library preparation that are critical for successful Myc ChIP assays. Key variables are discussed and include the starting quantity of cells or tissue, lysis and sonication conditions, the quantity and quality of antibody used, and the identification of reliable target genes for ChIP validation.

Histone modification profiling in normal and transformed human embryonic stem cells using micro chromatin immunoprecipitation, scalable to genome-wide microarray analyses

Comparing normal human embryonic stem cells (hESCs) to those that have acquired cellular properties of neoplasm provides a unique opportunity to study the distinguishing molecular features of human cellular transformation. As global alterations in the epigenetic landscape are a common feature of cancer, we sought to investigate the loci-specific and global differences between normal and transformed hESCs using ChIP-PCR and ChIP-microarray (also known as ChIP-chip). Here, specific emphasis was placed on optimizing ChIP for low cell numbers (termed micro-ChIP; μChIP) towards applications where the target population is rare, such as the case for somatic human tumors containing a low frequency of cancer stem cell populations and for single-colony analysis of embryonic and induced pluripotent stem cells emerging from initial derivation. Using these methods, we suggest that μChIP-PCR and microarray analysis is thus a powerful technology for epigenetic profiling of cell populations relevant to developmental biology, cancer, and regenerative medicine where target populations regulating the biological process can only be isolated in small numbers.

Evaluating the localization and DNA binding complexity of histones in mature sperm

The paternal genome in many animal taxa is efficiently packaged into the sperm nucleus by switching from a histone (nucleosome)-based chromatin configuration to one using predominantly protamines. Nonetheless, various studies have shown that some nucleosomes, often containing modified histones are retained in mature sperm and bind DNA with distinct sequence compositions. Considering the significance of histone modifications in epigenetic phenomena and the fact that sperm histones and their bound DNA must be carried into the oocyte, this chapter describes methods aimed at examining and analysing the histone composition of sperm chromatin. The focus is on both microscopic visualisation and evaluation of sequence composition of histones and histone-bound DNA in human and mouse spermatozoa. However, similar methods may be applicable to the sperm of other mammalian and even non-mammalian classes.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles

Genome-wide rhythms in RNA polymerase II loading and dynamic chromatin remodeling underlie periodic gene expression during diurnal cycles in the mouse liver.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

In mammalian organs such as the liver, many metabolic and physiological processes occur preferentially at specific times during the 24-hour daily cycle. The timing of these rhythmic functions depends on a complex interplay between the endogenous circadian clock and environmental timing cues relayed through the master circadian clock in the suprachiasmatic nucleus, or via feeding rhythms. These rhythms can be implemented on several regulatory levels, and here we aimed at a better understanding of the transcriptional and epigenetic changes that regulate diurnal rhythms. We performed genome-wide analysis of the locations of RNA polymerase II (Pol II) and the epigenetic histone modifications H3K4me3 and H3K36me3 at specific times of day, relating these data to mRNA expression levels. Our analyses show that Pol II transcriptional rhythms are biphasic in mouse liver, having predominant peak activities in the morning and evening. Moreover, dynamic changes in histone marks lag transcription rhythms genome-wide, indicating that the epigenetic landscape can be remodeled during the 24-hour cycle. Finally, a quantitative analysis of temporal Pol II and mRNA accumulation profiles indicates that posttranscriptional regulation significantly contributes to the amplitude and phase of mRNA accumulation profiles. While many studies have analyzed how transcription and chromatin states are modified during irreversible cell differentiation processes, our work highlights how these states can evolve reversibly in a system exhibiting periodicity in time.

The endoglycosidase heparanase enters the nucleus of T lymphocytes and modulates H3 methylation at actively transcribed genes via the interplay with key chromatin modifying enzymes

The methylation of histones is a fundamental epigenetic process regulating gene expression programs in mammalian cells. Dysregulated patterns of histone methylation are directly implicated in malignant transformation. Here, we report the unexpected finding that the invasive extracellular matrix degrading endoglycosidase heparanase enters the nucleus of activated human T lymphocytes and regulates the transcription of a cohort of inducible immune response genes by controlling histone H3 methylation patterns. It was found that nuclear heparanase preferentially associates with euchromatin. Genome-wide ChIP-on-chip analyses showed that heparanase is recruited to both the promoter and transcribed regions of a distinct cohort of transcriptionally active genes. Knockdown and overexpression of the heparanase gene also showed that chromatin-bound heparanase is a prerequisite for the transcription of a subset of inducible immune response genes in activated T cells. Furthermore, the actions of heparanase seem to influence gene transcription by associating with the demethylase LSD1, preventing recruitment of the methylase MLL and thereby modifying histone H3 methylation patterns. These data indicate that heparanase belongs to an emerging class of proteins that play an important role in regulating transcription in addition to their well-recognized extra-nuclear functions.

Hepatitis B viral core protein disrupts human host gene expression by binding to promoter regions

BACKGROUND

The core protein (HBc) of hepatitis B virus (HBV) has been implicated in the malignant transformation of chronically-infected hepatocytes and displays pleiotropic functions, including RNA- and DNA-binding activities. However, the mechanism by which HBc interacts with the human genome to exert effects on hepatocyte function remains unknown. This study investigated the distribution of HBc binding to promoters in the human genome and evaluated its effects on the related genes' expression.

RESULTS

Whole-genome chromatin immunoprecipitation microarray (ChIP-on-chip) analysis was used to identify HBc-bound human gene promoters. Gene Ontology and pathway analyses were performed on related genes. The quantitative polymerase chain reaction assay was used to verify ChIP-on-chip results. Five novel genes were selected for luciferase reporter assay evaluation to assess the influence of HBc promoter binding. The HBc antibody immunoprecipitated approximately 3100 human gene promoters. Among these, 1993 are associated with known biological processes, and 2208 regulate genes with defined molecular functions. In total, 1286 of the related genes mediate primary metabolic processes, and 1398 encode proteins with binding activity. Sixty-four of the promoters regulate genes related to the mitogen-activated protein kinase (MAPK) pathways, and 41 regulate Wnt/beta-catenin pathway genes. The reporter gene assay indicated that HBc binding up-regulates proto-oncogene tyrosine-protein kinase (SRC), type 1 insulin-like growth factor receptor (IGF1R), and neurotrophic tyrosine kinase receptor 2 (NTRK2), and down-regulates v-Ha-ras Harvey rat sarcoma viral oncogene (HRAS).

CONCLUSION

HBc has the ability to bind a large number of human gene promoters, and can disrupt normal host gene expression. Manipulation of the transcriptional profile in HBV-infected hepatocytes may represent a key pathogenic mechanism of HBV infection.

Role of MXD3 in proliferation of DAOY human medulloblastoma cells

A subset of medulloblastomas, the most common brain tumor in children, is hypothesized to originate from granule neuron precursors (GNPs) in which the sonic hedgehog (SHH) pathway is over-activated. MXD3, a basic helix-look-helix zipper transcription factor of the MAD family, has been reported to be upregulated during postnatal cerebellar development and to promote GNP proliferation and MYCN expression. Mxd3 is upregulated in mouse models of medulloblastoma as well as in human medulloblastomas. Therefore, we hypothesize that MXD3 plays a role in the cellular events that lead to medulloblastoma biogenesis. In agreement with its proliferative role in GNPs, MXD3 knock-down in DAOY cells resulted in decreased proliferation. Sustained overexpression of MXD3 resulted in decreased cell numbers due to increased apoptosis and cell cycle arrest. Structure-function analysis revealed that the Sin3 interacting domain, the basic domain, and binding to E-boxes are essential for this activity. Microarray-based expression analysis indicated up-regulation of 84 genes and down-regulation of 47 genes. Potential direct MXD3 target genes were identified by ChIP-chip. Our results suggest that MXD3 is necessary for DAOY medulloblastoma cell proliferation. However, increased level and/or duration of MXD3 expression ultimately reduces cell numbers via increased cell death and cell cycle arrest.

Genomic targets of Brachyury (T) in differentiating mouse embryonic stem cells

Background

The T-box transcription factor Brachyury (T) is essential for formation of the posterior mesoderm and the notochord in vertebrate embryos. Work in the frog and the zebrafish has identified some direct genomic targets of Brachyury, but little is known about Brachyury targets in the mouse.

Methodology/Principal Findings

Here we use chromatin immunoprecipitation and mouse promoter microarrays to identify targets of Brachyury in embryoid bodies formed from differentiating mouse ES cells. The targets we identify are enriched for sequence-specific DNA binding proteins and include components of signal transduction pathways that direct cell fate in the primitive streak and tailbud of the early embryo. Expression of some of these targets, such as Axin2, Fgf8 and Wnt3a, is down regulated in Brachyury mutant embryos and we demonstrate that they are also Brachyury targets in the human. Surprisingly, we do not observe enrichment of the canonical T-domain DNA binding sequence 5′-TCACACCT-3′ in the vicinity of most Brachyury target genes. Rather, we have identified an (AC)_n repeat sequence, which is conserved in the rat but not in human, zebrafish or Xenopus. We do not understand the significance of this sequence, but speculate that it enhances transcription factor binding in the regulatory regions of Brachyury target genes in rodents.

Conclusions/Significance

Our work identifies the genomic targets of a key regulator of mesoderm formation in the early mouse embryo, thereby providing insights into the Brachyury-driven genetic regulatory network and allowing us to compare the function of Brachyury in different species.

CREB critically regulates action potential shape and duration in the adult mouse ventricle

The cAMP response element binding protein (CREB) belongs to the CREB/cAMP response element binding modulator/activating transcription factor 1 family of cAMP-dependent transcription factors mediating a regulation of gene transcription in response to cAMP. Chronic stimulation of β-adrenergic receptors and the cAMP-dependent signal transduction pathway by elevated plasma catecholamines play a central role in the pathogenesis of heart failure. Ion channel remodeling, particularly a decreased transient outward current (I(to)), and subsequent action potential (AP) prolongation are hallmarks of the failing heart. Here, we studied the role of CREB for ion channel regulation in mice with a cardiomyocyte-specific knockout of CREB (CREB KO). APs of CREB KO cardiomyocytes were prolonged with increased AP duration at 50 and 70% repolarization and accompanied by a by 51% reduction of I(to) peak amplitude as detected in voltage-clamp measurements. We observed a 29% reduction of Kcnd2/Kv4.2 mRNA in CREB KO cardiomyocytes mice while the other I(to)-related channel subunits Kv4.3 and KChIP2 were not different between groups. Accordingly, Kv4.2 protein was reduced by 37% in CREB KO. However, we were not able to detect a direct regulation of Kv4.2 by CREB. The I(to)-dependent AP prolongation went along with an increase of I(Na) and a decrease of I(Ca,L) associated with an upregulation of Scn8a/Nav1.6 and downregulation of Cacna1c/Cav1.2 mRNA in CREB KO cardiomyocytes. Our results from mice with cardiomyocyte-specific inactivation of CREB definitively indicate that CREB critically regulates the AP shape and duration in the mouse ventricle, which might have an impact on ion channel remodeling in situations of altered cAMP-dependent signaling like heart failure.

Genome-wide in vivo cross-linking of sequence-specific transcription factors

Immunoprecipitation of cross-linked chromatin in combination with microarrays (ChIP-chip) or ultra high-throughput sequencing (ChIP-seq) is widely used to map genome-wide in vivo transcription factor binding. Both methods employ initial steps of in vivo cross-linking, chromatin isolation, DNA fragmentation, and immunoprecipitation. For ChIP-chip, the immunoprecipitated DNA samples are then amplified, labeled, and hybridized to DNA microarrays. For ChIP-seq, the immunoprecipitated DNA is prepared for a sequencing library, and then the library DNA fragments are sequenced using ultra high-throughput sequencing platform. The protocols described here have been developed for ChIP-chip and ChIP-seq analysis of sequence-specific transcription factor binding in Drosophila embryos. A series of controls establish that these protocols have high sensitivity and reproducibility and provide a quantitative measure of relative transcription factor occupancy. The quantitative nature of the assay is important because regulatory transcription factors bind to highly overlapping sets of thousands of genomic regions and the unique regulatory specificity of each factor is determined by relative moderate differences in occupancy between factors at commonly bound regions.

Genomic promoter occupancy of runt-related transcription factor RUNX2 in Osteosarcoma cells identifies genes involved in cell adhesion and motility

Runt-related transcription factors (RUNX1, RUNX2, and RUNX3) are key lineage-specific regulators of progenitor cell growth and differentiation but also function pathologically as cancer genes that contribute to tumorigenesis. RUNX2 attenuates growth and stimulates maturation of osteoblasts during bone formation but is also robustly expressed in a subset of osteosarcomas, as well as in metastatic breast and prostate tumors. To assess the biological function of RUNX2 in osteosarcoma cells, we examined human genomic promoter interactions for RUNX2 using chromatin immunoprecipitation (ChIP)-microarray analysis in SAOS-2 cells. Promoter binding of both RUNX2 and RNA polymerase II was compared with gene expression profiles of cells in which RUNX2 was depleted by RNA interference. Many RUNX2-bound loci (1550 of 2339 total) exhibit promoter occupancy by RNA polymerase II and contain the RUNX consensus motif 5'-((T/A/C)G(T/A/C)GG(T/G). Gene ontology analysis indicates that RUNX2 controls components of multiple signaling pathways (e.g. WNT, TGFβ, TNFα, and interleukins), as well as genes linked to cell motility and adhesion (e.g. the focal adhesion-related genes FAK/PTK2 and TLN1). Our results reveal that siRNA depletion of RUNX2, PTK2, or TLN1 diminishes motility of U2OS osteosarcoma cells. Thus, RUNX2 binding to diverse gene loci may support the biological properties of osteosarcoma cells.

Notch signaling inhibits hepatocellular carcinoma following inactivation of the RB pathway

Mice lacking all three Rb genes in the liver develop tumors resembling specific subgroups of human hepatocellular carcinomas, and Notch activity appears to suppress the growth and progression of these tumors.

Hepatocellular carcinoma (HCC) is the third cancer killer worldwide with >600,000 deaths every year. Although the major risk factors are known, therapeutic options in patients remain limited in part because of our incomplete understanding of the cellular and molecular mechanisms influencing HCC development. Evidence indicates that the retinoblastoma (RB) pathway is functionally inactivated in most cases of HCC by genetic, epigenetic, and/or viral mechanisms. To investigate the functional relevance of this observation, we inactivated the RB pathway in the liver of adult mice by deleting the three members of the Rb (Rb1) gene family: Rb, p107, and p130. Rb family triple knockout mice develop liver tumors with histopathological features and gene expression profiles similar to human HCC. In this mouse model, cancer initiation is associated with the specific expansion of populations of liver stem/progenitor cells, indicating that the RB pathway may prevent HCC development by maintaining the quiescence of adult liver progenitor cells. In addition, we show that during tumor progression, activation of the Notch pathway via E2F transcription factors serves as a negative feedback mechanism to slow HCC growth. The level of Notch activity is also able to predict survival of HCC patients, suggesting novel means to diagnose and treat HCC.

Lrig1 is an estrogen-regulated growth suppressor and correlates with longer relapse-free survival in ERα-positive breast cancer

Lrig1 is the founding member of the Lrig family and has been implicated in the negative regulation of several oncogenic receptor tyrosine kinases including ErbB2. Lrig1 is expressed at low levels in several cancer types but is overexpressed in some prostate and colorectal tumors. Given this heterogeneity, whether Lrig1 functions to suppress or promote tumor growth remains a critical question. Previously, we found that Lrig1 was poorly expressed in ErbB2-positive breast cancer, suggesting that Lrig1 has a growth-inhibitory role in this tumor type. However, breast cancer is a complex disease, with ErbB2-positive tumors accounting for just 25% of all breast cancers. To gain a better understanding of the role of Lrig1 in breast cancer, we examined its expression in estrogen receptor α (ERα)-positive disease which accounts for the majority of breast cancers. We find that Lrig1 is expressed at significantly higher levels in ERα-positive disease than in ERα-negative disease. Our study provides a molecular rationale for Lrig1 enrichment in ERα-positive disease by showing that Lrig1 is a target of ERα. Estrogen stimulates Lrig1 accumulation and disruption of this induction enhances estrogen-dependent tumor cell growth, suggesting that Lrig1 functions as an estrogen-regulated growth suppressor. In addition, we find that Lrig1 expression correlates with prolonged relapse-free survival in ERα-positive breast cancer, identifying Lrig1 as a new prognostic marker in this setting. Finally, we show that ErbB2 activation antagonizes ERα-driven Lrig1 expression, providing a mechanistic explanation for Lrig1 loss in ErbB2-positive breast cancer. This work provides strong evidence for a growth-inhibitory role for Lrig1 in breast cancer.

Long-range epigenetic silencing associates with deregulation of Ikaros targets in colorectal cancer cells

Transcription factors are common targets of epigenetic inactivation in human cancer. Promoter hypermethylation and subsequent silencing of transcription factors can lead to further deregulation of their targets. In this study, we explored the potential epigenetic deregulation in cancer of Ikaros family genes, which code for essential transcription factors in cell differentiation and exhibit genetic defects in hematologic neoplasias. Unexpectedly, our analysis revealed that Ikaros undergoes very specific promoter hypermethylation in colorectal cancer, including in all the cell lines studied and around 64% of primary colorectal adenocarcinomas, with increasing proportions in advanced Duke's stages. Ikaros hypermethylation occurred in the context of a novel long-range epigenetic silencing (LRES) region. Reintroduction of Ikaros in colorectal cancer cells, ChIP-chip analysis, and validation in primary samples led us to identify a number of direct targets that are possibly related with colorectal cancer progression. Our results not only provide the first evidence that LRES can have functional specific effects in cancer but also identify several deregulated Ikaros targets that may contribute to progression in colorectal adenocarcinoma.

Chromatin dynamics of gene activation and repression in response to interferon alpha (IFN(alpha)) reveal new roles for phosphorylated and unphosphorylated forms of the transcription factor STAT2

Signal transducer and activator of transcription 2 (STAT2), the critical component of type I interferons signaling, is a prototype latent cytoplasmic signal-dependent transcription factor. Activated tyrosine-phosphorylated STAT2 associates with STAT1 and IRF9 to bind the ISRE elements in the promoters of a subset of IFN-inducible genes (ISGs). In addition to activate hundreds of ISGs, IFNα also represses numerous target genes but the mechanistic basis for this dual effect and transcriptional repression is largely unknown. We investigated by ChIP-chip the binding dynamics of STAT2 and "active" phospho(P)-STAT2 on 113 putative IFNα direct target promoters before and after IFNα induction in Huh7 cells and primary human hepatocytes (PHH). STAT2 is already bound to 62% of our target promoters, including most "classical" ISGs, before IFNα treatment. 31% of STAT2 basally bound promoters also show P-STAT2 positivity. By correlating in vivo promoter occupancy with gene expression and changes in histone methylation marks we found that: 1) STAT2 plays a role in regulating ISGs expression, independently from its phosphorylation; 2) P-STAT2 is involved in ISGs repression; 3) "activated" ISGs are marked by H3K4me1 and H3K4me3 before IFNα; 4) "repressed" genes are marked by H3K27me3 and histone methylation plays a dominant role in driving IFNα-mediated ISGs repression.

Inhibitory mechanism of pure curcumin on Wilms' tumor 1 (WT1) gene expression through the PKCα signaling pathway in leukemic K562 cells

The aim of this study was to investigate the inhibitory mechanism of pure curcumin on WT1 expression in leukemic K562 cells. Pure curcumin suppressed WT1 expression, independent of effects on protein degradation or WT1 mRNA stability. Chromatin immunoprecipitation and reporter gene assays indicate that pure curcumin treatment attenuates WT1 auto-regulation. Interestingly, PKCα inhibition mimicks the repressive effects of pure curcumin in K562 cells. Conversely, myristoylated PKCα over-expression increased WT1 expression and reversed the inhibitory effect of pure curcumin. Our study indicates that pure curcumin attenuates WT1 auto-regulatory function through inhibition of PKCα signaling in K562 cells.

Single-tube linear DNA amplification (LinDA) for robust ChIP-seq

Genome-wide profiling of transcription factors based on massive parallel sequencing of immunoprecipitated chromatin (ChIP-seq) requires nanogram amounts of DNA. Here we describe a high-fidelity, single-tube linear DNA amplification method (LinDA) for ChIP-seq and reChIP-seq with picogram DNA amounts obtained from a few thousand cells. This amplification technology will facilitate global analyses of transcription-factor binding and chromatin with very small cell populations, such as stem or cancer-initiating cells.