Search

Gene regulatory networks (GRNs) serve as useful abstractions to understand transcriptional dynamics in developmental systems. Computational prediction of GRNs has been successfully applied to genome-wide gene expression measurements with the advent of microarrays and RNA-sequencing. However, these inferred networks are inaccurate and mostly based on correlative rather than causative interactions. In this review, we highlight three approaches that significantly impact GRN inference: (1) moving from one genome-wide functional modality, gene expression, to multi-omics, (2) single cell sequencing, to measure cell type-specific signals and predict context-specific GRNs, and (3) neural networks as flexible models. Together, these experimental and computational developments have the potential to significantly impact the quality of inferred GRNs. Ultimately, accurately modeling the regulatory interactions between transcription factors and their target genes will be essential to understand the role of transcription factors in driving developmental gene expression programs and to derive testable hypotheses for validation.

The transcription factor p63 shares a high sequence identity with the tumour suppressor p53 which manifests itself in high structural similarity and preference for DNA sequences. Mutations in the DNA binding domain (DBD) of p53 have been studied in great detail, enabling a general mechanism-based classification. In this study we provide a detailed investigation of all currently known mutations in the p63 DBD, which are associated with developmental syndromes, by measuring their impact on transcriptional activity, DNA binding affinity, zinc binding capacity and thermodynamic stability. Some of the mutations we have further characterized with respect to their ability to convert human dermal fibroblasts into induced keratinocytes. Here we propose a classification of the p63 DBD mutations based on the four different mechanisms of DNA binding impairment which we identified: direct DNA contact, zinc finger region, H2 region, and dimer interface mutations. The data also demonstrate that, in contrast to p53 cancer mutations, no p63 mutation induces global unfolding and subsequent aggregation of the domain. The dimer interface mutations that affect the DNA binding affinity by disturbing the interaction between the individual DBDs retain partial DNA binding capacity which correlates with a milder patient phenotype.

Introduction: Mutations in the FOXE1 gene are implicated in cleft palate and thyroid dysgenesis in humans.Methods: To investigate whether zebrafish could provide meaningful insights into the etiology of developmental defects in humans related to FOXE1, we generated a zebrafish mutant that has a disruption in the nuclear localization signal in the foxe1 gene, thereby restraining nuclear access of the transcription factor. We characterized skeletal development and thyroidogenesis in these mutants, focusing on embryonic and larval stages.Results: Mutant larvae showed aberrant skeletal phenotypes in the ceratohyal cartilage and had reduced whole body levels of Ca, Mg and P, indicating a critical role for foxe1 in early skeletal development. Markers of bone and cartilage (precursor) cells were differentially expressed in mutants in post-migratory cranial neural crest cells in the pharyngeal arch at 1 dpf, at induction of chondrogenesis at 3 dpf and at the start of endochondral bone formation at 6 dpf. Foxe1 protein was detected in differentiated thyroid follicles, suggesting a role for the transcription factor in thyroidogenesis, but thyroid follicle morphology or differentiation were unaffected in mutants.Discussion: Taken together, our findings highlight the conserved role of Foxe1 in skeletal development and thyroidogenesis, and show differential signaling of osteogenic and chondrogenic genes related to foxe1 mutation.

The High-Altitude Water Cherenkov (HAWC) observatory is a second-generation continuously operated, wide field-of-view, TeV gamma-ray observatory. The HAWC observatory and its analysis techniques build on experience of the Milagro experiment in using ground-based water Cherenkov detectors for gamma-ray astronomy. HAWC is located on the Sierra Negra volcano in M\'exico at an elevation of 4100 meters above sea level. The completed HAWC observatory principal detector (HAWC) consists of 300 closely spaced water Cherenkov detectors, each equipped with four photomultiplier tubes to provide timing and charge information to reconstruct the extensive air shower energy and arrival direction. The HAWC observatory has been optimized to observe transient and steady emission from sources of gamma rays within an energy range from several hundred GeV to several hundred TeV. However, most of the air showers detected are initiated by cosmic rays, allowing studies of cosmic rays also to be performed. This paper describes the characteristics of the HAWC main array and its hardware., Comment: Accepted for publications in Nuclear Inst. and Methods in Physics Research, A (2023) 168253 ( https://www.sciencedirect.com/science/article/abs/pii/S0168900223002437 ); 39 pages, 14 Figures

Therapeutic aryl hydrocarbon receptor (AHR) modulating agents gained attention in dermatology as non-steroidal anti-inflammatory drugs that improve skin barrier properties. By exploiting AHR's known ligand promiscuity, we generated novel AHR modulating agents by lead optimization of a selective AHR modulator (SAhRM; SGA360). Twenty-two newly synthesized compounds were screened yielding two novel derivatives, SGA360f and SGA388, in which agonist activity led to enhanced keratinocyte terminal differentiation. SGA388 showed the highest agonist activity with potent normalization of keratinocyte hyperproliferation, restored expression of skin barrier proteins and dampening of chemokine expression by keratinocytes upon Th2-mediated inflammation in vitro. The topical application of SGA360f and SGA388 reduced acute skin inflammation in vivo by reducing cyclooxygenase levels, resulting in less neutrophilic dermal infiltrates. The minimal induction of cytochrome P450 enzyme activity, lack of cellular toxicity and mutagenicity classifies SGA360f and SGA388 as novel potential therapeutic AHR ligands and illustrates the potential of medicinal chemistry to fine-tune AHR signaling for the development of targeted therapies in dermatology and beyond.

