- 1The First College of Clinical Medicine, Shandong University of Traditional Chinese Medicine, Jinan, China
- 2Department of Ophthalmology and Optometry, Shandong University of Traditional Chinese Medicine, Jinan, China
- 3Shandong Province Key Laboratory of Integrated Traditional Chinese and Western Medicine for Prevention and Therapy of Ocular Disease, Eye Institute of Shandong University of Traditional Chinese Medicine, Jinan, China
- 4Affiliated Eye Hospital of Shandong University of Traditional Chinese Medicine, Jinan, China
Excessive glucocorticoids (GC) may lead to the aggravation of several basic diseases including myopia, due to plasma hormone imbalances associated with the hypothalamic–pituitary–adrenal axis (HPAA). Electroacupuncture (EA) is an effective therapeutic method to treat many diseases, although it remains unclear whether EA at acupoints on the foot or back would be effective in treating eye diseases. It was recently found that visual cortex activity for responses to visual stimuli with spatial frequency and resting-state functional connectivity (FC) between the supramarginal gyrus and rostrolateral prefrontal cortex was significantly reduced in patients with high myopia. The present study aims to investigate the role of the alternation of resting-state FC among the bilateral visual cortex and hypothalamus in exerting anti-myopia effects of EA in GC-enhanced lens-induced myopic (LIM) guinea pigs such that the mechanisms of EA to treat GC-enhanced myopia at Shenshu (BL23) acupoints can be probed. To confirm the effects of EA, ocular parameters including axial length and GC-associated physiological parameters such as animal appearance, behavior, bodyweight, and levels of four HPAA-associated plasma hormones [free triiodothyronine (FT3), free thyroxine (FT4), estradiol (E2), and testosterone (T)] were also collected. Increased resting-state FC between the left and right visual cortex was detected in GC-enhanced lens-induced myopic guinea pigs with EA at BL23 acupoints (LIM+GC+EA) guinea pigs compared to GC-enhanced lens-induced myopic guinea pigs with EA at sham acupoints (LIM+GC+Sham) guinea pigs, as well as suppressed myopia and recovery of symptoms initially caused by overdose of GC. Recovered symptoms included improved animal appearance, behavior, bodyweight, and HPAA-associated plasma hormone levels were observed after 4 weeks of EA treatment. In contrast, the LIM+GC+Sham group showed decreased FC with elongation of axial length for myopization as compared to the control group and LIM group and exhibited a deterioration in physiological parameters including reduced body weight and balance disruption in the four measured HPAA-associated plasma hormones. Our findings suggest that EA could effectively treat GC-enhanced myopia by increasing resting-state FC between the left and right visual cortices, which may be pivotal to further understanding the application and mechanisms of EA in treating GC-enhanced myopia.
Introduction
High myopia is a major public health concern, often accompanied by several severe comorbidities, including retinal detachment, cataracts, and glaucoma due to the elongation of axial length (Rudnicka et al., 2016; Morgan et al., 2018). Currently, excessive myopic axial length elongation and increased risk of irreversible visual impairment have been found in experimental lens-induced myopia after intraperitoneal injection of the glucocorticoids (GC) (Ding et al., 2018). The imbalance of four plasma hormones associated with the hypothalamic–pituitary–adrenal axis (HPAA), including free triiodothyronine (FT3), free thyroxine (FT4), estradiol (E2), and testosterone (T), is often caused by excess GC. Excess GC also causes deteriorated physical conditions as well as a reduction in body weight and then resulted in enhancement of basic disorders, such as arthritis and diabetes (De Bosscher and Haegeman, 2009; Lu et al., 2011; Vieira et al., 2011; Wang et al., 2015; Ferreira et al., 2016; Oray et al., 2016; Xia et al., 2017; Yan et al., 2017; Hasona, 2018; Panettieri et al., 2019). Interestingly, electroacupuncture (EA) has been proven to be an effective therapeutic method to treat GC-induced diseases at the Shenshu (BL23) acupoint which was located adjacent to the second lumbar vertebra on the back (Wang et al., 2015; Feng et al., 2018).
Using resting-state functional MRI (rsfMRI), functional connectivity (FC) between brain regions can be assessed by analyzing the temporal relationships of blood oxygen level-dependent (BOLD) fluctuations between brain regions (Chong et al., 2019; O'Neill et al., 2019). In recent years, various studies utilized rsfMRI to investigate the underlying mechanisms of eye diseases including myopia and amblyopia (Hu et al., 2018; Dai et al., 2019). Previous studies reported that visual cortex activity for responses to visual stimuli with spatial frequency, and resting-state FC between the supramarginal gyrus and the rostrolateral prefrontal cortex, was significantly reduced in high-myopia patients (Zhai et al., 2016; Mirzajani et al., 2017). It was also reported that FC density significantly decreased in the posterior cingulate cortex/precuneus (PCC/preCun) (Zhai et al., 2016). However, the role of resting-state FC between the left and right visual cortex on the treatment of high myopia has not yet been explored.
Consequently, there is considerable interest in discovering means to explore whether EA would affect brain function in treating GC-enhanced eye diseases at BL23 acupoints located on the back. In the present study, we aimed to investigate alternation of resting-state FC between the visual cortex and hypothalamus to assess the effects of EA at BL23 acupoints on the treatment of GC enhanced myopia in guinea pigs. We also measured myopia-related ocular parameters including axial length, and GC-associated physiological parameters including animal appearance, behavior, body weight, and levels of four plasma hormones related to the HPAA (FT3, FT4, E2, and T). Our study may provide insights in deepening understanding of the mechanisms of acupuncture in the treatment of GC-enhanced eye diseases at acupoints far from the eyes.
Materials and Methods
Animals
Sixty male pigmented guinea pigs (Cavia porcellus) at the age of 2–3 weeks were obtained from the Jinan Xijueling Laboratory Animal Ltd. (Jinan, China) and raised in the animal lab center within the Eye Institute of Shandong University of Traditional Chinese Medicine. Food and water for the guinea pigs were available ad libitum, and the room temperature was maintained at 22°C. The guinea pigs were reared in plastic cages (15 cm × 26 cm × 32 cm) under a 12/12 h light–dark cycle. The average light in the cage was ~300 lux. All experimental protocols and animal handling procedures were approved by the ethics committee of the Eye Institute of Shandong University of Traditional Chinese Medicine (2017-002)s and were in accordance with the statement of the Association for Research in Vision and Ophthalmology for the use of animals in vision and ophthalmological research.
GC Administration and LIM Establishment
The guinea pigs were randomly divided into four groups: control, LIM, LIM+GC+Sham, and LIM+GC+EA. The control group includes animals with no treatment (n = 15), and the LIM group comprises of animals with lens-induced unilateral myopization by goggles with a refractive power of −10 diopters glued onto the orbital rim of right eyes (n = 15). The LIM+GC+Sham group includes animals with lens-induced unilateral myopization of the right eyes and intraperitoneal injection of hydrocortisone in a dose of 10 mg/kg once daily (8:00–10:00 a.m.) for 2 consecutive weeks and then followed in a dose of 5 mg·kg−1 for the next 4 consecutive weeks to maintain the treatment effect, with EA at sham acupoints (n = 15), and the LIM+GC+EA group includes animals with lens-induced unilateral myopization of the right eyes and intraperitoneal injection of hydrocortisone in a dose of 10 mg/kg once daily (8:00–10:00 a.m.) for 2 consecutive weeks and then followed in a dose of 5 mg·kg−1 for the next 4 consecutive weeks to maintain the treatment effect, with EA at bilateral BL23 acupoints (n = 15).
The animals underwent body weight measurement at baseline and at each follow-up examination. The sonographic ocular biometry for axial length measurement was also collected by A/B-mode scan (oscillator frequency: 11 MHz; Quantel Co., Les Ulis, France) at these time points. One drop of 1% cyclopentolate hydrochloride (Alcon, USA) was applied to both eyes to achieve a completely dilated pupil and cycloplegia.
Electroacupuncture
After the combined treatment of lens-induced myopia and intraperitoneal injection of hydrocortisone for 2 consecutive weeks, the guinea pigs in the LIM+GC+EA group received EA at the bilateral BL23 point for 30 min a day for 4 consecutive weeks. BL23 is located adjacent to the second lumbar vertebra on the back (Xiang et al., 2019). The guinea pigs in the LIM+GC+Sham group were treated with EA at a sham point, which was set to the “degenerated tail” on the gluteus muscle, a point further away from the traditional meridians (Wang et al., 2007). The animals were lightly immobilized using a manufactured apparatus to minimize restraint stress, and acupuncture needles (40 mm in length, 0.30 mm in diameter) were bilaterally inserted to a depth of 8 mm at BL23 once a day (2:00 p.m.). Acupuncture needles were stimulated with an electrical-stimulator (Suzhou Medical Appliance Factory of China, Model SDZ-V), and parameters were set as continuous wave electrical pulses (0.1 ms duration), with a frequency 2 Hz and an intensity of 2 mA.
Serum Collection and Radioimmunoassay
The blood was drawn by cardiac puncture under anesthesia at 2 p.m. at the 0, 2, and 6 week intervals. Plasma FT3, FT4, E2, and T concentrations were determined in duplicate using standard radioimmunoassay (RIA) techniques by means of 125I-RIA kits with detection limits of 5 × 10–13 M, 1 × 10–12 M, 7.7 × 10–12 M, and 6.6 × 10–11 M, respectively. The different concentrations of the various hormones were measured using a gamma counter (GC-911, Anhui Ustc Zonkia Scientific Instruments Co., Ltd, China). The experimental steps were performed according to the protocols of the kits: (1) 125I FT3 and FT4 RIA kits (North Biotechnology Research Institute, Beijing) and (2) 125I E2 and T RIA kits (Tianjin JiuDing Medicine Bio-Engineering Co., Ltd, Tianjin).
MRI
Each group included six randomly selected guinea pigs for MRI after 4 weeks of EA. MRI was performed on a BioSpec 70/20 animal MRI system (Bruker BioSpin) equipped with a 7.0-T magnet with a horizontal bore 20 cm in diameter. The operating system was ParaVision 6.0.1, and the maximum gradient strength of the gradient system was 100 mT/m, using a low-temperature phased array receiver coil. During MRI, low-dose isoflurane (Ruiward Life Technology Co., Ltd., Shenzhen, China) was used (3.5% for induction and 1.5% for maintenance), which was slightly adjusted throughout the experiment to maintain a stable breathing frequency of 90 bpm. The animal respiratory rate was monitored using a PC-SAM Small Animal Monitor (SA Instruments). The four groups of guinea pigs were each placed in a separate animal bed equipped with circulating warm water to ensure that body temperature was maintained at 37–38°C through the heated animal bed.
Anatomical images covering the entire guinea pig brain were acquired using a multislice rapid acquisition with a relaxation enhancement sequence with repetition time (TR) = 175 ms, echo time (TE) = 4.5 ms, effective echo time (TE eff) 36 ms, number of averages (NA) = 1, and number of repetitions (NR) = 4, matrix dimension (MD) = 256 * 256, pixel dimensions (V) = 50 * 50 mm2, slice thickness (STH) = 1 mm, interslice distance (ISD) = 1 mm, and number of slices (NSl) = 30.
For fMRI, gradient-echo echo-planar imaging (EPI) was used with TR = 1,500 ms, TE = 20 ms, number of repetition (NR) = 180, NA = 1, MD = 256* 256, pixel dimensions = 25 * 220 mm2, slice thickness = 1 mm, interslice distance = 1 mm, number of slices = 30.
MRI Data Analysis
We segmented the label of the visual cortex and hypothalamus manually according to the sixth edition of The Rat Brain in Stereotaxic Coordinates (Figure 1). For a subject's high-resolution T2 anatomical image, we manually performed skull-stripping to remove non-brain tissues. The next step was inhomogeneity correction. Subsequently, a control subject's brain was chosen as the brain template. We registered each subject's T2 volume to the brain template using non-linear registration. Using the transformation generated by the registration algorithm, the subject's visual cortex region (left and right) and the hypothalamus were segmented (Valdes-Hernandez et al., 2011).
Figure 1. Anatomical locations of the visual cortex ROI (A), based on Paxinos et al., 1980. Anatomical locations of the hypothalamus ROI (B), based on Papp et al., 2014 [control (n = 6), LIM (n = 5), LIM+GC+Sham (n = 4), and LIM+GC+EA (n = 6)].
The rsfMRI images were realigned and corrected for slice timing. Afterward, a base EPI volume was extracted and skull stripping was performed to remove non-brain tissues. Linear detrending was applied for the removal of a systematic linear trend (Zerbi et al., 2014). The data was then band-pass filtered using a range between 0.01 and 0.3 Hz (Zerbi et al., 2018), and the subject's T2 volume was co-registered to the subject's fMRI base EPI volume. The visual cortex and the hypothalamus in the native fMR space were labeled based on the transformation generated by the registration algorithm. The left and right visual cortex and the hypothalamus were region-of-interests (ROIs), and the mean time courses of these three ROIs were extracted. The pairwise correlation coefficient among the ROIs was calculated to assess the functional connectivity between a brain region pair.
Statistical Analysis
SPSS software (Version 21.0) was used to perform statistical analyses. All data were expressed as mean ± SEM. An independent sample t-test was used to detect differences between groups. Statistical significance was considered when P < 0.05.
Results
Changes in Morphological Behavior, Body Weight, and Hormone Levels
At baseline, all groups had similar morphological behavior, body weight, and hormone levels of FT3, FT4, T, and E2 (Figures 2B–F). After 2 weeks of treatment, GC-treated groups (LIM+GC+Sham and LIM+GC+EA) showed deteriorated animal appearance such as dull coats, shivering, and decreased activity when compared to the LIM group and the control group. In addition, the GC-treated groups (LIM+GC+Sham and LIM+GC+EA) showed significantly decreased body weight when compared to the LIM group and the control group (Figure 2B). No significant differences were observed between LIM and controls (Figure 2B). The GC-treated groups (LIM+GC+Sham and LIM+GC+EA) also had significantly decreased concentrations of FT3, FT4, and T and significantly increased concentrations of E2 compared to LIM and control groups (Figures 2C–F).
Figure 2. The change in morphological behavior, body weight, and hormone level (n = 15 for each group). The changes in morphological behavior after the treatment for 6 w (A), and alterations of body weight (B) and hormones, including FT3 (C), FT4 (D), T (E), and E2 (F) in serum in all time intervals. #P < 0.05 compared with the LIM group; &P < 0.05, compared with LIM+GC+Sham group.
After 4 weeks of EA treatment (at the time point of “6 weeks,” as EA treatment was conducted after 2 weeks of LIM+GC treatment), the LIM+GC+EA group showed significantly increased body weight compared to the LIM+GC+Sham group, suggesting the GC-induced symptoms were ameliorated by EA (Figures 2A,B). In addition, the LIM+GC+EA guinea pigs also showed significantly increased FT3, FT4, and T concentrations with significantly decreased E2 when compared to the LIM+GC+Sham group (Figures 2C–F).
Changes of Axial Length
At baseline, a similar mean axial length of right eyes (treated) and left eyes (untreated) was observed in all groups (Table 1, Figure 3). In the control group, mean axial length of right (treated) eyes increased from 7.77 ± 0.07 mm (mean ± standard deviation) at baseline to 8.62 ± 0.12 mm at the end of follow-up, and mean axial length of left (untreated) eyes increased from 7.79 ± 0.03 at baseline to 8.64 ± 0.07 mm at the end of follow-up. There was no significant difference between treated and untreated eyes (Table 1). In the LIM group, the mean axial length of right eyes (treated) was significantly longer after 2 weeks of treatment compared to the control group (LIM group vs. control group, 8.20 ± 0.04 mm vs. 8.11 ± 0.04 mm, P < 0.001) and this increase continued in a time-dependent manner (Table 1, Figure 3A).
Table 1. Sonographic biometric measurements (mean ± standard deviation; OD: right eyes; OS: right eyes).
Figure 3. Changes in axial length at various intervals (n = 15 for each group). The changes in axial length of right eyes (A) and left eyes (B) in control, LIM, LIM+GC+Sham, and LIM+GC+EA group at the 0, 2, 4, and 6 week interval (n = 15 per group). *P < 0.05 compared with the control group; #P < 0.05 compared with the LIM group; &P < 0.05, compared with LIM+GC+Sham group.
In groups with unilateral lens-induced myopization and intraperitoneal injection of hydrocortisone (LIM+GC+Sham and LIM+GC+EA), the mean axial lengths of the right eyes (treated eyes) were significantly longer compared to the LIM group after combined treatment of GC and lens-induced myopization for 2 weeks, thus suggesting that excess GC significantly increased the degree of lens-induced myopia [(LIM+GC+Sham group vs. LIM group, 8.26 ± 0.10 mm vs. 8.20 ± 0.04 mm, P < 0.05; LIM+GC+EA group vs. LIM group, 8.26 ± 0.07 mm vs. 8.20 ± 0.04 mm, P < 0.05] (Table 1, Figure 3A). However, there was no significant difference in axial length between LIM+GC+Sham and LIM+GC+EA at that time point (LIM+GC+EA vs. LIM+GC+Sham, 8.26 ± 0.10 mm vs. 8.26 ± 0.07 mm, P = 0.833) (Table 1, Figure 3A).
Interestingly, after the treatment of EA at BL23 for 4 weeks, the axial length of the right eyes was significantly shorter in the LIM+GC+EA group compared to the LIM+GC+Sham group (LIM+GC+EA vs. LIM+GC+Sham, 8.79 ± 0.07 mm vs. 8.87 ± 0.04 mm, P < 0.05) (Table 1, Figure 3A). No significant difference was found in the axial lengths of the left eyes (untreated eyes) among each group at any time interval (Table 1, Figure 3B).
Altered Resting-State FC
Among the four groups, control (n = 6), LIM (n = 5), LIM+GC+Sham (n = 4), and LIM+GC+EA (n = 6), there were significant differences across the groups, as shown through linear regression modeling (p = 0.05). As shown in Figure 4, we found that the FC between the left and right visual cortices of the LIM group was significantly lower than that of the control (p = 0.02). The FC of the visual cortex of the LIM+GC+Sham group was also lower than that of control (p = 0.03). Of note, the FC of visual cortex of the LIM+GC+EA group is higher than that of LIM and LIM+GC+Sham groups. The difference in FC between the visual cortices was not significant between control and LIM+GC+EA; no significant difference between the LIM group and the LIM+GC+Sham group was found either (p > 0.05). Meanwhile, we did not find any significant difference for FC between the hypothalamus and each side of the visual cortex among all the four groups.
Figure 4. The boxplot of functional connectivity in the visual cortex; *compared with the control group, P < 0.05; #compared with the LIM group, P < 0.05; &compared with LIM+GC+Sham group, P < 0.05 [control (n = 6), LIM (n = 5), LIM+GC+Sham (n = 4), and LIM+GC+EA (n = 6)].
To further elucidate the relationship of FC between the visual cortex and parameters of animals, we conducted regression analyses between them. Results from the analyses showed a positive association between FC of the visual cortex and bodyweight, FT3 and T, whereas a negative association was found between FC of the visual cortex and E2 or axial length (Figure 5). However, no significant relationship was found between FC of the visual cortex and FT4.
Figure 5. Regression analyses between the functional connectivity in the visual cortex, and body weight, axial length, and four plasma hormones (FT3, FT4, T, E2). (A) Regression analyses between the functional connectivity in the visual cortex and body weight (A) and axial length (B). (C–F) Regression analyses between the functional connectivity in the visual cortex and the concentrations of FT3 (C), FT4 (D), T (E), and E2 (F) [control (n = 6), LIM (n = 5), LIM+GC+Sham (n = 4), and LIM+GC+EA (n = 6)].
Discussion
The prevalence of myopia has markedly increased within the past three decades, especially in China. The main characteristic of myopia is the irreversible elongation of axial length, resulting in many severe visual impairments including myopic macular degeneration, glaucoma, and even blind (Rudnicka et al., 2016; Morgan et al., 2018). In addition, overdose of GC could also induce excessive myopic axial length elongation and imbalance of four plasma hormones associated with hypothalamic–pituitary–adrenal axis (HPAA), including E2, T, FT3, and FT4 (Ding et al., 2018). The receptors of E2 and T are widely expressed in various ocular tissues (Wickham et al., 1998, 2000; Suzuki et al., 2001).
The guinea pigs have been a classic model to evaluate the effects of various therapeutic methods on myopia as well as its complications. Similar to the process of eye development in human beings, guinea pigs show hyperopia at birth and then rapidly become emmetropia within the first 3 weeks (21 days) of age (Shan et al., 2018). Their eyes have a more similar structure and biometric changes in the development of myopia as human beings, compared to other experimental model animals such as chicks and mice as well (Wu et al., 2020). In the present study, a low dose of isoflurane was used to maintain stable breathing frequency in anesthetized animals, and low-dose isoflurane maintained the resting-state networks of anesthetized animals to be similar to wake ones in rats and mice. The protection in resting-state neuron activity may result from a moderate systemic vasodilator effect due to increasing resting blood flow from low-dose isoflurane (Iida et al., 1998; Guilfoyle et al., 2013; Zhou et al., 2014). Moreover, isoflurane was not reported to be a risk factor for GC signaling or myopiazation, suggesting that it would be suitable for the application in the present study. Due to the relative ease in signal acquisition and proficiency of rsfMRI technology to measure the functional connectivity (FC) between functional areas of the brain in different populations, we choose the axial length and the resting-state functional connectivity using rsfMRI technology to assess the effects of EA at BL23 on GC-enhanced myopia. It strongly suggested that FC between the visual cortices played an important role in myopization.
FC in the Myopic Guinea Pig's Brain
It has been demonstrated that the amplitude of low-frequency fluctuation values in high-myopia patients are reduced in the right cerebellum anterior lobe/calcarine/bilateral parahippocampal gyrus, bilateral posterior cingulate cortex, and bilateral middle cingulate cortex, while they are significantly increased in the left optic radiation, bilateral frontal parietal cortex, and left primary motor cortex (M1)/primary somatosensory cortex (S1) (Huang et al., 2016; Cheng et al., 2020). Meanwhile, it was found that the high myopia exhibited significantly decreased short- and long-range FC densities in the posterior cingulate cortex/precuneus (PCC/preCun), with a similar result reported in amblyopia patients, which showed decreased FC between the PCC/preCun and bilateral primary visual areas (Ding et al., 2013; Zhai et al., 2016).
Our results indicated that FC between the visual cortex in the LIM+GC+Sham group and LIM group was significantly lower compared to the control group with the elongation of axial length. It exhibited many similarities with previous studies. It was reported that expression of neurotransmitters and their receptors changed in the primary visual cortex during the development of myopia (Zhao et al., 2017). Researches found that induced high myopia caused a significant reduction in the visual cortex activity by presenting a high range of spatial frequencies using functional MRI compared to the normal vision state (Mirzajani et al., 2017). Meanwhile, a resting-state functional magnetic resonance imaging study also demonstrated that low/moderate myopia and high myopia will lead to decreased neuronal and physiological activities in the primary visual cortex by studying the amplitude of low-frequency fluctuations (Cheng et al., 2020).
Effect of Excess GC on HPAA and Myopia
Evidence has shown that excessive GC resulted in deterioration of various diseases including arthritis and myopia through affecting secretion of certain hormones including FT3, FT4, E2, and T secreted from the hypothalamic pituitary target gland (adrenal, thyroid, and gonad) axis into target gland (Yang et al., 2008; Pace et al., 2009; Ding et al., 2018). It was found that concentrations of FT3, FT4, and T decreased while E2 elevated after the treatment of GC in the present study. These findings are consistent with previous studies that intraperitoneal injection of hydrocortisone, a type of GC, results in deteriorated animal appearance and a reduced body weight, accompanied with the suppression of the HPAA function by affecting the HPAA plasma hormone expression (Yang et al., 2008; Zhao et al., 2013, 2016). Besides, E2, as a member of estrogen, was reported to be a modulating factor that maintains the biomechanical properties and stability of the cornea and upregulate MMP-2 activity and protein expression in human retinal pigment epithelium cells, whereas T acted as an androgen of the steroid family and was reported to be associated with the biochemical characteristics of the sclera, the aqueous outflow pathway, and the iris/ciliary body (Knepper et al., 1985; Marin-Castano et al., 2003; Song et al., 2014). Also, cortisol administration elevated both the default mode network and salience network activity to a normal level to treat traumatic stress disorder and anxiety (Soravia et al., 2018). Also, it was found that intraperitoneal injection of GC can enhance myopic shift and axial elongation in guinea pigs with lens-induced myopia (Ding et al., 2018). These symptoms were also defined as “kidney-yang deficiency” in traditional Chinese medicine, which could be effectively treated by EA (Shen, 1999).
Previous studies showed that rapid intravenous infusion of hydrocortisone significantly increased the fMRI BOLD signal within the hippocampus in a time-dependent manner (Symonds et al., 2012). Additionally, it was also found that increased endogenous GC can elevate the resting-state FC of brain regions highly expressing GC receptors, such as the medial prefrontal cortex and medial temporal lobe (Stomby et al., 2019). Nevertheless, only FC between visual cortices was found to be significantly associated with alteration of levels of HPAA-associated hormones including FT3, T, and E2 in the serum, instead of FC between the hypothalamus and each side of the visual cortex. It might also be suggested that neuron signals between the hypothalamus and visual cortex may be through some unknown intermediate medium. Moreover, FT3 rather than FT4 played an important role in the relationship between GC and FC of the visual cortex.
Effects of EA on the Brain
Increasing evidence demonstrates that acupuncture at acupoints located in the body such as the limbs and trunk could effectively treat nervous system diseases including stroke, migraines, motor system diseases, and other diseases such as functional dyspepsia hypertension, overweight and Crohn's disease (Cai et al., 2018). It was reported that acupuncture could effectively enhance the FC between left primary motor area and left inferior frontal gyrus to promote the compensatory response to treat refractory facial paralysis, and enhance the functional connectivity between the precentral gyrus and the hippocampus in the Alzheimer disease patients (Zheng et al., 2018; Ma et al., 2019). Acupuncture could effectively enhance the FC between the left primary motor area and left inferior frontal gyrus to promote the compensatory response, increase connectivity between the periaqueductal gray, anterior cingulate cortex, left posterior cingulate cortex, right anterior insula, limbic/paralimbic, and precuneus, and adjust the limbic-paralimbic-neocortical network, brainstem, cerebellum, and subcortical and hippocampus brain areas (Cai et al., 2018).
Researches have shown that EA at acupoints located near the eyes, including Hegu (LI4) and Taiyang (EX-HN5), was effective to improve myopia by downregulating the level of retinal GABA in a myopic guinea pig model (Sha et al., 2015). It was also found that the stimulation of acupoints in the body, such as LR3 located on the feet, activated some areas of the visual cortex (Liu et al., 2012). With the development of rsfMRI technology, more research has focused on exploring different mechanisms of treatment including acupuncture for eye diseases such as high myopia, amblyopia, and blindness, by measuring FC between brain regions (Huang et al., 2016; Mendola et al., 2018; Wen et al., 2018). It was well established that the development of myopia is highly associated with alteration of function of visual cortex (Mirzajani et al., 2017). In the present study, our data revealed that EA at BL23 acupoints, located adjacent to the second lumbar vertebra on the back, could effectively suppress the elongation of axial length induced by a combination of treatment of GC and negative lens through recovering the balance of HPAA-associated plasma hormones and effectively recovering FC between the visual cortex of LIM+GC animals to normal levels, providing strong support to the notion that FC is related to the mechanism of acupuncture. It was consistent with the previous studies that acupuncture at a group of acupoints including BL23 as one of the major acupoints could relieve the symptoms of many disorders by improving the cerebral hemodynamics and cognitive deficits in the hippocampal CA1 region and rebalancing HPA-associated plasma hormones including E2, T, CORT, LH, and GnRH, or inhibiting the expression of orexin in the lateral hypothalamus (Wang et al., 2017, 2020; Zhang et al., 2017; Ji et al., 2019; Jing et al., 2020). Despite the limited research conducted on elucidating the effects of stimulating the BL23 acupoint as a single point to cure Kidney Yang deficiency-associated disorders, it was once reported that acupuncture at BL23 could effectively treat senescence-accelerated mice by increasing levels of serum hormone T (Zhang et al., 2009). These results would provide further evidence for the hypothesis that acupuncture could treat visual impairments including high myopia, through the alteration of the function of the visual cortex at acupoints located far from the eyes.
Nevertheless, several limitations of the study should be mentioned. First, there have been no brain atlases on guinea pigs so far. We therefore performed the MRI data analysis according to the rat brain atlas. Second, the fMRI image in this study was relatively low, so a higher image resolution would definitely strengthen our conclusion. Third, the specific mechanism of glucocorticoids aggravating the development of myopia needs further exploration.
In summary, EA could effectively treat GC-enhanced myopia by increasing resting-state FC between the left and right visual cortices at BL23, which may be pivotal in understanding the underlying mechanisms of EA in the treatment of GC-enhanced myopia.
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 author/s.
Ethics Statement
The animal study was reviewed and approved by the Ethics Committee of Eye Institute of Shandong University of Traditional Chinese Medicine.
Author Contributions
HB and WJ conceived and formulated the research. TZ, QJ, FX, and RZ raised animals and performed treatment. WJ, QJ, FX, DL, and DG conducted rsfMRI measurement. WJ and QJ analyzed the data. WJ, QJ, and TZ wrote the paper. All authors read and approved the final manuscript.
Funding
This work was supported by National Key Research and Development Program (Nos. 2019YFC1710200 and 2019YFC1710204) and Key Research and Development Program of Shandong Province (2017CXGC1211, 2018JHZ005, and 2019GSF108252).
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.
Acknowledgments
We would like to thank Alice Chen, Schulich School of Medicine and Dentistry, Western University, London, ON, Canada, for her assistance in editing the manuscript. We also acknowledge Drs. Xianglin Li and Zhenbo Sun for their technical support in MRI acquisition.
References
Cai, R. L., Shen, G. M., Wang, H., and Guan, Y. Y. (2018). Brain functional connectivity network studies of acupuncture: a systematic review on resting-state fMRI. J. Integr. Med. 16, 26–33. doi: 10.1016/j.joim.2017.12.002
Cheng, Y., Huang, X., Hu, Y. X., Huang, M. H., Yang, B., Zhou, F. Q., et al. (2020). Comparison of intrinsic brain activity in individuals with low/moderate myopia versus high myopia revealed by the amplitude of low-frequency fluctuations. Acta Radiol. 61, 496–507. doi: 10.1177/0284185119867633
Chong, C. D., Schwedt, T. J., and Hougaard, A. (2019). Brain functional connectivity in headache disorders: a narrative review of MRI investigations. J. Cereb. Blood Flow Metab. 39, 650–669. doi: 10.1177/0271678X17740794
Dai, P., Zhang, J., Wu, J., Chen, Z., Zou, B., Wu, Y., et al. (2019). Altered spontaneous brain activity of children with unilateral amblyopia: a resting state fMRI study. Neural Plast 2019:3681430. doi: 10.1155/2019/3681430
De Bosscher, K., and Haegeman, G. (2009). Minireview: latest perspectives on antiinflammatory actions of glucocorticoids. Mol. Endocrinol. 23, 281–291. doi: 10.1210/me.2008-0283
Ding, K., Liu, Y., Yan, X., Lin, X., and Jiang, T. (2013). Altered functional connectivity of the primary visual cortex in subjects with amblyopia. Neural Plast 2013:612086. doi: 10.1155/2013/612086
Ding, M., Guo, D., Wu, J., Ye, X., Zhang, Y., Sha, F., et al. (2018). Effects of glucocorticoid on the eye development in guinea pigs. Steroids 139, 1–9. doi: 10.1016/j.steroids.2018.09.008
Feng, Y., Fang, Y., Wang, Y., and Hao, Y. (2018). Acupoint therapy on diabetes mellitus and its common chronic complications: a review of its mechanisms. Biomed. Res. Int. 2018:3128378. doi: 10.1155/2018/3128378
Ferreira, J. F., Ahmed Mohamed, A. A., and Emery, P. (2016). Glucocorticoids and rheumatoid arthritis. Rheum. Dis. Clin. North Am. 42, 33–46. doi: 10.1016/j.rdc.2015.08.006
Guilfoyle, D. N., Gerum, S. V., Sanchez, J. L., Balla, A., Sershen, H., Javitt, D. C., et al. (2013). Functional connectivity fMRI in mouse brain at 7T using isoflurane. J. Neurosci. Methods 214, 144–148. doi: 10.1016/j.jneumeth.2013.01.019
Hasona, N. A. (2018). Grape seed extract attenuates dexamethasone-induced testicular and thyroid dysfunction in male albino rats. Andrologia 50:e13002. doi: 10.1111/and.13002
Hu, Y. X., He, J. R., Yang, B., Huang, X., Li, Y. P., Zhou, F. Q., et al. (2018). Abnormal resting-state functional network centrality in patients with high myopia: evidence from a voxel-wise degree centrality analysis. Int. J. Ophthalmol. 11, 1814–1820. doi: 10.18240/ijo.2018.11.13
Huang, X., Zhou, F. Q., Hu, Y. X., Xu, X. X., Zhou, X., Zhong, Y. L., et al. (2016). Altered spontaneous brain activity pattern in patients with high myopia using amplitude of low-frequency fluctuation: a resting-state fMRI study. Neuropsychiatr. Dis. Treat. 12, 2949–2956. doi: 10.2147/NDT.S118326
Iida, H., Ohata, H., Iida, M., Watanabe, Y., and Dohi, S. (1998). Isoflurane and sevoflurane induce vasodilation of cerebral vessels via ATP-sensitive K+ channel activation. Anesthesiology 89, 954–960. doi: 10.1097/00000542-199810000-00020
Ji, Y. H., Ji, Y. H., and Sun, B. D. (2019). [Effect of acupuncture combined with repetitive transcranial magnetic stimulation on motor function and cerebral hemodynamics in children with spastic cerebral palsy with spleen-kidney deficiency]. Zhen Ci Yan Jiu 44, 757–761. doi: 10.13702/j.1000-0607.190154
Jing, Q., Ren, L., Deng, X., Zhang, N., Fu, M., Wang, G., et al. (2020). Electroacupuncture promotes neural proliferation in hippocampus of perimenopausal depression rats via wnt/β catenin signaling pathway. J. Acupunct. Meridian. Stud.. 13, 94–103. doi: 10.1016/j.jams.2020.03.065
Knepper, P. A., Collins, J. A., and Frederick, R. (1985). Effects of dexamethasone, progesterone, and testosterone on IOP and GAGs in the rabbit eye. Invest. Ophthalmol. Vis. Sci. 26, 1093–1100.
Liu, H., Xu, J., Shan, B., Li, Y., Li, L., Xue, J., et al. (2012). Determining the precise cerebral response to acupuncture: an improved FMRI study. PLoS ONE 7:e49154. doi: 10.1371/journal.pone.0049154
Lu, X., Xiong, Z., Li, J., Zheng, S., Huo, T., and Li, F. (2011). Metabonomic study on ‘Kidney-yang deficiency syndrome’ and intervention effects of rhizoma drynariae extracts in rats using ultra performance liquid chromatography coupled with mass spectrometry. Talanta 83, 700–708. doi: 10.1016/j.talanta.2010.09.026
Ma, B., Yang, J., Yuan, A. H., Li, C. F., Gao, T., and Song, S. T. (2019). [Effect of acupuncture and moxibustion on brain functional connectivity network in patients with refractory facial paralysis]. Zhongguo Zhen Jiu 39, 1321–1326. doi: 10.13703/j.0255-2930.2019.12.017
Marin-Castano, M. E., Elliot, S. J., Potier, M., Karl, M., Striker, L. J., Striker, G. E., et al. (2003). Regulation of estrogen receptors and MMP-2 expression by estrogens in human retinal pigment epithelium. Invest. Ophthalmol. Vis. Sci. 44, 50–59. doi: 10.1167/iovs.01-1276
Mendola, J. D., Lam, J., Rosenstein, M., Lewis, L. B., and Shmuel, A. (2018). Partial correlation analysis reveals abnormal retinotopically organized functional connectivity of visual areas in amblyopia. Neuroimage Clin. 18, 192–201. doi: 10.1016/j.nicl.2018.01.022
Mirzajani, A., Ghorbani, M., Rasuli, B., and Mahmoud-Pashazadeh, A. (2017). Effect of induced high myopia on functional MRI signal changes. Phys. Med. 37, 32–36. doi: 10.1016/j.ejmp.2017.04.004
Morgan, I. G., French, A. N., Ashby, R. S., Guo, X., Ding, X., He, M., et al. (2018). The epidemics of myopia: aetiology and prevention. Prog. Retin. Eye Res. 62, 134–149. doi: 10.1016/j.preteyeres.2017.09.004
O'Neill, J., Piacentini, J. C., and Peterson, B. S. (2019). Cingulate role in tourette syndrome. Handb. Clin. Neurol. 166, 165–221. doi: 10.1016/B978-0-444-64196-0.00011-X
Oray, M., Abu Samra, K., Ebrahimiadib, N., Meese, H., and Foster, C. S. (2016). Long-term side effects of glucocorticoids. Expert Opin. Drug Saf. 15, 457–465. doi: 10.1517/14740338.2016.1140743
Pace, T. W., Gaylord, R. I., Jarvis, E., Girotti, M., and Spencer, R. L. (2009). Differential glucocorticoid effects on stress-induced gene expression in the paraventricular nucleus of the hypothalamus and ACTH secretion in the rat. Stress 12, 400–411. doi: 10.1080/10253890802530730
Panettieri, R. A., Schaafsma, D., Amrani, Y., Koziol-White, C., Ostrom, R., and Tliba, O. (2019). Non-genomic effects of glucocorticoids: an updated view. Trends Pharmacol. Sci. 40, 38–49. doi: 10.1016/j.tips.2018.11.002
Papp, E. A., Leergaard, T. B., Calabrese, E., Johnson, G. A., and Bjaalie, J. G. (2014). Waxholm Space atlas of the Sprague Dawley rat brain. Neuroimage 97, 374–386. doi: 10.1016/j.neuroimage.2014.04.001.2014.04.001
Paxinos, G., Watson, C. R., and Emson, P. C. (1980). AChE-stained horizontal sections of the rat brain in stereotaxic coordinates. J. Neurosci. Methods 3, 129–149. doi: 10.1016/0165-0270(80)90021-7
Rudnicka, A. R., Kapetanakis, V. V., Wathern, A. K., Logan, N. S., Gilmartin, B., Whincup, P. H., et al. (2016). Global variations and time trends in the prevalence of childhood myopia, a systematic review and quantitative meta-analysis: implications for aetiology and early prevention. Br. J. Ophthalmol. 100, 882–890. doi: 10.1136/bjophthalmol-2015-307724
Sha, F., Ye, X., Zhao, W., Xu, C. L., Wang, L., Ding, M. H., et al. (2015). Effects of electroacupuncture on the levels of retinal gamma-aminobutyric acid and its receptors in a guinea pig model of lens-induced myopia. Neuroscience 287, 164–174. doi: 10.1016/j.neuroscience.2014.12.022
Shan, S. W., Tse, D. Y., Zuo, B., To, C. H., Liu, Q., McFadden, S. A., et al. (2018). Integrated SWATH-based and targeted-based proteomics provide insights into the retinal emmetropization process in guinea pig. J. Proteomics 181, 1–15. doi: 10.1016/j.jprot.2018.03.023
Song, X., Zhao, P., Wang, G., and Zhao, X. (2014). The effects of estrogen and androgen on tear secretion and matrix metalloproteinase-2 expression in lacrimal glands of ovariectomized rats. Invest. Ophthalmol. Vis. Sci. 55, 745–751. doi: 10.1167/iovs.12-10457
Soravia, L. M., Schwab, S., Weber, N., Nakataki, M., Wiest, R., Strik, W., et al. (2018). Glucocorticoid administration restores salience network activity in patients with spider phobia. Depress. Anxiety 35, 925–934. doi: 10.1002/da.22806
Stomby, A., Salami, A., Dahlqvist, P., Evang, J. A., Ryberg, M., Bollerslev, J., et al. (2019). Elevated resting-state connectivity in the medial temporal lobe and the prefrontal cortex among patients with Cushing's syndrome in remission. Eur. J. Endocrinol. 180, 329–338. doi: 10.1530/EJE-19-0028
Suzuki, T., Kinoshita, Y., Tachibana, M., Matsushima, Y., Kobayashi, Y., Adachi, W., et al. (2001). Expression of sex steroid hormone receptors in human cornea. Curr. Eye Res. 22, 28–33. doi: 10.1076/ceyr.22.1.28.6980
Symonds, C. S., McKie, S., Elliott, R., William Deakin, J. F., and Anderson, I. M. (2012). Detection of the acute effects of hydrocortisone in the hippocampus using pharmacological fMRI. Eur. Neuropsychopharmacol. 22, 867–874. doi: 10.1016/j.euroneuro.2012.03.008
Valdes-Hernandez, P. A., Sumiyoshi, A., Nonaka, H., Haga, R., Aubert-Vasquez, E., Ogawa, T., et al. (2011). An in vivo MRI template set for morphometry, tissue segmentation, and fMRI localization in rats. Front. Neuroinform. 5:26. doi: 10.3389/fninf.2011.00026
Vieira, J. S., Saraiva, K. L., Barbosa, M. C., Porto, R. C., Cresto, J. C., Peixoto, C. A., et al. (2011). Effect of dexamethasone and testosterone treatment on the regulation of insulin-degrading enzyme and cellular changes in ventral rat prostate after castration. Int. J. Exp. Pathol. 92, 272–280. doi: 10.1111/j.1365-2613.2011.00772.x
Wang, C., Zhou, D. F., Shuai, X. W., Liu, J. X., and Xie, P. Y. (2007). Effects and mechanisms of electroacupuncture at PC6 on frequency of transient lower esophageal sphincter relaxation in cats. World J. Gastroenterol. 13, 4873–4880. doi: 10.3748/wjg.v13.i36.4873
Wang, S. J., Zhang, J. J., Yang, H. Y., Wang, F., and Li, S. T. (2015). Acupoint specificity on acupuncture regulation of hypothalamic- pituitary-adrenal cortex axis function. BMC Complement. Altern. Med. 15:87. doi: 10.1186/s12906-015-0625-4
Wang, X., Zhang, B., Zhang, L., and Liu, S. (2017). Electroacupuncture suppresses morphine reward-seeking behavior: lateral hypothalamic orexin neurons implicated. Neurosci. Lett. 661, 84–89. doi: 10.1016/j.neulet.2017.09.057
Wang, Y. Y., Gao, M., Qiu, G. P., Xu, J., Wu, Y. J., Xu, Y., et al. (2020). [Effect of electroacupuncture on the P35/P25-cyclin-dependent kinase 5-Tau pathway in hippocampus of rats with Alzheimer's disease]. Zhen Ci Yan Jiu 45, 194–201. doi: 10.13702/j.1000-0607.190415
Wen, Z., Zhou, F. Q., Huang, X., Dan, H. D., Xie, B. J., and Shen, Y. (2018). Altered functional connectivity of primary visual cortex in late blindness. Neuropsychiatr. Dis. Treat. 14, 3317–3327. doi: 10.2147/NDT.S183751
Wickham, L. A., Gao, J., Toda, I., Rocha, E. M., Ono, M., and Sullivan, D. A. (2000). Identification of androgen, estrogen and progesterone receptor mRNAs in the eye. Acta Ophthalmol. Scand. 78, 146–153. doi: 10.1034/j.1600-0420.2000.078002146.x
Wickham, L. A., Rocha, E. M., Gao, J., Krenzer, K. L., da Silveira, L. A., Toda, I., et al. (1998). Identification and hormonal control of sex steroid receptors in the eye. Adv. Exp. Med. Biol. 438, 95–100. doi: 10.1007/978-1-4615-5359-5_12
Wu, S., Guo, D., Wei, H., Yin, X., Zhang, L., Guo, B., et al. (2020). Disrupted potassium ion homeostasis in ciliary muscle in negative lens-induced myopia in Guinea pigs. Arch. Biochem. Biophys. 688:108403. doi: 10.1016/j.abb.2020.108403
Xia, D., Chen, P., Du, P., Ding, L., and Liu, A. (2017). [Efficacy of acupoint catgut embedding combined with ginger-partitioned moxibustion on chronic fatigue syndrome of spleen-kidney yang deficiency syndrome and its effects on T lymphocyte subsets and activity of NK cell]. Zhongguo Zhen Jiu 37, 814–818. doi: 10.13703/j.0255-2930.2017.08.004
Xiang, S. J., Li, M. H., Chan, C. O., Shen, Q., Chen, S. B., An, B. C., et al. (2019). Altered metabolites in guinea pigs with allergic asthma after acupoint sticking therapy: new insights from a metabolomics approach. Phytomedicine 54, 182–194. doi: 10.1016/j.phymed.2018.09.021
Yan, Y., Pan, C., Wu, J., and Yang, L. (2017). [Efficacy of acupuncture-moxibustion periodic therapy on follicular maldevelopment differentiated as spleen and kidney yang deficiency]. Zhongguo Zhen Jiu 37, 39–44. doi: 10.13703/j.0255-2930.2017.01.009
Yang, J., Wang, Y., Bao, Y., and Guo, J. (2008). The total flavones from Semen cuscutae reverse the reduction of testosterone level and the expression of androgen receptor gene in kidney-yang deficient mice. J. Ethnopharmacol. 119, 166–171. doi: 10.1016/j.jep.2008.06.027
Zerbi, V., Ielacqua, G. D., Markicevic, M., Haberl, M. G., Ellisman, M. H., A, A. B., et al. (2018). Dysfunctional autism risk genes cause circuit-specific connectivity deficits with distinct developmental trajectories. Cereb. Cortex 28, 2495–2506. doi: 10.1093/cercor/bhy046
Zerbi, V., Wiesmann, M., Emmerzaal, T. L., Jansen, D., Van Beek, M., Mutsaers, M. P., et al. (2014). Resting-state functional connectivity changes in aging apoE4 and apoE-KO mice. J. Neurosci. 34, 13963–13975. doi: 10.1523/JNEUROSCI.0684-14.2014
Zhai, L., Li, Q., Wang, T., Dong, H., Peng, Y., Guo, M., et al. (2016). Altered functional connectivity density in high myopia. Behav. Brain Res. 303, 85–92. doi: 10.1016/j.bbr.2016.01.046
Zhang, M., Xv, G. H., Wang, W. X., Meng, D. J., and Ji, Y. (2017). Electroacupuncture improves cognitive deficits and activates PPAR-gamma in a rat model of Alzheimer's disease. Acupunct. Med. 35, 44–51. doi: 10.1136/acupmed-2015-010972
Zhang, X., Peng, Y., Yu, J., Liu, C., Cheng, H., Liu, L., et al. (2009). Changes in histomorphometric and mechanical properties of femurs induced by acupuncture at the Shenshu point in the SAMP6 mouse model of senile osteoporosis. Gerontology 55, 322–332. doi: 10.1159/000214845
Zhao, L., Wu, H., Qiu, M., Sun, W., Wei, R., Zheng, X., et al. (2013). Metabolic signatures of kidney yang deficiency syndrome and protective effects of two herbal extracts in rats using GC/TOF MS. Evid. Based Complement. Alternat. Med. 2013:540957. doi: 10.1155/2013/540957
Zhao, L., Zhao, A., Chen, T., Chen, W., Liu, J., Wei, R., et al. (2016). Global and targeted metabolomics evidence of the protective effect of chinese patent medicine jinkui shenqi pill on adrenal insufficiency after acute glucocorticoid withdrawal in rats. J. Proteome Res. 15, 2327–2336. doi: 10.1021/acs.jproteome.6b00409
Zhao, W., Bi, A. L., Xu, C. L., Ye, X., Chen, M. Q., Wang, X. T., et al. (2017). GABA and GABA receptors alterations in the primary visual cortex of concave lens-induced myopic model. Brain Res. Bull. 130, 173–179. doi: 10.1016/j.brainresbull.2017.01.017
Zheng, W., Su, Z., Liu, X., Zhang, H., Han, Y., Song, H., et al. (2018). Modulation of functional activity and connectivity by acupuncture in patients with Alzheimer disease as measured by resting-state fMRI. PLoS ONE 13:e0196933. doi: 10.1371/journal.pone.0196933
Keywords: functional connectivity, resting-state fMRI, neuroimaging, electroacupuncture, glucocorticoid, myopia
Citation: Zhang T, Jiang Q, Xu F, Zhang R, Liu D, Guo D, Wu J, Wen Y, Wang X, Jiang W and Bi H (2021) Alternation of Resting-State Functional Connectivity Between Visual Cortex and Hypothalamus in Guinea Pigs With Experimental Glucocorticoid Enhanced Myopia After the Treatment of Electroacupuncture. Front. Neuroinform. 14:579769. doi: 10.3389/fninf.2020.579769
Received: 03 July 2020; Accepted: 07 December 2020;
Published: 13 January 2021.
Edited by:
Jinghao Zhou, University of Maryland, Baltimore, United StatesCopyright © 2021 Zhang, Jiang, Xu, Zhang, Liu, Guo, Wu, Wen, Wang, Jiang and Bi. 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: Wenjun Jiang, bGl6X2ppYW5nJiN4MDAwNDA7MTM5LmNvbQ==; Hongsheng Bi, aG9uZ3NoZW5nYmkxJiN4MDAwNDA7MTYzLmNvbQ==
†These authors have contributed equally to this work and share first authorship