Abstract
Fibrosis occurs when collagen deposition and fibroblast proliferation replace healthy tissue. Red light (RL) may improve skin fibrosis via photobiomodulation, the process by which photosensitive chromophores in cells absorb visible or near-infrared light and undergo photophysical reactions. Our previous research demonstrated that high fluence RL reduces fibroblast proliferation, collagen deposition, and migration. Despite the identification of several cellular mechanisms underpinning RL phototherapy, little is known about the transcriptional changes that lead to anti-fibrotic cellular responses. Herein, RNA sequencing was performed on human dermal fibroblasts treated with RL phototherapy. Pathway enrichment and transcription factor analysis revealed regulation of extracellular matrices, proliferation, and cellular responses to oxygen-containing compounds following RL phototherapy. Specifically, RL phototherapy increased the expression of MMP1, which codes for matrix metalloproteinase-1 (MMP-1) and is responsible for remodeling extracellular collagen. Differential regulation of MMP1 was confirmed with RT-qPCR and ELISA. Additionally, RL upregulated PRSS35, which has not been previously associated with skin activity, but has known anti-fibrotic functions. Our results suggest that RL may benefit patients by altering fibrotic gene expression.
Similar content being viewed by others
Introduction
As part of the healing process, fibroblasts proliferate, differentiate, and increase collagen production1,2. Fibrosis occurs when collagen deposition and fibroblast proliferation replace healthy tissue due to impaired wound healing, immune dysfunction, or iatrogenic causes3. Fibrosis can affect most organs, including the skin, liver, ovaries, and lungs3,4,5,6,7. Keloids, hypertrophic scars, scleroderma, and radiation dermatitis are among the many manifestations of skin fibrosis that significantly burden patients' quality of life8,9. Current treatment options include 5-fluorouracil, radiation therapy, immunomodulators, and surgery, but fibrosis is often recalcitrant to existing approaches10. Ultraviolet-A1 (UVA1; 340–400 nm) phototherapy is an existing treatment for fibrosis but has a limited depth of penetration (less than 150 μm) and carries a risk for non-melanoma skin cancer and aging11,12,13,14. Red light (RL) may be a potentially safer alternative to UVA1 with a greater depth of penetration (6–50 mm)15,16,17,18. RL may improve skin fibrosis via photobiomodulation, the process by which photosensitive chromophores in skin cells absorb visible or near-infrared light and undergo photophysical reactions19,20. Photobiomodulation was previously known as low-level light therapy, but higher fluences (i.e., a higher amount of light energy delivered over the treatment period) may also induce cellular changes.
Photobiomodulation has a multiphasic response in which lower fluences are stimulatory, while higher fluences are inhibitory or cytotoxic21. Much of the existing research has focused on the effects of visible or near-infrared light at fluences lower than 150 J/cm2 to stimulate cell growth for skin rejuvenation or hair growth22. However, our laboratory team has studied the inhibitory properties of red light greater or equal to 320 J/cm2 to treat pathological hyperproliferative processes, including fibrosis23,24,25. Our previous research demonstrated that photobiomodulation using high fluence RL reduces human dermal fibroblast (HDF) proliferation, collagen deposition, and migration in vitro23,24,25. We sought to investigate changes in gene expression following RL phototherapy. We performed high-throughput RNA sequencing (RNA-Seq) to identify genes and pathways associated with RL. RNA-Seq allows for unbiased discovery of therapeutic targets and whole transcriptome gene expression analysis.
Results
Transcriptomic profiling of human dermal fibroblasts
Principal component analysis (PCA) demonstrated that samples segregate according to donors, which is characteristic of human subjects’ analysis (Table 1 and Fig. 1A). The heat map of the top 30 genes with maximum variance values, calculated for all samples, shows similar segregation by donor line (Fig. 1B). When PCA analysis was separately performed for each donor, samples collected at the 0-h timepoint clustered apart from the 4, 12, and 24 h (Fig. 1C,D). 640 J/cm2 treated samples clustered separately from control for 3 out of 4 donor lines (Fig. 1D). It should be noted that changes between treatment and control were found to be in the same direction of gene expression dimensionality reduced space.
To account for interpersonal variations in gene expression, DESeq2, capable of paired analysis, was used to identify differentially expressed genes (DEGs)26. Analysis of the samples revealed 859 DEGs following RL with a twofold change in expression and false discovery rate (FDR) < 0.05 across all timepoints. The complete list of DEGs from the analysis with fold change and FDR can be found in the supplemental dataset.
RL induces a temporal change in gene expression
At 0 h after RL, 147 out of 191 and 205 out of 239 DEGs were downregulated for both 320 J/cm2 and 640 J/cm2, respectively (Figs. 2A,B). Four hours post-treatment, the trend reversed for HDFs treated with 640 J/cm2 RL as 304 DEGs were upregulated (Fig. 2B). Four hours was the temporal apex of DEG upregulation (Fig. 2B). 640 J/cm2 RL may cause immediate downregulation of transcription and subsequent upregulatory cellular compensations. At 320 J/cm2, the number of downregulated DEGs decreased such that up and downregulation were approximately equal at 4 (28 up/17 down), 12 (19 up/18 down), and 24 h post-treatment (22 up/19 down) (Fig. 2A). These results indicate that the greatest changes in gene expression occur within 4 h and is dose-dependent. The differential expression of MMP1 (fold change: 2.36, p-value: 9.86 × 10–14, FDR: 5.02 × 10–11 at 24 h) is shown as a representative gene of interest involved in extracellular matrix organization (Fig. 2C). Volcano plots of fold-changes versus significance demonstrate that some highly expressed genes did not have an FDR greater than 0.05, likely due to a small sample size of fibroblast donors (Fig. 2D,E). 640 J/cm2 induced significant enrichment of the upregulation of cell migration (GO: 0030334, p-value: 4.39 × 10–7, q-value: 2.24 × 10–3), extracellular matrix organization (GO: 0030198, p-value: 1.96 × 10–5, q-value: 9.11 × 10–3), cellular response to oxygen-containing compounds (GO:1901701, p-value: 6.53 × 10–5, q-value: 1.67 × 10–2), and regulation of cell proliferation (GO: 0042127, p-value: 2.49 × 10–4, q-value: 4.23 × 10–2) by GO analysis via Fisher exact test (Fig. 2F). Complete GO and KEGG pathway analyses are provided in Figures S1-4. Transcription factor (TF) associated with oxidative stress (e.g., RELA) and fibrosis (SMAD3, FOSL1, FOSL2, JUN, and JUNB) are predicted to regulate the cellular response to RL phototherapy (Fig. 2G)27,28,29,30,31.
Identification of genes associated with fibrosis and oxidative stress
We aimed to investigate transcriptional changes induced by RL in HDFs. To determine if RL photobiomodulation induces differential expression beneficial for the treatment of fibrotic skin disease, candidate genes associated with fibrosis were further analyzed. The following DEGs were identified as related to fibrotic processes: DOCK4, WNT2, MMP1, PSAT1, MIR21, MIR145, and BTBD11 (Figure S5). PRSS35, RANBP3L, and MECOM had higher than tenfold changes in expression and are associated with fibrosis. The data was reanalyzed without a fold-change cutoff, and we observed differential expression of SMAD3, SMAD4, and SMAD7 (Figure S6). Additionally, we searched for candidate gene transcripts related to reactive oxygen species (ROS) regulation. ROS has been associated with photobiomodulation’s mechanism and changes in apoptotic, cytoprotective, and fibrotic pathways22,25. Identified ROS related DEGS were AKR1B1, NPAS2, AKR1C1, RCAN1, MRPS6, XDH, ADM2, NCF2, DDIT4L, SLC5A3, and GUCY1A2 (Figure S7).
Confirmation of differential expression and ontology
Verification of MMP1 protein and gene expression was performed (Figs. 3A–D). MMP1 produces matrix metalloproteinase-1 (MMP-1), also known as collagenase, and is capable of enzymatically degrading the collagen found in fibrosis32,33. Changes in gene expression between RT-qPCR and RNA-Seq were in the same direction and had similar fold changes at the 0, 4, 12, and 24 h time points (Figs. 3A). A Pearson correlation statistical test was performed to determine consistency between RT-qPCR and RNA-Seq found an R = 0.98 and a significant P-value of 0.015 (Fig. 3B). At 4, 12, and 24 post-RL irradiation, the supernatant was collected from the control and RL-treated sample. MMP-1 protein expression was quantified using ELISA. By 24 h post-treatment, RL treated samples released significantly more MMP-1 compared to control (Fig. 3C). A Pearson correlation statistical test to determine consistency between MMP-1 ELISA and RNA-Seq differential expression was performed. RNA-Seq and ELISA were highly (R = 0.98), but not significantly (P = 0.14), correlated (Fig. 3D).
We previously found that RL increased intracellular ROS generation and inhibited cell proliferation in HDF donor 125. GO pathway analysis revealed enrichment of cellular responses to oxygen containing compounds and regulation of cell proliferation in all 4 donors treated with 640 J/cm2 (Fig. 2F). We confirmed that ROS increases and cell count decreases in donors 1–4 when treated with RL. 320 and 640 J/cm2 RL increased ROS in a dose dependent manner at 0 h post-treatment (Fig. 3E). 640 J/cm2 RL decreased cell count by 48 h post irradiation (Fig. 3F).
Discussion
We present whole transcriptome gene expression analysis of HDFs at 0, 4, 12, and 24 h after treatment with 320 or 640 J/cm2 RL. There were more DEGs within 4 h of treatment than by 12 and 24 h after RL treatment (Fig. 2A). This is consistent with our previous finding that RL mediated effects on HDF migration dissipated within 12 h25. MMP1 expression was confirmed with RT-qPCR and ELISA for MMP1. Additionally, PRSS35, which produces a serine protease with collagen 1 degrading potential, was found to have greater than 30 fold increased expression at 4, 12, and 24 h after RL treatment, and thus may impart anti-fibrotic effects via collagen degradation34. This study provides foundational insights for future investigation into photobiomodulation and fibrosis.
We have previously shown that RL directly decreases collagen protein expression25. GO pathway analysis demonstrated enrichment of genes related to the organization of the extracellular matrix (Fig. 2F). As a result, we sought to investigate genes and pathways associated with anti-fibrotic activity. Other researchers have found that skin MMP-1 expression and secretion increases in response to UV and visible light35,36. Li et al. found that 3 J/cm2 RL increased the expression of multiple MMPs, including MMP-137. We confirmed similar increases in MMP-1 expression and secretion in HDF treated with high fluence RL (Fig. 3). This is significant as RL is not associated with skin cancer and aging like UVA phototherapy.
PRSS35 is a serine protease that may have collagen-1 degrading function34. PRSS35 has been previously linked to gonadal function, but human epididymis protein 4 (HE4), an inhibitor of PRSS35, was highly expressed in fibrotic kidneys34,38,39. We observed a greater than 30-fold increase of PRSS35 expression in the RL treated HDFs, suggesting that PRSS35 induction by RL may reduce collagen. The increased expression of collagen-1 degrading enzymes may highlight the anti-fibrotic mechanisms of RL through the degradation of extracellular collagen.
To further characterize the potential of RL in fibrosis, DEGs in the TGF-β signaling pathway were analyzed (Figure S3).40. We have previously shown that 640 J/cm2 RL decreases SMAD-2 phosphorylation in TGF-β1 induced HDFs within 4 h of irradiation41. Phosphorylated SMAD 2/3 translocate to the nucleus and increase the expression of collagen42,43,44. In our RNA-Seq analysis, HDFs were not TGF-β induced, but there was nevertheless downregulation of profibrotic SMAD3 at 4, 12, and 24 h post-RL treatment. TF analysis confirmed that SMAD3 was likely involved with regulation of cellular activity following RL phototherapy (Fig. 2G). SMAD-4 is pro-fibrotic, but its expression was increased at 4 h post-RL treatment45. SMAD-7 has anti-fibrotic properties and was slightly downregulated in our RNA-seq data46,47.
There were other highly expressed DEGs involved with SMAD regulation. RANBP3L produces a protein that acts as a nuclear factor that can export SMAD-1 and other proteins known to have anti-fibrotic properties as part of the TGF-β pathway48. MECOM produces a protein called EVI-1, which acts as a transcriptional regulator that can inhibit SMAD protein activity49,50.
We previously performed RNA-Seq for miRNA from HDF donor line 1 and found that miRNA-21, miRNA-23, and miRNA-31 were decreased, while miRNA-29, miRNA-196a, and let-7a were increased51. These microRNAs have been identified as mediators of skin fibrosis52. We repeated our miRNA analysis with HDF donor lines 1–4 and confirmed that miRNA-21 (Mir21) expression was significantly downregulated following RL irradiation (Figure S2). miRNA-21 regulates TGF-β/SMAD signaling, and decreased expression of miR-21 is anti-fibrotic52. Additionally, miRNA-145 expression was decreased. miRNA-145 is increased in hypertrophic scars and has been identified as a therapeutic target for anti-fibrotic therapies53.
FOS/JUN family of proteins (i.e., FOSL1, FOSL2, JUNB, and JUN) were identified as TF regulating the response of fibroblasts to RL phototherapy. JUN proteins form homodimers or heterodimers with FOS proteins to increase the expression of AP128,30. AP1 can regulate cell cycle progression and extracellular matrix organization28. c-Jun (JUN) is phosphorylated by c-Jun N-terminal kinases (JNKs) in response to cellular stress, growth factors, or cytokines28. c-Jun is important in IL-17 mediated production of MMP-1 in HDFs28,29. Similarly, inhibition of JNK in HDFs prevented the upregulation of MMP-3 and MMP-1 in response to UVB28,31. Fibrotic responses to FOS/JUN activity can change based on cell type and conditions28,29,30. c-Jun is upregulated in the skin of patients with systemic sclerosis and smooth muscle actin positive (SMA +) HDFs30. In fibrotic mouse models, phosphorylation of c-Jun is associated with profibrotic cellular responses via activation of AKT30. However, we previously found that increased AKT phosphorylation by RL in HDFs inhibited migration and was associated with decreased collagen deposition25,28.
RL and near-infrared radiation stimulate cytochrome C oxidase in the mitochondria, altering mitochondrial membrane potential and increasing intracellular ATP and free-radical ROS20,22. We confirmed that ROS increases following irradiation with 320 and 640 J/cm2 RL using flow cytometry (Fig. 3E). ROS may alter the activity of fibrotic pathways, including TGF-β, mTOR, and AKT22,25. GO analysis demonstrated the enrichment of cellular responses to oxygen-containing compounds (GO: 1901701, p-value: 6.53 × 10–5, q-value: 1.67 × 10–2).
Transcription factor co-regulatory network analysis indicated that RELA, which codes for the p65 subunit of NF-κB, is predicted to regulate gene expression in RL-treated HDFs (Fig. 2G)27,54. The heterodimer of RELA and p50 is the most abundant form of NF- κB27,54. NF-κB is involved with inflammation and cellular responses to stress, including ROS27,54. NF-κB has been previously implicated in photobiomodulatory mechanism since low fluences of near-infrared light-activated NF-κB in mouse embryonic fibroblasts55.
Our primary objective was the discovery of RL induced transcriptome modulation. Two previous studies by Kim et al. and Li et al. examined the effects of RL on HDF transcription using RNA-Seq and demonstrated similar regulation of genes/pathways related to MMPs, FOS, JUN, NF-kB, SMAD1/7, oxidative stress, and inflammation37,56. Our comprehensive data set is a strength of this study and may serve as a reference for future research. This objective was limited by sparse prior research on photobiomodulation transcriptomics. This limitation restricts our enrichment analysis because photobiomodulation pathways may be under-represented in the major databases such as KEGG and GO. Our and others’ research may contribute to the identification of photobiomodulation pathways37,56. Another limitation of our study was that we did not use HDFs isolated from fibrotic tissue. However, fibrotic HDFs may lose their fibrotic phenotype after being removed from their in vivo fibrotic niche. Non-fibrotic HDFs may be a good analog as RL has similar anti-fibrotic effects in normal and keloid-derived HDFs23,24. In future research, the transcriptomic effects on fibrotic skin from reconstructed three-dimensional, animal, and clinical models should be assessed57,58. In vivo and tissue culture models may respond to RL phototherapy differently.
In conclusion, we identified several genes that may contribute to the mechanism of RL treatment of fibrosis. MMP1 is a critical mediator of fibrotic disease that was modulated by RL treatment. We identified PRSS35 as a potential mechanism of RL anti-fibrosis due to its 30-fold increased expression. PRSS35 could be the focus of future photobiomodulation studies because of its previously limited characterization and profoundly differential expression in RL treated samples. Our results suggest that RL has the potential to benefit patients with fibrosis by altering gene expression.
Methods
Cell culture
Normal HDFs were obtained from the American Type Culture Collection (CRL-2617, CRL-2697, and CRL-2796) and Coriell Biorepository (AG13145). HDF Donor samples were used per relevant guidelines and regulations. HDFs were sub-cultured in DMEM (Invitrogen; Carlsbad, CA) supplemented with 10% FBS (R&D Systems; Minneapolis, MN) and 1% antibiotic/antimycotic (Invitrogen). Cells were maintained in a humidified incubator with 5% CO2 and 20% O2. RNA was collected from 35 mm tissue culture dishes (Corning, Corning, NY) that were initially seeded at low confluency (2 × 104 cells total; 4,000 cells per 1.77 cm2 surface area) between passages four and seven59. Twenty-four hours after seeding, samples were treated with RL, and RNA was collected at 0, 4, 12, and 24 h time points.
HDF donors
RNA-seq was performed with total RNA samples collected from four commercially available HDF cultures obtained from three different anatomical sites: two from the abdomen, one from the forearm, and one from the lower leg (Table 1).
RL treatment
HDFs were treated with RL as previously described23,24,25. Briefly, an LED unit (Omnilux; Globalmed technologies, Napa, CA) was utilized for all experiments. The LEDs have a rectangular aperture with dimensions 4.7 cm × 6.1 cm and emit visible red light at 633 ± 30 nm wavelength in the electromagnetic spectrum. The light has a power density of 872.6 W/m2 at room temperature and 10 mm from the bottom of the plastic culture dish. Cell cultures were treated to 320 J/cm2 or 640 J/cm2 (3667 s for 320 J/cm2 and 7334 s for 640 J/cm2) of RL at approximately 34 ˚C. During RL treatments, the cells were exposed to environmental 20% O2 and 412 parts per million CO2 concentrations outside of the incubator. Controls were placed on plate warmers set to 34 ˚C and protected from light with aluminum foil to match RL treated samples' environmental conditions.
RNA isolation
Total RNA from HDFs was collected at 0, 4, 12, and 24 h after RL treatment. The miRNeasy (Qiagen; Germantown, MD) kit was used to isolate RNA from cell cultures following the manufacturer’s suggested protocol. To briefly summarize, Qiazol reagent (Qiagen) was used to lyse cells, followed by chloroform extraction. The aqueous layer was obtained and mixed with 100% ethanol (Sigma; St. Louis). Spin columns further aided the separation of RNA and impurities from samples. Samples were treated with RNase-free DNase (Qiagen) to ensure no genomic DNA contamination. Finally, RNase free water was used to elute sample RNA. All samples had RNA quality assessed by Tapestation 2200 (Agilent Technologies; Santa Clara, CA). All samples had RNA integrity number values of 9.9 or 10.0.
Library preparation and sequencing
RNA integrity was measured using the RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent Technologies). Libraries were built with NEBNext Ultra Directional RNA Library Prep Kit for Illumina (New England Biolabs; Ipswich, MA). Sequencing of all samples was performed by pooling all of our indexed samples and putting equal amounts of the pooled sample into each lane of a flow cell performed on Hi Seq X (Illumina; San Diego, CA).
Mapping and identification of differentially expressed genes
Sequencing reads were mapped to the UCSC human reference genome (GRCh37/hg19), and the following read counts were evaluated by STAR (version 2.5.2)60. Gene expression level normalization and differential expression analysis were carried out by DESeq2 (version 1.6.3) bioconductor R package26. To compare samples before and after treatment for different cell lines, a multifactor design was used applying DESeq2 controlling for the effect of cell line difference. Differential expression p-values were corrected for multiple testing using the false discovery rate (FDR) method. Enrichment analysis was performed with Enrichr61,62.
RT-qPCR
RT-qPCR experiments used materials and equipment from Bio-Rad (Hercules, CA). 100 ng of RNA was synthesized into cDNA with the iScript reverse transcription kit using a C1000 thermal cycler. RT-qPCR was performed with 1 ng of cDNA on the BioRad CFX96 using SYBR green.
ELISA
At 4, 12, and 24 h following RL irradiation, we quantified total human MMP-1 in collected HDF culture media using ELISA (R&D Systems) according to the manufacturer's guidelines. For each sample, the concentration of released MMP-1 was indexed to the total intracellular protein. We quantified collected intracellular protein using Bradford reagent (Bio-Rad). Optical density was measured for ELISA and protein concentration using a 96-well plate reader (Synergy 2, Biotek; Winooski, VT). For each donor, the experiment was performed in technical triplicate. Relative MMP-1 expression (RL/control) was pooled from the 4 donor lines and compared to a hypothetical mean of 1, indicating no difference between RL and control, using a one sample T-Test. P < 0.05 was considered significant. ** indicated p < 0.01.
Cell count
Cell counts were assessed using crystal violet63. Following treatment with RL, experimental and control samples were placed in a humidified incubator for 48 h. Cells were fixed with 4% formaldehyde (Sigma) and stained with 0.1% crystal violet (Thermo-fisher Scientific; Waltham, MA). 10% acetic acid (Thermo-fisher Scientific) was used to elute the crystal violet. Optical density of eluted crystal violet was quantified with a plate reader at 595-nm. For each donor, the experiment was performed with a technical repeat of n = 3–5. Relative counts (RL/control) were pooled from the 4 donor lines and compared to a hypothetical mean of 1 (indicating no difference between RL and control), using a one sample T-Test. P < 0.05 was considered significant (*).
Free radical reactive oxygen species generation
For free radical ROS generation, HDFs were assayed using dihydrorhodamine-123 (DHR-123; Thermo-fisher Scientific). Cells were irradiated with RL and then treated with DHR-123 for 30 min. Non-fluorescent DHR-123 converts to fluorescent rhodamine-123 in the presence of ROS. RL treated and control cells were detached with 0.25% trypsin EDTA (Thermo-fisher Scientific), collected, and analyzed with flow cytometry (Fortessa; BD; San Jose, CA). Intracellular ROS generation was assessed immediately following irradiation (0 hours). Positive control cells were treated with 0.6 mM hydrogen peroxide (Thermo-fisher Scientific) for 30 min. ROS was quantified as the median fluorescent intensity (MFI) of rhodamine-123 using FlowJo Software (BD). For each donor, the experiment was performed with a technical repeat of n = 4 or 5. Relative MFIs of rhodamine-123 (RL/control) were pooled from the 4 donor lines and compared to a hypothetical mean of 1 (indicating no difference between RL and control), using a one sample T-Test. P < 0.05 was considered significant (*).
Data availability
The datasets generated during and/or analysed during the current study are included within manuscript or available from the corresponding author on reasonable request.
References
Darby, I. A. & Hewitson, T. D. Fibroblast differentiation in wound healing and fibrosis. Int. Rev. Cytol. 257, 143–179 (2007).
Werner, S., Krieg, T. & Smola, H. Keratinocyte–fibroblast interactions in wound healing. J. Investig. Dermatol. 127, 998–1008 (2007).
Andrews, J. P., Marttala, J., Macarak, E., Rosenbloom, J. & Uitto, J. Keloids: The paradigm of skin fibrosis—Pathomechanisms and treatment. Matrix Biol. J. Int. Soc. Matrix Biol. 51, 37–46. https://doi.org/10.1016/j.matbio.2016.01.013 (2016).
Nogueira, A., Pires, M. J. & Oliveira, P. A. Pathophysiological mechanisms of renal fibrosis: A review of animal models and therapeutic strategies. Vivo (Athens, Greece) 31, 1–22. https://doi.org/10.21873/invivo.11019 (2017).
Toosi, A. E. Liver fibrosis: Causes and methods of assessment, a review. Roman. J. Intern. Med. (Rev. Roumaine Med. Intern.) 53, 304–314, https://doi.org/10.1515/rjim-2015-0039 (2015).
Zhou, F., Shi, L. B. & Zhang, S. Y. Ovarian fibrosis: A phenomenon of concern. Chin. Med. J. 130, 365–371. https://doi.org/10.4103/0366-6999.198931 (2017).
Karsdal, M. A. et al. The good and the bad collagens of fibrosis—Their role in signaling and organ function. Adv. Drug Deliv. Rev. 121, 43–56. https://doi.org/10.1016/j.addr.2017.07.014 (2017).
Sierakowska, M. et al. Factors associated with quality of life in systemic sclerosis: a cross-sectional study. Qual. Life Res. 28, 3347–3354. https://doi.org/10.1007/s11136-019-02284-9 (2019).
Bock, O., Schmid-Ott, G., Malewski, P. & Mrowietz, U. Quality of life of patients with keloid and hypertrophic scarring. Arch. Dermatol. Res. 297, 433–438. https://doi.org/10.1007/s00403-006-0651-7 (2006).
Berman, B., Maderal, A. & Raphael, B. Keloids and hypertrophic scars: Pathophysiology, classification, and treatment. Dermatol. Surg. 43 Suppl 1, S3–S18, https://doi.org/10.1097/dss.0000000000000819 (2017).
Meinhardt, M., Krebs, R., Anders, A., Heinrich, U. & Tronnier, H. Wavelength-dependent penetration depths of ultraviolet radiation in human skin. J. Biomed. Opt. 13, 044030. https://doi.org/10.1117/1.2957970 (2008).
Anderson, R. R. & Parrish, J. A. The optics of human skin. J. Invest. Dermatol. 77, 13–19. https://doi.org/10.1111/1523-1747.ep12479191 (1981).
Gambichler, T. & Schmitz, L. Ultraviolet al phototherapy for fibrosing conditions. Front. Med. 5, 237. https://doi.org/10.3389/fmed.2018.00237 (2018).
Maiorino, A., De Simone, C., Perino, F., Caldarola, G. & Peris, K. Melanoma and non-melanoma skin cancer in psoriatic patients treated with high-dose phototherapy. J. Dermatol. Treat. 27, 443–447. https://doi.org/10.3109/09546634.2015.1133882 (2016).
Wang, E. B. et al. In Photobiomodulation in the Brain 49–66 (Elsevier, 2019).
Clement, M., Daniel, G. & Trelles, M. Optimising the design of a broad-band light source for the treatment of skin. J. Cosmet. Laser Ther. 7, 177–189 (2005).
Hu, D., van Zeyl, M., Valter, K. & Potas, J. R. Sex, but not skin tone affects penetration of red-light (660 nm) through sites susceptible to sports injury in lean live and cadaveric tissues. J. Biophoton. 12, e201900010 (2019).
Salehpour, F. et al. Penetration profiles of visible and near-infrared lasers and light-emitting diode light through the head tissues in animal and human species: A review of literature. Photobiomodul. Photomed. Laser Surg. 37, 581–595 (2019).
Chung, H. et al. The nuts and bolts of low-level laser (light) therapy. Ann. Biomed. Eng. 40, 516–533. https://doi.org/10.1007/s10439-011-0454-7 (2012).
Hamblin, M. R. Mechanisms and mitochondrial redox signaling in photobiomodulation. Photochem. Photobiol. 94, 199–212. https://doi.org/10.1111/php.12864 (2018).
Huang, Y.-Y., Sharma, S. K., Carroll, J. & Hamblin, M. R. Biphasic dose response in low level light therapy—An update. Dose-Response 9, 11-009 (Hamblin) (2011).
de Freitas, L. F. & Hamblin, M. R. Proposed mechanisms of photobiomodulation or low-level light therapy. IEEE J. Sel. Top. Quantum Electron. 22, https://doi.org/10.1109/jstqe.2016.2561201 (2016).
Lev-Tov, H., Mamalis, A., Brody, N., Siegel, D. & Jagdeo, J. Inhibition of fibroblast proliferation in vitro using red light-emitting diodes. Dermatol. Surg. 39, 1167–1170. https://doi.org/10.1111/dsu.12212 (2013).
Mamalis, A. & Jagdeo, J. Light-emitting diode-generated red light inhibits keloid fibroblast proliferation. Dermatol. Surg. 41, 35–39. https://doi.org/10.1097/01.Dss.0000452650.06765.51 (2015).
Mamalis, A. et al. High fluence light emitting diode-generated red light modulates characteristics associated with skin fibrosis. J. Biophoton. 9, 1167–1179. https://doi.org/10.1002/jbio.201600059 (2016).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. https://doi.org/10.1186/s13059-014-0550-8 (2014).
Lingappan, K. NF-κB in oxidative stress. Curr. Opin. Toxicol. 7, 81–86 (2018).
Hammouda, M. B., Ford, A. E., Liu, Y. & Zhang, J. Y. The JNK signaling pathway in inflammatory skin disorders and cancer. Cells 9, 857 (2020).
Brembilla, N. C. et al. Th17 cells favor inflammatory responses while inhibiting type I collagen deposition by dermal fibroblasts: differential effects in healthy and systemic sclerosis fibroblasts. Arthritis Res. Ther. 15, 1–13 (2013).
Wernig, G. et al. Unifying mechanism for different fibrotic diseases. Proc. Natl. Acad. Sci. 114, 4757–4762 (2017).
Murai, M. et al. Tryptophan photo-product FICZ upregulates AHR/MEK/ERK-mediated MMP1 expression: Implications in anti-fibrotic phototherapy. J. Dermatol. Sci. 91, 97–103. https://doi.org/10.1016/j.jdermsci.2018.04.010 (2018).
Aoki, M. et al. siRNA knockdown of tissue inhibitor of metalloproteinase-1 in keloid fibroblasts leads to degradation of collagen type I. J. Invest. Dermatol. 134, 818–826. https://doi.org/10.1038/jid.2013.396 (2014).
Iimuro, Y. et al. Delivery of matrix metalloproteinase-1 attenuates established liver fibrosis in the rat. Gastroenterology 124, 445–458. https://doi.org/10.1053/gast.2003.50063 (2003).
LeBleu, V. S. et al. Identification of human epididymis protein-4 as a fibroblast-derived mediator of fibrosis. Nat. Med. 19, 227–231. https://doi.org/10.1038/nm.2989 (2013).
Liebel, F., Kaur, S., Ruvolo, E., Kollias, N. & Southall, M. D. Irradiation of skin with visible light induces reactive oxygen species and matrix-degrading enzymes. J. Investig. Dermatol. 132, 1901–1907 (2012).
Buechner, N. et al. Changes of MMP-1 and collagen type Iα1 by UVA, UVB and IRA are differentially regulated by Trx-1. Exp. Gerontol. 43, 633–637 (2008).
Li, W. et al. RNA-Seq analysis revealed the molecular mechanisms of photobiomodulation effect on human fibroblasts. Photodermatol. Photoimmunol. Photomed. 36, 299–307 (2020).
Lee, R., Hwu, Y., Lu, C. & Li, S. P-139 Correlation of PRSS35 and SERPINE2 gene expression levels in cumulus cells with oocyte maturation and the potential as a biomarker to predict embryo quality. Hum. Reproduct. 28 (2013).
Yuda, E., Ogasawara, H., Yoshida, Y. & Hayano, J. Suppression of vagal cardiac modulation by blue light in healthy subjects. J. Physiol. Anthropol. 35, 24–24. https://doi.org/10.1186/s40101-016-0110-x (2016).
Biernacka, A., Dobaczewski, M. & Frangogiannis, N. G. TGF-β signaling in fibrosis. Growth Factors (Chur, Switzerland) 29, 196–202. https://doi.org/10.3109/08977194.2011.595714 (2011).
Mamalis, A. & Jagdeo, J. High-fluence light-emitting diode-generated red light modulates the transforming growth factor-beta pathway in human skin fibroblasts. Dermatol. Surg. 44, 1317–1322 (2018).
Flanders, K. C. Smad3 as a mediator of the fibrotic response. Int. J. Exp. Pathol. 85, 47–64. https://doi.org/10.1111/j.0959-9673.2004.00377.x (2004).
Ito, Y. et al. Overexpression of Smad2 reveals its concerted action with Smad4 in regulating TGF-beta-mediated epidermal homeostasis. Dev. Biol. 236, 181–194. https://doi.org/10.1006/dbio.2001.0332 (2001).
Lakos, G. et al. Targeted disruption of TGF-beta/Smad3 signaling modulates skin fibrosis in a mouse model of scleroderma. Am. J. Pathol. 165, 203–217. https://doi.org/10.1016/s0002-9440(10)63289-0 (2004).
Xue, M., Gong, S., Dai, J., Chen, G. & Hu, J. The treatment of fibrosis of joint synovium and frozen shoulder by Smad4 gene silencing in rats. PLoS ONE 11, e0158093 (2016).
Mori, Y., Chen, S. J. & Varga, J. Expression and regulation of intracellular SMAD signaling in scleroderma skin fibroblasts. Arthritis Rheum. 48, 1964–1978. https://doi.org/10.1002/art.11157 (2003).
Takagawa, S. et al. Sustained activation of fibroblast transforming growth factor-beta/Smad signaling in a murine model of scleroderma. J. Invest. Dermatol. 121, 41–50. https://doi.org/10.1046/j.1523-1747.2003.12308.x (2003).
Walton, K. L., Johnson, K. E. & Harrison, C. A. Targeting TGF-β mediated SMAD signaling for the prevention of fibrosis. Front. Pharmacol. 8, 461. https://doi.org/10.3389/fphar.2017.00461 (2017).
Alliston, T. et al. Repression of bone morphogenetic protein and activin-inducible transcription by Evi-1. J. Biol. Chem. 280, 24227–24237. https://doi.org/10.1074/jbc.M414305200 (2005).
Kurokawa, M. et al. The oncoprotein Evi-1 represses TGF-beta signalling by inhibiting Smad3. Nature 394, 92–96. https://doi.org/10.1038/27945 (1998).
Mamalis, A., Koo, E., Tepper, C. & Jagdeo, J. MicroRNA expression analysis of human skin fibroblasts treated with high-fluence light-emitting diode-red light. J. Biophoton. 12, e201800207. https://doi.org/10.1002/jbio.201800207 (2019).
Babalola, O., Mamalis, A., Lev-Tov, H. & Jagdeo, J. The role of microRNAs in skin fibrosis. Arch. Dermatol. Res. 305, 763–776. https://doi.org/10.1007/s00403-013-1410-1 (2013).
Gras, C. et al. miR-145 contributes to hypertrophic scarring of the skin by inducing myofibroblast activity. Mol. Med. (Cambridge, Mass.) 21, 296–304, https://doi.org/10.2119/molmed.2014.00172 (2015).
Serasanambati, M. & Chilakapati, S. R. Function of nuclear factor kappa B (NF-kB) in human diseases-a review. South Indian J. Biol. Sci. 2, 368–387 (2016).
Chen, A. C. et al. Low-level laser therapy activates NF-kB via generation of reactive oxygen species in mouse embryonic fibroblasts. PLoS ONE 6, e22453 (2011).
Kim, H. S. et al. Transcriptomic analysis of human dermal fibroblast cells reveals potential mechanisms underlying the protective effects of visible red light against damage from ultraviolet B light. J. Dermatol. Sci. 94, 276–283 (2019).
Do, N. & Eming, S. Skin fibrosis: Models and mechanisms. Curr. Res. Translat. Med. 64, 185–193 (2016).
Matei, A.-E. et al. Vascularised human skin equivalents as a novel in vitro model of skin fibrosis and platform for testing of antifibrotic drugs. Ann. Rheum. Dis. 78, 1686–1692 (2019).
Mignon, C., Uzunbajakava, N. E., Raafs, B., Botchkareva, N. V. & Tobin, D. J. Photobiomodulation of human dermal fibroblasts in vitro: decisive role of cell culture conditions and treatment protocols on experimental outcome. Sci. Rep. 7, 2797–2797. https://doi.org/10.1038/s41598-017-02802-0 (2017).
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics (Oxford, England) 29, 15–21. https://doi.org/10.1093/bioinformatics/bts635 (2013).
Chen, E. Y. et al. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 14, 128. https://doi.org/10.1186/1471-2105-14-128 (2013).
Kuleshov, M. V. et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90-97. https://doi.org/10.1093/nar/gkw377 (2016).
Feoktistova, M., Geserick, P. & Leverkus, M. Crystal violet assay for determining viability of cultured cells. Cold Spring Harbor Protoc. 2016, pdb. prot087379 (2016).
Acknowledgments
The work presented here was supported by the National Institute of General Medical Sciences of the NIH under Award No. K23GM117309 (Jagdeo), with additional support from the National Cancer Institute of the NIH under Award No. P30CA093373 (Maverakis/UC Davis Cancer Center), and UC Davis Department of Dermatology Funding (Jagdeo). The content is solely the responsibility of the authors and does not necessarily represent the views of the US Department of Veterans Affairs or the US government.
Funding
The work presented here was supported by the National Institute of General Medical Sciences of the NIH under Award No. K23GM117309 (Jagdeo), with additional support from the National Cancer Institute of the NIH under Award No. P30CA093373 (Maverakis/UC Davis Cancer Center), and UC Davis Department of Dermatology Funding (Jagdeo).
Author information
Authors and Affiliations
Contributions
E.A. designed studies, performed experiments, analyzed data, and co-wrote the manuscript. E.K. designed studies, performed experiments, analyzed data, and co-wrote the manuscript. A.Me. performed the statistical analysis and analyzed data. D.T. performed the statistical analysis and analyzed data. A.Mar. analyzed data. G.L. analyzed data. A.Mam. edited the manuscript. R.R.I. provided materials and equipment and edited the manuscript. A.Ma’a. provided materials and equipment and edited the manuscript. E.M. provided materials and equipment, obtained funding, and edited the manuscript. J.J. obtained funding, designed studies, and edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Austin, E., Koo, E., Merleev, A. et al. Transcriptome analysis of human dermal fibroblasts following red light phototherapy. Sci Rep 11, 7315 (2021). https://doi.org/10.1038/s41598-021-86623-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-86623-2
This article is cited by
-
High-throughput transcriptomics
Scientific Reports (2022)
-
Anti-inflammatory effect of green photobiomodulation in human adipose-derived mesenchymal stem cells
Lasers in Medical Science (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.