Abstract
Secondary benzyl alcohols (α-OH) in lignin β-O-4 linkages are reactive sites that can trigger undesirable condensation reactions during lignin valorization processes. The selective oxidation of α-OH groups can effectively inhibit lignin condensation reactions, thereby promoting lignin depolymerization. Meanwhile, this approach increases the functionalities, which is beneficial for the design of lignin-derived functional materials. In this study, we report a chemoselective oxidation method using 2-iodoxybenzoic acid (IBX) to selectively convert α-OH groups into ketones. The reaction proceeds under ambient conditions with green DMSO as the sole solvent. Oxidation of β–O–4 model compound (veratrylglycerol-β-guaiacyl ether, VG) revealed exceptional selectivity of this oxidation process, achieving 92% substrate conversion with 90% yield of the corresponding α-ketone product after 48 h. When applied to cellulolytic enzyme lignin, this protocol significantly increased the α-ketone content from 1.2% to 20.5% (relative to total interunit linkages), as identified and quantified by 2D HSQC NMR. The straightforward α-OH oxidation protocol demonstrates significant potential for precise lignin structure modification and rational design of lignin-derived functional materials in sustainable polymer systems.
1 Introduction
Lignocellulosic biomass, composed primarily of three polyol biopolymers—cellulose (40%–60%), hemicellulose (15%–35%), and lignin (20%–30%) (; ), has garnered extensive attention as a renewable resource for producing sustainable chemicals and functional materials (; ). Selective oxidation of lignocellulosic biopolymers to introduce highly reactive functional groups, such as ketone, into their polymer structures facilitates their downstream processing, selective modification, and grafting with other polymers (; ).
Lignin is a complex aromatic polymer composed of phenylpropyl structural units (; ; ; ; ). Among its inter-unit linkages, the β-aryl ether (β–O–4) bond is the most predominant, accounting for over 50%–80% of the total chemical connections (; ; ; ). The β-aryl ether structure features a primary alcohol at Cγ position (γ-OH) and a secondary benzyl alcohol at Cα position (α-OH). Selective oxidation of α-OH to carbonyl group in β–O–4 unit can effectively inhibit undesirable condensation reactions (; ; ) and promote the depolymerization of lignin by weakening the ether bond. Besides, the oxidation of the α-OH in β-aryl ether generates carbonyl group, thereby enhancing the functionality of lignin polymer and providing a specific reaction site for chemical modification, such as nucleophilic addition with various reagents (e.g., amines or alcohols). Existing oxidation systems, including the processes by using 2,3-dichloro-5,6-dicyano-1,4-benzoquinone (DDQ) (), TEMPO (), and N-hydroxyphthalimide (NHPI) () demonstrated effectiveness in the selective oxidation of α-OH in β–O–4 units. However, these reaction systems are often complex, typically requiring the use of co-catalysts or elevated reaction conditions ().
2-lodoxybenzoic acid (IBX) is a representative high-valence organic iodine reagent widely utilized for selective oxidation of reactive alcohols. The oxidation reaction is typically conducted in polar solvents such as dimethyl sulfoxide (DMSO) to stabilize the intermediates. Additionally, owing to its single-electron transfer mechanism, IBX-mediated oxidation generates fewer by-products and is easily operated (). These versatile characteristics prompted us to explore the potential of IBX for selective oxidation of polyol-containing lignin biopolymers.
In this work, IBX was employed as an oxidizing agent for the oxidation of β–O–4 model compound and lignin in DMSO under ambient conditions. The results demonstrated that IBX selectively oxidized the α-OH groups to α-ketone groups in β–O–4 units while maintained the structural integrity of γ–OH groups. This study presents a mild and facile method for the selective oxidation of lignin polymers, offering potential for precise lignin structure modification and rational design of lignin-derived functional materials.
2 Materials and methods
2.1 Chemicals and materials
2-lodoxybenzoic acid (97%), guaiacylglycerol-β-guaiacyl ether (≥97%), pyridine (≥99.5%), N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA, 99%), acetic anhydride (≥99%), dimethylsulfoxide-d6 were purchased from Shanghai Macklin Biochemical Technology Co., Ltd. Dichloromethane, sodium hydroxide (97%), veratrylglycerol-β-guaiacyl ether (VG, 95%), ethylacetate (99.5%), sodium acetate trihydrate (99%), dimethyl sulfoxide, n-phenylpentane (internal standard, 99%) were purchased from Shanghai Aladdin Biochemical Technology Co., Ltd. Sodium chloride (99.5%), was purchased from Beijing Inno Chem Science & Technology Co., Ltd. Viscozyme were purchased from Novozymes. Acetic acid, 1,4-dioxane, tetrahydrofuran (THF), and ethanol were all analytical reagent grade. The water used in all experiments was ultrapure (18.25 MΩ). All chemicals and reagents were used without further purification.
2.2 Oxidation of β-O-4 model compound
VG or GG (5 mg) and IBX (10 mg, 2.2–2.4 equivalent) were dissolved in 1 mL of DMSO. After reaction for a period of time, 1 mL of internal standard (3 mg/mL n-pentylbenzene) was added into the solution. The mixture was then transferred into a separating funnel charged with 10 mL of saturated NaCl solution. After acidifying with diluted hydrochloric acid (pH < 3), the mixture was extracted with ethyl acetate (3 × 5 mL). The combined organic phase was dried over anhydrous sodium sulphate, and concentrated under reduced pressure. The obtained crude product was dissolved in 5 mL of dichloromethane for GC-FID and GC-MS analysis. 10 μL of the crude product solution was placed into the sample bottle charged with 80 µL of BSTFA and 100 µL of pyridine. The mixture was then heated in an oven at 50 °C for 1 h. The sample is tested by a GC-MS and a GC-FID for qualitative and quantitative analysis, respectively. The product yield was calculated using the following
Equation 1(
):
Yproduct—the yield of the product (%);
nphenylpentane (mmol)—the molar amount of the internal standard (phenylpentane);
Sproduct—the peak area of product in the GC–FID chromatogram;
Cphenylpentane—the effective carbon number of phenylpentane;
nVG (mmol)—the molar amount of the internal standard (VG);
Sphenylpentane—the peak area of phenylpentane in the GC–FID chromatogram;
Cproduct—the effective carbon number of product.
After oxidation, the α-oxidized β-O-4 model compound (α-OVG) was separated using a preparative liquid chromatograph (Sepabean SP-2, Santai Technologies, Changzhou, China) equipped with a UV detector and a silica gel chromatographic column. An eluent consisting of ethyl acetate (EtOAc)/petroleum ether (1:3, v/v) was used for separation. Eventually, the crude α-OVG product was obtained with an isolated yield of 85%. The structure of the collected fraction was first confirmed by GC-MS and further characterized via NMR spectroscopy.
2.3 Cellulolytic enzyme lignin (CEL) isolation
Eucalyptus wood powder less than 40 mesh was milled by a Planetary Ball Mill for 48 h at 400 rpm. The ball milled wood powder (10 g) was then hydrolyzed in 400 mL of sodium acetate buffer (pH = 5.5) containing 300 mg of cellulosic enzyme (Celluelast 1.5L) and 400 mg of plant fiber complex enzyme (Viscozyme L) at 50 °C for 48 h (two times). After enzymatic hydrolysis, the suspension was centrifuged and washed with deionized water (2 × 60 mL). The solid residual obtained after centrifugation was extracted with 200 mL of dioxane/water (96:4, v/v) for 12 h. After centrifugation, the solution was concentrated at 45 °C under reduced pressure. The concentrated solution (∼3 mL) was droped into 50 mL of acidified deionized water (pH = 2), and the deposited lignin was obtained after centrifugation followed by vacuum drying.
2.4 Oxidation of lignin
100 mg of eucalyptus CEL and 100 mg of IBX were dissolved in 2 mL of DMSO and reacted at room temperature for 12–48 h. After reaction, the solution was added to 100 mL of DCM for deposition of the lignin. The mixture was then filtered and the solid was washed with DCM (3 × 10 mL). The oxidized lignin was finally obtained after drying in a vacuum oven.
2.5 Gel permeation chromatography (GPC) analysis
To a reaction vial, 5 mg of lignin sample was mixed with 0.5 mL of acetic anhydride and 0.5 mL of pyridine and reacted at room temperature for 24 h. The reaction mixture was then concentrated under reduced pressure at 45 °C, with anhydrous ethanol continuously added to remove residual reagents. The acetylated lignin was dissolved in 2 mL of tetrahydrofuran (THF) to achieve a concentration of 2.5 mg/mL, filtered through a 0.45 μm PTFE membrane, and subjected to GPC analysis by using an Agilent 1,260 high-performance liquid chromatography (HPLC) connecting with a refractive index detector. Elution was carried out at a column temperature of 35 °C using THF as the mobile phase at a flow rate of 1.0 mL/min. Molecular weight parameters were determined relative to polystyrene standards.
2.6 2D HSQC NMR characterization
The 2D HSQC NMR spectra of cellulolytic enzyme lignin (CEL) and oxidized CEL (OCEL) were analyzed by 600 MHz NMR spectrometer (AVANCE NEO 600 Bruker). CEL or OCEL (∼80 mg) was dissolved in 0.5 mL of DMSO-d6 in an NMR tube. The 2D HSQC NMR spectra were acquired from 10 to 0 ppm in F2 (1H) with 256 data points and a 1 s recycle delay, and 160–0 ppm in F1 (13C) with 256 increments. Total scan number of 32 was employed.
3 Results and discussion
3.1 Oxidation of β-O-4 model compounds by IBX
First, we evaluated the efficacy of IBX in oxidizing a β–O–4 model compound. The β–O–4 model compound, veratrylglycerol-β-guaiacyl ether (VG), was oxidized by IBX in DMSO solvent under ambient conditions. Post-oxidation, the resultant products were analyzed by GC-MS. The structure of main product was authenticated by mass spectrometry and 1H NMR after separation. In GC chromatograms, a typical product was identified at the retention time of 25.8 min (Figure 1A). Analyzing by mass spectrometry and 1H NMR (Figure 1B) confirmed its identity as the α-oxidized derivative of VG (α-OVG). γ-oxidized VG product (γ-OVG) was not identified from GC-MS analysis. This result demonstrates a high selectivity of IBX for oxidizing the second benzylic alcohol in the VG compound. The proposed mechanism for the oxidation of VG by IBX is depicted in Figure 1C. It is well established that the secondary alcohol at the benzylic carbon exhibits higher reactivity than the hydroxyl group at the γ-carbon, which accounts for the chemoselective oxidation occurring at the α-OH position.
FIGURE 1
To investigate the oxidation process of VG by IBX, 5 mg of substrate and 10 mg of IBX (∼2.3 equiv.) were mixed for reaction (Table 1, entry 1–4). As the reaction time was extended from 12 to 48 h, the conversion of VG increased from 81.6% to 92%, and the yield of α-OVG rose from 76.2% to 90.5%. However, when the amount of IBX was reduced to 5 mg (∼1.1 equiv.), both the VG conversion and the α-OVG yield remained low, even after a prolonged reaction time of 48 h (entry 5 in Table 1). These results suggest that the addition of approximately two equivalents of IBX is optimal for the efficient oxidation of VG. Furthermore, phenolic β–O–4 model compound (guaiacylglycerol-β-guaiacyl ether; GG) was employed as a substrate for oxidation. As shown in Table 2, only 3.3% yield of α-OGG was obtained. The major product observed was QGG. This result is attributed to the presence of a phenolic hydroxyl group and an ortho-methoxy group in GG. IBX preferentially oxidized the phenolic hydroxyl group and mediated the oxidative demethylation of the ortho-methoxy methyl group (-CH3). These results suggest that this oxidation process primarily targets the α-OH group within non-phenolic β-O-4 units.
TABLE 1
 | ||||||
