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

    Fas -670 A/G polymorphism predicts prognosis of hepatocellular carcinoma after curative resection in Chinese Han population

    2022-03-02 06:52:04XinYiGuoYiLinXiZhiYngLiRnLiYngFnWenMoYoFeiLiu

    Xin-Yi Guo , Yi-Lin Xi , Zhi-Yng Li , Rn Li , Yng Fn , Wen-Mo Yo , Fei Liu ,*

    a West China School of Medicine, West China Hospital, Sichuan University, Chengdu 610041, China

    b Department of Liver Surgery & Liver Transplantation Center, West China Hospital, Sichuan University, Chengdu 610041, China

    Keywords:Hepatocellular carcinoma Single nucleotide polymorphism Cell apoptosis Prognosis Resection

    ABSTRACT

    Background: Apoptosis, also called programmed cell death, is a genetically controlled process against hyperproliferation and malignancy. The Fas-Fas ligand (FasL) system is considered a major pathway for apoptosis in cells and tissues. Thus, this study aimed to investigate whether single nucleotide polymorphisms (SNPs) in Fas and FasL gene may have effects on the recurrence and survival of patients with hepatocellular carcinoma (HCC) after curative hepatectomy.

    Methods: We investigated the relationship between Fas rs1800682, rs2234767 and FasL rs763110 polymorphisms and recurrence-free survival (RFS) as well as overall survival (OS) in 117 Chinese Han patients with HCC who underwent hepatectomy.

    Results: In Kaplan-Meier survival analysis, only Fas rs1800682 (-670 A/G) was associated with RFS and OS. Compared with AA genotype, the AG/GG genotype was significantly associated with better RFS ( P = 0.008) and OS ( P = 0.020). Moreover, multivariate Cox regression analysis showed that Fas rs1800682 remained as a significant independent predictor of RFS for HCC patients with hepatectomy[AG/GG vs. AA: adjusted hazard ratio = 0.464, 95% confidence interval: 0.275-0.782, P = 0.004], but was not an independent predictor of OS ( P = 0.395).

    Conclusions: This study demonstrated that Fas -670 G allele may play a protective role in the recurrence and survival of HCC patients with hepatectomy. Furthermore, Fas rs1800682 polymorphism might be a promising biomarker for HCC patients after hepatectomy.

    Introduction

    Liver cancer is one of the most common malignant tumors worldwide. In China, liver cancer ranks the third in both most commonly diagnosed cancer and cancer mortality [1] . Of all the liver cancer cases, hepatocellular carcinoma (HCC) approximately comprises 75% [2] . Among the various therapeutic options of HCC,liver resection is one of the first-line treatments for HCC at early stage [3] . However, the prognosis after hepatectomy still remains unsatisfactory. The rate of cancer recurrence within five years after surgery keeps high to 60%-70% [4] . Postsurgical recurrence of HCC includes common intrahepatic recurrence (accounting for most of recurrences) and extrahepatic recurrence [5] . In recent years, various potential prognostic predictors of recurrences have been reported, including microvascular invasion (MVI) [ 2 , 4 , 6-8 ], tumor size [ 4 , 6-9 ], the number of tumors [ 2 , 4 , 7 , 8 ], and alpha-fetoprotein(AFP) [ 7 , 8 ]. However, none of them can meet the requirement of clinical accuracy. For example, several studies found that patients with tumors larger than 5 cm have a significantly higher risk of recurrence [ 4 , 6-9 ], while a few studies failed to find a significant effect between tumor size and recurrence rate, which means that patients with smaller tumor may also have high risk of recurrence [ 10, 11]. It was suggested that individual genetic factors may impact HCC recurrence. According to the development of individualized comprehensive therapy, accumulated evidences have indicated that single nucleotide polymorphisms (SNPs) were associated with prognoses of HCC [12] , non-small cell lung cancer(NSCLC) [ 13 , 14 ], and squamous cell carcinoma of the oropharynx(SCCOP) [15] .

    Apoptosis, also called programmed cell death, is a genetically controlled process against hyperproliferation and malignancy [16] .The Fas-Fas ligand (FasL) system is a major pathway for apoptosis in cells and tissues [17] . The Fas is mapped on the human chromosome 10q24.1 and widely expressed in normal and neoplastic cells, including normal hepatocytes and HCC cells [ 17 , 18 ].The FasL is mapped on the human chromosome 1q23 and expressed in activated T cells or natural killer (NK) cells [ 18 , 19 ].Combination of Fas and FasL transmits the signal for cell death to target cell and then triggers extrinsic and intrinsic pathways in human [ 16 , 20 ]. There is accumulating evidence that expression of Fas, FasL and related molecules plays an essential role in pathogenesis of many kinds of tumors including HCC [ 17 , 21 ]. Studies showed that down-regulation of theFasgene may help tumor cells escape from immune surveillance, and up-regulation of theFasLgene may increase the ability of tumor cells through inducing apoptosis of lymphocytes [ 22 , 23 ].

    Several potential functional SNPs ofFasandFasLgene were found to be associated with susceptibility, progression and prognosis of cancers. The Fas -670 A>G (rs1800682) have been proved to be related to higher risk of many kinds of cancers, including breast cancer, esophageal cancer, prostate cancer, HCC [24] , and prognoses of NSCLC [ 13 , 14 ], SCCOP [15] , and chronic myeloid leukemia [25] .One of these studies showed that the GG genotype and combined AG/GG genotypes are associated with a poor prognosis in NSCLC [14] . The FasL-844 C>T (rs763110) also was found to be associated with prognoses of chronic myeloid leukemia [25] , breast cancer [26] , and bladder cancer [27] . The CC genotype and C allele were suggested to be both possibly related to a better prognosis of breast cancer [26] . In addition, the Fas -1377 G>A (rs2234767)is related to prognosis of breast cancer [28] . Sunter et al. reported that patients with A allele have worse prognosis compared with GG genotype in acute promyelocytic leukemia [29] .

    The present study aimed to analyze the relationship between three key SNPs and the survival of patients with HCC and HCC recurrence after resection.

    Methods

    Study populations

    A total of 117 Chinese Han patients newly diagnosed with HCC from January 2011 to November 2014 were enrolled in this study.The diagnoses were established by postoperative histopathologic examination. All the patients who underwent liver resection in the Department of Liver Surgery, West China Hospital of Sichuan University (Chengdu, China), without any treatments before the resection, were included. The exclusion criteria were: 1) pregnancy; 2)having mental disorder; and 3) having cancer from other organs.

    Written informed consents were obtained from all patients.Clinical information included preoperative results of clinical laboratory tests, such as AFP level, tumor size, the number of tumor(s),existence of portal vein tumor thrombus and the other vascular invasion, distant metastasis, lymph node metastasis, liver cirrhosis,differentiation degree of cancer cell, and preoperative Child-Pugh grade. Identification of the clinical stage of tumor is based on the unified 8th edition TNM (tumor-node-metastasis) criteria [30] .

    We followed up the 117 patients to November 2019. The recurrence-free survival (RFS) time was calculated from the date of performing tumor resection to the date of diagnosing recurrence or the last follow-up. Tumor recurrence was considered as appearance of new intra- or extra-hepatic tumor nodule(s), with or without an increase in serum AFP level [5] . The overall survival (OS)time was defined from the date of the operation to the date of allcause death or the last follow-up. This study was coincident withtheDeclarationofHelsinki, and approved by the Ethics Committee of Sichuan University.

    Table 1 SNPs of cell death pathway genes evaluated in this study.

    Selection of candidate SNP

    The SNPs were selected based on previous articles reporting the association with development and prognosis of all kinds of cancers, especially HCC. At the same time, we focused on a series of comprehensive SNP-selecting principles to ensure that our research went on wheels as far as possible. The candidate SNPs were the SNPs genotyped in the Chinese Han populations (Han Chinese in Beijing and Southern Han Chinese), chosen from the 10 0 0 Genomes Project database using the Phase 3 data. Selection of tagging SNPs should satisfy the following requirements: the minor allele frequency (MAF) of selected SNPs ≥0.1 in Chinese Han people and only choosing one when several SNPs have a linkage coefficientr2 >0.8. Eventually, three tagging SNPs [31-33] in two candidate genes were chosen for this study ( Table 1 ), including two SNPs (rs1800682, rs2234767) from Fas and only one SNP (rs763110)from FasL.

    Genotype analysis

    Genomic DNA was extracted from whole blood samples of patients with the whole blood DNA kit (Biotake corporation, Beijing,China). According to the manufacturer’s protocol, DNA was isolated from 200 μL of the whole blood. And then, the concentration of isolated DNA was diluted to 20 ng/ μL for working solution, which was stored at -20 °C until analysis. SNP genotyping was conducted on MassARRAY system (Sequenom, San Diego, CA, USA) with matrix assisted laser desorption ionisation-time of flight mass spectrometry method (MALDI-TOF), according to the manufacturer’s protocol. Primers for polymerase chain reaction (PCR) and single base extension were both designed by Assay Design software package (Sequenom). More detailed process could be found in our previous study [34] and the supplementary materials (Table S1). For quality control, 5% of 117 samples were randomly chosen to analyze twice for each locus. Moreover, in order to further validate the genotyping result of MALDI-TOF, about 5% of samples were further verified by direct sequencing method.

    Statistical analysis

    The characteristics of patient data were expressed as numbers (percentages). The relations of the three tagging SNPs to RFS and OS were identified utilizing the Kaplan-Meier survival analysis with the log-rank test. The relations of the other variables to RFS and OS were estimated using the Cox univariate regression analysis. Those significant variables in univariate analysis were included in a Cox multivariate model, using the forward stepwise (likelihood ratio), to complete multivariate analysis. All statistical analyses were conducted utilizing SPSS version 23.0 statistical software(SPSS Inc., Chicago, IL, USA). AP<0.05 was considered statistically significant.

    Table 2 Baseline characteristics and clinical characteristics of the patients.

    Results

    Characteristics of populations and genotype frequency

    The baseline characteristics of the 117 HCC patients are summarized in Table 2 . Their age ranged from 24 to 71 years with a median age of 50 and 82.9% were male. Among the 117 patients, 68 (58.1%) had smoking history (all of them smoke more than 400 cigarettes one year), and 13 (11.1%) had family history of HCC, while nearly all the patients (99.1%) had an HBV infection history. Besides, only two patients had regional lymph node metastasis; one had distant metastasis at first diagnosis. One hundred(92.6%) patients had Child-Pugh grade A.

    Among the investigated SNPs, as shown in Table 3 , the AA, AG and GG genotypes in Fas -670 A/G (rs1800682) corresponded to 26.5%, 56.4% and 17.1% of patients, respectively; the frequencies of major allele (A) and minor allele (G) were 0.55 and 0.45, respectively. For Fas -1377 G/A (rs2234767), patients with genotypes AA,AG and GG account for 14.5%, 53.8% and 31.6%; the frequencies of major (G) and minor (A) allele were 0.59 and 0.41, respectively. The proportion of patients with CC, CT and TT genotypes in FasL -844 T/C (rs763110) were 50.4%, 40.2% and 9.4%; and the frequencies of major (C) and minor (T) allele were 0.71 and 0.29.

    Univariate analysis of genotype and clinicopathological variables

    The 1-, 3-, and 5-year RFS rates were 56.1%, 40.9%, 33.9%, respectively. The association between three genotypes (AA, AG, GG)and RFS were analyzed firstly. As shown in Fig. 1 A, the AG and GG genotypes in rs1800682 were both associated with better 5-year RFS after liver resection, compared with the AA genotype(P= 0.027). For the further exploration of the advantages of G allele, the analysis between AA homozygotes and the other two genotypes, as well as between GG homozygotes and the other two genotypes in rs1800682 were conducted. The RFS rate in the AG/GG genotypes was significantly better versus AA homozygotes(P= 0.008, Fig. 1 B). Conversely, there was no significant difference of RFS between the GG homozygotes and AA/AG genotypes(P= 0.837, Fig. 1 C). And Fig. 2 A, B show the RFS of the other two SNPs by every genotype.

    The 1-, 3-, and 5-year OS rates were 86.6%, 74.6%, and 70.1%, respectively. The association between three genotypes (AA, AG, GG)and OS were also analyzed. As shown in Fig. 3 A, no significant association among the OS of AA, AG and GG genotypes in rs1800682 were observed (P= 0.062), though an important phenomenon that AG and GG genotypes tend to manifest better 5-year OS compared with AA genotype was noticed. As the exploration in RFS,the analysis between AA homozygotes and the other two genotypes, as well as between GG homozygotes and the other two genotypes in rs1800682 were conducted, respectively. The OS rate in the AG/GG genotypes was significantly better than AA homozygotes (P= 0.020, Fig. 3 B). Conversely, there was no significant difference of OS between the GG homozygotes and AA/AG genotypes(P= 0.365, Fig. 3 C). And Fig. 2 C, D show the OS of the other two SNPs by every genotype.

    In addition, we detected the association between other clinicopathological factors and HCC prognosis after resection. By the Cox univariate regression analysis, a younger age, multiple tumors, tumor size larger than 5 cm, advanced tumor staging (TNM III/IV stages), and the presence of background cirrhosis were significantly associated with shorter RFS ( Table 4 ), and multiple tumors,advanced tumor staging (TNM III/IV stages), and the presence of background cirrhosis were significantly associated with worse OS( Table 5 ).

    Cox regression model in multivariate analysis

    As shown in Tables 4 and 5 , Fas -670 A/G (rs1800682) remained as a significant independent predictor of RFS for HCC patients with hepatectomy (AG/GG vs. AA: adjusted hazard ratio = 0.464, 95%confidence interval: 0.275-0.782,P= 0.004), but was not an independent predictor of OS (P= 0.395).

    Discussion

    To the best of our knowledge, this is the first study to investigate whether Fas -670 A/G polymorphism may have effects on the recurrence and survival of HCC patients undergoing curative resection. In this study, we investigated the relationship between three SNPs fromFasandFasLgene and clinical outcomes of HCC patients undergoing hepatic resection. The results suggested that AG and GG genotypes in rs1800682 were significantly associated with a better 5-year RFS and OS compared with AA genotype. Furthermore, the SNP, rs1800682 ofFasgene, could serve as an independent predictive factor of postoperative recurrence for resected HCC.

    There has been no recognized comprehensive prognostic evaluation system up to now, but in the past decades, numerous potentially prognostic predictors have been reported, not only including several common variables, like sex [ 4 , 8 ], age [ 2 , 35 ], tumor size [ 4 , 6-9 ], the number of tumors [ 2 , 4 , 7 , 8 ], presence of liver cirrhosis [ 4 , 36 ], presence of MVI [ 2 , 4 , 6-8 ], TNM stage [7] , BCLC stage [7] , Edmond-Steiner grade [2] , and Child-Pugh class [7] , but also including some newly discovered indices, such as neutrophilto-lymphocyte ratio [6] , Glasgow prognostic score [6] , serum level of protein induced by vitamin K absence or antagonist (PIVKA) [6]and AFP [ 7 , 8 ], albumin-bilirubin (ALBI) grade [8] , tumor expression of epithelial cell adhesion molecule (EpCAM) [9] , plateletcount [2] , and intraoperative estimated blood loss (EBL) [2] . However, as far as we know, no one has been universally recognized.Despite the existence of various staging systems for HCC, none of them could maintain enough accuracy for prognostic prediction. Although TNM stage, developed by the American Joint Committee on Cancer (AJCC) and Union for International Cancer Control (UICC), is used worldwide [37] and considered as a comparatively effective prognostic indicator in many kinds of malignant tumors, such as lung cancer [13] , it is not very satisfying in HCC.There was a study reporting that TNM stage was not an independent predictor of HCC recurrence [9] . Besides, BCLC stage was not always effective. Another study published in 2014 also indicated that both TNM and BCLC stage were not always associated with recurrence and survival when choosing different follow-up time [38] .

    Table 3 Allele frequency of Fas and FasL polymorphisms.

    Fig. 1. Kaplan-Meier curves of recurrence-free survival of patients with tumor resection by Fas -670 A/G (rs1800682).

    According to the development of individualized comprehensive therapy, genotype analysis plays an increasingly important role in prediction of clinical outcomes and guidance of clinical treatment. Therefore, discovering new genetic biomarkers could facilitate personalized treatment or be incorporated into certain staging systems to improve predictive accuracy of resected HCC, thereby improving the prognoses of patients. Until now, SNPs of gene in cell death pathway have achieved a lot of progress. In two metaanalysis published in 2014, Fas -670 A/G polymorphism was reported to be associated with susceptibility of HCC and many kinds of cancers [ 24 , 39 ]. In addition, a recent study from our team has shown that CASP8 -652 6N ins/del polymorphism in cell death pathway may play an important role in the development, progression, and survival of HBV-related HCC [16] . Nonetheless, related research on the relationship between Fas -670 A/G polymorphism and prognosis of resected HCC is still absent.

    Fig. 2. Recurrence-free survival and overall survival by genotype of rs2234767 and rs763110 in patients ( A-D ).

    Table 4 Univariate and multivariate analyses of genotypes and clinicopathological parameters associated with recurrence-free survival.

    Fig. 3. Kaplan-Meier curves of overall survival of patients with tumor resection by Fas -670 A/G (rs1800682).

    Table 5 Univariate and multivariate analyses of genotypes and clinicopathological parameters associated with overall survival .

    Fas, encoded byFasgene, is widely expressed on normal and neoplastic cells (including normal hepatocytes, HCC cells, immune cells and so on), and initiates death receptor pathway by the combination with FasL [17] . Our results illustrate that Fas rs1800682 A>G polymorphism was associated with better survival in patients with resected HCC. This SNP was an independent predictor of 5-year RFS for resected HCC patients. According to previous research data, rs1800682 was located in the enhancer region ofFasgene, which was combined with STAT1 to modulate the transcription activity ofFasgene; theA-to-Gsubstitution polymorphism can diminish promoter activity and downregulate gene expression [24] .The explanation of our result has three aspects. Firstly, activated T-lymphocytes, having a crucial role in immune surveillance, can be induced into apoptosis by activation induced cell death (AICD),mediated mainly by the death-receptor pathway, which could be beneficial to tumor cells escaping from immune surveillance [40] .Hence, in stimulated T-lymphocyte, this A>G polymorphism downregulatedFasgene expression [22] and thus reduced AICD;sequentially, patients with G allele are expected to have lower risk of HCC recurrence compared with patients without G allele, and it is biologically plausible that the Fas -670 A>G polymorphism may be associated with better prognosis of HCC. Secondly, if we hypothesize that no extra mutation appears at Fas 1800682 in tumor cells, and therebyFasgene expression was also downregulated in tumor tissues, which was similar to that in T lymphocytes, we could deduce that apoptosis of tumor cells might be reduced as well, and thus risk of recurrence and death might increase. However, the scenario in T lymphocytes is not applicable to the malignancy because malignant cells, including HCC cells, have acquired the capacity to resist apoptotic stimuli due to the aberrations of the death receptors pathway [16] . Furthermore, if there are exactly some mutations presenting at Fas 1800682 in tumor cells, it still could not function properly because of aforementioned reason.Therefore, Fas -670 A/G polymorphism might not significantly influence apoptosis of tumor cells, theoretically. From the above discussion, patients with G allele are highly possible to have longer RFS due to the reduced AICD of stimulated T lymphocytes. Thirdly,Fas protein is a double-edged sword, because not only can it mediate apoptosis by death receptor pathway, but also it has another form, circulating soluble Fas (sFas) that can inhibit Fas receptorligand binding in extracellular space to protect tumor cells from Fas-mediated apoptosis [41] . Therefore, we could deduce that patients with AG and GG genotypes might have lower level of sFas in serum compared with whom with AA genotype, which could protect tumor cells from Fas-mediated immune elimination, and thus patients with AG and GG genotypes have lower risk of recurrence.A study [41] in Chinese people showed that there was significantly higher level in sFas expression for patients with Fas -670 AA genotype compared with whom with the other two genotypes. Xu’s study was in consistent with our hypothesis. In addition, OS was influenced by many factors, such as post-operative surveillance for recurrence (regular vs. irregular) [4] , treatment modality for late recurrence (potentially curative vs. non-curative) [4] , the type of recurrence [36] (in this study, recurrence was divided into three types, marginal, nodular, and multinodular or diffuse; more specific information can be found in this article) and so on. With regular postoperative recurrence surveillance, recurrence could be discovered as soon as possible, which is beneficial for instant treatment [4] , thereby improving post-recurrence survival. This complexity might be the reason why Fas -670 A/G was not the independent predictor of OS in this study.

    In recent years, numerous studies including meta-analysis have discussed the association between individual Fas -670 A/G polymorphisms and many cancers; however, the results are controversial. Some studies showed that G allele was protective. A study published in 2013 [13] suggested that a favorable survival of NSCLC was significantly associated with genotypes of Fas rs1800682 GG,which was similar to our results. Another research on prostate cancer [42] reported that Fas -670 AG or GG genotypes were related to decreased prostate cancer risk. Furthermore, the present study reveals that Fas -670 AG and GG genotypes were associated with favorable prognosis. Some other studies showed that G allele was risky [ 14 , 15 ]. A study on NSCLC [14] indicated that Fas-670 GG genotype and the combined GG and AG genotypes exhibited a significantly decreased survival than AA genotype. Another study on SCCOP [15] reported that patients with Fas -670 AG/GG genotypes had higher risk of recurrence than those with AA genotype. Moreover, two meta-analyses [ 24 , 39 ] reported that Fas -670 G allele conferred an increased risk of HCC. However, some studies also found that no significant relationship was observed between Fas -670 A/G polymorphism and OS of malignant pleural mesothelioma (MPM) [43] , and nasopharyngeal carcinoma [44] . Some studies showed that there was no correlation between Fas rs1800682 polymorphism and risk of cervical cancer [45] , and breast cancer [46] . Cancer is a multifactorial disease resulting from complicated interactions between numerous genetic and environmental factors [24] . Explanation of the discrepancy might be attributable to the diversities of cancer types, ethnic backgrounds, treatment modalities, etc. for tumor [16] and living environment. For example, the frequency of Fas -670 A allele was significantly different among African, American, European and Han Chinese population (24.7% for African, 46.1% for American, 56.9% for European,58.6% for Southern Han Chinese, and 61.7% for Han Chinese in Beijing) [31] .

    One unexpected discovery was that age might be an independent predictor of recurrence for resected HCC. We speculate that it might not have practical and clinical value, though it is statistically significant, because adjusted HR (0.968) and 95% CI (0.947-0.990)were quite nearly 1. However, based on our relatively small sample size, the practical value might be underestimated. So, the predictive value of age might be explored by studies with larger sample size.

    The limitations of this study are the following. Firstly, potential selection bias might exist due to hardly fully controlled lost rate of follow-up. Secondly, the sample size of this study was relatively small, and thus more analyses utilizing a larger sample size would be required. Thirdly, our data were from Chinese population, so it is not sure whether these findings are applicable for other ethnic groups. Therefore, other validation enrolling other ethnics would be required. And fourth, the diversified tumor characteristics may affect the prognoses of patients undergoing HCC resection, though some of these factors were adjusted in multivariate analysis.

    In conclusion, the Fas -670 A/G polymorphism was associated with recurrence and survival in HCC patients who underwent hepatectomy; moreover, this polymorphism could serve as an independent predictive factor of postoperative recurrence for resected HCC among Chinese Han population. Our study showed that Fas -670Gallele may play a protective role in HCC patients with hepatectomy.

    Acknowledgments

    We gratefully thank the statistical guidance of Dr. Chao-Yue Chen from Department of Neurosurgery & State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, China.

    CRediT authorship contribution statement

    Xin-Yi Guo : Conceptualization, Methodology, Funding acquisition, Data curation, Investigation, Formal analysis, Project administration, Writing - original draft, Writing - review & editing. Yi-Lin Xia : Investigation, Formal analysis, Writing - original draft, Writing - review & editing. Zhi-Yang Li : Investigation. Ran Li : Investigation. Yang Fan : Investigation. Wen-Mo Yao : Investigation. Fei Liu : Conceptualization, Methodology, Investigation, Resources, Supervision, Validation, Writing - review & editing.

    Funding

    This study was supported by grants from the Sichuan Science and Technology Program (2019YFS0370) and National University Student Innovation and Entrepreneurship Training Programs(C2019104465).

    Ethical approval

    The study involving human participants was reviewed and approved by the Ethics Committee of Sichuan University. Written informed consent was obtained from all participants.

    Competing interest

    No benefits in any form have been received or will be received from a commercial party related directly or indirectly to the subject of this article.

    Supplementary materials

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.hbpd.2021.07.001 .

    亚洲aⅴ乱码一区二区在线播放| 午夜亚洲福利在线播放| 啦啦啦免费观看视频1| 婷婷精品国产亚洲av在线| 一个人看视频在线观看www免费 | 别揉我奶头~嗯~啊~动态视频| 久久天躁狠狠躁夜夜2o2o| 国产欧美日韩精品亚洲av| 日韩 欧美 亚洲 中文字幕| 啦啦啦观看免费观看视频高清| 别揉我奶头~嗯~啊~动态视频| 亚洲精品久久国产高清桃花| 黄频高清免费视频| 九九热线精品视视频播放| 欧美黑人欧美精品刺激| 成人三级黄色视频| 此物有八面人人有两片| 人妻夜夜爽99麻豆av| 精品国产亚洲在线| 白带黄色成豆腐渣| 成人高潮视频无遮挡免费网站| 精品国内亚洲2022精品成人| 亚洲第一电影网av| 九色国产91popny在线| 日韩欧美国产在线观看| aaaaa片日本免费| 国产亚洲av嫩草精品影院| 精品久久久久久久末码| 国产黄色小视频在线观看| 我要搜黄色片| 天天添夜夜摸| 国产av不卡久久| 亚洲av免费在线观看| 精品国产美女av久久久久小说| 一区福利在线观看| 亚洲 国产 在线| 老汉色∧v一级毛片| 中文亚洲av片在线观看爽| 成年版毛片免费区| 不卡av一区二区三区| 国产亚洲精品一区二区www| 亚洲一区二区三区不卡视频| 亚洲中文日韩欧美视频| 国产探花在线观看一区二区| 欧美一级a爱片免费观看看| 91老司机精品| 99热这里只有是精品50| 国产 一区 欧美 日韩| 欧美极品一区二区三区四区| 1024手机看黄色片| 亚洲欧美日韩高清专用| 99在线视频只有这里精品首页| 后天国语完整版免费观看| 国产精品自产拍在线观看55亚洲| 欧美一区二区国产精品久久精品| 亚洲国产欧洲综合997久久,| 国产精品av视频在线免费观看| 中文字幕人成人乱码亚洲影| 日韩欧美精品v在线| 美女大奶头视频| 无遮挡黄片免费观看| 亚洲成人中文字幕在线播放| 美女免费视频网站| www日本黄色视频网| 欧美av亚洲av综合av国产av| 美女免费视频网站| 国产蜜桃级精品一区二区三区| 午夜免费观看网址| 天堂动漫精品| 丰满的人妻完整版| 国产不卡一卡二| 久久99热这里只有精品18| 国产伦精品一区二区三区视频9 | 亚洲无线观看免费| 91九色精品人成在线观看| 一本综合久久免费| 日韩国内少妇激情av| 国产伦在线观看视频一区| 免费看光身美女| 在线十欧美十亚洲十日本专区| 不卡av一区二区三区| 757午夜福利合集在线观看| 国产精品国产高清国产av| 超碰成人久久| 在线十欧美十亚洲十日本专区| 亚洲最大成人中文| 老熟妇乱子伦视频在线观看| 老熟妇乱子伦视频在线观看| 婷婷精品国产亚洲av| 亚洲中文字幕日韩| 不卡一级毛片| 91字幕亚洲| 亚洲精品一区av在线观看| 97人妻精品一区二区三区麻豆| 少妇人妻一区二区三区视频| 91av网一区二区| 最新中文字幕久久久久 | 国产精品av久久久久免费| 色综合亚洲欧美另类图片| 97超视频在线观看视频| 男女床上黄色一级片免费看| 18禁国产床啪视频网站| 一二三四社区在线视频社区8| 欧美日韩中文字幕国产精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 99国产综合亚洲精品| 啪啪无遮挡十八禁网站| 丝袜人妻中文字幕| 别揉我奶头~嗯~啊~动态视频| 老司机午夜十八禁免费视频| 很黄的视频免费| 女人被狂操c到高潮| 老司机福利观看| 色尼玛亚洲综合影院| 亚洲欧美日韩高清在线视频| 亚洲av成人一区二区三| 免费看光身美女| av福利片在线观看| 日韩中文字幕欧美一区二区| 久久久久久九九精品二区国产| 波多野结衣高清无吗| 久久婷婷人人爽人人干人人爱| 在线免费观看的www视频| 亚洲精品美女久久av网站| 嫩草影院精品99| 国产成人精品久久二区二区91| 中国美女看黄片| 久久久国产成人免费| e午夜精品久久久久久久| 美女扒开内裤让男人捅视频| 中文字幕最新亚洲高清| 91老司机精品| 国产午夜精品论理片| 成人高潮视频无遮挡免费网站| 国产亚洲av高清不卡| 国产爱豆传媒在线观看| 国产三级在线视频| 久久亚洲精品不卡| 亚洲av免费在线观看| 90打野战视频偷拍视频| 国产成人精品久久二区二区91| 久9热在线精品视频| 99在线视频只有这里精品首页| xxxwww97欧美| 成年女人毛片免费观看观看9| 欧美日韩亚洲国产一区二区在线观看| 久久久久久久久免费视频了| 亚洲无线观看免费| 小说图片视频综合网站| 悠悠久久av| 久久人人精品亚洲av| 日本 av在线| 长腿黑丝高跟| 叶爱在线成人免费视频播放| 欧美zozozo另类| 亚洲色图 男人天堂 中文字幕| 成年免费大片在线观看| 国产久久久一区二区三区| 天堂网av新在线| 亚洲专区字幕在线| 亚洲成人免费电影在线观看| 日韩精品青青久久久久久| 少妇丰满av| 国产成人福利小说| 岛国在线免费视频观看| 一本久久中文字幕| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久久久久久久| 国产探花在线观看一区二区| 久久精品亚洲精品国产色婷小说| 网址你懂的国产日韩在线| 精品国产亚洲在线| 亚洲自拍偷在线| 欧美午夜高清在线| 国产精品自产拍在线观看55亚洲| 国模一区二区三区四区视频 | 久久国产精品影院| 午夜视频精品福利| 久久精品夜夜夜夜夜久久蜜豆| 美女黄网站色视频| 国产一区二区三区在线臀色熟女| 91在线精品国自产拍蜜月 | 亚洲欧美日韩卡通动漫| a在线观看视频网站| 久久这里只有精品19| 久久久久精品国产欧美久久久| 99国产精品99久久久久| 激情在线观看视频在线高清| 两个人的视频大全免费| 久久国产乱子伦精品免费另类| 亚洲,欧美精品.| 日本a在线网址| 午夜福利视频1000在线观看| 国产美女午夜福利| 美女扒开内裤让男人捅视频| 麻豆一二三区av精品| 亚洲精品一卡2卡三卡4卡5卡| 九九在线视频观看精品| 久久国产乱子伦精品免费另类| 国产一区二区激情短视频| 日韩av在线大香蕉| 国产美女午夜福利| 制服人妻中文乱码| 一本综合久久免费| 老汉色∧v一级毛片| 啦啦啦免费观看视频1| 欧美性猛交黑人性爽| 亚洲国产欧美一区二区综合| aaaaa片日本免费| 久久久精品大字幕| 国产亚洲精品一区二区www| 日韩中文字幕欧美一区二区| 欧美在线黄色| 久久久水蜜桃国产精品网| 亚洲七黄色美女视频| 久99久视频精品免费| svipshipincom国产片| 亚洲一区二区三区不卡视频| 色哟哟哟哟哟哟| 亚洲自拍偷在线| 亚洲熟妇中文字幕五十中出| av在线天堂中文字幕| 欧美一区二区国产精品久久精品| 性色av乱码一区二区三区2| 久久久国产精品麻豆| 国产精华一区二区三区| 欧美丝袜亚洲另类 | 精品日产1卡2卡| 怎么达到女性高潮| 成人高潮视频无遮挡免费网站| 999久久久国产精品视频| 国产高清激情床上av| 国产精品九九99| 午夜福利在线观看吧| 欧美丝袜亚洲另类 | 两性午夜刺激爽爽歪歪视频在线观看| 色哟哟哟哟哟哟| 在线免费观看的www视频| 亚洲av熟女| 国产高清有码在线观看视频| 两性夫妻黄色片| 成熟少妇高潮喷水视频| 欧美高清成人免费视频www| 成人av一区二区三区在线看| 久久99热这里只有精品18| 香蕉丝袜av| 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 国产高清视频在线播放一区| 黄色女人牲交| 国产高潮美女av| 一区福利在线观看| 亚洲精品中文字幕一二三四区| 国产成人一区二区三区免费视频网站| 午夜久久久久精精品| 国产成人精品久久二区二区免费| 精品电影一区二区在线| 日韩欧美一区二区三区在线观看| 免费人成视频x8x8入口观看| 免费在线观看成人毛片| 国产精品一区二区免费欧美| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观| 国产av不卡久久| 高清毛片免费观看视频网站| 亚洲国产中文字幕在线视频| 特级一级黄色大片| 国产爱豆传媒在线观看| АⅤ资源中文在线天堂| 成人无遮挡网站| 亚洲av电影不卡..在线观看| 亚洲无线观看免费| 欧美+亚洲+日韩+国产| 女生性感内裤真人,穿戴方法视频| 亚洲精华国产精华精| 欧美日韩亚洲国产一区二区在线观看| 亚洲无线观看免费| 日韩欧美国产一区二区入口| 亚洲av成人av| tocl精华| 免费在线观看日本一区| 亚洲国产高清在线一区二区三| 丰满的人妻完整版| 免费av毛片视频| 嫩草影院入口| 久久伊人香网站| 国产免费男女视频| 国产男靠女视频免费网站| 国产午夜精品久久久久久| 青草久久国产| 日韩欧美国产在线观看| 国产精品98久久久久久宅男小说| 亚洲国产看品久久| 俄罗斯特黄特色一大片| 黄色片一级片一级黄色片| 99久久成人亚洲精品观看| 国产精品98久久久久久宅男小说| 亚洲国产看品久久| 91久久精品国产一区二区成人 | 最近视频中文字幕2019在线8| 国产淫片久久久久久久久| 中文在线观看免费www的网站| av福利片在线观看| 99热网站在线观看| 毛片一级片免费看久久久久| 婷婷色av中文字幕| 国产伦理片在线播放av一区| 久久久午夜欧美精品| 偷拍熟女少妇极品色| 免费看av在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看| 只有这里有精品99| 国产黄片视频在线免费观看| 中文字幕制服av| 免费无遮挡裸体视频| 少妇被粗大猛烈的视频| 色播亚洲综合网| 高清毛片免费看| 99视频精品全部免费 在线| 久久久久久久久大av| 美女cb高潮喷水在线观看| 亚洲精品一区蜜桃| 精品国产一区二区三区久久久樱花 | 最近2019中文字幕mv第一页| 久久久久久久久中文| 午夜日本视频在线| 色综合站精品国产| 大香蕉97超碰在线| 国产激情偷乱视频一区二区| 日韩av不卡免费在线播放| a级毛片免费高清观看在线播放| 午夜爱爱视频在线播放| 日韩成人伦理影院| 免费观看a级毛片全部| 偷拍熟女少妇极品色| videossex国产| 亚洲久久久久久中文字幕| 亚洲人成网站高清观看| 亚洲av二区三区四区| 菩萨蛮人人尽说江南好唐韦庄 | 免费看av在线观看网站| 一个人看视频在线观看www免费| 深夜a级毛片| 欧美激情在线99| 国产在线一区二区三区精 | 少妇被粗大猛烈的视频| 男女视频在线观看网站免费| 晚上一个人看的免费电影| 精品酒店卫生间| 久久久久精品久久久久真实原创| 午夜老司机福利剧场| 日本一本二区三区精品| 国产三级在线视频| 嫩草影院新地址| 内地一区二区视频在线| 黄色日韩在线| 免费电影在线观看免费观看| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区四那| 男女啪啪激烈高潮av片| 亚洲最大成人av| 极品教师在线视频| 日本熟妇午夜| 菩萨蛮人人尽说江南好唐韦庄 | 欧美另类亚洲清纯唯美| 国产精品1区2区在线观看.| 爱豆传媒免费全集在线观看| 美女内射精品一级片tv| 国产高清三级在线| 久久热精品热| 视频中文字幕在线观看| 国产精品久久电影中文字幕| 国产老妇女一区| 一本一本综合久久| 亚洲精品乱码久久久v下载方式| 美女脱内裤让男人舔精品视频| 国产成年人精品一区二区| av女优亚洲男人天堂| 两性午夜刺激爽爽歪歪视频在线观看| 99久久精品热视频| 国产亚洲精品久久久com| 国产淫语在线视频| 亚洲激情五月婷婷啪啪| 视频中文字幕在线观看| 欧美高清性xxxxhd video| 九九在线视频观看精品| 欧美精品国产亚洲| 成人性生交大片免费视频hd| 91狼人影院| 亚洲美女视频黄频| 久久久久久久久中文| 日韩av在线免费看完整版不卡| 一级毛片久久久久久久久女| 亚洲国产精品久久男人天堂| 欧美区成人在线视频| 丰满乱子伦码专区| 亚洲一区高清亚洲精品| 精品不卡国产一区二区三区| 天堂√8在线中文| 欧美3d第一页| 人人妻人人澡人人爽人人夜夜 | 亚洲成色77777| 亚州av有码| 国产麻豆成人av免费视频| 亚洲精华国产精华液的使用体验| 少妇被粗大猛烈的视频| 国产精华一区二区三区| 欧美高清成人免费视频www| 26uuu在线亚洲综合色| 午夜免费激情av| 人人妻人人看人人澡| 日韩中字成人| 国产片特级美女逼逼视频| 国产精品99久久久久久久久| 日本午夜av视频| 亚洲av免费高清在线观看| 亚洲国产日韩欧美精品在线观看| 国产不卡一卡二| 日韩亚洲欧美综合| 国产午夜精品论理片| 日日撸夜夜添| 国产69精品久久久久777片| 一级黄色大片毛片| 淫秽高清视频在线观看| 久久精品国产亚洲av天美| 国产免费一级a男人的天堂| 搡老妇女老女人老熟妇| 国产乱人偷精品视频| 91久久精品电影网| 搡老妇女老女人老熟妇| 午夜精品一区二区三区免费看| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频| 久久亚洲国产成人精品v| 色综合亚洲欧美另类图片| 亚洲欧美一区二区三区国产| 欧美激情久久久久久爽电影| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 看十八女毛片水多多多| 国产亚洲最大av| 久久久久久久久久久丰满| 永久免费av网站大全| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产专区5o | 国产精品女同一区二区软件| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区| av视频在线观看入口| 十八禁国产超污无遮挡网站| 大话2 男鬼变身卡| 嫩草影院新地址| 国产欧美日韩精品一区二区| 国产一区二区在线观看日韩| 大香蕉97超碰在线| 日韩在线高清观看一区二区三区| 村上凉子中文字幕在线| 免费看日本二区| 亚洲中文字幕一区二区三区有码在线看| 国产精品日韩av在线免费观看| 成年av动漫网址| 国产精品综合久久久久久久免费| 综合色丁香网| 久久久a久久爽久久v久久| 亚洲精品成人久久久久久| 日本一二三区视频观看| kizo精华| 久久久午夜欧美精品| .国产精品久久| 婷婷六月久久综合丁香| 午夜老司机福利剧场| 亚洲欧美成人精品一区二区| 亚洲伊人久久精品综合 | 青春草国产在线视频| 精品国产一区二区三区久久久樱花 | 一级毛片我不卡| 欧美+日韩+精品| 男女边吃奶边做爰视频| 久久久久久久久久久丰满| 2022亚洲国产成人精品| 18禁在线无遮挡免费观看视频| 小蜜桃在线观看免费完整版高清| 日韩制服骚丝袜av| 美女高潮的动态| 久久精品人妻少妇| 好男人视频免费观看在线| 91狼人影院| 变态另类丝袜制服| 欧美另类亚洲清纯唯美| 日韩av在线免费看完整版不卡| 尤物成人国产欧美一区二区三区| 2022亚洲国产成人精品| 久久国内精品自在自线图片| 国产精品国产高清国产av| 男女国产视频网站| 国产成人a区在线观看| 国产精品av视频在线免费观看| 久久人人爽人人片av| 人妻夜夜爽99麻豆av| 男女视频在线观看网站免费| 91狼人影院| 久久精品影院6| 天堂影院成人在线观看| 亚洲国产精品sss在线观看| 国产精品一区www在线观看| 熟女电影av网| av视频在线观看入口| 伦精品一区二区三区| 色哟哟·www| 亚洲av成人精品一区久久| 久久亚洲精品不卡| 青春草国产在线视频| 精品人妻偷拍中文字幕| 久久6这里有精品| 亚洲性久久影院| 久久久精品94久久精品| 欧美一级a爱片免费观看看| 久久精品夜夜夜夜夜久久蜜豆| 三级毛片av免费| av.在线天堂| 亚洲五月天丁香| 人妻制服诱惑在线中文字幕| 久久久精品欧美日韩精品| 日韩视频在线欧美| 亚洲av.av天堂| 老司机福利观看| av播播在线观看一区| 国产麻豆成人av免费视频| 免费观看在线日韩| 久久午夜福利片| 亚洲色图av天堂| 两个人视频免费观看高清| 免费av观看视频| 久久精品综合一区二区三区| 国产色婷婷99| 又粗又硬又长又爽又黄的视频| 亚洲天堂国产精品一区在线| www日本黄色视频网| 国产一级毛片在线| 亚洲真实伦在线观看| 成人性生交大片免费视频hd| 精品一区二区三区人妻视频| 午夜免费激情av| 国产高清不卡午夜福利| 日韩av在线大香蕉| 不卡视频在线观看欧美| 18禁在线播放成人免费| 久久久久久久久中文| 一级毛片电影观看 | 在线a可以看的网站| 六月丁香七月| 看十八女毛片水多多多| 岛国在线免费视频观看| 麻豆成人av视频| 少妇裸体淫交视频免费看高清| 久久久久久久久中文| 男女国产视频网站| 99久久精品一区二区三区| av又黄又爽大尺度在线免费看 | 人人妻人人看人人澡| 男女国产视频网站| 秋霞伦理黄片| 亚洲内射少妇av| 午夜精品一区二区三区免费看| 精品午夜福利在线看| 一级爰片在线观看| 亚洲av免费在线观看| 国产成人a区在线观看| 亚洲熟妇中文字幕五十中出| 国产免费又黄又爽又色| 国产一区二区在线av高清观看| 啦啦啦观看免费观看视频高清| 久久精品久久久久久噜噜老黄 | 夜夜看夜夜爽夜夜摸| 亚洲国产最新在线播放| 九九在线视频观看精品| 插阴视频在线观看视频| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 日本-黄色视频高清免费观看| 免费观看的影片在线观看| 精品国产一区二区三区久久久樱花 | 欧美人与善性xxx| 久久久精品94久久精品| 看黄色毛片网站| 亚洲av电影在线观看一区二区三区 | 插逼视频在线观看| 春色校园在线视频观看| 亚洲欧美精品自产自拍| 久久久色成人| 91久久精品电影网| 国产精品一区www在线观看| 欧美精品一区二区大全| 欧美色视频一区免费| 国产精品综合久久久久久久免费| 久99久视频精品免费| 男女啪啪激烈高潮av片| 国产黄色视频一区二区在线观看 | 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 日日干狠狠操夜夜爽| 国产免费男女视频| 91久久精品国产一区二区成人| 麻豆成人av视频| 69人妻影院| 黄色日韩在线| 亚洲最大成人手机在线| 中文天堂在线官网| 国产精品久久久久久精品电影小说 | 日韩高清综合在线| 人人妻人人澡人人爽人人夜夜 | 久久这里只有精品中国| av黄色大香蕉| 亚洲熟妇中文字幕五十中出| 大又大粗又爽又黄少妇毛片口|