Abstract Neural activity-dependent gene regulation is central to the development of neural networks and neuronal plasticity. Induction of activity-dependent gene programs is equally important as repression of these programs, and both need to be balanced carefully. However, little is known about how repressive mechanisms modulate neuronal responsiveness across the genome before stimulation. Here, we identify nucleoporin-dependent regulation of neuronal responsiveness, in which Nup153 represses neuronal genes including activity-regulated genes (ARGs), in the basal state. By characterizing the genome-wide landscape of chromatin accessibility, histone modifications and Nup153 chromatin binding, we show that Nup153 bidirectionally regulates chromatin states through both basal activity-dependent and -independent mechanisms and influences associated genes. Mechanistically, Nup153 associates with HDAC1 to modulate histone acetylation and chromatin accessibility at target regulatory regions. Our data suggests Nup153 organizes chromatin states that regulate neuronal gene programs involved in maintaining neuronal responsiveness. Teaser Nup153 primes neuronal responsiveness through multi-layered epigenetic regulation.

Epigenetic mechanisms regulate gene expression programs during neurogenesis, but the extent of epigenetic remodeling during human cortical development remains unknown. Here, we characterize the epigenetic landscape of the human developing neocortex by leveraging Epi-CyTOF, a mass-cytometry-based approach for the simultaneous single-cell analysis of more than 30 epigenetic marks. We identify Polycomb repressive complex 2 (PRC2)-mediated H3K27me3 as the modification with the strongest cell-type-specific enrichment. Inhibition of PRC2 in human cortical organoids resulted in a shift of neural progenitor cell (NPC) proliferation toward differentiation. Cell-type-specific profiling of H3K27me3 identified neuronal differentiation and extracellular matrix (ECM) genes in the human neocortex. PRC2 inhibition resulted in increased production of the ECM proteins Syndecan 1 and laminin alpha 1. Overall, this study comprehensively characterizes the epigenetic state of specific neural cell types and highlights a novel role for H3K27me3 in regulating the ECM composition in the human developing neocortex.

While pancreatic beta-cell proliferation has been extensively studied, the role of cell death during islet development remains incompletely understood. Using a genetic model of caspase inhibition in beta cells coupled with mathematical modeling, we here discover an onset of beta-cell death in juvenile zebrafish, which regulates beta-cell mass. Histologically, this beta-cell death is underestimated due to phagocytosis by resident macrophages. To investigate beta-cell apoptosis at the molecular level, we implement a conditional model of beta-cell death linked to Ca2+ overload. Transcriptomic analysis reveals that metabolically-stressed beta cells follow paths to either de-differentiation or apoptosis. Beta cells destined to die activate inflammatory and immuno-regulatory pathways, suggesting that cell death regulates the crosstalk with immune cells. Consistently, inhibiting beta-cell death during development reduces pro-inflammatory resident macrophages and expands T-regulatory cells, the deficiency of which causes premature activation of NF-kB signaling in beta cells. Thus, developmental cell death not only shapes beta-cell mass but it also influences the islet's inflammatory milieu by shifting the immune-cell population towards pro-inflammatory.

A spatial gene expression pattern between the cortex (CTX) and inner medulla (IM) of the kidney has been observed, but the underlying mechanisms are unclear. Understanding these mechanisms is essential for elucidating renal function. Using the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) we analyzed the open chromatin structures and the involvement of epigenetic mechanisms in mediating gene expression differences between the renal CTX and IM. We also examined the role of the nuclear factor of activated T cells 5 (NFAT5), a key regulator of hypertonicity. ATAC-seq analysis was performed on CTX and IM samples from both wild-type (WT) and NFAT5 knockout (KO) mice.This work demonstrates for the first time that these differences in gene expression between renal CTX and IM are associated with an epigenetic mechanism driven by chromatin accessibility, which is partially modulated by the nuclear factor of activated T-cells 5 (NFAT5) in mice. Furthermore, spatial localization and NFAT5-promoted chromatin accessibility correlate with differential gene expression and altered promoter binding motif enrichment in CTX and IM.This study provides new insights into the spatial and NFAT5-mediated regulation of chromatin accessibility and gene expression in CTX and IM. This work advances our understanding of kidney physiology by uncovering previously unknown epigenetic factors influencing gene expression and provides a new perspective on renal adaptive mechanisms.TRANSLATIONAL STATEMENT The study reveals new insights into the spatial and epigenetic regulation of gene expression in the renal cortex (CTX) and inner medulla (IM) in the mouse kidney. We used the Assay for Transposase-Accessible Chromatin with High-Throughput Sequence Analysis (ATAC-seq) to identify a key role of NFAT5 in modulating chromatin accessibility and to uncover previously unknown epigenetic factors. This research enhances our understanding of renal physiology and has important implications for clinical care by providing insights into potential adaptive mechanisms in the kidney. These findings suggest future investigations targeting epigenetic signaling pathways for therapeutic intervention in renal diseases.### Competing Interest StatementThe authors have declared no competing interest.

