Deciphering the heterogeneity of differentiating hPSC-derived corneal limbal stem cells through single-cell RNA sequencing.

A comprehensive understanding of the human pluripotent stem cell (hPSC) differentiation process stands as a prerequisite for the development of hPSC-based therapeutics. In this study, single-cell RNA sequencing (scRNA-seq) was performed to decipher the heterogeneity during differentiation of three hPSC lines toward corneal limbal stem cells (LSCs). The scRNA-seq data revealed nine clusters encompassing the entire differentiation process, among which five followed the anticipated differentiation path of LSCs. The remaining four clusters were previously undescribed cell states that were annotated as either mesodermal-like or undifferentiated subpopulations, and their prevalence was hPSC line dependent. Distinct cluster-specific marker genes identified in this study were confirmed by immunofluorescence analysis and employed to purify hPSC-derived LSCs, which effectively minimized the variation in the line-dependent differentiation efficiency. In summary, scRNA-seq offered molecular insights into the heterogeneity of hPSC-LSC differentiation, allowing a data-driven strategy for consistent and robust generation of LSCs, essential for future advancement toward clinical translation.

CLDN1 knock out keratinocytes as a model to investigate multiple skin disorders.

The transmembrane protein claudin-1 is critical for formation of the epidermal barrier structure called tight junctions (TJ) and has been shown to be important in multiple disease states. These include neonatal ichthyosis and sclerosing cholangitis syndrome, atopic dermatitis and various viral infections. To develop a model to investigate the role of claudin-1 in different disease settings, we used CRISPR/Cas9 to generate human immortalized keratinocyte (KC) lines lacking claudin-1 (CLDN1 KO). We then determined whether loss of claudin-1 expression affects epidermal barrier formation/function and KC differentiation/stratification. The absence of claudin-1 resulted in significantly reduced barrier function in both monolayer and organotypic cultures. CLDN1 KO cells demonstrated decreases in gene transcripts encoding the barrier protein filaggrin and the differentiation marker cytokeratin-10. Marked morphological differences were also observed in CLDN1 KO organotypic cultures including diminished stratification and reduced formation of the stratum granulosum. We also detected increased proliferative KC in the basale layer of CLDN1 KO organotypic cultures. These results further support the role of claudin-1 in epidermal barrier and suggest an additional role of this protein in appropriate stratification of the epidermis.

Electrical Impedance Spectroscopy Quantifies Skin Barrier Function in Organotypic In Vitro Epidermis Models.

Three-dimensional human epidermal equivalents (HEEs) are a state-of-the-art organotypic culture model in preclinical investigative dermatology and regulatory toxicology. In this study, we investigated the utility of electrical impedance spectroscopy (EIS) for noninvasive measurement of HEE epidermal barrier function. Our setup comprised a custom-made lid fit with 12 electrode pairs aligned on the standard 24-transwell cell culture system. Serial EIS measurements for 7 consecutive days did not impact epidermal morphology, and readouts showed comparable trends with HEEs measured only once. We determined 2 frequency ranges in the resulting impedance spectra: a lower frequency range termed EISdiff correlated with keratinocyte terminal differentiation independent of epidermal thickness and a higher frequency range termed EISSC correlated with stratum corneum thickness. HEEs generated from CRISPR/Cas9-engineered keratinocytes that lack key differentiation genes FLG, TFAP2A, AHR, or CLDN1 confirmed that keratinocyte terminal differentiation is the major parameter defining EISdiff. Exposure to proinflammatory psoriasis- or atopic dermatitis-associated cytokine cocktails lowered the expression of keratinocyte differentiation markers and reduced EISdiff. This cytokine-associated decrease in EISdiff was normalized after stimulation with therapeutic molecules. In conclusion, EIS provides a noninvasive system to consecutively and quantitatively assess HEE barrier function and to sensitively and objectively measure barrier development, defects, and repair.

The Aryl Hydrocarbon Receptor Regulates Epidermal Differentiation through Transient Activation of TFAP2A.

