Literature DB >> 34909750

Deletion of TNFAIP6 Gene in Human Keratinocytes Demonstrates a Role for TSG-6 to Retain Hyaluronan Inside Epidermis.

Céline Evrard1, Emilie Faway1, Evelyne De Vuyst1, Olivier Svensek1, Valérie De Glas1, David Bergerat2, Michel Salmon3, Olivier De Backer1, Bruno Flamion1, Hélène Le-Buanec4, Catherine Lambert de Rouvroit1, Yves Poumay1.   

Abstract

TSG-6 is a soluble protein secreted in the extracellular matrix by various cell types in response to inflammatory stimuli. TSG-6 interacts with extracellular matrix molecules, particularly hyaluronan (HA), and promotes cutaneous wound closure in mice. Between epidermal cells, the discrete extracellular matrix contains HA and a tiny amount of TSG-6. However, challenges imposed to keratinocytes in reconstructed human epidermis revealed strong induction of TSG-6 expression, after exposure to T helper type 2 cytokines to recapitulate the atopic dermatitis phenotype or after fungal infection that causes secretion of cytokines and antimicrobial peptides. After both types of challenge, enhanced release of TSG-6 happens simultaneously with increased HA production. TSG-6 deficiency in N/TERT keratinocytes was created by inactivating TNFAIP6 using CRISPR/Cas9. Some TSG-6 -/- keratinocytes analyzed through scratch assays tend to migrate more slowly but produce reconstructed human epidermis that exhibits normal morphology and differentiation. Few significant alterations were noticed by transcriptomic analysis. Nevertheless, reduced HA content in TSG-6 -/- reconstructed human epidermis was observed, along with enhanced HA release into the culture medium, and this phenotype was even more pronounced after the challenging conditions. Reintroduction of cells producing TSG-6 in reconstructed human epidermis reduced HA leakage. Our results show a role for TSG-6 in sequestering HA between epidermal cells in response to inflammation.
© 2021 The Authors.

Entities:  

Keywords:  AD, atopic dermatitis; ECM, extracellular matrix; GEO, Gene Expression Omnibus; HA, hyaluronan; HAS, hyaluronan synthase; KC, keratinocyte; KLK, kallikrein; RHE, reconstructed human epidermis; RNA-seq, RNA sequencing; Th, T helper type; crRNA, CRISPR RNA

Year:  2021        PMID: 34909750      PMCID: PMC8659394          DOI: 10.1016/j.xjidi.2021.100054

Source DB:  PubMed          Journal:  JID Innov        ISSN: 2667-0267


Introduction

In the skin, although occupying most of the space in the dermis, the extracellular matrix (ECM) also plays crucial roles in the epidermis, where it is composed of basal lamina and 15–20 nm intercellular spaces around cell junctions. The main ECM components in the epidermis include various polysaccharides, particularly hyaluronan (HA), which is swollen by water. This environment creates a favorable milieu for the diffusion of nutrients, GFs, and cytokines (Sandjeu and Haftek, 2009; Toole, 2001). Langerhans cells and lymphocytes migrate through the epidermal ECM, bringing innate or adaptive immunity. HA is a linear glycosaminoglycan composed of repeated disaccharides, namely D-glucuronic acid, linked to N-acetyl-D-glucosamine, forming polymers up to 25 μm in length and 105–107 Da in mass (Sakai et al., 2000; Toole, 2004). Half of the total body HA resides in the skin (Fraser et al., 1997; Meyer and Stern, 1994), mostly in the dermis but also in the living layers of the epidermis (Pasonen-Seppänen et al., 2003). The functions of epidermal HA are not fully understood. Roles in proliferation and differentiation of epidermal keratinocytes (KCs) are unlikely, but altered regulation of HA synthases (HASs) and increased HA production were described in inflammatory conditions such as atopic dermatitis (AD) (Malaisse et al., 2014). Across the body, HA interacts with several partners, such as TSG-6 protein, secreted by fibroblasts on TNF-α stimulation (Jadin et al., 2014; Lee et al., 1990). TSG-6 is encoded by the TNFAIP6 gene located on human chromosome 2. The protein of 277 amino acids harbors a cleavable signal peptide, an HA-binding (LINK) domain typical of hyaladherins, and a CUB domain, typically part of ECM and plasma membrane proteins. Besides fibroblasts, several cell types secrete TSG-6 protein in response to proinflammatory stimuli (Baranova et al., 2011). Initially associated with inflammation, TSG-6 is also detected in conditions implying ECM remodeling (Fülöp et al., 2003) and in organs exerting a barrier function (Day and Milner, 2019; Tan et al., 2011). In healthy skin, TSG-6 is detected in the dermis and in living epidermal layers (Tan et al., 2011). Its expression is upregulated during skin injury (Shakya et al., 2020). Soluble TSG-6 protein dimerizes or interacts with ITI, and thereby becomes able to bind HA’s N-acetyl glucosamine unit and favor cross-linking of HA chains in ECM (Baranova et al., 2011; Milner and Day, 2003; Milner et al., 2007). Interaction between TSG-6 and HA increases the affinity of HA for the CD44 receptor involved in KCs adhesion, proliferation, differentiation, and migration (Baranova et al., 2013; Bourguignon, 2014; Bourguignon et al., 2004; Kavasi et al., 2017; Lesley et al., 2004). CD44 also contributes to maintaining a pericellular matrix around KCs (Muto et al., 2019; Pasonen-Seppänen et al., 2012). Similar to CD44, TSG-6 holds HA in the pericellular matrix (Selbi et al., 2006), thereby controlling differentiation of, for example, fibroblasts (Simpson et al., 2009). In this study, we investigated the human epidermis produced and challenged in vitro (Hennies and Poumay, 2021) to characterize the TSG-6 expression induced in such conditions. We then used genetic approaches to create TSG-6–/– in vitro reconstructed epidermis and analyze the potential role for TSG-6. Data show a crucial impact of TSG-6 deficiency on epidermal HA.

Results

TSG-6 is expressed and released at low levels in the culture medium of normal human epidermis

Reconstructed human epidermis (RHE) produced from primary KCs (Frankart et al., 2012) was analyzed for TSG-6-encoding mRNA expression and protein release in a culture medium during the whole reconstruction process (Figure 1). RT-qPCR showed that TSG-6-encoding mRNA expression remained very low (Cq > 32). No significant induction of TSG-6-encoding mRNA was observed (Figure 1a). Accordingly, in situ hybridization was unable to detect TSG-6-encoding mRNA in normal RHE and human epidermis (Figure 1b). Secreted TSG-6 protein levels were monitored, and small amounts were detected in the culture medium during the reconstruction process (Figure 1c).
Figure 1

Weak TSG-6-encoding mRNA expression and protein release in a culture medium during the reconstruction of the normal epidermis. (a) TSG-6-encoding mRNA expression of RHE at 1, 3, 5, 7, 9, and 11 days of growth assessed by RT-qPCR. Reference gene is RPLP0 (mean ± SD; n = 3; one-way ANOVA). (b) TSG-6, PPIB (CTL+), and DapB (CTL–) mRNA detection using in situ hybridization and hemalun counterstaining of day 11 RHE and normal human epidermis. Bar = 20 μm. (c) TSG-6 protein concentration assayed using ELISA (mean ± SD; n = 4; one-way ANOVA). CTL, control; RHE, reconstructed human epidermis.

Weak TSG-6-encoding mRNA expression and protein release in a culture medium during the reconstruction of the normal epidermis. (a) TSG-6-encoding mRNA expression of RHE at 1, 3, 5, 7, 9, and 11 days of growth assessed by RT-qPCR. Reference gene is RPLP0 (mean ± SD; n = 3; one-way ANOVA). (b) TSG-6, PPIB (CTL+), and DapB (CTL–) mRNA detection using in situ hybridization and hemalun counterstaining of day 11 RHE and normal human epidermis. Bar = 20 μm. (c) TSG-6 protein concentration assayed using ELISA (mean ± SD; n = 4; one-way ANOVA). CTL, control; RHE, reconstructed human epidermis.

TSG-6 expression and release strongly increase on exposure of RHE to challenging conditions

