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

    The use of widely targeted metabolite profiling to reveal the senescence changes in postharvest ‘Red Globe’ (Vitis vinifera)grape berries

    2022-03-16 03:05:22XUTengfeiYANGXinZHANGMengGUOShuihuanFUWenjingZHOUBijiangLlUYujiaMAHaijunFANGYulinYANGGangMENGJiangfei
    Journal of Integrative Agriculture 2022年4期

    XU Teng-fei,YANG Xin,ZHANG Meng,GUO Shui-huan,FU Wen-jing,ZHOU Bi-jiang,LlU Yu-jia,MA Hai-jun,FANG Yu-lin,YANG Gang,MENG Jiang-fei

    1 College of Horticulture,Northwest A&F University,Yangling 712100,P.R.China

    2 College of Enology,Northwest A&F University,Yangling 712100,P.R.China

    3 College of Biological Science and Engineering,North Minzu University,Yinchuan 750021,P.R.China

    4 Shenzhen Camartina Wine Co.,Ltd.,Shenzhen 518031,P.R.China

    Abstract Changes in the metabolites of table grapes (Vitis vinifera) reportedly occur during postharvest senescence.The aim of this study was to determine the metabolomic differences in postharvest table grapes (‘Red Globe’) after being subjected to different senescence periods.To this end,we used widely targeted metabolomics based on ultra-performance liquid chromatography and tandem mass spectrometry.A total of 135 differential metabolites were identified.During postharvest senescence,the levels of most differential flavonoids (e.g.,pelargonidin 3-O-glucoside,quercetin-3-O-glucoside,and cyanidin 3-O-glucoside) and L-aspartic acid decreased,while the levels of phenolic acids (e.g.,trans-4-hydroxycinnamic acid methyl ester) and pantothenol increased.During early and late senescence,the levels of most differential lipids,especially LysoPC,as well as those of nucleotides and their derivatives,such as uridine,decreased and increased,respectively.Collectively,the findings of this study provide fundamental insights into the reasonable control of table grape fruit postharvest senescence and lay a solid foundation for further research.

    Keywords:table grape, berry,postharvest senescence,widely targeted metabolites

    1.lntroduction

    Grapevines are fruit trees which are cultivated worldwide.The global grape plantation area was forecasted to reach 7.45 Mha in 2018 (OIV 2019).Currently,grapes are mainly consumed fresh or processed into raisins,wine,and juice.Most table grapes are not consumed immediately after harvest but are stored and made available in the market at the appropriate time.However,once grape clusters are harvested and separated from the plants,they enter senescence.Indeed,grapes are subject to postharvest senescence,exhibiting characteristics such as severe moisture loss,berry softening,off-flavor development,and decay that are caused mainly byBotrytiscinerea,which reduces the commodity value and suitability of grapes for consumption (Niet al.2016),thereby leading to reduced postharvest and shelf lives.

    Fruit ripening and senescence are inevitable and irreversible processes in the plant life cycle (Dinget al.2015).Fruits are generally divided into climacteric and nonclimacteric types based on the patterns of respiration and ethylene production.The ethylene content and respiration of climacteric fruit increase sharply during ripening as CO2evolution.In contrast,non-climacteric fruit do not undergo these changes and normally emit ethylene in markedly decreased levels (Bottonet al.2019).Furthermore,a series of metabolic changes related to carbohydrates,lipids,organic acids,and secondary metabolites takes place in fruit after harvesting (Pedreschiet al.2009).During fruit senescence,the sugar and organic acid contents decrease gradually with the extension of storage time because they are consumed during fruit respiration.The sucrose,fructose,and glucose contents of postharvest apples decline owing to the action of sucrose synthase and invertase (Zhuet al.2013).Sunet al.(2013) showed that the levels of succinic acid,γ-aminobutyric acid,and glutamine increase,but the level of 2-oxoglutaric acid decreases,in postharvestCitrusfruit.Additionally,the organic acid and hydrogen peroxide contents,and the superoxide dismutase and peroxidase activities,are obviously associated with senescencerelated physiological processes,indicating that organic acids could be used to measure the senescence level of postharvestCitrusfruit.During postharvest senescence,the phenolic and lignin contents of fruit decrease and increase,respectively (Caiet al.2006;Zenoniet al.2020).However,it is unclear whether any other metabolites are closely associated with postharvest senescence.

    Metabolomics refers to the analysis of all metabolites in an organism and allows for the simultaneous measurement of all metabolites in a given biological system (Dixon and Strack 2003).In recent years,metabolomic approaches have been widely used to explore aging-related mechanisms in fruit,as reported for melons,oranges,pears,cherries,bananas,and blackberries (Bernillonet al.2013;Dinget al.2015;Yuanet al.2017;Karagianniset al.2018;Xuet al.2018;Kimet al.2019;Sainiet al.2019).However,senescence differs depending on the fruit tree species because of chemical composition differences.Thus,the underlying mechanisms of senescence are unique among fruit types.To date,there have been no studies on the dynamic changes of endogenous metabolites within the integrated grape berry during postharvest senescence.

    Recently,widely targeted metabolomics technology based on ultra-performance liquid chromatography and tandem mass spectrometry (UPLC-MS/MS) has been used to quantify metabolites in a high-throughput manner.This technology combines the advantages of accurate measurements in terms of targeted metabolomics with the wide metabolite coverage conferred by the non-targeted metabolomics (Chenet al.2013).The aim of the present study was to characterize the metabolic changes during the natural senescence of postharvest grape berries by using a widely targeted UPLC-electrospray ionization (ESI)-MS/MS-based metabolomics approach with the application of multivariate statistical tools.The metabolomics approach used in this study will help in establishing a direct correlation between the senescence characteristics and metabolite expression during the postharvest senescence of table grapes,thereby contributing to the development of new strategies to slow the rate of grape senescence.

    2.Materials and methods

    2.1.Plant materials and fruit sampling process

    Grape clusters of ‘Red Globe’ (Vitisvinifera) table grapes were harvested at commercial maturity from an organic vineyard in Weinan City (Shaanxi Province,China).Grapes of uniform size,shape,and color,and free from disease and mechanical damage,were transported to the laboratory.Each grape cluster (approximately 40 berries) was weighed and labeled,and then stored at 20°C with a relative humidity of 80%.Six clusters were collected at 0,4,8,12,16,20,and 24 d postharvest and weighed every 4 d to calculate the water loss rate.Then,the berries of the six clusters were collected and mixed.Half of them were used immediately for the analysis of berry firmness,and the other half were frozen in liquid nitrogen and immediately stored at -80°C for further metabolome analysis.All experiments were performed in biological triplicates.

    2.2.Detection of grape berry firmness

    A total of 15 berries were randomly selected and cut along the latitude of the berry forming a kerf of approximately 1 cm2.Then,the probe of a fruit firmness meter (GY-4;Zhejiang Top Instrument Co.,Ltd.,Hangzhou,China) was vertically inserted into the berry kerf to a depth of 10 mm,while avoiding the seeds,to detect the hardness.The measurement results were recorded in N.

    2.3.Sample preparation and extraction for metabolome analysis

    A total of 30 grape berries randomly selected at 0,12,and 24 d postharvest (recorded as RG0,RG12,and RG24,respectively) were deseeded and ground in liquid nitrogen using a blender for 30 s.A 100-mg sample of the freezedried powder was incubated in a 1.5-mL centrifuge tube with 0.6 mL of 70% aqueous methanol at 4°C for 8 h.The mixture was centrifuged at 10 000×g for 10 min at 4°C and the supernatant liquor was absorbed and filtered prior to UPLC-MS/MS analysis.

    2.4.UPLC-ESl-MS/MS-based widely target metabolomics

    UPLC-ESI-MS/MS measurements were performed on the Shim-pack UFLC SHIMADZU CBM30A system(Shimadzu,Kyoto,Japan) coupled to the Applied Biosystems 4500 QTRAP (Applied Biosystems,Foster City,CA,USA).Chromatographic separation was achieved on a Waters ACQUITY UPLC HSS T3 C18 (1.8 μm,2.1 mm×100 mm) (Waters Corporation,Milford,MA,USA).Eluent A was ultrapure water with 0.04% acetic acid,and eluent B was acetonitrile with 0.04% acetic acid.Sample measurements were performed with a gradient program that employed starting conditions of 95% A and 5% B.Within 10 min,a linear gradient to 5% A and 95% B was run,and that composition was maintained for 1 min.Subsequently,the composition was adjusted to 95% A and 5% B within 0.10 min and maintained for 2.9 min.In addition,the column temperature was 40°C with a 4-μL injection volume.The metabolites separated by liquid chromatography were entered directly into the mass spectrometer for detection.

    The MS conditions were as follows:turbo spray;550°C for source temperature;5 500 V (positive ion mode)/-4 500 V(negative ion mode);50,60,and 30 psi for ion source gas I,gas II,and curtain gas,respectively;multiple reaction monitoring (MRM) mode;and 5 psi nitrogen for collision gas.A specific set of MRM transitions was monitored for each period according to the metabolites eluted within that period.

    2.5.Screening and function annotation of differential metabolites

    Orthogonal partial least squares-discriminant analysis(OPLS-DA) was used to initially screen for metabolites with differences during postharvest senescence.Differential metabolites could be further screened by combining theP-values or fold changes in the univariate analysis.Differentially expressed metabolites were filtered by variable importance in projection (VIP) ≥1 and absolute log2FC(fold change) ≥1 or ≤-1.The data were log transformed(log2) and their means were centered before OPLS-DA.To avoid overfitting,a permutation test (200 permutations)was performed.

    The identified metabolites were annotated using the KEGG Compound Database,and the annotated metabolites were mapped to the KEGG Pathway Database.Pathways with significantly expressed metabolites were mapped and subjected to metabolite set enrichment analysis (MSEA),and their significance was analyzed through hypergeometric testP-values.

    2.6.Statistical analysis

    Principal component analysis (PCA),OPLS-DA,and hierarchical cluster analysis (HCA) were performed using the R package (www.r-project.org).PCA was used to recombine the original variables into new,mutually independent variables through orthogonal transformations.HCA was performed on the accumulation patterns of metabolites between different samples.

    3.Results

    3.1.Water loss rate and berry firmness change

    The grape berry firmness and water loss rate changes during‘Red Globe’ postharvest senescence are shown in Fig.1.With the extension of senescence time,the water loss rate of ‘Red Globe’ grape clusters gradually increased,whereas grape berry firmness increased during the first 4 d,and then decreased during the final 20 d.The pattern between water loss and firmness fits the features.

    Fig.1 Changes in the water loss rate (A) and berry firmness (B) during ‘Red Globe’ grapes postharvest senescence.Error bars represent SD (n=3).Different letters indicate significant differences according to Duncan’s multiple range test (P<0.05).

    3.2.Data quality assessment of the UPLC-ESl-MS/MS results

    During analysis,quality control (QC) was performed by inserting the mixture of the sample extracts after each group of 10 test samples,to monitor the reproducibility of the unit operations during the analysis process.The total ion current (TIC) is the spectrum obtained by continuously summing the intensities of all ions in the mass spectrum at each time point.By using TIC,the accuracy and reproducibility of metabolite detection were confirmed according to Appendices A and B,satisfactory accuracy and reproducibility were achieved.In the MRM mode,the ion flux spectrum of each extracted substance can be determined by the multi-substance extracted ion chromatogram (XIC).The various metabolites were detected and shown by the variant colors of the mass spectral peaks.The peak area of each sample was calculated and represented the relative abundance of the corresponding substance (Appendix C).The MultiQuant Software was used to perform an integration and calibrationof the detected substance-specific peaks.

    3.3.Qualitative and quantitative metabolites

    We performed qualitative and quantitative analyses of the metabolites in the postharvest grape samples based on the KEGG database,MetWare database(MWDB),and MRM.A total of 424 metabolites were detected:57 amino acids and derivatives,31 nucleotides and derivatives,29 organic acids,17 free fatty acids,28 saccharides and alcohols,7 dihydroflavones,5 dihydroflavonols,3 phenolamines,38 phenolic acids,7 glycerol esters,12 anthocyanins,32 flavonoids,52 flavonols,11 flavanols,3 phosphatidylcholines,2 lignans,13 lysophosphatidylcholines,6 lysophosphatidylethanolamines,1 sphingolipid,2 tannins,9 triterpenes,8 alkaloids,9 vitamins,2 coumarins,2 isoflavones,1 plumerane,10 proanthocyanidins,6 stilbenes,and 21 other metabolites (Appendix D).

    3.4.Principal component analysis

    The PCA showed that the variability in each grape berry sample was small.The samples obtained during the same sampling session had similar metabolic characteristics.Therefore,our metabolomics results were stable and reproducible.Additionally,the separation trends among all treatments were obvious,suggesting significant metabolic differences during postharvest senescence (Fig.2-A).The three senescence periods were characterized by distinct metabolisms,and the effect of senescence time on the metabolism of postharvest ‘Red Globe’ berries was obvious(Fig.2-B-D).

    Fig.2 Principal component analysis (PCA) of all identified metabolites during ‘Red Globe’ grapes postharvest senescence.A,overall score scatter plot;B to D,PCA score maps of the RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.

    3.5.OPLS-DA

    To further maximize the distinctions between grape samples,we performed OPLS-DA,which is a more conclusive way to find differential metabolites.This method allows for theXmatrix information to decompose intoYcorrelation and irrelevance (OPLS-DA).In this method,the differential variables are extracted by eliminating the irrelevant differences.R2Yrepresents the interpretation rate of the model on theYmatrix,andQ2is the prediction ability of the model on theYmatrix.When the values ofR2YandQ2are closer to 1,the model is more stable and reliable.In general,the model is effective whenQ2>0.5.In this study,R2YandQ2were greater than 0.99 and 0.50,respectively,in each pair group (Fig.3-A-C),thereby demonstrating that this model was effective.The OPLS-DA model was confirmed through 200 alignment experiments.The horizontal line refers to theR2andQ2of the original model,whereas the red and blue dots represent theR2′ andQ2′ after replacement,respectively.These results showed that the current model was meaningful.Additionally,these differential metabolites could be screened on the basis of the VIP value in the subsequent analysis (Fig.3-D-F).

    Fig.3 Orthogonal partial least squares-discriminant analysis (OPLS-DA) scores and permutation verification of all identified metabolites during ‘Red Globe’ grapes postharvest senescence.A-C,scores of the OPLS-DA model with RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.D-F,OPLS-DA permutation analysis model verification charts of RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.

    3.6.HCA and volcano plot

    To evaluate differences in the characteristics of the ‘Red Globe’ berries postharvest senescence that led to the metabolite changes,HCA was used to analyze intratreatment homogeneity and inter-treatment variability.The extension of senescence increased the metabolite differences (Fig.4-A).The points in the volcano map represent all of the detected metabolites from two different senescence periods of the grape berries.Metabolites with a VIP of ≥1 and log2FC of ≥1 or ≤-1 were selected.Based on this method,no differences were detected for most metabolites.The metabolites that met the aforementioned conditions showed significant differences.Among these differential metabolites,the number of metabolites with decreasing levels was greater than the number with increasing levels for RG0vs.RG12 and RG0vs.RG24,while the number of metabolites with increasing levels was greater than the number with decreasing levels for RG12vs.RG24 (Fig.4-B-D).

    Fig.4 Hierarchical cluster analysis (HCA),volcano plot,and Venn diagram of differential metabolites during ‘Red Globe’ grapes postharvest senescence.A,hierarchical cluster of the differential metabolites among 0 (RG0),12 (RG12) and 24 (RG24) d postharvest.B-D,volcano plots of RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.E,Venn diagram of the differential metabolites.

    3.7.Statistical analysis of differential metabolites

    Among the 424 metabolites,135 were identified as differential metabolites during the postharvest senescence of ‘Red Globe’ grapes (Appendix E).These were separated into nine categories:amino acids and their derivatives,nucleotides and their derivatives,lipids,saccharides and alcohols,flavonoids,vitamins,phenolic acids,organic acids,and others (Table 1).A total of 50 differential metabolites were verified in the RG12 samples compared with those in the RG0 samples;the levels of 39 of them decreased,while those of the other 11 increased.A total of 82 differential metabolites were verified in the RG24 samples (relative to RG0);and the levels of 66 decreased while those of 16 increased.In addition,68 differential metabolites were confirmed in the RG24 samples (relative to RG12);and the levels of 23 decreased,while those of 45 increased.It is worth noting that,among all the differential metabolite classes,flavonoids were significantly more abundant during the entire senescence process,while the levels of most other metabolites decreased (Table 1;Appendix E).Furthermore,eight metabolites were common in all three postharvest periods examined,namely,xanthine,trans-4-hydroxycinnamic acid methyl ester,methylp-coumarate,O-phospho-L-serine,pantothenol,(R)-pantetheine,glutathione in its reduced form,and procyanidin A2(Fig.4-E).

    Table 1 Numbers of differential metabolites in ‘Red Globe’ grapes that were tested at 0,12,and 24 d postharvest,denoted as RG0,RG12,and RG24,respectively

    The metabolites that decreased the most in RG0 wereO-phospho-L-serine,allopurinol,glutathione in its reduced form,peonidin 3,5-O-diglucoside chloride,and 12,13-epoxy-9Z-octadecenoic acid (EODE) in sequence (log2FC<-2.5);while the metabolite that increased the most was methylp-coumarate (log2FC reached 9.34) (Fig.5-A and B).The log2FC values of 10 decreasing metabolites from RG24 were below -4.00 compared with those of RG0,including kaempferol-3-arabinopyranoside and gallate catechin gallate which reached -13.11 and -10.86,respectively(Fig.5-C and D).Between RG12 and RG24,the metabolites with increased levels exhibited relatively small fold changes(maximum log2FC only 2.66),although the number of metabolites with increased content was greater than the number with decreased content;and the metabolite that decreased the most was gallate catechin gallate (log2FC reached -12.26) (Fig.5-E and F).

    Fig.5 Fold change bar diagrams (A,C and E) and variable importance in projection (VIP) score plots (B,D and F) of the top ranked differential metabolites.A and B,RG0 vs.RG12.C and D,RG0 vs.RG24.E and F,RG12 vs.RG24.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.

    To further analyze the changes of these differential metabolites during postharvest senescence,we standardized and centralized the contents of all differential metabolites and then performed aK-means clustering analysis,which allowed us to divide the metabolites into six subclasses(Fig.6-A-F).Among them,subclass 5 contained the highest number of metabolites (42 differential metabolites),accounting for approximately one-third of the total.The metabolite content of this subclass decreased gradually as senescence increased (Fig.6-E).Next in order was subclass 4,which contained 28 metabolites.The content of this metabolite subclass mainly decreased from 0 to 12 d postharvest,while there were no significant changes over the next 12 d (Fig.6-E).It is also worth noting that subclass 6,which included pantothenol and 13 other metabolites,exhibited increased metabolite content during the senescence period (Fig.6-F).

    Fig.6 K-means clustering analysis of differential metabolites in the six subclasses.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.

    3.8.Functional annotation and enrichment analysis of the differential metabolites

    To screen and analyze the differential metabolic pathway of‘Red Globe’ grape berries during the fruit senescence periods,the metabolites that differed between RG0 and RG12,RG0 and RG24,and RG12 and RG24 were subjected to KEGG enrichment analysis.Three categories could be assigned to all the metabolites:genetic information processing,environmental information processing,and metabolism.According to the functional annotation results (Fig.7-A-C),the differential metabolites during the senescence of postharvest ‘Red Globe’ grapes were mainly involved in metabolism,especially metabolic pathways and the biosynthesis of secondary metabolites.Some metabolites were also involved in genetic information processing(aminoacyl-tRNA biosynthesis) and environmental information processing (ABC transporters),but both of these types of metabolites accounted for less than 20% of the total.Further enrichment analysis showed some differences in the metabolic pathways involved in the three different senescence phases.For the first 12 d,the metabolites were mainly enriched in purine metabolism;phenylalanine metabolism;pantothenate and CoA biosynthesis;and cysteine,methionine,and caffeine metabolism.Additionally,these metabolites were largely enriched in many other metabolic pathways,including flavonoid metabolism,biosynthesis of secondary metabolites,glycerophospholipid metabolism,and glycosylphosphatidylinositol (GPI)-anchor biosynthesis in the next 12 d (Fig.8-A-C).

    Fig.7 Functional annotation of differential metabolites during the postharvest senescence of ‘Red Globe’ grapes.A-C,RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.The numbers on column indicate the amount and proportion of the differential metabolites,respectively.

    Fig.8 Functional enrichment analysis of differential metabolites during the postharvest senescence of ‘Red Globe’ grapes.A-C,RG0 vs.RG12,RG0 vs.RG24,and RG12 vs.RG24,respectively.RG0,RG12,and RG24 indicate 0,12 and 24 d postharvest,respectively.

    4.Discussion

    Globally,grapes are among the most important fruit species,owing to their economic and nutritional values.Generally,table grape berries are not consumed immediately after harvest.Instead,there is usually a time lapse between harvest and consumption.The loss of firmness or softening of fruits is a very important quality parameter.Berry firmness is determined by the thickness of the outer mesocarp with the opaque colored flesh that contains turgid cells (Du Plessis 2008).The softening is due to cell wall degradation induced by several enzymes that are synergistically activated.Interestingly,grape berry firmness presents a slight rise in the initial stages of postharvest senescence,then drops steeply (Fig.1).The slight rise in berry firmness may be due to the structural changes of turgid cell walls caused by pectate lyase,etc.,which makes their arrangement higher in density and more compact within a certain time-frame at the beginning of harvesting.After that,the hydrolyzation of cell wall polysaccharides finally leads to the loosening of cell walls and loss of cell cohesion that are associated with fruit softening.

    During postharvest senescence,grape berries are not only subjected to water loss and physiological softening,but they also undergo a series of chemical metabolism changes,such as changes in sugars,organic acids,proteins,lipids,and phenolic and aroma compound.In the present study,135 differential metabolites were identified from a total of 424 metabolites quantified from ‘Red Globe’ grape berries at 0,12,and 24 d postharvest.These differential metabolites included polyphenolics,lipids,amino acids and their derivatives,and nucleotides and their derivatives.

    Polyphenolics constituted the class that included the most differential metabolites identified during the postharvest senescence of ‘Red Globe’ grapes.In grapes,these are generally classified as phenolic acids(e.g.,hydroxybenzoic acid and hydroxycinnamic acid),flavonoids (e.g.,flavanols,flavonols,and anthocyanins),and stilbenes (Menget al.2012).Polyphenols have high antioxidant activity and are associated with a reduced risk of a number of cardiovascular diseases and cancers (Oueslatiet al.2012).Therefore,they are used in functional foods,cosmetics,and medicines (Maisuthisakulet al.2007).In this study,the levels of most flavonoids decreased (Table 1),including those of anthocyanins (cyanidin 3-O-galactoside,peonidin 3-O-glucoside,cyanidin 3-O-rutinoside,and peonidin 3,5-O-diglucoside),flavonols (quercetin-3-Oα-L-rhamnopyranoside,quercetin 3-O-glucoside,and kaempferol 7-O-rhamnoside),and flavanols (gallate catechin gallate and gallocatechin 3-O-gallate).These compounds mainly exist in the grape skins and seeds and greatly affect grape color (Ferrandino and Guidoni 2010).Flavanols are the basic building blocks of proanthocyanidins,and their monomers and oligomers contribute bitterness to the astringency of grape berries and red wine (González-Manzanoet al.2004).Flavonols seem to contribute to the bitterness and color (Obreque-Slieret al.2010).Therefore,the decrease in the anthocyanin and flavonol contents may indicate a decrease in the intensity of grape skin color;however,flavanols may form proanthocyanidins after polymerization (the proanthocyanidin A1 and A2 contents increased during senescence),which could affect grape bitterness and the astringent taste during grape postharvest senescence.Furthermore,the nutritional value of grapes may gradually weaken owing to the decreases in the levels of these flavonoids.

    It is worth noting that,unlike the flavonoids,the number of phenolic acid differential metabolites (including methylp-coumarate,trans-4-hydroxycinnamic acid methyl ester,andp-coumaric acid-O-glycoside with increased levels was greater than the number of metabolites with decreased levels.Flavonoids participate in the construction of compounds with a C6-C3-C6carbon framework in their basic structure.These compounds are derivatives of chromanes or chromenes,such as flavans,flavones,flavonols,and anthocyanidins(Yonekura-Sakakibaraet al.2019),whereas phenolic acids are types of aromatic acid compounds containing a phenolic ring and an organic carboxylic acid (Neilson and Ferruzzi 2012).Therefore,during the postharvest senescence of grapes,these flavonoids may be degraded to phenylacetic and phenylpropionic acids (Hollman 2009;Murotaet al.2018) under the catalysis of related enzymes from the grape or some microbe,which leads to a decrease in flavonoid and an increase in phenolic acid contents.

    Cellular membrane systems play an important role in cellular and subcellular compartmentalization.In the process of enzymatic browning,cellular membrane systems can prevent contact between substrates and result in enzymatic alterations in the properties of the cell membrane composition during fruit ripening and postharvest senescence (Lurie and Ben-Arie 1983).During ripening and senescence,the desaturation levels of the fatty acyl chains of membrane phospholipids increase,which results in their peroxidation.Zhanget al.(2010) showed that the expression of a few desaturase isomers is closely related to the continuous flux of linoleate and linoleate substrates of the lipoxygenase pathway duringAmygdaluspersicafruit ripening.The initial excision of fatty acyl chains can result in the accumulation of lysophospholipids in dates,which display an early ripening profile.The late accumulation of lysophospholipid degradation products or their derivatives (phosphorylated head groups,free head groups,monoacylglycerols,other lysophosphatidic acids,and lysophospholipid derivatives) in dates results in a late ripening profile (Whitakeret al.2001;Coulonet al.2012;Dibounet al.2015).Therefore,in the present study,at the beginning of the 12 d postharvest period,the levels of all the differential lipids (most of which were unsaturated fatty acids,such as glyceryl linoleate,cis-10-heptadecenoic acid,and 12,13-EODE) decreased.On the contrary,during the second 12-d period,the levels of most differential lipids(most of which were lysophospholipid derivatives,such as LysoPE 18:2 (2n isomer),LysoPC 15:0,and LysoPE 16:0(2n isomer)) increased.In addition,pericarp browning and fruit disease appearance in postharvest longans are correlated with damage to the cellular membrane system,which entails cellular membrane phospholipid degradation and peroxidation of unsaturated fatty acids (Yiet al.2009;Linet al.2016;Wanget al.2018).The degradation and peroxidation of membrane lipids could cause damage to the structural integrity of cellular membranes and result in the disruption of cellular compartmentalization,thereby leading to the enzymatic browning process,owing to the contact of phenolic compounds with peroxidase and polyphenol oxidase (Linet al.2017;Wanget al.2018).Therefore,lipid metabolism also plays an important role in ‘Red Globe’ fruit postharvest softening and browning.

    Amino acids are not only the building blocks of proteins but also serve as precursors to some secondary metabolites,such as flavonoids,lignins,and phytoalexins (Yamakawa and Hakata 2010;Jinet al.2019).In the present study,the levels of most of the differential amino acids,including L-pyroglutamic acid,glutathione in its reduced form,oxidized glutathione,L-aspartic acid,and L-glutamic acid,but not 5-aminovaleric acid and L-asparagine anhydrous,decreased during postharvest senescence.Proteolysis has been shown to occur with a concomitant increase in free amino acid levels in the early phase of some postharvest crops(Matsumoto and Ikoma 2012;Yanet al.2019).However,in the present study,this phenomenon was not observed;free amino acid levels may increase within 12 d after harvest.In addition,both the oxidized and reduced glutathione levels decreased during the grape postharvest period.Glutathione and ascorbic acid are multifunctional metabolites playing important roles in the redox balance through the ascorbateglutathione cycle,suggesting that more reactive oxygen species were generated with aging,which would have caused further damage to the grape berry cells.

    Finally,the levels of pantothenol,a noteworthy metabolite identified in postharvest grapes,increased with senescence time.Therefore,it can be considered a characteristic component of grape berry senescence.Panthenol(pantothenol) is the alcohol analog of pantothenic acid(vitamin B5) and is,thus,a provitamin of B5.In organisms,it is quickly oxidized to pantothenate.Currently,it is mainly used in healthcare and cosmetic products,because it may suppress the growth of some harmful microbes (Chohnanet al.2014).

    Senescence of fruits involves several genetic networks where the phytohormone ethylene plays a key role,together with other hormones,integrating different signals and allowing the onset of conditions favorable for stage progression and organ longevity (Iqbalet al.2017).The free auxin increases during senescence and stimulates ethylene biosynthesis through its inductive action on the expression of the key enzyme 1-aminocyclopropane-1-carboxylic acid synthase (ACS).ABA can be a trigger for ethylene production and influence fruit ripening and senescence(Zhanget al.2009).Independently from ethylene,ABA can also increase all hydrolases to enhance softening (Lohaniet al.2004).Therefore,besides cell wall degradation,these hormones also enhance the fruit respiration rate,which may lead to the metabolite changes that were identified in the present study.However,further research should be undertaken to reveal the exact role of hormones in the postharvest senescence of grape fruit.

    5.Conclusion

    In the present study,we utilized widely targeted metabolomics using a UPLC-ESI-MS/MS detection platform to analyze the metabolic differences of ‘Red Globe’ grapes under postharvest senescence.Flavonoids included the highest number of differential metabolites,the content of which decreased as the postharvest senescence duration increased;and the same trend was observed for amino acids and their derivatives.Therefore,flavonoid and amino acid metabolism play key roles in grape fruit postharvest senescence.However,the levels of phenolic acids,procyanidins,and pantothenol,which are characteristic constituents of grapes,increased with increasing grape berry senescence duration.The lipid and nucleotide metabolite contents decreased during the early senescence phase but increased during the late senescence phase.A total of 27 flavonoids (e.g.,pelargonidin 3-O-glucoside,quercetin-3-O-glucoside,cyanidin 3-O-glucoside),4 phenolic acids (e.g.,trans-4-hydroxycinnamic acid methyl ester),L-aspartic acid and pantothenol can be considered as potential marker metabolites during the postharvest senescence of ‘Red Globe’ grape fruit.

    Acknowledgements

    The authors would like to thank Wuhan Metware Biotechnology Co.,Ltd.,China for providing the Metabolome Test Platform.This work was supported by the Natural Science Foundation of China (31801833 and 31801811),the Innovation Capability Support Programs of Shaanxi Province,China (2020KJXX-035),the China Postdoctoral Science Foundation (2019M653771 and 2019T120953),the Fundamental Research Funds for the Central Universities of China (2452019016),and the China Agriculture Research System of MOF and MARA (CARS-29-zp-6).

    Declaration of competing interest

    The authors declare that they have no conflict of interest.

    Appendicesassociated with this paper are available on http://www.ChinaAgriSci.com/V2/En/appendix.htm

    亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区| 桃红色精品国产亚洲av| 欧美国产精品一级二级三级| 岛国在线观看网站| 国产成人系列免费观看| 亚洲精品中文字幕一二三四区 | 999久久久精品免费观看国产| 国产主播在线观看一区二区| 纵有疾风起免费观看全集完整版| 午夜久久久在线观看| 老司机深夜福利视频在线观看| 男人舔女人的私密视频| 国产免费现黄频在线看| 操美女的视频在线观看| 中文字幕精品免费在线观看视频| 国产精品久久久av美女十八| 在线看a的网站| 亚洲成a人片在线一区二区| 亚洲色图综合在线观看| 自线自在国产av| 一二三四社区在线视频社区8| 久久精品91无色码中文字幕| √禁漫天堂资源中文www| 久久精品亚洲熟妇少妇任你| 久久精品熟女亚洲av麻豆精品| 男女之事视频高清在线观看| 精品午夜福利视频在线观看一区 | 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 精品少妇内射三级| 欧美av亚洲av综合av国产av| 欧美精品人与动牲交sv欧美| 久久久久精品国产欧美久久久| 国产精品久久久久久精品电影小说| 国产一区有黄有色的免费视频| 叶爱在线成人免费视频播放| 人妻久久中文字幕网| 两性午夜刺激爽爽歪歪视频在线观看 | 别揉我奶头~嗯~啊~动态视频| 国产成人av激情在线播放| 国产亚洲av高清不卡| 人妻 亚洲 视频| 久久久久久久精品吃奶| 中文字幕最新亚洲高清| 妹子高潮喷水视频| 黄色视频不卡| 菩萨蛮人人尽说江南好唐韦庄| 99国产精品99久久久久| 国产亚洲精品第一综合不卡| 天天操日日干夜夜撸| 欧美成人免费av一区二区三区 | 国产日韩欧美视频二区| 亚洲精品美女久久av网站| 香蕉国产在线看| 黑人欧美特级aaaaaa片| 免费在线观看影片大全网站| 久久精品人人爽人人爽视色| 亚洲成a人片在线一区二区| 搡老熟女国产l中国老女人| 极品少妇高潮喷水抽搐| 大型av网站在线播放| 91麻豆av在线| 女同久久另类99精品国产91| 蜜桃在线观看..| 高清毛片免费观看视频网站 | 国产男靠女视频免费网站| 12—13女人毛片做爰片一| 国产成人免费观看mmmm| 老司机午夜福利在线观看视频 | 国产高清videossex| 久久国产精品人妻蜜桃| 国产成人av激情在线播放| 在线观看免费高清a一片| 男男h啪啪无遮挡| av不卡在线播放| 另类精品久久| 人人妻,人人澡人人爽秒播| 亚洲成人免费av在线播放| 久久精品国产亚洲av香蕉五月 | 高清毛片免费观看视频网站 | 欧美人与性动交α欧美精品济南到| 50天的宝宝边吃奶边哭怎么回事| 日韩免费av在线播放| 国产亚洲欧美在线一区二区| av又黄又爽大尺度在线免费看| 欧美国产精品va在线观看不卡| 啦啦啦免费观看视频1| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频| 亚洲精品成人av观看孕妇| 国产91精品成人一区二区三区 | 国产亚洲精品第一综合不卡| 久久亚洲真实| 怎么达到女性高潮| 久久久国产欧美日韩av| 国产精品久久久久久人妻精品电影 | 欧美日韩亚洲高清精品| 俄罗斯特黄特色一大片| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品电影小说| av视频免费观看在线观看| 亚洲五月婷婷丁香| 9191精品国产免费久久| 成年动漫av网址| 汤姆久久久久久久影院中文字幕| 丝袜美足系列| 午夜福利乱码中文字幕| 麻豆国产av国片精品| 啦啦啦视频在线资源免费观看| 日本五十路高清| 国产免费福利视频在线观看| 精品第一国产精品| 97在线人人人人妻| 欧美精品一区二区大全| 久久青草综合色| 日韩欧美国产一区二区入口| 久久免费观看电影| 国产1区2区3区精品| 日韩欧美国产一区二区入口| 在线看a的网站| 国产高清videossex| 丝袜在线中文字幕| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| 黄网站色视频无遮挡免费观看| 久久久水蜜桃国产精品网| 亚洲伊人久久精品综合| 亚洲专区国产一区二区| 午夜福利在线免费观看网站| √禁漫天堂资源中文www| 精品一区二区三区视频在线观看免费 | 男男h啪啪无遮挡| 少妇精品久久久久久久| 日本av免费视频播放| 国产成人免费无遮挡视频| 男女午夜视频在线观看| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂| 精品亚洲成国产av| 久久性视频一级片| 99国产精品99久久久久| 一进一出好大好爽视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av成人一区二区三| a级片在线免费高清观看视频| 精品熟女少妇八av免费久了| 美女国产高潮福利片在线看| av有码第一页| 黑人猛操日本美女一级片| 妹子高潮喷水视频| 男女下面插进去视频免费观看| 国产精品一区二区在线不卡| 亚洲成人免费av在线播放| 最近最新免费中文字幕在线| 91老司机精品| 国产日韩欧美在线精品| 国产熟女午夜一区二区三区| 精品久久久精品久久久| 美女高潮喷水抽搐中文字幕| 免费不卡黄色视频| 建设人人有责人人尽责人人享有的| 国产午夜精品久久久久久| 久久久久久亚洲精品国产蜜桃av| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| 啦啦啦中文免费视频观看日本| 国产野战对白在线观看| 大香蕉久久成人网| 在线十欧美十亚洲十日本专区| 亚洲精品一二三| 黄色片一级片一级黄色片| av不卡在线播放| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 999精品在线视频| 精品国产乱码久久久久久男人| 免费av中文字幕在线| 丝袜在线中文字幕| 色精品久久人妻99蜜桃| 亚洲av成人一区二区三| 人妻一区二区av| 99热网站在线观看| 热re99久久精品国产66热6| 色94色欧美一区二区| 美女扒开内裤让男人捅视频| 欧美+亚洲+日韩+国产| 久久99一区二区三区| 国产高清videossex| 自线自在国产av| 啪啪无遮挡十八禁网站| 在线观看舔阴道视频| 欧美午夜高清在线| 亚洲国产av新网站| 一级黄色大片毛片| 国产亚洲精品久久久久5区| 丰满迷人的少妇在线观看| 丝瓜视频免费看黄片| 欧美精品人与动牲交sv欧美| 成人精品一区二区免费| a级片在线免费高清观看视频| 免费不卡黄色视频| 欧美日韩亚洲高清精品| 丝袜喷水一区| 十八禁网站免费在线| 俄罗斯特黄特色一大片| 一区二区av电影网| 国产成人av激情在线播放| 激情视频va一区二区三区| 在线观看一区二区三区激情| 国产男女超爽视频在线观看| 午夜精品久久久久久毛片777| 国产一区二区 视频在线| 亚洲精品成人av观看孕妇| 欧美变态另类bdsm刘玥| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| 亚洲va日本ⅴa欧美va伊人久久| tube8黄色片| 国产黄频视频在线观看| 激情在线观看视频在线高清 | 女警被强在线播放| 97人妻天天添夜夜摸| 性高湖久久久久久久久免费观看| 最新美女视频免费是黄的| 久久精品国产a三级三级三级| 亚洲精华国产精华精| 少妇猛男粗大的猛烈进出视频| 午夜福利一区二区在线看| 99久久国产精品久久久| 欧美亚洲日本最大视频资源| 久久精品成人免费网站| www日本在线高清视频| 99国产精品免费福利视频| 老司机午夜福利在线观看视频 | 高清毛片免费观看视频网站 | av又黄又爽大尺度在线免费看| 亚洲国产欧美日韩在线播放| 18在线观看网站| 国产在视频线精品| 午夜福利免费观看在线| 精品少妇久久久久久888优播| 日本欧美视频一区| 国产区一区二久久| 亚洲国产成人一精品久久久| 青草久久国产| 一级a爱视频在线免费观看| 精品国产一区二区三区久久久樱花| av天堂久久9| 免费黄频网站在线观看国产| 国产欧美日韩综合在线一区二区| 999久久久精品免费观看国产| 99国产精品一区二区蜜桃av | av天堂在线播放| 国产黄频视频在线观看| 黄色成人免费大全| 男女免费视频国产| 国产在视频线精品| 黄色视频,在线免费观看| 精品一区二区三区av网在线观看 | 一进一出抽搐动态| 12—13女人毛片做爰片一| 久久精品91无色码中文字幕| 一本大道久久a久久精品| 捣出白浆h1v1| 淫妇啪啪啪对白视频| 久久国产精品男人的天堂亚洲| 日韩有码中文字幕| 超碰97精品在线观看| 人成视频在线观看免费观看| 国产精品熟女久久久久浪| 久久精品人人爽人人爽视色| 国精品久久久久久国模美| 午夜福利一区二区在线看| 亚洲精品乱久久久久久| 精品人妻1区二区| 后天国语完整版免费观看| 黄色成人免费大全| 一级片免费观看大全| 丰满少妇做爰视频| 成人免费观看视频高清| 久久青草综合色| 一进一出好大好爽视频| 男女高潮啪啪啪动态图| 在线观看免费视频网站a站| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 国产色视频综合| 成年人午夜在线观看视频| 黄色视频,在线免费观看| 国产精品电影一区二区三区 | 嫁个100分男人电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 久久 成人 亚洲| 日韩欧美三级三区| 9色porny在线观看| 亚洲av美国av| 久久人人爽av亚洲精品天堂| 亚洲第一青青草原| 可以免费在线观看a视频的电影网站| 日韩一区二区三区影片| 18禁国产床啪视频网站| 男女之事视频高清在线观看| 波多野结衣av一区二区av| 性高湖久久久久久久久免费观看| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 国产精品.久久久| 亚洲精品乱久久久久久| a级毛片黄视频| 亚洲男人天堂网一区| 在线观看免费午夜福利视频| 免费一级毛片在线播放高清视频 | 成人国产一区最新在线观看| 日韩欧美免费精品| www.熟女人妻精品国产| 午夜精品国产一区二区电影| 免费观看av网站的网址| 国产aⅴ精品一区二区三区波| 亚洲av成人不卡在线观看播放网| 又大又爽又粗| 一区二区三区国产精品乱码| 一二三四在线观看免费中文在| 大片免费播放器 马上看| 亚洲美女黄片视频| 高清黄色对白视频在线免费看| 成年女人毛片免费观看观看9 | 久久精品91无色码中文字幕| 久久人妻av系列| 建设人人有责人人尽责人人享有的| 男女下面插进去视频免费观看| 无限看片的www在线观看| 国产精品自产拍在线观看55亚洲 | 丝袜在线中文字幕| 别揉我奶头~嗯~啊~动态视频| 精品国产超薄肉色丝袜足j| 亚洲av国产av综合av卡| 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 成人特级黄色片久久久久久久 | 欧美国产精品va在线观看不卡| 精品国产乱码久久久久久男人| 热re99久久国产66热| 黄色视频不卡| 亚洲av片天天在线观看| 叶爱在线成人免费视频播放| av不卡在线播放| 国产精品自产拍在线观看55亚洲 | 天天躁日日躁夜夜躁夜夜| 成人永久免费在线观看视频 | 80岁老熟妇乱子伦牲交| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 久久久久久久久久久久大奶| 97人妻天天添夜夜摸| 夜夜爽天天搞| 国产视频一区二区在线看| 国产精品 国内视频| 99久久99久久久精品蜜桃| 99国产极品粉嫩在线观看| 久久精品国产a三级三级三级| 久久99一区二区三区| 亚洲精品久久午夜乱码| 嫩草影视91久久| 亚洲av国产av综合av卡| 夜夜骑夜夜射夜夜干| 国产精品九九99| 国产精品 欧美亚洲| 丁香欧美五月| 日本一区二区免费在线视频| 亚洲午夜理论影院| 日本黄色日本黄色录像| 多毛熟女@视频| 亚洲国产成人一精品久久久| 最近最新中文字幕大全免费视频| 9色porny在线观看| 精品一区二区三卡| 两性夫妻黄色片| av不卡在线播放| 一二三四社区在线视频社区8| 亚洲欧洲精品一区二区精品久久久| 99久久人妻综合| 国产精品偷伦视频观看了| 国产精品98久久久久久宅男小说| 免费观看人在逋| 在线观看舔阴道视频| av有码第一页| 自拍欧美九色日韩亚洲蝌蚪91| 97在线人人人人妻| 在线观看免费日韩欧美大片| 在线 av 中文字幕| 久久中文字幕人妻熟女| 亚洲精品av麻豆狂野| 91麻豆精品激情在线观看国产 | 国产男靠女视频免费网站| 中文字幕人妻熟女乱码| 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| 亚洲av片天天在线观看| 亚洲精品一二三| √禁漫天堂资源中文www| 国产一区有黄有色的免费视频| 免费人妻精品一区二区三区视频| 天天躁日日躁夜夜躁夜夜| 亚洲精品一二三| 成年女人毛片免费观看观看9 | 国产精品98久久久久久宅男小说| 肉色欧美久久久久久久蜜桃| 黄色a级毛片大全视频| av欧美777| 成人国语在线视频| 久久精品国产a三级三级三级| 99精品久久久久人妻精品| 不卡av一区二区三区| 考比视频在线观看| 日本一区二区免费在线视频| 亚洲av电影在线进入| 18禁裸乳无遮挡动漫免费视频| 12—13女人毛片做爰片一| 狠狠婷婷综合久久久久久88av| 国产男女内射视频| 成人特级黄色片久久久久久久 | 男女免费视频国产| av国产精品久久久久影院| 涩涩av久久男人的天堂| 国产精品久久久久久人妻精品电影 | 色视频在线一区二区三区| 国产精品影院久久| 嫁个100分男人电影在线观看| 五月天丁香电影| 欧美日韩中文字幕国产精品一区二区三区 | 成年人黄色毛片网站| 国产精品99久久99久久久不卡| 午夜福利,免费看| 国产麻豆69| 成人特级黄色片久久久久久久 | 亚洲精品中文字幕一二三四区 | 国产精品九九99| 搡老乐熟女国产| 日韩制服丝袜自拍偷拍| 亚洲av欧美aⅴ国产| 淫妇啪啪啪对白视频| 亚洲成人免费av在线播放| 男女边摸边吃奶| 女人爽到高潮嗷嗷叫在线视频| 国产欧美亚洲国产| 丰满迷人的少妇在线观看| 精品国产一区二区三区四区第35| 99热网站在线观看| 大型黄色视频在线免费观看| 黄片大片在线免费观看| 精品久久久精品久久久| 久久av网站| 成人国产av品久久久| 国产欧美日韩一区二区三| 女警被强在线播放| 9热在线视频观看99| 久久中文字幕人妻熟女| 脱女人内裤的视频| 亚洲成人国产一区在线观看| 日韩中文字幕视频在线看片| 黑人巨大精品欧美一区二区mp4| 成人亚洲精品一区在线观看| 亚洲欧洲日产国产| 久久久久久久国产电影| 久久99一区二区三区| 精品国内亚洲2022精品成人 | 五月天丁香电影| 精品久久久久久久毛片微露脸| 精品少妇久久久久久888优播| 国产在线观看jvid| 久久九九热精品免费| 午夜福利影视在线免费观看| 久久中文看片网| 在线av久久热| 亚洲欧美激情在线| 欧美激情高清一区二区三区| 黑人欧美特级aaaaaa片| 日韩视频在线欧美| 国产成人欧美在线观看 | 久久久久久久大尺度免费视频| 国产精品 国内视频| 色播在线永久视频| 亚洲,欧美精品.| 女人精品久久久久毛片| 不卡一级毛片| 男女下面插进去视频免费观看| 精品国内亚洲2022精品成人 | 97人妻天天添夜夜摸| 黄色视频,在线免费观看| 久久国产亚洲av麻豆专区| 亚洲成人国产一区在线观看| 精品久久蜜臀av无| 老熟妇仑乱视频hdxx| 亚洲精华国产精华精| 欧美乱码精品一区二区三区| 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区 | 亚洲精品久久成人aⅴ小说| 无人区码免费观看不卡 | 在线看a的网站| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 在线av久久热| 久久久久久久大尺度免费视频| av超薄肉色丝袜交足视频| 一区二区三区乱码不卡18| 高清视频免费观看一区二区| 久久久国产一区二区| 成人国产一区最新在线观看| 桃红色精品国产亚洲av| 国产老妇伦熟女老妇高清| 美女高潮到喷水免费观看| 成人18禁高潮啪啪吃奶动态图| 欧美在线黄色| 国产精品亚洲av一区麻豆| 一个人免费看片子| 亚洲精品乱久久久久久| 男女午夜视频在线观看| 三级毛片av免费| 国产无遮挡羞羞视频在线观看| 久久亚洲真实| 国产成+人综合+亚洲专区| 国产亚洲欧美在线一区二区| 9色porny在线观看| 欧美黑人精品巨大| 首页视频小说图片口味搜索| 婷婷成人精品国产| 久久天躁狠狠躁夜夜2o2o| 757午夜福利合集在线观看| 午夜福利在线免费观看网站| 精品国产乱子伦一区二区三区| 欧美变态另类bdsm刘玥| 可以免费在线观看a视频的电影网站| 亚洲av成人不卡在线观看播放网| 亚洲五月婷婷丁香| 欧美日韩精品网址| 精品国产乱码久久久久久男人| 男男h啪啪无遮挡| 色94色欧美一区二区| 在线观看免费高清a一片| 女性生殖器流出的白浆| 精品少妇内射三级| 俄罗斯特黄特色一大片| 不卡av一区二区三区| 视频在线观看一区二区三区| 欧美日韩黄片免| 在线观看免费午夜福利视频| 成人av一区二区三区在线看| 国产精品美女特级片免费视频播放器 | 中文字幕色久视频| 国产日韩欧美在线精品| 热99re8久久精品国产| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 中文亚洲av片在线观看爽 | 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 国产午夜精品久久久久久| 亚洲欧美日韩高清在线视频 | 日韩中文字幕视频在线看片| 日韩制服丝袜自拍偷拍| 日本一区二区免费在线视频| 国产欧美日韩一区二区三区在线| 如日韩欧美国产精品一区二区三区| 国产视频一区二区在线看| 一级,二级,三级黄色视频| 久久久久网色| 亚洲一码二码三码区别大吗| 一本大道久久a久久精品| 大型av网站在线播放| 精品福利观看| 国产精品 国内视频| 午夜激情av网站| 一边摸一边抽搐一进一小说 | 超碰成人久久| 狠狠婷婷综合久久久久久88av| 国产亚洲精品一区二区www | 色94色欧美一区二区| 久久国产亚洲av麻豆专区| 欧美激情极品国产一区二区三区| 亚洲 欧美一区二区三区| 精品少妇黑人巨大在线播放| 国产成人av激情在线播放| 高清视频免费观看一区二区| 国产精品一区二区免费欧美| 黄色视频不卡| 热re99久久国产66热| 黄色视频不卡| 国产成人精品久久二区二区免费| 色94色欧美一区二区| 亚洲国产欧美网| 久久中文字幕人妻熟女| www日本在线高清视频| 欧美日韩成人在线一区二区| 好男人电影高清在线观看| 老司机午夜福利在线观看视频 | 久久天躁狠狠躁夜夜2o2o| 午夜两性在线视频| 日本一区二区免费在线视频| 日本欧美视频一区| 色视频在线一区二区三区| 免费观看人在逋| 建设人人有责人人尽责人人享有的| 一边摸一边抽搐一进一小说 | 99国产极品粉嫩在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产主播在线观看一区二区| 亚洲精品久久午夜乱码| 亚洲性夜色夜夜综合| 国产精品久久久久久精品电影小说| 国产精品久久久久久人妻精品电影 |