MotivationAnalyzing a functional genomics experiment, such as ATAC-, ChIP-, or RNA-sequencing, requires genomic resources such as a reference genome assembly and gene annotation. These data can generally be retrieved from different organizations and in different versions. Most bioinformatic workflows require the user to supply this genomic data manually, which can be a tedious and error-prone process.ResultsHere, we present genomepy, which can search, download, and preprocess the right genomic data for your analysis. Genomepy can search genomic data on NCBI, Ensembl, UCSC, and GENCODE, and inspect available gene annotations to enable an informed decision. The selected genome and gene annotation can be downloaded and preprocessed with sensible, yet controllable, defaults. Additional supporting data can be automatically generated or downloaded, such as aligner indexes, genome metadata, and blacklists.Availability and implementationGenomepy is freely available at https://github.com/vanheeringen-lab/genomepy under the MIT license and can be installed through pip or Bioconda.

Advanced computational methods exploit gene expression and epigenetic datasets to predict gene regulatory networks controlled by transcription factors (TFs). These methods have identified cell fate determining TFs but require large amounts of reference data and experimental expertise. Here, we present an easy to use network-based computational framework that exploits enhancers defined by bidirectional transcription, using as sole input CAGE sequencing data to correctly predict TFs key to various human cell types. Next, we applied this Analysis Algorithm for Networks Specified by Enhancers based on CAGE (ANANSE-CAGE) to predict TFs driving red and white blood cell development, and THP-1 leukemia cell immortalization. Further, we predicted TFs that are differentially important to either cell line- or primary- associated MLL-AF9-driven gene programs, and in primary MLL-AF9 acute leukemia. Our approach identified experimentally validated as well as thus far unexplored TFs in these processes. ANANSE-CAGE will be useful to identify transcription factors that are key to any cell fate change using only CAGE-seq data as input.

Notch1 is a crucial oncogenic driver in T-cell acute lymphoblastic leukemia (T-ALL), making it an attractive therapeutic target. However, the success of targeted therapy using γ-secretase inhibitors (GSIs), small molecules blocking Notch cleavage and subsequent activation, has been limited due to development of resistance, thus restricting its clinical efficacy. Here, we systematically compare GSI resistant and sensitive cell states by quantitative mass spectrometry-based phosphoproteomics, using complementary models of resistance, including T-ALL patient-derived xenografts (PDX) models. Our datasets reveal common mechanisms of GSI resistance, including a distinct kinase signature that involves protein kinase C delta. We demonstrate that the PKC inhibitor sotrastaurin enhances the anti-leukemic activity of GSI in PDX models and completely abrogates the development of acquired GSI resistance in vitro. Overall, we highlight the potential of proteomics to dissect alterations in cellular signaling and identify druggable pathways in cancer., NOTCH1 is a driver of T-cell acute lymphoblastic leukemia that can be inhibited by γ-secretase inhibitors (GSIs), but their clinical efficacy is limited. Here, the authors compare the phosphoproteomes of GSI resistant and sensitive models, and identify potential kinase targets to overcome GSI resistance.

Summary Polycomb Repressive Complex 2 (PRC2) plays an essential role in gene repression during development, catalyzing H3 lysine 27 trimethylation (H3K27me3). MTF2 in the PRC2.1 sub-complex, and JARID2 in PRC2.2, are central in core PRC2 recruitment to target genes in mouse embryonic stem cells (mESCs). To investigate how PRC2.1 and PRC2.2 cooperate, we combined Polycomb mutant mESCs with chemical inhibition of binding to H3K27me3. We find that PRC2.1 and PRC2.2 mediate two distinct paths for recruitment, which are mutually reinforced. Whereas PRC2.1 recruitment is mediated by MTF2 binding to DNA, JARID2-containing PRC2.2 recruitment is more dependent on PRC1. Both recruitment axes are supported by core subunit EED binding to H3K27me3, but EED inhibition exhibits a more pronounced effect in Jarid2 null cells. Finally, we show that PRC1 and PRC2 enhance reciprocal binding. Together, these data disentangle the interdependent interactions that are important for PRC2 recruitment., Graphical Abstract, Highlights • Systematic analysis of PRC2 genomic binding in mouse ESCs • Mutations and chemical inhibition uncover roles of EED and other interactions • Contributions of MTF2 to PRC2.1 binding and of JARID2 and PRC1 to PRC2.2 binding • PRC1, PRC2.1, and PRC2.2 are all mutually interdependent for binding to chromatin, Veenstra and Marks and colleagues systematically analyze the contributions of Polycomb Repressive Complexes PRC2.1, PRC2.2, and PRC1 in ESCs, using null mutations in combination with chemical inhibition of EED, the core subunit that interacts with H3K27me3. The mutations and the EED inhibition uncover the functional overlap, and the compensatory and highly interdependent binding of PRC2.1, PRC2.2, and PRC1.

