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

    Radiomics of rectal cancer for predicting distant metastasis and overall survival

    2020-10-29 02:04:54MouLiYuZhouZhuYongChangZhangYuFengYueHaoPengYuBinSong
    World Journal of Gastroenterology 2020年33期
    關(guān)鍵詞:閃光哲理古詩(shī)文

    Mou Li, Yu-Zhou Zhu, Yong-Chang Zhang, Yu-Feng Yue, Hao-Peng Yu, Bin Song

    Abstract

    Key words: Radiomics; Rectal cancer; Overall survival; Distant metastasis; Computed tomography

    INTRODUCTION

    More than 1.8 million cases of colorectal cancer were diagnosed in 2018 worldwide, making it the third most prevalent malignancy around the world. In terms of morbidity, colorectal cancer accounts for approximately 9.2% of cancer-related deaths worldwide, similar to that of stomach and liver cancer[1]. In up to 70% of patients with rectal cancer (RC), surgical removal of the primary tumor is successful. However, local recurrence and distant metastases are commonly detected in approximately 30% of RC patients, often within 3 years after surgery[2-4]. The overall prognosis becomes poor once distant metastases have developed, thus demonstrating the importance of prompt diagnosis and treatment of RC[5]. Some patients may be at a higher risk of developing adverse outcomes post-surgery. In these patients, alternative and adjunctive therapies, such as chemotherapy, radiotherapy, or other targeted therapies, may be needed to minimize the risk of developing distant metastases[6]. Hence, improved patient outcomes may be feasible by identifying unfavorable prognostic characteristics that could hinder the overall survival of patients. In return, personalized treatment strategies could be implemented to achieved improved outcomes in patients with RC[7]. Recently, the American Society of Clinical Oncology guidelines outlined the need for newer surveillance technologies to better characterize and detect rectal tumors during the early stages of the disease, which could aid in risk stratification and subsequent surveillance[8].

    Currently, there is no universal list of clinicoradiologic prognostic factors used for the detection of RC patients more likely to benefit from treatment. While the TNM classification system is widely utilized for staging cancer in the clinic, it has certain pitfalls that limit its clinical utility in RC[9]. While TNM considers the degree of cancer invasion, the involvement of surrounding lymph nodes, and metastatic spread of the tumor[10], it does not consider the importance of spatial heterogeneity. Spatial heterogeneity is an important characteristic of RC, as it indicates high cell density, hemorrhage, and necrosis[11].

    Radiomics can provide a comprehensive overview of intratumor heterogeneity by extracting multiple quantitative features from medical images. A wide array of parameters can be utilized to assess and quantify the degree of heterogeneity at relevant scales, such as skewness, entropy, kurtosis, and uniformity. Radiomic analysis has been successfully utilized to predict survival outcomes in other diseases[12-14]. When integrated with the clinicopathologic features, the radiomics signature is superior to that of a single biomarker, in terms of prognostic value[15-17]. Hence, in the current study, the aim was to compare the predictive abilities, in terms of distant metastases and 3-year OS, of a radiomics signature and clinicoradiologic risk model in patients with RC.

    MATERIALS AND METHODS

    Patients

    This study was approved by the medical ethics committee of our institution (No. 2019-1159; date: December 26, 2019). Patient approval or informed consent for the review of medical images was waived by the committee due to the retrospective nature of the study.

    The primary cohort of this study comprised an evaluation of the institutional database for medical records from October 2012 to December 2015 to identify patients with histologically-confirmed RC, who underwent curative surgery alone. A total of 148 patients (76 males and 72 females; mean age: 59.7 ± 11.7 years; age range: 27-94 years) were enrolled in our study according to the following inclusion criteria: (1) Patients with histologically-confirmed rectal adenocarcinoma who were treated with curative resection alone; (2) Preoperative enhanced computed tomography (CT) was performed within 1 mo before resection; (3) Follow-up was conducted for at least 3 years; and (4) T-stage, N-stage, and tumor grade were confirmed by pathology. The exclusion criteria were as follows: (1) Patients who received neoadjuvant or postoperative chemoradiotherapy; (2) Patients who presented distant metastasis at the time of diagnosis; (3) Patients who failed to undergo follow-up; (4) Lack of preoperative CT images; and (5) Poor CT image quality. The patients were allocated to a training or validation set at a ratio of 7:3, based on the scanning date. The early data before the 70thpercentile scanning date were allocated to the training set, while the other data were allocated to the validation set. The patient recruitment pathway is shown in Figure 1.

    The baseline clinical characteristics and pathological data of each patient, including age, sex, tumor size on CT, pathological TNM stage, histological grade, antigen Ki-67, circumferential resection margin, number of distant metastasis, and follow-up time were all derived from medical records, as shown in Table 1.

    這一期,我們來(lái)學(xué)習(xí)古詩(shī)文的第九個(gè)字——理。說(shuō)起古詩(shī)文,人們大多感受明顯的是兒女情長(zhǎng),其實(shí),在古詩(shī)文中也有很多哲理的閃光。只是其中需要更高的思想素養(yǎng)才能發(fā)掘,讓我們一起來(lái)發(fā)掘一下古詩(shī)文中的哲理。

    CT examinations

    Enhanced CT was performed using a 128-channel multidetector CT scanner (Simens, SOMATOM Definition AS+) with the following scanning parameters: Tube voltage of 120 kV, tube current of 200-210 mA, and slice thickness of 2.0 mm. With the mass injection (dose: 1.2 mL/kg; injection rate: 3 mL/s) of iodine contrast agent (300 mg/mL), the portal venous phase images (30 s after the trigger) were obtained.

    Reference standard for pathology

    All 148 patients had histopathologically-verified RC, which was diagnosed with resected surgical specimens. The pathological confirmatory report was acquired from electronic medical records. Samples were processed using standard procedures, fixed in formalin, embedded in paraffin, and stained with hematoxylin and eosin (HE).

    Follow-up

    All patients were followed for at least 3 years after surgery (3-mo intervals in the first year; every 6 mo in the following 2 years). The minimum postoperative follow-up period was 1 mo for metastasis (three cases had liver, lung, and bone metastases,respectively), and 11 mo for OS.

    Table 1 Patient and tumor characteristics in the training and validation sets

    Figure 1 Flow chart of the recruitment pathway for patients included in the study. CT: Computed tomography.

    The endpoint of this study was 3-year OS, defined as the time between the surgery and the date of the patient’s death. Local recurrence was defined as recurrence in the pelvis, and distant metastasis as recurrence at sites other than the pelvis. All distant metastatic cases were diagnosed by a multidisciplinary team based on clinical examinations, serum carcinoembryonic antigen levels, chest and abdominopelvic CT, abdominopelvic magnetic resonance imaging, endoscopy, and biopsy. Follow-up information was recorded in the database.

    Texture feature extraction

    A radiologist with 7 years of experience in RC imaging used the ITK-SNAP software (open source, www.itk-snap.org) for the three-dimensional (3D) manual segmentation of the primary tumor, as shown in Figure 2. The region of interest was carefully placed in an area that avoided the intestinal lumen, calcification, blood vessels, and necrosis. In all, 396 quantitative radiomic features of three types (first-order, second-order, and higher-order) were calculated and extracted automatically using the in-house Artificial Intelligence Kit software (GE Healthcare). First-order features are calculated from the pixel intensity histogram to quantify tumor intensity characteristics. Second-order features based on the co-occurrence matrix or run-length matrix account for the location of the pixels and analyze texture in a specific direction. Higher-order features, such as contrast, compare differences and relationships between multiple pixels.

    Model building and evaluation

    Figure 2 A 58-year-old woman with rectal cancer. A-C: Representative manual segmentation of the whole lesion in the axial, sagittal, and coronal enhanced computed tomography images. Red lines represent the delineations of the regions of interest used to derive the radiomic features; D: Three-dimensional volumetric reconstruction of the segmented lesion.

    Before selecting features, the redundant features needed to be eliminated. When the Pearson correlation coefficient of any two features was higher than 0.6, one of them was selected at random, and the rest could be treated in the same manner. Then, the least absolute shrinkage and selection operator (LASSO) regression analysis was performed to identify the most useful prognostic features. The LASSO regression is a variable selection method for linear regression models. The LASSO does this by imposing a constraint on the model parameter (λ) that causes regression coefficients for some variables to shrink toward zero. Variables with a regression coefficient equal to zero are excluded from the model. Variables with non-zero regression coefficients are selected by the LASSO. After that, z-score transformation was used to standardize all selected features. The CT images of 30 patients were segmented twice by the same radiologist to calculate intraclass correlation coefficients for assessing the stability of the features. Radiomic features with an intraclass correlation coefficient > 0.75 were included in the analysis. Each radiomic feature was assessed by univariate logistic regression, and statistical significance was assumed at a confidence level of 0.2 to avoid missing important features. Features withP< 0.2 in univariate logistic regression analysis were entered into the multivariate analysis. A nomogram was generated for model visualization. Receiver operating characteristic curves for each model were constructed to evaluate the performance for predicting distant metastasis of RC. The association of the predicting models with 3-year OS was investigated by Kaplan-Meier survival analysis.

    Statistical analysis

    All statistical analyses were performed using R software (version 3.6.1, R Foundation for Statistical Computing, Vienna, Austria), Statistic Package for Social Science version 21 (Statistic Package for Social Science Inc., Chicago, IL, United States), Stata 15.0 (Stata Corp LLC, Lakeway Drive College Station, United States), and Medcalc 15.2.2 (Medcalc Software, Ostend, Belgium). Differences between the training and validation sets were assessed by the chi-square test, Mann-Whitney test, andt-test. Survival curves were compared by the log-rank test.P< 0.05 was considered statistically significant.

    RESULTS

    A total of 148 patients were enrolled in our study, of whom 51 had distant metastases, 26 died, and 122 survived to the date of the last follow-up visit at 3 years post-surgery (Table 2).

    Table 2 Comparison of the models between the metastasis group and the non-metastasis group, and between the death group and the survival group

    Feature selection and model building

    A total of 17 radiomic features remained after removing the redundant features. Two features, “MaxIntensity” and “RelativeDeviation”, were selected by the LASSO method from the remaining features (Figure 3). To avoid missing important features, another feature “Inertia_AllDirection_offset7_SD” was selected by the univariate logistic regression (P= 0.192). These three features were presented in this calculation formula: Rad-score = 1/{1 + EXP [-(-0.649 - 0.786 × MaxIntensity - 0.868 × Inertia_AllDirection_offset7_SD - 2.743 × RelativeDeviation)]}. The radiomics signature (Rad-score), which consisted of three selected features, was statistically different between the metastasis and non-metastasis groups (Rad-score = 0.46 and 0.32, respectively;P< 0.001), as shown in Table 2.

    For clinical features, pathological T-stage and N-stage were selected by the univariate logistic regression analysis (P< 0.001 andP= 0.002, respectively). The clinicoradiologic risk model was built by the multivariate logistic regression analysis, using the formula 1/{1 + EXP [-(-6.199 + 1.368 × T-stage + 0.738 × N-stage + 4.673 × Rad-score)]}. The results of the combined model were significantly different between the metastasis and non-metastasis groups by 3 years after surgery (P< 0.001), as shown in Table 2. The Hosmer-Lemeshow test yielded no statistical significance (P= 0.280 and 0.151 for the Rad-score and combined model, respectively), which suggested that there was no evidence of lack of fit.

    A nomogram to predict the distant metastasis of RC was generated for the combined model visualization (Figure 4). The nomogram could be used by locating the score for each variable on the corresponding axis, adding the points together for all of the variables, and drawing a line from the total number of points to the risk axis, which allows for the determination of distant metastasis risk. Higher total scores were associated with a greater risk of distant metastasis. The contributions of these variables to the Rad-score and the combined model were measured by the value of the standardized logistic regression coefficient (Figure 5). The contribution of “MaxIntensity” to the Rad-score, and that of “T-stage” to the combined model were greater than those of the others.

    Classification results

    Figure 3 Feature selection using the least absolute shrinkage and selection operator regression model. Dotted vertical lines were drawn at the optimal values (two features were selected) by using the minimum criteria and one standard error of the minimum criteria (i.e., the 1-SE criteria).

    Figure 4 Developed radiomics nomogram. The radiomics nomogram was developed in the training cohort, with the Rad-score, T-stage, and N-stage incorporated.

    In the case of only considering the Rad-score, the resulting area under the curve (AUC) was 0.709 (95%CI: 0.612-0.795). Improved prediction of distant metastases could be achieved by combining the radiomic features from preoperative CT images with the clinical features. The AUCs of the clinical model (T-stage combined with N-stage) and the combined model were 0.782 (95%CI: 0.689-0.857) and 0.842 (95%CI: 0.757-0.906), respectively (Table 3 and Figure 6). There was a significant difference in AUCs between the Rad-score and combined model in the training cohort (AUC = 0.709vs0.842,P= 0.005), which was confirmed using the validation cohort (AUC = 0.687vs0.802,P= 0.020). We found that the clinical variables had higher classification contributions than Rad-score to build the combined model. This finding was consistent with the higher standardized logistic regression coefficient of T-stage than that of the Rad-score (Figure 5). Moreover, the AUC value of the clinical model was also higherthan that of the Rad-score. Even so, the Rad-score greatly helped to increase the AUC value from 0.782 to 0.842 (Table 3 and Figure 6). In the subgroups of the receiver operating characteristic analysis, the AUC of the combined model decreased in the 1-year (distant metastasis occurred in the first year after surgery) and 2-year (distant metastasis occurred within 2 years after surgery) subgroups (Table 3).

    Table 3 Diagnostic performance of the three models in both training and validation cohorts

    Figure 5 Histogram showing the contribution of each variable to the models. A: the Rad-score; B: the combined model. The contributions of the variables were measured by the values of the standardized logistic regression coefficients.

    Assessment of the combined model for OS

    Using the clinicoradiologic risk model, stratified analyses were performed for the whole set, training set, and validation set to evaluate the association with OS. As shown in Figure 7, the grouping results of the combined model were significantly associated with OS in the whole, training, and validation groups (P< 0.0001,P= 0.0001, andP= 0.0137, respectively). Although there was a lack of statistical significance for stage III in the stratified subgroup analysis, according to the overall pathological stage, the low-risk group displayed a longer OS than the high-risk group (P< 0.001,P= 0.0194, andP= 0.1401 for stages I, II, and III, respectively), which was significant in terms of individualized treatment (Figure 8)[18].

    Figure 6 Receiver operating characteristic curves of the models in the training and validation cohorts. A: In the training set, the combined model [area under the curve under the curve (AUC) = 0.842] achieved a better performance than the Rad-score (AUC = 0.709) and the clinical model containing Tstage and N-stage (AUC = 0.782); B: In the validation set, the AUCs of the Rad-score, the clinical model, and the combined model were 0.687, 0.766, and 0.802, respectively.

    Figure 7 Kaplan-Meier curves for overall survival stratified by risk group, as identified by the combined model. A: The whole set. P < 0.0001; B: The training set. P = 0.0001; C: The validation set. P = 0.0137.

    Figure 8 Overall survival curves for the low- and high-risk groups classified according to the combined model in the subgroups of the overall pathological stage. A: Stage I. P < 0.001; B: Stage II. P = 0.0194; C: Stage III. P = 0.1401.

    DISCUSSION

    In this study, we developed and validated a clinicoradiologic risk model that showed potential for predicting distant metastasis of RC within 3 years after surgery, which had a better prognostic performance than Rad-score (P= 0.005 and 0.020 for the training cohort and validation cohort, respectively). The combined model was used to stratify patients into low-risk and high-risk groups for the analysis of 3-year OS. The results showed that OS rates between low- and high-risk groups were significantly different in the training cohort, which was verified in the validation cohort.

    Compared with traditional image explanation, which is qualitative or subjective, radiomic analyses permit high-throughput extraction of radiomic features that can quantify differences between tissues invisible to human eyes. Recently, the use of radiomics has appeared as a potential technique for constructing decision-support models based on high-throughput quantificational characters extracted from medical images. Radiomics-based prognosis prediction models have been reported for advanced nasopharyngeal carcinoma[19], early-stage non-small cell lung cancer[20], and RC[12-14]. In their pioneering retrospective studies about the prognosis of RC[12-14], robust models and strong independent prognostic factors have been developed for the prediction of OS in patients with RC. Our results were consistent with previous studies, suggesting that radiomics could help predict the prognosis of patients with RC. The building methods between our study and those of previous studies were similar (machine learning), and all these studies lacked external validation. However, there were still some differences that need to be explained. For the follow-up time, Wanget al[13]followed the patients for 5 years, and Lovinfosseet al[12]for 4 years, which were longer than the time in the current study. These authors[12-14]focused on predicting the prognosis in locally advanced RC patients treated with neoadjuvant chemoradiation followed by surgery, which is different from our current study. Different chemoradiotherapy regimens might influence the prognostic evaluation of patients with RC. Moreover, the sample sizes of these two studies[12,14]were less than ours.

    In terms of feature selection, several relationships were uncovered between the specific features and their ability to predict distant metastasis. For the radiomic features included in this study, the first-order features show an excellent auxiliary classification effect, accounting for two of the three specified radiomic features. In the clinical features, we found that T-stage and N-stage had high standardized logistic regression coefficients, indicating the contribution to the combined model. This finding is consistent with the consensus that pathological T-stage and N-stage are very important for predicting the distant metastasis of RC. The histological grade, circumferential resection margin, Ki-67 score, and tumor diameter did not present enough predictive power for distant metastasis and prognosis. Subsequently, we integrated the Rad-score into a nomogram with clinical risk factors, and constructed a useful tool for individualized evaluation of distant metastasis and OS in patients with RC.

    In our study, only patients who did not receive chemoradiotherapy were selected for the analysis of prognosis. On the one hand, the variables caused by the different intensive chemoradiotherapy methods can be controlled. On the other hand, we could build a novel model to stratify the high-risk patients. This was in step with the current trends toward personalized medicine[21]. Considering individualized evaluation of patients with RC of different stages, the subgroup survival analysis was performed. The low-risk groups had longer OS than the high-risk groups, which was statically significant between stage I and II (P< 0.001 andP= 0.0194, respectively). However, there was no statistically significant difference for stage III (P= 0.1401). We speculate that this might be related to the small sample size of this study.

    There are several limitations to this study. The first limitation is the relatively small sample size. Second, even if we controlled therapeutic methods to surgery alone, inevitable bias may exist due to the retrospective design of this study. Therefore, prospective and external validation studies are required in future studies. Third, all of the CT images were obtained from a single institution. In the future, multicenter verification is necessary to extend the versatility of the experimental results.

    CONCLUSION

    In conclusion, this study describes a combined model that incorporates a radiomics signature and clinical risk factors. The model can aid in the individualized prediction of distant metastasis and prognosis in patients with RC.

    ARTICLE HIGHLIGHTS

    Research results

    A total of 51 (34.5%) patients had distant metastases, while 26 (17.6%) patients died, and 122 (82.4%) patients lived at least 3 years post-surgery. The values of both the Rad-score and the combined model were significantly different between the distant metastasis group and the non-metastasis group (0.46 ± 0.21vs0.32 ± 0.24 for the Radscore, 0.60 ± 0.23vs0.28 ± 0.26 for the combined model;P< 0.001 for both models). Predictors contained in the combined model included the Rad-score, pathological Nstage, and T-stage. The combined model showed good discrimination, with an area under the curve of 0.842 and 0.802 for the training set and validation set, respectively. For the survival analysis, the combined model was associated with an improved OS in the whole cohort and the respective subgroups.

    Research conclusions

    This study presents a novel model, visualized in a nomogram, that can be used to facilitate individualized prediction of distant metastasis and 3-year OS in patients with RC.

    Research perspectives

    Radiomics may change the practice of medicine, particularly for patients with RC. However, there are challenges to be overcome before its routine implementation including challenges related to sample size, model design, and the lack of robust multicenter validation set. Therefore, prospective multicenter studies of a larger size are needed to externally validate our proposed model in the future.

    猜你喜歡
    閃光哲理古詩(shī)文
    閃光
    閃光的枝條
    古詩(shī)文閱讀備考指津
    哲理漫畫
    特別文摘(2016年18期)2016-09-26 18:02:20
    漫畫哲理
    雜文選刊(2016年9期)2016-09-14 19:52:43
    漫畫哲理
    雜文選刊(2016年5期)2016-05-12 20:10:00
    哲理漫畫
    特別文摘(2016年5期)2016-05-04 22:14:40
    引領(lǐng)小學(xué)生誦讀古詩(shī)文之妙招
    八月,紀(jì)念碑在閃光
    海峽姐妹(2015年8期)2015-02-27 15:12:54
    對(duì)古詩(shī)文默寫有效性的思考
    啦啦啦 在线观看视频| 男女无遮挡免费网站观看| 可以免费在线观看a视频的电影网站| 欧美精品亚洲一区二区| 成年人黄色毛片网站| 免费女性裸体啪啪无遮挡网站| 亚洲av美国av| 日本91视频免费播放| 热re99久久国产66热| 久久久久国产一级毛片高清牌| 亚洲一卡2卡3卡4卡5卡精品中文| 纵有疾风起免费观看全集完整版| 一二三四社区在线视频社区8| 久久国产精品男人的天堂亚洲| 午夜老司机福利片| 操出白浆在线播放| 国产成人免费观看mmmm| 精品久久蜜臀av无| 99久久精品国产亚洲精品| 亚洲色图 男人天堂 中文字幕| 天堂中文最新版在线下载| 精品一区在线观看国产| 亚洲国产欧美在线一区| 亚洲成人国产一区在线观看 | 久久狼人影院| 在线观看www视频免费| 国产精品99久久99久久久不卡| 国产精品免费视频内射| 午夜免费成人在线视频| 亚洲成人免费av在线播放| 青春草视频在线免费观看| 国产真人三级小视频在线观看| av电影中文网址| 女性被躁到高潮视频| 国产精品熟女久久久久浪| 建设人人有责人人尽责人人享有的| 丝袜人妻中文字幕| 久久久久精品人妻al黑| 性高湖久久久久久久久免费观看| 精品久久久精品久久久| 亚洲国产最新在线播放| 一级毛片黄色毛片免费观看视频| 国产免费又黄又爽又色| 男人添女人高潮全过程视频| av视频免费观看在线观看| 久久久欧美国产精品| 两个人免费观看高清视频| av国产久精品久网站免费入址| 老汉色av国产亚洲站长工具| 欧美成人精品欧美一级黄| 免费日韩欧美在线观看| 国产男女超爽视频在线观看| 国产成人精品久久久久久| 欧美日韩综合久久久久久| 亚洲人成网站在线观看播放| 国产成人精品在线电影| 一二三四在线观看免费中文在| 纯流量卡能插随身wifi吗| 色婷婷av一区二区三区视频| 国产福利在线免费观看视频| 汤姆久久久久久久影院中文字幕| 极品人妻少妇av视频| 欧美人与性动交α欧美精品济南到| 老鸭窝网址在线观看| 亚洲欧美一区二区三区久久| 又粗又硬又长又爽又黄的视频| 激情五月婷婷亚洲| 欧美精品av麻豆av| 青春草亚洲视频在线观看| 久久久精品94久久精品| 亚洲av国产av综合av卡| 国产精品久久久久久人妻精品电影 | 亚洲精品国产av蜜桃| 啦啦啦视频在线资源免费观看| 热99国产精品久久久久久7| 另类精品久久| 欧美av亚洲av综合av国产av| 精品亚洲成国产av| 色视频在线一区二区三区| 欧美日韩视频高清一区二区三区二| svipshipincom国产片| 91精品伊人久久大香线蕉| 久久久久网色| 一级黄片播放器| 三上悠亚av全集在线观看| 久久人人爽人人片av| 999精品在线视频| 亚洲成人免费av在线播放| 欧美老熟妇乱子伦牲交| 看十八女毛片水多多多| 久久99一区二区三区| 久久九九热精品免费| 精品一区二区三区四区五区乱码 | 人人妻人人澡人人看| 色综合欧美亚洲国产小说| 亚洲国产精品一区三区| 人妻人人澡人人爽人人| 天堂俺去俺来也www色官网| 十八禁高潮呻吟视频| 久热这里只有精品99| 午夜免费男女啪啪视频观看| 亚洲国产精品999| 男女之事视频高清在线观看 | 99久久精品国产亚洲精品| 夫妻午夜视频| 两个人看的免费小视频| 91精品伊人久久大香线蕉| 亚洲自偷自拍图片 自拍| 中文字幕精品免费在线观看视频| 亚洲欧美一区二区三区久久| 国产成人91sexporn| 国产又爽黄色视频| 国产99久久九九免费精品| 日韩人妻精品一区2区三区| 色网站视频免费| 97在线人人人人妻| 一区二区三区精品91| 在线观看免费午夜福利视频| 操出白浆在线播放| 国产一区亚洲一区在线观看| 一本—道久久a久久精品蜜桃钙片| 男女边吃奶边做爰视频| 人妻人人澡人人爽人人| 国产精品一二三区在线看| 新久久久久国产一级毛片| 久久精品aⅴ一区二区三区四区| 成年人免费黄色播放视频| 精品少妇黑人巨大在线播放| 精品人妻在线不人妻| 国产97色在线日韩免费| 男女边摸边吃奶| 亚洲欧美精品自产自拍| 亚洲国产精品一区二区三区在线| 国产高清不卡午夜福利| 在线天堂中文资源库| 十八禁网站网址无遮挡| 亚洲中文av在线| 久久精品国产亚洲av高清一级| 日韩大码丰满熟妇| 香蕉丝袜av| 日韩视频在线欧美| 亚洲av美国av| 永久免费av网站大全| 一级片免费观看大全| 国产精品一国产av| 亚洲欧洲日产国产| av电影中文网址| 91字幕亚洲| 亚洲精品美女久久久久99蜜臀 | 欧美日韩亚洲国产一区二区在线观看 | 又大又黄又爽视频免费| 美女大奶头黄色视频| 深夜精品福利| 一级片免费观看大全| 欧美人与性动交α欧美精品济南到| 中文字幕色久视频| 免费在线观看日本一区| 少妇 在线观看| 久久久国产欧美日韩av| 99re6热这里在线精品视频| 欧美亚洲 丝袜 人妻 在线| 免费av中文字幕在线| 国产精品久久久久久精品电影小说| 久久狼人影院| 一级,二级,三级黄色视频| av一本久久久久| 午夜激情av网站| 久热这里只有精品99| 日韩一本色道免费dvd| 丰满人妻熟妇乱又伦精品不卡| 久热爱精品视频在线9| 1024香蕉在线观看| 欧美精品一区二区免费开放| 两个人看的免费小视频| 悠悠久久av| 久久久国产精品麻豆| 又大又黄又爽视频免费| av国产久精品久网站免费入址| 飞空精品影院首页| videosex国产| 人人妻人人澡人人看| 天天躁日日躁夜夜躁夜夜| 一个人免费看片子| 国产成人精品久久久久久| 亚洲国产av新网站| 久久国产精品人妻蜜桃| 亚洲熟女毛片儿| 一二三四社区在线视频社区8| 欧美日韩一级在线毛片| xxxhd国产人妻xxx| 999久久久国产精品视频| 女性生殖器流出的白浆| 99香蕉大伊视频| 丝袜美足系列| 自线自在国产av| 亚洲成色77777| 国产黄色免费在线视频| 午夜久久久在线观看| 亚洲一区中文字幕在线| 国产野战对白在线观看| 色网站视频免费| 亚洲国产av影院在线观看| av国产精品久久久久影院| 婷婷色av中文字幕| 男女床上黄色一级片免费看| 只有这里有精品99| 中国国产av一级| 99热全是精品| 欧美变态另类bdsm刘玥| 日韩av在线免费看完整版不卡| 亚洲综合色网址| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| 欧美日韩亚洲综合一区二区三区_| 日本vs欧美在线观看视频| 亚洲国产精品一区三区| 国产视频一区二区在线看| 少妇被粗大的猛进出69影院| 免费看十八禁软件| 久久国产精品人妻蜜桃| 国产精品一二三区在线看| 午夜福利,免费看| 香蕉国产在线看| 久久久久久久精品精品| 日韩中文字幕欧美一区二区 | 多毛熟女@视频| 无限看片的www在线观看| 欧美亚洲日本最大视频资源| 亚洲国产日韩一区二区| 国产一区二区三区av在线| 夫妻午夜视频| 最近手机中文字幕大全| 99久久人妻综合| 精品人妻在线不人妻| 19禁男女啪啪无遮挡网站| 久久免费观看电影| 各种免费的搞黄视频| 精品久久久精品久久久| 久久精品亚洲熟妇少妇任你| 亚洲图色成人| 亚洲精品美女久久av网站| 这个男人来自地球电影免费观看| 大码成人一级视频| 一本久久精品| 人人妻,人人澡人人爽秒播 | 成年人免费黄色播放视频| 国产黄频视频在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 国产免费福利视频在线观看| 精品久久久久久电影网| 最近手机中文字幕大全| 国产精品.久久久| 国产精品人妻久久久影院| 大片免费播放器 马上看| 免费黄频网站在线观看国产| 香蕉国产在线看| 女人爽到高潮嗷嗷叫在线视频| 男女午夜视频在线观看| 少妇裸体淫交视频免费看高清 | 人人妻人人爽人人添夜夜欢视频| 777米奇影视久久| 亚洲国产欧美日韩在线播放| 男女边摸边吃奶| 好男人视频免费观看在线| 国产精品香港三级国产av潘金莲 | 丝袜喷水一区| 两性夫妻黄色片| 国产日韩欧美亚洲二区| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久av美女十八| bbb黄色大片| 亚洲欧美日韩另类电影网站| 高清黄色对白视频在线免费看| 看免费av毛片| 亚洲,欧美精品.| 国产97色在线日韩免费| 91字幕亚洲| 亚洲av国产av综合av卡| 国产精品熟女久久久久浪| 19禁男女啪啪无遮挡网站| 国产精品 欧美亚洲| 国产在线一区二区三区精| 十八禁人妻一区二区| 亚洲自偷自拍图片 自拍| 久久久久精品人妻al黑| 成人18禁高潮啪啪吃奶动态图| 在线观看人妻少妇| 精品熟女少妇八av免费久了| 在线观看www视频免费| 成年美女黄网站色视频大全免费| 男男h啪啪无遮挡| 熟女av电影| 亚洲精品乱久久久久久| av天堂在线播放| 国产色视频综合| 一边摸一边抽搐一进一出视频| 一二三四社区在线视频社区8| 高清av免费在线| 黑人巨大精品欧美一区二区蜜桃| 午夜av观看不卡| 欧美日韩精品网址| 亚洲人成77777在线视频| 国产高清国产精品国产三级| 成年动漫av网址| 国产一区二区在线观看av| 激情视频va一区二区三区| 男女下面插进去视频免费观看| 99久久综合免费| 亚洲国产精品999| 午夜激情久久久久久久| 国产成人系列免费观看| 亚洲国产av影院在线观看| 超碰97精品在线观看| 母亲3免费完整高清在线观看| 日日爽夜夜爽网站| 99国产精品免费福利视频| 亚洲欧洲日产国产| 日韩电影二区| 精品久久蜜臀av无| 国产亚洲av片在线观看秒播厂| 国产精品成人在线| 777久久人妻少妇嫩草av网站| 亚洲av电影在线进入| 国产精品一区二区在线不卡| 国产真人三级小视频在线观看| 国产一区二区在线观看av| 蜜桃国产av成人99| 久久性视频一级片| 国产精品国产三级国产专区5o| 男女之事视频高清在线观看 | 蜜桃在线观看..| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕制服av| 欧美亚洲 丝袜 人妻 在线| 亚洲午夜精品一区,二区,三区| 最近中文字幕2019免费版| av国产精品久久久久影院| 91麻豆精品激情在线观看国产 | 日日爽夜夜爽网站| 在线天堂中文资源库| 国产日韩一区二区三区精品不卡| 国产av一区二区精品久久| 久久久久久亚洲精品国产蜜桃av| 午夜免费成人在线视频| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 制服诱惑二区| 国产亚洲一区二区精品| 免费高清在线观看视频在线观看| 久久久久国产精品人妻一区二区| a级毛片在线看网站| 2021少妇久久久久久久久久久| 女性生殖器流出的白浆| 免费在线观看日本一区| 国产精品二区激情视频| 在线观看国产h片| 美女午夜性视频免费| 十分钟在线观看高清视频www| 久久久久久人人人人人| 日韩伦理黄色片| 国产一级毛片在线| 观看av在线不卡| 亚洲精品国产区一区二| 如日韩欧美国产精品一区二区三区| 狂野欧美激情性bbbbbb| 亚洲国产欧美一区二区综合| 国产精品久久久久久精品电影小说| 国产亚洲精品久久久久5区| 99九九在线精品视频| 国产亚洲欧美在线一区二区| 精品高清国产在线一区| 老司机影院毛片| 热re99久久国产66热| 国产欧美亚洲国产| 啦啦啦中文免费视频观看日本| 深夜精品福利| 一级毛片女人18水好多 | 亚洲中文日韩欧美视频| 精品久久久久久久毛片微露脸 | 国语对白做爰xxxⅹ性视频网站| 视频区欧美日本亚洲| 啦啦啦视频在线资源免费观看| 日本黄色日本黄色录像| 黄色毛片三级朝国网站| 国产免费现黄频在线看| 亚洲av片天天在线观看| 黑丝袜美女国产一区| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 精品少妇内射三级| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区久久| 欧美日韩综合久久久久久| 精品国产国语对白av| 视频区欧美日本亚洲| 国产免费视频播放在线视频| 亚洲人成77777在线视频| 亚洲成人手机| 国产免费现黄频在线看| 天天躁日日躁夜夜躁夜夜| 国产高清不卡午夜福利| 亚洲av片天天在线观看| 精品国产乱码久久久久久男人| 咕卡用的链子| 国产成人精品久久久久久| 欧美国产精品va在线观看不卡| 欧美xxⅹ黑人| 9热在线视频观看99| 99国产精品一区二区蜜桃av | av网站免费在线观看视频| 男人爽女人下面视频在线观看| 亚洲国产精品一区二区三区在线| 精品福利观看| 只有这里有精品99| 欧美成人精品欧美一级黄| 99香蕉大伊视频| 午夜精品国产一区二区电影| 老司机午夜十八禁免费视频| 国产成人系列免费观看| www.精华液| 亚洲人成网站在线观看播放| 制服人妻中文乱码| 国产伦理片在线播放av一区| 精品少妇久久久久久888优播| 男女床上黄色一级片免费看| 免费一级毛片在线播放高清视频 | 精品亚洲成国产av| 午夜激情久久久久久久| 男人舔女人的私密视频| 国产97色在线日韩免费| 精品视频人人做人人爽| av天堂久久9| 一本一本久久a久久精品综合妖精| 精品一区二区三区av网在线观看 | 最近手机中文字幕大全| cao死你这个sao货| 麻豆乱淫一区二区| 如日韩欧美国产精品一区二区三区| 久久九九热精品免费| 午夜激情av网站| av网站免费在线观看视频| 九色亚洲精品在线播放| 久久99精品国语久久久| 中文字幕亚洲精品专区| 日本wwww免费看| 午夜福利视频在线观看免费| 99热国产这里只有精品6| 久久久精品区二区三区| 国产福利在线免费观看视频| 亚洲色图综合在线观看| 国产一区有黄有色的免费视频| 制服诱惑二区| 久久人人爽人人片av| 中文字幕人妻丝袜一区二区| 婷婷丁香在线五月| 黑人巨大精品欧美一区二区蜜桃| 成年人黄色毛片网站| 国产黄色视频一区二区在线观看| 男女免费视频国产| av国产久精品久网站免费入址| 亚洲,欧美,日韩| 制服人妻中文乱码| 人妻 亚洲 视频| 777米奇影视久久| av在线老鸭窝| 亚洲国产av新网站| 女人被躁到高潮嗷嗷叫费观| 久久精品亚洲av国产电影网| 久久精品亚洲熟妇少妇任你| 一级毛片黄色毛片免费观看视频| 9191精品国产免费久久| 久久久国产精品麻豆| 日日爽夜夜爽网站| 欧美黑人欧美精品刺激| 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 精品欧美一区二区三区在线| 纯流量卡能插随身wifi吗| 波多野结衣一区麻豆| 一级黄片播放器| 老鸭窝网址在线观看| 丰满少妇做爰视频| 90打野战视频偷拍视频| 久久99热这里只频精品6学生| 久久毛片免费看一区二区三区| 国产精品人妻久久久影院| 久久久国产欧美日韩av| 亚洲伊人久久精品综合| 丝瓜视频免费看黄片| 99精品久久久久人妻精品| 亚洲av综合色区一区| 91国产中文字幕| 一本久久精品| 亚洲国产日韩一区二区| www.精华液| 久久久国产欧美日韩av| 视频区图区小说| 精品少妇一区二区三区视频日本电影| 亚洲欧美精品自产自拍| 国产av国产精品国产| 免费看不卡的av| 各种免费的搞黄视频| 女人被躁到高潮嗷嗷叫费观| 国产一区二区三区av在线| www日本在线高清视频| 色94色欧美一区二区| xxxhd国产人妻xxx| 久久毛片免费看一区二区三区| 午夜激情久久久久久久| tube8黄色片| 国产黄色免费在线视频| 免费看十八禁软件| 国产成人免费无遮挡视频| 国产淫语在线视频| 狂野欧美激情性xxxx| av天堂在线播放| 国产国语露脸激情在线看| 水蜜桃什么品种好| 天堂中文最新版在线下载| av电影中文网址| 欧美日韩亚洲综合一区二区三区_| 热re99久久国产66热| 亚洲视频免费观看视频| 一边亲一边摸免费视频| 亚洲少妇的诱惑av| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 成人亚洲欧美一区二区av| 色视频在线一区二区三区| 久久毛片免费看一区二区三区| 啦啦啦在线观看免费高清www| 女人精品久久久久毛片| 老司机深夜福利视频在线观看 | 在线亚洲精品国产二区图片欧美| 麻豆国产av国片精品| 超碰成人久久| 亚洲七黄色美女视频| 亚洲一区中文字幕在线| 大型av网站在线播放| 国产片特级美女逼逼视频| cao死你这个sao货| 国产精品99久久99久久久不卡| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 汤姆久久久久久久影院中文字幕| 日韩一区二区三区影片| 免费看十八禁软件| 亚洲成人免费电影在线观看 | 国产精品一二三区在线看| 亚洲激情五月婷婷啪啪| 精品福利观看| 国产成人精品无人区| 热99国产精品久久久久久7| 久久国产精品大桥未久av| 美女福利国产在线| 熟女av电影| 成人18禁高潮啪啪吃奶动态图| a级片在线免费高清观看视频| 成年av动漫网址| 亚洲美女黄色视频免费看| 国产成人欧美| 欧美精品人与动牲交sv欧美| 日韩免费高清中文字幕av| 欧美日韩综合久久久久久| av国产精品久久久久影院| 性高湖久久久久久久久免费观看| 一区二区三区四区激情视频| 国产欧美亚洲国产| 大片免费播放器 马上看| 99精国产麻豆久久婷婷| 我要看黄色一级片免费的| 精品免费久久久久久久清纯 | 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 久久精品aⅴ一区二区三区四区| 国产一区亚洲一区在线观看| 亚洲欧美色中文字幕在线| 成人国产av品久久久| 亚洲中文字幕日韩| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av涩爱| www.熟女人妻精品国产| 亚洲欧美色中文字幕在线| 精品欧美一区二区三区在线| 国产淫语在线视频| www.精华液| 欧美日韩视频精品一区| 久久久久久久大尺度免费视频| 一区二区av电影网| 国产在线一区二区三区精| 国产一区亚洲一区在线观看| 精品视频人人做人人爽| 两个人看的免费小视频| 成年人免费黄色播放视频| 国产1区2区3区精品| 国产高清不卡午夜福利| 精品卡一卡二卡四卡免费| 男女之事视频高清在线观看 | 亚洲九九香蕉| 日韩中文字幕欧美一区二区 | 1024香蕉在线观看| 一本大道久久a久久精品| 日本午夜av视频| 欧美黑人精品巨大| 亚洲国产最新在线播放| 在线观看免费日韩欧美大片| 国产在线观看jvid| 啦啦啦在线免费观看视频4| 色精品久久人妻99蜜桃| av在线app专区| 亚洲精品久久午夜乱码| av网站免费在线观看视频|