• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Lithium promotes recovery of neurological function after spinal cord injury by inducing autophagy

    2018-10-22 03:33:24DuoZhangFangWangXuZhaiXiaoHuiLiXiJingHe

    Duo Zhang , Fang Wang, Xu Zhai, Xiao-Hui Li, Xi-Jing He,

    1 Department of Orthopedics, Beijing Tiantan Hospital, Capital Medical University, Beijing, China

    2 Department of Orthopedics, Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi Province, China

    3 Department of Emergency, Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi Province, China

    4 Department of Radiology, Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, Shaanxi Province, China

    Abstract

    Key Words: nerve regeneration; spinal cord injury; lithium; secondary injury; autophagy, diffusion tensor imaging; neuroprotection; functional recovery; immunohistochemistry; Beclin-1; light-chain 3B; neural regeneration

    Introduction

    Spinal cord injury (SCI) is the direct or indirect damage to any part of the spinal cord that results in permanent impairment in strength, sensation or other body function below the injury site (Wyndaele and Wyndaele, 2006; Krause et al.,2017). Secondary injury mechanisms play important roles in the acute, sub-acute and chronic phases and lead to vasospasm, ischemia, inflammation, and free radical production(Oyinbo, 2011). These events result in neuronal loss, which is the key cause of permanent neurological dysfunction(Kanno et al., 2009; Tang et al., 2014; Li et al., 2015a; Kwan et al., 2017).

    Cell autophagy or type II programmed cell death is an intracellular catabolic mechanism for recycling damaged organelles and senescent proteins and plays a very important role in cell survival, differentiation and homeostasis (Erlich et al., 2007; Mizushima et al., 2010). It was reported that enhancing autophagy promotes the recovery of neurological functions by inhibiting apoptosis (Sekiguchi et al., 2012; Liu et al., 2015; Colón and Miranda, 2016; Dai et al., 2017). Sev-eral studies have reported that autophagy is involved in SCI.Lysosomal dysfunction and the disruption of autophagy, as well as increased neuronal apoptosis, have been observed in SCI, suggesting that autophagy is involved in SCI (Silva et al., 2008; Liu et al., 2015). Accordingly, an increasing number of studies have focused on the therapeutic effect of modulating autophagy in SCI.

    Lithium is the first line drug for treating bipolar disorder, and provides neuroprotection in multiple neurologic diseases (Young, 2009; Huo et al., 2012). Accumulating evidence suggests that lithium has numerous actions, including neuroprotection, inflammation inhibition, induction of neurotrophic factor secretion, and the enhancement of neurogenesis (Son et al., 2003; Senatorov et al., 2004; Su et al., 2007; Yasuda et al., 2009; Yuskaitis and Jope, 2009; Chi-Tso and Chuang, 2011; De Meyer et al., 2011; Li et al., 2011;Lauterbach, 2016). However, it remains unclear whether autophagy plays a positive or negative role in SCI (O’Donovan et al., 2015; Del Grosso et al., 2016; Fabrizi et al., 2017).

    Multiple signaling pathways are involved in autophagy,including the PI3K/Akt/mTOR, AMPK/TSC/mTOR and eIF2α/Atg12 pathways (Periyasamy-Thandavan et al., 2009).It was reported that lithium affects several signaling pathways, including the PI3K/Akt and IP3 pathways, which are involved in autophagy. Therefore, it is reasonable to speculate that lithium promotes functional recovery by inducing autophagy in rat models of SCI, which has not yet been studied (Periyasamy-Thandavan et al., 2009; Chi-Tso and Chuang, 2011; Li et al., 2015b).

    In the present study, we investigate the neuroprotective effect of lithium and the role of autophagy in SCI using the autophagy inhibitor 3-methyladenine (3-MA). To objectively and accurately evaluate recovery following SCI, we performed diffusion tensor imaging (DTI), which is an effective method of assessing neurological recovery, in addition to the Basso, Beattie, Bresnahan (BBB) locomotor rating scale(Zhang et al., 2015).

    Materials and Methods

    Animal care and groups

    A total of 72 specific-pathogen-free adult healthy male Sprague-Dawley rats weighing 230–270 g were provided by the Experimental Animal Center of Xi’an Jiaotong University of China (production license No. SCXK (Shaan) 2007-001; user license No. SYXK (Shaan) 2007-003). All rats were housed under a 12-hour light/dark cycle. All experimental procedures were performed in accordance with the Guidelines for the Care and Use of Laboratory Animals published by the US National Institutes of Health. The protocols were approved by the Animal Ethics Committee of Xi’an Jiaotong University of China. All efforts were made to minimize distress to the rats.

    The rats were randomly separated into sham, SCI, lithium and 3-MA groups (n = 18 per group). Six rats from each group were randomly selected for BBB scoring and DTI examination, and three rats from each group were sacrificed for immunohistochemical staining at each time point.

    SCI model

    Rats in the SCI, lithium and 3-MA groups received SCI operation as previously described (Basso et al., 1995). Briefly, the rat was anesthetized with intraperitoneal injection of chloral hydrate (300 mg/kg) and placed in the prone position with a heating pad to maintain body temperature. After shaving and aseptic preparation, the spinal cord was exposed. Dorsal laminectomy was performed at T9–11. The T10 segment of the spinal cord was impacted with an NYU weight-drop impactor (10 g rod dropped from a height of 25 mm; Rutgers University, USA), which led to hemorrhage and edema at the site of impact, wagging tail reflex and lower limb spasm. All manifestations indicated the success of the injury model. The rats in the sham group underwent laminectomy alone. The tissue was sutured layer by layer, with a piece of fat sutured under the skin at the T10 level. After SCI surgery, manual bladder massage was performed three times, and intraperitoneal injection of penicillin 20 U/kg was given once daily until bladder function was reestablished.

    Lithium and 3-MA treatments

    Lithium chloride (LiCl; Kemiou, Tianjin, China) and 3-MA(Sigma-Aldrich, St. Louis, MO, USA) were dissolved in 0.9%NaCl. Rats in the 3-MA group were given intraperitoneal injection of 1 mL 3-MA (3 mg/kg) 2 hours before SCI (Chen et al., 2013; Tang et al., 2014). Rats in the lithium and 3-MA groups were administered 1 mL lithium by intraperitoneal injection (LiCl, 30 mg/kg) 6 hours after SCI and then once daily until sacrifice (Yick et al., 2004; Zakeri et al., 2014). The sham and SCI groups received 1 mL 0.9% NaCl via intraperitoneal injection.

    Neurological function assessment

    The BBB locomotor rating scale was used to assess neurological function after SCI (Zhang et al., 2015). The BBB score ranges from 0 (complete paralysis) to 21 (normal), based on the range and extent of motion, weight loading, coordination of the forelimbs and hindlimbs, and motion of the forepaw, hindpaw and tail. Three independent examiners blindly assessed the BBB score before operation, and at 6 hours and 1, 2, 3 and 4 weeks after SCI. The average score was taken as thefinal score for each rat at each time point.

    DTI examination

    A DTI scan was performed 24 hours before SCI and at 6 hours and 1, 2, 3 and 4 weeks after SCI using a 3.0 T SIGNA MRI scanner (GE Medical Systems, Milwaukee, WI,USA) at the same loci as the conventional MRI scan. The scanning parameters were as follows: diffusion-weighted coefficient (b-value) = 1000 s/mm2; diffusion-sensitive gradient = 15 different directions; repetition time = 3500 ms;echo time = 87.5 ms; thickness = 2.4 mm; space = 0;field of view = 10; acquisition matrix = 64 × 64. All data were input into a workstation running Advantage Windows 4.2 (GE Healthcare). The region of interest (ROI) was identified by the fat under the skin, which was displayed as a high signal on conventional T2WI MRI (Yan et al., 2007). Based on the fractional anisotropy (FA) map, the ROI was placed in the inferior medulla and the inferior oblongata. The ROI was selected by two independent testers, and apparent diffusion coefficient (ADC) and FA values were obtained. FA values reflect the degree of spatial displacement of water molecules,and higher FA values indicate stronger anisotropy. ADC values are independent of the diffusion directions, and indicate the diffusional displacement of water molecules.

    Immunohistochemical staining

    Three rats were randomly selected from the sham group,while three rats were randomly selected from the other groups at 6 hours and at 1, 3 and 7 days for immunohistochemistry.The rats were anesthetized with chloral hydrate (300 mg/kg)and received aortic cannulation through the apex of the leftventricle. The rat was then perfused with 4% paraformaldehyde. A 2.0-cm spinal cord tissue segment centered at the injury site (T10 segment) was harvested and immersed in 4%paraformaldehyde at 4°C for 12 hours, and thereafter, 6 μmthick coronal paraffin sections were prepared.

    Three slices from each rat were selected for immunohistochemical staining for the neuronal nuclear antigen NeuN and the autophagy markers Beclin-1 and microtubule-associated proteins 1A/1B light-chain 3B (LC3B). Briefly, the sections were deparaffinized, rehydrated in distilled water, placed in 3% H2O2to remove residual peroxidase, and then rinsed with phosphate-buffered saline (PBS). The slices were blocked with 10% normal goat serum for 2 hours following permeabilization with 0.1% Triton X-100. Afterwards, the sections were incubated with anti-NeuN antibody (rabbit anti-rat IgG; 1:1000; Abcam, Cambridge, UK), anti-Beclin-1 antibody(mouse anti-rat IgG; 1:400; Santa Cruz Biotechnology, Dallas,TX, USA) or anti-LC3B antibody (mouse anti-rat IgG; 1:400;Santa Cruz Biotechnology) at 4°C for 24 hours. The sections were then incubated with the corresponding horseradish peroxidase-labeled antibody (goat anti-mouse; 1:100; Beyotime,Shanghai, China) at 37°C for 30 minutes, followed by three washes with PBS. Specific staining was visualized with diaminobenzidine according to the supplier’s instructions (Beyotime), followed by counterstaining with hematoxylin. Finally,sections were washed with PBS, dehydrated through a graded alcohol series (50%, 75%, 95%, 100%), cleared with dimethylbenzene, and mounted using a coverslip.

    For analysis, four randomly selected fields were photographed at 400× magnification on a microscope (Olympus,Tokyo, Japan). NeuN-positive cells at the injury site were counted manually and blindly by three examiners. Images of Beclin-1 and LC3B-stained sections were imported into Image Plus Pro 6.0 software (Media Cybernetics, Rockville,MD, USA) to quantify the positively-stained area. Relative area, which was defined as the ratio of the average area in the experimental group to that in the sham group, was analyzed to compare autophagy levels among the groups.

    Statistical analysis

    Data are expressed as the mean ± SD. Statistical analysis was performed using SPSS 20.0 software (IBM Corporation, Armonk, NY, USA). One-way analysis of variance followed by the least significant difference post hoc test was used to compare differences in intergroup data at each time point. Pearson correlation analysis was used to analyze FA and ADC values and BBB scores. A value of P < 0.05 was considered statistically significant.

    Results

    General condition of the experimental animals

    All 72 rats recovered from anesthesia within 2 hours of surgery and survived over the course of the experimental period. All 72 rats were included in thefinal analysis.

    Neurological function assessment

    Lower hindlimb function was assessed with the BBB scale 24 hours before SCI, and at 6 hours and 1, 2, 3 and 4 weeks after SCI. All rats were evaluated on schedule and received 21 points before SCI. The rats in the SCI, lithium and 3-MA groups displayedflaccid paralysis and a failure of autonomic urination. Neurological function in the sham group was the same as in the pre-operative period at all time points after laminectomy. The graph shows the BBB scores at the different time points. Locomotor function was dramatically reduced after SCI and gradually improved with time in the SCI, lithium and 3-MA groups (P < 0.05; Figure 1). Recovery was significantly better in the lithium group than in the SCI and 3-MA groups at 1, 2, 3 and 4 weeks after SCI (P <0.05). There was no difference in BBB scores between the 3-MA and SCI groups at 1 and 2 weeks after SCI (P > 0.05),but BBB scores were higher in the 3-MA group than in the SCI group at 3 and 4 weeks (P < 0.05).

    Changes in FA and ADC values at the injury site

    FA values significantly decreased (P < 0.05), while ADC values increased significantly (P < 0.05) after SCI. FA values gradually increased over time in all groups (P < 0.05),and there was no difference among the three groups at 6 hours after SCI (P > 0.05). At 1 week after SCI, the FA value was higher in the lithium group than in the SCI and 3-MA groups (P < 0.05), and there was no difference between the SCI and 3-MA groups (P > 0.05). At 2, 3 and 4 weeks after SCI, the FA value was higher in the lithium group than in the SCI and 3-MA groups (P < 0.05), and higher in the 3-MA group than in the SCI group (P < 0.05; Table 1).

    Figure 1 Effects of lithium and 3-MA on motor function in rats with SCI.

    ADC values gradually decreased with time in all groups (P< 0.05), and there was no difference among the three groups at 6 hours or 1 week after SCI (P > 0.05). At 2 weeks after SCI, the ADC value was lower in the lithium group than in the SCI and 3-MA groups (P < 0.05), and there was no difference between the SCI and 3-MA groups (P > 0.05). At 3 and 4 weeks after SCI, the ADC value was again lower in the lithium group than in the SCI and 3-MA groups (P < 0.05),and was lower in the 3-MA group than in the SCI group (P< 0.05; Table 2).

    Correlation between DTI and neurological function

    Pearson correlation analysis showed that FA values were negatively correlated with ADC values in the rat model of spinal cord contusion injury (r = ?0.9537, P < 0.05; Figure 2A), consistent with our previous observations (Zhang et al., 2015; He, 2015). FA values were positively and linearly correlated with BBB scores (r = 0.9279, P < 0.05; Figure 2B).ADC values were negatively correlated with BBB scores, and the correlation was linear (r = ?0.9173, P < 0.05; Figure 2C).

    Immunolabeling for neurons

    Figure 2 Correlation between diffusion tensor imaging and neurological function assessment (Pearson correlation analysis).

    Table 1 FA values at the injury site at different time points after SCI

    Table 2 ADC values at the injury site at different time points after SCI

    Immunohistochemical staining showed that the number of neurons (NeuN+cells) at the site of injury was reduced in the SCI group at 6 hours after SCI, and continued to diminish with time compared with the sham group. In the lithium group, neurons were similarly reduced at 6 hours after SCI but were more numerous than in the SCI and 3-MA groups at 1 and 3 days and 1 week after SCI. In the 3-MA group,neurons were greatly reduced in number compared with the SCI group at 6 hours after SCI and compared with the SCI and lithium groups at 1 and 3 days and 1 week after SCI(Figure 3A). Furthermore, these neurons had an abnormal morphology.

    Cell counting showed that the number of neurons in all three experimental groups decreased significantly at 6 hours after SCI compared with the sham group (P < 0.05). The number of neurons in the experimental groups continued to decrease at 1 and 3 days and 1 week after SCI, compared with the previous time point (P < 0.05). More neuronal cells survived in the lithium group than in the SCI group, and more neuronal cells survived in the SCI group than in the 3-MA group (P < 0.05; Figure 3B).

    Beclin-1 immunohistochemistry

    The Beclin-1+area was larger and more strongly stained,and Beclin-1+cells were more numerous in the SCI group at 6 hours after SCI compared with the sham group. Further expansion of the Beclin-1+area was found at 1 day, but it diminished from 3 days after SCI. In the lithium group, the Beclin-1+area was larger and more intensely stained, and Beclin-1+cells were more numerous at 6 hours after SCI compared with the SCI group. The Beclin-1+area was even larger at 1 day, but it started to diminish from 3 days after SCI, although the staining was still more intense than in the SCI group at 1 week. In the 3-MA group, the staining was slightly more intense than in sham group 1 day after SCI,while it was weaker than in the SCI group at 3 days and 1 week after SCI (Figure 4A).

    The relative Beclin-1+area in the SCI group at 6 hours after SCI was significantly greater than 1 (P < 0.05), indicating that it was larger than in the sham group and that the level of autophagy increased after SCI. The relative Beclin-1+area reached a peak at 1 day after SCI and decreased from 3 days after SCI. There was a similar trend in the lithium group,with the relative area increasing from 6 hours after SCI,peaking at 1 day, and decreasing significantly from 3 days (P< 0.05). Compared with the SCI group, the relative Beclin-1+area in the lithium group was greater at 6 hours, 1 and 3 days and 1 week after SCI (P < 0.05). In comparison, the relative Beclin-1+area was significantly smaller in the 3-MA group than in the SCI group at 6 hours and 1 and 3 days after SCI (P < 0.05; Figure 4B).

    LC3B immunohistochemistry

    The LC3B+area was expanded, the staining intensity was greater, and positive cells were increased in the SCI group at 6 hours after SCI compared with the sham group. Further expansion of the LC3B+area was found at 1 day, but it decreased from 3 days after SCI. The LC3B+area was expanded, the staining intensity was greater, and positive cells were increased in the lithium group at 6 hours after SCI compared with the SCI group. Further expansion of the LC3B+area was found at 1 day, and it shrank from 3 days after SCI, although the staining was still more intense than in the SCI group at 1 week after SCI. In the 3-MA group, 1 day after SCI, the staining was slightly stronger than in the sham group, while it was weaker than in the SCI group at 3 days and 1 week after SCI (Figure 5A).

    The quantitative analysis revealed that the relative LC3B+area in the SCI group at 6 hours after SCI was significantly greater than 1 (P < 0.05), indicating that the LC3B+area was larger than in the sham group, and suggesting that the level of autophagy increased after SCI. The relative LC3B+area peaked at 1 day after SCI and decreased from 3 days after SCI. There was a similar trend in the lithium group, with the relative area increasing from 6 hours after SCI, peaking at 1 day, and significantly decreasing from 3 days (P < 0.05).Compared with the SCI group, the relative LC3B+area in the lithium group was greater at 6 hours, 1 and 3 days and 1 week after SCI (P < 0.05). However, the relative LC3B+area was significantly smaller in the 3-MA group than in the SCI group at 6 hours and 1 and 3 days after SCI (P < 0.05; Figure 5B).

    Discussion

    Advanced evaluation of SCI

    Conventional MRI is widely used for patients with SCI.However, it fails to clearly display the degree of injury or the recovery and regeneration of neuronal fibers in the spinal cord after injury. Therefore, in the present study, we used DTI for the three-dimensional reconstruction of white matterfiber bundles.

    Based on our previous study, DTI is an objective and accurate method for evaluating recovery following SCI and the effect of therapeutic interventions in complete transection SCI models (Zhang et al., 2015). SCI causes damage to cell membranes and myelin sheaths, leading to the destruction of the molecular diffusion barrier and the unrestricted movement of water (Li et al., 2016). This occurred immediately after SCI. Subsequently, along with glial scar formation,the displacement of water molecules was reduced, and the regeneration of axons forced the water molecules to diffuse primarily in one direction, which was reflected as a gradual increase in the FA value and a decrease in the ADC value.The DTI outcomes we observed in this study were consistent with the pathological changes. The FA and ADC values correlated well with the BBB scores. Lithium promoted recovery following SCI, while 3-MA reduced the therapeutic effectiveness of lithium. Therefore, DTI can accurately reflect axonal necrosis and degeneration, glial cell regeneration and demyelination after SCI, and display changes in the microstructure of the spinal cord in vivo (Zhang et al., 2015; Jirjis et al., 2017).

    Autophagy in lithium treatment for SCI

    Figure 3 Immunohistochemical staining and counting of neurons (NeuN+ cells) in the injured rat spinal cord at different time points.

    Figure 4 Immunohistochemical staining for Beclin-1 and the relative Beclin-1+ area in the injured spinal cord at different time points.

    Figure 5 Immunohistochemical staining for LC3B and the relative LC3B+ area in the injured spinal cord at different time points.

    Autophagy is an evolutionarily conserved process, and over 30 autophagy-related genes (Atgs) have been identified, of which LC3B (or LC3II) and Beclin-1 are standard markers(Kirisako et al., 1999; Ohsumi, 2001; Mizushima and Yoshimori, 2007; Periyasamy-Thandavan et al., 2009). Accumulating evidence suggests that lithium has neuroprotective properties, suggesting that it may have potential as a new therapy for SCI (Sarkar et al., 2005; Wada et al., 2005; Yan et al., 2007;Pasquali et al., 2009). Although lithium has been widely used for safely and effectively treating neuropsychiatric disorders,it is rarely used for acute SCI (Ohsumi, 2001; Chang et al.,2011; Chen et al., 2013; Kim et al., 2013; Duo and He, 2015;Hou et al., 2015; Seo et al., 2015; Quartini et al., 2016; Zhou et al., 2016; Wu et al., 2018). Therefore, the effectiveness of lithium treatment for acute SCI remained unclear.

    The role of autophagy in recovery following injury is still controversial. While some studies have reported that enhanced autophagy improves neuroprotection, others have suggested that the suppression of autophagy is beneficial to recovery (Li et al., 2010; Shimada et al., 2012; O’Donovan et al., 2015; Del Grosso et al., 2016). Previous studies in otherfields have demonstrated that lithium can enhance autophagy, or in contrast, reduce apoptosis and autophagy (Wong et al., 2011; Raja et al., 2015; Guttuso, 2016). Our results show that lithium promotes neurological functional recovery and neural cell survival, which supports the notion that lithium has neuroprotective properties. Furthermore, we observed that these neuroprotective effects were inhibited by 3-MA,which downregulated the autophagy induced by lithium.This implies that lithium reduces neuronal damage and promotes functional recovery by inducing autophagy. Nevertheless, BBB scores were still higher in the 3-MA group than in the SCI group from 3 weeks after SCI, in accordance with the neural cell counting results. This suggests that the mechanisms of autophagy are complex and that other signaling pathways that are not inhibited by 3-MA or activated by lithium play a role in the process (Galluzzi et al., 2017)Furthermore, lithium may also promote neurotrophin secretion, inhibit inflammation or enhance the proliferation of neural progenitor cells (Son et al., 2003; Senatorov et al.,2004; Su et al., 2007; Yasuda et al., 2009; Chi-Tso and Chuang, 2011; Li et al., 2011).

    The opposing concept that autophagy aggravates injury may be explained by differences in lithium concentrations and target cells in previous studies. Discrepancies may also be caused by differences in animal models, the therapeutic window and the treatment period. Fang et al. (2016) found that early activated autophagy alleviates spinal cord ischemia/reperfusion injury, while later excessively elevated au-tophagy aggravates the injury. Therefore, autophagy appears to be a dynamic process with differential effects, depending on the time frame and model. Further study is needed to examine the signaling pathways affected by lithium and the dynamic changes in the autophagy pathway.

    Summary

    Further clinical trials are required to explore the effect of lithium therapy in acute SCI patients. In addition, studies are needed to optimize the time window of treatment, the treatment dose and protocol and to reduce the side effects of lithium.

    In conclusion, our findings demonstrate that lithium protects neurons and promotes autophagy in a rat model of acute SCI. DTI is an effective method for evaluating recovery following SCI, and correlates well with neurological functional scores in our rat model of spinal cord contusion injury. The dynamic changes in autophagy after SCI and the effects of lithium on this process need to be investigated further.

    Acknowledgments:We are very grateful to Jia-Lin Zhu from Xi’an Jiaotong University of China for technical support.

    Author contributions:Data provision and integration: DZ and XHL;study concept and design: DZ and XJH; experimental data analysis: DZ,XZ and XHL; manuscript drafting: DZ; statistical analysis: FW; paper supervisor and instructor: XJH; technical or information support: DZ,XHL, FW and XJH. All authors approved thefinal version of the paper.

    Conflicts of interest:None declared.

    Financial support:This study was supported by the Beijing Excellent Talent Training Funding in China, No. 2017000021469G215 (to DZ);the Youth Science Foundation of Beijing Tiantan Hospital of China, No.2016-YQN-14 (to DZ); the Natural Science Foundation of Capital Medical University of China, No. PYZ2017082 (to DZ); the Xi’an Science and Technology Project in China, No. 2016048SF/YX04(3) (to XHL). The funding sources had no role in study conception and design, data analysis or interpretation, paper writing or deciding to submit this paper for publication.

    Institutional review board statement:All experimental procedures and protocols were approved by the Animal Ethics Committee of Xi’an Jiaotong University of China. All experimental procedures described here were in accordance with the National Institutes of Health (NIH)guidelines for the Care and Use of Laboratory Animals.

    Copyright license agreement:The Copyright License Agreement has been signed by all authors before publication.

    Data sharing statement:Datasets analyzed during the current study are available from the corresponding author on reasonable request.

    Plagiarism check:Checked twice by iThenticate.

    Peer review:Externally peer reviewed.

    Open access statement:This is an open access journal, and articles are distributed under the terms of the Creative Commons Attribution-Non-Commercial-ShareAlike 4.0 License, which allows others to remix, tweak,and build upon the work non-commercially, as long as appropriate credit is given and the new creations are licensed under the identical terms.

    Open peer reviewer:Chang Ho Hwang, University of Ulsan, Republic of Korea.

    Additionalfile:Open peer review report 1.

    欧美一区二区国产精品久久精品| 国产v大片淫在线免费观看| 亚洲欧洲日产国产| 只有这里有精品99| av在线蜜桃| 人妻制服诱惑在线中文字幕| 国产精品久久电影中文字幕| 久久久久性生活片| 好男人在线观看高清免费视频| 国产黄片视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲va在线va天堂va国产| 久久久久久大精品| 久久精品国产清高在天天线| 国产 一区精品| 22中文网久久字幕| 国产高潮美女av| 伦理电影大哥的女人| 一本精品99久久精品77| 国产亚洲欧美98| 99热全是精品| 国产亚洲精品久久久com| 国产视频首页在线观看| 日韩制服骚丝袜av| 久久精品91蜜桃| 在线播放无遮挡| 久久精品国产鲁丝片午夜精品| 国产精品伦人一区二区| 日本色播在线视频| 午夜爱爱视频在线播放| 99国产精品一区二区蜜桃av| 久久久久久国产a免费观看| 国产精品一及| 亚洲真实伦在线观看| 国产精品久久久久久久电影| av在线亚洲专区| 国产大屁股一区二区在线视频| 在线免费观看不下载黄p国产| 性插视频无遮挡在线免费观看| 国产精品三级大全| eeuss影院久久| 午夜福利在线观看免费完整高清在 | 成年av动漫网址| 热99re8久久精品国产| 久久午夜亚洲精品久久| 国产黄片美女视频| 亚洲av中文字字幕乱码综合| 精品人妻一区二区三区麻豆| 亚洲一区高清亚洲精品| 精品人妻熟女av久视频| 国产成人影院久久av| 日本欧美国产在线视频| 国产成人aa在线观看| 午夜激情欧美在线| 国产亚洲av嫩草精品影院| 中国国产av一级| 久久久久久伊人网av| 国产高清激情床上av| 色噜噜av男人的天堂激情| 国产一区二区三区av在线 | 欧美成人a在线观看| 美女内射精品一级片tv| 欧美在线一区亚洲| 综合色丁香网| 国产精品三级大全| 亚洲最大成人av| 亚洲av成人av| 99久久中文字幕三级久久日本| 黄片wwwwww| 亚洲国产精品sss在线观看| 日韩av不卡免费在线播放| 午夜福利在线在线| 最好的美女福利视频网| 久久精品人妻少妇| 在线天堂最新版资源| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久久成人| 欧美变态另类bdsm刘玥| 欧美高清成人免费视频www| 日韩欧美在线乱码| 亚洲精华国产精华液的使用体验 | 国产伦精品一区二区三区视频9| 少妇被粗大猛烈的视频| 久久精品久久久久久噜噜老黄 | 亚洲精品亚洲一区二区| 一进一出抽搐gif免费好疼| 国内少妇人妻偷人精品xxx网站| 国产精品综合久久久久久久免费| 国内精品美女久久久久久| 成人毛片60女人毛片免费| 不卡一级毛片| 国产成人精品久久久久久| 美女脱内裤让男人舔精品视频 | 99久久久亚洲精品蜜臀av| 小说图片视频综合网站| 99久久九九国产精品国产免费| 亚洲欧美成人精品一区二区| 国产在线男女| 午夜视频国产福利| 国产成人午夜福利电影在线观看| 国产又黄又爽又无遮挡在线| 又黄又爽又刺激的免费视频.| 高清在线视频一区二区三区 | 18禁在线播放成人免费| 有码 亚洲区| 久久久久久久午夜电影| 精品久久久久久久久av| 国产男人的电影天堂91| 成年女人永久免费观看视频| 欧美性猛交黑人性爽| 国产精品爽爽va在线观看网站| 国产精品女同一区二区软件| 亚洲av熟女| 免费人成视频x8x8入口观看| av视频在线观看入口| 一本久久中文字幕| 一本一本综合久久| 亚洲欧美日韩高清专用| 亚洲最大成人手机在线| 1024手机看黄色片| 久99久视频精品免费| 亚洲va在线va天堂va国产| 亚洲精品影视一区二区三区av| 毛片一级片免费看久久久久| 性插视频无遮挡在线免费观看| 搡老妇女老女人老熟妇| 亚洲五月天丁香| 日韩欧美精品免费久久| 插阴视频在线观看视频| 日本欧美国产在线视频| 亚洲欧美日韩高清专用| 一级黄片播放器| 我要看日韩黄色一级片| av在线天堂中文字幕| 国产精品嫩草影院av在线观看| 好男人视频免费观看在线| 亚洲精品影视一区二区三区av| 欧美xxxx黑人xx丫x性爽| 亚洲欧美中文字幕日韩二区| 国产单亲对白刺激| 中文字幕免费在线视频6| 99热这里只有精品一区| 亚洲天堂国产精品一区在线| 国产日韩欧美在线精品| 啦啦啦韩国在线观看视频| 99精品在免费线老司机午夜| 欧美日韩国产亚洲二区| 欧美+日韩+精品| 国产精品一区二区在线观看99 | 97人妻精品一区二区三区麻豆| 亚洲美女视频黄频| 久久精品国产清高在天天线| 97在线视频观看| 18禁在线播放成人免费| 欧美xxxx黑人xx丫x性爽| 国产精品嫩草影院av在线观看| 麻豆国产av国片精品| 国产精品久久久久久精品电影小说 | 热99re8久久精品国产| 久久午夜福利片| 成人毛片60女人毛片免费| 内射极品少妇av片p| 国产大屁股一区二区在线视频| 国产成人a区在线观看| 91精品国产九色| 乱码一卡2卡4卡精品| 午夜免费男女啪啪视频观看| 国产精品99久久久久久久久| av福利片在线观看| 亚洲av中文av极速乱| 少妇被粗大猛烈的视频| 日韩制服骚丝袜av| 淫秽高清视频在线观看| 岛国在线免费视频观看| 精品欧美国产一区二区三| 三级毛片av免费| 精品日产1卡2卡| 中文字幕久久专区| 日日摸夜夜添夜夜添av毛片| 色综合色国产| 精品人妻偷拍中文字幕| 国产精品福利在线免费观看| 国产精品一及| 乱人视频在线观看| 日韩高清综合在线| 国产精品一二三区在线看| 高清在线视频一区二区三区 | 91久久精品国产一区二区成人| 亚洲精品影视一区二区三区av| 国产精品美女特级片免费视频播放器| 如何舔出高潮| 国产在线精品亚洲第一网站| 国产精品野战在线观看| 岛国毛片在线播放| 边亲边吃奶的免费视频| av国产免费在线观看| 日韩三级伦理在线观看| av在线天堂中文字幕| 岛国在线免费视频观看| 国产真实乱freesex| 欧美高清性xxxxhd video| 91精品一卡2卡3卡4卡| av.在线天堂| 久久国产乱子免费精品| kizo精华| a级一级毛片免费在线观看| 亚洲最大成人中文| 久久久久久九九精品二区国产| 久久久成人免费电影| 亚洲精品国产成人久久av| 色播亚洲综合网| 秋霞在线观看毛片| 国产av一区在线观看免费| 一级毛片aaaaaa免费看小| 色综合色国产| 国产免费一级a男人的天堂| 99国产极品粉嫩在线观看| 美女黄网站色视频| 日韩欧美 国产精品| 久久久午夜欧美精品| 亚洲成人中文字幕在线播放| 99在线人妻在线中文字幕| 国产黄片视频在线免费观看| 在现免费观看毛片| 亚洲欧美精品综合久久99| 高清在线视频一区二区三区 | 久久久精品94久久精品| 在线观看美女被高潮喷水网站| 久久韩国三级中文字幕| 久久久久九九精品影院| av女优亚洲男人天堂| 啦啦啦啦在线视频资源| 最近最新中文字幕大全电影3| 三级男女做爰猛烈吃奶摸视频| 插逼视频在线观看| av天堂在线播放| 久久婷婷人人爽人人干人人爱| 免费电影在线观看免费观看| 国产精品日韩av在线免费观看| 国产高清有码在线观看视频| 能在线免费观看的黄片| 精品久久久久久成人av| 成人午夜精彩视频在线观看| 小说图片视频综合网站| 蜜臀久久99精品久久宅男| 国产成人91sexporn| 最近视频中文字幕2019在线8| 国产成人a∨麻豆精品| 好男人在线观看高清免费视频| 一级黄片播放器| 国产久久久一区二区三区| 18+在线观看网站| 熟女电影av网| 亚洲真实伦在线观看| 日韩欧美三级三区| 亚洲国产欧美人成| 亚洲五月天丁香| 国产亚洲欧美98| 国产又黄又爽又无遮挡在线| 一个人看的www免费观看视频| 黄片无遮挡物在线观看| 久久久精品94久久精品| 精品人妻熟女av久视频| 日韩一区二区视频免费看| 亚洲成人av在线免费| 草草在线视频免费看| 观看美女的网站| 日韩 亚洲 欧美在线| 啦啦啦观看免费观看视频高清| av福利片在线观看| 国产一区二区在线观看日韩| 舔av片在线| 最近2019中文字幕mv第一页| 国产av麻豆久久久久久久| 日本-黄色视频高清免费观看| 国产精品av视频在线免费观看| 亚洲av一区综合| 成人二区视频| 久久久久免费精品人妻一区二区| 亚洲精品乱码久久久久久按摩| 国产亚洲91精品色在线| 国内精品宾馆在线| 天堂影院成人在线观看| 精品少妇黑人巨大在线播放 | 亚洲成人av在线免费| 亚洲国产精品sss在线观看| 免费观看人在逋| 欧美3d第一页| 久久精品国产亚洲av涩爱 | 在线免费十八禁| 亚洲欧美日韩高清在线视频| 夫妻性生交免费视频一级片| 一区福利在线观看| 99国产精品一区二区蜜桃av| 国产v大片淫在线免费观看| 亚洲在久久综合| 91av网一区二区| 91在线精品国自产拍蜜月| 啦啦啦啦在线视频资源| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区| 一级毛片电影观看 | 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 成熟少妇高潮喷水视频| 看免费成人av毛片| avwww免费| 久久6这里有精品| 18禁裸乳无遮挡免费网站照片| 欧美高清性xxxxhd video| 国产毛片a区久久久久| 最新中文字幕久久久久| 男女视频在线观看网站免费| 亚洲av电影不卡..在线观看| 欧美日韩国产亚洲二区| 少妇人妻一区二区三区视频| 欧美潮喷喷水| 亚洲中文字幕日韩| 丰满的人妻完整版| 嫩草影院新地址| 欧洲精品卡2卡3卡4卡5卡区| 国产精品99久久久久久久久| 国产精品久久久久久亚洲av鲁大| 久久久精品欧美日韩精品| 一级黄片播放器| 亚洲国产日韩欧美精品在线观看| 人人妻人人澡欧美一区二区| 一级二级三级毛片免费看| 少妇熟女欧美另类| 色5月婷婷丁香| 哪个播放器可以免费观看大片| 岛国在线免费视频观看| 成人亚洲精品av一区二区| 婷婷精品国产亚洲av| 久久这里只有精品中国| 亚洲高清免费不卡视频| 尤物成人国产欧美一区二区三区| 女同久久另类99精品国产91| 美女国产视频在线观看| 精品国产三级普通话版| 久久久精品94久久精品| 狂野欧美白嫩少妇大欣赏| 久久久欧美国产精品| 日本五十路高清| 超碰av人人做人人爽久久| 国产精品蜜桃在线观看 | 男人舔奶头视频| 一级二级三级毛片免费看| 亚洲国产精品久久男人天堂| 日本黄大片高清| АⅤ资源中文在线天堂| 成人av在线播放网站| 在线国产一区二区在线| 只有这里有精品99| 亚洲国产欧美人成| 少妇熟女欧美另类| 日本熟妇午夜| 亚洲美女视频黄频| 最近手机中文字幕大全| 黄色欧美视频在线观看| 亚洲av中文av极速乱| 精品久久久久久久末码| av专区在线播放| 男女做爰动态图高潮gif福利片| 能在线免费观看的黄片| 99久久九九国产精品国产免费| 亚洲综合色惰| 久99久视频精品免费| 亚洲无线观看免费| 99久久中文字幕三级久久日本| 三级毛片av免费| 麻豆成人av视频| 色综合亚洲欧美另类图片| 国产麻豆成人av免费视频| 伦理电影大哥的女人| 国产精品,欧美在线| 中文字幕精品亚洲无线码一区| 成人高潮视频无遮挡免费网站| 97超视频在线观看视频| 亚洲精品456在线播放app| 国产精品1区2区在线观看.| 九色成人免费人妻av| 欧美色欧美亚洲另类二区| 日日摸夜夜添夜夜添av毛片| 久久精品国产自在天天线| 国产亚洲av嫩草精品影院| 国产精品野战在线观看| av免费观看日本| 老司机影院成人| 两性午夜刺激爽爽歪歪视频在线观看| 色噜噜av男人的天堂激情| 男人和女人高潮做爰伦理| 在线国产一区二区在线| 亚洲精品成人久久久久久| 最近手机中文字幕大全| 看黄色毛片网站| 国产精品av视频在线免费观看| 简卡轻食公司| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 99热这里只有是精品在线观看| 国语自产精品视频在线第100页| 一边亲一边摸免费视频| 久久久国产成人精品二区| 熟女电影av网| 亚洲最大成人av| 亚洲国产精品久久男人天堂| 国产乱人视频| 久久韩国三级中文字幕| 亚洲人成网站高清观看| 日韩欧美精品v在线| 久久久久九九精品影院| 色5月婷婷丁香| av免费在线看不卡| 精品一区二区三区视频在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚州av有码| 国产一级毛片在线| 寂寞人妻少妇视频99o| 乱人视频在线观看| 亚洲国产精品久久男人天堂| 日韩精品青青久久久久久| 亚洲五月天丁香| 日韩在线高清观看一区二区三区| 99久久精品热视频| 亚洲人与动物交配视频| 看免费成人av毛片| 国产精品一区二区三区四区久久| 啦啦啦啦在线视频资源| 国产亚洲精品久久久久久毛片| 青春草亚洲视频在线观看| 国产又黄又爽又无遮挡在线| 神马国产精品三级电影在线观看| videossex国产| 国产精品人妻久久久久久| 热99在线观看视频| 少妇熟女欧美另类| 欧美最新免费一区二区三区| 亚洲自偷自拍三级| 嘟嘟电影网在线观看| 国产在线男女| 亚洲精品成人久久久久久| 久久久久网色| 特级一级黄色大片| 哪里可以看免费的av片| 免费黄网站久久成人精品| 男的添女的下面高潮视频| 国产精品久久久久久av不卡| 尾随美女入室| 我的女老师完整版在线观看| 秋霞在线观看毛片| 欧美日韩一区二区视频在线观看视频在线 | 好男人视频免费观看在线| 亚洲在线自拍视频| 成人国产麻豆网| 午夜激情欧美在线| 在线天堂最新版资源| 免费大片18禁| 99久久精品一区二区三区| 日韩精品有码人妻一区| 韩国av在线不卡| 人妻夜夜爽99麻豆av| 性插视频无遮挡在线免费观看| 亚洲无线观看免费| 成年免费大片在线观看| 日日干狠狠操夜夜爽| 国产色爽女视频免费观看| 亚洲国产精品sss在线观看| 人妻少妇偷人精品九色| 人妻制服诱惑在线中文字幕| 亚洲精品粉嫩美女一区| 国产久久久一区二区三区| 国产精品福利在线免费观看| 乱系列少妇在线播放| 亚洲人成网站在线播放欧美日韩| 在线观看美女被高潮喷水网站| 国产伦精品一区二区三区四那| 国产精品人妻久久久影院| 中文字幕人妻熟人妻熟丝袜美| 国产精品福利在线免费观看| 简卡轻食公司| 欧美三级亚洲精品| 日韩成人av中文字幕在线观看| 亚洲欧美日韩高清在线视频| avwww免费| 国产精品一二三区在线看| 女同久久另类99精品国产91| 观看美女的网站| 日韩在线高清观看一区二区三区| 亚洲国产欧洲综合997久久,| 直男gayav资源| 国产午夜精品久久久久久一区二区三区| 日本欧美国产在线视频| 成人亚洲欧美一区二区av| 一夜夜www| 亚洲精华国产精华液的使用体验 | 嘟嘟电影网在线观看| 亚洲精品久久久久久婷婷小说 | 成人特级av手机在线观看| 日本一本二区三区精品| 久久99热这里只有精品18| 特大巨黑吊av在线直播| 久久人人精品亚洲av| 少妇熟女欧美另类| 三级国产精品欧美在线观看| 欧美zozozo另类| 亚洲性久久影院| 淫秽高清视频在线观看| 午夜久久久久精精品| 一区福利在线观看| 国产v大片淫在线免费观看| 深夜a级毛片| a级毛色黄片| av天堂在线播放| 乱系列少妇在线播放| 精品一区二区三区人妻视频| 美女cb高潮喷水在线观看| 成人亚洲欧美一区二区av| 亚洲婷婷狠狠爱综合网| www.av在线官网国产| 日本免费a在线| 久久久久免费精品人妻一区二区| 九九热线精品视视频播放| 亚洲丝袜综合中文字幕| 国产成人精品婷婷| 老女人水多毛片| 免费不卡的大黄色大毛片视频在线观看 | 国产精品1区2区在线观看.| 亚洲av一区综合| 国产伦在线观看视频一区| 深爱激情五月婷婷| 国产亚洲精品av在线| 91久久精品国产一区二区三区| 啦啦啦韩国在线观看视频| 99在线人妻在线中文字幕| 日韩大尺度精品在线看网址| 欧美精品国产亚洲| 久久久a久久爽久久v久久| 久久99热这里只有精品18| 直男gayav资源| 男女做爰动态图高潮gif福利片| 国产精品久久电影中文字幕| 最新中文字幕久久久久| 特大巨黑吊av在线直播| avwww免费| 亚洲图色成人| avwww免费| 亚洲第一电影网av| 99久久中文字幕三级久久日本| 中文字幕人妻熟人妻熟丝袜美| 久久久久久伊人网av| 丝袜喷水一区| 国产精品一区二区三区四区久久| 99热这里只有是精品50| 少妇的逼好多水| 精品久久久久久久人妻蜜臀av| 欧美一区二区国产精品久久精品| 欧美潮喷喷水| 亚洲欧洲日产国产| 日韩在线高清观看一区二区三区| 久久久精品欧美日韩精品| 亚洲性久久影院| 亚洲av一区综合| 日韩国内少妇激情av| 岛国在线免费视频观看| 白带黄色成豆腐渣| 亚洲五月天丁香| 中文亚洲av片在线观看爽| 日韩亚洲欧美综合| 欧美一区二区精品小视频在线| 久99久视频精品免费| 大型黄色视频在线免费观看| 国产麻豆成人av免费视频| 亚洲精品久久久久久婷婷小说 | 国产精品女同一区二区软件| 97超碰精品成人国产| 免费看美女性在线毛片视频| 特级一级黄色大片| 一夜夜www| 亚洲色图av天堂| 在线观看午夜福利视频| 国产老妇女一区| 日韩亚洲欧美综合| 欧美一区二区精品小视频在线| 熟女人妻精品中文字幕| 精品人妻偷拍中文字幕| 热99re8久久精品国产| 亚洲真实伦在线观看| 久久久国产成人精品二区| 看非洲黑人一级黄片| 欧洲精品卡2卡3卡4卡5卡区| av免费在线看不卡| 成人综合一区亚洲| 成年女人永久免费观看视频| 久久久久久久久中文| 欧美性感艳星| 国产精品乱码一区二三区的特点| 欧美zozozo另类| 亚洲人成网站在线播放欧美日韩| 在现免费观看毛片| 色哟哟哟哟哟哟| 最近的中文字幕免费完整| av在线观看视频网站免费| 欧美丝袜亚洲另类| 日韩人妻高清精品专区| 人妻少妇偷人精品九色| 亚洲真实伦在线观看| 中出人妻视频一区二区| 成人国产麻豆网| 国产 一区 欧美 日韩| 长腿黑丝高跟| 日韩在线高清观看一区二区三区| 日本爱情动作片www.在线观看| 看十八女毛片水多多多| 男人的好看免费观看在线视频|