SummaryMesendodermal specification is one of the earliest events in embryogenesis, where cells first acquire distinct identities. Cell differentiation is a highly regulated process that involves the function of numerous transcription factors (TFs) and signaling molecules, which can be described with gene regulatory networks (GRNs). Cell differentiation GRNs are difficult to build because existing mechanistic methods are low-throughput, and high-throughput methods tend to be non-mechanistic. Additionally, integrating highly dimensional data comprised of more than two data types is challenging. Here, we use linked self-organizing maps to combine ChIP-seq/ATAC-seq with temporal, spatial and perturbation RNA-seq data fromXenopus tropicalismesendoderm development to build a high resolution genome scale mechanistic GRN. We recovered both known and previously unsuspected TF-DNA/TF-TF interactions and validated through reporter assays. Our analysis provides new insights into transcriptional regulation of early cell fate decisions and provides a general approach to building GRNs using highly-dimensional multi-omic data sets.HighlightsBuilt a generally applicable pipeline to creating GRNs using highly-dimensional multi-omic data setsPredicted new TF-DNA/TF-TF interactions during mesendoderm developmentGenerate the first genome scale GRN for vertebrate mesendoderm and expanded the core mesendodermal developmental network with high fidelityDeveloped a resource to visualize hundreds of RNA-seq and ChIP-seq data using 2D SOM metaclusters.

Concurrent mutation of a RAS oncogene and the tumor suppressor p53 is common in tumorigenesis, and inflammation can promote RAS-driven tumorigenesis without the need to mutate p53. Here, we show, using a well-established mutant RAS and an inflammation-driven mouse skin tumor model, that loss of the p53 inhibitor iASPP facilitates tumorigenesis. Specifically, iASPP regulates expression of a subset of p63 and AP1 targets, including genes involved in skin differentiation and inflammation, suggesting that loss of iASPP in keratinocytes supports a tumor-promoting inflammatory microenvironment. Mechanistically, JNK-mediated phosphorylation regulates iASPP function and inhibits iASPP binding with AP1 components, such as JUND, via PXXP/SH3 domain-mediated interaction. Our results uncover a JNK-iASPP-AP1 regulatory axis that is crucial for tissue homeostasis. We show that iASPP is a tumor suppressor and an AP1 coregulator.

Fibroblast growth factor-2 (FGF2) has multiple roles in cutaneous wound healing but its natural low stability prevents the development of its use in skin repair therapies. Here we show that FGF2 binds the outer surface of dermal fibroblast (DF)-derived extracellular vesicles (EVs) and this association protects FGF2 from fast degradation. EVs isolated from DF cultured in the presence of FGF2 harbor FGF2 on their surface and FGF2 can bind purified EVs in absence of cells. Remarkably, FGF2 binding to EVs is restricted to a specific subpopulation of EVs, which do not express CD63 and CD81 markers. Treatment of DF with FGF2-EVs activated ERK and STAT signaling pathways and increased cell proliferation and migration. Local injection of FGF2-EVs improved wound healing in mice. We further demonstrated that binding to EVs protects FGF2 from both thermal and proteolytic degradation, thus maintaining FGF2 function. This suggests that EVs protect soluble factors from degradation and increase their stability and half-life. These results reveal a novel aspect of EV function and suggest EVs as a potential tool for delivering FGF2 in skin healing therapies.

Functional human atrial myocytes derived from foetal atrial tissue can be massively expanded via a conditional cell-immortalization method, and used in in vitro models of atrial fibrillation.The lack of a scalable and robust source of well-differentiated human atrial myocytes constrains the development of in vitro models of atrial fibrillation (AF). Here we show that fully functional atrial myocytes can be generated and expanded one-quadrillion-fold via a conditional cell-immortalization method relying on lentiviral vectors and the doxycycline-controlled expression of a recombinant viral oncogene in human foetal atrial myocytes, and that the immortalized cells can be used to generate in vitro models of AF. The method generated 15 monoclonal cell lines with molecular, cellular and electrophysiological properties resembling those of primary atrial myocytes. Multicellular in vitro models of AF generated using the immortalized atrial myocytes displayed fibrillatory activity (with activation frequencies of 6-8 Hz, consistent with the clinical manifestation of AF), which could be terminated by the administration of clinically approved antiarrhythmic drugs. The conditional cell-immortalization method could be used to generate functional cell lines from other human parenchymal cells, for the development of in vitro models of human disease.

Contains fulltext : 247923.pdf (Publisher’s version ) (Open Access)

Contains fulltext : 249539.pdf (Publisher’s version ) (Open Access) Radboud University, 16 mei 2022 Promotor : Veenstra, G.J.C. Co-promotores : Zhou, H., Heeringen, S.J. van 233 p.

Embryonic development is a highly intricate and complex process. Different regulatory mechanisms cooperatively dictate the fate of cells as they progress from pluripotent stem cells to terminally differentiated cell types in tissues. A crucial regulator of these processes is the Polycomb Repressive Complex 2 (PRC2). By catalyzing the mono-, di-, and tri-methylation of lysine residues on histone H3 tails (H3K27me3), PRC2 compacts chromatin by cooperating with Polycomb Repressive Complex 1 (PRC1) and represses transcription of target genes. Proteomic and biochemical studies have revealed two variant complexes of PRC2, namely PRC2.1 which consists of the core proteins (EZH2, SUZ12, EED, and RBBP4/7) interacting with one of the Polycomb-like proteins (MTF2, PHF1, PHF19), and EPOP or PALI1/2, and PRC2.2 which contains JARID2 and AEBP2 proteins. MTF2 and JARID2 have been discovered to have crucial roles in directing and recruiting PRC2 to target genes for repression in embryonic stem cells (ESCs). Following these findings, recent work in the field has begun to explore the roles of different PRC2 variant complexes during different stages of embryonic development, by examining molecular phenotypes of PRC2 mutants in both in vitro (2D and 3D differentiation) and in vivo (knock-out mice) assays, analyzed with modern single-cell omics and biochemical assays. In this review, we discuss the latest findings that uncovered the roles of different PRC2 proteins during cell-fate and lineage specification and extrapolate these findings to define a developmental roadmap for different flavors of PRC2 regulation during mammalian embryonic development.