To monitor TSG-6 regulation in inflammatory settings, two very different challenging conditions were imposed on RHE to analyze the tissue response. A first challenge, mimicking AD by deprivation of plasma membrane cholesterol followed by incubation with T helper type (Th) 2 ILs and exhibiting barrier alterations (De Vuyst et al., 2018) were selected on the basis of previous transcriptomic analysis (Hubaux et al., 2018; Park et al., 2021). In response to IL-4, IL-13, and IL-25 incubation for 48 hours, TNFAIP6 appeared as the sixth most upregulated gene, with a 27-fold increase in expression and a P-value of 1.37 × 10–8. Results are accessible at Gene Expression Omnibus (GEO) Series (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE149638; GSE149639). Unlike control tissues, AD-like RHE exhibited spongiosis and hypogranulosis (Figure 2a) and brighter detection of HA around KCs (Figure 2b). Analysis of TSG-6-encoding mRNA expression revealed an approximately 90-fold induction in AD-like RHE compared with that of the controls (Figure 2c). The amount of TSG-6 protein released in the culture medium corroborated increased epidermal expression (Figure 2d). TSG-6-encoding mRNA expression in RHE incubated for 48 hours with either IL-4, IL-13, or IL-25 or with methyl-β-cyclodextrin identified IL-4 as the stronger inducer of TSG-6 expression, whereas IL-25 and methyl-β-cyclodextrin had no effect (Figure 2e). Thus, IL-4 and IL-13 were chosen to stimulate TSG-6 expression. In accordance with RT-qPCR data, TSG-6-encoding mRNA were detected using in situ hybridization in living layers of RHE challenged with IL-4 and IL-13 (Figure 2f). Conversely, addition of TNF-α in the culture medium of RHE did not induce TSG-6 expression (unpublished observations). In skin biopsies from lesional areas of patients with AD and from healthy controls, in situ hybridization was not sensitive enough to detect TSG-6-encoding mRNA (Figure 3a). However, RT-qPCR analysis of biopsies revealed increased TSG-6-encoding mRNA levels in lesional AD compared with that in the normal skin with a P-value as low as 0.0506 (Figure 3b). This increase in TSG-6 expression in the skin of patients suffering from AD is also detected in large-scale analyses. A search of GEO Profiles databases (Edgar et al., 2002) highlighted an increase in TSG-6 expression (GEO accession GSE5667 [Plager et al., 2010]) (Figure 3c). Likewise, transcriptomic analysis of a single-cell RNA sequencing (RNA-seq) dataset allowed a precise study of mRNA expression in undifferentiated and differentiated KCs of healthy and lesional AD skins (Reynolds et al., 2021; https://developmentcellatlas.ncl.ac.uk/datasets/ hca_skin_portal). As expected, expression of CA2 and NELL2, two typical hallmarks of AD lesions, was detected in more AD KCs than in normal KCs. In parallel, TSG-6-encoding mRNA was present in some undifferentiated AD KCs, whereas the signal was barely detectable in healthy cells (Figure 3d).
Figure 2

Increase of TSG-6-encoding mRNA and protein in challenged RHE (AD-like RHE model). RHEs were treated with MβCD and IL-4, IL-13, and IL-25 for 48 hours. (a) HE staining and (b) HA detection. Bar = 20 μm. (c) TSG-6-encoding mRNA expression of AD-like RHE assessed by RT-qPCR (mean ± SD; n = 3; Student t-test; ∗∗P < 0.01). (d) TSG-6 protein concentration of challenged RHE assayed using ELISA (mean ± SD; n = 3; Student t-test; ∗P < 0.05). (e) TSG-6-encoding mRNA expression of RHE treated for 48 hours by ILs or MβCD separately assessed by RT-qPCR (mean ± SD; n = 3; one-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. normal RHE). (f) TSG-6 and PPIB (CTL+) mRNA detection using in situ hybridization and hemalun counterstaining of RHE challenged or not challenged with IL-4 and IL-13. AD, atopic dermatitis; HA, hyaluronan; HE, hemalun/eosin; MβCD, methyl-β-cyclodextrin; RHE, reconstructed human epidermis.

Figure 3

Increase of TSG-6-encoding mRNA in lesional AD skin. (a) TSG-6 and PPIB (CTL+) mRNA detection using in situ hybridization and hemalun counterstaining of healthy and lesional AD skin. Bar = 20 μm. (b) TSG-6-encoding mRNA expression in normal and lesional AD skin assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 5 for normal skin and n = 10 for AD skin; unpaired t-test with Welch’s correction). (c) TSG-6-encoding mRNA expression in normal, nonlesional AD, and lesional AD skins analyzed from GEO Profiles database (accession GSE5667; GPL96; 206026_s_at) (mean ± SD; n = 5 or 6; one-way ANOVA). (d) UMAP visualization of undifferentiated and differentiated KC in healthy and lesional AD skins found in single-cell RNA-seq dataset (https://developmentcellatlas.ncl.ac.uk/datasets/ hca_skin_portal) (n = 5). Normalized gene expressions of K14, K10, CA2, NELL2, and TNFAIP6 are represented in green for a low expression and in red for a strong expression. AD, atopic dermatitis; CTL, control; GEO, Gene Expression Omnibus; K, keratin; KC, keratinocyte; RNA-seq, RNA sequencing; UMAP, Uniform Manifold Approximation and Projection.

Increase of TSG-6-encoding mRNA and protein in challenged RHE (AD-like RHE model). RHEs were treated with MβCD and IL-4, IL-13, and IL-25 for 48 hours. (a) HE staining and (b) HA detection. Bar = 20 μm. (c) TSG-6-encoding mRNA expression of AD-like RHE assessed by RT-qPCR (mean ± SD; n = 3; Student t-test; ∗∗P < 0.01). (d) TSG-6 protein concentration of challenged RHE assayed using ELISA (mean ± SD; n = 3; Student t-test; ∗P < 0.05). (e) TSG-6-encoding mRNA expression of RHE treated for 48 hours by ILs or MβCD separately assessed by RT-qPCR (mean ± SD; n = 3; one-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. normal RHE). (f) TSG-6 and PPIB (CTL+) mRNA detection using in situ hybridization and hemalun counterstaining of RHE challenged or not challenged with IL-4 and IL-13. AD, atopic dermatitis; HA, hyaluronan; HE, hemalun/eosin; MβCD, methyl-β-cyclodextrin; RHE, reconstructed human epidermis. Increase of TSG-6-encoding mRNA in lesional AD skin. (a) TSG-6 and PPIB (CTL+) mRNA detection using in situ hybridization and hemalun counterstaining of healthy and lesional AD skin. Bar = 20 μm. (b) TSG-6-encoding mRNA expression in normal and lesional AD skin assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 5 for normal skin and n = 10 for AD skin; unpaired t-test with Welch’s correction). (c) TSG-6-encoding mRNA expression in normal, nonlesional AD, and lesional AD skins analyzed from GEO Profiles database (accession GSE5667; GPL96; 206026_s_at) (mean ± SD; n = 5 or 6; one-way ANOVA). (d) UMAP visualization of undifferentiated and differentiated KC in healthy and lesional AD skins found in single-cell RNA-seq dataset (https://developmentcellatlas.ncl.ac.uk/datasets/ hca_skin_portal) (n = 5). Normalized gene expressions of K14, K10, CA2, NELL2, and TNFAIP6 are represented in green for a low expression and in red for a strong expression. AD, atopic dermatitis; CTL, control; GEO, Gene Expression Omnibus; K, keratin; KC, keratinocyte; RNA-seq, RNA sequencing; UMAP, Uniform Manifold Approximation and Projection. A second challenge imposed on RHE consisted in superficial infection by Trichophyton rubrum, which induces cytokine release in concomitance with the significant induction of expression and release of TSG-6 protein (Faway et al., 2019, 2017). After 4 days of infection, RHE becomes invaded by fungal elements only in the stratum corneum (Figure 4a). At the same time, the signal of HA seems higher around the KCs of infected RHE than around the control (Figure 4b). In parallel, TSG-6-encoding mRNA can be detected through in situ hybridization solely in infected RHE (Figure 4c) when its mRNA expression is significantly increased (induction of 46-fold) according to RT-qPCR analysis. Alongside, we found an interesting significant increase in the mRNA encoding LL37, an antimicrobial peptide that blocks the proinflammatory HA‒toll-like receptor 4‒CD44 axis (Morioka et al., 2008) (Figure 4d).
Figure 4

Elevated TSG-6-encoding mRNA expression and protein release in RHE infected by RHEs were infected or not infected by T. rubrum arthroconidia and maintained in culture for 4 additional days. (a) PAS staining with α-amylase pretreatment and hemalun counterstaining and (b) HA detection. Bar = 20 μm. (c) TSG-6-encoding mRNA detection using in situ hybridization and hemalun counterstaining. Bar = 20 μm. (d) LL37 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; one-way ANOVA; ∗∗P < 0.01 vs. day 0). HA, hyaluronan; PAS, periodic-acid schiff; RHE, reconstructed human epidermis.

Elevated TSG-6-encoding mRNA expression and protein release in RHE infected by RHEs were infected or not infected by T. rubrum arthroconidia and maintained in culture for 4 additional days. (a) PAS staining with α-amylase pretreatment and hemalun counterstaining and (b) HA detection. Bar = 20 μm. (c) TSG-6-encoding mRNA detection using in situ hybridization and hemalun counterstaining. Bar = 20 μm. (d) LL37 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; one-way ANOVA; ∗∗P < 0.01 vs. day 0). HA, hyaluronan; PAS, periodic-acid schiff; RHE, reconstructed human epidermis.

Generation, validation, and characterization of TSG-6–/– N/TERT clones

