The NSL complex-mediated nucleosome landscape is required to maintain transcription fidelity and suppression of transcription noise.

Nucleosomal organization at gene promoters is critical for transcription, with a nucleosome-depleted region (NDR) at transcription start sites (TSSs) being required for transcription initiation. How NDRs and the precise positioning of the +1 nucleosomes are maintained on active genes remains unclear. Here, we report that the Drosophila nonspecific lethal (NSL) complex is necessary to maintain this stereotypical nucleosomal organization at promoters. Upon NSL1 depletion, nucleosomes invade the NDRs at TSSs of NSL-bound genes. NSL complex member NSL3 binds to TATA-less promoters in a sequence-dependent manner. The NSL complex interacts with the NURF chromatin remodeling complex and is necessary and sufficient to recruit NURF to target promoters. Not only is the NSL complex essential for transcription, but it is required for accurate TSS selection for genes with multiple TSSs. Furthermore, loss of the NSL complex leads to an increase in transcriptional noise. Thus, the NSL complex establishes a canonical nucleosomal organization that enables transcription and determines TSS fidelity.

DYRK1B blockade promotes tumoricidal macrophage activity in pancreatic cancer.

OBJECTIVE: Highly malignant pancreatic ductal adenocarcinoma (PDAC) is characterised by an abundant immunosuppressive and fibrotic tumour microenvironment (TME). Future therapeutic attempts will therefore demand the targeting of tumours and stromal compartments in order to be effective. Here we investigate whether dual specificity and tyrosine phosphorylation-regulated kinase 1B (DYRK1B) fulfil these criteria and represent a promising anticancer target in PDAC. DESIGN: We used transplantation and autochthonous mouse models of PDAC with either genetic Dyrk1b loss or pharmacological DYRK1B inhibition, respectively. Mechanistic interactions between tumour cells and macrophages were studied in direct or indirect co-culture experiments. Histological analyses used tissue microarrays from patients with PDAC. Additional methodological approaches included bulk mRNA sequencing (transcriptomics) and proteomics (secretomics). RESULTS: We found that DYRK1B is mainly expressed by pancreatic epithelial cancer cells and modulates the influx and activity of TME-associated macrophages through effects on the cancer cells themselves as well as through the tumour secretome. Mechanistically, genetic ablation or pharmacological inhibition of DYRK1B strongly attracts tumoricidal macrophages and, in addition, downregulates the phagocytosis checkpoint and 'don't eat me' signal CD24 on cancer cells, resulting in enhanced tumour cell phagocytosis. Consequently, tumour cells lacking DYRK1B hardly expand in transplantation experiments, despite their rapid growth in culture. Furthermore, combining a small-molecule DYRK1B-directed therapy with mammalian target of rapamycin inhibition and conventional chemotherapy stalls the growth of established tumours and results in a significant extension of life span in a highly aggressive autochthonous model of PDAC. CONCLUSION: In light of DYRK inhibitors currently entering clinical phase testing, our data thus provide a novel and clinically translatable approach targeting both the cancer cell compartment and its microenvironment.

Lactate induces metabolic and epigenetic reprogramming of pro-inflammatory Th17 cells.

Increased lactate levels in the tissue microenvironment are a well-known feature of chronic inflammation. However, the role of lactate in regulating T cell function remains controversial. Here, we demonstrate that extracellular lactate predominantly induces deregulation of the Th17-specific gene expression program by modulating the metabolic and epigenetic status of Th17 cells. Following lactate treatment, Th17 cells significantly reduced their IL-17A production and upregulated Foxp3 expression through ROS-driven IL-2 secretion. Moreover, we observed increased levels of genome-wide histone H3K18 lactylation, a recently described marker for active chromatin in macrophages, in lactate-treated Th17 cells. In addition, we show that high lactate concentrations suppress Th17 pathogenicity during intestinal inflammation in mice. These results indicate that lactate is capable of reprogramming pro-inflammatory T cell phenotypes into regulatory T cells.

