- 1Department of Orthopaedic Surgery, Wuyi County First People's Hospital, Zhejiang, China
- 2Department of Orthopaedic Surgery, Changxing People's Hospital of Chongming District, Shanghai, China
- 3Spine Center, Xinhua Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China
Backgrounds: The incidence and characteristics of postoperative dysesthesia (POD) have not been reported for posterior vertebral column resection (PVCR) in the treatment of severe spinal kyphoscoliosis.
Objective: The objective of the study is to investigate the incidence and characteristics of POD in PVCR without anterior support applied in Yang's type A severe spinal kyphoscoliosis.
Material and methods: From August 2010 to December 2019, 167 patients diagnosed with Yang's type A severe spinal kyphoscoliosis who underwent PVCR without anterior support applied were retrospectively reviewed. All the patients were monitored using five modes of intraoperative multimodal neurophysiological monitoring. Neuromonitoring data, radiographic parameters, and neurological complications were reviewed and analyzed. The incidence and characteristics of POD were further summarized. POD was defined as dysesthetic pain or burning dysesthesia which could be caused by spinal cord kinking or dorsal root ganglion (DRG) injury but with no motor deficits.
Results: PVCR without anterior support was successfully conducted in all 167 patients. Intraoperative monitoring events occurred in five patients. One out of these five patients showed postoperative spinal cord injury (Frankel level C) but completely recovered within 9 months postoperation (Frankel level E). The number of levels and osteotomy space for vertebra resection were 1.28 and 3.6 cm, respectively. POD was confirmed in three patients (3/167, 1.8%), characterized as kyphosis with the apex vertebrae in T12 with the kyphotic Cobb angles of 100°, 115°, and 122°, respectively. The osteotomy space of vertebra resection in these three patients were 3.9, 3.8, and 4.2 cm, respectively. After the treatment by drug administration, they reported pain relief for 12–36 days. The pain gradually moved to the distal end of a proper DRG innervated region near the end.
Conclusions: In this study, the incidence rate of POD in Yang's type A severe spinal kyphoscoliosis patients who underwent PVCR without anterior support applied was 1.8% (3/167). Evoked potential monitoring could not detect the occurrence of POD. POD in Yang's type A severe spinal kyphoscoliosis after PVCR could be ascribed to spinal cord kinking and DRG injury.
Introduction
In 2016, Yang et al. first reported on a proposed classification system to guide surgical strategy in severe spinal deformity based on spinal cord function (1). The classification system defined type A as patients with normal spinal cord and evoked potentials and no neurological symptoms. Posterior vertebral column resection (PVCR) without anterior support applied was routinely used in our center to treat Yang's type A severe kyphoscoliosis, and similar clinical outcomes were observed as those in previous studies (2–4).
Dysesthesia is commonly reported in diabetes, Guillain–Barre syndrome, Lyme disease, multiple sclerosis, percutaneous endoscopic lumbar discectomy, and oral surgery patients (5–8). Dorsal root ganglion (DRG) injury–related postoperative dysesthesia (POD) has been reported in 0.9% patients with minimally invasive transforaminal lumbar interbody fusion, less than 5% patients with percutaneous endoscopic lumbar discectomy, and 7.7% patients with lumbar spine osteotomies (9–11). Although PVCR-related neurological deficit risks have been repeatedly emphasized (3, 12, 13), the POD has not been reported for PVCR in the treatment of severe spinal kyphoscoliosis (9–11), especially for Yang's type A kyphoscoliosis. Thus, the purpose of this study was to investigate the incidence and clinical characteristics of POD in our institution and present the outcomes of these patients.
Materials and methods
Patients
From August 2010 to December 2019, 167 consecutive Yang's type A kyphoscoliosis patients treated with PVCR without anterior support were retrospectively reviewed. The inclusion criteria included (1) scoliosis curve >80° or kyphosis >80°; (2) complete medical record data including radiography, intraoperative neurological monitoring, and spinal cord function; and (3) deformity diagnosed as type A according to Yang's spinal cord classification (normal spinal cord and evoked potentials, no neurological symptoms).
Radiographic evaluation
Standing anteroposterior and lateral radiographs of the entire spine were conducted before and after surgery. Left and right lateral flexion radiographs, whole-spine computed tomography (CT), and magnetic resonance imaging (MRI) were also performed before surgery. The radiographic variables evaluated in this study were the Cobb angle of major curve and its corrective rate, flexibility, and angle of kyphosis and its corrective rate.
Surgical technique
The indications for PVCR in the current patients included (1) severe scoliosis (>80°), (2) severe angular or tuberculous scoliosis (>80°), and (3) severity based on quantitative strategy on deformity angular ratio (14). The same senior doctor (JY) performed all surgeries under neurological monitoring. General anesthesia administration and endotracheal intubation were first performed, and then, a skin incision was made over the spinous processes to expose the intended fusion levels. Following pedicle screw placement using the free-hand technique, provisional instrumentation was then placed at the concave side using pedicle screws above and below the level of resection prior to the corpectomy and discectomy. Osteotomies were performed at the apex of kyphosis to increase the effectiveness of correction. Laminectomy was performed to expose the spinal cord. The corresponding rib heads were removed to expose the lateral walls of the pedicles. Then, compression force was applied to shorten the spinal column, and the osteotomy site was closed gradually. A wide laminectomy is often performed to extend one level above and below the PVCR level for better orientation and safety during compression. This approach can prevent shortening of the osteotomy space, which might result in kinking and bulking of the spinal cord. Under neuromonitoring, the osteotomy site must be completely closed so that the bony endplates could contact with each other. More attention should be given to prevent excessive shortening and interference with the intervertebral foramen, especially during osteotomy and correction. The osteotomy distance at the convex lateral border of apical vertebra was then measured and recorded.
POD
POD was defined as dysesthetic pain or burning dysesthesia which could be caused by spinal cord kinking or DRG injury but with no motor deficits. The visual analog scale/score (VAS) was assessed in patients with POD. Non-steroidal anti-inflammatory drugs (NSAIDs), central non-opioid analgesic agents (tramadol), and neuropathic pain medications (pregabalin) were selectively used to treat POD based on VAS.
Statistical analysis
Statistical analysis was performed using SPSS version 19.0 for Windows (SPSS Inc., Chicago, IL, USA). The normality of continuous data was determined using the Shapiro–Wilk test. Analysis of variance was performed to compare parameters among three vertebral column resection (VCR) locations. A P-value less than 0.05 was considered significant. Data were presented as mean ± standard deviation.
Results
Patients
A total of 167 patients (99 males and 68 females), who met the inclusion criteria and received PVCR without anterior support, were included in the analysis. All patients were successfully followed up, with an average duration of 5 years. The mean age at the time of surgery was 25.7 years. The etiological diagnoses included 96 idiopathic cases, 51 congenital cases, and 20 neurofibromatosis cases. Table 1 summarizes the preoperative demographics.
Content of PVCR surgery
No structural anterior support was applied in all 167 PVCR patients (Figure 1). The number of levels for vertebra resection was 1.3 (range: 1–2 levels), with an average osteotomy distance of 3.6 cm (range: 2.7–5.0 cm) (Figure 2). There were 145 one-level VCRs and 22 two-level contiguous VCRs. According to the location of VCR, there were 22 high thoracic cases (T1–T4) with 3.6 ± 0.9 cm osteotomy distance, 112 middle thoracic cases (T5–T8) with 3.8 ± 1.2 cm osteotomy distance, 33 lower thoracic cases (T9–T12) with 3.9 ± 0.7 cm osteotomy distance. The average intraoperative blood loss and surgery time were 2,831.2 ± 965.4 ml and 479.2 ± 68.3 min, respectively. There were no significant statistical differences in the number of levels for vertebra resection, osteotomy distance, blood loss, and surgery time among the three VCR locations (P-value > 0.05).
Figure 1. A 21-year-old male with congenital scoliosis with Yang's type A kyphoscoliosis treated with PVCR without anterior support. The patient was diagnosed with congenital scoliosis with the apex located at T12 preoperatively. (A–C) The preoperative X-ray and CT of the patient; (D,E) postoperative X-ray of the patient; (F,G) state of the spinal cord before and after the closure of the osteotomy gap during the operation. It can be clearly seen that there are bulking/kinking in the spinal cord after the closure of the osteotomy gap.
Figure 2. Neuroelectrophysiological monitoring outcome of the patient with postoperative dysesthesia. (A) Intraoperative SEP monitoring for bilateral lower extremities; (B) intraoperative left-MEP monitoring; (C) intraoperative right-MEP monitoring; (D) postoperative SEP monitoring 3 day later. There were no intraoperative evoked potential changes (A–D).
Scoliosis and kyphosis correction rates
The preoperative coronal plane major curve (with a mean flexibility of 17.6%) was corrected from 126.4 ± 28.3° to 55.0 ± 16.5°, showing a mean of 58.2% correction. Preoperative thoracic kyphosis was corrected from 128.4 ± 24.3° to 53.3 ± 16.5°, showing a mean of 62.1% correction. There were no significant statistical differences in all these radiographic parameters among the three VCR locations (P-value > 0.05) (Table 2 and Figure 1).
Intraoperative and postoperative neural event
In this study, intraoperative monitoring events occurred in five patients (5/167). One out of the five patients showed transient spinal cord injury (a decrease from Frankel level E to Frankel level C) after the operation, whose cord function recovered to Frankel level E after 9 months of follow-up.
POD
Immediately after operation, three cases of POD were confirmed (3/167, 1.8%), which were characterized as kyphosis with the apex vertebrae in T12, with kyphotic Cobb angles of 100°, 115°, and 122°, respectively. The osteotomy space of vertebra resection in these three patients were 3.9, 3.8, and 4.2 cm, respectively. After the patients woke up from anesthesia, POD was characterized by persistent lateral or bilateral lower limb hyperalgesia without motor impairment. The degree of pain was intense, with a VAS score of more than seven points. Slight skin touch can induce or aggravate the degree of pain. After treatment by drug administration (NSAIDs, tramadol, or pregabalin), all three patients reported pain relief with a duration of 12–36 days (Table 3). We routinely adopted five evoked potential monitoring modes, and no intraoperative evoked potential changes occurred in the three patients. Somatosensory evoked potential (SEP) was re-monitored for these patients 3 days after operation, and the results were also normal, indicating that the occurrence of POD could not be detected by neuroelectrophysiological monitoring (Figure 2). In the process of rehabilitation, the POD region of three patients gradually moved to the distal lower extremity, which finally disappeared. None complained of any symptoms of POD at the final follow-up. Thus, there was a favorable prognosis for POD.
Discussion
This study found that the incidence rate of POD in Yang's type A severe spinal kyphoscoliosis patients who underwent PVCR without anterior support applied was 1.8% (3/167). Although the occurrence of POD is not associated with motor dysfunction, its impact on the patients’ quality of life is quite significant and should be taken seriously by surgeons.
It has been reported that intraoperative DRG injury was responsible for POD, and it was caused by intraoperative pulling of nerve tissues, mechanical stimulation, bipolar electrical coagulation disturbance, and accumulation of closed spinal cord in osteotomy (15, 16). In our series, 1.8% (3 out of 167) patients underwent POD, while all of them were characterized as severe kyphosis with apex vertebrae in T12. The osteotomy space of vertebra resection in the three patients were 3.9, 3.8, and 4.2 cm. Dural kinking and bucking due to the bone-on-bone technique after PVCR were obvious in these three patients. Thus, we speculated that, in addition to related causes reported in literature, the magnitude of kyphosis angle, kyphosis apex position, osteotomy closure distance, and spinal cord kinking and bulking might be the main risk factors for POD in PVCR surgery in the treatment of Yang's type A severe spinal kyphoscoliosis.
The main manifestation of POD in this study was hyperalgesia in the lateral or bilateral lower limbs, with severe pain (VAS scores > 7), which is more serious than that of previous investigations (9–11). The pain can be induced or aggravated by slight skin touch, similar to the characteristics in published literatures (17, 18). However, all the three patients showed no lower limb motor dysfunction, and the manifestations of their nervous system were mainly hyperalgesia. The deep reflex, deep sensation, and pathology of lower limbs were negative. One patient initially presented with impaired left lower limb movement caused by severe pain (VAS 9), and the pain relief by analgesic drugs (celecoxib and pregabalin) proved that motor function was barrier-free.
There is no unified standard for the diagnosis of POD, which mainly depends on clinical symptoms and clinical experiences (19, 20). In our series, two patients (2/3) showed dysesthesia in the lateral or bilateral limbs, one patient showed anaphylaxis in the L1 and L2 nerve root innervation regions, and another showed anaphylaxis in the L3 and L4 regions without motor dysfunction. Another patient (1/3) showed hyperalgesia in areas L4 and L5 combined with left lower limb extensor knee movement disorder. After the use of analgesics, motor dysfunction was confirmed to be caused by pain. POD is easily misdiagnosed as nerve root injury, making it essential to exercise caution when diagnosing POD, especially when the nerve injury and POD occur at the same time (21, 22). Our experience is that isolated POD cannot be diagnosed until motor dysfunction is completely excluded. We routinely adopted five evoked potential monitoring modes throughout the whole operation, and no intraoperative evoked potential changes occurred in the three patients. We re-monitored the three patients with SEP 3 days after operation, and the results were still normal, indicating that the occurrence of POD could not be recognized by neuroelectrophysiological monitoring. In all the three patients, the pain gradually moved to the distal end of a proper DRG innervated region during the recovering period.
The treatment of POD aimed to relieve pain and shorten the recovery period. NSAIDS, opioids, glucocorticoid drugs, antidepressants, and neuralgia drugs have been reported to show certain clinical effects (9–11). All the three patients presented initial VAS scores of more than seven, and favorable outcomes were obtained with the treatment of pregabalin combined with celecoxib. The pain of three patients was relieved by about 50% the first day after medication treatment, and the lower limb pain of three patients was completely relieved with medication treatment for 12–36 days. Previous literature recommended antidepressant drugs for the treatment of POD (23). In addition, intervertebral foramen block, with analgesic drug and steroid administration, also has been shown effective (24, 25). However, we did not use this treatment strategy in this study because the combination of celecoxib and pregabalin already achieved significant pain relief.
Although it has been reported that the prognosis of POD after routine spinal surgery is generally satisfactory, the prognosis of POD after PVCR in severe spinal kyphoscoliosis has not been well described. According to the condition of our three patients, the paresthesia could be completely relieved after 12–36 days of treatment, which is similar to previous studies (23–25). In addition, we believe that the gradual migration of the anaphylaxis plane to the distal extremity during the treatment process could be a significant feature indicating the potential for a cure for POD.
This study has several limitations. First, this was a retrospective study performed in a single center by a single surgeon using one technique. Therefore, it is challenging to generalize the incidence and prognosis of POD to other patients and other techniques of PVCR. In addition, there were insufficient patients (three patients) diagnosed with POD in the current study, so the symptom, pathogenic mechanism, and treatment strategy should be further investigated and discussed.
In conclusion, in the current study, the incidence of POD in Yang's type A severe spinal kyphoscoliosis patients who underwent PVCR without anterior support applied was 1.8% (3/167). Evoked potential monitoring could not detect the occurrence of POD. POD in Yang's type A severe spinal kyphoscoliosis after PVCR could be ascribed to spinal cord kinking and DRG injury.
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 Ethics Committee of the Xinhua Hospital Affiliated to Shanghai Jiaotong University School of Medicine reviewed and approved the studies involving human participants. The patients/participants provided their written informed consent to participate in this study.
Author contributions
HW and JD designed the study, analyzed and interpreted the data, and wrote and edited the manuscript. JY designed the study, analyzed and interpreted the data, and wrote and edited the manuscript. WS designed the study, analyzed and interpreted the data, and wrote and edited the manuscript. HW and JD contributed equally and were considered co-first authors. JY and WS contributed equally and were considered co-corresponding authors. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Medical-engineering Cross Fund of Rising Star Plan in Shanghai Jiao Tong University (YG2021QN42), the National Key Research and Development Program (2018YFC0116500), and the Key Project of Transformational Medicine Cross-Research Fund of Shanghai Jiaotong University (ZH2018ZDB04).
Acknowledgment
We acknowledge all the participants enrolled in the current study.
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.
Abbreviations
POD, postoperative dysesthesia; PVCR, posterior vertebral column resection; SEP, somatosensory evoked potential; MEP, motor evoked potential; DNEP, descending neurogenic evoked potential; DRG, dorsal root ganglion.
References
1. Yang JL, Huang ZF, Yin JQ, Deng YL, Xie XB, Li FB, et al. A proposed classification system for guiding surgical strategy in cases of severe spinal deformity based on spinal cord function. Eur Spine J. (2016) 25:1821–9. doi: 10.1007/s00586-015-4367-2
2. Suk SI, Chung ER, Kim JH, Kim SS, Lee JS, Choi WK. Posterior vertebral column resection for severe rigid scoliosis. Spine (Phila Pa 1976). (2005) 30:1682–7. doi: 10.1097/01.brs.0000170590.21071.c1
3. Lenke LG, O'leary PT, Bridwell KH, Sides BA, Koester LA, Blanke KM. Posterior vertebral column resection for severe pediatric deformity: minimum two-year follow-up of thirty-five consecutive patients. Spine (Phila Pa 1976). (2009) 34:2213–21. doi: 10.1097/BRS.0b013e3181b53cba
4. Lee BH, Hyun SJ, Han S, Kim KJ, Jahng TA, Kim YJ, et al. Surgical and radiological outcomes after posterior vertebral column resection according to the surgeon's experience. Medicine (Baltimore). (2018) 97:e11660. doi: 10.1097/MD.0000000000011660
5. Baba K, Aridome K, Haketa T, Kino K, Ohyama T. Sensory perceptive and discriminative abilities of patients with occlusal dysesthesia. Nihon Hotetsu Shika Gakkai Zasshi. (2005) 49:599–607. doi: 10.2186/jjps.49.599
6. Toyofuku A, Kikuta T. Treatment of phantom bite syndrome with milnacipran—a case series. Neuropsychiatr Dis Treat. (2006) 2:387–90. doi: 10.2147/nedt.2006.2.3.387
7. Reeves JL, Merrill RL. Diagnostic and treatment challenges in occlusal dysesthesia. J Calif Dent Assoc. (2007) 35:198–207. doi: 10.1080/19424396.2007.12221218
8. Hara ES, Matsuka Y, Minakuchi H, Clark GT, Kuboki T. Occlusal dysesthesia: a qualitative systematic review of the epidemiology, aetiology and management. J Oral Rehabil. (2012) 39:630–8. doi: 10.1111/j.1365-2842.2012.02300.x
9. Yeung AT, Tsou PM. Posterolateral endoscopic excision for lumbar disc herniation: surgical technique, outcome, and complications in 307 consecutive cases. Spine (Phila Pa 1976). (2002) 27:722–31. doi: 10.1097/00007632-200204010-00009
10. Wang H, Zhou Y, Zhang Z. Postoperative dysesthesia in minimally invasive transforaminal lumbar interbody fusion: a report of five cases. Eur Spine J. (2016) 25:1595–600. doi: 10.1007/s00586-015-4365-4
11. Zhang Z, Wang H, Zheng W. Postoperative dysesthesia in lumbar three-column resection osteotomies. Eur Spine J. (2016) 25:2622–8. doi: 10.1007/s00586-014-3574-6
12. Wang XB, Lenke LG, Thuet E, Blanke K, Koester LA, Roth M. Deformity angular ratio describes the severity of spinal deformity and predicts the risk of neurologic deficit in posterior vertebral column resection surgery. Spine (Phila Pa 1976). (2016) 41:1447–55. doi: 10.1097/BRS.0000000000001547
13. Huang ZF, Chen L, Yang JF, Deng YL, Sui WY, Yang JL. Multimodality intraoperative neuromonitoring in severe thoracic deformity posterior vertebral column resection correction. World Neurosurg. (2019) 127:e416–26. doi: 10.1016/j.wneu.2019.03.140
14. Huang Z, Sui W, Huang H, Deng Y, Li J, Liu L, et al. Quantitative determining of pre-operative osteotomy plan for severe spinal deformity: an analysis of 131 consecutive Yang's A type cases from single center. Eur Spine J. (2021) 30:3200–8. doi: 10.1007/s00586-021-06840-1
15. Moon HS, Kim YD, Song BH, Cha YD, Song JH, Lee MH. Position of dorsal root ganglia in the lumbosacral region in patients with radiculopathy. Korean J Anesthesiol. (2010) 59:398–402. doi: 10.4097/kjae.2010.59.6.398
16. Roychoudhury S, Nagori SA, Roychoudhury A. Neurosensory disturbance after bilateral sagittal split osteotomy: a retrospective study. J Oral Biol Craniofac Res. (2015) 5:65–8. doi: 10.1016/j.jobcr.2015.04.006
17. Hoss D, Segal S. Scalp dysesthesia. Arch Dermatol. (1998) 134:327–30. doi: 10.1001/archderm.134.3.327
18. Thoppay JR, De Rossi SS, Ciarrocca KN. Burning mouth syndrome. Dent Clin North Am. (2013) 57:497–512. doi: 10.1016/j.cden.2013.04.010
19. Bruckmoser E, Bulla M, Alacamlioglu Y, Steiner I, Watzke IM. Factors influencing neurosensory disturbance after bilateral sagittal split osteotomy: retrospective analysis after 6 and 12 months. Oral Surg Oral Med Oral Pathol Oral Radiol. (2013) 115:473–82. doi: 10.1016/j.oooo.2012.08.454
20. Kuroyanagi N, Miyachi H, Ochiai S, Kamiya N, Kanazawa T, Nagao T, et al. Prediction of neurosensory alterations after sagittal split ramus osteotomy. Int J Oral Maxillofac Surg. (2013) 42:814–22. doi: 10.1016/j.ijom.2012.11.016
21. Auerbach JD, Lenke LG, Bridwell KH, Sehn JK, Milby AH, Bumpass D, et al. Major complications and comparison between 3-column osteotomy techniques in 105 consecutive spinal deformity procedures. Spine (Phila Pa 1976). (2012) 37:1198–210. doi: 10.1097/BRS.0b013e31824fffde
22. Norton RP, Bianco K 2nd, Lafage V, Schwab FJ. Complications and intercenter variability of three-column resection osteotomies for spinal deformity surgery: a retrospective review of 423 patients. Evid Based Spine Care J. (2013) 4:157–9. doi: 10.1055/s-0033-1357364
23. Kalia LV, O'connor PW. Severity of chronic pain and its relationship to quality of life in multiple sclerosis. Mult Scler. (2005) 11:322–7. doi: 10.1191/1352458505ms1168oa
24. Stromberg L, Toohey JS, Neidre A, Ramsey M, Brantigan JW. Complications and surgical considerations in posterior lumbar interbody fusion with carbon fiber interbody cages and Steffee pedicle screws and plates. Orthopedics. (2003) 26:1039–43. doi: 10.3928/0147-7447-20031001-13
Keywords: postoperative dysesthesia, posterior vertebral column resection, severe spinal deformity, kyphoscoliosis, neurological monitoring
Citation: Wu H, Ding J, Yang J and Sui W (2023) Postoperative dysesthesia after PVCR without anterior support applied in Yang's type A severe spinal kyphoscoliosis. Front. Surg. 10:1222520. doi: 10.3389/fsurg.2023.1222520
Received: 14 May 2023; Accepted: 3 July 2023;
Published: 19 July 2023.
Edited by:
Yu Song, Huazhong University of Science and Technology, ChinaReviewed by:
Rami Rachkidi, Saint Joseph University, LebanonXueshi Li, First Affiliated Hospital of Jinan University, China
Shengru Wang, Peking Union Medical College Hospital (CAMS), China
Shifeng Wen, Guangzhou First People's Hospital, China
© 2023 Wu, Ding, Yang and Sui. 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: Junlin Yang eWFuZ2p1bmxpbkB4aW5odWFtZWQuY29tLmNu Wenyuan Sui c3Vpd2VueXVhbkBob3RtYWlsLmNvbQ==
†These authors share first authorship