Ewastools: Infinium Human Methylation BeadChip pipeline for population epigenetics integrated into Galaxy.

BACKGROUND: Infinium Human Methylation BeadChip is an array platform for complex evaluation of DNA methylation at an individual CpG locus in the human genome based on Illumina's bead technology and is one of the most common techniques used in epigenome-wide association studies. Finding associations between epigenetic variation and phenotype is a significant challenge in biomedical research. The newest version, HumanMethylationEPIC, quantifies the DNA methylation level of 850,000 CpG sites, while the previous versions, HumanMethylation450 and HumanMethylation27, measured >450,000 and 27,000 loci, respectively. Although a number of bioinformatics tools have been developed to analyse this assay, they require some programming skills and experience in order to be usable. RESULTS: We have developed a pipeline for the Galaxy platform for those without experience aimed at DNA methylation analysis using the Infinium Human Methylation BeadChip. Our tool is integrated into Galaxy (http://galaxyproject.org), a web-based platform. This allows users to analyse data from the Infinium Human Methylation BeadChip in the easiest possible way. CONCLUSIONS: The pipeline provides a group of integrated analytical methods wrapped into an easy-to-use interface. Our tool is available from the Galaxy ToolShed, GitHub repository, and also as a Docker image. The aim of this project is to make Infinium Human Methylation BeadChip analysis more flexible and accessible to everyone.

The ELIXIR Human Copy Number Variations Community: building bioinformatics infrastructure for research.

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While "High-Throughput" sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR's recently established human CNV Community, with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context.

Cellular sheddases are induced by Merkel cell polyomavirus small tumour antigen to mediate cell dissociation and invasiveness.

Merkel cell carcinoma (MCC) is an aggressive skin cancer with a high propensity for recurrence and metastasis. Merkel cell polyomavirus (MCPyV) is recognised as the causative factor in the majority of MCC cases. The MCPyV small tumour antigen (ST) is considered to be the main viral transforming factor, however potential mechanisms linking ST expression to the highly metastatic nature of MCC are yet to be fully elucidated. Metastasis is a complex process, with several discrete steps required for the formation of secondary tumour sites. One essential trait that underpins the ability of cancer cells to metastasise is how they interact with adjoining tumour cells and the surrounding extracellular matrix. Here we demonstrate that MCPyV ST expression disrupts the integrity of cell-cell junctions, thereby enhancing cell dissociation and implicate the cellular sheddases, A disintegrin and metalloproteinase (ADAM) 10 and 17 proteins in this process. Inhibition of ADAM 10 and 17 activity reduced MCPyV ST-induced cell dissociation and motility, attributing their function as critical to the MCPyV-induced metastatic processes. Consistent with these data, we confirm that ADAM 10 and 17 are upregulated in MCPyV-positive primary MCC tumours. These novel findings implicate cellular sheddases as key host cell factors contributing to virus-mediated cellular transformation and metastasis. Notably, ADAM protein expression may be a novel biomarker of MCC prognosis and given the current interest in cellular sheddase inhibitors for cancer therapeutics, it highlights ADAM 10 and 17 activity as a novel opportunity for targeted interventions for disseminated MCC.

MCV-miR-M1 Targets the Host-Cell Immune Response Resulting in the Attenuation of Neutrophil Chemotaxis.

Virus-encoded microRNAs are emerging as key regulators of persistent infection and host-cell immune evasion. Merkel cell polyomavirus, the predominant etiological agent of Merkel cell carcinoma, encodes a single microRNA, MCV-miR-M1, which targets the oncogenic Merkel cell polyomavirus large T antigen. MCV-miR-M1 has previously been shown to play an important role in the establishment of long-term infection, however, the underlying mechanism is not fully understood. A key unanswered question is whether, in addition to autoregulating large T antigen, MCV-miR-M1 also targets cellular transcripts to orchestrate an environment conducive to persistent infection. To address this, we adopted an RNA sequencing-based approach to identify cellular targets of MCV-miR-M1. Intriguingly, bioinformatics analysis of transcripts that are differentially expressed in cells expressing MCV-miR-M1 revealed several genes implicated in immune evasion. Subsequent target validation led to the identification of the innate immunity protein, SP100, as a direct target of MCV-miR-M1. Moreover, MCV-miR-M1-mediated modulation of SP100 was associated with a significant decrease in CXCL8 secretion, resulting in the attenuation of neutrophil chemotaxis toward Merkel cells harboring synthetic Merkel cell polyomavirus. Based on these observations, we propose that MCV-miR-M1 targets key immune response regulators to help facilitate persistent infection, which is a prerequisite for cellular transformation in Merkel cell carcinoma.

