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

    Simultaneously quantifying hundreds of acylcarnitines in multiple biological matrices within ten minutes using ultrahigh-performance liquid-chromatography and tandem mass spectrometry

    2024-03-21 05:51:14JingxinZhngQinshengChenLinglongZhngBiruShiMenYuQingxiHungHuiruTng
    Journal of Pharmaceutical Analysis 2024年1期

    Jingxin Zhng , Qinsheng Chen , Linglong Zhng , Biru Shi , Men Yu ,Qingxi Hung ,**, Huiru Tng ,*

    a State Key Laboratory of Genetic Engineering, School of Life Sciences, Human Phenome Institute, Zhangjiang Fudan International Innovation Center,Metabonomics and Systems Biology Laboratory at Shanghai International Centre for Molecular Phenomics, Zhongshan Hospital, Fudan University,Shanghai, 200032, China

    b Wuhan Laboratory for Shanghai Metabolome Institute (SMI) Ltd., Wuhan, 430000, China

    Keywords:Acylcarnitine UPLC-MS/MS Quantitative structure-retention relationship Molecular phenotype

    A B S T R A C T Acylcarnitines are metabolic intermediates of fatty acids and branched-chain amino acids having vital biofunctions and pathophysiological significances.Here, we developed a high-throughput method for quantifying hundreds of acylcarnitines in one run using ultrahigh performance liquid chromatography and tandem mass spectrometry (UPLC-MS/MS).This enabled simultaneous quantification of 1136 acylcarnitines (C0-C26) within 10-min with good sensitivity (limit of detection <0.7 fmol), linearity (correlation coefficient >0.992), accuracy (relative error <20%), precision (coefficient of variation (CV),CV <15%),stability(CV <15%),and inter-technician consistency(CV <20%,n=6).We also established a quantitative structure-retention relationship(goodness of fit >0.998)for predicting retention time(tR)of acylcarnitines with no standards and built a database of their multiple reaction monitoring parameters(tR,ion-pairs,and collision energy).Furthermore,we quantified 514 acylcarnitines in human plasma and urine, mouse kidney, liver, heart, lung, and muscle.This provides a rapid method for quantifying acylcarnitines in multiple biological matrices.

    1.Introduction

    Acylcarnitines are esters of carnitine (3-hydroxy-4-trimethylaminobutyric acid) derived from fatty-acids and branched-chain amino acids in mitochondria, peroxisome, and endoplasmic reticulum [1-10].Although several thousand acylcarnitine species were predicted to be possibly present in the human samples [9], only a few hundred acylcarnitines were collectively reported with their acyl-chains containing 2-28 carbons(C2-C28)[11,12].Based on the acyl-chain length,they can be categorized as short-chain (C2-5), medium-chain (C6-12), longchain (C13-20), and very-long chain (>C20) acylcarnitines [8-10](Fig.1).The acyl-chains may also contain hydroxyl,carboxylic,and side-chain methyl groups so that these subgroups can be termed as hydroxylated, carboxylated, and branched-chain acylcarnitines,respectively [9,10].More recently, some acylcarnitines conjugated with numerous hydrophilic compounds including amino acids,glucuronic acid,and sulfoacid(Fig.1)were discovered though only in urine so far which were probably derived from phase II biotransformation [12].

    These acylcarnitines play many vital roles including metabolic homeostasis of lipids and amino acids, fatty acid transportation,and host-microbial symbiosis [3-10].Numerous acylcarnitines have been reported as biomarkers for type-2 diabetes, cardiovascular diseases, Alzheimer's disease, non-alcoholic fatty liver disease, tumors, and inborn errors of metabolism [13-21].However,the functions of branched-chain, hydroxylated, carboxylated, and conjugated acylcarnitines are far from fully understood.Conceivably, simultaneous quantification of all these acylcarnitines is the prerequisite to effectively investigate their diverse functions thus reveal acylcarnitine-related molecular phenotypes of biological systems.Unfortunately, such analysis is nontrivial with huge diversity for the structural and physiochemical properties of acylcarnitines.Concurrent high-coverage and high-throughput quantification of acylcarnitines is still needed especially for large cohort population studies.

    Fig.1.Structure of carnitines detected in biological samples.(A) Carnitine (CA).(B)Typical acylcarnitines whose acyl-group R (C2-C26) may contain straight chains and branched-chains, double bonds, hydroxyl (OH), and carboxyl (DC) groups.(C) Some typical conjugated acylcarnitines.

    Ultrahigh performance liquid chromatography coupled with mass spectrometry (UPLC-MS) has shown great potentials for simultaneous detection and quantification of acylcarnitines[11,12,22-24], whilst flow injection MS methods [25] cannot distinguish isomers.For example, 733 acylcarnitines in the pooled samples from rat liver tissue, human plasma, and urine were reportedly detected with UPLC-coupled high-resolution mass spectrometry(HRMS)although the method required more than 10 runs (with 19 min each) [11].A recent HRMS study also reported coverage of 586 acylcarnitines and 125 conjugated acylcarnitines from human urine[12]with 108 individual runs(22 min each)but without quantitative data.For simultaneous quantifications,HRMS is limited by its scanning-speed and dynamic range hence accuracy for low abundant acylcarnitines such as those having polyunsaturated, hydroxylated, and carboxylated acyl-chains [24].UPLC-MS/MS in scheduled multiple reaction monitoring (MRM)mode should be the best choice.However, the reported UPLC-MS/MS methods either had thermal instability problems or limited coverage, throughput or isomer differentiability [26-30].Even with the combination of reversed- and normal-phase liquid chromatography (LC), such LC-MS/MS method could only quantify 117 acylcarnitines in biofluids with a 26-min run[29].All the reported LC-MS/MS methods covered only a few hydroxylated and carboxylated acylcarnitines with no conjugated acylcarnitines [9].

    One of the bottlenecks for high-coverage LC-MS/MS quantification of acylcarnitines is the lack of standards to reliably obtain the retention time(tR)of all possible acylcarnitines for scheduled MRM analysis.For hundreds of acylcarnitines reportedly detected, only dozens of them have commercial standards.Under such circumstances, reliable modelling for quantitative structure-retention relationship (QSRR) becomes helpful to predict their tR[31,32] for setting up the MRM-needed time-windows and to assist identification.Previous studies have already shown potential linear relationship between tRand some structural descriptors for acylcarnitines [11,30,33].Unfortunately, these required massive calculation resources and gave no information for the tRdependence on functional groups.Understandably, the retention behaviors of acylcarnitines collectively depend on structural characteristics of their acyl-chains including the number of carbons, double-bonds, side-chains, hydroxyl and carboxyl groups,and their positions.So far, however, the best such study has considered only the number of carbon atoms and double bonds in acyl-chains with several linear models constructed for different classes of acylcarnitines[34].A unified QSRR model for demanded structural features has not been reported so far, to the best of our knowledge.

    Here, we developed an MRM-based sensitive UPLC-MS/MS method for simultaneous quantification of all identifiable acylcarnitines in biofluids and tissues with a single 10-min run.We also built a QSRR model for predicting acylcarnitine tRvalues with their seven structural features and constructed a database of MRM parameters (tR, MRM ion-pairs, and collision energy) for 1136 acylcarnitines detected in literatures and with UPLC-based quadrupole time-of-flight MS (UPLC-QTOFMS) in this study.We further quantified acylcarnitines in two typical human biofluids (plasma and urine) and five typical tissues of animal model to reveal their molecular phenotypic characteristics.

    2.Materials and methods

    2.1.Chemicals

    HPLC grade methanol (MeOH), isopropyl alcohol (IPA), and dichloromethane were obtained from Merck(Darmstadt,Germany)whereas acetonitrile (ACN), formic acid, and acetic acid were purchased from Sigma-Aldrich (St.Louis, MO, USA).Ultrapure water(H2O)was obtained from a Milli-Q system(Millipore,Billerica,MA,USA).49 acylcarnitine standards and 26 deuterated ones as internal standards (IS) were purchased from commercial sources with detailed information listed in Table S1.

    2.2.Collection of biological samples

    Human plasma and urine were obtained from Chinese adult volunteers recruited for the Human Phenome Project which was approved by the Ethic Committee of Fudan University (Approval No.:FE21087)with informed consent from all participants.Human lung cancer cell lines (A549, 16HBE) were obtained from China Center for Type Culture Collection.The New Zealand rabbit liver and five tissue samples of C57BL/6 mice were collected according to the National Guidelines for Experimental Animal Welfare(Ministry of Science and Technology of China, 2006).All biological samples were snap-frozen with liquid-nitrogen and stored at -80°C until further analysis.

    2.3.Preparation of stock and working solutions

    Each acylcarnitine standard was dissolved in MeOH or CH2Cl2:MeOH (3:2, V/V) to prepared individual stock solution(1-15 mM).These solutions were mixed, dried with nitrogen-gas,and reconstituted into aqueous MeOH (90%, V/V) to obtain the stock solution of acylcarnitine mixture.This mixture solution was then sequentially diluted to obtain 13 solutions for constructing calibration curves (Levels 1-13, Tables S2-S4).The mixture of 26 deuterated acylcarnitines in 90% (V/V) aqueous MeOH was used as IS.

    2.4.Preparation of biological samples for analysis of acylcarnitines

    For UPLC-QTOFMS analysis, acylcarnitines were extracted from multiple biological matrices including rabbit liver tissue, human plasma,urine,and cells as described in the Supplementary data.For quantification of acylcarnitines in biofluids and tissues, analytes were extracted with an optimized procedure.Briefly,20 μL of each human plasma or urine was mixed with 500 μL pre-cooled IPA containing 0.5%(V/V)acetic acid,20 μL mixed IS solution,and 30 μL water, and sonicated in an ice-bath for 5 min.After centrifugation(4°C,16,000 g)for 10 min,520 μL of supernatant was obtained.The above procedure was repeated once more and resultant 2 supernatants were combined and dried with nitrogen-gas.Each tissue sample (10 mg) was added with 400 μL pre-cooled IPA containing 0.5% (V/V) acetic acid, 20 μL mixed IS, and 30 μL water, and then homogenized with a tissuelyser (20 Hz) for 3 min to obtain supernatant after 10-min centrifugation (4°C, 16,000 g).This extraction was repeated two more times and three resultant supernatants were pooled and dried.These dried extracts were individually reconstituted into 90% (V/V) MeOH (40 μL) for UPLCMS/MS analysis.

    2.5.UPLC-MS analysis

    A ZenoTOF 7600 and TripleTOF 5600 plus (SCIEX, Chromos,Singapore) coupled to a Shimazu UPLC system were used in the information-dependent acquisition (IDA) mode (Supplementary data) to detect acylcarnitines and their fragment ions in multiple biological samples.A QTRAP 6500 Plus (SCIEX) coupled to two Shimadzu Nexera UPLC systems (Kyoto, Japan) was used for simultaneous quantification of acylcarnitines with scheduled MRM.Data were acquired and processed with Analyst and OS (v1.7,SCIEX).After comprehensive optimization of various C18columns and chromatographic conditions, the best performer Agilent ZORBAX Eclipse Plus C18column (2.1 mm × 100 mm, 1.8 μm; Santa Clara,CA,USA)was chosen at 40°C with 0.5 μL of sample injection.H2O and ACN(both containing 0.1%(V/V)formic acid)were used as mobile phase A and B,respectively,with elution gradients listed in Table S5.The dual-column-switching mode (Fig.S1) with Multiplexing (MPX) software was used to double throughput.By using acylcarnitine standards, mass spectrometric parameters were optimized including ionspray voltage, curtain gas, ion source temperature, collision energy, and declustering potential(Supplementary data).Collision energy values were further adjusted to ensure simultaneous quantification of all acylcarnitines in a single run for all biological matrices (Table S6).

    2.6.Method validation

    Method validation was conducted for sensitivity, linearity, accuracy, precision, stability, and inter-technician consistency.Linearity and sensitivity were assessed using acylcarnitine standards with variable concentration and IS.The limit of detection(LOD)was obtained as the amount of analytes on column at the signal-tonoise ratio (S/N) of 3 whereas the low limit of quantification(LLOQ) was obtained from the lowest amount in calibration curve according to the U.S.Food and Drug Administration guideline[35].Accuracy, precision, and stability were evaluated with plasma,urine,and liver samples at three different concentration levels(L3,L5,and L7).Accuracy was assessed with relative error(RE)at high,medium, and low levels whilst precision was assessed with the coefficient of variation (CV) for the intra- and inter-day measurements.Stability was assessed with extracted samples stored at -4°C for 24 h and at -80°C for 72 h.Inter-technician consistency was evaluated using the same batch of plasma samples(n = 6) with the independent results from three different technicians.

    2.7.QSRR and MRM-parameter database

    The QSRR model was built with multivariate polynomial regression approach [36,37] using experimental data from 45 known acylcarnitines (Table S1) containing saturated (18), unsaturated (12), hydroxylated (9), and dicarboxylic (6) acyl-chains.LM functions in R software package (v4.1.0) was used to optimize models with tRdependence on various structural features for acyl chains with stepwise regression and full subset regression.These features included the number of carbon atoms(c),branched-chain methyl groups (b), carbon-carbon double bonds (d), hydroxyl groups(h),carboxyl groups(j),position of branched-chain(p),and double bonds(q).10-fold cross-validation was conducted to obtain the optimized model with minimum prediction error and favorable adjusted-goodness of fit (R2a, Eq.(1)).Based on the experimental and calculated tRresults, mean absolute error (MAE, Eq.(2)) and Pearson correlation coefficient (r2) were used to assess the overall performance of model.

    3.Results and discussion

    3.1.Detection of acylcarnitines in multiple biological matrices with UPLC-QTOFMS

    Although thousands of acylcarnitines have been suggested in the Human Metabolome Database,so far,only dozens of them have experimental mass spectrometric data due to limited standards available.To obtain such data,we employed UPLC-QTOFMS in fullscan and IDA modes to detect all possible acylcarnitines present in multiple biological matrices including biofluids, cells, and tissues.Chromatographic parameters were systematically optimized including column, temperature, elution gradients, flow-rate, and injection volume to achieve maximum separation (data not shown).Since carnitine moiety has a quaternary ammonium cation,(CH3)3N+R,with permanent charge, acylcarnitines are expected to have some unique mass spectrometric features.In the negative ion mode, all acylcarnitines gave parent ions as [M+HCOO--H]-,fragment ions for carnitine moiety (m/z 188.0934 and 144.1042)and [FA-H]-from fatty-acid chains (Fig.S2).In the positive ion mode,in contrast,much more intense signals including parent ions[M]+due to the carnitine quaternary ammonium cations,[M-59.0735]+due to neutral loss of N(CH3)3, and three fragment ions from carnitine moiety (m/z 144.1019, 85.0284, and 60.0808)were observable for all acylcarnitines(Figs.2A and S2)as reported[38-40].All conjugated acylcarnitines gave parent ions [M]+,[M-conjugator-59.0735]+,and three fragment ions from carnitine moiety (Figs.2B and S2) as reported [12,38].Therefore, acylcarnitines and conjugated acylcarnitines were both readily identifiable collectively with their aforementioned five ions, respectively.Acylcarnitines gave the most intense fragment ion m/z 85.0284,thus MRM ion-pair [M]+/85.0284 was useable for quantifying acylcarnitines(Fig.S3).In contrast,[M-conjugator-59.0735]+was a unique intense fragment ion for conjugated acylcarnitines(Fig.2B),thus[M]+/[M-conjugator-59.0735]+was employable for their quantification.Typical ions were also recorded for 26 deuterated acylcarnitines, which were commercially available and used as internal standards for quantification.

    By concurrently extracting the above ions, we identified 471 acylcarnitines including 395 acylcarnitines and 76 conjugated ones from above biological matrices (rabbit liver tissue, human plasma,urine, and cells) and their pooled samples (Fig.S4).After taking structural isomers having different retention time into consideration,the coverage of acylcarnitines having different retention time was expanded to 1136 species (Table S7) including 1124 possible species derived from these collectively reported in literatures[11,12,29,41,42].Unambiguous identification was achieved for acylcarnitine isomers having commercial standards (Tables S7 and S8)such as butyryl-carnitine(C4:0)and isobutyry-carnitine(2-M-C3:0),succinyl-carnitine (C4-DC), and methylmalonyl-carnitine (2-M-C3-DC) as well as adipoyl-carnitine (C6-DC) and 3-methylglutarylcarnitine (3-M-C5-DC).For acylcarnitines with saturated fivecarbon chains (C12H23NO4), for instance, four isomers were detected with clearly different retention times.By using the commercially available standards,three of them were unambiguously assigned as valerylcarnitine (C5:0), 2-methylbutyryl-carnitine (2-M-C4:0), and isovalerylcarnitine (3-M-C4:0), respectively; pivaloylcarntine (2-2M-C3:0)became the only possibility for the remaining isomer and was further confirmed using its standard.Furthermore,the isomers without commercial standards were only considered as present here when detected with clearly different retention times(Table S8).They can be unambiguously identified with extensive syntheses of their standards, which is beyond the scope of this study.Nevertheless,we identified these isomers as analyte-a,analyte-b,etc.here to ensure their separate identities.For instance, five isomers were detected for acylcarnitines with saturated six-carbon chains giving same precursor and fragment ions in both Q-TOFMS and MS/MS.Hexanoyl-carnitine was unambiguously assigned (C6:0) with its standard and the other four were tentatively assigned to be C6:0-a,C6:0-b, C6:0-c, and C6:0-d (Table S8) respectively.It is conceivable that the identification of isomers can also be assisted to some extent by reliable prediction of their tR.

    Fig.2.Fragment assignments for (A) decanoylcarnitine (C10:0) and (B) N-acetylcysteine-conjugated-hydroxydecenoylcarnitine (N-acetylcys-C10:1-OH-CA) in the positive ion mode.* Fragment ions for carnitine moiety.Conj: conjugator.

    3.2.QSRR and MRM parameter database

    The parameters for tRand precursor-product ion pairs for MRM transitions (MRM ion-pairs) are necessary for simultaneous quantification of all acylcarnitines in the scheduled MRM mode using UPLC-MS/MS [43].However, only less than 50 acylcarnitines have standards available currently.Nevertheless, our tRresults of these available acylcarnitines clearly showed some dependence on their structural features(Fig.S5A)as indicated previously[44,45].These structural features included the number of carbons, branchedchains, double-bonds, hydroxyl, carboxylic groups, and positions in their acyl-chains.With the same number of acyl-chain carbons,the presence of these groups resulted in earlier elution.Amongst four isomers of acylcarnitines with five acyl-chain carbons, for example, C5:0 had the greatest tRfollowed with 3-M-C4:0, 2-MC4:0, and then 2-2M-C3:0.

    By using 45 acylcarnitine standards (Table S1) excluding carnitine and acetylcarnitine for their close-to-void elution, a unified QSRR model was built for tRvalues (i.e., the model-calculated retention times (tRC)) obtained from our chromatographic method(including column and elution schemes)as a function of acyl-chain structural characteristics (Eq.(3)) with the fitted k0-k10data obtained (i.e., k0= -0.252, k1= 8.54 × 10-4, k2= -4.82 × 10-2,k3= 1.06, k4= -0.487, k5= -0.172, k6= -1.65, k7= -1.55,k8=0.152,k9=-1.59×10-2,and k10=7.20×10-2).In the model,variables c, b, d, h, and j denote the number of acyl-chain carbon,branched-chain methyl group, carbon-carbon double-bond, hydroxyl, and carboxylic groups, respectively, whilst p and q denote the positions of the branched-chains and double-bonds.

    Fig.3.Comparison of the model-calculated retention times (tRC) and experimental results (tRE) for (A) 45 acylcarnitine standards and (B) acylcarnitines detected in biological samples.

    Such a unified model could predict tRvalues for all different classes of acylcarnitines with good performance (R2a>0.998) for chain length of C3-C26,which were only describable with multiple linear models previously [11,30,34].The tRCand experimental results (tRE) showed excellent linear correlation (r2, 0.998) and accuracy (ΔtR<0.20 min and MAE <0.07 min) for 45 acylcarnitine standards (Fig.3A).Such tRdifferences between the experimental and predicted values were 0.46 and 0.01 min for 2-2M-C3:0 and C4-DC,respectively,when using the present chromatographic method.

    The model also showed good tRCand tREcorrelation(r2~0.984)and prediction accuracy(MAE <0.21 min and ΔtR<0.9 min)for all 211 acylcarnitines detected here from the pooled sample having variable acyl-chain carbons, branch-chains, double-bonds, or hydroxyl or dicarboxylic groups (Fig.3B).Even with concurrent presence of these structure features,the model showed reasonable tRCand tREcorrelation (r2~0.956) and accuracy (MAE <0.25 min,ΔtR<0.9 min) for 97% acylcarnitines identified from the pooled sample(Fig.S5B).By setting the tRrange of 2 min for MRM window,even the acylcarnitine isomers without commercial standards can be reliably covered in UPLC-MS/MS quantification.Taking all this together, we built a database of MRM parameters (tR, MRM ionpairs, and collision energy) to cover 1136 acylcarnitine species with 398 formulae including acylcarnitines, conjugated ones and their isomers (Table S7).Among them, hydroxylated, unsaturated,dicarboxylated,conjugated,and saturated acylcarnitines accounted for 31.6%, 24.9%,18.5%,11.5%, and 7.9%, respectively; acylcarnitines with dicarboxylated chains and hydroxyl group accounted for 5.5%.

    3.3.Development of quantitative method for acylcarnitines using UPLC-MS/MS

    Based on the above investigations, we established an efficient method for simultaneous quantification of all possible carnitines within 10 min using UPLC-MS/MS in MRM mode.With all standards available commercially, the method showed good separation powers for acylcarnitine isomers with acyl-chains containing C0-C26(Fig.4),which was exemplified with resolution of 1.14-1.24 for acylcarnitine isomers with C5:0 (peaks 16-18).All deuterated standards with N-CD3(peaks 50-53 and 55-75) had similar retention times to the corresponding acylcarnitines(Fig.4).This is in contrast to methyl-D3-malonylcarnitine(peak 54)and suggests the N-CD3-labeled acylcarnitines as good IS for quantifications.

    Fig.4.Ultrahigh performance liquid chromatography and tandem mass spectrometry (UPLC-MS/MS) chromatogram (up: 0-6.0 min; bottom: 6.0-10.0 min) for 49 acylcarnitine standards (solid line) and 26 deuterated internal standards (dotted line) with keys listed in Table S1.

    Fragment ion m/z 85.0284 reached maximum with collision energy above 30 consistently for short chain(C4:0),medium chain(C10:0), long chain (C16:0), and very-long chain (C22:0) acylcarnitines hence being employable as a quantitative ion(Fig.S3).Other MS parameters were further optimized using all acylcarnitine standards available(Fig.S6),which were used with the deuterated standards to generate calibration curves.For these without authentic standards, such MS parameters for analytes in the same subclass with the closest retention time were used and the results were considered as relative quantification.For instance, MS parameters for C10:0 were used for quantification of C9:0, whose standard was not available.

    The method was further validated in terms of linearity, sensitivity, accuracy, precision, stability, and inter-technician consistency.For all acylcarnitines having authentic standards, excellent linearity (r2>0.992) was observed with concentration over 3-4 orders of magnitude and the LOD below 0.7 fmol on column(Table S9), with the exception for only carnitine (<3.5 fmol).The method had good accuracy (RE <20%), intra- and inter-day precision (CV <15%), as well as stability (CV <15%) at 4°C (24 h)and-80°C(72 h)for all spiked standards in plasma,liver,and urine samples (Figs.S7-S9).Acceptable inter-technician consistency(CV <20%) was further confirmed from independent experiments conducted by three different technicians (Fig.S10).With dualcolumn switching (MPX mode), finally, the method was capable of quantifying 1136 carnitines within 10 min;both throughput and coverage were better than the existing methods (Table S10)[11,12,26-30].

    3.4.Quantification of acylcarnitines in multiple biological matrices to characterize their molecular phenotypes

    By using this method, we managed to quantify 514 acylcarnitines collectively in two typical human biofluids(plasma and urine)and five mouse tissues (kidney, liver, heart, lung, and muscle) to characterize their metabolome-based molecular phenotypes(Fig.5 and Table S8).The results showed carnitine itself as the most abundant in all samples.However, clear acylcarnitine phenotypic differences were observable between biofluids and tissues, between two biofluids and between different tissues from the same animal(Fig.5 and Table S8).Human urine samples contained more acylcarnitine species (421) than human plasma (149) and mouse tissues including kidney (211), liver (195), heart (189), lung (188),and muscle (169), especially conjugated acylcarnitines and these with hydroxyl and carboxylic groups (Fig.5 and Table S8).Noticeably, some carnitine species were only detected in urine including 3-dehydrocarnitine, pivaloylcarnitine, and isomers of C6:0, C8:0,C10:0,C11:0,C12:0,C7:1,C9:1,C12:1,C8:2,C10:2,C11:2,C12:2,and C14:2.Although some of these urinary acylcarnitines (e.g., Phe-C4:1, Phe-C4:0, and odd-carbon chain acylcarnitines) were probably from diets and gut microbiota[46],their sources and functions clearly warrant further investigations.

    Fig.5.Concentration (c) in multiple biological matrices for acylcarnitine species and different types of acylcarnitines with saturated (sAC), unsaturated (uAC), hydroxylated (ACOH), carboxylated (AC-DC) chains, and conjugated acylcarnitine (CAC) with bold line (sum) and thin lines indicating total concentration of each class and each species.

    Additionally,acylcarnitines had a huge concentration range from 0.31 nM(butenoylcarnitine)to 58.73 μM(carnitine)in human plasma and even greater range in other samples (Fig.5 and Table S8).The abundances of acylcarnitines also differed greatly in two human biofluids(urine and plasma).For example, the concentration of isobutyrylcanitine (2-M-C3:0) in urine was 164 times higher than in plasma.In contrast, butyrylcarnitine (C4:0) and valerylcarnitines were not detected in urine but abundant in human plasma and all mouse tissues.The acylcarnitine compositions further showed obvious diversities for these mice tissues in terms of different acylcarnitines with the levelsranging from 0.14 pmol/mg(C9:0 and C11:0)to 1.66 μmol/mg (carnitine).Kidney had more acylcarnitine species especially conjugated ones with unsaturated and hydroxylated acylchains than the other tissues.In heart and lung tissues,the levels of carnitine (1.3-1.6 μmol/mg) and γ-butyrobetaine (21-24 nmol/mg)were about twice as high as in the other three tissues, whilst acetylcarnitine levels in heart,lung,and muscle tissues(0.16-0.29 μmol/mg) were three times higher than in the other two tissues(0.04-0.06 μmol/mg).These acylcarnitines are broadly attributable to oxidation of fatty acids and branched-chain amino acids in mitochondria and peroxisomes, although the functions of all these detected acylcarnitines warrant further investigations.

    Numerous acylcarnitine isomers were concurrently detected in this study although they could be derived from completely different metabolic pathways.Amongst isomers for acylcarnitines with five acyl-carbons, for example, 2-methylbutyrylcarnitine and isovalerylcarnitine are from isoleucine (Ile) and leucine (Leu) metabolism,respectively,whereas pivaloylcarnitine and valerylcarnitine are adducts of exogenous metabolites pivalic acid and valeric acid,respectively, originated from antibodies [47], microbes [48], and plants [49].Interestingly, both 2-methylbutyrylcarnitine and isovalerylcarnitine from two essential amino acids (Ile and Leu) had low levels in all mouse tissues.Hydroxyisobutyrylcarnitine and 3-hydroxybutyrylcarnitine are two isomers but derived from valine and β-hydroxybutyrate, respectively [4,50].The former had high concentration in heart, lung, and muscle, whereas the latter had high concentration in heart, lung, kidney, and liver (Fig.S11).Apparently, new strategies are required to ambiguously assign these acylcarnitines detected but with no standards.

    Moreover, obvious level differences were observable for acylcarnitines in human blood plasma and serum from the same blood sample (Fig.S12).The concentrations of most acylcarnitines were lower in serum than in plasma,including about 80%acylcarnitines with short chains and very long chains as well as 60%-70% acylcarnitines with medium chains and long chains.Therefore,plasma samples are recommended for metabolomic studies of acylcarnitines.

    4.Conclusions

    We developed a high-throughput UPLC-MS method and detected several hundred acylcarnitines and conjugated acylcarnitines in tested multiple biological matrices with their typical fragment ions recorded.Based on the experimental data from our method, we built a good quality model of QSSR for predicting the tRvalues of acylcarnitines with no standards.We also established MRM parameters (tR, MRM ion pairs, and collision energy) to enable simultaneous quantification of 1136 acylcarnitines using a 10-min UPLC-MS/MS method with good linearity, accuracy, precision, stability, and inter-technician consistency.To the best of our knowledge, this method has the higher coverage and throughput than reported methods.By using this method,moreover,we managed to quantify more than 500 acylcarnitine species in two human biofluids and five mouse tissues offering data for molecular phenotyping these biological matrices.This provides some essential basic biochemical data for these biological matrices and a highthroughput method to quantify acylcarnitines for large cohort studies, which are already ongoing in our lab for various diseases.

    CRediT author statement

    Jingxian Zhang:Conceptualization, Methodology, Validation,Formal analysis,Investigation,Data curation,Visualization,Writing- Original draft preparation;Qinsheng Chen:Validation, Formal analysis, Visualization, Writing - Original draft preparation;Lianglong Zhang:Validation,Formal analysis;Biru Shi:Validation;Men Yu:Conceptualization, Validation;Qingxia Huang:Conceptualization, Methodology, Validation, Formal analysis, Investigation, Data curation, Visualization, Writing - Original draft preparation;Huiru Tang:Resources, Conceptualization, Methodology, Validation, Formal analysis, Investigation, Data curation,Visualization, Writing - Reviewing and Editing, Supervision,Funding acquisition, Project administration.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    We acknowledge financial supports from the National Key R&D Program of China(Grant Nos.:2022YFC3400700,2022YFA0806400,and 2020YFE0201600),Shanghai Municipal Science and Technology Major Project(Grant No.:2017SHZDZX01),and the National Natural Science Foundation of China(Grant No.:31821002).

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jpha.2023.10.004.

    av天堂久久9| 欧美日本中文国产一区发布| 新久久久久国产一级毛片| 久久精品国产a三级三级三级| 午夜福利一区二区在线看| 色94色欧美一区二区| 亚洲五月婷婷丁香| 午夜免费鲁丝| 亚洲五月色婷婷综合| 天天影视国产精品| 欧美xxⅹ黑人| 国产男人的电影天堂91| 不卡av一区二区三区| 男女午夜视频在线观看| 国产亚洲精品久久久久5区| 精品国产一区二区久久| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 久久精品国产亚洲av香蕉五月 | 午夜免费观看性视频| 欧美激情高清一区二区三区| 在线观看免费午夜福利视频| 少妇人妻久久综合中文| 久久精品国产综合久久久| 天天躁夜夜躁狠狠躁躁| 国产成人精品在线电影| 国产区一区二久久| 欧美国产精品va在线观看不卡| av在线app专区| 亚洲av日韩精品久久久久久密| 亚洲熟女精品中文字幕| 日日摸夜夜添夜夜添小说| 999久久久精品免费观看国产| 他把我摸到了高潮在线观看 | 亚洲avbb在线观看| 亚洲精品第二区| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 正在播放国产对白刺激| 国产精品免费视频内射| 久久狼人影院| 最新的欧美精品一区二区| √禁漫天堂资源中文www| 一边摸一边抽搐一进一出视频| 国产成人系列免费观看| 黄频高清免费视频| 夫妻午夜视频| 亚洲五月色婷婷综合| 他把我摸到了高潮在线观看 | 精品第一国产精品| 欧美 日韩 精品 国产| 91精品三级在线观看| 久久精品亚洲av国产电影网| 日日爽夜夜爽网站| 极品少妇高潮喷水抽搐| 大陆偷拍与自拍| 日本vs欧美在线观看视频| 亚洲av成人一区二区三| 天天影视国产精品| 18在线观看网站| 正在播放国产对白刺激| 一区二区三区激情视频| 十八禁网站网址无遮挡| 国产极品粉嫩免费观看在线| 女人爽到高潮嗷嗷叫在线视频| av线在线观看网站| 亚洲欧美一区二区三区黑人| 欧美少妇被猛烈插入视频| 九色亚洲精品在线播放| 国产精品二区激情视频| 国产精品成人在线| 亚洲欧美日韩另类电影网站| 午夜日韩欧美国产| 精品国产国语对白av| 欧美性长视频在线观看| 精品高清国产在线一区| 97精品久久久久久久久久精品| 亚洲专区中文字幕在线| 国产精品久久久av美女十八| 在线观看免费视频网站a站| 日本黄色日本黄色录像| 肉色欧美久久久久久久蜜桃| 久久免费观看电影| 18禁裸乳无遮挡动漫免费视频| 一区二区三区乱码不卡18| 亚洲三区欧美一区| 久久久久久免费高清国产稀缺| 亚洲中文字幕日韩| 亚洲av成人不卡在线观看播放网 | 天天躁日日躁夜夜躁夜夜| 丁香六月天网| 亚洲七黄色美女视频| 国产男人的电影天堂91| 国产精品一二三区在线看| 免费在线观看完整版高清| 波多野结衣一区麻豆| 国产免费视频播放在线视频| 美女中出高潮动态图| 成年美女黄网站色视频大全免费| 夫妻午夜视频| 三上悠亚av全集在线观看| 国产成人精品在线电影| 国产亚洲午夜精品一区二区久久| 日韩中文字幕欧美一区二区| 深夜精品福利| 夜夜夜夜夜久久久久| 中文字幕精品免费在线观看视频| 80岁老熟妇乱子伦牲交| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 性色av乱码一区二区三区2| 免费久久久久久久精品成人欧美视频| 激情视频va一区二区三区| 久久久水蜜桃国产精品网| 久久久精品国产亚洲av高清涩受| 国产老妇伦熟女老妇高清| 国产精品香港三级国产av潘金莲| 成年人黄色毛片网站| 中文字幕高清在线视频| 91字幕亚洲| 亚洲,欧美精品.| 大香蕉久久网| 在线天堂中文资源库| 亚洲国产精品一区二区三区在线| 日韩人妻精品一区2区三区| 亚洲五月婷婷丁香| 国产成人免费观看mmmm| 多毛熟女@视频| 午夜福利影视在线免费观看| av视频免费观看在线观看| 免费人妻精品一区二区三区视频| 国产成人精品久久二区二区免费| 国产xxxxx性猛交| 亚洲精品自拍成人| 丝袜美足系列| 久久毛片免费看一区二区三区| 如日韩欧美国产精品一区二区三区| 真人做人爱边吃奶动态| 一区二区日韩欧美中文字幕| av福利片在线| 夜夜骑夜夜射夜夜干| av欧美777| 免费高清在线观看视频在线观看| 9191精品国产免费久久| av视频免费观看在线观看| 麻豆国产av国片精品| 91精品伊人久久大香线蕉| 在线观看人妻少妇| 青草久久国产| 我要看黄色一级片免费的| 国产精品香港三级国产av潘金莲| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 少妇被粗大的猛进出69影院| 中文字幕人妻熟女乱码| 亚洲欧美清纯卡通| 高清视频免费观看一区二区| 亚洲精品乱久久久久久| 国产无遮挡羞羞视频在线观看| 国产三级黄色录像| 男女高潮啪啪啪动态图| 亚洲av成人不卡在线观看播放网 | 亚洲精品在线美女| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| 欧美97在线视频| 亚洲人成电影观看| 日韩免费高清中文字幕av| h视频一区二区三区| 欧美精品一区二区免费开放| 欧美 亚洲 国产 日韩一| 亚洲成人国产一区在线观看| 操美女的视频在线观看| 国产精品免费视频内射| 亚洲中文av在线| 老司机午夜福利在线观看视频 | 伊人久久大香线蕉亚洲五| 精品免费久久久久久久清纯 | 2018国产大陆天天弄谢| 自线自在国产av| 满18在线观看网站| 久久国产亚洲av麻豆专区| 久久国产精品男人的天堂亚洲| 黄片大片在线免费观看| 在线观看舔阴道视频| 一区福利在线观看| 日本wwww免费看| 水蜜桃什么品种好| 久久国产精品大桥未久av| 法律面前人人平等表现在哪些方面 | 999久久久精品免费观看国产| 国产精品.久久久| 亚洲国产中文字幕在线视频| 99久久人妻综合| 啪啪无遮挡十八禁网站| 91九色精品人成在线观看| 国产主播在线观看一区二区| 国产一区二区在线观看av| 国产精品免费视频内射| 国产伦人伦偷精品视频| 国精品久久久久久国模美| 国产成人系列免费观看| 美女主播在线视频| 一边摸一边做爽爽视频免费| 高清黄色对白视频在线免费看| xxxhd国产人妻xxx| 国产精品国产av在线观看| 免费在线观看视频国产中文字幕亚洲 | 精品人妻在线不人妻| 亚洲国产精品一区二区三区在线| 18禁黄网站禁片午夜丰满| 久久ye,这里只有精品| 下体分泌物呈黄色| 波多野结衣av一区二区av| 90打野战视频偷拍视频| 午夜精品久久久久久毛片777| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品电影小说| 男女午夜视频在线观看| 精品国产乱码久久久久久男人| 亚洲欧洲日产国产| 91成人精品电影| 在线观看免费日韩欧美大片| 欧美激情久久久久久爽电影 | 操出白浆在线播放| 午夜两性在线视频| 亚洲欧美精品自产自拍| 国精品久久久久久国模美| 一区福利在线观看| 国产成人欧美在线观看 | 国产精品99久久99久久久不卡| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 搡老岳熟女国产| 国产有黄有色有爽视频| 免费不卡黄色视频| 手机成人av网站| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 搡老乐熟女国产| 少妇 在线观看| 精品国产一区二区三区四区第35| 国产伦人伦偷精品视频| 久久久久视频综合| 99国产精品99久久久久| 12—13女人毛片做爰片一| 亚洲成人国产一区在线观看| 狠狠狠狠99中文字幕| 老司机午夜福利在线观看视频 | 国产亚洲一区二区精品| 亚洲av欧美aⅴ国产| 精品高清国产在线一区| 欧美日本中文国产一区发布| 国产精品一区二区免费欧美 | 国产片内射在线| 亚洲精品在线美女| 青青草视频在线视频观看| 精品国产国语对白av| 国产又爽黄色视频| 精品免费久久久久久久清纯 | 精品一区二区三卡| 成年美女黄网站色视频大全免费| www.精华液| 18禁观看日本| 亚洲天堂av无毛| 一边摸一边抽搐一进一出视频| 在线观看免费视频网站a站| 久久精品亚洲av国产电影网| 日日摸夜夜添夜夜添小说| 亚洲五月色婷婷综合| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美精品济南到| 啦啦啦视频在线资源免费观看| 国产精品久久久久久精品古装| 两个人免费观看高清视频| 电影成人av| 国产高清videossex| 老司机在亚洲福利影院| 午夜福利视频在线观看免费| 国产av又大| 国产精品一区二区精品视频观看| 俄罗斯特黄特色一大片| 国产成人a∨麻豆精品| 欧美激情 高清一区二区三区| 可以免费在线观看a视频的电影网站| 亚洲专区中文字幕在线| 日韩有码中文字幕| 黄色 视频免费看| 久久天躁狠狠躁夜夜2o2o| 99国产精品一区二区三区| 午夜福利免费观看在线| 在线观看免费高清a一片| 久久亚洲国产成人精品v| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻,人人澡人人爽秒播| 香蕉国产在线看| 午夜激情av网站| 成年美女黄网站色视频大全免费| 又黄又粗又硬又大视频| 最近中文字幕2019免费版| 免费在线观看日本一区| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 黄片小视频在线播放| 亚洲成国产人片在线观看| 国产精品国产三级国产专区5o| 秋霞在线观看毛片| av在线老鸭窝| 最近最新中文字幕大全免费视频| www.熟女人妻精品国产| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精| 最新在线观看一区二区三区| 99国产精品一区二区三区| 国产免费av片在线观看野外av| 日韩有码中文字幕| 秋霞在线观看毛片| 国产伦人伦偷精品视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲 国产 在线| 欧美在线一区亚洲| xxxhd国产人妻xxx| 欧美av亚洲av综合av国产av| 制服诱惑二区| 日韩 亚洲 欧美在线| 亚洲精品美女久久久久99蜜臀| 亚洲欧洲日产国产| 午夜福利视频在线观看免费| 日本av免费视频播放| 亚洲自偷自拍图片 自拍| 在线天堂中文资源库| 大片电影免费在线观看免费| 婷婷成人精品国产| 久久人妻福利社区极品人妻图片| 制服人妻中文乱码| 一区福利在线观看| 午夜精品久久久久久毛片777| 水蜜桃什么品种好| 国产有黄有色有爽视频| 亚洲欧美日韩高清在线视频 | 精品少妇一区二区三区视频日本电影| 亚洲精品美女久久av网站| 国产高清videossex| 两性午夜刺激爽爽歪歪视频在线观看 | 免费少妇av软件| 久久99一区二区三区| 亚洲精品国产一区二区精华液| 日韩一区二区三区影片| 午夜福利免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 少妇精品久久久久久久| 日韩有码中文字幕| 777久久人妻少妇嫩草av网站| 捣出白浆h1v1| 一本—道久久a久久精品蜜桃钙片| 中文字幕制服av| 欧美亚洲 丝袜 人妻 在线| 色老头精品视频在线观看| 国产一区二区三区在线臀色熟女 | 日韩电影二区| 国产成人免费观看mmmm| 91字幕亚洲| 久9热在线精品视频| 日韩中文字幕欧美一区二区| 久久久久久久久久久久大奶| av不卡在线播放| 国产精品一区二区在线观看99| 色视频在线一区二区三区| 亚洲视频免费观看视频| av网站在线播放免费| 一个人免费看片子| 日本撒尿小便嘘嘘汇集6| 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻丝袜制服| 夜夜夜夜夜久久久久| 如日韩欧美国产精品一区二区三区| 久久国产亚洲av麻豆专区| 搡老岳熟女国产| 国产伦理片在线播放av一区| 国产亚洲精品久久久久5区| 国产免费视频播放在线视频| 日本精品一区二区三区蜜桃| 中国国产av一级| 另类亚洲欧美激情| 欧美一级毛片孕妇| 18禁观看日本| 亚洲一区二区三区欧美精品| 国产99久久九九免费精品| 丝瓜视频免费看黄片| 久久99一区二区三区| av天堂久久9| 在线观看一区二区三区激情| 一级,二级,三级黄色视频| 日韩 亚洲 欧美在线| 国产黄色免费在线视频| 国产精品国产三级国产专区5o| 91国产中文字幕| 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 成人国产一区最新在线观看| 亚洲成国产人片在线观看| 欧美97在线视频| svipshipincom国产片| 国产一卡二卡三卡精品| 黄色怎么调成土黄色| 19禁男女啪啪无遮挡网站| 久久 成人 亚洲| 国产成人精品在线电影| 热99国产精品久久久久久7| 人成视频在线观看免费观看| 亚洲欧美激情在线| 久久香蕉激情| 国产野战对白在线观看| av网站在线播放免费| 人妻 亚洲 视频| 成人18禁高潮啪啪吃奶动态图| 99国产精品免费福利视频| 免费在线观看日本一区| 国产成人精品在线电影| av不卡在线播放| 一级片免费观看大全| 亚洲成人免费电影在线观看| 别揉我奶头~嗯~啊~动态视频 | 午夜免费观看性视频| 国产成人欧美在线观看 | 男人舔女人的私密视频| 一本色道久久久久久精品综合| 日韩 欧美 亚洲 中文字幕| 欧美另类一区| 日本黄色日本黄色录像| 人妻人人澡人人爽人人| 欧美精品一区二区大全| 动漫黄色视频在线观看| 一区二区av电影网| 亚洲第一av免费看| 国产精品一区二区在线不卡| 日本黄色日本黄色录像| 亚洲精品国产av成人精品| 国产成人精品久久二区二区91| 一级黄色大片毛片| 亚洲人成电影免费在线| 国产精品麻豆人妻色哟哟久久| 国产精品秋霞免费鲁丝片| av在线app专区| 国产亚洲欧美精品永久| 他把我摸到了高潮在线观看 | 另类精品久久| 国产麻豆69| 久久国产精品男人的天堂亚洲| 久久久久精品国产欧美久久久 | 国产区一区二久久| 午夜福利乱码中文字幕| 天堂8中文在线网| 午夜福利免费观看在线| 亚洲午夜精品一区,二区,三区| 欧美亚洲 丝袜 人妻 在线| 在线亚洲精品国产二区图片欧美| 最新在线观看一区二区三区| 久久精品aⅴ一区二区三区四区| 久久久久久久久免费视频了| 天天添夜夜摸| 一二三四在线观看免费中文在| 青草久久国产| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 蜜桃在线观看..| 日日摸夜夜添夜夜添小说| 午夜福利在线免费观看网站| 在线 av 中文字幕| 亚洲色图 男人天堂 中文字幕| 久久人人97超碰香蕉20202| 日韩大片免费观看网站| 欧美变态另类bdsm刘玥| 国产av国产精品国产| 午夜福利在线免费观看网站| 欧美日韩精品网址| 国产欧美日韩一区二区精品| tube8黄色片| 久久ye,这里只有精品| 国产一区二区激情短视频 | 国产成人影院久久av| 午夜精品国产一区二区电影| 亚洲综合色网址| 高清视频免费观看一区二区| 女性被躁到高潮视频| 午夜福利在线免费观看网站| 亚洲精品一卡2卡三卡4卡5卡 | 午夜激情av网站| 久久久国产成人免费| tube8黄色片| 搡老乐熟女国产| 国产在线免费精品| 欧美性长视频在线观看| 欧美成狂野欧美在线观看| 成人国产一区最新在线观看| 国产成人啪精品午夜网站| 日韩三级视频一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品熟女久久久久浪| 超色免费av| 男女高潮啪啪啪动态图| 国产成人系列免费观看| 国产精品亚洲av一区麻豆| 欧美av亚洲av综合av国产av| 爱豆传媒免费全集在线观看| 丁香六月天网| 国产成人欧美在线观看 | 91精品国产国语对白视频| 日韩三级视频一区二区三区| 亚洲精品av麻豆狂野| 久久精品人人爽人人爽视色| 亚洲中文日韩欧美视频| 成在线人永久免费视频| 成年女人毛片免费观看观看9 | 国产精品 欧美亚洲| 国产免费av片在线观看野外av| av一本久久久久| 欧美黄色片欧美黄色片| 久久99热这里只频精品6学生| 宅男免费午夜| 精品一区在线观看国产| 一区二区日韩欧美中文字幕| 高清av免费在线| 亚洲国产看品久久| 国产激情久久老熟女| 两个人免费观看高清视频| 久久午夜综合久久蜜桃| 一二三四社区在线视频社区8| 精品国内亚洲2022精品成人 | 精品一区在线观看国产| 91麻豆精品激情在线观看国产 | 妹子高潮喷水视频| 日本av免费视频播放| 色婷婷久久久亚洲欧美| 久久久精品区二区三区| 色精品久久人妻99蜜桃| 久久久久网色| 69av精品久久久久久 | 欧美日韩成人在线一区二区| 水蜜桃什么品种好| av福利片在线| 午夜福利一区二区在线看| 午夜日韩欧美国产| 免费高清在线观看视频在线观看| 亚洲成人免费av在线播放| 人人澡人人妻人| 亚洲七黄色美女视频| 午夜视频精品福利| 亚洲欧洲日产国产| 国产男女超爽视频在线观看| 黄网站色视频无遮挡免费观看| 日韩大片免费观看网站| 十八禁网站网址无遮挡| 亚洲午夜精品一区,二区,三区| 一进一出抽搐动态| 国产欧美日韩一区二区三 | 欧美激情极品国产一区二区三区| 精品福利观看| 黄色片一级片一级黄色片| 亚洲成人免费av在线播放| 老司机福利观看| 新久久久久国产一级毛片| www.999成人在线观看| 亚洲情色 制服丝袜| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 亚洲精品一区蜜桃| 成年女人毛片免费观看观看9 | 免费看十八禁软件| 国产亚洲一区二区精品| 美女主播在线视频| av在线app专区| 欧美一级毛片孕妇| 精品少妇黑人巨大在线播放| 下体分泌物呈黄色| 日韩一卡2卡3卡4卡2021年| 黄片小视频在线播放| 老司机亚洲免费影院| 美女主播在线视频| 另类亚洲欧美激情| 国产精品久久久久成人av| 中文字幕人妻熟女乱码| 他把我摸到了高潮在线观看 | 日韩视频在线欧美| 一级毛片精品| 男女边摸边吃奶| 久久久国产欧美日韩av| 免费在线观看日本一区| 亚洲激情五月婷婷啪啪| 天天躁狠狠躁夜夜躁狠狠躁| 国产av国产精品国产| 免费看十八禁软件| 在线观看一区二区三区激情| 午夜福利视频精品| 最新的欧美精品一区二区| 国产在线观看jvid| 91精品伊人久久大香线蕉| 国产精品偷伦视频观看了| 免费少妇av软件| 大片免费播放器 马上看| 成年人黄色毛片网站| 久久精品aⅴ一区二区三区四区| 搡老熟女国产l中国老女人| 69精品国产乱码久久久| 精品一品国产午夜福利视频| 99久久99久久久精品蜜桃| 欧美日韩视频精品一区| 国产精品久久久人人做人人爽| 免费在线观看视频国产中文字幕亚洲 | 美女国产高潮福利片在线看| 日韩一卡2卡3卡4卡2021年| 精品人妻熟女毛片av久久网站| 欧美精品亚洲一区二区| 99热网站在线观看|