Next Article in Journal
Genome Wide Transcriptome Analysis Reveals Complex Regulatory Mechanisms Underlying Phosphate Homeostasis in Soybean Nodules
Next Article in Special Issue
Advances in Understanding the Immunological Pathways in Psoriasis
Previous Article in Journal
Microparticles Produced by Activated Platelets Carry a Potent and Functionally Active Angiogenic Signal in Subjects with Crohn’s Disease
Previous Article in Special Issue
Oxidized Low-Density Lipoprotein-Deteriorated Psoriasis Is Associated with the Upregulation of Lox-1 Receptor and Il-23 Expression In Vivo and In Vitro
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 10
10.3390/ijms19102923
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Tissue-Engineered Human Psoriatic Skin Substitute: A Valuable In Vitro Model to Identify Genes with Altered Expression in Lesional Psoriasis

by
Geneviève Rioux
1,2,†,
Claudia Pouliot-Bérubé
1,2,†,
Mélissa Simard
1,2,
Manel Benhassine
3,4,
Jacques Soucy
5,
Sylvain L. Guérin
3,4 and
Roxane Pouliot
1,2,*
1
Centre LOEX de l’Université Laval, Génie Tissulaire et Régénération, Centre de Recherche FRQS du CHU de Québec, Axe Médecine Régénératrice, Québec, QC G1J 1Z4, Canada
2
Faculté de Pharmacie, Université Laval, Québec, QC G1V 0A6, Canada
3
Centre Universitaire d’Ophtalmologie-Recherche, Centre de Recherche FRQS du CHU de Québec, Axe Médecine Régénératrice, Québec, QC G1S4L8, Canada
4
Département d’Ophtalmologie, Université Laval, Québec, QC G1V 0A6, Canada
5
Département de Dermatologie, Hôpital de l’Enfant-Jésus, Québec, QC G1J 1Z4, Canada
*
Author to whom correspondence should be addressed.
†
These two authors contributed equally to the work and therefore should be considered as equivalent first authors.
Int. J. Mol. Sci. 2018, 19(10), 2923; https://doi.org/10.3390/ijms19102923
Submission received: 4 September 2018 / Accepted: 16 September 2018 / Published: 26 September 2018
(This article belongs to the Special Issue Molecular Mechanisms and Pathology of Psoriasis)
Browse Figures
Versions Notes

Abstract

:
Psoriasis is a chronic inflammatory skin disease for which no cure has emerged. Its complex etiology requires the development of an in vitro model representative of the pathology. In this study, we exploited gene profiling analyses on microarray in order to characterize and further optimize the production of a human psoriatic skin model representative of this in vivo skin disease. Various skin substitutes were produced by tissue-engineering using biopsies from normal, healthy donors, or from lesional or non-lesional skin samples from patients with psoriasis, and their gene expression profiles were examined by DNA microarray. We demonstrated that more than 3540 and 1088 genes (two-fold change) were deregulated between healthy/lesional and lesional/non-lesional psoriatic substitutes, respectively. Moreover, several genes related to lipid metabolism, such as PLA2G4E and PLA2G4C, were identified as repressed in the lesional substitutes. In conclusion, gene profiling analyses identified a list of deregulated candidate genes associated with various metabolic pathways that may contribute to the progression of psoriasis.

1. Introduction

Psoriasis is a chronic skin affection characterized by well-demarcated red scaly plaques that appear most of the time on elbows, knees, the sacroiliac region, nails, and the scalp [1,2]. The disease affects both men and women in the same proportion, and can be classified into two subtypes, according to the onset of the disease [3,4]. Plaque psoriasis is the most common, but other forms, such as erythrodermic psoriasis, inverse, pustular, and guttate, are frequently observed [5,6,7]. The disease has a different incidence depending on ethnicity and geographical distribution [8]. The latest statistics indicate that 7.4 million people in the United States, and 125 million individuals worldwide are affected by psoriasis [9]. The physical consequences of plaque psoriasis are more important than the simple fact that they provide pain and pruritus. Indeed, the impacts on general health go beyond what we previously believed as there is a substantial list of psoriasis-associated comorbidities, including Crohn’s disease, psoriatic arthritis, atherogenic dyslipidemia, hypertension, diabetes, and several more [10]. The psychological consequences are worth mentioning as many patients will suffer from anxiety and depression [11].
Despite being an incurable disease, several treatments have emerged. They are essentially based on the deep knowledge we have about immunology and genetics, the two main components of the pathogenesis that makes psoriasis a very heterogeneous disease [12,13]. In this sense, the quest for the miracle drug to treat psoriasis is not a simple task. In addition, a better psoriasis model could help reduce the huge costs associated with the development of new drugs aimed at reducing the skin lesions caused by this disease [14]. Then, a more accurate model for an early stage of preclinical research would perhaps eliminate molecules earlier in the developing process. Indeed, research in this area is hampered by a lack of a model that adequately represents the pathology. The models used for psoriasis research can be classified into two broad, well-described categories: in vivo and in vitro [15,16]. In vivo models, that mostly exploit murine models, are divided into three sub-categories: spontaneous, genetically engineered, and xenotransplantation [17,18]. These models have several strengths and weaknesses, but none accurately represent human skin. In vitro models, in turn, can be valuable research tools to study specific pathways [19]. The general classification for in vitro psoriatic models relies on whether they represent 2D or 3D reconstructed skin substitutes. The monolayer model (2D) was largely used despite the fact that no interactions can be observed between the different cell subtypes that normally constitute the human skin, which also caused an increased interest in the 3D model [20]. Our team has been working for several years on a reconstructed human psoriatic skin that uses fibroblasts and keratinocytes from psoriatic plaques [21]. Our psoriatic skin model has been extensively characterized for its biochemical and immunohistochemical properties [22,23,24]. However, the genetic component, very important in the pathology of this disease, has not yet been explored. The present study was aimed at comparing the pattern of genes expressed in our model of human psoriatic lesional skin with that of healthy and non-lesional skin substitutes produced using the same technique.

2. Results

2.1. Macroscopic and Histological Analyses

An examination of reconstructed skin substitutes under phase contrast microscopy showed that the healthy substitutes (Figure 1A–C) all had the same macroscopic aspect without any differences related to cell populations. Indeed, the appearance of the skin substitutes was regular and uniform. The macroscopic aspect of the non-lesional substitute showed two different phenotypes, one closer to the healthy skin substitutes (Figure 1D) and the other being more similar to the lesional skin substitutes (Figure 1E,F). Skin substitutes reconstructed with lesional cells isolated from psoriatic plaques were very similar to one another, all being irregular with thinner and thicker skin areas (Figure 1G–I).
The histological appearance of healthy skin substitutes reconstructed from different human cell populations appeared similar to one another (Figure 2A–C). Indeed, the thickness of the dermis (in blue) did not show any significant difference, and the epidermis (in purple) maintained a uniform thickness. On the other hand, the non-lesional skin substitutes had a thicker epidermis than that seen for healthy substitutes (Figure 2D–F). The lesional skin substitutes demonstrated a thicker epidermis than the healthy skin substitutes and showed a more disorganized structure (Figure 2G–I).

2.2. Gene Profiling Analysis of the Most Deregulated Genes between Healthy, Lesional and Non-Lesional Skin Substitutes

We next conducted a gene profiling analysis on the microarray in order to compare the gene expression profile between healthy, lesional, and non-lesional skin substitutes. The scatter plot analysis of the 60,000 probes loaded on the chip showed important changes in the pattern of genes expressed by the lesional (L) skin substitutes against the healthy (H) skin substitutes (Figure 3A, middle panel; R2 = 0.9525). The comparison between the healthy and the non-lesional (NL) skin substitutes revealed that the number of genes deregulated between these two conditions decreased (Figure 3A, left panel; R2 = 0.9686). Very much the same was observed when deregulated genes from both the lesional and non-lesional skin substitutes were analyzed against each other (Figure 3A, right panel; R2 = 0.9739).
A heatmap representation for all the genes showing a two-fold or more expression variation unique to the non-lesional or lesional substitutes relative to the healthy substitutes was generated (Figure 3B). An analysis of the heatmap from the healthy against the lesional condition indicates that 3540 candidate genes were deregulated by more than two-fold (Figure 3B, left; Figure 3C). The number of genes that fit into this category dropped to 1088 (69% reduction) when non-lesional and lesional conditions were compared to each other (Figure 3B, right; Figure 3C).
We next analyzed the data files from the microarray to sort out genes whose expression was the most greatly deregulated between the lesional and either the non-lesional or healthy conditions (Figure 3D). The Arraystar program was used to restrict the search to the 55 most deregulated genes in the healthy against lesional condition and in the non-lesional against lesional condition. However, in order to avoid inter-individual variations that may occur in the pattern of expressed genes (see Figure S1), data shown in Figure 3D (left and right columns) were collected from the comparison between pools of reconstructed skin substitutes produced from two healthy cell populations and from lesional and non-lesional reconstructed skin substitutes produced from three different cell populations. Among these highly deregulated genes, 12 (GUCA2B, LCOR, LRG1, RDH12, NIPAL4, FGFR3, C10orf99, FGFBP1, KRT1, KREMEN2, AKR1B10, and PRSS56) were similarly deregulated when comparing the healthy and lesional, and the non-lesional and lesional, conditions (genes in red on Figure 3D).

