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

    Urdnet: A Cryo-EM Particle Automatic Picking Method

    2022-08-24 12:56:42JianquanOuyangYueZhangKunFangTianmingLiuandXiangyuPan
    Computers Materials&Continua 2022年7期

    Jianquan Ouyang, Yue Zhang, Kun Fang,,*, Tianming Liuand Xiangyu Pan

    1School of Computer Science & School of Cyberspace Science, Xiangtan University, Xiangtan, 411105, China

    2Hunan Meteorological Information Center, Hunan Meteorological Bureau, Changsha, Hunan, 410118, China

    3Department of Computer Science, The University of Georgia, Athens, Georgia, USA

    Abstract: Cryo-Electron Microscopy (Cryo-EM) images are characterized by the low signal-to-noise ratio, low contrast, serious background noise, more impurities, less data, difficult data labeling, simpler image semantics, and relatively fixed structure, while U-Net obtains low resolution when downsampling rate information to complete object category recognition, obtains highresolution information during upsampling to complete precise segmentation and positioning, fills in the underlying information through skip connection to improve the accuracy of image segmentation, and has advantages in biological image processing like Cryo-EM image.This article proposes A U-Net based residual intensive neural network (Urdnet), which combines point-level and pixel-level tags, used to accurately and automatically locate particles from cryo-electron microscopy images, and solve the bottleneck that cryo-EM Single-particle biological macromolecule reconstruction requires tens of thousands of automatically picked particles.The 80S ribosome, HCN1 channel and TcdA1 toxin subunits, and other public protein datasets have been trained and tested on Urdnet.The experimental results show that Urdnet could reach the same excellent particle picking performances as the mainstream methods of RELION,DeepPicker, and acquire the 3Dstructure of picked particles with higher resolution.

    Keywords: Deep learning; convolutional neural network; particle picking;cryo-electron microscopy; single-particle reconstruction

    1 Introduction

    Cryo-EM has become an essential structural biology technology.It freezes the sample and keeps it in the microscope at a low temperature.Subsequently, the highly coherent electrons, used as a light source, illuminate from above and are scattered by the sample and the nearby ice layer.The scatter signal is then imaged and recorded using a detector and a lens system.Finally, signal processing is performed to obtain the structure of the sample, which is a valuable means of understanding the mechanism of biochemical reactions.Just as the catalytic sites of some proteins are known through structures,inhibitors can be designed to deactivate these proteins, and drugs can be screened with these proteins as targets.

    The other two leading technologies are X-ray crystallography and nuclear magnetic resonance(NMR), which can analyze the structure of a biomacromolecule to understand its function.For decades, X-ray crystallography has been the dominant technology for obtaining high-resolution structures of biomacromolecule.In recent years, with the improvement of cryo-EM [1,2] and the latest technological advances in sample preparation, calculation, and instrumentation, the structural analysis of protein complexes has undergone a considerable breakthrough [3] that the resolution of the large protein 3D structure increased to 3 ?[4,5].

    The cryo-EM image contains 2D projections of the particles at different angles.Due to the similarity of the sample density of the protein to its surrounding solution, and the limited electron dose used in data collection, the cryo-EM image has a low SNR, low contrast, uneven background intensity, and irregular internal grain texture.Besides, bad particles that do not meet the requirements for 3D reconstruction emerged, which is caused by image impurities such as frozen liquids, carbon film, and stacked or dissociated particles.To obtain a high-resolution 3D structure, tens of thousands of projection images are typically required.Thus, particle picking is the first step toward the 3D reconstruction of macromolecules, including the particle recognition in micrograph and localization of the particle region, rather than the noise region and impurity region.

    Existing methods for particle picking often use low-resolution 2D particle templates for template matching, where templates are generated fromthe clustering projections of manually selected particles.The basic idea of it is that the micrograph region has high cross-correlations with popular particle templates [6].RELION [7] is a cryo-EM structure determination software that uses a template matching framework to select particles [8].In the software, the user manually picks approximately 1000 particles from a fewmicrographs.These picked particle images are two-dimensionally classified to produce a small number of template images for automatically selecting particles from all micrographs.One problem with the approach is that it is susceptible to reference dependency biases resulting in high error detection rates.

    In addition to RELION, many tools include automatic or semi-automated particle picking steps, such as PICKER [9], EMAN2 [10], XMIPP [11], cryoSPARC [12], most of which are based on traditional computational vision algorithms, such as edge detection, feature recognition, and template matching aforementioned.These methods are not entirely suitable for processing cryo-EM images with poor contrast and low SNR, for they do not take full advantage of the inherent and unique particle features.Moreover, their performance degrades significantly as the quality of the microscopic image decreases.Consequently, it is significant to develop an efficient, fully-automated, template-free method for particle picking.

    In cryo-EM image analysis, many machine learning algorithms are widely used from support vector machines to convolutional neural networks.In the past fewyears, deep learning has proliferated.It can be superior to lots of traditional algorithms in computer vision by generating hierarchical features from big data analysis with deep neural networks [13,14].Besides, some deep learning applications are robust to low SNR images [15].As the cryo-EM image set continues to grow, and the SNR of the microscopic image is still low, deep learning seems to be very suitable for processing cryo-EM data.Inspired by some successful applications, we seek to apply deep learning methods to cryo-EM particle picking.

    In this paper, an improved U-Net-based residual dense neural convolutional network, Urdnet,is proposed to select the cryo-EM particles accurately and automatically.The method introduces a learning method combining point-level and pixel-level labels, which significantly saves the time of manual labeling and improves the efficiency and accuracy of particle picking.

    2 Related Work

    DeepPicker [16], a fully automatic particle selection method based on deep learning, applied the convolutional neural network to the particle picking task for the first time.Before the advent of DeepPicker, semi-automatic solutions such as RELION and EMAN2 were used to pick particles.DeepPicker converts particle picking into an image classification problem; it crops the microscopic images through a sliding window and classifies the sub-images into particles or backgrounds.Considering the lack of training data, DeepPicker uses a novel cross-molecular training strategy to train networks.However,the disadvantage is the long processing time, and the average particle picking speed is 1.5min/mrc (mrc is a format for cryo-EM images).

    DeepEM [17] uses an eight-layer convolutional neural network to recognize single particles of noisy cryo-EM micrographs and to achieve automated particle picking, selection, and verification in an integrated manner.DeepEM increases particle images through image rotation to augment the training set, but still requires manually selecting thousands of particles to train the data.

    FastParticlePicker [18] is based on the object detection algorithm, Fast R-CNN, which includes a“region-of-interest proposal”network and a classification network.However, the FastParticlePicker crops the microscopic image with a sliding window instead of selecting the area of interest in the microscopic image.Therefore, its performance mainly depends on the classification network.The classification network has three types of objects, particles, ice, and background, which reduces the false-positive results caused by ice.

    FCRN [19] is also an automatic particle picking method based on deep learning.It proposes a Fully Convolution Regression Network (FCRN) mapping particle images to continuous distance maps to recognize particles from different data sets.Experimental results on EMPIAR data show that FCRN achieves better particle picking performance than Faster-RCNN and RELION.

    crYOLO [20] is based on the YOLO9000 algorithm to achieve particle detection.For small particles, crYOLO achieves higher precision than the original YOLO network.It achieves a processing speed of up to 6 mrcs/s and can be extended to other biomacromolecules outside the training set.Despite the modifications to the original YOLO model, the experiment did not mention how crYOLO detects particles of different sizes and aspect ratios.

    Topaz [21] is an efficient and accurate particle selection method based on deep learning.Unlike other methods, it uses a positive unlabeled (PU) learning framework to reuse the remaining unlabeled particles to train neural networks of a small number of labeled particles.Even for the challenging dataset of rod-like particles with low SNR, the experimental results are superior to the general PU learning method.

    PIXER [22] is an automatic particle picking method based on image segmentation using deep neural networks.To adapt to the low SNR, it uses a segmentation network to convert the microscopic image into a probability density map to detect the particles.A grid-based local maximum method is proposed to locate the particles from the probability density map.Compared to mainstream methods,PIXER can achieve as good results as the semi-automatic methods RELION and DeepEM.

    Compared with natural images, cryo-EM images have the following characteristics: (1) The semantics of the images are relatively simple, and the structure is fixed.Therefore, high-level semantic information and low-level features are both important.U-Net’s skip connection and U-shaped structure are quite suitable to deal with our scenario.(2) Few sources of data and the difficulty of image labeling lead us to design a model that should not be too large.If parameters of the model are excessive, it will easily lead to overfitting.(3) Compared with natural images, cryo-EM images have a higher signal-to-noise ratio (SNR), higher complexity, a larger grayscale, and more unclear boundaries.U-Net obtains low-resolution information during down-sampling to recognize different objects and obtains high-resolution information during up-sampling to achieve segmentation and positioning of objects.In them eantime, it fills in the underlying information through skip-connections to improve segmentation accuracy.Simple CNN cannot achieve the global feature fusion of U-Net,thatis, it cannot connect shallow features and deep features.Therefore, direct training of CNNis often suitable for a large number of datasets, while U-Net with a U-shaped network uses skip-connection to reduce the time of feature extraction so that the entire network can better remember the overall picture information.In addition, the goal of this article is not only to use CNN feature extraction to detect cryo-EM particles but to obtain high-quality particle images for the three-dimensional reconstruction of biological macromolecules.In addition to identifying particles, it is necessary to accurately segment artifacts to avoid picking out defective low-mass particles.However, the efficiency of U-Net in biomedical image segmentation is well known [23].

    For the shortcomings of the above methods, Combine the advantages of U-Net, we developed an accurate and fully automatic particle picking approach, Urdnet, based on deep learning.Among it, the pertinent data preprocessing raised the quality of the raw Cryo-EM images.The image annotation method that combines with multiple labels significantly reduced time and effort for manually building a training dataset.Urdnet can thoroughly learn the features of particles with different particle sizes and aspect ratios via a small amount of training data.The connected component analysis is performed on the pixel map predicted by Urdnet to locate the candidate particles, which avoided selecting most of the false-positive particles.Moreover, experimental results show our method can apply to multiple types of micrographs from different cryo-EM detectors and has superior particle detection precision compared to the DeepPicker method [24-26].Compared with the semi-automatic picking ofRELION, our method resulted in a 3D single-particle structure with a higher resolution.More details of our proposed method will be introduced next.

    3 Proposed Method

    3.1 Data Preprocessing

    3.1.1 Intensity Adjustment

    Particle picking is significantly dependent on the intensity of the grayscale cryo-EM image.Due to the low-dose electron microscopy imaging mode, the randomly distributed particle and thin ice in high defocus areas are exposed to extremely the low-intensity beam, and low-contrast cryo-EM images are collected in that imaging mode.To alleviate this problem, we applied a contrast enhancement method to adjust the global image intensity and improve the SNR of cryo-EM images [27,28].The steps are as follows: (1) calculate two segmentation thresholds xland xhaccording to a percentage of pixel saturation of 1%, respectively indicating that the number of pixels less than xland more than xhboth accounts for 1%; (2) the gray value in [xl, xh] of the original image is linearly mapped to [0, 1], as shown in Eq.(1), where x is a pixel value of the original image, and xfis the new pixel value after the mapping.

    Fig.1 shows the raw cryo-EM image, the intensity-adjusted image, and their corresponding histograms.(c) has higher contrast than (a), and correspondingly, (d) has a wider range of pixel distribution than (b).

    Figure 1: (a) is a cryo-EM image of the TcdA1 toxin subunit; (b) is the histogram of (a); (c) is the intensity-adjusted image of (a); (d) is the histogram of (c)

    3.1.2 Noise Suppressing

    Cryo-EM images tend to be so noisy and blurred that proteins and solvents are visually similar since low-electron dose imaging only produces a small number of available signals of the taken images.Hence, we apply image restoration to improve image quality.For the Gaussian noise in the cryo-EM image, we choose the Wiener filter [29] to reduce the noise of the blur area and improve the sharpness of the defocused phase-plate cryo-EM image.Wiener filtering is a method of filtering the noise-mixed signal by minimizing the overall mean square error at the inverse filtering and noise smoothing.Its mathematical expression is as in Eq.(2) Where H(f) is the Fourier transform of h in the frequency domain f, SNR(f) is the SNR, and * represents the complex conjugate.

    Fig.2 shows the intensity adjusted cryo-EM image in Section 3.1.1 and its Wiener filtered image,where PSNR stands for Peak Signal-to-Noise Ratio, which is a widely used index for evaluating image quality and signal restoration.The index is calculated by the mean square error (MSE).The MSE is calculated by Eq.(3).I and K are images of size m * n, which represent respectively the raw cryo-EM image and the processed cryo-EM image.The peak signal-to-noise ratio is defined as Eq.(4), where MAXIrepresents the maximum pixel value of the image.

    Figure 2: (a) is a cryo-EM image of the TcdA1 toxin subunit (EMPIAR-10089); (b) is Wiener filtered image of (a)

    As shown in Fig.2, the PSNRs of the intensity-adjusted image aforementioned, and its noise suppressed image were calculated.Comparing the two PSNR values, 22.60 of Fig.2a is higher than 17.00of Fig.2b,which indicates that the Wiener filter has a direct noise suppression effect on the cryo--EM images.Compared with Figs.2a and 2b makes it easier to distinguish the particle boundaries.

    3.2 Combination of Point Annotation and Per-Pixel Annotation

    Various weak annotation forms have been explored in weakly supervised annotation methods for semantic segmentation of natural scene images, such as image-level labels (like classification labels)[30-32], extreme points [33], graffiti [34] and bounding boxes (like detection labels) [35].The point form is the most straightforward and simplest of all weak annotation forms [36].Since a cryo-EM image often contains a large number of particles, we proposed an annotation method of point labeling of particle center for cryo-EM particle picking to reduce the burden of manual annotation greatly.We divide objects into two classes: particles and artifacts.Artifacts are usually areas of cryo-EM impurities.We annotate the particles as points so that they are only associated with several pixels.The number of artifacts is much smaller than particles, and artifacts are particularly large in the cryo-EM images.To prevent the bad particles in the artifact from being selected, the artifacts are labeled at the pixel level like fully supervised segmentation labels, and the remaining unlabeled pixels are automatically classified as the background.Hence, we combine different forms of annotation to label particles and artifacts in a cryo-electron image.Fig.3b shows the label image of an 80S ribosome cryo-EM image.

    Figure 3: (a) is a cryo-EM image of the 80S ribosome; (b) is a label image of (a), the red region is the artifact’s label, and the green points are the particles’labels

    3.3 Convolutional Neural Network Design

    U-Net [37] has achieved excellent results in medical image segmentation.Its most significant advantage is that it combines shallow and in-depth features and can reach high segmentation precision with a small training dataset.In the particle picking task, the precision of artifact segmentation affects the quality of picked particles.Imprecise artifact segmentation can easily lead to selecting false-positive particles.To improve the precision of artifact segmentation, we embedded the residual dense block (RDB) [38] in the overall model of improved U-Net, as shown in Fig.4 when designing the Urdnet model.As shown in Fig.5, we add two convolution layers as the input port to the basic RDB structure for local shallow feature extraction.RDB mainly integrates the residual block (Eq.(5)) [39] and the dense block (Eq.(6)) [40].The residual block short circuits the output layer with the previous input layer by adding the feature vectors.The dense block concatenates each middle layer with all previous layers in the channel dimension to implement feature multiplexing.Combining the advantages of both, RDB performs not only local residual learning but also strengthens local feature fusion.Accordingly, it can effectively alleviate gradient vanishing, enhance feature propagation and reduce the number of parameters.RDB was proposed for processing image restoration, image super-resolution, and denoising tasks at first.Consequently, RDB is well suited for processing low-SNR and low-resolution cryo-EM images.To the best of our knowledge, our work is the first to combine U-Net with RDB to solve cryo-EM particle picking tasks.

    As shown in formula (5), the outputxlof the residual block is defined as the addition of f(xl-1) ,the output ofxl-1after multiple convolutions and activations andxl-1, the output of thel- 1 layer.

    As shown in formula (6), the lth layer receives the feature matrixx1,x2...xl-1from all previous layers as input,and[x1,x2...xl-1]is the connection feature vector of all previous layers,Hl(.)represents the amount of a single vector after convolution, activation, and other operations.

    Figure 4: Structure of the Urdnet.The hyperparameter f1:64 indicates the vector channels number of the input and output of the RDB is 64, f2, f3, and f4 are the same.Different colored arrows indicate different operations on the feature vector, and the numbers above the blue vectors indicate the number of channels

    Figure 5: RDB of Urdnet

    The network of Fig.4 consists of 49 convolution layers, four max-pooling layers, four up-sampling layers, and four feature concatenations.In the encoder, i.e., feature extraction part, the first convolution layer learns the shallow global features, then followed by four RDBs for deep feature fusion.Each RDB consists of 8 standard convolution layers and one 1*1 convolution layer.In an RDB, the number of filters in the first two convolution layers and the last 1*1 convolution layer is determined by the external hyperparameters, such as 128 for the second RDB, and the number of filters of the remaining six convolution layers is fixed at 64, and the dashed arrows in Fig.5 refer to the six standard convolution layers.The decoder combines the global context features with four times of upsampling and four feature concatenations, and finally restores the feature vector to the input vector size.We input a 512 × 512 2D cryo-EM image into Urdnet and apply three 1*1 convolution filters to output the feature map onto three channels, to obtain a softmax function (Eq.(7)) score map of 512×512×3 where 3 indicates the number of object classes.While testing, the Arg max function gets the class of the highest score at each coordinate of the score map [40] and finally outputs the predicted pixel map(background: black, grain: green, artifact: red).

    The softmax function is shown in formula (7),xlis the output feature map of thelth network layer, e is the natural constant,represents the output value on the ith channel, C is the number of network categories, and the softmax score vector for each pixel is calculated by the formula (7).

    3.4 Particle Picking Via Connected Component Analysis

    We input the micrographs into the trained Urdnet model to output the prediction results of particles and artifacts.The prediction map is then binarized as shown in Figs.6a and 6b, and the foreground data only contains pixels labeled as “particles,” i.e., green pixels of Fig.6a.Next, the foreground data is analyzed as the connected component [41], and the positions of the centroids of each connected component are taken as the particle positions.Using the connected component analysis, we can reduce interference from false-positive particles and locate the particles more accurately where the particles are close to each other.Toweaken the impact caused by false positives,we remove the particles whose connected components are smaller than 30% of the mean component area.As shown in case 1 of Fig.6, the dissociated particle in the red circle in Fig.6c will not be picked for its small connected components.Where the bounding boxes of two particles overlap in a range of 30% area or more, the particles are considered to be stacked particles.In case 2 of Fig.6, we reject the two particles in the red circles, since the overlapping area of them exceeds the prescribed threshold.

    Figure 6: Particle picking workflow.(a) represents the prediction map, (b) is the corresponding binary image of (a), and (c) represents the particle picking result, where the green circle indicates the selected particles and the red circle indicates the rejected particles

    4 Results and Discussion

    4.1 Datasets

    EMPIAR, the Electron Microscopy Public Image Archive, is a public resource for raw 2D electron microscopy images.Our experimental data includes cryo-EM images of three proteins from EMPIAR, 80S ribosomes [42] (EMPIAR-10153), HCN1 channel [43] (EMPIAR-10081), and TcdA1 toxin subunit [44] (EMPIAR-10089).The three protein datasets from EMPIAR are new datasets rreleased in the last three years.Besides, we test the common benchmark dataset,KLH (Keyhole limpet hemocyanin) [45] collected in 2004 to validate the performance of our automatic particle picking algorithm.Since structural analyses of viruses and large, well-ordered molecules with high pointgroup symmetry have been well known, and particle picking of them is relatively simple tasks, our experimental data mainly focuses on the new three datasets of small and medium protein molecules with no or low symmetry.

    Human 80S ribosomes are used to synthesize proteins in cells and have complex molecular structures.The diameter of the 80S ribosomes is between 25nm and 30 (250-300 ?), and its molecular mass is 3.9-4.5 MDa.The micrographs of EMPIAR-10153 were collected with low defocus and volta phase plate (VPP) [46].VPP improves the phase contrast and SNR in the low-frequency range by introducing an additional phase shift in the unscattered beam.However, VPP not only enhances the contrast of the particles of interest but also the contrast of all weak phase objects,including contamination, ice dregs, carbon film, which increases the difficulty of accurately selecting the particles.

    HCN1 channel is a hyperpolarization-activated cyclic nucleotide-gated ion channel that underlies the control of rhythmic activity in cardiac and neuronal pacemaker cells.Moreover, it forms a structure of the channel tetramer in a ligand-free state with a molecular weight of ~74.6 kDa.Due to its small size, the number and quality of picked particles will affect the final resolution of 3D reconstruction.

    KLH, a highly immunogenic protein macromolecule, is used as a carrier protein for the preparation of immunogens, which has a molecular weight of 7.9Mda and a D5 point-group symmetry with a size of ~40 nm.There are two main types of projection views of the KLH particles, the side view, and the top view.We select both views when picking particles.In addition to the particles of both views, the micrographs consist of clearly overlapping KLH particles and broken particles that make it difficult to extract unbroken single KLH particles.

    Each protein dataset has different parameters such as electronic dose, defocus value, pixel size,and particle size.The main difference between the four datasets is that they have different defocus ranges and electronic dose.The SNRs of the 80S ribosome and KLH are higher for their lowest or none electron dose and low defocus value, and those of the TcdA1 and HCN1 channel are quite low for their higher electron dose and broad defocus range.The specific parameters of these data sets are listed in Tab.1.

    Table 1: Experimental datasets

    Table 1: Continued

    4.2 Model Training

    Our experimental hardware for model training is equipped with NVIDIA GeForce GTX 1080Ti graphics, 64 GB RAM, and Intel Core i7 8700K CPU.We trained the original U-Net and Urdnet models on the EMPIAR datasets of three proteins using the deep learning library Keras based on the Tensorflow backend.The data used for training and testing included 71 human 80S ribosomal micrographs, 30 HCN1 micrographs, and 24 TcdA1micrographs, and each protein contains approximately 2,000-4,000 particles.We set 20% of the training set for validation.To reduce the risk of overfitting training, we augmented the training images by horizontal flip, rotation, width shift, height change, cropping, zooming, and filling, expanding one image to 32 images.The batch size was set to 2, the model optimizer is Adam [47], and the loss function is cross-entropy.The initial learning rate was set to 1E-4; as training progresses, the learning rate gradually decreases.To obtain a generic model for particles of different scales and aspect ratios, we added 16 KLH micrographs to the training dataset and separately trained Urdnet with the raw data and preprocessed data of the four proteins.The loss and accuracy of all trained models are shown in Tab.2.

    Table 2: Model training results

    After being trained 40 epochs, all Urdnet models’losses are less than 0.1, which are lower than those of U-Net, and their accuracies are higher than those of U-Net.It sufficiently proves that the Urdnet model has better classification accuracy than the U-Net model when processing the same protein data.Furthermore, the training loss of the Urdnet generic model is shown in Fig.7.The loss of preprocessed images is always lower than the raw images, indicating the importance of image preprocessing in Section 3.1.For the preprocessed data, the Urdnet generic model has a model loss of less than 0.1 and an accuracy of higher than 0.9, as shown in Tab.2, which indicates that the Urdnet generic model also has excellent performance when dealing with different biomacromolecules.To a certain extent, it reveals the possibility of generalizing the Urdnet generic model to unknown biomacromolecules if more protein data of different particle shapes and sizes are trained.

    Figure 7: The training loss of raw micrographs and preprocessed micrographs

    4.3 Particle Picking

    To evaluate the particle picking performance of Urdnet, we calculated the precision, recall,precision-recall curve, and IoU (Intersection over Union) for the three proteins’test data, and compared them with the DeepPicker method.The results showed in Fig.7 and Tab.3.True Positive means a positive sample predicted by the model as positive, and True Negative means a negative sample predicted as negative.False Positive means a negative sample predicted to be positive, and False Negative means a positive sample predicted as negative.Precision is the proportion of the samples that are correctly predicted to be positive in all predicted positive samples (Eq.(8)).The recall is the proportion of the correctly predicted positive samples in all real positive samples (Eq.(9)),which indicates the ability of Urdnet to detect positive.Different precisions and recalls make up the precision-recall curve by changing the threshold of the particle prediction.As the threshold increases,the accuracy also increases, and the recall rate decreases.For the excellent performance of particle picking, we have to pick out the threshold with the highest recall before the precision drops sharply.We regard manually selected particles as Ground Truth.The mean IoU represents the accuracy of particle position, which defined the ratio of the intersection area of ground truth and the testing result with their united area (Eq.(10)).

    Table 3: Precision, recall, and IOU

    As shown in Fig.8, we find out Urdnet outperformed the DeepPicker at 80S ribosomes, HCN1 channels, and TcdA1 toxin subunits datasets.The 80S ribosome performs best on precision, recall,and IoU for its highest SNR.The precisions of the three datasets are higher than 0.75, the recalls of them are higher than 0.83, and the IoUs of them are higher than 0.81, which fully demonstrates the high performance of our method in particle picking.

    Figure 8: (Continued)

    Figure 8: The precision-recall curve of three protein EMPIAR datasets

    To further evaluate the quality of the particles selected by Urdnet, we compared our method to the semi-automatic selectionmethod of RELION via the entire public datasets.80S ribosomes (EMPIAR- 10153) contain 318 micrographs,HCN1channel (EMPIAR-10081) is the largest dataset, including 997 micrographs, TcdA1 toxin subunit (EMPIAR-10089) contains 97 micrographs, and KLH consists of 82 micrographs from the US National Resource for Automated Molecular Microscopy (NRAMM).2D classification in the RELION software worked on the picked particles to identify suitable 2D average classes.Then, 2D average classes were further 3D classified using a low-resolution 3d map as the initial model, and good templates after 3D classification were kept and refined to construct the final 3Dstructure.We recorded the number of picked particles by the two methods, the “good”particle number related to the 2D class average template and the final resolution of the 3D reconstruction computed by the “gold standard” Fourier shell correlation [48].The results are listed in Tab.4.

    Table 4: Particle auto-picking and 3D reconstruction results

    As shown in Tab.4, the particles’ number of 80S ribosomes, HCN1 channels, and TcdA1 extracted by RELION was approximately 16.8%, 37.1%, and 12.1% more than those ofUrdnet.Due to the small number of the three protein particles used to train Urdnet, it picked fewer particles than RELION at the start.Nonetheless, the number of “good”particles after 2Dclassification is close to that of particles picked by RELION, i.e., the particles picked by two methods ultimately used for 3D reconstruction are quite close.It suggests that Urdnet can identify more true-positive particles and pick fewer falsepositive particles, which makes the final resolution of 3D reconstruction is slightly higher.We think that theUrdnet, different fromRELION’s semi auto-pick job by using themanually selected templates, picked particles with broader angular coverage, which reflects that Urdnet can choose particles with more angles and views.Especially for processing low-SNR and low-contrast datasets such as the HCN1 channel and TcDA1, Urdnet can thoroughly learn the inherent and unique particle features and avoid the dependency of users’ intervention.Besides, Urdnet can significantly improve the accuracy of the picked particle.Since artifacts are labeled at the pixel level as a separate class, Urdnet avoids selecting “bad” particles in large-area artifacts, substantially eliminating the infuences of ice dregs, carbon films, and background noise.The gap between the reconstruction resolutions acquired by Urdnet and RELION is small, which is probably because the reconstructions reached resolutions close to those proteins’ theoretical resolution limit.Nevertheless, the particles picked by Urdnet achieved higher resolutions of 3D reconstruction of the three protein data than RELION.

    For KLH benchmark data, we find that the particles selected by RELION’s semi auto-picking and the particles after 2D classification are slightly more than Urdnet, and its resolution of the final 3d reconstruction is higher than Urdnet.On the one hand, it is because a small number of KLH particles are used in the experiment.The dataset contains 82 micrographs, and each micrograph only has an average of about 25 particles, which results in little difference between the particles selected by the two methods.On the other hand, in recent years, significant changes have taken place in the imaging system of cryo-EM.Micrographs collected today (after 2012) are quite different from the older ones.Old datasets like KLH usually have higher pixel size, higher SNR, and the molecular weights of their particles are often enormous.Traditional methods, such as template matching, can easily detect these large particles.However, as the quality of microscopic images decreases, their performance dramatically decreases.Therefore, there is usually no one method to invariably get the best performance in different data sets.We think our method is quite effective when the image quality is mainly limited by low SNR, or the training data is insufficient.

    4.4 Discussion

    This paper proposes a residual dense convolutional network model based on multiple annotations and improved U-Net for automatic particle picking of cryo-EM biomacromolecules.It can automatically and accurately select particles from cryo-EM images and can reconstruct high-resolution 3D structures of particles from different proteins.But meanwhile, there is still room for improvement in our solution:

    (1) Compared with other methods, it has little advantage in the resolution of 3D reconstruction and is not highly applicable to proteins of different shapes.It is necessary to improve the network model by increasing the types of protein training data.

    (2) The image restoration performance of the residual intensive module is not further analyzed.Setting up experiments to calculate image noise indicators is requisite to evaluate whether the network has additional denoising effects after adding the RDB module.

    Future work will focus on improving the performance of the Urdnet general model to achieve the goal of automatically picking a variety of challenging biological macromolecule cryo-EM image particles without retraining the model, deepening the research on image restoration issues to strive to eliminate data pre-processing steps and enabling the network to autonomously achieve image denoising to facilitate particle feature extraction.

    5 Conclusion

    In this paper, we proposed an automatic particle picking method, Urdnet, based on the U-Net architecture and residual dense block.We introduced a method of combining multiple forms of labels to build training data, which significantly releases the burden of manual labeling and improves annotation efficiency.Urdnet demonstrates its excellent performance on the public data of 80S ribosomal, HCN1 channel, TcdA1 toxin subunit, and KLH that it can effortlessly process cryo-EM particles from different proteins.Compared with the DeepPicker method, Urdnet achieves higher precision, recall, and IoU in particle picking.Compared with RELION’s semi-automatic selection,our method can achieve higher resolutions of 3D reconstruction at most data.In future work, we will focus on improving the performance of the Urdnet generic model to pick new protein particles of multiple morphologies without the need to retrain the model.

    Acknowledgement:The authors thank Yi He for his valuable suggestions on cryo-EM data preprocessing and Xiaoran Yu for her assistance with training data annotation.

    Funding Statement:This research was supported by Key Projects of the Ministry of Science and Technologyofthe People’s Republicof China(2018AAA0102301);the Open Research Fundof Hunan Provincial Key Laboratory of Network Investigational Technology, Grant No.2018WLZC001.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    成年女人看的毛片在线观看| 边亲边吃奶的免费视频| 久久鲁丝午夜福利片| 不卡视频在线观看欧美| 亚洲在线观看片| 国产69精品久久久久777片| 禁无遮挡网站| 综合色丁香网| 欧美 日韩 精品 国产| 欧美zozozo另类| 欧美性猛交╳xxx乱大交人| 免费电影在线观看免费观看| 久久精品久久久久久噜噜老黄| 亚洲欧美精品自产自拍| 国产黄片美女视频| 亚洲在久久综合| 亚洲丝袜综合中文字幕| 99热这里只有是精品在线观看| 欧美成人a在线观看| 国产精品爽爽va在线观看网站| 97超视频在线观看视频| 久久久久久久久久黄片| 亚洲av福利一区| 淫秽高清视频在线观看| 欧美另类一区| 一级片'在线观看视频| 亚洲av男天堂| 久久精品久久久久久噜噜老黄| 成人国产麻豆网| 亚洲欧美成人综合另类久久久| 看非洲黑人一级黄片| 亚洲精品自拍成人| ponron亚洲| 两个人的视频大全免费| 大话2 男鬼变身卡| 精品熟女少妇av免费看| 高清午夜精品一区二区三区| 夫妻性生交免费视频一级片| 在线a可以看的网站| 全区人妻精品视频| 国产有黄有色有爽视频| 最近中文字幕2019免费版| 成人亚洲精品一区在线观看 | 免费观看a级毛片全部| 亚洲精品自拍成人| 我的老师免费观看完整版| 永久免费av网站大全| 国产精品综合久久久久久久免费| 熟女电影av网| 特级一级黄色大片| 亚洲无线观看免费| 天天一区二区日本电影三级| 国产一区二区三区综合在线观看 | 日韩av在线免费看完整版不卡| 寂寞人妻少妇视频99o| 亚洲,欧美,日韩| 国产综合精华液| 免费人成在线观看视频色| 亚洲电影在线观看av| 99九九线精品视频在线观看视频| 91狼人影院| 嫩草影院入口| 麻豆精品久久久久久蜜桃| 欧美高清性xxxxhd video| 亚洲精品中文字幕在线视频 | 欧美日韩精品成人综合77777| 久久热精品热| h日本视频在线播放| 少妇的逼水好多| 国产乱人偷精品视频| 色视频www国产| 韩国av在线不卡| 国产精品嫩草影院av在线观看| 激情 狠狠 欧美| 亚洲内射少妇av| 亚洲美女视频黄频| av播播在线观看一区| 大话2 男鬼变身卡| 人人妻人人澡人人爽人人夜夜 | 99视频精品全部免费 在线| 天天一区二区日本电影三级| 国产亚洲精品av在线| 精品国内亚洲2022精品成人| 免费无遮挡裸体视频| 精品熟女少妇av免费看| 成年版毛片免费区| 免费看光身美女| 18+在线观看网站| 国产色婷婷99| 国产乱人偷精品视频| 18禁动态无遮挡网站| 国产成人aa在线观看| 国产精品熟女久久久久浪| 好男人视频免费观看在线| 日韩欧美精品免费久久| 日本与韩国留学比较| 永久网站在线| 国产老妇女一区| www.av在线官网国产| 高清视频免费观看一区二区 | 久久人人爽人人爽人人片va| 又粗又硬又长又爽又黄的视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 超碰97精品在线观看| 最近的中文字幕免费完整| 久久久久精品性色| 午夜免费激情av| 国模一区二区三区四区视频| 久久6这里有精品| 亚洲欧美日韩卡通动漫| 我要看日韩黄色一级片| 精品午夜福利在线看| 免费电影在线观看免费观看| 又大又黄又爽视频免费| 日本与韩国留学比较| 国产亚洲5aaaaa淫片| 国产精品国产三级专区第一集| 日韩精品青青久久久久久| 1000部很黄的大片| 欧美xxⅹ黑人| 天堂中文最新版在线下载 | 99热6这里只有精品| 久热久热在线精品观看| 男人狂女人下面高潮的视频| 嘟嘟电影网在线观看| 精品99又大又爽又粗少妇毛片| 内地一区二区视频在线| 精品亚洲乱码少妇综合久久| 国产色爽女视频免费观看| 亚洲在线观看片| 国产伦精品一区二区三区四那| 在线播放无遮挡| 熟妇人妻久久中文字幕3abv| 天天一区二区日本电影三级| 九草在线视频观看| 黄色配什么色好看| 三级男女做爰猛烈吃奶摸视频| 99视频精品全部免费 在线| 国产免费一级a男人的天堂| 男人舔奶头视频| 五月伊人婷婷丁香| 少妇丰满av| 草草在线视频免费看| 国产精品久久久久久精品电影| 一级毛片 在线播放| av又黄又爽大尺度在线免费看| 亚州av有码| 久久精品综合一区二区三区| 欧美一区二区亚洲| 国产高清有码在线观看视频| 午夜亚洲福利在线播放| 亚洲aⅴ乱码一区二区在线播放| 天天躁日日操中文字幕| 国产亚洲av嫩草精品影院| 中文乱码字字幕精品一区二区三区 | 性插视频无遮挡在线免费观看| 欧美xxxx黑人xx丫x性爽| 亚洲欧美中文字幕日韩二区| 51国产日韩欧美| 亚洲成人一二三区av| 九九爱精品视频在线观看| 天堂影院成人在线观看| 国产 一区 欧美 日韩| av在线老鸭窝| 99热这里只有是精品50| 欧美变态另类bdsm刘玥| av在线亚洲专区| 夫妻性生交免费视频一级片| 两个人的视频大全免费| 人人妻人人澡欧美一区二区| 黄片wwwwww| 五月玫瑰六月丁香| 18禁在线播放成人免费| 一区二区三区免费毛片| 成人无遮挡网站| 麻豆成人午夜福利视频| 婷婷色av中文字幕| 中文在线观看免费www的网站| 亚洲在线自拍视频| 女人久久www免费人成看片| 国产老妇女一区| 国产精品一区二区性色av| 亚洲乱码一区二区免费版| 亚洲怡红院男人天堂| 精品久久国产蜜桃| 国产黄色小视频在线观看| 亚洲国产最新在线播放| 亚洲成人久久爱视频| 午夜福利成人在线免费观看| 亚洲在线观看片| 不卡视频在线观看欧美| 天堂中文最新版在线下载 | 全区人妻精品视频| 日产精品乱码卡一卡2卡三| 99久久精品热视频| 国产综合精华液| 亚洲精品国产av成人精品| 国产精品美女特级片免费视频播放器| 男人狂女人下面高潮的视频| 欧美极品一区二区三区四区| 国产色婷婷99| 精品少妇黑人巨大在线播放| 久久久久久久久久久免费av| 大陆偷拍与自拍| 熟女电影av网| 日韩精品青青久久久久久| 日韩中字成人| 国产一区亚洲一区在线观看| or卡值多少钱| 成人无遮挡网站| 精品一区二区三区视频在线| 能在线免费观看的黄片| 欧美日韩亚洲高清精品| 丝袜喷水一区| 亚洲国产日韩欧美精品在线观看| 成人亚洲精品一区在线观看 | 中文欧美无线码| 在线观看一区二区三区| 黑人高潮一二区| 亚洲av电影在线观看一区二区三区 | 人人妻人人澡人人爽人人夜夜 | 亚洲欧洲国产日韩| 中文精品一卡2卡3卡4更新| 日韩成人av中文字幕在线观看| 91精品伊人久久大香线蕉| 亚洲欧美中文字幕日韩二区| 99热这里只有是精品在线观看| 18禁动态无遮挡网站| a级一级毛片免费在线观看| 白带黄色成豆腐渣| 亚洲精品视频女| 99久久中文字幕三级久久日本| 日韩av在线大香蕉| 青春草国产在线视频| 国产色婷婷99| 国产老妇女一区| 日韩精品有码人妻一区| 青春草亚洲视频在线观看| 国产精品一区www在线观看| 国产色爽女视频免费观看| h日本视频在线播放| 国产 一区精品| 国产亚洲午夜精品一区二区久久 | 毛片一级片免费看久久久久| 青春草国产在线视频| 日韩电影二区| 肉色欧美久久久久久久蜜桃 | 少妇人妻一区二区三区视频| av一本久久久久| 久久久色成人| 中文字幕人妻熟人妻熟丝袜美| 成年版毛片免费区| 日韩欧美精品v在线| or卡值多少钱| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 人人妻人人澡人人爽人人夜夜 | 久久6这里有精品| 超碰97精品在线观看| 国产又色又爽无遮挡免| 午夜福利网站1000一区二区三区| 高清日韩中文字幕在线| 一级毛片久久久久久久久女| 午夜激情福利司机影院| 精品久久久噜噜| 成年女人看的毛片在线观看| 色综合亚洲欧美另类图片| 卡戴珊不雅视频在线播放| 晚上一个人看的免费电影| 亚洲三级黄色毛片| 精品国产露脸久久av麻豆 | 欧美一区二区亚洲| 国产精品熟女久久久久浪| 伊人久久国产一区二区| 国产黄片美女视频| 国产一区二区三区av在线| 国国产精品蜜臀av免费| 日韩欧美精品免费久久| 亚洲精品国产av蜜桃| 国产精品人妻久久久久久| 国产女主播在线喷水免费视频网站 | 夫妻午夜视频| 18禁裸乳无遮挡免费网站照片| 国产精品伦人一区二区| 天堂av国产一区二区熟女人妻| 精品熟女少妇av免费看| 九色成人免费人妻av| 欧美日韩国产mv在线观看视频 | 亚洲国产精品sss在线观看| 国产精品日韩av在线免费观看| 国产亚洲av嫩草精品影院| 国产在线一区二区三区精| 天天一区二区日本电影三级| 日产精品乱码卡一卡2卡三| 久久久国产一区二区| 久久午夜福利片| 亚洲av不卡在线观看| 欧美日本视频| 麻豆久久精品国产亚洲av| 国产有黄有色有爽视频| 少妇丰满av| 国产高清国产精品国产三级 | 三级国产精品欧美在线观看| 国产成人免费观看mmmm| 嫩草影院入口| 国产大屁股一区二区在线视频| 看十八女毛片水多多多| 日韩中字成人| 国产精品精品国产色婷婷| 99热6这里只有精品| 国产在线男女| 纵有疾风起免费观看全集完整版 | 一区二区三区免费毛片| 亚洲人成网站在线播| 搞女人的毛片| 久久99蜜桃精品久久| 中国国产av一级| 麻豆精品久久久久久蜜桃| 免费看不卡的av| 国产爱豆传媒在线观看| 亚洲av成人av| 一本久久精品| 国产亚洲av片在线观看秒播厂 | 国内揄拍国产精品人妻在线| 日本免费a在线| 又爽又黄无遮挡网站| 麻豆成人av视频| 黄片无遮挡物在线观看| 亚洲精品国产成人久久av| 51国产日韩欧美| 日韩欧美一区视频在线观看 | 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va| 国产亚洲精品久久久com| 亚洲av成人精品一区久久| 一级毛片黄色毛片免费观看视频| 菩萨蛮人人尽说江南好唐韦庄| av卡一久久| 成年av动漫网址| 国产毛片a区久久久久| 3wmmmm亚洲av在线观看| 日日干狠狠操夜夜爽| 亚洲精品456在线播放app| 91狼人影院| ponron亚洲| 精品一区二区三区人妻视频| 一级av片app| 亚洲av国产av综合av卡| 国内精品宾馆在线| 日韩精品有码人妻一区| 亚洲内射少妇av| 男人爽女人下面视频在线观看| 成人美女网站在线观看视频| 永久免费av网站大全| 国产一区亚洲一区在线观看| 99久国产av精品国产电影| 真实男女啪啪啪动态图| 亚洲第一区二区三区不卡| 又粗又硬又长又爽又黄的视频| 在线免费十八禁| 亚洲精品一二三| av在线观看视频网站免费| 国产精品精品国产色婷婷| 久久久久性生活片| 天堂av国产一区二区熟女人妻| 国产成人精品久久久久久| 99久久精品热视频| 国产亚洲精品久久久com| 亚洲在线观看片| 国产精品一区二区在线观看99 | or卡值多少钱| 水蜜桃什么品种好| 国产亚洲精品av在线| 久久久午夜欧美精品| 亚洲最大成人av| 国产69精品久久久久777片| 高清视频免费观看一区二区 | 春色校园在线视频观看| 亚洲av成人精品一区久久| 日本-黄色视频高清免费观看| xxx大片免费视频| 国产精品久久久久久精品电影小说 | 国产一区二区三区综合在线观看 | 国产高潮美女av| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产av成人精品| 天堂影院成人在线观看| av.在线天堂| 九色成人免费人妻av| 午夜激情欧美在线| 中文资源天堂在线| 亚洲熟女精品中文字幕| 日韩av在线大香蕉| 又爽又黄a免费视频| 亚洲av日韩在线播放| 欧美精品国产亚洲| 最近的中文字幕免费完整| 日韩成人伦理影院| 日韩精品青青久久久久久| 啦啦啦中文免费视频观看日本| 精品99又大又爽又粗少妇毛片| h日本视频在线播放| 国产精品一及| 毛片女人毛片| 国产伦精品一区二区三区四那| 2022亚洲国产成人精品| 国产av国产精品国产| 别揉我奶头 嗯啊视频| 最近中文字幕2019免费版| 国产精品一区二区三区四区久久| 能在线免费观看的黄片| 日韩欧美一区视频在线观看 | 中文精品一卡2卡3卡4更新| 能在线免费观看的黄片| 欧美人与善性xxx| 日韩精品有码人妻一区| 欧美成人一区二区免费高清观看| 国产乱人偷精品视频| 国产免费福利视频在线观看| 亚洲精品色激情综合| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看| 看非洲黑人一级黄片| 男人和女人高潮做爰伦理| 国产亚洲精品av在线| 婷婷六月久久综合丁香| 男人和女人高潮做爰伦理| 欧美日本视频| 亚洲伊人久久精品综合| 亚洲国产精品成人久久小说| 十八禁网站网址无遮挡 | 精品午夜福利在线看| 人妻少妇偷人精品九色| 国产精品99久久久久久久久| 日本wwww免费看| 中文乱码字字幕精品一区二区三区 | 精品酒店卫生间| 一级毛片 在线播放| 在线免费十八禁| 国产伦理片在线播放av一区| 亚洲熟女精品中文字幕| 免费观看无遮挡的男女| 中文字幕人妻熟人妻熟丝袜美| 国产成人a区在线观看| 色网站视频免费| 亚洲精品乱久久久久久| 麻豆乱淫一区二区| 黄色配什么色好看| 色吧在线观看| 亚洲欧美一区二区三区国产| 内射极品少妇av片p| 亚洲伊人久久精品综合| 久久精品夜色国产| 伦理电影大哥的女人| 日韩三级伦理在线观看| 免费在线观看成人毛片| 欧美不卡视频在线免费观看| 国产精品人妻久久久影院| 97精品久久久久久久久久精品| 免费av观看视频| 日日撸夜夜添| 国产精品一区www在线观看| 人妻制服诱惑在线中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜日本视频在线| 两个人视频免费观看高清| 日韩视频在线欧美| 亚洲精品乱码久久久久久按摩| 两个人视频免费观看高清| 国产在线一区二区三区精| ponron亚洲| 久久97久久精品| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 久久久精品欧美日韩精品| 婷婷色av中文字幕| 熟女电影av网| 国产精品爽爽va在线观看网站| 高清毛片免费看| 日韩欧美 国产精品| 亚洲av一区综合| 尤物成人国产欧美一区二区三区| 国产精品综合久久久久久久免费| 亚洲欧洲日产国产| 精品久久久久久成人av| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站| 成人鲁丝片一二三区免费| 国产精品一区二区在线观看99 | 大片免费播放器 马上看| 80岁老熟妇乱子伦牲交| 91aial.com中文字幕在线观看| 大片免费播放器 马上看| 日韩视频在线欧美| 午夜福利在线在线| 日本免费a在线| 亚洲在线自拍视频| 91av网一区二区| 国产亚洲精品久久久com| 尾随美女入室| 国产 亚洲一区二区三区 | 99热6这里只有精品| 看免费成人av毛片| 国产亚洲午夜精品一区二区久久 | 亚洲一区高清亚洲精品| 看十八女毛片水多多多| 亚洲国产欧美人成| 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| 欧美成人a在线观看| 少妇的逼水好多| 菩萨蛮人人尽说江南好唐韦庄| 午夜激情福利司机影院| 国产高清不卡午夜福利| 青春草亚洲视频在线观看| 直男gayav资源| 亚洲精品视频女| 午夜老司机福利剧场| 亚洲性久久影院| 黑人高潮一二区| 高清毛片免费看| 色综合亚洲欧美另类图片| 国产69精品久久久久777片| 日本免费在线观看一区| 卡戴珊不雅视频在线播放| 波多野结衣巨乳人妻| 国产一区亚洲一区在线观看| 晚上一个人看的免费电影| 亚洲精品亚洲一区二区| 久久久久网色| 全区人妻精品视频| 亚洲国产精品成人综合色| 人妻一区二区av| 一级a做视频免费观看| 亚洲美女搞黄在线观看| 亚洲精品,欧美精品| 久久久久网色| 在线a可以看的网站| av在线播放精品| 特级一级黄色大片| 噜噜噜噜噜久久久久久91| 免费人成在线观看视频色| 亚洲av成人精品一区久久| 看非洲黑人一级黄片| 色视频www国产| 国产乱人视频| 国产淫片久久久久久久久| 99视频精品全部免费 在线| 久久韩国三级中文字幕| 伊人久久精品亚洲午夜| av在线天堂中文字幕| 亚洲丝袜综合中文字幕| xxx大片免费视频| 观看免费一级毛片| 国产成人免费观看mmmm| 丝瓜视频免费看黄片| 男人舔女人下体高潮全视频| 亚洲一区高清亚洲精品| 汤姆久久久久久久影院中文字幕 | 亚洲国产色片| kizo精华| 最近的中文字幕免费完整| 色综合色国产| 超碰av人人做人人爽久久| 97超碰精品成人国产| 男女边吃奶边做爰视频| 欧美97在线视频| 亚洲精品乱久久久久久| av线在线观看网站| 男女下面进入的视频免费午夜| 久久草成人影院| 日本色播在线视频| 国产v大片淫在线免费观看| 亚洲自拍偷在线| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区| 三级毛片av免费| 日韩一本色道免费dvd| 久久久久网色| 亚洲成色77777| 天天躁日日操中文字幕| 国产一区二区三区综合在线观看 | 亚州av有码| 国产色爽女视频免费观看| 97超视频在线观看视频| 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 国产成人a区在线观看| av又黄又爽大尺度在线免费看| 观看免费一级毛片| 日韩三级伦理在线观看| 午夜精品一区二区三区免费看| 国产成年人精品一区二区| 一级爰片在线观看| 亚洲18禁久久av| 国产日韩欧美在线精品| 国产成人精品婷婷| 国产高清三级在线| 中文精品一卡2卡3卡4更新| 联通29元200g的流量卡| a级毛片免费高清观看在线播放| 国产探花在线观看一区二区| 国产av不卡久久| 成年女人在线观看亚洲视频 | 男人舔女人下体高潮全视频| 蜜桃亚洲精品一区二区三区| 男女下面进入的视频免费午夜| 午夜精品国产一区二区电影 | 内射极品少妇av片p| 久久综合国产亚洲精品| 高清视频免费观看一区二区 | 国产成人精品一,二区| 九色成人免费人妻av| 成人午夜精彩视频在线观看| 淫秽高清视频在线观看|