Quantitative Analysis of Cellular Morphology During In Vitro Decidualization

Decidualization is a differentiation process involving shape reorganization from a fibroblast to an epithelioid-like appearance characteristic of endometrial stromal cells. For the study of in vitro decidualization, one needs to check that the cells have undergone this process effectively. Verification is usually done by analyzing the expression of decidual markers, but changes in morphology are a more comprehensive feature. However, morphological specificities (i.e., flatness) of endometrial cells prevent the use of existing automated tools. A simple and accurate methodology was developed to quantify the phenotypic changes that occur in an in vitro decidualization system. This approach analyzes cell circularity directly from light microscopy images to follow the effects of progesterone or progestin R5020 in combination with estradiol (E2) and cAMP in inducing the decidualization of human endometrial cells. A statistical model to detect the differences in the kinetics of decidualization of the two hormonal stimuli before all the cell population acquire the decidual phenotype was implemented. It was found that statistical differences in morphology between decidualized and control cells could be detected 2 days after the treatments. Here we detail the model applied, scripts, and input files in order to provide a useful, practical, and low-cost tool to evaluate morphological aspects of endometrial stromal differentiation. This method allows the verification of the effectiveness of the decidualization process of the stromal endometrial cells without having to use cell replicates, as other methods such as immunofluorescence and RT-qPCR assays require. Consequently, this approach can follow the kinetics of a living single replicate throughout the experiment. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Cell circularity quantification of human stromal endometrial cells using ImageJ Basic Protocol 2: Statistical analysis of cell circularity of human stromal endometrial cells.

Coordinated changes in gene expression, H1 variant distribution and genome 3D conformation in response to H1 depletion

Up to seven members of the histone H1 family may contribute to chromatin compaction and its regulation in human somatic cells. In breast cancer cells, knock-down of multiple H1 variants deregulates many genes, promotes the appearance of genome-wide accessibility sites and triggers an interferon response via activation of heterochromatic repeats. However, how these changes in the expression profile relate to the re-distribution of H1 variants as well as to genome conformational changes have not been yet studied. Here, we combined ChIP-seq of five endogenous H1 variants with Chromosome Conformation Capture analysis in wild-type and H1.2/H1.4 knock-down T47D cells. The results indicate that H1 variants coexist in the genome in two large groups depending on the local GC content and that their distribution is robust with respect to H1 depletion. Despite the small changes in H1 variants distribution, knock-down of H1 translated into more isolated but de-compacted chromatin structures at the scale of topologically associating domains (TADs). Such changes in TAD structure correlated with a coordinated gene expression response of their resident genes. This is the first report describing simultaneous profiling of five endogenous H1 variants and giving functional evidence of genome topology alterations upon H1 depletion in human cancer cells.

In vivo temporal resolution of acute promyelocytic leukemia progression reveals a role of Klf4 in suppressing early leukemic transformation

In this study, Mas et al. used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. Their multiomics-integrated analysis identified Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis, and they characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα.

Genome organization plays a pivotal role in transcription, but how transcription factors (TFs) rewire the structure of the genome to initiate and maintain the programs that lead to oncogenic transformation remains poorly understood. Acute promyelocytic leukemia (APL) is a fatal subtype of leukemia driven by a chromosomal translocation between the promyelocytic leukemia (PML) and retinoic acid receptor α (RARα) genes. We used primary hematopoietic stem and progenitor cells (HSPCs) and leukemic blasts that express the fusion protein PML-RARα as a paradigm to temporally dissect the dynamic changes in the epigenome, transcriptome, and genome architecture induced during oncogenic transformation. We found that PML-RARα initiates a continuum of topologic alterations, including switches from A to B compartments, transcriptional repression, loss of active histone marks, and gain of repressive histone marks. Our multiomics-integrated analysis identifies Klf4 as an early down-regulated gene in PML-RARα-driven leukemogenesis. Furthermore, we characterized the dynamic alterations in the Klf4 cis-regulatory network during APL progression and demonstrated that ectopic Klf4 overexpression can suppress self-renewal and reverse the differentiation block induced by PML-RARα. Our study provides a comprehensive in vivo temporal dissection of the epigenomic and topological reprogramming induced by an oncogenic TF and illustrates how topological architecture can be used to identify new drivers of malignant transformation.

Chromatin topology defines estradiol-primed progesterone receptor and PAX2 binding in endometrial cancer cells

