scANANSE gene regulatory network and motif analysis of single-cell clusters.

The recent development of single-cell techniques is essential to unravel complex biological systems. By measuring the transcriptome and the accessible genome on a single-cell level, cellular heterogeneity in a biological environment can be deciphered. Transcription factors act as key regulators activating and repressing downstream target genes, and together they constitute gene regulatory networks that govern cell morphology and identity. Dissecting these gene regulatory networks is crucial for understanding molecular mechanisms and disease, especially within highly complex biological systems. The gene regulatory network analysis software ANANSE and the motif enrichment software GimmeMotifs were both developed to analyse bulk datasets. We developed scANANSE, a software pipeline for gene regulatory network analysis and motif enrichment using single-cell RNA and ATAC datasets. The scANANSE pipeline can be run from either R or Python. First, it exports data from standard single-cell objects. Next, it automatically runs multiple comparisons of cell cluster data. Finally, it imports the results back to the single-cell object, where the result can be further visualised, integrated, and interpreted. Here, we demonstrate our scANANSE pipeline on a publicly available PBMC multi-omics dataset. It identifies well-known cell type-specific hematopoietic factors. Importantly, we also demonstrated that scANANSE combined with GimmeMotifs is able to predict transcription factors with both activating and repressing roles in gene regulation.

Seq2science: an end-to-end workflow for functional genomics analysis.

Sequencing databases contain enormous amounts of functional genomics data, making them an extensive resource for genome-scale analysis. Reanalyzing publicly available data, and integrating it with new, project-specific data sets, can be invaluable. With current technologies, genomic experiments have become feasible for virtually any species of interest. However, using and integrating this data comes with its challenges, such as standardized and reproducible analysis. Seq2science is a multi-purpose workflow that covers preprocessing, quality control, visualization, and analysis of functional genomics sequencing data. It facilitates the downloading of sequencing data from all major databases, including NCBI SRA, EBI ENA, DDBJ, GSA, and ENCODE. Furthermore, it automates the retrieval of any genome assembly available from Ensembl, NCBI, and UCSC. It has been tested on a variety of species, and includes diverse workflows such as ATAC-, RNA-, and ChIP-seq. It consists of both generic as well as advanced steps, such as differential gene expression or peak accessibility analysis and differential motif analysis. Seq2science is built on the Snakemake workflow language and thus can be run on a range of computing infrastructures. It is available at https://github.com/vanheeringen-lab/seq2science.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease.

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Lead optimization of aryl hydrocarbon receptor ligands for treatment of inflammatory skin disorders.

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.

Enhanced pro-apoptosis gene signature following the activation of TAp63α in oocytes upon γ irradiation.

Specialized surveillance mechanisms are essential to maintain the genetic integrity of germ cells, which are not only the source of all somatic cells but also of the germ cells of the next generation. DNA damage and chromosomal aberrations are, therefore, not only detrimental for the individual but affect the entire species. In oocytes, the surveillance of the structural integrity of the DNA is maintained by the p53 family member TAp63α. The TAp63α protein is highly expressed in a closed and inactive state and gets activated to the open conformation upon the detection of DNA damage, in particular DNA double-strand breaks. To understand the cellular response to DNA damage that leads to the TAp63α triggered oocyte death we have investigated the RNA transcriptome of oocytes following irradiation at different time points. The analysis shows enhanced expression of pro-apoptotic and typical p53 target genes such as CDKn1a or Mdm2, concomitant with the activation of TAp63α. While DNA repair genes are not upregulated, inflammation-related genes become transcribed when apoptosis is initiated by activation of STAT transcription factors. Furthermore, comparison with the transcriptional profile of the ΔNp63α isoform from other studies shows only a minimal overlap, suggesting distinct regulatory programs of different p63 isoforms.

Proteomics of resistance to Notch1 inhibition in acute lymphoblastic leukemia reveals targetable kinase signatures.

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.

Characterization of In Vitro Differentiation of Human Primary Keratinocytes by RNA-Seq Analysis.

Human primary keratinocytes are often used as in vitro models for studies on epidermal differentiation and related diseases. Methods have been reported for in vitro differentiation of keratinocytes cultured in two-dimensional (2D) submerged manners using various induction conditions. Described here is a procedure for 2D in vitro keratinocyte differentiation method by contact inhibition and subsequent molecular characterization by RNA-seq. In brief, keratinocytes are grown in defined keratinocyte medium supplemented with growth factors until they are fully confluent. Differentiation is induced by close contacts between the keratinocytes and further stimulated by excluding growth factors in the medium. Using RNA-seq analyses, it is shown that both 1) differentiated keratinocytes exhibit distinct molecular signatures during differentiation and 2) the dynamic gene expression pattern largely resembles cells during epidermal stratification. As for comparison to normal keratinocyte differentiation, keratinocytes carrying mutations of the transcription factor p63 exhibit altered morphology and molecular signatures, consistent with their differentiation defects. In conclusion, this protocol details the steps for 2D in vitro keratinocyte differentiation and its molecular characterization, with an emphasis on bioinformatic analysis of RNA-seq data. Because RNA extraction and RNA-seq procedures have been well-documented, it is not the focus of this protocol. The experimental procedure of in vitro keratinocyte differentiation and bioinformatic analysis pipeline can be used to study molecular events during epidermal differentiation in healthy and diseased keratinocytes.

The acute myeloid leukemia associated AML1-ETO fusion protein alters the transcriptome and cellular progression in a single-oncogene expressing in vitro induced pluripotent stem cell based granulocyte differentiation model.

Acute myeloid leukemia (AML) is characterized by recurrent mutations that affect normal hematopoiesis. The analysis of human AMLs has mostly been performed using end-point materials, such as cell lines and patient derived AMLs that also carry additional contributing mutations. The molecular effects of a single oncogenic hit, such as expression of the AML associated oncoprotein AML1-ETO on hematopoietic development and transformation into a (pre-) leukemic state still needs further investigation. Here we describe the development and characterization of an induced pluripotent stem cell (iPSC) system that allows in vitro differentiation towards different mature myeloid cell types such as monocytes and granulocytes. During in vitro differentiation we expressed the AML1-ETO fusion protein and examined the effects of the oncoprotein on differentiation and the underlying alterations in the gene program at 8 different time points. Our analysis revealed that AML1-ETO as a single oncogenic hit in a non-mutated background blocks granulocytic differentiation, deregulates the gene program via altering the acetylome of the differentiating granulocytic cells, and induces t(8;21) AML associated leukemic characteristics. Together, these results reveal that inducible oncogene expression during in vitro differentiation of iPS cells provides a valuable platform for analysis of aberrant regulation in disease.