• <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ī)文默寫有效性的思考
    变态另类丝袜制服| 国产v大片淫在线免费观看| 亚洲成人精品中文字幕电影| 自拍偷自拍亚洲精品老妇| 亚洲真实伦在线观看| 日本黄大片高清| 国产极品精品免费视频能看的| 色综合色国产| 搞女人的毛片| 欧美不卡视频在线免费观看| 天天一区二区日本电影三级| 国产精品一区二区免费欧美| 国产 一区 欧美 日韩| 亚洲人与动物交配视频| 精品一区二区三区视频在线| 亚洲国产色片| 中文字幕精品亚洲无线码一区| 亚洲高清免费不卡视频| 观看美女的网站| 亚洲av成人av| 九九久久精品国产亚洲av麻豆| 日韩精品有码人妻一区| 精品一区二区免费观看| 亚洲av美国av| 亚洲精品国产成人久久av| 日韩成人av中文字幕在线观看 | 美女黄网站色视频| 亚洲熟妇熟女久久| 亚洲内射少妇av| 亚洲精品粉嫩美女一区| 噜噜噜噜噜久久久久久91| 国产精品爽爽va在线观看网站| av在线播放精品| 男女视频在线观看网站免费| 国产精品嫩草影院av在线观看| 97热精品久久久久久| 亚洲av免费高清在线观看| 最近在线观看免费完整版| 人妻少妇偷人精品九色| 午夜激情福利司机影院| 久久精品国产亚洲网站| 成年av动漫网址| 观看免费一级毛片| 亚洲成av人片在线播放无| 黄色一级大片看看| 伊人久久精品亚洲午夜| 午夜福利在线观看免费完整高清在 | 有码 亚洲区| 免费在线观看影片大全网站| 国产精品久久久久久亚洲av鲁大| 国产精品不卡视频一区二区| 国产在线男女| 精品久久久久久久久久免费视频| 国产色婷婷99| 欧美日本视频| 亚洲人成网站在线观看播放| 日本一本二区三区精品| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 赤兔流量卡办理| 久久久久免费精品人妻一区二区| 黄色日韩在线| 亚洲中文字幕日韩| 一区二区三区高清视频在线| 日韩av不卡免费在线播放| 身体一侧抽搐| 午夜福利高清视频| 久久久午夜欧美精品| 免费看日本二区| 搡女人真爽免费视频火全软件 | 搡老妇女老女人老熟妇| 日韩成人av中文字幕在线观看 | 亚洲av中文字字幕乱码综合| 欧美人与善性xxx| 人人妻人人澡人人爽人人夜夜 | 特级一级黄色大片| 国产毛片a区久久久久| 久久久久久久久中文| 国产伦在线观看视频一区| 在线a可以看的网站| 亚洲电影在线观看av| 欧美一区二区精品小视频在线| 亚洲国产精品sss在线观看| 亚洲av.av天堂| 深夜精品福利| 狂野欧美激情性xxxx在线观看| 91在线精品国自产拍蜜月| 亚洲精品影视一区二区三区av| 精品一区二区三区视频在线| 亚洲欧美精品综合久久99| 乱码一卡2卡4卡精品| 在线播放无遮挡| 一级黄色大片毛片| 日本免费一区二区三区高清不卡| av在线蜜桃| 此物有八面人人有两片| 少妇丰满av| 久久久成人免费电影| 免费高清视频大片| a级毛片免费高清观看在线播放| 日本爱情动作片www.在线观看 | 欧美激情国产日韩精品一区| 高清午夜精品一区二区三区 | 亚洲乱码一区二区免费版| 久久亚洲精品不卡| 久久中文看片网| 日韩强制内射视频| 毛片一级片免费看久久久久| 国产欧美日韩一区二区精品| 高清毛片免费观看视频网站| 精品一区二区三区视频在线| 成人鲁丝片一二三区免费| 国产大屁股一区二区在线视频| 亚洲中文日韩欧美视频| 波野结衣二区三区在线| 国产精品久久久久久亚洲av鲁大| 少妇丰满av| 色哟哟·www| 最近中文字幕高清免费大全6| 午夜福利视频1000在线观看| 午夜爱爱视频在线播放| 特大巨黑吊av在线直播| 亚洲久久久久久中文字幕| 真人做人爱边吃奶动态| 婷婷亚洲欧美| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩无卡精品| 国产精品人妻久久久影院| 黄色视频,在线免费观看| av视频在线观看入口| 国产真实乱freesex| 人妻久久中文字幕网| 99热精品在线国产| 精品不卡国产一区二区三区| 国产男靠女视频免费网站| 免费av毛片视频| а√天堂www在线а√下载| 午夜老司机福利剧场| 精品人妻熟女av久视频| 精品人妻一区二区三区麻豆 | 此物有八面人人有两片| 国产亚洲精品久久久久久毛片| 最近视频中文字幕2019在线8| 岛国在线免费视频观看| 国产精品一区www在线观看| 天堂影院成人在线观看| 久久精品国产亚洲av涩爱 | 亚洲最大成人中文| 国产女主播在线喷水免费视频网站 | 国产av麻豆久久久久久久| 亚洲专区国产一区二区| 日本免费一区二区三区高清不卡| 亚洲成人精品中文字幕电影| 日韩欧美三级三区| 亚洲欧美日韩卡通动漫| 国产视频内射| 天堂av国产一区二区熟女人妻| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 成人特级av手机在线观看| 一本久久中文字幕| 色尼玛亚洲综合影院| 亚洲欧美成人精品一区二区| 夜夜看夜夜爽夜夜摸| 无遮挡黄片免费观看| 男女边吃奶边做爰视频| 男人的好看免费观看在线视频| 国产久久久一区二区三区| 久久久久久久久久成人| 91久久精品电影网| 亚洲av美国av| videossex国产| 亚洲av电影不卡..在线观看| 久久久久免费精品人妻一区二区| 中文字幕av成人在线电影| 国产精品,欧美在线| 六月丁香七月| 晚上一个人看的免费电影| 精品久久久久久久久久久久久| 国产精品乱码一区二三区的特点| av中文乱码字幕在线| 少妇的逼好多水| 日本一二三区视频观看| 久久精品国产亚洲av天美| 精品少妇黑人巨大在线播放 | 国产久久久一区二区三区| 观看免费一级毛片| 国产伦一二天堂av在线观看| 22中文网久久字幕| 久久中文看片网| videossex国产| 色av中文字幕| 国产激情偷乱视频一区二区| 身体一侧抽搐| 婷婷六月久久综合丁香| 在线观看美女被高潮喷水网站| 熟女人妻精品中文字幕| 老师上课跳d突然被开到最大视频| 亚洲色图av天堂| 99国产极品粉嫩在线观看| 国产一区二区激情短视频| 精品欧美国产一区二区三| 国产高清三级在线| 国产精品久久久久久久电影| 丰满人妻一区二区三区视频av| 久久韩国三级中文字幕| 99久国产av精品| 免费看光身美女| 国产三级在线视频| 成年女人永久免费观看视频| 国产亚洲av嫩草精品影院| 特大巨黑吊av在线直播| 91在线观看av| 天天躁夜夜躁狠狠久久av| 18禁在线播放成人免费| 国产69精品久久久久777片| 伦理电影大哥的女人| 国产伦精品一区二区三区视频9| 18禁在线无遮挡免费观看视频 | 99久久久亚洲精品蜜臀av| 免费不卡的大黄色大毛片视频在线观看 | 自拍偷自拍亚洲精品老妇| 午夜亚洲福利在线播放| 国产精品伦人一区二区| 国产精品99久久久久久久久| 欧美又色又爽又黄视频| 黄色视频,在线免费观看| 麻豆乱淫一区二区| 亚洲人成网站在线观看播放| 蜜臀久久99精品久久宅男| 成年女人永久免费观看视频| 久久精品国产亚洲网站| 波多野结衣巨乳人妻| 男女啪啪激烈高潮av片| 国产精品日韩av在线免费观看| 一进一出抽搐动态| 国产69精品久久久久777片| 亚洲不卡免费看| 久久鲁丝午夜福利片| 成人三级黄色视频| 亚洲欧美日韩无卡精品| 国产精品三级大全| 免费看美女性在线毛片视频| 成人av在线播放网站| 午夜精品一区二区三区免费看| 1000部很黄的大片| 有码 亚洲区| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 99热网站在线观看| 日韩亚洲欧美综合| 青春草视频在线免费观看| 中文亚洲av片在线观看爽| 99久久精品热视频| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 欧美一区二区精品小视频在线| 久久久久久大精品| 中出人妻视频一区二区| 看十八女毛片水多多多| 免费电影在线观看免费观看| 免费人成在线观看视频色| 日本黄色视频三级网站网址| 精品一区二区免费观看| 男人和女人高潮做爰伦理| 看黄色毛片网站| 亚洲精品在线观看二区| 欧美性猛交╳xxx乱大交人| 国产69精品久久久久777片| 亚洲最大成人中文| 成人国产麻豆网| 插阴视频在线观看视频| av黄色大香蕉| 亚洲第一区二区三区不卡| 少妇的逼好多水| 国产中年淑女户外野战色| 中国美白少妇内射xxxbb| 国产精品无大码| 欧美激情在线99| 亚洲久久久久久中文字幕| 国产亚洲精品久久久久久毛片| 国产男人的电影天堂91| 男插女下体视频免费在线播放| 可以在线观看毛片的网站| 日韩av不卡免费在线播放| 狂野欧美激情性xxxx在线观看| 麻豆av噜噜一区二区三区| 内射极品少妇av片p| 最新中文字幕久久久久| 国产成人影院久久av| 国产又黄又爽又无遮挡在线| 你懂的网址亚洲精品在线观看 | 又黄又爽又刺激的免费视频.| 九九在线视频观看精品| 丰满的人妻完整版| 在线观看午夜福利视频| 级片在线观看| 久久精品91蜜桃| 国产精品久久久久久精品电影| videossex国产| 午夜精品国产一区二区电影 | 精品久久久久久久末码| 99热这里只有精品一区| 不卡一级毛片| 丰满乱子伦码专区| 久久久久性生活片| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 夜夜爽天天搞| 好男人在线观看高清免费视频| 久久人妻av系列| 亚洲国产高清在线一区二区三| 日韩三级伦理在线观看| 午夜爱爱视频在线播放| 国产成人一区二区在线| 午夜精品在线福利| 午夜日韩欧美国产| 天天躁夜夜躁狠狠久久av| 久久人人爽人人爽人人片va| 国产亚洲av嫩草精品影院| 欧美潮喷喷水| 最好的美女福利视频网| 亚洲欧美精品自产自拍| 丰满人妻一区二区三区视频av| 国产亚洲精品久久久com| 亚洲人成网站在线观看播放| 色在线成人网| 嫩草影院入口| 老师上课跳d突然被开到最大视频| 老熟妇仑乱视频hdxx| 亚洲欧美日韩卡通动漫| 十八禁国产超污无遮挡网站| ponron亚洲| 国产私拍福利视频在线观看| 久久久成人免费电影| 简卡轻食公司| 欧美潮喷喷水| 精华霜和精华液先用哪个| 俺也久久电影网| 人人妻人人看人人澡| 日韩精品有码人妻一区| 性色avwww在线观看| 国产精品野战在线观看| 久久人人精品亚洲av| 校园人妻丝袜中文字幕| 亚洲第一电影网av| 日本成人三级电影网站| 国产av不卡久久| 极品教师在线视频| 亚洲中文字幕一区二区三区有码在线看| 激情 狠狠 欧美| 亚洲图色成人| 国产免费一级a男人的天堂| 听说在线观看完整版免费高清| 国产视频内射| 两性午夜刺激爽爽歪歪视频在线观看| 欧美bdsm另类| 日韩高清综合在线| 麻豆一二三区av精品| 日本免费一区二区三区高清不卡| 天美传媒精品一区二区| 特大巨黑吊av在线直播| 中文亚洲av片在线观看爽| 亚洲精品国产成人久久av| 在线看三级毛片| 伦理电影大哥的女人| 亚洲熟妇熟女久久| 三级毛片av免费| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品成人久久久久久| 午夜老司机福利剧场| 99热这里只有精品一区| 免费观看在线日韩| 免费不卡的大黄色大毛片视频在线观看 | 欧美激情在线99| 最新在线观看一区二区三区| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 欧美xxxx性猛交bbbb| 欧美潮喷喷水| 欧美bdsm另类| 一本久久中文字幕| 赤兔流量卡办理| av天堂中文字幕网| 成人特级av手机在线观看| 男人狂女人下面高潮的视频| 欧美3d第一页| 欧美一区二区精品小视频在线| 亚洲精品影视一区二区三区av| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 给我免费播放毛片高清在线观看| 亚洲第一电影网av| 12—13女人毛片做爰片一| 国产黄色小视频在线观看| 可以在线观看毛片的网站| 日韩欧美一区二区三区在线观看| 午夜爱爱视频在线播放| 国产视频一区二区在线看| 一夜夜www| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 亚洲婷婷狠狠爱综合网| 国产精品电影一区二区三区| 成人三级黄色视频| 午夜福利在线观看吧| 国产精品不卡视频一区二区| 国产精品免费一区二区三区在线| 男插女下体视频免费在线播放| 搞女人的毛片| 丰满人妻一区二区三区视频av| 精品午夜福利视频在线观看一区| 成人高潮视频无遮挡免费网站| 亚洲熟妇熟女久久| 国产精品伦人一区二区| 亚洲国产高清在线一区二区三| 欧美性猛交黑人性爽| 国产精品久久电影中文字幕| 身体一侧抽搐| 天堂动漫精品| 蜜桃亚洲精品一区二区三区| 熟女电影av网| 最近视频中文字幕2019在线8| 色av中文字幕| 男人舔奶头视频| 观看美女的网站| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 搡老妇女老女人老熟妇| 午夜精品国产一区二区电影 | 成年女人看的毛片在线观看| 日韩成人av中文字幕在线观看 | 精品久久久久久成人av| 国产色爽女视频免费观看| 嫩草影视91久久| 嫩草影院新地址| 国产高清三级在线| 男人舔奶头视频| 久久久色成人| 两性午夜刺激爽爽歪歪视频在线观看| www.色视频.com| 男插女下体视频免费在线播放| 午夜福利成人在线免费观看| 在线观看av片永久免费下载| 乱人视频在线观看| 国产乱人偷精品视频| 日本熟妇午夜| 久久久久久久午夜电影| 国产成人影院久久av| 精品国内亚洲2022精品成人| 一级黄色大片毛片| 国产精品乱码一区二三区的特点| 久久99热6这里只有精品| 女同久久另类99精品国产91| 搡老妇女老女人老熟妇| 黄色欧美视频在线观看| 日韩欧美一区二区三区在线观看| 人妻少妇偷人精品九色| 亚洲人成网站在线播放欧美日韩| 国产一区二区在线av高清观看| 国产av一区在线观看免费| 身体一侧抽搐| 岛国在线免费视频观看| 午夜日韩欧美国产| 亚洲乱码一区二区免费版| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 免费看美女性在线毛片视频| 国产又黄又爽又无遮挡在线| 黄片wwwwww| 亚洲人与动物交配视频| 亚洲人成网站在线播放欧美日韩| av在线天堂中文字幕| 一级av片app| 在线观看av片永久免费下载| 超碰av人人做人人爽久久| 麻豆一二三区av精品| 国产不卡一卡二| 色视频www国产| 国语自产精品视频在线第100页| 在线观看一区二区三区| 欧美成人精品欧美一级黄| 日本黄大片高清| 最近手机中文字幕大全| 国产成人一区二区在线| 听说在线观看完整版免费高清| 乱码一卡2卡4卡精品| 我的女老师完整版在线观看| 美女 人体艺术 gogo| 国产女主播在线喷水免费视频网站 | 搡老熟女国产l中国老女人| 舔av片在线| 精品一区二区三区av网在线观看| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 老熟妇仑乱视频hdxx| 人妻少妇偷人精品九色| 国产高清视频在线播放一区| 一级毛片我不卡| 欧美日韩一区二区视频在线观看视频在线 | 看十八女毛片水多多多| 综合色丁香网| 精品久久久久久久久久久久久| 91av网一区二区| 国产亚洲av嫩草精品影院| 九九在线视频观看精品| 99热只有精品国产| 在线播放国产精品三级| 此物有八面人人有两片| 久久精品国产亚洲网站| 你懂的网址亚洲精品在线观看 | av黄色大香蕉| 免费av毛片视频| 99在线视频只有这里精品首页| 99久久久亚洲精品蜜臀av| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 狠狠狠狠99中文字幕| 亚洲久久久久久中文字幕| 一区福利在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品亚洲一区二区| 久99久视频精品免费| 夜夜看夜夜爽夜夜摸| 一夜夜www| 国产高清三级在线| 日韩三级伦理在线观看| 欧美极品一区二区三区四区| 深夜a级毛片| 一边摸一边抽搐一进一小说| 精品福利观看| 一级毛片电影观看 | 日本免费一区二区三区高清不卡| 日韩av在线大香蕉| 少妇高潮的动态图| 免费观看人在逋| 亚洲人成网站在线观看播放| 国产高清视频在线观看网站| 成人综合一区亚洲| 97人妻精品一区二区三区麻豆| 女人十人毛片免费观看3o分钟| 精品久久久久久久久久免费视频| 亚洲国产欧洲综合997久久,| 国产欧美日韩精品一区二区| 亚州av有码| 美女高潮的动态| av在线亚洲专区| 免费看美女性在线毛片视频| 久久久久国产网址| 国产精品三级大全| 丰满少妇做爰视频| 日本黄色日本黄色录像| 伦理电影免费视频| 熟女av电影| 日韩av免费高清视频| 久热这里只有精品99| 国产熟女午夜一区二区三区 | 男人和女人高潮做爰伦理| 久久99一区二区三区| 亚洲婷婷狠狠爱综合网| 久久久久国产网址| av一本久久久久| 人人妻人人爽人人添夜夜欢视频 | 在线观看三级黄色| 在线观看人妻少妇| 大片电影免费在线观看免费| 亚洲精品乱码久久久v下载方式| 视频区图区小说| 乱码一卡2卡4卡精品| 久久久久精品久久久久真实原创| 乱人伦中国视频| 黄片无遮挡物在线观看| 婷婷色综合www| 久久99一区二区三区| 免费观看无遮挡的男女| 亚洲成人av在线免费| 国产高清国产精品国产三级| 久久av网站| 黑人猛操日本美女一级片| 国产伦理片在线播放av一区| 色视频在线一区二区三区| h视频一区二区三区| 伦理电影免费视频| 嫩草影院新地址| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 免费看不卡的av| 亚洲av福利一区| 日韩一区二区视频免费看| 国产视频首页在线观看| 国产免费又黄又爽又色| 国产精品三级大全| 老女人水多毛片| 中文精品一卡2卡3卡4更新| 99久久人妻综合| 爱豆传媒免费全集在线观看| 欧美丝袜亚洲另类| 欧美日韩视频精品一区| 日日啪夜夜爽| 亚洲美女黄色视频免费看| 一级av片app| 深夜a级毛片| 香蕉精品网在线| 久久热精品热| 最近的中文字幕免费完整| 两个人的视频大全免费| 国产中年淑女户外野战色| 国产精品福利在线免费观看| 久久精品久久久久久久性| 人人妻人人澡人人看| 日本91视频免费播放| 一级毛片 在线播放| h视频一区二区三区| 国产国拍精品亚洲av在线观看| 国产毛片在线视频|