|---|---|---|---|---|---|---|
| Entry | VG (mg) | IBX (mg) | Time (h) | Conversion (%) | α-OVG (%) | γ-OVG (%) |
| 1 | 5 | 10 | 12 | 81.6 | 76.2 | ND |
| 2 | 5 | 10 | 24 | 84.6 | 79.8 | ND |
| 3 | 5 | 10 | 36 | 90.8 | 88.2 | ND |
| 4 | 5 | 10 | 48 | 92.0 | 90.5 | ND |
| 5 | 5 | 5 | 48 | 48.6 | 44.5 | ND |
Oxidation of nonphenolic β–O–4 model compound (VG) by IBX.
TABLE 2
 | |||||
|---|---|---|---|---|---|
| GG (mg) | IBX (mg) | Time (h) | Conversion (%) | QGG (%) | OGG (%) |
| 5 | 10 | 48 | 89.7 | 56.2 | 3.3 |
Oxidation of phenolic β–O–4 model compound (GG) by IBX.
3.2 Oxidation of lignin polymer
Based on the successful oxidation of α-OH in β–O–4 model compound by IBX, we further investigate the effectivity of IBX for oxidation of lignin polymer. Eucalyptus cellulolytic enzyme lignin (CEL) was isolated and subjected to IBX oxidation with an equal mass (lignin/IBX: 1:1, w/w). The CEL and oxidized CEL (OCEL) were characterized by 2D HSQC NMR and GPC to elucidate the change in chemical structure and molecular weight.
As shown in 2D HSQC NMR spectra (Figure 2), CEL consists of high content of β–O–4 linkage (56.8%), lower amounts of β–5 (3.9%) and β–β (10.5%) linkages, and trace of α-oxidized β–O–4 linkage (1.2%). After IBX oxidation, the content of β–5 and β–β remained nearly unchanged. However, the content of the native β–O–4 linkage significantly decreased. In the side chain region of the OCEL spectrum, a substantial number of α-oxidized β–O–4 structures (A′) were identified, with their Cβ–Hβ and Cγ–Hγ crosslinked signals appearing at 81/5.3–5.7 ppm and 62/3.7 ppm, respectively. Additionally, α-ketoyl aromatic ring structures (G′ and S′) were detected in the conjugated aromatic ring region of the OCEL spectrum. Quantitative analysis via 2D HSQC NMR revealed that the content of α-oxidized β–O–4 structures in CEL increased markedly to 20.8% after IBX oxidation (Figure 3A), with an oxidation yield of 34.5%. These results suggests that IBX can selectively oxidize the native β–O–4 units in lignin to form α-ketoyl β–O–4 structures. The lower oxidation efficiency in lignin compared to that in VG compound can be attributed to the presence of other linkages and free phenolic hydroxyl groups in lignin polymer, which can consume a portion of IBX.
FIGURE 2
FIGURE 3
The molecular weight distributions of CEL before and after oxidation were analyzed using gel permeation chromatography (GPC) to investigate the potential degradation of lignin during the oxidation process. As shown in Figure 3B, the oxidized product (OCEL), exhibited a lower molecular weight (weight-average molecular weight of 6,582 Da), compared to the unoxidized CEL (weight-average molecular weight of 15,420 Da). This reduction in molecular weight indicates that the IBX-mediated oxidation induces the degradation of lignin polymer. This degradation can be attributed to the formation of α-ketone groups during oxidation, which destabilize the β-aryl ether linkages in the lignin structure (; ; ; ), thereby rendering the polymer more susceptible to cleavage under the IBX oxidation conditions.
Overall, our findings demonstrate that IBX serves as a versatile oxidizing agent, capable of effectively oxidizing secondary benzyl alcohols to ketone groups in the β–O–4 units of lignin. Although the relatively high cost of IBX and stoichiometric oxidation of this method may limit its large-scale application, this method offers a mild and straightforward approach for the modification or functionalization of lignin materials, particularly in scenarios where precise chemical transformations are required.
4 Conclusion
In summary, we have successfully developed a mild and facile IBX-based oxidation method capable of selectively oxidizing of secondary benzyl alcohols in the β–O–4 units of lignin polymers without the use of toxic metallic catalysts. The method provides a mild and efficient strategy for the modification and functionalization of lignin materials. Its operational simplicity and selectivity make it a promising tool for advanced biomass valorization and functional material design.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
FJ: Formal Analysis, Writing – original draft, Methodology, Writing – review and editing, Investigation. ZD: Formal Analysis, Writing – review and editing, Methodology. XL: Methodology, Writing – review and editing, Formal Analysis. ZZ: Investigation, Writing – review and editing. LS: Funding acquisition, Project administration, Writing – review and editing, Supervision. CZ: Conceptualization, Writing – review and editing, Supervision, Funding acquisition, Project administration.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by National Natural Science Foundation of China (Grant No.: 32301543), National Natural Science Foundation of China for Distinguished Young Scholars (Grant No.: 32425039).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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
1
DengW.FengY.FuJ.GuoH.GuoY.HanB.et al (2023). Catalytic conversion of lignocellulosic biomass into chemicals and fuels. Green Energy and Environ.8 (1), 10–114. 10.1016/j.gee.2022.07.003
2
DessbesellL.PaleologouM.LeitchM.PulkkiR.XuC. (2020). Global lignin supply overview and kraft lignin potential as an alternative for petroleum-based polymers. Renew. Sustain. Energy Rev.123, 109768. 10.1016/j.rser.2020.109768
3
HuangD.LiR.XuP.LiT.DengR.ChenS.et al (2020). The cornerstone of realizing lignin value-addition: exploiting the native structure and properties of lignin by extraction methods. Chem. Eng. J.402, 126237. 10.1016/j.cej.2020.126237
4
JiQ.YuX.ChenL.MustaphaA. T.OkonkwoC. E.ZhouC.et al (2023). Comprehensive depolymerization of lignin from lignocellulosic biomass: a review. Crit. Rev. Environ. Sci. Technol.53 (21), 1866–1887. 10.1080/10643389.2023.2190314
5
LanW.AmiriM. T.HunstonC. M.LuterbacherJ. S. (2018). Protection group Effects during α,γ-Diol lignin Stabilization promote high-selectivity Monomer production. Angew. Chem. Int. Ed.57 (5), 1356–1360. 10.1002/anie.201710838
6
LanW.de BuerenJ. B.LuterbacherJ. S. (2019). Highly selective oxidation and depolymerization of α,γ-Diol-Protected lignin. Angew. Chem. Int. Ed.58 (9), 2649–2654. 10.1002/anie.201811630
7
LancefieldC. S.OjoO. S.TranF.WestwoodN. J. (2015). Isolation of functionalized phenolic Monomers through selective oxidation and CO bond cleavage of the β-O-4 linkages in lignin. Angew. Chem. Int. Ed.54 (1), 258–262. 10.1002/anie.201409408
8
LiT.MaH.WuS.YinY. (2020). Effect of highly selective oxypropylation of phenolic hydroxyl groups on subsequent lignin pyrolysis: toward the lignin valorization. Energy Convers. Manag.207, 112551. 10.1016/j.enconman.2020.112551
9
LiuM.DysonP. J. (2023). Direct conversion of lignin to functionalized diaryl ethers via oxidative cross-coupling. Nat. Commun.14 (1), 2830. 10.1038/s41467-023-38534-1
10
MaityA.HyunS.-M.WortmanA. K.PowersD. C. (2018). Oxidation Catalysis by an Aerobically generated Dess–Martin Periodinane Analogue. Angew. Chem. Int. Ed.57 (24), 7205–7209. 10.1002/anie.201804159
11
PanovicI.Miles-BarrettD. M.LancefieldC. S.WestwoodN. J. (2019). Preparation and reaction of β-O-4 γ-Aldehyde-containing Butanosolv lignins. ACS Sustain. Chem. and Eng.7 (14), acssuschemeng.9b01188–12104. 10.1021/acssuschemeng.9b01188
12
RahimiA.AzarpiraA.KimH.RalphJ.StahlS. S. (2013). Chemoselective Metal-free Aerobic alcohol oxidation in lignin. J. Am. Chem. Soc.135 (17), 6415–6418. 10.1021/ja401793n
13
RenX.WangP.HanX.ZhangG.GuJ.DingC.et al (2017). Depolymerization of lignin to aromatics by selectively oxidizing cleavage of C–C and C–O bonds using CuCl2/Polybenzoxazine catalysts at room temperature. ACS Sustain. Chem. and Eng.5 (8), 6548–6556. 10.1021/acssuschemeng.7b00732
14
SternbergJ.SequerthO.PillaS. (2021). Green chemistry design in polymers derived from lignin: review and perspective. Prog. Polym. Sci.113, 101344. 10.1016/j.progpolymsci.2020.101344
15
SudarsanamP.RuijtenD.LiaoY.RendersT.KoelewijnS.-F.SelsB. F. (2020). Towards lignin-derived chemicals using Atom-efficient Catalytic Routes. Trends Chem.2 (10), 898–913. 10.1016/j.trechm.2020.07.011
16
SunZ.FridrichB.de SantiA.ElangovanS.BartaK. (2018). Bright side of lignin depolymerization: toward New platform chemicals. Chem. Rev.118 (2), 614–678. 10.1021/acs.chemrev.7b00588
17
VangeelT.SchutyserW.RendersT.SelsB. F. (2018). Perspective on lignin oxidation: advances, challenges, and future directions. Top. Curr. Chem.376 (4), 30. 10.1007/s41061-018-0207-2
18
XieS.KumagaiS.KimY.-M.SaitoY.YoshiokaT. (2024). Combining hierarchical clustering analysis with a simplex lattice mixture design in rapidly identifying the pyrolytic interactions and predicting the product yields during the co-pyrolysis of cellulose, xylan, and milled wood lignin. Chem. Eng. J.493, 152434. 10.1016/j.cej.2024.152434
19
YamadaR.HasunumaT.KondoA. (2013). Endowing non-cellulolytic microorganisms with cellulolytic activity aiming for consolidated bioprocessing. Biotechnol. Adv.31 (6), 754–763. 10.1016/j.biotechadv.2013.02.007
20
YangG.GongZ.LuoX.ChenL.ShuaiL. (2023). Bonding wood with uncondensed lignins as adhesives. Nature621 (7979), 511–515. 10.1038/s41586-023-06507-5
21
YuX.WeiZ.LuZ.PeiH.WangH. (2019). Activation of lignin by selective oxidation: an emerging strategy for boosting lignin depolymerization to aromatics. Bioresour. Technol.291, 121885. 10.1016/j.biortech.2019.121885
22
YuanZ.DaiW.ZhangS.WangF.JianJ.ZengJ.et al (2022). Heterogeneous strategies for selective conversion of lignocellulosic polysaccharides. Cellulose29 (6), 3059–3077. 10.1007/s10570-022-04434-8
23
ZhangC.WangF. (2020). Catalytic lignin depolymerization to aromatic chemicals. Accounts Chem. Res.53 (2), 470–484. 10.1021/acs.accounts.9b00573
24
ZhangC.LiH.LuJ.ZhangX.MacArthurK. E.HeggenM.et al (2017). Promoting lignin depolymerization and restraining the condensation via an Oxidation−Hydrogenation strategy. ACS Catal.7 (5), 3419–3429. 10.1021/acscatal.7b00148
25
ZhuY.LiaoY.LuL.LvW.LiuJ.SongX.et al (2023). Oxidative catalytic fractionation of lignocellulose to high-yield aromatic aldehyde monomers and pure cellulose. ACS Catal.13 (12), 7929–7941. 10.1021/acscatal.3c01309
Summary
Keywords
biomass, lignin, selective oxidation, structural modification, lignin condensation, lignin depolymerization
Citation
Jiang F, Dong Z, Lou X, Zhang Z, Shuai L and Zhao C (2025) Selective and mild oxidation of α-hydroxyl groups in lignin by 2-iodoxybenzoic acid. Front. Bioeng. Biotechnol. 13:1654739. doi: 10.3389/fbioe.2025.1654739
Received
26 June 2025
Accepted
22 August 2025
Published
09 September 2025
Volume
13 - 2025
Edited by
Xiaojun Shen, Beijing Forestry University, China
Reviewed by
Wu Lan, South China University of Technology, China
Yongzhuang Liu, Northeast Forestry University, China
Updates
Copyright
© 2025 Jiang, Dong, Lou, Zhang, Shuai and Zhao.
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: Li Shuai, lishuai@fafu.edu.cn; Chengke Zhao, zhaock@fafu.edu.cn
Disclaimer
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.