Contains fulltext : 285239.pdf (Publisher’s version ) (Open Access)

Single-cell omics such as single-cell RNA-sequencing (RNA-seq) have been used extensively to obtain single-cell genome-wide expression data. This technique can be used to compare mutant and wild-type embryos at predifferentiation stages when individual tissues are not yet formed (therefore requiring genotyping to distinguish among embryos), for example, to determine effects of mutations on developmental trajectories or congenital disease phenotypes. It is, however, hard to use single cells for this technique, because such embryos or cells would need to be frozen until genotyping is complete to capture a given developmental stage precisely, but intact cells cannot be isolated from frozen samples. We developed a protocol in which high-quality nuclei are isolated from frozen cell suspensions, allowing for genotyping individual embryos based on a small fraction of a single embryo suspension. The remaining suspension is frozen. After genotyping is complete, nuclei are isolated from embryo suspensions with the desired genotype and encapsulated in 10× Genomics barcoded gel beads for single-nucleus RNA-seq. We provide a step-by-step protocol that can be used for single transcriptomic analysis as well as single-nucleus chromatin accessibility assays such as ATAC-seq. This technique allows for high-quality high-throughput single-nucleus analysis of gene expression in genotyped embryos. This approach may also be valuable for collection of wild-type embryonic material, for example, when collecting tissue from a particular developmental stage. In addition, freezing of tissue suspensions allows precise staging of collected embryos or tissue that may be difficult to manage when collecting and processing cells from living embryos for single-cell RNA-seq.

Contains fulltext : 284055.pdf (Publisher’s version ) (Open Access)

Polycomb Repressive Complex 2 (PRC2) is crucial for the coordinated expression of genes during early embryonic development, catalyzing histone H3 lysine 27 trimethylation. Two distinct PRC2 complexes, PRC2.1 and PRC2.2, contain respectively MTF2 and JARID2 in embryonic stem cells (ESCs). In this study, we explored their roles in lineage specification and commitment, using single-cell transcriptomics and mouse embryoid bodies derived from Mtf2 and Jarid2 null ESCs. We observe that the loss of Mtf2 results in enhanced and faster differentiation towards cell fates from all germ layers, while the Jarid2 null cells are predominantly directed towards early differentiating precursors, with reduced efficiency towards mesendodermal lineages. These effects are caused by derepression of developmental regulators that are poised for activation in pluripotent cells and gain H3K4me3 at their promoters in the absence of PRC2 repression. Upon lineage commitment, the differentiation trajectories are relatively similar to those of wild-type cells. Together, our results uncover a major role for MTF2-containing PRC2.1 in balancing poised lineage-specific gene activation, whereas the contribution of JARID2-containing PRC2 is more selective in nature compared to MTF2. These data explain how PRC2 imposes thresholds for lineage choice during the exit of pluripotency., Polycomb Repressive Complex 2, an important regulator of embryonic development, exists in two configurations, PRC2.1 and PRC2.2. Here the authors dissect the functional contributions of these two PRC2 subunits and observed complex-specific alterations in the cell state of pluripotent and early differentiating cells.

Contains fulltext : 248251.pdf (Publisher’s version ) (Open Access) CRISPR-Cas9 is the most straightforward genome-editing tool to date. However, its implementation across disciplines is hampered by variable genome-editing efficiencies, reduced cell viability, and low success rates in obtaining clonal cell lines. This review aims to recognize all CRISPR-Cas9‒related work within the experimental dermatology field to identify key factors for successful strategies in the different keratinocyte (KC) cell sources available. On the basis of these findings, we conclude that most groups use immortalized KCs for generating knockout KCs. Our critical considerations for future studies using CRISPR-Cas9, both for fundamental and clinical applications, may guide implementation strategies of CRISPR-Cas9 technologies in the (experimental) dermatology field.

Contains fulltext : 204872.pdf (Publisher’s version ) (Open Access) We aimed to identify novel deletions and variants of TP63 associated with orofacial clefting (OFC). Copy number variants were assessed in three OFC families using microarray analysis. Subsequently, we analyzed TP63 in a cohort of 1072 individuals affected with OFC and 706 population-based controls using molecular inversion probes (MIPs). We identified partial deletions of TP63 in individuals from three families affected with OFC. In the OFC cohort, we identified several TP63 variants predicting to cause loss-of-function alleles, including a frameshift variant c.569_576del (p.(Ala190Aspfs*5)) and a nonsense variant c.997C>T (p.(Gln333*)) that introduces a premature stop codon in the DNA-binding domain. In addition, we identified the first missense variants in the oligomerization domain c.1213G>A (p.(Val405Met)), which occurred in individuals with OFC. This variant was shown to abrogate oligomerization of mutant p63 protein into oligomeric complexes, and therefore likely represents a loss-of-function allele rather than a dominant-negative. All of these variants were inherited from an unaffected parent, suggesting reduced penetrance of such loss-of-function alleles. Our data indicate that loss-of-function alleles in TP63 can also give rise to OFC as the main phenotype. We have uncovered the dosage-dependent functions of p63, which were previously rejected.

Summary: Regenerative responses predispose tissues to tumor formation by largely unknown mechanisms. High-mobility group box 1 (HMGB1) is a danger-associated molecular pattern contributing to inflammatory pathologies. We show that HMGB1 derived from keratinocytes, but not myeloid cells, delays cutaneous wound healing and drives tumor formation. In wounds of mice lacking HMGB1 selectively in keratinocytes, a marked reduction in neutrophil extracellular trap (NET) formation is observed. Pharmacological targeting of HMGB1 or NETs prevents skin tumorigenesis and accelerates wound regeneration. HMGB1-dependent NET formation and skin tumorigenesis is orchestrated by tumor necrosis factor (TNF) and requires RIPK1 kinase activity. NETs are present in the microenvironment of keratinocyte-derived tumors in mice and lesional and tumor skin of patients suffering from recessive dystrophic epidermolysis bullosa, a disease in which skin blistering predisposes to tumorigenesis. We conclude that tumorigenicity of the wound microenvironment depends on epithelial-derived HMGB1 regulating NET formation, thereby establishing a mechanism linking reparative inflammation to tumor initiation. : Inflammation mediates tissue repair but can be hijacked to promote tumorigenesis. Hoste et al. demonstrate that HMGB1 delays regeneration and drives tumor formation in skin by recruitment and priming of neutrophils. These data indicate that therapies targeting HMGB1 or NET formation might be relevant in chronic and diabetic wound treatment. Keywords: tumor microenvironment, skin inflammation, innate immunity, neutrophil extracellular traps, HMGB1, TNF, wound healing, diabetes, epidermolysis bullosa

Contains fulltext : 231651.pdf (Publisher’s version ) (Open Access) In vitro neuronal models are essential for studying neurological physiology, disease mechanisms and potential treatments. Most in vitro models lack controlled vasculature, despite its necessity in brain physiology and disease. Organ-on-chip models offer microfluidic culture systems with dedicated micro-compartments for neurons and vascular cells. Such multi-cell type organs-on-chips can emulate neurovascular unit (NVU) physiology, however there is a lack of systematic data on how individual cell types are affected by culturing on microfluidic systems versus conventional culture plates. This information can provide perspective on initial findings of studies using organs-on-chip models, and further optimizes these models in terms of cellular maturity and neurovascular physiology. Here, we analysed the transcriptomic profiles of co-cultures of human induced pluripotent stem cell (hiPSC)-derived neurons and rat astrocytes, as well as one-day monocultures of human endothelial cells, cultured on microfluidic chips. For each cell type, large gene expression changes were observed when cultured on microfluidic chips compared to conventional culture plates. Endothelial cells showed decreased cell division, neurons and astrocytes exhibited increased cell adhesion, and neurons showed increased maturity when cultured on a microfluidic chip. Our results demonstrate that culturing NVU cell types on microfluidic chips changes their gene expression profiles, presumably due to distinct surface-to-volume ratios and substrate materials. These findings inform further NVU organ-on-chip model optimization and support their future application in disease studies and drug testing.

Summary To further our understanding of how biochemical information flows through cells upon external stimulation, we require single-cell multi-omics methods that concurrently map changes in (phospho)protein levels across signaling networks and the associated gene expression profiles. Here, we present quantification of RNA and intracellular epitopes by sequencing (QuRIE-seq), a droplet-based platform for single-cell RNA and intra- and extracellular (phospho)protein quantification through sequencing. We applied QuRIE-seq to quantify cell-state changes at both the signaling and the transcriptome level after 2-, 4-, 6-, 60-, and 180-min stimulation of the B cell receptor pathway in Burkitt lymphoma cells. Using the multi-omics factor analysis (MOFA+) framework, we delineated changes in single-cell (phospho)protein and gene expression patterns over multiple timescales and revealed the effect of an inhibitory drug (ibrutinib) on signaling and gene expression landscapes.

Summary Retinoic acid (RA) signaling is an important and conserved pathway that regulates cellular proliferation and differentiation. Furthermore, perturbed RA signaling is implicated in cancer initiation and progression. However, the mechanisms by which RA signaling contributes to homeostasis, malignant transformation, and disease progression in the intestine remain incompletely understood. Here, we report, in agreement with previous findings, that activation of the Retinoic Acid Receptor and the Retinoid X Receptor results in enhanced transcription of enterocyte-specific genes in mouse small intestinal organoids. Conversely, inhibition of this pathway results in reduced expression of genes associated with the absorptive lineage. Strikingly, this latter effect is conserved in a human organoid model for colorectal cancer (CRC) progression. We further show that RXR motif accessibility depends on progression state of CRC organoids. Finally, we show that reduced RXR target gene expression correlates with worse CRC prognosis, implying RA signaling as a putative therapeutic target in CRC., Graphical abstract, Highlights • RA signaling contributes to enterocyte differentiation in murine intestinal organoids • Inhibition of RXR decreases enterocyte gene expression in colon cancer organoids • Accessibility of RXR motifs correlates with RXRi susceptibility • High expression of RA signaling targets correlates with higher CRC patient survival, Biological sciences; Cancer; Cell biology; Functional aspects of cell biology; Oncology

Contains fulltext : 236639.pdf (Publisher’s version ) (Open Access) Radboud University, 08 september 2021 Promotor : Veenstra, G.J.C. 141 p.

Simple Summary The protein p63 belongs to the family of the p53 tumor suppressor. Mouse models have, however, shown that it is not a classical tumor suppressor but instead involved in developmental processes. Mutations in the p63 gene cause several developmental defects in human patients characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia due to p63’s role as a master regulator of epidermal development. In addition, p63 plays a key role as a quality control factor in oocytes and p63 mutations can result either in compromised genetic quality control or premature cell death of all oocytes. Abstract The p63 gene encodes a master regulator of epidermal commitment, development, and differentiation. Heterozygous mutations in the DNA binding domain cause Ectrodactyly, Ectodermal Dysplasia, characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia while mutations in in the C-terminal domain of the α-isoform cause Ankyloblepharon-Ectodermal defects-Cleft lip/palate (AEC) syndrome, a life-threatening disorder characterized by skin fragility, severe, long-lasting skin erosions, and cleft lip/palate. The molecular disease mechanisms of these syndromes have recently become elucidated and have enhanced our understanding of the role of p63 in epidermal development. Here we review the molecular cause and functional consequences of these p63-mutations for skin development and discuss the consequences of p63 mutations for female fertility.

During vertebrate gastrulation, mesoderm is induced in pluripotent cells, concomitant with dorsal-ventral patterning and establishing of the dorsal axis. We applied single-cell chromatin accessibility and transcriptome analyses to explore the emergence of cellular heterogeneity during gastrulation in Xenopus tropicalis. Transcriptionally inactive lineage-restricted genes exhibit relatively open chromatin in animal caps, whereas chromatin accessibility in dorsal marginal zone cells more closely reflects transcriptional activity. We characterized single-cell trajectories and identified head and trunk organizer cell clusters in early gastrulae. By integrating chromatin accessibility and transcriptome data, we inferred the activity of transcription factors in single-cell clusters and tested the activity of organizer-expressed transcription factors in animal caps, alone or in combination. The expression profile induced by a combination of Foxb1 and Eomes most closely resembles that observed in the head organizer. Genes induced by Eomes, Otx2, or the Irx3-Otx2 combination are enriched for maternally regulated H3K4me3 modifications, whereas Lhx8-induced genes are marked more frequently by zygotically controlled H3K4me3. Taken together, our results show that transcription factors cooperate in a combinatorial fashion in generally open chromatin to orchestrate zygotic gene expression.

Contains fulltext : 244255.pdf (Publisher’s version ) (Open Access) Tissue engineering (TE) of long tracheal segments is conceptually appealing for patients with inoperable tracheal pathology. In tracheal TE, stem cells isolated from bone marrow or adipose tissue have been employed, but the ideal cell source has yet to be determined. When considering the origin of stem cells, cells isolated from a source embryonically related to the trachea may be more similar. In this study, we investigated the feasibility of isolating progenitor cells from pleura and pericard as an alternative cells source for tracheal tissue engineering. Porcine progenitor cells were isolated from pleura, pericard, trachea and adipose tissue and expanded in culture. Isolated cells were characterized by PCR, RNA sequencing, differentiation assays and cell survival assays and were compared to trachea and adipose-derived progenitor cells. Progenitor-like cells were successfully isolated and expanded from pericard and pleura as indicated by gene expression and functional analyses. Gene expression analysis and RNA sequencing showed a stem cell signature indicating multipotency, albeit that subtle differences between different cell sources were visible. Functional analysis revealed that these cells were able to differentiate towards chondrogenic, osteogenic and adipogenic lineages. Isolation of progenitor cells from pericard and pleura with stem cell features is feasible. Although functional differences with adipose-derived stem cells were limited, based on their gene expression, pericard- and pleura-derived stem cells may represent a superior autologous cell source for cell seeding in tracheal tissue engineering.

The epidermal compartment of the skin is regenerated constantly by proliferation of epidermal keratinocytes. Differentiation of a subset of these keratinocytes allows the epidermis to retain its barrier properties. Regulation of keratinocyte fate—whether to remain proliferative or terminally differentiate—is complex and not fully understood. The objective of our study was to assess if DNA methylation changes contribute to the regulation of keratinocyte fate. We employed genome‐wide MethylationEPIC beadchip array measuring approximately 850 000 probes combined with RNA sequencing of in vitro cultured non‐differentiated and terminally differentiated adult human primary keratinocytes. We did not observe a correlation between methylation status and transcriptome changes. Moreover, only two differentially methylated probes were detected, of which one was located in the TRIM29 gene. Although TRIM29 knock‐down resulted in lower expression levels of terminal differentiation genes, these changes were minor. From these results, we conclude that—in our in vitro experimental setup—it is unlikely that changes in DNA methylation have an important regulatory role in terminal keratinocyte differentiation.

Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.

Contains fulltext : 215040.pdf (Publisher’s version ) (Open Access) Radboud University, 09 januari 2020 Promotor : Veenstra, G.J.C. Co-promotor : Mulder, K.W. 179 p.

Motivation Cellular identity and behavior is controlled by complex gene regulatory networks. Transcription factors (TFs) bind to specific DNA sequences to regulate the transcription of their target genes. On the basis of these TF motifs in cis-regulatory elements we can model the influence of TFs on gene expression. In such models of TF motif activity the data is usually modeled assuming a linear relationship between the motif activity and the gene expression level. A commonly used method to model motif influence is based on Ridge Regression. One important assumption of linear regression is the independence between samples. However, if samples are generated from the same cell line, tissue, or other biological source, this assumption may be invalid. This same assumption of independence is also applied to different yet similar experimental conditions, which may also be inappropriate. In theory, the independence assumption between samples could lead to loss in signal detection. Here we investigate whether a Bayesian model that allows for correlations results in more accurate inference of motif activities. Results We extend the Ridge Regression to a Bayesian Linear Mixed Model, which allows us to model dependence between different samples. In a simulation study, we investigate the differences between the two model assumptions. We show that our Bayesian Linear Mixed Model implementation outperforms Ridge Regression in a simulation scenario where the noise, which is the signal that can not be explained by TF motifs, is uncorrelated. However, we demonstrate that there is no such gain in performance if the noise has a similar covariance structure over samples as the signal that can be explained by motifs. We give a mathematical explanation to why this is the case. Using four representative real datasets we show that at most ∼​40% of the signal is explained by motifs using the linear model. With these data there is no advantage to using the Bayesian Linear Mixed Model, due to the similarity of the covariance structure. Availability & implementation The project implementation is available at https://github.com/Sim19/SimGEXPwMotifs.

High-throughput sequencing methods have created exciting opportunities to explore the regulatory landscape of the entire genome. Here we introduce methods to characterize the genomic locations of bound proteins, open chromatin, and sites of DNA–DNA contact in Xenopus embryos. These methods include chromatin immunoprecipitation followed by sequencing (ChIP-seq), a combination of DNase I digestion and sequencing (DNase-seq), the assay for transposase-accessible chromatin and sequencing (ATAC-seq), and the use of proximity-based DNA ligation followed by sequencing (Hi-C).

Background Estimating the genetic component of a complex phenotype is a complicated problem, mainly because there are many allele effects to estimate from a limited number of phenotypes. In spite of this difficulty, linear methods with variable selection have been able to give good predictions of additive effects of individuals. However, prediction of non-additive genetic effects is challenging with the usual prediction methods. In machine learning, non-additive relations between inputs can be modeled with neural networks. We developed a novel method (NetSparse) that uses Bayesian neural networks with variable selection for the prediction of genotypic values of individuals, including non-additive genetic effects. Results We simulated several populations with different phenotypic models and compared NetSparse to genomic best linear unbiased prediction (GBLUP), BayesB, their dominance variants, and an additive by additive method. We found that when the number of QTL was relatively small (10 or 100), NetSparse had 2 to 28 percentage points higher accuracy than the reference methods. For scenarios that included dominance or epistatic effects, NetSparse had 0.0 to 3.9 percentage points higher accuracy for predicting phenotypes than the reference methods, except in scenarios with extreme overdominance, for which reference methods that explicitly model dominance had 6 percentage points higher accuracy than NetSparse. Conclusions Bayesian neural networks with variable selection are promising for prediction of the genetic component of complex traits in animal breeding, and their performance is robust across different genetic models. However, their large computational costs can hinder their use in practice.

Contains fulltext : 225108.pdf (Publisher’s version ) (Open Access) Radboud University, 29 oktober 2020 Promotor : Veenstra, G.J.C. Co-promotor : Zhou, H. 274 p.

Summary Epidermal homeostasis requires balanced and coordinated adult stem cell renewal and differentiation. These processes are controlled by both extracellular signaling and by cell intrinsic transcription regulatory networks, yet how these control mechanisms are integrated to achieve this is unclear. Here, we developed single-cell Immuno-Detection by sequencing (scID-seq) and simultaneously measured 69 proteins (including 34 phosphorylated epitopes) at single-cell resolution to study the activation state of signaling pathways during human epidermal differentiation. Computational pseudo-timing inference revealed dynamic activation of the JAK-STAT, WNT, and BMP pathways along the epidermal differentiation trajectory. We found that during differentiation, cells start producing BMP2-ligands and activate the canonical intracellular effectors SMAD1/5/9. Mechanistically, the BMP pathway is responsible for activating the MAF/MAFB/ZNF750 transcription factor network to drive late-stage epidermal differentiation. Our work indicates that incorporating signaling pathway activation into this transcription regulatory network enables coordination of transcription programs during epidermal differentiation., Graphical Abstract, Highlights • scID-seq allows quantification of 70 (phospho-)proteins at single-cell level • Pseudo-time inference reveals signaling dynamics during epidermal differentiation • BMP signaling drives a late differentiation transcription program • BMP signaling activates the MAF/MAFB/ZNF750 transcription factor network, Stem Cells Research; Systems Biology; Proteomics

Contains fulltext : 190857.pdf (Publisher’s version ) (Closed access) 23 p.

Hybridization of eggs and sperm from closely related species can give rise to genetic diversity, or can lead to embryo inviability owing to incompatibility. Although central to evolution, the cellular and molecular mechanisms underlying post-zygotic barriers that drive reproductive isolation and speciation remain largely unknown. Species of the African clawed frog Xenopus provide an ideal system to study hybridization and genome evolution. Xenopus laevis is an allotetraploid with 36 chromosomes that arose through interspecific hybridization of diploid progenitors, whereas Xenopus tropicalis is a diploid with 20 chromosomes that diverged from a common ancestor approximately 48 million years ago. Differences in genome size between the two species are accompanied by organism size differences, and size scaling of the egg and subcellular structures such as nuclei and spindles formed in egg extracts. Nevertheless, early development transcriptional programs, gene expression patterns, and protein sequences are generally conserved. Whereas the hybrid produced when X. laevis eggs are fertilized by X. tropicalis sperm is viable, the reverse hybrid dies before gastrulation. Here we apply cell biological tools and high-throughput methods to study the mechanisms underlying hybrid inviability. We reveal that two specific X. laevis chromosomes are incompatible with the X. tropicalis cytoplasm and are mis-segregated during mitosis, leading to unbalanced gene expression at the maternal to zygotic transition, followed by cell-autonomous catastrophic embryo death. These results reveal a cellular mechanism underlying hybrid incompatibility that is driven by genome evolution and contributes to the process by which biological populations become distinct species.

Contains fulltext : 195277.pdf (Publisher’s version ) (Open Access)

International audience; The chromatin state is finely tuned to regulate function and specificity for transcription factors such as oestrogen receptor alpha (ER), which contributes to cell growth in breast cancer. ER transcriptional potential is mediated, in large part, by the specific associated proteins and co-factors that interact with it. Despite the identification and characterisation of several ER coregulators, a complete and systematic view of ER-regulating chromatin modifiers is lacking. By exploiting a focused siRNA screen that investigated the requirement for a library of 330 chromatin regulators in ER-mediated cell growth, we find that the NuRD and coREST histone deacetylation complexes are critical for breast cancer cell proliferation. Further, by proteomic and genomics approaches, we discover the transcription factor TRPS1 to be a key interactor of the NuRD and coREST complexes. Interestingly, TRPS1 gene amplification occurs in 28% of human breast tumours and is associated with poor prognosis. We propose that TRPS1 is required to repress spurious binding of ER, where it contributes to the removal of histone acetylation. Our data suggest that TRPS1 is an important ER-associated transcriptional repressor that regulates cell proliferation, chromatin acetylation and ER binding at the chromatin of cis-regulatory elements.

The prefrontal cortex (PFC) is one of the latest brain regions to mature, which allows the acquisition of complex cognitive abilities through experience. To unravel the underlying gene expression changes during postnatal development, we performed RNA-sequencing (RNA-seq) in the rat medial PFC (mPFC) at five developmental time points from infancy to adulthood, and analyzed the differential expression of protein-coding genes, long intergenic noncoding RNAs (lincRNAs), and alternative exons. We showed that most expression changes occur in infancy, and that the number of differentially expressed genes reduces toward adulthood. We observed 137 differentially expressed lincRNAs and 796 genes showing alternative exon usage during postnatal development. Importantly, we detected a genetic switch from neuronal network establishment in infancy to maintenance of neural networks in adulthood based on gene expression dynamics, involving changes in protein-coding and lincRNA gene expression as well as alternative exon usage. Our gene expression datasets provide insights into the multifaceted transcriptional regulation of the developing PFC. They can be used to study the basic developmental processes of the mPFC and to understand the mechanisms of neurodevelopmental and neuropsychiatric disorders. Our study provides an important contribution to the ongoing efforts to complete the "brain map", and to the understanding of PFC development.

Contains fulltext : 174194.pdf (Publisher’s version ) (Open Access) Keratolytic winter erythema (KWE) is a rare autosomal-dominant skin disorder characterized by recurrent episodes of palmoplantar erythema and epidermal peeling. KWE was previously mapped to 8p23.1-p22 (KWE critical region) in South African families. Using targeted resequencing of the KWE critical region in five South African families and SNP array and whole-genome sequencing in two Norwegian families, we identified two overlapping tandem duplications of 7.67 kb (South Africans) and 15.93 kb (Norwegians). The duplications segregated with the disease and were located upstream of CTSB, a gene encoding cathepsin B, a cysteine protease involved in keratinocyte homeostasis. Included in the 2.62 kb overlapping region of these duplications is an enhancer element that is active in epidermal keratinocytes. The activity of this enhancer correlated with CTSB expression in normal differentiating keratinocytes and other cell lines, but not with FDFT1 or NEIL2 expression. Gene expression (qPCR) analysis and immunohistochemistry of the palmar epidermis demonstrated significantly increased expression of CTSB, as well as stronger staining of cathepsin B in the stratum granulosum of affected individuals than in that of control individuals. Analysis of higher-order chromatin structure data and RNA polymerase II ChIA-PET data from MCF-7 cells did not suggest remote effects of the enhancer. In conclusion, KWE in South African and Norwegian families is caused by tandem duplications in a non-coding genomic region containing an active enhancer element for CTSB, resulting in upregulation of this gene in affected individuals.

Whole-transcriptome or RNA sequencing (RNA-Seq) is a powerful and versatile tool for functional analysis of different types of RNA molecules, but sample reagent and sequencing cost can be prohibitive for hypothesis-driven studies where the aim is to quantify differential expression of a limited number of genes. Here we present an approach for quantification of differential mRNA expression by targeted resequencing of complementary DNA using single-molecule molecular inversion probes (cDNA-smMIPs) that enable highly multiplexed resequencing of cDNA target regions of ∼100 nucleotides and counting of individual molecules. We show that accurate estimates of differential expression can be obtained from molecule counts for hundreds of smMIPs per reaction and that smMIPs are also suitable for quantification of relative gene expression and allele-specific expression. Compared with low-coverage RNA-Seq and a hybridization-based targeted RNA-Seq method, cDNA-smMIPs are a cost-effective high-throughput tool for hypothesis-driven expression analysis in large numbers of genes (10 to 500) and samples (hundreds to thousands)., Transcriptome sequencing is a powerful tool for functional analysis of different types of RNA molecules in a wide range of applications. Here the authors use targeted resequencing of cDNA with single-molecule molecular inversion probes as a cost-effective, high-throughput tool for mRNA quantification.