• <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
    亚洲片人在线观看| 在线播放国产精品三级| 91麻豆精品激情在线观看国产| 亚洲中文字幕一区二区三区有码在线看| 欧美成人a在线观看| 欧美一区二区亚洲| 亚洲黑人精品在线| 久久香蕉国产精品| 国产欧美日韩一区二区三| 成人三级黄色视频| 国产一区在线观看成人免费| 国产精品亚洲一级av第二区| 精品久久久久久久久久免费视频| 夜夜夜夜夜久久久久| 最后的刺客免费高清国语| 有码 亚洲区| 无遮挡黄片免费观看| 动漫黄色视频在线观看| 成人精品一区二区免费| 亚洲avbb在线观看| 亚洲成人久久爱视频| 亚洲欧美一区二区三区黑人| 丰满人妻一区二区三区视频av | 波多野结衣高清无吗| 婷婷亚洲欧美| 99国产综合亚洲精品| 欧美又色又爽又黄视频| 欧美日韩瑟瑟在线播放| 久久久精品大字幕| 97超级碰碰碰精品色视频在线观看| 欧美黑人巨大hd| 天天添夜夜摸| 欧美黄色片欧美黄色片| 欧美日本亚洲视频在线播放| 精品福利观看| 精品久久久久久久人妻蜜臀av| 免费看日本二区| 一本一本综合久久| а√天堂www在线а√下载| 国产私拍福利视频在线观看| 午夜免费成人在线视频| 亚洲中文字幕日韩| 国产免费男女视频| 国产精品一区二区免费欧美| aaaaa片日本免费| 午夜精品在线福利| 嫩草影院入口| 日韩亚洲欧美综合| 高清毛片免费观看视频网站| 免费搜索国产男女视频| 极品教师在线免费播放| 搡女人真爽免费视频火全软件 | 啦啦啦观看免费观看视频高清| 亚洲不卡免费看| 性欧美人与动物交配| 最后的刺客免费高清国语| 欧美日韩乱码在线| 精品乱码久久久久久99久播| 欧美日本亚洲视频在线播放| 精品一区二区三区人妻视频| 午夜两性在线视频| 老司机深夜福利视频在线观看| 久久精品国产亚洲av涩爱 | 有码 亚洲区| 香蕉av资源在线| 日本与韩国留学比较| 在线观看免费视频日本深夜| 香蕉久久夜色| 精品国内亚洲2022精品成人| 色尼玛亚洲综合影院| 欧美丝袜亚洲另类 | 国产精品亚洲一级av第二区| 欧美精品啪啪一区二区三区| 欧美精品啪啪一区二区三区| netflix在线观看网站| 看黄色毛片网站| 中文字幕人妻熟人妻熟丝袜美 | 欧美又色又爽又黄视频| 久久香蕉精品热| 亚洲精品一卡2卡三卡4卡5卡| 国产精品乱码一区二三区的特点| 久久久久免费精品人妻一区二区| 一区二区三区激情视频| 中文资源天堂在线| 窝窝影院91人妻| 精品乱码久久久久久99久播| avwww免费| 午夜福利18| 亚洲美女黄片视频| 亚洲精品日韩av片在线观看 | 听说在线观看完整版免费高清| 亚洲国产欧洲综合997久久,| 国内揄拍国产精品人妻在线| 精品久久久久久,| 大型黄色视频在线免费观看| 中文亚洲av片在线观看爽| 成人国产综合亚洲| 中文字幕高清在线视频| 国产高清videossex| 免费在线观看亚洲国产| 精品熟女少妇八av免费久了| 一边摸一边抽搐一进一小说| 成人午夜高清在线视频| 国产三级黄色录像| 亚洲18禁久久av| 国产精品免费一区二区三区在线| av国产免费在线观看| 99国产极品粉嫩在线观看| 日韩人妻高清精品专区| 身体一侧抽搐| 精品国产超薄肉色丝袜足j| 黄色日韩在线| 久久精品国产亚洲av涩爱 | 九色国产91popny在线| 在线播放国产精品三级| 色av中文字幕| 久久伊人香网站| 99热这里只有精品一区| 男女那种视频在线观看| 99热这里只有是精品50| www.熟女人妻精品国产| 久久久久精品国产欧美久久久| 日韩欧美国产在线观看| 国产成人a区在线观看| 国产综合懂色| 日本三级黄在线观看| 午夜福利成人在线免费观看| 蜜桃久久精品国产亚洲av| 国产在视频线在精品| 丰满人妻熟妇乱又伦精品不卡| 高清在线国产一区| 国产av一区在线观看免费| 校园春色视频在线观看| 久久精品国产99精品国产亚洲性色| 狂野欧美激情性xxxx| 尤物成人国产欧美一区二区三区| 级片在线观看| 美女 人体艺术 gogo| 国产精品嫩草影院av在线观看 | 国产成人欧美在线观看| 国产97色在线日韩免费| 91麻豆av在线| 国产高清videossex| 国产精品一区二区免费欧美| 日韩欧美国产在线观看| 亚洲av成人精品一区久久| 亚洲不卡免费看| 国产精品 欧美亚洲| svipshipincom国产片| 久久精品国产99精品国产亚洲性色| 国内久久婷婷六月综合欲色啪| 级片在线观看| 亚洲午夜理论影院| 日日夜夜操网爽| 欧美乱色亚洲激情| 久久精品亚洲精品国产色婷小说| 香蕉久久夜色| 人人妻人人看人人澡| 久久精品国产综合久久久| 亚洲最大成人中文| 国产精品一区二区三区四区免费观看 | 亚洲无线在线观看| 深爱激情五月婷婷| a在线观看视频网站| 欧美日韩福利视频一区二区| 精品日产1卡2卡| 两个人看的免费小视频| 最新美女视频免费是黄的| 制服人妻中文乱码| 天堂网av新在线| 国产色婷婷99| 啦啦啦韩国在线观看视频| 久久人妻av系列| 亚洲中文日韩欧美视频| 国产熟女xx| 国产精品亚洲一级av第二区| 日韩欧美精品v在线| 国内精品美女久久久久久| 岛国在线观看网站| 又紧又爽又黄一区二区| 精品国产亚洲在线| 国产乱人视频| 波多野结衣高清无吗| 五月伊人婷婷丁香| 婷婷六月久久综合丁香| www.999成人在线观看| 国产黄a三级三级三级人| 天天添夜夜摸| а√天堂www在线а√下载| 亚洲五月天丁香| 俺也久久电影网| 久9热在线精品视频| 18禁黄网站禁片午夜丰满| 性欧美人与动物交配| 国产 一区 欧美 日韩| 在线国产一区二区在线| 久久精品亚洲精品国产色婷小说| 全区人妻精品视频| 每晚都被弄得嗷嗷叫到高潮| 人人妻人人澡欧美一区二区| 少妇丰满av| 久久久久久久久久黄片| www.999成人在线观看| 亚洲五月天丁香| 国产精品嫩草影院av在线观看 | 99热6这里只有精品| 亚洲欧美精品综合久久99| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美黑人欧美精品刺激| 国产成人福利小说| 嫁个100分男人电影在线观看| 国产成人啪精品午夜网站| 啦啦啦韩国在线观看视频| 欧美高清成人免费视频www| 欧美日韩瑟瑟在线播放| 午夜两性在线视频| 国内少妇人妻偷人精品xxx网站| 亚洲自拍偷在线| 一进一出好大好爽视频| 麻豆国产97在线/欧美| 国产一区二区激情短视频| а√天堂www在线а√下载| 一本久久中文字幕| 久久久久性生活片| 在线播放国产精品三级| 国产激情偷乱视频一区二区| 国产精品1区2区在线观看.| 午夜福利在线在线| 黄片小视频在线播放| 亚洲av日韩精品久久久久久密| 一本久久中文字幕| 欧美3d第一页| 国产久久久一区二区三区| 美女被艹到高潮喷水动态| 岛国在线观看网站| 日韩中文字幕欧美一区二区| 人人妻人人澡欧美一区二区| 在线免费观看的www视频| 偷拍熟女少妇极品色| 亚洲av第一区精品v没综合| 热99re8久久精品国产| 国产成人av教育| 精品久久久久久久久久久久久| 欧美绝顶高潮抽搐喷水| 亚洲成人免费电影在线观看| 日韩精品青青久久久久久| 男女做爰动态图高潮gif福利片| 欧美一区二区国产精品久久精品| 夜夜爽天天搞| 久9热在线精品视频| 国内久久婷婷六月综合欲色啪| 国产精品亚洲美女久久久| 国产高潮美女av| 最新中文字幕久久久久| 国产一区二区三区在线臀色熟女| 男人的好看免费观看在线视频| www国产在线视频色| 免费一级毛片在线播放高清视频| 变态另类成人亚洲欧美熟女| www日本黄色视频网| 综合色av麻豆| 日本一二三区视频观看| 久久久国产成人精品二区| 夜夜躁狠狠躁天天躁| 在线十欧美十亚洲十日本专区| 欧美一区二区精品小视频在线| 3wmmmm亚洲av在线观看| 午夜亚洲福利在线播放| 日本撒尿小便嘘嘘汇集6| 精品一区二区三区av网在线观看| 欧美乱码精品一区二区三区| 一区二区三区激情视频| 亚洲精品日韩av片在线观看 | 国产一区二区在线av高清观看| 一本久久中文字幕| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| 国产一级毛片七仙女欲春2| 国产成人系列免费观看| 99热只有精品国产| 免费观看人在逋| 好男人在线观看高清免费视频| 国内久久婷婷六月综合欲色啪| 激情在线观看视频在线高清| 尤物成人国产欧美一区二区三区| 免费在线观看日本一区| 久久香蕉精品热| 久久这里只有精品中国| 18禁在线播放成人免费| 国产视频内射| 成人精品一区二区免费| 老汉色∧v一级毛片| 国产欧美日韩一区二区三| 国产伦人伦偷精品视频| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品一区二区| 久久精品91蜜桃| avwww免费| 亚洲av一区综合| 欧美色欧美亚洲另类二区| 丁香六月欧美| 日韩欧美 国产精品| 岛国视频午夜一区免费看| 久久久久九九精品影院| 男女那种视频在线观看| 99视频精品全部免费 在线| 九九久久精品国产亚洲av麻豆| 99国产精品一区二区蜜桃av| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看| 国产三级中文精品| 国产一区二区三区在线臀色熟女| 精品一区二区三区视频在线 | 在线a可以看的网站| 男女午夜视频在线观看| 九九热线精品视视频播放| 国产亚洲精品久久久com| 欧美极品一区二区三区四区| 欧美日韩国产亚洲二区| 日本免费一区二区三区高清不卡| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 夜夜躁狠狠躁天天躁| 午夜两性在线视频| 一a级毛片在线观看| 亚洲av电影在线进入| 狂野欧美白嫩少妇大欣赏| 日本三级黄在线观看| 亚洲欧美日韩东京热| 一级a爱片免费观看的视频| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 黄色成人免费大全| 美女免费视频网站| 90打野战视频偷拍视频| 欧美日韩一级在线毛片| 色综合婷婷激情| 在线视频色国产色| 女人十人毛片免费观看3o分钟| 欧美日韩亚洲国产一区二区在线观看| 国产精品国产高清国产av| 欧美在线黄色| 51午夜福利影视在线观看| 熟女少妇亚洲综合色aaa.| 久久精品国产综合久久久| 久久欧美精品欧美久久欧美| 天天添夜夜摸| 午夜精品在线福利| 蜜桃亚洲精品一区二区三区| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 国产精品自产拍在线观看55亚洲| 99久久综合精品五月天人人| 成人特级黄色片久久久久久久| 偷拍熟女少妇极品色| 久久久久久久亚洲中文字幕 | 一级黄色大片毛片| 黄色女人牲交| 成人无遮挡网站| 欧美国产日韩亚洲一区| 久久精品91蜜桃| 中文字幕人妻丝袜一区二区| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 熟女电影av网| 亚洲av美国av| 成年人黄色毛片网站| 免费观看精品视频网站| 欧美一区二区精品小视频在线| 亚洲性夜色夜夜综合| 麻豆国产av国片精品| 美女被艹到高潮喷水动态| 国产精品久久久久久亚洲av鲁大| 噜噜噜噜噜久久久久久91| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 久久婷婷人人爽人人干人人爱| 日本 欧美在线| xxx96com| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 日本撒尿小便嘘嘘汇集6| 无限看片的www在线观看| 国产成人系列免费观看| 欧美激情在线99| 国产精品 欧美亚洲| 午夜两性在线视频| 99精品欧美一区二区三区四区| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻| 亚洲精品在线观看二区| 偷拍熟女少妇极品色| 国产v大片淫在线免费观看| 亚洲色图av天堂| 丰满乱子伦码专区| 精品无人区乱码1区二区| 国产乱人伦免费视频| 男女床上黄色一级片免费看| 俺也久久电影网| 国产在线精品亚洲第一网站| 国产爱豆传媒在线观看| 午夜福利免费观看在线| 国产69精品久久久久777片| 亚洲最大成人中文| www.色视频.com| 男插女下体视频免费在线播放| 国产成人a区在线观看| 国产高清videossex| 久久香蕉精品热| 免费高清视频大片| 国产麻豆成人av免费视频| 伊人久久精品亚洲午夜| 热99在线观看视频| 亚洲成av人片免费观看| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 日韩高清综合在线| 可以在线观看毛片的网站| 欧美日本亚洲视频在线播放| 国产精品久久久久久人妻精品电影| 国产综合懂色| 精品一区二区三区人妻视频| 美女高潮喷水抽搐中文字幕| 可以在线观看的亚洲视频| 日韩大尺度精品在线看网址| 久久香蕉精品热| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播| 久久精品人妻少妇| 欧美黑人巨大hd| 人人妻人人看人人澡| 99久久精品热视频| 熟女人妻精品中文字幕| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 哪里可以看免费的av片| 欧美日韩亚洲国产一区二区在线观看| 久久精品91蜜桃| 深爱激情五月婷婷| 国内精品久久久久精免费| 日韩欧美精品v在线| 毛片女人毛片| 亚洲国产日韩欧美精品在线观看 | av女优亚洲男人天堂| 成人高潮视频无遮挡免费网站| 免费搜索国产男女视频| 最新美女视频免费是黄的| 久久国产乱子伦精品免费另类| 毛片女人毛片| 亚洲最大成人中文| 国产精品三级大全| 日本在线视频免费播放| 国产三级在线视频| 男人的好看免费观看在线视频| 国产午夜精品论理片| 免费高清视频大片| 欧美日韩亚洲国产一区二区在线观看| 日韩人妻高清精品专区| 国产成人影院久久av| 欧美激情久久久久久爽电影| 欧美性猛交黑人性爽| 真实男女啪啪啪动态图| 一a级毛片在线观看| 久99久视频精品免费| 国产精品久久久久久久电影 | 亚洲国产欧美网| 日韩欧美三级三区| 欧美成人免费av一区二区三区| 亚洲精品亚洲一区二区| 午夜福利视频1000在线观看| 757午夜福利合集在线观看| 一个人免费在线观看的高清视频| 观看免费一级毛片| 18禁黄网站禁片免费观看直播| aaaaa片日本免费| 成人特级黄色片久久久久久久| 国产黄色小视频在线观看| 国产色婷婷99| 久久这里只有精品中国| 亚洲成a人片在线一区二区| 亚洲 欧美 日韩 在线 免费| 性色av乱码一区二区三区2| 久久久久久久亚洲中文字幕 | 日韩 欧美 亚洲 中文字幕| 午夜精品久久久久久毛片777| 成年版毛片免费区| 欧美性猛交黑人性爽| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 在线免费观看不下载黄p国产 | 色精品久久人妻99蜜桃| 国产午夜精品论理片| 天堂动漫精品| 亚洲avbb在线观看| 色在线成人网| 亚洲人成电影免费在线| 不卡一级毛片| 日韩高清综合在线| 婷婷精品国产亚洲av| 免费看美女性在线毛片视频| 午夜久久久久精精品| 久久国产乱子伦精品免费另类| 99在线视频只有这里精品首页| 亚洲av五月六月丁香网| 亚洲无线在线观看| 高潮久久久久久久久久久不卡| avwww免费| 亚洲第一欧美日韩一区二区三区| 精品人妻偷拍中文字幕| 日本一二三区视频观看| 国产乱人伦免费视频| 免费大片18禁| 国产精品电影一区二区三区| 国产精品香港三级国产av潘金莲| 国产成人aa在线观看| 午夜两性在线视频| 搡老妇女老女人老熟妇| 在线国产一区二区在线| 男女床上黄色一级片免费看| 亚洲avbb在线观看| 国产精品影院久久| 1024手机看黄色片| 又紧又爽又黄一区二区| 亚洲成人免费电影在线观看| а√天堂www在线а√下载| 日韩欧美在线乱码| 亚洲av成人av| 久久人妻av系列| 搡女人真爽免费视频火全软件 | avwww免费| 亚洲一区高清亚洲精品| 日韩免费av在线播放| 免费搜索国产男女视频| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 亚洲国产欧美网| 亚洲国产欧美人成| 久久久久久久亚洲中文字幕 | 九九久久精品国产亚洲av麻豆| 国产成人啪精品午夜网站| 国产一区二区在线观看日韩 | 99国产精品一区二区三区| 听说在线观看完整版免费高清| 丁香六月欧美| 有码 亚洲区| 国产亚洲精品综合一区在线观看| 嫩草影院精品99| 看片在线看免费视频| 少妇的逼好多水| 成人欧美大片| 久久久久久久午夜电影| 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 欧美最黄视频在线播放免费| 欧美黄色片欧美黄色片| 日本一本二区三区精品| 国产精品免费一区二区三区在线| 国产高清三级在线| 国产av不卡久久| 一区二区三区国产精品乱码| 亚洲精品国产精品久久久不卡| 日日摸夜夜添夜夜添小说| 色综合站精品国产| 中文资源天堂在线| 少妇裸体淫交视频免费看高清| 亚洲一区二区三区不卡视频| 我要搜黄色片| 好看av亚洲va欧美ⅴa在| 国内毛片毛片毛片毛片毛片| 国产伦在线观看视频一区| 国产一区二区三区在线臀色熟女| 午夜视频国产福利| 国产三级黄色录像| 一级毛片女人18水好多| 成年女人永久免费观看视频| 99久久精品一区二区三区| 变态另类丝袜制服| 中文字幕久久专区| 亚洲av电影在线进入| 亚洲成av人片在线播放无| 极品教师在线免费播放| 亚洲一区二区三区不卡视频| 免费大片18禁| 免费看光身美女| 少妇裸体淫交视频免费看高清| 欧美高清成人免费视频www| 欧美性猛交黑人性爽| 日日夜夜操网爽| 国产成人福利小说| 亚洲成av人片在线播放无| 黄色视频,在线免费观看| 久久精品亚洲精品国产色婷小说| 国产精品精品国产色婷婷| 国产午夜精品论理片| 在线观看美女被高潮喷水网站 | 美女cb高潮喷水在线观看| 在线天堂最新版资源| 国产真实乱freesex| 91字幕亚洲| 成年免费大片在线观看| 亚洲精品日韩av片在线观看 | 国产探花极品一区二区| 欧美在线一区亚洲| 国产黄片美女视频| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 欧美乱妇无乱码| 亚洲 国产 在线| 国内精品久久久久久久电影| 国产野战对白在线观看| 精品日产1卡2卡| 久久精品国产清高在天天线| 国产在线精品亚洲第一网站|