The cellular chloride channels CLIC1 and CLIC4 contribute to virus-mediated cell motility.

Ion channels regulate many aspects of cell physiology, including cell proliferation, motility, and migration, and aberrant expression and activity of ion channels is associated with various stages of tumor development, with K+ and Cl- channels now being considered the most active during tumorigenesis. Accordingly, emerging in vitro and preclinical studies have revealed that pharmacological manipulation of ion channel activity offers protection against several cancers. Merkel cell polyomavirus (MCPyV) is a major cause of Merkel cell carcinoma (MCC), primarily because of the expression of two early regulatory proteins termed small and large tumor antigens (ST and LT, respectively). Several molecular mechanisms have been attributed to MCPyV-mediated cancer formation but, thus far, no studies have investigated any potential link to cellular ion channels. Here we demonstrate that Cl- channel modulation can reduce MCPyV ST-induced cell motility and invasiveness. Proteomic analysis revealed that MCPyV ST up-regulates two Cl- channels, CLIC1 and CLIC4, which when silenced, inhibit MCPyV ST-induced motility and invasiveness, implicating their function as critical to MCPyV-induced metastatic processes. Consistent with these data, we confirmed that CLIC1 and CLIC4 are up-regulated in primary MCPyV-positive MCC patient samples. We therefore, for the first time, implicate cellular ion channels as a key host cell factor contributing to virus-mediated cellular transformation. Given the intense interest in ion channel modulating drugs for human disease. This highlights CLIC1 and CLIC4 activity as potential targets for MCPyV-induced MCC.

Inhibition of HOX/PBX dimer formation leads to necroptosis in acute myeloid leukemia cells.

The HOX genes encode a family of transcription factors that have key roles in both development and malignancy. Disrupting the interaction between HOX proteins and their binding partner, PBX, has been shown to cause apoptotic cell death in a range of solid tumors. However, despite HOX proteins playing a particularly significant role in acute myeloid leukemia (AML), the relationship between HOX gene expression and patient survival has not been evaluated (with the exception of HOXA9), and the mechanism by which HOX/PBX inhibition induces cell death in this malignancy is not well understood. In this study, we show that the expression of HOXA5, HOXB2, HOXB4, HOXB9, and HOXC9, but not HOXA9, in primary AML samples is significantly related to survival. Furthermore, the previously described inhibitor of HOX/PBX dimerization, HXR9, is cytotoxic to both AML-derived cell lines and primary AML cells from patients. The mechanism of cell death is not dependent on apoptosis but instead involves a regulated form of necrosis referred to as necroptosis. HXR9-induced necroptosis is enhanced by inhibitors of protein kinase C (PKC) signaling, and HXR9 combined with the PKC inhibitor Ro31 causes a significantly greater reduction in tumor growth compared to either reagent alone.

5C analysis of the Epidermal Differentiation Complex locus reveals distinct chromatin interaction networks between gene-rich and gene-poor TADs in skin epithelial cells.