Estrogen (E2) and Progesterone (Pg), via their specific receptors (ERalpha and PR), are major determinants in the development and progression of endometrial carcinomas, However, their precise mechanism of action and the role of other transcription factors involved are not entirely clear. Using Ishikawa endometrial cancer cells, we report that E2 treatment exposes a set of progestin-dependent PR binding sites which include both E2 and progestin target genes. ChIP-seq results from hormone-treated cells revealed a non-random distribution of PAX2 binding in the vicinity of these estrogen-promoted PR sites. Altered expression of hormone regulated genes in PAX2 knockdown cells suggests a role for PAX2 in fine-tuning ERalpha and PR interplay in transcriptional regulation. Analysis of long-range interactions by Hi-C coupled with ATAC-seq data showed that these regions, that we call ‘progestin control regions’ (PgCRs), exhibited an open chromatin state even before hormone exposure and were non-randomly associated with regulated genes. Nearly 20% of genes potentially influenced by PgCRs were found to be altered during progression of endometrial cancer. Our findings suggest that endometrial response to progestins in differentiated endometrial tumor cells results in part from binding of PR together with PAX2 to accessible chromatin regions. What maintains these regions open remains to be studied.

STAG2 loss-of-function affects short-range genomic contacts and modulates the basal-luminal transcriptional program of bladder cancer cells

Cohesin exists in two variants containing STAG1 or STAG2. STAG2 is one of the most mutated genes in cancer and a major bladder tumor suppressor. Little is known about how its inactivation contributes to tumorigenesis. Here, we analyze the genomic distribution of STAG1 and STAG2 and perform STAG2 loss-of-function experiments using RT112 bladder cancer cells; we then analyze the genomic effects by integrating gene expression and chromatin interaction data. Functional compartmentalization exists between the cohesin complexes: cohesin-STAG2 displays a distinctive genomic distribution and mediates short and mid-ranged interactions that engage genes at higher frequency than those established by cohesin-STAG1. STAG2 knockdown results in down-regulation of the luminal urothelial signature and up-regulation of the basal transcriptional program, mirroring differences between STAG2-high and STAG2-low human bladder tumors. This is accompanied by rewiring of DNA contacts within topological domains, while compartments and domain boundaries remain refractive. Contacts lost upon depletion of STAG2 are assortative, preferentially occur within silent chromatin domains, and are associated with de-repression of lineage-specifying genes. Our findings indicate that STAG2 participates in the DNA looping that keeps the basal transcriptional program silent and thus sustains the luminal program. This mechanism may contribute to the tumor suppressor function of STAG2 in the urothelium.

MyoD induces ARTD1 and nucleoplasmic poly-ADP-ribosylation during fibroblast to myoblast transdifferentiation

While protein ADP-ribosylation was reported to regulate differentiation and dedifferentiation, it has so far not been studied during transdifferentiation. Here, we found that MyoD-induced transdifferentiation of fibroblasts to myoblasts promotes the expression of the ADP-ribosyltransferase ARTD1. Comprehensive analysis of the genome architecture by Hi-C and RNA-seq analysis during transdifferentiation indicated that ARTD1 locally contributed to A/B compartmentalization and coregulated a subset of MyoD target genes that were however not sufficient to alter transdifferentiation. Surprisingly, the expression of ARTD1 was accompanied by the continuous synthesis of nuclear ADP ribosylation that was neither dependent on the cell cycle nor induced by DNA damage. Conversely to the H₂O₂-induced ADP-ribosylation, the MyoD-dependent ADP-ribosylation was not associated to chromatin but rather localized to the nucleoplasm. Together, these data describe a MyoD-induced nucleoplasmic ADP-ribosylation that is observed particularly during transdifferentiation and thus potentially expands the plethora of cellular processes associated with ADP-ribosylation.

•

MyoD-dependent transdifferentiation of IMR90 to myoblasts induces ARTD1 expression

•

Transdifferentiation induces nuclear ARTD1-dependent ADP-ribosylation in myoblasts

•

This ADP-ribosylation is induced independent of cell cycle and of DNA damage

•

ARTD1-mediated poly-ADP-ribosylation localizes to the nucleoplasm in myoblasts

90 YEARS OF PROGESTERONE: Molecular mechanisms of progesterone receptor action on the breast cancer genome

