• <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)

    久久欧美精品欧美久久欧美| 国产一区二区激情短视频| 久久香蕉国产精品| 国产精品精品国产色婷婷| 99国产精品一区二区蜜桃av| 欧美一级毛片孕妇| 激情在线观看视频在线高清| 国产高清激情床上av| 天美传媒精品一区二区| 午夜福利成人在线免费观看| 亚洲av五月六月丁香网| 午夜福利欧美成人| 观看美女的网站| 国产精品美女特级片免费视频播放器| 怎么达到女性高潮| 老熟妇仑乱视频hdxx| 亚洲av免费在线观看| 亚洲精品久久国产高清桃花| 亚洲人成网站高清观看| 日本黄色视频三级网站网址| 五月玫瑰六月丁香| 成年人黄色毛片网站| 久久久久精品国产欧美久久久| 国产毛片a区久久久久| 亚洲成人中文字幕在线播放| 亚洲五月天丁香| 91麻豆av在线| 成人国产综合亚洲| 99久久精品热视频| 久久性视频一级片| 国产综合懂色| 亚洲精品一卡2卡三卡4卡5卡| 色在线成人网| 香蕉丝袜av| 免费在线观看亚洲国产| 搞女人的毛片| 欧美bdsm另类| 最近最新免费中文字幕在线| 久久国产精品人妻蜜桃| 夜夜看夜夜爽夜夜摸| 国产欧美日韩一区二区三| 99久久无色码亚洲精品果冻| 日本在线视频免费播放| 精品一区二区三区人妻视频| 亚洲精品美女久久久久99蜜臀| 国模一区二区三区四区视频| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 国产成人福利小说| 在线看三级毛片| 看免费av毛片| 欧美日韩综合久久久久久 | 亚洲av二区三区四区| 又紧又爽又黄一区二区| 午夜福利高清视频| 亚洲精品成人久久久久久| 国产成年人精品一区二区| 最近视频中文字幕2019在线8| 亚洲精华国产精华精| 少妇人妻精品综合一区二区 | 男女床上黄色一级片免费看| 嫩草影院精品99| 国模一区二区三区四区视频| 女人被狂操c到高潮| 国产高清videossex| 国产黄片美女视频| 久久久久久久久久黄片| 亚洲国产色片| 国产激情偷乱视频一区二区| 亚洲专区中文字幕在线| 噜噜噜噜噜久久久久久91| 99热6这里只有精品| 韩国av一区二区三区四区| av中文乱码字幕在线| 99国产综合亚洲精品| 国产视频内射| 在线十欧美十亚洲十日本专区| 天堂影院成人在线观看| 一个人看的www免费观看视频| 国产主播在线观看一区二区| 久久九九热精品免费| 脱女人内裤的视频| av中文乱码字幕在线| 国产精品女同一区二区软件 | 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 日韩欧美国产一区二区入口| 69av精品久久久久久| 精品久久久久久成人av| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久久电影 | 国产熟女xx| 亚洲乱码一区二区免费版| 我要搜黄色片| 精品人妻一区二区三区麻豆 | 国产精品久久久久久人妻精品电影| 成熟少妇高潮喷水视频| 国产精品98久久久久久宅男小说| 十八禁网站免费在线| av女优亚洲男人天堂| 一级毛片高清免费大全| 精华霜和精华液先用哪个| 啪啪无遮挡十八禁网站| 日日夜夜操网爽| 伊人久久大香线蕉亚洲五| 国产一区二区三区在线臀色熟女| 人妻久久中文字幕网| 三级男女做爰猛烈吃奶摸视频| 欧美黄色片欧美黄色片| 国产日本99.免费观看| 国产精品亚洲美女久久久| 美女大奶头视频| 69av精品久久久久久| 精品久久久久久久久久久久久| 久久精品91蜜桃| 变态另类丝袜制服| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品av在线| 亚洲精品成人久久久久久| 亚洲五月天丁香| 制服人妻中文乱码| 身体一侧抽搐| 免费观看人在逋| 亚洲av电影不卡..在线观看| 亚洲专区国产一区二区| 午夜福利免费观看在线| 欧美另类亚洲清纯唯美| 亚洲 欧美 日韩 在线 免费| 叶爱在线成人免费视频播放| 亚洲天堂国产精品一区在线| 每晚都被弄得嗷嗷叫到高潮| 日本免费a在线| 99热这里只有精品一区| 国内少妇人妻偷人精品xxx网站| 国产一区二区三区在线臀色熟女| 久久久久精品国产欧美久久久| 国产一区二区三区视频了| 国产一区二区激情短视频| 欧美丝袜亚洲另类 | 在线国产一区二区在线| 搡女人真爽免费视频火全软件 | 黄色日韩在线| 有码 亚洲区| 97超级碰碰碰精品色视频在线观看| 伊人久久精品亚洲午夜| 美女黄网站色视频| 中文字幕精品亚洲无线码一区| 国产精品亚洲美女久久久| 搡女人真爽免费视频火全软件 | 日本成人三级电影网站| 久久久久久大精品| 亚洲第一电影网av| 日韩欧美国产一区二区入口| www.www免费av| 中文字幕久久专区| 欧美在线黄色| 别揉我奶头~嗯~啊~动态视频| 91九色精品人成在线观看| 国产一区二区在线观看日韩 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧洲综合997久久,| 色av中文字幕| 老熟妇仑乱视频hdxx| 午夜亚洲福利在线播放| 免费观看精品视频网站| 一级a爱片免费观看的视频| 国产精品久久视频播放| 色综合婷婷激情| 欧美在线黄色| 男人和女人高潮做爰伦理| 亚洲国产欧美网| eeuss影院久久| 免费观看的影片在线观看| 尤物成人国产欧美一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| www日本在线高清视频| 中文字幕av成人在线电影| 亚洲av免费在线观看| 19禁男女啪啪无遮挡网站| 中亚洲国语对白在线视频| 少妇的丰满在线观看| 久久精品国产综合久久久| 97碰自拍视频| 亚洲精品国产精品久久久不卡| 最近最新免费中文字幕在线| 欧美极品一区二区三区四区| 欧美激情久久久久久爽电影| 成人三级黄色视频| 90打野战视频偷拍视频| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品综合一区在线观看| 91久久精品电影网| 亚洲精品色激情综合| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 亚洲性夜色夜夜综合| 国产精品 欧美亚洲| ponron亚洲| 亚洲七黄色美女视频| 午夜久久久久精精品| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 少妇的逼水好多| 一区福利在线观看| 国产69精品久久久久777片| 麻豆成人av在线观看| 女生性感内裤真人,穿戴方法视频| 成人亚洲精品av一区二区| 好男人电影高清在线观看| 叶爱在线成人免费视频播放| 脱女人内裤的视频| 91麻豆精品激情在线观看国产| 无人区码免费观看不卡| 久久精品91蜜桃| 精品久久久久久久毛片微露脸| 欧美成人一区二区免费高清观看| 在线播放国产精品三级| 亚洲成人免费电影在线观看| 国产黄色小视频在线观看| 精品一区二区三区av网在线观看| 亚洲国产精品999在线| 真人一进一出gif抽搐免费| 国产精品99久久久久久久久| 九色成人免费人妻av| 成熟少妇高潮喷水视频| 国产男靠女视频免费网站| 免费看十八禁软件| 欧美激情久久久久久爽电影| 俄罗斯特黄特色一大片| 脱女人内裤的视频| 婷婷六月久久综合丁香| 99久国产av精品| 精品一区二区三区视频在线观看免费| 18禁美女被吸乳视频| www国产在线视频色| 欧美日韩精品网址| 深爱激情五月婷婷| 天堂网av新在线| 午夜福利欧美成人| 成人鲁丝片一二三区免费| 身体一侧抽搐| 桃红色精品国产亚洲av| 国产精品国产高清国产av| 19禁男女啪啪无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 97人妻精品一区二区三区麻豆| 免费观看的影片在线观看| 亚洲欧美日韩东京热| 国产一区二区激情短视频| www.熟女人妻精品国产| 精华霜和精华液先用哪个| 18禁黄网站禁片午夜丰满| av福利片在线观看| 亚洲七黄色美女视频| 看片在线看免费视频| 国产免费av片在线观看野外av| 久久久久性生活片| 国产爱豆传媒在线观看| 俺也久久电影网| 少妇熟女aⅴ在线视频| 国产成人福利小说| 亚洲欧美日韩高清在线视频| 国产精品一及| 亚洲国产中文字幕在线视频| www日本黄色视频网| 51国产日韩欧美| 日本黄色视频三级网站网址| 久久久久国内视频| 国语自产精品视频在线第100页| 中文字幕熟女人妻在线| av专区在线播放| 99久久精品热视频| 久久久国产成人精品二区| 成年女人看的毛片在线观看| 十八禁人妻一区二区| 天天一区二区日本电影三级| 99国产综合亚洲精品| 日日夜夜操网爽| 国产精品香港三级国产av潘金莲| 国产亚洲精品久久久com| 久久亚洲精品不卡| 男女床上黄色一级片免费看| 可以在线观看的亚洲视频| 91在线观看av| 亚洲成人中文字幕在线播放| 俺也久久电影网| 一区二区三区国产精品乱码| 搡老岳熟女国产| 久久久久九九精品影院| 成人一区二区视频在线观看| 婷婷精品国产亚洲av在线| 久久久久久久精品吃奶| 尤物成人国产欧美一区二区三区| 91久久精品电影网| 在线观看美女被高潮喷水网站 | 久久性视频一级片| 国产激情偷乱视频一区二区| 欧美bdsm另类| xxxwww97欧美| 亚洲精品亚洲一区二区| 在线观看免费午夜福利视频| 国产高清视频在线播放一区| 亚洲av免费高清在线观看| h日本视频在线播放| 深爱激情五月婷婷| 欧美日韩国产亚洲二区| 欧美性感艳星| 国产伦在线观看视频一区| 国产亚洲欧美在线一区二区| 午夜视频国产福利| 国产私拍福利视频在线观看| 一级黄片播放器| 久久久久国产精品人妻aⅴ院| 搡女人真爽免费视频火全软件 | 久久午夜亚洲精品久久| 欧美日本视频| 最新中文字幕久久久久| 精品欧美国产一区二区三| 久久久久免费精品人妻一区二区| 一a级毛片在线观看| 禁无遮挡网站| 在线观看舔阴道视频| 一区福利在线观看| 亚洲欧美日韩高清专用| 午夜日韩欧美国产| 久久精品夜夜夜夜夜久久蜜豆| 久久九九热精品免费| 搡老熟女国产l中国老女人| 国产精品一及| 中出人妻视频一区二区| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| ponron亚洲| 1024手机看黄色片| 一区二区三区高清视频在线| 亚洲国产色片| 大型黄色视频在线免费观看| 精品无人区乱码1区二区| 老熟妇仑乱视频hdxx| 精品欧美国产一区二区三| 久久精品91蜜桃| 久久久久九九精品影院| 国产午夜福利久久久久久| 国产精品爽爽va在线观看网站| 国产高清videossex| 2021天堂中文幕一二区在线观| 欧美日韩乱码在线| 国产精品一区二区免费欧美| 人人妻人人澡欧美一区二区| 18美女黄网站色大片免费观看| 色综合婷婷激情| 色综合亚洲欧美另类图片| 女警被强在线播放| 少妇高潮的动态图| 国产成人福利小说| 99热这里只有精品一区| 欧美一区二区精品小视频在线| 亚洲av熟女| 免费观看精品视频网站| 久久精品国产自在天天线| 麻豆成人午夜福利视频| 成人国产综合亚洲| 搡老妇女老女人老熟妇| 成人国产综合亚洲| 久久精品91无色码中文字幕| 精品不卡国产一区二区三区| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 搡老岳熟女国产| 最新在线观看一区二区三区| 美女大奶头视频| 老司机午夜福利在线观看视频| 亚洲乱码一区二区免费版| 1024手机看黄色片| 国产黄片美女视频| 亚洲内射少妇av| 久久久久性生活片| 亚洲电影在线观看av| 亚洲欧美一区二区三区黑人| 夜夜躁狠狠躁天天躁| 少妇丰满av| 少妇人妻精品综合一区二区 | 一级黄片播放器| 亚洲成人免费电影在线观看| 国产黄色小视频在线观看| 脱女人内裤的视频| 搡女人真爽免费视频火全软件 | 一区二区三区高清视频在线| 国产成人av教育| 国产精品影院久久| 成人鲁丝片一二三区免费| 一级a爱片免费观看的视频| 精品日产1卡2卡| 久久亚洲精品不卡| 丰满人妻熟妇乱又伦精品不卡| 欧美bdsm另类| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美一区二区三区在线观看| 亚洲精品456在线播放app | 亚洲国产高清在线一区二区三| 精品久久久久久久毛片微露脸| 国产色爽女视频免费观看| 在线视频色国产色| 成人性生交大片免费视频hd| 久久香蕉精品热| 超碰av人人做人人爽久久 | 久久欧美精品欧美久久欧美| 亚洲av免费高清在线观看| 白带黄色成豆腐渣| 亚洲,欧美精品.| 性欧美人与动物交配| 99久久久亚洲精品蜜臀av| 18禁美女被吸乳视频| 看片在线看免费视频| 亚洲,欧美精品.| 亚洲人成网站在线播| 最好的美女福利视频网| 真人一进一出gif抽搐免费| 一个人看的www免费观看视频| 制服丝袜大香蕉在线| 成人精品一区二区免费| 岛国视频午夜一区免费看| 亚洲久久久久久中文字幕| 丁香欧美五月| 一级黄片播放器| 一区二区三区激情视频| 国产伦人伦偷精品视频| 可以在线观看毛片的网站| 午夜免费观看网址| 欧美一级毛片孕妇| 人人妻,人人澡人人爽秒播| 亚洲最大成人手机在线| 三级男女做爰猛烈吃奶摸视频| 91在线精品国自产拍蜜月 | 老汉色av国产亚洲站长工具| 日韩欧美在线二视频| 亚洲成人中文字幕在线播放| 很黄的视频免费| 日韩欧美 国产精品| 亚洲av熟女| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 黄色女人牲交| 国产高清videossex| 高清日韩中文字幕在线| 国产色爽女视频免费观看| 日本精品一区二区三区蜜桃| 午夜激情欧美在线| 啦啦啦观看免费观看视频高清| 国产91精品成人一区二区三区| 蜜桃久久精品国产亚洲av| 伊人久久大香线蕉亚洲五| 久久久色成人| 亚洲成人久久性| 最近在线观看免费完整版| 中文亚洲av片在线观看爽| 亚洲国产精品成人综合色| 久久精品国产综合久久久| 麻豆国产97在线/欧美| 亚洲欧美日韩高清在线视频| 久久这里只有精品中国| 欧美3d第一页| 神马国产精品三级电影在线观看| 成年人黄色毛片网站| 熟女电影av网| 国产精品久久视频播放| 老司机在亚洲福利影院| 免费大片18禁| 看片在线看免费视频| 天堂动漫精品| 日本三级黄在线观看| 国产真实伦视频高清在线观看 | 亚洲片人在线观看| 99久久精品国产亚洲精品| 国产一区二区三区视频了| 一个人观看的视频www高清免费观看| 亚洲成av人片在线播放无| 国产黄a三级三级三级人| 国产色爽女视频免费观看| 精品国产三级普通话版| 欧美中文日本在线观看视频| 国产精品久久久久久亚洲av鲁大| 在线观看免费午夜福利视频| 色在线成人网| 最近最新中文字幕大全电影3| 色综合站精品国产| 欧美激情久久久久久爽电影| 1024手机看黄色片| 成人三级黄色视频| 黑人欧美特级aaaaaa片| 国产三级黄色录像| 欧美激情久久久久久爽电影| 久9热在线精品视频| 精品久久久久久,| 天堂影院成人在线观看| 国产精品爽爽va在线观看网站| 亚洲av五月六月丁香网| 精品熟女少妇八av免费久了| 国产真人三级小视频在线观看| 真实男女啪啪啪动态图| 叶爱在线成人免费视频播放| 熟女电影av网| 日本黄色片子视频| 丰满人妻一区二区三区视频av | 变态另类丝袜制服| 久久精品国产亚洲av涩爱 | 啦啦啦韩国在线观看视频| 观看美女的网站| 国产精品久久久久久久电影 | 国产伦精品一区二区三区视频9 | 欧美黄色片欧美黄色片| 香蕉久久夜色| 日本一二三区视频观看| 国产av在哪里看| 国产高清videossex| 国产精品亚洲av一区麻豆| 久久久成人免费电影| 国产激情欧美一区二区| 观看美女的网站| 成年女人看的毛片在线观看| 久久久国产成人精品二区| 99热只有精品国产| 欧美+日韩+精品| 久久性视频一级片| 人人妻,人人澡人人爽秒播| 熟女少妇亚洲综合色aaa.| 免费高清视频大片| 少妇熟女aⅴ在线视频| 真人一进一出gif抽搐免费| 国产成人系列免费观看| 亚洲欧美日韩高清在线视频| 日日干狠狠操夜夜爽| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 日韩欧美三级三区| 内射极品少妇av片p| 男人的好看免费观看在线视频| 亚洲精品在线观看二区| 在线播放国产精品三级| 黄色视频,在线免费观看| 国产亚洲精品久久久久久毛片| 很黄的视频免费| 91麻豆精品激情在线观看国产| 亚洲av中文字字幕乱码综合| 久久精品91无色码中文字幕| 日本五十路高清| 午夜久久久久精精品| 偷拍熟女少妇极品色| 亚洲精品在线美女| 成年女人毛片免费观看观看9| 美女高潮喷水抽搐中文字幕| 欧美绝顶高潮抽搐喷水| 青草久久国产| 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜添小说| 99国产精品一区二区蜜桃av| 国产乱人视频| 日韩欧美国产在线观看| 超碰av人人做人人爽久久 | 午夜福利在线观看吧| 国产三级中文精品| 精品久久久久久久久久久久久| 极品教师在线免费播放| 国产一级毛片七仙女欲春2| 婷婷丁香在线五月| 天堂动漫精品| 白带黄色成豆腐渣| 国内精品久久久久久久电影| 日韩欧美三级三区| 国内精品一区二区在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲激情在线av| 久久精品国产综合久久久| 久久九九热精品免费| 午夜福利欧美成人| 看黄色毛片网站| 午夜老司机福利剧场| 久久人人精品亚洲av| 中文字幕久久专区| 日本 欧美在线| 最好的美女福利视频网| 日韩大尺度精品在线看网址| 母亲3免费完整高清在线观看| 三级国产精品欧美在线观看| 99久久精品国产亚洲精品| 国产伦精品一区二区三区视频9 | 亚洲乱码一区二区免费版| 日本在线视频免费播放| 国产黄a三级三级三级人| 搞女人的毛片| 啪啪无遮挡十八禁网站| av在线蜜桃| 欧美黄色淫秽网站| 国产91精品成人一区二区三区| 国产av麻豆久久久久久久| 观看美女的网站| 成人无遮挡网站| 国产成人av教育| 成人三级黄色视频| 欧美极品一区二区三区四区| 久久久成人免费电影| 国产精品永久免费网站| 亚洲国产欧洲综合997久久,| 免费av观看视频| 国产高清有码在线观看视频| 免费无遮挡裸体视频| 国产精品 国内视频| 午夜福利在线观看免费完整高清在 | 亚洲国产欧美网| 免费搜索国产男女视频| 国产高潮美女av| 国产成人欧美在线观看| 香蕉av资源在线| 国产亚洲精品综合一区在线观看|