The related coactivators SRC-2 and SRC-3 interact with peroxisome proliferator activated receptor γ (PPARγ) to coordinate transcriptional circuits to promote adipogenesis. To identify potential coactivator redundancy during human adipogenesis at single cell resolution, we used high content analysis to quantify links between PPARγ, SRC-2, SRC-3, and lipogenesis. Because we detected robust increases and significant cell-cell heterogeneity in PPARγ and lipogenesis, without changes in SRC-2 or SRC-3, we hypothesized that permissive coregulator levels comprise a necessary adipogenic equilibrium. We probed this equilibrium by down-regulating SRC-2 and SRC-3 while simultaneously quantifying PPARγ. Individual or joint knockdown equally inhibits lipid accumulation by preventing lipogenic gene engagement, without affecting PPARγ protein levels. Supporting dominant, pro-adipogenic roles for SRC-2 and SRC-3, SRC-1 knockdown does not affect adipogenesis. SRC-2 and SRC-3 knockdown increases the proportion of cells in a PPARγ(hi)/lipid(lo) state while increasing phospho-PPARγ-S114, an inhibitor of PPARγ transcriptional activity and adipogenesis. Together, we demonstrate that SRC-2 and SRC-3 concomitantly promote human adipocyte differentiation by attenuating phospho-PPARγ-S114 and modulating PPARγ cellular heterogeneity.

Estrogen receptor-alpha (ER) transcription function is regulated in a ligand-dependent (e.g., estradiol, E2) or ligand-independent (e.g., growth factors) manner. Our laboratory seeks to understand these two modes of action. Using a cell line that contains a visible prolactin enhancer/promoter array (PRL-HeLa) regulated by ER, we analyzed ER response to E2 and EGF by quantifying image-based results. Data show differential recruitment of GFP-ER to the array, with the AF1 domain playing a vital role in EGF-mediated responsiveness. Temporal analyses of large-scale chromatin dynamics, and accumulation of array-localized reporter mRNA over 24 hours showed that the EGF response consists of a single pulse of reporter mRNA accumulation concomitant with transient increase in array decondensation. Estradiol induced a novel cyclical pattern of mRNA accumulation with a sustained increase in array decondensation. Collectively, our work shows that there is a stimuli-specific pattern of large-scale chromatin modification and transcript levels by ER.

Annexin A1 (ANXA1) is an anti-inflammatory protein reported to play a role in cell proliferation and apoptosis, and to be deregulated in breast cancer. The exact role of annexin A1 in the biology of breast cancer remains unclear. We hypothesized that the annexin A1 plays an oncogenic role in basal subtype of breast cancer by modulating key growth pathway(s).

The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage.

TP53 is commonly altered in human cancer, and Tp53 reactivation suppresses tumours in vivo in mice (TP53 and Tp53 are also known as p53). This strategy has proven difficult to implement therapeutically, and here we examine an alternative strategy by manipulating the p53 family members, Tp63 and Tp73 (also known as p63 and p73, respectively). The acidic transactivation-domain-bearing (TA) isoforms of p63 and p73 structurally and functionally resemble p53, whereas the ΔN isoforms (lacking the acidic transactivation domain) of p63 and p73 are frequently overexpressed in cancer and act primarily in a dominant-negative fashion against p53, TAp63 and TAp73 to inhibit their tumour-suppressive functions. The p53 family interacts extensively in cellular processes that promote tumour suppression, such as apoptosis and autophagy, thus a clear understanding of this interplay in cancer is needed to treat tumours with alterations in the p53 pathway. Here we show that deletion of the ΔN isoforms of p63 or p73 leads to metabolic reprogramming and regression of p53-deficient tumours through upregulation of IAPP, the gene that encodes amylin, a 37-amino-acid peptide co-secreted with insulin by the β cells of the pancreas. We found that IAPP is causally involved in this tumour regression and that amylin functions through the calcitonin receptor (CalcR) and receptor activity modifying protein 3 (RAMP3) to inhibit glycolysis and induce reactive oxygen species and apoptosis. Pramlintide, a synthetic analogue of amylin that is currently used to treat type 1 and type 2 diabetes, caused rapid tumour regression in p53-deficient thymic lymphomas, representing a novel strategy to target p53-deficient cancers.

