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

    Prediction of meat quality traits in the abattoir using portable near-infrared spectrometers:heritability of predicted traits and genetic correlations with laboratory-measured traits

    2021-09-19 13:14:34SimoneSavoiaAndreaAlberaAlbertoBrugiapagliaLilianaDiStasioAlessioCecchinatoandGiovanniBittante

    Simone Savoia,Andrea Albera ,Alberto Brugiapaglia ,Liliana Di Stasio ,Alessio Cecchinato and Giovanni Bittante

    Abstract Background:The possibility of assessing meat quality traits over the meat chain is strongly limited,especially in the context of selective breeding which requires a large number of phenotypes.The main objective of this study was to investigate the suitability of portable infrared spectrometers for phenotyping beef cattle aiming to genetically improving the quality of their meat.Meat quality traits (pH,color,water holding capacity,tenderness) were appraised on rib eye muscle samples of 1,327 Piemontese young bulls using traditional (i.e.,reference/gold standard) laboratory analyses;the same traits were also predicted from spectra acquired at the abattoir on the intact muscle surface of the same animals 1 d after slaughtering.Genetic parameters were estimated for both laboratory measures of meat quality traits and their spectra-based predictions.Results: The prediction performances of the calibration equations,assessed through external validation,were satisfactory for color traits (R2 from 0.52 to 0.80),low for pH and purge losses (R2 around 0.30),and very poor for cooking losses and tenderness (R2 below 0.20).Except for lightness and purge losses,the heritability estimates of most of the predicted traits were lower than those of the measured traits while the genetic correlations between measured and predicted traits were high (average value 0.81).Conclusions:Results showed that NIRS predictions of color traits,pH,and purge losses could be used as indicator traits for the indirect genetic selection of the reference quality phenotypes.Results for cooking losses were less effective,while the NIR predictions of tenderness were affected by a relatively high uncertainty of estimate.Overall,genetic selection of some meat quality traits,whose direct phenotyping is difficult,can benefit of the application of infrared spectrometers technology.

    Keywords:Genetic parameters,Meat quality,Near-infrared spectroscopy,Piemontese

    Background

    Improving meat quality attributes through genetic selection proved to be theoretically feasible as many quality traits display moderate to medium heritability values [1].However,establishing a selection procedure depends on the availability of phenotypes collected as part of a routine recording scheme.This is a serious limitation when it comes to meat quality traits as it currently requires meat samples to be collected at the slaughterhouse,which depreciates the carcasses,while the subsequent laboratory analyses are expensive,and time consuming.

    Visible and near-infrared spectroscopy (Vis-NIRS),which is based on the principle that different chemical bonds in organic matter absorb or emit light of different wavelengths when the sample is irradiated,offers a number of important advantages over conventional methods:ability to take rapid and frequent measurements,fast and simple or no sample preparation,suitability for online use,ability to simultaneously determine different attributes [2].

    Several studies have assessed the use of reflectance spectroscopy to accurately predict the chemical composition of beef [3–5] and,with lower accuracy,different meat quality attributes [6–8].

    However,with regards to genetic improvement,scientific knowledge is very scarce in beef cattle although genetic parameters of NIRS predictions of meat quality traits have been estimated in pigs [9].In a previous study by our group [10],the only one that we are aware of to have made a genetic comparison between laboratorymeasured and laboratory infrared-predicted meat quality traits of beef samples,we found medium-high genetic relationships between some of the measured and corresponding predicted meat quality traits.The genetic correlations for all the color traits and purge losses were high,and were greater than the corresponding phenotypic correlations,whereas both the phenotypic and genetic correlations for tenderness and cooking losses were negligible.These findings suggest the feasibility of genetically improving some meat quality traits using their NIR spectrometry predictions from meat samples.However,a selection program for meat quality traits could be more easily established if it were possible to routinely record phenotypes at the slaughterhouse without having to collect samples.In our previous study [10],meat samples were taken in the abattoir and after aging were transported to the laboratory where the muscle portions were dissected and minced;spectra were then acquired using a bench-top NIR spectrometer and laboratory analyses were carried out on the same sample,on the same day,and in the same laboratory.

    The availability of new,portable NIR and Vis-NIR spectrometers able to collect spectra directly from the muscle surface at the slaughterhouse [7] means that we now have more efficient phenotyping tools for use in selection programs to improve meat quality traits.

    In the absence of large-scale studies on the prediction of meat quality using portable NIR spectrometers,a research project (Qualipiem project) was set up with the aims of analyzing meat quality traits in Piedmontese young bulls,and of proposing innovative selection strategies for their improvement.The first steps taken were:to evaluate beef farming systems and other phenotypic sources of variation in carcass and meat quality traits using gold standard laboratory analyses [11];to estimate their quantitative genetic parameters [12];and to study their genome-wide associations and carry out pathway analyses [13].In addition,two spectrometers very different in size,ease of use,and cost were compared for their ability to predict meat quality on the muscle surface in the abattoir without the need for meat sampling [14].The obtained prediction accuracies showed a large variation,ranging from high values for most of the color traits to low values for meat tenderness and cooking losses.However,the prediction performance of calibration models was penalized because meat samples were affected by many processing factors prior to laboratory analyses (i.e.carcass aging and dissection,samples uptaking,transport and processing),which could not be predicted by spectra taken at the slaughterhouse.The conclusion of that study was that NIRS predictions of meat quality traits could be useful to capture the animal’s“native”characteristics,which is the case for the genetic improvement [14].

    The main objective of this study was,therefore,to investigate the suitability of portable infrared spectrometers for phenotyping beef cattle as a basis for genetically improving the quality of their meat.The specific aims were:1) to analyze the genetic variation in meat quality traits as predicted by two very different portable spectrometers on the intact cross-sectional muscle surface in the abattoir;2) to compare this with the genetic variation in traits measured in the laboratory after sampling,aging,transport and analysis;3) to assess the genetic relationships between laboratory measures of meat quality traits and their spectrum-based predictions.

    Methods

    Animals

    The study was carried out on samples from 1,327 Piemontese young bulls slaughtered at the same commercial abattoir over 106 slaughter days.The young bulls were progeny of 204 A.I.purebred sires and 1,286 dams,all registered in the Italian Piemontese Herd Book.

    The animals were fattened on 115 farms representative of the beef production systems in the Piemonte region(north-west Italy).The beef farming systems,feeding regimes,fattening conditions and slaughter performances of the young bulls are described in detail by Savoia et al.[11].In brief,the young bulls were reared on farms operating one of the following beef production systems:traditional with restricted feeding and animals either kept in tie stalls or loose housed;modern breeders and fatteners or specialized fatteners using ad libitum feeding and loose housing (the last two systems were further subdivided according to whether or not they used total mixed rations).

    The sampled young bulls had an average carcass weight of 438.1±43.6 kg,and an average age at slaughter of 541±63 days,giving an average daily carcass gain of 0.818±0.107 kg/d.Average carcass conformation score(using the SEUROP classification system with each category divided into 3 subclasses giving a scale of 18 points) was 14.66,corresponding to an average evaluation close to“E+”in the EU linear grading system;the average rib eye area measured at the 5th rib was 92±14.3 cm2.

    Spectra collection

    Spectra collection and the technical characteristics of the instruments are described in detail by Savoia et al.[14].Briefly,the spectra were acquired with the following spectrometers:

    · Vis-NIRS:LabSpec 2500(ASD Inc.,Boulder,CO,USA),which has a spectral range in the visible and near-infrared sections of electromagnetic radiation(wavelengths 350 to 1830 nm),measurements taken every 1 nm producing 1481 data points per sample;the instrument’s dimensions are 12.7 cm×36.8 cm×29.2 cm,and it weighs 5,600 g;the spectra are collected with a probe (26 cm ×10 cm ×5 cm) connected to the instrument with an optical fiber;

    · Micro-NIRS:Micro NIR Pro (JDSU San Jose,CA,USA),which has a spectral range in the nearinfrared region (wavelengths 905 to 1,649 nm),measurements taken every 6 nm producing 125 data points per sample;the instrument’s dimensions are 4.5 cm ×4.4 cm ×4.0 cm,and it weighs 60 g;the spectra are collected directly by the instrument,which should be connected to a lap-top or tablet via a USB cable.

    The right side of each carcass was divided into 2 quarters between the 5th and 6th ribs(pistol cut) in the abattoir the day after slaughter (about 24 h post-mortem).The spectra were collected on the cross-sectional surface of the longissimus thoracis muscle using the scanning head of the fiber-optic contact probe (10 mm in diameter) of the Vis-NIRS or by applying the Micro-NIRS to the surface of the muscle.Five spectra were acquired with each instrument from different positions on the cut surface of the same muscle.

    Beef sample collection and meat quality analyses

    Twenty-four hours after slaughter and immediately following spectra collection,individual samples (4.0 cm thick) of the longissimus thoracis muscle were collected from between the 5th and 6th ribs,then were individually vacuum packed and transferred under refrigerated conditions to the laboratory,where they were stored in a chilling room at 4°C for 6 d (meat aging),after which meat quality traits were measured on all samples.

    Meat quality was assessed 7 d after slaughter by measuring:muscle pH in triplicate using a digital pH-meter;lightness (L*),redness (a*),yellowness (b*),hue angle(h*),chroma (C*) in triplicate after 1 h blooming using a Minolta CR-331C colorimeter;purge losses (PL,%) by difference between the weight of the sample cut at carcass dissection (day 1 after slaughtering) and before analyses (day 7);cooking losses (CL,%) by difference between of the meat subsample before and after cooking meat sample in a water-bath till an internal temperature of 70°C;tenderness(Warner Bratzler shear force;WBSF,N) on 6 cylindrical 1.27 cm-diameter cores of cooked meat.Details of the procedures used to assess the meat quality traits can be found in Savoia et al.[11].

    Statistical analyses

    Spectral data editing and validation procedure

    The spectral data were edited and processed according to the model described in detail in the previous study by Savoia et al.[14].In brief,records with errors (e.g.,individually identified spectra not matching the reference samples) and outlier spectra identified by Mahalanobis distance were discarded from the two original spectral datasets obtained with the Vis-NIR and Micro-NIR portable spectrometers.Before developing the calibration equations,the spectral data were centered and standardized to improve the goodness of fit of the chemometric modeling.

    A Bayesian model (Bayes B) implemented with the BGLR library [15] of the R-software was used to develop calibration equations for each beef quality trait,as described in detail by Ferragina et al.[16].

    In order to reproduce operational conditions,an external validation procedure was carried out to predict meat quality traits from the calibration equations based on the Vis-NIR and Micro-NIR spectra.As the most important source of variation was the slaughter batch (all animals slaughtered on the same day),external validation consisted in predicting the observations of a given batch from the calibration equations developed using the measured meat quality data of animals slaughtered on all the other days (“l(fā)eave-one-batch-out”procedure).This procedure has been repeated for each batch accumulating the predicted values.Therefore,the predictions of meat quality data for each animal were obtained using prediction equations developed without the observations measured on that animal and on all the animals slaughtered in the same batch.As a consequence,the entire set of predicted data could be used for the estimation of genetic parameters without any risk of inflation of the estimates.The R2values of external validation have been used to evaluate the prediction accuracy of calibration models.The“l(fā)eave-one-batch-out”procedure’s details are described in Savoia et al.[14].

    The final dataset used to estimate (co)variance components and to evaluate the magnitude of the genetic correlations contained the measured and predicted observations on a number of animals ranging from 1,117 to 1,134 depending on the trait.

    Estimates of (co)variance components and genetic parameters

    (Co)variance components were estimated by REML procedures using the VCE software [17].For each of the meat quality traits,(co)variance components were estimated through separate bivariate analyses including the measured trait and its prediction obtained with either the Vis-NIR or Micro-NIR spectrometer.

    The general model can be written in matrix notation as:

    where y contains the observations for traits 1 and 2,β is the vector of non-genetic fixed effects,c is the vector of random herd effects (98 levels),q is the vector of the random batch effect (106 levels),u is the vector of animal additive genetic effects,e is the vector of random residual effects,and X,W1,W2 and Z are incidence matrices of proper dimensions.Fixed effects of parity of the dam (4 classes:1st,2nd,3rd–8th,>8th) and age at slaughter (5 classes:<450 d,450–510 d,511–570 d,571–630 d,>630 d) have been included in the model for pH and L*,respectively.Effects of different herds were assumed to be normally and independently distributed,c~N(0,C ?I);effects of batch were assumed to be normally and independently distributed,q~N(0,Q ?I).A minimum cell size of 3 observations was required for both herd and batch effects.Animal additive genetic effects were assumed to be normally distributed,u~N(0,G ?A),where G is the (co)variance matrix between the animal effects,and A is the numerator of the Wright’s relationship matrix.Additive relationships were computed using a pedigree file that included the phenotyped animals and all their known ancestors (13,122 animals).Residuals were assumed to have a normal distribution,e~N(0,R ?I).

    To facilitate comparison with literature estimates,we estimated intra-herd heritability defined as:

    Results

    Descriptive statistics of the meat quality traits measured in the laboratory on aged samples according to gold standard methods,and their predictions obtained from the spectra acquired in the abattoir the day after slaughter are reported in Table 1.

    Large differences in variability across the traits measured in the laboratory were observed.Water loss traits(PL and CL) and WBSF showed the highest variability followed by color traits,whereas the SD of the pH measurements was very low.For all the traits considered,the average values of the predictions obtained with both instruments were very similar to the corresponding laboratory measurements.On the other hand,the variability in the predicted traits was always much lower than in the measured traits.This was particularly marked for PL,CL,and WBSF,where the standard deviation of the predictions was 50% to 78% lower than that of the measured traits.The reduction was less pronounced (?10%to ?27%) for color traits.Loss of variability was in general more marked in the predictions obtained with the Micro-NIRS spectrometer (?40%) than with the Vis-NIRS instrument (?30%).The externally validated prediction performance of the calibration equations was satisfactory for all color traits (0.52 to 0.80),low for pH and PL (around 0.30),and very poor for CL and WBSF (below 0.20).No relevant differences were observed between the two spectrometers in terms of the magnitude ofAcross traits,there was a clear relationship between loss of variability in the predictions with respect to the measurements and the quality of the prediction performance (R20.96 for Vis-NIRS,R20.98 for Micro-NIRS,Fig.1).

    Fig. 1 Relationship between R2 of external validations of predicted meat quality traits and the decrease in phenotypic variance(Δ)of predictions compared with the laboratory measured traits

    Table 2compares laboratory-measured and spectrapredicted traits with respect to the variance components and heritabilities of color traits.The batch effect was the most important source of variation for all traits,with the exception of L*,accounting for 15% to 30% of the total variance.The incidence of this effect was always lower in the Micro-NIRS predictions than in the Vis-NIRS predictions and the laboratory-measured traits.The effect of the fattening herd was small (5% to 10% of totalvariance according to the trait) and relatively homogeneous for the laboratory-measured and spectra-predicted traits.

    Table 1 Descriptive statistics of meat quality traitsa and of their predictions,and calibration equations prediction performances

    Table 2 Variance components and intraherd heritability of color traitsa and of Vis-NIRS and Micro-NIRS predictions

    The animal additive genetic effect explained 20% to 30% of the total variance of L*,according to the instrument,with a higher incidence in the predicted than in the measured traits for both instruments.As a consequence,the intraherd heritability values were relatively high,ranging from 0.30 in the lab measurements to 0.41 in the Vis-NIRS predictions.All the other color traits behaved differently:the proportion of variance in the animal effect was around 10% and consistent across the measured traits,and was much lower with both prediction techniques.Heritability values for the predictions were therefore quite low compared with the heritability of 0.14 for the laboratory measured traits,and fairly homogeneous,with values ranging from 0.04 to 0.08 with the only exception of h* from the Micro-NIRS spectrometer.

    The variance of the batch effect was very high for measured pH,CL and WBSF,and for most of the corresponding predictions,accounting for between 40 and 60% of total variance (Table 3).Only the Micro-NIRS predictions of CL exhibited a small variance in batch(13%).The variances in batch of both the measured and predicted PL were similar to those of color traits.Likewise,there was little variability in the meat quality traits due to the herd effect,in most of the cases not exceeding 7% of the total variance.

    Table 3 Variance components and intraherd heritability of meat quality traitsa and of Vis-NIRS and Micro-NIRS predictions

    The proportion of variance explained by the additive genetic effect was much higher in the measured CL and WBSF than in their predictions.The heritabilities of the measured traits were moderate for CL (0.19) and relatively high for WBSF (0.31),but those of the predicted traits were considerably lower.In particular,the predictions of CL obtained with Micro-NIRS and WBSF obtained with the Vis-NIRS instrument showed an almost null incidence of additive genetic variance,and theresulting heritabilities were close to zero.For both pH and PL,the estimated heritabilities of the predictions obtained from the Vis-NIRS instrument were higher than those from the Micro-NIRS (0.18 vs.0.13 for pH,0.22 vs.0.13 for PL).The heritability of measured pH was higher than that of the corresponding predictions,whereas the heritability of measured PL was similar to that of the Micro-NIRS predictions and markedly lower than that of the Vis-NIRS predictions.

    Estimates of the genetic and residual correlations obtained by bivariate analyses of color and meat quality traits measured on aged meat samples in the laboratory,and their predictions obtained from meat spectra acquired in the abattoir by both the Vis-NIRS and Micro-NIRS are presented in Table 4.The values of the residual correlations reflect the prediction performance of the calibration equations.Genetic correlations between the lab-measured and spectra-predicted traits were always higher than the corresponding residual correlations.Their average value across traits was 0.81 compared with values of around 0.50 for the residual correlations for both instruments.Large differences in the genetic correlations among traits were observed.These were extremely high for color traits and PL,almost always 1.0 for the Vis-NIRS,and on average 0.9 for the Micro-NIRS.Among the other traits,the genetic correlations were of a lower magnitude and differed between the two spectrometers.The estimated genetic correlations obtained from the Vis-NIRS were greater than those from the Micro-NIRS,particularly for pH (0.70 vs.0.45),CL (0.70 vs.0.25) and WBSF (0.81 vs.0.42).

    Table 4 Additive genetic(ra) and residual correlations (re) of color and meat quality traits with their predictions

    Overall,the prediction performances of the two spectrometers were quite similar for all the traits considered.The genetic parameters of the predictions for color traits were comparable,whereas the Vis-NIRS performed better than the Micro-NIRS in the other meat quality traits.

    Discussion

    The main objective of this research was to evaluate the possible use of NIRS technology to predict phenotypes for meat quality traits for genetic evaluation purposes.The use of portable instruments at the slaughterhouse,and spectra acquisition from the muscle surface naturally exposed during routine subdivision of half-carcasses into quarters could eliminate the need for meat sample collection,which depreciates the carcasses,and for the subsequent transport,aging and laboratory analyses.In a recent work from the same project,we described,compared and discussed in detail the Vis-NIR and Micro-NIR spectroscopic techniques in terms of instrument characteristics,repeatability,calibration,cross-validation,and external validation of the predictions,and their ability to capture the main phenotypic sources of variation in meat quality traits [14].In this study,discussion of the two spectroscopic techniques was focused on estimates of genetic parameters and their use for the genetic improvement of meat quality in beef cattle populations.In our investigation the predictive performance of NIRS was evaluated using external validation,which led to R2values lower than those in most of the published studies using random cross-validation.In other foods,Bittante et al.[18] and Eskildsen et al.[19] showed that the random cross-validation procedure tended to overestimate the predictive ability of FTIR predicted traits compared with external validation.This over-fitting is particularly large for traits affected by high environmental variation related to farms,seasons,or batches.Calibration model parameters,such as cross or external validation (R2or RMSE),are not sufficient to establish the usefulness of spectral predictions for the purposes of genetic improvement.Although Soyeurt et al.[20] consider an R2of cross-validation exceeding 0.75 to be necessary in order to use predictions of dairy traits for selection,some authors report satisfactory results with moderate or even low prediction performances for milk[18,21] and meat quality traits [10].Aside from R2,the suitability of infrared predictions as indicator traits for selection relies especially on a combination of their heritability and loss of additive genetic variance with respect to measured traits and their genetic correlations with the corresponding measurements [21,22].

    Heritability of measured and predicted meat quality traits

    The heritabilities of the laboratory-measured meat quality traits were in the range of most literature reports [23,24],and were extensively discussed in our previous study using this dataset [12].In a study on the genetic basis of meat quality traits in the Piemontese breed,Boukha et al.[1] found similar heritabilities to this study for L*,b* and C*,but their estimates were considerably higher for a* and h*.However,different tools were used to measure color traits in that study.The heritabilities of the other measured meat quality traits were in some cases lower (CL and WBSF),and in other cases higher(pH and PL) than in our study.The theoretical basis of the heritability of infrared predictions lies not only in the phenotypic correlations with meat quality traits (prediction accuracy),but also in the fact that the absorbance of electromagnetic radiation by the meat surface may be affected by the animal’s genetics,i.e.,it may be heritable.Different sections of the milk infrared spectrum have been shown to have variable degrees of heritability [25,26],but this has not so far been shown for meat.

    In general,the predictions of meat quality traits obtained from spectra acquired in the abattoir yielded lower heritability values than the corresponding traits measured on aged meat samples in the laboratory,with the exception of L* and PL.The difference in heritabilities between the predicted and measured traits was on average 0.08 for both spectrometers,but there was variability across traits.Although the heritabilities of the meat quality predictions were lower than the heritabilities of the measured traits,in most cases they were large enough to be exploited for selection.Only the predictions of CL obtained with Micro-NIRS,and the predictions of WBSF obtained with Vis-NIRS were of no use,as their additive genetic variance was almost null.So far,our previous survey has been the only one to address the estimation of genetic parameters for meat quality traits predicted by NIR spectroscopy [10].That investigation was carried out on cattle of the same breed and sex,and similar age to those in this study.However,the spectra were acquired from minced,aged samples with a bench-top spectrometer in laboratory conditions,and the calibration equations were developed using a different procedure based on partial least squares regression and random assignment of samples to the calibration and testing sets.Unlike the present study,the predictions of all the meat quality traits,except L*,had similar or greater heritabilities compared with the measured traits [10].In general,the estimates of heritability of the predicted traits were higher than in the present study for most of the traits,and no null heritabilities were found also for CL and WBSF predictions.It should be borne in mind that in the previous study the predictions were obtained from spectra acquired on the same day and on the same material used for the meat analyses,and that the validation data were not entirely independent of the calibration data.In the present study,the spectra were acquired in the abattoir 6 d before the laboratory analyses,so these calibrations predict the quality of the meat after sampling,aging,transport,and analysis.

    In this investigation,the lower heritability of the predictions compared with the measurements was due to the additive genetic variance being on average 70%lower,which exceeded the decrease observed in the phenotypic variance,ranging from 50% to 60%,depending on the instrument.These results are consistent with those of Bonfatti et al.[22],who found a similar pattern in the infrared predictions of a large number of traits related to milk composition and technological properties at the population level.On the other hand,other studies on the same traits reported that the lower additive genetic variance in the infrared predictions compared with the measured traits was also associated with a lower residual variance leading to increased heritability values,particularly for poorly predicted traits [18,27].

    The lower phenotypic variability in the meat quality predictions observed in this study was strongly related to the accuracy of the calibrations measured by the R2of the external validations (Fig.1).As expected,using predictions instead of original traits results in lower variability,which is directly related to the predictive performance of the model adopted.The results were almost identical for the two spectrometers,and were in agreement with the findings of Cecchinato et al.[10].On the other hand,there was not such a close association between the decrease in additive genetic variability and the R2of the external validations of the calibration models,particularly with the Vis-NIR spectrometer(Fig.2).As a consequence,neither the heritabilities of the predictions nor their loss in additive genetic variation compared with the measured traits displayed a consistent relationship with the predictive performance of the calibration models.These findings,which are in agreement with Cecchinato et al.[10],confirm that meat quality traits predicted from infrared spectra with moderate or even low accuracy may also be heritable,and that there is potential to exploit the genetic variability.

    Fig. 2 Relationship between R2 of external validations of predicted meat quality traits and the decrease in additive genetic variance(Δσ2a)of predictions compared with the laboratory measured traits

    Correlations between laboratory-measured and abattoirpredicted meat quality traits,and their possible use for genetic improvement

    The genetic correlations between the infrared-predicted and laboratory-measured traits are important in determining the effectiveness of their use as indicator traits for selective breeding [21,27].The degree of genetic gain achievable with indirect selection is affected by the genetic correlations between the desired and indicator traits [28].For all the color traits,the genetic correlations between the measured traits and the Vis-NIRS predictions were extremely high,as were those between the measured traits and the Micro-NIRS predictions,despite the absence of the visible wavelengths of the spectrum.These results are consistent with the findings of the previous study [10],where the genetic correlations between the bench-top NIRS predictions of meat color traits and the corresponding measurements ranged from 0.85 to 0.99.The genetic correlations between measured and Vis-NIRS predictions for a*,b*,C* and h* were equal to 1 despite the reduced σ2afor the predictions.

    For PL,the estimated genetic correlations obtained in the present study with both instruments were substantial,and larger than in the previous study.Similarly,the Vis-NIRS predictions of pH had rather high genetic correlations with the measurements,whereas the corresponding correlations with the Micro-NIRS predictions were only moderate.In the previous study [10],CL and WBSF also proved to be difficult traits to predict from NIR spectra from the perspective of genetic improvement,as both the estimated phenotypic and genetic correlations with their laboratory measurements were inconsistent,even though the spectra were acquired from the same material,at the same location and on the same day as the laboratory analyses.These traits exhibited a partially different pattern in the present study.Although the predictive abilities of the calibration models were low,the Vis-NIRS predictions correlated well with the measured traits from the genetic standpoint,whereas the corresponding correlations with the Micro-NIRS were weaker.However,given the considerable decrease in the additive genetic variance in the Vis-NIRS predictions of WBSF,only the predictions of CL obtained from this instrument seem to be useful for selection purposes.

    As with the earlier findings [10],the genetic correlations between the predicted and measured traits were positively associated with the predictive ability of the calibration models (Fig.3),but to a lesser extent than the phenotypic or residual correlations.Traits that are very accurately predicted by calibration models always show high genetic correlations with their measurements,while there is greater variability in the genetic correlations when the prediction performance is moderate or poor [22].

    Fig. 3 Relationship between R2 of external validations of predicted meat quality traits and the genetic correlation(ra) of predictions with the laboratory measured traits

    The use of infrared predicted traits at the population level for genetic purposes has been shown to be possible for some milk traits [29],but the use of phenotypic predictions obtained with imprecise methods is the subject of much debate.It worth noting that,in response to criticisms of a technique for predicting enteric methane emissions,Bovenhuis,van Engelen,and Visker [30] recently stated that“… even if measurements are inaccurate,imprecise,or biased,they might provide valuable information for selective breeding”,and that“When given the choice,accurate and unbiased measurements are preferred.However,such measurements are seldom available on a large scale and at reasonable cost.…However,inaccurate and biased sniffer methane phenotypes do not automatically imply inaccurate and biased methane breeding values”.When evaluating the usefulness of infrared predictions for genetically improving the quality of animal products,it should be considered that the heritability of the predictions and particularly their genetic correlations with measured traits are probably more important than phenotypic accuracy,the precision of the technique,and the lack of bias in it.A poor prediction performance of calibration model certainly affects phenotypic accuracy of NIRS predictions,but not necessarily their usefulness for the prediction of the genetic component of the same traits.In this regard,both infrared spectrometers proved to be useful tools for establishing programs for indirect genetic improvement of most of the meat quality traits,that cannot directly be measured and improved at population level.

    Conclusions

    This study investigated the feasibility of selection for meat quality traits using NIRS predictions obtained from spectra acquired in the abattoir from the intact muscle surface using portable instruments.The accuracy of the predictions was good for color traits,but low for the other traits investigated.Nevertheless,the estimated genetic parameters showed that NIRS predictions of color traits,pH and PL can be used as indicator traits of the corresponding measurements for selection purposes.The results for CL were more doubtful,while the estimates for WBSF predictions were unreliable.

    The two spectrometers compared in this study yielded similar results for the prediction of color traits and their relative genetic parameters.However,the Vis-NIRS instrument produced better estimates of the genetic parameters of the other predicted meat quality traits.

    On the whole,the selection of complex traits such as those related to meat quality,which are difficult to phenotype directly,can benefit from the application of NIRS technology.

    Acknowledgements

    Not applicable.

    Authors’contributions

    Conceptualization:A.A.,A.B.,L.D.S.,A.C.and G.B.;Data curation:S.S.and A.A.;Formal analysis:S.S.,A.A.,A.B.and L.D.S.;Funding acquisition:A.A.;Investigation:S.S.,A.A.and G.B.;Methodology:A.A.,A.C.and G.B.;Project administration:A.A.;Resources:A.A.,A.B.,L.D.S.and G.B.;Supervision:A.A.and G.B.;Writing,review &editing:S.S.,A.A.,A.C.and G.B.All authors read and approved the final manuscript.

    Funding

    This research was funded by FONDAZIONE CASSA DI RISPARMIO DI CUNEO and is part of the project“QUALIPIEM -Innovative tools for the selection of meat quality in the Piemontese breed”,project number 2014/0249 coordinator Andrea Albera.

    Availability of data and materials

    The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

    Ethics approval and consent to participate

    Animal Care and Use Committee approval was not obtained for this study because the data were obtained from an existing database;the analysed records were collected after slaughtering of animals in a commercial abattoir(Operti,Centallo– CN,Italy) from April 2015 to February 2017.The authors did not have direct control over the care of the animals included in this study.

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1Associazione Nazionale Allevatori Bovini di Razza Piemontese,strada provinciale Trinita’32/A,12061 Carrù,CN,Italy.2Department of Agronomy,Food,Natural Resources,Animals and Environment (DAFNAE),University of Padova (Padua),viale dell’Università 16,35020 Legnaro,PD,Italy.3Department of Agricultural,Forest and Food Sciences,University of Torino,Via L.Da Vinci 44,10095 Grugliasco,TO,Italy.

    国产午夜精品久久久久久一区二区三区| 少妇的逼好多水| av免费在线看不卡| 日本免费a在线| 人妻少妇偷人精品九色| 中文字幕人妻熟人妻熟丝袜美| 麻豆成人av视频| 国产极品精品免费视频能看的| av免费观看日本| 中文字幕熟女人妻在线| 日韩一区二区三区影片| 内地一区二区视频在线| 99久久九九国产精品国产免费| 免费看a级黄色片| 久久久久久久久大av| 人人妻人人澡欧美一区二区| 亚洲真实伦在线观看| 午夜福利高清视频| 国产精品一区二区三区四区免费观看| 国产老妇伦熟女老妇高清| 最近2019中文字幕mv第一页| 国产精品伦人一区二区| 精品日产1卡2卡| 99在线视频只有这里精品首页| 久久久久性生活片| 国产精品一区www在线观看| 青春草亚洲视频在线观看| 又粗又硬又长又爽又黄的视频 | 亚洲精品粉嫩美女一区| 老司机福利观看| 自拍偷自拍亚洲精品老妇| 可以在线观看毛片的网站| 综合色丁香网| 色播亚洲综合网| 午夜福利视频1000在线观看| 久久久久九九精品影院| 日韩亚洲欧美综合| 午夜爱爱视频在线播放| 欧美xxxx性猛交bbbb| 久久精品国产亚洲网站| 能在线免费看毛片的网站| 熟女人妻精品中文字幕| 成熟少妇高潮喷水视频| 看免费成人av毛片| 婷婷六月久久综合丁香| 午夜免费激情av| 老女人水多毛片| 成人亚洲欧美一区二区av| 精品国产三级普通话版| 18禁裸乳无遮挡免费网站照片| 别揉我奶头 嗯啊视频| 国产精品国产高清国产av| 亚洲欧美日韩无卡精品| 中文字幕熟女人妻在线| 色视频www国产| 少妇的逼水好多| 一区二区三区高清视频在线| 狂野欧美激情性xxxx在线观看| 熟女人妻精品中文字幕| 亚洲,欧美,日韩| 深夜a级毛片| 一级毛片久久久久久久久女| 中出人妻视频一区二区| 国产精品一区二区在线观看99 | 99热精品在线国产| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 日韩一区二区三区影片| 少妇高潮的动态图| 欧美日韩精品成人综合77777| 午夜亚洲福利在线播放| 亚洲图色成人| 色哟哟哟哟哟哟| 少妇的逼好多水| 日韩 亚洲 欧美在线| 亚洲av中文字字幕乱码综合| 少妇人妻一区二区三区视频| 久久韩国三级中文字幕| 老女人水多毛片| 亚洲中文字幕一区二区三区有码在线看| 黑人高潮一二区| 在现免费观看毛片| 在线国产一区二区在线| 久久精品国产自在天天线| 国产成人a∨麻豆精品| 亚洲久久久久久中文字幕| 久久久久久久久久久免费av| 国产高清视频在线观看网站| 熟女人妻精品中文字幕| 国内精品一区二区在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜免费男女啪啪视频观看| 精品人妻视频免费看| 国产精品女同一区二区软件| 午夜福利成人在线免费观看| 国产白丝娇喘喷水9色精品| 熟女电影av网| 三级经典国产精品| 一级黄色大片毛片| 国产69精品久久久久777片| 精品久久国产蜜桃| 日本免费a在线| 99热全是精品| 韩国av在线不卡| 最后的刺客免费高清国语| 久久久久九九精品影院| 麻豆乱淫一区二区| 国产美女午夜福利| 欧美一区二区国产精品久久精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 特级一级黄色大片| 成人午夜精彩视频在线观看| 久久99蜜桃精品久久| 国产精品一及| 人妻系列 视频| 国产精品乱码一区二三区的特点| 少妇被粗大猛烈的视频| 亚洲精品影视一区二区三区av| 国产在视频线在精品| 亚洲第一区二区三区不卡| 长腿黑丝高跟| 干丝袜人妻中文字幕| 99久久九九国产精品国产免费| 99热网站在线观看| 最近2019中文字幕mv第一页| 精品久久久久久久人妻蜜臀av| 亚洲无线观看免费| 国产乱人视频| 日韩欧美三级三区| 麻豆乱淫一区二区| 国产精品久久视频播放| 在线免费十八禁| 插阴视频在线观看视频| 久99久视频精品免费| 99热精品在线国产| 欧美+亚洲+日韩+国产| 欧美性感艳星| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 有码 亚洲区| 青春草国产在线视频 | 又粗又硬又长又爽又黄的视频 | 人妻夜夜爽99麻豆av| 一级毛片电影观看 | 亚洲不卡免费看| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 日韩在线高清观看一区二区三区| 亚洲人成网站高清观看| 免费人成在线观看视频色| 色5月婷婷丁香| 国产熟女欧美一区二区| 人人妻人人澡欧美一区二区| 国产男人的电影天堂91| 日本欧美国产在线视频| 国产成人午夜福利电影在线观看| 亚洲国产精品sss在线观看| 免费人成视频x8x8入口观看| 日韩大尺度精品在线看网址| 日韩中字成人| 成人av在线播放网站| 久久久久免费精品人妻一区二区| 97超视频在线观看视频| 男人狂女人下面高潮的视频| 麻豆久久精品国产亚洲av| 中国国产av一级| 欧美另类亚洲清纯唯美| 欧美区成人在线视频| 亚洲精品粉嫩美女一区| 久久精品国产鲁丝片午夜精品| 国内精品美女久久久久久| 欧美性感艳星| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久 | 欧美成人免费av一区二区三区| 国产精品久久久久久久电影| 观看美女的网站| 99久国产av精品| 99视频精品全部免费 在线| 精品久久久久久久久av| 久久草成人影院| 免费观看在线日韩| 一区二区三区四区激情视频 | 男女边吃奶边做爰视频| 最近2019中文字幕mv第一页| 成人欧美大片| 欧美激情国产日韩精品一区| 岛国毛片在线播放| 边亲边吃奶的免费视频| 欧美精品一区二区大全| 亚洲中文字幕一区二区三区有码在线看| 亚洲av中文av极速乱| 黄片wwwwww| 国产高清有码在线观看视频| 亚洲欧美中文字幕日韩二区| 九九爱精品视频在线观看| 国产在线精品亚洲第一网站| 高清午夜精品一区二区三区 | 免费黄网站久久成人精品| 免费一级毛片在线播放高清视频| 国内精品久久久久精免费| 亚洲欧美日韩东京热| 99视频精品全部免费 在线| 久久精品国产自在天天线| 精品欧美国产一区二区三| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄 | 久久九九热精品免费| 久久精品国产亚洲av天美| 蜜桃亚洲精品一区二区三区| 日韩三级伦理在线观看| 日日啪夜夜撸| 中出人妻视频一区二区| 成人欧美大片| 欧美区成人在线视频| 亚洲精品亚洲一区二区| 久久久久久国产a免费观看| 性色avwww在线观看| 亚洲丝袜综合中文字幕| 日韩精品青青久久久久久| 18禁在线播放成人免费| 亚洲成a人片在线一区二区| 国产精品人妻久久久久久| 少妇人妻一区二区三区视频| 国产黄色小视频在线观看| 美女内射精品一级片tv| 国产高清三级在线| 女同久久另类99精品国产91| 色噜噜av男人的天堂激情| 能在线免费看毛片的网站| 国产在视频线在精品| 亚洲精品乱码久久久久久按摩| 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕| 亚洲人成网站在线观看播放| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站| 中文字幕制服av| 看片在线看免费视频| 只有这里有精品99| 男女视频在线观看网站免费| 美女被艹到高潮喷水动态| 欧美激情久久久久久爽电影| 日产精品乱码卡一卡2卡三| or卡值多少钱| 亚洲成a人片在线一区二区| 免费看日本二区| 亚洲成av人片在线播放无| 国产精品,欧美在线| 午夜老司机福利剧场| 真实男女啪啪啪动态图| 久久热精品热| 观看美女的网站| 亚洲美女搞黄在线观看| 色哟哟哟哟哟哟| 最近中文字幕高清免费大全6| 在线播放国产精品三级| ponron亚洲| 国产一区二区亚洲精品在线观看| 国产精品蜜桃在线观看 | 国产在线男女| 欧美又色又爽又黄视频| 国国产精品蜜臀av免费| 中出人妻视频一区二区| 在线国产一区二区在线| 中文字幕av成人在线电影| 成年女人永久免费观看视频| 欧美zozozo另类| 麻豆乱淫一区二区| 亚洲三级黄色毛片| 国产精品1区2区在线观看.| 国产精品久久久久久精品电影| 91麻豆精品激情在线观看国产| av视频在线观看入口| 特大巨黑吊av在线直播| 国产成人freesex在线| 一本一本综合久久| 99热6这里只有精品| 黄色欧美视频在线观看| 久久人人爽人人片av| 国模一区二区三区四区视频| 欧美潮喷喷水| 老熟妇乱子伦视频在线观看| 久久99热这里只有精品18| av在线播放精品| 国产精品美女特级片免费视频播放器| 国产 一区 欧美 日韩| 欧美日韩一区二区视频在线观看视频在线 | 伦精品一区二区三区| 18禁黄网站禁片免费观看直播| 久久6这里有精品| 国产精品一二三区在线看| 观看免费一级毛片| 日韩成人伦理影院| 大又大粗又爽又黄少妇毛片口| 少妇人妻精品综合一区二区 | 婷婷色av中文字幕| 一本精品99久久精品77| 欧美最新免费一区二区三区| 亚洲av男天堂| 久久午夜福利片| 人妻系列 视频| 欧美日本视频| 日韩一本色道免费dvd| 日日干狠狠操夜夜爽| 国产高清三级在线| 国产爱豆传媒在线观看| 女人十人毛片免费观看3o分钟| 国产探花在线观看一区二区| 欧美日韩在线观看h| 免费搜索国产男女视频| 日本三级黄在线观看| 国产精品一区二区在线观看99 | 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 日本黄大片高清| 国产成人福利小说| 美女 人体艺术 gogo| 免费无遮挡裸体视频| av又黄又爽大尺度在线免费看 | 深夜精品福利| 嫩草影院新地址| 最近手机中文字幕大全| 久久精品久久久久久久性| 99在线视频只有这里精品首页| 老司机影院成人| 国产精品野战在线观看| 在线天堂最新版资源| 99久久中文字幕三级久久日本| 在线免费观看不下载黄p国产| 久久亚洲国产成人精品v| 国产不卡一卡二| 久久韩国三级中文字幕| 欧美人与善性xxx| 最近视频中文字幕2019在线8| 成人午夜精彩视频在线观看| 国产私拍福利视频在线观看| 亚洲国产精品国产精品| 精华霜和精华液先用哪个| 国模一区二区三区四区视频| 十八禁国产超污无遮挡网站| 亚洲精品久久国产高清桃花| 欧美色欧美亚洲另类二区| 亚洲欧美成人精品一区二区| 小说图片视频综合网站| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲网站| 婷婷精品国产亚洲av| 波多野结衣高清作品| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| 亚洲av成人av| 午夜精品一区二区三区免费看| 午夜a级毛片| 成人亚洲精品av一区二区| 亚洲精品成人久久久久久| 嘟嘟电影网在线观看| 久久99蜜桃精品久久| 青春草国产在线视频 | 最近视频中文字幕2019在线8| 国产午夜精品一二区理论片| 亚洲av免费在线观看| 最近最新中文字幕大全电影3| 在线a可以看的网站| 神马国产精品三级电影在线观看| 国产一区二区三区在线臀色熟女| 熟女人妻精品中文字幕| 亚洲经典国产精华液单| 久久久久久久午夜电影| 久久草成人影院| 国产午夜福利久久久久久| 91久久精品国产一区二区三区| 国产单亲对白刺激| 欧美日韩一区二区视频在线观看视频在线 | av在线亚洲专区| eeuss影院久久| 1024手机看黄色片| 老熟妇乱子伦视频在线观看| 国产精品免费一区二区三区在线| 18禁在线无遮挡免费观看视频| 日韩成人伦理影院| 日韩欧美一区二区三区在线观看| 美女高潮的动态| 亚洲经典国产精华液单| 噜噜噜噜噜久久久久久91| 男女那种视频在线观看| 日本免费一区二区三区高清不卡| 一区二区三区四区激情视频 | 国产av不卡久久| 亚洲国产精品久久男人天堂| 2021天堂中文幕一二区在线观| 亚洲高清免费不卡视频| 国产精品99久久久久久久久| 亚洲国产欧美在线一区| 久久久国产成人精品二区| 久久久久久九九精品二区国产| 能在线免费看毛片的网站| 国产精品99久久久久久久久| 色哟哟·www| 老女人水多毛片| 少妇丰满av| 国产精品一二三区在线看| 国产免费一级a男人的天堂| 99久久无色码亚洲精品果冻| av在线亚洲专区| 黑人高潮一二区| 别揉我奶头 嗯啊视频| 精品人妻视频免费看| 最近2019中文字幕mv第一页| 久久精品91蜜桃| 天天躁日日操中文字幕| 国产伦理片在线播放av一区 | 国产日本99.免费观看| 国产精品一区二区性色av| 美女内射精品一级片tv| 搞女人的毛片| 淫秽高清视频在线观看| 精品久久久噜噜| 国产蜜桃级精品一区二区三区| 免费看av在线观看网站| 边亲边吃奶的免费视频| 给我免费播放毛片高清在线观看| 国产伦一二天堂av在线观看| 亚洲国产欧美人成| 午夜久久久久精精品| 国产高清激情床上av| 中文字幕制服av| 日韩av不卡免费在线播放| 欧美色视频一区免费| 天堂中文最新版在线下载 | 卡戴珊不雅视频在线播放| 欧美激情国产日韩精品一区| 国产高清视频在线观看网站| 日本三级黄在线观看| 亚洲欧美日韩无卡精品| 日本撒尿小便嘘嘘汇集6| 女同久久另类99精品国产91| 免费黄网站久久成人精品| 亚洲第一区二区三区不卡| 亚洲高清免费不卡视频| 看十八女毛片水多多多| 69人妻影院| 久久久久国产网址| 青春草亚洲视频在线观看| 成人欧美大片| 亚洲精品色激情综合| 高清在线视频一区二区三区 | 欧美三级亚洲精品| 亚洲精品影视一区二区三区av| 国产色婷婷99| 国产精品久久久久久精品电影| 国产精品久久久久久久久免| 国产精品综合久久久久久久免费| av卡一久久| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 亚洲av不卡在线观看| 我的老师免费观看完整版| 亚洲欧美清纯卡通| 在线观看66精品国产| 国产精品女同一区二区软件| 成人一区二区视频在线观看| 亚洲一级一片aⅴ在线观看| 精品久久久久久成人av| 91久久精品电影网| 久久久久网色| av在线蜜桃| 日韩成人av中文字幕在线观看| 人人妻人人看人人澡| 中文精品一卡2卡3卡4更新| 级片在线观看| 91aial.com中文字幕在线观看| 可以在线观看的亚洲视频| 日本黄大片高清| 在现免费观看毛片| 边亲边吃奶的免费视频| 国产淫片久久久久久久久| 美女脱内裤让男人舔精品视频 | 成人美女网站在线观看视频| 热99在线观看视频| 午夜激情欧美在线| 国产 一区精品| 久久久久网色| 国产午夜福利久久久久久| 岛国在线免费视频观看| 在线国产一区二区在线| 欧美精品国产亚洲| 极品教师在线视频| 日本-黄色视频高清免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美成人精品一区二区| 国产一区二区亚洲精品在线观看| 亚洲av成人精品一区久久| 亚洲av免费在线观看| 国产成人a区在线观看| 国产精品一及| 国产久久久一区二区三区| www日本黄色视频网| 美女高潮的动态| 你懂的网址亚洲精品在线观看 | 99久国产av精品国产电影| 色5月婷婷丁香| 九九热线精品视视频播放| 赤兔流量卡办理| 亚洲精品粉嫩美女一区| 成人欧美大片| 看非洲黑人一级黄片| 欧美成人免费av一区二区三区| 国产极品天堂在线| 在现免费观看毛片| 身体一侧抽搐| 亚洲欧美中文字幕日韩二区| 悠悠久久av| 亚洲人与动物交配视频| 在线天堂最新版资源| 麻豆成人午夜福利视频| 在线观看免费视频日本深夜| 国模一区二区三区四区视频| 插逼视频在线观看| 99国产精品一区二区蜜桃av| 麻豆成人av视频| 成人午夜高清在线视频| 欧美不卡视频在线免费观看| 免费搜索国产男女视频| 波野结衣二区三区在线| 好男人视频免费观看在线| 国产伦在线观看视频一区| 国产视频首页在线观看| 乱码一卡2卡4卡精品| 午夜精品一区二区三区免费看| 深夜精品福利| 天堂网av新在线| 亚州av有码| 桃色一区二区三区在线观看| 亚洲精品自拍成人| 亚洲av免费高清在线观看| 永久网站在线| 亚洲人与动物交配视频| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩高清专用| 亚洲欧美精品专区久久| 国产乱人偷精品视频| 久久久精品94久久精品| 国产黄片美女视频| 少妇人妻精品综合一区二区 | 成人二区视频| 国产高清视频在线观看网站| 国产探花极品一区二区| 在线免费十八禁| 3wmmmm亚洲av在线观看| 成人一区二区视频在线观看| 男女视频在线观看网站免费| 国产精品久久久久久亚洲av鲁大| 国内久久婷婷六月综合欲色啪| 国产一区二区三区在线臀色熟女| 波野结衣二区三区在线| 麻豆国产av国片精品| 欧美又色又爽又黄视频| 观看免费一级毛片| 亚洲av第一区精品v没综合| 日本黄大片高清| 午夜福利成人在线免费观看| 一级毛片电影观看 | 神马国产精品三级电影在线观看| 高清毛片免费看| 国产一区二区在线av高清观看| 国语自产精品视频在线第100页| 99热6这里只有精品| 国产私拍福利视频在线观看| 国产探花极品一区二区| 三级经典国产精品| 亚洲av免费高清在线观看| 婷婷色综合大香蕉| 精品久久久久久久人妻蜜臀av| 国产精品久久视频播放| 久久99蜜桃精品久久| 欧美一级a爱片免费观看看| 欧美激情在线99| 午夜福利高清视频| 尤物成人国产欧美一区二区三区| 国产精品女同一区二区软件| 精品一区二区三区人妻视频| 精品久久久久久成人av| 日韩强制内射视频| 亚洲欧美精品自产自拍| 亚洲久久久久久中文字幕| or卡值多少钱| 天堂中文最新版在线下载 | 国产精品无大码| 99精品在免费线老司机午夜| 六月丁香七月| 中文欧美无线码| 99精品在免费线老司机午夜| 午夜福利在线观看吧| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 99九九线精品视频在线观看视频| 身体一侧抽搐| 日本在线视频免费播放| 青青草视频在线视频观看| 亚洲成人久久性| 99热这里只有精品一区| 成人漫画全彩无遮挡| 成人毛片60女人毛片免费| 别揉我奶头 嗯啊视频| 在线观看美女被高潮喷水网站| 国产精品一区二区性色av| 日本三级黄在线观看| 少妇裸体淫交视频免费看高清| 欧美一区二区精品小视频在线| 午夜久久久久精精品| av在线播放精品| 亚洲精品久久久久久婷婷小说 | 国产伦在线观看视频一区|