2.3. Alteration in Cytokines, Chemokines and Growth Factors Gene Expression in Healthy vs. Lesional Skin Substitutes

As the immune system is expected to have a significant impact on the pattern of genes expressed by the pathologic skin of patients with psoriasis, we therefore searched our gene profiling data files to sort out cytokines, chemokines, growth factors, and some of their receptors encoding genes whose expression is altered between healthy and lesional skin substitutes. From the 257 cytokine/growth factor genes present on the array, only 14 had their expression altered by more than two-fold (Table 1), with the remaining genes being either not expressed or having their expression level unaltered in both conditions (healthy vs. lesional). Among the affected genes, IL1R2, CXCL13, CXCL14, EPO, INHBA, and TNFSF9 had their expression increased between 2.1 and 6.8-fold, whereas the expression of IL15, CXCL2, CCL20, CXCL1, LTB, TNFRSF10A, IL24, and CCL27 was reduced by a factor of 2.1 to 6.6-fold in lesional compared with healthy skin substitutes. Expression of the gene encoding cathelicidin antimicrobial peptide (CAMP), whose expression has been reported to increase in psoriasis [25], was also found to increase by 2.4-fold in our lesional skin substitutes (Table 1).

2.4. Gene Ontology Analysis

Gene ontology (GO) assignments were conducted on the 55 most deregulated genes (Figure 3D), using the amigo gene ontology database, to illustrate which biological processes are altered between the different conditions (Table 2) [26]. The most deregulated pathways in the healthy against the lesional condition include ‘keratinization’ (p = 3.625 × 10−4), ‘isoprenoid metabolic process’ (p = 5.922 × 10−4), and ‘retinoid metabolic process’ (p = 2.321 × 10−3). The products encoded by the most strongly deregulated genes between the non-lesional against the lesional condition are particularly involved in ‘skin development’ (p = 1.151 × 10−11), ‘keratinization’ (p = 7.952 × 10−6), and ‘epidermis development’ (p = 1.088 × 10−5). A closer and more extended examination of the genes identified as deregulated in our lesional substitutes within each of these biological processes is presented in Table S1. Of the 28 genes selected within these biological functions and identified as deregulated in psoriasis in both the present study and that reported by Gudjonsson et al. [27], a total of 19 (68%) were similarly deregulated in both analyses, therefore further validating the reliability of our model.
The results presented in Table 3 illustrate the differences (both the linear signals and fold-changes are indicated) in the expression of genes involved in lipid synthesis and the establishment of the skin barrier. The conditions compared are: healthy and non-lesional substitutes (H vs. NL), healthy and lesional substitutes (H vs. L), and non-lesional and lesional substitutes (NL vs. L). Several genes related to collagen expression were more strongly deregulated between healthy and lesional substitutes (H vs. NL). Also, several enzymes related to lipid metabolism such as phospholipases, protein kinases, and lipoprotein lipase had their expression altered in the lesional (H vs. KL) and non-lesional (H vs. NL) substitutes compared with healthy ones.

3. Discussion