MicroRNAs modulate tumorigenesis through suppression of specific genes. As many tumour types rely on overlapping oncogenic pathways, a core set of microRNAs may exist, which consistently drives or suppresses tumorigenesis in many cancer types. Here we integrate The Cancer Genome Atlas (TCGA) pan-cancer data set with a microRNA target atlas composed of publicly available Argonaute Crosslinking Immunoprecipitation (AGO-CLIP) data to identify pan-tumour microRNA drivers of cancer. Through this analysis, we show a pan-cancer, coregulated oncogenic microRNA 'superfamily' consisting of the miR-17, miR-19, miR-130, miR-93, miR-18, miR-455 and miR-210 seed families, which cotargets critical tumour suppressors via a central GUGC core motif. We subsequently define mutations in microRNA target sites using the AGO-CLIP microRNA target atlas and TCGA exome-sequencing data. These combined analyses identify pan-cancer oncogenic cotargeting of the phosphoinositide 3-kinase, TGFβ and p53 pathways by the miR-17-19-130 superfamily members.

Prostate cancer remains a deadly disease especially when patients become resistant to drugs that target the Androgen Receptor (AR) ligand binding domain. At this stage, patients develop recurring castrate-resistant prostate cancers (CRPCs). Interestingly, CRPC tumors maintain dependency on AR for growth; moreover, in CRPCs, constitutively active AR splice variants (e.g., AR-V7) begin to be expressed at higher levels. These splice variants lack the ligand binding domain and are rendered insensitive to current endocrine therapies. Thus, it is of paramount importance to understand what regulates the expression of AR and its splice variants to identify new therapeutic strategies in CRPCs. Here, we used high throughput microscopy and quantitative image analysis to evaluate effects of selected endocrine disruptors on AR levels in multiple breast and prostate cancer cell lines. Bisphenol AP (BPAP), which is used in chemical and medical industries, was identified as a down-regulator of both full length AR and the AR-V7 splice variant. We validated its activity by performing time-course, dose-response, Western blot and qPCR analyses. BPAP also reduced the percent of cells in S phase, which was accompanied by a ~60% loss in cell numbers and colony formation in anchorage-independent growth assays. Moreover, it affected mitochondria size and cell metabolism. In conclusion, our high content analysis-based screening platform was used to classify the effect of compounds on endogenous ARs, and identified BPAP as being capable of causing AR (both full-length and variants) down-regulation, cell cycle arrest and metabolic alterations in CRPC cell lines.

Loss of Dicer, an enzyme critical for microRNA biogenesis, results in lethality due to a block in mouse embryonic stem cell (mES) differentiation. Using ChIP-Seq we found increased H3K9me2 at over 900 CpG islands in the Dicer(-/-)ES epigenome. Gene ontology analysis revealed that promoters of chromatin regulators to be among the most impacted by increased CpG island H3K9me2 in ES (Dicer(-/-)). We therefore, extended the study to include H3K4me3 and H3K27me3 marks for selected genes. We found that the ES (Dicer(-/-)) mutant epigenome was characterized by a shift in the overall balance between transcriptionally favorable (H3K4me3) and unfavorable (H3K27me3) marks at key genes regulating ES cell differentiation. Pluripotency genes Oct4, Sox2 and Nanog were not impacted in relation to patterns of H3K27me3 and H3K4me3 and showed no changes in the rates of transcript down-regulation in response to RA. The most striking changes were observed in regards to genes regulating differentiation and the transition from self-renewal to differentiation. An increase in H3K4me3 at the promoter of Lin28b was associated with the down-regulation of this gene at a lower rate in Dicer(-/-)ES as compared to wild type ES. An increase in H3K27me3 in the promoters of differentiation genes Hoxa1 and Cdx2 in Dicer(-/-)ES cells was coincident with an inability to up-regulate these genes at the same rate as ES upon retinoic acid (RA)-induced differentiation. We found that siRNAs Ezh2 and post-transcriptional silencing of Ezh2 by let-7 g rescued this effect suggesting that Ezh2 up-regulation is in part responsible for increased H3K27me3 and decreased rates of up-regulation of differentiation genes in Dicer(-/-)ES.