Gene regulation by steroid hormones has been at the forefront in elucidating the intricacies of transcriptional regulation in eukaryotes ever since the discovery by Karlson and Clever that the insect steroid hormone ecdysone induces chromatin puffs in giant chromosomes. After the successful cloning of the hormone receptors toward the end of the past century, detailed mechanistic insight emerged in some model systems, in particular the MMTV provirus. With the arrival of next generation DNA sequencing and the omics techniques, we have gained even further insight into the global cellular response to steroid hormones that in the past decades also extended to the function of the 3D genome topology. More recently, advances in high resolution microcopy, single cell genomics and the new vision of liquid-liquid phase transitions in the context of nuclear space bring us closer than ever to unravelling the logic of gene regulation and its complex integration of global cellular signaling networks. Using the function of progesterone and its cellular receptor in breast cancer cells, we will briefly summarize the history and describe the present extent of our knowledge on how regulatory proteins deal with the chromatin structure to gain access to DNA sequences and interpret the genomic instructions that enable cells to respond selectively to external signals by reshaping their gene regulatory networks.

Three-Dimensional Genomic Structure and Cohesin Occupancy Correlate with Transcriptional Activity during Spermatogenesis

Mammalian gametogenesis involves dramatic and tightly regulated chromatin remodeling, whose regulatory pathways remain largely unexplored. Here, we generate a comprehensive high-resolution structural and functional atlas of mouse spermatogenesis by combining in situ chromosome conformation capture sequencing (Hi-C), RNA sequencing (RNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) of CCCTC-binding factor (CTCF) and meiotic cohesins, coupled with confocal and super-resolution microscopy. Spermatogonia presents well-defined compartment patterns and topological domains. However, chromosome occupancy and compartmentalization are highly re-arranged during prophase I, with cohesins bound to active promoters in DNA loops out of the chromosomal axes. Compartment patterns re-emerge in round spermatids, where cohesin occupancy correlates with transcriptional activity of key developmental genes. The compact sperm genome contains compartments with actively transcribed genes but no fine-scale topological domains, concomitant with the presence of protamines. Overall, we demonstrate how genome-wide cohesin occupancy and transcriptional activity is associated with three-dimensional (3D) remodeling during spermatogenesis, ultimately reprogramming the genome for the next generation.

Hormone-control regions mediate steroid receptor-dependent genome organization

In breast cancer cells, some topologically associating domains (TADs) behave as hormonal gene regulation units, within which gene transcription is coordinately regulated in response to steroid hormones. Here we further describe that responsive TADs contain 20- to 100-kb-long clusters of intermingled estrogen receptor (ESR1) and progesterone receptor (PGR) binding sites, hereafter called hormone-control regions (HCRs). In T47D cells, we identified more than 200 HCRs, which are frequently bound by unliganded ESR1 and PGR. These HCRs establish steady long-distance inter-TAD interactions between them and organize characteristic looping structures with promoters in their TADs even in the absence of hormones in ESR1⁺-PGR⁺ cells. This organization is dependent on the expression of the receptors and is further dynamically modulated in response to steroid hormones. HCRs function as platforms that integrate different signals, resulting in some cases in opposite transcriptional responses to estrogens or progestins. Altogether, these results suggest that steroid hormone receptors act not only as hormone-regulated sequence-specific transcription factors but also as local and global genome organizers.

Rapid reversible changes in compartments and local chromatin organization revealed by hyperosmotic shock

Nuclear architecture is decisive for the assembly of transcriptional responses. However, how chromosome organization is dynamically modulated to permit rapid and transient transcriptional changes in response to environmental challenges remains unclear. Here we show that hyperosmotic stress disrupts different levels of chromosome organization, ranging from A/B compartment changes to reduction in the number and insulation of topologically associating domains (TADs). Concomitantly, transcription is greatly affected, TAD borders weaken, and RNA Polymerase II runs off from hundreds of transcription end sites. Stress alters the binding profiles of architectural proteins, which explains the disappearance of local chromatin organization. These processes are dynamic, and cells rapidly reconstitute their default chromatin conformation after stress removal, uncovering an intrinsic organization. Transcription is not required for local chromatin reorganization, while compartment recovery is partially transcription-dependent. Thus, nuclear organization in mammalian cells can be rapidly modulated by environmental changes in a reversible manner.

Arginine Citrullination at the C-Terminal Domain Controls RNA Polymerase II Transcription