Despite the fact that the cellular and immunohistological aspects of our in vitro psoriatic model were largely characterized [21,24,28], they missed the detailed transcriptomic analyses that could help verify the integrity of our reconstructed psoriatic skin. To fulfill this need, we produced a large number of healthy and both lesional and non-lesional psoriatic substitutes in order to demonstrate whether the progression of this disease might result from alterations in the pattern of genes expressed by the cells from these reconstructed tissues. Both macroscopic and immunohistochemical analyses revealed that the healthy substitutes presented in Figure 1A–C have a very similar phenotype to one another. A well-differentiated stratum corneum that contributes to maintaining an appropriate barrier function is also observed. Previous studies have determined that the lipid domain of the stratum corneum from healthy substitutes was comparable to that found in the in vivo healthy skin [28].
Unlike healthy substitutes, psoriatic substitutes made from non-lesional cells presented phenotypic variations between the different cell populations tested (Figure 1D–F). First, we observed that one of them (Figure 1D) had a phenotype similar to that observed with the healthy substitutes, which means that the appearance was regular and opaque. Nonetheless, others (such as Figure 1E,F) appeared similar to the psoriatic substitutes by their irregular and contracted appearance. As a result, two different profiles could be proposed for these substitutes: one related to the healthy and the other to the psoriatic phenotype. These results support previous research conducted in our laboratory where two different expression profiles were confirmed by histological, immunohistochemical, and macroscopic analyses for substitutes reconstructed using non-lesional psoriatic cells [22]. Finally, the macroscopic and histological data support the psoriatic phenotype of our lesional psoriatic substitutes, as also described in previous work [21]. However, it is important to point out that the histological features of the cell population used in Figure 2I appeared to be less pronounced than the other two cell populations used for the reconstruction of lesional psoriatic skin substitutes. Indeed, several studies on psoriasis characterized it as a cyclic disease, with periods of remission and exacerbation of plaques, meaning that the characteristics of the plaques change over time, sometimes appearing to be exacerbated or healing [29,30]. The state wherein the psoriatic plaque was at the time of the biopsy may be an intrinsic component of the cells used to produce the reconstructed psoriatic substitutes.
A recent review by Niehues et al. has listed the most complete and best-validated models that include major psoriasis hallmarks with regard to the gene and protein expression profile and epidermal morphology [31]. Our model appears to be the only one that fulfills the complete expression signatures related to acanthosis, parakeratosis, and hyperproliferation. Although tissue-engineered, 3D models are structures that advantageously represent the in vivo condition, the fact remains that the in vitro results are not necessarily the perfect mirror of the in vivo observations [32]. Therefore, to extensively characterize and validate the psoriasis-like features of the 3D models, Niehues et al. suggested the use of not just one marker, but rather of a panel of disease-associated genes and proteins for the validation of novel psoriasis models [31].
The data files from the gene profiling studies were first examined using scatter plot analyses and have provided particularly interesting details on important changes occurring in the pattern of genes expressed in our three studied conditions (healthy, non-lesional, and lesional). A total of five cell populations were pooled to represent the healthy condition (H), whereas both the lesional (L) and non-lesional (NL) conditions were each constituted of three and four pooled cell populations, respectively. Three scatter plots have been generated from these data. First, the graph of the healthy against the lesional condition identified a significant amount of deregulated genes, which was supported by the slope of the linear regression of 0.9525. Thereafter, the graph with the most reduced number of deregulated genes was that which compared the non-lesional against the lesional condition, with an R2 value of 0.9739. Based on these two scatter plots, it can be concluded that there are many deregulated genes between healthy and lesional substitutes, but much less between non-lesional and lesional cells from the same patient. This graph is particularly important since the comparison that is made here did not take into account the inter-individual differences since the samples (lesional and non-lesional) came from the same four patients. Moreover, one of the most important points to emerge from this study was the scatter plot comparing the healthy against the non-lesional substitutes. The latter had a linear regression slope of 0.9686, demonstrating that the healthy skin of an individual who has never been affected by psoriasis is different from healthy skin of a patient that suffers from psoriasis. To our knowledge, there are no non-lesional psoriatic skin models available on the market to which we can refer in order to establish any possible correlation with our results, despite the fact that the non-lesional skin has been increasingly characterized [33,34]. In regard to the in vivo skin, in all studies encountered, morphological and histological characteristics were similar between the healthy skin of a patient that had never suffered from psoriasis and the healthy skin of a psoriatic patient. Nevertheless, a few studies have reported that many genes, whose encoded protein products are associated with the processes of angiogenesis and lymphogenesis, are overexpressed in non-lesional skin [35,36]. Our gene profiling data also show that despite similar histological characteristics, non-lesional psoriatic cells demonstrate a disruption in their gene expression pattern that is less important than the one observed for lesional tissues, but more important than the healthy one.
The Venn diagram analysis was then used to define the number of deregulated genes specific to each experimental condition. An examination of the diagram revealed that a large number of genes are expressed in very different ways between healthy and lesional substitutes. Indeed, of the previous 3540 genes identified as deregulated between healthy and lesional substitutes, 2850 genes were found to be specific to that condition (healthy vs. lesional), whereas 690 genes were also found to be deregulated between non-lesional and lesional substitutes. In the same sense, 398 genes out of the 1088 genes previously identified for non-lesional and lesional substitutes were also found to be specific to this condition (non-lesional vs. lesional). It was no surprise to note a greater number of deregulated genes between healthy and lesional than non-lesional and lesional substitutes. Interestingly, and based on these results, we can assume that the 398 genes specifically deregulated between the non-lesional and lesional conditions most likely account for the progression of the non-lesional state to the lesional one.
The type of array used and the nature of the samples, for instance, tissue-engineered psoriatic skin substitutes, make our work even more innovative, but also very difficult to compare to those published in the literature. Indeed, several groups have studied the pattern of genes expressed by lesional and non-lesional skin in vivo [33,34]. Among them, Gudjonsson et al. identified an impressive amount of deregulated genes between psoriatic and healthy skin (1326 genes), while 1085 genes were identified as deregulated between lesional and non-lesional skin [27]. Oestreicher et al. also generated a complete list of 159 genes that define psoriasis in molecular terms and whose expression was compared between lesional and non-lesional psoriatic skin [34]. These data are in accordance with those presented in the present study. Indeed, from the 41 randomly selected genes among those identified as deregulated in psoriasis by Gudjonsson et al. and Oestreicher et al., a total of 22 (54%) were found to be similarly deregulated (repressed or activated) in our psoriatic skin substitutes, therefore demonstrating the reliability of our in vitro model (Table S2). In spite of this, the difference in the expression we observed for some of these genes in regard to other studies may rely on differences such as the type of array used and the nature of the samples. It is important to stress that the RNAs extracted directly from biopsies do not go through the same steps as those extracted from our reconstructed tissues. Indeed, and unlike RNA extracted from biopsies, the cells of the reconstructed psoriatic tissues from which RNA was extracted were exposed to additional stresses that can also impact the pattern of expressed genes. Our reconstructed psoriatic skin only contains keratinocytes and fibroblasts, whereas a skin biopsy contains all the different cell types from the skin. Some cell types may have a gene expression profile that stands out in the dataset obtained in the ex vivo studies, whereas the profile in our study only involves skin fibroblasts and keratinocytes. This aspect may also explain the difference observed between the two studies.
The mystery of whether the primary abnormality of psoriasis resides in epidermal keratinocytes or dermal immunocytes remains unsolved. Conventional wisdom portrays psoriasis as an autoimmune disease, where misguided T lymphocyte activities cause secondary epithelial abnormalities despite the fact that a few researchers believe that psoriasis is a genetically determined, abnormal epithelial response to infectious and/or physicochemical skin insults [37]. There is a cumulative body of evidence that the psoriatic epidermis has both structural and functional abnormalities [38]. It is therefore plausible that psoriatic keratinocytes tease the immune system, waking it from its slumber to secrete a barrage of cytokines and factors in a misguided attempt to restore keratinocytes to their normal proliferative state. Consistent with this hypothesis, we observed variations in the expression of a few cytokines or growth factors by our lesional skin substitutes, despite the fact that they are devoid of any immune cells, suggesting that their expression is likely keratinocyte-specific. In agreement with the results generated using a mouse model of skin psoriasis [39], a near seven-fold increase in the expression of the CXCL13 chemokine was observed in our lesional skin substitutes relative to the level observed in the healthy substitutes. It is also particularly interesting to point out that the near seven-fold repression we observed for CCL27 gene expression in our lesional substitutes was also reported in punch biopsies from human psoriatic lesions [40]. Consequently, our in vitro model can, indeed, be used to identify genes with altered expression in psoriasis, and maybe new, yet unidentified genes as well, such as CXCL2, CXCL14, INHBA, EPO, TNFSF9, and IL1R2. Therefore, our results confirmed and extended previous studies highlighting the remarkable potential of keratinocytes to produce pro-inflammatory cytokines, adhesion molecules, growth factors, and chemotactic polypeptides [23].
Several research teams around the world are working on the gene expression of healthy, non-lesional and lesional skin, allowing us to compare our in vitro model to theirs. A few studies reported that increased gene expression in psoriatic skin is partly attributable to keratinocytes, T lymphocytes, and macrophages, but also to stimulation by pro-inflammatory cytokines [41,42]. In contrast, the repression of some specific genes in the psoriatic plaques would be mediated by subcutaneous adipose tissues and the epidermis [43]. These findings were made possible through ontological analyses that assign genes differentially expressed to biological processes or cell types. Since only two cell types are present in our model (fibroblasts and keratinocytes), an analysis of the biological processes has been achieved using the 55 most deregulated genes identified in our study. Among these biological activities, the keratinization process was greatly altered across all the deregulated genes, as well as that of isoprenoid metabolism and retinoids in lesional and healthy substitutes. As for non-lesional and lesional substitutes, the most altered biological processes were more related to skin and epidermis development. Although our model is not entirely complete as it lacks vascularization, innervation, and immune cells, it nevertheless demonstrates similarities with the metabolic pathways affected in psoriasis [2]. So far, there is no doubt that the keratinization process is greatly affected in psoriasis, also leading to disturbances in the organization and composition of lipids of various epidermal layers, thus explaining the ontology results [44,45]. In contrast to the various lipid changes orchestrated into the psoriatic plaque itself, healthy skin areas in psoriatic patients seem to have very little lipid changes [46].
Table 3 helped identify several key players in the lipid metabolism of the stratum corneum and skin in general. These genes are repressed in the lesional psoriatic substitutes compared to the healthy substitutes. Indeed, the phospholipase A2 epsilon and gamma genes (PLA2G4E and PLA2G4C) were downregulated by 3.7- and 3.9-fold, respectively. Type A2 phospholipase catalyzes the hydrolysis of glycerophospholipids, thereby releasing arachidonic acid and lysophospholipids [47]. It is interesting to note that lysophospholipids are mediators of cell communication that can affect various cellular processes, including proliferation/apoptosis and cellular remodeling [48]. Specifically, phospholipase A2 epsilon (PLA2G4E) is involved in the input mechanisms of the cell, or the clathrin-independent mechanism of tubule formation involved in compound transportation to the cell surface [49]. The phospholipase A2 gamma (encoded by the gene PLA2G4C) participates in the constitutive release of various lipids, including arachidonic acid, oleic acid, and some fatty acids [50]. The phospholipase A2 gamma also plays a role in oxidative stress mechanisms in order to fix lipid damage [51]. Psoriasis is known and recognized as having disturbances in lipid, as well as oxidant-antioxidant, metabolism [52,53]. These imbalances may be linked to the products of genes for which expression is altered, such as is the case with PLA2G4E and PLA2G4C. Unfortunately, no direct link has been reported between the expression of these genes and the imbalances mentioned above. The gene for dystonin (DST), a member of the plakines family associated with adhesion and cellular junction through intermediate filament, was significantly down-regulated between healthy and involved substitutes (four-fold change) [54]. No direct link has been mentioned for any modification of its expression in psoriasis; however one of the isoforms causes the expression of a protein called BPAG1-e (bullous pemphigoid antigen 1), which is found in bullous pemphigoid skin disease [55]. This protein would affect both the migration and adhesion of keratinocytes [56]. In addition, besides its interaction with both clathrin and actin, dystonin participates in the formation of the hemidesmosomes through its interaction with the β4 integrin subunit [57,58]. Interestingly, a few studies reported that some patients with psoriasis also developed bullous dermatosis spontaneously or following medical treatments [59,60,61], thereby suggesting that a relationship between the autoimmune disease of bullous pemphigoid and psoriasis exists.
Detailed analysis of the genes shown in Figure 3D allowed us to make multiple comparisons with a number of genes associated with the late cornified envelope family, as is the case for LCE2A, LCE2B, LCE2C, LCE2D, and LCE6A. Indeed, these genes were down-regulated in lesional compared to non-lesional substitutes. It is particularly interesting to point out that the protein products of these genes are associated with the stratification of the epidermis. Consistent with our observations, deletion within the LCE3B and LCE3C genes (LCE3B_LCE3C) was shown to be associated with an increased risk of developing psoriasis [62,63]. Similarly, genes associated with loricrin (LOR), various keratins (KRT1, KRT2, KRT5, KRT31, and KRT77), and interleukins (IL1A and IL1B) were also found to be strongly deregulated in our human psoriatic skin model compared to the reconstructed non-lesional skin [64,65].
Recent findings have established the skin as a peripheral neuroendocrine organ that is tightly networked to central stress axes. This capability contributes to the maintenance of the skin’s and body’s homeostasis. Specifically, epidermal and dermal cells produce and respond to classical stress neurotransmitters, neuropeptides, and hormones, and this production is modified by ultraviolet radiation and biological, chemical, and physical factors [66,67,68]. Therefore, skin removed from the natural neuro-immuno-endocrine connections in the body is certainly deprived of the local neuroendocrine capabilities that skin can normally compensate for. Although this withdrawal from central communication may affect the results of our research, it remains interesting to think that in view of local neuro-endocrine-cytokine activities, organ culture models can provide clinically useful information [66,68,69,70,71]. For instance, it has been demonstrated that the corticotrophin-releasing factor, the pro-opiomelanocortin, and corresponding receptors were co-expressed in cultured keratinocytes, melanocytes, or dermal fibroblasts [72,73]. Moreover, the expression of the executive arm of the cutaneous hypothalamo-pituitary adrenal axis, i.e., the production of cortisol and corticosterone, has been clearly demonstrated in cultured epidermal keratinocytes and melanocytes, as well as in dermal fibroblasts [73,74,75,76,77]. Nevertheless, it is important to stress that the model described in our study, even if more complete than cultured cells grown as monolayers, still remains a simplified representation of one of the most important and complex organs of the human body: the skin.
In conclusion, transcriptome analysis of our tissue-engineered lesional and non-lesional substitutes allowed us to establish a list of the most deregulated genes whose encoded products can contribute to the induction and/or the progression of this disease. Moreover, we have been able to associate these genes with altered biological processes and make correlations with other pathologies. The changes in the expression noted for the genes herein identified may have resulted from processes within the cell itself, but also from feedback activation mechanisms and repression in various metabolic pathways. Plaque psoriasis is a disease for which causes seem to be multifactorial, making the creation of a unique gene profile implausible. The more there will be analyses of gene expression, methylation status, and post-transcriptional processes, the clearer the explanatory picture of the disease will be.

