File Name : supp.figure_1.pdf Caption : supplementary figure 1. normalization of signal to exclude kcs. a) sample correlation before kcs normalization. b) principal component analysis indicating the kcs contamination score. c) heatmap showing the differential expression level of kcs top marker genes across samples. untreated donors 1 and 4 showed a substantial contamination. d) volcano plot showing the differential gene expression of untreated control vs tgfb1-treated before and after normalization. File Name : supp.figure_2.pdf Caption : supplementary figure 2. additional pathway analysis. go (a) and kegg (b) show similar pathway enrichment, mostly associated with extracellular organization, adhesion, collagen metabolism, integrin signaling and organ development. File Name : supp.figure_3.pdf Caption : supplementary figure 3. paired t-test of fibrosis-associated genes (tgfb vs tv, donors 1, 3, and 4). col6a3 was significantly downregulated under tv, while other genes showed non-significant downward trends. donor 2 was excluded due to lack of response to tgfb. File Name : supp.figure_4.pdf Caption : supplementary figure 4. bulk rna-sequencing analysis comparing selected tgfb1-treated vs tgfb1+vtp-treated donors. heatmap showing differentially expressed fibrosis-associated genes between tgfb and tv (verteporfin) treatment across three donors. the data highlight donor-specific transcriptional responses to tgfb-induced fibrosis and verteporfin treatment. File Name : supp.figure_5.pdf Caption : supplementary figure 5. quantification of mechanosignaling mediators localized in the nucleus. if was used to assess the nuclear localization of yap, pyap, psmad2, pstat3, pfak, psting, and to calculate the yap/pyap ratio (e-j) (n>1700 per condition). the statistical analysis was performed using one-way anova. File Name : supp.figure_6.pdf Caption : supplementary figure 6. representative load and indentation depth versus time curves with association load relaxation experimental and poroviscoelastic model fit. solid black line represents the time-dependent force response of the sample during an indentation test. dotted blue line represents the indentation depth of the probe tip. force data in the green boxed region was used to compute material properties with the combined poroelastic viscoelastic (pve) model, which is represented by the red solid line. File Name : supp.figure_7.pdf Caption : supplementary figure 7. benchmarking of microindentation using neonatal foreskin. a) elastic modulus and b) permeability were tested on foreskins from four different donors.