Genome-wide location analysis reveals distinct transcriptional circuitry by paralogous regulators Foxa1 and Foxa2

Interrogating Estrogen Signaling Pathways in Human ER-Positive Breast Cancer Cells Forming Bone Metastases in Mice

Breast cancer bone metastases (BMET) are incurable, primarily osteolytic, and occur most commonly in estrogen receptor-α positive (ER+) breast cancer. ER+ human breast cancer BMET modeling in mice has demonstrated an estrogen (E2)-dependent increase in tumor-associated osteolysis and bone-resorbing osteoclasts, independent of estrogenic effects on tumor proliferation or bone turnover, suggesting a possible mechanistic link between tumoral ERα-driven osteolysis and ER+ bone progression. To explore this question, inducible secretion of the osteolytic factor, parathyroid hormone-related protein (PTHrP), was utilized as an in vitro screening bioassay to query the osteolytic potential of estrogen receptor- and signaling pathway-specific ligands in BMET-forming ER+ human breast cancer cells expressing ERα, ERß, and G protein-coupled ER. After identifying genomic ERα signaling, also responsibility for estrogen's proliferative effects, as necessary and sufficient for osteolytic PTHrP secretion, in vivo effects of a genomic-only ER agonist, estetrol (E4), on osteolytic ER+ BMET progression were examined. Surprisingly, while pharmacologic effects of E4 on estrogen-dependent tissues, including bone, were evident, E4 did not support osteolytic BMET progression (vs robust E2 effects), suggesting an important role for nongenomic ER signaling in ER+ metastatic progression at this site. Because bone effects of E4 did not completely recapitulate those of E2, the relative importance of nongenomic ER signaling in tumor vs bone cannot be ascertained here. Nonetheless, these intriguing findings suggest that targeted manipulation of estrogen signaling to mitigate ER+ metastatic progression in bone may require a nuanced approach, considering genomic and nongenomic effects of ER signaling on both sides of the tumor/bone interface.

FOXA1 forms biomolecular condensates that unpack condensed chromatin to function as a pioneer factor

Eukaryotic DNA is packaged into chromatin in the nucleus, restricting the binding of transcription factors (TFs) to their target DNA sites. FOXA1 functions as a pioneer TF to bind condensed chromatin and initiate the opening of local chromatin for gene expression. However, the principles of FOXA1 recruitment and how it subsequently unpacks the condensed chromatin remain elusive. Here, we revealed that FOXA1 intrinsically forms submicron-sized condensates through its N- and C-terminal intrinsically disordered regions (IDRs). Notably, both IDRs enable FOXA1 to dissolve the condensed chromatin. In addition, the DNA-binding capacity of FOXA1 contributes to its ability to both form condensates and dissolve condensed chromatin. Further genome-wide investigation showed that IDRs enable FOXA1 to bind and unpack the condensed chromatin to regulate the proliferation and migration of breast cancer cells. This work provides a principle of how pioneer TFs function to initiate competent chromatin states using their IDRs.

Sequential ChIP-Seq

ChIP-Seq has been used extensively to profile genome-wide transcription factor binding and post-translational histone modifications. A sequential ChIP assay determines the in vivo co-localization of two proteins to the same genomic locus. In this chapter, we combine the two protocols in Sequential ChIP-Seq, a method for identifying genome-wide sites of in vivo protein co-occupancy.

regioneReloaded: evaluating the association of multiple genomic region sets

MOTIVATION

Next-generation sequencing methods continue improving the annotation of genomes in part by determining the distribution of features such as epigenetic marks. Evaluating and interpreting the association between genomic regions and their features has become a common and challenging analysis in genomic and epigenomic studies.

RESULTS

With regioneR we provided an R package allowing to assess the statistical significance of pairwise associations between genomic region sets using permutation tests. We now present the R package regioneReloaded that builds upon regioneR's statistical foundation and extends the functionality for the simultaneous analysis and visualization of the associations between multiple genomic region sets. Thus, we provide a novel discovery tool for the identification of significant associations that warrant to be tested for functional interdependence.

AVAILABILITY AND IMPLEMENTATION

regioneReloaded is an R package released under an Artistic-2.0 License. The source code and documentation are freely available through Bioconductor: http://www.bioconductor.org/packages/regioneReloaded.

Joint sequence & chromatin neural networks characterize the differential abilities of Forkhead transcription factors to engage inaccessible chromatin

The DNA-binding activities of transcription factors (TFs) are influenced by both intrinsic sequence preferences and extrinsic interactions with cell-specific chromatin landscapes and other regulatory proteins. Disentangling the roles of these binding determinants remains challenging. For example, the FoxA subfamily of Forkhead domain (Fox) TFs are known pioneer factors that can bind to relatively inaccessible sites during development. Yet FoxA TF binding also varies across cell types, pointing to a combination of intrinsic and extrinsic forces guiding their binding. While other Forkhead domain TFs are often assumed to have pioneering abilities, how sequence and chromatin features influence the binding of related Fox TFs has not been systematically characterized. Here, we present a principled approach to compare the relative contributions of intrinsic DNA sequence preference and cell-specific chromatin environments to a TF's DNA-binding activities. We apply our approach to investigate how a selection of Fox TFs (FoxA1, FoxC1, FoxG1, FoxL2, and FoxP3) vary in their binding specificity. We over-express the selected Fox TFs in mouse embryonic stem cells, which offer a platform to contrast each TF's binding activity within the same preexisting chromatin background. By applying a convolutional neural network to interpret the Fox TF binding patterns, we evaluate how sequence and preexisting chromatin features jointly contribute to induced TF binding. We demonstrate that Fox TFs bind different DNA targets, and drive differential gene expression patterns, even when induced in identical chromatin settings. Despite the association between Forkhead domains and pioneering activities, the selected Fox TFs display a wide range of affinities for preexiting chromatin states. Using sequence and chromatin feature attribution techniques to interpret the neural network predictions, we show that differential sequence preferences combined with differential abilities to engage relatively inaccessible chromatin together explain Fox TF binding patterns at individual sites and genome-wide.

Pioneer Factor Foxa2 Mediates Chromatin Conformation Changes for Activation of Bile Acid Targets of FXR

BACKGROUND & AIMS

Transcription factors regulate gene expression that orchestrates liver physiology. Many bind at distal enhancers and chromatin looping is required to activate their targets. Chromatin architecture has been linked to essential functions of the liver, including metabolism and sexually dimorphic gene expression. We have previously shown that pioneer factor Foxa2 opens chromatin for binding of nuclear receptors farnesoid X receptor (FXR) and liver X receptor-α during acute ligand activation. FXR is activated by bile acids and deletion of Foxa2 in the liver results in intrahepatic cholestasis. We hypothesized that Foxa2 also enables chromatin conformational changes during ligand activation and performed genome-wide studies to test this hypothesis.

METHODS

We performed Foxa2 HiChIP (Hi-C and ChIP) to assess Foxa2-dependent long-range interactions in mouse livers treated with either vehicle control or FXR agonist GW4064.

RESULTS

HiChIP contact analysis shows that global chromatin interactions are dramatically increased during FXR activation. Ligand-treated livers exhibit extensive redistribution of topological associated domains and substantial increase in Foxa2-anchored loops, suggesting Foxa2 is involved in dynamic chromatin conformational changes. We demonstrate that chromatin conformation, including genome-wide interactions, topological associated domains, and intrachromosomal and interchromosomal Foxa2-anchored loops, drastically changes on addition of FXR agonist. Additional Foxa2 binding in ligand-activated state leads to formation of Foxa2-anchored loops, leading to distal interactions and activation of gene expression of FXR targets.

CONCLUSIONS

Ligand activation of FXR, and likely of related receptors, requires global changes in chromatin architecture. We determine a novel role for Foxa2 in enabling these conformational changes, extending its function in bile acid metabolism.

FOXA1/2 depletion drives global reprogramming of differentiation state and metabolism in a human liver cell line and inhibits differentiation of human stem cell-derived hepatic progenitor cells

FOXA factors are critical members of the developmental gene regulatory network (GRN) composed of master transcription factors (TF) which regulate murine cell fate and metabolism in the gut and liver. How FOXA factors dictate human liver cell fate, differentiation, and simultaneously regulate metabolic pathways is poorly understood. Here, we aimed to determine the role of FOXA2 (and FOXA1 which is believed to compensate for FOXA2) in controlling hepatic differentiation and cell metabolism in a human hepatic cell line (HepG2). siRNA mediated knockdown of FOXA1/2 in HepG2 cells significantly downregulated albumin (p < .05) and GRN TF gene expression (HNF4α, HEX, HNF1ß, TBX3) (p < .05) and significantly upregulated endoderm/gut/hepatic endoderm markers (goosecoid [GSC], FOXA3, and GATA4), gut TF (CDX2), pluripotent TF (NANOG), and neuroectodermal TF (PAX6) (p < .05), all consistent with partial/transient reprograming. shFOXA1/2 targeting resulted in similar findings and demonstrated evidence of reversibility of phenotype. RNA-seq followed by bioinformatic analysis of shFOXA1/2 knockdown HepG2 cells demonstrated 235 significant downregulated genes and 448 upregulated genes, including upregulation of markers for alternate germ layers lineages (cardiac, endothelial, muscle) and neurectoderm (eye, neural). We found widespread downregulation of glycolysis, citric acid cycle, mitochondrial genes, and alterations in lipid metabolism, pentose phosphate pathway, and ketogenesis. Functional metabolic analysis agreed with these findings, demonstrating significantly diminished glycolysis and mitochondrial respiration, with concomitant accumulation of lipid droplets. We hypothesized that FOXA1/2 inhibit the initiation of human liver differentiation in vitro. During human pluripotent stem cells (hPSC)-hepatic differentiation, siRNA knockdown demonstrated de-differentiation and unexpectedly, activation of pluripotency factors and neuroectoderm. shRNA knockdown demonstrated similar results and activation of SOX9 (hepatobiliary). These results demonstrate that FOXA1/2 controls hepatic and developmental GRN, and their knockdown leads to reprogramming of both differentiation and metabolism, with applications in studies of cancer, differentiation, and organogenesis.

FOXA1 is a transcriptional activator of Odf2/Cenexin and regulates primary ciliation

Primary cilia are sensory organelles essential for embryonic and postnatal development, and tissue homeostasis in adulthood. They are generated in a cell cycle-dependent manner and found on most cells of the body. Although cilia formation is intensively investigated virtually nothing is known about the transcriptional regulation of primary ciliation. We used here Odf2/Cenexin, encoding a protein of the mother centriole and the basal body that is mandatory for primary cilia formation, as the target gene for the identification of transcriptional activators. We identified a consensus binding site for Fox transcription factors (TFs) in its promoter region and focused here on the Fox family. We found transcriptional activation of Odf2 neither by FOXO TFs nor by the core TF for multiciliation, FOXJ1. However, we identified FOXA1 as a transcriptional activator of Odf2 by reporter gene assays and qRT-PCR, and showed by qWB that Foxa1 knockdown caused a decrease in ODF2 and CP110 proteins. We verified the binding sequence of FOXA1 in the Odf2 promoter by ChIP. Finally, we demonstrated that knockdown of FOXA1 affected primary cilia formation. We, thus, showed for the first time, that FOXA1 regulates primary ciliation by transcriptional activation of ciliary genes.

Paralogous Genes Involved in Embryonic Development: Lessons from the Eye and Other Tissues

During evolution, gene duplications lead to a naturally increased gene dosage. Duplicated genes can be further retained or eliminated over time by purifying selection pressure. The retention probability is increased by functional diversification and by the acquisition of novel functions. Interestingly, functionally diverged paralogous genes can maintain a certain level of functional redundancy and at least a partial ability to replace each other. In such cases, diversification probably occurred at the level of transcriptional regulation. Nevertheless, some duplicated genes can maintain functional redundancy after duplication and the ability to functionally compensate for the loss of each other. Many of them are involved in proper embryonic development. The development of particular tissues/organs and developmental processes can be more or less sensitive to the overall gene dosage. Alterations in the gene dosage or a decrease below a threshold level may have dramatic phenotypic consequences or even lead to embryonic lethality. The number of functional alleles of particular paralogous genes and their mutual cooperation and interactions influence the gene dosage, and therefore, these factors play a crucial role in development. This review will discuss individual interactions between paralogous genes and gene dosage sensitivity during development. The eye was used as a model system, but other tissues are also included.

Redistribution of lamina-associated domains reshapes binding of pioneer factor FOXA2 in development of nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is highly prevalent in type 2 diabetes mellitus and the elderly, impacting 40% of individuals over 70. Regulation of heterochromatin at the nuclear lamina has been associated with aging and age-dependent metabolic changes. We previously showed that changes at the lamina in aged hepatocytes and laminopathy models lead to redistribution of lamina-associated domains (LADs), opening of repressed chromatin, and up-regulation of genes regulating lipid synthesis and storage, culminating in fatty liver. Here, we test the hypothesis that change in the expression of lamina-associated proteins and nuclear shape leads to redistribution of LADs, followed by altered binding of pioneer factor FOXA2 and by up-regulation of lipid synthesis and storage, culminating in steatosis in younger NAFLD patients (aged 21-51). Changes in nuclear morphology alter LAD partitioning and reduced lamin B1 signal correlate with increased FOXA2 binding before severe steatosis in young mice placed on a western diet. Nuclear shape is also changed in younger NAFLD patients. LADs are redistrubted and lamin B1 signal decreases similarly in mild and severe steatosis. In contrast, FOXA2 binding is similar in normal and NAFLD patients with moderate steatosis and is repositioned only in NAFLD patients with more severe lipid accumulation. Hence, changes at the nuclear lamina reshape FOXA2 binding with progression of the disease. Our results suggest a role for nuclear lamina in etiology of NAFLD, irrespective of aging, with potential for improved stratification of patients and novel treatments aimed at restoring nuclear lamina function.