4. Materials and Methods

This study was conducted in agreement with the Helsinki declaration and was performed under the guidelines of the research ethics committee of the ‘CHU de Québec’ (ethic code: DR-002-1387, protocol renewal approved on 6 June 2018). All patients were given adequate information to provide written consents.

4.1. Cell Culture

Fibroblasts were cultured in Dulbecco–Vogt modification of Eagle’s medium (DMEM) supplemented with 10% bovine growth serum (HyClone®, Thermo Fisher Scientific, Ottawa, ON, Canada), 100 UI/mL penicillin G (Sigma, Oakville, ON, Canada), and 25 µg/mL gentamicin (Schering, Pointe-Claire, QC, Canada). Keratinocytes were cultured in a combination of DMEM with Ham’s F12 (3:1), supplemented with 5% Fetal Clone II serum (Hyclone), 5 µg/mL insulin (Sigma), 0.4 µg/mL hydrocortisone (Calbiochem, EMD Biosciences, Gibbstown, NJ, USA), 10−10 M cholera toxin (MP Biomedicals, Montreal, QC, Canada), 10 ng/mL human epidermal growth factor (EGF; Austral Biological, San Ramon, CA, USA), 100 UI/mL penicillin G (Sigma), and 25 µg/mL gentamicin (Schering).

4.2. Production of Tissue-Engineered Substitutes

All skin substitutes (healthy, lesional, and non-lesional) were produced using the self-assembly method with some modifications [21,78]. Different substitutes were produced using nine different cell populations; five from healthy patients (18-, 22-, 23-, 38-, and 46-year old) and four from psoriatic patients (47-, 49-, 64-, and 70-year old). The psoriatic cells were used to produce lesional and non-lesional substitutes. Briefly, fibroblasts (passage 6) were cultured in the presence of ascorbic acid (Sigma) at a concentration of 50 µg/mL. After 28 days of culture, these cells formed sheets that could be superimposed and incubated further for seven days to form a new dermal layer before keratinocytes (passage 2) were seeded on top of it to form a new epidermal layer. After one week of submerged culture in DMEM-HAM, the substitutes were raised to the air–liquid interface. Biopsies of skin substitutes were taken after 21 days of culture at the air–liquid interface and analyzed by gene profiling on DNA microarrays [79].

4.3. Gene Expression Profiling