Mammalian genomes contain several dozens of large (>0.5 Mbp) lineage-specific gene loci harbouring functionally related genes. However, spatial chromatin folding, organization of the enhancer-promoter networks and their relevance to Topologically Associating Domains (TADs) in these loci remain poorly understood. TADs are principle units of the genome folding and represents the DNA regions within which DNA interacts more frequently and less frequently across the TAD boundary. Here, we used Chromatin Conformation Capture Carbon Copy (5C) technology to characterize spatial chromatin interaction network in the 3.1 Mb Epidermal Differentiation Complex (EDC) locus harbouring 61 functionally related genes that show lineage-specific activation during terminal keratinocyte differentiation in the epidermis. 5C data validated by 3D-FISH demonstrate that the EDC locus is organized into several TADs showing distinct lineage-specific chromatin interaction networks based on their transcription activity and the gene-rich or gene-poor status. Correlation of the 5C results with genome-wide studies for enhancer-specific histone modifications (H3K4me1 and H3K27ac) revealed that the majority of spatial chromatin interactions that involves the gene-rich TADs at the EDC locus in keratinocytes include both intra- and inter-TAD interaction networks, connecting gene promoters and enhancers. Compared to thymocytes in which the EDC locus is mostly transcriptionally inactive, these interactions were found to be keratinocyte-specific. In keratinocytes, the promoter-enhancer anchoring regions in the gene-rich transcriptionally active TADs are enriched for the binding of chromatin architectural proteins CTCF, Rad21 and chromatin remodeler Brg1. In contrast to gene-rich TADs, gene-poor TADs show preferential spatial contacts with each other, do not contain active enhancers and show decreased binding of CTCF, Rad21 and Brg1 in keratinocytes. Thus, spatial interactions between gene promoters and enhancers at the multi-TAD EDC locus in skin epithelial cells are cell type-specific and involve extensive contacts within TADs as well as between different gene-rich TADs, forming the framework for lineage-specific transcription.

p63 Transcription Factor Regulates Nuclear Shape and Expression of Nuclear Envelope-Associated Genes in Epidermal Keratinocytes.

The maintenance of a proper nuclear architecture and three-dimensional organization of the genes, enhancer elements, and transcription machinery plays an essential role in tissue development and regeneration. Here we show that in the developing skin, epidermal progenitor cells of mice lacking p63 transcription factor display alterations in the nuclear shape accompanied by a marked decrease in expression of several nuclear envelope-associated components (Lamin B1, Lamin A/C, Sun1, Nesprin-3, Plectin) compared with controls. Furthermore, chromatin immunoprecipitation-quantitative PCR assay showed enrichment of p63 on Sun1, Syne3, and Plec promoters, suggesting them as p63 targets. Alterations in the nuclei shape and expression of nuclear envelope-associated proteins were accompanied by altered distribution patterns of the repressive histone marks trimethylation on lysine 27 of histone H3, trimethylation on lysine 9 of histone H3, and heterochromatin protein 1-alpha in p63-null keratinocytes. These changes were also accompanied by downregulation of the transcriptional activity and relocation of the keratinocyte-specific gene loci away from the sites of active transcription toward the heterochromatin-enriched repressive nuclear compartments in p63-null cells. These data demonstrate functional links between the nuclear envelope organization, chromatin architecture, and gene expression in keratinocytes and suggest nuclear envelope-associated genes as important targets mediating p63-regulated gene expression program in the epidermis.

Cbx4 maintains the epithelial lineage identity and cell proliferation in the developing stratified epithelium.

During development, multipotent progenitor cells establish lineage-specific programmers of gene activation and silencing underlying their differentiation into specialized cell types. We show that the Polycomb component Cbx4 serves as a critical determinant that maintains the epithelial identity in the developing epidermis by repressing nonepidermal gene expression programs. Cbx4 ablation in mice results in a marked decrease of the epidermal thickness and keratinocyte (KC) proliferation associated with activation of numerous neuronal genes and genes encoding cyclin-dependent kinase inhibitors (p16/p19 and p57). Furthermore, the chromodomain- and SUMO E3 ligase-dependent Cbx4 activities differentially regulate proliferation, differentiation, and expression of nonepidermal genes in KCs. Finally, Cbx4 expression in KCs is directly regulated by p63 transcription factor, whereas Cbx4 overexpression is capable of partially rescuing the effects of p63 ablation on epidermal development. These data demonstrate that Cbx4 plays a crucial role in the p63-regulated program of epidermal differentiation, maintaining the epithelial identity and proliferative activity in KCs via repression of the selected nonepidermal lineage and cell cycle inhibitor genes.

MicroRNA-214 controls skin and hair follicle development by modulating the activity of the Wnt pathway.

Skin development is governed by complex programs of gene activation and silencing, including microRNA-dependent modulation of gene expression. Here, we show that miR-214 regulates skin morphogenesis and hair follicle (HF) cycling by targeting ß-catenin, a key component of the Wnt signaling pathway. miR-214 exhibits differential expression patterns in the skin epithelium, and its inducible overexpression in keratinocytes inhibited proliferation, which resulted in formation of fewer HFs with decreased hair bulb size and thinner hair production. The inhibitory effects of miR-214 on HF development and cycling were associated with altered activities of multiple signaling pathways, including decreased expression of key Wnt signaling mediators ß-catenin and Lef-1, and were rescued by treatment with pharmacological Wnt activators. Finally, we identify ß-catenin as one of the conserved miR-214 targets in keratinocytes. These data provide an important foundation for further analyses of miR-214 as a key regulator of Wnt pathway activity and stem cell functions during normal tissue homeostasis, regeneration, and aging.