The circRNA CNEACR regulates necroptosis of cardiomyocytes through Foxa2 suppression

Circular RNAs (circRNAs) are differentially expressed in various cardiovascular disease including myocardial ischemia-reperfusion (I/R) injury. However, their functional impact on cardiomyocyte cell death, in particular, in necrotic forms of death remains elusive. In this study, we found that the level of mmu_circ_000338, a cardiac- necroptosis-associated circRNA (CNEACR), was reduced in hypoxia-reoxygenation (H/R) exposed cardiomyocytes and I/R-injured mice hearts. The enforced expression of CNEACR attenuated the necrotic form of cardiomyocyte death caused by H/R and suppressed of myocardial necrosis in I/R injured mouse heart, which was accompanied by a marked reduction of myocardial infarction size and improved cardiac function. Mechanistically, CNEACR directly binds to histone deacetylase (HDAC7) in the cytoplasm and interferes its nuclear entry. This leads to attenuation of HDAC7-dependent suppression of forkhead box protein A2 (Foxa2) transcription, which can repress receptor-interacting protein kinase 3 (Ripk3) gene by binding to its promoter region. In addition, CNEACR-mediated upregulation of FOXA2 inhibited RIPK3-dependent necrotic/necroptotic death of cardiomyocytes. Our study reveals that circRNAs such as CNEACR can regulate the cardiomyocyte necroptosis associated activity of HDACs, promotes cell survival and improves cardiac function in I/R-injured heart. Hence, the CNEACR/HDAC7/Foxa2/ RIPK3 axis could be an efficient target for alleviating myocardial damage caused by necroptotic death in ischemia heart diseases.

Comprehensive epigenomic profiling of human alveolar epithelial differentiation identifies key epigenetic states and transcription factor co-regulatory networks for maintenance of distal lung identity

BACKGROUND

Disruption of alveolar epithelial cell (AEC) differentiation is implicated in distal lung diseases such as chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, and lung adenocarcinoma that impact morbidity and mortality worldwide. Elucidating underlying disease pathogenesis requires a mechanistic molecular understanding of AEC differentiation. Previous studies have focused on changes of individual transcription factors, and to date no study has comprehensively characterized the dynamic, global epigenomic alterations that facilitate this critical differentiation process in humans.

RESULTS

We comprehensively profiled the epigenomic states of human AECs during type 2 to type 1-like cell differentiation, including the methylome and chromatin functional domains, and integrated this with transcriptome-wide RNA expression data. Enhancer regions were drastically altered during AEC differentiation. Transcription factor binding analysis within enhancer regions revealed diverse interactive networks with enrichment for many transcription factors, including NKX2-1 and FOXA family members, as well as transcription factors with less well characterized roles in AEC differentiation, such as members of the MEF2, TEAD, and AP1 families. Additionally, associations among transcription factors changed during differentiation, implicating a complex network of heterotrimeric complex switching in driving differentiation. Integration of AEC enhancer states with the catalog of enhancer elements in the Roadmap Epigenomics Mapping Consortium and Encyclopedia of DNA Elements (ENCODE) revealed that AECs have similar epigenomic structures to other profiled epithelial cell types, including human mammary epithelial cells (HMECs), with NKX2-1 serving as a distinguishing feature of distal lung differentiation.

CONCLUSIONS

Enhancer regions are hotspots of epigenomic alteration that regulate AEC differentiation. Furthermore, the differentiation process is regulated by dynamic networks of transcription factors acting in concert, rather than individually. These findings provide a roadmap for understanding the relationship between disruption of the epigenetic state during AEC differentiation and development of lung diseases that may be therapeutically amenable.

Suppression of Cell Tumorigenicity by Non-neural Pro-differentiation Factors via Inhibition of Neural Property in Tumorigenic Cells

Our studies have demonstrated that cell tumorigenicity and pluripotent differentiation potential stem from neural stemness or a neural ground state, which is defined by a regulatory network of higher levels of machineries for basic cell physiological functions, including cell cycle, ribosome biogenesis, protein translation, spliceosome, epigenetic modification factors, reprogramming factors, etc., in addition to the neural stemness specific factors. These machineries and neural stemness factors mostly play cancer-promoting roles. It can be deduced that differentiation requires the repression of neural ground state and causes the reduction or loss of neural ground state and thus tumorigenicity in tumorigenic cells. Formerly, we showed that neuronal differentiation led to reduced tumorigenicity in tumorigenic cells. In the present study, we show that non-neural pro-differentiation factors, such as GATA3, HNF4A, HHEX, and FOXA3 that specify mesodermal or/and endodermal tissues during vertebrate embryogenesis, suppress tumorigenicity via repression of neural stemness and promotion of non-neural property in tumorigenic cells. Mechanistically, these transcription factors repress the transcription of neural enriched genes and meanwhile activate genes that specify non-neural properties via direct binding to the promoters of these genes. We also show that combined expression of HHEX and FOXA3 suppresses tumorigenesis effectively in the AOM/DSS model of colitis-associated cancer. We suggest that targeting the property of neural stemness could be an effective strategy for cancer therapy.

Expression of cell adhesion molecule 1 in human and murine endometrial glandular cells and its increase during the proliferative phase by estrogen and cell density

AIMS

Cell adhesion molecule 1 (CADM1) mediates interepithelial adhesion and is upregulated in crowded epithelial monolayers. This study aimed to examine CADM1 expression in the human endometrium of proliferative and secretory phases, and its transcriptional regulation in terms of estrogen stimuli and higher cellularity.

MAIN METHODS

CADM1 immunohistochemistry was conducted on endometrial tissues from women in their 40s and adult mice subcutaneously injected with estradiol following ovariectomy. Dual-luciferase reporter assays were conducted using human endometrial HEC-50B and HEC-1B cells and reporter plasmids harboring the human CADM1 3.4-kb promoter and its deleted and mutated forms. Cells were transfected with estrogen receptor α cDNA and reporter plasmids, and treated with estradiol before luciferase activity measurement.

KEY FINDINGS

Immunohistochemistry revealed that CADM1 was clearly expressed on the lateral membranes of the simple columnar glandular cells in the proliferative phase, but not in the secretory phase, from both women and the mouse model. The glandular cell density increased two-fold in the proliferative phase. Reporter assays identified three Sp1-binding sites as estradiol-responsive elements in the proximal region (from -223 to -84) of the transcription start site (+1) in HEC-50B cells. When the cell culture was started at eight-fold higher cell density, the CADM1 3.4-kb promoter was transactivated at a two-fold higher level in HEC-50B cells. This cell density effect was not detected for the CADM1 2.3-kb or 1.6-kb promoter.

SIGNIFICANCE

Two (proximal and distal) promoter regions are suggested to function additively to transactivate CADM1 in endometrial glandular cells that crowd in the proliferative phase.

Resolving diverse protein-DNA footprints from exonuclease-based ChIP experiments

MOTIVATION

High-throughput chromatin immunoprecipitation (ChIP) sequencing-based assays capture genomic regions associated with the profiled transcription factor (TF). ChIP-exo is a modified protocol, which uses lambda exonuclease to digest DNA close to the TF-DNA complex, in order to improve on the positional resolution of the TF-DNA contact. Because the digestion occurs in the 5'-3' orientation, the protocol produces directional footprints close to the complex, on both sides of the double stranded DNA. Like all ChIP-based methods, ChIP-exo reports a mixture of different regions associated with the TF: those bound directly to the TF as well as via intermediaries. However, the distribution of footprints are likely to be indicative of the complex forming at the DNA.

RESULTS

We present ExoDiversity, which uses a model-based framework to learn a joint distribution over footprints and motifs, thus resolving the mixture of ChIP-exo footprints into diverse binding modes. It uses no prior motif or TF information and automatically learns the number of different modes from the data. We show its application on a wide range of TFs and organisms/cell-types. Because its goal is to explain the complete set of reported regions, it is able to identify co-factor TF motifs that appear in a small fraction of the dataset. Further, ExoDiversity discovers small nucleotide variations within and outside canonical motifs, which co-occur with variations in footprints, suggesting that the TF-DNA structural configuration at those regions is likely to be different. Finally, we show that detected modes have specific DNA shape features and conservation signals, giving insights into the structure and function of the putative TF-DNA complexes.

AVAILABILITY AND IMPLEMENTATION

The code for ExoDiversity is available on https://github.com/NarlikarLab/exoDIVERSITY.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

m6A RNA methylation-mediated HNF3γ reduction renders hepatocellular carcinoma dedifferentiation and sorafenib resistance

Hepatocyte nuclear factor 3γ (HNF3γ) is a hepatocyte nuclear factor, but its role and clinical significance in hepatocellular carcinoma (HCC) remain unclear. Herein, we report that HNF3γ expression is downregulated in patient HCC and inversely correlated with HCC malignancy and patient survival. Moreover, our data suggested that the HNF3γ reduction in HCC could be mediated by METTL14-dependent m6A methylation of HNF3γ mRNA. HNF3γ expression was increased during hepatic differentiation and decreased in dedifferentiated HCC cells. Interestingly, HNF3γ delivery promoted differentiation of not only HCC cells but also liver CSCs, which led to suppression of HCC growth. Mechanistic analysis suggested an HNF3γ-centered regulatory network that includes essential liver differentiation-associated transcription factors and functional molecules, which could synergistically facilitate HCC cell differentiation. More importantly, enforced HNF3γ expression sensitized HCC cells to sorafenib-induced growth inhibition and cell apoptosis through transactivation of OATP1B1 and OATP1B3 expression, which are major membrane transporters for sorafenib uptake. Clinical investigation showed that patient-derived HCC xenografts with high HNF3γ expression exhibited a sorafenib response and patients with high HCC HNF3γ levels benefited from sorafenib therapy. Together, these results suggest that HNF3γ plays an essential role in HCC differentiation and may serve as a therapeutic target and predictor of sorafenib benefit in patients.

Genome-Wide Analysis of the FOXA1 Transcriptional Network Identifies Novel Protein-Coding and Long Noncoding RNA Targets in Colorectal Cancer Cells

Differentiation status of tumors is correlated with metastatic potential and malignancy. FOXA1 (forkhead box A1) is a transcription factor known to regulate differentiation in certain tissues. Here, we investigate FOXA1 function in human colorectal cancer (CRC). We found that FOXA1 is robustly expressed in the normal human colon but significantly downregulated in colon adenocarcinoma. Applying FOXA1 chromatin immunoprecipitation coupled with deep sequencing and transcriptome analysis upon FOXA1 knockdown in well-differentiated CRC cells and FOXA1 overexpression in poorly differentiated CRC cells, we identified novel protein-coding and lncRNA genes regulated by FOXA1. Among the numerous novel FOXA1 targets we identified, we focused on CEACAM5, a tumor marker and facilitator of cell adhesion. We show that FOXA1 binds to a distal enhancer downstream of CEACAM5 and strongly activates its expression. Consistent with these data, we show that FOXA1 inhibits anoikis in CRC cells. Collectively, our results uncover novel protein-coding and noncoding targets of FOXA1 and suggest a vital role of FOXA1 in enhancing CEACAM5 expression and anoikis resistance in CRC cells.

Functional informed genome-wide interaction analysis of body mass index, diabetes and colorectal cancer risk

BACKGROUND

Body mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology.

METHODS

To improve statistical power and interpretation for gene-environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2 ) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype-Tissue Expression Project for all genes with heritability ≥1%. We used a mixed-effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI-CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2.

RESULTS

Among 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10-5 ), PSMC5 (P = 4.51 × 10-4 ) and CD33 (P = 2.71 × 10-4 ) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10-5 ) and SCN1B (P = 2.76 × 10-4 ) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10-5 ).

CONCLUSIONS

Aggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk.

FOXA2 controls the cis-regulatory networks of pancreatic cancer cells in a differentiation grade-specific manner

Differentiation of normal and tumor cells is controlled by regulatory networks enforced by lineage-determining transcription factors (TFs). Among them, TFs such as FOXA1/2 bind naïve chromatin and induce its accessibility, thus establishing new gene regulatory networks. Pancreatic ductal adenocarcinoma (PDAC) is characterized by the coexistence of well- and poorly differentiated cells at all stages of disease. How the transcriptional networks determining such massive cellular heterogeneity are established remains to be determined. We found that FOXA2, a TF controlling pancreas specification, broadly contributed to the cis-regulatory networks of PDACs. Despite being expressed in both well- and poorly differentiated PDAC cells, FOXA2 displayed extensively different genomic distributions and controlled distinct gene expression programs. Grade-specific functions of FOXA2 depended on its partnership with TFs whose expression varied depending on the differentiation grade. These data suggest that FOXA2 contributes to the regulatory networks of heterogeneous PDAC cells via interactions with alternative partner TFs.

LINC00261 Is an Epigenetically Regulated Tumor Suppressor Essential for Activation of the DNA Damage Response

Lung cancer is the leading cause of cancer-related death in the United States. Long noncoding RNAs (lncRNA) are a class of regulatory molecules whose role in lung carcinogenesis is poorly understood. In this study, we profiled lncRNA expression in lung adenocarcinoma (LUAD) cell lines, compared their expression with that of purified alveolar epithelial type II cells (the purported cell of origin for LUAD), cross-referenced these with lncRNAs altered in the primary human tumors, and interrogated for lncRNAs whose expression correlated with patient survival. We identified LINC00261, a lncRNA with unknown function in LUAD, adjacent to the pioneering transcription factor FOXA2. Loss of LINC00261 was observed in multiple tumor types, including liver, breast, and gastric cancer. Reintroduction of LINC00261 into human LUAD cell lines inhibited cell migration and slowed proliferation by inducing G₂-M cell-cycle arrest, while upregulating DNA damage pathway genes and inducing phosphorylation-mediated activation of components of the DNA damage pathway. FOXA2 was able to induce LINC00261 expression, and the entire locus underwent hypermethylation in LUAD, leading to loss of expression. We have thus identified an epigenetically deregulated lncRNA, whose loss of expression in LUAD promotes the malignant phenotype and blocks activation of the DNA damage machinery, predisposing lung cells to cancer development. SIGNIFICANCE: These findings identify LINC00261 as a tumor suppressor that blocks cellular proliferation by activating the DNA damage response and suggest that epigenetic therapy to inhibit DNA methylation may enhance treatment of LUAD. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/79/12/3050/F1.large.jpg.See related commentary by Davalos and Esteller, p. 3028.

TRAP-seq identifies cystine/glutamate antiporter as a driver of recovery from liver injury

Understanding the molecular basis of the regenerative response following hepatic injury holds promise for improved treatment of liver diseases. Here, we report an innovative method to profile gene expression specifically in the hepatocytes that regenerate the liver following toxic injury. We used the Fah-/- mouse, a model of hereditary tyrosinemia, which conditionally undergoes severe liver injury unless fumarylacetoacetate hydrolase (FAH) expression is reconstituted ectopically. We used translating ribosome affinity purification followed by high-throughput RNA sequencing (TRAP-seq) to isolate mRNAs specific to repopulating hepatocytes. We uncovered upstream regulators and important signaling pathways that are highly enriched in genes changed in regenerating hepatocytes. Specifically, we found that glutathione metabolism, particularly the gene Slc7a11 encoding the cystine/glutamate antiporter (xCT), is massively upregulated during liver regeneration. Furthermore, we show that Slc7a11 overexpression in hepatocytes enhances, and its suppression inhibits, repopulation following toxic injury. TRAP-seq allows cell type-specific expression profiling in repopulating hepatocytes and identified xCT, a factor that supports antioxidant responses during liver regeneration. xCT has potential as a therapeutic target for enhancing liver regeneration in response to liver injury.

Gene regulatory networks in differentiation and direct reprogramming of hepatic cells

Liver development proceeds by sequential steps during which gene regulatory networks (GRNs) determine differentiation and maturation of hepatic cells. Characterizing the architecture and dynamics of these networks is essential for understanding how cell fate decisions are made during development, and for recapitulating these processes during in vitro production of liver cells for toxicology studies, disease modelling and regenerative therapy. Here we review the GRNs that control key steps of liver development and lead to differentiation of hepatocytes and cholangiocytes in mammals. We focus on GRNs determining cell fate decisions and analyse subcircuitry motifs that may confer specific dynamic properties to the networks. Finally, we put our analysis in the perspective of recent attempts to directly reprogram cells to hepatocytes by forced expression of transcription factors.

DGCA: A comprehensive R package for Differential Gene Correlation Analysis

BACKGROUND

Dissecting the regulatory relationships between genes is a critical step towards building accurate predictive models of biological systems. A powerful approach towards this end is to systematically study the differences in correlation between gene pairs in more than one distinct condition.

RESULTS

In this study we develop an R package, DGCA (for Differential Gene Correlation Analysis), which offers a suite of tools for computing and analyzing differential correlations between gene pairs across multiple conditions. To minimize parametric assumptions, DGCA computes empirical p-values via permutation testing. To understand differential correlations at a systems level, DGCA performs higher-order analyses such as measuring the average difference in correlation and multiscale clustering analysis of differential correlation networks. Through a simulation study, we show that the straightforward z-score based method that DGCA employs significantly outperforms the existing alternative methods for calculating differential correlation. Application of DGCA to the TCGA RNA-seq data in breast cancer not only identifies key changes in the regulatory relationships between TP53 and PTEN and their target genes in the presence of inactivating mutations, but also reveals an immune-related differential correlation module that is specific to triple negative breast cancer (TNBC).

CONCLUSIONS

DGCA is an R package for systematically assessing the difference in gene-gene regulatory relationships under different conditions. This user-friendly, effective, and comprehensive software tool will greatly facilitate the application of differential correlation analysis in many biological studies and thus will help identification of novel signaling pathways, biomarkers, and targets in complex biological systems and diseases.

Changes in Gene Expression and Estrogen Receptor Cistrome in Mouse Liver Upon Acute E2 Treatment

Transcriptional regulation by the estrogen receptor-α (ER) has been investigated mainly in breast cancer cell lines, but estrogens such as 17β-estradiol (E2) exert numerous extrareproductive effects, particularly in the liver, where E2 exhibits both protective metabolic and deleterious thrombotic actions. To analyze the direct and early transcriptional effects of estrogens in the liver, we determined the E2-sensitive transcriptome and ER cistrome in mice after acute administration of E2 or placebo. These analyses revealed the early induction of genes involved in lipid metabolism, which fits with the crucial role of ER in the prevention of liver steatosis. Characterization of the chromatin state of ER binding sites (BSs) in mice expressing or not ER demonstrated that ER is not required per se for the establishment and/or maintenance of chromatin modifications at the majority of its BSs. This is presumably a consequence of a strong overlap between ER and hepatocyte nuclear factor 4α BSs. In contrast, 40% of the BSs of the pioneer factor forkhead box protein a (Foxa2) were dependent upon ER expression, and ER expression also affected the distribution of nucleosomes harboring dimethylated lysine 4 of Histone H3 around Foxa2 BSs. We finally show that, in addition to a network of liver-specific transcription factors including CCAAT/enhancer-binding protein and hepatocyte nuclear factor 4α, ER might be required for proper Foxa2 function in this tissue.

The Pioneer Transcription Factor FoxA Maintains an Accessible Nucleosome Configuration at Enhancers for Tissue-Specific Gene Activation

