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

    High-Throughput Profiling of Serological Immunoglobulin G N-Glycome as a Noninvasive Biomarker of Gastrointestinal Cancers

    2024-01-22 10:32:42PngchngLiuXiobingWngAishDunYutongLiHouqingLiLuWngYichunZhngCncnLiJinxiZhngXioyuZhngLixingHingHou
    Engineering 2023年7期

    Pngchng Liu, Xiobing Wng, Aish Dun, Yutong Li, Houqing Li, Lu Wng, Yichun Zhng,Cncn Li, Jinxi Zhng, Xioyu Zhng, Lixing M, Hing Hou,g,*

    a School of Public Health, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan 250117, China

    b State Key Laboratory of Molecular Oncology, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100021, China

    c School of Stomatology, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan 250117, China

    d Beijing Key Laboratory of Clinical Epidemiology, School of Public Health, Capital Medical University, Beijing 100069, China

    e Beijing Sanbo Brain Hospital, Capital Medical University, Beijing 100093, China

    f Department of Gastroenterology, The Second Affiliated Hospital of Shandong First Medical University, Taian 271000, China

    g Department of Epidemiology, School of Public Health, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan 250117, China

    Keywords:Gastrointestinal cancer Glycosylation Immunoglobulin G Diagnostic biomarker

    ABSTRACT Immunoglobulin G (IgG) N-glycosylation plays a crucial role in the development of inflammatory diseases.This study aimed to evaluate the diagnostic performance of IgG for gastrointestinal(GI)cancer subtypes.A total of 749 GI cancer patients were enrolled from the Cancer Hospital, Chinese Academy of Medical Sciences,including esophageal cancer(EC),gastric cancer(GC),colorectal cancer(CRC),and pancreatic cancer(PC)patients.Hydrophilic interaction liquid chromatography using ultra-performance liquid chromatography (HILIC-UPLC) was employed to analyze the composition of the plasma IgG Nglycome.The levels of circulating inflammatory cytokines were detected by means of a Bio-Plex Pro Human Th17 Cytokine Assay.Canonical correlation analysis (CCA) was used to explore the correlation between IgG N-glycosylation patterns and inflammatory cytokines.A Lasso algorithm, accompanied by a logistic regression model, was used to develop a glycan-based model for differentiating GI cancer patients from healthy individuals.The levels of sialylation and galactosylation were significantly decreased among EC, GC, CRC, and PC patients, whereas the abundance of glycans with bisecting Nacetylglucosamine(GlcNAc)was increased in GI cancer patients in comparison with the healthy controls.Moreover, only PC patients had a decreased level of fucosylation.The levels of interleukin 1β (IL-1β), IL-31,and soluble CD40 ligand(sCD40L)were significantly higher in GI cancer patients than in the controls.In addition, the composition of IgG N-glycans was correlated with that of inflammatory cytokines(r=0.556).The glycan-based models for diagnosing GI cancers exhibited an excellent performance,with areas under the receiver operating characteristic curves (AUCs) of 0.972 for EC, 0.871 for GC, 0.867 for CRC, and 0.907 for PC.Our findings demonstrate that IgG N-glycosylation plays an important role in modulating the pathogenesis of GI cancers.Serological IgG N-glycosylation is thus a potential candidate for noninvasively assisting in the clinical diagnosis of GI cancer subtypes.

    1.Introduction

    Gastrointestinal (GI) cancers, including digestive tract (mouth,throat,esophagus,stomach, colorectum, and appendix) and digestive organ(pancreatic and hepatobiliary)carcinomas,are the most common malignancies, accounting for approximately 39% of all cancer deaths worldwide [1].In China, among the new cancer patients reported in 2020 (nearly 4.5 million), 478 508 (10.47%)had gastric cancer (GC), 555 477 (12.16%) had colorectal cancer(CRC), and 324 422 (7.10%) had esophageal cancer (EC) [2].The main health challenges caused by GI cancers have imposed a heavy burden on the healthcare system [3].Early diagnosis has been shown to be of great clinical significance for cancer management;however, commonly used noninvasive diagnostic approaches lack sufficient accuracy.

    Glycosylation is an enzyme-catalyzed post-translational reaction by which oligosaccharides are transferred to specific sites on biological macromolecules (e.g., proteins and lipids) and together form glycoconjugates (i.e., glycoproteins and glycolipids) [4–6].Glycosylation aberrations can be detected in the tumor tissues and biological fluids of cancer patients.The glycosylation signatures change dramatically along with the process of malignant transformation and tumor progression [7,8].As such, manifold serological glycoproteins have been widely employed as tumor markers, such as carbohydrate antigen 19–9 (CA19-9) for GI cancers, carcinoembryonic antigen (CEA) for colon cancer, alphafetoprotein (AFP) for liver cancer, carbohydrate antigen 125(CA125) for ovarian cancer, and prostate-specific antigen (PSA)for prostate cancer[9].However,these biomarkers possess insufficient sensitivity (Se) and specificity (Sp) [10].In addition, cancerspecific cell-surface glycophenotypes have been investigated in clinical settings and serve as treatment targets rather than as diagnostic biomarkers for high-throughput screening at the population level.Serological glycoproteins that carry out regulating functions in tumorigenesis,as well as relevant detection techniques,must be identified for the diagnosis and subtype differentiation of GI cancers.

    As a major component of humoral immunity and the most abundant immunoglobulin in human plasma, immunoglobulin G(IgG) plays a vital role in the inflammatory response [11,12].Each IgG contains two N-linked glycosylation sites at asparagine 297(Asn-297) of the fragment crystallizable (Fc) segment in heavy chains [13].N-glycans in the Fc region regulate the function of IgG effectors by affecting the binding affinity of Fcγ receptors[4,14,15].IgG glycosylation responds quickly to physiological or pathological changes but remains relatively stable within healthy individuals [16].By detecting glycoconjugates, the features of cancer-associated IgG glycosylation can be displayed.Studies have explored the potential role of the glycol phenotypes of IgG in hepatic, colorectal, gastric, lung, and ovarian cancers [17–20].

    Taking advantage of the progress that has been achieved in glycomics analyses, we conducted this investigation to comprehensively analyze the features of IgG N-glycosylation in four types of GI cancer—that is, EC, GC, CRC, and pancreatic cancer (PC)—and to evaluate the performance of glycan-based models for diagnosing GI cancers and subtypes.

    2.Materials and methods

    2.1.Study participants

    This case-control study included 749 GI cancer patients who were hospitalized at the Cancer Hospital, Chinese Academy of Medical Sciences, Beijing, China, between October 2018 and January 2020 (Fig.S1 in Appendix A).Among them, 100 patients had EC, 121 had GC, 328 had CRC, and 200 had PC.In addition,112 healthy participants of similar ages were recruited as controls.

    The inclusion criteria for cancer patients were as follows:①patients diagnosed with EC (International Classification of Diseases Tenth Revision (ICD-10) code, C15), GC (C16), CRC (C18), or PC(C25)through histopathological examination of tissue biopsies;②histological diagnostic criteria conforming to the World Health Organization (WHO) Classification of Tumours (4th edition), Volume 1: Digestive System Tumours [21]; ③no other severe somatic or psychiatric disorders; and ④not enrolled in other clinical trials.The exclusion criteria were:①patients with severe infectious diseases and cardiovascular diseases, such as coronary heart disease(I25.1) and stroke (I64); ②patients with psychiatric disorders(F99); and ③patients who declined to participate in this study.

    The study protocol was approved by the ethics committee of the National Cancer Center/Cancer Hospital, Chinese Academy of Medical Sciences (No.NCC1839).Written informed consent was obtained from each participant.

    2.2.Collection and storage of blood samples

    Blood samples were collected by venipuncture after 12 h of overnight fasting.The separated plasma and serum were used for the detection of IgG N-glycans and routine biochemical tests,respectively.All collected blood samples were processed within 8 h and stored at -80 °C until further measurement.

    2.3.Measurements of demographic and clinical variables

    Demographic information was obtained through in-person interviews with participants and/or their families.Hematology and biochemical parameters, including fasting blood glucose(FBG), total cholesterol (TC), high-density lipoprotein cholesterol(HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides(TG),were determined with an automatic analyzer(Hitachi,Japan) (Table 1).

    2.4.Analysis of cancer biomarkers and inflammatory cytokines

    The levels of circulating CA19-9 and CEA were measured by means of immunochemiluminometric assays (AutoLumo A2000Plus, China).The normal range of CA19-9 was 0–37 U?mL-1,and participants with CA19-9 > 37 U?mL-1were defined as abnormal [22].Serological CEA levels were classified into two groups:normal (≤5.0 ng?mL-1) and elevated (> 5.0 ng?mL-1) [23].

    To identify the levels of inflammatory cytokines among GI cancers, a subset of the participants (45 EC, 45 GC, 45 CRC, and 45 PC patients, and 48 healthy controls) were randomly selected to test against nine inflammatory cytokines—that is, interleukin-1β (IL-1β), IL-4, IL-6, IL-17A, IL-17F, IL-31, IL-33, soluble CD40 ligand(sCD40L), and tumor necrosis factor-α (TNF-α), using the Bio-Plex 200 suspension array system (Bio-Rad, USA) with the Bio-Plex Pro Human Th17 Cytokine 15-plex Assay (Bio-Rad) [24].

    2.5.Analysis of IgG glycans

    Isolation of plasma IgG and release of IgG Fc N-linked glycans were performed in accordance with the approaches established previously [25,26].In brief, a 100 μL plasma sample was added to a 96-well protein G monolithic plate(BIA Separations,Slovenia)for IgG isolation.Isolated IgG samples were denaturized with 1 mL of 0.1 mol?L-1formic acid and immediately neutralized with 1 mol?L-1ammonium bicarbonate.For N-glycan release, 4 μL of PNGase F enzyme was added to the IgG samples and incubated in a 37 °C water bath for 18 h.

    The released N-glycans were labeled with 2-aminobenzamide(2-AB) and then transferred into an oven at 65 °C for 3 h.The 2-AB labeled glycans were analyzed with hydrophilic interaction liquid chromatography using an ultra-performance liquid chromatography (HILIC-UPLC) instrument (Walters Corporation, USA), by which 24 IgG glycan peaks (GPs) were detected.The structures of these GPs were identified by reference to a dataset established by mass spectrometry [27].The level of each glycan was quantitated by the percentage of the area of the relevant chromatogram peak to the integrated area of all glycans.Furthermore, we calculated the 54 derived glycan traits using the measurements of the initial glycans, which consisted of sialylation, bisectingN-acetylglucosamine (GlcNAc), galactosylation, and fucosylation traits (Table S1 in Appendix A) [28].

    Table 1 Characteristics of the study participants.

    2.6.Statistical analysis

    The assumption of normal distribution of continuous data was tested using the Kolmogorov–Smirnov test.The mean and standard deviation (SD) were used to represent normally distributed data;otherwise, the median (M) and interquartile range (IQR) were employed when a normal distribution was absent.One-way analysis of variance (ANOVA) was performed to analyze the betweengroup difference in normally distributed variables; otherwise, the nonparametric test (i.e., Kruskal–Wallis test) was carried out.Moreover, multiple comparison correction was performed with the Benjamini–Hochberg method of false discovery rate(FDR)correction.Categorical variables were expressed by frequency and percentage (%) and analyzed by the chi-squared test (χ2test),accompanied by a post hoc Bonferroni correction.Canonical correlation analysis (CCA) was used to explore the overall correlation between the 24 glycan structures (x) and inflammatory cytokines(y).Identified variables with a statistically significant impact on the canonical variables were judged by the canonical loadings.In general,an absolute value greater than 0.30 was used to define significant loadings [29].A clustering analysis and heatmap were undertaken to further explore the inflammatory cytokines significantly associated with GI cancers [30].

    We established glycan-based diagnostic models for GI cancers.Internal validation was performed by randomly splitting the participants into training and validation sets at a certain ratio (7:3)[31].First, the least absolute shrinkage and selection operator(Lasso) algorithm was used to screen initial glycans that were significantly associated with GI cancers,thus avoiding overfitting and effectively controlling the model’s complexity.Significant glycans selected with Lasso were then introduced into multivariate logistic regression analysis, with the forward maximum likelihood ratio test being taken as the stopping rule during the model building;in this way, adjusted odds ratios (ORs) and 95% confidence intervals (CIs) were then obtained.A diagnostic model for each GI cancer was developed using the significant glycans.Then, a receiver operator characteristic (ROC) curve was created to evaluate the performance of this model, by which the area under the curve(AUC) value, Se, and Sp were calculated.R software version 4.1.1(R Core Team, New Zealand) and SPSS software version 25.0(IBM,USA) were used for statistical analyses.The R package‘‘glmnet”was used to run the Lasso algorithm.ROC curves and AUC values were generated using the ‘‘pROC” R package and GraphPad Prism(GraphPad Software,Inc.,USA).A P value<0.05 was considered statistically significant.

    3.Results

    3.1.Characteristics of study participants

    The demographics,biochemical indicators,and clinical features of the 749 GI cancer patients(mean age 58 years,484 men and 265 women)and 112 healthy controls are listed in Table 1,as well as in Tables S2–S11 and Figs.S2 and S3 in Appendix A.No statistically significant difference in age was detected between each cancer group and the controls.Significantly lower levels of serum TC and HDL-C were observed among the EC, GC,CRC,and PC patients compared with the healthy control group.The FBG in the EC, GC,and CRC groups was also significantly lower than that in the control group.The prevalence of dyslipidemia was significantly higher in the GC,CRC,and PC patients than in the controls.In addition,the PC patients had a higher prevalence of hyperglycemia than the GC,CRC, and control groups.

    3.2.Description of clinically used cancer biomarkers

    The level of CA19-9 in the healthy control group was 7.63 U?mL-1(IQR:5.11,10.70),which was significantly lower than that in the GC (M = 9.23; IQR: 5.15, 18.95), CRC (M = 12.56; IQR: 7.50,25.31), and PC (M = 66.88; IQR: 11.89, 529.83) groups.The level of CEA in the control group was 1.53 ng?mL-1(IQR: 1.21, 2.00),which was significantly lower than that in the EC (M = 2.04;IQR: 1.19, 2.97), GC (M = 1.78; IQR: 1.04, 3.73), CRC (M = 3.07;IQR: 1.48, 8.37), and PC (M = 2.76; IQR: 1.80, 6.40) groups.

    3.3.Inflammatory cytokines

    As shown in Tables S12–S22 in Appendix A,the circulating concentrations of IL-1β, IL-31, sCD40L, and TNF-α were significantly higher in the four cancer groups than in the control group.Compared with the EC group, PC patients had significantly higher plasma levels of IL-31, sCD40L, and TNF-α.No significant differences were observed for IL-4,IL-6,IL-17A,IL-17F,or IL-33 between the groups.

    A clustering analysis was conducted to identify the patterns of circulating inflammatory cytokines in GI cancers.As shown in Fig.1, individuals with GI cancers shared similar inflammatory cytokines, which were significantly different from those of the healthy controls (e.g., sCD40L and IL-31).

    Fig.1.Clustering analysis of the patterns of inflammatory cytokines in GI cancers.

    3.4.The IgG glycome composition in GI cancer patients

    Our HILIC-UPLC assay initially identified 24 oligosaccharide chains that were linked to the IgG Fc segment.Of these, there were significant differences in the abundance of 21 glycans between the EC, GC, CRC, PC, and control groups (Table 2;Figs.S4–S7 in Appendix A).We further calculated 54 derived glycan traits using the measurements of the initial glycans [28].As indicated in Tables S23–S33 in Appendix A, 48 of the 54 derived glycan traits differed between the four types of GI cancer groups and the controls.

    Moreover, to explore the modulatory role of IgG N-glycans in the development of GI cancers, we calculated the summary abundances of the four major glycosylation features—that is,sialylation,bisecting GlcNAc, galactosylation, and fucosylation (Table S34 in Appendix A) [13].

    3.4.1.Galactosylation

    As depicted in Table 3, the percentages of glycan structures without a terminal galactose (agalactosylation, G0) were 37.37%in EC, 36.00% in GC, 36.78% in CRC, and 36.85% in PC patients,which were significantly higher than those in the control group(31.42%).The galactosylation traits include glycans with one terminal galactose (monogalactosylation, G1) and two galactoses(digalactosylation, G2) [6].The abundances of G1 in the EC, GC,CRC, and PC groups were 42.22%, 42.90%, 42.98%, and 43.05%,respectively, which were lower than those in the control group(44.47%).In addition, the abundances of G2 in the EC, GC, CRC,and PC patients (19.75%,20.40%,19.30%, and 18.99%,respectively)were also lower than those in the control group (22.93%).

    3.4.2.Sialylation

    In total, a lower level of sialylation was identified in the EC(20.70%), GC (20.93%), CRC (20.22%), and PC (20.86%) groups than in the controls (21.63%) (Table 3).

    Table 2 A comparison of abundance (%) between 24 initial glycans in the studied groups.

    Table 3 A comparison of the relative abundance (%) of four IgG glycome features.

    3.4.3.Bisecting GlcNAc

    The percentage of glycans with bisecting GlcNAc was significantly higher in the EC (17.68%), GC (17.10%), CRC (17.31%), and PC (17.09%) patients than in the controls (15.92%) (Table 3).

    3.4.4.Fucosylation

    The level of fucosylation was significantly lower in the PC patients (94.96%) than in the controls (95.49%) (Table 3).

    These findings indicated a significant decrease in galactosylated, sialylated, and fucosylated glycans among GI cancers,while bisecting GlcNAc was significantly higher among these glycans.Details on galactosylation, sialylation, bisecting GlcNAc,and fucosylation are provided in Tables S35-S38 and Box S1 in Appendix A.

    3.5.Association of IgG N-glycans with GI cancers

    As listed in Table S39 in Appendix A, after being adjusted for sex, FBG, TC, HDL-C, and the prevalence of dyslipidemia, four glycans were significantly associated with EC: GP10 (OR: 7.371; 95%CI: 2.301, 23.608; P = 0.001); GP11 (OR: 5.943 × 10-6; 95% CI:3.061 × 10-8, 0.001; P < 0.001); GP14 (OR: 0.758; 95% CI: 0.607,0.948; P = 0.015); and GP23 (OR: 0.195; 95% CI: 0.055, 0.697;P = 0.012).Meanwhile, sex, FBG, TC, HDL-C, and the prevalence of dyslipidemia showed no significance.

    For the logistic regression analysis of GC, GP5 (OR: 4.489 ×10-4; 95% CI: 1.817 × 10-5, 0.011; P < 0.001), GP6 (OR: 2.652;95% CI: 1.625, 4.326; P < 0.001), sex (OR: 0.190; 95% CI: 0.070,0.519; P = 0.001), FBG (OR: 0.585; 95% CI: 0.435, 0.787;P <0.001), and TC(OR:0.438;95%CI:0.25,0.745;P =0.002) were of significance (Table S40 in Appendix A).

    With regard to CRC,two glycans and three covariates were statistically significant: GP2 (OR: 12.941; 95% CI: 3.840, 43.616;P < 0.001), GP20 (OR: 1.305 × 10-10; 95% CI: 3.535 × 10-14,4.818 × 10-7; P < 0.001), sex (OR: 0.428; 95% CI: 0.203, 0.901;P = 0.026), FBG (OR: 0.647; 95% CI: 0.496, 0.844; P = 0.001), and HDL-C (OR: 0.222; 95% CI: 0.087, 0.566; P = 0.002) (Table S41 in Appendix A).

    Four glycans and one covariate were significantly associated with PC: GP17 (OR: 53.358; 95% CI: 9.165, 310.659; P < 0.001), GP19(OR: 3.855; 95% CI: 1.506, 9.869; P = 0.005), GP20 (OR: 3.076 ×10-5; 95% CI: 1.251 × 10-7, 0.008; P < 0.001), GP23 (OR: 0.115;95% CI: 0.046, 0.286; P < 0.001), and HDL-C (OR: 0.061; 95% CI:0.022,0.175;P<0.001)(Table S42 in Appendix A).

    3.6.Association between IgG N-glycans and inflammatory cytokines

    Fig.2.Canonical structures of the IgG N-glycans and inflammatory cytokines.An absolute value of canonical loadings greater than 0.30 was considered significant.All variables are sorted by the absolute value of their canonical loadings.Positive relationships are represented in red boxes, while negative relationships are shown in black boxes.

    The compositions of IgG N-glycans were significantly correlated with inflammatory cytokines, with a canonical correlation coefficient (r) of 0.556 (P = 0.005) (Fig.2).Six glycans (GP1, GP2, GP4,GP14, GP15, and GP18) were associated with sCD40L, IL-31, TNFα, IL-17A, and IL-6.Moreover, a strong association was observed between GP18 and the canonical variables, with a loading of 0.497,and the response variable with the highest canonical loading was -0.669 (sCD40L).

    3.7.Discrimination of EC,GC,CRC,and PC from healthy controls based on IgG N-glycosylation

    Fig.3.Logistic regression analyses of the associations between the initial glycans and GI cancers.(a) EC; (b) GC; (c) CRC; (d) PC.

    Fig.4.Receiver operating characteristic curves (ROC) with AUC values for glycan-based prediction models.Glycan model of (a) EC; (b) GC; (c) CRC; (d) PC.

    As shown in Figs.3 and 4, to establish a diagnostic model of GI cancers on the basis of the 24 initial glycans,a Lasso algorithm was used and a logistic regression analysis was carried out,according to which four glycans (GP10, GP11, GP14, and GP23) were employed for the diagnosis of EC(Fig.3(a));a ROC curve analysis was used to evaluate the performance of this diagnostic model (Fig.4(a)).The AUC value was 0.972 (95% CI: 0.945, 0.998), with an Se of 95.3%and an Sp of 95.0%, indicating that the glycan-based model had a significantly better performance than CA19-9 (AUC: 0.693; 95%CI: 0.603, 0.783) and CEA (AUC: 0.797; 95% CI: 0.723, 0.871).

    GP5,GP6,and GP18 were employed for the diagnosis of GC,and also demonstrated a better performance (AUC: 0.871; 95% CI:0.815, 0.926) compared with CA19-9 (AUC: 0.662; 95% CI: 0.577,0.747) and CEA (AUC: 0.731; 95% CI: 0.653, 0.810) (Figs.3(b) and 4(b)).

    In addition,the diagnostic model of CRC,which consisted of GP2 and GP20,showed a significantly better performance(AUC:0.867;95% CI: 0.816, 0.917) than CA19-9 (AUC: 0.521; 95% CI: 0.453,0.590) and CEA (AUC: 0.530; 95% CI: 0.466, 0.594) (Figs.3(c) and 4(c)).

    Furthermore, GP6, GP8, GP10, GP15, GP17, GP19, GP20, GP21,and GP23 were selected for the diagnosis of PC, and also demonstrated a significantly better performance (AUC: 0.907; 95% CI:0.865, 0.949) than CA19-9 (AUC: 0.690; 95% CI: 0.620, 0.759) and CEA (AUC: 0.562; 95% CI: 0.486, 0.639) (Figs.3(d) and 4(d)).

    In the validation dataset, the glycan-based models had AUCs of 0.991 (95% CI: 0.977, 0.999), 0.874 (95% CI: 0.791, 0.957), 0.856(95% CI: 0.785, 0.928), and 0.872 (95% CI: 0.800, 0.944) for EC,GC, CRC, and PC, respectively (Fig.5).

    3.8.Comparison between glycans in early and advanced GI cancers

    As shown in Tables S43–S46 in Appendix A, the glycosylation patterns were similar in early-and advanced-stage cancer patients.

    3.9.Diagnostic model for GI cancer patients in early and advanced stages

    In this study,clinical stages I–II were defined as early stage,and stages III–IV were defined as advanced stage (Fig.S3 in Appendix A).Taking 24 glycans as a starting point, a diagnostic model for early-stage EC was established, with four glycans (GP11, GP21,GP23,and GP24)being included in the model(Table S47 in Appendix A).Three glycans (GP11, GP17, and GP23) were selected to establish a diagnostic model for advanced EC(Table S47 in Appendix A).

    Fig.5.ROC curves in the internal validation set.(a) EC; (b) GC; (c) CRC; (d) PC.

    Regarding a subgroup analysis of GC, two glycans (GP5 and GP6)were employed for the diagnosis of early GC,and two glycans(GP4 and GP5)were selected for developing a diagnostic model for advanced GC (Table S48 in Appendix A).With regard to CRC, four glycans (GP8, GP10, GP11, and GP20) were employed for the diagnosis of early-stage patients,and two glycans (GP2 and GP18)were selected for the diagnosis of advanced patients (Table S49 in Appendix A).Furthermore, two glycans (GP2 and GP13) were employed for the diagnosis of early PC, and six glycans (GP7,GP8,GP17,GP18,GP20, and GP23)were selected for the diagnosis of advanced PC (Table S50 in Appendix A).Details of the model validation are provided in Fig.S8 in Appendix A.

    4.Discussion

    The massive amount of data generated by high-throughput genomics and proteomics experiments has advanced our understanding of cancer; nevertheless, there is still a surprising scarcity of well-validated,clinically useful biomarkers.In addition to genomics and proteomics, the emerging field of glycomics is gaining importance in glycomedicine and cancer research [32].Glycan moieties of IgG Fc segments (i.e., galactose, sialic acid, bisecting GlcNAc, and fucose) participate in the pathophysiological process of inflammatory disorders and provide potential for the early diagnosis and targeted treatment of inflammatory diseases, including cancers (Table S51 in Appendix A) [9,25,28,33,34].

    This study comprehensively analyzed serological IgG N-glycan profiles in four types of GI cancer (EC, GC, CRC, and PC) by means of an HILIC-UPLC-based high-throughput method.We demonstrated a decrease in IgG sialylation and galactosylation and an increase in bisecting GlcNAc among GI patients.These alterations in IgG glycosylation induced a pathological increase in inflammatory cytokines, such as IL-1β, IL-31, sCD40L, and TNF-α, and contributed to the development of GI cancers.In addition, the IgG N-glycosylation trait has the capacity to serve as a potential biomarker for GI cancer diagnosis and subtype differentiation from healthy individuals, and demonstrates a significantly higher performance than clinically used serological biomarkers such as CA19-9 and CEA.Consequently, high-throughput assay of IgG N-glycans is expected to be exploited for the early diagnosis and targeted treatment of GI cancers.

    The mechanism of tumorigenesis is regulated by both acute and chronic inflammation [35,36].A variety of molecules, especially pro-inflammatory cytokines (e.g., IL-1β, IL-6, and TNF-α), are engaged in the development of cancer [37–40].Studies have demonstrated that IgG N-glycans modulate the pro- and antiinflammatory balance in both inflammatory and autoimmune diseases [15,41–43].In the current study, an overall correlation between GI cancer-related IgG N-glycans and inflammatory cytokines was identified, supporting previous findings in the investigation of patients with primary tumors of the GI tract [44].These findings emphasize the role of IgG N-glycosylation in the regulation of inflammatory response.

    IgG sialylation—also known as sialic acid or N-acetylneuraminic acid (Neu5Ac)—plays a role in the anti-inflammatory properties of IgG [45–47].During the development of IgG N-glycans, the terminal sialic acid residue covalently attaches to galactose and leads to a reduced capacity to ligate FcγRIIIa on natural killer(NK) cells, thereby lowering inflammatory activity via antibodydependent cell cytotoxicity (ADCC) [48,49].The absence of sialic acid dramatically impacts the function of IgG, shifting from anti- to pro-inflammatory.Consistent with our findings, a decrease in sialylation has also been associated with GI cancers[7,8,50,51].

    Galactosylation can promote IgG of its affinity for the inhibitory FcγRIIb and thereby increase its anti-inflammatory activities [52].A decrease in IgG galactosylation—termed agalactosylation—can expose GlcNAc residues; this increases IgG binding with mannose-binding lectin, which results in the upregulation of CDC activity and leads to the initiation of inflammation[5,52–54].Studies have reported that the complement effector system is extensively activated in various tumor microenvironments, leading to tumor initiation and growth [55].Consistent with our findings,an elevated level of agalactosylated IgG has been detected in GC,CRC,and PC patients,indicating its role in the development of cancers [8,51,54,56–58].Studies have also suggested that agalactosylation is strongly correlated with an increase in proinflammatory cytokines such as TNF-α and C-reactive protein(CRP) [59,60].Our findings showed significantly higher TNF-α levels in EC,GC,CRC,and PC patients,evidencing that the agalactosylation of IgG plays a crucial role in the upregulation of the inflammatory process in cancer.

    Fucosylation has been intensively studied due to its mechanisms of antibody-based therapeutics against tumors via ADCC.Fucosylation on IgG inhibits ADCC via decreased FcγRIIIa binding on NK cells, downregulating the pro-inflammatory activity of IgG[61,62].The ADCC activity of IgG without a core fucose on the N-glycan is enhanced by up to 100-fold [63].A decreased level of fucosylation was observed in PC patients in our study, supporting the findings in another Chinese PC patients, among which significantly lower fucosylation, sialylation and galactosylation were demonstrated[7].Moreover,lower fucosylated IgG promotes the production of pro-inflammatory cytokines that are secreted by monocytes,including IL-1β,IL-6,and TNF-α[64].Similarly,fucoserich IgG reduces the production of TNF-α and IL-6, resulting in anti-inflammatory activity [65].Our findings showed that TNF-α,IL-1β,IL-31,and sCD40L levels were significantly higher in GI cancer patients than in the controls, supporting the abovementioned findings [66–68].

    The bisecting GlcNAc of IgG regulates ADCC by causing an increase in affinity to FcγRIIIa, leading to the pro-inflammatory function of IgG[47,69].This disruption in the balance of inflammation progression is considered to play an important role in the development of carcinomas [62,70].In the present study, we observed an increase in bisecting GlcNAc among EC, GC, CRC, and PC patients, which was consistent with analogous findings in previous reports [7,50,54,57].An in vitro experiment showed that a decrease in bisecting GlcNAc was associated with increasing levels of IL-21[71].However,no IL-21 values were detected in the serum of our study participants.

    Although ADCC is a crucial pathway against tumors,the upregulation of ADCC induced by abnormal IgG N-glycans has not been demonstrated to be associated with more accurate diagnosis.Studies have reported that decreased galactosylation, decreased sialylation, and increased bisecting GlcNAc in IgG N-glycans are associated with poorer prognosis among CRC patients [51].

    Massive quantities of serum biomarkers have been investigated to monitor the clinical progression of cancer or to define treatment regimens.Although some of these tumor biomarkers play an important role in disease management, they are of limited use as diagnostic tools due to insufficient Sp or Se.At present,the noninvasive IgG N-glycan assay provides a high level of accuracy in the diagnosis of inflammatory and immune diseases [72].The combination of the HILIC analysis of fluorescently tagged N-glycans with UPLC analysis has been demonstrated to be a valuable approach[27].In the current study, the high-throughput analysis of IgG Nglycosylation presented significantly high AUC values of 0.972,0.871, 0.867, and 0.907 for distinguishing EC, GC, CRC, and PC,respectively.Consistent with findings in prostate cancer,lung cancer, and bladder cancer[73–75], the present findings demonstrate the potential of IgG N-glycosylation for the early diagnosis of GI cancers.

    This study has several limitations that should be interpreted.First, the case-control study design limited the inference of the causal effect between IgG N-glycosylation and diseases.Second,even though glycosylation can have major etiological effects in cancers, environmental and lifestyle factors were not measured in this study,although they play important roles in the development of cancers.

    In conclusion, our findings indicate that decreases in IgG sialylation and galactosylation, together with an increase in bisecting GlcNAc,might play important roles in the carcinogenesis and progression of GI cancers.Serological IgG N-glycosylation may serve as a potential candidate for noninvasively assisting the clinical diagnosis of GI cancer subtypes.To the best of our knowledge, this is the first report on integrative comparisons of IgG N-glycosylation profiles in EC, GC, CRC, and PC patients; it offers opportunities for the identification of new diagnostic biomarkers for GI cancers.

    Acknowledgments

    This work was supported by the Shandong Provincial Natural Science Foundation (ZR2022MH082) and the National Natural Science Foundation of China (81872682).We thank Professor Youxin Wang at Capital Medical University for his helpful suggestions on this work.

    Compliance with ethics guidelines

    Pengcheng Liu,Xiaobing Wang,Aishe Dun,Yutong Li,Houqiang Li, Lu Wang, Yichun Zhang, Cancan Li, Jinxia Zhang, Xiaoyu Zhang,Lixing Ma, and Haifeng Hou declare that they have no conflicts of interest or financial conflicts to disclose.

    Appendix A.Supplementary data

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

    欧美激情国产日韩精品一区| 免费av不卡在线播放| 少妇被粗大猛烈的视频| 亚洲人成网站高清观看| 久久久久久大精品| 久久久久国产网址| 观看免费一级毛片| 六月丁香七月| 一二三四中文在线观看免费高清| 成人午夜高清在线视频| 国产亚洲午夜精品一区二区久久 | 精品熟女少妇av免费看| 日日摸夜夜添夜夜添av毛片| 午夜福利视频1000在线观看| 亚洲av一区综合| 狠狠狠狠99中文字幕| 长腿黑丝高跟| 亚洲国产欧洲综合997久久,| 久久久久久久久大av| 色噜噜av男人的天堂激情| 男女下面进入的视频免费午夜| 亚洲激情五月婷婷啪啪| 午夜精品国产一区二区电影 | 人妻少妇偷人精品九色| 18禁在线无遮挡免费观看视频| 日韩强制内射视频| 免费不卡的大黄色大毛片视频在线观看 | 一边亲一边摸免费视频| 男人舔女人下体高潮全视频| 国内少妇人妻偷人精品xxx网站| 亚洲经典国产精华液单| 国产精品不卡视频一区二区| 国产亚洲av片在线观看秒播厂 | ponron亚洲| 免费无遮挡裸体视频| 中国美白少妇内射xxxbb| 最近最新中文字幕大全电影3| 一级av片app| 中文乱码字字幕精品一区二区三区 | 国产免费福利视频在线观看| 搞女人的毛片| 国产又色又爽无遮挡免| 亚洲人成网站高清观看| 国产乱人偷精品视频| 桃色一区二区三区在线观看| 午夜福利视频1000在线观看| 国产免费视频播放在线视频 | 久久久久久久午夜电影| 日本wwww免费看| 国产午夜精品久久久久久一区二区三区| 日本-黄色视频高清免费观看| 女人久久www免费人成看片 | АⅤ资源中文在线天堂| 男女下面进入的视频免费午夜| 男人舔奶头视频| 精品国产露脸久久av麻豆 | 免费看日本二区| 国产精品嫩草影院av在线观看| 丝袜美腿在线中文| 色哟哟·www| 国产成人精品婷婷| 精品久久久久久电影网 | 亚洲在线观看片| 亚洲人与动物交配视频| 亚洲va在线va天堂va国产| 少妇裸体淫交视频免费看高清| 国产一区亚洲一区在线观看| 看非洲黑人一级黄片| 神马国产精品三级电影在线观看| 亚洲国产高清在线一区二区三| 成人高潮视频无遮挡免费网站| 联通29元200g的流量卡| 国产综合懂色| 一卡2卡三卡四卡精品乱码亚洲| 国产黄色视频一区二区在线观看 | 国语自产精品视频在线第100页| 日韩欧美在线乱码| 中文字幕久久专区| 国产综合懂色| 高清毛片免费看| 国产精品伦人一区二区| 午夜激情福利司机影院| av女优亚洲男人天堂| 午夜福利在线观看吧| 日韩成人av中文字幕在线观看| 有码 亚洲区| 免费不卡的大黄色大毛片视频在线观看 | 麻豆一二三区av精品| 观看免费一级毛片| 国产精品永久免费网站| 欧美极品一区二区三区四区| 超碰97精品在线观看| 韩国av在线不卡| 亚洲欧美一区二区三区国产| 国产亚洲av嫩草精品影院| 亚洲高清免费不卡视频| 一本久久精品| 丰满少妇做爰视频| 日本免费a在线| 欧美性猛交╳xxx乱大交人| 国产一区二区亚洲精品在线观看| 九九久久精品国产亚洲av麻豆| 亚洲自拍偷在线| av国产免费在线观看| 久久99蜜桃精品久久| 波多野结衣高清无吗| 日韩欧美 国产精品| 国产一级毛片七仙女欲春2| h日本视频在线播放| 看免费成人av毛片| 免费看光身美女| 国产精品人妻久久久久久| 国产欧美日韩精品一区二区| 丝袜喷水一区| 国产高清三级在线| 99久久成人亚洲精品观看| 亚洲精品一区蜜桃| 国产人妻一区二区三区在| 嫩草影院新地址| 综合色丁香网| 国内揄拍国产精品人妻在线| 国产女主播在线喷水免费视频网站 | 国产视频首页在线观看| 国产一区二区亚洲精品在线观看| 国产高清视频在线观看网站| 国产男人的电影天堂91| 亚洲精品色激情综合| 久久人人爽人人爽人人片va| 又粗又硬又长又爽又黄的视频| 国产成人精品婷婷| 亚洲av福利一区| 亚洲精品色激情综合| 国产成人精品婷婷| 亚洲国产欧美在线一区| 国产精品人妻久久久影院| 国产欧美另类精品又又久久亚洲欧美| 日日撸夜夜添| 99热网站在线观看| 亚洲国产精品专区欧美| 色综合亚洲欧美另类图片| 国产私拍福利视频在线观看| 美女国产视频在线观看| 我的女老师完整版在线观看| 国产精品,欧美在线| www.色视频.com| 色吧在线观看| 国产精品国产三级国产av玫瑰| 最近最新中文字幕免费大全7| 日韩制服骚丝袜av| 免费不卡的大黄色大毛片视频在线观看 | 亚洲在久久综合| 看免费成人av毛片| 女人久久www免费人成看片 | 久久久久久久久久久免费av| 秋霞伦理黄片| 天堂av国产一区二区熟女人妻| 久久这里有精品视频免费| 中文资源天堂在线| 欧美不卡视频在线免费观看| 伦理电影大哥的女人| 国产白丝娇喘喷水9色精品| 国产精品1区2区在线观看.| 色网站视频免费| 寂寞人妻少妇视频99o| 亚洲精品aⅴ在线观看| 欧美精品一区二区大全| 日韩av不卡免费在线播放| 亚洲av中文av极速乱| 亚洲真实伦在线观看| 国产一区二区在线观看日韩| 精品一区二区免费观看| 亚洲18禁久久av| 男女国产视频网站| 综合色丁香网| 卡戴珊不雅视频在线播放| 午夜福利在线观看吧| 免费观看的影片在线观看| 国产淫片久久久久久久久| 日韩一区二区三区影片| 国产一级毛片在线| 国产午夜精品论理片| 99久久九九国产精品国产免费| 人体艺术视频欧美日本| 国产乱人视频| 亚洲一级一片aⅴ在线观看| 日日撸夜夜添| 亚洲av男天堂| 国产精品久久久久久精品电影| 欧美性猛交╳xxx乱大交人| 国产免费视频播放在线视频 | 毛片一级片免费看久久久久| 免费观看性生交大片5| 性插视频无遮挡在线免费观看| 春色校园在线视频观看| 晚上一个人看的免费电影| 乱码一卡2卡4卡精品| 亚洲丝袜综合中文字幕| 亚洲成av人片在线播放无| 中国美白少妇内射xxxbb| 黄色配什么色好看| 午夜福利在线观看免费完整高清在| 午夜福利网站1000一区二区三区| 日韩精品青青久久久久久| 免费观看a级毛片全部| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久精品电影| 国产精品,欧美在线| 久久久久久久午夜电影| 嫩草影院新地址| 日本欧美国产在线视频| 免费黄色在线免费观看| 国产黄色视频一区二区在线观看 | 国产伦一二天堂av在线观看| 亚洲欧美成人精品一区二区| 国产成人午夜福利电影在线观看| 日韩欧美在线乱码| 波野结衣二区三区在线| 国产不卡一卡二| 综合色丁香网| 日本一二三区视频观看| 春色校园在线视频观看| 麻豆成人午夜福利视频| 在线免费观看不下载黄p国产| 国产69精品久久久久777片| 国产精品伦人一区二区| 午夜精品国产一区二区电影 | 日韩国内少妇激情av| 精品国产一区二区三区久久久樱花 | 综合色丁香网| av卡一久久| 国产精品熟女久久久久浪| 国产精品一区www在线观看| 青春草国产在线视频| 国产91av在线免费观看| 真实男女啪啪啪动态图| 欧美另类亚洲清纯唯美| 国产免费又黄又爽又色| 看片在线看免费视频| 国产探花极品一区二区| 免费av不卡在线播放| 永久免费av网站大全| 国产三级在线视频| 国内揄拍国产精品人妻在线| 国产精品人妻久久久影院| 噜噜噜噜噜久久久久久91| 欧美激情在线99| 日本一二三区视频观看| 少妇人妻精品综合一区二区| 亚洲国产成人一精品久久久| 又黄又爽又刺激的免费视频.| 赤兔流量卡办理| 免费看光身美女| 国产午夜福利久久久久久| 日日啪夜夜撸| 非洲黑人性xxxx精品又粗又长| 午夜精品一区二区三区免费看| 国产一级毛片七仙女欲春2| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 国产老妇女一区| 亚洲国产精品专区欧美| 国产在视频线精品| 亚洲天堂国产精品一区在线| 99视频精品全部免费 在线| 99久久中文字幕三级久久日本| 男人的好看免费观看在线视频| 天堂影院成人在线观看| 高清午夜精品一区二区三区| 天堂√8在线中文| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 嫩草影院新地址| 欧美变态另类bdsm刘玥| 少妇的逼好多水| 国产久久久一区二区三区| 亚洲真实伦在线观看| 久久久久网色| 黄色一级大片看看| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 看非洲黑人一级黄片| 亚洲五月天丁香| 久久久久久久久久成人| 乱人视频在线观看| 久久久欧美国产精品| 久久99热这里只有精品18| 国产亚洲一区二区精品| 久久欧美精品欧美久久欧美| 一本一本综合久久| 日本-黄色视频高清免费观看| 91精品伊人久久大香线蕉| 91狼人影院| 有码 亚洲区| 最近最新中文字幕大全电影3| 日韩精品有码人妻一区| 欧美成人精品欧美一级黄| 日韩av在线免费看完整版不卡| 大香蕉97超碰在线| 亚洲国产高清在线一区二区三| 日韩av在线大香蕉| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| 国产亚洲91精品色在线| 国产免费又黄又爽又色| av又黄又爽大尺度在线免费看 | 国产精品国产三级专区第一集| 男人和女人高潮做爰伦理| 伦理电影大哥的女人| 我要搜黄色片| 内射极品少妇av片p| 亚洲精品国产成人久久av| 啦啦啦观看免费观看视频高清| 国产成人免费观看mmmm| 乱码一卡2卡4卡精品| 免费播放大片免费观看视频在线观看 | 99热这里只有精品一区| 看非洲黑人一级黄片| 国产精品一二三区在线看| 日韩一本色道免费dvd| 精品欧美国产一区二区三| 国产老妇女一区| 一个人观看的视频www高清免费观看| 高清视频免费观看一区二区 | 国产乱来视频区| 两个人的视频大全免费| 亚洲国产色片| ponron亚洲| 亚洲欧美日韩高清专用| 我的老师免费观看完整版| 禁无遮挡网站| 午夜福利在线观看免费完整高清在| 午夜福利在线观看吧| 国产精品.久久久| 日本av手机在线免费观看| 国产成人一区二区在线| 综合色丁香网| 久久鲁丝午夜福利片| 免费大片18禁| 天天一区二区日本电影三级| 亚洲欧美日韩卡通动漫| 免费无遮挡裸体视频| 99国产精品一区二区蜜桃av| 国产精品.久久久| 一区二区三区四区激情视频| 大香蕉久久网| 男女视频在线观看网站免费| 波野结衣二区三区在线| 日韩高清综合在线| 在线天堂最新版资源| 国模一区二区三区四区视频| 美女脱内裤让男人舔精品视频| 欧美xxxx黑人xx丫x性爽| 国产白丝娇喘喷水9色精品| 女人十人毛片免费观看3o分钟| 亚洲最大成人中文| av在线播放精品| 不卡视频在线观看欧美| 女人十人毛片免费观看3o分钟| 少妇熟女aⅴ在线视频| 日本av手机在线免费观看| 欧美人与善性xxx| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 久久精品久久精品一区二区三区| 99久久九九国产精品国产免费| 69av精品久久久久久| 久久鲁丝午夜福利片| 国产黄片视频在线免费观看| 欧美不卡视频在线免费观看| 国模一区二区三区四区视频| 久久精品91蜜桃| 99久国产av精品国产电影| 天堂中文最新版在线下载 | 国产欧美另类精品又又久久亚洲欧美| 黄色欧美视频在线观看| 精品久久国产蜜桃| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 人人妻人人看人人澡| 嫩草影院精品99| 国产精品不卡视频一区二区| 午夜亚洲福利在线播放| 三级经典国产精品| 18禁裸乳无遮挡免费网站照片| 欧美区成人在线视频| 性色avwww在线观看| 日本免费在线观看一区| 亚洲国产色片| 亚洲av中文av极速乱| 青春草国产在线视频| 中文欧美无线码| 男人狂女人下面高潮的视频| 男女下面进入的视频免费午夜| 九九爱精品视频在线观看| 婷婷六月久久综合丁香| 十八禁国产超污无遮挡网站| 美女大奶头视频| 精品人妻熟女av久视频| 亚洲国产日韩欧美精品在线观看| 中文字幕av在线有码专区| 日本五十路高清| 国产精品久久久久久精品电影小说 | 久久精品国产鲁丝片午夜精品| av播播在线观看一区| 国产在视频线精品| 亚洲成色77777| 国产精品女同一区二区软件| 国产成人午夜福利电影在线观看| 色吧在线观看| 国产v大片淫在线免费观看| 国产高清不卡午夜福利| 一夜夜www| 亚洲欧美成人精品一区二区| 变态另类丝袜制服| 国产成人精品一,二区| 久久久成人免费电影| 国产精品久久电影中文字幕| 91aial.com中文字幕在线观看| 国产精品一区www在线观看| 国内精品美女久久久久久| 一区二区三区乱码不卡18| 日日摸夜夜添夜夜添av毛片| АⅤ资源中文在线天堂| 国产 一区 欧美 日韩| 黄片wwwwww| 伊人久久精品亚洲午夜| 亚洲内射少妇av| 尤物成人国产欧美一区二区三区| 26uuu在线亚洲综合色| 人妻少妇偷人精品九色| 99久久精品一区二区三区| 国产免费一级a男人的天堂| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 日韩三级伦理在线观看| 成年av动漫网址| 中文欧美无线码| 精品人妻熟女av久视频| 日本黄色视频三级网站网址| 视频中文字幕在线观看| 免费一级毛片在线播放高清视频| 日韩欧美三级三区| 亚洲av电影不卡..在线观看| 人人妻人人澡欧美一区二区| 水蜜桃什么品种好| 日本黄色视频三级网站网址| 国产精品久久久久久av不卡| 91狼人影院| 久久久久久大精品| 亚洲最大成人中文| 插阴视频在线观看视频| 日本免费一区二区三区高清不卡| 久久久国产成人精品二区| 亚洲av成人精品一二三区| 最近最新中文字幕大全电影3| 一边亲一边摸免费视频| 国产精品99久久久久久久久| av国产久精品久网站免费入址| 精品国产一区二区三区久久久樱花 | 久久久久久久国产电影| 亚洲精品色激情综合| 中文字幕av在线有码专区| 校园人妻丝袜中文字幕| 全区人妻精品视频| 在线免费观看不下载黄p国产| 免费观看a级毛片全部| 国产在线一区二区三区精 | 在线播放国产精品三级| 高清av免费在线| 亚洲无线观看免费| av国产久精品久网站免费入址| 三级国产精品欧美在线观看| 啦啦啦啦在线视频资源| 美女国产视频在线观看| 赤兔流量卡办理| 69人妻影院| 九九在线视频观看精品| 麻豆成人午夜福利视频| 深夜a级毛片| 久久精品国产鲁丝片午夜精品| 又粗又硬又长又爽又黄的视频| 国产精品不卡视频一区二区| 亚洲五月天丁香| 成人毛片60女人毛片免费| 久久99热6这里只有精品| 久久久久精品久久久久真实原创| 国产色婷婷99| 免费观看人在逋| 在线a可以看的网站| 熟女人妻精品中文字幕| 婷婷色综合大香蕉| 真实男女啪啪啪动态图| 一级二级三级毛片免费看| 国产精品1区2区在线观看.| 男人和女人高潮做爰伦理| av在线观看视频网站免费| 高清毛片免费看| АⅤ资源中文在线天堂| 日韩国内少妇激情av| 男人的好看免费观看在线视频| 日韩av在线免费看完整版不卡| 亚洲av不卡在线观看| 亚洲精品乱久久久久久| 色5月婷婷丁香| 午夜福利在线观看吧| 天天躁夜夜躁狠狠久久av| 国产精品蜜桃在线观看| 精品午夜福利在线看| 欧美日韩国产亚洲二区| 日韩高清综合在线| 日产精品乱码卡一卡2卡三| 蜜桃久久精品国产亚洲av| 欧美精品国产亚洲| 国产探花极品一区二区| 国产精品熟女久久久久浪| 又粗又爽又猛毛片免费看| 秋霞在线观看毛片| 亚洲中文字幕一区二区三区有码在线看| 精品人妻视频免费看| 最近的中文字幕免费完整| 亚洲成人av在线免费| 午夜免费激情av| h日本视频在线播放| 嫩草影院入口| 白带黄色成豆腐渣| 国产免费视频播放在线视频 | 久久久a久久爽久久v久久| 日韩大片免费观看网站 | 黄色日韩在线| 91精品伊人久久大香线蕉| 白带黄色成豆腐渣| 丝袜美腿在线中文| 最近中文字幕2019免费版| 最近手机中文字幕大全| 全区人妻精品视频| 亚洲av免费高清在线观看| 亚洲av日韩在线播放| 男人和女人高潮做爰伦理| 在线a可以看的网站| 午夜老司机福利剧场| 国产老妇女一区| 99久国产av精品国产电影| 亚洲精品自拍成人| 男人狂女人下面高潮的视频| www日本黄色视频网| 国产乱人偷精品视频| 久久久久久久久久黄片| 精品久久久噜噜| 久久久精品94久久精品| 联通29元200g的流量卡| 亚洲精品乱久久久久久| 国产一区有黄有色的免费视频 | 麻豆久久精品国产亚洲av| 国产精品一区www在线观看| 亚洲精品久久久久久婷婷小说 | 国产精品99久久久久久久久| 久久99热这里只频精品6学生 | 爱豆传媒免费全集在线观看| 免费观看精品视频网站| 国产亚洲最大av| 欧美人与善性xxx| 少妇猛男粗大的猛烈进出视频 | 亚洲美女搞黄在线观看| 最新中文字幕久久久久| 色尼玛亚洲综合影院| 99热网站在线观看| 免费看日本二区| 热99在线观看视频| 一级爰片在线观看| 床上黄色一级片| 99久久精品一区二区三区| 国产精品日韩av在线免费观看| 色5月婷婷丁香| 两个人视频免费观看高清| 三级经典国产精品| 久久午夜福利片| 小说图片视频综合网站| 国产真实乱freesex| 日日啪夜夜撸| 久久久国产成人精品二区| 最新中文字幕久久久久| 国产精品伦人一区二区| 免费播放大片免费观看视频在线观看 | 99久国产av精品| 午夜a级毛片| 国产在线一区二区三区精 | 亚洲国产精品成人综合色| 久久久久久久久久久免费av| 久久久久性生活片| 国产私拍福利视频在线观看| 内射极品少妇av片p| 精品人妻熟女av久视频| 久久久久久久久大av| 亚洲,欧美,日韩| 亚洲国产色片| 久久久成人免费电影| 欧美激情国产日韩精品一区| 乱系列少妇在线播放| 亚洲av.av天堂| 日本熟妇午夜| 国产中年淑女户外野战色| 最近最新中文字幕免费大全7| 日本午夜av视频| 国产精品野战在线观看| 成人无遮挡网站| 91av网一区二区| 国产午夜精品一二区理论片| 国产精华一区二区三区| 色噜噜av男人的天堂激情| 亚洲欧美成人精品一区二区| 一本一本综合久久| 草草在线视频免费看| 国产成人一区二区在线|