Total RNA was isolated from five different healthy (18-, 22-, 23-, 38-, and 46-year old) and eight different psoriatic reconstructed skin substitutes (four lesional and four non-lesional (47-, 49-, 64-, and 70-year old)) using the RNeasy Mini Kit (QIAGEN, Toronto, ON, Canada). Cyanine 3-CTP labeled cRNA targets were prepared from 50 ng of total RNA using the Agilent One-Color Microarray-Based Gene Expression Analysis kit (Agilent Technologies, Santa-Clara, CA, USA). Then, 600 ng cRNA was incubated on a G4851A SurePrint G3 Human GE 8 × 60 K array slide (60,000 probes, Agilent Technologies). Slides were then hybridized and scanned on an Agilent SureScan Scanner. Data were analyzed using ArrayStar V4.1 (DNASTAR, Madison, WI, USA) software for the generation of heat maps and scatter plots. Biological replicates were as follows: for the experiment conducted on the healthy tissue-engineered skin, total RNA was obtained from two different reconstructed skin substitutes produced using keratinocytes and fibroblasts cultured from the skin of two different healthy donors (22- and 46-year old); for the experiments that used reconstructed psoriatic skin substitutes, RNA was obtained from four different reconstructed non-lesional (NL) and lesional (L) skin substitutes produced using keratinocytes and fibroblasts cultured from the skin of four different psoriatic donors (47-, 49-, 64-, and 70-year old). All data generated from the arrays were also analyzed by robust multiarray analysis (RMA) for background correction of the raw values. They were then transformed in Log2 base and quantile normalized before a linear model was fitted to the normalized data to obtain an expression measure for each probe set on each array. All microarray data presented in this study comply with the Minimum Information about a Microarray Experiment (MIAME) requirements (GEO # (GSE120464, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE120464)).

4.4. Histological Analysis

Three biopsies of each reconstructed tissue combination were fixed with HistoChoice’s solution and embedded in paraffin wax. Five-micrometer-thick sections were cut and stained with Masson’s Trichrome (Zeiss, Axio Imager, North York, ON, Canada).

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/10/2923/s1.

Author Contributions

Conceptualization, R.P. and S.L.G.; Data curation, G.R., C.P.-B., M.S., and M.B.; Formal analysis, C.P.-B.; Funding acquisition, R.P.; Investigation, M.S. and J.S.; Writing—original draft, C.P.-B.; Writing—review & editing, R.P. and S.L.G.

Acknowledgments

The authors are grateful to Gaëtan Le-Bel (CUO-Recherche, Centre de Recherche FRQS du CHU de Québec) for his help with the analysis of the microarray data files. This study was supported by the Canadian Institutes of Health Research (CIHR) to R. Pouliot (MOP-311262). The ‘Fonds d’Enseignement et de Recherche’ (FER) of the Faculté de Pharmacie, Université Laval, Québec, Canada (scholarship to G Rioux). R Pouliot is an FRQS career award scholar.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lowes, M.A.; Bowcock, A.M.; Krueger, J.G. Pathogenesis and therapy of psoriasis. Nature 2007, 445, 866–873. [Google Scholar] [CrossRef] [PubMed]
  2. Stern, R.S. Psoriasis. Lancet 1997, 350, 349–353. [Google Scholar] [CrossRef]
  3. Smith, R.L.; Warren, R.B.; Griffiths, C.E.; Worthington, J. Genetic susceptibility to psoriasis: An emerging picture. Genome Med. 2009, 1, 72. [Google Scholar] [CrossRef] [PubMed]
  4. Henseler, T.; Christophers, E. Psoriasis of early and late onset: Characterization of two types of psoriasis vulgaris. J. Am. Acad. Dermatol. 1985, 13, 450–456. [Google Scholar] [CrossRef]
  5. Boyd, A.S.; Menter, A. Erythrodermic psoriasis. J. Am. Acad. Dermatol. 1989, 21, 985–991. [Google Scholar] [CrossRef]
  6. Menter, A.; Barker, J.N. Psoriasis in practice. Lancet 1991, 338, 231–234. [Google Scholar] [CrossRef]
  7. Naldi, L.; Gambini, D. The clinical spectrum of psoriasis. Clin. Dermatol. 2007, 25, 510–518. [Google Scholar] [CrossRef] [PubMed]
  8. Ryan, C.; Korman, N.J.; Gelfand, J.M.; Lim, H.W.; Elmets, C.A.; Feldman, S.R.; Gottlieb, A.B.; Koo, J.Y.; Lebwohl, M.; Leonardi, C.L.; et al. Research gaps in psoriasis: Opportunities for future studies. J. Am. Acad. Dermatol. 2014, 70, 146–167. [Google Scholar] [CrossRef] [PubMed]
  9. Rachakonda, T.D.; Schupp, C.W.; Armstrong, A.W. Psoriasis prevalence among adults in the United States. J. Am. Acad. Dermatol. 2014, 70, 512–516. [Google Scholar] [CrossRef] [PubMed]
  10. Grozdev, I.; Korman, N.; Tsankov, N. Psoriasis as a systemic disease. Clin. Dermatol. 2014, 32, 343–350. [Google Scholar] [CrossRef] [PubMed]
  11. Gupta, M.A.; Gupta, A.K. Psychiatric and Psychological Co-Morbidity in Patients with Dermatologic Disorders. J. Am. Acad. Dermatol. 2003, 4, 833–842. [Google Scholar] [CrossRef]
  12. Weigle, N.; Mcbane, S. Psoriasis. Am. Fam. Physician 2013, 87, 626–633. [Google Scholar] [PubMed]
  13. Nickoloff, B.J.; Qin, J.Z.; Nestle, F.O. Immunopathogenesis of psoriasis. Clin. Rev. Allergy Immunol. 2007, 33, 45–56. [Google Scholar] [CrossRef] [PubMed]
  14. Zollner, T.M.; Renz, H.; Igney, F.H.; Asadullah, K. Animal models of T-cell-mediated skin diseases. Bioessays 2004, 26, 693–696. [Google Scholar] [CrossRef] [PubMed]
  15. Schon, M.P. Animal Models of Psoriasis—What Can We Learn from Them? J. Investig. Dermatol. 1999, 112, 405–410. [Google Scholar] [CrossRef] [PubMed]
  16. Jean, J.; Pouliot, R. In vivo and In vitro Models of Psoriasis. In Tissue Engineering; Eberli, D., Ed.; InTech Publishers: Novi Sad, Serbia, 2010; pp. 359–382. ISBN 978-953-307-079-7. [Google Scholar]
  17. Danilenko, D.M. Review paper: Preclinical models of psoriasis. Vet. Pathol. 2008, 45, 563–575. [Google Scholar] [CrossRef] [PubMed]
  18. Mizutani, H.; Yamanaka, K.; Konishi, H.; Murakami, T. Animal models of psoriasis and pustular psoriasis. Arch. Dermatol. Res. 2003, 295 (Suppl. 1), S67–S68. [Google Scholar] [CrossRef] [PubMed]
  19. Benam, K.H.; Dauth, S.; Hassell, B.; Herland, A.; Jain, A.; Jang, K.J.; Karalis, K.; Kim, H.J.; MacQueen, L.; Mahmoodian, R.; et al. Engineered in vitro disease models. Annu. Rev. Pathol. 2015, 10, 195–262. [Google Scholar] [CrossRef] [PubMed]
  20. Jean, J.; Garcia-Perez, M.E.; Pouliot, R. Psoriatic Skin Models: A Need for the Pharmaceutical Industry. In Psoriasis; Soung, J., Ed.; InTech Publishers: Novi Sad, Serbia, 2011; pp. 47–62. ISBN 978-953-307-878-6. [Google Scholar]
  21. Jean, J.; Lapointe, M.; Soucy, J.; Pouliot, R. Development of an in vitro psoriatic skin model by tissue engineering. J. Dermatol. Sci. 2009, 53, 19–25. [Google Scholar] [CrossRef] [PubMed]
  22. Jean, J.; Leroy, M.; Duque-Fernandez, A.; Bernard, G.; Soucy, J.; Pouliot, R. Characterization of a psoriatic skin model produced with involved or uninvolved cells. J. Tissue Eng. Regen. Med. 2012, 9, 789–798. [Google Scholar] [CrossRef] [PubMed]
  23. Martin, G.; Guerard, S.; Fortin, M.M.; Rusu, D.; Soucy, J.; Poubelle, P.E.; Pouliot, R. Pathological crosstalk in vitro between T lymphocytes and lesional keratinocytes in psoriasis: Necessity of direct cell-to-cell contact. Lab. Investig. 2012, 92, 1058–1070. [Google Scholar] [CrossRef] [PubMed]
  24. Ayata, R.E.; Bouhout, S.; Auger, M.; Pouliot, R. Study of in vitro capillary-like structures in psoriatic skin substitutes. Biores. Open Access 2014, 3, 197–205. [Google Scholar] [CrossRef] [PubMed]
  25. Morizane, S.; Yamasaki, K.; Muhleisen, B.; Kotol, P.F.; Murakami, M.; Aoyama, Y.; Iwatsuki, K.; Hata, T.; Gallo, R.L. Cathelicidin antimicrobial peptide LL-37 in psoriasis enables keratinocyte reactivity against TLR9 ligands. J. Investig. Dermatol. 2012, 132, 135–143. [Google Scholar] [CrossRef] [PubMed]
  26. Carbon, S.; Ireland, A.; Mungall, C.J.; Shu, S.; Marshall, B.; Lewis, S. AmiGO: Online access to ontology and annotation data. Bioinformatics 2009, 25, 288–289. [Google Scholar] [CrossRef] [PubMed]
  27. Gudjonsson, J.E.; Ding, J.; Johnston, A.; Tejasvi, T.; Guzman, A.M.; Nair, R.P.; Voorhees, J.J.; Abecasis, G.R.; Elder, J.T. Assessment of the psoriatic transcriptome in a large sample: Additional regulated genes and comparisons with in vitro models. J. Investig. Dermatol. 2010, 130, 1829–1840. [Google Scholar] [CrossRef] [PubMed]
  28. Bernard, G.; Auger, M.; Soucy, J.; Pouliot, R. Physical characterization of the stratum corneum of an in vitro psoriatic skin model by ATR-FTIR and Raman spectroscopies. Biochim. Biophys. Acta 2007, 1770, 1317–1323. [Google Scholar] [CrossRef] [PubMed]
  29. Koo, J.; Lebwohl, M. Duration of remission of psoriasis therapies. J. Am. Acad. Dermatol. 1999, 41, 51–59. [Google Scholar] [CrossRef]
  30. Ragaz, A.; Ackerman, A.B. Evolution, maturation, and regression of lesions of psoriasis. Am. J. Dermatopathol. 1979, 1, 199–214. [Google Scholar] [CrossRef] [PubMed]
  31. Niehues, H.; van den Bogaard, E.H. Past, present and future of in vitro 3D reconstructed inflammatory skin models to study psoriasis. Exp. Dermatol. 2018, 27, 512–519. [Google Scholar] [CrossRef] [PubMed]
  32. Chiricozzi, A.; Romanelli, M.; Panduri, S.; Donetti, E.; Prignano, F. Relevance of in vitro 3-D skin models in dissecting cytokine contribution to psoriasis pathogenesis. Histol. Histopathol. 2017, 32, 893–898. [Google Scholar] [PubMed]
  33. Zhou, X.; Krueger, J.G.; Kao, M.C.; Lee, E.; Du, F.; Menter, A.; Wong, W.H.; Bowcock, A.M. Novel mechanisms of T-cell and dendritic cell activation revealed by profiling of psoriasis on the 63,100-element oligonucleotide array. Physiol. Genom. 2003, 13, 69–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Oestreicher, J.; Walters, I.B.; Kikuchi, T.; Gilleaudeau, P.; Surette, J.; Schwertschlag, U.; Dorner, A.J.; Krueger, J.G.; Trepicchio, W.L. Molecular classification of psoriasis disease-associated genes through pharmacogenomic expression profiling. Pharmacogenomics 2001, 1, 272–287. [Google Scholar] [CrossRef] [Green Version]
  35. Henno, A.; Blacher, S.; Lambert, C.; Colige, A.; Seidel, L.; Noel, A.; Lapiere, C.; de la Brassinne, M.; Nusgens, B.V. Altered expression of angiogenesis and lymphangiogenesis markers in the uninvolved skin of plaque-type psoriasis. Br. J. Dermatol. 2009, 160, 581–590. [Google Scholar] [CrossRef] [PubMed]
  36. Henno, A.; Blacher, S.; Lambert, C.A.; Deroanne, C.; Noel, A.; Lapiere, C.; de la Brassinne, M.; Nusgens, B.V.; Colige, A. Histological and transcriptional study of angiogenesis and lymphangiogenesis in uninvolved skin, acute pinpoint lesions and established psoriasis plaques: An approach of vascular development chronology in psoriasis. J. Dermatol. Sci. 2010, 57, 162–169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Nickoloff, B.J.; Schroder, J.M.; von den Driesch, P.; Raychaudhuri, S.P.; Farber, E.M.; Boehncke, W.H.; Morhenn, V.B.; Rosenberg, E.W.; Schon, M.P.; Holick, M.F. Is psoriasis a T-cell disease? Exp. Dermatol. 2000, 9, 359–375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Krueger, J.G.; Krane, J.F.; Carter, D.M.; Gottlieb, A.B. Role of growth factors, cytokines, and their receptors in the pathogenesis of psoriasis. J. Investig. Dermatol. 1990, 94 (Suppl. 6), 135s–140s. [Google Scholar] [CrossRef] [PubMed]
  39. McColl, A.; Thomson, C.A.; Nerurkar, L.; Graham, G.J.; Cavanagh, J. TLR7-mediated skin inflammation remotely triggers chemokine expression and leukocyte accumulation in the brain. J. Neuroinflamm. 2016, 13, 102. [Google Scholar] [CrossRef] [PubMed]
  40. Riis, J.L.; Johansen, C.; Vestergaard, C.; Bech, R.; Kragballe, K.; Iversen, L. Kinetics and differential expression of the skin-related chemokines CCL27 and CCL17 in psoriasis, atopic dermatitis and allergic contact dermatitis. Exp. Dermatol. 2011, 20, 789–794. [Google Scholar] [CrossRef] [PubMed]
  41. Tian, S.; Krueger, J.G.; Li, K.; Jabbari, A.; Brodmerkel, C.; Lowes, M.A.; Suarez-Farinas, M. Meta-analysis derived (MAD) transcriptome of psoriasis defines the “core” pathogenesis of disease. PLoS ONE 2012, 7, e44274. [Google Scholar] [CrossRef] [PubMed]
  42. Lowes, M.A.; Russell, C.B.; Martin, D.A.; Towne, J.E.; Krueger, J.G. The IL-23/T17 pathogenic axis in psoriasis is amplified by keratinocyte responses. Trends Immunol. 2013, 34, 174–181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Swindell, W.R.; Johnston, A.; Voorhees, J.J.; Elder, J.T.; Gudjonsson, J.E. Dissecting the psoriasis transcriptome: Inflammatory- and cytokine-driven gene expression in lesions from 163 patients. Genomics 2013, 14, 527. [Google Scholar] [CrossRef] [PubMed]
  44. Motta, S.; Monti, M.; Sesana, S.; Caputo, R.; Carelli, S.; Ghidoni, R. Ceramide composition of the psoriatic scale. Biochim. Biophys. Acta 1993, 1182, 147–151. [Google Scholar] [CrossRef]
  45. Van Smeden, J.; Janssens, M.; Gooris, G.S.; Bouwstra, J.A. The important role of stratum corneum lipids for the cutaneous barrier function. Biochim. Biophys. Acta 2014, 1841, 295–313. [Google Scholar] [CrossRef] [PubMed]
  46. Farwanah, H.; Raith, K.; Neubert, R.H.; Wohlrab, J. Ceramide profiles of the uninvolved skin in atopic dermatitis and psoriasis are comparable to those of healthy skin. Arch Dermatol. Res. 2005, 296, 514–521. [Google Scholar] [CrossRef] [PubMed]
  47. Murakami, M.; Kudo, I. Phospholipase A2. J. Biochem. 2002, 131, 285–292. [Google Scholar] [CrossRef] [PubMed]
  48. Ye, X.; Ishii, I.; Kingsbury, M.A.; Chun, J. Lysophosphatidic acid as a novel cell survival/apoptotic factor. Biochim. Biophys. Acta 2002, 1585, 108–113. [Google Scholar] [CrossRef]
  49. Capestrano, M.; Mariggio, S.; Perinetti, G.; Egorova, A.V.; Iacobacci, S.; Santoro, M.; Di Pentima, A.; Iurisci, C.; Egorov, M.V.; Di Tullio, G.; et al. Cytosolic phospholipase A(2)epsilon drives recycling through the clathrin-independent endocytic route. J. Cell Sci. 2014, 127, 977–993. [Google Scholar] [CrossRef] [PubMed]
  50. Stewart, A.; Ghosh, M.; Spencer, D.M.; Leslie, C.C. Enzymatic properties of human cytosolic phospholipase A(2)gamma. J. Biol. Chem. 2002, 277, 29526–29536. [Google Scholar] [CrossRef] [PubMed]
  51. Ghosh, M.; Tucker, D.E.; Burchett, S.A.; Leslie, C.C. Properties of the Group IV phospholipase A2 family. Prog. Lipid Res. 2006, 45, 487–510. [Google Scholar] [CrossRef] [PubMed]
  52. Hong, K.K.; Cho, H.R.; Ju, W.C.; Cho, Y.; Kim, N.L. A Study on Altered Expression of Serine Palmitoyltransferase and Ceramidase in Psoriatic Skin Lesion. J. Korean Med. Sci. 2007, 22, 862–867. [Google Scholar] [CrossRef] [PubMed]
  53. Pietrzak, A.; Michalak-Stoma, A.; Chodorowska, G.; Szepietowski, J.C. Lipid disturbances in psoriasis: An update. Mediat. Inflamm. 2010, 2010, 535612. [Google Scholar] [CrossRef] [PubMed]
  54. National Center for Biotechnology Information. Available online: http://www.ncbi.nlm.nih.gov/gene/667 (accessed on 13 March 2018).
  55. Sonnenberg, A.; Liem, R.K. Plakins in development and disease. Exp. Cell Res. 2007, 313, 2189–2203. [Google Scholar] [CrossRef] [PubMed]
  56. Michael, M.; Begum, R.; Fong, K.; Pourreyron, C.; South, A.P.; McGrath, J.A.; Parsons, M. BPAG1-e restricts keratinocyte migration through control of adhesion stability. J. Investig. Dermatol. 2014, 134, 773–782. [Google Scholar] [CrossRef] [PubMed]
  57. Ferrier, A.; Boyer, J.G.; Kothary, R. Cellular and Molecular Biology of Neuronal Dystonin. In International Review of Cell and Molecular Biology; Jeon, K.W., Ed.; Elsevier: New York, NY, USA, 2013; pp. 88–95. [Google Scholar]
  58. Litjens, S.H.; Koster, J.; Kuikman, I.; van Wilpe, S.; de Pereda, J.M.; Sonnenberg, A. Specificity of binding of the plectin actin-binding domain to beta 4 integrin. Mol. Biol. Cell 2003, 14, 4039–4050. [Google Scholar] [CrossRef] [PubMed]
  59. Ishida, S.; Takahashi, K.; Kanaoka, M.; Okawa, T.; Tateishi, C.; Yasukochi, A.; Ishii, N.; Li, X.; Hashimoto, T.; Aihara, M. Case of subepidermal autoimmune bullous disease with psoriasis vulgaris reacting to both BP180 C-terminal domain and laminin gamma-1. J. Dermatol. 2015, 42, 391–393. [Google Scholar] [CrossRef] [PubMed]
  60. Stoica, L.E.; Patrascu, V.; Dascalu, R.C.; Ciurea, M.E. Bullous pemphigoid associated with psoriasis, breast cancer and Parkinson’s disease. Curr. Health Sci. J. 2014, 40, 62–66. [Google Scholar] [PubMed]
  61. Nakayama, C.; Iwata, H.; Haga, N.; Hamade, Y.; Mizuno, O.; Nishie, W.; Shimizu, H. The different intensity of autoantibody deposits in bullous pemphigoid associated with psoriasis vulgaris. Eur. J. Dermatol. 2014, 25, 70–71. [Google Scholar]
  62. Bergboer, J.G.; Dulak, M.G.; van Vlijmen-Willems, I.M.; Jonca, N.; van Wijk, E.; Hendriks, W.J.; Zeeuwen, P.L.; Schalkwijk, J. Analysis of protein-protein interaction between late cornified envelope proteins and corneodesmosin. Exp. Dermatol. 2014, 23, 769–771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. De Cid, R.; Riveira-Munoz, E.; Zeeuwen, P.L.; Robarge, J.; Liao, W.; Dannhauser, E.N.; Giardina, E.; Stuart, P.E.; Nair, R.; Helms, C.; et al. Deletion of the late cornified envelope LCE3B and LCE3C genes as a susceptibility factor for psoriasis. Nat. Genet. 2009, 41, 211–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Hohl, D. Expression Patterns of Loricrin in Dermatological Disorders. Am. J. Dermatopathol. 1993, 15, 20–27. [Google Scholar] [CrossRef] [PubMed]
  65. Rashmi, R.; Rao, K.S.; Basavaraj, K.H. A comprehensive review of biomarkers in psoriasis. Clin. Exp. Dermatol. 2009, 34, 658–663. [Google Scholar] [CrossRef] [PubMed]
  66. Slominski, A.T.; Zmijewski, M.A.; Zbytek, B.; Tobin, D.J.; Theoharides, T.C.; Rivier, J. Key role of CRF in the skin stress response system. Endocr. Rev. 2013, 34, 827–884. [Google Scholar] [CrossRef] [PubMed]
  67. Slominski, A.T.; Zmijewski, M.A.; Skobowiat, C.; Zbytek, B.; Slominski, R.M.; Steketee, J.D. Sensing the environment: Regulation of local and global homeostasis by the skin’s neuroendocrine system. Adv. Anat. Embryol. Cell Biol. 2012, 212, 1–115. [Google Scholar]
  68. Slominski, A.T.; Manna, P.R.; Tuckey, R.C. On the role of skin in the regulation of local and systemic steroidogenic activities. Steroids 2015, 103, 72–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Slominski, A.T.; Zmijewski, M.A.; Plonka, P.M.; Szaflarski, J.P.; Paus, R. How UV Light Touches the Brain and Endocrine System Through Skin, and Why. Endocrinology 2018, 159, 1992–2007. [Google Scholar] [CrossRef] [PubMed]
  70. Skobowiat, C.; Postlethwaite, A.E.; Slominski, A.T. Skin Exposure to Ultraviolet B Rapidly Activates Systemic Neuroendocrine and Immunosuppressive Responses. Photochem. Photobiol. 2017, 93, 1008–1015. [Google Scholar] [CrossRef] [PubMed]
  71. Slominski, A.T.; Brozyna, A.A.; Tuckey, R.C. Cutaneous Glucocorticoidogenesis and Cortisol Signaling Are Defective in Psoriasis. J. Investig. Dermatol. 2017, 137, 1609–1611. [Google Scholar] [CrossRef] [PubMed]
  72. Slominski, A.; Wortsman, J.; Luger, T.; Paus, R.; Solomon, S. Corticotropin releasing hormone and proopiomelanocortin involvement in the cutaneous response to stress. Physiol. Rev. 2000, 80, 979–1020. [Google Scholar] [CrossRef] [PubMed]
  73. Slominski, A.; Zbytek, B.; Szczesniewski, A.; Wortsman, J. Cultured human dermal fibroblasts do produce cortisol. J. Investig. Dermatol. 2006, 126, 1177–1178. [Google Scholar] [CrossRef] [PubMed]
  74. Hannen, R.F.; Michael, A.E.; Jaulim, A.; Bhogal, R.; Burrin, J.M.; Philpott, M.P. Steroid synthesis by primary human keratinocytes; implications for skin disease. Biochem. Biophys. Res. Commun. 2011, 404, 62–67. [Google Scholar] [CrossRef] [PubMed]
  75. Slominski, A.; Zbytek, B.; Semak, I.; Sweatman, T.; Wortsman, J. CRH stimulates POMC activity and corticosterone production in dermal fibroblasts. J. Neuroimmunol. 2005, 162, 97–102. [Google Scholar] [CrossRef] [PubMed]
  76. Slominski, A.; Zbytek, B.; Szczesniewski, A.; Semak, I.; Kaminski, J.; Sweatman, T.; Wortsman, J. CRH stimulation of corticosteroids production in melanocytes is mediated by ACTH. Am. J. Physiol. Endocrinol. Metab. 2005, 288, E701–E706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Vukelic, S.; Stojadinovic, O.; Pastar, I.; Rabach, M.; Krzyzanowska, A.; Lebrun, E.; Davis, S.C.; Resnik, S.; Brem, H.; Tomic-Canic, M. Cortisol synthesis in epidermis is induced by IL-1 and tissue injury. J. Biol. Chem. 2011, 286, 10265–10275. [Google Scholar] [CrossRef] [PubMed]
  78. Jean, J.; Garci-Perez, M.E.; Pouliot, R. Bioengineered Skin: The Self-Assembly Approach. J. Tissue Sci. Eng. 2011, 3, 001. [Google Scholar] [CrossRef]
  79. Pouliot, R.; Larouche, D.; Auger, F.A.; Juhasz, J.; Xu, W.; Li, H.; Germain, L. Reconstructed human skin produced in vitro and grafted on athymic mice. Transplantation 2002, 73, 1751–1757. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02923 g001 550
Figure 1. Macroscopic analysis of the reconstructed skin substitutes. For each group, tissue-engineered skin substitutes were produced with three different cell populations. (A–C) Tissue-engineered skin substitutes produced with healthy fibroblasts and keratinocytes. (D–I) Tissue-engineered skin substitutes produced with fibroblasts and keratinocytes isolated from either non-lesional (D–F) or lesional (G–I) psoriatic skin. Scale bar = 1 cm. Black arrowheads indicate the position of protuberant regions, whereas white arrowheads position thinner regions of the reconstructed skin substitutes.
Figure 1. Macroscopic analysis of the reconstructed skin substitutes. For each group, tissue-engineered skin substitutes were produced with three different cell populations. (A–C) Tissue-engineered skin substitutes produced with healthy fibroblasts and keratinocytes. (D–I) Tissue-engineered skin substitutes produced with fibroblasts and keratinocytes isolated from either non-lesional (D–F) or lesional (G–I) psoriatic skin. Scale bar = 1 cm. Black arrowheads indicate the position of protuberant regions, whereas white arrowheads position thinner regions of the reconstructed skin substitutes.
Ijms 19 02923 g001
Ijms 19 02923 g002 550
Figure 2. Histological analysis of the reconstructed skin substitutes. Mason’s trichrome staining after 21 days of culture at the air-liquid interface. For each group, tissue-engineered skin substitutes were produced with three different cell populations. (A–C) Healthy skin substitutes. (D–F) Non-lesional psoriatic skin substitutes. (G–I) Lesional psoriatic skin substitutes. Scale bar = 100 µm. C: Stratum corneum; E: Epidermis; D: Dermis.
Figure 2. Histological analysis of the reconstructed skin substitutes. Mason’s trichrome staining after 21 days of culture at the air-liquid interface. For each group, tissue-engineered skin substitutes were produced with three different cell populations. (A–C) Healthy skin substitutes. (D–F) Non-lesional psoriatic skin substitutes. (G–I) Lesional psoriatic skin substitutes. Scale bar = 100 µm. C: Stratum corneum; E: Epidermis; D: Dermis.
Ijms 19 02923 g002
Ijms 19 02923 g003 550
Figure 3. Microarray analysis. (A) Scatter plot of log2 of signal intensity from 60,000 different targets covering the entire human transcriptome of healthy or non-lesional tissue-engineered skin substitutes (in the y-axis) against non-lesional (first graph) or lesional (last two graphs) (in the x-axis). (B) Heatmap representation of genes whose expression is differentially regulated by at least two-fold in healthy against lesional substitutes and in non-lesional against lesional substitutes. (C) Venn diagram that indicates the number of deregulated genes in healthy against lesional substitutes (red circle) and in non-lesional against lesional substitutes (green circle). Genes that are commonly deregulated between these two groups are indicated in yellow. (D) Heatmap representation of the 55 most deregulated genes expressed by lesional substitutes relative to healthy substitutes (left) and non-lesional substitutes against lesional substitutes (right). The most highly expressed genes are shown in red, while the most weakly expressed are in blue. The genes with red writing are those identified as similarly deregulated between the two heatmaps. H: healthy; L: lesional; NL: non-lesional.
Figure 3. Microarray analysis. (A) Scatter plot of log2 of signal intensity from 60,000 different targets covering the entire human transcriptome of healthy or non-lesional tissue-engineered skin substitutes (in the y-axis) against non-lesional (first graph) or lesional (last two graphs) (in the x-axis). (B) Heatmap representation of genes whose expression is differentially regulated by at least two-fold in healthy against lesional substitutes and in non-lesional against lesional substitutes. (C) Venn diagram that indicates the number of deregulated genes in healthy against lesional substitutes (red circle) and in non-lesional against lesional substitutes (green circle). Genes that are commonly deregulated between these two groups are indicated in yellow. (D) Heatmap representation of the 55 most deregulated genes expressed by lesional substitutes relative to healthy substitutes (left) and non-lesional substitutes against lesional substitutes (right). The most highly expressed genes are shown in red, while the most weakly expressed are in blue. The genes with red writing are those identified as similarly deregulated between the two heatmaps. H: healthy; L: lesional; NL: non-lesional.
Ijms 19 02923 g003
Table
Table 1. Deregulated cytokines and growth factors.
Table 1. Deregulated cytokines and growth factors.
Gene SymbolGene NameLinear SignalFold Change (L/H)
Healthy (H)Lesional (L)
IL1R2Interleukin-1 receptor type 2, soluble form21.11152.4872.486 up
CXCL13C-X-C motif chemokine 138.95660.8446.793 up
CXCL14C-X-C motif chemokine 145475.35526,650.3084.867 up
CAMPCathelicidin antimicrobial peptide26.90465.1272.416 up
EPOErythropoietin23.08194.424.090 up
INHBAInhibin beta A chain132.158293.072.217 up
TNFSF9Tumor necrosis factor ligand superfamily member 9110.858332.4362.998 up
CXCL2C-X-C motif chemokine 2334.366158.7472.106 down
CCL20C-C motif chemokine 20135.67148.32.808 down
CXCL1Growth-regulated alpha protein1739.554656.3182.650 down
LTBLymphotoxin-beta57.05116.1613.530 down
TNFRSF10ATumor necrosis factor receptor superfamily member 10A179.23284.8482.112 down
IL24Interleukin-24164.62567.682.432 down
CCL27C-C motif chemokine 274802.716731.6516.564 down
IL15Interleukin-15368.781174.5542.112 down
Table
Table 2. Gene ontology analysis.
Table 2. Gene ontology analysis.
PathwaysSample Frequency (n)p-Value
Healthy against lesional condition
Keratinization (GO: 0031424)43.625 × 10−4
Isoprenoid metabolic process (GO: 0006720)55.922 × 10−4
Retinoid metabolic process (GO: 0001523)42.321 × 10−3
Diterpenoid metabolic process (GO: 0016101)43.077 × 10−3
Biological process (GO: 0008150)504.007 × 10−3
Keratinocyte differentiation (GO: 0030216)45.312 × 10−3
Non-lesional against lesional condition
Skin development (GO: 0043588)131.151 × 10−11
Keratinization (GO: 0031424)57.952 × 10−6
Epidermis development (GO: 0008544)81.088 × 10−5
Single-organism developmental process (GO: 0044767)272.897 × 10−5
Organ development (GO: 0048513)202.984 × 10−5
Developmental process (GO: 0032502)273.619 × 10−5
Table
Table 3. Selection of genes whose expression is altered through the conditions.
Table 3. Selection of genes whose expression is altered through the conditions.
GeneProtein NameLinear SignalsFold Change
HNLLH vs. NLH vs. LNL vs. L
AREGAmphiregulin4556.1143649.7481740.8641.248 down2.617 down2.096 down
CCL27C-C motif chemokine 274802.7161196.792327.9914.012 down14.646 down3.649 down
CERS3Ceramide synthase 31205.2981076.517545.5221.119 down2.209 down1.973 down
COL10A1Collagen alpha-1(X) chain, Collagen type X alpha 140.35153.551223.5813.805 up5.540 up1.456 up
COL4A1Collagen alpha-1(IV) chain624.3841325.6451377.7162.123 up2.206 up1.039 up
COL4A2Collagen type IV alpha 2, Collagen alpha-2(IV) chain3008.0667008.37538.5322.329 up2.506 up1.075 up
COL5A3Collagen alpha-3(V) chain81.195201.437272.1052.480 up3.351 up1.350 up
COL6A3Collagen alpha-3(VI) chain, Uncharacterized protein216.945264.225666.8491.217 up3.073 up2.523 up
COL8A1Collagen alpha-1(VIII) chain179.866373.282419.7552.075 up2.333 up1.124 up
COL9A3Collagen alpha-3(IX) chain; Collagen, type IX, alpha 3672.004561.211177.5321.197 down3.785 down3.161 down
CXCL10C-X-C motif chemokine 1018.22316.4716.6641.106 down2.734 down2.471 down
DMRTA1Doublesex- and mab-3-related transcription factor A1149.17172.60929.6442.054 down5.032 down2.449 down
DSTDystonin2666.3232193.123654.8541.215 down4.071 down3.349 down
FABP6Gastrotropin396.131270.887105.9031.462 down3.740 down2.557 down
LPLLipoprotein lipase354.755560.3411138.3061.579 up3.208 up2.031 up
NOD2Nucleotide-binding oligomerization domain-containing protein 21002.491893.358445.8881.122 down2.248 down2.003 down
PIK3R2Phosphatidylinositol 3-kinase regulatory subunit beta2620.7555460.4446354.3032.083 up2.424 up1.163 up
PLA2G4CCytosolic phospholipase A2 gamma287.49273.72473.5533.899 down3.908 down1.002 down
PLA2G4ECytosolic phospholipase A2 epsilon260.291205.14271.1761.268 down3.657 down2.882 down
PNPLA5Patatin-like phospholipase domain-containing protein 515.7881.34798.8065.155 up6.261 up1.214 up
POSTNPeriostin85.28478.011044.0075.605 up12.242 up2.184 up

Share and Cite

MDPI and ACS Style

Rioux, G.; Pouliot-Bérubé, C.; Simard, M.; Benhassine, M.; Soucy, J.; Guérin, S.L.; Pouliot, R. The Tissue-Engineered Human Psoriatic Skin Substitute: A Valuable In Vitro Model to Identify Genes with Altered Expression in Lesional Psoriasis. Int. J. Mol. Sci. 2018, 19, 2923. https://doi.org/10.3390/ijms19102923

AMA Style

Rioux G, Pouliot-Bérubé C, Simard M, Benhassine M, Soucy J, Guérin SL, Pouliot R. The Tissue-Engineered Human Psoriatic Skin Substitute: A Valuable In Vitro Model to Identify Genes with Altered Expression in Lesional Psoriasis. International Journal of Molecular Sciences. 2018; 19(10):2923. https://doi.org/10.3390/ijms19102923

Chicago/Turabian Style

Rioux, Geneviève, Claudia Pouliot-Bérubé, Mélissa Simard, Manel Benhassine, Jacques Soucy, Sylvain L. Guérin, and Roxane Pouliot. 2018. "The Tissue-Engineered Human Psoriatic Skin Substitute: A Valuable In Vitro Model to Identify Genes with Altered Expression in Lesional Psoriasis" International Journal of Molecular Sciences 19, no. 10: 2923. https://doi.org/10.3390/ijms19102923

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop