Signalling by senescent melanocytes hyperactivates hair growth.

Niche signals maintain stem cells in a prolonged quiescence or transiently activate them for proper regeneration1. Altering balanced niche signalling can lead to regenerative disorders. Melanocytic skin nevi in human often display excessive hair growth, suggesting hair stem cell hyperactivity. Here, using genetic mouse models of nevi2,3, we show that dermal clusters of senescent melanocytes drive epithelial hair stem cells to exit quiescence and change their transcriptome and composition, potently enhancing hair renewal. Nevus melanocytes activate a distinct secretome, enriched for signalling factors. Osteopontin, the leading nevus signalling factor, is both necessary and sufficient to induce hair growth. Injection of osteopontin or its genetic overexpression is sufficient to induce robust hair growth in mice, whereas germline and conditional deletions of either osteopontin or CD44, its cognate receptor on epithelial hair cells, rescue enhanced hair growth induced by dermal nevus melanocytes. Osteopontin is overexpressed in human hairy nevi, and it stimulates new growth of human hair follicles. Although broad accumulation of senescent cells, such as upon ageing or genotoxic stress, is detrimental for the regenerative capacity of tissue4, we show that signalling by senescent cell clusters can potently enhance the activity of adjacent intact stem cells and stimulate tissue renewal. This finding identifies senescent cells and their secretome as an attractive therapeutic target in regenerative disorders.

How and Why the Circadian Clock Regulates Proliferation of Adult Epithelial Stem Cells.

First described in the early 20th century, diurnal oscillations in stem cell proliferation exist in multiple internal epithelia, including in the gastrointestinal tract, and in the epidermis. In the mouse epidermis, 3- to 4-fold more stem cells are in S-phase during the night than during the day. More recent work showed that an intact circadian clock intrinsic to keratinocytes is required for these oscillations in epidermal stem cell proliferation. The circadian clock also regulates DNA excision repair and DNA damage in epidermal stem cells in response to ultraviolet B radiation. During skin inflammation, epidermal stem cell proliferation is increased and diurnal oscillations are suspended. Here we discuss possible reasons for the evolution of this stem cell phenomenon. We argue that the circadian clock coordinates intermediary metabolism and the cell cycle in epidermal stem cells to minimize the accumulation of DNA damage from metabolism-generated reactive oxygen species. Circadian disruption, common in modern society, leads to asynchrony between metabolism and the cell cycle, and we speculate this will lead to oxidative DNA damage, dysfunction of epidermal stem cells, and skin aging.

Circadian control of interferon-sensitive gene expression in murine skin.

The circadian clock coordinates a variety of immune responses with signals from the external environment to promote survival. We investigated the potential reciprocal relationship between the circadian clock and skin inflammation. We treated mice topically with the Toll-like receptor 7 (TLR7) agonist imiquimod (IMQ) to activate IFN-sensitive gene (ISG) pathways and induce psoriasiform inflammation. IMQ transiently altered core clock gene expression, an effect mirrored in human patient psoriatic lesions. In mouse skin 1 d after IMQ treatment, ISGs, including the key ISG transcription factor IFN regulatory factor 7 (Irf7), were more highly induced after treatment during the day than the night. Nuclear localization of phosphorylated-IRF7 was most prominently time-of-day dependent in epidermal leukocytes, suggesting that these cell types play an important role in the diurnal ISG response to IMQ. Mice lacking Bmal1 systemically had exacerbated and arrhythmic ISG/Irf7 expression after IMQ. Furthermore, daytime-restricted feeding, which affects the phase of the skin circadian clock, reverses the diurnal rhythm of IMQ-induced ISG expression in the skin. These results suggest a role for the circadian clock, driven by BMAL1, as a negative regulator of the ISG response, and highlight the finding that feeding time can modulate the skin immune response. Since the IFN response is essential for the antiviral and antitumor effects of TLR activation, these findings are consistent with the time-of-day-dependent variability in the ability to fight microbial pathogens and tumor initiation and offer support for the use of chronotherapy for their treatment.

Overexpression of Grainyhead-like 3 causes spina bifida and interacts genetically with mutant alleles of Grhl2 and Vangl2 in mice.

The genetic basis of human neural tube defects (NTDs), such as anencephaly and spina bifida (SB), is complex and heterogeneous. Grainyhead-like genes represent candidates for involvement in NTDs based on the presence of SB and exencephaly in mice carrying loss-of-function alleles of Grhl2 or Grhl3. We found that reinstatement of Grhl3 expression, by bacterial artificial chromosome (BAC)-mediated transgenesis, prevents SB in Grhl3-null embryos, as in the Grhl3 hypomorphic curly tail strain. Notably, however, further increase in expression of Grhl3 causes highly penetrant SB. Grhl3 overexpression recapitulates the spinal NTD phenotype of loss-of-function embryos, although the underlying mechanism differs. However, it does not phenocopy other defects of Grhl3-null embryos such as abnormal axial curvature, cranial NTDs (exencephaly) or skin barrier defects, the latter being rescued by the Grhl3-transgene. Grhl2 and Grhl3 can form homodimers and heterodimers, suggesting a possible model in which defects arising from overexpression of Grhl3 result from sequestration of Grhl2 in heterodimers, mimicking Grhl2 loss of function. This hypothesis predicts that increased abundance of Grhl2 would have an ameliorating effect in Grhl3 overexpressing embryo. Instead, we observed a striking additive genetic interaction between Grhl2 and Grhl3 gain-of-function alleles. Severe SB arose in embryos in which both genes were expressed at moderately elevated levels that individually do not cause NTDs. Furthermore, moderate Grhl3 overexpression also interacted with the Vangl2Lp allele to cause SB, demonstrating genetic interaction with the planar cell polarity signalling pathway that is implicated in mouse and human NTDs.

GRHL3 binding and enhancers rearrange as epidermal keratinocytes transition between functional states.

Transcription factor binding, chromatin modifications and large scale chromatin re-organization underlie progressive, irreversible cell lineage commitments and differentiation. We know little, however, about chromatin changes as cells enter transient, reversible states such as migration. Here we demonstrate that when human progenitor keratinocytes either differentiate or migrate they form complements of typical enhancers and super-enhancers that are unique for each state. Unique super-enhancers for each cellular state link to gene expression that confers functions associated with the respective cell state. These super-enhancers are also enriched for skin disease sequence variants. GRHL3, a transcription factor that promotes both differentiation and migration, binds preferentially to super-enhancers in differentiating keratinocytes, while during migration, it binds preferentially to promoters along with REST, repressing the expression of migration inhibitors. Key epidermal differentiation transcription factor genes, including GRHL3, are located within super-enhancers, and many of these transcription factors in turn bind to and regulate super-enhancers. Furthermore, GRHL3 represses the formation of a number of progenitor and non-keratinocyte super-enhancers in differentiating keratinocytes. Hence, chromatin relocates GRHL3 binding and enhancers to regulate both the irreversible commitment of progenitor keratinocytes to differentiation and their reversible transition to migration.

Neural tube closure depends on expression of Grainyhead-like 3 in multiple tissues.

Failure of neural tube closure leads to neural tube defects (NTDs), common congenital abnormalities in humans. Among the genes whose loss of function causes NTDs in mice, Grainyhead-like3 (Grhl3) is essential for spinal neural tube closure, with null mutants exhibiting fully penetrant spina bifida. During spinal neurulation Grhl3 is initially expressed in the surface (non-neural) ectoderm, subsequently in the neuroepithelial component of the neural folds and at the node-streak border, and finally in the hindgut endoderm. Here, we show that endoderm-specific knockout of Grhl3 causes late-arising spinal NTDs, preceded by increased ventral curvature of the caudal region which was shown previously to suppress closure of the spinal neural folds. This finding supports the hypothesis that diminished Grhl3 expression in the hindgut is the cause of spinal NTDs in the curly tail, carrying a hypomorphic Grhl3 allele. Complete loss of Grhl3 function produces a more severe phenotype in which closure fails earlier in neurulation, before the stage of onset of expression in the hindgut of wild-type embryos. This implicates additional tissues and NTD mechanisms in Grhl3 null embryos. Conditional knockout of Grhl3 in the neural plate and node-streak border has minimal effect on closure, suggesting that abnormal function of surface ectoderm, where Grhl3 transcripts are first detected, is primarily responsible for early failure of spinal neurulation in Grhl3 null embryos.

Characterization of enhancers and the role of the transcription factor KLF7 in regulating corneal epithelial differentiation.

During tissue development, transcription factors bind regulatory DNA regions called enhancers, often located at great distances from the genes they regulate, to control gene expression. The enhancer landscape during embryonic stem cell differentiation has been well characterized. By contrast, little is known about the shared and unique enhancer regulatory mechanisms in different ectodermally derived epithelial cells. Here we use ChIP sequencing (ChIP-seq) to identify domains enriched for the histone marks histone H3 lysine 4 trimethylation, histone H3 lysine 4 monomethylation, and histone H3 lysine 27 acetylation (H3K4me3, H3K4me1, and H3K27ac) and define, for the first time, the super enhancers and typical enhancers active in primary human corneal epithelial cells. We show that regulatory regions are often shared between cell types of the ectodermal lineage and that corneal epithelial super enhancers are already marked as potential regulatory domains in embryonic stem cells. Kruppel-like factor (KLF) motifs were enriched in corneal epithelial enhancers, consistent with the important roles of KLF4 and KLF5 in promoting corneal epithelial differentiation. We now show that the Kruppel family member KLF7 promotes the corneal progenitor cell state; on many genes, KLF7 antagonized the corneal differentiation-promoting KLF4. Furthermore, we found that two SNPs linked previously to corneal diseases, astigmatism, and Stevens-Johnson syndrome fall within corneal epithelial enhancers and alter their activity by disrupting transcription factor motifs that overlap these SNPs. Taken together, our work defines regulatory enhancers in corneal epithelial cells, highlights global gene-regulatory relationships shared among different epithelial cells, identifies a role for KLF7 as a KLF4 antagonist in corneal epithelial cell differentiation, and explains how two SNPs may contribute to corneal diseases.

Cofactors of LIM Domains Associate with Estrogen Receptor α to Regulate the Expression of Noncoding RNA H19 and Corneal Epithelial Progenitor Cell Function.

Cofactors of LIM domain proteins, CLIM1 and CLIM2, are widely expressed transcriptional cofactors that are recruited to gene regulatory regions by DNA-binding proteins, including LIM domain transcription factors. In the cornea, epithelium-specific expression of a dominant negative (DN) CLIM under the keratin 14 (K14) promoter causes blistering, wounding, inflammation, epithelial hyperplasia, and neovascularization followed by epithelial thinning and subsequent epidermal-like differentiation of the corneal epithelium. The defects in corneal epithelial differentiation and cell fate determination suggest that CLIM may regulate corneal progenitor cells and the transition to differentiation. Consistent with this notion, the K14-DN-Clim corneal epithelium first exhibits increased proliferation followed by fewer progenitor cells with decreased proliferative potential. In vivo ChIP-sequencing experiments with corneal epithelium show that CLIM binds to and regulates numerous genes involved in cell adhesion and proliferation, including limbally enriched genes. Intriguingly, CLIM associates primarily with non-LIM homeodomain motifs in corneal epithelial cells, including that of estrogen receptor α. Among CLIM targets is the noncoding RNA H19 whose deregulation is associated with Silver-Russell and Beckwith-Wiedemann syndromes. We demonstrate here that H19 negatively regulates corneal epithelial proliferation. In addition to cell cycle regulators, H19 affects the expression of multiple cell adhesion genes. CLIM interacts with estrogen receptor α at the H19 locus, potentially explaining the higher expression of H19 in female than male corneas. Together, our results demonstrate an important role for CLIM in regulating the proliferative potential of corneal epithelial progenitors and identify CLIM downstream target H19 as a regulator of corneal epithelial proliferation and adhesion.

Dominant mutations in GRHL3 cause Van der Woude Syndrome and disrupt oral periderm development.

Mutations in interferon regulatory factor 6 (IRF6) account for ∼70% of cases of Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate. In 8 of 45 VWS-affected families lacking a mutation in IRF6, we found coding mutations in grainyhead-like 3 (GRHL3). According to a zebrafish-based assay, the disease-associated GRHL3 mutations abrogated periderm development and were consistent with a dominant-negative effect, in contrast to haploinsufficiency seen in most VWS cases caused by IRF6 mutations. In mouse, all embryos lacking Grhl3 exhibited abnormal oral periderm and 17% developed a cleft palate. Analysis of the oral phenotype of double heterozygote (Irf6(+/-);Grhl3(+/-)) murine embryos failed to detect epistasis between the two genes, suggesting that they function in separate but convergent pathways during palatogenesis. Taken together, our data demonstrated that mutations in two genes, IRF6 and GRHL3, can lead to nearly identical phenotypes of orofacial cleft. They supported the hypotheses that both genes are essential for the presence of a functional oral periderm and that failure of this process contributes to VWS.

The co-factor of LIM domains (CLIM/LDB/NLI) maintains basal mammary epithelial stem cells and promotes breast tumorigenesis.

Mammary gland branching morphogenesis and ductal homeostasis relies on mammary stem cell function for the maintenance of basal and luminal cell compartments. The mechanisms of transcriptional regulation of the basal cell compartment are currently unknown. We explored these mechanisms in the basal cell compartment and identified the Co-factor of LIM domains (CLIM/LDB/NLI) as a transcriptional regulator that maintains these cells. Clims act within the basal cell compartment to promote branching morphogenesis by maintaining the number and proliferative potential of basal mammary epithelial stem cells. Clim2, in a complex with LMO4, supports mammary stem cells by directly targeting the Fgfr2 promoter in basal cells to increase its expression. Strikingly, Clims also coordinate basal-specific transcriptional programs to preserve luminal cell identity. These basal-derived cues inhibit epidermis-like differentiation of the luminal cell compartment and enhance the expression of luminal cell-specific oncogenes ErbB2 and ErbB3. Consistently, basal-expressed Clims promote the initiation and progression of breast cancer in the MMTV-PyMT tumor model, and the Clim-regulated branching morphogenesis gene network is a prognostic indicator of poor breast cancer outcome in humans.

GRHL3/GET1 and trithorax group members collaborate to activate the epidermal progenitor differentiation program.

The antagonistic actions of Polycomb and Trithorax are responsible for proper cell fate determination in mammalian tissues. In the epidermis, a self-renewing epithelium, previous work has shown that release from Polycomb repression only partially explains differentiation gene activation. We now show that Trithorax is also a key regulator of epidermal differentiation, not only through activation of genes repressed by Polycomb in progenitor cells, but also through activation of genes independent of regulation by Polycomb. The differentiation associated transcription factor GRHL3/GET1 recruits the ubiquitously expressed Trithorax complex to a subset of differentiation genes.

Brain and muscle Arnt-like protein-1 (BMAL1) controls circadian cell proliferation and susceptibility to UVB-induced DNA damage in the epidermis.

The role of the circadian clock in skin and the identity of genes participating in its chronobiology remain largely unknown, leading us to define the circadian transcriptome of mouse skin at two different stages of the hair cycle, telogen and anagen. The circadian transcriptomes of telogen and anagen skin are largely distinct, with the former dominated by genes involved in cell proliferation and metabolism. The expression of many metabolic genes is antiphasic to cell cycle-related genes, the former peaking during the day and the latter at night. Consistently, accumulation of reactive oxygen species, a byproduct of oxidative phosphorylation, and S-phase are antiphasic to each other in telogen skin. Furthermore, the circadian variation in S-phase is controlled by BMAL1 intrinsic to keratinocytes, because keratinocyte-specific deletion of Bmal1 obliterates time-of-day-dependent synchronicity of cell division in the epidermis leading to a constitutively elevated cell proliferation. In agreement with higher cellular susceptibility to UV-induced DNA damage during S-phase, we found that mice are most sensitive to UVB-induced DNA damage in the epidermis at night. Because in the human epidermis maximum numbers of keratinocytes go through S-phase in the late afternoon, we speculate that in humans the circadian clock imposes regulation of epidermal cell proliferation so that skin is at a particularly vulnerable stage during times of maximum UV exposure, thus contributing to the high incidence of human skin cancers.

The Ets transcription factor EHF as a regulator of cornea epithelial cell identity.

The cornea is the clear, outermost portion of the eye composed of three layers: an epithelium that provides a protective barrier while allowing transmission of light into the eye, a collagen-rich stroma, and an endothelium monolayer. How cornea development and aging is controlled is poorly understood. Here we characterize the mouse cornea transcriptome from early embryogenesis through aging and compare it with transcriptomes of other epithelial tissues, identifying cornea-enriched genes, pathways, and transcriptional regulators. Additionally, we profiled cornea epithelium and stroma, defining genes enriched in these layers. Over 10,000 genes are differentially regulated in the mouse cornea across the time course, showing dynamic expression during development and modest expression changes in fewer genes during aging. A striking transition time point for gene expression between postnatal days 14 and 28 corresponds with completion of cornea development at the transcriptional level. Clustering classifies co-expressed, and potentially co-regulated, genes into biologically informative categories, including groups that exhibit epithelial or stromal enriched expression. Based on these findings, and through loss of function studies and ChIP-seq, we show that the Ets transcription factor EHF promotes cornea epithelial fate through complementary gene activating and repressing activities. Furthermore, we identify potential interactions between EHF, KLF4, and KLF5 in promoting cornea epithelial differentiation. These data provide insights into the mechanisms underlying epithelial development and aging, identifying EHF as a regulator of cornea epithelial identity and pointing to interactions between Ets and KLF factors in promoting epithelial fate. Furthermore, this comprehensive gene expression data set for the cornea is a powerful tool for discovery of novel cornea regulators and pathways.

The effects of age and dysfunction on meibomian gland population dynamics.

PURPOSE: While meibomian gland dysfunction (MGD) is widely recognized as a major cause of evaporative dry eye disease, little is known about normal gland differentiation and lipid synthesis or the mechanism underlying gland atrophy and abnormal lipid secretion. The purpose of this study was to use single-cell and spatial transcriptomics to probe changes in cell composition, differentiation, and gene expression associated with two murine models of MGD: age-related gland atrophy in wild-type mice and altered meibum quality in acyl-CoA wax alcohol acyltransferase 2 (Awat2) knockout (KO) mice. METHODS: Young (6 month) and old (22 month) wild type, C57Bl/6 mice and young (3 month) and old (13 month) Awat2 KO mice were used in these studies. For single-cell analysis, the tarsal plate was dissected from the upper and lower eyelids, and single cells isolated and submitted to the UCI Genomic Core, while for the spatial analysis frozen tissue sections were shipped to Resolve Biosciences on dry ice and sections probed in duplicate using a meibomian gland specific, 100 gene Molecular Chartography panel. RESULTS: Analysis of gene expression patterns identified the stratified expression of lipogenic genes during meibocyte differentiation, which may control the progressive synthesis of meibum lipids; an age-related decrease in meibocytes; and increased immune cell infiltration. Additionally, we detected unique immune cell populations in the Awat2 KO mouse suggesting activation of psoriasis-like, inflammatory pathways perhaps caused by ductal dilation and hyperplasia. CONCLUSION: Together these findings support novel mechanism controlling gland function and dysfunction.

tauFisher predicts circadian time from a single sample of bulk and single-cell pseudobulk transcriptomic data.

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.

The estrogen-responsive Agr2 gene regulates mammary epithelial proliferation and facilitates lobuloalveolar development.

Agr2 is a putative protein disulfide isomerase (PDI) initially identified as an estrogen-responsive gene in breast cancer cell lines. While Agr2 expression in breast cancer is positively correlated with estrogen receptor (ER) expression, it is upregulated in both hormone dependent and independent carcinomas. Several in vitro and xenograft studies have implicated Agr2 in different oncogenic features of breast cancer; however, the physiological role of Agr2 in normal mammary gland development remains to be defined. Agr2 expression is developmentally regulated in the mammary gland, with maximum expression during late pregnancy and lactation. Using a mammary gland specific knockout mouse model, we show that Agr2 facilitates normal lobuloalveolar development by regulating mammary epithelial cell proliferation; we found no effects on apoptosis in Agr2(-/-) mammary epithelial cells. Consequently, mammary glands of Agr2(-/-) females exhibit reduced expression of milk proteins, and by two weeks post-partum their pups are smaller in size. Utilizing a conditional mouse model, we show that Agr2 constitutive expression drives precocious lobuloalveolar development and increased milk protein expression in the virgin mammary gland. In vitro studies using knock down and overexpression strategies in estrogen receptor positive and negative mammary epithelial cell lines demonstrate a role for Agr2 in estradiol-induced cell proliferation. In conclusion, the estrogen-responsive Agr2, a candidate breast cancer oncogene, regulates epithelial cell proliferation and lobuloalveolar development in the mammary gland. The pro-proliferative effects of Agr2 may explain its actions in early tumorigenesis.

The estrogen-regulated anterior gradient 2 (AGR2) protein in breast cancer: a potential drug target and biomarker.

Initially discovered as an estrogen-responsive gene in breast cancer cell lines, anterior gradient 2 (AGR2) is a developmentally regulated gene belonging to the protein disulfide isomerase (PDI) gene family. Developmentally, AGR2 is expressed in the mammary gland in an estrogen-dependent manner, and AGR2 knockout and overexpression mouse models indicate that the gene promotes lobuloalveolar development by stimulating cell proliferation. Although AGR2 overexpression alone seems insufficient for breast tumorigenesis in mice, several lines of investigations suggest that AGR2 promotes breast tumorigenesis. Overexpression of AGR2 in several breast cancer cell lines increases cell survival in clonogenic assays and cell proliferation, whereas AGR2 loss of function leads to decreased cell cycle progression and cell death. In addition, AGR2 was shown to promote metastasis of breast epithelial cells in an in vivo metastasis assay. As a PDI, AGR2 is thought to be involved in the unfolded protein response that alleviates endoplasmic reticulum stress. Since cancer has to overcome proteotoxic stress due to excess protein production, AGR2 may be one of many pro-survival factors recruited to assist in protein folding or degradation or both. When AGR2 is secreted, it plays a role in cellular adhesion and dissemination of metastatic tumor cells. In breast cancer, AGR2 expression is associated with estrogen receptor (ER)-positive tumors; its overexpression is a predictor of poor prognosis. The AGR2 gene is directly targeted by ER-alpha, which is preferentially bound in tumors with poor outcome. Whereas aromatase inhibitor therapy decreases AGR2 expression, tamoxifen acts as an agonist of AGR2 expression in ER-positive tumors, perhaps contributing to tamoxifen resistance. AGR2 is also overexpressed in a subset of ER-negative tumors. Furthermore, AGR2 expression is associated with the dissemination of metastatic breast cancer cells and can be used as a marker to identify circulating tumor cells and metastatic cells in sentinel lymph nodes. In conclusion, AGR2 is a promising drug target in breast cancer and may serve as a useful prognostic indicator as well as a marker of breast cancer metastasis.

The epidermal differentiation-associated Grainyhead gene Get1/Grhl3 also regulates urothelial differentiation.

Skin and bladder epithelia form effective permeability barriers through the activation of distinct differentiation gene programs. Using a genome-wide gene-expression study, we identified transcriptional regulators whose expression correlates highly with that of differentiation markers in both the bladder and skin, including the Grainyhead factor Get1/Grhl3, which is already known to be important for epidermal barrier formation. In the bladder, Get1 is most highly expressed in the differentiated umbrella cells and its mutation in mice leads to a defective bladder epithelial barrier formation due to the failure of apical membrane specialization. Genes encoding components of the specialized urothelial membrane, the uroplakins, were downregulated in Get1(-/-) mice. At least one of these genes, uroplakin II, is a direct target of Get1. The urothelial-specific activation of the uroplakin II gene is due to selective binding of Get1 to the uroplakin II promoter in urothelial cells, which is most likely regulated by histone modifications. These results show a crucial role for Get1 in urothelial differentiation and barrier formation.