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

    Quantitative Secretome Analysis Reveals Clinical Values of Carbonic Anhydrase II in Hepatocellular Carcinoma

    2021-12-03 09:02:02XiaohuaXingHuiYuanHongzhiLiuXionghongTanBixingZhao
    Genomics,Proteomics & Bioinformatics 2021年1期

    Xiaohua Xing, Hui Yuan, Hongzhi Liu, Xionghong Tan, Bixing Zhao,

    Yingchao Wang2, Jiahe Ouyang2, Minjie Lin2, Xiaolong Liu1,2,*, Aimin Huang1,*

    1 Department of Pathology and Institute of Oncology,School of Basic Medical Sciences of Fujian Medical University,Fuzhou 350004,China

    2 The United Innovation of Mengchao Hepatobiliary Technology Key Laboratory of Fujian Province, Mengchao Hepatobiliary Hospital of Fujian Medical University, Fuzhou 350025, China

    KEYWORDS

    Abstract Early detection and intervention are key strategies to reduce mortality, increase longterm survival, and improve the therapeutic effects of hepatocellular carcinoma (HCC) patients.Herein, the isobaric tag for relative and absolute quantitation (iTRAQ)-based quantitative proteomic strategy was used to study the secretomes in conditioned media from HCC cancerous tissues,surrounding noncancerous tissues,and distal noncancerous tissues to identify diagnostic and prognostic biomarkers for HCC. In total, 22 and 49 dysregulated secretory proteins were identified in the cancerous and surrounding noncancerous tissues, respectively, compared with the distal noncancerous tissues. Among these proteins, carbonic anhydrase II (CA2) was identified to be significantly upregulated in the secretome of cancerous tissues; correspondingly, the serum concentrations of CA2 were remarkably increased in HCC patients compared with that in normal populations.Interestingly,a significant increase of serum CA2 in recurrent HCC patients after radical resection was also confirmed compared with HCC patients without recurrence, and the serum level of CA2 could act as an independent prognostic factor for time to recurrence and overall survival. Regarding the mechanism, the secreted CA2 enhances the migration and invasion of HCC cells by activating the epithelial mesenchymal transition pathway.Taken together,this study identified a novel biomarker for HCC diagnosis and prognosis,and provided a valuable resource of HCC secretome for investigating serological biomarkers.

    Introduction

    Hepatocellular carcinoma (HCC) has a high incidence and mortality,making it the sixth most common malignant cancer worldwide[1].Although surgical strategy was proven to be the most suitable option for HCC therapy, most HCC patients were not diagnosed or intervened until the advanced stage,rendering them unsuitable for surgical treatments[2]or resulting in a poor prognosis after surgical excision. Currently, the five-year survival rate of HCC patients is less than 20%,while the five-year recurrence and metastasis rates are even higher than 80% [3–6]. So far, alpha fetoprotein (AFP) and desgamma-carboxy prothrombin (DCP) are the most widely accepted and clinically applied biomarkers for HCC diagnosis and monitoring. However, the sensitivity and specificity of AFP and DCP in the early diagnosis as well as prognosis evaluation of HCC are insufficient [7–9]. Therefore, it is urgently required to find novel biomarkers that are highly specific and sensitive for early diagnosis and prognosis evaluation of HCC.

    Proteomics-based technology has become a very useful and powerful analytical tool for biomarker screening [10–13]. A desirable biomarker for HCC diagnosis or monitoring should allow measurement using body fluid samples such as serum and plasma [14], since these samples are low-cost and easy to collect and process, and are amenable to repeat sampling whenever it is necessary. Therefore, the serum and plasma are also the ideal targets for proteomic studies that aim to identify diagnostic or prognostic biomarkers for HCC[15,16]. However, the complex nature of serum and plasma,as well as their large dynamic concentration range of different proteins, significantly hinders the progress of proteomicsbased biomarker screening.

    Secretory proteins play important roles in signal transduction, cell growth, proliferation, and apoptosis, and even in tumorigenesis, development, invasion, and metastasis, and are ideal sources for biomarker screening [17]. Investigating the secretome of HCC tissues or cells may provide valuable information for identifying potential biomarkers. Recently,the application of secretomics in screening diagnostic or prognostic protein biomarkers in HCC cell lines has been reported by many groups [18–20]. However, these results still must be clinically validated [21]. Therefore, it would be more straightforward and convincing to analyze the secretome of primary tumor tissue cultures to identify the diagnostic or prognostic biomarkers for HCC.For example,Yang et al.[22]have established an in vitro tissue culture system for HCC and identified matrix metalloproteinase 1(MMP1)as a diagnostic biomarker for HCC; however, the influence of hepatitis B virus (HBV)infection was not analyzed in this study.

    In the present study, we collected serum-free conditioned media (CM) from the tissue cultures of cancerous (C) tissues,surrounding noncancerous (SN)tissues, and distal noncancerous (DN) tissues of HCC patients, respectively, and analyzed their secretomes to identify potential diagnostic and prognostic biomarkers for HCC via an isobaric tag for relative and absolute quantitation (iTRAQ)-based quantitative proteomic approach. Meanwhile, the sensitivity, specificity, and clinical significance of the identified biomarkers were also carefully validated in a large-scale HCC patient cohort by enzyme linked immunosorbent assay(ELISA)and targeted proteomics of parallel reaction monitoring (PRM). Furthermore, the molecular mechanisms of the identified biomarker, carbonic anhydrase II (CA2), were also carefully explored.

    Results and discussion

    Cells in the in vitro-cultured tissues were alive and secretory

    To ensure that the cells in the in vitro-cultured tissues were alive and the secretome was not contaminated by intracellular proteins, a series of analyses were performed (Figure 1A). We used haematoxylin and eosin (H&E) staining to evaluate the cell morphology changes of tissues cultured for 0 day, 1 day,and 2 days. As revealed in Figure 1B, the H&E stained tissue sections showed corresponding characteristic anatomical details of HCC C, SN, and DN tissues in all respective cultures. With the extension of incubation time, cells in the cultured tissues were starved and showed necrosis or apoptosis due to lack of nutrients. The morphology of tissues cultured for 1 day was still very similar to that of the fresh tissues (cultured for 0 day). By contrast, the number of cell nuclei in tissues cultured for 2 days significantly decreased compared with the fresh tissues due to cell necrosis or apoptosis during the culture process. TdT-mediated dUTP nick-end labeling(TUNEL) staining was further used to evaluate the apoptosis rates of tissues cultured for 0 day, 1 day, and 2 days. As revealed in Figure 1C,the tissues cultured for 2 days had a significantly higher apoptosis rate compared with tissues cultured for 0 day or 1 day.In addition,the proteins extracted from culture supernatants with different culture times were examined by SDS-PAGE. As shown in Figure 1D,the molecular weight distribution of the extracted proteins was significantly changed with the increase of culture time, especially when the culture time was more than 1 day. Furthermore, Western blot assays clearly demonstrated the prevention of contamination by intracellular proteins(Figure 1E).Taken together,these results suggested that the 1-day culture time was the optimal time point for collecting culture supernatants.

    The overall characteristics and Genetic Ontology enrichment of the secretory proteins identified from tissue secretomes

    iTRAQ labeling combined with mass spectrometry was applied to investigate the secretome difference of the three types of tissues from primary HCC patients. The features of the ten HCC patients in current study are listed in Table S1.Total proteins separately extracted from the supernatants collected from the CM of HCC C,SN,and DN tissues were analyzed using two-dimensional liquid chromatography–tandem mass spectrometry (2D LC–MS/MS), as shown in Figure 2.We quantified 2388 proteins in total using Scaffold_4.3.2, of which 1312 proteins were annotated or predicted as secretory proteins. This result covered 75.7% of the secretory proteins previously reported by Yang and colleagues[22](Figure S1A),and the percentage of secretory proteins in our dataset(54.9%;1312/2388) was much higher than that in the human protein database (23%) (Figure S1B). The complete list of quantified secretory proteins, as well as their detailed features including the NN-score, MetazSecKB characteristic, isoelectric point(pI), molecular weight, hydrophobicity, and quantification results, is shown in Table S2. Among these secretory proteins,936 were commonly found in 5 biological replicates, which accounted for 71.3% of the quantified secretory proteins.

    Figure 1 Strict quality control of the secretory proteins from the cultured HCC-related tissuesA. The quality control workflow for the cultured HCC-related tissues and the secretory proteins in the supernatant. B. H&E staining of cultured HCC-related tissues.Scale bar,1 mm.C.TUNEL staining of cultured HCC-related tissues.We examined the densities of DAPI(blue) and FITC (green) for each tissue. Scale bar, 150 μm. D. The molecular weight distribution of secretory proteins determined by SDS-PAGE.E.Western blot demonstrating the effective prevention of intracellular protein contamination in extracted secretory proteins.Actin was used as a control. HCC, hepatocellular carcinoma; C, cancerous tissue; SN, surrounding noncancerous tissue; DN, distal noncancerous tissue; CM, conditioned medium; H&E, haematoxylin and eosin; TUNEL, TdT-mediated dUTP nick-end labeling.

    The molecular weight distribution of the secretory proteins ranged from 7396 Da to 628,685 Da, with a primary range of 10–40 kDa, indicating smaller molecular weight of the secretory proteins (Figure 3A). The pI values ranged from 3.67 to 12.56 and were mainly in the range of 4.4–7.6,which coincided with the microenvironment of liver tissues (Figure 3B). The hydrophobicity of the secretory proteins ranged from 1.2 to 3.3 and was mainly in the range of 1.8–2.2, which implied the enrichment of membrane or transmembrane proteins(Figure S1C).

    We also investigated the Genetic Ontology (GO) annotations containing cell components, molecular functions, and biological processes of the secretory proteins. Cell component category showed that the secretory proteins were primarily extracellular, which suggested excellent purity of the HCCrelated tissue secretomes (Figure 3C). The biological process category indicated that these secretory proteins were mainly involved in biosynthetic process, signal transduction, and transport process (Figure 3D). The molecular function category indicated that these secretory proteins played major roles in RNA binding, enzyme binding, DNA binding, and transmembrane transporter activity (Figure 3E).These results were consistent with those of previously reported secretory proteins from HCC tissues [22].

    Figure 2 Experimental workflow for the secretome analysis of HCC-related tissuesThe three types of tissues (C, SN, and DN) from HCC patients were cultured in vitro, and the CM was collected to extract secretory proteins. Secretory proteins were digested with trypsin, directly labeled using iTRAQ-8plex, and analyzed through 2D LC–MS/MS. The target proteins screened by bioinformatics were then verified in vitro to find potential biomarkers of HCC and to investigate the molecular mechanisms of HCC recurrence. iTRAQ, isobaric tag for relative and absolute quantitation; 2D LC–MS/MS, two-dimensional liquid chromatography–tandem mass spectrometry; ELISA, enzyme linked immunosorbent assay; PRM, parallel reaction monitoring.

    Various functions of secretory proteins in different HCC-related tissues

    Next, we identified the differentially abundant secretory proteins in C and SN tissues by comparing with the DN tissue,which should have the same change tendency in all 10 biological replicates and presented fold change>1.5 or<0.67 in at least 5 biological replicates.Under this standard,there were 22 differentially abundant secretory proteins in the C/DN group(Table S3), and 49 differentially abundant secretory proteins in the SN/DN group (Table S4). Among these differentially abundant secretory proteins, 13 were shared between the two comparison groups(Figure 3F).GO biological process enrichment analysis indicated that these shared proteins mainly participated in extracellular matrix organization, extracellular structure organization,and tissue morphogenesis(Figure S1D),which might be contributing factors for HCC development.In addition,9 differentially abundant secretory proteins that were only found in the C/DN group mainly participated in the disorder of primary metabolism (Figure S1E), which is tightly linked to the development, and even invasion and metastasis of HCC. Meanwhile, 36 differentially abundant secretory proteins that were only found in the SN/DN group mainly participated in acute-phase response, acute inflammatory response, and post-transcriptional regulation on gene expression(Figure S1F),suggesting that the SN tissues are distinctly different from DN tissues. Interestingly, these processes were all related to primary metabolism, suggesting that changes in primary metabolism might play a crucial role in the occurrence of HCC.

    Figure 3 The overall characteristics and GO enrichment analyses for the secretory proteins identified from HCC tissue secretomesA. and B. The distribution of molecular weight (A) and isoelectric point (B) of the identified secretory proteins. C.–E. GO enrichment analyses of the cell components (C), biological processes (D), and molecular functions (E) of the secretory proteins. F. Venn diagrams showing the numbers of differentially abundant secretory proteins in the two comparison groups (the C/DN group and the SN/DN group).G.and H.The signaling pathway networks in which the differentially abundant secretory proteins of the C/DN group(G)and the SN/DN group (H) involved. GO, Genetic Ontology.

    To further study the potential molecular mechanisms of the occurrence and development of HCC, Ingenuity Pathways Analysis (IPA) was applied to investigate the signaling pathways in which the differentially abundant secretory proteins participated. The results showed that the two comparison groups (the C/DN group and the SN/DN group) indeed had unique dysregulated signaling pathways, although they also had shared signaling pathways.As shown by IPA,the differentially abundant secretory proteins in the C/DN group were mostly involved in ERK/MAPK signaling, while the differentially abundant secretory proteins in the SN/DN group mainly participated in PI3K-Akt signaling (Figure 3G and H). There were 11 differentially abundant secretory proteins (2 upregulated and 9 down-regulated)participating in ERK/MAPK signaling (Figure 3G). The ERK/MAPK pathway can transduce extracellular signals through intracellular signal transduction cascades to control the expression of proteins that regulate tumorigenesis and aggressive behaviors [23,24].The dysregulation of ERK/MAPK signaling pathway in the C/DN group revealed that the secreted factors from HCC C tissue might modulate tumor microenvironment to exert important roles in tumorigenesis.In the SN/DN group,20 differentially abundant secretory proteins (4 up-regulated and 16 down-regulated) participated in PI3K-Akt signaling(Figure 3H).PI3K-Akt signaling plays a crucial role in the regulation of inflammation and metabolism[25–27].Hence,alterations in the PI3K-Akt pathway might be closely linked to the occurrence and development of tumors. The dysregulation of PI3K-Akt pathway in the secretory environment of HCC SN tissues compared with that of DN tissues suggested the importance of changing microenvironment in tumorigenesis.Targeting these important effectors in tumor microenvironment might be a promising therapeutic strategy.

    CA2 might be a valuable biomarker for HCC diagnosis

    According to the IPA results, we found that CA2 was upregulated in both the C/DN and SN/DN groups. CA2, a zinc metal enzyme, carries out the reversible hydration of carbon dioxide, and plays a key role in adjusting the pH of tumor microenvironment [28–33]. It is frequently abnormally expressed in different cancers[34–40].Here,the relative intensity of the iTRAQ reporter ions related to CA2 in the MS/MS spectra was in good agreement with the protein level of CA2. As shown in Figure S2,the relative intensities of the iTRAQ reporter ions related to CA2 were remarkably increased in C and SN tissues compared with those in DN tissues(iTRAQ-114> iTR AQ-115 > iTRAQ-116 and iTRAQ-117 > iTRAQ-118 >iTRAQ-119). Therefore, CA2 might be a potential interesting biomarker for HCC diagnosis and prognosis prediction.

    Next, we verified the clinical significance of CA2 on HCC diagnosis and prognosis in two additional serum cohorts via the PRM targeted proteomics method and traditional ELISA.The CA2 serum concentration was monitored in 49 HCC patients and 23 healthy volunteers through PRM targeted proteomics. A total of nine peptides from CA2 were identified,seven of which were commonly found in the five iTRAQ-8plex experiments with one carrying a miscleavage site(Table S5). Therefore, only six peptides from CA2 could be used for PRM analysis. The annotated spectra and detailed peptide spectrum match (PSM) information for the six identified CA2 peptides were provided in Table S5 and Figure S3,respectively.Using Skyline,the peak contributions of the individual fragment ions from the unique peptide were determined(Figure S4A), and the representative quantification information based on peak areas of the peptides including endogenic and synthetic heavy peptides is displayed in Figure S4B. The average serum concentration of CA2 in HCC patients was remarkably higher than that in healthy volunteers (HCC patients 13.72 pg/ml vs. healthy volunteers 7.01 pg/ml;P < 0.01). The receiver operating characteristic (ROC) curve analysis of CA2 revealed that the area under the curve(AUC)was 0.715 for HCC patients relative to the healthy volunteers(Figure 4A).These results were consistent with the data from the proteomic studies, indicating that CA2 might be a valuable biomarker for HCC diagnosis.

    CA2 is a novel prognostic biomarker for HCC

    By analyzing our previous protemoics data of a series of serum samples collected at different time points from four HCC patients after radical resection [41], we found that the serum levels of CA2 in the patients with short-term recurrence were steadily increased from 3 months to 9 months after radical resection compared with the patients with long-term recurrence (Figure S5A). Thus, we cellected the serum samples of 159 HCC patients, containing 94 relapsed patients and 65 relapse-free patients, at 5 months after surgery to analyze the association between CA2 serum concentration and HCC prognosis by ELISA. As shown in Figure 4B, the average serum concentration of CA2 was remarkably upregulated in HCC patients with recurrence compared with HCC patients without recurrence (P < 0.0001). According to the ROC curve, the AUC reached 0.708 in the validation cohort, indicating that CA2 might be a HCC prognostic biomarker. To evaluate the prognostic significance of serum CA2 in HCC,we further studied the cumulative recurrence rate and overall survival rate of patients with low and high serum levels of CA2 using the 159 samples.The median concentration(105.3 ng/mL)was used as the optimal cutoff value to stratify patients into low-CA2 (≤105.3 ng/ml)and high-CA2(>105.3 ng/ml)groups.As shown by the Kaplan-Meier analysis,HCC patients with a high serum CA2 level had a higher recurrence rate than those with a low serum CA2 level. Moreover, the overall survival rate of HCC patients who had a high serum CA2 level was remarkably lower than those with a low CA2 level (Figure 4C). Furthermore, we repeated this study using a training cohort of 49 HCC patients and obtained similar results (Figure S5B). In addition, we found a strong negative linear correlation between CA2 serum concentration and recurrence time, i.e.,the higher the CA2 serum concentration,the shorter the recurrence time(Figure S5C). These results suggested that the CA2 serum concentration might be a novel prognostic biomarker for HCC.

    CA2 is a novel prognostic biomarker for AFP-negative HCC patients

    As the gold standard for clinical diagnosis and monitoring in HCC, the average serum concentration of AFP was remarkably higher in relapsed HCC patients than that in relapse-free patients (P < 0.001) (Figure S5D). Although the cumulative recurrence rate was substantially lower and the overall survival rate was remarkably higher in AFPnegative patients(AFP serum concentration<20 ng/ml)than those in AFP-positive patients (AFP serum concentration≥ 20 ng/ml) in general, some AFP-negative patients still experienced rapid recurrence or metastasis(Figure S5E).However, there was no significant difference in the average serum concentration of AFP between the relapsed and relapse-free AFP-negative HCC patients (Figure S5F). Therefore, there is an urgent need to identify novel biomarkers that can predict the recurrence/metastasis in AFP-negative HCC patients.Interestingly, among the AFP-negative HCC patients, the average serum concentration of CA2 in relapsed patients was remarkably higher than that in relapse-free patients(P < 0.001) (Figure 4D). Furthermore, in the AFP-negative HCC patients, the cumulative recurrence rate was also higher in patients with a high serum CA2 level, suggesting the prognostic value of CA2 in AFP-negative patients (Figure 4E).However,the overall survival rate did not show significant differences between the high- and low-CA2 groups in the AFPnegative HCC patients (Figure S5G).

    Figure 4 Clinical values of serum CA2 in HCC diagnosis and prognosisA.Distribution(upper panel)and ROC curve(lower panel)of serum CA2 levels in HCC patients and healthy volunteers from the training cohort through PRM targeted proteomics. B. Distribution (upper panel) and ROC curve (lower panel) of serum CA2 levels in HCC patients, with and without recurrence, from the validation cohort assessed by ELISA. C. Comparison of cumulative recurrence rate and overall survival rate between the high-and low-CA2 groups in the validation cohort by Kaplan-Meier analysis.D.Distribution of serum CA2 levels in AFP-negative HCC patients.E.Comparison of cumulative recurrence rate between the high-and low-CA2 groups in AFPnegative HCC patients by Kaplan-Meier analysis. F. AUC of the CA2/AFP combination in HCC patients. G. Prognostic value of the CA2/AFP combination in HCC patients. Two-tailed unpaired Student’s t-test (**, P < 0.01; ***, P < 0.001; ****, P < 0.0001). CA2,carbonic anhydrase II; ROC, receiver operating characteristic; AUC, area under the curve; CI, confidence interval; AFP, alpha fetoprotein.

    Combination of serum CA2 and AFP improves prognostic performance in HCC

    Considering the prognostic value of CA2 in AFP-negative patients, we then jointly considered the serum levels of CA2 and AFP. As shown in Figure 4F, when serum levels of CA2 and AFP were considered together, the AUC reached to 0.803 for the combination (AUC of 0.708 for CA2 and of 0.765 for AFP),suggesting the clinical values for the combination of these two markers. Next, we evaluated the prognostic value of the CA2/AFP combination for HCC patients in the validation cohort. As shown in Figure 4G, the cumulative recurrence rate of patients in the AFP-/low-CA2 group was remarkably lower than that of patients in the AFP+/high-CA2 group.Correspondingly,the overall survival rate of patients in the AFP-/low-CA2 group was substantially higher than that of patients in the AFP+/high-CA2 group.Patients of the AFP-/high-CA2 and AFP+/low-CA2 groups had moderate cumulative recurrence rate and overall survival rate in the validation cohort. Taken together, these results suggested an improved prognostic value when using the serum levels of CA2 and AFP in combination for HCC patients.

    Serum CA2 significantly correlates with tumor size and microvascular invasion

    To further study molecular mechanisms of CA2 on the prognosis of HCC, we checked the association of serum CA2 with various clinicopathologic features of HCC patients. As revealed in Table S6, the pearson’s Chi-square test indicated that higher serum CA2 level was significantly correlated with larger tumor size (P = 0.006), more advanced Barcelona Clinic Liver Cancer (BCLC) stage (P = 0.011), and higher post-operation level of AFP (P = 0.004), as well as microvascular invasion (MVI) (P = 0.008) and cirrhosis (P = 0.003),but was not associated with other clinicopathologic features.These results suggested that CA2 might be a metastasis/recurrence-related protein in HCC.

    Secreted CA2 increases the migration and invasion abilities of HCC cells by activating epithelial-mesenchymal transition signaling

    Figure 5 Extracellular CA2 promotes cell migration and invasion through activating EMTA. Representative images and quantification results of migration and invasion of MHCC97L cells treated with exogenous recombinant CA2. B. Cellular morphology of MHCC97L cells treated with exogenous recombinant CA2. Scale bar, 50 μm. C. Western blot showing the expression levels of E-cadherin, N-cadherin, and Zeb1 in MHCC97L cells treated with exogenous recombinant CA2. Two-tailed unpaired Student’s t-test (*, P < 0.05; **, P < 0.01).

    Figure 6 Expression pattern of intracellular CA2 is opposite to its extracellular formA.Western blot showing the intracellular CA2 expression in C tissues and their paired SN tissues in HCC patients.B.Quantification of the intracellular CA2 expression shown in (A). C. Statistical analysis of the intracellular CA2 expression in C tissues and their paired SN tissues on in-house TMAs determined by IHC.D.Representative images of CA2 IHC in a C tissue(score 1)and its paired SN tissue(score 3). Two-tailed paired Student’s t-test (****, P < 0.0001). TMA, tissue microarray; IHC, immunohistochemistry.

    To investigate the molecular mechanism of extracellular CA2 in HCC recurrence/metastasis, exogenous recombinant CA2 protein (obtained from Abcam) was used to examine its influences on the migration and invasion of HCC cell line MHCC97L using a trans-well strategy.As shown in Figure 5A,exogenous recombinant CA2 treatment significantly promoted the migration and invasion of MHCC97L cells compared to the control group exposed to the same concentration of BSA(P < 0.05), suggesting that extracellular CA2 enhances cell migration and invasion in HCC, which could explain the observed clinical data. Epithelial-mesenchymal transition(EMT) is well known to be responsible for tumor metastasis.As shown in Figure 5B, the MHCC97L cells treated with exogenous recombinant CA2 exhibited a spindle-like fibroblastic morphology,while the control cells were round with a more epithelial morphology. Such morphological changes indicated that secreted CA2 might be involved in the EMT process. To further investigate the involved molecular mechanisms, the related key markers of EMT were also analyzed. As shown in Figure 5C, the addition of exogenous recombinant CA2 up-regulated N-cadherin (a mesenchymal marker) and downregulated E-cadherin (an epithelial maker); meanwhile, the expression level of Zeb1 (EMT-promoting transcription factor) was also significantly up-regulated by the addition of exogenous recombinant CA2. Taken together, these results demonstrated that secreted CA2 promoted EMT to further modulate the migration and invasion of HCC cells, in turn affecting HCC metastasis.

    Intracellular CA2 may perform opposite functions from its extracellular form

    We also analyzed the intracellular levels of CA2 in HCC patients. In this experiment, 28 pairs of C and SN tissues of HCC patients were detected by Western blot. As shown in Figure 6A and B,the intracellular CA2 expression was remarkably decreased in C tissues when compared with their paired SN tissues. Meanwhile, the down-regulated expression of intracellular CA2 in C tissues was further confirmed by immunohistochemistry (IHC) using in-house tissue microarrays (TMAs) containing 75 C tissues and their paired SN tissues (Figure 6C and D, Figure S6). These results were exactly opposite from those of the extracellular CA2 in the serum of HCC patients,suggesting that intracellular and extracellular CA2 might carry out opposite functions [42] and the involved molecular mechanisms should be further explored.

    Conclusion

    Herein, we applied the iTRAQ-based quantitative proteomic approach to investigate the secretome of primary cultures of HCC tissues, and identified secreted CA2 as a diagnostic and prognostic biomarker for HCC. In particular, CA2 showed good predicative performance in AFP-negative HCC patients,and the combination of CA2 and AFP improved the sensitivity and specificity of HCC prognosis. Regarding the mechanism,extracellular CA2 might regulate HCC cell migration and invasion by targeting the EMT signaling pathway to affect HCC patients’ prognosis. This secretome investigation enabled us to identify a novel HCC diagnostic and prognostic biomarker.The information from this study constitutes a valuable resource for further HCC investigation and for identifying potential serological biomarkers of HCC.

    Materials and methods

    Patients and follow-up

    In total, 293 HCC patients and 23 healthy volunteers were enrolled in the current study.All the HCC patients underwent surgical procedures at the Mengchao Hepatobiliary Hospital of Fujian Medical University (Fuzhou, China). Computer tomography (CT) scanning and ultrasonography (US) were used to monitor the absence of intrahepatic recurrence and metastasis in the residual liver. All patients met the enrolment eligibility criteria as follows: 1) The patients were diagnosed with HCC through pathological examination after operation.2) Serum hepatitis B surface antigen (HBs Ag) was positive,and hepatitis B surface antibody (HBs Ab) and hepatitis C virus(HCV)were negative before surgery.3)Standard radical resection was performed:no distal metastasis was found before or during the operation; intraoperative US examination revealed no other liver lesions; no obvious tumor thrombus was found in the hepatic portal vein or primary venous branch;postoperative pathological examination showed no cancer cells at the cutting edge;and no recurrence/metastasis was found in the US and CT examination two months after the operation.4)Serum AFP was increased before operation and decreased to normal level two months after surgical operation. 5) The patient had not received any other intervention or treatment before surgery.

    For the proteomic analysis, C, SN, and DN tissues were obtained from ten HCC patients who underwent surgical operation. Serum from 49 HCC patients before hepatectomy and 23 healthy volunteers in a health screening program at the Mengchao Hepatobiliary Hospital of Fujian Medical University were collected for PRM analysis. Another validation serum cohort from 159 HCC patients after hepatectomy with long-term follow-up that ended with death was used for ELISA.The collection of serum samples was strictly controlled following the in-house standard operating procedure, which was established on the basis of the previous study[43].In addition,a total of 75 formalin-fixed and paraffin-embedded C tissues and their paired SN tissues from HCC patients who underwent surgical resection were collected to generate the TMAs for IHC investigation.

    In vitro tissue culture and quality control

    Following surgical operation, the primary tissue specimens were immediately transferred into PBS on ice and sent to the laboratory within 30 min. The tissues were then rinsed and cut into 2–3 mm3pieces, and extensively washed several times at room temperature (RT) with PBS to eliminate major blood and serum contaminants. Subsequently, the samples were transferred to 10-cm cell culture dishes and incubated in serum-free DMEM medium supplemented with 1% penicillin and streptomycin at 37 °C. For protein extraction, the supernatants were centrifuged at 100 g for 2 min at 4 °C to further remove all remaining cells and debris; the samples were then concentrated with a 3 K cutoff centrifugal filter device (Catalog No. UFC500396, Millipore, Billerica, MA) and stored at-80 °C until further analysis.

    Culture medium was collected every 4 h in following 72-h period (a total of 18 time points). The proteins extracted from every time point were subjected to 10%SDS-PAGE to analyze the molecular weight distribution, and Western blot was performed to determine any contamination by intracellular proteins. At the same time, the cultured tissues collected at every time point were evaluated by histological observations.The cultured tissues were fixed in 10% neutral formaldehyde,and paraffin sections were made in the conventional way.One slide was stained with H&E (Catalog No. D006, Nanjing Jiancheng Bioengineering Institute, Nanjing, China) and observed under an optical microscope (AE2000, Motic, Xiamen, China); another slide was biotin-labeled and stained by TUNEL to detect cell apoptosis during the culture period.Finally, the slides were observed under a confocal microscope(LSM 780, Carl Zeiss, Jena, Germany).

    Bottom-up proteomics and data analysis

    The proteomic studies and data analyses were modified from our previously reported protocols[44].Briefly,the thiol groups of the proteins extracted from the aforementioned culture supernatants of three groups (C, SN, and DN) were reducted with 8 mM DTT at 55 °C for 1 h, and alkylated with 50 mM iodoacetamide for 30 min in dark at RT. Then, the proteins were precipitated by ice-cold acetone (5× volume),and re-dissolved in 100 mM tetraethyl-ammonium bromide(TEAB). Subsequently, 100 μg of each protein sample was digested by trypsin (Catalog No. V511, Promega, Madison,WI) using filter aided sample preparation (FASP), and the peptides were labeled with the 8-plex iTRAQ reagent(Catalog No. 4381663, AB SCIEX, Boston, MA) as follows: C group,SN group, and DN group were labeled with 114, 115, and 116,respectively;and one biological replicate of the aforementioned three groups was labeled with 117,118,and 119,respectively. A, B, C, D, and E were defined as the five independent replicates of iTRAQ 8-plex labeling. In addition, ten DN group samples mixed in equal amounts were labeled with 113 and included in every 8-plex labeling reaction as an internal standard to balance each 8-plex labeling. The labeled peptides were mixed in equal amounts in every 8-plex labeling and desalted by sep-pak vac C18cartridges (Catalog No.WAT023590, Waters, Milford, MA). The samples were then dried by the vacuum centrifuge (Catalog No. 7310038,LABCONCO, Kansas, MO) for further use.

    The peptide mixture was separated using an offline LC system (Acquity UPLC, Waters) via high-pH separation.High-pH (pH = 10) separation was performed in a reversephase column (C18, 2.1 mm × 50 mm, 1.7 μm; catalog No.186002350, Waters) using a 20-min linear gradient from 5%B buffer to 35% B buffer [A buffer: 20 mM ammonium formate in water; B buffer: 20 mM ammonium formate added in 90% acetonitrile (ACN), ammonium hydroxide was used to adjust the pH].Finally,a total of 30 fractions were collected and 2 equal-interval fractions were combined to reduce the MS running time, such as 1 and 16, and 2 and 17 [45]. In total, 15 fractions were dried and subsequently separated on a nano-LC system (Nano-Aquity UPLC, Waters) with a 75-min linear gradient from 2% D buffer to 40% D buffer (C buffer:0.1% formic acid in water; D buffer: 0.1% formic acid in ACN), which was performed on an analytical column (C18,75 μm×15 cm,3 μm;catalog No.164534,ThermoFisher Scientific, Waltham, MA). Next, the peptides were detected by MS(Q-Exactive,ThermoFisher Scientific)with 2.1 kV electrospray voltage at the mass spectrometer inlet.In addition,70 K mass resolution was applied in the full-scan of MS spectra processing (m/z 350–1200), and 17.5 K resolution was used in the following 15 sequential MS/MS scans of high energy collisional dissociation. In all studies, 1 microscan was recorded by using a dynamic exclusion of 30 s.

    The data processing strategy was also carried out according to our previous publication,with certain modification[44].The LC–MS/MS data were processed using Proteome Discoverer(version 1.4, ThermoFisher Scientific) and searched using the Sequest HT (ThermoFisher Scientific) search algorithms against the human database(20,264 entries,UniProt,released at April 10, 2014). Proteome Discoverer was searched with trypsin for protease digestion and maximally we only allowed two missed cleavages. Ten parts per million (ppm) parent ion tolerance and 0.02 Da fragment ion mass tolerance were set according to MS precision. Fixed modifications included the iTRAQ labeling of lysine and the peptide N-terminus and the carbamidomethylation of cysteine,while the variable modifications included the iTRAQ labeling of tyrosine and the oxidation of methionine. In addition, Scaffold (version 4.3.2,Proteome Software, Portland, OR) was used to validate MS/MS-based peptide and protein identifications. Peptide identifications were accepted if they could be established at greater than 90.0%probability to achieve a false discovery rate(FDR) less than 1.0% by the Scaffold local FDR algorithm.Protein identifications were accepted if they could be established at greater than 92.0%probability to achieve a FDR less than 1.0%and contained at least two identified peptides. Protein probabilities were assigned by the protein prophet algorithm. Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony. For protein quantitation, we only considered peptides that were unique to a certain given protein. The fold change between different samples was calculated by the ratio between iTRAQ reporter ion intensity and MS/MS spectra (m/z 113–119).

    Two complementary methods were combined to further characterize the secretory proteins. First, MetazSecKB, a secretome proteome knowledgebase of metazoan, was performed to screen secretory proteins, which is the most direct source to characterize the secretory proteins. Second, SecretomeP (version 2.0, DTU Health Tech, Lyngby, Denmark),a sequence-based prediction strategy for mammalian secretory proteins, was used to classify the secretory proteins:the classical secretory proteins could be correctly predicted with an N-terminal signal peptide; non-classical secretory proteins without an N-terminal signal peptide could be correctly predicted as secretory according to an NN-score> 0.6. The proteins without an N-terminal signal peptide but having a low NN-score would not be identified as secretory proteins.

    The GO annotation and signaling pathway investigation of differentially abundant secretory proteins were performed using the free online tool DAVID (http://david.abcc.ncifcrf.-gov/). The enriched signaling pathways of the differentially abundant secretory proteins were analyzed using IPA (version 7.5, Ingenuity Systems, Redwood, CA).

    Targeted proteomics and data analysis

    Before applying the targeted proteomic method of PRM analysis, the 14 most abundant proteins in serum samples were depleted by an immunoaffinity column IgY14 LC20 (Catalog No. 5188-6557, Agilent, Santa Clara, CA). The depletion was performed on an Agilent 1260 HPLC system (1260 Infinity, Agilent) following the manufacturer’s protocol. The depleted serum was concentrated and the buffer was exchanged to 100 mM TEAB using amicon 3 K concentrators(Catalog No. UFC500396, Millipore). The procedures for serum protein denaturation, reduction, alkylation, and digestion were described above.

    For PRM analysis, the unique peptides for CA2 were synthesized by fmoc solid-phase synthesis with isotope-labeled on the carboxyl side of the amino acid lysine (13C6,15N2) and purified by HPLC with purity > 99% (Anhui Guoping Pharmaceutical, Hefei, China). The heavy-labeled peptides were mixed and spiked into the tryptic digests of serum proteins at a concentration of 30 fg/μl. Detailed information for the unique CA2 peptides is displayed in Table S5.The unique peptides for CA2 should preferably have a narrow, symmetrical chromatographic peak,be 8–25 amino acids in length and ionize efficiently, provide a stable and intense signal without any modification, and not be eluted at the beginning or end time.

    The PRM analyses were also performed on a mass spectrometer(Q Exactive Plus,ThermoFisher Scientific).LC separation was executed with a trap column (C18, 75 μm × 2 cm,3 μm; catalog No. 164946, ThermoFisher Scientific) and an analytical C18column (C18, 75 μm × 15 cm, 3 μm; catalog No. 164534, ThermoFisher Scientific) on a nano-LC system(EASY-nLC 1000,ThermoFisher Scientific).2 μl of the tryptic digests and depleted serum samples were injected,and a gradient of 2% D buffer up to 35% D buffer over 30 min was applied. In all experiments, a PRM scan was performed at m/z 200,with a resolution of 70 K,an automated gain control target of 1×106,a maximum injection time of 200 ms,and an isolation window of ± 2. In addition, a normalized collision energy of 27% was used for ion dissociation, and a fixed first mass at m/z 120 was set. The inclusion list including the m/z and corresponding retention time of precursor peptides of interest is displayed in Table S5.

    The PRM raw data were analyzed with Skyline software(version 4.2, MacCoss Lab Software, Washington, WA), and a 5-min window was used.The sum of peak areas of the 5 most intense product ions was considered for protein quantification.Some ions needed to be excluded, e.g., the ions that did not match the retention time of other monitored ions, or that showed interference signals, or that gave intense signals at other retention time.

    ELISA

    The serum levels of CA2 in HCC patients were analyzed using an ELISA assay kit (Catalog No. LS-F29508, LifeSpan BioScience,Seattle,WA)following the manufacturer’s instructions.Briefly, the standard proteins and patient serum samples were diluted with the sample dilution buffer and then 50 μl of the diluted standard proteins and serum samples were added to 96-well plates. Next, 100 μl of horseradish peroxidase(HRP)-conjugate solution was added to each well and incubated for 1 h at 37 °C. Then, the plate was washed four times with PBS,and 100 μl of chromogen solution(1:1 of solution A and solution B) was added to each well and further incubated for 15 min at 37 °C. Finally, 50 μl of stop solution was added to stop the reaction, and the optical density at 450 nm was measured by a spectrophotometer(SpectraMax M5e,Molecular Devices, San Jose, CA).

    Western blot

    Western blot was performed according to a previous publication [44]. Briefly, proteins were separated on 12%SDS-PAGE and then transferred to a nitrocellulose membrane. The membrane was blocked in 5% BSA at RT for 2 h, and then incubated with the primary antibody against CA2(1:1000;catalog No.ab226987,Abcam,Cambridge,England) at 4 °C overnight. Subsequently, the mambrane was washed four times with TBST buffer,and then incubated with the corresponding secondary antibody(Catalog No.HS101-01,TransGen Biotech,Beijing,China)at RT for 1 h.Next,the blot was washed four times with TBST buffer, revealed using enhanced chemiluminescence reagents (Catalog No. 34080,ThermoFisher Scientific), and visualized by autoradiography.

    Immunohistochemistry

    IHC was carried out on HCC TMAs according to a previous publication [46]. Briefly, after pre-treating at pH = 6 and blocking with peroxidase, the sections were further incubated with the primary antibody against CA2 (1:50; catalog No.ab226987, Abcam) for 0.5 h. The sections were then treated with EnVision FLEX/HRP reagent (Catalog No. K8000,Dako, Glostrup, Denmark) for 20 min, washed three times,treated with the EnVision FLEX-DAB chromogen (Catalog No. DM827, Dako) for 3 min, and stained by Mayer’s hematoxylin (Lille’s Modification) for 3 min. They were then washed in distilled water for 5 min. All pathological sections underwent double-blind scoring as follows: negative (0), weak(1),strong (2),or very strong(3)by two different pathologists double-blindly.A score of 0 or 1 indicates low expression,and a score of 2 or 3 indicates high expression.

    Cell migration and invasion assays

    Cell migration and invasion assays were carried out according to our previously published protocols [44]. The cell migration ability was investigated by using transwell units with 8-μm pores (Catalog No. 3428, Corning Incorporated, Corning,NY), and cell invasion was investigated by using 24-well transwell inserts that had been precoated with matrigel with 8-μm pores (Catalog No. 354480, BD Biosciences, Franklin Lake,NY).A total of 1×105cells treated with exogenous recombinant CA2 (Catalog No. ab85484, Abcam) or BSA were cultured in serum-free medium in the upper chamber. To induce cell invasion and migration, DMEM supplemented with 10%FBS was placed in the lower chamber. The cells that adhered to the lower surface were then fixed with paraformaldehyde after 18 h of incubation and then stained by crystal violet(0.1%). The cells that adhered to the bottom surface were counted within 5 different views under a microscope using 20× magnification.

    Statistical analysis

    A threshold for the iTRAQ ratio was set to screen secretory proteins whose abundance was substantially changed in the C or SN tissue group compared with its corresponding DN tissue group.The proteins were defined as differentially dysregulated if the iTRAQ ratio was higher than 1.5 or lower than 0.67 in at least 5 patients, and they also must have the same direction of alteration in all 10 biological replicates. The iTRAQ ratio was based on a comparison of the reporter ion intensities between the C or SN tissue group and the corresponding DN tissue group.

    SPSS 19.0 was used for statistical analysis. Two-tailed unpaired Student’s t-test was used to compare quantitative data of CA2 from PRM and ELISA between two groups.Fisher’s exact test was applied to analyze the relationships between CA2 and clinical-pathological features. The Kaplan-Meier method was applied to calculate survival curves, while the differences were determined by using a log-rank test.Two-tailed paired Student’s t-test was used to compare quantitative data of CA2 from IHC between C and SN groups. In all analyses, P < 0.05 was considered statistically significant.

    Ethical statement

    This project was approved by the Institution Review Board of Mengchao Hepatobiliary Hospital of Fujian Medical University, Fuzhou, China. Informed consent was obtained from each participant before the operation.The use of clinical specimens was completely in compliance with the ‘‘Declaration of Helsinki”.

    Data availability

    The MS-based proteomic data have been deposited to the integrated proteome resources (iProX: IPX0001425001) [47], and are publicly accessible at https://www.iprox.cn/.

    CRediT author statement

    Xiaohua Xing: Conceptualization, Methodology, Software,Validation, Formal analysis, Investigation, Data curation,Writing - original draft, Funding acquisition. Hui Yuan:Methodology,Validation.Hongzhi Liu:Methodology,Validation. Xionghong Tan: Methodology, Validation. Bixing Zhao:Methodology, Funding acquisition. Yingchao Wang: Methodology, Resources. Jiahe Ouyang: Validation, Data curation.Minjie Lin: Resources. Xiaolong Liu: Conceptualization,Methodology, Supervision, Writing - original draft, Writing -review & editing, Project administration. Aimin Huang: Conceptualization, Supervision, Funding acquisition. All authors read and approved the final manuscript.

    Competing interests

    The authors have declared no competing interests.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China (Grant Nos. 81702910 and 81672376),the Educational Commission of Fujian Province, China(Grant No. 2018B013), and the Natural Science Foundation of Fujian Province, China (Grant Nos. 2019J01298,2017J01159, and 2016J01417). We thank Prof. Ping Xu and Dr. Yang Zhang for their help in proteomics data analysis.

    Supplementary material

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

    ORCID

    0000-0001-8971-8577 (Xiaohua Xing)

    0000-0002-7858-3539 (Hui Yuan)

    0000-0001-7102-6855 (Hongzhi Liu)

    0000-0001-9113-4159 (Xionghong Tan)

    0000-0002-0055-3989 (Bixing Zhao)

    0000-0003-3590-9233 (Yingchao Wang)

    0000-0001-9248-8805 (Jiahe Ouyang)

    0000-0002-0523-2922 (Minjie Lin)

    0000-0002-3096-4981 (Xiaolong Liu)

    0000-0002-3141-0315 (Aimin Huang)

    亚洲图色成人| 麻豆av在线久日| 成年女人毛片免费观看观看9 | 飞空精品影院首页| 婷婷色综合www| 欧美 亚洲 国产 日韩一| 久久国产亚洲av麻豆专区| 国产亚洲最大av| 国产精品 国内视频| 日本欧美视频一区| 18禁观看日本| 伊人久久国产一区二区| 久久久久久人妻| 妹子高潮喷水视频| 久久99热这里只频精品6学生| 青草久久国产| 欧美日韩一区二区视频在线观看视频在线| 宅男免费午夜| 老女人水多毛片| 久久精品久久久久久久性| 免费日韩欧美在线观看| 国产成人精品婷婷| av又黄又爽大尺度在线免费看| 亚洲国产欧美在线一区| av.在线天堂| 久热这里只有精品99| 韩国高清视频一区二区三区| 美女福利国产在线| 国产97色在线日韩免费| 18在线观看网站| 90打野战视频偷拍视频| 久久这里有精品视频免费| 亚洲国产欧美日韩在线播放| 欧美人与性动交α欧美软件| 青春草亚洲视频在线观看| 精品一区二区免费观看| 欧美精品人与动牲交sv欧美| 免费大片黄手机在线观看| 欧美97在线视频| 欧美亚洲日本最大视频资源| 亚洲精品国产av成人精品| 最近最新中文字幕免费大全7| 免费人妻精品一区二区三区视频| 欧美精品一区二区免费开放| 国产熟女午夜一区二区三区| 婷婷色综合大香蕉| 国产成人精品久久二区二区91 | 最近中文字幕高清免费大全6| 最近中文字幕高清免费大全6| 国产免费一区二区三区四区乱码| 我要看黄色一级片免费的| 久久精品国产综合久久久| 久久午夜综合久久蜜桃| 黑人欧美特级aaaaaa片| 亚洲色图综合在线观看| 国产免费一区二区三区四区乱码| 免费黄色在线免费观看| 免费看av在线观看网站| 国产免费视频播放在线视频| 丝袜人妻中文字幕| 国产在线免费精品| 亚洲人成77777在线视频| 纵有疾风起免费观看全集完整版| 亚洲精品日韩在线中文字幕| 黄片小视频在线播放| 亚洲三级黄色毛片| 成人午夜精彩视频在线观看| 欧美激情 高清一区二区三区| 熟女av电影| 最近2019中文字幕mv第一页| 美国免费a级毛片| 精品少妇久久久久久888优播| 亚洲国产精品一区三区| 可以免费在线观看a视频的电影网站 | 国产爽快片一区二区三区| 久久人人爽人人片av| 99久久精品国产国产毛片| 亚洲av免费高清在线观看| 一边摸一边做爽爽视频免费| 又黄又粗又硬又大视频| 久久人人爽av亚洲精品天堂| 91久久精品国产一区二区三区| 久久久精品94久久精品| 制服人妻中文乱码| 黑人欧美特级aaaaaa片| 大片电影免费在线观看免费| 国产白丝娇喘喷水9色精品| 午夜福利一区二区在线看| 国产毛片在线视频| 久久精品人人爽人人爽视色| 99国产综合亚洲精品| 麻豆精品久久久久久蜜桃| 在线观看美女被高潮喷水网站| 成年女人毛片免费观看观看9 | 成人亚洲精品一区在线观看| 久久亚洲国产成人精品v| av有码第一页| 电影成人av| 寂寞人妻少妇视频99o| 久久久久精品人妻al黑| 秋霞伦理黄片| videos熟女内射| 欧美成人精品欧美一级黄| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区久久| 国精品久久久久久国模美| www.熟女人妻精品国产| 亚洲成国产人片在线观看| 亚洲熟女精品中文字幕| a级毛片在线看网站| 亚洲欧美一区二区三区黑人 | 国产成人精品婷婷| 女人精品久久久久毛片| 99香蕉大伊视频| 久久99精品国语久久久| 天堂中文最新版在线下载| 2022亚洲国产成人精品| 精品一区二区免费观看| av一本久久久久| 亚洲四区av| 久久久久网色| 99久久综合免费| 男男h啪啪无遮挡| 国产日韩一区二区三区精品不卡| 国产免费一区二区三区四区乱码| 9191精品国产免费久久| 麻豆精品久久久久久蜜桃| 男人添女人高潮全过程视频| 国产高清不卡午夜福利| 亚洲精品,欧美精品| av福利片在线| 国产成人91sexporn| 黑人猛操日本美女一级片| 中文字幕av电影在线播放| 一级黄片播放器| 在线观看人妻少妇| 免费大片黄手机在线观看| 日韩av不卡免费在线播放| 日韩人妻精品一区2区三区| 亚洲av电影在线进入| 深夜精品福利| 波野结衣二区三区在线| 男女啪啪激烈高潮av片| 极品人妻少妇av视频| 欧美日韩视频高清一区二区三区二| 亚洲精品国产av成人精品| 国产成人欧美| 少妇精品久久久久久久| 香蕉精品网在线| 美女视频免费永久观看网站| 女人被躁到高潮嗷嗷叫费观| 国产成人欧美| 午夜福利一区二区在线看| 天天操日日干夜夜撸| 日本vs欧美在线观看视频| 国产日韩一区二区三区精品不卡| 中文乱码字字幕精品一区二区三区| 少妇熟女欧美另类| 亚洲国产精品国产精品| 精品一区二区三卡| 欧美在线黄色| 亚洲成人av在线免费| 最黄视频免费看| 国产一区亚洲一区在线观看| 一区二区三区四区激情视频| 久久久久久免费高清国产稀缺| 亚洲第一av免费看| 日韩在线高清观看一区二区三区| 丝袜美腿诱惑在线| 只有这里有精品99| 2018国产大陆天天弄谢| 国产成人免费无遮挡视频| 欧美精品亚洲一区二区| 老司机亚洲免费影院| 精品卡一卡二卡四卡免费| 好男人视频免费观看在线| 亚洲国产色片| 免费黄网站久久成人精品| 伦精品一区二区三区| 美女大奶头黄色视频| 天堂中文最新版在线下载| 99精国产麻豆久久婷婷| 国产熟女欧美一区二区| 永久免费av网站大全| 亚洲人成电影观看| 国产熟女午夜一区二区三区| 国产成人91sexporn| 日韩av免费高清视频| 丰满乱子伦码专区| 91在线精品国自产拍蜜月| 曰老女人黄片| 国产成人aa在线观看| 日韩人妻精品一区2区三区| 国产成人精品福利久久| av.在线天堂| 狠狠婷婷综合久久久久久88av| 人成视频在线观看免费观看| 亚洲国产av影院在线观看| 精品第一国产精品| 曰老女人黄片| 国产欧美日韩一区二区三区在线| 欧美 日韩 精品 国产| 熟妇人妻不卡中文字幕| 亚洲图色成人| 九草在线视频观看| 91aial.com中文字幕在线观看| 一级毛片电影观看| a 毛片基地| 亚洲精品在线美女| 久久久久久久国产电影| 人人妻人人爽人人添夜夜欢视频| 永久免费av网站大全| 亚洲成国产人片在线观看| av国产久精品久网站免费入址| 日韩中文字幕欧美一区二区 | 精品一品国产午夜福利视频| 丝袜喷水一区| 人人妻人人澡人人爽人人夜夜| 欧美精品一区二区免费开放| 伦理电影免费视频| 最近2019中文字幕mv第一页| 丁香六月天网| 久久狼人影院| 99re6热这里在线精品视频| 免费高清在线观看日韩| 欧美日韩视频精品一区| 人人妻人人澡人人看| 亚洲,一卡二卡三卡| 飞空精品影院首页| 捣出白浆h1v1| 亚洲美女黄色视频免费看| 搡老乐熟女国产| 少妇 在线观看| 亚洲欧美色中文字幕在线| 亚洲久久久国产精品| 国产精品国产三级专区第一集| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 亚洲精品第二区| 一个人免费看片子| 99国产综合亚洲精品| 亚洲国产av影院在线观看| 久久午夜福利片| 亚洲国产毛片av蜜桃av| 成年美女黄网站色视频大全免费| av网站在线播放免费| 春色校园在线视频观看| 最近2019中文字幕mv第一页| 777久久人妻少妇嫩草av网站| 日韩制服丝袜自拍偷拍| 久久久久国产网址| 最新中文字幕久久久久| 精品少妇内射三级| 亚洲av在线观看美女高潮| 看十八女毛片水多多多| av在线观看视频网站免费| 亚洲av男天堂| 久久午夜综合久久蜜桃| 亚洲精品国产av成人精品| 蜜桃国产av成人99| 天天影视国产精品| 亚洲精品美女久久久久99蜜臀 | 久久青草综合色| a 毛片基地| 午夜免费男女啪啪视频观看| freevideosex欧美| 国产激情久久老熟女| 国产成人a∨麻豆精品| 曰老女人黄片| 久久人人97超碰香蕉20202| 人妻 亚洲 视频| 免费高清在线观看日韩| 各种免费的搞黄视频| 亚洲综合精品二区| 曰老女人黄片| 成年av动漫网址| 国产午夜精品一二区理论片| 精品酒店卫生间| 亚洲一区二区三区欧美精品| 精品一区二区三区四区五区乱码 | 2021少妇久久久久久久久久久| 最近的中文字幕免费完整| 女人高潮潮喷娇喘18禁视频| 黄片小视频在线播放| 亚洲第一区二区三区不卡| 七月丁香在线播放| 人人妻人人添人人爽欧美一区卜| 亚洲av免费高清在线观看| 亚洲国产毛片av蜜桃av| 国产成人精品久久二区二区91 | 国产成人精品久久久久久| 各种免费的搞黄视频| 黄片播放在线免费| 亚洲国产精品成人久久小说| av在线老鸭窝| 成年美女黄网站色视频大全免费| www.熟女人妻精品国产| 天堂俺去俺来也www色官网| 欧美激情极品国产一区二区三区| 少妇人妻 视频| 咕卡用的链子| 十八禁高潮呻吟视频| 国产免费现黄频在线看| 欧美人与性动交α欧美精品济南到 | 永久网站在线| 日本午夜av视频| 亚洲国产欧美网| 男女啪啪激烈高潮av片| 午夜福利在线观看免费完整高清在| 国产亚洲最大av| 99久久人妻综合| av网站免费在线观看视频| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 国产成人精品无人区| 久久久久久久久免费视频了| 久久精品国产鲁丝片午夜精品| 成人毛片a级毛片在线播放| 日本av手机在线免费观看| 午夜福利视频在线观看免费| 午夜日本视频在线| 美女福利国产在线| 日韩精品免费视频一区二区三区| 午夜av观看不卡| 高清欧美精品videossex| 国产成人午夜福利电影在线观看| 久久精品久久久久久噜噜老黄| 99久久精品国产国产毛片| 亚洲欧美成人精品一区二区| 国产国语露脸激情在线看| 男人爽女人下面视频在线观看| 在线精品无人区一区二区三| 亚洲一级一片aⅴ在线观看| 亚洲精品在线美女| 亚洲国产精品一区二区三区在线| 久久99热这里只频精品6学生| 一区二区日韩欧美中文字幕| 制服丝袜香蕉在线| 青青草视频在线视频观看| 久久久久久久亚洲中文字幕| 日韩中字成人| 亚洲欧美精品综合一区二区三区 | 中文天堂在线官网| 久久久久视频综合| 在线观看国产h片| 亚洲av福利一区| 久久人人爽人人片av| 在线精品无人区一区二区三| 十分钟在线观看高清视频www| 人妻少妇偷人精品九色| 亚洲精品美女久久久久99蜜臀 | 天天影视国产精品| 国产极品天堂在线| 久久久久久久精品精品| 一级爰片在线观看| 欧美精品一区二区免费开放| 亚洲av免费高清在线观看| 黄色配什么色好看| 欧美日韩亚洲国产一区二区在线观看 | 中文乱码字字幕精品一区二区三区| 捣出白浆h1v1| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 久久99一区二区三区| 校园人妻丝袜中文字幕| 精品一区二区三卡| 哪个播放器可以免费观看大片| 精品卡一卡二卡四卡免费| 精品视频人人做人人爽| 亚洲天堂av无毛| 我要看黄色一级片免费的| 人人妻人人爽人人添夜夜欢视频| 亚洲伊人久久精品综合| 9热在线视频观看99| 久久久精品免费免费高清| 这个男人来自地球电影免费观看 | 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 亚洲男人天堂网一区| 免费播放大片免费观看视频在线观看| 亚洲欧美中文字幕日韩二区| 老汉色av国产亚洲站长工具| 日韩大片免费观看网站| 美女视频免费永久观看网站| 国产成人精品无人区| 国产一区二区 视频在线| 免费黄网站久久成人精品| 亚洲国产精品999| 少妇人妻 视频| 欧美日韩精品网址| 尾随美女入室| 桃花免费在线播放| 热re99久久精品国产66热6| 黄频高清免费视频| 啦啦啦中文免费视频观看日本| 日韩电影二区| 美女国产视频在线观看| 蜜桃国产av成人99| 亚洲精品日韩在线中文字幕| 欧美日韩国产mv在线观看视频| 国产伦理片在线播放av一区| 色婷婷av一区二区三区视频| 秋霞在线观看毛片| 亚洲人成77777在线视频| 黄片播放在线免费| 国产免费又黄又爽又色| 观看av在线不卡| 国产精品嫩草影院av在线观看| 亚洲精品一二三| 成年女人毛片免费观看观看9 | 国产精品女同一区二区软件| 中文字幕制服av| 免费观看在线日韩| √禁漫天堂资源中文www| 精品国产超薄肉色丝袜足j| 女人久久www免费人成看片| 可以免费在线观看a视频的电影网站 | 欧美bdsm另类| 久久午夜福利片| 伦精品一区二区三区| 国产精品秋霞免费鲁丝片| 2021少妇久久久久久久久久久| 熟女av电影| 日韩中字成人| 天天躁日日躁夜夜躁夜夜| 亚洲av男天堂| 亚洲一区中文字幕在线| 高清视频免费观看一区二区| 又大又黄又爽视频免费| 人妻人人澡人人爽人人| 精品国产一区二区三区四区第35| 如日韩欧美国产精品一区二区三区| av电影中文网址| 国产在线视频一区二区| 久久ye,这里只有精品| 自线自在国产av| 波多野结衣一区麻豆| 99热全是精品| 最近最新中文字幕大全免费视频 | 国产免费福利视频在线观看| 91精品国产国语对白视频| 人妻 亚洲 视频| 成人国语在线视频| 欧美日韩亚洲高清精品| a级毛片在线看网站| 啦啦啦视频在线资源免费观看| 两个人免费观看高清视频| 国产男人的电影天堂91| 捣出白浆h1v1| 在线精品无人区一区二区三| 桃花免费在线播放| 秋霞在线观看毛片| 久久女婷五月综合色啪小说| www.精华液| 蜜桃国产av成人99| 你懂的网址亚洲精品在线观看| 校园人妻丝袜中文字幕| 亚洲三级黄色毛片| 精品一区二区免费观看| 超碰97精品在线观看| 麻豆av在线久日| 亚洲国产色片| 国产精品av久久久久免费| 国产av国产精品国产| 免费看不卡的av| 日本vs欧美在线观看视频| 亚洲三区欧美一区| 欧美激情高清一区二区三区 | 久久久久久久久免费视频了| 午夜影院在线不卡| 日本av免费视频播放| 高清欧美精品videossex| 看免费成人av毛片| 亚洲精品国产色婷婷电影| 国产精品久久久av美女十八| 欧美bdsm另类| 青春草国产在线视频| 视频在线观看一区二区三区| 久久久久精品性色| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 国产成人精品无人区| 亚洲精品日本国产第一区| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久av不卡| 欧美黄色片欧美黄色片| 又粗又硬又长又爽又黄的视频| 婷婷色av中文字幕| 亚洲精品aⅴ在线观看| 日日摸夜夜添夜夜爱| 国产av一区二区精品久久| 亚洲国产av新网站| 激情视频va一区二区三区| 在线看a的网站| 人妻人人澡人人爽人人| 九草在线视频观看| 青青草视频在线视频观看| 热99久久久久精品小说推荐| 18禁裸乳无遮挡动漫免费视频| 青春草国产在线视频| 天堂中文最新版在线下载| 啦啦啦视频在线资源免费观看| 国产精品国产av在线观看| 日日撸夜夜添| 亚洲成人一二三区av| 美女大奶头黄色视频| av电影中文网址| 国产一区二区激情短视频 | 少妇人妻久久综合中文| 久久久久久免费高清国产稀缺| 老鸭窝网址在线观看| 在线 av 中文字幕| 黑丝袜美女国产一区| av卡一久久| 亚洲中文av在线| 久久精品久久精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲国产av影院在线观看| 国产精品一二三区在线看| 男女下面插进去视频免费观看| av卡一久久| 国产精品二区激情视频| 久久久国产一区二区| 99久久人妻综合| 欧美日韩视频精品一区| 亚洲欧美精品综合一区二区三区 | 91aial.com中文字幕在线观看| 中文乱码字字幕精品一区二区三区| 波野结衣二区三区在线| xxxhd国产人妻xxx| 国产午夜精品一二区理论片| 精品午夜福利在线看| 午夜福利视频在线观看免费| 桃花免费在线播放| 久久精品国产a三级三级三级| 欧美日韩精品成人综合77777| 久久久精品区二区三区| 免费高清在线观看视频在线观看| 免费黄色在线免费观看| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 亚洲男人天堂网一区| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 日韩不卡一区二区三区视频在线| 久久久久久久大尺度免费视频| 青草久久国产| 久久久久久免费高清国产稀缺| 久久99蜜桃精品久久| 性色av一级| 亚洲色图 男人天堂 中文字幕| 少妇人妻精品综合一区二区| 中文天堂在线官网| 五月天丁香电影| 美女脱内裤让男人舔精品视频| 最近手机中文字幕大全| 肉色欧美久久久久久久蜜桃| 高清欧美精品videossex| 中文字幕人妻丝袜制服| 国产精品熟女久久久久浪| 亚洲欧美精品自产自拍| 在线 av 中文字幕| 交换朋友夫妻互换小说| 卡戴珊不雅视频在线播放| 久久久久国产网址| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 精品一区在线观看国产| 一级爰片在线观看| 精品午夜福利在线看| 美女中出高潮动态图| 免费av中文字幕在线| 另类亚洲欧美激情| 久久人人爽av亚洲精品天堂| 精品少妇黑人巨大在线播放| 一区二区三区四区激情视频| 久久久久久久久久久免费av| 午夜91福利影院| 中文字幕人妻丝袜制服| 老司机亚洲免费影院| 免费观看在线日韩| 我要看黄色一级片免费的| 你懂的网址亚洲精品在线观看| 亚洲av免费高清在线观看| 交换朋友夫妻互换小说| 亚洲精品第二区| 波多野结衣一区麻豆| 丝袜在线中文字幕| 69精品国产乱码久久久| 亚洲av综合色区一区| 国产精品一国产av| 久久久久久人人人人人| 搡女人真爽免费视频火全软件| 精品亚洲成a人片在线观看| 大码成人一级视频| 国产av一区二区精品久久| 国产成人av激情在线播放| 亚洲中文av在线| 午夜老司机福利剧场| 满18在线观看网站| 最近2019中文字幕mv第一页| freevideosex欧美| 99热全是精品| 亚洲人成电影观看| 国产乱来视频区| 国产一区二区三区av在线| 精品99又大又爽又粗少妇毛片| 亚洲综合色惰| 男女无遮挡免费网站观看| 国产片内射在线| 久久久久精品人妻al黑| 亚洲成人手机| 亚洲av日韩在线播放| 狂野欧美激情性bbbbbb| 大香蕉久久成人网|