Complex changes in the apoptotic and cell differentiation programs during initiation of the hair follicle response to chemotherapy.

Chemotherapy has severe side effects in normal rapidly proliferating organs, such as hair follicles, and causes massive apoptosis in hair matrix keratinocytes followed by hair loss. To define the molecular signature of hair follicle response to chemotherapy, human scalp hair follicles cultured ex vivo were treated with doxorubicin (DXR), and global microarray analysis was performed 3 hours after treatment. Microarray data revealed changes in expression of 504 genes in DXR-treated hair follicles versus controls. Among these genes, upregulations of several tumor necrosis factor family of apoptotic receptors (FAS, TRAIL (tumor necrosis factor-related apoptosis-inducing ligand) receptors 1/2), as well as of a large number of keratin-associated protein genes, were seen after DXR treatment. Hair follicle apoptosis induced by DXR was significantly inhibited by either TRAIL-neutralizing antibody or caspase-8 inhibitor, thus suggesting a previously unreported role for TRAIL receptor signaling in mediating DXR-induced hair loss. These data demonstrate that the early phase of the hair follicle response to DXR includes upregulation of apoptosis-associated markers, as well as substantial reorganization of the terminal differentiation programs in hair follicle keratinocytes. These data provide an important platform for further studies toward the design of effective approaches for the management of chemotherapy-induced hair loss.

Bone morphogenetic protein signaling suppresses wound-induced skin repair by inhibiting keratinocyte proliferation and migration.

Bone morphogenetic protein (BMP) signaling plays a key role in the control of skin development and postnatal remodeling by regulating keratinocyte proliferation, differentiation, and apoptosis. To study the role of BMPs in wound-induced epidermal repair, we used transgenic mice overexpressing the BMP downstream component Smad1 under the control of a K14 promoter as an in vivo model, as well as ex vivo and in vitro assays. K14-caSmad1 (transgenic mice overexpressing a constitutively active form of Smad1 under K14 promoter) mice exhibited retarded wound healing associated with significant inhibition of proliferation and increased apoptosis in healing wound epithelium. Furthermore, microarray and quantitative real-time reverse-transcriptase-PCR (qRT-PCR) analyses revealed decreased expression of a number of cytoskeletal/cell motility-associated genes including wound-associated keratins (Krt16, Krt17) and Myosin VA (Myo5a), in the epidermis of K14-caSmad1 mice versus wild-type (WT) controls during wound healing. BMP treatment significantly inhibited keratinocyte migration ex vivo, and primary keratinocytes of K14-caSmad1 mice showed retarded migration compared with WT controls. Finally, small interfering RNA (siRNA)-mediated silencing of BMPR-1B in primary mouse keratinocytes accelerated cell migration and was associated with increased expression of Krt16, Krt17, and Myo5a compared with controls. Thus, this study demonstrates that BMPs inhibit keratinocyte proliferation, cytoskeletal organization, and migration in regenerating skin epithelium during wound healing, and raises a possibility for using BMP antagonists for the management of chronic wounds.

p63 and Brg1 control developmentally regulated higher-order chromatin remodelling at the epidermal differentiation complex locus in epidermal progenitor cells.

Chromatin structural states and their remodelling, including higher-order chromatin folding and three-dimensional (3D) genome organisation, play an important role in the control of gene expression. The role of 3D genome organisation in the control and execution of lineage-specific transcription programmes during the development and differentiation of multipotent stem cells into specialised cell types remains poorly understood. Here, we show that substantial remodelling of the higher-order chromatin structure of the epidermal differentiation complex (EDC), a keratinocyte lineage-specific gene locus on mouse chromosome 3, occurs during epidermal morphogenesis. During epidermal development, the locus relocates away from the nuclear periphery towards the nuclear interior into a compartment enriched in SC35-positive nuclear speckles. Relocation of the EDC locus occurs prior to the full activation of EDC genes involved in controlling terminal keratinocyte differentiation and is a lineage-specific, developmentally regulated event controlled by transcription factor p63, a master regulator of epidermal development. We also show that, in epidermal progenitor cells, p63 directly regulates the expression of the ATP-dependent chromatin remodeller Brg1, which binds to distinct domains within the EDC and is required for relocation of the EDC towards the nuclear interior. Furthermore, Brg1 also regulates gene expression within the EDC locus during epidermal morphogenesis. Thus, p63 and its direct target Brg1 play an essential role in remodelling the higher-order chromatin structure of the EDC and in the specific positioning of this locus within the landscape of the 3D nuclear space, as required for the efficient expression of EDC genes in epidermal progenitor cells during skin development.

Remodeling of three-dimensional organization of the nucleus during terminal keratinocyte differentiation in the epidermis.

The nucleus of epidermal keratinocytes (KCs) is a complex and highly compartmentalized organelle, whose structure is markedly changed during terminal differentiation and transition of the genome from a transcriptionally active state seen in the basal and spinous epidermal cells to a fully inactive state in the keratinized cells of the cornified layer. Here, using multicolor confocal microscopy, followed by computational image analysis and mathematical modeling, we demonstrate that in normal mouse footpad epidermis, transition of KCs from basal epidermal layer to the granular layer is accompanied by marked differences in nuclear architecture and microenvironment including the following: (i) decrease in the nuclear volume; (ii) decrease in expression of the markers of transcriptionally active chromatin; (iii) internalization and decrease in the number of nucleoli; (iv) increase in the number of pericentromeric heterochromatic clusters; and (v) increase in the frequency of associations between the pericentromeric clusters, chromosomal territory 3, and nucleoli. These data suggest a role for nucleoli and pericentromeric heterochromatin clusters as organizers of nuclear microenvironment required for proper execution of gene expression programs in differentiating KCs, and provide important background information for further analyses of alterations in the topological genome organization seen in pathological skin conditions, including disorders of epidermal differentiation and epidermal tumors.

Genome organizing function of SATB1 in tumor progression.

When cells change functions or activities (such as during differentiation, response to extracellular stimuli, or migration), gene expression undergoes large-scale reprogramming, in cell type- and function-specific manners. Large changes in gene regulation require changes in chromatin architecture, which involve recruitment of chromatin remodeling enzymes and epigenomic modification enzymes to specific genomic loci. Transcription factors must also be accurately assembled at these loci. SATB1 is a genome organizer protein that facilitates these processes, providing a nuclear architectural platform that anchors hundreds of genes, through its interaction with specific genomic sequences; this activity allows expression of all these genes to be regulated in parallel, and enables cells to thereby alter their function. We review and describe future perspectives on SATB1 function in higher-order chromatin structure and gene regulation, and its role in metastasis of breast cancer and other tumor types.

Cigarette smoke-induced transgenerational alterations in genome stability in cord blood of human F1 offspring.

The relevance of preconceptional and prenatal toxicant exposures for genomic stability in offspring is difficult to analyze in human populations, because gestational exposures usually cannot be separated from preconceptional exposures. To analyze the roles of exposures during gestation and conception on genomic stability in the offspring, stability was assessed via the Comet assay and highly sensitive, semiautomated confocal laser scans of γH2AX foci in cord, maternal, and paternal blood as well as spermatozoa from 39 families in Crete, Greece, and the United Kingdom. With use of multivariate linear regression analysis with backward selection, preconceptional paternal smoking (% tail DNA: P>0.032; γH2AX foci: P>0.018) and gestational maternal (% tail DNA: P>0.033) smoking were found to statistically significantly predict DNA damage in the cord blood of F1 offspring. Maternal passive smoke exposure was not identified as a predictor of DNA damage in cord blood, indicating that the effect of paternal smoking may be transmitted via the spermatozoal genome. Taken together, these studies reveal a role for cigarette smoke in the induction of DNA alterations in human F1 offspring via exposures of the fetus in utero or the paternal germline. Moreover, the identification of transgenerational DNA alterations in the unexposed F1 offspring of smoking-exposed fathers supports the claim that cigarette smoke is a human germ cell mutagen.

Lhx2 differentially regulates Sox9, Tcf4 and Lgr5 in hair follicle stem cells to promote epidermal regeneration after injury.

The Lhx2 transcription factor plays essential roles in morphogenesis and patterning of ectodermal derivatives as well as in controlling stem cell activity. Here, we show that during murine skin morphogenesis, Lhx2 is expressed in the hair follicle (HF) buds, whereas in postnatal telogen HFs Lhx2(+) cells reside in the stem cell-enriched epithelial compartments (bulge, secondary hair germ) and co-express selected stem cell markers (Sox9, Tcf4 and Lgr5). Remarkably, Lhx2(+) cells represent the vast majority of cells in the bulge and secondary hair germ that proliferate in response to skin injury. This is functionally important, as wound re-epithelization is significantly retarded in heterozygous Lhx2 knockout (+/-) mice, whereas anagen onset in the HFs located closely to the wound is accelerated compared with wild-type mice. Cell proliferation in the bulge and the number of Sox9(+) and Tcf4(+) cells in the HFs closely adjacent to the wound in Lhx2(+/-) mice are decreased in comparison with wild-type controls, whereas expression of Lgr5 and cell proliferation in the secondary hair germ are increased. Furthermore, acceleration of wound-induced anagen development in Lhx2(+/-) mice is inhibited by administration of Lgr5 siRNA. Finally, Chip-on-chip/ChIP-qPCR and reporter assay analyses identified Sox9, Tcf4 and Lgr5 as direct Lhx2 targets in keratinocytes. These data strongly suggest that Lhx2 positively regulates Sox9 and Tcf4 in the bulge cells, and promotes wound re-epithelization, whereas it simultaneously negatively regulates Lgr5 in the secondary hair germ and inhibits HF cycling. Thus, Lhx2 operates as an important regulator of epithelial stem cell activity in the skin response to injury.

p63 regulates Satb1 to control tissue-specific chromatin remodeling during development of the epidermis.

During development, multipotent progenitor cells establish tissue-specific programs of gene expression. In this paper, we show that p63 transcription factor, a master regulator of epidermal morphogenesis, executes its function in part by directly regulating expression of the genome organizer Satb1 in progenitor cells. p63 binds to a proximal regulatory region of the Satb1 gene, and p63 ablation results in marked reduction in the Satb1 expression levels in the epidermis. Satb1(-/-) mice show impaired epidermal morphology. In Satb1-null epidermis, chromatin architecture of the epidermal differentiation complex locus containing genes associated with epidermal differentiation is altered primarily at its central domain, where Satb1 binding was confirmed by chromatin immunoprecipitation-on-chip analysis. Furthermore, genes within this domain fail to be properly activated upon terminal differentiation. Satb1 expression in p63(+/-) skin explants treated with p63 small interfering ribonucleic acid partially restored the epidermal phenotype of p63-deficient mice. These data provide a novel mechanism by which Satb1, a direct downstream target of p63, contributes in epidermal morphogenesis via establishing tissue-specific chromatin organization and gene expression in epidermal progenitor cells.

BMP signaling induces cell-type-specific changes in gene expression programs of human keratinocytes and fibroblasts.

BMP signaling has a crucial role in skin development and homeostasis, whereas molecular mechanisms underlying its involvement in regulating gene expression programs in keratinocytes and fibroblasts remain largely unknown. We show here that several BMP ligands, all BMP receptors, and BMP-associated Smad1/5/8 are expressed in human primary epidermal keratinocytes and dermal fibroblasts. Treatment of both cell types by BMP-4 resulted in the activation of the BMP-Smad, but not BMP-MAPK pathways. Global microarray analysis revealed that BMP-4 treatment induces distinct and cell type-specific changes in gene expression programs in keratinocytes and fibroblasts, which are far more complex than the effects of BMPs on cell proliferation/differentiation described earlier. Furthermore, our data suggest that the potential modulation of cell adhesion, extracellular matrix remodeling, motility, metabolism, signaling, and transcription by BMP-4 in keratinocytes and fibroblasts is likely to be achieved by the distinct and cell-type-specific sets of molecules. Thus, these data provide an important basis for delineating mechanisms that underlie the distinct effects of the BMP pathway on different cell populations in the skin, and will be helpful in further establishing molecular signaling networks regulating skin homeostasis in health and disease.