The post-translational modification of key residues at the C-terminal domain of RNA polymerase II (RNAP2-CTD) coordinates transcription, splicing, and RNA processing by modulating its capacity to act as a landing platform for a variety of protein complexes. Here, we identify a new modification at the CTD, the deimination of arginine and its conversion to citrulline by peptidyl arginine deiminase 2 (PADI2), an enzyme that has been associated with several diseases, including cancer. We show that, among PADI family members, only PADI2 citrullinates R1810 (Cit1810) at repeat 31 of the CTD. Depletion of PADI2 or loss of R1810 results in accumulation of RNAP2 at transcription start sites, reduced gene expression, and inhibition of cell proliferation. Cit1810 is needed for interaction with the P-TEFb (positive transcription elongation factor b) kinase complex and for its recruitment to chromatin. In this way, CTD-Cit1810 favors RNAP2 pause release and efficient transcription in breast cancer cells.

Lamin B1 mapping reveals the existence of dynamic and functional euchromatin lamin B1 domains

Lamins (A/C and B) are major constituents of the nuclear lamina (NL). Structurally conserved lamina-associated domains (LADs) are formed by genomic regions that contact the NL. Lamins are also found in the nucleoplasm, with a yet unknown function. Here we map the genome-wide localization of lamin B1 in an euchromatin-enriched fraction of the mouse genome and follow its dynamics during the epithelial-to-mesenchymal transition (EMT). Lamin B1 associates with actively expressed and open euchromatin regions, forming dynamic euchromatin lamin B1-associated domains (eLADs) of about 0.3 Mb. Hi-C data link eLADs to the 3D organization of the mouse genome during EMT and correlate lamin B1 enrichment at topologically associating domain (TAD) borders with increased border strength. Having reduced levels of lamin B1 alters the EMT transcriptional signature and compromises the acquisition of mesenchymal traits. Thus, during EMT, the process of genome reorganization in mouse involves dynamic changes in eLADs.

Parallel sequencing lives, or what makes large sequencing projects successful

T47D_rep2 and b1913e6c1_51720e9cf were 2 Hi-C samples. They were born and processed at the same time, yet their fates were very different. The life of b1913e6c1_51720e9cf was simple and fruitful, while that of T47D_rep2 was full of accidents and sorrow. At the heart of these differences lies the fact that b1913e6c1_51720e9cf was born under a lab culture of Documentation, Automation, Traceability, and Autonomy and compliance with the FAIR Principles. Their lives are a lesson for those who wish to embark on the journey of managing high-throughput sequencing data.

Distinct roles of cohesin-SA1 and cohesin-SA2 in 3D chromosome organization

Two variant cohesin complexes containing SMC1, SMC3, RAD21 and either SA1 (also known as STAG1) or SA2 (also known as STAG2) are present in all cell types. We report here their genomic distribution and specific contributions to genome organization in human cells. Although both variants are found at CCCTC-binding factor (CTCF) sites, a distinct population of the SA2-containing cohesin complexes (hereafter referred to as cohesin-SA2) localize to enhancers lacking CTCF, are linked to tissue-specific transcription and cannot be replaced by the SA1-containing cohesin complex (cohesin-SA1) when SA2 is absent, a condition that has been observed in several tumors. Downregulation of each of these variants has different consequences for gene expression and genome architecture. Our results suggest that cohesin-SA1 preferentially contributes to the stabilization of topologically associating domain boundaries together with CTCF, whereas cohesin-SA2 promotes cell-type-specific contacts between enhancers and promoters independently of CTCF. Loss of cohesin-SA2 rewires local chromatin contacts and alters gene expression. These findings provide insights into how cohesin mediates chromosome folding and establish a novel framework to address the consequences of mutations in cohesin genes in cancer.

Transcription factors orchestrate dynamic interplay between genome topology and gene regulation during cell reprogramming

Chromosomal architecture is known to influence gene expression, yet its role in controlling cell fate remains poorly understood. Reprogramming of somatic cells into pluripotent stem cells (PSCs) by the transcription factors (TFs) OCT4, SOX2, KLF4 and MYC offers an opportunity to address this question but is severely limited by the low proportion of responding cells. We have recently developed a highly efficient reprogramming protocol that synchronously converts somatic into pluripotent stem cells. Here, we used this system to integrate time-resolved changes in genome topology with gene expression, TF binding and chromatin-state dynamics. The results showed that TFs drive topological genome reorganization at multiple architectural levels, often before changes in gene expression. Removal of locus-specific topological barriers can explain why pluripotency genes are activated sequentially, instead of simultaneously, during reprogramming. Together, our results implicate genome topology as an instructive force for implementing transcriptional programs and cell fate in mammals.

Signaling by Steroid Hormones in the 3D Nuclear Space

Initial studies showed that ligand-activated hormone receptors act by binding to the proximal promoters of individual target genes. Genome-wide studies have now revealed that regulation of transcription by steroid hormones mainly depends on binding of the receptors to distal regulatory elements. Those distal elements, either enhancers or silencers, act on the regulation of target genes by chromatin looping to the gene promoters. In the nucleus, this level of chromatin folding is integrated within dynamic higher orders of genome structures, which are organized in a non-random fashion. Terminally differentiated cells exhibit a tissue-specific three-dimensional (3D) organization of the genome that favors or restrains the activity of transcription factors and modulates the function of steroid hormone receptors, which are transiently activated upon hormone exposure. Conversely, integration of the hormones signal may require modifications of the 3D organization to allow appropriate transcriptional outcomes. In this review, we summarize the main levels of organization of the genome, review how they can modulate the response to steroids in a cell specific manner and discuss the role of receptors in shaping and rewiring the structure in response to hormone. Taking into account the dynamics of 3D genome organization will contribute to a better understanding of the pleiotropic effects of steroid hormones in normal and cancer cells.

Distinct structural transitions of chromatin topological domains correlate with coordinated hormone-induced gene regulation

The human genome is segmented into topologically associating domains (TADs), but the role of this conserved organization during transient changes in gene expression is not known. Here we describe the distribution of progestin-induced chromatin modifications and changes in transcriptional activity over TADs in T47D breast cancer cells. Using ChIP-seq (chromatin immunoprecipitation combined with high-throughput sequencing), Hi-C (chromosome capture followed by high-throughput sequencing), and three-dimensional (3D) modeling techniques, we found that the borders of the ∼ 2000 TADs in these cells are largely maintained after hormone treatment and that up to 20% of the TADs could be considered as discrete regulatory units where the majority of the genes are either transcriptionally activated or repressed in a coordinated fashion. The epigenetic signatures of the TADs are homogeneously modified by hormones in correlation with the transcriptional changes. Hormone-induced changes in gene activity and chromatin remodeling are accompanied by differential structural changes for activated and repressed TADs, as reflected by specific and opposite changes in the strength of intra-TAD interactions within responsive TADs. Indeed, 3D modeling of the Hi-C data suggested that the structure of TADs was modified upon treatment. The differential responses of TADs to progestins and estrogens suggest that TADs could function as "regulons" to enable spatially proximal genes to be coordinately transcribed in response to hormones.

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.

COUP-TFI modulates estrogen signaling and influences proliferation, survival and migration of breast cancer cells

We previously showed that COUP-TFI interacts with the Estrogen Receptor alpha (ER alpha) to recruit Extracellular signal Regulated Kinases (ERKs) in an Estradiol (E2)-independent manner, resulting in an enhancement of ER alpha transcriptional activity. However, the involvement of COUP-TFI in physiologically relevant functions of ER alpha, such as the mitogenic activity that E2 has on breast cancer cells, remains poorly understood. Here, we first showed that the amounts of COUP-TFI protein are higher in dedifferentiated mammary cell lines (MDA-MB-231) and tumor breast cells as compared to the differentiated MCF-7 cell line and normal breast cells. To evaluate the functional relevance of the COUP-TFI/ER alpha interplay in mammary cells, we generated MCF-7 cells that stably over-express COUP-TFI. We found that the over-expression of COUP-TFI enhances motility and invasiveness of MCF-7 cells. COUP-TFI also promotes the proliferation of MCF-7 cells through ER alpha-dependent mechanisms that target cell cycle progression and cell survival. To further investigate the mechanisms underlying these effects of COUP-TFI, we evaluated the expression of known E2-target genes in breast cancer, and found that COUP-TFI differentially regulated genes involved in cell proliferation, apoptosis, and migration/invasion. Notably, Cathepsin D (CTSD) transcript and protein levels were significantly higher in presence and absence of E2 in MCF-7 over-expressing COUP-TFI. Chromatin Immunoprecipitation assays showed that ER alpha, phospho-RNA Polymerase II, as well as p68 RNA Helicase, a phospho-Serine 118 dependent co-activator of ER alpha, were preferentially recruited onto the CTSD gene proximal promoter in COUP-TFI over-expressing cells. These results suggest that COUP-TFI selectively regulates the expression of endogenous E2-target genes and consequently modifies ER alpha positive mammary cells response to E2.