- 1Department of Radiation Oncology, The First Affiliated Hospital of Xi’an Jiao Tong University, Xi’an, China
- 2Department of Radiation Oncology, Hanzhong Center Hospital, Hanzhong, China
- 3Department of Radiation Oncology, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 4Department of Radiation Oncology, The Second Affiliated hospital of Xi’an Jiao Tong University, Xi’an, China
- 5Department of Radiation Oncology, First Affiliated Hospital of Air Force Medical University, Xi’an, China
- 6Department of Radiation Oncology, The 940th Hospital of Joint Logistics Support force of Chinese People’s Liberation Army, Lanzhou, China
- 7Department of Radiation Oncology, The Second Affiliated Hospital of Shaanxi University of Chinese Medicine, Xianyang, China
Objective: The study aims to investigate if a relationship exists between vaginal doses and vaginal stenosis (VS) using posterior–inferior border of symphysis (PIBS) points and the International Commission on Radiation Units-Rectum (ICRU-R) point evaluation system for definitive radio(chemo)therapy in locally advanced cervical cancer.
Methods and Materials: From a vaginal dose study in China, 351 patients were prospectively assessed. For every reference point of the PIBS system and ICRU-R point was calculated for all BT and summed with EBRT. Pearson’s chi-square test and Student’s unpaired t-test compared variables with and without vaginal stenosis (VS) G ≥2. The risk factors were assessed for VS G ≥2 in multi- and univariate analyses through Cox proportional hazards model followed by a dose–effect curve construction. The VS morbidity rate was compared via the log-rank test using the median vaginal reference length (VRL).
Results: The patients (38-month median follow-up) had 21.3% three-year actuarial estimate for VS G ≥2. Compared to G <2 patients, VS G ≥2 patients received higher doses to PIBS points except for PIBS − 2 and had significantly shorter VRL. VRL (HR = 1.765, P = 0.038), total EBRT and BT ICRU-R point dose (HR = 1.017, p = 0.003) were risk factors for VS. With VRL >4.6 cm, the 3-year actuarial estimate was 12.8% vs. 29.6% for VRL ≤4.6 cm. According to the model curve, the risks were 21, 30, and 39% at 75, 85, and 95 Gy, respectively (ICRU-R point dose).
Conclusions: PIBS system point doses correlated with late vaginal toxicity. VRL combined with both EBRT and BT dose to the ICRU-R point contribute to VS risk.
Introduction
The vagina is a target organ and organ at risk (OAR) in radiation therapy for cervical cancer (CC). To date, morbidity rates due to vaginal conditions for grades II and III or higher vaginal stenosis (VS) in patients of CC treated with radio(chemo)therapy were 17–41.7 and 1–17.6% respectively (1–5). VS is commonly defined as the loss of elasticity and mucosal atrophy, shortening and narrowing of the vaginal canal, and pain during exam/intercourse (6, 7). It can seriously worsen the quality of life of patients and even recurrent tumors cannot be found due to inability to open a vaginal speculum. A few studies have revealed that the dose of ICRU-R point and D2cc in the vagina as the predictors of VS occurring after pelvic radiotherapy, besides, is no clear specification regarding accepted dosimetry, vaginal radiation tolerance, and the limitations of the vagina as an OAR (8).
The PIBS point system proposed for the first time at the Medical University, Vienna is a newly defined method of vaginal dose-reporting for both EBRT (external beam radiotherapy) and BT (brachytherapy) (9). This method contains different points total (BT + EBRT) dose administered to the lower, upper, and middle vaginal parts. This system of evaluation can be applied for 2D/3D BT. We earlier studied whether this method can be applied in patients, and found that Chinese patients with a shorter median vaginal reference length (VRL) of <4.5 cm received relatively higher doses of radiation at the PIBS − 2 cm, PIBS, and PIBS + 2 cm points than those of American and European patients (10). In this study, we further investigated if a relationship exists between vaginal doses and VS using PIBS points and ICRU-R point evaluation system in locally advanced CC treated with definitive radio(chemo)therapy.
Methods
Patient Characteristics, Treatment Modality, and Endpoint Assessment
Three hundred fifty-one patients with cervical cancer FIGO stage IB–IVB except IIIA were recruited between December 2016 and June 2018 from six radiation centers in China with locally advanced CC participated in the clinical trial (NCT03257475) and were treated with 45–50 Gy 3D-conformal radiotherapy or IMRT (intensity-modulated radiotherapy), with weekly simultaneous administration of 25 mg/m2 cisplatin (11), or a combination of 25 mg/m2 cisplatin (days 1–3) with 135 mg/m2 liposome paclitaxel (day 1). High dose rate BT with two T&O implants were conducted separately, was administered doses of 24 or 30 Gy in four to five fractions. For treatment, the method of point A dose-evaluation was applied, and the dose of loading was according to the Manchester system of BT. Doses were reported in total EBRT and BT EQD2 dose using the linear quadratic model (α/β = 3 Gy for vaginal, α/β = 10 Gy for tumor) (12).
VS was scored in accordance with the CTCAE v3.0 (Common Terminology Criteria for Adverse Events v3.0) and prospectively assessed at baseline by the physician, at three months interval after the treatment ended in the first two years and at six months interval in the years three to five, and thereafter, annually. As per CTCAE v3.0, the definition of G1 is narrowing and/or shortening but not causing functional interference, G2 as narrowing and/or shortening along with functional impairment, and G3 as absolute vaginal obliteration and unmenable to correction by surgery (5). Here, the definition of primary outcome was the time from the completion of radiotherapy to first G ≥2 VS occurrence before local recurrence or death. Censored patients were those with no documented G ≥2 VS before recurrence for the first time or death on the date of the first relapse or last visit.
Data Collection and Definitions
All reference points of the PIBS system proposed by the Westerveld group have been defined in our previous study (8). ICRU-R and ICRU-B points were also done based on the ICRU guidelines (Report 89) and the ICRU (Report 38) (13, 14). The ICRU-R point is associated with the applicator and positioned 5 mm on an anterior–posterior line at the back of the posterior vaginal wall drawn from the center of the vaginal source. The dose of A point, PIBS, ICRU-R, PIBS ± 2 cm, and ICRU-B points were recorded. Meanwhile, the value of VRL and other clinical data of patients were all obtained. A CT scan of patients was done at the time of the first and third fractions (BT1, BT3). The parameters of BT2, 4, and 5 were recorded by averaging the BT1 and BT3 dose values.
Statistical Analysis
This study identified a set of variables in the dosimetric and clinical data for predicting of VS G ≥2 occurrence. Based on the distribution, data were presented as mean and SD (standard deviation) or median and IQR (interquartile range). Proportions were presented as the percentage of patients without and with the characteristic. The differentiation of continuous variables was done by student’s (unpaired) t-test and the comparison of categorical variables was done using Pearson’s chi-square test. Several a priori chosen, clinically relevant characteristics of disease-, patient-, and treatment were assessed as risk factors for VS G ≥2 in uni- and multivariate analyses by Cox regression model. Ninety-five percent CI (confidence intervals) and HR (Hazard Ratios) were estimated. Analysis of data at the time of the event was done by Kaplan–Meier analyses. The level of two-sided significance was set at 5%. Analyses of statistical data were conducted using SPSS v 23 (SPSS Inc., USA).
Results
For this study, the time for a median follow-up was three years and two months (IQR 36–42 months). The disease, demography, and treatment characteristics of 351 patients are shown in Table 1. The median age was 50 years (range 31–60 years). Most patients (64.7%) had stages IIB (227/351) or IIIB (63/351, 17.9%) lesions. At the time of diagnosis, among 351 patients, 52 (14.8%) had a vaginal involvement according to clinical assessment and imaging tests. Most lesions were observed to be confined to the fornix of the vagina. The median VRL was 4.6 (IQR 3.8–5.2) cm. The doses of median total (BT + EBRT) at PIBS, PIBS − 2 cm, and PIBS + 2 cm were 55.8 (IQR 54.2–59.3) Gy3, 2.4 (IQR 2–3.1) Gy3, and 80.9 (IQR 65.8–102.8) Gy3, respectively. At point A, the median dose administered was 84.1 (IQR 81–87.5) Gy10. In 44 patients (12.5%), vaginal stenosis was not found (G0); 213 patients were observed to have G1 (60.6%); 90 G2 (25.6%) and 4 G3 (1.1%). The 3-year actual estimate VS G ≥2 was 21.31% (Figure 1).
An analysis of factors as per the morbidity groups revealed different vaginal doses and VRL in two groups (Table 2). The results show that doses delivered to PIBS + 2, PIBS, and PIBS − 2 were higher (108.1 ± 55.4 Gy vs. 90.2 ± 40.8 Gy, P = 0.0059, 59.49 ± 7.1 6 vs. 57.8 ± 7. P = 0.067 and 2.8 ± 1.2vs. 2.6 ± 1.0 Gy P = 0.099) for patients having VS G ≥2 compared to those delivered in patients with G <2. Likewise, doses delivered to ICRU-R and ICRU-B points were significantly larger (85.8 ± 18.1 Gy vs. 75.7 ± 12.47Gy, P <0.0001 and 83.2 ± 13.6Gy vs. 76.1 ± 12.4 Gy, P = 0.007) for patients with VS G ≥2 compared with patients with G <2. However, values of VRL were significantly shorter for patients with VS G ≥2 compared to those in patients having G <2 (4.25 ± 0.99 cm vs. 4.69 ± 0.97 cm, P = 0.0003).
Furthermore, in the multivariable model, VS risk factors related to an increase in the risk of G ≥2 VS were the following: values of VRL (binary, ≤4.6 cm vs. >4.6 cm, HR = 1.765, 95% CI 1.033–3.016, P = 0.038), and total EBRT and BT ICRU-R point dose (continuous EQD2 in Gy, HR = 1.017, 95% CI 1.006–1.028, p = 0.003). An overview of the uni- and multivariate Cox regression model is presented in Table 3.
In VS G ≥2 group, most patients (69.7%) had VRL of ≤4.6 cm, while only 30.3% had VRL >4.6 cm. The probability of occurrence of VS G ≥2 according to this length cut-off is mentioned in Figure 2. With longer VRL (>4.6 cm), the three-year actuarial estimate was 12.8%, while with shorter VRL (≤4.6 cm), the value was 29.6%.
Figure 2 Actuarial estimates for VS G ≥2 in patients VRL ≤4.6 cm in comparison to those with VRL >4.6 cm.
In the relationship between dose–effect of VS with an increase in the dose to the ICRU-R point, there was a significant increase in the probability of VS G ≥2 to increase (p <0.001). A dose–effect curve according to the risk determined by univariate Cox regression test is shown in Figure 3 with estimation for VS G ≥2 and ICRU-R point dose in EQD2 (HR = 1.021 per 1 Gy with 95% CI 1.012–1.032). As per the model for dose–effect, the probability to develop VS G ≥2 were 21, 30, and 39% with ICRU-R point doses of 75, 85, and 95 Gy, respectively.
Figure 3 The relationship of dose–effect of the combined brachytherapy and EBRT dose to ICRU-R point in EQD2 and VS G≥2 in (N = 351) patients. Three black points represent ICRU-R point doses of 75, 85, and 95 Gy, respectively.
Discussion
In this multi-center, prospective study using the PIBS points and ICRU-R point evaluation system, we found VRL and ICRU-R point-dose (total EBRT + BT) as two risk factors for VS G ≥2 in CC post-definitive radio(chemo)therapy. With shorter VRL and an increase in dose to the ICRU-R point, there was a significant increase in the probability of the occurrence of VS.
In this study, most patients received IMRT + 2D BT and the median dose of point A in EQD2 was 84.1 Gy. Although three-dimensional and IGABT (image-guided adaptive brachytherapy) has been recommended by EMBRACE study (15), many hospitals around the world cannot promote this technology due to the limitations of equipment and personnel, so the two-dimensional and A point-dose evaluation system is still an important method in BT of cervical cancer. Furthermore, the rate of morbidity due to VS after definitive radiotherapy was similar to that obtained in other studies that applied IGABT. Severe or no morbidity VS is rare, and most studies have reported mild-to-moderate morbidity (3, 16–19). In a few earlier retrospective studies, patterns of standard loading directed to point A was used and attempted to define a maximum vaginal mucosa tolerance dose. Most of these reports on morbidity due to vaginal side effects are from mono-institutional, retrospective studies and have used varying systems for side-effect grading, and therefore, cannot be used for comparing with the results of this study.
As far as we are aware, this report is the first on the use of the PIBS system to evaluate the dose–effect of VS in patients with CC. The PIBS system is a new, more detailed assessment method of vaginal radiation dose, which was recently proposed by the EMBRACE study and ICRU-89 document (12, 20). This point dosimetry evaluation system has many advantages compared to traditional ICRU-R point, such as the inclusion of dose points in different vaginal parts, both in the region of BT high dose and in the lower and mid vaginal points, so, it can show the dose administration along the cranio-caudal region of the vagina. Besides, this system also includes dose information of both BT and EBRT (8). In the PIBS system, the dose of PIBS series points is an important index to evaluate the vaginal dose. However, so far there is no study on the relationship between the dose of PIBS points and vaginal toxicity in pelvic radiotherapy. In this study, we observed that the dose delivered to PIBS + 2 was 108.1 ± 55.4 Gy3 for the patients with VS G ≥2, but only 90.2 ± 40.8 Gy3 with patients with VS G = 0/1. Meanwhile, doses delivered to PIBS and PIBS − 2 for patients with VS G ≥2 were also higher than those in patients with VS G = 0/1, although the differences were not statistically significant.
VRL is another important parameter for the PIBS evaluation system. It was found to be a strong risk factor for the morbidity of VS in this study. The cut-off value was set at 4.6 cm, which is the median of VRL, similar to the results of a previous dose study on many patients (10). The incidence rate of VS G ≥2 was 2.3 times higher in patients with VRL ≤4.6 cm compared to VRL >4.6 cm. We found that because of having shorter VRL, Chinese patients with CC received higher doses of radiation for total EBRT and BT in the vaginal upper two-thirds region. In univariate Cox regression, the dose of PIBS + 2 was a high-risk factor for VS. However, doses of PIBS points were not the risk factors in the further analysis using the multi-Cox regression model. A few facts can explain the results. First, PIBS points dose, especially PIBS + 2 and PIBS are significantly affected by VRL, so there is strong correlation of each other. Second, VRL in Chinese patients is shorter than that in patients of European and American. Therefore, part of vaginal that received higher dose is large. It may be a risk factor for VS. However, the PIBS point system only involves the dosimetric evaluation, so it cannot assess the dose-volume of the vagina, although, this still needs further study to verify.
Currently, the ICRU-R point is designated as a representative dose point to assess the intermediate dose for the upper portion of the vagina (13). In this study, the dose of ICRU-R point was found to be higher in VS G ≥2 group than in the G = 0/1 group. It is further confirmed that the dose of ICRU-R point is a strong risk factor for morbidity due to VS in CC patients who received radiotherapy treatment by uni- and multivariate Cox regression models. The association between dose and effect indicated that the possibility to develop VS G ≥2 is 21% with an ICRU-R point dose of 75 Gy, which is higher than the findings of the EMBRACE study (7). The reason for the difference in results may be because we used the A points planning not the three-dimensional planning.
In conclusion, doses of PIBS system points were associated with late vaginal toxicity. VRL and a dose of the combination of EBRT and BT to the ICRU-R point added to the VS risk. Physicians should pay attention to VS when the VRL of the patient is ≤4.6 cm. Besides, using a dose of ICRU-R point ≤75 Gy EQD2 (EBRT + BT dose) could reduce the risk of severe VS.
Data Availability Statement
The raw data supporting the conclusions 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 First Affiliated Hospital of Xi’an Jiaotong University (No. XJTU1AF2017LSK-11). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
Data collection, statistical analysis, data analysis and interpretation, manuscript preparation: JW, K-sZ, and TW. Conception and design of study, data analysis and interpretation, senior reviewer: ZL. Data analysis and interpretation: R-hW and WY. Data collection and interpretation, reviewer: F-qZ, LY, Y-lW, L-cW, MS, SL, B-gL, FS, JS, Q-yZ, and JZ. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
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
1. Susko M, Craciunescu O, Meltsner S, Yun Y, Beverly S, Jing C, et al. Vaginal Dose is Associated With Toxicity in Image Guided Tandem Ring or Ovoid-Based Brachytherapy. Int J Radiat Oncol Biol Phys (2016) 94(5):1099–105. doi: 10.1016/j.ijrobp.2015.12.360
2. Brand AH, Bull CA, Cakir B. Vaginal Stenosis in Patients Treated With Radiotherapy for Carcinoma of the Cervix. Int J Gynecol Cancer (2006) 16(1):288–93. doi: 10.1111/j.1525-1438.2006.00348.x
3. Kirchheiner K, Nout RA, Tanderup K, Lindegaard JC, Westerveld H, Haie-Meder C, et al. Manifestation Pattern of Early-Late Vaginal Morbidity After Definitive Radiation (Chemo) Therapy and Image-Guided Adaptive Brachytherapy for Locally Advanced Cervical Cancer: An Analysis From the EMBRACE Study. Int J Radiat Oncol Biol Phys (2014) 89(1):88–95. doi: 10.1016/j.ijrobp.2014.01.032
4. Zhang KS, Liu Z, Wang T, Wang J, Shi F, Wang RH, et al. Occurrence of Radiation Induced Injury in Vagina After Radical Radiotherapy of Cervical Cancer and its Affecting Factors. Zhonghua Fu Chan Ke Za Zhi (2018) 53(4):257–62. doi: 10.3760/cma.j.issn.0529-567x.2018.04.009
5. Cancer Therapy Evaluation Program. Common Terminology Criteria for Adverse Events (CTCAE). Version 3.0, DCTD, NCI, NIH, DHHS. August 9 (2006). Available at: http://ctep.cancer.gov/protocol.Development/electronic_applications/docs/ctcaev3.pdf.
6. Bergmark K, Avall-Lundqvist E, Dickman PW, Henningsohn L, Steineck G. Vaginal Changes and Sexuality in Women With a History of Cervical Cancer. N Eng J Med (1999) 340(18):1383–9. doi: 10.1056/NEJM199905063401802
7. Bruner DW, Lanciano R, Keegan M, Corn B, Martin E, Hanks GE. Vaginal Stenosis and Sexual Function Following Intracavitary Radiation for the Treatment of Cervical and Endometrial Cancer. Int J Radiat Oncol Biol Phys (1993) 27(4):825–30. doi: 10.1016/0360-3016(93)90455-5
8. Kathrin K, Remi A, Jacob C, Haie-Meder C, Mahantshetty U, Segedin B, et al. Dose–effect Relationship and Risk Factors for Vaginal Stenosis After Definitive Radio(Chemo)Therapy With Image-Guided Brachytherapy for Locally Advanced Cervical Cancer in the EMBRACE Study. Radiother Oncol (2016) 118(1):160–6. doi: 10.1016/j.radonc.2015.12.025
9. Westerveld H, Pötter R, Berger D, Dankulchai P, Dörr W, Sora MC, et al. Vaginal Dose Point Reporting in Cervical Cancer Patients Treated With Combined 2D/3D External Beam Radiotherapy and 2D/3D Brachytherapy. Radiother Oncol (2013) 107(1):99–105. doi: 10.1016/j.radonc.2013.04.009
10. Wang J, Zhang K, Wang T, Liu Z, Wang RH, Zhang FQ, et al. Vaginal Dose of Radical Radiotherapy for Cervical Cancer in China: A Multicenter Study. BMC Cancer (2019) 19(1):1219. doi: 10.1186/s12885-019-6423-5
11. Zahra MA, Taylor A, Mould G, Coles C, Crawford R, Tan LT. Concurrent Weekly Cisplatin Chemotherapy and Radiotherapy in a Haemodialysis Patient With Locally Advanced Cervix Cancer. Clin Oncol (2008) 20(1):6–11. doi: 10.1016/j.clon.2007.10.007
12. Pötter R, Haie-Mede C, Limbergen E. Recommendations From Gynaecological (GYN) GEC ESTRO Working Group (II): Concepts and Terms in 3D Image-Based Treatment Planning in Cervix Cancer Brachytherapy-3D Dose Volume Parameters and Aspects of 3D Image-Based Anatomy, Radiation Physics, Radiobiology. Radiother Oncol (2006) 78(1):67–77. doi: 10.1016/j.radonc.2005.11.014
13. Potter R, Kirisits C, Erickson B, Haie-Meder C, Van Limbergen E, Lindegaard JC, et al. Prescribing, Recording, and Reporting Brachytherapy for Cancer of the Cervix. J ICRU (2013) 13(1–2). doi: 10.1093/jicru/ndw027
14. Chassagne D, Dutreix A, Almond P, Burgers J, Busch M, Joslin C eds. ICRU Report No. 38, Dose and Volume Specification for Reporting Intracavitary Therapy in Gynaecology. In: International Commissioning on Radiation Units and Measurements. Bethesda, MD: ICRU (1985).
15. Kenji Y, Noha J, Alina S, Lindegaard J, Segedin B, Mahantshetty U, et al. Assessment of Parametrial Response by Growth Pattern in Patients With International Federation of Gynecology and Obstetrics Stage IIB and IIIB Cervical Cancer: Analysis of Patients From a Prospective, Multicenter Trial (EMBRACE). Int J Radiat Oncol Biol Phys (2015) 93(4):788–96. doi: 10.1016/j.ijrobp.2015.08.007
16. Pötter R, Dimopoulos J, Georg P, Lang S, Waldhäusl C, Wachter-Gerstner N, et al. Clinical Impact of MRI Assisted Dose Volume Adaptation and Dose Escalation in Brachytherapy of Locally Advanced Cervix Cancer. Radiother Oncol (2007) 83(2):148–55. doi: 10.1016/j.radonc.2007.04.012
17. Pötter R, Georg P, Dimopoulos JC, Grimm M, Berger D, Nesvacil N, et al. Clinical Outcome of Protocol Based Image (MRI) Guided Adaptive Brachytherapy Combined With 3D Conformal Radiotherapy With or Without Chemotherapy in Patients With Locally Advanced Cervical Cancer. Radiother Oncol (2011) 100(1):116–23. doi: 10.1016/j.radonc.2011.07.012
18. Charra-Brunaud C, Harter V, Delannes M, Haie-Meder C, Quetin P, Kerr C, et al. Impact of 3D Image-Based PDR Brachytherapy on Outcome of Patients Treated for Cervix Carcinoma in France: Results of the French STIC Prospective Study. Radiother Oncol (2012) 103(3):305–13. doi: 10.1016/j.radonc.2012.04.007
19. Lindegaard JC, Fokdal LU, Nielsen SK, Juul-Christensen J, Tanderup K. MRI Guided Adaptive Radiotherapy in Locally Advanced Cervical Cancer From a Nordic Perspective. Acta Oncol (2013) 52(7):1510–9. doi: 10.3109/0284186X.2013.818253
20. EMBRACE, International MRI-Guided Brachytherapy in Cervical Cancer. Available at: https://www.Emb-racestudy.dk.
Keywords: radiotherapy, PIBS points system, ICRU rectum point, vaginal stenosis, cervical cancer
Citation: Wang J, Zhang K-s, Liu Z, Wang T, Wang R-h, Zhang F-q, Yu L, Wang Y-l, Wei L-c, Shi M, Li S, Liu B-g, Shi F, Su J, Yuan W, Zhang Qy and Zhang J (2022) Using New Vaginal Doses Evaluation System to Assess the Dose–Effect Relationship for Vaginal Stenosis After Definitive Radio(Chemo)Therapy for Cervical Cancer. Front. Oncol. 12:840144. doi: 10.3389/fonc.2022.840144
Received: 20 December 2021; Accepted: 21 March 2022;
Published: 19 April 2022.
Edited by:
Valerio Gallotta, Division of Gynecologic Oncology (IRCCS), ItalyReviewed by:
Alex Federico, Agostino Gemelli University Polyclinic (IRCCS), ItalyNicolò Bizzarri, Department of Women’s Health, Children’s Health and Public Health (IRCCS), Italy
Copyright © 2022 Wang, Zhang, Liu, Wang, Wang, Zhang, Yu, Wang, Wei, Shi, Li, Liu, Shi, Su, Yuan, Zhang and Zhang. 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: Zi Liu, bGl1em1haUAxNjMuY29t
†These authors share first authorship