Despite advances in early detection and adjuvant targeted therapies, breast cancer is still the second most common cause of cancer mortality among women. Tumor recurrence is one of the major contributors to breast cancer mortality. However, the mechanisms underlying this process are not completely understood. In this study, we investigated the mechanisms of tumor dormancy and recurrence in a preclinical mouse model of breast cancer.

We used targeted mass spectrometry to study the metabolic fingerprint of urothelial cancer and determine whether the biochemical pathway analysis gene signature would have a predictive value in independent cohorts of patients with bladder cancer.

Metastatic prostate cancer (PCa) is primarily an androgen-dependent disease, which is treated with androgen deprivation therapy (ADT). Tumors usually develop resistance (castration-resistant PCa [CRPC]), but remain androgen receptor (AR) dependent. Numerous mechanisms for AR-dependent resistance have been identified including expression of constitutively active AR splice variants lacking the hormone-binding domain. Recent clinical studies show that expression of the best-characterized AR variant, AR-V7, correlates with resistance to ADT and poor outcome. Whether AR-V7 is simply a constitutively active substitute for AR or has novel gene targets that cause unique downstream changes is unresolved. Several studies have shown that AR activation alters cell metabolism. Using LNCaP cells with inducible expression of AR-V7 as a model system, we found that AR-V7 stimulated growth, migration, and glycolysis measured by ECAR (extracellular acidification rate) similar to AR. However, further analyses using metabolomics and metabolic flux assays revealed several differences. Whereas AR increased citrate levels, AR-V7 reduced citrate mirroring metabolic shifts observed in CRPC patients. Flux analyses indicate that the low citrate is a result of enhanced utilization rather than a failure to synthesize citrate. Moreover, flux assays suggested that compared to AR, AR-V7 exhibits increased dependence on glutaminolysis and reductive carboxylation to produce some of the TCA (tricarboxylic acid cycle) metabolites. These findings suggest that these unique actions represent potential therapeutic targets.

Osteosarcoma (OS) is a highly aggressive pediatric bone cancer in which most tumor cells remain immature and fail to differentiate into bone-forming osteoblasts. However, OS cells readily respond to adipogenic stimuli suggesting they retain mesenchymal stem cell-like properties. Here we demonstrate that nuclear receptor PPARγ agonists such as the anti-diabetic, thiazolidinedione (TZD) drugs induce growth arrest and cause adipogenic differentiation in human, mouse and canine OS cells as well as in tumors in mice. Gene expression analysis reveals that TZDs induce lipid metabolism pathways while suppressing targets of the Hippo-YAP pathway, Wnt signaling and cancer-related proliferation pathways. Significantly, TZD action appears to be restricted to the high Sox2 expressing cancer stem cell population and is dependent on PPARγ expression. TZDs also affect growth and cell fate by causing the cytoplasmic sequestration of the transcription factors SOX2 and YAP that are required for tumorigenicity. Finally, we identify a TZD-regulated gene signature based on Wnt/Hippo target genes and PPARγ that predicts patient outcomes. Together, this work highlights a novel connection between PPARγ agonist in inducing adipogenesis and mimicking the tumor suppressive hippo pathway. It also illustrates the potential of drug repurposing for TZD-based differentiation therapy for osteosarcoma.

Cutaneous squamous cell carcinoma (cuSCC) comprises 15-20% of all skin cancers, accounting for over 700,000 cases in USA annually. Most cuSCC arise in association with a distinct precancerous lesion, the actinic keratosis (AK). To identify potential targets for molecularly targeted chemoprevention, here we perform integrated cross-species genomic analysis of cuSCC development through the preneoplastic AK stage using matched human samples and a solar ultraviolet radiation-driven Hairless mouse model. We identify the major transcriptional drivers of this progression sequence, showing that the key genomic changes in cuSCC development occur in the normal skin to AK transition. Our data validate the use of this ultraviolet radiation-driven mouse cuSCC model for cross-species analysis and demonstrate that cuSCC bears deep molecular similarities to multiple carcinogen-driven SCCs from diverse sites, suggesting that cuSCC may serve as an effective, accessible model for multiple SCC types and that common treatment and prevention strategies may be feasible.