Comorbidity defines asthmatic patients' risk of COVID-19 hospitalization: A global perspective.

BACKGROUND: The global epidemiology of asthma among patients with coronavirus disease 2019 (COVID-19) presents striking geographic differences, defining prevalence zones of high and low co-occurrence of asthma and COVID-19. OBJECTIVE: We aimed to compare asthma prevalence among hospitalized patients with COVID-19 in major global hubs across the world by applying common inclusion criteria and definitions. METHODS: We built a network of 6 academic hospitals in Stanford (Stanford University)/the United States; Frankfurt (Goethe University), Giessen (Justus Liebig University), and Marburg (Philipps University)/Germany; and Moscow (Clinical Hospital 52 in collaboration with Sechenov University)/Russia. We collected clinical and laboratory data for patients hospitalized due to COVID-19. RESULTS: Asthmatic individuals were overrepresented among hospitalized patients with COVID-19 in Stanford and underrepresented in Moscow and Germany as compared with their prevalence among adults in the local community. Asthma prevalence was similar among patients hospitalized in an intensive care unit and patients hospitalized in other than an intensive care unit, which implied that the risk for development of severe COVID-19 was not higher among asthmatic patients. The numbers of males and comorbidities were higher among patients with COVID-19 in the Stanford cohort, and the most frequent comorbidities among these patients with asthma were other chronic inflammatory airway disorders such as chronic obstructive pulmonary disease. CONCLUSION: The observed disparity in COVID-19-associated risk among asthmatic patients across countries and continents is connected to the varying prevalence of underlying comorbidities, particularly chronic obstructive pulmonary disease.

IL-17A-producing CD8+ T cells promote PDAC via induction of inflammatory cancer-associated fibroblasts.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDAC) is characterised by an abundant desmoplastic stroma composed of cancer-associated fibroblasts (CAF) and interspersed immune cells. A non-canonical CD8+ T-cell subpopulation producing IL-17A (Tc17) promotes autoimmunity and has been identified in tumours. Here, we evaluated the Tc17 role in PDAC. DESIGN: Infiltration of Tc17 cells in PDAC tissue was correlated with patient overall survival and tumour stage. Wild-type (WT) or Il17ra-/- quiescent pancreatic stellate cells (qPSC) were exposed to conditional media obtained from Tc17 cells (Tc17-CM); moreover, co-culture of Tc17-CM-induced inflammatory (i)CAF (Tc17-iCAF) with tumour cells was performed. IL-17A/F-, IL-17RA-, RAG1-deficient and Foxn1nu/nu mice were used to study the Tc17 role in subcutaneous and orthotopic PDAC mouse models. RESULTS: Increased abundance of Tc17 cells highly correlated with reduced survival and advanced tumour stage in PDAC. Tc17-CM induced iCAF differentiation as assessed by the expression of iCAF-associated genes via synergism of IL-17A and TNF. Accordingly, IL-17RA controlled the responsiveness of qPSC to Tc17-CM. Pancreatic tumour cells co-cultured with Tc17-iCAF displayed enhanced proliferation and increased expression of genes implicated in proliferation, metabolism and protection from apoptosis. Tc17-iCAF accelerated growth of mouse and human tumours in Rag1-/- and Foxn1nu/nu mice, respectively. Finally, Il17ra-expressed by fibroblasts was required for Tc17-driven tumour growth in vivo. CONCLUSIONS: We identified Tc17 as a novel protumourigenic CD8+ T-cell subtype in PDAC, which accelerated tumour growth via IL-17RA-dependent stroma modification. We described a crosstalk between three cell types, Tc17, fibroblasts and tumour cells, promoting PDAC progression, which resulted in poor prognosis for patients.

High-Affinity Sites Form an Interaction Network to Facilitate Spreading of the MSL Complex across the X Chromosome in Drosophila.

Dosage compensation mechanisms provide a paradigm to study the contribution of chromosomal conformation toward targeting and spreading of epigenetic regulators over a specific chromosome. By using Hi-C and 4C analyses, we show that high-affinity sites (HAS), landing platforms of the male-specific lethal (MSL) complex, are enriched around topologically associating domain (TAD) boundaries on the X chromosome and harbor more long-range contacts in a sex-independent manner. Ectopically expressed roX1 and roX2 RNAs target HAS on the X chromosome in trans and, via spatial proximity, induce spreading of the MSL complex in cis, leading to increased expression of neighboring autosomal genes. We show that the MSL complex regulates nucleosome positioning at HAS, therefore acting locally rather than influencing the overall chromosomal architecture. We propose that the sex-independent, three-dimensional conformation of the X chromosome poises it for exploitation by the MSL complex, thereby facilitating spreading in males.

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq.

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites.

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

C1q/TNF-Related Protein 3 (CTRP-3) Deficiency of Adipocytes Affects White Adipose Tissue Mass but Not Systemic CTRP-3 Concentrations.

CTRP-3 (C1q/TNF-related protein-3) is an adipokine with endocrine and immunological function. The impact of adipocyte CTRP-3 production on systemic CTRP-3 concentrations and on adipocyte biology is unknown. A murine model of adipocyte CTRP-3 knockout (KO) was established (via the Cre/loxP system). Serum adipokine levels were quantified by ELISA and adipose tissue (AT) gene expression by real-time PCR. Preadipocytes were isolated from AT and differentiated into adipocytes. Comparative transcriptome analysis was applied in adipocytes and liver tissue. Body weight and AT mass were reduced in CTRP-3 KO mice together with decreased serum leptin. In primary cells from visceral AT of KO mice, expression of adiponectin, progranulin, and resistin was induced, while peroxisome proliferator activated receptor γ (PPARγ) was decreased. M1/M2 macrophage polarization markers were shifted to a more anti-inflammatory phenotype. CTRP-3 expression in AT did not contribute to serum concentrations. AT and liver morphology remained unaffected by CTRP-3 KO. Myelin transcription factor 1-like (Myt1l) was identified as a highly upregulated gene. In conclusion, adipocyte CTRP-3 has a role in adipogenesis and AT weight gain whereas adipocyte differentiation is not impaired by CTRP-3 deficiency. Since no effects on circulating CTRP-3 levels were observed, the impact of adipocyte CTRP-3 KO is limited to adipose tissue. Modified AT gene expression indicates a rather anti-inflammatory phenotype.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

Genomic dissection of enhancers uncovers principles of combinatorial regulation and cell type-specific wiring of enhancer-promoter contacts.

Genomic binding of transcription factors, like the glucocorticoid receptor (GR), is linked to the regulation of genes. However, as we show here, GR binding is a poor predictor of GR-dependent gene regulation even when taking the 3D organization of the genome into account. To connect GR binding sites to the regulation of genes in the endogenous genomic context, we turned to genome editing. By deleting GR binding sites, individually or in combination, we uncovered how cooperative interactions between binding sites contribute to the regulation of genes. Specifically, for the GR target gene GILZ, we show that the simultaneous presence of a cluster of GR binding sites is required for the activity of an individual enhancer and that the GR-dependent regulation of GILZ depends on multiple GR-bound enhancers. Further, by deleting GR binding sites that are shared between different cell types, we show how cell type-specific genome organization and enhancer-blocking can result in cell type-specific wiring of promoter-enhancer contacts. This rewiring allows an individual GR binding site shared between different cell types to direct the expression of distinct transcripts and thereby contributes to the cell type-specific consequences of glucocorticoid signaling.

coTRaCTE predicts co-occurring transcription factors within cell-type specific enhancers.

Cell-type specific gene expression is regulated by the combinatorial action of transcription factors (TFs). In this study, we predict transcription factor (TF) combinations that cooperatively bind in a cell-type specific manner. We first divide DNase hypersensitive sites into cell-type specifically open vs. ubiquitously open sites in 64 cell types to describe possible cell-type specific enhancers. Based on the pattern contrast between these two groups of sequences we develop "co-occurring TF predictor on Cell-Type specific Enhancers" (coTRaCTE) - a novel statistical method to determine regulatory TF co-occurrences. Contrasting the co-binding of TF pairs between cell-type specific and ubiquitously open chromatin guarantees the high cell-type specificity of the predictions. coTRaCTE predicts more than 2000 co-occurring TF pairs in 64 cell types. The large majority (70%) of these TF pairs is highly cell-type specific and overlaps in TF pair co-occurrence are highly consistent among related cell types. Furthermore, independently validated co-occurring and directly interacting TFs are significantly enriched in our predictions. Focusing on the regulatory network derived from the predicted co-occurring TF pairs in embryonic stem cells (ESCs) we find that it consists of three subnetworks with distinct functions: maintenance of pluripotency governed by OCT4, SOX2 and NANOG, regulation of early development governed by KLF4, STAT3, ZIC3 and ZNF148 and general functions governed by MYC, TCF3 and YY1. In summary, coTRaCTE predicts highly cell-type specific co-occurring TFs which reveal new insights into transcriptional regulatory mechanisms.

ChIP-exo signal associated with DNA-binding motifs provides insight into the genomic binding of the glucocorticoid receptor and cooperating transcription factors.

The classical DNA recognition sequence of the glucocorticoid receptor (GR) appears to be present at only a fraction of bound genomic regions. To identify sequences responsible for recruitment of this transcription factor (TF) to individual loci, we turned to the high-resolution ChIP-exo approach. We exploited this signal by determining footprint profiles of TF binding at single-base-pair resolution using ExoProfiler, a computational pipeline based on DNA binding motifs. When applied to our GR and the few available public ChIP-exo data sets, we find that ChIP-exo footprints are protein- and recognition sequence-specific signatures of genomic TF association. Furthermore, we show that ChIP-exo captures information about TFs other than the one directly targeted by the antibody in the ChIP procedure. Consequently, the shape of the ChIP-exo footprint can be used to discriminate between direct and indirect (tethering to other DNA-bound proteins) DNA association of GR. Together, our findings indicate that the absence of classical recognition sequences can be explained by direct GR binding to a broader spectrum of sequences than previously known, either as a homodimer or as a heterodimer binding together with a member of the ETS or TEAD families of TFs, or alternatively by indirect recruitment via FOX or STAT proteins. ChIP-exo footprints also bring structural insights and locate DNA:protein cross-link points that are compatible with crystal structures of the studied TFs. Overall, our generically applicable footprint-based approach uncovers new structural and functional insights into the diverse ways of genomic cooperation and association of TFs.

Inference of interactions between chromatin modifiers and histone modifications: from ChIP-Seq data to chromatin-signaling.

Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro. Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

A link between Sas2-mediated H4 K16 acetylation, chromatin assembly in S-phase by CAF-I and Asf1, and nucleosome assembly by Spt6 during transcription.

The histone acetyltransferase Sas2 is part of the SAS-I complex and acetylates lysine 16 of histone H4 (H4 K16Ac) in the genome of Saccharomyces cerevisiae. Sas2-mediated H4 K16Ac is strongest over the coding region of genes with low expression. However, it is unclear how Sas2-mediated acetylation is incorporated into chromatin. Our previous work has shown physical interactions of SAS-I with the histone chaperones CAF-I and Asf1, suggesting a link between SAS-I-mediated acetylation and chromatin assembly. Here, we find that Sas2-dependent H4 K16Ac in bulk histones requires passage of the cells through the S-phase of the cell cycle, and the rate of increase in H4 K16Ac depends on both CAF-I and Asf1, whereas steady-state levels and genome-wide distribution of H4 K16Ac show only mild changes in their absence. Furthermore, H4 K16Ac is deposited in chromatin at genes upon repression, and this deposition requires the histone chaperone Spt6, but not CAF-I, Asf1, HIR or Rtt106. Altogether, our data indicate that Spt6 controls H4 K16Ac levels by incorporating K16-unacetylated H4 in strongly transcribed genes. Upon repression, Spt6 association is decreased, resulting in less deposition of K16-unacetylated H4 and therefore in a concomitant increase of H4 K16Ac that is recycled during transcription.

Inferring nucleosome positions with their histone mark annotation from ChIP data.

MOTIVATION: The nucleosome is the basic repeating unit of chromatin. It contains two copies each of the four core histones H2A, H2B, H3 and H4 and about 147 bp of DNA. The residues of the histone proteins are subject to numerous post-translational modifications, such as methylation or acetylation. Chromatin immunoprecipitiation followed by sequencing (ChIP-seq) is a technique that provides genome-wide occupancy data of these modified histone proteins, and it requires appropriate computational methods. RESULTS: We present NucHunter, an algorithm that uses the data from ChIP-seq experiments directed against many histone modifications to infer positioned nucleosomes. NucHunter annotates each of these nucleosomes with the intensities of the histone modifications. We demonstrate that these annotations can be used to infer nucleosomal states with distinct correlations to underlying genomic features and chromatin-related processes, such as transcriptional start sites, enhancers, elongation by RNA polymerase II and chromatin-mediated repression. Thus, NucHunter is a versatile tool that can be used to predict positioned nucleosomes from a panel of histone modification ChIP-seq experiments and infer distinct histone modification patterns associated to different chromatin states. AVAILABILITY: The software is available at http://epigen.molgen.mpg.de/nuchunter/.

Finding associations among histone modifications using sparse partial correlation networks.

Histone modifications are known to play an important role in the regulation of transcription. While individual modifications have received much attention in genome-wide analyses, little is known about their relationships. Some authors have built Bayesian networks of modifications, however most often they have used discretized data, and relied on unrealistic assumptions such as the absence of feedback mechanisms or hidden confounding factors. Here, we propose to infer undirected networks based on partial correlations between histone modifications. Within the partial correlation framework, correlations among two variables are controlled for associations induced by the other variables. Partial correlation networks thus focus on direct associations of histone modifications. We apply this methodology to data in CD4+ cells. The resulting network is well supported by common knowledge. When pairs of modifications show a large difference between their correlation and their partial correlation, a potential confounding factor is identified and provided as explanation. Data from different cell types (IMR90, H1) is also exploited in the analysis to assess the stability of the networks. The results are remarkably similar across cell types. Based on this observation, the networks from the three cell types are integrated into a consensus network to increase robustness. The data and the results discussed in the manuscript can be found, together with code, on http://spcn.molgen.mpg.de/index.html.

Genome-wide H4 K16 acetylation by SAS-I is deposited independently of transcription and histone exchange.

The MYST HAT Sas2 is part of the SAS-I complex that acetylates histone H4 lysine 16 (H4 K16Ac) and blocks the propagation of heterochromatin at the telomeres of Saccharomyces cerevisiae. In this study, we investigated Sas2-mediated H4 K16Ac on a genome-wide scale. Interestingly, H4 K16Ac loss in sas2Δ cells outside of the telomeric regions showed a distinctive pattern in that there was a pronounced decrease of H4 K16Ac within the majority of open reading frames (ORFs), but little change in intergenic regions. Furthermore, regions of low histone H3 exchange and low H3 K56 acetylation showed the most pronounced loss of H4 K16Ac in sas2Δ, indicating that Sas2 deposited this modification on chromatin independently of histone exchange. In agreement with the effect of Sas2 within ORFs, sas2Δ caused resistance to 6-azauracil, indicating a positive effect on transcription elongation in the absence of H4 K16Ac. In summary, our data suggest that Sas2-dependent H4 K16Ac is deposited into chromatin independently of transcription and histone exchange, and that it has an inhibitory effect on the ability of PolII to travel through the body of the gene.