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

    LncRNA Airsci increases the inflammatory response after spinal cord injury in rats through the nuclear factor kappa B signaling pathway

    2021-11-02 02:59:10TaoZhangKangLiZiLuZhangKaiGaoChaoLiangLv

    Tao Zhang , Kang Li, Zi-Lu Zhang Kai Gao, , Chao-Liang Lv,

    Abstract Spinal cord injury (SCI) is a serious traumatic event to the central nervous system. Studies show that long non-coding RNAs (lncRNAs)play an important role in regulating the inflammatory response in the acute stage of SCI. Here, we investigated a new lncRNA related to spinal cord injury and acute inflammation. We analyzed the expression profile of lncRNAs after SCI, and explored the role of lncRNA Airsci(acute inflammatory response in SCI) on recovery following acute SCI. The rats were divided into the control group, SCI group, and SCI +lncRNA Airsci-siRNA group. The expression of inflammatory factors, including nuclear factor kappa B [NF-κB (p65)], NF-κB inhibitor IκBα and phosphorylated IκBα (p-IκBα), and the p-IκBα/IκBα ratio were examined 1–28 days after SCI in rats by western blot assay. The differential lncRNA expression profile after SCI was assessed by RNA sequencing. The differentially expressed lncRNAs were analyzed by bioinformatics technology. The differentially expressed lncRNA Airsci, which is involved in NF-κB signaling and associated with the acute inflammatory response, was verified by quantitative real-time PCR. Interleukin (IL-1β), IL-6 and tumor necrosis factor (TNF-α) at 3 days after SCI were measured by western blot assay and quantitative real-time PCR. The histopathology of the spinal cord was evaluated by hematoxylin-eosin and Nissl staining. Motor function was assessed with the Basso, Beattie and Bresnahan Locomotor Rating Scale. Numerous differentially expressed lncRNAs were detected after SCI, including 151 that were upregulated and 186 that were downregulated in the SCI 3 d group compared with the control group. LncRNA Airsci was the most significantly expressed among the five lncRNAs involved in the NF-κB signaling pathway. LncRNA Airsci-siRNA reduced the inflammatory response by inhibiting the NF-κB signaling pathway, alleviated spinal cord tissue injury, and promoted the recovery of motor function in SCI rats. These findings show that numerous lncRNAs are differentially expressed following SCI, and that inhibiting lncRNA Airsci reduces the inflammatory response through the NF-κB signaling pathway, thereby promoting functional recovery. All experimental procedures and protocols were approved by the approved by the Animal Ethics Committee of Jining Medical University (approval No. JNMC-2020-DW-RM-003) on January 18, 2020.

    Key Words: apoptosis; functional recovery; inflammatory response; long non-coding RNA; neuroprotection; NF-κB signaling pathway; RNA sequencing; spinal cord injury

    Introduction

    Spinal cord injury (SCI) can lead to severe dysfunction of the limbs below the damaged segment, which not only affects the physical and mental health and quality of life of patients, but also imparts a huge economic burden on families and society(Kumar et al., 2018; Jung et al., 2019). The pathogenesis of SCI is complex, involving two stages of primary and secondary injury (Cao and Krause, 2020). After the occurrence of primary injuries such as fracture and compression, a series of secondary injuries such as inflammation, hypoxia, neuronal necrosis and apoptosis, and local microenvironmental changes further aggravate the SCI, resulting in severe loss of sensory and motor functions (Dai et al., 2018; Li et al., 2019b; Lv et al.,2019). At present, clinical therapeutics for SCI are limited and the pathogenesis of SCI still needs to be further studied.

    It has been shown that the inflammatory response plays a key role in the pathogenesis of SCI (do Espírito Santo et al.,2019; He et al., 2019; Wang et al., 2019; Oh et al., 2020).Dysregulation of the inflammatory response after SCI causes a cascade effect, inducing apoptosis, edema, oxidative stress,and other reactions that aggravate tissue damage (Gao et al.,2019; Huang et al., 2019). Early anti-inflammatory treatment can improve the clinical prognosis of patients with SCI (Gensel and Zhang, 2015). It has also been found that stem cell transplantation in the acute inflammatory phase can provide better functional recovery (Nishimura et al., 2013; Cheng et al., 2017). Therefore, an in-depth study of the molecular mechanisms of SCI, including the dysregulation of the inflammatory response, may provide new molecular targets for the diagnosis and treatment of SCI.

    Long non-coding RNA (lncRNA) is a type of RNA with a transcript length of over 200 nt and no protein-coding function (Rui et al., 2018). LncRNAs are widely expressed in organisms and have important regulatory functions(McDonel and Guttman, 2019; Qian et al., 2019). The main functions of lncRNAs include participating with competitive endogenous RNAs, interacting with proteins, regulating mRNA shearing, chromatin and histone remodeling, and transcriptional regulation (Shi et al., 2018). Preliminary studies show that lncRNAs play a critical role in the regulation of the molecular mechanisms involved in the acute phase of inflammatory dysregulation in SCI (Shao et al., 2020). LncRNA LOC685699-OT1 is a new transcript of gene LOC685699(ENSRNOG00000043358), located on chromosome X:23419209–23467865. No studies on lncRNA LOC685699-OT1 in SCI have been found to date. Therefore, we named this lncRNA Airsci (acute inflammatory response in SCI) according to the rules of the Hugo Gene Naming Committee (HGNG)naming guidelines. Here, we found that lncRNA Airsci had the most significantly differential expression after SCI in rats, and its expression was significantly up-regulated during the acute inflammatory phase. Hence, exploring the molecular functions of lncRNA Airsci in SCI may provide new therapeutic targets and intervention strategies for treatment and rehabilitation following SCI.

    Materials and Methods

    Animals

    A total of 105 healthy and specific pathogen-free adult(6–7-week-old) male Sprague-Dawley (SD) rats, weighing 250–280 g, were purchased from Jinan PengYue Experimental Animal Breeding Co., Ltd (Jinan, Shandong Province, China).The rats were kept at 25 ± 0.5°C in the SPF Experimental Animal Center of Jining Medical University. The light/dark cycle was altered for 12 hours. To examine the regulatory function of lncRNA Airsci in SCI rats, we divided the rats into the control group, SCI group, and SCI + lncRNA Airsci-siRNA group. Control group: only laminectomy, no SCI; SCI group: no intervention after SCI; SCI + lncRNA Airsci-siRNA group: The pLKO.1-sh Airsci small interfering RNA gene silencing plasmid was constructed. After transformation and identification,the plasmid was amplified using the plasmid extraction kit(Beyotime, Shanghai, China). The plasmid was packaged with lentivirus and transfected in 293T cells according to 750 mL OPTI medium + 1 μg pMD2.G (envelope plasmid)+ 2 μg pCMV-r8.74 (packaging plasmid) + 3 μg pLKO.1-sh Airsci to collect the lentivirus solution (Genechem, Shanghai,China). The filtered and purified lentivirus solution (0.5 mL)was injected into the SCI site to inhibit the expression of lncRNA Airsci. All procedures involving animals concurred with the Animal Research Reporting ofIn VivoExperiments(ARRIVE) guidelines, and were approved by the Animal Ethics Committee of Jining Medical University (approval No. JNMC-2020-DW-RM-003) on January 18, 2020.

    Establishment of the rat SCI model

    The SCI model was established using the modified Allen technique (Lu et al., 2016). Pentobarbital sodium (Germany,Solarbio) was used for intraperitoneal injection anesthesia at a dose of 65 mg/kg. After anesthesia, the skin was prepared and disinfected in the back operation area. A 2–3 cm skin incision was made to separate the fascia and muscle, exposing the spinous process and lamina of T9–11. The ligaments and muscles on both sides of the T9–10 spinous process were detached to expose the T9–10 spinous process and transverse process, and the T9–10 spinous process and lamina were removed to expose the spinal cord fully. Using a spinal cord strike machine, a small 15 g blunt metal stick was allowed to fall freely from a height of 30 mm to produce blunt spinal cord contusion. After successful modeling, the wound was sewn up layer by layer. The postoperative manifestations of the rats were paralysis of the lower limbs and bladder distension. The rats were given bladder massage once in the morning and once in the evening for artificial urination.

    Sample collection

    Rats were sacrificed with 25% CO2asphyxiation on days 1, 2, 3,5, 7, 14, 21 and 28 after SCI. Cardiac perfusion was performed using 10% formaldehyde solution for fixation. The surrounding segmental lamina were removed to fully expose the injured spinal cord tissues. The lesioned spinal cord tissues were cut 6 mm from the top to the bottom of the lesion site, and the tissue was preserved in liquid nitrogen for subsequent western blot assay, qPCR, and RNA-sequencing.

    Western blot assay

    Three days after SCI, the T9–11 spinal cord tissues were dissected out and homogenized in radioimmunoprecipitation assay (RIPA) lysis buffer. The protein concentration of the supernatant was measured, and electrophoreses and membrane transfer were carried out. Membranes were blocked with 5% skim milk for 2 hours at room temperature and incubated overnight at 4°C with the following primary antibodies: anti-nuclear factor kappa-B (NF-κB (p65),1:500, Cat# ab16502, Abcam, Cambridge, MA, USA), antiphosphorylated NF-κB inhibitor (p-IκBα, 1:1000, Cat#ab109300, Abcam), anti-IκBα (1:1000, Cat# ab109300,Abcam), anti-interleukin (IL-1β, 1:500, Cat# ab200478,Abcam), anti-IL-6 (1:500, Cat# ab9324, Abcam), anti-tumor necrosis factor (TNF-α, 1:1000, Cat# ab109322, Abcam),anti-β-actin (1:1000, Cat# ab8226, Abcam). Subsequently,membranes were incubated with peroxidase-labeled secondary antibody (goat anti-rabbit IgG and goat antimouse IgG, 1:2000, Cat# ab6721 and ab6789, Abcam) at room temperature for 2 hours. Images were obtained by continuous exposures with the enhanced chemiluminescence(ECL) imaging system (Tanon Imager, Shanghai, China). The final results were semi-quantitatively analyzed relative to the optical density of β-actin.

    RNA-sequencing

    From the western blot results, we found that 3 days after SCI was the critical time point during the acute inflammatory response period (Figure 1). Therefore, rats were divided into the control group and SCI 3d group for RNA-sequencing. After pentobarbital sodium anesthesia, spinal cord tissues of 3 mm length at each end of the head and tail of the rats were excised(the same parts were taken from the control group) with the injury point as the center, and total RNA was extracted with the Trizol method (Invitrogen, Carlsbad, CA, USA). The RNA samples were tested, and then, a specific library was constructed after removing ribosomal RNA. The library was initially quantified by Qubit 2.0 (Shanghai, China) and diluted to 1.5 ng/μL. Then, the insert size (length of the sequence between the pair adapters) was analyzed on an Agilent 2100(Shanghai, China). Once the insert size was satisfactory, the effective concentration of the library (greater than 3 nM) was accurately quantified by qRT-PCR to ensure library quality.Then, Illumina PE150 sequencing was conducted after pooling of the library according to the coverage required (sequencing was conducted by Shanghai Jikai Gene Chemical Technology Co., Ltd., Shanghai, China).

    Screening of lncRNAs by bioinformatics technology

    We used RNA-sequencing to detect changes in lncRNA expression in the spinal cord tissues of rats in the control group and the SCI 3d group, and obtained the expression profile heat map of lncRNAs. Then, bioinformatics technology was used to screen the differential lncRNAs in the SCI 3d group that were related to the inflammatory response.Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.kegg.jp/) is a database for the systematic analysis of relationships, gene function and gene information between genes and their encoded products. To identify lncRNAs associated with acute inflammation, we first screened lncRNAs according to the set fold change (over 10 times) and the significance level of expression difference (< 0.05) between the SCI 3 days and control groups. Next, target gene prediction and KEGG function enrichment analysis were performed.Finally, the lincRNA Airsci, which is involved in the NF-κB signaling pathway and associated with the acute inflammatory response, was selected for further analysis.

    Quantitative real-time PCR

    Quantitative real-time PCR (qRT-PCR) was used to measure the expression levels of related genes in each group 3 days after SCI. Total RNA in tissues (40–50 mg) was extracted by the Trizol

    method (Invitrogen). The purity and quantity of RNA were determined with an ultraviolet spectrophotometer (SYNERGY H1, BioTek, USA) and then diluted with diethylpyrocarbonate(DEPC) water (Beyotime, Shanghai, China). cDNA was then synthesized according to the instructions in the TaKaRa RNA PCR Kit (TaKaRa, Otsu, Japan). qRT-PCR parameters were as follows: 95°C for 5 minutes, followed by 40 cycles

    (amplification) of 95°C for 15 seconds (denaturation), 60°C for

    40 seconds (annealing), and 72°C for 20 seconds (extension).Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a reference gene, and the relative quantification of gene expression levels was performed using the 2–??Ctmethod. The sequences of all primers used for the qRT-PCR experiments are presented in Table 1.

    Hematoxylin-eosin and Nissl staining

    Hematoxylin-eosin (HE) and Nissl staining were used to evaluate the histopathology of the spinal cord in each group at 3 days after SCI. The collected tissues were paraffinembedded, sectioned into 10 μm sections, and then HE and Nissl staining were performed. HE staining: The sections were stained with hematoxylin for 1 minute, and washed with water. Then, they were immersed in 1% hydrochloric acid/ethanol for 5 seconds and rinsed with water. Sections were then treated with ammonia water for 5 minutes and rinsed with water. Eosin staining was conducted for 1 min and rinsed with water. Sections were then dehydrated sequentially with 80%, 95% and anhydrous ethanol for 10 seconds each.Nissl: The sections were immersed overnight in an equal volume of anhydrous ethanol and chloroform solution. On the next day, they were immersed in anhydrous ethanol, 95%ethanol, and distilled water. The sections were then placed in 0.1% methylphenol solution at 40°C for 10 minutes and washed quickly with distilled water. We observed the tissue morphology of HE staining sections with an optical microscope(OLYMPUS, Jinan, China) and measured the proportional spinal cord tissue lesion size with ImageJ software (NIH,Bethesda, USA). We observed the neurons in the anterior horn of the spinal cord and performed count analysis with the optical microscope (five animals were included; counting was performed at 400× magnification; results were expressed as averages).

    Basso, Beattie and Bresnahan Locomotor Rating Scale assessment

    The motor function of rat hind limbs was measured by the Basso, Beattie and Bresnahan Locomotor Rating Scale (BBB)score (All et al., 2020) at different time points (0, 1, 3, 7,14, 21, 28 days) after SCI. BBB scores ranged from 0 to 21 points. The minimum number of points (0) indicates complete paralysis and the maximum number of points (21) represents normal function. Three uninformed researchers assessed the motor function of rats. As the score increased, the motor function improved.

    Statistical analysis

    Data were analyzed using GraphPad Prism version 5.0 for Windows (GraphPad Software, La Jolla, CA, USA) and expressed as the mean ± SD. Western blot, qPCR, HE and Nissl staining differences among multiple groups were analyzed with oneway analysis of variance and the least significant difference test. BBB scores were analyzed using the Mann-WhitneyU-test.P< 0.05 was considered statistically significant.

    Results

    The NF-κB signaling pathway is activated after SCI in rats

    Western blot results showed that NF-κB (p65) and p-IκBα protein expression increased initially and then decreased between 1 and 28 days after SCI (Figure 1). The p-IκBα/IκBα ratio initially increased and reached a peak at 3 days after SCI,and thereafter decreased (Figure 1).

    lncRNAs are differentially expressed after SCI in rats

    RNA-sequencing results showed that lncRNAs were differentially expressed at 3 days after SCI (Figure 2A).Compared with the control group, 151 lncRNAs were upregulated and 186 lncRNAs were downregulated in the SCI 3d group (Figure 2B).

    Screening of lncRNA Airsci and examining its effect on the NF-κB signaling pathway

    KEGG functional enrichment analysis of lncRNA expression profiles at 3 days after SCI revealed a variety of signaling pathways, including the NF-κB signaling pathway (Figure 3A).Bioinformatics was used to screen for lncRNAs involved in the NF-κB signaling pathway and associated with the acute inflammatory response. The screen identified the following five lncRNAs: LOC685699-OT1 (Airsci), AABR07015057.1-OT11,LINC5298, AABR07000398.1-OT18 and AABR07055878.1-OT1. All five lncRNAs were verified by qRT-PCR, and the differences were found to be statistically significant (P< 0.05),consistent with the results of RNA-sequencing. Among these,the expression of lncRNA Airsci increased significantly after SCI (Figure 3B). Accordingly, we divided the rats into the control group, SCI group and SCI + lncRNA Airsci-siRNA group,and western blot was used to detect the expression levels of NF-κB (p65), phospho-IκBα and IκBα (Figure 3C). Inhibition of lncRNA Airsci significantly reduced the protein expression levels of NF-κB (p65) and p-IκBα, and reduced the p-IκBα/IκBα ratio (Figure 3D–G).

    Inhibiting lncRNA Airsci reduces the inflammatory response after SCI in rats

    Western blot and qRT-PCR results showed that the expression levels of IL-1β, IL-6 and TNF-α were dramatically increased in the SCI group at 3 days after SCI compared with the control group. Inhibition of lncRNA Airsci significantly reduced the protein and mRNA expression levels of the inflammatory factors IL-1β, IL-6 and TNF-α (Figure 4).

    Inhibiting lncRNA Airsci ameliorates histopathological changes in the spinal cord and promotes the recovery of motor function in rats

    HE staining was performed to assess the histopathology of the spinal cord tissues 3 days after SCI (Figure 5A). The proportional lesion size was larger in the SCI group compared with the control group. In comparison, the proportional lesion size was smaller in the Airsci-siRNA group than in the SCI group (Figure 5B). Nissl staining was used for counting spinal cord motor neurons and for detecting neuronal injury (Figure 5A). The number of Nissl-positive cells in the anterior horn of the spinal cord was less in the SCI group compared with the control group. Airsci-siRNA dramatically increased the number of Nissl-positive cells compared with the SCI group (Figure 5C). The BBB score was used to assess motor function after SCI (Figure 5D). The BBB score in the SCI group decreased significantly compared with the control group. Notably, the score was markedly higher in the Airsci-siRNA group than in the SCI group at 7 days after SCI.

    Figure 1 |NF-κB signaling pathway expression before and after SCI.

    Figure 2 |LncRNA expression profile after SCI (RNA-sequencing).

    Figure 3 |Screening of lncRNA Airsci 3 days after SCI and regulation of the NF-κB signaling pathway.

    Figure 4|LncRNA Airsci regulates the inflammatory response at 3 days after SCI.

    Figure 5 |LncRNA Airsci regulates histopathological changes in spinal cord tissue and motor function in rats.

    Table 1 |Quantitative real-time PCR primers

    Discussion

    SCI is a serious traumatic event to the central nervous system(Kleene et al., 2019). The inflammatory response plays an important role in the pathophysiology of SCI and affects recovery and functional outcome (Lemmens et al., 2019;?aker et al., 2019). Studies have shown that lncRNAs play an important role in regulating the inflammatory response in the acute stage of SCI (Zhou et al., 2018; Jia et al., 2019).Therefore, further study of the molecular mechanisms underpinning lncRNA function may provide new ideas and methods for the treatment of SCI.

    In this study, our major findings are as follows: (1) The NFκB signaling pathway is significantly activated after SCI, and the inflammatory response is significantly increased, with the most significant increase at 3 days after SCI; (2) The expression of lncRNA after SCI is significantly different and plays a pivotal role in the pathological mechanism of SCI; (3) LncRNA Airsci is involved in the regulation of the NF-κB inflammatory signaling pathway; (4) Inhibiting lncRNA Airsci reduces the inflammatory response, lessens the area of spinal cord tissue damage,reduces the death of neurons, and promotes the recovery of motor function.

    Current studies demonstrate that lncRNAs are significantly differentially expressed after SCI, and play an important regulatory role in the process of tissue repair (Zhang et al.,2018; Ren et al., 2019). For example, lncRNA MALAT1 inhibits autophagy and apoptosis of neurons and promotes SCI repair by regulating the Wnt/β-catenin and PI3K/Akt signaling pathways (Liu et al., 2019).In vitrostudies of PC12 cells show that lncRNA Mirt2 inhibits the p38MAPK signaling pathway by downregulating miRNA-429 and reducing apoptosis (Li et al., 2019a). Unfortunately, these studies on lncRNA in SCI are mainly focused on neuronal proliferation, apoptosis and autophagy, and few studies have been conducted on the regulatory function of lncRNAs in the acute phase of SCI (Ren et al., 2019).

    Preliminary research has shown that lncRNA MALAT1 promotes the inflammatory response in the SCI acute phase by regulating miRNA-199b, demonstrating that lncRNA plays an important role in the regulation of the molecular mechanisms of the acute inflammatory response (Zhou et al., 2018). However, indepth studies of the expression profile of lncRNAs in the SCI acute phase are lacking. Further study, including additional screening of lncRNAs, may clarify their important regulatory functions in the SCI acute inflammatory response, and may provide new ideas and methods for SCI treatment.

    A modified Allen’s apparatus was used to produce the rat SCI model, and we assessed changes in NF-κB(p65), p-IκBα, IκBα and the p-IκBα/IκBα ratio at 1, 2, 3, 5, 7, 14, 21 and 28 days after SCI. We found that the NF-κB signaling pathway reached peak activation on the third day after SCI. We used RNAsequencing to detect changes in lncRNA expression in spinal cord tissues after SCI, and found that different lincRNAs were expressed. Compared with the control group, 151 lincRNAs were upregulated and 186 lincRNAs were downregulated in the SCI 3d group. To identify lincRNAs involved in the NF-κB signaling pathway and associated with the acute inflammatory response, we screened the lncRNAs through bioinformatics analysis and then verified them by qRT-PCR. Among the five lncRNAs of interest, lncRNA Airsci was significantly upregulated during the SCI acute inflammatory phase. We hypothesized that exploring the role and mechanism of lncRNA Airsci in SCI may provide new therapeutic targets and intervention strategies for the treatment and rehabilitation of SCI. We found that inhibiting lncRNA Airsci expression suppressed the NF-κB signaling pathway, decreased the inflammatory factors IL-1β, IL-6 and TNF-α, and ameliorated spinal cord tissue damage and improved motor function.Therefore, we speculated that inhibiting the expression of lncRNA Airsci could reduce the inflammatory response through NF-κB signaling and promote functional recovery in SCI rats. However, the underlying mechanisms still need to be elucidated, and will be the focus of upcoming studies.

    In summary, the pathological progression of SCI is closely related to the regulation of the inflammatory response, as well as the differential expression of lncRNAs (Ding et al., 2020). We screened lncRNAs by bioinformatics and verified candidates by qRT-PCR. LncRNA Airsci was significantly upregulated during the SCI acute inflammatory phase. Notably, inhibition of lncRNA Airsci reduced the inflammatory response through the NF-κB signaling pathway and promoted functional recovery. Therefore,we speculate that lncRNA Airsci mediates the inflammatory response after SCI in rats through the NF-κB signaling pathway.Our further studies will perform expression analysis on the brain and spinal cord by tissue hybridization to provide insight into the expression profile of Airsci in the CNS in the future.

    Acknowledgments:We thank Dong-Mei Shi from Department of the Key Laboratory of Fungi in the Jining First People’s Hospital, China for providing valuable technical assistance in this work.

    Author contributions:TZ implemented the experiment and wrote the manuscript. KL and ZLZ collected and analyzed the data. KG and CLL designed and reviewed this manuscript. All authors approved the final version of the manuscript.

    Conflicts of interest:The authors declare that they have no conflicts of interest.

    Financial support:This work was supported by the National Natural Science Foundation of China (No. 81801906; to KG) and the Natural Science Foundation of Shandong Province of China (No. ZR2018PH024; to KG).

    Institutional review board statement:All experimental procedures and protocols were approved by the approved by the Animal Ethics Committee of Jining Medical University (approval No. JNMC-2020-DW-RM-003) on January 18, 2020. All experimental procedures described here were in accordancewith the National Institutes of Health (NIH) guidelines for the Care and Use of Laboratory Animals (NIH Publication No. 85-23, revised 1996).

    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-NonCommercial-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.

    午夜福利在线观看免费完整高清在| 日韩精品有码人妻一区| 在线观看免费高清a一片| 久久久亚洲精品成人影院| 久久鲁丝午夜福利片| 特大巨黑吊av在线直播| 99久久综合免费| 国产伦精品一区二区三区视频9| 制服丝袜香蕉在线| 亚洲人成网站在线观看播放| √禁漫天堂资源中文www| 久久狼人影院| 中文字幕精品免费在线观看视频 | 大香蕉久久网| av国产久精品久网站免费入址| 黑人高潮一二区| 久久99一区二区三区| 日韩,欧美,国产一区二区三区| 丝袜在线中文字幕| 最近的中文字幕免费完整| 综合色丁香网| 极品人妻少妇av视频| 欧美日韩成人在线一区二区| 国产高清三级在线| 狂野欧美激情性bbbbbb| 日韩制服骚丝袜av| 夫妻性生交免费视频一级片| 成人二区视频| 国语对白做爰xxxⅹ性视频网站| 老司机影院毛片| 男女高潮啪啪啪动态图| 国产男女内射视频| 精品熟女少妇av免费看| 五月开心婷婷网| 蜜桃久久精品国产亚洲av| 久久久久国产精品人妻一区二区| av一本久久久久| 日韩中文字幕视频在线看片| 亚洲天堂av无毛| 久久鲁丝午夜福利片| 中文字幕久久专区| 精品久久久久久久久亚洲| 一区二区三区精品91| 国产在线一区二区三区精| 天堂俺去俺来也www色官网| 欧美少妇被猛烈插入视频| 天天躁夜夜躁狠狠久久av| 欧美成人午夜免费资源| 日本黄色片子视频| 在线观看国产h片| 中文字幕最新亚洲高清| 爱豆传媒免费全集在线观看| 18在线观看网站| 国产欧美亚洲国产| 乱人伦中国视频| 另类精品久久| 九色亚洲精品在线播放| 一区二区三区乱码不卡18| xxxhd国产人妻xxx| 久久青草综合色| 亚洲色图综合在线观看| 精品少妇久久久久久888优播| 国产一区亚洲一区在线观看| 国产精品无大码| 日本黄大片高清| 国产欧美亚洲国产| 99国产综合亚洲精品| 欧美最新免费一区二区三区| 国产白丝娇喘喷水9色精品| 九九在线视频观看精品| 国产在线视频一区二区| 精品视频人人做人人爽| 国产av一区二区精品久久| 日本91视频免费播放| 久久精品国产亚洲av天美| 久久精品国产亚洲av天美| 五月玫瑰六月丁香| 色婷婷av一区二区三区视频| 亚洲精品久久久久久婷婷小说| 精品酒店卫生间| 看十八女毛片水多多多| 99国产精品免费福利视频| 欧美97在线视频| 欧美一级a爱片免费观看看| 91精品国产九色| 91精品国产九色| 熟妇人妻不卡中文字幕| 特大巨黑吊av在线直播| a级毛色黄片| 婷婷成人精品国产| 婷婷成人精品国产| av黄色大香蕉| 80岁老熟妇乱子伦牲交| 亚洲少妇的诱惑av| 精品久久久噜噜| 欧美一级a爱片免费观看看| 最后的刺客免费高清国语| 美女大奶头黄色视频| 免费黄频网站在线观看国产| 一区二区三区乱码不卡18| 99久国产av精品国产电影| 99久久精品国产国产毛片| 久久 成人 亚洲| 大片免费播放器 马上看| 久久99热这里只频精品6学生| 久久青草综合色| 久久影院123| 好男人视频免费观看在线| 久久综合国产亚洲精品| 久久韩国三级中文字幕| 在现免费观看毛片| 在线观看www视频免费| 午夜精品国产一区二区电影| 成年av动漫网址| 亚洲不卡免费看| 高清在线视频一区二区三区| 午夜福利视频在线观看免费| 国产亚洲av片在线观看秒播厂| 国产淫语在线视频| 日本欧美视频一区| 免费看不卡的av| 日韩伦理黄色片| 欧美日韩精品成人综合77777| 亚洲综合色惰| av又黄又爽大尺度在线免费看| 国产精品免费大片| 免费大片18禁| 欧美性感艳星| av在线播放精品| 热99国产精品久久久久久7| 综合色丁香网| 欧美日韩国产mv在线观看视频| 亚洲国产最新在线播放| 日韩欧美精品免费久久| 亚洲国产色片| 亚洲成人一二三区av| 亚洲人成网站在线播| 乱人伦中国视频| 精品人妻在线不人妻| 在线观看免费视频网站a站| av视频免费观看在线观看| 免费看光身美女| 免费高清在线观看日韩| 男女边摸边吃奶| 国产熟女午夜一区二区三区 | 人妻人人澡人人爽人人| 久久久久国产精品人妻一区二区| 国产视频内射| 久热久热在线精品观看| 国产成人精品一,二区| 熟女av电影| 国产探花极品一区二区| 久久人妻熟女aⅴ| 久久久亚洲精品成人影院| 亚洲天堂av无毛| 爱豆传媒免费全集在线观看| 观看美女的网站| 能在线免费看毛片的网站| 美女国产高潮福利片在线看| 我要看黄色一级片免费的| 最后的刺客免费高清国语| 夜夜骑夜夜射夜夜干| 国产无遮挡羞羞视频在线观看| 精品久久久噜噜| 久久鲁丝午夜福利片| 亚洲av日韩在线播放| 亚洲精品,欧美精品| 久久久久久久久大av| 青春草亚洲视频在线观看| 人妻人人澡人人爽人人| 妹子高潮喷水视频| 成年女人在线观看亚洲视频| 国产探花极品一区二区| 三级国产精品欧美在线观看| 国产一区二区在线观看日韩| 久久久国产一区二区| 在线观看免费高清a一片| 三级国产精品片| 欧美97在线视频| 成人亚洲精品一区在线观看| 蜜臀久久99精品久久宅男| 日韩中字成人| 9色porny在线观看| 亚洲欧美中文字幕日韩二区| 纯流量卡能插随身wifi吗| 亚洲av成人精品一区久久| 亚洲经典国产精华液单| 亚洲美女黄色视频免费看| xxx大片免费视频| 久久国产亚洲av麻豆专区| 一级黄片播放器| 国产无遮挡羞羞视频在线观看| 如何舔出高潮| 秋霞伦理黄片| 黑丝袜美女国产一区| 看十八女毛片水多多多| 久久久久久久久久久久大奶| 免费av不卡在线播放| 18在线观看网站| 五月玫瑰六月丁香| 美女内射精品一级片tv| 亚洲欧美中文字幕日韩二区| 亚洲不卡免费看| 国产成人一区二区在线| 欧美精品人与动牲交sv欧美| 日日撸夜夜添| 日本av手机在线免费观看| 韩国av在线不卡| 99久久精品一区二区三区| 伊人亚洲综合成人网| 五月玫瑰六月丁香| 一区二区三区乱码不卡18| 欧美 日韩 精品 国产| 国产伦理片在线播放av一区| 满18在线观看网站| 哪个播放器可以免费观看大片| 色哟哟·www| 亚洲国产最新在线播放| 热re99久久精品国产66热6| 日韩一区二区视频免费看| 高清午夜精品一区二区三区| 制服诱惑二区| 黑人猛操日本美女一级片| 色婷婷久久久亚洲欧美| 国产深夜福利视频在线观看| 国产欧美亚洲国产| 永久免费av网站大全| 亚洲精品色激情综合| 亚洲国产av新网站| tube8黄色片| 51国产日韩欧美| 欧美日韩成人在线一区二区| 亚洲欧洲日产国产| 2021少妇久久久久久久久久久| av卡一久久| 精品人妻熟女av久视频| 又黄又爽又刺激的免费视频.| 亚洲在久久综合| 一级,二级,三级黄色视频| 亚洲av福利一区| 精品少妇内射三级| a级毛色黄片| av视频免费观看在线观看| 精品熟女少妇av免费看| 99热网站在线观看| 免费大片18禁| 国产成人aa在线观看| 国产日韩欧美在线精品| 久久女婷五月综合色啪小说| 国产午夜精品久久久久久一区二区三区| 国产国拍精品亚洲av在线观看| 高清欧美精品videossex| 日日啪夜夜爽| 国产男女超爽视频在线观看| 街头女战士在线观看网站| 成人影院久久| 欧美日韩av久久| 精品99又大又爽又粗少妇毛片| 岛国毛片在线播放| 国产精品国产三级专区第一集| 国产精品人妻久久久影院| 一区二区三区四区激情视频| 婷婷色av中文字幕| 少妇的逼水好多| 亚洲精品aⅴ在线观看| 午夜激情久久久久久久| 精品国产国语对白av| 亚洲av综合色区一区| 国产成人午夜福利电影在线观看| 我的老师免费观看完整版| 桃花免费在线播放| 性色avwww在线观看| 久久久久久久亚洲中文字幕| 大陆偷拍与自拍| 久久久久久伊人网av| 免费观看的影片在线观看| 高清在线视频一区二区三区| 国产精品国产av在线观看| 夜夜看夜夜爽夜夜摸| 国产精品成人在线| 久久久精品94久久精品| 亚洲欧美日韩另类电影网站| 国产毛片在线视频| 亚洲精品日韩av片在线观看| 亚洲情色 制服丝袜| 久久人人爽人人片av| av专区在线播放| 国产午夜精品一二区理论片| 日韩av不卡免费在线播放| 精品午夜福利在线看| 日韩欧美精品免费久久| 午夜福利视频在线观看免费| 寂寞人妻少妇视频99o| 亚洲av成人精品一二三区| av国产精品久久久久影院| 91精品国产九色| 久久久久国产网址| av福利片在线| 91久久精品国产一区二区成人| 欧美 日韩 精品 国产| 在线观看免费视频网站a站| 久久久久久久久久久久大奶| 一级毛片电影观看| 国产亚洲一区二区精品| 久久久久国产网址| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 黄色欧美视频在线观看| 亚洲国产精品一区二区三区在线| 色哟哟·www| 丰满饥渴人妻一区二区三| 国产成人freesex在线| 国产 一区精品| 日本av手机在线免费观看| 999精品在线视频| 女人久久www免费人成看片| 亚洲精品视频女| 最近最新中文字幕免费大全7| 欧美另类一区| 久久精品久久久久久久性| 日韩成人伦理影院| 成人亚洲欧美一区二区av| 成人黄色视频免费在线看| 成人影院久久| 晚上一个人看的免费电影| 亚洲经典国产精华液单| 五月玫瑰六月丁香| 韩国av在线不卡| 九草在线视频观看| 日韩制服骚丝袜av| a 毛片基地| a级片在线免费高清观看视频| 大话2 男鬼变身卡| 一级a做视频免费观看| 51国产日韩欧美| 丁香六月天网| 全区人妻精品视频| 美女国产视频在线观看| 国产免费一区二区三区四区乱码| 性色av一级| 岛国毛片在线播放| 极品人妻少妇av视频| 亚洲av免费高清在线观看| 日韩在线高清观看一区二区三区| 免费观看av网站的网址| 男女边吃奶边做爰视频| 久热久热在线精品观看| 国产高清不卡午夜福利| 午夜久久久在线观看| 免费黄色在线免费观看| 麻豆成人av视频| 少妇被粗大猛烈的视频| 美女福利国产在线| 国产日韩欧美在线精品| 美女福利国产在线| 一区二区三区免费毛片| 亚洲国产精品国产精品| 在线天堂最新版资源| 多毛熟女@视频| 韩国高清视频一区二区三区| 欧美日本中文国产一区发布| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久com| 热99久久久久精品小说推荐| 日韩精品有码人妻一区| 免费日韩欧美在线观看| 中国美白少妇内射xxxbb| 男人操女人黄网站| 午夜福利影视在线免费观看| 亚洲人成网站在线观看播放| 超色免费av| 亚洲,一卡二卡三卡| 亚洲国产毛片av蜜桃av| 夜夜看夜夜爽夜夜摸| 免费久久久久久久精品成人欧美视频 | 熟女人妻精品中文字幕| 亚洲精品乱码久久久v下载方式| 国产免费一级a男人的天堂| 欧美性感艳星| 精品亚洲成a人片在线观看| 亚洲av国产av综合av卡| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡 | 成人毛片a级毛片在线播放| 观看美女的网站| 中文字幕人妻熟人妻熟丝袜美| 日韩一本色道免费dvd| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 欧美日韩一区二区视频在线观看视频在线| 一级a做视频免费观看| 国产无遮挡羞羞视频在线观看| 国产精品一区二区三区四区免费观看| 亚洲综合色惰| 性色av一级| 欧美最新免费一区二区三区| 午夜久久久在线观看| 在线亚洲精品国产二区图片欧美 | 一级毛片我不卡| 亚洲精华国产精华液的使用体验| 国产欧美日韩综合在线一区二区| 高清欧美精品videossex| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 国产片特级美女逼逼视频| 在线精品无人区一区二区三| 日韩三级伦理在线观看| 欧美3d第一页| 国产免费福利视频在线观看| 99九九线精品视频在线观看视频| 内地一区二区视频在线| 久久热精品热| 精品少妇内射三级| 亚洲精品一二三| 亚洲综合色惰| 韩国av在线不卡| 日本免费在线观看一区| 国产精品一区二区在线不卡| 999精品在线视频| 日韩人妻高清精品专区| 少妇猛男粗大的猛烈进出视频| 青春草视频在线免费观看| 日本午夜av视频| 国产精品久久久久成人av| 久久亚洲国产成人精品v| 日韩三级伦理在线观看| 人妻系列 视频| 免费高清在线观看日韩| 成人漫画全彩无遮挡| 一级爰片在线观看| 免费观看a级毛片全部| 亚洲欧美中文字幕日韩二区| 国产免费现黄频在线看| 一边亲一边摸免费视频| 亚洲图色成人| 中文欧美无线码| 大香蕉97超碰在线| 大香蕉久久网| 激情五月婷婷亚洲| 老女人水多毛片| 26uuu在线亚洲综合色| 久久久久国产精品人妻一区二区| 久久女婷五月综合色啪小说| 国产日韩一区二区三区精品不卡 | 黄色配什么色好看| xxx大片免费视频| 免费黄网站久久成人精品| 亚洲av欧美aⅴ国产| av电影中文网址| 亚洲国产日韩一区二区| 91久久精品国产一区二区三区| 亚洲国产毛片av蜜桃av| 精品人妻偷拍中文字幕| 免费av中文字幕在线| 国产在线一区二区三区精| 久久99精品国语久久久| 日本欧美国产在线视频| 99久久精品一区二区三区| 亚洲av中文av极速乱| 国模一区二区三区四区视频| 国产精品麻豆人妻色哟哟久久| 久久免费观看电影| 亚洲欧美清纯卡通| 欧美最新免费一区二区三区| 精品久久国产蜜桃| 久久久久久人妻| 2022亚洲国产成人精品| 亚洲av免费高清在线观看| 精品国产一区二区三区久久久樱花| av在线播放精品| 色5月婷婷丁香| √禁漫天堂资源中文www| 在线观看一区二区三区激情| 久久久久精品久久久久真实原创| 日本黄色片子视频| a级毛片在线看网站| 老司机亚洲免费影院| 欧美最新免费一区二区三区| av国产精品久久久久影院| 久久亚洲国产成人精品v| 少妇的逼好多水| 在线观看一区二区三区激情| 国产伦精品一区二区三区视频9| 丝袜在线中文字幕| 亚洲av国产av综合av卡| 只有这里有精品99| 啦啦啦中文免费视频观看日本| 国产成人午夜福利电影在线观看| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 亚洲国产av新网站| 日韩中字成人| 久久久久精品性色| 一级毛片黄色毛片免费观看视频| 人人澡人人妻人| 熟女av电影| 黑人猛操日本美女一级片| 欧美一级a爱片免费观看看| 我的女老师完整版在线观看| 久久久久久人妻| 欧美人与性动交α欧美精品济南到 | 国产精品99久久99久久久不卡 | 亚洲色图 男人天堂 中文字幕 | 国产av国产精品国产| 中国美白少妇内射xxxbb| 黑人猛操日本美女一级片| 欧美一级a爱片免费观看看| 国产成人91sexporn| 免费黄网站久久成人精品| 如日韩欧美国产精品一区二区三区 | 精品人妻熟女av久视频| 寂寞人妻少妇视频99o| 高清欧美精品videossex| 色网站视频免费| 日日摸夜夜添夜夜爱| 伊人久久精品亚洲午夜| 91久久精品国产一区二区三区| 日韩精品免费视频一区二区三区 | 日韩一本色道免费dvd| 精品一区二区三区视频在线| av.在线天堂| 国产精品久久久久成人av| 国产av码专区亚洲av| 国产免费一区二区三区四区乱码| 黄色一级大片看看| 天堂8中文在线网| 亚洲精品久久成人aⅴ小说 | a级毛片黄视频| 永久网站在线| 国产在视频线精品| 中文字幕人妻丝袜制服| 一本大道久久a久久精品| 狠狠婷婷综合久久久久久88av| 亚洲综合精品二区| 插逼视频在线观看| 婷婷色综合www| 黑丝袜美女国产一区| av卡一久久| 人人妻人人澡人人爽人人夜夜| 中文字幕av电影在线播放| 国产欧美亚洲国产| 高清不卡的av网站| 伦理电影大哥的女人| 亚洲成人手机| 欧美最新免费一区二区三区| 99久久综合免费| 国产高清不卡午夜福利| 性色avwww在线观看| 国产精品久久久久久久电影| 精品卡一卡二卡四卡免费| 大香蕉久久网| 啦啦啦视频在线资源免费观看| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂| 亚洲丝袜综合中文字幕| 多毛熟女@视频| 欧美变态另类bdsm刘玥| av又黄又爽大尺度在线免费看| 黄色欧美视频在线观看| 有码 亚洲区| 秋霞伦理黄片| 男女高潮啪啪啪动态图| 18+在线观看网站| 国产精品一区二区在线不卡| 91精品三级在线观看| 久久韩国三级中文字幕| .国产精品久久| 18禁在线无遮挡免费观看视频| 成人亚洲欧美一区二区av| 亚洲欧美清纯卡通| 高清黄色对白视频在线免费看| 熟女人妻精品中文字幕| 亚洲一区二区三区欧美精品| 国产精品99久久99久久久不卡 | 色婷婷av一区二区三区视频| 91成人精品电影| 国产精品人妻久久久影院| 免费人成在线观看视频色| 免费av中文字幕在线| 97超碰精品成人国产| 在线观看美女被高潮喷水网站| 99re6热这里在线精品视频| 国精品久久久久久国模美| 国产成人精品在线电影| 午夜福利影视在线免费观看| 国产成人精品久久久久久| 黄色视频在线播放观看不卡| 婷婷色麻豆天堂久久| 欧美亚洲日本最大视频资源| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 日韩大片免费观看网站| 大又大粗又爽又黄少妇毛片口| 免费av中文字幕在线| 麻豆成人av视频| 高清不卡的av网站| 国产av精品麻豆| 久久久久网色| 欧美精品人与动牲交sv欧美| 亚洲精品第二区| 哪个播放器可以免费观看大片| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看| 久久婷婷青草| 青青草视频在线视频观看| 国产白丝娇喘喷水9色精品| 国语对白做爰xxxⅹ性视频网站| 人妻夜夜爽99麻豆av| 人妻系列 视频| 亚洲av成人精品一二三区| 免费少妇av软件| 丰满乱子伦码专区| 蜜桃在线观看..| 亚洲欧美日韩卡通动漫| 国产av精品麻豆|