EMBO/SourceData huggingface dataset:
Figure 3A. Human GSDMD at a concentration of 2 μM was incubated at room temperature with 5 nM caspase-1. The protein was cleaved in a time-dependent manner into two bands of 31 kDa (GSDMDNterm) and 22 kDa (GSDMDCterm).B. Cross-linking experiment of full-length and cleaved GSDMD. GSDMD at a concentration of 2 μM was incubated at room temperature with 5 nM caspase-1. After enzymatic cleavage, GSDMDNterm is highly cross-linked by DSS, resulting in the gel-impenetrating species highlighted by the arrow. GSDMDCterm is not cross-linked.C. GSDMD at a concentration of 1 μM was incubated at room temperature with 5 nM caspase-1 and liposomes composed of 4 mM DMPC or polar lipid extract derived from E. coli. After 2 h, the lipid fraction (L) was separated from the supernatant (S) by ultracentrifugation at 4C for 1 h at 120,000g.
Figure 4A-I. Dye release time courses from liposomes as a percentage of maximal release. A. Five different reactions, where 5 nM caspase-1 and 400 μM 6-carboxifluorescein-loaded liposomes prepared with E. Coli polar lipid extract were incubated with GSDMD concentrations of (nM): 520, 260, 130, 100, 65 (colored dark to light orange). The time point of 20 min is highlighted by a vertical dashed line.A-I. Dye release time courses from liposomes as a percentage of maximal release. B. Five different reactions, where 130 nM of GSDMD and 400μM 6-carboxifluorescein-loaded liposomes prepared with E. Coli extract polar lipid were incubated with caspase-1 concentration of (nM): 15, 8, 5, 2.5, 1.2 (colored dark to light blue). The time point of 20 min is highlighted by a vertical dashed line.A-I. Dye release time courses from liposomes as a percentage of maximal release. C. Three different reactions, where 5 nM caspase-1 and 400 μM 6-carboxifluorescein-loaded liposomes prepared with porcine brain total lipid extract were incubated with GSDMD concentrations of (nM): 520, 260,100 (colored dark to light green).A-I. Dye release time courses from liposomes as a percentage of maximal release. D. Dye release at 20 minutes reaction as a function of GSDMD (dark orange) and caspase-1 (dark blue) concentrations. Errors bars for three independent experiments are shown.A-I. Dye release time courses from liposomes as a percentage of maximal release. E. Two different sets of reactions, where wild type GSDMD (dark to light orange) and the mutant GSDMDI104N (dark to light blue) were independently incubated at the concentrations of 260, 130 and 65 nM with 5 nM caspase-1 and 400 μM 6-carboxifluorescein-loaded liposomes. The time point of 60 minutes is highlighted by a vertical dashed line.A-I. Dye release time courses from liposomes as a percentage of maximal release. F. Dye release at 60 minutes of reaction as a function of GSDMD wild type (dark orange) and GSDMDI104N (dark blue) concentration. Error bars for three independent experiments are shown.A-I. Dye release time courses from liposomes as a percentage of maximal release. G-I. Dye release from 400 μM liposomes loaded with the 6-carboxifluorescein-derivates FD-20, FD-40 and FD-150, with variable Stokes diameters, as indicated. 130 nM of GSDMD and 5 nM caspase-1 were incubated with the liposomes.
Figure 5A-C. Cryo-EM micrographs of GSDMDNterm pores in E. coli extract polar lipids liposomes. The micrographs were acquired at protein/lipid molar ratios of 1/1000, 1/500 and 1/100, respectively. Black arrows indicate ring-shaped structures corresponding to oligomeric GSDMDNterm pore forms. Scale bars = 80 nm.D. Proteoliposome with protein/lipid molar ratio of 1/100 at higher magnification. Scale bar = 80 nm.
Figure 6A. AFM topograph of GSDMDNterm bound to lipid membranes of E. coli extract polar lipids. Overview topograph showing arc-, slit- and ring-like GSDMD oligomers.B-E. High-resolution AFM topographs and height profiles of GSDMDNterm oligomers (red curves). Height profiles were measured along the red lines in the topographs. Dashed zero lines (0 nm) indicate the surface of the lipid membrane; heights of <−4 nm indicate that the oligomers formed transmembrane pores. Topographs were taken in buffer solution (150 NaCl, 20 mM Hepes, pH 7.8) at room temperature. The experiment was reproduced at least three times using independent liposome and GSDMD preparations.F-H. Analysis of GSDMDNterm oligomers bound to lipid membranes of E. coli extract polar lipids. F. GSDMDNterm oligomers assembled into arcs, slits and rings leading to transmembrane (TM) pores or not transmembrane aggregate (n = 218). G. Height of GSDMDNterm oligomers protruding from the lipid membrane. Bars represent mean and error bars represent s.e.m. (n = 218). H. Distribution of the diameters of rings formed by GSDMDNterm oligomers. The average distribution was 21.2 ± 5.6 nm (n = 164; average ± s.d.) and the bin size was 3 nm.
Figure 1(E) Expression levels of LSD1 in MGC-803, MKN-45, HGC-27, BGC-823, and NCI-N87 cell lines and their corresponding sEVs. The samples with equal amounts of proteins were loaded (n = 3 biological replicates; mean ± SEM; compared with the NCI-N87; *P = 0.0432 (cell/MGC-803), *P = 0.0300 (cell/HGC-27), **P = 0.0306 (cell/BGC-823), **P = 0.0018 (sEVs/MGC-803), **P = 0.0029 (sEVs/HGC-27), and ***P = 0.0003 (sEVs/BGC-823); two-tailed unpaired Student's t-test; GAPDH was used as a loading control for cell lysis; CD9 was used as a loading control for sEV lysis).(J) Expression levels of LSD1, CD63, and CD9 in sEVs with indicated treatment.
Figure 2(A) Confocal microscopy image analysis of sEV fusion to MGC-803 cells. The MGC-803 (left side) and MKN-45 (right side) cells were treated with sEVs derived from MGC-803 cells (Con sEVs) or LSD1 knockout (KO) MGC-803 cells (KO sEVs) and stained with PKH26 for 12 h. Additionally, the cell membrane was stained with Dio, while the nuclei were stained with 4',6-diamidino-2-phenylindole (DAPI). Scale bar = 100 µm.(B) Immunofluorescence confocal microscopy analysis of LSD1 (green) in LSD1 KO MGC-803 cells (left panel) and LSD1 KO MKN-45 cells (right panel) incubated with 20 μg/mL Con sEVs and KO sEVs for 12 h. The cell membrane was stained with PKH26, while the nuclei were stained with DAPI. Scale bar = 100 µm.
Figure 3(F) Expression levels of LSD1, Nanog, OCT4, SOX2, and CD44 in LSD1 KO MGC-803 cells.(N) Immunoprecipitation of Kme1/2 on SOX2 (WT indicates HEK293T cells co-transfected with LSD1-WT and SOX2-WT; Mut indicates HEK293T cells co-transfected with LSD1-WT and SOX2-Mut; WT, wild type; Mut, mutant).
Figure 4(C) Expression levels of LSD1, OCT4, and SOX2 in MGC-803 and MKN-45 cells incubated with 20 μg/mL sEVs from MGC-803 or LSD1 KO MGC-803 cells for 48 h. (D) Quantification of the results of (C) (n = 3 biological replicates; mean ± SEM; *P = 0.0385 (LSD1), *P = 0.0160 (OCT4), and ***P = 0.0006 (SOX2) for MGC-803; *P = 0.0299 (LSD1), *P = 0.0147 (OCT4), and *P = 0.0258 (SOX2) for MKN-45; two-tailed unpaired Student's t-test). (F) Sphere formation assay results of MGC-803 (upper panel) and MKN-45 (bottom panel) cells incubated with 20 μg/mL sEVs from HEK293T cells (left panel), WT-LSD1 sEVs (middle panel), and LSD1 K661A sEVs (right panel) for 7 days. Scale bar = 100 µm (n = 3 biological replicates; mean ± SEM; **P = 0.0021 (MGC-803) and **P = 0.0020 (MKN-45); two-tailed unpaired Student's t-test).
Figure 5(F) Immunofluorescence image and quantification of the expression of SOX2 in tumor tissues subjected to in vivo limiting dilution assay. Scale bar = 50 µm (n = 3 biological replicates; mean ± SEM; **P = 0.0024, **P = 0.0088, and *P = 0.0342; two-tailed unpaired Student's t-test).(I and J) Expression levels of LSD1 in sEVs isolated from the plasma. CD9 and CD63 were used as markers of sEVs. CD9 was used as a loading control for sEV lysis (n = 3 biological replicates; mean ± SEM; *P = 0.0256 (patient 1), **P = 0.0068 (patient 2), **P = 0.0013 (patient 3), **P = 0.0031 (patient 4), ***P = 0.0009 (patient 5), **P = 0.0012 (patient 6), ***P = 0.0007 (patient 7), ***P = 0.0010 (patient 8), ***P = 0.0001 (patient 9), and *P = 0.0461 (patient 10); two-tailed unpaired Student's t-test).(K) Sphere formation assay results of MGC-803 cells treated with sEVs as indicated for 7 days. Scale bar = 100 µm (n = 3 biological replicates; mean ± SEM; **P = 0.0019 and **P = 0.0050; two-tailed unpaired Student's t-test).
Figure 6Proliferation assay results of MGC-803 (A), MKN-45 (B), NCI-N87 (C) cells treated with oxaliplatin along with indicated treatments (n = 3 biological replicates, mean ± standard error of mean (SEM); *P = 0.0204 (MGC-803), *P = 0.0498 (MKN-45), *P = 0.0473 (NCI-N87), and *P = 0.0019 (LSD1 KO MGC-803); two-tailed unpaired Student's t-test).Representative images of tumors (E) and the tumor weight (F) of mice treated with oxaliplatin in the presence or absence of sEVs (n = 3 biological replicates; mean ± SEM; *P = 0.0488; two-tailed unpaired Student's t-test).
Figure 1Loss of LRRC8D causes resistance to carboplatin and cisplatin, but not to oxaliplatin. Survival of parental, vector‐transduced, or LRRC8D‐deficient GT1 and GT2 KBM7 cells exposed for 96 h to increasing concentrations of cisplatin, carboplatin, and oxaliplatin. The corresponding IC50 values and 95% confidence interval (CI) are given in Appendix Table S2. Data are presented as mean ± SEM.
Figure 3A, BVRAC currents (ICl,vol) of the HAP1 (A) and KBM7 (B) haploid cell lines. Left panels, example current traces of ICl,vol fully activated by hypotonic cell swelling measured with the voltage‐clamp protocol shown in (A). Dashed lines indicate zero current. Right panels, averaged current/voltage relationships of maximally activated ICl,vol. Consistent with VRAC currents, they needed hypotonic swelling for activation, displayed an I− > Cl− permeability sequence, and were blocked by DCPIB (Appendix Fig S2A-H). The difference in current inactivation between HAP1 and KBM7 cells can be explained by the fact that KBM7 cells hardly express LRRC8E (Fig EV1) which accelerates VRAC inactivation (Voss et al, 2014). At potentials > +100 mV, also KBM7 currents inactivated (Appendix Fig S2I). Data are presented as mean ± SEM; n = 5-10.CCDependence of regulatory volume decrease (RVD) of HEK cells on LRRC8 genes. Cells were exposed to hypotonic medium starting at t = 0, and intracellular calcein fluorescence was followed over ˜1 h as semiquantitative measure of cell volume. Data are presented as mean values ± SEM from sixteen wells.
Figure 4LRRC8 subunit‐ and osmolarity‐dependent caspase induction in HCT116 cellsA, BCisplatin‐induced caspase activity in the continuous presence of 200 μM cisplatin under isotonic conditions (A), or after 1 h exposure to 200 μM cisplatin under iso‐ and hypotonic conditions (B), was followed over time in WT, HAP1−/−, LRRC8D−/−, and LRRC8 −/−HCT116 cells. Results from HAP1−/− and LRRC8D−/− were obtained with two different clonal cell lines each and averaged.CCaspase activation after 1‐h exposure to 4 μM staurosporine under iso‐ or hypotonic conditions of WT, HAP1−/−, LRRC8D−/−, and LRRC8 −/−HCT116 cells.Data information: Data are presented as mean ± SEM, n = 3-6. *P < 0.05; **P < 0.01; and ***P < 0.001. Similar results were obtained in three independent experiments. Fold change in (A) refers to t = 0. Control experiments indicated that hypotonicity per se had no effect (Appendix Fig S4).
Figure 5Activation of LRRC8 channels by pro‐apoptotic stimuliA-CCisplatin‐induced iodide influx into WT (black ■), but not HAP1−/− (green ▼) HEK cells indicates VRAC halide current activation during apoptosis. Cells expressing an iodide‐sensitive YFP variant were exposed to 200 μM cisplatin for periods of 0.5 h (A), 4.5 h (B), or 8.5 h (C) before adding extracellular I− (50 mM final). The difference in slopes of YFPfluorescence quenching between control and cisplatin‐treated cells semiquantitatively reflects VRAC current activation. Note that increased YFP quenching with cisplatin preincubation is not due to large non‐specific leaks resulting from cell morbidity. Such leaks should lead to a fast component of YFP quenching in WT, but not HAP1−/− cells after the pipetting artifact that immediately follows addition of iodide (indicated by arrows).DSwelling‐induced iodide influx into WT (black ■) and HAP1−/− (green ▼) HEK cells for comparison. Iodide (50 mM final) was added in isotonic or hypotonic (230 mOsm final) solution at the time indicated by arrow.E, FTime course of VRAC activation by 200 μM cisplatin (E) or 4 μM staurosporine (F) determined as in (A-C). Averaged maximal slopes of YFP quenching from eight wells (E) or 16 wells (F) each were evaluated to estimate iodide influx rates. WT (black ■) and HAP1−/− (green ▼). Data are presented as mean ± SEM.
Figure 6Carboplatin uptake into control KBM7 cells and two LRRC8D‐deficient clones (n = 6).Hypotonicity‐stimulated cisplatin uptake in HEK WT cells was selectively blocked by 100 μM carbenoxolone (CBX), a non‐specific blocker of VRAC (n = 3).Cisplatin uptake into HEK cells of indicated genotypes using 40 μM cisplatin under long‐term isotonic (C) or 200 μM cisplatin in short‐term hypo‐ and isotonic (D) conditions as function of time (n = 3). Similar results were obtained in HCT116 cells (Fig EV4).Cisplatin uptake (200 μM) into HEK cells of indicated genotypes. LRRC8(B,C,E) −/− and LRRC8(B,D,E) −/− cells express only LRRC8A and LRRC8D, and LRRC8A and LRRC8C, respectively (n = 3 for WT, 6 for LRRC8D −/−, 9 for LRRC8A −/− in E; n = 3 in F).Mean current densities of maximally activated ICl,vol at −80 mV. The number of cells is indicated for each column.Ratio of LRRC8‐dependent swelling‐activated cisplatin uptake (60 min) to mean ICl,vol (as in G) as function of genotype.Data information: Data are presented as mean ± SEM. *P < 0.05; **P < 0.01; and ***P < 0.001 (for (C, E) compared to HAP1−/− cells). Dotted lines in (D) and (F) highlight that there is no significant difference in isotonic 60‐min uptake between the genotypes.
Figure 7Cisplatin-DMSO inhibition of ICl,vol depends on the LRRC8D subunitUpper panel, example current traces (as in Fig 3A) of fully activated ICl,vol in HEK cells exposed to hypotonic solution containing vehicle (0.3% DMSO) or 200 μM cisplatin in 0.3% DMSO. Dashed lines indicate zero current. Lower panel, ICl,vol current densities (at −100 mV and 100 mV) of WT HEK cells treated with different cisplatin concentrations.No effect of 200 μM cisplatin/DMSO on ICl,vol in LRRC8D−/− HEK cells.Data information: Data are presented as mean ± SEM; the number of experiments is given for each bar; *P < 0.05; **P < 0.01.
Figure 8Swelling‐induced efflux of 3[H]‐taurine from WT and LRRC8D −/− HEK cells (A, B) and partial rescue by transient transfection of LRRC8D (B). Rescue is incomplete due to low transfection/expression efficiency of LRRC8D (Voss et al, 2014).C-ESwelling‐induced efflux of 3[H]‐taurine from LRRC8(B,C,E) −/− (C), LRRC8(B,D,E) −/− (D), or LRRC8(B,C,D) −/− (E) HEK cells compared to WT cells.FCisplatin‐induced taurine efflux (over 30 min) from WT and HAP1−/−HEK cells after preincubation with 200 μM cisplatin for 4 or 12 h, or without cisplatin (0 h). Carbenoxolone (CBX; 100 μM) blocks taurine efflux from WT cells treated for 12 h with cisplatin, excluding taurine flux through unspecific leaks.Data information: Data are presented as mean ± SEM. Dashed lines in (A-E) represent isotonic efflux from the same individual experiments. Red arrows, change to hypotonic solution. Each panel shows the average of two independent experiments (n = 4 for each experiments, i.e., n = 8 in total).
Figure 1(C). The table summarises the number and the sex of the embryos of the indicated genotypes, recovered from timed matings of Rif1+/- x Rif1+/- mice, in a C57BL/6J genetic background. The day of gestation (E) is indicated.
Figure 2(C). Time course analysis of Xist RNA expression by RT-qPCR during EB differentiation of Rif1+/+ (Rif1+/+ +OHT) and Rif1-/- (Rif1F/F +OHT) cells at the indicated timepoints. Rif1+/+ (solid line) and Rif1-/- (dashed line), female (black) and male (grey). Data are presented as mean ± standard deviation from three (female lines) or two (male lines) independent experiments. Statistical significance was determined using 2-way ANOVA comparing female Rif1+/+ to female Rif1-/- cell lines (****p ≤ 0.0001). Xist RT-primers Xist ex7 F and R were used. Values are normalised to a geometric mean consisting of the expression of Gapdh, Ubiquitin and β-Actin.
Figure 3RIF1 association with the Xist promoter assessed by ChIP-qPCR in two independent Rif1+/+ (Rif1+/+ +OHT, black) and two Rif1-/- (Rif1F/F +OHT, grey) female cell lines, in ESCs (A) P1 and P2 indicate the two Xist promoters, 5' indicates a region 2 kb upstream of Xist TSS. Inter1 and 2 are two intergenic regions that serve as negative controls. Peak and cRAD represent two previously identified regions of RIF1 association (positive control). See Appendix Fig. S2C for primer positions within Xist. Mean ±standard deviation from 3 independent experiments (A) and 2 independent experiments (B). p calculated by Student's two-tailed, paired t test comparing RIF1 association in Rif1+/+ cells on Xist P2 and P1 versus 5'. *p ≤ 0.05, ****p ≤ 0.0001 and ns=not significant.(D). Association of RIF1 with Xist P2 in the Fa2L cells (black) and a wild type female mESC line (grey), also harbouring one castaneus and one 129 X chromosome. Allele-specific ChIP-qPCR primers were used, cast indicates association with the castaneus Xist P2 and 129 indicates association with the 129 Xist P2. Enrichments are presented relative to input DNA. Mean ±standard deviation from 3 independent experiments. p calculated by Student's two-tailed, paired t test comparing RIF1 association with the castaneus and with the 129 X chromosome Xist P2, ***p ≤ 0.001. Below, schematic of the Xist/Tsix alleles in the Fa2L undifferentiated cells. (E). Association of RIF1 with Xist P2 in the Fa2L cells (black) upon differentiation. The analysis was performed as in (D). Mean ±standard deviation from 3 independent experiments. p calculated by Student's two-tailed, paired t test comparing RIF1 association with the castaneus and with the 129 X chromosome Xist P2. *p ≤ 0.05. Below, schematic of the Xist/Tsix alleles in the Fa2L differentiated cells. ((F). Quantification by ChIP-qPCR of RIF1 association with the indicated regions in the Fa2L cell line, following treatment with DMSO only (black) or flavopiridol (grey). Primers as in (E). (G). Same as in (F) but for a wild type female mESC line. All enrichments are presented relative to input DNA. Mean ±standard deviation from three (F) and two (G) independent experiments are presented. Statistical significance was determined using Student's two-tailed, paired t test (*p ≤ 0.05, **p ≤ 0.01 and ns=not significant).
Figure 4(A). Time course analysis of Xist expression by RT-qPCR during EB differentiation of female mESCs following knock down of Luciferase (Control, black) and Kap1 (Kap1 KD, grey), at the indicated timepoints. Data are presented as mean ±standard deviation from three independent experiments. Statistical significance was determined using 2-way ANOVA. Xist primers Xist ex3 F and Xist ex4 R were used. Normalisation was performed using a geometric mean consisting of the expression of Rplp0, Ubiquitin and Sdha. (*p ≤ 0.05). (B). RT-qPCR analysis of Xist expression levels during differentiation of the Fa2L cells, following expression of shRNA against Luciferase (Control, black) and Kap1 (Kap1 KD, grey), at the indicated timepoints. Mean ±standard deviation from a minimum of three independent experiments is presented. 2-way ANOVA was used to determine statistical significance. ns=not significant. ((C). Tsix RNA levels in female mESCs infected with shRNA directed against Luciferase (Control, black) and KAP1 (Kap1 KD, grey), during differentiation. Mean ±standard deviation from four independent experiments are shown. Statistical significance was determined using 2-way ANOVA. (**p ≤ 0.01). Values have first been normalised to a geometric mean consisting of the expression of Rplp0, Ubiquitin and Sdha.(D). RNA FISH analysis of Tsix expression during differentiation of female mESCs expressing an shRNA directed against Luciferase (control) or against Kap1 (Kap1 KD). Left: Cells with no (0, red), one (1, grey) and two or more (2, black) Tsix foci were counted in two independent experiments, shown averaged. Statistical significance was determined by χ2. A minimum of 110 cells were counted per each time point for each line (**p ≤ 0.01, ****p ≤ 0.0001). Right: examples of cells with one (top) or two (central and bottom) Tsix FISH signals. Scale bars: 5μm.
Figure 5(A). ChIP-qPCR analysis of KAP1 association with the indicated sites in wild type female mESCs (Rif1+/+, same as used in Fig. 3B but without OHT) and during early differentiation. ZFP629 is a well-characterised KAP1 associated region (positive control). Xite A and C indicate two regions within the Tsix enhancer Xite, Tsix region 1 indicates Tsix major promoter, Tsix region 2 indicates the Dxpas34 region, Tsix region 3 indicates a region slightly downstream of the Dxpas34 region. P1 and P2 indicate the two Xist promoters, 5' indicates a region 2 kb upstream of Xist TSS. Inter1 is an intergenic region. See Appendix Fig. S2C for the positions of the primers within Xist and Tsix. The data are presented as mean ±standard deviation from three (2d EB and 1d EB) and two (ESCs) independent experiments. Statistical significance was calculated by Student's two-tailed unpaired t test comparing RIF1 association to Xist P2 and P1 in 2d EB versus 1d EB (*p ≤ 0.05 and ns=not significant). In the inset, Tsix RNA levels were quantified by RT-qPCR during the differentiation of wild type female ESCs shown in Fig. 4F. The average of two experiments is shown. Tsix values are normalised to a geometric mean consisting of the expression of Rplp0, Ubiquitin and Sdha. Error bars indicate standard deviations.(C). KAP1 association with Xist promoter in two independent Rif1+/+ (Rif1+/+ +OHT, black) and two Rif1-/- (Rif1F/F +OHT, grey) female mESC cell lines. Ezr is an additional region known to be associated with KAP1 in mESCs. Enrichments are presented relative to input DNA. Mean ±standard deviation from a minimum of three independent experiments per cell line are displayed. Statistical significance was determined using Student's two-tailed, unpaired t test comparing the KAP1 association with Xist P2 and P1 in Rif1+/+ versus Rif1-/- cells (**p ≤ 0.01).
Figure 6(C). HEXIM1 association with the Xist promoter (P1 and P2) in two independent Rif1+/+ and two Rif1-/- mESC lines, analysed by ChIP-qPCR. Tkbp1 and Myf5 are two control regions. As in the case of KAP1 association, deletion of Rif1 induces accumulation of HEXIM1 on Xist P1 and P2. Average ±standard deviation of two independent experiments. p values were calculated by two-tailed, unpaired, equal variance t-test. (**p ≤ 0.01, ***p ≤ 0.001).(E). Upon infection of Fa2L cells with shRNA against Kap1 or control, against Luciferase, HEXIM1 association with P2 was analysed, on both alleles, by ChIP-qPCR. As in the case of KAP1, HEXIM1 shows preferential association with Xist P2 on the Cast allele (future Xa, Control). The association is lost upon knock down of Kap1 (Kap1 KD). Myf5 serves as a negative and Neat1 as a positive control regions.
Figure 1(B) Whole mount low magnification images show marking in interscale (K10+) vs scale regions,C) interscale (K10+) vs scale regions, quantified in (C). Error bars are SDs. The data shown is average of 4 different timepoints (2 week, 3month, 6month and 1 year) each consisting of 2 mice and 5 images per mice. p-value is calculated for difference between scale vs interscale using Student's t-test (mixed Anova).(D) Tail skin sections immunofluorescence stained for interscale (K10) or scale (K31) supra-basal markers, and basal markers. Scale bar =20 µm or 50 µm as indicated on the panels. Asterisks indicate that photoshop was used to erase the autoflourescent cornified layer. Dashed lines indicate ROI in basement membrane.(E) Quantification of images like those in (D) with n=3-4 mice and 6-7 images per mouse showing preferred-gene enrichment. Error bars are SDs. P-values calculated by Student's t-test (mixed effect model)