Nuclear DNA wraps around core histones to form nucleosomes, which restricts the binding of transcription factors to gene regulatory sequences. Pioneer transcription factors can bind DNA sites on nucleosomes and initiate gene regulatory events, often leading to the local opening of chromatin. However, the nucleosomal configuration of open chromatin and the basis for its regulation is unclear. We combined low and high levels of micrococcal nuclease (MNase) digestion along with core histone mapping to assess the nucleosomal configuration at enhancers and promoters in mouse liver. We find that MNase-accessible nucleosomes, bound by transcription factors, are retained more at liver-specific enhancers than at promoters and ubiquitous enhancers. The pioneer factor FoxA displaces linker histone H1, thereby keeping enhancer nucleosomes accessible in chromatin and allowing other liver-specific transcription factors to bind and stimulate transcription. Thus, nucleosomes are not exclusively repressive to gene regulation when they are retained with, and exposed by, pioneer factors.

Lhx1 functions together with Otx2, Foxa2, and Ldb1 to govern anterior mesendoderm, node, and midline development

Costello et al. demonstrate that Smad4/Eomes-dependent Lhx1 expression in the epiblast marks the entire definitive endoderm lineage, the anterior mesendoderm, and midline progenitors. In proteomic experiments, they characterize a complex comprised of Lhx1, Otx2, and Foxa2 as well as the chromatin-looping protein Ldb1.

Gene regulatory networks controlling functional activities of spatially and temporally distinct endodermal cell populations in the early mouse embryo remain ill defined. The T-box transcription factor Eomes, acting downstream from Nodal/Smad signals, directly activates the LIM domain homeobox transcription factor Lhx1 in the visceral endoderm. Here we demonstrate Smad4/Eomes-dependent Lhx1 expression in the epiblast marks the entire definitive endoderm lineage, the anterior mesendoderm, and midline progenitors. Conditional inactivation of Lhx1 disrupts anterior definitive endoderm development and impedes node and midline morphogenesis in part due to severe disturbances in visceral endoderm displacement. Transcriptional profiling and ChIP-seq (chromatin immunoprecipitation [ChIP] followed by high-throughput sequencing) experiments identified Lhx1 target genes, including numerous anterior definitive endoderm markers and components of the Wnt signaling pathway. Interestingly, Lhx1-binding sites were enriched at enhancers, including the Nodal-proximal epiblast enhancer element and enhancer regions controlling Otx2 and Foxa2 expression. Moreover, in proteomic experiments, we characterized a complex comprised of Lhx1, Otx2, and Foxa2 as well as the chromatin-looping protein Ldb1. These partnerships cooperatively regulate development of the anterior mesendoderm, node, and midline cell populations responsible for establishment of the left–right body axis and head formation.

The ratio of FoxA1 to FoxA2 in lung adenocarcinoma is regulated by LncRNA HOTAIR and chromatin remodeling factor LSH

The lncRNA HOTAIR is a critical regulator of cancer progression. Chromatin remodeling factor LSH is critical for normal development of plants and mammals. However, the underlying mechanisms causing this in cancer are not entirely clear. The functional diversification of the FOXA1 and FOXA2 contributes to the target genes during evolution and carcinogenesis. Little is known about the ratio of FOXA1 to FOXA2 in cancer. We here found that both HOTAIR and LSH overexpression was significantly correlated with poor survival in patients with lung adenocarcinoma cancer (ADC). Also, the ratio of FOXA1 and FOXA2 is linked with poor survival in patients with lung ADC. HOTAIR regulates the ratio of FOXA1 to FOXA2 and migration and invasion. HOTAIR and the ratio of FOXA1 to FOXA2 are negatively correlated. HOTAIR knockdown inhibits migration and invasion. HOTAIR is associated with LSH, and this association linked with the binding of LSH in the promoter of FOXA1, not FOXA2. Targeted inhibition of HOTAIR suppresses the migratory and invasive properties. These data suggest that HOTAIR is an important mediator of the ratio of FOXA1 and FOXA2 and LSH involves in, and suggest that HOTAIR inhibition may represent a promising therapeutic option for suppressing lung ADC progression.

Revising the embryonic origin of thyroid C cells in mice and humans

Current understanding infers a neural crest origin of thyroid C cells, the major source of calcitonin in mammals and ancestors to neuroendocrine thyroid tumors. The concept is primarily based on investigations in quail–chick chimeras involving fate mapping of neural crest cells to the ultimobranchial glands that regulate Ca²⁺ homeostasis in birds, reptiles, amphibians and fishes, but whether mammalian C cell development involves a homologous ontogenetic trajectory has not been experimentally verified. With lineage tracing, we now provide direct evidence that Sox17+ anterior endoderm is the only source of differentiated C cells and their progenitors in mice. Like many gut endoderm derivatives, embryonic C cells were found to coexpress pioneer factors forkhead box (Fox) a1 and Foxa2 before neuroendocrine differentiation takes place. In the ultimobranchial body epithelium emerging from pharyngeal pouch endoderm in early organogenesis, differential Foxa1/Foxa2 expression distinguished two spatially separated pools of C cell precursors with different growth properties. A similar expression pattern was recapitulated in medullary thyroid carcinoma cells in vivo, consistent with a growth-promoting role of Foxa1. In contrast to embryonic precursor cells, C cell-derived tumor cells invading the stromal compartment downregulated Foxa2, foregoing epithelial-to-mesenchymal transition designated by loss of E-cadherin; both Foxa2 and E-cadherin were re-expressed at metastatic sites. These findings revise mammalian C cell ontogeny, expand the neuroendocrine repertoire of endoderm and redefine the boundaries of neural crest diversification. The data further underpin distinct functions of Foxa1 and Foxa2 in both embryonic and tumor development.

Highlighted article: Mouse thyroid C cell precursors arise in foregut endoderm, and not the neural crest, disproving the current concept of a neural crest origin of thyroid neuroendocrine cells.

Chimeric MicroRNA-1291 Biosynthesized Efficiently in Escherichia coli Is Effective to Reduce Target Gene Expression in Human Carcinoma Cells and Improve Chemosensitivity

In contrast to the growing interests in studying noncoding RNAs (ncRNAs) such as microRNA (miRNA or miR) pharmacoepigenetics, there is a lack of efficient means to cost effectively produce large quantities of natural miRNA agents. Our recent efforts led to a successful production of chimeric pre-miR-27b in bacteria using a transfer RNA (tRNA)-based recombinant RNA technology, but at very low expression levels. Herein, we present a high-yield expression of chimeric pre-miR-1291 in common Escherichia coli strains using the same tRNA scaffold. The tRNA fusion pre-miR-1291 (tRNA/mir-1291) was then purified to high homogeneity using affinity chromatography, whose primary sequence and post-transcriptional modifications were directly characterized by mass spectrometric analyses. Chimeric tRNA/mir-1291 was readily processed to mature miR-1291 in human carcinoma MCF-7 and PANC-1 cells. Consequently, recombinant tRNA/mir-1291 reduced the protein levels of miR-1291 target genes, including ABCC1, FOXA2, and MeCP2, as compared with cells transfected with the same doses of control methionyl-tRNA scaffold with a sephadex aptamer (tRNA/MSA). In addition, tRNA-carried pre-miR-1291 suppressed the growth of MCF-7 and PANC-1 cells in a dose-dependent manner, and significantly enhanced the sensitivity of ABCC1-overexpressing PANC-1 cells to doxorubicin. These results indicate that recombinant miR-1291 agent is effective in the modulation of target gene expression and chemosensitivity, which may provide insights into high-yield bioengineering of new ncRNA agents for pharmacoepigenetics research.

FOXA2 regulates a network of genes involved in critical functions of human intestinal epithelial cells

The forkhead box A (FOXA) family of pioneer transcription factors is critical for the development of many endoderm-derived tissues. Their importance in regulating biological processes in the lung and liver is extensively characterized, though much less is known about their role in intestine. Here we investigate the contribution of FOXA2 to coordinating intestinal epithelial cell function using postconfluent Caco2 cells, differentiated into an enterocyte-like model. FOXA2 binding sites genome-wide were determined by ChIP-seq and direct targets of the factor were validated by ChIP-qPCR and siRNA-mediated depletion of FOXA1/2 followed by RT-qPCR. Peaks of FOXA2 occupancy were frequent at loci contributing to gene ontology pathways of regulation of cell migration, cell motion, and plasma membrane function. Depletion of both FOXA1 and FOXA2 led to a significant reduction in the expression of multiple transmembrane proteins including ion channels and transporters, which form a network that is essential for maintaining normal ion and solute transport. One of the targets was the adenosine A2B receptor, and reduced receptor mRNA levels were associated with a functional decrease in intracellular cyclic AMP. We also observed that 30% of FOXA2 binding sites contained a GATA motif and that FOXA1/A2 depletion reduced GATA-4, but not GATA-6 protein levels. These data show that FOXA2 plays a pivotal role in regulating intestinal epithelial cell function. Moreover, that the FOXA and GATA families of transcription factors may work cooperatively to regulate gene expression genome-wide in the intestinal epithelium.

YAP Inhibition Restores Hepatocyte Differentiation in Advanced HCC, Leading to Tumor Regression

Defective Hippo/YAP signaling in the liver results in tissue overgrowth and development of hepatocellular carcinoma (HCC). Here, we uncover mechanisms of YAP-mediated hepatocyte reprogramming and HCC pathogenesis. YAP functions as a rheostat in maintaining metabolic specialization, differentiation, and quiescence within the hepatocyte compartment. Increased or decreased YAP activity reprograms subsets of hepatocytes to different fates associated with deregulation of the HNF4A, CTNNB1, and E2F transcriptional programs that control hepatocyte quiescence and differentiation. Importantly, treatment with small interfering RNA-lipid nanoparticles (siRNA-LNPs) targeting YAP restores hepatocyte differentiation and causes pronounced tumor regression in a genetically engineered mouse HCC model. Furthermore, YAP targets are enriched in an aggressive human HCC subtype characterized by a proliferative signature and absence of CTNNB1 mutations. Thus, our work reveals Hippo signaling as a key regulator of the positional identity of hepatocytes, supports targeting of YAP using siRNA-LNPs as a paradigm of differentiation-based therapy, and identifies an HCC subtype that is potentially responsive to this approach.

Decoupling of evolutionary changes in transcription factor binding and gene expression in mammals

To understand the evolutionary dynamics between transcription factor (TF) binding and gene expression in mammals, we compared transcriptional output and the binding intensities for three tissue-specific TFs in livers from four closely related mouse species. For each transcription factor, TF-dependent genes and the TF binding sites most likely to influence mRNA expression were identified by comparing mRNA expression levels between wild-type and TF knockout mice. Independent evolution was observed genome-wide between the rate of change in TF binding and the rate of change in mRNA expression across taxa, with the exception of a small number of TF-dependent genes. We also found that binding intensities are preferentially conserved near genes whose expression is dependent on the TF, and the conservation is shared among binding peaks in close proximity to each other near the TSS. Expression of TF-dependent genes typically showed an increased sensitivity to changes in binding levels as measured by mRNA abundance. Taken together, these results highlight a significant tolerance to evolutionary changes in TF binding intensity in mammalian transcriptional networks and suggest that some TF-dependent genes may be largely regulated by a single TF across evolution.

The Populus ARBORKNOX1 homeodomain transcription factor regulates woody growth through binding to evolutionarily conserved target genes of diverse function

The class I KNOX homeodomain transcription factor ARBORKNOX1 (ARK1) is a key regulator of vascular cambium maintenance and cell differentiation in Populus. Currently, basic information is lacking concerning the distribution, functional characteristics, and evolution of ARK1 binding in the Populus genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) technology to identify ARK1 binding loci genome-wide in Populus. Computational analyses evaluated the distribution of ARK1 binding loci, the function of genes associated with bound loci, the effect of ARK1 binding on transcript levels, and evolutionary conservation of ARK1 binding loci. ARK1 binds to thousands of loci which are highly enriched proximal to the transcriptional start sites of genes of diverse functions. ARK1 target genes are significantly enriched in paralogs derived from the whole-genome salicoid duplication event. Both ARK1 and a maize (Zea mays) homolog, KNOTTED1, preferentially target evolutionarily conserved genes. However, only a small portion of ARK1 target genes are significantly differentially expressed in an ARK1 over-expression mutant. This study describes the functional characteristics and evolution of DNA binding by a transcription factor in an undomesticated tree, revealing complexities similar to those shown for transcription factors in model animal species.

Changes in nucleosome occupancy associated with metabolic alterations in aged mammalian liver

Aging is accompanied by physiological impairments, which, in insulin-responsive tissues, including the liver, predispose individuals to metabolic disease. However, the molecular mechanisms underlying these changes remain largely unknown. Here, we analyze genome-wide profiles of RNA and chromatin organization in the liver of young (3 months) and old (21 months) mice. Transcriptional changes suggest that derepression of the nuclear receptors PPARα, PPARγ, and LXRα in aged mouse liver leads to activation of targets regulating lipid synthesis and storage, whereas age-dependent changes in nucleosome occupancy are associated with binding sites for both known regulators (forkhead factors and nuclear receptors) and candidates associated with nuclear lamina (Hdac3 and Srf) implicated to govern metabolic function of aging liver. Winged-helix transcription factor Foxa2 and nuclear receptor corepressor Hdac3 exhibit a reciprocal binding pattern at PPARα targets contributing to gene expression changes that lead to steatosis in aged liver.

Strategies to fine-map genetic associations with lipid levels by combining epigenomic annotations and liver-specific transcription profiles

