Cactus: A user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis.

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner.

Histone Purification Combined with High-Resolution Mass Spectrometry to Examine Histone Post-Translational Modifications and Histone Variants in Caenorhabditis elegans.

Histones are the major proteinaceous component of chromatin in eukaryotic cells and an important part of the epigenome, affecting most DNA-related events, including transcription, DNA replication, and chromosome segregation. The properties of histones are greatly influenced by their post-translational modifications (PTMs), over 200 of which are known today. Given this large number, researchers need sophisticated methods to study histone PTMs comprehensively. In particular, mass spectrometry (MS)-based approaches have gained popularity, allowing for the quantification of dozens of histone PTMs at once. Using these approaches, even the study of co-occurring PTMs and the discovery of novel PTMs become feasible. The success of MS-based approaches relies substantially on obtaining pure and well-preserved histones for analysis, which can be difficult depending on the source material. Caenorhabditis elegans has been a popular model organism to study the epigenome, but isolation of pure histones from these animals has been challenging. Here, we address this issue, presenting a method for efficient isolation of pure histone proteins from C. elegans at good yield. Further, we describe an MS pipeline optimized for accurate relative quantification of histone PTMs from C. elegans. We alkylate and tryptically digest the histones, analyze them by bottom-up MS, and then evaluate the resulting data by a C. elegans-adapted version of the software EpiProfile 2.0. Finally, we show the utility of this pipeline by determining differences in histone PTMs between C. elegans strains that age at different rates and thereby achieve very different lifespans. © 2020 The Authors. Basic Protocol 1: Large-scale growth and harvesting of synchronized C. elegans Basic Protocol 2: Nuclear preparation, histone extraction, and histone purification Basic Protocol 3: Bottom-up mass spectrometry analysis of histone PTMs and histone variants.

Transcriptomics-Based Screening Identifies Pharmacological Inhibition of Hsp90 as a Means to Defer Aging.

Aging strongly influences human morbidity and mortality. Thus, aging-preventive compounds could greatly improve our health and lifespan. Here we screened for such compounds, known as geroprotectors, employing the power of transcriptomics to predict biological age. Using age-stratified human tissue transcriptomes and machine learning, we generated age classifiers and applied these to transcriptomic changes induced by 1,309 different compounds in human cells, ranking these compounds by their ability to induce a "youthful" transcriptional state. Testing the top candidates in C. elegans, we identified two Hsp90 inhibitors, monorden and tanespimycin, which extended the animals' lifespan and improved their health. Hsp90 inhibition induces expression of heat shock proteins known to improve protein homeostasis. Consistently, monorden treatment improved the survival of C. elegans under proteotoxic stress, and its benefits depended on the cytosolic unfolded protein response-inducing transcription factor HSF-1. Taken together, our method represents an innovative geroprotector screening approach and was able to identify a class that acts by improving protein homeostasis.

Histone H1 depletion triggers an interferon response in cancer cells via activation of heterochromatic repeats.

Histone H1 has seven variants in human somatic cells and contributes to chromatin compaction and transcriptional regulation. Knock-down (KD) of each H1 variant in breast cancer cells results in altered gene expression and proliferation differently in a variant specific manner with H1.2 and H1.4 KDs being most deleterious. Here we show combined depletion of H1.2 and H1.4 has a strong deleterious effect resulting in a strong interferon (IFN) response, as evidenced by an up-regulation of many IFN-stimulated genes (ISGs) not seen in individual nor in other combinations of H1 variant KDs. Although H1 participates to repress ISG promoters, IFN activation upon H1.2 and H1.4 KD is mainly generated through the activation of the IFN response by cytosolic nucleic acid receptors and IFN synthesis, and without changes in histone modifications at induced ISG promoters. H1.2 and H1.4 co-KD also promotes the appearance of accessibility sites genome wide and, particularly, at satellites and other repeats. The IFN response may be triggered by the expression of noncoding RNA generated from heterochromatic repeats or endogenous retroviruses upon H1 KD. In conclusion, redundant H1-mediated silencing of heterochromatin is important to maintain cell homeostasis and to avoid an unspecific IFN response.

Specificities and genomic distribution of somatic mammalian histone H1 subtypes.

Histone H1 is a structural component of chromatin that may have a role in the regulation of chromatin dynamics. Unlike core histones, the linker histone H1 family is evolutionarily diverse and many organisms have multiple H1 variants or subtypes, distinguishable between germ-line and somatic cells. In mammals, the H1 family includes seven somatic H1 variants with a prevalence that varies between cell types and over the course of differentiation, H1.1 to H1.5 being expressed in a replication-dependent manner, whereas H1.0 and H1X are replication-independent. Until recently, it has not been known whether the different variants had specific roles in the regulation of nuclear processes or were differentially distributed across the genome. To address this, an increasing effort has been made to investigate divergent features among H1 variants, regarding their structure, expression patterns, chromatin dynamics, post-translational modifications and genome-wide distribution. Although H1 subtypes seem to have redundant functions, several reports point to the idea that they are also differently involved in specific cellular processes. Initial studies investigating the genomic distribution of H1 variants have started to suggest that despite a wide overlap, different variants may be enriched or preferentially located at different chromatin types, but this may depend on the cell type, the relative abundance of the variants, the differentiation state of the cell, or whether cells are derived from a neoplastic process. Understanding the heterogeneity of the histone H1 family is crucial to elucidate their role in chromatin organization, gene expression regulation and other cellular processes.

PARP1- and CTCF-Mediated Interactions between Active and Repressed Chromatin at the Lamina Promote Oscillating Transcription.

Transcriptionally active and inactive chromatin domains tend to segregate into separate sub-nuclear compartments to maintain stable expression patterns. However, here we uncovered an inter-chromosomal network connecting active loci enriched in circadian genes to repressed lamina-associated domains (LADs). The interactome is regulated by PARP1 and its co-factor CTCF. They not only mediate chromatin fiber interactions but also promote the recruitment of circadian genes to the lamina. Synchronization of the circadian rhythm by serum shock induces oscillations in PARP1-CTCF interactions, which is accompanied by oscillating recruitment of circadian loci to the lamina, followed by the acquisition of repressive H3K9me2 marks and transcriptional attenuation. Furthermore, depletion of H3K9me2/3, inhibition of PARP activity by olaparib, or downregulation of PARP1 or CTCF expression counteracts both recruitment to the envelope and circadian transcription. PARP1- and CTCF-regulated contacts between circadian loci and the repressive chromatin environment at the lamina therefore mediate circadian transcriptional plasticity.

Genome distribution of replication-independent histone H1 variants shows H1.0 associated with nucleolar domains and H1X associated with RNA polymerase II-enriched regions.

Unlike core histones, the linker histone H1 family is more evolutionarily diverse, and many organisms have multiple H1 variants or subtypes. In mammals, the H1 family includes seven somatic H1 variants; H1.1 to H1.5 are expressed in a replication-dependent manner, whereas H1.0 and H1X are replication-independent. Using ChIP-sequencing data and cell fractionation, we have compared the genomic distribution of H1.0 and H1X in human breast cancer cells, in which we previously observed differential distribution of H1.2 compared with the other subtypes. We have found H1.0 to be enriched at nucleolus-associated DNA repeats and chromatin domains, whereas H1X is associated with coding regions, RNA polymerase II-enriched regions, and hypomethylated CpG islands. Further, H1X accumulates within constitutive or included exons and retained introns and toward the 3' end of expressed genes. Inducible H1X knockdown does not affect cell proliferation but dysregulates a subset of genes related to cell movement and transport. In H1X-depleted cells, the promoters of up-regulated genes are not occupied specifically by this variant, have a lower than average H1 content, and, unexpectedly, do not form an H1 valley upon induction. We conclude that H1 variants are not distributed evenly across the genome and may participate with some specificity in chromatin domain organization or gene regulation.

Dynamics and dispensability of variant-specific histone H1 Lys-26/Ser-27 and Thr-165 post-translational modifications.

In mammals, the linker histone H1, involved in DNA packaging into chromatin, is represented by a family of variants. H1 tails undergo post-translational modifications (PTMs) that can be detected by mass spectrometry. We developed antibodies to analyze several of these as yet unexplored PTMs including the combination of H1.4 K26 acetylation or trimethylation and S27 phosphorylation. H1.2-T165 phosphorylation was detected at S and G2/M phases of the cell cycle and was dispensable for chromatin binding and cell proliferation; while the H1.4-K26 residue was essential for proper cell cycle progression. We conclude that histone H1 PTMs are dynamic over the cell cycle and that the recognition of modified lysines may be affected by phosphorylation of adjacent residues.

Mapping of six somatic linker histone H1 variants in human breast cancer cells uncovers specific features of H1.2.

Seven linker histone H1 variants are present in human somatic cells with distinct prevalence across cell types. Despite being key structural components of chromatin, it is not known whether the different variants have specific roles in the regulation of nuclear processes or are differentially distributed throughout the genome. Using variant-specific antibodies to H1 and hemagglutinin (HA)-tagged recombinant H1 variants expressed in breast cancer cells, we have investigated the distribution of six H1 variants in promoters and genome-wide. H1 is depleted at promoters depending on its transcriptional status and differs between variants. Notably, H1.2 is less abundant than other variants at the transcription start sites of inactive genes, and promoters enriched in H1.2 are different from those enriched in other variants and tend to be repressed. Additionally, H1.2 is enriched at chromosomal domains characterized by low guanine-cytosine (GC) content and is associated with lamina-associated domains. Meanwhile, other variants are associated with higher GC content, CpG islands and gene-rich domains. For instance, H1.0 and H1X are enriched at gene-rich chromosomes, whereas H1.2 is depleted. In short, histone H1 is not uniformly distributed along the genome and there are differences between variants, H1.2 being the one showing the most specific pattern and strongest correlation with low gene expression.

Progestins activate 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) in breast cancer cells.

PFKFB (6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase) catalyses the synthesis and degradation of Fru-2,6-P2 (fructose-2,6-bisphosphate), a key modulator of glycolysis and gluconeogenesis. The PFKFB3 gene is extensively involved in cell proliferation owing to its key role in carbohydrate metabolism. In the present study we analyse its mechanism of regulation by progestins in breast cancer cells. We report that exposure of T47D cells to synthetic progestins (ORG2058 or norgestrel) leads to a rapid increase in Fru-2,6-P2 concentration. Our Western blot results are compatible with a short-term activation due to PFKFB3 isoenzyme phosphorylation and a long-term sustained action due to increased PFKFB3 protein levels. Transient transfection of T47D cells with deleted gene promoter constructs allowed us to identify a PRE (progesterone-response element) to which PR (progesterone receptor) binds and thus transactivates PFKFB3 gene transcription. PR expression in the PR-negative cell line MDA-MB-231 induces endogenous PFKFB3 expression in response to norgestrel. Direct binding of PR to the PRE box (-3490 nt) was confirmed by ChIP (chromatin immunoprecipiation) experiments. A dual mechanism affecting PFKFB3 protein and gene regulation operates in order to assure glycolysis in breast cancer cells. An immediate early response through the ERK (extracellular-signal-regulated kinase)/RSK (ribosomal S6 kinase) pathway leading to phosphorylation of PFKFB3 on Ser461 is followed by activation of mRNA transcription via cis-acting sequences on the PFKFB3 promoter.

Histone H1 variants are differentially expressed and incorporated into chromatin during differentiation and reprogramming to pluripotency.

There are seven linker histone variants in human somatic cells (H1.0 to H1.5 and H1X), and their prevalence varies as a function of cell type and differentiation stage, suggesting that the different variants may have distinct roles. We have revisited this notion by using new methodologies to study pluripotency and differentiation, including the in vitro differentiation of human embryonic stem (ES) and teratocarcinoma cells and the reprogramming of keratinocytes to induced pluripotent stem cells. Our results show that pluripotent cells (PCs) have decreased levels of H1.0 and increased levels of H1.1, H1.3, and H1.5 compared with differentiated cells. PCs have a more diverse repertoire of H1 variants, whereas in differentiated cells, H1.0 expression represents ∼80% of the H1 transcripts. In agreement with their prevalent expression in ES cells, the regulatory regions of H1.3 and H1.5 genes were found to be occupied by pluripotency factors. Moreover, the H1.0 gene promoter contains bivalent domains (H3K4me2 and H3K27me3) in PCs, suggesting that this variant is likely to have an important role during differentiation. Indeed, the knockdown of H1.0 in human ES did not affect self-renewal but impaired differentiation. Accordingly, H1.0 was recruited to the regulatory regions of differentiation and pluripotency genes during differentiation, confirming that this histone variant plays a critical role in the regulation of these genes. Thus, histone H1 variant expression is controlled by a variety of mechanisms that produce distinct but consistent H1 repertoires in pluripotent and differentiated cells that appear critical to maintain the functionality of such cells.

Mutational analysis of progesterone receptor functional domains in stable cell lines delineates sets of genes regulated by different mechanisms.

Steroid hormone receptors act directly in the nucleus on the chromatin organization and transcriptional activity of several promoters. Furthermore, they have an indirect effect on cytoplasmic signal transduction pathways, including MAPK, impacting ultimately on gene expression. We are interested in distinguishing between the two modes of action of progesterone receptor (PR) on the control of gene expression and cell proliferation. For this, we have stably expressed, in PR-negative breast cancer cells, tagged forms of the PR isoform B mutated at regions involved either in DNA binding (DNA-binding domain) or in its ability to interact with the estrogen receptor and to activate the c-Src/MAPK/Erk/Msk cascade (estrogen receptor-interacting domain). Both mutants impair PR-mediated activation of a well-understood model promoter in response to progestin, as well as hormone-induced cell proliferation. Additional mutants affecting transactivation activity of PR (activation function 2) or a zinc-finger implicated in dimerization (D-box) have also been tested. Microarrays and gene expression experiments on these cell lines define the subsets of hormone-responsive genes regulated by different modes of action of PR isoform B, as well as genes in which the nuclear and nongenomic pathways cooperate. Correlation between CCND1 expression in the different cell lines and their ability to support cell proliferation confirms CCND1 as a key controller gene.

Progesterone induction of the 11beta-hydroxysteroid dehydrogenase type 2 promoter in breast cancer cells involves coordinated recruitment of STAT5A and progesterone receptor to a distal enhancer and polymerase tracking.

Steroid hormone receptors regulate gene expression, interacting with target DNA sequences but also activating cytoplasmic signaling pathways. Using the human 11beta-hydroxysteroid dehydrogenase type 2 (11beta-HSD2) gene as a model, we have investigated the contributions of both effects on a human progesterone-responsive promoter in breast cancer cells. Chromatin immunoprecipitation has identified two different mechanisms of hormone-induced progesterone receptor (PR) recruitment to the 11beta-HSD2 promoter: (i) direct PR binding to DNA at the proximal promoter, abrogated when PR contains a mutated DNA binding domain (DBD), and (ii) STAT5A (signal transducer and activator of transcription 5A)-mediated recruitment of PR to an upstream distal region, impaired by AG490, a JAK/STAT pathway inhibitor. The JAK/STAT inhibitor, as well as expression of dominant-negative STAT5A, impairs hormone induction of 11beta-HSD2. On the other hand, the DBD-mutated PR fully supports 11beta-HSD2 expression. These results, along with data from a deletion analysis, indicate that the distal region is crucial for hormone regulation of 11beta-HSD2. We show active RNA polymerase II tracking from the distal region upon PR and STAT5A binding, concomitant with synthesis of noncoding, hormone-dependent RNAs, suggesting that this region works as a hormone-dependent transcriptional enhancer. In conclusion, coordination of PR transcriptional effects and cytoplasmic signaling activation, in particular the JAK/STAT pathway, are critical in regulating progestin-induced endogenous 11beta-HSD2 gene expression in breast cancer cells. This is not unique to this promoter, as AG490 also alters the expression of other progesterone-regulated genes.