• <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成人av| 男女视频在线观看网站免费| 亚洲五月婷婷丁香| 99热这里只有精品一区 | 亚洲国产精品sss在线观看| 1024香蕉在线观看| 天天躁日日操中文字幕| 国产成人欧美在线观看| 欧美日本视频| 每晚都被弄得嗷嗷叫到高潮| 99热只有精品国产| 久久性视频一级片| 精品久久久久久久久久免费视频| 亚洲熟女毛片儿| 热99在线观看视频| 亚洲精华国产精华精| 首页视频小说图片口味搜索| 色av中文字幕| 青草久久国产| 精品乱码久久久久久99久播| 亚洲中文av在线| 波多野结衣高清作品| 国产成人福利小说| a级毛片a级免费在线| 欧美成人免费av一区二区三区| 久久精品人妻少妇| 听说在线观看完整版免费高清| 久久亚洲真实| 99riav亚洲国产免费| x7x7x7水蜜桃| 男插女下体视频免费在线播放| 麻豆成人午夜福利视频| 欧美成人性av电影在线观看| 波多野结衣高清作品| 我要搜黄色片| 国产av不卡久久| 久久精品国产99精品国产亚洲性色| 色av中文字幕| 国产精品野战在线观看| 757午夜福利合集在线观看| 欧美日韩综合久久久久久 | 91麻豆av在线| 午夜福利成人在线免费观看| 亚洲欧洲精品一区二区精品久久久| 亚洲国产欧美网| 亚洲国产精品合色在线| 免费在线观看影片大全网站| 村上凉子中文字幕在线| 亚洲美女视频黄频| 天天添夜夜摸| 成人无遮挡网站| 国产日本99.免费观看| 欧美另类亚洲清纯唯美| 九色国产91popny在线| 一区二区三区国产精品乱码| 国产精品1区2区在线观看.| 1000部很黄的大片| 一级a爱片免费观看的视频| 久久久久久久久久黄片| 香蕉丝袜av| 青草久久国产| 91麻豆av在线| 国产成人精品久久二区二区免费| 欧美日韩精品网址| 又粗又爽又猛毛片免费看| 热99在线观看视频| 成在线人永久免费视频| 在线观看午夜福利视频| 国内精品久久久久久久电影| 亚洲天堂国产精品一区在线| 国产极品精品免费视频能看的| 精品乱码久久久久久99久播| 美女 人体艺术 gogo| 国产1区2区3区精品| 国产爱豆传媒在线观看| 51午夜福利影视在线观看| 中出人妻视频一区二区| 国产视频一区二区在线看| 他把我摸到了高潮在线观看| 欧美成狂野欧美在线观看| 91老司机精品| 一区二区三区高清视频在线| 欧美日韩中文字幕国产精品一区二区三区| 久久久国产成人精品二区| 欧美绝顶高潮抽搐喷水| 草草在线视频免费看| a级毛片在线看网站| 级片在线观看| 午夜福利高清视频| 国产成人福利小说| 女人被狂操c到高潮| 少妇熟女aⅴ在线视频| 免费看a级黄色片| 91av网站免费观看| 黄色丝袜av网址大全| 亚洲人成网站在线播放欧美日韩| 美女免费视频网站| 精品久久久久久久末码| 欧美一级a爱片免费观看看| 老汉色∧v一级毛片| 久久久水蜜桃国产精品网| 欧美成狂野欧美在线观看| 美女被艹到高潮喷水动态| 欧美+亚洲+日韩+国产| 九九热线精品视视频播放| xxxwww97欧美| 一二三四在线观看免费中文在| 国产私拍福利视频在线观看| 亚洲成人久久爱视频| 欧美高清成人免费视频www| 久久国产乱子伦精品免费另类| 亚洲aⅴ乱码一区二区在线播放| 国产久久久一区二区三区| 亚洲av五月六月丁香网| 国产精品一区二区三区四区久久| 日本黄色片子视频| 老司机在亚洲福利影院| 国产高清视频在线观看网站| 色视频www国产| 国产免费男女视频| 在线永久观看黄色视频| 欧美色视频一区免费| 一区二区三区国产精品乱码| 少妇人妻一区二区三区视频| 嫩草影院入口| 真人做人爱边吃奶动态| 国产成人系列免费观看| 午夜精品在线福利| 一本综合久久免费| 色吧在线观看| 99热6这里只有精品| 一个人看视频在线观看www免费 | 嫩草影院入口| 国产欧美日韩精品一区二区| 国产一区在线观看成人免费| 好男人在线观看高清免费视频| 亚洲在线观看片| 婷婷精品国产亚洲av| 精品一区二区三区视频在线观看免费| 成人特级av手机在线观看| 毛片女人毛片| 麻豆一二三区av精品| 网址你懂的国产日韩在线| 午夜福利在线在线| 18禁国产床啪视频网站| 亚洲欧美一区二区三区黑人| 国产97色在线日韩免费| 国产熟女xx| 一进一出抽搐gif免费好疼| 婷婷精品国产亚洲av在线| 宅男免费午夜| 一个人免费在线观看的高清视频| 99riav亚洲国产免费| 亚洲国产中文字幕在线视频| 亚洲专区字幕在线| www.自偷自拍.com| 亚洲精品美女久久av网站| 免费看十八禁软件| 色av中文字幕| 一本综合久久免费| 久久久久亚洲av毛片大全| 一二三四社区在线视频社区8| 麻豆av在线久日| 精品午夜福利视频在线观看一区| 看黄色毛片网站| 欧美色视频一区免费| 最近最新中文字幕大全电影3| 国产精品国产高清国产av| 午夜福利在线在线| АⅤ资源中文在线天堂| 狠狠狠狠99中文字幕| www.熟女人妻精品国产| 午夜影院日韩av| 成人无遮挡网站| 国产伦精品一区二区三区四那| 亚洲国产精品sss在线观看| 性色avwww在线观看| 国产伦精品一区二区三区四那| 成人三级黄色视频| 久久性视频一级片| 色老头精品视频在线观看| 色综合亚洲欧美另类图片| 欧美激情在线99| 日韩高清综合在线| 国产一区二区在线观看日韩 | 中文字幕最新亚洲高清| 无限看片的www在线观看| 日韩中文字幕欧美一区二区| 欧美日本视频| 两个人的视频大全免费| 淫秽高清视频在线观看| 久久久国产成人精品二区| 中文字幕av在线有码专区| 国产成人系列免费观看| 最近最新中文字幕大全免费视频| 欧美国产日韩亚洲一区| 欧美三级亚洲精品| 亚洲中文av在线| 国产成人一区二区三区免费视频网站| 亚洲国产欧美网| 女人高潮潮喷娇喘18禁视频| 午夜亚洲福利在线播放| 精品99又大又爽又粗少妇毛片 | 悠悠久久av| 亚洲av免费在线观看| 1024手机看黄色片| 伦理电影免费视频| 不卡一级毛片| 丰满的人妻完整版| 午夜福利视频1000在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产真实乱freesex| 精品久久久久久久久久久久久| 免费在线观看视频国产中文字幕亚洲| 两个人的视频大全免费| 亚洲av电影在线进入| 香蕉国产在线看| 国产在线精品亚洲第一网站| 亚洲精品一卡2卡三卡4卡5卡| 91av网一区二区| 久久中文字幕人妻熟女| 久久九九热精品免费| 久久久久国产一级毛片高清牌| 亚洲av成人一区二区三| 午夜福利在线观看免费完整高清在 | 91在线观看av| 国产av在哪里看| 黄色视频,在线免费观看| 欧美黄色淫秽网站| 无限看片的www在线观看| 亚洲成人精品中文字幕电影| 午夜亚洲福利在线播放| 欧美日韩亚洲国产一区二区在线观看| 精品福利观看| 九九久久精品国产亚洲av麻豆 | 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区免费欧美| 手机成人av网站| 国内精品美女久久久久久| www.自偷自拍.com| 国产黄片美女视频| 午夜福利在线观看免费完整高清在 | 国产久久久一区二区三区| 亚洲熟女毛片儿| 91麻豆精品激情在线观看国产| 亚洲,欧美精品.| 99国产精品一区二区蜜桃av| 美女扒开内裤让男人捅视频| 国产精品1区2区在线观看.| 国产精品av视频在线免费观看| 99久久无色码亚洲精品果冻| 成人18禁在线播放| 精品国产亚洲在线| 久久香蕉国产精品| 久久这里只有精品19| 亚洲成av人片在线播放无| 亚洲欧美日韩无卡精品| 亚洲精品中文字幕一二三四区| 欧洲精品卡2卡3卡4卡5卡区| 小说图片视频综合网站| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 国产成人精品久久二区二区免费| 美女大奶头视频| 亚洲五月天丁香| 九色国产91popny在线| 成人av一区二区三区在线看| 90打野战视频偷拍视频| 女人被狂操c到高潮| 黄色日韩在线| 中文字幕人妻丝袜一区二区| 热99re8久久精品国产| 又黄又粗又硬又大视频| 国产私拍福利视频在线观看| 日本 欧美在线| 一级毛片高清免费大全| 亚洲av成人av| 国产av一区在线观看免费| 色在线成人网| 热99re8久久精品国产| 欧美三级亚洲精品| 丝袜人妻中文字幕| 日韩精品中文字幕看吧| 观看美女的网站| 国产激情偷乱视频一区二区| aaaaa片日本免费| 精品国产亚洲在线| 十八禁人妻一区二区| 亚洲最大成人中文| 午夜福利免费观看在线| 一级毛片高清免费大全| 日韩三级视频一区二区三区| 99久国产av精品| 岛国在线观看网站| 欧美日本亚洲视频在线播放| 日本 av在线| 窝窝影院91人妻| 好男人电影高清在线观看| av在线天堂中文字幕| 在线免费观看不下载黄p国产 | 一级毛片精品| 精品久久久久久,| 好男人在线观看高清免费视频| 国产亚洲精品综合一区在线观看| 色噜噜av男人的天堂激情| 99精品久久久久人妻精品| 少妇熟女aⅴ在线视频| 黄片小视频在线播放| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 88av欧美| 久久久久国产一级毛片高清牌| 999精品在线视频| 免费大片18禁| av黄色大香蕉| 国产爱豆传媒在线观看| 亚洲精华国产精华精| 亚洲成av人片在线播放无| 五月伊人婷婷丁香| 国产亚洲av嫩草精品影院| 麻豆一二三区av精品| 午夜a级毛片| 亚洲精品久久国产高清桃花| 欧美成人免费av一区二区三区| 老熟妇仑乱视频hdxx| 国产精品久久视频播放| 国产亚洲精品综合一区在线观看| 香蕉av资源在线| 三级毛片av免费| 九九在线视频观看精品| 久久香蕉国产精品| 欧美大码av| 美女大奶头视频| 俺也久久电影网| 欧美日本视频| 草草在线视频免费看| 一进一出抽搐动态| 国产乱人伦免费视频| 综合色av麻豆| 成熟少妇高潮喷水视频| 一区二区三区高清视频在线| 精品免费久久久久久久清纯| 精品国产乱子伦一区二区三区| 精品久久久久久,| 一本一本综合久久| 亚洲午夜精品一区,二区,三区| 国产三级黄色录像| 国产久久久一区二区三区| 午夜免费成人在线视频| 18禁美女被吸乳视频| 亚洲第一电影网av| 久久人人精品亚洲av| 一本一本综合久久| 国内精品一区二区在线观看| 国产高清视频在线观看网站| 欧美极品一区二区三区四区| ponron亚洲| 久久精品91无色码中文字幕| 欧美一级a爱片免费观看看| 2021天堂中文幕一二区在线观| 久久九九热精品免费| 成人国产综合亚洲| 最近视频中文字幕2019在线8| 九色国产91popny在线| 亚洲自拍偷在线| 18禁裸乳无遮挡免费网站照片| 少妇的逼水好多| 国产一区在线观看成人免费| 1024手机看黄色片| 麻豆成人av在线观看| 亚洲性夜色夜夜综合| 日韩欧美在线二视频| 国产精品久久久久久亚洲av鲁大| 禁无遮挡网站| 超碰成人久久| 色综合亚洲欧美另类图片| 久久草成人影院| 午夜福利高清视频| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 亚洲精品色激情综合| 99久久综合精品五月天人人| 1000部很黄的大片| 日本 欧美在线| 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 久久这里只有精品19| 国产精品亚洲一级av第二区| 午夜亚洲福利在线播放| 精品国产乱子伦一区二区三区| 99re在线观看精品视频| 后天国语完整版免费观看| 又黄又爽又免费观看的视频| 久久久久久久午夜电影| 欧美大码av| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 国语自产精品视频在线第100页| 国产三级黄色录像| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 午夜免费激情av| 亚洲最大成人中文| 色综合站精品国产| 国产aⅴ精品一区二区三区波| www.自偷自拍.com| 一本久久中文字幕| 特级一级黄色大片| 久久久久久久午夜电影| 黑人欧美特级aaaaaa片| 欧美av亚洲av综合av国产av| 国产成人一区二区三区免费视频网站| 国产亚洲av高清不卡| 男女之事视频高清在线观看| 亚洲自偷自拍图片 自拍| 欧美成人性av电影在线观看| 日本成人三级电影网站| 久久九九热精品免费| 一进一出抽搐gif免费好疼| 三级毛片av免费| 国产三级在线视频| 亚洲国产中文字幕在线视频| 亚洲av美国av| 成人精品一区二区免费| 国产黄色小视频在线观看| 欧美乱色亚洲激情| 精品久久久久久久末码| 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 99热6这里只有精品| 亚洲美女黄片视频| 亚洲精品乱码久久久v下载方式 | 亚洲性夜色夜夜综合| 国产精品乱码一区二三区的特点| 嫩草影院入口| 麻豆久久精品国产亚洲av| 亚洲成人久久爱视频| 亚洲精品一区av在线观看| 91麻豆av在线| 黄片小视频在线播放| 精品国内亚洲2022精品成人| 国产精品av久久久久免费| 亚洲国产精品合色在线| 午夜免费激情av| ponron亚洲| e午夜精品久久久久久久| 97人妻精品一区二区三区麻豆| 欧美一级毛片孕妇| 色吧在线观看| 美女被艹到高潮喷水动态| av在线天堂中文字幕| 国产精品国产高清国产av| 日韩高清综合在线| 好男人电影高清在线观看| 国产成人欧美在线观看| 日韩欧美一区二区三区在线观看| 亚洲 欧美一区二区三区| 色在线成人网| 国内久久婷婷六月综合欲色啪| 岛国视频午夜一区免费看| 美女被艹到高潮喷水动态| 一区二区三区高清视频在线| 亚洲五月婷婷丁香| 九九在线视频观看精品| 国产麻豆成人av免费视频| 日韩欧美精品v在线| 亚洲avbb在线观看| 男女做爰动态图高潮gif福利片| 国内精品久久久久久久电影| 男女下面进入的视频免费午夜| 天堂影院成人在线观看| 欧美一区二区国产精品久久精品| 成人高潮视频无遮挡免费网站| 啪啪无遮挡十八禁网站| 两个人的视频大全免费| www国产在线视频色| 欧美午夜高清在线| 看片在线看免费视频| 少妇裸体淫交视频免费看高清| 欧美日韩精品网址| 成人av在线播放网站| 村上凉子中文字幕在线| 日韩欧美国产在线观看| 色综合站精品国产| 成年版毛片免费区| 久9热在线精品视频| 村上凉子中文字幕在线| 中文在线观看免费www的网站| 国产精品九九99| 又黄又粗又硬又大视频| 桃红色精品国产亚洲av| 日本五十路高清| 欧美中文综合在线视频| 亚洲精品粉嫩美女一区| 日韩国内少妇激情av| 日本 欧美在线| 一进一出抽搐gif免费好疼| 国产成人欧美在线观看| 亚洲av美国av| 国产精品亚洲美女久久久| 午夜福利在线观看吧| 麻豆久久精品国产亚洲av| 欧美3d第一页| 久久这里只有精品中国| 日日摸夜夜添夜夜添小说| 51午夜福利影视在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美日韩一级在线毛片| 在线观看一区二区三区| 欧美最黄视频在线播放免费| 51午夜福利影视在线观看| 精品国产超薄肉色丝袜足j| 精品无人区乱码1区二区| 国产高清激情床上av| 可以在线观看毛片的网站| 搡老妇女老女人老熟妇| 亚洲 欧美 日韩 在线 免费| 日本黄色片子视频| 美女黄网站色视频| 国产伦在线观看视频一区| 又大又爽又粗| 日韩中文字幕欧美一区二区| 免费搜索国产男女视频| 欧美日韩综合久久久久久 | 九色成人免费人妻av| 亚洲av电影在线进入| 美女被艹到高潮喷水动态| 国产亚洲欧美在线一区二区| 精品久久久久久久人妻蜜臀av| 极品教师在线免费播放| 日本熟妇午夜| 91久久精品国产一区二区成人 | 桃色一区二区三区在线观看| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 国产真人三级小视频在线观看| 国产成人啪精品午夜网站| 亚洲乱码一区二区免费版| 国产97色在线日韩免费| 亚洲国产看品久久| 欧洲精品卡2卡3卡4卡5卡区| 宅男免费午夜| 亚洲av成人精品一区久久| 岛国视频午夜一区免费看| 亚洲欧美日韩高清在线视频| 19禁男女啪啪无遮挡网站| 99国产精品99久久久久| 最近视频中文字幕2019在线8| 亚洲中文日韩欧美视频| a级毛片a级免费在线| 久久精品国产综合久久久| 成年女人永久免费观看视频| 精品久久久久久久人妻蜜臀av| or卡值多少钱| 老司机深夜福利视频在线观看| 天天添夜夜摸| 搞女人的毛片| 亚洲av熟女| 成人国产综合亚洲| 91字幕亚洲| 成年免费大片在线观看| 午夜福利18| 欧美成人免费av一区二区三区| 国产伦精品一区二区三区视频9 | 嫁个100分男人电影在线观看| 国产99白浆流出| 永久免费av网站大全| 毛片一级片免费看久久久久| 麻豆一二三区av精品| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品.久久久| 一区二区三区乱码不卡18| 亚洲国产高清在线一区二区三| 国产精品美女特级片免费视频播放器| 村上凉子中文字幕在线| 亚洲国产精品国产精品| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 一级毛片久久久久久久久女| 国产成人freesex在线| 嫩草影院精品99| 黄色配什么色好看| 成人无遮挡网站| 边亲边吃奶的免费视频| 女人被狂操c到高潮| 热99在线观看视频| 日韩视频在线欧美| 久久久精品94久久精品| 久久久久免费精品人妻一区二区| 岛国在线免费视频观看| 老司机影院毛片| 禁无遮挡网站| 国产乱人视频| 久久久国产成人精品二区| 久久精品国产自在天天线| 真实男女啪啪啪动态图| 男女国产视频网站| 少妇熟女aⅴ在线视频| av国产久精品久网站免费入址| 精品熟女少妇av免费看| 欧美xxxx黑人xx丫x性爽| 欧美成人精品欧美一级黄| 日韩制服骚丝袜av| 午夜精品一区二区三区免费看| 亚洲精品国产成人久久av| 日韩,欧美,国产一区二区三区 | 在线播放无遮挡| 国产精品久久久久久av不卡| 在线天堂最新版资源| 亚洲精品国产av成人精品| 国产精品国产三级专区第一集|