The hotspots of structural polymorphisms and structural mutability in the human genome remain to be explained mechanistically. We examine associations of structural mutability with germline DNA methylation and with non-allelic homologous recombination (NAHR) mediated by low-copy repeats (LCRs). Combined evidence from four human sperm methylome maps, human genome evolution, structural polymorphisms in the human population, and previous genomic and disease studies consistently points to a strong association of germline hypomethylation and genomic instability. Specifically, methylation deserts, the ~1% fraction of the human genome with the lowest methylation in the germline, show a tenfold enrichment for structural rearrangements that occurred in the human genome since the branching of chimpanzee and are highly enriched for fast-evolving loci that regulate tissue-specific gene expression. Analysis of copy number variants (CNVs) from 400 human samples identified using a custom-designed array comparative genomic hybridization (aCGH) chip, combined with publicly available structural variation data, indicates that association of structural mutability with germline hypomethylation is comparable in magnitude to the association of structural mutability with LCR-mediated NAHR. Moreover, rare CNVs occurring in the genomes of individuals diagnosed with schizophrenia, bipolar disorder, and developmental delay and de novo CNVs occurring in those diagnosed with autism are significantly more concentrated within hypomethylated regions. These findings suggest a new connection between the epigenome, selective mutability, evolution, and human disease.

Brown adipose tissue (BAT) dissipates chemical energy as heat via thermogenesis and protects against obesity by increasing energy expenditure. However, regulation of BAT by dietary factors remains largely unexplored at the mechanistic level. We investigated the effect of eicosapentaenoic acid (EPA) on BAT metabolism. Male C57BL/6J (B6) mice were fed either a high-fat diet (HF, 45% kcal fat) or HF diet supplemented with EPA (HF-EPA, 6.75% kcal EPA) for 11 weeks. RNA sequencing (RNA-Seq) and microRNA (miRNA) profiling were performed on RNA from BAT using Illumina HiSeq and Illumina Genome Analyzer NextSeq, respectively. We conducted pathway analyses using ingenuity pathway analysis software (IPA®) and validated some genes and miRNAs using qPCR. We identified 479 genes that were differentially expressed (2-fold change, n = 3, P ≤ 0.05) in BAT from HF compared to HF-EPA. Genes negatively correlated with thermogenesis such as hypoxia inducible factor 1 alpha subunit inhibitor (Hif1an), were downregulated by EPA. Pathways related to thermogenesis such as peroxisome proliferator-activated receptor (PPAR) were upregulated by EPA while pathways associated with obesity and inflammation such as nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) were downregulated by EPA. MiRNA profiling identified nine and six miRNAs that were upregulated and downregulated by EPA, respectively (log2 fold change > 1.25, n = 3, P ≤ 0.05). Key regulatory miRNAs which were involved in thermogenesis, such as miR-455-3p and miR-129-5p were validated using qPCR. In conclusion, the depth of transcriptomic and miRNA profiling revealed novel mRNA-miRNA interaction networks in BAT which are involved in thermogenesis, and regulated by EPA.

The architectural organization of chromatin can play an important role in genome regulation by affecting the mobility of molecules within its surroundings via binding interactions and molecular crowding. The diffusion of molecules at specific locations in the nucleus can be studied by fluorescence correlation spectroscopy (FCS), a well-established technique based on the analysis of fluorescence intensity fluctuations detected in a confocal observation volume. However, detecting subtle variations of mobility between different chromatin regions remains challenging with currently available FCS methods. Here, we introduce a method that samples multiple positions by slowly scanning the FCS observation volume across the nucleus. Analyzing the data in short time segments, we preserve the high temporal resolution of single-point FCS while probing different nuclear regions in the same cell. Using the intensity level of the probe (or a DNA marker) as a reference, we efficiently sort the FCS segments into different populations and obtain average correlation functions that are associated to different chromatin regions. This sorting and averaging strategy renders the method statistically robust while preserving the observation of intranuclear variations of mobility. Using this approach, we quantified diffusion of monomeric GFP in high versus low chromatin density regions. We found that GFP mobility was reduced in heterochromatin, especially within perinucleolar heterochromatin. Moreover, we found that modulation of chromatin compaction by ATP depletion, or treatment with solutions of different osmolarity, differentially affected the ratio of diffusion in both regions. Then, we used the approach to probe the mobility of estrogen receptor-α in the vicinity of an integrated multicopy prolactin gene array. Finally, we discussed the coupling of this method with stimulated emission depletion FCS for performing FCS at subdiffraction spatial scales.

