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

    Biochemical and Metabolomics on Rice Cultivars

    2019-05-23 01:13:02HurugalavadiPuttaswamaiahRanjithaRameGowdaNagarajappaNethraNarayanappaAmrutaHemanthKumarKandikattu
    Rice Science 2019年3期

    Hurugalavadi Puttaswamaiah Ranjitha, Rame Gowda, Nagarajappa Nethra, Narayanappa Amruta, Hemanth Kumar Kandikattu

    ?

    Biochemical and Metabolomics on Rice Cultivars

    Hurugalavadi Puttaswamaiah Ranjitha1, Rame Gowda1, Nagarajappa Nethra1, Narayanappa Amruta1, Hemanth Kumar Kandikattu2

    ()

    Plant metabolites are important for both plant life and human nutrition. However, the genetic control of plant metabolome remains largely unknown. Here, we performed a genetic analysis of the different rice metabolome and isozymes which are highly versatile and non-destructive as bio-markers. Five isozymes peroxidase, catalase, malate dehydrogenase, alcohol dehydrogenase (ADH), polyphenol oxidase were studied to characterize the thirty rice cultivars and two hybrids KRH-2 and KRH-4 along with their parental lines. Among the zymograms, ADH was found to be useful for the detection of cultivars, like CTH1, IR64, IR30864, with anRvalue of 0.549. The metabolomics of rice cultivars by using gas chromatography coupled with mass spectrometry instrument with selected reaction monitoring mode software identified the 66 metabolites in the rice samples, including amino, organic, fatty acids, alcohols and sugars (mono-/dis-accharides). All metabolites investigated varied significantly among rice samples. Jaya had the higher number of metabolites (15) with a peak for each metabolite, followed by Jyothi (13). This study demonstrated a powerful tool and provided a high-quality data for understanding the plant metabolome and isozymes, which may help bridge the gap between the genome and phenome.

    gas chromatography; mass spectrometry; isozyme; spectroscopy; metabolomics; rice

    Plants are highly enriched in specific metabolites with extensive quantitative and qualitative variation both among and within different plant species. Understanding the genes involved in metabolism and dissection of the metabolic pathway is essential to improve plant adaptation to environmental stress, food quality and crop yield. Electrophoresis is a technique used to separate charged particles under the influence of the electric field. It is most useful for separating mixtures of large molecules such as proteins and nucleic acids, which in solution behave as charged particles. The use of proteins and isozymes as genetic markers has been accepted as a reliable tool. Glaszman et al (1988) studied the electrophoretic variation of isozymes in plumules and coleoptile of 4- to 10-day-old seedlings of rice genotypes using starch gel electrophoresis and developed a key for identification of 76 alleles of 24 loci. Hussain et al(1989) characterized 11 rice genotypes for specific protein pattern using polyacrylamide gel electrophoresis (PAGE). Band variation is conspicuous in slow migration regions on the gel. Yu and Wang (1989) conducted a study on esterase isozyme pattern of 84 rice varieties and reported that typical sickle-shaped rice mostly belonged to thegroup. Isozyme patterns of esterase, malate dehydrogenase, alcohol dehydrogenase, peroxidase, catalase and acid phosphatase are compared to distinguish cotton cultivars. The esterase zymogram reveals a unique banding pattern for each genotype. Peroxidase and catalase enzymes are useful in distinguishing the species group and interspecific hybrids, whereas alcohol dehydrogenase polymorphism is not sufficient to aid in cultivar identification (Nerkar and Rao, 1993). Shion et al (1994) revealed the genetic nature of several isozymes segregating in Mendelian inheritance, which showed have simple genetic control, therefore, isozyme variation can be used as genetic markers.

    Metabolomics is the scientific study of chemical processes involving metabolites. Specifically, metabolomics is the ‘systematic study of the unique chemical fingerprints that specific cellular processes leave behind’, and the study of their small-molecule metabolite profiles (Daviss, 2005; Gong et al, 2013). Metabolomics is a high-throughput platform which is used to profile the small molecules or metabolites involved in metabolic processes. Mass spectrometry (MS) is used to identify and quantify metabolites after optional separation by Gas chromatography (GC), high performance liquid chromatography (HPLC), liquid chromatography-mass spectrometry (LC-MS), and is the first hyphenated technique to be developed. Identification leverages the distinct patterns in which analytes fragment can be as a mass spectral fingerprint; libraries exist to allow identification of a metabolite according to this fragmentation pattern. MS is sensitive and can be very specific as well. There are also a number of techniques which use MS as a stand-alone technology: the sample is infused directly into the mass spectrometer with no prior separation, and MS provides sufficient selectivity to separate and detect metabolites. The metabolomics study of rice may be used to ‘fingerprint’ different rice samples in order to identify strains that hold significant promise for development, like rice quality improvement and increased nutritional content. As crude extracts tend to be used for metabolite profiling in a nontargeted manner, chromatographic separation is usually required (Fukushima and Kusano, 2014), and fractionated compounds were analyzed. Gas chromatography coupled with mass spectrometry (GC-MS) is particularly well suited to the study of low molecular-weight metabolites that can be made amenable to gas chromatography by chemical derivitization. The main advantage of GC-MS lies in its capability for both metabolite profiling and targeted metabolite quantitation for the analysis of hundreds of metabolites in one GC-MS/MS multiple reaction monitoring (MRM) run. The MRM acquisition mode of unique and timed precursor ion-product, ion transitions for metabolites of interest eliminates interference from matrix or noise. Likewise, co-eluting compounds can be quantified separately, forgoing the need for de-convolution by mathematical methods.

    MATERIALS AND METHODS

    Isozyme analysis

    Pure rice seeds of thirty cultivars and two hybrids (KRH-2 and KRH-4) along with their parental lines (IR58025 and KMR-3 for KRH-2, and CMRS32A and MSN36 for KRH-4) were obtained from Rice Breeding Unit, Zonal Agricultural Research Station, VC Farm, Mandya, Karnataka, India. The seeds were cleaned, dried to the safe moisture level and preserved for further use. Five isozymes peroxidase (POD), catalase (CAT), malate dehydrogenase (MDH), alcohol dehydrogenase (ADH), polyphenol oxidase (PPO) were studied.

    Five seeds were soaked in water under an aerobic condition for 48 h and then for germination. Four-day-old single seedling was ground thoroughly in pestle and mortar along with 200 μL extraction buffer, and then a pinch of poly vinyl pyrrolidone (PVP) was added. The extract was taken in 2 mL centrifuge tube and incubated at 40 oC for extraction of the enzyme for 2 h. After 2 h, sample was centrifuge at 40 oC at 12 000 ×for 15 min. The supernatant was collected and 10 μL of tracking dye (1% bromophenol blue and drop of glycerol) was added and again centrifuged at 40 oC at 12 000 ×for 5 min, and 30 μL supernatant was loaded into the gel for electrophoresis. Resolving stacking gel and staining solutions were prepared by using the protocol of Glaszman et al (1988). Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was performed according to Laemmli (1970).

    TheRvalue was calculated asR=D/D, whereDis the distance traveled by the band, andDis the distance traveled by the tracking dye. Bands were numbered on the basis of increasingRvalues, the presence or absence of specific bands or combination of different bands.

    Metabolomics analysis

    Twelve cultivars BR2655, Thanu, Tunga, BPT5204, IR64, Raksha, Rasi, Jaya, Jyothi, CTH1, KRH-4 and JGL1798 were selected to metabolomics analysis based on their popularity and usage by the public.

    Derivatives of rice extracts were prepared for metabolic profiling according to Zhou et al (2009). Freshly dehulled grains were immediately dipped in liquid nitrogen prior to grinding the rice samples using a 60-mesh sieve. A total of 300 mg rice samples were weighed, and mixed with 5 mL methanol/water (4:1) with 120 μL cupric acid solution (0.30 mg/mL) as an internal standard. The sample vial was vortexed for 5 min and let stand for 45 min followed by 60 min of sonication before centrifuging at 12 000 ×for 10 min. Then, 7 mL supernatant was removed, freeze-dried and subsequently derivatized using 100 μL N,O-bis (trimethylsilyl) tri-fluoro- acetamide (BSTFA) with 3% trimethyl-chlorosilane (TMCS) and 80 μL pyridine. The mixture was heated in a water bath at 75 oC for 45 min before being transferred to 2 mL GC vials with 200 μL inserts, ready for analysis by GC-MS/MS. Steps were followed to ensure that the derivatized samples were analyzed within 24 h after derivatization.

    Fig. 1. Alcohol dehydrogenase activity of rice cultivars.

    1, Jaya; 2, Mandya Sona 1; 3, CTH1; 4, BI33; 5, Ratna Choodi; 6, KMP200; 7, MAS946-1; 8, Rajamudi; 9, JGL1798; 10, KMO201; 11, Mandya Sona 2; 12, MAS26; 13, BR2655; 14, BPT5204; 15, KMP175; 16, Mandya Sona 3; 17, Thanu; 18, KMP128; 19, MTU1001; 20, Gangavathi Sona; 21, Rasi; 22, Raksha; 23, Tella Hamsa; 24, MTU1001; 25, IR30864; 26, IR64; 27, Jyothi; 28, CTH3; 29, KMP153; 30, Tunga; 2A, IR58025; 2H, KRH-2; 2R, KMR-3; 4A, CRMS32A; 4H, KRH-4; 4R, MSN36.

    Phase-I requires full scan GC-MS analysis for identification of compounds by searching library databases. Phase-II requires analysis in the selected reaction monitoring mode (SRM) in order to accurately identify the compounds at trace levels in the heavy matrix background. The ability to perform simultaneously full scan and SRM analyses in a single sample run allows the users to quantitate the known metabolites (targeted) while allowing identification of unknowns (discovery). One of these samples was chosen for the development and optimization of complete SRM. The compound retention time was constant for both workflow stages, served as an additional qualifier for identification, and greatly facilitated the setup of the targeted MRM method.

    RESULTS AND DISCUSSION

    Isozyme analysis

    We selected the five isozymes mainly because they are the major enzymes found in almost all the rice cultivars. And among the isozymes, ADH was effective in distinguishing the cultivars CTH1, IR30864 and IR64 with anRvalue of 0.549, while the other isozymes did not show any significant differences between the cultivars.

    Among the different isozymes (Fig. 1 and Supplemental Fig. 1), ADH showed different banding pattern or difference in the Rm value. The zymograms of ADH are presented in Fig. 1. The electrophoretic analysis of rice cultivars for ADH revealed the presence of only two bands. Band number one with medium intensity was present in most of the cultivars, except Raksha, Tella Hamsa, MTU1010, IR30864, Jyothi, CTH3, KMP153, and Tunga showing very high-intensity bands, with anRvalue of 0.427. Whereas, this band was absent in CTH1 and IR64. Band number two of high intensity was absent in all the cultivars studied except CTH1, IR30864, and IR64 with anRvalue of 0.549. Therefore, by using ADH enzyme, these cultivars can be differentiated and the zymogram profile was showed in the Fig. 1-A.

    The electrophoretic analysis of rice parental lines and hybrid for ADH revealed the presence of only four bands (Fig. 1-B). Band number one with high intensity was present in the hybrid KRH-2 and its parental lines with anRvalue of 0.398. The difference in terms of relative intensity among CRMS32A, MSN36 and KRH-4 was observed in the band number one. KRH-4 exhibited dark band intensity whereas light intensity band was observed in both CRMS32A and MSN36. It was also present with very high intensity in KRH-4 and medium intensity in its parental lines. Band number two was present in KRH-2 with low intensity and in its parental lines with very low intensity (Rof 0.439). In case of KRH-4 and its parent lines, the band was low intensity. Band number three was present in hybrid KRH-2 and its parental lines with medium intensity. In KRH-4, high-intensity band was present and whereas in its parental lines medium intensity band was present with anRvalue 0.464. Band number four was present in KRH2 hybrid and its parental lines but it was very low in intensity. KRH-4 and its parental lines also showed a very low-intensity band with anRvalue of 0.509. Band number five was low intensity in IR58025 line and it was very low in hybrid KRH-2 and KMR-3. KRH-4 and CRMS32A hybrid showed low-intensity band and MSN36 had shown very low-intensity band with anRvalue of 0.578 (Rohini, 2000; Nethra, 2007; Bimb et al,2010; Rao et al, 2012; Pushpa et al, 2014; Venuprasad et al, 2014; Galal et al, 2015).

    Metabolomics analysis

    Rice produced various types of metabolites which may be major bioactive (Supplemental Table 1). Phytohormones were secondary metabolites that can coordinate cellular activities for germination, growth, protection against stress, and other physiological activities in the life of plants. As crude extracts tend to be used for metabolite profiling in a non-targeted manner, chromatographic separation is usually required (Fukushima and Kusano, 2014), and fractionated compounds are then analyzed by MS. These analytical methods are most commonly used in not only metabolomics (Kim et al, 2011), but also natural product chemistry. We analyzed the metabolite composition in the rice cultivars and also discussed current knowledge on the secondary metabolites produced by rice and their biological activities as well as the strategies and challenges for using constituents. It captured differences in the metabolite abundance and correlations between rice cultivars, suggesting that these differences may reflect the extent of metabolic adaptation and the specific phenotypes in rice genotypes.

    During phase I of the workflow, five different rice samples KRH-4, Jyothi, CTH1, Raksha, Jaya were analyzed. These samples were chosen for their diverse characteristics and the likelihood they would display significant differences in the amounts of various metabolites.

    Precursor ion study from full scan data, 66 interest metabolites of and 3–4 precursor ions were chosen for each compound. This step required a single sample injection GC-MS/MS run and the selected precursor ions were then taken to the product ion study (Kandikattu et al, 2014, 2015). Product ions were selected for each of the precursor ions. Several product ions were chosen for each compound, resulting in 3–7 precursor-product ion pairs per compound, which were then taken for the SRM optimization study.

    Fig. 2. Selected reaction monitoring mode comparison heat map for 66 compounds in 12 tested rice cultivars.

    The colors in the heat map represent the relative concentration with green and red corresponding to 0 and 100%, respectively. The heat map allows direct visual comparison of the relative amounts of these compounds in the 12 tested samples.

    The SRM transitions from the Auto SRM optimization described above were exported to the Trace Finder quantitation software and a processing method was automatically created. KRH-4, Jyothi, CTH1, Raksha, Jaya as then run again in SRM mode to check the validity of the SRM optimization. The optimized SRM method used to analyze the twelve different rice strains to identify the potential biomarker, compounds that could be used to discriminate between rice accessions. The relative number of individual metabolites can be represented in the form of the heat map.

    The heat map showed the relative concentration of the 66 analysts in the 12 rice samples (Fig. 2). The peak areas were internal standard-corrected, and the final peak areas were normalized. In this way, the sample with the highest concentration of a particular analyst was assigned a relative value of 100% and all other sample concentrations were normalized as a fraction of this. The colors in the heat map represent the relative concentration with green and red corresponding to 0 and 100%, respectively (Fig. 2). The heat map allows direct visual comparison of the relative amounts of these compounds in the 12 tested samples.

    The heat map showed that the relative concentration of the 66 metabolites analyzed varied widely among the 12 tested samples. Jaya showed the highest traces (15) the most of metabolites detected, followed by Jyothi (13) and KRH-4 (12), while the remaining cultivars showed low to medium traces for the metabolites (Fig. 2 and Supplemental Table 1). Based on the peak values, retention time and amount of traces of metabolites in particular cultivar, the cultivars can be easily identified, which is helpful to distinguish the cultivars from one other. In fact, all the metabolites studied are potential biomarkers and may be used to distinguish between different rice accessions. This quantitative metabolic data can be used with principal component analysis to associate metabolic biomarkers with grain quality features of the contrasting rice, and National Institute of Standards and Technologylibrary search was used to identify unknowns with aim of identifying potential biomarkers for characterization in support of grain quality research. The identification of unknown compounds could be achieved by library search with high confidence. Isomeric compounds have been distinguished on the basis of their elution at differing retention times and could be selectively quantified based on their specific fragmentation profiles. The 66 metabolites in the rice samples includes amino, organic and fatty acids, alcohols and sugars (mono-/dis- accharides). The optimized SRM method was completed in less than 24 h and used to analyze these metabolites in 12 diverse rice samples. All metabolites investigated varied significantly among rice samples and therefore have the potential for use as biomarkers for the characterization of diverse rice accessions. This work forms a basis for future studies to tap into the value of diversity to improve grain quality and nutritional value of rice in the larger panels (Wechwerth et al, 2004; Zhou, 2009; Kusano et al, 2015; Nam et al, 2015).

    Acknowledgements

    This study was supported by the Rajiv Gandhi National Fellowship of the Ministry of Science and Technology, India.

    SUPPLEMENTAL DATA

    The following materials are available in the online version of this article at http://www.sciencedirect.com/science/journal/ 16726308; http://www.ricescience.org.

    Supplemental Table 1. Metabolites obtained from GC/MS SRM transitions method for 12 diverse rice cultivars selected from Phase-I discovery.

    Supplemental Fig. 1. Malate dehydrogenase (A), catalase (B), peroxidase (C), esterase (D), and phenylalanine ammonia lyase (E) activities of thirty rice cultivars.

    Bimb H P, Sah R P, Karn N L. 2010. Isozyme variations in fine and aromatic rice genotypes.,5: 59–66.

    Daviss B. 2005. Growing pains for metabolomics., 19(8): 25–28.

    Fukushima A, Kusano M. 2014. A network perspective on nitrogen metabolism from model to crop plants using integrated ‘omics’ approaches., 65: 5619–5630.

    Galal O A, Abo-Youssef M I, Abdelaziz M, Gharib A T, Dora S A. 2015. Assessment of genetic purity of some hybrid rice parental lines using protein profile and fertility restorer gene linked markers.,2(6):75–88.

    Glaszman J C, de Los Reyes B G, Khush G S. 1988. Electrophoretic variation of isozymes in plumules of rice (L.) a key to the identification of 76 alleles at 24 loci., 134: 1–3.

    Gong L, Chen W, Gao Y Q, Liu X Q, Zhang H Y, Xu C G, Yu S B, Zhang Q F, Luo J. 2013. Genetic analysis of the metabolome exemplified using a rice population., 110: 20320–20325.

    Hussain A, Scanlon M G, Juliano B O, Bushu K W. 1989. Discrimination of rice cultivars by PAGE and HPLC.,66(4): 353–356.

    Kandikattu H K, Venuprasad M P, Pal A, Khanum F. 2014. Phytochemical analysis and exercise enhancing effects of hydroalcoholic extract of., 55: 217–224.

    Kandikattu H K, Rachitha P, Krupashree K, Jayashree G V, Abhishek V, Khanum F. 2015. LC-ESI-MS/MS analysis of total oligomeric flavonoid fraction ofand its antioxidant, macromolecule damage protective and antihemolytic effects., 22(4): 165–173.

    Kim H K, Choi Y H, Verpoorte R. 2011. NMR-based plant metabolomics: Where do we stand, where do we go?, 29(6): 267–275.

    Kusano M, Yang Z, Okazaki Y, Nakabayashi R, Fukushima A, Saito K. 2015. Using metabolomic approaches to explore chemical diversity in rice., 8(1): 58–67.

    Laemmli U K. 1970. Cleavage of structural proteins during the assembly of the head of bacteriophage., 227: 680–685.

    Nam M H, Bang E, Kwon T Y, Kim Y, Kim E H, Cho K, Park W J, Kim B G, Yoon I S. 2015. Metabolite profiling of diverse rice germplasm and identification of conserved metabolic markers of rice roots in response to long-term mild salinity stress., 16: 21959–21974.

    Nerkar M, Rao G N. 1993. Studies on isozyme variability of rice (L.)., 19(3): 477–480.

    Nethra N. 2007 Studies on varietal characterization based on morphological, biochemical and molecular markers in rice. [Master Thesis]. Bangalore, India: University Agriculture Science of Bangalore.

    Pushpa C, Gowda R, Nethra N, Nataraj K, Uma Rani K, Gangaraju N. 2014. Identification of parental lines and rice hybrid (KRH-4) using protein and isozyme electrophoresis.: 2328–3491.

    Rao P, Bharathi M, Reddy K, Keshavulu K, Rao L, Neeraja C. 2012. Varietal identification in rice () through chemical tests and gel electrophoresis of soluble seed proteins., 82(4): 304–311.

    Rohini D D. 2000. Studies on the characterization of varieties based on morphological and biochemical traits in rice (L.). [Master Thesis]. Bangalore, India: University Agriculture Science of Bangalore.

    Shion B B, Morishma H, Oka H I. 1994. A survey of variations in peroxidase, acid phosphatase and esterase isozymes of wild and cultivatedsp.,44(5): 303–319.

    Venuprasad M P, Kandikattu H K, Razack S, Khanum F. 2014. Phytochemical analysis ofby LC-ESI-MS/MS and its antioxidant and anxiolytic effects., 92: 151–158.

    Weckwerth W, Wenzel K, Fiehn O. 2004. Process for the integrated extraction identification, and quantification of metabolites, proteins and RNA to reveal their co-regulation in biochemical networks., 4(1): 78–83.

    Yu W J, Wang X. 1989. Classification of sickle-shaped rice (L.) based on esterase isoenzymes.,3: 4–7. (in Chinese)

    Zhou J, Ma C F, Xu H L, Yuan K L, Lu X, Zhu Z, Wu Y N, Xu G W. 2009. Metabolic profiling of transgenic rice withandgenes: An evaluation of unintended effects at metabolic level by using GC-FID and GC-MS., 877: 725–732.

    Copyright ? 2019, China National Rice Research Institute. Hosting by Elsevier B V

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

    Peer review under responsibility of China National Rice Research Institute

    http://dx.doi.org/10.1016/j.rsci.2018.08.007

    8 March 2018;

    28 August 2018

    Hurugalavadi Puttaswamaiah Ranjitha (ranjuranjitha91@gmail.com)

    (Managing Editor: Wang Caihong)

    日韩大尺度精品在线看网址| 成人性生交大片免费视频hd| 色视频www国产| 首页视频小说图片口味搜索| 国产高清三级在线| 久久久久久久久久成人| 免费av不卡在线播放| av天堂在线播放| 日本 av在线| 99热只有精品国产| 久久精品综合一区二区三区| 国产黄色小视频在线观看| 男人的好看免费观看在线视频| 欧美精品国产亚洲| 亚洲性夜色夜夜综合| 日韩高清综合在线| 成年女人毛片免费观看观看9| 国产成人啪精品午夜网站| 国内精品久久久久精免费| 伦理电影大哥的女人| 嫩草影视91久久| 国产精品免费一区二区三区在线| 露出奶头的视频| a在线观看视频网站| 免费看日本二区| 国产欧美日韩精品一区二区| 国产乱人伦免费视频| 最好的美女福利视频网| 一二三四社区在线视频社区8| 欧美日韩乱码在线| 麻豆国产av国片精品| 亚洲avbb在线观看| 在线国产一区二区在线| 在线观看66精品国产| 波多野结衣高清无吗| 欧美黄色片欧美黄色片| 国产精品久久久久久人妻精品电影| 麻豆国产av国片精品| 蜜桃久久精品国产亚洲av| 成人精品一区二区免费| 国产精品亚洲一级av第二区| 亚洲成人精品中文字幕电影| 亚洲熟妇熟女久久| 欧美成人a在线观看| 热99re8久久精品国产| 精品久久久久久久久久久久久| 色综合婷婷激情| 在线观看免费视频日本深夜| 国产av在哪里看| 五月玫瑰六月丁香| 一个人观看的视频www高清免费观看| 亚洲av免费高清在线观看| 亚洲av不卡在线观看| 99久久99久久久精品蜜桃| 波野结衣二区三区在线| 国内少妇人妻偷人精品xxx网站| 日韩精品青青久久久久久| 国产主播在线观看一区二区| 有码 亚洲区| 美女xxoo啪啪120秒动态图 | 最近最新免费中文字幕在线| 国产亚洲av嫩草精品影院| 狂野欧美白嫩少妇大欣赏| 国产中年淑女户外野战色| 国内久久婷婷六月综合欲色啪| 99热精品在线国产| .国产精品久久| 一级a爱片免费观看的视频| 欧美黑人欧美精品刺激| 亚洲国产精品成人综合色| x7x7x7水蜜桃| 久久国产精品人妻蜜桃| 激情在线观看视频在线高清| 自拍偷自拍亚洲精品老妇| АⅤ资源中文在线天堂| 欧美成人一区二区免费高清观看| 少妇人妻精品综合一区二区 | 男女那种视频在线观看| 国产精品国产高清国产av| 精品人妻1区二区| 十八禁网站免费在线| 女人十人毛片免费观看3o分钟| 一级作爱视频免费观看| 一本精品99久久精品77| 婷婷亚洲欧美| 村上凉子中文字幕在线| 有码 亚洲区| 欧美成人一区二区免费高清观看| 亚洲人成伊人成综合网2020| 99热这里只有是精品50| 三级男女做爰猛烈吃奶摸视频| av在线观看视频网站免费| 99在线人妻在线中文字幕| 97热精品久久久久久| 精品午夜福利视频在线观看一区| 亚洲人成网站高清观看| 国产探花极品一区二区| netflix在线观看网站| 免费无遮挡裸体视频| 日本一二三区视频观看| 亚洲成av人片免费观看| 一卡2卡三卡四卡精品乱码亚洲| av在线老鸭窝| 少妇高潮的动态图| 国产野战对白在线观看| 国产精品三级大全| 高清毛片免费观看视频网站| 别揉我奶头 嗯啊视频| 欧美精品啪啪一区二区三区| 麻豆成人av在线观看| 午夜福利18| 长腿黑丝高跟| 色播亚洲综合网| 欧美黄色淫秽网站| 国产亚洲精品av在线| 国产乱人视频| 欧美黑人欧美精品刺激| 精品久久久久久久末码| 老司机福利观看| 成人永久免费在线观看视频| 婷婷精品国产亚洲av在线| 九色成人免费人妻av| 毛片女人毛片| 亚洲在线自拍视频| 一区二区三区激情视频| 2021天堂中文幕一二区在线观| av国产免费在线观看| 成年女人永久免费观看视频| 亚洲国产色片| 人妻丰满熟妇av一区二区三区| 在线十欧美十亚洲十日本专区| 此物有八面人人有两片| 国产亚洲欧美在线一区二区| 一级作爱视频免费观看| 怎么达到女性高潮| 999久久久精品免费观看国产| www日本黄色视频网| 搡老妇女老女人老熟妇| 性欧美人与动物交配| 国产伦在线观看视频一区| 国产精品影院久久| a级毛片a级免费在线| 国产成人av教育| a级一级毛片免费在线观看| 不卡一级毛片| 一二三四社区在线视频社区8| 午夜老司机福利剧场| 国产在视频线在精品| 小说图片视频综合网站| 高清日韩中文字幕在线| 日本 欧美在线| 日本在线视频免费播放| 亚洲经典国产精华液单 | 99久久精品国产亚洲精品| 国产激情偷乱视频一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产国拍精品亚洲av在线观看| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 亚洲成人免费电影在线观看| 国产精品,欧美在线| 俺也久久电影网| 国产高清三级在线| av在线蜜桃| 亚洲熟妇熟女久久| 国产成人aa在线观看| 天堂av国产一区二区熟女人妻| 国产高清视频在线播放一区| 久久久久久国产a免费观看| 一级黄片播放器| 欧美激情国产日韩精品一区| 国内揄拍国产精品人妻在线| 99riav亚洲国产免费| 午夜精品一区二区三区免费看| 伦理电影大哥的女人| 精品久久久久久,| 国产黄a三级三级三级人| 97热精品久久久久久| 久久人妻av系列| 日韩欧美 国产精品| 国产三级黄色录像| 757午夜福利合集在线观看| 国产午夜精品论理片| 欧美bdsm另类| 18禁裸乳无遮挡免费网站照片| 色尼玛亚洲综合影院| 欧美区成人在线视频| 免费高清视频大片| 亚洲人与动物交配视频| 国产精品99久久久久久久久| 国产主播在线观看一区二区| 禁无遮挡网站| 人人妻人人看人人澡| 日韩欧美精品v在线| 夜夜爽天天搞| 又黄又爽又刺激的免费视频.| 国产一级毛片七仙女欲春2| 国产三级黄色录像| 国产69精品久久久久777片| 伦理电影大哥的女人| 日韩欧美三级三区| 三级毛片av免费| 熟妇人妻久久中文字幕3abv| av国产免费在线观看| 在线观看午夜福利视频| 成人欧美大片| 91av网一区二区| a在线观看视频网站| 天堂av国产一区二区熟女人妻| 精品人妻偷拍中文字幕| 俄罗斯特黄特色一大片| 亚洲久久久久久中文字幕| 久久精品91蜜桃| 久久精品综合一区二区三区| 国产精品久久久久久久久免 | 我的女老师完整版在线观看| 能在线免费观看的黄片| 搞女人的毛片| 亚洲激情在线av| 亚洲国产色片| 国产午夜精品论理片| 一区福利在线观看| 男人的好看免费观看在线视频| 级片在线观看| 免费av不卡在线播放| 男女床上黄色一级片免费看| 日韩 亚洲 欧美在线| av国产免费在线观看| 床上黄色一级片| 午夜激情福利司机影院| 欧美激情国产日韩精品一区| 丰满的人妻完整版| 久久久久九九精品影院| 成人一区二区视频在线观看| 97碰自拍视频| 一个人免费在线观看电影| 蜜桃久久精品国产亚洲av| 成人特级黄色片久久久久久久| 久久久成人免费电影| 男人和女人高潮做爰伦理| 国产在线男女| 精品久久久久久,| 毛片一级片免费看久久久久 | 在线看三级毛片| 特级一级黄色大片| 亚洲美女视频黄频| 给我免费播放毛片高清在线观看| 久久久久久久亚洲中文字幕 | 亚洲av免费高清在线观看| 丁香六月欧美| av在线观看视频网站免费| 99国产精品一区二区蜜桃av| 级片在线观看| 99久久久亚洲精品蜜臀av| 亚洲性夜色夜夜综合| 午夜福利成人在线免费观看| 中文亚洲av片在线观看爽| 国产三级黄色录像| 国产成人欧美在线观看| 国产精品98久久久久久宅男小说| 国产免费av片在线观看野外av| 成人无遮挡网站| 免费人成在线观看视频色| 2021天堂中文幕一二区在线观| 制服丝袜大香蕉在线| 国产久久久一区二区三区| 国内精品美女久久久久久| 亚洲最大成人av| 免费黄网站久久成人精品 | 在线观看午夜福利视频| 美女被艹到高潮喷水动态| 久久国产乱子伦精品免费另类| 岛国在线免费视频观看| 最近最新中文字幕大全电影3| 精品国产亚洲在线| 亚洲无线在线观看| 99久国产av精品| 在线观看免费视频日本深夜| 琪琪午夜伦伦电影理论片6080| 桃色一区二区三区在线观看| 精品福利观看| 国产高清视频在线观看网站| 美女高潮的动态| 一级黄色大片毛片| 国产黄a三级三级三级人| 男人狂女人下面高潮的视频| 小说图片视频综合网站| 久久久久久国产a免费观看| 精品午夜福利视频在线观看一区| 国产精品98久久久久久宅男小说| 丰满人妻一区二区三区视频av| 欧美色视频一区免费| 在线国产一区二区在线| 热99在线观看视频| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 禁无遮挡网站| 国产 一区 欧美 日韩| 91字幕亚洲| 岛国在线免费视频观看| 亚洲欧美日韩高清在线视频| 日本与韩国留学比较| 亚洲精品久久国产高清桃花| 成人精品一区二区免费| 亚洲片人在线观看| 一级作爱视频免费观看| 少妇熟女aⅴ在线视频| 97超视频在线观看视频| 99久久精品国产亚洲精品| 久久国产精品影院| 亚洲狠狠婷婷综合久久图片| 国内久久婷婷六月综合欲色啪| 日本成人三级电影网站| 国产精品爽爽va在线观看网站| 三级毛片av免费| 看十八女毛片水多多多| 麻豆国产97在线/欧美| 成熟少妇高潮喷水视频| 12—13女人毛片做爰片一| 两个人视频免费观看高清| 亚洲欧美日韩高清在线视频| 亚洲精品日韩av片在线观看| 国产亚洲精品久久久com| 欧美日韩国产亚洲二区| 精品乱码久久久久久99久播| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 国产伦人伦偷精品视频| 亚洲欧美日韩高清在线视频| 久久国产乱子免费精品| 久久九九热精品免费| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区久久| 久久久久性生活片| 国产精品久久久久久久电影| 18禁在线播放成人免费| 亚洲av成人精品一区久久| 美女大奶头视频| 一本一本综合久久| 欧美激情久久久久久爽电影| 国产亚洲欧美在线一区二区| 久久伊人香网站| 亚洲真实伦在线观看| 俺也久久电影网| 免费看光身美女| 九色成人免费人妻av| 亚洲久久久久久中文字幕| 久久欧美精品欧美久久欧美| 俺也久久电影网| 91九色精品人成在线观看| 日韩欧美三级三区| 中出人妻视频一区二区| 1000部很黄的大片| 成人午夜高清在线视频| 最近最新中文字幕大全电影3| 亚洲精品在线美女| 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 欧美乱色亚洲激情| 天堂√8在线中文| 久久伊人香网站| av在线观看视频网站免费| 午夜精品在线福利| 婷婷亚洲欧美| 露出奶头的视频| 老女人水多毛片| 91av网一区二区| 亚洲成人久久爱视频| 黄色日韩在线| 一级毛片久久久久久久久女| 怎么达到女性高潮| 一区福利在线观看| 欧美色欧美亚洲另类二区| 成年女人看的毛片在线观看| 国产高清有码在线观看视频| 中文字幕高清在线视频| 亚洲av.av天堂| 国产精品亚洲av一区麻豆| 99久久久亚洲精品蜜臀av| 精品久久久久久久末码| 人人妻人人看人人澡| 国产欧美日韩一区二区三| 国内精品久久久久久久电影| 国产亚洲精品久久久久久毛片| 欧美黄色片欧美黄色片| 亚洲精品456在线播放app | 国产色婷婷99| 色噜噜av男人的天堂激情| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产成人福利小说| 色av中文字幕| 日本免费一区二区三区高清不卡| 日韩av在线大香蕉| 亚洲成av人片在线播放无| 亚洲18禁久久av| 欧美成人免费av一区二区三区| 12—13女人毛片做爰片一| 欧美xxxx黑人xx丫x性爽| 搞女人的毛片| 国产精品三级大全| 国产在视频线在精品| 亚洲自拍偷在线| 国产淫片久久久久久久久 | 久久国产乱子伦精品免费另类| 俺也久久电影网| 国产精品乱码一区二三区的特点| 丰满人妻一区二区三区视频av| 国产中年淑女户外野战色| 国产毛片a区久久久久| 九色成人免费人妻av| 国产久久久一区二区三区| 免费看a级黄色片| 亚洲av免费高清在线观看| 白带黄色成豆腐渣| 嫩草影院精品99| 国产色婷婷99| 国产免费av片在线观看野外av| 美女 人体艺术 gogo| 国产老妇女一区| 搡老妇女老女人老熟妇| 少妇的逼水好多| 99在线人妻在线中文字幕| 又黄又爽又刺激的免费视频.| 精品久久久久久,| 老司机午夜十八禁免费视频| 亚洲成人久久性| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 国产国拍精品亚洲av在线观看| 欧美高清性xxxxhd video| 无遮挡黄片免费观看| 热99re8久久精品国产| 我的老师免费观看完整版| 最后的刺客免费高清国语| 欧美不卡视频在线免费观看| 我要搜黄色片| 国产精品不卡视频一区二区 | 亚洲天堂国产精品一区在线| 久久精品久久久久久噜噜老黄 | 成年人黄色毛片网站| 久久国产精品影院| 一个人观看的视频www高清免费观看| 欧美日本视频| a级毛片免费高清观看在线播放| 中文在线观看免费www的网站| 美女免费视频网站| 亚洲欧美激情综合另类| ponron亚洲| 久久久国产成人免费| 亚洲电影在线观看av| 中国美女看黄片| 国产精品98久久久久久宅男小说| 亚洲av一区综合| 国产亚洲精品av在线| 淫秽高清视频在线观看| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 国产aⅴ精品一区二区三区波| 色噜噜av男人的天堂激情| 国产aⅴ精品一区二区三区波| 国产精品人妻久久久久久| 国产欧美日韩精品亚洲av| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 欧美中文日本在线观看视频| 人人妻,人人澡人人爽秒播| 国产麻豆成人av免费视频| 我的老师免费观看完整版| 午夜精品在线福利| 日本黄色视频三级网站网址| 国模一区二区三区四区视频| 成年版毛片免费区| 女人十人毛片免费观看3o分钟| 免费av观看视频| 国产精品爽爽va在线观看网站| 国内精品久久久久精免费| 日本精品一区二区三区蜜桃| 男女那种视频在线观看| 88av欧美| 精品一区二区三区av网在线观看| 亚洲激情在线av| 国产又黄又爽又无遮挡在线| 色综合亚洲欧美另类图片| 亚洲美女搞黄在线观看 | 色综合婷婷激情| 90打野战视频偷拍视频| 中文字幕高清在线视频| 国产精品爽爽va在线观看网站| 国产探花在线观看一区二区| 免费高清视频大片| 欧美最黄视频在线播放免费| 国产欧美日韩精品亚洲av| 国产色婷婷99| 亚洲熟妇熟女久久| 国产精品野战在线观看| 床上黄色一级片| 精品一区二区三区av网在线观看| 久久99热这里只有精品18| 亚洲欧美精品综合久久99| 午夜福利高清视频| 国产精品嫩草影院av在线观看 | 中文资源天堂在线| 亚洲最大成人av| 国产伦在线观看视频一区| 欧美一区二区亚洲| 成人永久免费在线观看视频| а√天堂www在线а√下载| 757午夜福利合集在线观看| 国产高清有码在线观看视频| 少妇高潮的动态图| 毛片一级片免费看久久久久 | 国产高清三级在线| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 九色成人免费人妻av| 精品久久久久久久人妻蜜臀av| 久久人人爽人人爽人人片va | 国产亚洲av嫩草精品影院| 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 成人av一区二区三区在线看| 国产成人aa在线观看| 久久久久亚洲av毛片大全| 国产单亲对白刺激| av欧美777| 国产在视频线在精品| 波多野结衣高清作品| 美女被艹到高潮喷水动态| 露出奶头的视频| 熟妇人妻久久中文字幕3abv| 亚洲av熟女| 欧美最新免费一区二区三区 | 国内精品久久久久精免费| 女生性感内裤真人,穿戴方法视频| 怎么达到女性高潮| 色综合站精品国产| 韩国av一区二区三区四区| 51午夜福利影视在线观看| 白带黄色成豆腐渣| 俺也久久电影网| 一进一出抽搐gif免费好疼| 国产精品美女特级片免费视频播放器| 午夜福利高清视频| 夜夜夜夜夜久久久久| 亚洲国产精品sss在线观看| 午夜福利成人在线免费观看| 丰满的人妻完整版| 国内精品久久久久精免费| 国产精品爽爽va在线观看网站| 亚洲成av人片在线播放无| 日日干狠狠操夜夜爽| 麻豆成人av在线观看| 亚洲欧美日韩高清专用| 97碰自拍视频| 国产成人欧美在线观看| 久久亚洲精品不卡| av在线老鸭窝| 欧美日韩中文字幕国产精品一区二区三区| 波多野结衣巨乳人妻| 欧美日韩乱码在线| 亚洲经典国产精华液单 | 国产精品久久久久久久电影| 亚洲人成电影免费在线| av中文乱码字幕在线| 午夜福利成人在线免费观看| 丰满人妻一区二区三区视频av| 天堂网av新在线| 我要看日韩黄色一级片| 少妇高潮的动态图| 亚洲五月婷婷丁香| 91字幕亚洲| 欧美黄色淫秽网站| 久久99热这里只有精品18| 日韩大尺度精品在线看网址| 极品教师在线免费播放| 婷婷亚洲欧美| 毛片女人毛片| 日本五十路高清| 国产久久久一区二区三区| 波多野结衣高清作品| 91九色精品人成在线观看| 99久久无色码亚洲精品果冻| 精品一区二区三区人妻视频| 麻豆av噜噜一区二区三区| 亚洲综合色惰| 赤兔流量卡办理| 51午夜福利影视在线观看| 中文亚洲av片在线观看爽| 午夜激情欧美在线| 一区福利在线观看| 免费观看精品视频网站| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 久久午夜亚洲精品久久| 人人妻,人人澡人人爽秒播| 精品久久久久久成人av| 男女床上黄色一级片免费看| netflix在线观看网站| 久久亚洲真实| ponron亚洲| 色av中文字幕| АⅤ资源中文在线天堂| 免费看光身美女| 成年人黄色毛片网站| 午夜影院日韩av| 天堂影院成人在线观看| 90打野战视频偷拍视频| 国语自产精品视频在线第100页| www.www免费av| 神马国产精品三级电影在线观看| 一级av片app| 久久久久久久久久黄片| 亚洲无线观看免费| 久久久久九九精品影院|