The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as an indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multiomics analyses using human skin keratinocytes revealed that upon ligand activation, AHR binds open chromatin to induce expression of transcription factors, for example, TFAP2A, as a swift response to environmental stimuli. The terminal differentiation program, including upregulation of barrier genes, FLG and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides additional insights into the molecular mechanism behind AHR-mediated barrier function and identifies potential targets for the treatment of skin barrier diseases.

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.

Novel methodologies for host-microbe interactions and microbiome-targeted therapeutics in 3D organotypic skin models.

BACKGROUND: Following descriptive studies on skin microbiota in health and disease, mechanistic studies on the interplay between skin and microbes are on the rise, for which experimental models are in great demand. Here, we present a novel methodology for microbial colonization of organotypic skin and analysis thereof. RESULTS: An inoculation device ensured a standardized application area on the stratum corneum and a homogenous distribution of bacteria, while preventing infection of the basolateral culture medium even during prolonged culture periods for up to 2 weeks at a specific culture temperature and humidity. Hereby, host-microbe interactions and antibiotic interventions could be studied, revealing diverse host responses to various skin-related bacteria and pathogens. CONCLUSIONS: Our methodology is easily transferable to a wide variety of organotypic skin or mucosal models and different microbes at every cell culture facility at low costs. We envision that this study will kick-start skin microbiome studies using human organotypic skin cultures, providing a powerful alternative to experimental animal models in pre-clinical research. Video Abstract.

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.

The aryl hydrocarbon receptor regulates epidermal differentiation through transient activation of TFAP2A.

The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multi-omics analyses using human skin keratinocytes revealed that, upon ligand activation, AHR binds open chromatin to induce expression of transcription factors (TFs), e.g., Transcription Factor AP-2α (TFAP2A), as a swift response to environmental stimuli. The terminal differentiation program including upregulation of barrier genes, filaggrin and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides novel insights into the molecular mechanism behind AHR-mediated barrier function and potential novel targets for the treatment of skin barrier diseases.

Investigations into the FLG Null Phenotype: Showcasing the Methodology for CRISPR/Cas9 Editing of Human Keratinocytes.

Ever since the association between FLG loss-of-function variants and ichthyosis vulgaris and atopic dermatitis disease onset was identified, FLGs function has been under investigation. Intraindividual genomic predisposition, immunological confounders, and environmental interactions complicate the comparison between FLG genotypes and related causal effects. Using CRISPR/Cas9, we generated human FLG-knockout (ΔFLG) N/TERT-2G keratinocytes. FLG deficiency was shown by immunohistochemistry of human epidermal equivalent cultures. Next to (partial) loss of structural proteins (involucrin, hornerin, keratin 2, and transglutaminase 1), the stratum corneum was denser and lacked the typical basket weave appearance. In addition, electrical impedance spectroscopy and transepidermal water loss analyses highlighted a compromised epidermal barrier in ΔFLG human epidermal equivalents. Correction of FLG reinstated the presence of keratohyalin granules in the stratum granulosum, FLG protein expression, and expression of the proteins mentioned earlier. The beneficial effects on stratum corneum formation were reflected by the normalization of electrical impedance spectroscopy and transepidermal water loss. This study shows the causal phenotypical and functional consequences of FLG deficiency, indicating that FLG is not only central in epidermal barrier function but also vital for epidermal differentiation by orchestrating the expression of other important epidermal proteins. These observations pave the way to fundamental investigations into the exact role of FLG in skin biology and disease.

Targeting Skin Barrier Function in Atopic Dermatitis.

Atopic dermatitis (AD) is the most common chronic inflammatory skin disease in the general population. Skin barrier dysfunction is the central abnormality leading to AD. The cause of skin barrier dysfunction is complex and rooted in genetic mutations, interactions between the immune pathway activation and epithelial cells, altered host defense mechanisms, as well as environmental influences that cause epithelial cell activation and release of alarmins (such as thymic stromal lymphopoietin) that can activate the type 2 immune pathway, including generation of interleukins 4 and 13, which induces defects in the skin barrier and increased allergic inflammation. These inflammatory pathways are further influenced by environmental factors including the microbiome (especially Staphylococcus aureus), air pollution, stress, and other factors. As such, AD is a syndrome involving multiple phenotypes, all of which have in common skin barrier dysfunction as a key contributing factor. Understanding mechanisms leading to skin barrier dysfunction in AD is pointing to the development of new topical and systemic treatments in AD that helps keep skin borders secure and effectively treat the disease.

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.

CRISPR-Cas9‒Based Genomic Engineering in Keratinocytes: From Technology to Application.

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.

Carboxamide Derivatives Are Potential Therapeutic AHR Ligands for Restoring IL-4 Mediated Repression of Epidermal Differentiation Proteins.

Atopic dermatitis (AD) is a common T-helper 2 (Th2) lymphocyte-mediated chronic inflammatory skin disease characterized by disturbed epidermal differentiation (e.g., filaggrin (FLG) expression) and diminished skin barrier function. Therapeutics targeting the aryl hydrocarbon receptor (AHR), such as coal tar and tapinarof, are effective in AD, yet new receptor ligands with improved potency or bioavailability are in demand to expand the AHR-targeting therapeutic arsenal. We found that carboxamide derivatives from laquinimod, tasquinimod, and roquinimex can activate AHR signaling at low nanomolar concentrations. Tasquinimod derivative (IMA-06504) and its prodrug (IMA-07101) provided full agonist activity and were most effective to induce FLG and other epidermal differentiation proteins, and counteracted IL-4 mediated repression of terminal differentiation. Partial agonist activity by other derivatives was less efficacious. The previously reported beneficial safety profile of these novel small molecules, and the herein reported therapeutic potential of specific carboxamide derivatives, provides a solid rationale for further preclinical assertation.

Antimicrobial Late Cornified Envelope Proteins: The Psoriasis Risk Factor Deletion of LCE3B/C Genes Affects Microbiota Composition.

Late cornified envelope proteins are predominantly expressed in the skin and other cornified epithelia. On the basis of sequence similarity, this 18-member homologous gene family has been subdivided into six groups. The LCE3 proteins have been the focus of dermatological research because the combined deletion of LCE3B and LCE3C genes (LCE3B/C-del) is a risk factor for psoriasis. We previously reported that LCE3B/C-del increases the expression of the LCE3A gene and that LCE3 proteins exert antibacterial activity. In this study, we analyzed the antimicrobial properties of other family members and the role of LCE3B/C-del in the modulation of microbiota composition of the skin and oral cavity. Differences in killing efficiency and specificity between the late cornified envelope proteins and their target microbes were found, and the amino acid content rather than the order of the well-conserved central domain of the LCE3A protein was found responsible for its antibacterial activity. In vivo, LCE3B/C-del correlated with a higher beta-diversity in the skin and oral microbiota. From these results, we conclude that all late cornified envelope proteins possess antimicrobial activity. Tissue-specific and genotype-dependent antimicrobial protein profiles impact skin and oral microbiota composition, which could direct toward LCE3B/C-del‒associated dysbiosis and a possible role for microbiota in the pathophysiology of psoriasis.

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.

Research Techniques Made Simple: Delivery of the CRISPR/Cas9 Components into Epidermal Cells.

CRISPR/Cas9 technology is a powerful tool used to alter the genetic landscape of various hosts. This has been exemplified by its success in the transgenic animal world where it has been utilized to develop novel mouse lines modeling numerous disease states. The technology has helped to develop both in vitro and in vivo systems that simulate diseases within the fields of epithelial biology, skin cancer biology, dermatology, and beyond. Importantly, the delivery of the single-guide RNA/Cas9 editing complex to the host cell is key for its success. In this paper, we discuss the various methods that have been utilized as delivery techniques for CRISPR/Cas9 components, the benefits and pitfalls of each, and how successful they have been at genetically modifying epidermal cells. In addition, we acknowledge recent advances in the field of dermatology that have harnessed these methods to better understand epidermal biology, identify potential therapeutic targets, or serve as novel methods to treat disease states.

Terminal keratinocyte differentiation in vitro is associated with a stable DNA methylome.

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.