This preclinical study in mice shows that the β3-adrenergic receptor can be a potential therapeutic approach to counteracting glucocorticoid (GC)-induced obesity and metabolic dysfunction. Both cold acclimation and β3-adrenergic receptor stimulation in a mouse model of excess glucocorticoids were adequate in not only preventing obesity, adiposity, and adipose tissue dysfunction but also correcting hyperinsulinemia, hyperleptinemia, and dyslipidemia.

[This corrects the article DOI: 10.3389/fimmu.2023.1157373.].

Allergic inflammation of the airways such as allergic asthma is a major health problem with growing incidence world-wide. One cardinal feature in severe type 2-dominated airway inflammation is the release of lipid mediators of the eicosanoid family that can either promote or dampen allergic inflammation. Macrophages are key producers of prostaglandins and leukotrienes which play diverse roles in allergic airway inflammation and thus require tight control. Using RNA- and ATAC-sequencing, liquid chromatography coupled to mass spectrometry (LC-MS/MS), enzyme immunoassays (EIA), gene expression analysis and in vivo models, we show that the aryl hydrocarbon receptor (AhR) contributes to this control via transcriptional regulation of lipid mediator synthesis enzymes in bone marrow-derived as well as in primary alveolar macrophages. In the absence or inhibition of AhR activity, multiple genes of both the prostaglandin and the leukotriene pathway were downregulated, resulting in lower synthesis of prostanoids, such as prostaglandin E2 (PGE2), and cysteinyl leukotrienes, e.g., Leukotriene C4 (LTC4). These AhR-dependent genes include PTGS1 encoding for the enzyme cyclooxygenase 1 (COX1) and ALOX5 encoding for the arachidonate 5-lipoxygenase (5-LO) both of which major upstream regulators of the prostanoid and leukotriene pathway, respectively. This regulation is independent of the activation stimulus and partially also detectable in unstimulated macrophages suggesting an important role of basal AhR activity for eicosanoid production in steady state macrophages. Lastly, we demonstrate that AhR deficiency in hematopoietic but not epithelial cells aggravates house dust mite induced allergic airway inflammation. These results suggest an essential role for AhR-dependent eicosanoid regulation in macrophages during homeostasis and inflammation.

Cohesin SMC1β; ensures silencing at chromosome ends, restricts TERRA expression and RNA-DNA hybrid formation, and prevents telomer damage. Previous data showed that meiotic cohesin SMC1β protects spermatocyte telomeres from damage. The underlying reason, however, remained unknown as the expressions of telomerase and shelterin components were normal in Smc1β−/− spermatocytes. Here. we report that SMC1β restricts expression of the long noncoding RNA TERRA (telomeric repeat containing RNA) in spermatocytes. In somatic cell lines increased TERRA was reported to cause telomere damage through altering telomere chromatin structure. In Smc1β−/− spermatocytes, we observed strongly increased levels of TERRA which accumulate on damaged chromosomal ends, where enhanced R-loop formation was found. This suggested a more open chromatin configuration near telomeres in Smc1β−/− spermatocytes, which was confirmed by ATAC-seq. Telomere-distal regions were not affected by the absence of SMC1β but RNA-seq revealed increased transcriptional activity in telomere-proximal regions. Thus, SMC1β promotes closed chromatin specifically near telomeres and limits TERRA expression in spermatocytes.

Loss of limbs evolved many times in squamate reptiles. Here we investigated the genomic basis of convergent limb loss in reptiles. We sequenced the genomes of a closely related pair of limbless-limbed gymnophthalmid lizards and performed a comparative genomic analysis including five snakes and the limbless glass lizard. Our analysis of these three independent limbless lineages revealed that signatures of shared sequence or transcription factor binding site divergence in individual limb regulatory elements are generally rare. Instead, shared divergence occurs more often at the level of signaling pathways, involving different regulatory elements associated with the same limb genes (such as Hand2 or Hox) and/or patterning mechanisms (such as Shh signaling). Interestingly, although snakes are known to have mutations in the Shh ZRS limb enhancer, this enhancer lacks relevant mutations in limbless lizards. Thus, different mechanisms could contribute to limb loss, and there are likely multiple evolutionary paths to limblessness in reptiles.

Hematopoietic stem cells (HSCs) produce highly diverse cell lineages. Here, we chart native lineage pathways emanating from HSCs and define their physiological regulation by computationally integrating experimental approaches for fate mapping, mitotic tracking, and single-cell RNA sequencing. We find that lineages begin to split when cells leave the tip HSC population, marked by high Sca-1 and CD201 expression. Downstream, HSCs either retain high Sca-1 expression and the ability to generate lymphocytes, or irreversibly reduce Sca-1 level and enter into erythro-myelopoiesis or thrombopoiesis. Thrombopoiesis is the sum of two pathways that make comparable contributions in steady state, a long route via multipotent progenitors and CD48 hi megakaryocyte progenitors (MkPs), and a short route from HSCs to developmentally distinct CD48 −/lo MkPs. Enhanced thrombopoietin signaling differentially accelerates the short pathway, enabling a rapid response to increasing demand. In sum, we provide a blueprint for mapping physiological differentiation fluxes from HSCs and decipher two functionally distinct pathways of native thrombopoiesis.

Hematopoietic stem cells (HSCs) are the ultimate source of a large variety of mature blood and immune cells. While the potential of HSCs has been well studied by perturbative transplantation assays, their properties and differentiation pathways remain poorly defined in physiological hematopoiesis. Here, we chart native lineage pathways emanating from HSCs and define their physiological regulation by integrating non-perturbative experimental approaches for fate mapping, mitotic tracking and single-cell RNA sequencing. We find that lineages begin to split when cells leave the tip HSC population, marked by high Sca-1 and CD201 surface expression. Downstream, HSCs either retain high Sca-1 expression and the ability to generate lymphocytes, or irreversibly reduce Sca-1 level and enter into erythro-myelopoiesis or thrombopoiesis. Thrombopoiesis is the sum of two pathways that make comparable contributions in steady state, a long route via multipotent progenitors and CD48hi megakaryocyte progenitors (MkPs), and a short route from HSCs to developmentally distinct CD48-/lo MkPs. This newly defined progenitor represents the missing link between HSCs and megakaryocytes. Enhanced thrombopoietin signaling differentially accelerates this short pathway, enabling a rapid response to increasing demand. In sum, we provide a blueprint for mapping physiological differentiation fluxes from HSCs and decipher two functionally distinct pathways of native thrombopoiesis. Hematopoietic stem cells (HSCs) are the ultimate source of a large variety of mature blood and immune cells. While the potential of HSCs has been well studied by perturbative transplantation assays, their properties and differentiation pathways remain poorly defined in physiological hematopoiesis. Here, we chart native lineage pathways emanating from HSCs and define their physiological regulation by integrating non-perturbative experimental approaches for fate mapping, mitotic tracking and single-cell RNA sequencing. We find that lineages begin to split when cells leave the tip HSC population, marked by high Sca-1 and CD201 surface expression. Downstream, HSCs either retain high Sca-1 expression and the ability to generate lymphocytes, or irreversibly reduce Sca-1 level and enter into erythro-myelopoiesis or thrombopoiesis. Thrombopoiesis is the sum of two pathways that make comparable contributions in steady state, a long route via multipotent progenitors and CD48hi megakaryocyte progenitors (MkPs), and a short route from HSCs to developmentally distinct CD48-/lo MkPs. This newly defined progenitor represents the missing link between HSCs and megakaryocytes. Enhanced thrombopoietin signaling differentially accelerates this short pathway, enabling a rapid response to increasing demand. In sum, we provide a blueprint for mapping physiological differentiation fluxes from HSCs and decipher two functionally distinct pathways of native thrombopoiesis.

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cisregulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Background:Reptiles are a species-rich group with great phenotypic and life history diversity but are highly underrepresented among the vertebrate species with sequenced genomes.Results:Here, we report a high-quality genome assembly of the tegu lizard, Salvator merianae, the first lacertoid with a sequenced genome. We combined 74X Illumina short-read, 29.8X Pacific Biosciences long-read, and optical mapping data to generate a high-quality assembly with a scaffold N50 value of 55.4 Mb. The contig N50 value of this assembly is 521 Kb, making it the most contiguous reptile assembly so far. We show that the tegu assembly has the highest completeness of coding genes and conserved non-exonic elements (CNEs) compared to other reptiles. Furthermore, the tegu assembly has the highest number of evolutionarily conserved CNE pairs, corroborating a high assembly contiguity in intergenic regions. As in other reptiles, long interspersed nuclear elements comprise the most abundant transposon class. We used transcriptomic data, homology- and de novo gene predictions to annotate 22,413 coding genes, of which 16,995 (76%) likely have human orthologs as inferred by CESAR-derived gene mappings. Finally, we generated a multiple genome alignment comprising 10 squamates and 7 other amniote species and identified conserved regions that are under evolutionary constraint. CNEs cover 38 Mb (1.8%) of the tegu genome, with 3.3 Mb in these elements being squamate specific. In contrast to placental mammal-specific CNEs, very few of these squamate-specific CNEs (<20 Kb) overlap transposons, highlighting a difference in how lineage-specific CNEs originated in these two clades.Conclusions:The tegu lizard genome together with the multiple genome alignment and comprehensive conserved element datasets provide a valuable resource for comparative genomic studies of reptiles and other amniotes.

Next-generation sequencers such as Illumina can now produce reads up to 300 bp with high throughput, which is attractive for genome assembly. A first step in genome assembly is to computationally correct sequencing errors. However, correcting all errors in these longer reads is challenging. Here, we show that reads with remaining errors after correction often overlap repeats, where short erroneous k-mers occur in other copies of the repeat. We developed an iterative error correction pipeline that runs the previously published String Graph Assembler (SGA) in multiple rounds of k-mer-based correction with an increasing k-mer size, followed by a final round of overlap-based correction. By combining the advantages of small and large k-mers, this approach corrects more errors in repeats and minimizes the total amount of erroneous reads. We show that higher read accuracy increases contig lengths two to three times. We provide SGA-Iteratively Correcting Errors (https://github.com/hillerlab/IterativeErrorCorrection/) that implements iterative error correction by using modules from SGA.

BACKGROUND:Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease.RESULTS:To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression".CONCLUSIONS:These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

To understand regulatory systems, it would be useful to uniformly determine how different components contribute to the expression of all other genes. We therefore monitored mRNA expression genome-wide, for individual deletions of one-quarter of yeast genes, focusing on (putative) regulators. The resulting genetic perturbation signatures reflect many different properties. These include the architecture of protein complexes and pathways, identification of expression changes compatible with viability, and the varying responsiveness to genetic perturbation. The data are assembled into a genetic perturbation network that shows different connectivities for different classes of regulators. Four feed-forward loop (FFL) types are overrepresented, including incoherent type 2 FFLs that likely represent feedback. Systematic transcription factor classification shows a surprisingly high abundance of gene-specific repressors, suggesting that yeast chromatin is not as generally restrictive to transcription as is often assumed. The data set is useful for studying individual genes and for discovering properties of an entire regulatory system.

Transcription plays a key role in cellular processes and its regulation is of paramount importance. The aim of the work described in this thesis is to study the transcription regulatory network of Saccharomyces cerevisiae, employing genome-wide approaches. All the three presented research studies have in common that individual genes are deleted and resulting gene expression changes are monitored by DNA microarrays. It is first described how gene expression changes can be used as detailed molecular phenotypes to study the transcription regulatory network underlying a signalling pathway. Genome-wide expression changes of 91 viable deletion mutants of different glucose signalling and metabolic pathways are analysed. A gene signature is used to group pathway members with similar effects on transcription. A new network approach is developed that is designed to explain gene expression changes upon deletion of one pathway member through the transcriptional regulation of another pathway member. This new approach reveals hierarchy and feedback in the transcription regulatory network. In particular, it predicts that the different glucose pathways converge on the transcriptional regulation of the biosynthesis of storage carbohydrates. For understanding the transcription regulatory network, it is important to discern target genes of gene-specific transcription factors (GSTFs). Genome-wide gene expression changes of 183 viable GSTF deletion mutants are compared with available DNA binding data, re-evaluating the overlap between both data types as well as discerning direct target genes of GSTFs. Besides determining roles of previously uncharacterized GSTFs, for example Stp3, this comparison has led to the first systematic classification of GSTFs into activators and/or repressors. Of all surveyed gene-specific transcription factors that could be classified, activators account for less than 54%. The remaining 46% of gene-specific transcription factors are repressors (37%) or have a dual function (9%). The unanticipated high number of gene-specific repressors indicates that in the yeast S. cerevisiae, chromatin is not as restrictive to transcription as has previously been thought and indicates that a considerable part of the gene-specific machinery is aimed at restricting unwanted transcription. Recent studies have systematically exposed large numbers of non-additive genetic interactions, the majority of which are functionally uncharacterized. To investigate such genetic interactions between GSTFs, we systematically analysed 72 viable double deletion mutants and compared them to the respective single deletion mutants. By generating a high-resolution gene expression atlas, epistatic effects of GSTF pairs on the expression of individual genes are investigated. Known genetic interactions are confirmed, and new ones are revealed. The analysis also provides evidence for two previously uncharacterized mechanisms, one for a negative (“buffering by induced dependency” between Hac1 and Rpn4) and one for a positive genetic interaction (“alleviation by derepression” between Gln3 and Gzf3). The study provides general insight into the complex nature of epistasis and proposes new models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

High levels of TFIID, a basal transcription factor, are found to be essential to induce and maintain the transcriptional program of pluripotent cells.