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

    Prediction of clinically actionable genetic alterations from colorectal cancer histopathology images using deep learning

    2021-01-13 09:34:14HyunJongJangAhwonLeeKangInHyeSongSungHakLee
    World Journal of Gastroenterology 2020年40期

    Hyun-Jong Jang, Ahwon Lee, J Kang, In Hye Song, Sung Hak Lee

    Abstract

    Key Words: Colorectal cancer; Mutation; Deep learning; Computational pathology; Computer-aided diagnosis; Digital pathology

    INTRODUCTION

    Identifying genetic mutations in cancer patients has been increasingly important because mutational status can be very informative to determine the optimal therapeutic strategy[1]. However, molecular analysis is not performed routinely in every cancer patient, since it is not time and cost effective[2]. Thus, cost-effective alternatives for current molecular tests can be helpful in making appropriate treatment decisions. It has long been recognized that the histologic phenotypes reflect the genetic alterations in cancer tissues[3]. Since hematoxylin and eosin (H&E)-stained tissue slides are produced for almost every cancer patient, mutation prediction from the tissue slides can be a time- and cost-effective alternative method for individualized treatment. Thus, researchers attempted to examine the genotype–phenotype relationship in the H&E-stained tissue slides, and some gross tissue patterns related to specific molecular aberrations have been reported[4-9]. However, it remains largely unknown how specific molecular abnormalities are related to the specific histomorphologic findings, as it is not easy to capture the subtle features underlying the specific molecular alterations with the naked eye. To overcome the limitation of visual inspection of tissue structures by pathologists, various image analysis techniques have been applied for many decades to detect the subvisual characteristics of tissue patterns, not discernible to the unaided eyes[1]. Particularly, deep learning has been successfully applied to perform tasks considered too challenging for conventional image analysis techniques because it learns discriminative features directly from the large training dataset for any given task[10]. Therefore, deep learning is increasingly applied for tissue analysis tasks[11]. With the approval to use the digitized whole-slide images (WSIs) for diagnostic purposes, the digitization of tissue slides has been explosively increasing, providing huge digitized tissue data[12]. Combining the routine digitization of tissue slides with deep learning, the computer-aided analysis of WSIs could be adopted to support the evaluation of molecular alterations in H&E-stained cancer tissues in the near future. Although deep learning-based tissue analysis is still in its early phase, few promising results have been published. For example, a recent study reported that deep learning-based molecular cancer subtyping can be performed directly from the standard H&E sections obtained from patients with colorectal cancers (CRCs)[13]. Microsatellite instability can also be predicted from the tissue slides[14]. Furthermore, positive results for the mutation prediction of specific genes from histopathologic images have been reported in patients with various cancer types[3,15-17].

    Motivated by these recent studies, we tried to predict the frequently occurring and clinically meaningful mutations from the H&E-stained CRC tissue WSIs with deep learning-based classifiers. Based on the frequency of mutation and prognostic values of the genes, we choseAPC,KRAS,PIK3CA,SMAD4, andTP53genes for the current study. The area under the curves (AUCs) for the receiver operating characteristic (ROC) curves ranged from 0.645 to 0.809 for The Cancer Genome Atlas (TCGA) datasets, showing the potential for deep learning-based mutation prediction in the CRC tissue slides. By combining two different datasets for training, the prediction performance can be enhanced with the expansion of datasets.

    MATERIALS AND METHODS

    Tests with TCGA WSI dataset

    TCGA program offers the opportunity to reveal the genotype-phenotype relationship because it provides extensive archives of digital pathology slides with multi-omics test results[18]. Both frozen section tissue slides and formalin-fixed paraffin-embedded (FFPE) diagnostic slides were provided by the program. The WSIs from the TCGACOAD (colon cancer) and TCGA-READ (rectal cancer) projects were combined in this study because colonic and rectal adenocarcinoma share similar molecular and histological features[18]. After removing the WSIs with poor quality, 629 patients were included in the present study. We chose to include the genetic alteration including frame shift insertion and deletion, missense mutations, and nonsense mutation. ForAPC,KRAS,PIK3CA,SMAD4, andTP53genes, 436, 249, 133, 74, and 340 patients were confirmed to have the mutations, respectively. Deep learning did not perform optimally when there was a huge imbalance between classes[19]. In a previous study, we failed to obtain the balanced performance in tissue classification tasks unless the dataset itself was forced to have similar numbers between the classes[20]. Thus, we limited the difference in patient numbers between the mutation group and wild-type group by less than 1.4 fold through a random sampling. To match this limitation, we selected 263 patients withAPCmutation as there were only 188 patients with theAPCwild-type gene in the cohorts. The final patient IDs with their respective mutations are listed in Supplementary Table 1.

    Various artifacts including air bubbles, compression artifacts, out-of-focus blur, pen markings, tissue folding, and white background are unavoidable in the WSIs. To make the prediction process fully automated, these artifacts should be automatically removed. Because it is impractical to analyze a WSI as a whole, small image patches are often sliced from a WSI and used for the analysis. Thus, we built a deep learningbased tissue/non-tissue classifier for 360 × 360 pixel image patches at 20 × magnification to remove all of these artifacts at once (Figure 1A). The classifier was a simple convolutional neural network (CNN) with 12 (5 × 5), 24 (5 × 5), and 24 (5 × 5) convolutional filters, each followed by a (2 × 2) max pooling layer. The tissue/nontissue classifier could filter out more than 99.9% of improper patches. Next, tumor tissues should be delineated to predict the mutational status of cancer cells. Because of the freezing process for frozen tissue preparation, the frozen and FFPE tissue WSIs can differ in their morphologic features. Thus, we built separate normal/tumor classifiers for the frozen and FFPE WSIs based on the 360 × 360 pixel tissue image patches using the Inception-v3 model, a widely used CNN architecture. To train the wildtype/mutation classifiers for each gene, frozen and FFPE tissue patches with tumor probability higher than 0.9 by each tumor classifier were collected (Figure 1B). We arbitrarily chose the tumor probability as 0.9 because we decided to only include tissues with prominent tumor features. Although each slide may contain mixed regions of wild-type and mutated tissues considering the tumor heterogeneity, we assigned the same label for all tumor tissue patches in a WSI based on the mutational status of the patients. This labeling strategy was inevitable since we had no methods to delineate the wild-type and mutated regions before the classifiers could be built. The classifiers for the five genes were separately trained and validated with a patient-level ten-fold cross-validation scheme for frozen and FFPE WSIs. The slide-level mutation probability was calculated as the average of the probabilities of all the tumor patches in the WSI. For the training of the Inception-v3 models, we used a mini-batch size of 128, and the cross entropy loss function was adopted as a loss function. Deep neural networks were implemented using the TensorFlow deep learning library (http://tensorflow.org). To minimize overfitting, data augmentation techniques, including random rotations by 90°, random horizontal/vertical flipping, and random perturbation of the contrast and brightness, were applied to the tissue patches during training. In addition, 10% of the training slides were used as a validation dataset for the early stopping of the training. At least five separate classifiers were trained for each gene and tissue modality, and the classifier with the best AUC on the test dataset was included in the results.

    Figure 1 Fully automated prediction of mutation with three consecutive classifiers. A: Proper tissue patches can be selected by the tissue/non-tissue classifier. The four insets in the middle panel demonstrated the tissue patches representing pen marking, blurry scanned area, background rich region, and tissue folding, clockwise from top left, all removed by the tissue/non-tissue classifiers. Then, the normal/tumor classifier delineates the tumor patches among the proper tissue patches; B: The wild-type/mutation classifiers are applied only for patches with tumor probability higher than 0.9. The patch-level probabilities of mutation are averaged to yield the slide-level probability.

    Tests on the external cohorts

    Patient cohort:A total of 142 patients with CRC who previously underwent surgical resection in Seoul St. Mary’s hospital between 2017 and 2019 were enrolled (SMH dataset). All cases were sporadic, without any familial history of CRCs. The clinicopathological parameters including age, sex, and tumor location were retrospectively reviewed from the medical records. The study was approved by the Institutional Review Board of the College of Medicine at the Catholic University of Korea, No. KC19SESI0787.

    Mutation prediction on SMH dataset:ForAPC,KRAS,PIK3CA,SMAD4, andTP53genes, 66, 75, 31, 23, and 98 patients were confirmed to have the mutations, respectively. The sequencing methods are described in Supplementary Methods. Because the SMH dataset was originally collected to extra-validate the model trained on the TCGA datasets, we did not adjust the patient numbers between the classes. The normal/tumor classifier for TCGA FFPE tissues was also used to discriminate the tumor tissue patches of SMH WSIs. The normal/tumor classification accuracy was reviewed by Lee SH and Song IH and was confirmed to be valid. Again, patches with tumor probability higher than 0.9 were collected for mutational status classification. Then, the SMH data were split into ten folds, and each training fold was mixed with TCGA training fold to build new classifiers trained on both datasets. The classification results of the new classifiers on TCGA or SMH datasets were compared with the TCGA-based classifiers to investigate the effects of the expanded training dataset.

    Statistical analysis

    The ROC curves and their AUCs for all classifiers were presented to demonstrate the performance of each classifier. We used a permutation test with 1000 iterations to compare the differences between the two paired or unpaired ROC curves when necessary[21]. APvalue of < 0.05 was considered significant.

    RESULTS

    This study aimed to investigate the feasibility of mutation prediction for the frequently occurring mutations in the CRC tissue WSIs. Since only tumor tissues would be meaningful for the prediction of the mutational status in the tissue slides, three different tissue patch classifiers were sequentially applied to discriminate between tissue/non-tissue, normal/tumor, and wild-type/mutation in order (Figure 1). Only proper tissue patches with high tumor probabilities were used to determine the mutational status (Figure 1B). Patient-level ten-fold cross validation was applied for both frozen and FFPE datasets to fully evaluate the properties of the TCGA CRC WSIs.

    From Figures 2 to 6, the classification results forAPC,KRAS,PIK3CA,SMAD4, andTP53genes are presented for both frozen (upper panels) and FFPE (lower panels) TCGA WSIs. In A and C of every figure, the representative binary heatmaps demonstrating the distribution of tissue patches classified as wild-type or mutation are presented. From left to right, WSIs with gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation are presented, which were determined by the probability threshold set to 0.5. The sensitivity and specificity of a classifier can be much improved by setting the threshold appropriately. However, we set the threshold to 0.5 in the figures for simplicity because every classifier for different folds had different optimal thresholds. To demonstrate the differences in the performance between folds, slide-level ROC curves for folds with the lowest and highest AUCs were presented (left and middle ROC curves in the figures). Finally, the overall performance was inferred based on the slidelevel ROC curves drawn for the concatenated results from all ten folds (right ROC curves). For theAPCgene (Figure 2), the AUCs per fold ranged from 0.648 to 0.819 for the frozen tissues and from 0.655 to 0.880 for the FFPE tissues. The concatenated AUCs were 0.771 and 0.742 for the frozen and FFPE tissues, respectively. For theKRASgene (Figure 3), the performance was much better for the frozen tissues than for the FFPE tissues with a per fold AUC for the frozen tissues of 0.675-0.937 and a concatenated AUC of 0.778. For the FFPE tissues, the concatenated AUC was only 0.645, while the per fold AUCs ranged from 0.594 to 0.736. With regard to thePIK3CAgene (Figure 4), the lowest and highest AUCs per fold were 0.669 and 0.775 for the frozen tissues and 0.597 and 0.857 for the FFPE tissues. The concatenated AUCs were 0.713 and 0.690, respectively. For theSMAD4gene (Figure 5), AUCs per fold ranged from 0.619 to 0.849 for the frozen tissues and from 0.587 to 0.926 for the FFPE tissues, while the concatenated AUCs were 0.693 and 0.763, respectively. With regard to theTP53gene (Figure 6), the lowest and highest AUCs per fold were 0.707 and 0.963 for the frozen tissues and 0.737 and 0.805 for the FFPE tissues. The concatenated AUCs were 0.809 and 0.783, respectively. Overall, the wild-type/mutation classifiers for theTP53gene yielded the highest AUCs for both frozen and FFPE tissues of the TCGA datasets. Between the ROC curves of the frozen and FFPE tissues, classifiers for the frozen tissues yielded better results for theAPCandKRASgenes (P< 0.05,P< 0.001,P= 0.068,P= 0.057, andP= 0.115 between the frozen and FFPE classifiers forAPC,KRAS,PIK3CA,SMAD4, andTP53genes, respectively, by Venkatraman’s permutation test for unpaired ROC curves).

    The generalizability of a deep learning model for the external dataset is an important issue to be validated. Thus, we collected our own CRC FFPE WSIs with information on genetic mutation. The normal/tumor classifier for the TCGA FFPE tissues was applied to collect the tissue patches with high tumor probabilities. Then, the mutation classifiers for each gene trained on the TCGA FFPE tissues were applied to the tumor patches. The slide-level ROC curves for the five genes are presented in Supplementary Figure 1. The AUCs were 0.654, 0.581, 0.570, 0.652, and 0.775 forAPC,KRAS,PIK3CA,SMAD4, andTP53genes, respectively. For theAPC,KRAS, andPIK3CAgenes, the performance of the TCGA-based mutation classifiers on the SMH dataset were worse than that on the TCGA dataset (P< 0.01,P< 0.05,P< 0.05,P= 0.107, andP= 0.263 forAPC,KRAS,PIK3CA,SMAD4, andTP53genes, respectively, by Venkatraman’s permutation test for unpaired ROC curves). These results indicated that the mutation classifiers did not have an excellent generalizability when they were trained only with the TCGA WSI datasets. It remains unclear whether the performance could be improved when more data are used for the training. Thus, we combined the TCGA and SMH datasets to train new sets of mutation classifiers. Patient-level tenfold cross validation schemes were also used for the mixed dataset. The performance of the SMH dataset showed an obvious improvement, since the SMH data were included in the training data in this setting. The AUCs forAPCandKRASgenes increased to 0.812 and 0.832 (Figure 7,P< 0.01 andP< 0.001 compared with the TCGA-trained classifiers by Venkatraman’s permutation test for paired ROC curves). Improved results were also obtained forPIK3CA,SMAD4, andTP53with AUCs of 0.769, 0.782, and 0.845, respectively (Figure 8,P< 0.05,P< 0.01, andP< 0.05 by Venkatraman’s permutation test for paired ROC curves). More importantly, the performance of the TCGA data was also generally improved by the classifiers trained on both datasets (Supplementary Figure 2). The AUCs were 0.766, 0.694, 0.708, 0.791, and 0.822 for theAPC,KRAS,PIK3CA,SMAD4, andTP53genes, respectively (P= 0.072,P< 0.01,P= 0.091,P= 0.074, andP< 0.05 compared with the TCGA-trained classifiers). These results indicated that the deep learning-based classifiers for mutation prediction in tissue slides can yield better performance when more data are collected from various sources.

    Figure 2 Classifiers to predict APC gene mutation for the Cancer Genome Atlas colorectal cancer tissue slides. A: Representative whole slide images (WSIs) of the frozen slides with APC gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right, obtained with the classifiers trained with the frozen tissues; C and D: Same as A and B, but the results were for the formalin-fixed paraffin-embedded WSIs. APC-M: APC mutated; APC-W: APC wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    Figure 3 Classifiers to predict KRAS gene mutation for the Cancer Genome Atlas colorectal cancer tissue slides. A: Representative whole slide images (WSIs) of the frozen slides with KRAS gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right, obtained with the classifiers trained with the frozen tissues; C and D: Same as A and B, but the results were for the formalin-fixed paraffin-embedded WSIs. KRAS-M: KRAS mutated; KRAS-W: KRAS wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    DISCUSSION

    Figure 4 Classifiers to predict PIK3CA gene mutation for the Cancer Genome Atlas colorectal cancer tissue slides. A: Representative whole slide images (WSIs) of the frozen slides with PIK3CA gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right, obtained with the classifiers trained with the frozen tissues; C and D: Same as A and B, but the results were for the formalin-fixed paraffin-embedded WSIs. PIK3CA-M: PIK3CA mutated; PIK3CA-W: PIK3CA wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    In the present study, we selected theAPC,KRAS,PIK3CA,SMAD4, andTP53genes because they were frequently occurring in both TCGA and SMH CRC datasets and had prognostic values.APCis an important tumor suppressor known to play a role in CRC development. DeactivatingAPCleads to the constitutive activation of the Wnt signaling pathway, which may contribute to tumor progression[22]. The frequency ofAPCmutations was 47% for the SMH dataset, which is a slightly higher mutational rate compared with that in previous studies (24.2%-44.8%). The RAS proto-oncogenes (HRAS,KRAS, andNRAS) play a pivotal role in numerous basic cellular functions, such as control of cell growth, differentiation, and apoptosis, and regulate key signaling cascades including phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) pathways[23,24]. Mutations in RAS family members are found in 20% of all human cancers, of whichKRASmutations account for 85%[25].KRASmutated in 30% to 50% of patients with CRCs[25]. In the SMH dataset, the frequency was 53%.KRASis a critical oncogene involved in the MAPK signaling pathway, andKRASmutations promote colorectal adenoma growth in the early phase of carcinogenesis[26]. The presence of activatingKRASandNRASmutations is a predictor of resistance to epidermal growth factor receptor (EGFR) inhibitors, such as cetuximab or panitumumab[27,28]. ThePIK3CAgene is responsible for coordinating various cellular processes, including proliferation, migration, and survival. ThePIK3CAmutation is associated with the activation of downstream PI3K/Akt signaling, which in turn deregulates other signaling pathways that contribute to oncogenic transformations[29]. ThePIK3CAmutation occurs in 10%-30% of patients with CRCs[30]. In the present study, the frequency of thePIK3CAmutation was observed to be 22%. Recent studies have shown thatPIK3CAmutations are associated with a worse clinical outcome and with a negative prediction for anti-EGFR targeted therapy[31].SMAD4is an essential intermediator in the TGFβ signaling pathway, exhibiting a pivotal role as a tumor suppressor gene in CRC[32].SMAD4mutations occur in 10%-20% of patients with CRC[32,33]. In the SMH dataset, the rate of theSMAD4mutation was 16%. Recent studies have demonstrated that somaticSMAD4mutations are more common in patients with advanced stages, and a decrease in the level ofSMAD4expression is associated with worse recurrence-free and overall survival in patients with CRC[32]. The tumor suppressor geneTP53regulates DNA repair mechanism and apoptosis. Loss ofTP53function is one of the major events in the development of CRC, which is thought to occur in the later stages of colon cancer progression[34]. TheTP53mutation rate in the SMH dataset was 69%, which is consistent with the frequencies reported in various studies (45%-84%)[35].

    Figure 5 Classifiers to predict SMAD4 gene mutation for the Cancer Genome Atlas colorectal cancer tissue slides. A: Representative whole slide images (WSIs) of the frozen slides with SMAD4 gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right, obtained with the classifiers trained with the frozen tissues; C and D: Same as A and B, but the results were for the formalin-fixed paraffin-embedded WSIs. SMAD4-M: SMAD4 mutated; SMAD4-W: SMAD4 wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    Figure 6 Classifiers to predict TP53 gene mutation for the Cancer Genome Atlas colorectal cancer tissue slides. A: Representative whole slide images (WSIs) of the frozen slides with TP53 gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right, obtained with the classifiers trained with the frozen tissues; C and D: Same as A and B, but the results were for the formalin-fixed paraffin-embedded WSIs. TP53-M: TP53 mutated; TP53-W: TP53 wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    Figure 7 Mutation prediction of APC and KRAS genes for the Seoul St. Mary Hospital colorectal cancer tissue slides by the classifiers trained with both The Cancer Genome Atlas and Seoul St. Mary Hospital data. A: Representative whole slide images of the slides with APC gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right; C and D: Same as A and B, but the results were for the KRAS gene. SMH: Seoul St. Mary Hospital; APC-M: APC mutated; APC-W: APC wild-type; KRAS-M: KRAS mutated; KRAS-W: KRAS wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    In general, theAPCmutation is thought to have no prognostic significance[36]. However, in a specific situation such as in a microsatellite stable proximal colon cancer, wild-typeAPChas been associated with poorer survival[37]. On the contrary,KRAS,PIK3CA,SMAD4, andTP53gene mutations were associated with poorer prognosis in CRCs[34,38-40]. Thus, information on the mutational status of these genes can be useful in making therapeutic decisions for CRC patients. On occasion, a specific gene mutation can be related to a specific visual characteristic in tissue histology. For example, thePIK3CAmutation often coincides with lymphovascular invasion, tumor budding, and a high number of poorly differentiated clusters in CRC tissues[39]. However, it is not always possible to discover the visually discernible features reflecting the mutation of a specific gene. Therefore, we adopted deep learning to predict the mutational status of the five genes because the discriminative features of the mutations can be automatically learned directly from the large training data of tissue images. To our knowledge, this is the first study to evaluate the mutation prediction capabilities of deep learning models for the frequently occurring mutations in the pathologic tissue slides of CRC patients.

    Figure 8 Mutation prediction of PIK3CA, SMAD4, and TP53 genes for the Seoul St. Mary Hospital colorectal cancer tissue slides by the classifiers trained with both The Cancer Genome Atlas and Seoul St. Mary Hospital data. A: Representative whole slide images of the slides with PIK3CA gene mutation correctly classified as mutation, with wild-type gene correctly classified as wild-type, with gene mutation falsely classified as wild-type, and with wild-type gene falsely classified as mutation, from left to right; B: Receiver operating characteristic curves for the fold with lowest area under the curve (AUC), for the fold with highest AUC, and for the concatenated results of all ten folds, from left to right; C and D: Same as A and B, but the results were for the SMAD4 gene; E and F: Same as A and B, but the results were for the TP53 gene. PIK3CA-M: PIK3CA mutated; PIK3CA-W: PIK3CA wild-type; SMAD4-M: SMAD4 mutated; SMAD4-W: SMAD4 wild-type; TP53-M: TP53 mutated; TP53-W: TP53 wild-type; AUC: Area under the curve; FFPE: Formalin-fixed paraffin-embedded.

    In all the mutation classifiers applied to the TCGA frozen and FFPE tissues, the slide-level discrimination capabilities were much better against chance performance (P< 0.001 for all five genes by permutation test). These results indicated that the Inception-v3 model learned valid features to discriminate the mutated tissue phenotypes of each gene. In the case ofAPCandKRASgenes, the classifiers for the frozen tissues yielded better results compared with the FFPE tissues, although the frozen sections generally showed poorer tissue quality than did the FFPE sections. It can be explained by the fact that the frozen sections provided the best representation of the tissue contents on which the genomic signatures were tested[18]. Since the FFPE sections can be taken far from the frozen tissue sections, the mutational status can be different between them, considering the heterogeneity of large tumors. When we validated the classifiers trained with the TCGA FFPE tissues on the SMH WSIs, the performance was generally poorer (Supplementary Figure 1). Deep learning operates well under a condition where both the training and test datasets come from the same distribution[41]. For the H&E-stained tissue slides, the quality may vary because they undergo multiple processes for preparation including formalin fixation, paraffin embedding, sectioning, and staining, which can be slightly different between institutes[42]. Furthermore, the ethnic difference between the TCGA and SMH datasets may also contribute to the difference in the performance. Although the difference can be negligible to human eye, deep learning can be very sensitive to the subtle difference in tissue conditions. Therefore, many researchers insisted on the necessity of using large multi-national and multi-institutional datasets to enhance the generalizability of the deep learning model[2,12]. Thus, we combined the two datasets to build new classifiers trained on both TCGA and SMH datasets. Naturally, the performance for the SMH data was greatly enhanced because the tissue features of the data were exposed to the classifiers in this setting. More importantly, the performance of the TCGA data was also enhanced by adding the WSIs from the SMH dataset for training. These results clearly demonstrated that multi-national and multi-institutional datasets can improve the performance of the mutation classifiers. However, it remains unclear how far the performance can be improved if much more data are supplied.

    When we scrutinized the binary heatmaps of falsely classified WSIs, we recognized that the wild-type and mutated patches were generally aggregated rather than dispersed. The patterns implied the possibility that the tumor tissues in a tissue slide may have different mutational statuses between different regions. Large tumors can be molecularly heterogeneous, and the tumor heterogeneity can contribute to the resistance to treatment[43]. Therefore, tumor heterogeneity has been an important issue for both researchers and clinicians. To elucidate the spatial heterogeneity of a tumor, molecular methods with high spatial specificity such as multi-region sequencing and single-cell sequencing can be applied to examine a tissue sample. However, a random sampling of tissues for these molecular tests would be very inefficient. If possible regions of molecular heterogeneity in a tissue slide could be identified before the tests, molecular testing can be more specific and efficient. Furthermore, there are possibilities of false negative molecular tests because of the imprecise delineation of target regions in a tissue block[12]. Therefore, it is very important to objectively discriminate the tumor regions for the molecular evaluation of the tumor tissues. Thus, both normal/tumor and wild-type/mutation classifiers can be used to delineate the appropriate target sites for various molecular tests in cancer tissues. For example, Supplementary Figure 3 presents the heatmaps for the mutational status of all five genes in a TCGA frozen tissue slide, demonstrating how different regions of a slide can have different mutational statuses. When an overlaid probability map of mutation was drawn, areas with low and high mutational statuses can be recognized. It may not be easy to obtain this kind of information without the help of deep learning. Hence, molecular tests with high spatial specificity can be targeted to specific regions depending on the purpose of the tests. Therefore, these classifiers can make the selection of lesional regions for relevant multi-omics testing fully automated in the near future[2].

    Limitations also exist for the deep learning-based tissue classifiers. One of the limitations is the sensitive nature of deep learning to minute differences in the datasets. Because of the sensitive nature, classifiers applied to very subtly different conditions should be separately built. For example, classifiers for the frozen and FFPE tissues should be separately trained for the same tasks. It requires additional data collection and training overload. In clinical practice, pathologists should take an additional step to determine the kind of classifiers that should be applied for a specific specimen. It is currently inevitable to separately build classifiers to support various real-world tasks in the pathology laboratories. Therefore, manual selection of appropriate classifiers for target tasks is a necessary step that can limit the fully automated adoption of deep learning-based classifiers in the pathology workflows.

    In the current study, we used the high-throughput cancer panel to identify mutations in CRC tissues of the SMH dataset. This panel test approach makes it possible to identify diverse clinically actionable mutations in a single assay. However, it is quite expensive to prepare the equipment necessary to perform the test and to save a large number of data generated. This study demonstrated that a deep learningbased method could be a useful and effective tool for the prediction of actionable mutations from CRC WSIs. However, the interpretation of decision made by the deep learning-based classifier is unclear because of the black box nature of deep learning and should be further studied. Besides this aspect, the advantages and disadvantages between the mutation panel test (molecular test) and deep learning method were described in Table 1.

    Despite the limitation, with the increasing digitization of tissue slides, various computer-assisted methods will be introduced for histopathologic interpretation and clinical care. In the present study, we demonstrated the potential of deep learningbased classifiers to predict mutations in the CRC WSIs. Although the classifiers in this study are not yet enough to be used for predicting the genetic mutations in the clinic, deep learning-based methods have the potential to learn features for discriminating the wild-type tissues from the mutated tissues, which are not easily discernible to the human eye. Thus, deep learning will be increasingly adopted to discover new tissuebased biomarkers, which provide fundamental information for personalized medicine. With the accumulation of large sets of WSI data, deep learning-based tissue analyses will play important roles in the better characterization of cancer patients and will be an essential part of digital pathology in the era of precision medicine.

    CONCLUSION

    In the present study, we demonstrated that theAPC,KRAS,PIK3CA,SMAD4andTP53mutation can be predicted from H&E pathology images using the deep learningbased classifiers. Furthermore, by combining the TCGA and our datasets for training, the prediction performance was enhanced. Therefore, with the accumulation of tissue image data for training, deep learning can be used to supplement current molecular testing methods in the near future.

    Table 1 The advantages and disadvantages between the mutation panel test and deep learning-based method

    ARTICLE HIGHLIGHTS

    日韩有码中文字幕| 麻豆成人av在线观看| av黄色大香蕉| 首页视频小说图片口味搜索| 免费观看人在逋| 久久精品国产清高在天天线| 精华霜和精华液先用哪个| 2021天堂中文幕一二区在线观| 成人av一区二区三区在线看| 久久精品91蜜桃| 国产乱人伦免费视频| 亚洲自偷自拍图片 自拍| 9191精品国产免费久久| 久久久精品大字幕| 亚洲第一电影网av| 欧洲精品卡2卡3卡4卡5卡区| 国产激情久久老熟女| 午夜精品久久久久久毛片777| 99精品久久久久人妻精品| 俺也久久电影网| 日本a在线网址| 日日干狠狠操夜夜爽| 国产精品久久久久久人妻精品电影| 一个人免费在线观看的高清视频| 五月伊人婷婷丁香| 午夜福利18| av中文乱码字幕在线| 成人av一区二区三区在线看| 久久久久久九九精品二区国产| 国产爱豆传媒在线观看| 91字幕亚洲| 久久精品人妻少妇| 日韩欧美 国产精品| 久久精品aⅴ一区二区三区四区| 国产精品av视频在线免费观看| 香蕉av资源在线| 美女黄网站色视频| 亚洲国产欧洲综合997久久,| 美女黄网站色视频| 亚洲国产欧洲综合997久久,| 亚洲精品粉嫩美女一区| 亚洲精品国产精品久久久不卡| 天天添夜夜摸| 制服丝袜大香蕉在线| 嫩草影院精品99| 国产精品精品国产色婷婷| 免费看美女性在线毛片视频| 亚洲欧美一区二区三区黑人| 国产精品av久久久久免费| or卡值多少钱| 老鸭窝网址在线观看| 久久性视频一级片| 在线观看免费午夜福利视频| 国产美女午夜福利| 成人性生交大片免费视频hd| 在线观看免费午夜福利视频| 午夜福利成人在线免费观看| 看片在线看免费视频| 好看av亚洲va欧美ⅴa在| 欧美日韩综合久久久久久 | 国产精品野战在线观看| 国内精品一区二区在线观看| av福利片在线观看| 国产av在哪里看| 久久99热这里只有精品18| 91麻豆精品激情在线观看国产| 日本免费a在线| 男人的好看免费观看在线视频| 老司机午夜十八禁免费视频| 校园春色视频在线观看| 久久热在线av| 91九色精品人成在线观看| 国产熟女xx| 日本黄色视频三级网站网址| 日韩成人在线观看一区二区三区| 18禁黄网站禁片免费观看直播| www.自偷自拍.com| 三级男女做爰猛烈吃奶摸视频| 女同久久另类99精品国产91| 1000部很黄的大片| 亚洲欧美日韩高清在线视频| 一本精品99久久精品77| 小蜜桃在线观看免费完整版高清| 国产午夜精品论理片| 老熟妇乱子伦视频在线观看| 一级黄色大片毛片| 精品人妻1区二区| 日本 av在线| 十八禁网站免费在线| 国产精品女同一区二区软件 | 99久久精品国产亚洲精品| 欧美国产日韩亚洲一区| 窝窝影院91人妻| 亚洲专区字幕在线| 欧美三级亚洲精品| 村上凉子中文字幕在线| 天堂√8在线中文| 久久天躁狠狠躁夜夜2o2o| x7x7x7水蜜桃| 亚洲人成网站高清观看| 亚洲成av人片在线播放无| 国产精品99久久99久久久不卡| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久亚洲av鲁大| 中文字幕人成人乱码亚洲影| 久久久水蜜桃国产精品网| 黄频高清免费视频| 嫩草影院入口| 久久久久性生活片| 丰满人妻熟妇乱又伦精品不卡| 91老司机精品| 国产精品亚洲美女久久久| 一二三四在线观看免费中文在| 中文亚洲av片在线观看爽| 男女午夜视频在线观看| 99精品欧美一区二区三区四区| 在线观看美女被高潮喷水网站 | 黑人巨大精品欧美一区二区mp4| av天堂在线播放| 丰满的人妻完整版| 国产黄色小视频在线观看| 18禁裸乳无遮挡免费网站照片| 91麻豆av在线| 看免费av毛片| 一级黄色大片毛片| 国产真人三级小视频在线观看| 亚洲精品美女久久av网站| 久久天躁狠狠躁夜夜2o2o| 亚洲国产中文字幕在线视频| 两个人看的免费小视频| 三级男女做爰猛烈吃奶摸视频| 91字幕亚洲| 草草在线视频免费看| 性色av乱码一区二区三区2| 精品一区二区三区视频在线 | 女警被强在线播放| 可以在线观看毛片的网站| 欧美高清成人免费视频www| 在线a可以看的网站| 久久久久久九九精品二区国产| 日韩精品青青久久久久久| 亚洲国产看品久久| 亚洲狠狠婷婷综合久久图片| 啪啪无遮挡十八禁网站| 色视频www国产| 欧美乱色亚洲激情| 两个人的视频大全免费| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av熟女| 一个人观看的视频www高清免费观看 | 欧美最黄视频在线播放免费| 曰老女人黄片| 日本黄大片高清| 久久人人精品亚洲av| 午夜a级毛片| 国产一区二区激情短视频| a级毛片在线看网站| 美女 人体艺术 gogo| 欧洲精品卡2卡3卡4卡5卡区| 搡老熟女国产l中国老女人| 女同久久另类99精品国产91| 日本与韩国留学比较| 久久精品综合一区二区三区| 色精品久久人妻99蜜桃| 久久久国产成人免费| 亚洲国产精品久久男人天堂| 国产熟女xx| 亚洲人成网站在线播放欧美日韩| 巨乳人妻的诱惑在线观看| 亚洲熟妇中文字幕五十中出| 久久久久久久精品吃奶| 欧美色视频一区免费| 蜜桃久久精品国产亚洲av| 中文资源天堂在线| 成人av一区二区三区在线看| 五月玫瑰六月丁香| 好男人在线观看高清免费视频| 国产美女午夜福利| 欧美中文综合在线视频| 看黄色毛片网站| 亚洲中文字幕日韩| 中文字幕人成人乱码亚洲影| 极品教师在线免费播放| 熟女人妻精品中文字幕| 天堂影院成人在线观看| 国产毛片a区久久久久| 欧美又色又爽又黄视频| 91麻豆精品激情在线观看国产| 亚洲七黄色美女视频| www.www免费av| 日本三级黄在线观看| av福利片在线观看| 听说在线观看完整版免费高清| svipshipincom国产片| 12—13女人毛片做爰片一| 亚洲七黄色美女视频| 啦啦啦观看免费观看视频高清| 午夜影院日韩av| 十八禁人妻一区二区| 亚洲av成人一区二区三| 亚洲av美国av| 亚洲国产中文字幕在线视频| 久久人人精品亚洲av| 美女cb高潮喷水在线观看 | 一个人观看的视频www高清免费观看 | 国产高潮美女av| 亚洲精品中文字幕一二三四区| 夜夜躁狠狠躁天天躁| 美女大奶头视频| 男人和女人高潮做爰伦理| 特级一级黄色大片| 97超视频在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲真实伦在线观看| 色精品久久人妻99蜜桃| 免费看光身美女| 一二三四社区在线视频社区8| 亚洲第一欧美日韩一区二区三区| www.自偷自拍.com| 精品欧美国产一区二区三| 搞女人的毛片| 久久久精品欧美日韩精品| 最新美女视频免费是黄的| 亚洲自拍偷在线| 成年人黄色毛片网站| 国产免费男女视频| 变态另类丝袜制服| 观看免费一级毛片| 国产伦精品一区二区三区视频9 | 精品久久久久久久久久久久久| 此物有八面人人有两片| 麻豆国产av国片精品| 激情在线观看视频在线高清| 噜噜噜噜噜久久久久久91| www.熟女人妻精品国产| 亚洲精品一区av在线观看| 俺也久久电影网| 国产伦在线观看视频一区| 深夜精品福利| 天堂影院成人在线观看| 最新中文字幕久久久久 | 国产午夜福利久久久久久| 国产伦人伦偷精品视频| 国产精品电影一区二区三区| 国产欧美日韩一区二区三| 国产69精品久久久久777片 | 日韩欧美在线二视频| 久久久久亚洲av毛片大全| 禁无遮挡网站| 麻豆成人av在线观看| 在线观看美女被高潮喷水网站 | 国产麻豆成人av免费视频| 日本成人三级电影网站| 亚洲aⅴ乱码一区二区在线播放| 中文字幕高清在线视频| 哪里可以看免费的av片| 日本一本二区三区精品| 香蕉av资源在线| 久久草成人影院| 一进一出抽搐gif免费好疼| 麻豆av在线久日| 久久精品亚洲精品国产色婷小说| or卡值多少钱| 99国产精品99久久久久| 国产成人福利小说| 日韩大尺度精品在线看网址| 午夜免费观看网址| 丰满人妻熟妇乱又伦精品不卡| 无遮挡黄片免费观看| 国产淫片久久久久久久久 | 日韩欧美国产在线观看| 婷婷六月久久综合丁香| 久久草成人影院| 嫩草影视91久久| 欧美乱妇无乱码| 成人国产综合亚洲| 欧美成人性av电影在线观看| 国产亚洲精品综合一区在线观看| 国产av不卡久久| 欧美国产日韩亚洲一区| av欧美777| 欧美色欧美亚洲另类二区| 999精品在线视频| 中国美女看黄片| 不卡av一区二区三区| 亚洲av日韩精品久久久久久密| 黑人欧美特级aaaaaa片| 九九久久精品国产亚洲av麻豆 | 此物有八面人人有两片| 国产美女午夜福利| 午夜福利高清视频| 99精品久久久久人妻精品| 九九在线视频观看精品| 国模一区二区三区四区视频 | 香蕉国产在线看| av福利片在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人中文| 少妇裸体淫交视频免费看高清| 成人欧美大片| 老汉色av国产亚洲站长工具| 久久午夜综合久久蜜桃| 免费电影在线观看免费观看| 中文在线观看免费www的网站| 男女下面进入的视频免费午夜| 国产视频内射| 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 久久国产精品影院| 又大又爽又粗| 国语自产精品视频在线第100页| 免费观看人在逋| 国产亚洲精品综合一区在线观看| 噜噜噜噜噜久久久久久91| 国产激情欧美一区二区| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 国内精品久久久久精免费| 1024香蕉在线观看| 国产高清videossex| 国产又色又爽无遮挡免费看| 欧美一级毛片孕妇| 美女免费视频网站| 精品99又大又爽又粗少妇毛片 | 久久久国产成人精品二区| 波多野结衣高清无吗| 国产淫片久久久久久久久 | 啦啦啦观看免费观看视频高清| 色吧在线观看| 久久热在线av| 国产精品免费一区二区三区在线| 日本三级黄在线观看| 99re在线观看精品视频| 日韩大尺度精品在线看网址| 黄色 视频免费看| 俄罗斯特黄特色一大片| 日韩欧美国产在线观看| 久久欧美精品欧美久久欧美| 久久天躁狠狠躁夜夜2o2o| 国产成人精品久久二区二区免费| 久久久成人免费电影| 美女 人体艺术 gogo| 在线观看66精品国产| 午夜福利免费观看在线| 久久久久精品国产欧美久久久| 成人av一区二区三区在线看| 性色av乱码一区二区三区2| 国产高清三级在线| 亚洲成人精品中文字幕电影| 成人无遮挡网站| 黑人操中国人逼视频| 在线免费观看的www视频| 色av中文字幕| 伊人久久大香线蕉亚洲五| 免费看十八禁软件| 女生性感内裤真人,穿戴方法视频| 男人和女人高潮做爰伦理| 亚洲自偷自拍图片 自拍| 国产欧美日韩精品亚洲av| 免费在线观看成人毛片| 又紧又爽又黄一区二区| 欧美性猛交╳xxx乱大交人| 最近最新中文字幕大全电影3| 人人妻人人澡欧美一区二区| 黄片大片在线免费观看| 熟妇人妻久久中文字幕3abv| 97超视频在线观看视频| 老司机午夜福利在线观看视频| 一区二区三区激情视频| 日本免费a在线| 日韩大尺度精品在线看网址| 国产一区二区在线观看日韩 | 久久亚洲精品不卡| bbb黄色大片| 国产精品99久久99久久久不卡| 亚洲狠狠婷婷综合久久图片| 久久热在线av| 在线观看一区二区三区| 99热6这里只有精品| 99久久久亚洲精品蜜臀av| 色精品久久人妻99蜜桃| 亚洲自偷自拍图片 自拍| 国产三级黄色录像| 国产成人精品久久二区二区免费| 91av网站免费观看| 天天一区二区日本电影三级| 国产精品久久久久久人妻精品电影| 国产精品一区二区三区四区免费观看 | cao死你这个sao货| 日本成人三级电影网站| 后天国语完整版免费观看| 亚洲成人久久性| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 国产97色在线日韩免费| 色噜噜av男人的天堂激情| 国产99白浆流出| 高清毛片免费观看视频网站| 免费在线观看影片大全网站| x7x7x7水蜜桃| 五月伊人婷婷丁香| 国产高清videossex| 99热只有精品国产| 日日夜夜操网爽| 琪琪午夜伦伦电影理论片6080| 日韩欧美在线乱码| 两个人视频免费观看高清| 一个人看的www免费观看视频| 中亚洲国语对白在线视频| 欧美日韩综合久久久久久 | 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久精品吃奶| 变态另类成人亚洲欧美熟女| 久久性视频一级片| 成人国产一区最新在线观看| 免费看日本二区| 午夜福利欧美成人| 国产高清videossex| 亚洲av免费在线观看| 男女那种视频在线观看| 国产精品日韩av在线免费观看| 看黄色毛片网站| 我要搜黄色片| 黄色日韩在线| 亚洲五月婷婷丁香| 欧美+亚洲+日韩+国产| 精品无人区乱码1区二区| 久久人人精品亚洲av| 欧美极品一区二区三区四区| 宅男免费午夜| 日日摸夜夜添夜夜添小说| 午夜福利免费观看在线| 国内精品一区二区在线观看| 日韩高清综合在线| 午夜视频精品福利| 亚洲一区高清亚洲精品| 国产免费av片在线观看野外av| 国产aⅴ精品一区二区三区波| 免费高清视频大片| 亚洲av片天天在线观看| 亚洲av熟女| 在线a可以看的网站| 真人一进一出gif抽搐免费| 999久久久国产精品视频| 亚洲成a人片在线一区二区| 国产一区在线观看成人免费| 高清在线国产一区| 国产精品久久久久久精品电影| 国产成人av激情在线播放| 99国产极品粉嫩在线观看| e午夜精品久久久久久久| 国产一区二区在线观看日韩 | 日日夜夜操网爽| 99国产综合亚洲精品| 国产毛片a区久久久久| 男插女下体视频免费在线播放| 久久久久久九九精品二区国产| 色视频www国产| 精品国产三级普通话版| 亚洲在线自拍视频| 99国产综合亚洲精品| 精品福利观看| 久久午夜亚洲精品久久| av中文乱码字幕在线| 精品99又大又爽又粗少妇毛片 | 熟妇人妻久久中文字幕3abv| 精品久久久久久久末码| 99久久久亚洲精品蜜臀av| 国产亚洲欧美98| 草草在线视频免费看| 怎么达到女性高潮| 久久中文字幕人妻熟女| 亚洲精品一卡2卡三卡4卡5卡| 精品欧美国产一区二区三| 国产av在哪里看| 亚洲精品乱码久久久v下载方式 | 精品一区二区三区视频在线 | 日本免费a在线| 在线a可以看的网站| 精品熟女少妇八av免费久了| 日韩欧美精品v在线| 高清在线国产一区| 视频区欧美日本亚洲| 亚洲国产精品合色在线| 真人一进一出gif抽搐免费| 女人高潮潮喷娇喘18禁视频| 精品午夜福利视频在线观看一区| 亚洲av美国av| 一本综合久久免费| 黄片大片在线免费观看| 亚洲精品一区av在线观看| 日本三级黄在线观看| 国产成人福利小说| 日韩精品中文字幕看吧| 精品一区二区三区视频在线 | cao死你这个sao货| 每晚都被弄得嗷嗷叫到高潮| 国产成人欧美在线观看| 1024香蕉在线观看| 在线观看舔阴道视频| 久久精品91无色码中文字幕| 国产三级中文精品| 国产激情久久老熟女| 亚洲熟妇熟女久久| 少妇熟女aⅴ在线视频| 青草久久国产| 亚洲自拍偷在线| 国产又黄又爽又无遮挡在线| 欧美性猛交╳xxx乱大交人| av天堂中文字幕网| 国产麻豆成人av免费视频| 欧美在线黄色| 久久久精品大字幕| 亚洲精品中文字幕一二三四区| 精品国产亚洲在线| 国产成人欧美在线观看| 五月玫瑰六月丁香| 天天一区二区日本电影三级| 真人一进一出gif抽搐免费| 欧美av亚洲av综合av国产av| 亚洲在线自拍视频| 亚洲精品中文字幕一二三四区| 国产一级毛片七仙女欲春2| 久久久国产成人精品二区| 日韩欧美三级三区| 亚洲自偷自拍图片 自拍| 久久九九热精品免费| 亚洲美女黄片视频| 色尼玛亚洲综合影院| 香蕉av资源在线| 亚洲人成网站高清观看| 天堂动漫精品| www.熟女人妻精品国产| 国产亚洲精品一区二区www| 丰满人妻一区二区三区视频av | 欧美成人性av电影在线观看| 色老头精品视频在线观看| 欧美激情在线99| 亚洲专区国产一区二区| a在线观看视频网站| 国内精品美女久久久久久| 男插女下体视频免费在线播放| 身体一侧抽搐| 可以在线观看的亚洲视频| 神马国产精品三级电影在线观看| 99久久精品国产亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av电影不卡..在线观看| 国产精品av视频在线免费观看| 级片在线观看| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合| 香蕉丝袜av| 手机成人av网站| 一本久久中文字幕| 国产精品综合久久久久久久免费| 亚洲人成伊人成综合网2020| 亚洲自拍偷在线| 特大巨黑吊av在线直播| 丰满的人妻完整版| 国产成人av激情在线播放| 久久久久久久午夜电影| 亚洲自偷自拍图片 自拍| 亚洲精华国产精华精| 欧美一区二区精品小视频在线| www国产在线视频色| 国产精品久久久人人做人人爽| 精品久久久久久久久久免费视频| 久久香蕉国产精品| 一本综合久久免费| 91久久精品国产一区二区成人 | 国产一区二区三区在线臀色熟女| 亚洲国产欧美网| 免费搜索国产男女视频| 老汉色∧v一级毛片| 91久久精品国产一区二区成人 | 夜夜躁狠狠躁天天躁| 亚洲精品粉嫩美女一区| 国产亚洲精品一区二区www| 夜夜爽天天搞| 三级国产精品欧美在线观看 | 人人妻,人人澡人人爽秒播| 久久久国产成人精品二区| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品av一区二区| 丰满人妻一区二区三区视频av | 日韩欧美在线二视频| 国产高清视频在线观看网站| 久久久久久久久中文| 成人午夜高清在线视频| 女同久久另类99精品国产91| 欧美日韩一级在线毛片| 国产成人系列免费观看| 亚洲第一电影网av| www.999成人在线观看| 国产欧美日韩精品一区二区| 亚洲九九香蕉| 亚洲欧美日韩东京热| 俺也久久电影网| 男人舔女人的私密视频| 亚洲一区高清亚洲精品| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添小说| 18禁美女被吸乳视频| 亚洲成人免费电影在线观看| 99热6这里只有精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲在线观看片| 国产精品女同一区二区软件 | 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 69av精品久久久久久| 麻豆成人av在线观看| 精品久久久久久久末码|