Background: Eosinophilic esophagitis (EoE) is a disorder of the esophagus that has become increasingly recognized in children. Because these children undergo multiple endoscopies, discovering a non-invasive biomarker of disease activity is highly desirable. The aim of this study was to use targeted plasma metabolomics to identify potential biomarker candidates for EoE in a discovery phase. Methods: A prospective, single-center clinical trial was performed on 24 children ages 2-18 years with and without EoE undergoing upper endoscopy for any indication. Blood samples were collected for metabolomics profiling using the subclasses: amino acids, tricarboxylic acid cycle, acetylation, and methylation. Using mass spectrometry and systematic bioinformatics analysis, 48 metabolites were measured and compared between children with active EoE (+EoE) and controls (-EoE). To investigate the effect of proton pump inhibitor (PPI) use on metabolites, patients were also stratified based on PPI use (+PPI, -PPI). Results: Seven children had active EoE at the time of endoscopy. Eleven children were on PPI (4 with EoE). Of the 48 metabolites measured, 8 plasma metabolites showed statistically significant differences (p < 0.05) comparing +EoE -PPI to -EoE -PPI, a few of which were upregulated metabolites involved in the urea cycle. There were 14 significant differences comparing +EoE +PPI to +EoE -PPI. This demonstrated that in EoE patients, PPI use upregulated metabolites involved in the urea cycle, while it downregulated metabolites involved in methylation. Comparison among all four groups, +EoE +PPI, +EoE -PPI, -EoE +PPI, and -EoE -PPI, revealed 27 significantly different metabolites. +EoE +PPI had downregulated methionine and N-acetyl methionine, while both +EoE groups and -EoE +PPI had upregulated homocysteine, N-acetylputrescine, N-acetylornithine, arginine, and ornithine. Conclusion: The present study revealed key plasma metabolite differences in children with EoE compared to unaffected controls. Notable candidate biomarkers include dimethylarginine, putrescine, and N-acetylputrescine. PPI use was shown to influence these urea cycle metabolites, regardless of EoE presence. Therefore, future studies should distinguish patients based on PPI use or determine metabolites while not on treatment. These findings will be confirmed in a larger validation phase, as this may represent a significant discovery in the search for a non-invasive biomarker for EoE. Clinical Trial Registration: This clinical trial was registered with ClinicalTrials.gov, identifier: NCT 03107819.

Monozygotic twins have long been studied to estimate heritability and explore epigenetic influences on phenotypic variation. The phenotypic and epigenetic similarities of monozygotic twins have been assumed to be largely due to their genetic identity.

Estrogen receptor α (ERα) is the major driving transcription factor in the mammary gland development as well as breast cancer initiation and progression. However, the genomic landscape of ERα binding sites in the normal mouse mammary gland has not been completely elucidated. Here, we mapped genome-wide ERα binding events by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) in the mouse mammary gland in response to estradiol. We identified 6237 high confidence ERα binding sites in two biological replicates and showed that many of these were located at distal enhancer regions. Furthermore, we discovered 3686 unique genes in the mouse genome that recruit ER in response to estradiol. Interrogation of ER-DNA binding sites in ER-positive luminal epithelial cells showed that the ERE, PAX2, SF1, and AP1 motifs were highly enriched at distal enhancer regions. In addition, comprehensive transcriptome analysis by RNA-seq revealed that 493 genes are differentially regulated by acute treatment with estradiol in the mouse mammary gland in vivo. Through integration of RNA-seq and ERα ChIP-seq data, we uncovered a novel ERα targetome in mouse mammary epithelial cells. Taken together, our study has identified the genomic landscape of ERα binding events in mouse mammary epithelial cells. Furthermore, our study also highlights the cis-regulatory elements and cofactors that are involved in estrogen signaling and may contribute to ductal elongation in the normal mouse mammary gland.

DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific.