- 1School of Microeletronics, Xidian University, Xi’an, China
- 2Guangdong Institute of Intelligence Science and Technology, Zhuhai, China
Transcranial ultrasound stimulation (TUS) is a young neuromodulation technology, which uses ultrasound to achieve non-invasive stimulation or inhibition of deep intracranial brain regions, with the advantages of non-invasive, deep penetration, and high resolution. It is widely considered to be one of the most promising techniques for probing brain function and treating brain diseases. In preclinical studies, developing miniaturized transducers to facilitate neuromodulation in freely moving small animals is critical for understanding the mechanism and exploring potential applications. In this article, a miniaturized transducer with a half-concave structure is proposed. Based on the finite element simulation models established by PZFlex software, several ultrasound transducers with different concave curvatures were designed and analyzed. Based on the simulation results, half-concave focused ultrasonic transducers with curvature radii of 5 mm and 7.5 mm were fabricated. Additionally, the emission acoustic fields of the ultrasonic transducers with different structures were characterized at their thickness resonance frequencies of 1 MHz using a multifunctional ultrasonic test platform built in the laboratory. To verify the practical ability for neuromodulation, different ultrasound transducers were used to induce muscle activity in mice. As a result, the stimulation success rates were (32 ± 10)%, (65 ± 8)%, and (84 ± 7)%, respectively, by using flat, #7, and #5 transducers, which shows the simulation and experimental results have a good agreement and that the miniaturized half-concave transducer could effectively converge the acoustic energy and achieve precise and effective ultrasonic neuromodulation.
Introduction
Neuromodulation technology is fundamental to applications such as the probe of the nervous system and treating neurological disorders (Miniussi and Vallar, 2011; Kar and Sarkar, 2016). The available neuromodulation techniques such as deep brain stimulation, transcranial magnetic stimulation, and the newly developed optogenetics for preclinical use are either invasive or lack satisfying spatiotemporal resolution. Alternatively, transcranial ultrasound stimulation (TUS), enabling non-invasive neuromodulation with the advantages of non-invasive (Hariz, 2013; Panczykowski et al., 2014; Baek et al., 2017; Li et al., 2018), high spatial resolution, and high penetration depth, has been considered a next generation technology for probing human brain functions and treating dysfunctions (Yamaoka et al., 2011; Mead et al., 2017; Guo et al., 2021). However, since it is still in its infant stages, preclinical studies in small animals are required to figure out the fundamental mechanisms and optimal parameters for treating brain dysfunction models.
In preclinical models, it has been shown that TUS can be used to treat various brain diseases including Parkinson’s diseases, Alzheimer’s diseases, and addiction (Yang Y. et al., 2021). On the other hand, sonogenetics facilitates the non-invasive dissection of neural circuits in mouse brains. In small animal studies, anesthesia level is hard to control, which could bring additional confound to the studies as several results demonstrated that the effects of ultrasound stimulation are strongly dependent on the level of anesthesia, which can affect the reliability of experimental results (King et al., 2013; Lee et al., 2016). The use of more miniature transducers mounting on the skull is a potential solution to avoid the operation of anesthesia and immobilization. However, a small size transducer (sub-MHz) is always accompanied by low acoustic energy and spatial resolution (Zhang et al., 2020), so appropriate stimulation conditions and methods for converging acoustic energy become the critical factor of neuromodulation.
With the development of neural regulation technology, commercial planar ultrasonic transducers gradually cannot meet the precise neural regulation. In recent years, some studies have used focused ultrasonic transducers to improve accuracy (Yoo et al., 2011; Deffieux et al., 2013; Kim et al., 2014). In this study, we proposed a transducer that both has a small wearable size and a half-concave structure to converge acoustic energy, as shown in Figure 1. The half-concave structure was made by the dimple grinder. The finite element simulation software PZFlex was used to design and analyze the single element ultrasonic transducer with different dimpling curvatures. Muscle activity in mice was induced by ultrasonic transducers with different structures. The results indicate that the half-concave ultrasonic transducer can achieve more effective convergence of sound energy and make more acoustic energy transmit through the skull when compared to the plane transducer, which leads to a higher success rate.
Methods and Results
The characterization of material properties is the first and critical step in the design of devices because the piezoelectric dielectric properties of materials will directly affect the final performance of the device, which is why scholars have studied it extensively over the past years (Kuzenko, 2021; Yang L. et al., 2021). Several kinds of literature have pointed out that the energy transmission efficiency is higher when the ultrasonic stimulation frequency is lower than 1 MHz (Hayner and Hynynen, 2001; White et al., 2006), so it is suitable to study the acoustic characteristics of ultrasonic transducers with different structures at 400 kHz. Due to the impedance matching, the diameter of the piezoelectric element of a conventional 400 kHz transducer is much higher than 85 mm. Inversely, the diameter of the half-concave transducer is just 5 mm, and the piezoelectric element is 2 mm in thickness. Therefore, a miniaturized device is designed through half-concave structural optimization. The electrical impedance and the phase plots of the piezoelectric element were measured using the WK6500B 1J65120B impedance analyzer (Wayne Kerr Electronics, United Kingdom), as shown in Figure 2.
As shown in Figure 2, the measurement plot of the planar piezoelectric indicates that there are a lot of resonance peaks, and in order to further determine the excitation frequency for subsequent experiments, the spectrum of different curvatures was simulated with the help of finite element simulation. In the simulation results of the piezoelectric element in Figure 3, it is seen that there is a resonance frequency near 1 MHz, but with the change of curvature, the resonant peak has a certain deviation, but at 400 kHz, resonance coincidence is better.
The finite element simulation has been proved to be an efficient simulation method by many previous studies (Hou et al., 2021; Li et al., 2021b; Lin et al., 2021). In this study, the finite element simulation model was established to design the half-concave structure and simulate the performance of the half-concave ultrasonic transducer with the software PZFlex and Comsol. All materials used in this simulation are listed in Table 1.
The simulation parameters are as follows. The thickness of the backing material is set to 6 mm in case of backpropagation. Combined with the limitation of the actual manufacturing process, the dimpling curvature of the half-concave ultrasonic transducer was increased from 5 to 7.5 mm at equal intervals of 0.5 mm. Besides, an insulating epoxy that surrounds the piezoelectric material is 0.4 mm in width to make the condition of simulation consistent with the experiment. During the simulation, the ultrasonic transducer is excited by one cycle of sinusoidal signal with a frequency of 400 kHz and a driving voltage of 1 V peak-to-peak. Figures 4A,B show the structure diagram of the plane ultrasonic transducer and the half-concave ultrasonic transducer in simulation, respectively. The blue part represents the epoxy that acted as the backing layer and the insulating, and the purple part represents the piezoelectric material, and a reflecting plate is placed at the end of the water body, which is considered as the focal area.
Figure 4. Simulation structure: (A) The half-concave ultrasonic transducer and (B). The plane ultrasonic transducer.
Under the above parameters, a normalized emission sound field is obtained through MATLAB processing, as shown in Figure 5. Due to the excitation frequency, the wavelength of the emitted sound wave is similar to the size of the transducer, so the transducer can be regarded as a point source, and the acoustic wavefront is similar to the spherical wave (Crew et al., 1900; Simon et al., 2015). The maximum acoustic pressure data obtained from the emission acoustic field are listed in Table 2. It can be seen that the difference in the maximum acoustic pressure between the transducers with different radii of dimpling curvature is slight, but the maximum acoustic pressure of the half-concave ultrasonic transducer is slightly larger than that of the plane transducer, and the area with higher energy in the acoustic field is also larger than that of the plane transducer, which demonstrates that half-concave transducer can effectively converge the sound energy and is more beneficial to the neuromodulation of small animals.
Based on the simulation results, subsequent experiments were carried out. The manufacturing process of the half-concave structure and the structure diagram of the dimpled ultrasonic transducer are shown in Figures 6A,B. As the piezoelectric material, the piezoelectric ceramic pic255 with double-side silver electrodes was prepared in advance. Then one side of the pic255 piezoelectric ceramics was dimpled to a half-concave structure by the dimple grinder (model656, Gatan), and a 0.9 mm diameter lead was attached to the backside of the dimpled using conductive epoxy e-solder 3022 (VonRoll Isola, New Haven, CT). After that, the piezoelectric material was placed at the center of the brass housing, and then the housing was filled with insulating epoxy (EPO-TEK 301) to act as both a backing and an insulating layer. After epoxy curing, a layer of gold (∼200 μm) was sputtered on the dimpled side to act as the front electrode, and subsequently, a thin layer of parylene (∼1.5 μm) was vapor-deposited by a PDS 2010 Labcoator (Specialty Coating Systems, Indianapolis, IN) to protect the ultrasonic transducer. Finally, an SMA connector was attached to the ultrasonic transducer. The dimpled piezoelectric ceramic pic255 and the packaged ultrasonic transducer are shown in Figures 6C,D, respectively. Besides, a plane transducer was fabricated as the control group.
Figure 6. (A) Manufacture process of dimpling, (B) the structure diagram of the half-concave ultrasonic transducer, (C) dimpled pic255 piezoelectric ceramic, and (D) the packaged ultrasonic transducers.
The laboratory multi-function ultrasonic testing system based on LabView software was a high precision test platform (Li et al., 2021a,c). Based on the platform, the emission acoustic fields of ultrasonic transducers with different structures were measured at 400 kHz. In deionized water, the transducer was placed opposite a needle hydrophone (NH0200, Precision Acoustics, Dorchester, United Kingdom). Additionally, the acoustic pressure was measured in the range of 4 mm2 × 5 mm2 in the axial plane of the transducer. Then it was used to get the normalized emission acoustic field with MATLAB software, as shown in Figure 7. It can be seen that the emission acoustic fields of different structures of ultrasonic transducers are similar, and the difference in the maximum acoustic pressure is slight in the experiment, which is basically consistent with the simulation results.
However, in biological applications, it is not comprehensive to consider only the emission acoustic field of ultrasonic transducers in water. In order to further prove the ability of the half-concave ultrasonic transducer to carry out biological experiments, a skull with a thickness of 0.5 mm was placed in front of the transducer, and the finite element simulation software PZFlex was also used to study the acoustic field distribution after transmitting through the skull. All models were set in water. The simulation model of the skull with a radius of curvature of 9 mm appressed to the surface of the transducer was established, and the relative position of the skull and the transducer is shown in Figure 8D, in which the right inset shows the skull appressed to piezoelectric material of different transducers. MATLAB software was used to process the acoustic field data after transmitting it through the skull. The normalized acoustic field was obtained as shown in Figures 8A,B. The red dash line shows the –6 dB line in the acoustic pressure distribution and the normalized lateral acoustic field distribution of different transducers at maximum acoustic pressure point is shown in Figure 8C, which shows the area with higher acoustic field energy can fall in the brain area that needs neuromodulation as far as possible. The maximum acoustic pressure is shown in Table 3.
Figure 8. (A) Normalized acoustic field after transmitting through the skull by simulation. (B) Normalized lateral acoustic field distribution of different transducers by simulation at a maximum acoustic pressure point. (C) The normalized acoustic pressure distribution of the focal area in the X-Y plane. (D) The relative position of the skull and different transducers.
It can be seen that the dimpling structure allows more energy to transmit through the skull, and the maximum acoustic pressure of the half-concave ultrasonic transducer is increased by about 20% compared with the plane transducer. The dimpling depth is greater, and the maximum acoustic pressure is greater. The greater dimpling depth achieves the greater focus effect which shows the higher resolution and near focus depth. So, the effective focus and short transmission distance can achieve higher sound energy. In summary, the half-concave structure can converge acoustic energy, and the curvature of the half-concave ultrasonic transducer is more suitable for the head of small animals. Different ultrasonic transducers were used to induce muscle activity in mice to verify the biological application value of the half-concave structure.
To further demonstrate the neuromodulation effects by using the miniaturized transducers, we stimulated the mouse brain with these transducers and monitored the motor responses determined by electromyography (EMG). The EMG recording procedures were similar to previous work (Qiu et al., 2021). Briefly, mice underwent anesthetic induction with 3% isoflurane in 100% oxygen and were maintained on 0.6% isoflurane for TUS. Eye ointments were applied to both eyes. Two 32-gauge enamel-coated copper electrodes were inserted into the triceps muscles of both forelimbs. After the attachment of all leads, the US transducer was fixed to a three-axis positioning system. Ultrasound gel was used to couple the transducer to the head of the mouse. For all experiments, the transducer was positioned approximately 1 mm from the surface of the animal’s head, targeting the primary motor cortex as shown in Figure 9A. A heat lamp distant from the experiment table was used to keep the animal warm. All the stimulation and recordings were finished within 5 min. After each mouse experiment, the EMG leads were removed, and the mouse was transferred to an induction chamber for recovery and then returned to its cage. The EMG signals were amplified with a gain of 1,000, and the bandpass was filtered between 10 Hz and 1 kHz with a preamplifier (World Precision Instruments, Sarasota, FL, United States). Data were acquired at a 2 kHz sampling rate (Lab-Jack U3, LabJack, Lakewood, CO, United States). The sonication parameters (PRF 1 kHz, 40% duty cycle, 200 ms) were controlled by a computer running software written in MATLAB (Mathworks, Natick, MA, United States). As shown in Figure 9B, The miniaturized transducers can induce EMG responses in mice experiments. The focused transducers resulted in a higher success rate compared to the flat transducer. The EMG signals were post-processed and analyzed using additional software written in MATLAB. For each trial, the DC drift was removed with a 10 Hz Butterworth IIR filter, followed by a notch filter for line frequency (60 Hz) removal. The signal was then full-wave rectified and smoothed with a 15-point moving average filter. A signal was considered representative of a muscle contraction if two conditions were met: First, the smoothed EMG signal exceeded a contraction threshold, defined as 6 SDs of the signal 100 ms before sonication plus the mean of the signal during this period. Second, the smoothed EMG signal exceeded the contraction threshold within a temporal latency of less than 200 ms. The representative responses was shown in Figure 9C. The response was then calculated for each set of ultrasound signals as a success rate defined as the number of responses divided by the total number of sonications. As shown in Figure 9D, the stimulation success rates were (32 ± 10)%, (65 ± 8)%, and (84 ± 7)%, respectively, by using flat, #7, and #5 transducers. The response latencies (Figure 9E) showed a decreasing trend, but there is no significant difference. These results were expected as these transducers have different ultrasound intensities transmitted through the skull, as shown in Figure 9.
Figure 9. (A) Schematic diagram of in vivo mouse EMG recordings using ultrasound neuromodulation to the motor cortex. (B) Example sync pulses (orange) and EMG responses (blue) resulting from sonication with sham, flat transducer, focused transducer #7.5, and #5 at 1.6 MPa. (C) Representative EMG signals induced by a single trial, zooming from (B). (D) The success rate of ultrasound evoked responses and the response latencies (E) sonicated by sham, flat, focused transducer #7.5, and #5 (N = 6, 20 trials/animal, *p < 0.05, **p < 0.01, one-way ANOVA with post-hoc Tukey test. All statistically significant differences are shown).
Conclusion
In this study, a new structure of the miniaturized ultrasonic transducer was proposed for neuromodulation, and one side of the piezoelectric material was dimpled to form the half-concave structure. Both simulated and experimental emission acoustic fields show that the half-concave structure converges the acoustic energy so that the higher energy areas of the sound field fall in brain regions where neuromodulation is needed, which would lead to more effective neuromodulation. Besides, the size of the half-concave ultrasonic transducer is smaller than most transducers in neuromodulation-related articles, which makes the biological experiment more flexible. To verify the validity of the design, different transducers were used to induce muscle activity in mice, and the result shows the half-concave ultrasonic transducer led to a higher success rate. Using the fabricated transducer with flat, #7, and #5, the stimulation success rates were (32 ± 10)%, (65 ± 8)%, and (84 ± 7)%, respectively. All results indicate the research on the half-concave structure could promote the development of neuromodulation.
Data Availability Statement
The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
All animal experiments were approved by the Animal Subjects Ethics 7 Sub-Committee (ASESC) of the Hong Kong Polytechnic University. Animal use and 8 care were performed following the guidelines of the Department of Health—Animals 9 (Control of Experiments) of the Hong Kong SAR government.
Author Contributions
CH, YW, and XS: drafting and refining the manuscript. ZL and CZ: conducting of the experiment and simulation. CF, ZQ, and YY: critical reading of the manuscript. All authors have read and approved the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (No. 61974110), the Key Projects of National Defense Basic Research Program of National Defense Science and Industry Administration (No. JCKY2019210B003), the Shaanxi Provincial Association of Science and Technology Young Talents Support Project (No. 20190105), the Natural Science Foundation of Shaanxi Province (No. 2021JM-253), and Special Fund for Science and Technology Innovation Strategy of Guangdong Province (No. 2021B090905000).
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
Baek, H., Pahk, K., and Kim, H. (2017). A review of low-intensity focused ultrasound for neuromodulation. Biomed. Eng. Lett. 7, 135–142. doi: 10.1007/s13534-016-0007-y
Crew, H., Huygens, C., Young, T., Fresnel, A., and Arago, F. (1900). The Wave Theory Of Light. New York: American Book Company.
Deffieux, T., Wattiez, N., Tanter, M., Pouget, P., Aubry, J., and Younan, Y. (2013). Low intensity focused ultrasound modulates monkey visuomotor behavior. J. Ther. Ultra. 23, 2430–2433. doi: 10.1016/j.cub.2013.10.029
Guo, J., Song, X., Chen, X., Xu, M., and Ming, D. (2021). Mathematical model of ultrasound attenuation with skull thickness for transcranial-focused ultrasound. Front. Neurosci 15:778616. doi: 10.3389/fnins.2021.778616
Hariz, M. (2013). Focused ultrasound thalamotomy improves essential tremor. Mov. Dis. 28, 1803–1803. doi: 10.1002/mds.25695
Hayner, M., and Hynynen, K. (2001). Numerical analysis of ultrasonic transmission and absorption of oblique plane waves through the human skull. J. Acoust. Soc. Am. 110, 3319–3330. doi: 10.1121/1.1410964
Hou, C., Fei, C., Li, Z., Zhang, S., Man, J., Chen, D., et al. (2021). Optimized backing layers design for high frequency broad bandwidth ultrasonic transducer. IEEE Trans. Biomed. Eng. 69, 475–481. doi: 10.1109/tbme.2021.3098567
Kar, S., and Sarkar, S. (2016). Neuro-stimulation techniques for the management of anxiety disorders: an update. Clin. Psychopharmacol. Neurosci. 14, 330–337. doi: 10.9758/cpn.2016.14.4.330
Kim, H., Chiu, A., Lee, S., Fischer, K., and Yoo, S. (2014). Focused ultrasound-mediated non-invasive brain stimulation: examination of sonication parameters. Brain Stimulat. 7, 748–756. doi: 10.1016/j.brs.2014.06.011
King, R., Brown, J., Newsome, W., and Pauly, K. (2013). Effective parameters for ultrasound-induced in vivo neurostimulation. Ultra. Med. Biol. 39, 312–331. doi: 10.1016/j.ultrasmedbio.2012.09.009
Kuzenko, D. V. (2021). Critical temperature below the curie temperature of ferroelectric ceramics PZT. J. Adv. Dielect. 11:2150006. doi: 10.1142/S2010135X21500065
Lee, W., Lee, S., Park, M., Foley, L., Purcell-Estabrook, E., Kim, H., et al. (2016). Image-guided focused ultrasound-mediated regional brain stimulation in sheep. Ultra. Med. Biol. 42, 459–470. doi: 10.1016/j.ultrasmedbio.2015.10.001
Li, G., Qiu, W., Zhang, Z., Jiang, Q., Su, M., Cai, R., et al. (2018). Noninvasive ultrasonic neuromodulation in freely moving mice. IEEE Trans. Biomed. Eng. 66, 217–224. doi: 10.1109/TBME.2018.2821201
Li, Z., Guo, R., Fei, C., Li, D., Chen, D., Zheng, C., et al. (2021a). Liquid lens with adjustable focus for ultrasonic imaging. Appl. Acoust. 175:107787. doi: 10.1016/j.apacoust.2020.107787
Li, Z., Wang, D., Fei, C., Qiu, Z., Hou, C., Wu, R., et al. (2021b). The forbidden band and size selectivity of acoustic radiation force trapping. iScience 24:101988. doi: 10.1016/j.isci.2020.101988
Li, Z., Yang, S., Fei, C., Guo, R., Chen, D., Zheng, C., et al. (2021c). Adjustable acoustic field controlled by “ultrasonic projector” on ultrasound application. IEEE Trans. Ultra. Ferroelect. Freq. Control 2021:3106712. doi: 10.1109/tuffc.2021.3106712
Lin, P., Zhu, Y., Fei, C., Chen, D., Chen, Z., Zheng, C., et al. (2021). Multilayer stairstep piezoelectric structure design for ultrabroad-bandwidth ultrasonic transducer. IEEE Sen. J. 21, 19889–19895. doi: 10.1109/jsen.2021.3100126
Mead, B. P., Kim, N., Miller, G. W., and Price, R. J. (2017). Novel focused ultrasound gene therapy approach noninvasively restores dopaminergic neuron function in a rat parkinson’s disease model. Nano Lett. 17, 3533–3542. doi: 10.1021/acs.nanolett.7b00616
Miniussi, C., and Vallar, G. (2011). Brain stimulation and behavioural cognitive rehabilitation: a new tool for neurorehabilitation. Neuropsychol. Rehabili. 21, 553–559. doi: 10.1080/09602011.2011.622435
Panczykowski, D., Monaco, E., and Friedlander, R. (2014). Transcranial focused ultrasound modulates the activity of primary somatosensory cortex in humans. Neurosurgery 74, N8–N9. doi: 10.1227/neu.0000000000000365
Qiu, Z., Kala, S., Guo, J., Xian, Q., Zhu, J., Zhu, T., et al. (2021). Targeted neurostimulation in mouse brains with non-invasive ultrasound. Cell Rep. 34:108595. doi: 10.1016/j.celrep.2020.108595
Simon, J., Sapozhnikov, O., Khokhlova, V., Crum, L., and Bailey, M. (2015). Ultrasonic atomization of liquids in drop-chain acoustic fountains. J. Fluid Mech. 766, 129–146. doi: 10.1017/jfm.2015.11
White, P. J., Clement, G. T., and Hynynen, K. (2006). Local frequency dependence in transcranial ultrasound transmission. Phys. Med. Biol. 51, 2293–2305. doi: 10.1088/0031-9155/51/9/013
Yamaoka, Y., Nambu, M., and Takamatsu, T. (2011). Fine depth resolution of two-photon absorption-induced photoacoustic microscopy using low-frequency bandpass filtering. Opt. Exp. 19, 13365–13377. doi: 10.1364/OE.19.013365
Yang, L., Song, L., Li, Q., and Zhang, T. (2021). Dielectric properties and electrical response of yttrium-doped Bi2/3Cu3Ti4 O12 ceramics. J. Adv. Dielect. 11:2150007. doi: 10.1142/S2010135X21500077
Yang, Y., Pacia, C. P., Ye, D., Zhu, L., Baek, H., Yue, Y., et al. (2021). Sonothermogenetics for noninvasive and cell-type specific deep brain neuromodulation. Brain Stimulat. 14, 790–800. doi: 10.1016/j.brs.2021.04.021
Yoo, S., Bystritsky, A., Lee, J., Zhang, Y., Fischer, K., Min, B., et al. (2011). Focused ultrasound modulates region-specific brain activity. NeuroImage 56, 1267–1275.
Keywords: half-concave structure, ultrasonic transducer, converge acoustic energy, neuromodulation, FEM simulation
Citation: Hou C, Wu Y, Fei C, Qiu Z, Li Z, Sun X, Zheng C and Yang Y (2022) An Optimized Miniaturized Ultrasound Transducer for Transcranial Neuromodulation. Front. Neurosci. 16:893108. doi: 10.3389/fnins.2022.893108
Received: 10 March 2022; Accepted: 26 May 2022;
Published: 21 June 2022.
Edited by:
Yaoheng Yang, Washington University in St. Louis, United StatesCopyright © 2022 Hou, Wu, Fei, Qiu, Li, Sun, Zheng and Yang. 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: Chunlong Fei, Y2xmZWlAeGlkaWFuLmVkdS5jbg==; Zhihai Qiu, cWl1emhpaGFpQGdkaWlzdC5jbg==
†These authors have contributed equally to this work