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

    新疆一枝蒿多糖的提取、純化及其結構分析

    2018-08-30 11:41:46王旭輝涂振東晁群芳
    新疆農業(yè)科學 2018年6期
    關鍵詞:群芳新疆大學科學院

    王旭輝,徐 鑫,王 卉,葉 凱,涂振東,晁群芳

    (1.新疆大學生命科學與技術學院,烏魯木齊 830046;2.新疆農業(yè)科學院生物質能源研究所,烏魯木齊 830091;3. 新疆農業(yè)科學院糧食作物研究所,烏魯木齊 830091)

    0 引 言

    【Research significance】ArtemisiarupestrisL., belonging to theCompositaefamily, is a well-known traditional Uygur medicinal herb of Xinjiang province in China, and is mainly distributed in Kazakhstan, Mongolia and Russia. This plant is commonly used for the treatment of many diseases, such as detoxication, anti-hypersusceptibility, antitumor, antibacterial, antivirus and protecting liver, and it contains many kinds of bioactive chemicals, primarily includingavonoids, terpenoids, alkaloid, rupestonic acid[1]. The recent pharmacology research indicated that it has the functions of antitumor, antimicrobial, antioxidant, and improving cellular immunity.【Previous research progress】Polysaccharides from plant, epiphyte and animals extracts are an interesting source of additives for several industries, in particular food and drug industry[2]. They play an important role in the growth and development of living organisms and have been widely studied in recent years due to their unique biological, chemical and physical properties[3]. Published literatures have indicated that plant polysaccharides in general have strong antioxidant activities and can be explored as novel potential antioxidants[4-6]. Hot water technology is the main and most conventional extraction method for polysaccharides mentioned in previous studies[7, 8]. However, it is usually associated with longer extraction time and higher temperature but lower extraction efciency. Therefore, it is desirable tond an effective and economical method for bioactive compounds extraction. Ultrasonic has been used to increase extraction yield of bioactive substances from natural products, which is mainly attributed to disruption of cell walls, particle-size reduction, and enhanced mass transfer to the cell contents as a result of cavitation bubble collapse[9-12]. It offers high reproducibility at shorter times, simplied manipulation, and lowered energy input, as well as solvent consumption[13]. Optimizing refers to improving the performance of a system, a process, or a product in order to obtain the maximum benet from it. When many factors and interactions affect desired response, response surface methodology (RSM) is an effective tool for optimizing the process, which was originally developed by Box and Wilson in the 50s[14]. The main advantage of RSM is the reduced number of experimental trials needed to evaluate multiple parameters and their interactions. Therefore, it is less laborious and time-consuming than other approaches required to optimize a process[15], and it has successfully been applied in order to optimize the conditions in food and pharmaceutical research[16-20]. 【Entry point】However, so far there is not any information published about the optimization, purication and structural characterization of polysaccharides fromA.rupestrisL.. For that reason, extraction of polysaccharides is an important step for its application or further research and development. 【The key problems intended to be solved】Therefore, the aim of this work was mainly to report on the extraction and purication of polysaccharides fromA.rupestrisL.. Firstly, single-factor experimental designs (ultrasonic power, ratio of water to raw material, ultrasonic time and extracting times) were carried out before RSM experiments. Secondly, three factors (ultrasonic power, ratio of water to raw material and ultrasonic time) were chosen based on single-factor designs for further optimization by employing a three-level, three-variable BBD from RSM. Furthermore, after purication by gelltration chromatography, ultraviolet spectrum (UV) and Fourier transform-infrared spectroscopy (FT-IR) were employed to identify the chemical structure.

    1 Materials and methods

    1.1 Materials

    DriedA.rupestrisL. was obtained from Urumqi, Xinjiang province, China in September, 2010, the whole plants were cut into smaller pieces and further ground into ane powder in a high speed disintegrator (160C, Fei Da Chinese Traditional Medicine Machine Co., Ltd, Zhejiang, China), and passed through a 40 mesh sieve. The materials were stored at room temperature in a desiccator until use(less than two weeks). All other reagents used in this study were of analytical grade.

    1.2 Ultrasonic extraction

    The extraction of ARP was done according to the reported method with some modications[21-23]. Briey, the extraction process was performed using an ultrasonic cell disrupter (JY92-2D, Ningbo Xinzhi Biological Technology Co., Ltd., Zhejiang, China), with different ultrasonic power, ratios of water to raw material and the ultrasonic time.

    DriedA.rupestrisL. powder was defatted in a Soxhlet apparatus with Petroleum ether (boiling point: 60-90℃) and then were extracted for three times with 80% ethanol at 60℃ and 2 h each time to defat and remove some colored materials, oligosaccharides, and some small molecule materials[24, 25]. The pretreated samples were separated from the organic solvent by centrifugation (3,500 g for 15 min). Each dried pretreated sample was extracted by water in a designed ultrasonic power, ratios of water to raw material, the ultrasonic time and number. The water extraction solutions were separated from insoluble residue through the nylon cloth (Pore diameter: 38 nm), was then concentrated to 25 mL with a rotary evaporator (RE-5299, Gongyi Yuhua Instrument Co., Ltd., He’nan, China) at 55℃ under vacuum[6]. The proteins in the extract were removed using the Sevag reagent[26]. After removal of the Sevag reagent, the supernatant was precipitated by the addition of 95% ethanol to anal concentration of 80% (v/v) to obtain polysaccharides, and kept overnight at 4℃. Then the solution was centrifuged at 4,500 rpm/min for 20 min, washed three times with dehydrated ethanol, and then the precipitate was solubilized in deionized water and lyophilized by vacuum freeze drying machine (D-37520, Ostevode, telefon, Germany). The content of the polysaccharides was measured by anthrone-sulfuric method. The percentage polysaccharides extraction yield (%) is calculated as the polysaccharides content of extraction divided by dried sample weight. All experiments were performed at least in duplicate.

    1.3 Experimental design and statistical analysis

    To explore the effects of independent variables on the response within the range of investigation, proper ranges of ultrasonic power, ratio of water to raw material, ultrasonic time, and extraction number were preliminarily determined. A three-level-three-factor, Box-Behnken factorial design (BBD) employed in this optimization study. Based on the investigations on single-factor experiments, ultrasonic powers (X1), ratio of water to raw material (X2), ultrasonic time (X3) were chosen for independent variables to be optimized for the extraction yield of ARP. Yield of polysaccharides (Y) was taken as the response of the design experiments. For statistical calculation, the variables were coded according to

    (1)

    Wherexiwas a coded value of the variable;Xiwas the actual value of variable;X0was the actual value of theXion the center point; and △Xwas the step change value. The range of design factors and their levels is presented in Table 1, which was based on the results of preliminary experiments. As seen from Table 2, the whole design consisted of 17 experimental points carried out in random order. Five replicates at the center of the design were used for the estimation of a pure error sum of squares. The response value in each trial was average of duplicates. Table 1, Table 2

    Table 1 Response surface methodology design factors and levels

    Independent variablesFactor level-101X1: Ultrasonic power (w)600700800X2: Ratio of water to raw material303540X3: Ultrasonic time (min)202530

    Table 2 Box-Behnken design and the response value for yield of polysaccharides

    RunX1Ultrasonic power(w)X2Ratio of water to raw materialX3Ultrasonic time (min)YPolysaccharides yield (%)1-1-102.6120-1-13.3230003.9240003.9650003.9560-113.3470113.5781103.2890003.93100003.931101-13.5112-1102.721310-13.3114-10-12.73151-102.8516-1012.89171013.25

    Based on the experimental data, regression analysis was performed and wastted into an empirical second-order polynomial model:

    (2)

    WhereYis the response function,A0was constant.Ai,AiiandAijwere the coefcients of the linear, quadratic and interactive terms, respectively. AndXiandXjrepresented the coded independent variables. The coefcients of the second polynomial model and the responses obtained from each set of experimental design were subjected to multiple nonlinear regressions using software Design-Expert 8.0.5b (Trial Version, State-Ease Inc., Minneapolis, MN, USA). Thetness of the polynomial model equation was expressed by the coefcient of determinationR2, and its statistical signicance was checked by F-test at a probability (P) of 0.05. The signicances of the regression coefcients were also tested by F-test. Additional conrmation experiments were subsequently conducted to verify the validity of the experimental design.

    1.4 Purication of the polysaccharide

    The crude polysaccharide was re-dissolved in 5 ml distilled water,ltered through 4.5 × 10-4mmlters and applied to a DEAE-52 cellulose column (3.0 × 60 cm) equilibrated with distilled water. The polysaccharide was fractionated and eluted with distilled water and different concentrations of stepwise NaCl solution (0,0.1,0.2 ,0.3 and 0.4 M NaCl) at aow rate of 2.0 mL/min. The elutes were concentrated to obtain the main fractions, which were then fractionated by size-exclusion chromatography on a Sephadex G-100 column (2.5 × 40 cm) eluted with 0.2 M NaCl at aow rate of 0.5 mL/min. The relevant fractions were collected, concentrated, dialyzed and lyophilized.

    1.5 Ultraviolet analysis

    ARP-2 was dissolved and diluted to 2 mg/mL respectively, and the solutions of the polysaccharide fraction were scanned from 200 to 400 nm with a UV-VIS-NIR spectrophotometer (UV3600, Shimadzu Corporation, Japan).

    1.6 FT-IR spectrometric analysis

    The FT-IR spectra were recorded using the KBr-disk method with a Nicolet Fourier transform infrared (FTIR) spectrometer (EQINOX 55, BRUKER Corporation, Germany) in the range 400-4,000 cm-1.

    1.7 Statistical analysis

    Comparison of means was performed by one-way analysis of variance (ANOVA) followed by LSD (Least significant difference method). Statistical analyses (P< 0.05) were performed using SPSS 16.0 software. Design Expert 8.0.5b (Trial Version, State-Ease Inc., Minneapolis, MN, USA) was employed for experimental design, analysis of variance (ANOVA), and model building. All analyses were performed in triplicate.

    2 Results and discussion

    2.1 Effect of different ultrasonic power on extraction yield of ARP

    The effect of different ultrasonic power on extraction yield of ARP is shown in Fig. 1(a). Extraction was carried out at different ultrasonic power (200-1,000 W) conditions while other extraction parameters were as follows: ratios of water to raw material of 30, ultrasonic time of 25 min, and number of extraction of 2. The extraction yield of ARP signicantly increased from 2.08% to 3.71% as ultrasonic power increased from 200 to 1,000 W. It can be found that signicant differences were existing among 200 W, 400 W and 600 W, and among 200 W, 400 W, 600 W and 1,000 W (P< 0.05), but there was no signicant difference between 600 W and 800 W, and among 800 W and 1,000 W (P> 0.05). Although the extraction yield of ARP was also high at 1,000 W, increasing ultrasonic power will bring about the increase in cost for the extraction process from the industrialization point of view. Therefore, ultrasonic power range of 600-800 W was considered to be optimal in the present experiment. Fig.1

    Fig.1 Effort of ultrasonic power (a), ratio of water to raw material (b), ultrasonic time(c) and number of extraction (d) on the yield of polysaccharides.Values have no significant statistical difference at 0.05 with the same superscript letters

    2.2 Effect of different ratio of water to raw material on extraction yield of ARP

    The effect of different ratio of water to raw material on extraction yield of ARP is shown in Fig. 1(b). Extraction was carried out at different ratio of water to material (10-50) conditions while other extraction parameters were as follows: ultrasonic power of 700 W, ultrasonic time of 25 min, and number of extraction of 2. The extraction yield of ARP signicantly increased from 1.92% to 3.86%. Results indicated that signicant differences were existing among 10, 20, 30 and 40, and among 10, 20, 30 and 50 (P< 0.05), but there was no signicant difference between 30 and 40 (P>0.05). Therefore, ratio of water to raw material range of 30-40 was considered to be optimal in the present experiment.

    2.3 Effect of different ultrasonic time on extraction yield of ARP

    Ultrasonic time is a factor that would inuence the extraction efciency and selectivity of theuid. This might be due to the time requirement of the exposure ofA.rupestrisL. to the release medium where the liquid penetrated into the driedA.rupestrisL. powdered, dissolved it and subsequently diffused out from the material. A longer ultrasonic time presents a positive effect on the extraction yield of ARP. The effect of different ultrasonic time on the extraction yield of ARP is shown in Fig. 1(c). Extraction was carried out at different ultrasonic time conditions while other extraction variables were set as follows: ultrasonic power of 700 W, numbers of extraction of 2, and ratio of water to raw material of 30. When ultrasonic time varied from 10 to 30 min, the variance of extraction yield was relatively rapid, and the extraction yield of ARP reached a maximum at 20-30 min, and then no longer changed as the extraction proceeded. Results showed that signicant differences were existing among 10 min, 15 min, 20 min and 25 min, and among 10 min, 15 min, 20 min and 30 min (P< 0.05), but there was no signicant difference between 25 min and 30 min (P> 0.05). This indicated that ultrasonic time of 25 min was sufcient to obtain the polysaccharide production. Thus, extraction of 20-30 min was favorable for producing the polysaccharides.

    2.4 Effect of different number of extraction on extraction yield of ARP

    The effect of number of extraction on extraction yield of ARP is shown in Fig. 1(d). While other extraction parameters weretted as follows: ultrasonic power 700 W, ultrasonic time 25 min, and ratio of water to raw material 30. From Fig. 1 (d), it can be found that there was an increasing trend in the yield of ARP accompanying the increase of extracting times, but there was not signicant difference (P>0.05) between 2 times and 3 times. Taking the yield and processing cost into consideration, 2 times was sufcient for the extraction of ARP. Thus, 2 times was selected as the extracting times in the next experiments.

    According to the single-parameter study, ultrasonic power of 600-800 W, ratio of water to raw material of 30-40 and ultrasonic time of 20-30 min for RSM experiments were adopted.

    2.5 Optimization of the procedure

    2.5.1 Statistical analysis and the modeltting

    The value of responses (the extraction yield of ARP) at different experimental combination for coded variables is given in Table 2. The extraction yield of ARP ranged from 2.61% to 3.96%. By employing multiple regression analysis on the experimental data, the predicted response for the yield of ARP can be obtained by the following second-order polynomial equation:

    Y=3.94+0.22×X1+0.12×X2+0.023×X3+0.080×X1×X2-0.055×X1×X3+0.010×X2×X3-0.73×X1×X1-0.34×X2×X2-0.16×X3×X3

    (3)

    Table 3 Analysis of variance for thetted quadratic polynomial model of extraction of polysaccharides

    ItemStd.devMeanC.V.%PressR2R2AdjR2PredAdeq precisionValue0.0283.360.830.0700.998 50.996 60.980 562.574

    The Lack of Fit is an indication of the failure for a model representing the experimental data at which points were not included in the regression or variations in the models cannot be accounted for random error[27]. If there is a signicant Lack of Fit which could be indicated by a low probability value, the response predictor is discarded. The Lack of Fit did not result in a signicant P-value for selected variables, meaning that these models were sufciently accurate for predicting the relevant responses. The Lack of Fit F-value of 5.31 implies that the Lack of Fit is signicant. There is only a 7.03% chance that a “Lack of Fit F-value” could occur due to noise, which indicated that the model equation was adequate for predicting the yield of ARP under any combination of values of the variables. Table 4

    Table 4 Regression coefcient estimation and their signicance test for the quadratic polynomial model

    VariablesSun of squaresDFMean squareF valueP value prob.>FModel3.6190.04521.25<0.000 1X10.3810.38492.41<0.000 1X20.1210.12149.89<0.000 1X34.05×10314.05×1035.270.055 4X1×X12.2512.252 931.44<0.000 1X2×X20.4910.49638.9<0.000 1X3×X30.1110.11142.89<0.000 1X1×X20.02610.02633.310.000 7X1×X30.01210.01215.740.005 4X2×X34×10414×1040.520.494 0Residual5.38×10377.686×104--Lack of fit4.3×10331.433×1035.310.070 3Pure error1.08×10342.7×104--Correlation total3.6116---

    The model was found to be adequate for prediction within the range of experimental variables. The regression coefcient values of Eq. (3) were listed in Table 4. The P-values were used as a tool to check the signicance of each coefcient, which in turn may indicate the pattern of the interactions between the variables. The smaller the value of P was, the more signicant the corresponding coefcient was[28]. It can be seen from this table that the linear coefcients (X1,X2), a quadratic term coefcient (X12,X22,X32) and cross product coefcients (X1×X2) were signicant, with very small P-values (P< 0.05). The other term coefcients were not signicant (P> 0.05). Therefore,X1,X2,X12,X22,X32andX1×X2were important factors in the extraction process of the polysaccharides. The full modellled Eq. (3) was made the 3-D response surface plots and contour plots to predict the relationships between the independent and dependent variables.

    2.5.2 Optimization of extraction conditions of ARP

    The response surface curves were plotted using Design-Expert to study the effects of parameters and their interactions on polysaccharides yield. The 3-D response surface plots were drawn to illustrate the main and interactive effects of the independent variables on the dependent one. They provide a method to visualize the relationship between responses and experimental levels of each variable and the type of interactions between two test variables. These types of plot (circular or elliptical) indicated the mutual interactions between the variables are signicant or not. Circular contour plot indicates that the interactions between the corresponding variables are negligible, while elliptical contour plot indicates that the interactions between the corresponding variables are signicant[29].In the present study, the 3-D response surface plots and the contour plots, as presented in Figs. 2. It is clear that the yield of ARP is sensitive to minor alterations of the test variables (ultrasonic power, ratio of water to raw material and ultrasonic time). Through these 3-D response surface plots and their respective contour plots, it is very easy and convenient to understand the interactions between two variables and to locate their optimum ranges.

    The 3-D response surface plot and the contour plot in Figs. 2 (a) and (d), which give the extraction yield of ARP as a function of ultrasonic power and ratio of water to raw material atxed ultrasonic time (0 level), indicated that the extraction yield of ARP increased with increase of ultrasonic power from 600 to 715.78 W, then decreased slightly from 715.78 to 800 W, and increased rapidly with increase of ratio of water to raw material from 30 to 35.98, then dropped from 35.98 to 40. Figs. 2 (b) and (e) shows the 3-D response surface plot and the contour plot at varying ultrasonic power and ultrasonic time atxed ratio of water to raw material (0 level). From twogures, it indicated that the extraction yield of ARP increased with increase of ultrasonic power from 600 to 715.78 W, then decreased from 715.78 to 800 W, and the extraction yield of ARP was found to increase rapidly with increase of ultrasonic time from 20 to 25.25 min, but beyond 25.25 min, the extraction yield of ARP reached the plateau region where the yield was maximized and did not further increase the yield. The extraction yield of ARP affected by different ultrasonic time and ratio of water to raw material was seen in Figs. 2 (c) and (f), when ultrasonic power wasxed at 0 level. It can be seen that the extraction yield of ARP increased with the increase of ratio of water to raw material from 30 to 35.98, then further increased with the increase of ultrasonic time, and reached the maximum value when ultrasonic time at 25.25 min, and beyond this level, the extraction yield of ARP did not further increase.

    From Figs. 2, it can be concluded that the optimum extraction conditions of polysaccharides fromA.rupestrisL. are ultrasonic power 715.78 W, ratio of water to material 35.98 and ultrasonic time 25.25 min. The theoretical extraction yield of ARP that was predicted under the above conditions was 3.97%.

    Fig.2

    The suitability of the model equations for predicting optimum response values was tested under the conditions: ultrasonic power 715.78 W, ratio of water to material 35.98 and ultrasonic time 25.25 min. A predicted value of 3.97% was obtained for yield of polysaccharides under the optimal conditions. In order to facilitate the practical extraction process of ARP, the optimal conditions were modied as follows: ultrasonic power 716 W, ratio of water to material 36 and ultrasonic time 25 min. A predicted value of 3.94% was obtained under the modied conditions. The modied conditions were used to validate the suitability of thetted model equation for accurately predicting the responses values. The results showed that the actual values of polysaccharides yield were 3.92%± 0.065% under the modied conditions (Table 5), which were in agreement with the predict values signicantly (P>0.05). Thus, the model can be used to optimize the process of polysaccharides extraction fromA.rupestrisL.. Table 5

    Fig.2 Response surface plots (a-c) and contour plots (d-f) showing the effects of the variables on the yield of polysaccharide

    Table 5 Predicted and experimental values of the responses at optimum conditions

    Ultrasonic power(W)Ratio of water to raw material(mL/g)Extraction time(min)Extraction yield of ARP(%)(Experimental)Extraction yield of ARP(%)(Predicted)Optimum conditions715.7835.9825.253.963.97Modi ed conditions71636253.923.94

    The representative values were obtained from three independent experiments, and the results are presented as mean value ± SD

    2.6 Isolation and purication of polysaccharides

    GroundA.rupestrisL. was reuxed with ethanol to deactivate the endogenous enzymes and remove some soluble materials, including free sugars, amino acids and some phenols. Then the dried ethanol-extracted residue was extracted with distilled water at 80℃. After dialysis and precipitation with ethanol, a crude water-soluble polysaccharide (ARP) was obtained as a brownish powder. The ARP was separated and puried by DEAE-52 cellulose column. Four fractions of ARP-1, ARP-2, ARP-3 and ARP-4 were eluted. ARP-2 was further puried through Sephadex G-100 column (Fig. 3(a).and (b).). Fig.3

    (a) Elution prole of ARP on a DEAE-cellulose ion-exchange column (at 2 mL/min mobile phaseow rate); (b) Elution prole of ARP-2 (from a) on a Sephadex G-100 column (eluted with 0.2 M mobile phase at 0.5 mL/min). Absorbance at 620 nm (A620) represents the total carbohydrate content (relative) determined by the anthrone-sulfuric method. Spectral Analysis of ARP-2:(c)The UV spectra of ARP-2;(d)FT-IR analysis of ARP-2

    Fig.3 Fractionation of ARP by column chromatography

    2.7 Ultraviolet analysis

    ARP-2 had no absorption at 280 and 260 nm in the UV spectrum (Fig.3(c).), indicating the absence of protein and nucleic acid.

    2.8 FTIR analysis

    The FTIR spectrum of the polysaccharide in the wavelength range of 4,000-400 cm-1was shown in (Fig.3(d).). In the spectrum of polysaccharide, the bands in the region of 3,410 cm-1were due to the hydroxyl stretching vibration of the polysaccharides. The bands in the region of 2,938-1were due to C-H stretching vibration[5,31]. The characteristic absorptions at 917 cm-1in the IR spectra indicated that was β-congurations. Stronger bands occurring between 1,742 cm-1and 1,611 cm-1are derived from the ester carbonyl (-COOR) groups and carboxylate ion stretching band (-COO-)[30].

    3 Conclusions

    Compared to conventional extraction techniques, ultrasonic technology improves the extraction efciency within a shortened period of time. Ultrasonic technology was performed for the polysaccharides extraction fromA.rupestrisL. in order to increase the yield extraction. Based on the single-factor experiments, Box-Behnken design (BBD) was used to optimizes process parameters and improves polysaccharide extraction. ANOVA showed that the effects of the variables are signicant and quadratic models are obtained for predicting responses. The optimum set of the independent variables was obtained graphically in order to obtain the desired levels of polysaccharides extraction. The optimum extraction conditions were modified as follows: ultrasonic power 716 W, ratio of water to raw material 36 and ultrasonic time 25 min. Under these optimal conditions, the mean polysaccharides extraction yield was 3.92%± 0.065%, which corresponded well with the predicted yield of 3.94%. Then, the crude polysaccharides were puried byltration, DEAE-52 cellulose column and Sephadex G-100 size-exclusion chromatography in that order. A main fraction, ARP-2 was obtained through the extraction and purication steps. Finally, The absorption spectra of the polysaccharides showed that the polysaccharides-related absorption peaks at 917, 1,100, 1,147, 1,332, 1,423, 1,611, 1,742, 2,939 and 3,410 cm-1. This study could become the basis for later research according to the optimization predictions of the extraction. Further study should be carried out to study bioactivity in vitro and in vivo, with the purpose of applying the polysaccharides as a potential natural antioxidant functional ingredient in the food industry.

    猜你喜歡
    群芳新疆大學科學院
    新疆大學紡織與服裝學院攝影作品選登
    詠 梅
    咸寧市農業(yè)科學院情況簡介
    咸寧市農業(yè)科學院農機所簡介
    湖北農機化(2020年4期)2020-07-24 09:06:58
    A Corpus-Based Study on Linguistic Variables of CET Writings
    《河北省科學院學報》稿約
    鴨蘭
    詠梅詩
    評《農產品物流》
    一位科學院院士的文學人生
    海峽姐妹(2018年8期)2018-09-08 07:58:52
    色5月婷婷丁香| 麻豆精品久久久久久蜜桃| 精品久久久久久久久亚洲| 少妇被粗大猛烈的视频| 久久婷婷人人爽人人干人人爱| 看十八女毛片水多多多| av在线观看视频网站免费| 少妇丰满av| 成年女人永久免费观看视频| 99国产精品一区二区蜜桃av| 久久久国产成人精品二区| 亚洲av成人精品一区久久| 久久久久国产精品人妻aⅴ院| 国产亚洲精品久久久久久毛片| 日韩一本色道免费dvd| av免费在线看不卡| 午夜福利在线观看吧| 可以在线观看的亚洲视频| 色尼玛亚洲综合影院| 国产亚洲精品综合一区在线观看| 久久热精品热| 十八禁网站免费在线| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 老熟妇乱子伦视频在线观看| 热99在线观看视频| 六月丁香七月| 日本在线视频免费播放| 成人亚洲欧美一区二区av| 搡老岳熟女国产| 久久久精品欧美日韩精品| 麻豆国产av国片精品| 久久久a久久爽久久v久久| 日韩精品中文字幕看吧| 国产高清激情床上av| 欧美日韩乱码在线| 老司机福利观看| 国产午夜精品论理片| 国产精品久久久久久久久免| 久久久久国产网址| 久久精品国产自在天天线| 成人漫画全彩无遮挡| 国产黄色视频一区二区在线观看 | 舔av片在线| 99久国产av精品国产电影| 亚洲精品乱码久久久v下载方式| 伦精品一区二区三区| 国产三级中文精品| 免费观看人在逋| 麻豆精品久久久久久蜜桃| av在线老鸭窝| 悠悠久久av| 国产欧美日韩精品亚洲av| 美女cb高潮喷水在线观看| 日本三级黄在线观看| 看免费成人av毛片| 精品国产三级普通话版| 在线观看午夜福利视频| 欧美国产日韩亚洲一区| 国产精品久久久久久久电影| 久久精品影院6| 久久午夜亚洲精品久久| 国产激情偷乱视频一区二区| 男人舔奶头视频| 亚洲内射少妇av| av女优亚洲男人天堂| 亚洲欧美日韩高清在线视频| 少妇熟女欧美另类| 18+在线观看网站| 成人毛片a级毛片在线播放| 国产一区亚洲一区在线观看| 午夜爱爱视频在线播放| 欧美日韩国产亚洲二区| 看非洲黑人一级黄片| 免费大片18禁| 亚洲久久久久久中文字幕| 欧美不卡视频在线免费观看| 亚洲综合色惰| 最近在线观看免费完整版| 国内精品久久久久精免费| 午夜福利在线观看免费完整高清在 | 人人妻人人看人人澡| 22中文网久久字幕| 俄罗斯特黄特色一大片| 校园春色视频在线观看| 亚洲专区国产一区二区| 国产一区亚洲一区在线观看| 午夜激情欧美在线| 国产精品99久久久久久久久| 午夜福利18| 少妇裸体淫交视频免费看高清| a级毛色黄片| 久久亚洲国产成人精品v| 免费看a级黄色片| 国产精品精品国产色婷婷| a级毛色黄片| 亚洲,欧美,日韩| 国产午夜精品久久久久久一区二区三区 | 日韩亚洲欧美综合| 欧美高清成人免费视频www| 日本五十路高清| 欧美性感艳星| 日本在线视频免费播放| 亚洲国产高清在线一区二区三| 亚洲精品色激情综合| 亚洲av免费高清在线观看| 国产色婷婷99| 99国产极品粉嫩在线观看| 亚洲无线在线观看| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 日韩中字成人| 午夜福利视频1000在线观看| 真人做人爱边吃奶动态| 日本与韩国留学比较| 我要搜黄色片| 中文字幕人妻熟人妻熟丝袜美| 又黄又爽又免费观看的视频| 亚洲精品色激情综合| 亚洲国产精品成人久久小说 | 国产高清视频在线播放一区| 少妇裸体淫交视频免费看高清| 欧美成人精品欧美一级黄| 色综合色国产| 国产成年人精品一区二区| 老司机影院成人| 国产精品电影一区二区三区| 久久国产乱子免费精品| 日本一二三区视频观看| 国产亚洲精品av在线| av中文乱码字幕在线| 搡老妇女老女人老熟妇| 观看免费一级毛片| 午夜福利在线在线| 久久韩国三级中文字幕| 免费人成在线观看视频色| 嫩草影视91久久| 美女xxoo啪啪120秒动态图| 国产精品伦人一区二区| 中文亚洲av片在线观看爽| 免费高清视频大片| 91精品国产九色| 看十八女毛片水多多多| 女生性感内裤真人,穿戴方法视频| 少妇人妻精品综合一区二区 | 欧美日韩一区二区视频在线观看视频在线 | 色噜噜av男人的天堂激情| 性插视频无遮挡在线免费观看| 国产又黄又爽又无遮挡在线| 国产av在哪里看| 久久6这里有精品| 久久精品国产鲁丝片午夜精品| 插逼视频在线观看| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 永久网站在线| 日韩成人av中文字幕在线观看 | 国产av不卡久久| 国产中年淑女户外野战色| 日韩,欧美,国产一区二区三区 | eeuss影院久久| 国产成人影院久久av| 日韩一区二区视频免费看| 久久人人精品亚洲av| 亚洲婷婷狠狠爱综合网| 97超碰精品成人国产| 久久精品久久久久久噜噜老黄 | 国产精品久久视频播放| 国产精品野战在线观看| 国产伦精品一区二区三区视频9| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看| 日韩欧美三级三区| 美女高潮的动态| 人人妻人人澡人人爽人人夜夜 | av黄色大香蕉| av在线观看视频网站免费| 天天躁日日操中文字幕| 国产精品不卡视频一区二区| 日韩欧美三级三区| 色av中文字幕| 成人av一区二区三区在线看| 国产成人a∨麻豆精品| 国内精品久久久久精免费| 亚洲乱码一区二区免费版| 一边摸一边抽搐一进一小说| 特级一级黄色大片| 国产av不卡久久| 老司机午夜福利在线观看视频| 在线观看免费视频日本深夜| 日韩成人av中文字幕在线观看 | 听说在线观看完整版免费高清| 蜜桃久久精品国产亚洲av| av在线蜜桃| 在线观看美女被高潮喷水网站| av黄色大香蕉| 五月伊人婷婷丁香| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 丰满的人妻完整版| 午夜激情欧美在线| 国产成人精品久久久久久| 听说在线观看完整版免费高清| 尤物成人国产欧美一区二区三区| 亚洲精品日韩av片在线观看| 一进一出抽搐动态| 欧美日韩一区二区视频在线观看视频在线 | 日韩制服骚丝袜av| 日本黄色片子视频| 精品人妻偷拍中文字幕| 最近手机中文字幕大全| 成年版毛片免费区| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 看非洲黑人一级黄片| 成人av在线播放网站| 看免费成人av毛片| av在线天堂中文字幕| 精品久久久久久久久久久久久| 国产精品久久久久久亚洲av鲁大| 一本一本综合久久| 日本成人三级电影网站| 国产精品女同一区二区软件| 国产av麻豆久久久久久久| 久久精品国产自在天天线| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| 麻豆国产97在线/欧美| 亚洲av不卡在线观看| ponron亚洲| 少妇猛男粗大的猛烈进出视频 | 午夜影院日韩av| 99久久精品国产国产毛片| 国产探花极品一区二区| 嫩草影院精品99| 国产在视频线在精品| 婷婷精品国产亚洲av在线| 成人无遮挡网站| 国产精品不卡视频一区二区| 有码 亚洲区| 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| 国产精品久久久久久av不卡| 久久精品国产鲁丝片午夜精品| 欧美+日韩+精品| 国产精品久久久久久av不卡| 国产麻豆成人av免费视频| 级片在线观看| 内射极品少妇av片p| 久久久午夜欧美精品| 亚洲一级一片aⅴ在线观看| 久久韩国三级中文字幕| 午夜爱爱视频在线播放| 99精品在免费线老司机午夜| 人妻夜夜爽99麻豆av| 一夜夜www| 国产精品久久久久久久久免| 全区人妻精品视频| 十八禁网站免费在线| 免费看av在线观看网站| 欧美不卡视频在线免费观看| 亚洲美女黄片视频| 美女xxoo啪啪120秒动态图| 国产av一区在线观看免费| 成人av一区二区三区在线看| 狠狠狠狠99中文字幕| 天堂网av新在线| 在线观看av片永久免费下载| 精品免费久久久久久久清纯| 乱码一卡2卡4卡精品| 悠悠久久av| 最新中文字幕久久久久| а√天堂www在线а√下载| 男人狂女人下面高潮的视频| 99久久精品一区二区三区| 精品福利观看| 香蕉av资源在线| 高清毛片免费观看视频网站| 观看免费一级毛片| 熟妇人妻久久中文字幕3abv| 亚洲中文字幕日韩| 精品国产三级普通话版| 丝袜美腿在线中文| 中文字幕免费在线视频6| 亚洲中文日韩欧美视频| 精品人妻偷拍中文字幕| 久久精品国产清高在天天线| 少妇猛男粗大的猛烈进出视频 | 天堂√8在线中文| 99riav亚洲国产免费| 中出人妻视频一区二区| 老熟妇乱子伦视频在线观看| 91在线精品国自产拍蜜月| 日韩欧美在线乱码| 欧美精品国产亚洲| 亚洲av成人精品一区久久| 婷婷六月久久综合丁香| 精品99又大又爽又粗少妇毛片| 日本五十路高清| 成人av在线播放网站| 欧美日韩一区二区视频在线观看视频在线 | 长腿黑丝高跟| 日本一本二区三区精品| 久久久久久久久久黄片| 国产 一区 欧美 日韩| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 91久久精品电影网| 简卡轻食公司| 国产爱豆传媒在线观看| 简卡轻食公司| 日本免费一区二区三区高清不卡| 亚洲欧美中文字幕日韩二区| 欧美一区二区精品小视频在线| 春色校园在线视频观看| 久久精品91蜜桃| 亚洲国产精品成人综合色| 国产免费男女视频| 嫩草影院新地址| 日韩在线高清观看一区二区三区| 国产精品爽爽va在线观看网站| 国内精品久久久久精免费| 久久中文看片网| 久久精品国产亚洲av天美| 亚洲国产精品sss在线观看| 秋霞在线观看毛片| 午夜福利在线在线| 97超级碰碰碰精品色视频在线观看| 国产爱豆传媒在线观看| 最近2019中文字幕mv第一页| 日本成人三级电影网站| 成人特级av手机在线观看| 夜夜夜夜夜久久久久| 亚洲国产日韩欧美精品在线观看| 搡老岳熟女国产| 国产人妻一区二区三区在| 少妇高潮的动态图| 我的女老师完整版在线观看| 美女 人体艺术 gogo| 亚洲国产日韩欧美精品在线观看| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 亚洲国产精品sss在线观看| 久久久久国内视频| 精品国内亚洲2022精品成人| 一进一出抽搐gif免费好疼| 一夜夜www| 男女那种视频在线观看| 22中文网久久字幕| 两个人视频免费观看高清| 免费黄网站久久成人精品| 一区二区三区高清视频在线| 免费大片18禁| 一本精品99久久精品77| 午夜福利高清视频| 夜夜爽天天搞| 青春草视频在线免费观看| 男女下面进入的视频免费午夜| 成人性生交大片免费视频hd| 人人妻人人澡人人爽人人夜夜 | 成人欧美大片| 国产精品三级大全| 看黄色毛片网站| 欧美性猛交黑人性爽| 一级毛片aaaaaa免费看小| 天堂动漫精品| 99热精品在线国产| 亚洲精品国产成人久久av| 看片在线看免费视频| 国产成人福利小说| 久久久久久久久中文| 综合色av麻豆| 久久久久久久久久成人| 亚洲aⅴ乱码一区二区在线播放| 黄色一级大片看看| 亚洲婷婷狠狠爱综合网| 你懂的网址亚洲精品在线观看 | 成人午夜高清在线视频| 婷婷精品国产亚洲av| 精品久久久久久久末码| 夜夜看夜夜爽夜夜摸| 国产亚洲av嫩草精品影院| 成人特级黄色片久久久久久久| 最近手机中文字幕大全| 亚洲成人久久性| 最近在线观看免费完整版| 国产高潮美女av| 国产真实乱freesex| 久久热精品热| 麻豆av噜噜一区二区三区| 国产精品无大码| 男人狂女人下面高潮的视频| 国内揄拍国产精品人妻在线| 久久久久精品国产欧美久久久| av免费在线看不卡| 国产探花在线观看一区二区| 国国产精品蜜臀av免费| 日韩在线高清观看一区二区三区| 日本爱情动作片www.在线观看 | 毛片女人毛片| 久久热精品热| 日本黄色片子视频| avwww免费| 日本在线视频免费播放| 国产精品人妻久久久影院| 精品国内亚洲2022精品成人| 噜噜噜噜噜久久久久久91| 亚洲,欧美,日韩| av在线播放精品| 国产黄片美女视频| 男插女下体视频免费在线播放| av.在线天堂| 精品一区二区三区视频在线| av天堂在线播放| 成人鲁丝片一二三区免费| 亚洲精品一区av在线观看| 最新在线观看一区二区三区| 国产不卡一卡二| 亚洲国产精品久久男人天堂| 精品午夜福利在线看| 国产探花在线观看一区二区| 欧美性猛交黑人性爽| 伦精品一区二区三区| 欧美激情在线99| 成人亚洲欧美一区二区av| 免费观看的影片在线观看| av在线蜜桃| 少妇人妻精品综合一区二区 | 香蕉av资源在线| 亚洲在线观看片| 一级毛片电影观看 | 美女高潮的动态| 欧美bdsm另类| 国产高清三级在线| 九九热线精品视视频播放| 国产午夜精品论理片| 一区二区三区四区激情视频 | 亚洲综合色惰| 中文字幕久久专区| 欧美一区二区亚洲| .国产精品久久| 99久国产av精品国产电影| 中文字幕av成人在线电影| 夜夜夜夜夜久久久久| 久久精品久久久久久噜噜老黄 | 午夜福利成人在线免费观看| 精品乱码久久久久久99久播| a级一级毛片免费在线观看| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 国产乱人偷精品视频| 国产精品人妻久久久影院| 小蜜桃在线观看免费完整版高清| 一夜夜www| 变态另类成人亚洲欧美熟女| 别揉我奶头~嗯~啊~动态视频| av在线观看视频网站免费| 成人亚洲精品av一区二区| 免费电影在线观看免费观看| av专区在线播放| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 国模一区二区三区四区视频| ponron亚洲| 久久久国产成人精品二区| 国产伦精品一区二区三区视频9| 淫妇啪啪啪对白视频| a级毛片免费高清观看在线播放| 国产午夜精品论理片| 99久久无色码亚洲精品果冻| 国产免费一级a男人的天堂| a级毛色黄片| 久久亚洲精品不卡| 啦啦啦韩国在线观看视频| 俄罗斯特黄特色一大片| 国产高清不卡午夜福利| 亚洲人与动物交配视频| 日韩中字成人| 黄色欧美视频在线观看| 成人漫画全彩无遮挡| 网址你懂的国产日韩在线| 免费观看的影片在线观看| 日本免费一区二区三区高清不卡| 国产成人aa在线观看| 日韩人妻高清精品专区| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人精品一区久久| 精品一区二区三区av网在线观看| 日韩欧美精品v在线| 色av中文字幕| 在线观看一区二区三区| 在线播放国产精品三级| 久久久久久久久久黄片| 联通29元200g的流量卡| 床上黄色一级片| 日日摸夜夜添夜夜添小说| 精品久久国产蜜桃| 亚洲欧美成人精品一区二区| 免费观看人在逋| 搡老熟女国产l中国老女人| 综合色av麻豆| 久久久久久久久久成人| 99国产精品一区二区蜜桃av| 国产免费男女视频| 一级毛片久久久久久久久女| 亚洲av成人av| 99热这里只有精品一区| avwww免费| 国产91av在线免费观看| 男女下面进入的视频免费午夜| 极品教师在线视频| www日本黄色视频网| 特大巨黑吊av在线直播| www日本黄色视频网| 精品99又大又爽又粗少妇毛片| avwww免费| 尤物成人国产欧美一区二区三区| 九九久久精品国产亚洲av麻豆| 九九爱精品视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产极品精品免费视频能看的| 给我免费播放毛片高清在线观看| 成年女人看的毛片在线观看| 国产精华一区二区三区| 国产午夜福利久久久久久| 国产片特级美女逼逼视频| 国产精品综合久久久久久久免费| 久99久视频精品免费| 日韩av不卡免费在线播放| 一卡2卡三卡四卡精品乱码亚洲| 哪里可以看免费的av片| 日本熟妇午夜| 国产片特级美女逼逼视频| 精华霜和精华液先用哪个| 黄色欧美视频在线观看| 一区二区三区四区激情视频 | 女生性感内裤真人,穿戴方法视频| 久久久久国产精品人妻aⅴ院| 99在线视频只有这里精品首页| 97人妻精品一区二区三区麻豆| 精品人妻一区二区三区麻豆 | 国产成人a区在线观看| 欧美性感艳星| 精品熟女少妇av免费看| 蜜桃亚洲精品一区二区三区| 国产精品人妻久久久影院| 天天一区二区日本电影三级| 老司机影院成人| 18禁黄网站禁片免费观看直播| 日韩大尺度精品在线看网址| 日日摸夜夜添夜夜爱| 国产乱人视频| 级片在线观看| 在线观看美女被高潮喷水网站| 国产极品精品免费视频能看的| 99久国产av精品| 床上黄色一级片| 国语自产精品视频在线第100页| 在线免费观看不下载黄p国产| 日韩欧美一区二区三区在线观看| 欧美高清成人免费视频www| 99riav亚洲国产免费| 午夜福利18| 成人综合一区亚洲| 一个人观看的视频www高清免费观看| 久久热精品热| 精品人妻一区二区三区麻豆 | 精品久久久久久久人妻蜜臀av| 中国美女看黄片| 蜜臀久久99精品久久宅男| 十八禁国产超污无遮挡网站| 天美传媒精品一区二区| 一区二区三区高清视频在线| 国产久久久一区二区三区| 中文字幕熟女人妻在线| 国产视频一区二区在线看| 女人被狂操c到高潮| 哪里可以看免费的av片| .国产精品久久| 桃色一区二区三区在线观看| 波多野结衣高清作品| 免费无遮挡裸体视频| 非洲黑人性xxxx精品又粗又长| 亚洲一级一片aⅴ在线观看| 天天一区二区日本电影三级| 一进一出抽搐动态| 免费电影在线观看免费观看| 天天一区二区日本电影三级| 最近2019中文字幕mv第一页| av免费在线看不卡| 日韩,欧美,国产一区二区三区 | 国产成人影院久久av| 久久国产乱子免费精品| 免费一级毛片在线播放高清视频| 成人无遮挡网站| 日本色播在线视频| 天堂av国产一区二区熟女人妻| 99riav亚洲国产免费| 18禁黄网站禁片免费观看直播| 精品久久久久久久人妻蜜臀av| 亚洲精品国产av成人精品 | 精品人妻一区二区三区麻豆 | 亚洲精华国产精华液的使用体验 | 可以在线观看毛片的网站| 亚洲精品影视一区二区三区av| 2021天堂中文幕一二区在线观| 欧美一级a爱片免费观看看| 丝袜美腿在线中文| 国产男靠女视频免费网站| 精品国内亚洲2022精品成人| 国模一区二区三区四区视频| 美女大奶头视频| 精品久久久久久久人妻蜜臀av| 国产精华一区二区三区| 最新中文字幕久久久久|