To investigate the potential roles of epidermal TSG-6, we used CRISPR/Cas9 to disrupt the TNFAIP6 gene in N/TERT KCs. RNA-seq analysis of epidermis reconstructed from N/TERT cells and stimulated with IL-4 and IL-13 showed a transcriptional profile similar to that of AD-like models made of primary KCs (Table 1). RNA-seq data are accessible in ArrayExpress (http://www.ebi.ac.uk/arrayexpresse; E-MTAB-10095).
Table 1

Transcriptomic Analysis of Epidermis Reconstructed from N/TERT Cells and Stimulated or Not with IL-4 and IL-13 through RNA Sequencing

GeneNamelog2FCPadj
CISHCytokine inducible SH2 containing protein3.076840250.000E+00
ANO1Anoctamin 12.805917401.042E–219
SLC16A14Solute carrier family 16 member 144.475507068.627E–192
CA2Carbonic anhydrase 22.186545821.411E–188
CTSCCathepsin C1.544622455.680E–123
GLDCGlycine decarboxylase6.484258337.039E–122
SLC26A2Solute carrier family 26 member 21.421196757.896E–113
AC100801.22.675824844.386E–86
C1QTNF1C1q and TNF related 12.765586252.474E–85
TNFAIP6TNF alpha-induced protein 64.614100571.402E–80
BCL2L15BCL2 like 154.386938211.861E–75
DIO2Iodothyronine deiodinase 2–1.386209738.042E–75
ST6GAL1ST6 beta-galactoside alpha-2,6-sialyltransferase 11.949049865.983E–73
EGR3Early growth response 3–1.438240801.248E–70
KRT16Keratin 161.766154761.724E–67
HRH1Histamine receptor H12.045225441.800E–60
POSTNPeriostin2.965690639.221E–60
ETS2ETS proto-oncogene 2, transcription factor–0.979909963.121E–59
EGLN3Egl-9 family hypoxia inducible factor 31.902227187.334E–58
APOL6Apolipoprotein L61.111877341.488E–57
LBHLBH regulator Of WNT signaling pathway1.631856913.760E–57
MXRA5Matrix remodeling associated 5–1.756464507.512E–56
RGMBRepulsive guidance molecule BMP coreceptor B1.270979078.026E–56
NEK6NIMA related kinase 61.469070884.158E–55
C1orf68Chromosome 1 open reading frame 68–2.175358435.709E–54
NELL2Neural EGFL like 22.776698731.809E–52
PTGS2Prostaglandin-endoperoxide synthase 2–2.485319232.546E–51
GCNT3Glucosaminyl (N-acetyl) transferase 3, mucin type2.139743161.036E–50
MAP3K14Mitogen-activated protein kinase 141.622403561.088E–50
DUOX2Dual oxidase 22.121968563.027E–50
PADI3Peptidyl arginine deiminase 31.839770155.216E–50
CD59CD59 molecule (CD59 blood group)0.753275731.087E–49
KRT4Keratin 41.999771311.734E–49
BMP4Bone morphogenetic protein 42.536704451.151E–47
BTBD3BTB domain containing 30.780793741.162E–47
AREGAmphiregulin–1.474073427.747E–47
PDE6APhosphodiesterase 6A2.222773441.398E–46
DENND4ADENN domain containing 4A0.832935011.443E–44
CCN1Cellular communication network factor 1–1.814120092.467E–44
SORL1Sortilin related receptor 10.968621925.845E–44
BCL10BCL10 immune signaling adaptor–1.018148509.519E–44
RASGRP1RAS guanyl releasing protein 11.158676821.093E–43
CXCL2C-X-C motif chemokine ligand 2 (IL-2)–3.207123717.322E–43
HBEGFHeparin-binding EGF like growth factor–1.921969761.046E–42
FLT1Fms related receptor tyrosine kinase 12.523707481.357E–42
LGALSLGalectin like–0.935018393.556E–42
CASP14Caspase 14–1.201336516.805E–41
CTSHCathepsin H–1.485722798.798E–40
TXNThioredoxin1.094196412.633E–39

Abbreviations: Padj, adjusted P; FC, fold change; RHE, reconstructed human epidermis.

List of 50 more regulated genes in IL-4‒ and IL-13‒challenged RHE compared with nonstimulated RHE based on Padj-value.

Transcriptomic Analysis of Epidermis Reconstructed from N/TERT Cells and Stimulated or Not with IL-4 and IL-13 through RNA Sequencing Abbreviations: Padj, adjusted P; FC, fold change; RHE, reconstructed human epidermis. List of 50 more regulated genes in IL-4‒ and IL-13‒challenged RHE compared with nonstimulated RHE based on Padj-value. To generate deletions in the TNFAIP6 gene, five targets were chosen in the first, third, and sixth exons (Figure 5a). DNA extracted from the pool of N/TERT KCs electroporated with Cas9/ribonucleoprotein was analyzed by PCR to characterize deletions. Figure 5b illustrates the PCR data that revealed deletions at various sites. Among 100 clones sequenced and presenting many different mutations of several base pairs (not shown), two were selected for large deletions in both alleles of the TNFAIP6 gene leading to probable null alleles. Genomic positions were defined by numbering nucleotide A of the methionine initiation codon as 1. The first clone TNFAIP6UNΔ1/UNΔ2 named N/TERT TSG-6 (a) exhibits a UNΔ1 deletion in allele 1 that encompasses nucleotide ‒120 (i.e., 44 nucleotides upstream of exon 1) to +21,696 and a UNΔ2 deletion in allele 2 encompassing nucleotide +44 (in exon 1) to +21,635 (62 bases upstream from exon 6). The second clone TNFAIP6UNΔ3/UNΔ3 named N/TERT TSG-6 (b) exhibits identical deletions in both alleles encompassing nucleotide +33 (in exon 1) to +21,695 (two nucleotides upstream from exon 6) (Figure 5c). Both clones encode putative peptides lacking LINK and CUB domains.
Figure 5

Deletions in (a) Genomic organization of TNFAIP6 gene using the localization of five sites targeted by crRNA (red bars). (b) PCR screening of deletions in DNA from a pool of N/TERT cells electroporated with a mixture of five crRNA (red bars) and Cas9. Primers are indicated by arrowheads. Obtained amplicon sizes of WT and deletant alleles are shown on the right. (c) Genomic deletions in TSG-6–/– clone (a) and clone (b) using primer combinations PCR A and PCR E illustrated in panel b. All PCR products have been sequenced. Only the sequences of the junctions are shown. bp, base pair; crRNA, CRISPR RNA; KC, keratinocyte; WT, wild-type

Deletions in (a) Genomic organization of TNFAIP6 gene using the localization of five sites targeted by crRNA (red bars). (b) PCR screening of deletions in DNA from a pool of N/TERT cells electroporated with a mixture of five crRNA (red bars) and Cas9. Primers are indicated by arrowheads. Obtained amplicon sizes of WT and deletant alleles are shown on the right. (c) Genomic deletions in TSG-6–/– clone (a) and clone (b) using primer combinations PCR A and PCR E illustrated in panel b. All PCR products have been sequenced. Only the sequences of the junctions are shown. bp, base pair; crRNA, CRISPR RNA; KC, keratinocyte; WT, wild-type RHEs were produced using unedited TSG-6+/+ or deletants TSG-6 (a) and TSG-6 (b) N/TERT cells and were treated or not treated for 48 hours with IL-4 and IL-13. Both mutant TSG-6 (a) and TSG-6 (b) clones produced RHE with a functional barrier (not shown) and exhibited identical morphological features compared with unedited N/TERT cells, including localization of keratin 10, involucrin, loricrin, and FLG (Figure 6a). As expected, TSG-6-encoding mRNA was only overexpressed in the RHE produced with TSG-6+/+ N/TERT cells challenged with IL-4 and IL-13. Conversely, TSG-6-encoding mRNA was undetectable in the RHE produced with TSG-6 (a) or TSG-6 (b) N/TERT cells (Figure 6b). Using in situ hybridization, TSG-6-encoding mRNA were solely detected in living layers of TSG-6+/+ N/TERT RHE challenged with IL-4 and IL-13 and were absent in knockout RHE (Figure 6c). A strong increase in TSG-6 protein level was found in the medium of challenged TSG-6+/+ N/TERT but not in mutant RHE (Figure 6d). Transcriptomes of TSG-6 (a) and TSG-6 (b) RHE challenged with IL-4 and IL-13 were compared with those of TSG-6+/+ RHE (http://www.ebi.ac.uk/arrayexpresse; E-MTAB-10095), confirming suppressed TNFAIP6 expression. Deletion of TSG-6 resulted in weak upregulation and downregulation in a few genes (Table 2). Among them, suppressed expression of kallikreins (KLKs) 5 and 8, confirmed by RT-qPCR analysis in clone TSG-6–/– (b) (Figure 7), could indicate consequences on epidermal homeostasis and migration processes. Weak alterations in the gene expression observed after knocking out TSG-6 concurred with the normal morphology of TSG-6 RHE.
Figure 6

Validation and characterization of the two clones exhibiting large deletions in the TSG-6+/+ and TSG-6–/– (a and b) RHE stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HE staining and immunofluorescence (in green) of K10, IVL, LOR, and FLG. Bar = 20 μm. (b) TSG-6-encoding mRNA expression of TSG-6+/+ and TSG-6–/– RHE assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; two-way ANOVA; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (c) TSG-6-encoding mRNA detection using in situ hybridization and hemalun counterstaining of TSG-6+/+ and TSG-6–/– RHE. Bar = 20 μm. (d) TSG-6 protein concentration of TSG-6+/+ and TSG-6–/– RHE assayed using ELISA (mean ± SD; n = 3; two-way ANOVA; ∗∗∗P < 0.001 vs. TSG-6+/+). CTL, control; HE, hemalun/eosin; IVL, involucrin; K, keratin; LOR, loricrin; RHE, reconstructed human epidermis.

Table 2

Transcriptomic Analysis of TSG-6–/– Epidermis through RNA Sequencing



TSG-6–/– (a)
TSG-6–/– (b)
IDGenelog2FCPadjlog2FCPadj
ENSG00000123610TNFAIP6–7.415032.96E–32–7.327363E–30
ENSG00000234444ZNF736–4.446362.81E–22–4.693333.14E–18
ENSG00000170430MGMT–4.337128.06E–16–4.188493.18E–12
ENSG00000105750ZNF85–3.875943.01E–37–2.405979.47E–27
ENSG00000144115THNSL2–3.860522.91E–15–3.785591.02E–13
ENSG00000187210GCNT1–3.04386E–50–2.614446.09E–38
ENSG00000011465DCN–2.821651.1E–84–1.54913.33E–29
ENSG00000188257PLA2G2A–2.74611.93E–24–2.009171.94E–14
ENSG00000177990DPY19L2–2.371623.21E–10–2.424156.58E–09
ENSG00000153802TMPRSS11D–2.352789.53E–15–1.245440.00015
ENSG00000129455KLK8–2.276283.22E–18–1.709971.05E–10
ENSG00000122691TWIST1–2.263962.75E–06–2.355148.97E–08
ENSG00000103546SLC6A2–2.230511.73E–11–1.015280.002253
ENSG00000102837OLFM4–2.108581.07E–09–2.294763.6E–09
ENSG00000188868ZNF563–2.010961.75E–17–2.3947.67E–21
ENSG00000244242IFITM10–1.939996.41E–12–1.93491.62E–10
ENSG00000153292ADGRF1–1.936596.88E–08–1.085170.008469
ENSG00000206195DUXAP8–1.880833.76E–05–2.08416.63E–07
ENSG00000213096ZNF254–1.86031.41E–38–1.781.78E–33
ENSG00000260027HOXB7–1.720592.43E–05–2.214491.82E–05
ENSG00000081803CADPS2–1.649892.05E–05–1.675671.61E–05
ENSG00000237440AC008554.1–1.619360.000136–2.200675.76E–08
ENSG00000105357MYH14–1.489541.63E–16–1.00353.51E–07
ENSG00000101017CD40–1.477995.85E–06–2.084131.64E–08
ENSG00000172478MAB21L4–1.447839.05E–36–1.824771.8E–52
ENSG00000197647ZNF433–1.362085.15E–08–2.484963.31E–17
ENSG00000167754KLK5–1.347852.03E–16–1.922991.2E–30
ENSG00000198182ZNF607–1.3443.92E–17–1.094627.89E–11
ENSG00000250644AC068580.4–1.331368.21E–07–1.143640.000121
ENSG00000160284SPATC1L–1.275790.001851–1.651450.000263
ENSG00000196263ZNF471–1.207799.05E–05–1.206740.000435
ENSG00000077063CTTNBP2–1.165977.2E–06–2.105725.59E–12
ENSG00000235584AC008268.1–1.158262.67E–07–1.345381.76E–08
ENSG00000091656ZFHX4–1.143860.00455–1.276650.005149
ENSG00000196172ZNF681–1.132073.24E–16–1.082221.35E–13
ENSG00000150687PRSS23–1.045038.83E–14–1.572558.84E–28
ENSG00000139445FOXN4–1.031120.000567–1.483523.45E–06
ENSG00000219665ZNF433-AS1–1.01840.003964–1.370530.000331
ENSG00000130038CRACR2A4.1555044.42E–371.3396323.03E–05
ENSG00000007038PRSS212.8256154.82E–271.2299793.46E–05
ENSG00000248370LINC024342.611156.57E–431.6374528.7E–16
ENSG00000150722PPP1R1C2.3947061.72E–101.0899670.003993
ENSG00000182389CACNB42.1797911.85E–231.1750441.83E–06
ENSG00000131771PPP1R1B2.17771.76E–111.0579330.004363
ENSG00000141753IGFBP42.0836091.45E–421.2196613.22E–13
ENSG00000108176DNAJC122.0573472.89E–151.8009884.46E–11
ENSG00000234929AC018685.21.777641.76E–101.3251831.51E–05
ENSG00000166897ELFN21.6606263.69E–061.2685150.000785
ENSG00000153064BANK11.4817945.99E–081.3061591.21E–05
ENSG00000100968NFATC41.4717112.94E–111.0933999.51E–06
ENSG00000251003ZFPM2-AS11.3541863.59E–201.060912.67E–11
ENSG00000131634TMEM2041.3058411.52E–102.3889954.19E–33
ENSG00000164684ZNF7041.2153934.55E–091.0859851.47E–06
ENSG00000141622RNF1651.1622477.5E–051.0460220.001585
ENSG00000179914ITLN11.1440140.0001951.1045320.001204
ENSG00000279108AC008537.31.1147010.0010841.411620.000283
ENSG00000225968ELFN11.049660.0023121.3566170.000548
ENSG00000172005MAL1.006040.0008871.8814283.6E–10

Abbreviations: Padj, adjusted P; FC, fold change; RHE, reconstructed human epidermis.

List of commonly downregulated and upregulated genes in both IL-4‒ and IL-13‒challenged TSG-6–/– RHE compared with challenged TSG-6+/+ RHE (FC ≥ 2; P < 0.01).

Figure 7

mRNA levels of TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) KLK5 and (b) KLK8 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 4; two-way ANOVA; ∗P < 0.05 vs. TSG-6+/+). CTL, control; KLK, kallikrein; RHE, reconstructed human epidermis.

Validation and characterization of the two clones exhibiting large deletions in the TSG-6+/+ and TSG-6–/– (a and b) RHE stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HE staining and immunofluorescence (in green) of K10, IVL, LOR, and FLG. Bar = 20 μm. (b) TSG-6-encoding mRNA expression of TSG-6+/+ and TSG-6–/– RHE assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; two-way ANOVA; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (c) TSG-6-encoding mRNA detection using in situ hybridization and hemalun counterstaining of TSG-6+/+ and TSG-6–/– RHE. Bar = 20 μm. (d) TSG-6 protein concentration of TSG-6+/+ and TSG-6–/– RHE assayed using ELISA (mean ± SD; n = 3; two-way ANOVA; ∗∗∗P < 0.001 vs. TSG-6+/+). CTL, control; HE, hemalun/eosin; IVL, involucrin; K, keratin; LOR, loricrin; RHE, reconstructed human epidermis. Transcriptomic Analysis of TSG-6–/– Epidermis through RNA Sequencing Abbreviations: Padj, adjusted P; FC, fold change; RHE, reconstructed human epidermis. List of commonly downregulated and upregulated genes in both IL-4‒ and IL-13‒challenged TSG-6–/– RHE compared with challenged TSG-6+/+ RHE (FC ≥ 2; P < 0.01). mRNA levels of TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) KLK5 and (b) KLK8 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 4; two-way ANOVA; ∗P < 0.05 vs. TSG-6+/+). CTL, control; KLK, kallikrein; RHE, reconstructed human epidermis.

Deficiency in TSG-6 protein slows down KC migration

Because some genes implicated in migration (KLK5 and 8) were suppressed in TSG-6–/– RHE and because Tsg-6–/– mice show delayed wound closure (Shakya et al., 2020), the migrative ability of knocked-out TSG-6 KCs was monitored in a scratch assay. TSG-6 (a) N/TERT KCs migrated more slowly than TSG-6+/+ KCs, being unable to close the scratch within 24 hours. Migration of TSG-6 (b) N/TERT KCs was also slower than that of nonedited N/TERT cells (Figure 8a). The migration defect was partially restored in the presence of recombinant TSG-6 (Figure 8b).
Figure 8

Lack of TSG-6 slows down KCs migration. Monolayers cultures of TSG-6+/+ or TSG-6–/– (a and b) N/TERT KCs at confluence stage treated with IL-4 and IL-13 for 24 h. KCs were growth arrested by mitomycin C before performing a scratch. KCs migration was monitored for 24 h in the presence of IL-4 and IL-13. (a) Scratch closure was analyzed every 6 h by Cytonote and was quantified by ImageJ software (mean ± SD; n = 5; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (b) Recombinant human TSG-6 protein (50 ng/ml) was added in the culture medium just after scratching monolayers cultures of TSG-6–/– (a) and (b) cells (mean ± SD; n = 3 for TSG-6–/– (a) and n = 4 for TSG-6–/– (b); two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. 0 ng/ml rhTSG-6). Scratch closure was monitored by Cytonote and quantified by ImageJ software. h, hour; KC, keratinocyte; rhTSG-6, recombinant human TSG-6.

Lack of TSG-6 slows down KCs migration. Monolayers cultures of TSG-6+/+ or TSG-6–/– (a and b) N/TERT KCs at confluence stage treated with IL-4 and IL-13 for 24 h. KCs were growth arrested by mitomycin C before performing a scratch. KCs migration was monitored for 24 h in the presence of IL-4 and IL-13. (a) Scratch closure was analyzed every 6 h by Cytonote and was quantified by ImageJ software (mean ± SD; n = 5; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (b) Recombinant human TSG-6 protein (50 ng/ml) was added in the culture medium just after scratching monolayers cultures of TSG-6–/– (a) and (b) cells (mean ± SD; n = 3 for TSG-6–/– (a) and n = 4 for TSG-6–/– (b); two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. 0 ng/ml rhTSG-6). Scratch closure was monitored by Cytonote and quantified by ImageJ software. h, hour; KC, keratinocyte; rhTSG-6, recombinant human TSG-6.

Deficiency in TSG-6 protein is responsible for enhanced leakage of HA from the epidermis

The effect of TSG-6 deficiency on the synthesis and degradation of HA as well as on HA partners (CD44 and ITIH5) was investigated at the mRNA level. No alteration in knock-out tissues was observed (Figure 9). CD44 and HAS3 appeared overexpressed in both TSG-6+/+ and TSG-6–/– tissues stimulated with IL-4 and IL-13, in accordance with previous results (De Vuyst et al., 2018; do Nascimento Pedrosa et al., 2017; Malaisse et al., 2014).
Figure 9

mRNA levels of genes related to HA are unchanged in TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HAS1, (b) HAS2, (c) HAS3, (d) Hyal1, (e) Hyal2, (f) Hyal3, (g) Kiaa1199, (h) CD44, and (i) ITIH5 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05). CTL, control; HA, hyaluronan; RHE, reconstructed human epidermis.

mRNA levels of genes related to HA are unchanged in TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HAS1, (b) HAS2, (c) HAS3, (d) Hyal1, (e) Hyal2, (f) Hyal3, (g) Kiaa1199, (h) CD44, and (i) ITIH5 mRNA expression assessed by RT-qPCR. The reference gene is RPLP0 (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05). CTL, control; HA, hyaluronan; RHE, reconstructed human epidermis. Because HAS3 is enhanced in RHE exposed to IL-4 and IL-13 and produces proinflammatory low-molecular-weight HA fragments (Weigel et al., 1997), HA released in the culture medium was analyzed using size exclusion chromatography (Figure 10). In control conditions, both TSG-6+/+ and TSG-6–/– molecular-weight distribution of HA revealed similar profiles with a peak around 500 kDa (Figure 10a). When RHEs were stimulated with IL-4 and IL-13, the size of HA slightly shifted to lower molecular weights in both nonedited and knock-out tissues (Figure 10b).
Figure 10

Similar molecular-weight distribution of HA released in the culture medium of normal and challenged The culture medium of TSG-6+/+ and TSG-6–/– (a and b) RHE (a) control or (b) treated with IL-4 and IL-13 for 48 hours was collected and chromatographed on a Sephacryl S-1000 column. HA fragments of 2,500, 500, and 150 kDa were used as standards. Graphs compared HA size distribution in each fraction (mean ± SD; n = 2). #, number of; HA, hyaluronan; K, thousand; RHE, reconstructed human epidermis.

Similar molecular-weight distribution of HA released in the culture medium of normal and challenged The culture medium of TSG-6+/+ and TSG-6–/– (a and b) RHE (a) control or (b) treated with IL-4 and IL-13 for 48 hours was collected and chromatographed on a Sephacryl S-1000 column. HA fragments of 2,500, 500, and 150 kDa were used as standards. Graphs compared HA size distribution in each fraction (mean ± SD; n = 2). #, number of; HA, hyaluronan; K, thousand; RHE, reconstructed human epidermis. Morphologically, HA distribution inside TSG-6+/+ and TSG-6–/– tissues appeared similar. However, the intensity of the HA fluorescent signal was weaker in TSG-6–/– than in TSG-6+/+ RHE (Figure 11a). This concurred with a larger HA amount released into the culture medium from TSG-6–/– than that from TSG-6+/+ RHE, especially after IL-4 and IL-13 treatment (Figure 11b). Simultaneously, the amount of HA inside TSG-6–/– RHE was reduced, although it was nonsignificant (Figure 11c). The distribution of HA was similar in TSG-6+/+ and TSG-6–/– RHE infected with dermatophytes, and HA leakage was increased from infected TSG-6–/– RHE, concomitantly with a reduced content inside tissues. However, the differences were nonsignificant (Figure 12).
Figure 11

The change of HA localization in TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HE staining and HA detection (green). Bar = 20 μm. (b) HA released in culture medium measured by ELISA-like assay (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (c) HA remaining inside RHE (mean ± SD; n = 3; two-way ANOVA). CTL, control; HA, hyaluronan; HE, hemalun/eosin; RHE, reconstructed human epidermis; Th, T helper type.

Figure 12

The change of HA localization in TSG-6+/+ and TSG-6–/– (a and b) RHE infected or not infected with TR and recovered 4 days after infection. (a) PAS staining with α-amylase pretreatment and hemalun counterstaining and HA detection. Bar = 20 μm. (b) HA released in culture medium measured by ELISA-like assay (mean ± SD; n = 4; two-way ANOVA). (c) HA remaining inside RHE (mean ± SD; n = 3; two-way ANOVA). CTL, control; HA, hyaluronan; PAS, periodic-acid schiff; RHE, reconstructed human epidermis; TR, Trichophyton rubrum.

The change of HA localization in TSG-6+/+ and TSG-6–/– (a and b) RHEs were stimulated or not stimulated for 48 hours with IL-4 and IL-13. (a) HE staining and HA detection (green). Bar = 20 μm. (b) HA released in culture medium measured by ELISA-like assay (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+). (c) HA remaining inside RHE (mean ± SD; n = 3; two-way ANOVA). CTL, control; HA, hyaluronan; HE, hemalun/eosin; RHE, reconstructed human epidermis; Th, T helper type. The change of HA localization in TSG-6+/+ and TSG-6–/– (a and b) RHE infected or not infected with TR and recovered 4 days after infection. (a) PAS staining with α-amylase pretreatment and hemalun counterstaining and HA detection. Bar = 20 μm. (b) HA released in culture medium measured by ELISA-like assay (mean ± SD; n = 4; two-way ANOVA). (c) HA remaining inside RHE (mean ± SD; n = 3; two-way ANOVA). CTL, control; HA, hyaluronan; PAS, periodic-acid schiff; RHE, reconstructed human epidermis; TR, Trichophyton rubrum. To reintroduce TSG-6 in knocked-out tissues, knocked-out cells were mixed with nonedited N/TERT cells in varying proportions while reconstructing tissues (Figure 13). Interestingly, HA leakage from TSG-6–/– tissues was rescued depending on the proportion of TSG-6‒producing cells (Figure 13).
Figure 13

Introduction of Schematic representation of TSG-6–/– rescue test: epidermis were reconstructed from different ratios of TSG-6+/+ (green) and TSG-6–/– (red) KCs. HA released in the culture medium of epidermis reconstructed from mixtures of KCs TSG-6+/+ and TSG-6–/– (a) or TSG-6–/– (b) was measured by ELISA-like assay (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+ 100%). CTL, control; HA, hyaluronan; KC, kerationcyte; RHE, reconstructed human epidermis.

Introduction of Schematic representation of TSG-6–/– rescue test: epidermis were reconstructed from different ratios of TSG-6+/+ (green) and TSG-6–/– (red) KCs. HA released in the culture medium of epidermis reconstructed from mixtures of KCs TSG-6+/+ and TSG-6–/– (a) or TSG-6–/– (b) was measured by ELISA-like assay (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001 vs. TSG-6+/+ 100%). CTL, control; HA, hyaluronan; KC, kerationcyte; RHE, reconstructed human epidermis. For yet unidentified reason, clone TSG-6 (a) constantly releases a higher amount of HA than clone TSG-6 (b) as shown in Figures 11b and 14a despite similar expression levels of HASs and hyaluronidases (Figure 9).
Figure 14

Change of HA localization in Autocrine monolayers cultures of TSG-6+/+ or TSG-6–/– (a and b) N/TERT KCs at confluence stage treated or not treated with IL-4 and IL-13 for 48 hours. HA amount (a) released in the culture medium present (b) in pericellular spaces and (c) in intercellular compartment measured by ELISA-like assay and relative to the total amount of protein (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. TSG-6+/+). CTL, control; HA, hyaluronan; KC, keratinocyte.

Change of HA localization in Autocrine monolayers cultures of TSG-6+/+ or TSG-6–/– (a and b) N/TERT KCs at confluence stage treated or not treated with IL-4 and IL-13 for 48 hours. HA amount (a) released in the culture medium present (b) in pericellular spaces and (c) in intercellular compartment measured by ELISA-like assay and relative to the total amount of protein (mean ± SD; n = 3; two-way ANOVA; ∗P < 0.05; ∗∗P < 0.01 vs. TSG-6+/+). CTL, control; HA, hyaluronan; KC, keratinocyte. In addition, pericellular HA was quantified. The RHE model does not allow to distinguish pericellular from intracellular compartment. Therefore, HA distribution was studied using monolayer cultures stimulated or not with IL-4 and IL-13 for 48 hours at the confluence. Whereas HA released in the culture medium of TSG-6–/– monolayers (Figure 14a) was increased in response to treatment, pericellular (Figure 14b) or intracellular (Figure 14c) HA was significantly reduced.

Discussion

In normal conditions, TSG-6 is weakly expressed in mouse and human skin (Shakya et al., 2020; Tan et al., 2011). During in vitro reconstruction of the human epidermis, TSG-6 expression remains very weak as well. In infected RHE, TSG-6 expression is induced in response to recognition of fungal elements by toll-like receptors, leading to production of antimicrobial peptides and cytokines such as TNF-α (Faway et al., 2019; Lee et al., 1990). TSG-6 expression is also induced in RHE exposed to Th2 cytokines to mimic typical features of AD epidermis (do Nascimento Pedrosa et al., 2017), and HA production simultaneously increases as a result of HAS3 overexpression (De Vuyst et al., 2018; Malaisse et al., 2014). In vivo, TSG-6 expression is increased in lesional AD skin and in contact dermatitis (GEO accession GSE57225 [Quaranta et al., 2014] and GEO accession GSE711996 [Tan et al., 2017]) and is upregulated in Th1/Th17 diseases such as psoriasis (GEO accession GSE166388) and wound healing (GEO accession GSE50425). Altogether, these observations show the involvement of TSG-6 in many inflammatory conditions. To understand the role of epidermal TSG-6, deletions in the TNFAIP6 gene were carried out by CRISPR/Cas9 in N/TERT KCs. This immortalized human cell line was chosen because it allows for proper reconstruction of epidermis able to respond to Th2 cytokines (Dickson et al., 2000; Moran et al., 2021; Rheinwald et al., 2002; Smits et al., 2017). Recently, genome editing was performed in N/TERT cells using Cas9 plasmids for studies of Harlequin ichthyosis (Enjalbert et al., 2020) to target the genes linked to glycan metabolism (Dabelsteen et al., 2020) and to knockout claudin 1 and CD40 (Moran et al., 2021). In this study, ribonucleoprotein/Cas9 complexes were preferred to plasmid to reach higher efficiency of gene editing and benefit from a lower occurrence of off targets (Farboud et al., 2018). Indeed, all analyzed N/TERT cell clones harbored mutations. Among them, two clones exhibiting major deletions that conduct to null alleles were chosen for further studies. No off-target effect was found by screening RNA-seq results. TSG-6–/– KCs were used for RHE production. These tissues exhibit normal morphology in accordance with the unaltered epidermal phenotype of Tsg-6-knockout mice (Fülöp et al., 2003). In control and challenging conditions, deficiency of TSG-6 has no severe impact on gene expression levels as revealed by transcriptomic analysis. We noticed that the expression of decorin, a proteoglycan playing a role in cutaneous inflammatory reactions (Bocian et al., 2013), is downregulated. Of interest, KLK5 and KLK8 involved in skin homeostasis, KCs migration, and wound remodeling (Kishibe et al., 2012; Nauroy and Nyström, 2020) are also downregulated. In the scratch assay, whereas TSG-6–/– N/TERT clone (b) shows a very discreet slowdown, the second clone (a) migrates more slowly than control cells, reminiscent of the delayed wound closure in Tsg-6–/– mouse skin (Shakya et al., 2020). The cell migration rate is partially restored by reintroduction of soluble recombinant TSG-6 in the medium, as reported with epithelial corneal cells (Lee et al., 2015). Because the amount of HA is decreased in the pericellular matrix of TSG-6–/– KCs, we hypothesize that the defect in migration could result from poor pericellular HA organization. A disorganized ECM was described around preovulatory follicles in Tsg-6‒knockout mice (Fülöp et al., 2003). The HA interaction with its CD44 membrane receptor is facilitated by TSG-6 (Baranova et al., 2013), and this interaction promotes HA internalization together with proliferation and migration (Bourguignon, 2014; Csoka et al., 2001; Pasonen-Seppänen et al., 2012). TSG-6 deficiency may therefore influence the HA‒CD44 signaling axis and pathways located downstream of other receptors such as EGFR (Perez et al., 2013). In TSG-6–/– RHE, a change in the localization of HA is observed. A large amount of HA diffuses into the subepidermal compartment, suggesting a role for TSG-6 in maintaining HA inside the tissue. HA leakage from the epidermis was already reported in skin equivalent models and was attributed to an improper organization of the basement membrane (Iriyama et al., 2021; Tammi et al., 2000). Nonetheless, Rosdy et al. (1993) had illustrated typical basal lamina elements in RHE. Because TSG-6 is not a constituent of the basement membrane, enhanced release of HA from RHE made of TSG-6–/– cells cannot be attributed to a lack of integrity or to an altered composition of the basement membrane. We propose that the lack of HA retention in ECM is a more likely explanation. Future studies are required to complement the data obtained with our single-cell type in vitro model. Understanding the consequences of massive HA release due to the absence of TSG-6 on underlying dermal fibroblasts and immune cells is of interest. For instance, recruitment and migration of immune cells to an inflammatory site depend on the CD44-mediated organization of HA into cables (Jokela et al., 2008; Kobayashi et al., 2020). Our serum-free culture conditions do not allow to assess the enzymatic activity of TSG-6. Adding one source of ITI is required to investigate HA covalent cross-linking through a heavy chain. More precisely, the involvement of heavy chain 5 in ECM stabilization was shown in AD and psoriasis and should therefore be further analyzed (Huth et al., 2020). In conclusion, our data suggest a role for TSG-6 in retaining HA in inflamed epidermis by interaction and cross-linking in the ECM, as shown in other tissues (Day and Milner, 2019). In AD lesions, Th2-stimulated production of TSG-6 could participate in HA accumulation between KCs (Malaisse et al., 2014) (Figure 15). More broadly, when interacting with ECM components, TSG-6 potentially modulates the presentation of GFs and cytokines to neighboring receptors. In the epidermis, it could influence the signaling pathways related to HA (Kavasi et al., 2017) and play roles in inflammatory skin conditions.
Figure 15

TSG-6 retains HA in the ECM of the epidermis. (a) Schematic representations of in vitro observations in TSG-6+/+ and TSG-6–/– RHE challenged or not challenged. In control conditions, weak amounts of HA (green) and TSG-6 (red) are released in the culture medium. In inflammatory conditions, increased TSG-6 secretion goes along with enhanced HA release. In the absence of TSG-6, HA release is significantly increased, especially after the challenge. (b) Hypothetical role of TSG-6 in HA distribution and organization in healthy and lesional atopic dermatitis epidermis. In the healthy epidermis, the weak amount of secreted TSG-6 is involved in ECM organization by HA cross-linking. In atopic dermatitis lesions, enlarged intercellular spaces go along with enhanced HA production. Secreted TSG-6 links HA cables to form a dense matrix. ECM, extracellular matrix; HA, hyaluronan; RHE, reconstructed human epidermis.

TSG-6 retains HA in the ECM of the epidermis. (a) Schematic representations of in vitro observations in TSG-6+/+ and TSG-6–/– RHE challenged or not challenged. In control conditions, weak amounts of HA (green) and TSG-6 (red) are released in the culture medium. In inflammatory conditions, increased TSG-6 secretion goes along with enhanced HA release. In the absence of TSG-6, HA release is significantly increased, especially after the challenge. (b) Hypothetical role of TSG-6 in HA distribution and organization in healthy and lesional atopic dermatitis epidermis. In the healthy epidermis, the weak amount of secreted TSG-6 is involved in ECM organization by HA cross-linking. In atopic dermatitis lesions, enlarged intercellular spaces go along with enhanced HA production. Secreted TSG-6 links HA cables to form a dense matrix. ECM, extracellular matrix; HA, hyaluronan; RHE, reconstructed human epidermis.

Materials and Methods

KCs origins, cultures, and reconstruction of the epidermis

Skin biopsies (University Hospital in Liège, Belgium) from patients with AD were collected on lesional zones as described (Mathay et al., 2011). Healthy skin biopsies and normal human KCs were isolated from abdominoplasties (Clinique St Luc, Namur, Belgium) as described (Minner et al., 2010). All samples were obtained after written informed consent, in accordance with Medical Ethical Committees of the University Hospital in Liège, Clinique St-Luc in Namur, and Saint-Louis Hospital in Paris (France) and according to the Declaration of Helsinki Principles. Immortalized human KCs (N/TERT-1) were obtained from J.G. Rheinwald’s laboratory (Dickson et al., 2000). Cells authenticity was assessed at ATCC (Manassas, VA) (http://www.atcc.org) by a comparison of an original frozen vial of cells (Bertrand-Vallery et al., 2010) with current cultured N/TERT cells. Normal and immortalized KCs were cultured as autocrine monolayer cultures (Minner et al., 2010) or were grown on a polycarbonate filter and exposed to the air‒liquid interface for 11 days to produce a reconstructed epidermis (De Vuyst et al., 2014; Frankart et al., 2012). Whereas autocrine monolayers were cultured in EpiLife culture medium supplemented by amino acids (Sigma-Aldrich, St. Louis, MO) and hydrocortisone (0.18 μg/ml; Sigma-Aldrich), RHEs were grown in EpiLife medium containing Human Keratinocytes Growth Supplement (Cascade Biologics, Mansfield, United Kingdom) and calcium ion (1.5 mM), KC GF (10 ng/ml) (R&D Systems, Abingdon, United Kingdom), and vitamin C (50 mg/ml). Both cultures were maintained at 37 °C and 5% carbon dioxide.

Challenges applied to RHE to mimic AD or dermatophytosis

A described procedure was used to mimic AD alterations in RHE (De Vuyst et al., 2018). Briefly, on the 11th day of reconstruction and after cholesterol depletion with methyl-β-cyclodextrin (7.5 mM; Sigma-Aldrich, Diegem, Belgium) for 2 hours, RHEs were incubated for 48 hours with IL-4 (50 ng/ml), IL-13 (50 ng/ml), and IL-25 (20 ng/ml) (PeproTech, London, United Kingdom). RHEs were eventually challenged with IL-4 and IL-13 only for 48 hours to produce AD-like alterations in a simpler way (do Nascimento Pedrosa et al., 2017). Another challenge was applied using infection with T. rubrum arthroconidia (Faway et al., 2017). Briefly, topical application of dermatophyte arthroconidia in suspension in PBS over RHE was performed on the 11th day of reconstruction. After 4 hours of exposure to these dermatophytes, RHEs were washed with PBS and then re-exposed to air for 4 additional days of culturing. This protocol was adapted for infection of RHE using N/TERT KCs. Owing to a weaker barrier, an inoculum of 850 arthroconidia per cm2 was used instead of 1,700 per cm2.

Scratch assay

KCs migration was studied using a scratch assay in culture (Mound et al., 2017). Confluent monolayers were treated 1 hour with 10 μg/ml mitomycin C (#sc 3514B; Santa Cruz Biotechnology, Dallas, TX) before scraping using a p200 micropipette tip and then washed and incubated with fresh medium for 24 hours. Rescue was performed by adding 50 ng/ml recombinant human TSG-6 (#2104-TS; R&D Systems) in the culture medium during the 24 hours of scratch closure. The scratch closure was monitored by time-lapse imaging using an Iprasense Cytonote device (Iprasense, Montpellier, France) and analyzed by ImageJ software (National Institutes of Health, Bethesda, MD).

Histological staining, immunofluorescence, and HA labeling

For morphological analysis, RHE and skin biopsies were fixed in a 4% formaldehyde, 70% ethanol, and 5% acetic acid solution (Lin et al., 1997) and then dehydrated and embedded in paraffin to prepare 7 μm tissue sections using a microtome. Tissue sections were stained with hemalun-eosin(-saffron). Infected RHEs were stained by periodic-acid schiff after incubation of sections with α-amylase (Faway et al., 2017). For immunolabeling, tissue sections were incubated 1 hour at room temperature with primary anti-keratin 10 (1:100 dilution, Dako, Glostrup, Denmark), anti-involucrin (1:200 dilution, Sigma-Aldrich), and anti-loricrin (1:100 dilution, Abcam. Cambridge, United Kingdom) or overnight at 4 °C with primary anti-FLG (1:75 dilution, Thermo Fisher Scientific, Waltham, MA) antibodies and then incubated 1 hour at room temperature with secondary antibodies conjugated with Alexa fluorochrome A488: anti-mouse diluted 1:200 for keratin 10, involucrin, and FLG and anti-rabbit diluted at 1:500 for loricrin (Life Technologies, Carlsbad, CA). HA was localized on tissue sections using an overnight incubation at 4 °C with a biotinylated HA-binding protein at 1:100 dilution (Calbiochem-Merck, Darmstadt, Germany), followed by 1-hour incubation at room temperature with fluorescein-labeled streptavidin (Vector Laboratories, Burlingame, CA). Sections were mounted with coverslips using Mowiol (Sigma-Aldrich) and analyzed under an Olympus BX63 microscope.

RNA extraction and RT-qPCR

Total RNA was extracted using a High Pure RNA isolation kit (Roche, Basel, Switzerland) for monolayers and an RNeasy kit (Qiagen, Hilden, Germany) for RHE. First-strand cDNA was synthesized using SuperScript III (Invitrogen, Carlsbad, CA) and amplified by FastStart Universal SYBR Green Master (Roche) in the Light Cycler (Light Cycler 96 Instrument; Roche). Relative mRNA expression level was normalized to RPLP0 reference gene (Minner and Poumay, 2009) and calculated using the ΔΔCq method (Livak and Schmittgen, 2001). Specifics primers (Eurogentec, Liège, Belgium) were used at 300 nM and are listed in Table 3.
Table 3

The Sequence of Primer Pairs Used for RT-qPCR

Gene SymbolForward PrimerReverse Primer
RPLP0ATCAACGGGTACAAACGAGTCCAGATGGATCAGCCAAGAAGG
TSG-6GCTAGAGGCAGCCAGAAAAAATCCAACTCTGCCCTTAGCC
HAS1CACCTCACCAACCGCATGCTGTCTGCTGGCTCAGCCACCG
HAS2CCTCATCATCCAAAGCCTGTGATGCAAAGGGCAACTGTTT
HAS3GCGATTCGGTGGACTACATGGATCCTCCTCCAGGACTC
Hyal1TCCAGATCTTCTATGACACGACAAATCCTTGATGGCCTGACATGA
Hyal2GCACTCCCAGTCTACGTCTTCAGCACTCTCGCCAATGGTAGAG
Hyal3GATCTGGGAGGTTCCTGTCCAGAGCTGGAGAGGCTCAGGT
Kiaa1199TACTGGGACGAGGACTCAGGTCACAGCCTTTCATGGAGCA
CD44CATTGCAGTCAACAGTCGAAGAAATTGCCACTGTTGATCACTAGCTT
ITIH5ACTGTCGCTGGAGAACTGTGCGGGGTCCTGATTTCATCGT
LL37CCAGGACGACACAGCAGTCACTTCACCAGCCCGTCCTTC
The Sequence of Primer Pairs Used for RT-qPCR

RNA-seq

Total RNA was extracted using the ReliaPrep RNA Tissue Miniprep System (Promega, Madison, WI) before checking integrity through an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA). RNA-seq and standard analysis were performed by the GIGA-Genomics platform (University of Liège, Belgium). RNA-seq data are available in the ArrayExpress database at EMBL’s European Bioinformatics Institute (http://www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-10095.

ELISA assay

TSG-6 protein levels were quantified in a culture medium using a homemade ELISA assay. Plates were coated overnight with an anti‒TSG-6 capture antibody (#A38.1.20; Santa Cruz Biotechnology) and were blocked for 1 hour with 0.5% BSA in PBS. After 2 hours of incubation in the presence of samples of culture medium or in the presence of recombinant human TSG-6 protein (#2104-TS; R&D Systems) samples at different concentrations chosen to calibrate the assay, TSG-6 was secondarily detected using an anti‒TSG-6 biotinylated antibody (#BAF2104; R&D Systems), horseradish peroxidase-streptavidin, and chromogenic substrate (R&D Systems). The concentration of TSG-6 protein is proportional to absorbance (λ: 450–570 nm) measured by spectrophotometer (Molecular Device, Berkshire, United Kingdom).

Procedures for HA analysis

The amount of HA present in the culture medium, in pericellular, and in intracellular compartments was analyzed using Hyaluronan DuoSet kit (R&D Systems) as described by Malaisse et al. (2014). To measure the size of HA fragments, the total amount of HA recovered in the RHE culture medium was assessed using the Hyaluronan DuoSet kit. Samples were then prepared to load an identical amount of HA (10 ng of HA per sample) on a Sephacryl S-1000 column (Amersham Biosciences, Uppsala, Sweden). Samples were eluted at 70 μl/min with 0.15 M sodium acetate and 0.1% CHAPS (Sigma-Aldrich, Steinheim, Germany) to collect 30 fractions from each sample. Standard HA fragments were chromatographed in identical conditions.

In situ hybridization

RNAscope kit was used (Advanced Cell Diagnostics, Silicon Valley, CA) in accordance with the manufacturer’s instructions (RNAscope 2.5 HD assay Red). Tissue sections were targeted by specific probes Homo sapiens-TNFAIP6 (reference 444421). Homo sapiens-PPIB (Peptidyl-prolyl cis-trans isomerase B) probe (reference 313901) was used as the positive control, and DapB (Bacillus subtilis dihydrodipicolinate reductase) probe (reference 310043) as the negative control. Fast Red staining in sections counterstained by hemalun (blue) reveals sequence hybridization as described by Wang et al. (2012).

Genome editing by CRISPR/Cas9

N/TERT KCs genome was edited to generate TNFAIP6-knockout cells. The procedure is a modified version of what is described in the user guide provided with the Alt-R CRISPR-Cas9 system (Integrated DNA Technologies, Coralville, IA). Practically, CRISPR RNA (crRNA) trans-activating crRNA duplexes were formed by mixing crRNA (200 μM) and trans-activating crRNA (200 μM, #1072532; Integrated DNA Technologies), hybridizing them by heating up to 95 °C for 5 minutes, and then cooling them down to room temperature. Five crRNAs were designed using http://crispor.tefor.net/, https://zlab.bio/guide-design-resources and eu.idtdna.com website and are listed as follows: #Exon1 crRNA Up 5′-TTGCTATGGGAAGACACTCA AGG-3′, #Exon1 crRNA Down 5′-AAGGTTCTTACCAAGCCATA TGG-3′, #Exon3 crRNA Up 5′-ATTATGGAATCCGTCTCAAT AGG-3′, #Exon6 crRNA Up 5′-GAGAATCGACCTCTGAATTC TGG-3′, and #Exon6 crRNA Down 5′-GGTCATGACATTTCCTGTTG GGG-3' (Figure 5a). Duplexes were then mixed to Cas9 enzyme (#1081058; Integrated DNA Technologies) to form ribonucleoprotein complexes, which were electroporated (Amaxa Nucleofector I device, #AAD-1001; Lonza, Bâle, Switzerland) in proliferative N/TERT KCs according to the manufacturer’s instructions (Amaxa Human Keratinocyte Nucleofector Kit of Lonza, #VAPD-1002). After 48 hours of recovery, edited cells were seeded in 96-well plates at low density to obtain clonal populations (limiting dilution). The mutations in the TNFAIP6 gene were amplified by PCR. PCR products corresponding to homologous alleles of TSG-6 were isolated by TA cloning and sequenced by Sanger sequencing (Genewiz company, Takeley, United Kingdom). PCR primers (Eurogentec) used at 10 μM were exon 1 forward 5′-ACAGCGGAAAATGCCATTGG-3′, exon 1 reverse 5′-TCTCAACAAGAGGCTCATGAGT-3′, exon 3 forward 5′-TTGGAAGTGGCTATCTTTGCT-3′, exon 3 reverse 5′-GGAACACAACTCAATGTCGGG-3′, exon 6 forward 5′-CATGCGGCAATGAAACCACA-3′, and exon 6 reverse 5′-AAGTGGCTAAATCTTCCAGCT-3' (Figure 5b).

Statistical analyses

All data were analyzed using SigmaPlot software, version 11.0 (Systat Software, Chicago, IL). Statistical differences were evaluated by Student t-test, one-way ANOVA, and two-way ANOVA repeated measure tests, and error bars represented SD. All graphs were performed using GraphPad Prism 5 software (GraphPad Software, San Diego, CA).

Data availability statement

Dataset related to this article can be found at https://data.mendeley.com/datasets/4fbgg2syzk/1, an open-source online data repository hosted at Mendeley Data.

ORCIDs

Céline Evrard: http://orcid.org/0000-0002-7450-0759 Emilie Faway: http://orcid.org/0000-0003-3406-1103 Evelyne De Vuyst: http://orcid.org/0000-0001-7152-1370 Olivier Svensek: http://orcid.org/0000-0001-9254-0060 Valérie De Glas: http://orcid.org/0000-0001-7366-2601 David Bergerat: http://orcid.org/0000-0003-0013-6621 Michel Salmon: http://orcid.org/0000-0002-1728-8448 Olivier De Backer: http://orcid.org/0000-0002-1490-4508 Bruno Flamion: http://orcid.org/0000-0003-4458-7627 Hélène Le-Buanec: http://orcid.org/0000-0002-2139-0357 Catherine Lambert de Rouvroit: http://orcid.org/0000-0002-0273-8995 Yves Poumay: http://orcid.org/0000-0001-5200-3367

Author Contributions

Conceptualization: CE, CLDR, BF, YP; Data Curation: CE, CLDR, YP; Formal Analysis: CE; Funding Acquisition: CLDR, YP; Investigation: CE, VDG, EF, MS, DB; Methodology: CE, EF, EDV, OS, ODB, DB, HLB, CLDR, YP; Project Administration: CE, CLDR, BF, YP; Supervision: CE, CLDR, YP; Validation: CE, CLDR, YP; Visualization: CE, CLDR, YP; Writing - Original Draft Preparation: CE, CLDR, YP; Writing - Review and Editing: CE, CLDR, YP
  68 in total

Review 1.  Hyaluronan: from extracellular glue to pericellular cue.

Authors:  Bryan P Toole
Journal:  Nat Rev Cancer       Date:  2004-07       Impact factor: 60.716

2.  Study of epidermal differentiation in human keratinocytes cultured in autocrine conditions.

Authors:  Frédéric Minner; Françoise Herphelin; Yves Poumay
Journal:  Methods Mol Biol       Date:  2010

3.  Cutaneous Wounds in Mice Lacking TSG-6 Exhibit Delayed Closure and an Abnormal Inflammatory Response.

Authors:  Sajina Shakya; Judith A Mack; Minou Alipour; Edward V Maytin
Journal:  J Invest Dermatol       Date:  2020-05-15       Impact factor: 8.551

Review 4.  The six hyaluronidase-like genes in the human and mouse genomes.

Authors:  A B Csoka; G I Frost; R Stern
Journal:  Matrix Biol       Date:  2001-12       Impact factor: 11.583

5.  Reconstruction of normal and pathological human epidermis on polycarbonate filter.

Authors:  Evelyne De Vuyst; Céline Charlier; Séverine Giltaire; Valérie De Glas; Catherine Lambert de Rouvroit; Yves Poumay
Journal:  Methods Mol Biol       Date:  2014

6.  Intraindividual genome expression analysis reveals a specific molecular signature of psoriasis and eczema.

Authors:  Maria Quaranta; Bettina Knapp; Natalie Garzorz; Martina Mattii; Venu Pullabhatla; Davide Pennino; Christian Andres; Claudia Traidl-Hoffmann; Andrea Cavani; Fabian J Theis; Johannes Ring; Carsten B Schmidt-Weber; Stefanie Eyerich; Kilian Eyerich
Journal:  Sci Transl Med       Date:  2014-07-09       Impact factor: 17.956

7.  Inter-α-Trypsin Inhibitor Heavy Chain 5 (ITIH5) Is a Natural Stabilizer of Hyaluronan That Modulates Biological Processes in the Skin.

Authors:  Sebastian Huth; Laura Huth; Yvonne Marquardt; Katharina Fietkau; Edgar Dahl; Philipp R Esser; Stefan F Martin; Ruth Heise; Hans F Merk; Jens Malte Baron
Journal:  Skin Pharmacol Physiol       Date:  2020-08-14       Impact factor: 3.479

8.  Generalized Correlation Coefficient for Non-Parametric Analysis of Microarray Time-Course Data.

Authors:  Qihua Tan; Mads Thomassen; Mark Burton; Kristian Fredløv Mose; Klaus Ejner Andersen; Jacob Hjelmborg; Torben Kruse
Journal:  J Integr Bioinform       Date:  2017-06-06

9.  Developmental cell programs are co-opted in inflammatory skin disease.

Authors:  Gary Reynolds; Peter Vegh; James Fletcher; Elizabeth F M Poyner; Emily Stephenson; Issac Goh; Rachel A Botting; Ni Huang; Bayanne Olabi; Anna Dubois; David Dixon; Kile Green; Daniel Maunder; Justin Engelbert; Mirjana Efremova; Krzysztof Polański; Laura Jardine; Claire Jones; Thomas Ness; Dave Horsfall; Jim McGrath; Christopher Carey; Dorin-Mirel Popescu; Simone Webb; Xiao-Nong Wang; Ben Sayer; Jong-Eun Park; Victor A Negri; Daria Belokhvostova; Magnus D Lynch; David McDonald; Andrew Filby; Tzachi Hagai; Kerstin B Meyer; Akhtar Husain; Jonathan Coxhead; Roser Vento-Tormo; Sam Behjati; Steven Lisgo; Alexandra-Chloé Villani; Jaume Bacardit; Philip H Jones; Edel A O'Toole; Graham S Ogg; Neil Rajan; Nick J Reynolds; Sarah A Teichmann; Fiona M Watt; Muzlifah Haniffa
Journal:  Science       Date:  2021-01-22       Impact factor: 47.728

Review 10.  Hyaluronan: Metabolism and Function.

Authors:  Takashi Kobayashi; Theerawut Chanmee; Naoki Itano
Journal:  Biomolecules       Date:  2020-11-07
View more
  1 in total

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

Authors:  Jos P H Smits; Luca D Meesters; Berber G W Maste; Huiqing Zhou; Patrick L J M Zeeuwen; Ellen H van den Bogaard
Journal:  JID Innov       Date:  2021-12-01
  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.