Characterization of the epigenome promises to yield the functional elements buried in the human genome sequence, thus helping to annotate non-coding DNA polymorphisms with regulatory functions. Here, we develop two novel strategies to combine epigenomic data with transcriptomic profiles in humans or mice to prioritize potential candidate SNPs associated with lipid levels by genome-wide association study (GWAS). First, after confirming that lipid-associated loci that are also expression quantitative trait loci (eQTL) in human livers are enriched for ENCODE regulatory marks in the human hepatocellular HepG2 cell line, we prioritize candidate SNPs based on the number of these marks that overlap the variant position. This method recognized the known SORT1 rs12740374 regulatory SNP associated with LDL-cholesterol, and highlighted candidate functional SNPs at 15 additional lipid loci. In the second strategy, we combine ENCODE chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-seq) data and liver expression datasets from knockout mice lacking specific transcription factors. This approach identified SNPs in specific transcription factor binding sites that are located near target genes of these transcription factors. We show that FOXA2 transcription factor binding sites are enriched at lipid-associated loci and experimentally validate that alleles of one such proxy SNP located near the FOXA2 target gene BIRC5 show allelic differences in FOXA2-DNA binding and enhancer activity. These methods can be used to generate testable hypotheses for many non-coding SNPs associated with complex diseases or traits.

Identification of transcription factors and single nucleotide polymorphisms of Lrh1 and its homologous genes in Lrh1-knockout pancreas of mice

BACKGROUND

To identify transcription factors (TFs) and single nucleotide polymorphisms (SNPs) of Lrh1 (also named Nr5a2) and its homologous genes in Lrh1-knockout pancreas of mice.

METHODS

The RNA-Seq data GSE34030 were downloaded from Gene Expression Omnibus (GEO) database, including 2 Lrh1 pancreas knockout samples and 2 wild type samples. All reads were processed through TopHat and Cufflinks package to calculate gene-expression level. Then, the differentially expressed genes (DEGs) were identified via non-parametric algorithm (NOISeq) methods in R package, of which the homology genes of Lrh1 were identified via BLASTN analysis. Furthermore, the TFs of Lrh1 and its homologous genes were selected based on TRANSFAC database. Additionally, the SNPs were analyzed via SAM tool to record the locations of mutant sites.

RESULTS

Total 15683 DEGs were identified, of which 23 was Lrh1 homology genes (3 up-regulated and 20 down-regulated). Fetoprotein TF (FTF) was the only TF of Lrh1 identified and the promoter-binding factor of FTF was CYP7A. The SNP annotations of Lrh1 homologous genes showed that 92% of the mutation sites were occurred in intron and upstream. Three SNPs of Lrh1 were located in intron, while 1819 SNPs of Phkb were located in intron and 1343 SNPs were located in the upstream region.

CONCLUSION

FTF combined with CYP7A might play an important role in Lrh1 regulated pancreas-specific transcriptional network. Furthermore, the SNPs analysis of Lrh1 and its homology genes provided the candidate mutant sites that might affect the Lrh1-related production and secretion of pancreatic fluid.

ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

Leveraging biological replicates to improve analysis in ChIP-seq experiments

ChIP-seq experiments identify genome-wide profiles of DNA-binding molecules including transcription factors, enzymes and epigenetic marks. Biological replicates are critical for reliable site discovery and are required for the deposition of data in the ENCODE and modENCODE projects. While early reports suggested two replicates were sufficient, the widespread application of the technique has led to emerging consensus that the technique is noisy and that increasing replication may be worthwhile. Additional biological replicates also allow for quantitative assessment of differences between conditions. To date it has remained controversial about how to confirm peak identification and to determine signal strength across biological replicates, particularly when the number of replicates is greater than two. Using objective metrics, we evaluate the consistency of biological replicates in ChIP-seq experiments with more than two replicates. We compare several approaches for binding site determination, including two popular but disparate peak callers, CisGenome and MACS2. Here we propose read coverage as a quantitative measurement of signal strength for estimating sample concordance. Determining binding based on genomic features, such as promoters, is also examined. We find that increasing the number of biological replicates increases the reliability of peak identification. Critically, binding sites with strong biological evidence may be missed if researchers rely on only two biological replicates. When more than two replicates are performed, a simple majority rule (>50% of samples identify a peak) identifies peaks more reliably in all biological replicates than the absolute concordance of peak identification between any two replicates, further demonstrating the utility of increasing replicate numbers in ChIP-seq experiments.

FOXA1 mediates p16(INK4a) activation during cellular senescence

Mechanisms governing the transcription of p16(INK4a), one of the master regulators of cellular senescence, have been extensively studied. However, little is known about chromatin dynamics taking place at its promoter and distal enhancer. Here, we report that Forkhead box A1 protein (FOXA1) is significantly upregulated in both replicative and oncogene-induced senescence, and in turn activates transcription of p16(INK4a) through multiple mechanisms. In addition to acting as a classic sequence-specific transcriptional activator, FOXA1 binding leads to a decrease in nucleosome density at the p16(INK4a) promoter in senescent fibroblasts. Moreover, FOXA1, itself a direct target of Polycomb-mediated repression, antagonizes Polycomb function at the p16(INK4a) locus. Finally, a systematic survey of putative FOXA1 binding sites in the p16(INK4a) genomic region revealed an ∼150 kb distal element that could loop back to the promoter and potentiate p16(INK4a) expression. Overall, our findings establish several mechanisms by which FOXA1 controls p16(INK4a) expression during cellular senescence.

Comprehensive meta-analysis of Signal Transducers and Activators of Transcription (STAT) genomic binding patterns discerns cell-specific cis-regulatory modules

BACKGROUND

Cytokine-activated transcription factors from the STAT (Signal Transducers and Activators of Transcription) family control common and context-specific genetic programs. It is not clear to what extent cell-specific features determine the binding capacity of seven STAT members and to what degree they share genetic targets. Molecular insight into the biology of STATs was gained from a meta-analysis of 29 available ChIP-seq data sets covering genome-wide occupancy of STATs 1, 3, 4, 5A, 5B and 6 in several cell types.

RESULTS

We determined that the genomic binding capacity of STATs is primarily defined by the cell type and to a lesser extent by individual family members. For example, the overlap of shared binding sites between STATs 3 and 5 in T cells is greater than that between STAT5 in T cells and non-T cells. Even for the top 1,000 highly enriched STAT binding sites, ~15% of STAT5 binding sites in mouse female liver are shared by other STATs in different cell types while in T cells ~90% of STAT5 binding sites are co-occupied by STAT3, STAT4 and STAT6. In addition, we identified 116 cis-regulatory modules (CRM), which are recognized by all STAT members across cell types defining a common JAK-STAT signature. Lastly, in liver STAT5 binding significantly coincides with binding of the cell-specific transcription factors HNF4A, FOXA1 and FOXA2 and is associated with cell-type specific gene transcription.

CONCLUSIONS

Our results suggest that genomic binding of STATs is primarily determined by the cell type and further specificity is achieved in part by juxtaposed binding of cell-specific transcription factors.