- 1Clinical College of Ophthalmology, Tianjin Medical University, Tianjin, China
- 2Tianjin Key Lab of Ophthalmology and Visual Science, Tianjin Eye Hospital, Tianjin Eye Institute, Nankai University Affiliated Eye Hospital, Tianjin, China
- 3The Second Affiliated Hospital of Anhui Medical University, Hefei, China
Background: To investigate the correlation between corneal biomechanical characteristics (in vitro and in vivo) and corneal densitometry in myopia.
Methods: The Pentacam (Oculus, Wetzlar, Germany) corneal densitometry (CD) and Corvis ST (Oculus, Wetzlar, Germany) exams were conducted prior to surgery for myopic patients who were intended to undergo small-incision lenticule extraction (SMILE). CD values (grayscale units, GSUs), and in vivo biomechanical parameters were obtained. The stromal lenticule was subjected to a uniaxial tensile test to obtain the elastic modulus E in vitro. We exam the correlations among in vivo, in vitro biomechanical characteristics and CD values.
Results: In this study, 37 myopic patients (63 eyes) were included. The mean age of participants was 25.14 ± 6.74 years (range:16–39 years). The mean CD values of the total cornea, anterior layer, intermediate layer, posterior layer, 0–2 mm region and 2–6 mm region were 15.03 ± 1.23 GSU, 20.35 ± 1.98 GSU, 11.76 ± 1.01 GSU, 10.95 ± 0.83 GSU, 15.57 ± 1.12 GSU and 11.94 ± 1.77 GSU, respectively. Elastic modulus E (in vitro biomechanical indicator) was negatively correlated with intermediate layer CD (r = −0.35, p = 0.01) and 2–6 mm region CD (r = −0.39, p = 0.00). A negative correlation was also found between 0-2 mm central region CD and in vivo biomechanical indicator SP-HC (r = −0.29, p = 0.02).
Conclusion: In myopic patients, densitometry is negatively correlated with biomechanical properties both in vivo and in vitro. With an increase in CD, the cornea deformed more easily.
Introduction
With the expansion of the myopia population in China, the number of refractive surgeries is also increasing year by year (Wang et al., 2020). Corneal biomechanical properties are closely related to multiple aspects of refractive surgery (Ma et al., 2018). Biomechanical inspection can help clinicians exclude patients with corneal ectasia tendencies before surgery and choose personalized surgical methods. For patients who are prone to undergo refractive regression, postoperative medications should be adjusted in a timely fashion (Qi et al., 2017). This means that a deep understanding of corneal biomechanical properties will be very helpful in clinical diagnosis and treatment (Raevdal et al., 2019). At present, corneal biomechanical characteristics can be obtained in vitro and in vivo, while in vivo measurements are more common. Corvis ST (Oculus, Wetzlar, Germany), a biomechanical instrument that has good accuracy and repeatability, is increasingly being used.
As an index of corneal morphology, corneal densitometry (CD, also known as optical density) (Charpentier et al., 2021) was obtained by a high-speed Scheimpflug Pentacam camera (Oculus, Wetzlar, Germany). By quantifying the backscatter light from the cornea, the transparency of each layer, region, and total cornea was measured. The stroma accounts for approximately 90% of the human cornea. The diameter and arrangement of the stromal collagen fibers directly affects the corneal transparency (Espana and Birk, 2020), which can be discribed as CD value (grayscale units, GSUs). The spacing, distribution and direction of the collagen fibers also impacts the biomechanical indicators. At present, knowledge of the relationship between corneal optical indicators and biomechanics is deficient. SMILE surgery provides stromal tissue that can be used for mechanical tensile testing. Therefore, we studied the correlation between biomechanical properties and CD.
Limited research has been performed on the correlation between corneal biomechanics and optical densitometry. The results of Li et al. (2020) showed that the central region CD had a significant correlation with Corvis ST SP-A1. To the best of our knowledge, there have been no studies on the correlation between the in vitro elastic modulus and optical densitometry of human corneas. We hope that our study will provide some references for clinical work through analyzing the correlation between biomechanical indicators (both in vivo and in vitro) and corneal densitometry.
Materials and methods
Participants
A retrospective study was conducted of 37 patients (63 eyes) with myopia treated in the refractive center of Tianjin Eye Hospital, Tianjin Medical University from 1 March 2017 to 1 September 2018, including 29 eyes of 16 men and 34 eyes of 21 women. All participants were between 16 and 39 years old, with an average age of 25.14 ± 6.74 years, and the equivalent sphericity (SE) was between −1.5D and −10.0D. All participants underwent routine preoperative examinations, including vision, optometry, slit-lamp microscopy, fundus examination and intraocular pressure (IOP, by Topcon). After SMILE surgery, the corneal lenticule obtained during the operation was subjected to a uniaxial tensile test to measure its biomechanical properties in vitro. The exclusion criteria were as follows: a previous history of eye surgery, previous wearing of contact lenses, a history of eye trauma, glaucoma, pterygium and other eye diseases, pregnancy or lactation, rheumatic immune system diseases and abnormal thyroid function. This study followed the principles of the Declaration of Helsinki and was approved by the Ethics Committee of the Tianjin Eye Hospital. All enrolled participants signed an informed consent form.
Methods
Corneal densitometry
The examination was performed by the same skilled technician in one dark room. Participants were told not to move or blink their eyes during the inspection. Twenty-five images for each patient were taken with a Scheimpflug camera, and the machine automatically analyzed the CD values of four circular or concentric ring areas (0–2 mm, 2–6 mm, 6–10 mm, and 10–12 mm diameter). Meanwhile, the machine will layer the cornea according to the depth and measure the CD values of each layer. The anterior layer is 120 μm thick. The posterior layer is 60 μm near the atria, while the intermediate layer is the section between the anterior and posterior layers. The instrument provided the central corneal thickness (CCT) of the participants at the same time.
Corvis ST biomechanical analyzer
The examination was performed by the same skilled technician. Participants were told to keep their eyes on the indicator light and not move or blink during the examination. The air pulse will start automatically, and the cornea will go through the process from stillness to depression and then rebound to the physiological state. The Scheimpflug camera automatically took 140 images to record the deformation process. The results of the quality index “OK” were included in the analysis, and the average value was taken from three consecutive measurements. After the test, the instrument automatically analyzed the data and calculated the corneal deformation indicators. Eight commonly used mechanical indices of Corvis ST were selected for research. The indicators are shown as follows.
Eight commonly used Corvis ST indicators.
Small incision femtosecond laser lenticule extraction (SMILE) surgery
All surgeries were successfully completed by the same experienced chief physician. Levofloxacin 0.5% (Santen, Japan) eye drops were routinely used 4 times a day for 3 days before the operation. Eye drops of obukaine hydrochloride were used for surface anesthesia prior to the operation. A VisuMax femtosecond laser (Zesis Company, Germany) was used for the surgery at frequency 500 kHz, energy 115–130 NJ. The thickness of the corneal cap was 110 μm. The diameter of the stromal lenticule was 6–6.5 mm, and the corneal cap diameter was the lenticule’s diameter plus 1 mm. Laser scanning was first performed under the lenticule, then lateral cutting, and finally scanning was performed above the lenticule. A 3 mm micro incision was made at 12 o’clock, and the base thickness was 10–15 μm. After scanning, the matrix tissue was separated and removed from the cornea and then carefully checked to determine whether the lenticule was integrated.
Uniaxial tensile test
After removing the stromal lenticule during the surgery, the 12 o’clock position was marked with gentian violet, and the tissue was immediately placed in preservation solution (Eusol-C, Alchima, Padova). The lenticule was stored in a refrigerator at 4°C for less than 24 h. The stromal tissue was placed on a rubber pad and smoothed out, then we cut the tissue from the top (12 o’clock position) to the bottom into a long strip with a 1 mm width, 6–6.5 mm length (equals to the diameter of stromal lenticule). The test piece was held on the fixture and kept wet with a normal saline water bath. Both ends were fixed with screws. First, the lenticule was preadjusted, and then the tensile force was applied using a mechanical test system (ITBC-50, Kyle Measurement and Control Test System Co., Ltd.) with a speed of 0.01 mm/s uniaxial tension. The stress corresponding to less than 5% of the strain produced was analyzed, and Young’s modulus E was obtained by linear regression analysis. This study used σ = A (eBe-1) to describe the relationship between stress and strain and E = B (σ+A) to deduce the linear relationship of the elastic modulus. The results of the elastic modulus corresponded to a stress of 0.02 MPa which resembles the intraocular pressure.
Statistical methods
This study used SPSS 26.0 statistical software for data analysis. Because some densitometry and biomechanical values were not normally distrubuted, Pearson and Spearman correlation analysis were used to exam the correlations between variables. Multiple liner regression analysis was used to analysis the relationship between the Young’s modulus E and Corvis ST parameters. p < 0.05 was considered statistically significant.
Results
A total of 37 participants (63 eyes) were included in this study, including 29 eyes of 16 men and 34 eyes of 21 women. The descriptive characteristics are shown in Table 1.
Generally speaking, the corneal densitometry was higher at the anterior layer. The CD values of each layer and region were shown in Table 2.
The stromal lenticule was removed from the central cornea, so we analyzed the in vitro elastic modulus E with three layers and the 0–6 mm central region CD values. The correlations between elastic modulus E and CD values are shown in Table 3 and Figure 1 and Figure 2. There was a negative correlation between the elastic modulus E and intermediate layer CD value (r = −0.35, p = 0.01). A significant correlation also observed between elastic modulus E and 2–6 mm annulus region CD value (r = −0.39, p = 0.00).
FIGURE 1. Correlation between the corneal densitometry of intermediate layer and Elastic modulus E (R = −0.35,p = 0.01).
FIGURE 2. Correlation between the corneal densitometry of 2–6 mm region and Elastic modulus E (R = −0.39,p = 0.00).
Correlation existed between multiple biomechanical parameters of Corvis ST and the CD values. A negative correlation was found between 0-2 mm central region CD and in vivo biomechanical indicator SP-HC (r = −0.29, p = 0.02). The results are shown in Table 4. Most corneal densitometry values which had a correlation were distributed at the 0–6 mm central region of cornea.
TABLE 4. Correlation analysis between the biomechanical parameters of Corvis ST and the corneal densitometry values.
There was a statistically significant positive correlation between age and in vitro elastic modulus E, and A1 V of Corvis ST was also positively correlated with age. CCT had a statistically significant correlation with the many Corvis ST indicators, especially with the corneal hardness parameter SP-HC. Details are shown in Table 5.
Correlation between in vitro elastic modulus E and the in vivo biomechanical indicators: there was a correlation between the elastic modulus E and the Corvis ST biomechanical indicators. Our team had previously studied and equated the two (Xue et al., 2021). The elastic modulus E of myopia can be calculated from the in vivo biomechanical parameters as follows: E = −49.47 + 9.188A1T + 9.946A1 V-0.195DA +0.077SP-HC. Refer to our previous study for the details. Our research had good consistency with the previous results (R2 = 0.543). The regression model is statistically significant (F = 2.90, p = 0.00).
Discussion
Previous researcher (Atalay et al., 2021) thought that changes in the corneal mechanical properties preceded any structural lesions. Some corneal diseases, such as keratoconus, have abnormal biomechanical properties prior to morphological corneal dilation. However, with the development of optical technology, the priority of morphological and mechanical changes has become a problem worthy of reconsideration. Recently, using a Pentacam optical densitometry to measure corneal transparency has become increasingly common. Pentacam is an anterior segment analyzer based on a Scheimpflug camera, which has the advantages of fast detection, non-invasiveness and good repeatability. Its CD measurement has been widely used in the early diagnosis of keratitis, corneal dystrophy, keratoconus, Fuchs endothelial dystrophy and other eye diseases in recent years (Guo et al., 2021). The measurement of the CD is of great value for the evaluation of corneal health.
As the main component of the corneal stroma, collagen fiber plays a decisive role in corneal transparency and load-bearing ability (Doughty and Jonuscheit, 2019). Maurice’s lattice theory (Meek and Knupp, 2015) hypothesized that the matrix collagen had the same diameter and that the spacing between fibers was less than one wavelength. This structure cancelled the scattered light, while the transmitted light was not disturbed, making the tissue transparent. The reason for corneal opacity was that the arrangement of the collagen fibers changed (Wilson et al., 2022), and the secondary reflected waves interfered with each other. This hypothesis has some defects, but it can still explain many phenomena. The corneal structure and arrangement of the fibers also determines the mechanical properties of the cornea (Wang et al., 2020). As a heterogeneous and viscoelastic tissue, the cornea has complex mechanical characteristics. The viscosity is closely related to the water content of the stroma (Dackowski et al., 2020), and the elasticity of the cornea is determined by fibers (Goh et al., 2007).
Many diseases, especially in the early stage, have an abnormal CD and biomechanical properties at the same time. Shen et al. (2019) found that the CD value may be related to the severity of keratoconus and negatively related to the SP-A1 of Corvis ST. Patients with Fabry disease have increased CD and abnormal biomechanical properties (Guo et al., 2021), but there has been no report of any correlation analysis between the two. At present, the finite element model that simulates the mechanics and optics of the human cornea needs to be improved. The refractive index of the cornea and various mechanical properties have not been integrated into the model, so the predictability of the results needs to be improved. By discussing the relationship between CD and biomechanical characteristics, we hope that our study will provide references for the optimization of corneal mechanics and optical models.
The development of SMILE surgery not only provides a better refractive correction method for myopia patients but also provides an opportunity for in vitro biomechanical research into the cornea. As a viscoelastic material, the Young’s modulus of the cornea is not invariable, and the values reported in previous studies were quite different. This difference may be partly due to the non-linear elasticity of the cornea. The materials used in previous studies were cadaver eyes or diseased corneal tissue (Zeng et al., 2001), which can affect the accuracy of the results.
The materials and experiments on the matrix lens in our study were completed within 24 h after the surgery. The elastic modulus corresponding to 0.02 MPa was selected as the result of the elastic modulus E because the stress was similar to the value of the physiological intraocular pressure. Our results showed that there was a negative correlation between the elastic modulus E and CD of the intermediate layer (r = −0.35, p = 0.01). The CD of the anterior layer measured by Pentacam corresponds to the uppermost 120 μm. The posterior layer is the 60 μm cornea next to the anterior chamber, and the intermediate layer is the part between the two layers. This layer was exactly the area corresponding to the matrix tissue removed by surgery. The higher the CD in this layer is, the more irregular the arrangement of the collagen fibers, and a certain level of decline in mechanical properties may also occur. A negative correlation also existed between the 2–6 mm region CD and the elastic modulus E.
An in vitro tensile test can obtain the stress–strain curve of the cornea, which is considered to be the gold standard of mechanical measurements (Spiru et al., 2017). However, this test is destructive and cannot be widely applied, so in vivo biomechanical measurements are more feasible in the clinic.
At present, the Corvis ST has been widely used in the preoperative examination before refractive surgery. In the correlation analysis between CD and the in vivo biomechanical parameters, CD was related to many in vivo biomechanical indicators, including A1 V, PD, integrated radius, SP-HC, and DA ratio 2 mm. The related areas were mostly concentrated in the central 0–6 mm region. SP-HC is the index of corneal hardness, which means the pressure at the maximum corneal compress divided by the displacement. PD refers to the distance between the two highest points on corneal rhinometry and the temporal side at the maximum depression. A stiffer cornea has a larger SP-HC and a shorter PD. In this study, SP-HC was negatively correlated with the 0–2 mm central region CD, and PD was positively correlated. This suggested that there is a negative correlation between corneal hardness and CD. This relationship also exists in keratoconus patients (Shen et al., 2019). In participants with myopia, we assumed that the irregular arrangement of collagen fibers may be the reason for the CD increase. Meanwhile, this irregular structure may also reduce the mechanical properties of the cornea.
CCT did not correlate with the elastic modulus E in vitro but was correlated with the Corvis ST indicators A1 V, integrated radius, SP-A1, and SP-HC. This is consistent with previous studies (Wang et al., 2017). The elastic modulus E refers to the stress applied on the specimen, so it is not affected by CCT naturally. However, the Corvis ST indicators was greatly affected by CCT. The thicker the cornea is, the higher the corneal stiffness and the smaller the shape variability. Therefore, we can conclude that the elastic modulus E in vitro is not affected by corneal thickness. But the influence of CCT should be considered when referring to in vivo biomechanical parameters.
The results of our study suggested that there was a statistically significant positive correlation between age and in vitro elastic modulus E and a positive correlation between age and the Corvis ST indicator A1 V. With increasing age, the corneal collagen fibers become thicker, the fiber spacing decreases, and the glycosaminoglycan content in the extracellular matrix decreases, resulting in an increase in corneal stiffness and a decrease in viscosity (Blackburn et al., 2019). The increase in collagen cross-linking caused by aging may be the reason for the improvement in corneal mechanical properties.
A limitation of this study was that the data included were limited; only 63 eyes were enrolled from young people aged 16-39. Because it was a retrospective study, comprehensive corneal biomechanical (CBI) data were lacking. The physical examination of this study was carried out between 8 a.m. and 5 p.m., and there was no fixed examination time, so the results may have systematic variations. However, considering that the research target was the correlation analysis, our results still have good reliability.
In conclusion, there is a certain degree of correlation between corneal densitometry and biomechanical parameters. The higher CD of myopic participants indicates a more irregular arrangement of collagen fibers, which means that their corneal mechanical properties may decline. Hospitals lacking biomechanical instruments can pay appropriate attention to the changes in CD. Research on CD has been increasing recently. However, there are still many problems that need to be studied further and explored in more depth, such as whether the content of glycosaminoglycan in the cornea affects CD value, will CD value of myopic patients differ from normal people, and so on.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the Tianjin Eye Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
Study concept and design YZ; Data collection YZ, CX, and XC; Analysis and interpretation of data YZ, JW, and XW; Writing the manuscript YZ; Critical revision of the manuscript YZ and YW; Administrative, technical, or material support YW; Supervision YW.
Funding
Neither author has a financial or proprietary interest in any material or method mentioned. This study was supported by the National Natural Science Foundation of China (NO.82271118) and Tianjin Health Science and Technology Project (TJWJ2021QN072).
Acknowledgments
The authors are grateful for the technical support of the Tianjin Eye Hospital and Eye Institute, Tianjin Key Lab of Ophthalmology and Visual Science.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Atalay, E., Özalp, O., and Yıldırım, N. (2021). Advances in the diagnosis and treatment of keratoconus. Ther. Adv. Ophthalmol. 13, 251584142110127. doi:10.1177/25158414211012796
Blackburn, B. J., Jenkins, M. W., Rollins, A. M., and Dupps, W. J. (2019). A review of structural and biomechanical changes in the cornea in aging, disease, and photochemical crosslinking. Front. Bioeng. Biotechnol. 7, 66. doi:10.3389/fbioe.2019.00066
Charpentier, S., Keilani, C., Maréchal, M., Friang, C., De Faria, A., Froussart-Maille, F., et al. (2021). Corneal haze post photorefractive keratectomy. J. Fr. Ophtalmol. 44, 1425–1438. doi:10.1016/j.jfo.2021.05.006
Dackowski, E. K., Lopath, P. D., and Chuck, R. S. (2020). Preoperative, intraoperative, and postoperative assessment of corneal biomechanics in refractive surgery. Curr. Opin. Ophthalmol. 31, 234–240. doi:10.1097/ICU.0000000000000663
Doughty, M. J., and Jonuscheit, S. (2019). Corneal structure, transparency, thickness and optical density (densitometry), especially as relevant to contact lens wear-a review. Cont. Lens Anterior Eye 42, 238–245. doi:10.1016/j.clae.2018.11.014
Espana, E. M., and Birk, D. E. (2020). Composition, structure and function of the corneal stroma. Exp. Eye Res. 198, 108137. doi:10.1016/j.exer.2020.108137
Goh, K. L., Meakin, J. R., Aspden, R. M., and Hukins, D. W. (2007). Stress transfer in collagen fibrils reinforcing connective tissues: Effects of collagen fibril slenderness and relative stiffness. J. Theor. Biol. 21, 305–311. doi:10.1016/j.jtbi.2006.10.008
Guo, H., Ren, S., and Zhao, D. (2021). Application of corneal densitometry in clinical practices. Chin. J. Exp. Ophthalmol. 39, 454–458. doi:10.3760/cma.j.cn115989-20200723-00523
Li, Y. Z., Liang, G., and Zhang, J. Y. (2020). Correlation between corneal biomechanics and corneal densitometry. Int. Eye Sci. 20, 328–331. doi:10.3980/j.issn.1672-5123.2020.2.32
Ma, J., Wang, Y., Wei, P., and Jhanji, V. (2018). Biomechanics and structure of the cornea: Implications and association with corneal disorders. Surv. Ophthalmol. 63, 851–861. doi:10.1016/j.survophthal.2018.05.004
Meek, K. M., and Knupp, C. (2015). Corneal structure and transparency. Prog. Retin. Eye Res. 49, 1–16. doi:10.1016/j.preteyeres.2015.07.001
Qi, H., Gao, C., Li, Y., Feng, X., Wang, M., Zhang, Y., et al. (2017). The effect of timolol 0.5% on the correction of myopic regression after LASIK. Med. Baltim. 96, e6782. doi:10.1097/MD.0000000000006782
Raevdal, P., Grauslund, J., and Vestergaard, A. H. (2019). Comparison of corneal biomechanical changes after refractive surgery by noncontact tonometry: Small-incision lenticule extraction versus flap-based refractive surgery - a systematic review. Acta Ophthalmol. 97, 127–136. doi:10.1111/aos.13906
Shen, Y., Han, T., Jhanji, V., Shang, J., Zhao, J., Li, M., et al. (2019). Correlation between corneal topographic, densitometry, and biomechanical parameters in keratoconus eyes. Transl. Vis. Sci. Technol. 8, 12. doi:10.1167/tvst.8.3.12
Spiru, B., Kling, S., Hafezi, F., and Sekundo, W. (2017). Biomechanical differences between femtosecond lenticule extraction (FLEx) and small incision lenticule extraction (SmILE) tested by 2D-extensometry in ex vivo porcine eyes. Invest Ophthalmol. Vis. Sci. 58, 2591–2595. doi:10.1167/iovs.16-20211
Wang, W., He, M., He, H., Zhang, C., Jin, H., and Zhong, X. (2017). Corneal biomechanical metrics of healthy Chinese adults using Corvis ST. Cont. Lens Anterior Eye 40, 97–103. doi:10.1016/j.clae.2016.12.003
Wang, Y., Shi, W., and Li, Y. (2020). The rapid development and changes of corneal refractive surgery in China. Zhonghua Yan Ke Za Zhi 56, 81–85. doi:10.3760/cma.j.issn.0412-4081.2020.02.001
Wilson, S. E., Sampaio, L. P., Shiju, T. M., Hilgert, G. S. L., and de Oliveira, R. C. (2022). Corneal opacity: Cell biological determinants of the transition from transparency to transient haze to scarring fibrosis, and resolution, after injury. Invest Ophthalmol. Vis. Sci. 63, 22. doi:10.1167/iovs.63.1.22
Xue, C., Wang, Y., and Xiang, R. Q. (2021). Correlation of corneal biomechanical properties of myopic patients in vivo and in vitro. Recent Adv. Ophthalmol. 41, 943–947. doi:10.13389/j.cnki.rao.2021.0198
Keywords: corneal densitometry, Corvis ST, small incision lenticule extraction (SMILE), uniaxial tensile test, biomechanic
Citation: Zheng Y, Xue C, Wang J, Chen X, Wang X and Wang Y (2023) Analysis of the correlation between biomechanical properties and corneal densitometry in myopic eyes. Front. Bioeng. Biotechnol. 11:1182372. doi: 10.3389/fbioe.2023.1182372
Received: 08 March 2023; Accepted: 04 April 2023;
Published: 26 April 2023.
Edited by:
Chiara Giulia Fontanella, University of Padua, ItalyCopyright © 2023 Zheng, Xue, Wang, Chen, Wang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yan Wang, wangyan7143@vip.sina.com