Research Article | Open Access
Dianbao Zhang, Jing Wang, Zhe Wang, Tao Zhang, Ping Shi, Xiliang Wang, Feng Zhao, Xiaoyu Liu, Xuewen Lin, Xining Pang, "miR-136 Modulates TGF-β1-Induced Proliferation Arrest by Targeting PPP2R2A in Keratinocytes", BioMed Research International, vol. 2015, Article ID 453518, 8 pages, 2015. https://doi.org/10.1155/2015/453518
miR-136 Modulates TGF-β1-Induced Proliferation Arrest by Targeting PPP2R2A in Keratinocytes
Keratinocytes proliferation is critical for the capacity to heal wounds and accumulating evidences have proved that dysregulation of microRNAs is involved in proliferation of keratinocytes. However, the molecular mechanisms remain to be completely elucidated. Here, we show that miR-136 was significantly decreased by TGF-β1 treatment in HaCaT cells and normal human epidermal keratinocytes (NHEK), and it was a Smad3-dependent manner. By cell proliferation assay and cell cycle analysis, we found that reintroduction of miR-136 by transfection, as well as PPP2R2A silencing, counteracted TGF-β-induced proliferation arrest in HaCaT cells. Further, PPP2R2A was verified as a direct target of miR-136 by dual-luciferase reporter assays and Western blotting. These data suggest that miR-136 may play an important role during TGF-β1-induced proliferation arrest by targeting PPP2R2A in keratinocytes, which might represent a potential target for improving skin wound healing.
Keratinocytes proliferation is essential for the formation of skin appendages and reepithelialization after wounding , and epidermal homeostasis depends on a balance of keratinocytes between proliferation in the basal layer and differentiation/apoptosis in the superbasal layer [2–4]. So far, these functions have been shown to be stimulated by many growth factors, including transforming growth factor-β1 (TGF-β1) [5, 6]. Moreover, TGF-β1 secreted by fibroblasts and keratinocytes would cause proliferation arrest of keratinocytes through activating the signaling pathways by binding a receptor complex comprising two transmembrane serine/threonine kinases known as the type I and type II receptors . Thus, exploring the regulation mechanisms of TGF-β1-induced keratinocytes proliferation arrest is necessary to understand skin development and wound repair and to design novel therapeutic strategies for skin wound healing.
MicroRNAs (miRNAs) represent a large family of conserved noncoding small RNAs (~21 nucleotides long), which have emerged as key regulators of gene expression at both the transcriptional and translational levels. Recently, increasing evidences suggest that miRNAs are responsible for regulating proliferation and self-renewal of keratinocytes . Depletion of the miRNA processing enzyme Dicer in the epidermis results in hyperproliferation , and several miRNAs have been shown to modulate keratinocytes proliferation, such as miR-34 family , miR-125b , miR-200 family , miR-203 , and miR-210 . Nonetheless, the crucial role of miR-136 in the TGF-β1-induced proliferation arrest is largely unknown.
PPP2R2A, also known as B55α, is a regulatory subunit of the protein phosphatase 2A (PP2A), which is one of the four major serine/threonine (Ser/Thr) phosphatases, and is implicated in the negative control of cell growth and division in mammalian cells [15, 16]. The crystal structure of PPP2R2A has provided insight on how the B subunit forms a pocket with the catalytic subunit that may confer substrate specificity . PPP2R2A plays an important role in the modulation of the phosphorylation state of pocket proteins and mitotic CDK substrates throughout the cell cycle and in quiescent cells .
In the present study, we identified and characterized PPP2R2A as a functional downstream target of miR-136 and examined the roles of miR-136 and PPP2R2A in TGF-β1-induced proliferation arrest in keratinocytes.
2. Materials and Methods
2.1. Cell Isolation and Culture
Normal human epidermal keratinocytes (NHEK) were isolated from foreskins (18–26 years old, ) using Dispase II (Roche, IN, USA) to remove the dermis from the epidermis followed by trypsinization and culture in EpiLife serum-free medium with 60 μM calcium supplemented with HKGS (Gibco, NY, USA). Written informed consent was obtained from all patients and approval of the China Medical University Ethics Committee was obtained.
The immortalized human keratinocyte cell line HaCaT was kindly provided by Chundi He. Cos-7 cells were purchased from ATCC. HaCaT cells and Cos-7 cells were cultured in DMEM (Hyclone, MA, USA) supplemented with 10% fetal bovine serum. All the cells were maintained under standard conditions at 37°C and 5% CO2 in a humid atmosphere. For TGF-β1 treatment, NHEK or HaCaT cells were plated at 3,000 cells/cm2 and maintained for 24 h, after which they were serum-starved for another 24 h. Then, the cells were cultured in the absence or presence of TGF-β1 (R&D Systems, MN, USA) for 48 h and 72 h.
2.2. Quantitative Real-Time PCR (qRT-PCR)
Total RNA was extracted using TRIzol reagents (Invitrogen, CA, USA) according to the manufacturer’s protocol. To analyze miR-136 expression, reverse transcription PCR was performed using specific stem-loop RT primers from a Hairpin-it miRNAs qPCR Quantitation Kit (GenePharma, Shanghai, China), and quantitative real-time PCR was performed using the same kit on Applied Biosystems 7500 system. U6 was used as an internal control. To quantify mRNA level of PPP2R2A and Smad3, reverse transcription PCR was applied using PrimeScript RT Reagent Kit with gDNA Eraser (Takara, Dalian, China) and quantitative real-time PCR was performed using SYBR Premix Ex Taq (Takara, Dalian, China). GAPDH was used as an internal control. Fold changes of both miRNA and mRNA were calculated using the method. Each experiment was performed in triplicate. The primer sequences were as follows: 5′-AACCCACTTCCTGCTTAGTTGAG-3′ (forward) and 5′-GAACACCACAGTGATGAATCCAC-3′ (reverse) for human PPP2R2A; 5′-ACTACATCGGAGGGGAGGTC-3′ (forward) and 5′-TGCATCCTGGTGGGATCTTG-3′ (reverse) for human Smad3; 5′-GCACCGTCAAGGCTGAGAAC-3′ (forward) and 5′-TGGTGAAGACGCCAGTGGA-3′ (reverse) for human GAPDH.
2.3. miRNA and siRNA Transfection
Cells cultured in 6-well plates were transfected at 30–50% confluence using Lipofectamine RNAiMAX Transfection Reagent (Invitrogen, CA, USA) with Opti-MEM I (Gibco, NY, USA) according to the manufacturer’s instructions. Briefly, 25 pmol of miR-136 mimics, PPP2R2A siRNA (Santa Cruz, CA, USA), Smad3 siRNA, or negative control (NC) was used with 7.5 μL of the transfection reagent and the medium was replaced with fresh medium 24 h after transfection. Smad3 siRNA, miR-136 mimics, and NC were chemical synthesized by GenePharma (Shanghai, China) and the sequences were as follows: 5′-GAGCCTGGTCAAGAAACTCAATT-3′ for Smad3 siRNA, 5′-ACUCCAUUUGUUUUGAUGAUGGA-3′ for miR-136 mimics, and 5′-UUCUCCGAACGUGUCACGUTT-3′ for NC.
2.4. Cell Proliferation Assay
Cell proliferation was measured using Cell Counting Kit-8 (CCK-8, Dojindo, Kumamoto, Japan). 10 μL CCK-8 was added to cells maintained in 96-well plates and the cells were subsequently incubated for 2 h at 37°C. The absorbance of each sample at 450 nm was measured with a reference wavelength of 630 nm. All experiments were performed in triplicate and repeated three times.
2.5. Cell Cycle Analysis
Cell Cycle and Apoptosis Analysis Kit (Beyotime, Nantong, China) was used for cell cycle analysis. Cells were fixed in 70% ethanol in PBS at −20°C for 12 h after washing three times with cold PBS; then the cells were washed again and stained with 0.5 mL of propidium iodide (PI) staining buffer (it contains 200 mg/mL RNase A and 50 μg/mL PI) at 37°C for 30 min in the dark. Analyses were performed on BD LSR flow cytometer. The experiments were repeated three times.
2.6. Dual-Luciferase Reporter Assays
Bioinformatic analysis of microRNA target sites was performed using TargetScan Human Release 6.2 (http://www.targetscan.org/). Oligonucleotide pairs were designed to contain wild-type or mutant putative miRNA target sequences in 3′-UTR of PPP2R2A, with proper overhangs and an internal KpnI site for clone confirmation. Annealing and ligating these pairs into the pmirGLO Dual-Luciferase miRNA Target Expression Vector maintained the miRNA target region in the correct 5′ to 3′ orientation. Oligonucleotide pairs were annealed (95°C for 3 min, cooled to 37°C over 60 min and kept at 4°C for 30 min), ligated with the SacI-XbaI digested and gel-purified pmirGLO vector, and transformed into E. coli DH5α. The positive clones were identified by digesting miniprep-purified DNA with KpnI and verified by sequencing. Cos-7 cells plated on 24-well plates were cotransfected with 50 nmol of either miR-136 mimics or NC oligos and with 200 ng reporter plasmid containing either wild-type or mutant target site using Lipofectamine 2000 Transfection Reagent according to the manufacturer’s instructions. Cells were lysed and the relative luciferase activities were quantified at 48 h after transfection using Dual Glo Luciferase Assay System (Promega, WI, USA). All transfections were carried out in triplicate. The oligonucleotide sequences used for each reporter plasmid construction are listed in Table 1.
2.7. Western Blot
Protein was harvested using RIPA lysis buffer (Dingguo, Beijing, China) supplemented with 0.1 mmol PMSF, diluted with 5 × RSB, and denatured at 95°C for 10 min. Protein lysates were separated by 10% SDS-PAGE and transferred to PVDF membranes; the membranes were blocked in 5% nonfat milk for 1 h and then incubated with mouse anti-PPP2R2A (1 : 1000, CST, MA, USA) and mouse anti-GAPDH (1 : 5000, KangChen, Shanghai, China) antibodies overnight at 4°C. After incubation with HRP-conjugated secondary antibody (1 : 5000, KangChen, Shanghai, China), protein bands were visualized using Amersham ECL Prime Western Blotting Detection Reagent (GE Healthcare, NJ, USA) on Fujifilm LAS-3000 Imager and quantified by Image J.
2.8. Statistical Analysis
Data were expressed as the mean ± SD from at least three independent experiments. Student’s t-test and one-way ANOVA were used to compare the differences between groups. All statistical analyses were performed using SPSS 16.0 software (SPSS Inc., USA). was considered statistically significant.
3.1. miR-136 Was Downregulated by TGF-β1 in HaCaT and NHEK Cells
To investigate miR-136 expression changes during proliferation regulation of TGF-β1, we performed quantitative real-time PCR assay to measure miR-136 in keratinocytes treated with different concentration of TGF-β1. The result showed that the expression level of miR-136 in HaCaT decreased to 0.19-fold at 48 h after TGF-β1 treatment at 2 ng/mL and 0.13-fold for 5 ng/mL TGF-β1 (Figure 1(a)). To further confirm this result, we isolated and cultured primary epidermal keratinocytes from normal human skin. The results showed that miR-136 was also significantly downregulated by TGF-β1 treatment in NHEK (Figure 1(b)).
3.2. miR-136 Suppression by TGF-β1 Was Smad3-Dependent
It was reported that Smad2/3 signaling as well as downstream gene targets such as PAI-1 expression is required for TGF-β mediated cytostasis in HaCaT cells [19, 20]. To investigate whether canonical Smad2/3 signaling is required for miR-136 repression during TGF-β1-induced proliferation arrest in keratinocytes, HaCaT cells were transfected with Smad3 siRNA for 24 h and then treated with 2 ng/mL TGF-β1 for 48 h. By quantitative real-time PCR, the mRNA level of Smad3 was decreased to 0.2-fold at 48 h after transfection (Figure 1(c)), and in Smad3-silenced HaCaT cells, TGF-β1 treatment showed little effect on miR-136 expression (Figure 1(d)). These results indicated that miR-136 was suppressed by TGF-β1 in a Smad3-dependent manner in HaCaT cells.
3.3. Overexpression of miR-136 Overcomes TGF-β1-Induced Proliferation Arrest
Previous studies showed that TGF-β1 inhibited cell proliferation in keratinocytes ; to investigate whether forced expression of miR-136 is able to modulate TGF-β1-induced proliferation arrest of keratinocytes, HaCaT cells were transfected with miR-136 mimics or miR-NC and subsequently stimulated with 2 ng/mL TGF-β1 for 72 h. Successful increased expression of miR-136 was confirmed by qRT-PCR (Figure 2(a)). The results of cell proliferation assay showed that TGF-β1 treatment inhibited HaCaT cells proliferation effectively and led to >20% proliferation suppression (Figure 2(b)). Moreover, transfection of miR-136 mimics could negatively regulate the proliferation inhibition of TGF-β1 and caused 29% proliferation enhancement. As shown in Figure 2(c), the cell subpopulation in S phase was obviously decreased in TGF-β1 treated cells compared with the control group, while the cell number in S phase was significantly increased in cells transfected with miR-136 mimics versus miR-NC group. Taken together, these findings indicated that miR-136 might act as a modulator of TGF-β1-induced proliferation arrest in keratinocytes.
3.4. PPP2R2A Was a Direct Target of miR-136
To elucidate the underlying mechanisms by which miR-136 exerts its function, we explored miR-136 targets using the TargetScan bioinformatics algorithm. Our analysis revealed that PPP2R2A was a potential target of miR-136 based on putative conserved target sequences at positions 149–155, 712–719, and 1471–1478 of the PPP2R2A 3′-UTR (Figure 3(a)). To further examine whether miR-136 directly targets PPP2R2A, the luciferase reporters containing wild-type or mutant predicted miR-136 binding sites were cotransfected with miR-136 mimics or NC into Cos-7 cells. Luciferase assays were applied 48 h after transfection and the results showed that, compared to NC, transfection with miR-136 resulted in a significant decrease in renilla/firefly luciferase activity of wild-type site 1 and site 2 reporter (Figures 3(b) and 3(c)), while there was no significant decrease of wild-type site 3 reporter (Figure 3(d)). These results suggested that miR-136 repressed PPP2R2A through 2 specific 3′-UTR binding sites at positions 149–155 and 712–719. Notably, the expression of PPP2R2A in HaCaT cells substantially decreased at 48 h after miR-136 transfection (Figure 3(e)). Taken together, these results indicated that miR-136 negatively regulated PPP2R2A in a posttranscriptional manner in HaCaT cells.
3.5. PPP2R2A Was Involved in TGF-β1-Induced Proliferation Arrest
To further establish whether the counteraction of miR-136 overexpression against TGF-β1-induced proliferation arrest is mediated by repression of PPP2R2A, the expression of PPP2R2A in HaCaT cells treated with 2 ng/mL TGF-β1 for 48 h was analyzed by Western blot and the proliferation of PPP2R2A knockdown HaCaT cells stimulated with TGF-β1 was assessed by proliferation assay and flow cytometry. As shown in Figures 4(a) and 4(b), PPP2R2A in HaCaT cells was upregulated to about 2-fold by 2 ng/mL TGF-β1 treatment for 48 hours, and in Figures 4(c) and 4(d), knockdown of PPP2R2A clearly counteracted the proliferation arrest induced by TGF-β1. In addition, the reduction of cell cycle progression at G1/S transition induced by TGF-β1 was abrogated in PPP2R2A knockdown HaCaT cells. The effects of PPP2R2A knockdown were similar to those induced by miR-136 overexpression. Taken together, these findings indicated that PPP2R2A was a functionally important target of miR-136 and was involved in the TGF-β1 regulated proliferation of HaCaT cells.
Wound healing is a complex biological process, during which keratinocyte proliferation and migration are crucial steps for the rapid closure of the epidermis. TGF-β1 causes wound margin contraction at the early stage of wound healing and it is responsible for scar formation . These processes are controlled by a network of biomolecules in a spatiotemporal manner. Several recent reports have indicated that miRNAs are involved in regulating keratinocyte proliferation during wound healing [10–14]. Here, we aimed to clarify the biological role of miR-136 in keratinocytes proliferation regulation by TGF-β1. The experiments showed significant reduction of miR-136 in keratinocytes treated with TGF-β1 and canonical Smad2/3 signaling pathway was involved. Reintroduction of miR-136 by transient transfection, as well as silencing by siRNA of target PPP2R2A, blocked TGF-β1-induced proliferation arrest and increased the percentage of keratinocytes in the S phase of the cell cycle, while reducing the percentage of those in the G0/G1 phase. Our results supported the notion that TGF-β1-induced proliferation arrest was partially mediated by miR-136 reduction in HaCaT cells.
There are several reports that miR-136 was implicated in cell proliferation and played different roles in different types of cells. miR-136 is proposed to be a tumor suppressor in glioma and is capable of targeting the antiapoptosis genes AEG-1 and BCL-2 . However, miR-136 was found to target tumor suppressor PTEN in breast cancer cells . Recently, results of others indicated that miR-136 enhanced phosphorylation of Erk1/2 through inhibition of PPP2R2A expression to promoted cell proliferation in human non-small cell lung cancer, and the sequence at position 149–155 of the PPP2R2A 3′-UTR was determined to be the target site of miR-136 . In this study, we clarified miR-136 suppressed PPP2R2A expression by directly targeting two conserved regions at positions 149–155 and 712–719 in its 3′-UTR. Furthermore, the protein level of PPP2R2A was significantly reduced by miR-136 mimics transfection in HaCaT cells. Therefore, we proposed that, in the presence of TGF-β1, miR-136 was a positive modulator of keratinocytes proliferation and PPP2R2A was a functional target of miR-136 via binding to its 3′-UTR region.
In conclusion, our results indicated that the miR-136 treatment bypassed TGF-β1-induced proliferation arrest, at least partially through downregulation of PPP2R2A expression (Figure 5). Our data suggested that modulation of keratinocytes proliferation by miR-136 depends on the cellular environment, especially in the presence of TGF-β1, which likely has major effects on the expression of miR-136 target gene PPP2R2A. Further studies are needed to determine whether miR-136 based therapeutic strategies could be exploited to protect keratinocytes proliferation and function in the condition of enhanced TGF-β1 during wound healing.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Dianbao Zhang and Jing Wang contributed equally to this work.
This work was funded by the National Basic Research Program of China (2012CB518103), National Natural Science Foundation of China (81370883), Doctoral Fund of Ministry of Education of China (20132104110020), Science and Technology Planning Project of Liaoning Province (2012225080), and Science and Technology Planning Project of Shenyang (F11-262-9-01). The authors thank Dr. Hongchuan Xin and Dr. Ziwei Miao for critical reading and editing of this paper.
- I. Haase, R. Evans, R. Pofahl, and F. M. Watt, “Regulation of keratinocyte shape, migration and wound epithelialization by IGF-1 and EGF-dependent signalling pathways,” Journal of Cell Science, vol. 116, no. 15, pp. 3227–3238, 2003.
- E. Fuchs, “Skin stem cells: rising to the surface,” Journal of Cell Biology, vol. 180, no. 2, pp. 273–284, 2008.
- C. Pincelli and A. Marconi, “Keratinocyte stem cells: friends and foes,” Journal of Cellular Physiology, vol. 225, no. 2, pp. 310–315, 2010.
- J. Nie, X. Fu, and W. Han, “Microenvironment-dependent homeostasis and differentiation of epidermal basal undifferentiated keratinocytes and their clinical applications in skin repair,” Journal of the European Academy of Dermatology and Venereology, vol. 27, no. 5, pp. 531–535, 2013.
- S. Guo and L. A. Dipietro, “Factors affecting wound healing,” Journal of Dental Research, vol. 89, no. 3, pp. 219–229, 2010.
- G. S. Schultz and A. Wysocki, “Interactions between extracellular matrix and growth factors in wound healing,” Wound Repair and Regeneration, vol. 17, no. 2, pp. 153–162, 2009.
- J. S. Kang, C. Liu, and R. Derynck, “New regulatory mechanisms of TGF-β receptor function,” Trends in Cell Biology, vol. 19, no. 8, pp. 385–394, 2009.
- R. Yi and E. Fuchs, “MicroRNA-mediated control in the skin,” Cell Death & Differentiation, vol. 17, no. 2, pp. 229–235, 2010.
- M. Teta, Y. S. Choi, T. Okegbe et al., “Inducible deletion of epidermal Dicer and Drosha reveals multiple functions for miRNAs in postnatal skin,” Development, vol. 139, no. 8, pp. 1405–1416, 2012.
- D. Antonini, M. T. Russo, L. de Rosa, M. Gorrese, L. del Vecchio, and C. Missero, “Transcriptional repression of miR-34 family contributes to p63-mediated cell cycle progression in epidermal cells,” Journal of Investigative Dermatology, vol. 130, no. 5, pp. 1249–1257, 2010.
- N. Xu, P. Brodin, T. Wei et al., “MiR-125b, a MicroRNA downregulated in psoriasis, modulates keratinocyte proliferation by targeting FGFR2,” Journal of Investigative Dermatology, vol. 131, no. 7, pp. 1521–1529, 2011.
- Z. Lin, X. Wang, C. Fewell, J. Cameron, Q. Yin, and E. K. Flemington, “Differential expression of the miR-200 family microRNAs in epithelial and B cells and regulation of Epstein-Barr virus reactivation by the miR-200 family member miR-429,” Journal of Virology, vol. 84, no. 15, pp. 7892–7897, 2010.
- A. M. Lena, R. Shalom-Feuerstein, P. R. di Val Cervo et al., “miR-203 represses “stemness” by repressing ΔNp63,” Cell Death & Differentiation, vol. 15, no. 7, pp. 1187–1195, 2008.
- S. Biswas, S. Roy, J. Banerjee et al., “Hypoxia inducible microRNA 210 attenuates keratinocyte proliferation and impairs closure in a murine model of ischemic wounds,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 15, pp. 6976–6981, 2010.
- M. H. Schmitz, M. Held, V. Janssens et al., “Live-cell imaging RNAi screen identifies PP2A-B55α and importin-β 21 as key mitotic exit regulators in human cells,” Nature Cell Biology, vol. 12, no. 9, pp. 886–893, 2010.
- J. Batut, B. Schmierer, J. Cao, L. A. Raftery, C. S. Hill, and M. Howell, “Two highly related regulatory subunits of PP2A exert opposite effects on TGF-β/Activin/Nodal signalling,” Development, vol. 135, no. 17, pp. 2927–2937, 2008.
- A. Kurimchak and X. Graña, “PP2A counterbalances phosphorylation of pRB and mitotic proteins by multiple CDKs: potential implications for PP2A disruption in cancer,” Genes & Cancer, vol. 3, no. 11-12, pp. 739–748, 2012.
- G. Jayadeva, A. Kurimchak, J. Garriga et al., “B55α PP2A holoenzymes modulate the phosphorylation status of the retinoblastoma-related protein p107 and its activation,” The Journal of Biological Chemistry, vol. 285, no. 39, pp. 29863–29873, 2010.
- J. M. Overstreet, R. Samarakoon, K. K. Meldrum, and P. J. Higgins, “Redox control of p53 in the transcriptional regulation of TGF-β1 target genes through SMAD cooperativity,” Cellular Signalling, vol. 26, no. 7, pp. 1427–1436, 2014.
- M. Motizuki, K. Isogaya, K. Miyake et al., “Oligodendrocyte transcription factor 1 (Olig1) is a smad cofactor involved in cell motility induced by transforming growth factor-β,” The Journal of Biological Chemistry, vol. 288, no. 26, pp. 18911–18922, 2013.
- K. Räsänen and A. Vaheri, “TGF-beta1 causes epithelial-mesenchymal transition in HaCaT derivatives, but induces expression of COX-2 and migration only in benign, not in malignant keratinocytes,” Journal of Dermatological Science, vol. 58, no. 2, pp. 97–104, 2010.
- G. Viticchiè, A. M. Lena, F. Cianfarani et al., “MicroRNA-203 contributes to skin re-epithelialization,” Cell Death & Disease, vol. 3, no. 11, article e435, 2012.
- Y. Yang, J. Wu, H. Guan et al., “MiR-136 promotes apoptosis of glioma cells by targeting AEG-1 and Bcl-2,” FEBS Letters, vol. 586, no. 20, pp. 3608–3612, 2012.
- D. Y. Lee, Z. Jeyapalan, L. Fang et al., “Expression of versican 3′-untranslated region modulates endogenous microrna functions,” PLoS ONE, vol. 5, no. 10, Article ID e13599, 2010.
- S. Shen, H. Yue, Y. Li et al., “Upregulation of miR-136 in human non-small cell lung cancer cells promotes Erk1/2 activation by targeting PPP2R2A,” Tumor Biology, vol. 35, no. 1, pp. 631–640, 2014.
Copyright © 2015 Dianbao Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.