CHIPIN: ChIP-seq inter-sample normalization based on signal invariance across transcriptionally constant genes.

BACKGROUND: Multiple studies rely on ChIP-seq experiments to assess the effect of gene modulation and drug treatments on protein binding and chromatin structure. However, most methods commonly used for the normalization of ChIP-seq binding intensity signals across conditions, e.g., the normalization to the same number of reads, either assume a constant signal-to-noise ratio across conditions or base the estimates of correction factors on genomic regions with intrinsically different signals between conditions. Inaccurate normalization of ChIP-seq signal may, in turn, lead to erroneous biological conclusions. RESULTS: We developed a new R package, CHIPIN, that allows normalizing ChIP-seq signals across different conditions/samples when spike-in information is not available, but gene expression data are at hand. Our normalization technique is based on the assumption that, on average, no differences in ChIP-seq signals should be observed in the regulatory regions of genes whose expression levels are constant across samples/conditions. In addition to normalizing ChIP-seq signals, CHIPIN provides as output a number of graphs and calculates statistics allowing the user to assess the efficiency of the normalization and qualify the specificity of the antibody used. In addition to ChIP-seq, CHIPIN can be used without restriction on open chromatin ATAC-seq or DNase hypersensitivity data. We validated the CHIPIN method on several ChIP-seq data sets and documented its superior performance in comparison to several commonly used normalization techniques. CONCLUSIONS: The CHIPIN method provides a new way for ChIP-seq signal normalization across conditions when spike-in experiments are not available. The method is implemented in a user-friendly R package available on GitHub: https://github.com/BoevaLab/CHIPIN.

Steroidogenic differentiation and PKA signaling are programmed by histone methyltransferase EZH2 in the adrenal cortex.

Adrenal cortex steroids are essential for body homeostasis, and adrenal insufficiency is a life-threatening condition. Adrenal endocrine activity is maintained through recruitment of subcapsular progenitor cells that follow a unidirectional differentiation path from zona glomerulosa to zona fasciculata (zF). Here, we show that this unidirectionality is ensured by the histone methyltransferase EZH2. Indeed, we demonstrate that EZH2 maintains adrenal steroidogenic cell differentiation by preventing expression of GATA4 and WT1 that cause abnormal dedifferentiation to a progenitor-like state in Ezh2 KO adrenals. EZH2 further ensures normal cortical differentiation by programming cells for optimal response to adrenocorticotrophic hormone (ACTH)/PKA signaling. This is achieved by repression of phosphodiesterases PDE1B, 3A, and 7A and of PRKAR1B. Consequently, EZH2 ablation results in blunted zF differentiation and primary glucocorticoid insufficiency. These data demonstrate an all-encompassing role for EZH2 in programming steroidogenic cells for optimal response to differentiation signals and in maintaining their differentiated state.

COUP-TFI modifies CXCL12 and CXCR4 expression by activating EGF signaling and stimulates breast cancer cell migration.

BACKGROUND: The orphan receptors COUP-TF (chicken ovalbumin upstream promoter transcription factor) I and II are members of the nuclear receptor superfamily that play distinct and critical roles in vertebrate organogenesis. The involvement of COUP-TFs in cancer development has recently been suggested by several studies but remains poorly understood. METHODS: MCF-7 breast cancer cells overexpressing COUP-TFI and human breast tumors were used to investigate the role of COUP-TFI in the regulation of CXCL12/CXCR4 signaling axis in relation to cell growth and migration. We used Immunofluorescence, western-blot, RT-PCR, Formaldehyde-assisted Isolation of Regulatory Elements (FAIRE) assays, as well as cell proliferation and migration assays. RESULTS: Previously, we showed that COUP-TFI expression is enhanced in breast cancer compared to normal tissue. Here, we report that the CXCL12/CXCR4 signaling pathway, a crucial pathway in cell growth and migration, is an endogenous target of COUP-TFI in breast cancer cells. The overexpression of COUP-TFI in MCF-7 cells inhibits the expression of the chemokine CXCL12 and markedly enhances the expression of its receptor, CXCR4. Our results demonstrate that the modification of CXCL12/CXCR4 expression by COUP-TFI is mediated by the activation of epithelial growth factor (EGF) and the EGF receptor. Furthermore, we provide evidence that these effects of COUP-TFI increase the growth and motility of MCF-7 cells in response to CXCL12. Cell migration toward a CXCL12 gradient was inhibited by AMD3100, a specific antagonist of CXCR4, or in the presence of excess CXCL12 in the cell culture medium. The expression profiles of CXCR4, CXCR7, CXCL12, and COUP-TFI mRNA in 82 breast tumors and control non-tumor samples were measured using real-time PCR. CXCR4 expression was found to be significantly increased in the tumors and correlated with the tumor grade, whereas the expression of CXCL12 was significantly decreased in the tumors compared with the healthy samples. Significantly higher COUP-TFI mRNA expression was also detected in grade 1 tumors. CONCLUSIONS: Together, our mechanistic in vitro assays and in vivo results suggest that a reduction in chemokine CXCL12 expression, with an enhancement of CXCR4 expression, provoked by COUP-TFI, could be associated with an increase in the invasive potential of breast cancer cells.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs.

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Estrogen represses CXCR7 gene expression by inhibiting the recruitment of NFκB transcription factor at the CXCR7 promoter in breast cancer cells.

Although many studies reported mechanisms involved in the positive regulation of estrogens (E2) target genes, very little is known concerning the repressive effect of E2. In this study, we explored the molecular mechanisms by which E2 regulates CXCR7 expression in breast cancer cells. Our results show that E2-mediated down-regulation of CXCR7 occurs at the transcriptional level as demonstrated using actinomycin D and requires estrogen receptor alpha (ERα). In addition, CXCR7 is a primary ERα-target gene because the effect of E2 does not require the synthesis of an intermediary protein as revealed by the translational inhibitor cycloheximide treatment. Using an inhibitor of the NFκB pathway and chromatin immunoprecipitation assays, we demonstrated that NFκB is necessary for the high expression of CXCR7 gene and is recruited to the proximal promoter of the CXCR7 gene. Interestingly, the chromatin immunoprecipitation analyses also showed that E2-treatment significantly prevented the recruitment of NFκB to the promoter. Altogether, our results demonstrate that E2, through ERα, directly down-regulates CXCR7 expression by interfering with NFκB transcription factor at the promoter level.

DNA hypermethylation driven by DNMT1 and DNMT3A favors tumor immune escape contributing to the aggressiveness of adrenocortical carcinoma.

BACKGROUND: Adrenocortical carcinoma is rare and aggressive endocrine cancer of the adrenal gland. Within adrenocortical carcinoma, a recently described subtype characterized by a CpG island methylator phenotype (CIMP) has been associated with an especially poor prognosis. However, the drivers of CIMP remain unknown. Furthermore, the functional relation between CIMP and poor clinical outcomes of patients with adrenocortical carcinoma stays elusive. RESULTS: Here, we show that CIMP in adrenocortical carcinoma is linked to the increased expression of DNA methyltransferases DNMT1 and DNMT3A driven by a gain of gene copy number and cell hyperproliferation. Importantly, we demonstrate that CIMP contributes to tumor aggressiveness by favoring tumor immune escape. This effect could be at least partially reversed by treatment with the demethylating agent 5-azacytidine. CONCLUSIONS: In sum, our findings suggest that co-treatment with demethylating agents might enhance the efficacy of immunotherapy and could represent a novel therapeutic approach for patients with high CIMP adrenocortical carcinoma.

Chromatin Immunoprecipitation Followed by Next-Generation Sequencing (ChIP-Seq) Analysis in Ewing Sarcoma.

ChIP-seq is the method of choice for profiling protein-DNA interactions, and notably for characterizing the landscape of transcription factor binding and histone modifications. This technique has been widely used to study numerous aspects of tumor biology and led to the development of several promising cancer therapies. In Ewing sarcoma research, ChIP-seq provided important insights into the mechanism of action of the major oncogenic fusion protein EWSR1-FLI1 and related epigenetic and transcriptional changes. In this chapter, we provide a detailed pipeline to analyze ChIP-seq experiments from the preprocessing of raw data to tertiary analysis of detected binding sites. We also advise on best practice to prepare tumor samples prior to sequencing.

Opposite T3 Response of ACTG1-FOS Subnetwork Differentiate Tailfin Fate in Xenopus Tadpole and Post-hatching Axolotl.

Amphibian post-embryonic development and Thyroid Hormones (TH) signaling are deeply and intimately connected. In anuran amphibians, TH induce the spectacular and complex process known as metamorphosis. In paedomorphic salamanders, at similar development time, raising levels of TH fail to induce proper metamorphosis, as many "larval" tissues (e.g., gills, tailfin) are maintained. Why does the same evolutionary conserved signaling pathway leads to alternative phenotypes? We used a combination of developmental endocrinology, functional genomics and network biology to compare the transcriptional response of tailfin to TH, in the post-hatching paedormorphic Axolotl salamander and Xenopus tadpoles. We also provide a technological framework that efficiently reduces large lists of regulated genes down to a few genes of interest, which is well-suited to dissect endocrine regulations. We first show that Axolotl tailfin undergoes a strong and robust TH-dependent transcriptional response at post embryonic transition, despite the lack of visible anatomical changes. We next show that Fos and Actg1, which structure a single and dense subnetwork of cellular sensors and regulators, display opposite regulation between the two species. We finally show that TH treatments and natural variations of TH levels follow similar transcriptional dynamics. We suggest that, at the molecular level, tailfin fate correlates with the alternative transcriptional states of an fos-actg1 sub-network, which also includes transcription factors and regulators of cell fate. We propose that this subnetwork is one of the molecular switches governing the initiation of distinct TH responses, with transcriptional programs conducting alternative tailfin fate (maintenance vs. resorption) 2 weeks post-hatching.

Significant random signatures reveals new biomarker for breast cancer.

BACKGROUND: In 2012, Venet et al. proposed that at least in the case of breast cancer, most published signatures are not significantly more associated with outcome than randomly generated signatures. They suggested that nominal p-value is not a good estimator to show the significance of a signature. Therefore, one can reasonably postulate that some information might be present in such significant random signatures. METHODS: In this research, first we show that, using an empirical p-value, these published signatures are more significant than their nominal p-values. In other words, the proposed empirical p-value can be considered as a complimentary criterion for nominal p-value to distinguish random signatures from significant ones. Secondly, we develop a novel computational method to extract information that are embedded within significant random signatures. In our method, a score is assigned to each gene based on the number of times it appears in significant random signatures. Then, these scores are diffused through a protein-protein interaction network and a permutation procedure is used to determine the genes with significant scores. The genes with significant scores are considered as the set of significant genes. RESULTS: First, we applied our method on the breast cancer dataset NKI to achieve a set of significant genes in breast cancer considering significant random signatures. Secondly, prognostic performance of the computed set of significant genes is evaluated using DMFS and RFS datasets. We have observed that the top ranked genes from this set can successfully separate patients with poor prognosis from those with good prognosis. Finally, we investigated the expression pattern of TAT, the first gene reported in our set, in malignant breast cancer vs. adjacent normal tissue and mammospheres. CONCLUSION: Applying the method, we found a set of significant genes in breast cancer, including TAT, a gene that has never been reported as an important gene in breast cancer. Our results show that the expression of TAT is repressed in tumors suggesting that this gene could act as a tumor suppressor in breast cancer and could be used as a new biomarker.

De Novo Transcriptomic Approach to Study Thyroid Hormone Receptor Action in Non-mammalian Models.

Thyroid hormones are pleiotropic hormones involved in chordates physiology. Understanding their functions and mechanisms is also instrumental to diagnose dys-regulations and get a predictive power that can applied to medicine, ecology, etc. Today, high-throughput sequencing technologies offer the opportunity to address this issue not only in model organisms but also in non-model organisms. Here, we describe a method that makes use of RNA-seq to address differential expression analysis in non-model organism.

Rapid assessment of estrogenic compounds by CXCL-test illustrated by the screening of the UV-filter derivative benzophenones.

CXCL-test is a method that uses the estrogen-dependent secretion of the natural endogenous chemokine CXCL12 to evaluate the estrogenic activity of molecules. CXCL12 chemokine is involved in the estrogen dependent proliferation of breast cancer cells. Its measure is an indicator of cell proliferation and is used as an alternative test to classical proliferation test. Here we aimed to optimize this test, first to increase the number of tested molecules in a single assay and then to decrease the number of intermediate steps. The optimized CXCL-test was finally used for the evaluation of the estrogenic potency of emerging chemical pollutants: the UV filter benzophenones (BPs). The effect of BPs on CXCL12 secretion was also validated by real time quantitative RT-PCR. The optimized CXCL-test allowed a fast and direct assessment of estrogenic potency of molecules. The estrogenic activities of benzophenones were characterized and divided in two groups. The first one contains weak estrogenic compounds (BP, BP1, BP2, BP3, 234BP and 2344'BP). The second one contains medium estrogenic compounds (4BP, 44'BP, BP8, THB).

Activation of the MKL1/actin signaling pathway induces hormonal escape in estrogen-responsive breast cancer cell lines.

Estrogen receptor alpha (ERα) is generally considered to be a good prognostic marker because almost 70% of ERα-positive tumors respond to anti-hormone therapies. Unfortunately, during cancer progression, mammary tumors can escape from estrogen control, resulting in resistance to treatment. In this study, we demonstrate that activation of the actin/megakaryoblastic leukemia 1 (MKL1) signaling pathway promotes the hormonal escape of estrogen-sensitive breast cancer cell lines. The actin/MKL1 signaling pathway is silenced in differentiated ERα-positive breast cancer MCF-7 and T47D cell lines and active in ERα-negative HMT-3522 T4-2 and MDA-MB-231 breast cancer cells, which have undergone epithelial-mesenchymal transition. We showed that MKL1 activation in MCF-7 cells, either by modulating actin dynamics or using MKL1 mutants, down-regulates ERα expression and abolishes E2-dependent cell growth. Interestingly, the constitutively active form of MKL1 represses PR and HER2 expression in these cells and increases the expression of HB-EGF, TGFß, and amphiregulin growth factors in an E2-independent manner. The resulting expression profile (ER-, PR-, HER2-) typically corresponds to the triple-negative breast cancer expression profile.

Modulation of estrogen receptor alpha activity and expression during breast cancer progression.

Seventy percent of breast tumors express the estrogen receptor (ER), which is generally considered to predict a better outcome relative to ER-negative tumors, as they often respond to antiestrogen therapies. During cancer progression, mammary tumors can escape from estrogen control, resulting in the acquisition of invasive properties and resistance to treatment. ER expression is a dynamic phenomenon and is finely regulated at numerous levels, including the gene, mRNA, and protein levels. As a consequence, many molecular mechanisms have been implicated in modulating ER activity and estrogen signaling in mammary cancer. In fact, one-third of ER-positive breast cancer cells do not respond to first-line endocrine therapies, and a large subset of relapsing tumors retain ER expression. Increased knowledge of these mechanisms has led to the development of better prognostic methods and targeted therapies for patients; however, additional research is still needed to improve patient survival. In this chapter, we focus on the signaling pathways leading to changes in or loss of ER activity in breast cancer progression.

Assessment and molecular actions of endocrine-disrupting chemicals that interfere with estrogen receptor pathways.

In all vertebrate species, estrogens play a crucial role in the development, growth, and function of reproductive and nonreproductive tissues. A large number of natural or synthetic chemicals present in the environment and diet can interfere with estrogen signaling; these chemicals are called endocrine disrupting chemicals (EDCs) or xenoestrogens. Some of these compounds have been shown to induce adverse effects on human and animal health, and some compounds are suspected to contribute to diverse disease development. Because xenoestrogens have varying sources and structures and could act in additive or synergistic effects when combined, they have multiple mechanisms of action. Consequently, an important panel of in vivo and in vitro bioassays and chemical analytical tools was used to screen, evaluate, and characterize the potential impacts of these compounds on humans and animals. In this paper, we discuss different molecular actions of some of the major xenoestrogens found in food or the environment, and we summarize the current models used to evaluate environmental estrogens.

Estrogenic potency of benzophenone UV filters in breast cancer cells: proliferative and transcriptional activity substantiated by docking analysis.

The results from recent studies show that some benzophenones (BPs) and their hydroxylated metabolites can function as weak estrogens (E2) in the environment. However, little is known about the structure-activity relationship of these molecules. We have examined the effects of exposure to ten different BPs on the proliferation of estrogen receptor (ER)-positive breast cancer cells and on the transcriptional activity of E2-target genes. We analyzed two genes that are tightly linked with estrogen-mediated proliferation, the CXCL12 and amphiregulin genes and two classical estrogen-responsive genes, the pS2 and progesterone receptor. Significant differences in the BPs efficiency to induce cell proliferation and endogenous E2-target gene expressions were observed. Using ERE-, Sp1-, AP1- and C3-reporter genes that contain different ER-binding sites in their promoter, we also showed significant differences in the BPs efficiency in activation of the ER transactivation. Together, our analyzes showed that the most active molecule is 4-hydroxy-BP. Docking analysis of the interaction of BPs in the ligand-binding pocket of ERα suggests that the minimum structural requirement for the estrogenic activity of BPs is a hydroxyl (OH) group in the phenyl A-ring that allows interaction with Glu-353, Arg-394 or Phe-404, which enhances the stability between BPs and ERα. Our modeling also indicates a loss of interaction between the OH groups of the phenyl B-ring and His-524. In addition, the presence of some OH groups in the phenyl B-ring can create repulsion forces, which may constrain helix 12 in an unfavorable position, explaining the differential estrogenic effects of BPs. These results, together with our analysis of BPs for their potency in activation of cell proliferation and ER-mediated transcription, report an improved understanding of the mechanism and structure-activity relationship of BPs.

Differential estrogen-regulation of CXCL12 chemokine receptors, CXCR4 and CXCR7, contributes to the growth effect of estrogens in breast cancer cells.

CXCR4 and CXCR7 are the two receptors for the chemokine CXCL12, a key mediator of the growth effect of estrogens (E2) in estrogen receptor (ER)-positive breast cancers. In this study we examined E2-regulation of the CXCL12 axis components and their involvement in the growth of breast cancer cells. CXCR4 and CXCR7 were differentially regulated by E2 which enhanced the expression of both CXCL12 and CXCR4 but repressed the expression of CXCR7. Formaldehyde-associated isolation of regulatory elements (FAIRE) revealed that E2-mediated transcriptional regulation of these genes is linked to the control of the compaction state of chromatin at their promoters. This effect could be accomplished via several distal ER-binding sites in the regions surrounding these genes, all of which are located 20-250 kb from the transcription start site. Furthermore, individual down-regulation of CXCL12, CXCR4 or CXCR7 expression as well as the inhibition of their activity significantly decreases the rate of basal cell growth. In contrast, E2-induced cell growth was differentially affected. Unlike CXCR7, the inhibition of the expression or activity of either CXCL12 or CXCR4 significantly blunted the E2-mediated stimulation of cellular growth. Besides, CXCR7 over-expression increased the basal MCF-7 cell growth rate and decreased the growth effect of E2. These findings indicate that E2 regulation of the CXCL12 signaling axis is important for the E2-mediated growth effect of breast cancer cells. These data also provide support for distinct biological functions of CXCR4 and CXCR7 and suggest that targeting CXCR4 and/or CXCR7 would have distinct molecular effects on ER-positive breast tumors.