{'fraction_non_alphanumeric': 0.04106761565836299, 'fraction_numerical': 0.016298932384341637, 'mean_word_length': 4.995519522082356, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 6, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Ionizing radiation exposure of stem cell-derived chondrocytes affects their gene and microRNA expression profiles and cytokine production

Human induced pluripotent stem cells (hiPSCs) can be differentiated into chondrocyte-like cells. However, implantation of these cells is not without risk given that those transplanted cells may one day undergo ionizing radiation (IR) in patients who develop cancer. We aimed to evaluate the effect of IR on chondrocyte-like cells differentiated from hiPSCs by determining their gene and microRNA expression profile and proteomic analysis. Chondrocyte-like cells differentiated from hiPSCs were placed in a purpose-designed phantom to model laryngeal cancer and irradiated with 1, 2, or 3 Gy. High-throughput analyses were performed to determine the gene and microRNA expression profile based on microarrays. The composition of the medium was also analyzed. The following essential biological processes were activated in these hiPSC-derived chondrocytes after IR: "apoptotic process", "cellular response to DNA damage stimulus", and "regulation of programmed cell death". These findings show the microRNAs that are primarily responsible for controlling the genes of the biological processes described above. We also detected changes in the secretion level of specific cytokines. This study demonstrates that IR activates DNA damage response mechanisms in differentiated cells and that the level of activation is a function of the radiation dose.

Upper airway reconstruction, including tracheal and laryngotracheal resection and reconstruction, is a common surgical procedure after resection of tumours or stenotic inflammatory lesions 1 . However, restoring proper airway function in these cases is difficult 2 . In recent years, the use of exogenous autologous stem or progenitor cells, which include autologous adipose-derived stem cells, mesenchymal stromal cells, and induced pluripotent stem cells (iPSC), has been proposed to improve airway reconstruction 3 . While the individual use of chondrocytes or mesenchymal stem cells (MSC) is, at present, not considered ideal for cell-based cartilage repair in the head and neck (H&N) area 4 , the combination of these two cell sources holds great promise for cartilage tissue engineering, as it reduces the required number of chondrocytes and attenuates most of the disadvantages of the individual cell types. Mixed cell cultures of chondrocytes and MSCs have also been shown to improve chondrogenesis and to reduce hypertrophy and tissue mineralization 5 .
Human iPSCs (hiPSCs) possess the ability to differentiate into chondrocytes to produce a cartilaginous matrix, thus obviating the need to biopsy healthy cartilage 6 . The most significant benefit of hiPSC-derived progenitor cells for cartilage regeneration is that they may provide an unlimited supply of chondrogenic cells for in vivo applications. Consequently, cell-based therapy has been recognized as a potentially effective approach to healing tracheal cartilage 7 . Imaizumi et al. 8 suggested that hiPSCs could be a new cell source to regenerate tracheal cartilage. However, this strategy is not without risk given that those transplanted cells may eventually be subjected to ionizing radiation (IR) in patients who develop cancer 9 . At present, published data on irradiated

Results
Differentiated cells possess features of chondrocytes. We first formed EBs from hiPSCs in a culture medium deprived of FGF-2 for 7 days. Those EBs were differentiated towards chondrocytes in the presence of TGF-β3 for 21 days (Fig. 1A). These chondrocyte-like cells exhibited the characteristics of chondrocytes, including expression of desirable markers such as COMP, type II collagen, and aggrecan ( Fig. 1B and Supplementary Fig. 4). Glycosaminoglycan (GAG) deposition was marked by alcian blue staining. Blue staining and condensed cartilage-like nodules were evident in EBs formed during chondrogenic differentiation (Fig. 1B). After differentiation, the chondrocyte-like cells were irradiated according to conventional therapy protocols (1, 2, 3 Gy) (Fig. 2). After IR, no notable changes in cell morphology were observed (Fig. 1C); however, we clearly We assumed the following selection criteria for changed gene expression: an expression fold difference > abs. 1.5 and an adjusted p value ≤ 0.05 with 20% false discovery rate (FDR) correction. According to these criteria, the number of downregulated or upregulated genes in the cells 1 h after IR (1, 2 and 3 Gy, respectively) versus non-irradiated controls was as follows: downregulated: 10, 0, 1 genes; upregulated: 6, 0, 1 genes. At 24 h, the numbers of significantly down-or upregulated genes were as follows: downregulated: 28, 40, and 70, respectively and upregulated: 7, 40 and 24, respectively (Supplementary Fig. 1). Given that changes at 1 h post-IR were barely perceptible, we decided to analyze only the material from cells at 24 h after irradiation.
Due to the structure of the GO database, single genes can often be assigned to many ontological terms. For this reason, the relationship between genes and GO terms were mapped with circos plots, with visualization of logFC values and gene symbols. All of those genes were either upregulated (Fig. 4A) or downregulated (Fig. 4B) in the irradiated cells compared to controls. The most upregulated or downregulated genes (CCNE2, DLGAP5, KIF20B, FAS, TP53INP1, BTG2, DDB2, ZMAT3 and CENPE) were then subjected to further qPCR evaluation (Fig. 5). All these genes are involved in several different pathways activated after IR, including reprogrammed cell death, response to radiation, and mitotic cell cycle processes.
The expression of genes engaged in DDR in the irradiated hiPSC-derived chondrocytes was high when compared to non-irradiated cells. Genes involved in apoptosis were also upregulated in the irradiated chondrocytelike cells. By contrast, expression of genes involved in cell division was strongly downregulated in irradiated cells.
We compared expression of the most downregulated genes engaged in DDR in the hiPSC-derived chondrocytes (1 Gy vs. control; 2 Gy vs. control; 3 Gy vs. control) to those obtained from laryngeal tumor tissues compared to normal samples (data from the Cancer Genome Atlas [TCGA] were obtained through the FireBrowse server (http:// gdac. broad insti tute. org/) 13 (Fig. 6). After IR, the expression of genes engaged in the cell cycle and cell division (e.g., CENPE, KIF11, NUSAP1 and MCM6) was strongly downregulated in hiPSC-derived cells. By contrast, the gene expression profile of tumor samples were strongly upregulated compared to normal samples. www.nature.com/scientificreports/ In addition, we observed a radiation dose-dependent upregulation of some genes (e.g. SPATA18, DDB2 and EDA2R) involved in DDR in the differentiated cells. By contrast, the comparison between laryngeal tumor samples and normal samples did not find any upregulation of expression of genes involved in DDR. The differing gene expression profile between (irradiated vs. non-irradiated) chondrocyte-like cells and (tumor vs. normal) laryngeal samples may indicate some degree of safety (i.e., low tumorigenic potential) of using differentiated cells that may be subjected to IR. 3 Gy vs. controls) assigned to gene ontology (GO) terms where the red color represents GO terms whose genes are downregulated while green corresponds to GO terms of upregulated genes.  Fig. 2). The smallest change in microRNA expression was observed in cells that received the hyperfractionated dose of 1 Gy.

Analysis of the composition of the medium in chondrocyte-like cells differentiated from hiP-SCs after IR treatment.
One of the main aims of the present study was to assess the secretory activity of cells treated with IR and compare their activity level with non-irradiated cells. Thus, media from the cells were collected 24 h after IR. We found significant alterations in the following 38 pro-and anti-inflammatory cytokines: basic FGF; CTACK; Eotaxin; G-CSF; GM-CSF; GRO-alpha; HGF; IFN-alpha 2; IFN-gamma; IL-1 alpha; IL-1 beta; IL-2; IL-2R alpha; IL-4; IL-6; IL-7; IL-8; IL-9; IL-12; IL-13; IL-17; IP-10; LIF; MCP-1; MCP-3; M-CSF; MIG; MIP-1 alpha; MIP-1 beta; PDGF-BB; NGF-beta; RANTES; SCF; SCGF-beta; SDF-1 alpha; TNFalpha; TNF-beta; and TRAIL) in the medium of hiPSC-derived cells after IR. By contrast, we did not observe any notable changes in the level of IL-3 and MIF. These results are shown in Table 1 as mean values (pg/ml of secreted cytokine) with standard deviation. The mean values ranged from 0.223 pg/ml (IL-3 in controls) to 3448.34 pg/ ml (SCGF-beta in controls). In the graphs (Supplementary Fig. 3) we highlight the most statistically significant relationships between particular variants (controls vs. 1 Gy; controls vs. 2 Gy; controls vs. 3 Gy; 1 Gy vs. 2 Gy; 1 vs. 3 Gy and 2 Gy vs. 3 Gy). The statistical significance for each comparison is shown in Supplementary Table 1. Interestingly, in most cytokines, there was a notable decrease in secretion (e.g., G-CSF, GRO-alpha, HGF, IP-10, LIF, SCGF-beta, TNF-alpha) after 2 Gy irradiation. The cytokine secretion levels after 1 and 3 Gy were, to a large extent, more similar to that of non-irradiated cells than the 2 Gy dose (conventional therapy). These findings demonstrate that the specific cytokine secretion profile depends on the applied dose (and treatment protocol).

Discussion
Findings from many recent studies suggest that hiPSCs can differentiate into various types of functional cells such as cardiomyocytes and osteocytes 7 . However, transplantation of these cells into humans presents a major risk, as the data indicate that a small number of hiPSCs remain undifferentiated at transplantation and can spontaneously transform into rapidly proliferating tumors 14 .
Our group previously applied molecular methods to analyze DDR mechanisms in chondrocyte-like cells after IR 10 . In that study, we assessed, inter alia, the kinetics of DNA DSB formation, activation of DNA repair mechanisms, cell cycle changes, and apoptosis. Our findings showed that although hiPSC-derived chondrocyte cells readily formed DSBs and accumulated in the G2 phase, these cells were also characterized by highly efficient DNA repair mechanisms. Moreover, those cells appeared to be more likely to undergo senescence rather than apoptosis after IR exposure. However, in that study, we did not use conventional radiotherapy regimens/dose for  Fig. 7) in irradiated cells. Moreover, in the present study, we used a specially modified anthropomorphic phantom to ensure similar conditions to those present in the patient's larynx during treatment. This is an important advantage of the present study.
In this study, we focused on the high-throughput analysis of the response of chondrocyte-like cells differentiated from hiPSCs ( Fig. 1 and Supplementary Fig. 4) to IR treatment based on hypofractionated and hyperfractionated protocols (Fig. 2). Our findings suggest that those cells are characterized by specific DDR mechanisms. Despite the immense promise of hiPSCs, data on the response of those cells to IR are limited, which is why we decided to investigate the gene and microRNA expression profile in hiPSC-derived chondrocytes after IR (Supplementary Fig. 1 and Supplementary Fig. 2). As noted above, since the gene expression profile of irradiated cells after 1 h did not reveal any significant changes ( Supplementary Fig. 1), we only used the material collected 24 h after irradiation for all further analyses.
Below we discuss our findings in the context of the limited published data on the response of iPSC-derived cells and chondrocytes to IR.
Becker et al. 15 irradiated hiPSC-derived cardiomyocytes with 5 Gy of x-ray irradiation and then performed deep-sequencing of RNA to investigate gene expression changes in these cells. Those authors identified genes implicated in cardiac calcium homeostasis (PDE3B), oxidative stress response (FDXR and SPATA18), and the etiology of cardiomyopathy (SGCD, BBC3 and GDF15). The large fraction of downregulated genes was dominated by genes involved in cell cycle processes 15 . In our study, we also observed decreased expression of genes engaged in certain cell cycle processes in irradiated cells, including "cell cycle", "mitotic cell cycle process", "cell division" (Figs. 3, 4). Furthermore, we have demonstrated the close relationship between the most upregulated and downregulated genes (e.g., TP53INP1, DDB2, ZMAT3, SESN1) and the microRNAs that control their expression, such as hsa-miR-17, hsa-miR-5252, hsa-miR-19a-3p, and hsa-miR-551b-5p (Figs. 5, 7). www.nature.com/scientificreports/ Miyake et al. 16 successfully derived keratinocytes from iPSCs and characterized the differentiated states and DDR of these cells. The derived keratinocytes showed progenitor-like characteristics (e.g., high expression of integrin a6/CD71 and resistance to IR) as a result of DDR. After IR exposure, g-H2AX foci decreased more slowly in passages 2 and 3 keratinocytes compared to passage 1 keratinocytes, which suggests that matured cells attenuate DNA repair activity 17 .
Hong et al. showed that low-dose γ-radiation (LDR) 19 (0.5-1 cGy) suppresses IL-1β-induced dedifferentiation and inflammation of chondrocytes without any noticeable side effects. In that study, LDR attenuated signaling pathways both upstream and downstream of catenin, blocked dedifferentiation process throughout preventing IL-1β-induced Sox-9 suppression, and inhibited inflammation by counteracting IL-1β-induced activation of NF-κB. These anti-disease effects establish LDR as a potentially valuable therapeutic tool for patients with cytokine-mediated cartilage disorders 19,20 .
Lepleux et al. 18 found that the radiation-induced bystander response can be transferred from chondrosarcoma cells treated with IR to non-irradiated chondrocyte cells. A great impact on chondrocyte survival was observed at low doses. Besides, this effect was more notable after low LET (X-rays) in comparison with high LET irradiation (C-ions). Bystander factors secreted in the conditioned medium contributed to the reduced proliferation and elevated DNA damage level at low doses (X-rays and C-ions). Bystander biological activity was absent after dilution and heat treatment of the conditioned medium, leading the authors to suggest that factors such as TNF-α and IL-6 may have contributed to this effect 18 .
The findings of the present study, based on TCGA data, indicate that the downregulated expression of genes engaged in the cell cycle and cell division in hiPSC-derived chondrocytes differed from the gene expression profile of tumor and normal tissue samples from laryngeal tissues, in which those genes were strongly upregulated. Similarly, the differentiated cells presented a dose-dependent upregulation of some genes involved in DDR compared to laryngeal tumor cells (Fig. 6). This negative correlation suggests that DDR mechanisms in hiPSC-derived chondrocytes may be more similar to differentiated mature cells than to potentially tumorigenic hiPSCs or partially differentiated cells. We assume that irradiation of these cells (hiPSC-derived chondrocytes) is unlikely to induce malignant transformation; consequently, it should be safe to use such cells. However, more research is needed to elucidate this poorly understood area. www.nature.com/scientificreports/ In general, inflammatory cytokines, GFs, and proteases can affect cancer cell invasion, the bystander effect, and radiation-derived tissue complications including fibrosis and genomic instability. As a result, these factors can greatly impact intrinsic cellular radiosensitivity 21 . IR can provoke a multi-layered signaling response by activating many pro-survival pathways accompanied by transient activation of key transcription factors such as Nuclear Factor kappa B (NF-ĸB) as well as signal transducers and activators of transcription members (STATs) 22 . NF-ĸB plays a crucial role in immune and inflammatory responses due to its ability to regulate the expression of proinflammatory cytokines and chemokines such as tumor necrosis factor-alpha (TNF-α), interleukin-1 (IL-1), IL-2, IL-6, and monocyte chemoattractant protein-1 for inflammatory cells (MCP-1) 23 .  www.nature.com/scientificreports/ There is evidence 24,25 that exposure of cells and tissues to IR stimulates the expression of many different cytokines and GFs, including the following: TNF-α, interleukin-1 alpha (IL-1α), interleukin-1 beta (IL-1β), IL-6, type I interferon (type I IFN), granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-4, IL-5, IL-10, IL-12, IL-18, and TGF-β. The inflammatory reaction induced by IR is mediated by inflammation-related cytokine genes such as TNF-a, IL-1, IL-6, IL-8, interferon-gamma (IFN-γ), granulocyte-macrophage colony-stimulating factor (G-CSF), vascular endothelial growth factor (VEGF), and epidermal growth factor receptor. Additionally, downregulation of the existing inflammatory response depends largely on the half-life of the proinflammatory cytokines and the production of the antiinflammatory cytokines, such as IL-4, IL-10, IL-13, and TGF-β 24,25 . Similarly, Schaue et al. 26 found that in vitro and in vivo exposure of cells and tissues to IR treatment induces expression of many cytokines and GFs, including the following: TNF-a, IL-1aa, IL-1b, type I IFN, GM-CSF, IL-4, IL-5, IL-6, IL-10, IL-12, IL-18, VEGF, bFGF, and TGF-βa. Many of these appear as immediate early genes and thus qualify as "radiation-inducible". In tissues, radiation damage ultimately manifests as fibrosis. Macrophages accumulate in the damaged tissue and subsequently switch from being proinflammatory to producing fibrosispromoting cytokines (PDGF and TGF-β) 26 . The radioprotective effect of IL-12 in hematopoietic tissue depends on its interaction with IL-1 and SCF. Whereas, the sensitizing effect on gastrointestinal tissue is mediated by IFN-y and TNF 27 .

Mean (pg/ml) SD Mean (pg/ml) SD Mean (pg/ml) SD Mean (pg/ml) SD
Schrӧder et al. 28 found that IR causes nonlinear dose-dependent effects on the secretion of proinflammatory cytokines. Among the 23 inflammatory markers assessed in that study, four of them: KC, MCP-1, RANTES, and G-CSF were found to be significantly altered. The findings of that study indicate that LDR has an immune modifying capacity on the response of endothelial cells, both in a conventional 2D culture system and an ECMbased 3D model 28 .
In a mouse model, Shan et al. 29 described the effect of various doses (range, 0.05-6 Gy) of whole-body irradiation (WBI) on peritoneal macrophages. Those authors found a marked stimulation of IL-12 and IL-18 secretion correlated with the nuclear translocation of NF-ĸB. Moreover, upregulated expression of cytoplasmic MyD88 and the surface molecules CD14 and TLR4-MD2 was also proved. The authors also demonstrated stimulation of IL-12 and IL-18 secretion by mouse peritoneal macrophages after WBI. They explored the possible mechanisms and implications of these findings in cancer radiotherapy. Both low (0.075 Gy) and high (2 Gy) dose IR led to sustained stimulation of IL-12 and IL-18 secretion by mouse macrophages. This finding was paralleled by activation of NF-ĸB as well as upregulated expression of CD14 and TLR4-MD2 on the macrophage surface and MyD88 in the cytoplasm 29 .
Kiang et al., based on cytokine profiling of blood from B6D2F1 female mice irradiated at low dose rates, selected G-CSF and IL-18 for the extended experiment involving higher dose rates and large dose range. Those authors found that radiation dose rates do not influence the cytokine response to IR. Moreover, G-CSF responded to all of the tested radiation dose rates, which ranged from 0.04 Gy/min to 1.94 Gy/min. The increases were biphasic, occurring on day 2. In contrast to G-CSF, circulating IL-18 increased later and in a persisted manner 30 .
Our findings also show that cytokines play a key role in the DDR of irradiated cells. We found 38 cytokines that differed significantly after IR (Table 1 and Supplementary Fig. 3). We observed altered secretion of numerous relevant cytokines, including IL-12, IL-6, IL-8, HGF, MCP-3, TNF-beta, TRAILGRO-alpha, and IFN-gamma. Interestingly, after irradiating cells with 2 Gy of IR, the secretion level of most cytokines decreased, which may partly explain why 2 Gy is commonly used in conventional radiotherapy: it provides an effective, therapeutic fractional dose but does not induce an excessive immune response.

Conclusions
In the present study, we have demonstrated, in a purpose-designed phantom, the broad context of activated DDR mechanisms of irradiated hiPSC-derived chondrocytes. The activation of crucial signaling pathways following IR was notable, with both upregulation and downregulation of a wide range of genes involved in DDR. Our analysis revealed that expression of those genes is regulated by specific microRNAs. Based on proteomics, we identified the profile of 38 cytokines secreted into the medium after IR. We believe the findings reported here help explain how stem cell-derived cells respond to irradiation, findings that may have clinically significant implications. Our results show that changes in both the mRNA and microRNA expression profiles become increasingly evident as the fractional radiation dose increases (hyperfractionation < conventional therapy < hypofractionation). However, a fractional dose of 2 Gy-the conventional radiotherapy dose-does not cause excessive activation of immune response components. Finally, based on TCGA data, it appears that irradiation of chondrocyte-like cells is unlikely to lead to malignant transformation, which suggests that they can be safely used for transplantation purposes.

Materials and methods
Cell culture. The commercially available hiPSC cell line ND41658*H (Coriell Cell Repository, NY, USA) was cultured as described elsewhere 31 .
HiPSCs were used to form EBs. After 7 days, the EBs were differentiated into a chondrogenic lineage in a defined medium supplemented with transforming growth factor-beta (TGF-β) (10 ng/ml) according to a previously established protocol 32 (Fig. 1).

Irradiation and dosimetry.
To irradiate the chondrocyte-like cells, the specially modified and redesigned Alderson Radiation Therapy phantom (ART) was used to obtain conditions similar to those present in the larynx of a patient with laryngeal cancer. It is crucial to establish proper attenuation of radiation (total mass attenuation coefficients similar to those present in the patient's tissue, adequate energy spectrum, absorbed and scatter radiation to the cell layer). The ART phantom was immobilized with the QFixEncompassTM system (Avondale, PA, www.nature.com/scientificreports/ USA) to ensure positioning reproducibility during preparation and irradiation procedure. Tomographic scans of the phantom were performed using Siemens SOMATOM Definition AS (120 kV X-Ray Tube Voltage, 1 mm slice thickness, head, and neck protocol). An Eppendorf tube with 1 × 10 6 chondrocyte-like cells was used to irradiate the cells. The quasi-gross tumor volume (GTV) was contoured. Geometrical margins of 1 cm around the GTV were applied to create the quasiclinical target volume (CTV). Following typical procedures for H&N cancer, a 0.5 cm margin around the CTV was added to ensure that the quasi-PTV received the prescribed doses. In the phantom, we also outlined quasi-OARs representing the spinal cord, parotid glands, and mandible. Treatment plans were prepared using the volumetric modulated arc therapy (VMAT) technique with two co-planar arcs (Fig. 2). The goal of the plans was to cover at least 95% of the PTV with the planned prescription dose (1, 2, or 3 Gy). Optimization and calculations were done on the Eclipse Planning System, v. 15.6 (Varian Medical Systems, USA) using the anisotropic analytical algorithm (AAA).
Irradiation was performed with a 6 MV photon beam produced by the Unique linac (Varian Medical Systems, USA) with a 120 leaf multileaf collimator and a maximal nominal dose rate of 600 MU/min. Phantom positioning was done using the X-ray 6D IGRT (orthogonal kV-kV) imaging system with the BrainLabExacTrac (BrainLab AG, Feldkirchen, Germany). Before irradiation, the electronic portal imaging device (EPID) system was used to perform dosimetric pre-verification. For the gamma evaluation criteria 3%/3 mm, the following results were obtained for the plans (%):  (Fig. 2).
Immunofluorescence analysis and alcian blue staining. HiPSC-derived chondrocytes were stained with the use of following antibodies: anti-type II collagen, aggrecan and anti-cartilage oligomerix matrix protein (COMP) according to a previously established protocol 32 .
Alcian blue staining was performed as previously described 33 to assess the deposition of proteoglycans in the differentiated cells.
Microarray study. Total RNA and miRNA-enriched fraction from hiPSCs after chondrogenic differentiation was purified with the miRNeasy Kit and RNeasy MinElute according to the manufacturer's instructions (217004, Qiagen, Hilden, Germany). The microarray study was conducted according to procedures described elsewhere [34][35][36] . The obtained CEL files were further analyzed using the R statistical language and Bioconductor package, including the relevant Bioconductor libraries. For the normalization, background correction, and calculation of the expression values of the examined genes or microRNAs, the Robust Multiarray Average (RMA) normalization algorithm implemented in the "Affy" library was applied. A complete gene/microRNA data table involving normalized gene/microRNA expression values, gene/microRNA symbols, gene/microRNA names and Entrez IDs was generated on the basis of assigned biological annotations taken from the "pd.hugene.2.1.st" and "pd.mirna3.1" libraries. For the expression and statistical assessment, the linear models for microarray data included in the "limma" library was applied. The established cut-off criteria were based on both differences in expression fold change (FC) greater than abs. 1.5 and an adjusted p value of ≤ 0.05. Genes/microRNAs fulfilling the aforementioned selection criteria were considered to be significantly different and therefore were subjected to the further analyses. www.nature.com/scientificreports/ lyze cytokines, chemokines, and growth factors (GF) potentially involved in the inflammatory process and in the initiation and progression of cancer. The concentrations (pg/ml) of these molecules were compared to the standard curve of serial dilutions of standard molecules. The analysis was performed according to the manufacturer's instructions. Each sample was analyzed in duplicate.

RTCA cell proliferation assay.
To verify the effect of IR on the proliferation rate of chondrocyte-like cells differentiated from hiPSCs, we applied an electrical impedance-based cell proliferation assay, named Real-Time Cell Analyser (RTCA, Roche Applied Science, GmbH, Penzberg, Germany) 38 . The RTCA system detects fluctuations in electrical impedance on the integrated sensory electrodes located at the bottom of the chamber's 16-hole slide plates (E-Plate 16), which are covered by dividing cells. Electrical impedance is measured at 15-min intervals throughout the cultivation period. The main RTCA readout is the "Cell Index"-a measurable parameter corresponding to the relative change in electrical impedance depending on the rate of proliferation or apoptosis of the cultivated cells. Each experimental group was seeded in the eight E-plate wells to a final volume of 200 μl per well. The cell index was normalized (normalized cell index) at time of examined substances administration using software from the same manufacturer (RTCA Software, v. 1.2, November 2009) 38 . Irradiated cells were cultivated per approximately 48 h. Each experiment was repeated ≥ three times. We decided to present results involving proliferation completed in 39 h because after this time, the control (0 Gy) reached a plateu.