• <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.

    久久久久久久午夜电影| 一级爰片在线观看| 小蜜桃在线观看免费完整版高清| 高清视频免费观看一区二区 | 一级a做视频免费观看| or卡值多少钱| 国模一区二区三区四区视频| 性插视频无遮挡在线免费观看| 亚洲四区av| 成人av在线播放网站| 少妇熟女欧美另类| 亚洲国产成人一精品久久久| av线在线观看网站| 亚洲精品一二三| 99久久九九国产精品国产免费| 精品熟女少妇av免费看| 日本午夜av视频| 日韩欧美三级三区| 久久久久久久久大av| 亚洲图色成人| 久久久久精品性色| 久久久久久久大尺度免费视频| 国国产精品蜜臀av免费| 校园人妻丝袜中文字幕| 啦啦啦啦在线视频资源| 成人一区二区视频在线观看| 人妻制服诱惑在线中文字幕| 久久久久久九九精品二区国产| freevideosex欧美| 亚洲精品亚洲一区二区| 国产成年人精品一区二区| 亚洲av电影在线观看一区二区三区 | 久久热精品热| 亚洲伊人久久精品综合| 亚洲av一区综合| 美女内射精品一级片tv| 麻豆精品久久久久久蜜桃| 禁无遮挡网站| 欧美 日韩 精品 国产| 夜夜爽夜夜爽视频| 亚洲av二区三区四区| 男的添女的下面高潮视频| 美女主播在线视频| 最近最新中文字幕免费大全7| eeuss影院久久| 日韩 亚洲 欧美在线| 两个人视频免费观看高清| 国产高清国产精品国产三级 | av国产免费在线观看| 寂寞人妻少妇视频99o| 尾随美女入室| 99久久人妻综合| 一级黄片播放器| 久久久色成人| 久久久久久久久久久丰满| 成年女人在线观看亚洲视频 | 高清在线视频一区二区三区| 亚洲av中文字字幕乱码综合| 一个人看的www免费观看视频| 国内少妇人妻偷人精品xxx网站| 欧美bdsm另类| 搡老妇女老女人老熟妇| 国产av国产精品国产| 午夜激情福利司机影院| 亚洲久久久久久中文字幕| 中文字幕av在线有码专区| 亚洲成色77777| 男人舔奶头视频| 国产亚洲5aaaaa淫片| 欧美三级亚洲精品| 午夜福利视频1000在线观看| 欧美3d第一页| 日韩在线高清观看一区二区三区| 观看免费一级毛片| ponron亚洲| 亚洲一级一片aⅴ在线观看| 免费看美女性在线毛片视频| 国产精品熟女久久久久浪| 秋霞在线观看毛片| 天堂网av新在线| 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 国产男人的电影天堂91| 午夜老司机福利剧场| 亚洲国产精品成人久久小说| 久热久热在线精品观看| 天堂av国产一区二区熟女人妻| 日韩,欧美,国产一区二区三区| 两个人的视频大全免费| 高清日韩中文字幕在线| 七月丁香在线播放| 成人国产麻豆网| 亚洲最大成人av| 一级毛片我不卡| 亚洲成人一二三区av| 国产伦精品一区二区三区四那| 日韩人妻高清精品专区| 日韩成人av中文字幕在线观看| 蜜桃亚洲精品一区二区三区| 国产伦一二天堂av在线观看| 亚洲精品久久久久久婷婷小说| 97在线视频观看| 男插女下体视频免费在线播放| 免费看日本二区| 可以在线观看毛片的网站| 人人妻人人澡人人爽人人夜夜 | 久久久久久伊人网av| 舔av片在线| 午夜福利在线观看吧| 嫩草影院入口| 十八禁国产超污无遮挡网站| 国产精品久久久久久久电影| www.色视频.com| 久久97久久精品| 亚洲精品成人av观看孕妇| 夜夜爽夜夜爽视频| 三级经典国产精品| 七月丁香在线播放| 夫妻午夜视频| 国精品久久久久久国模美| 99热全是精品| 国产在线男女| 亚洲欧美精品自产自拍| 男女下面进入的视频免费午夜| 日韩av在线免费看完整版不卡| 亚洲av成人av| 亚洲精品日本国产第一区| av线在线观看网站| 黄色一级大片看看| 免费高清在线观看视频在线观看| 婷婷六月久久综合丁香| 亚洲精品国产av成人精品| 日本三级黄在线观看| 欧美激情国产日韩精品一区| 黄片无遮挡物在线观看| 少妇的逼好多水| 卡戴珊不雅视频在线播放| 国产精品久久久久久久久免| 日本熟妇午夜| av国产免费在线观看| 偷拍熟女少妇极品色| 在线 av 中文字幕| 久久亚洲国产成人精品v| 一区二区三区高清视频在线| 国产亚洲午夜精品一区二区久久 | 一区二区三区乱码不卡18| 亚洲精品aⅴ在线观看| 亚洲综合色惰| 99久久精品国产国产毛片| a级毛色黄片| 中文字幕av在线有码专区| 国产高潮美女av| 国产精品一及| 一区二区三区乱码不卡18| 国产亚洲精品久久久com| 成人午夜精彩视频在线观看| 三级国产精品欧美在线观看| 国产男人的电影天堂91| 免费观看的影片在线观看| 午夜福利成人在线免费观看| 亚洲自偷自拍三级| 色播亚洲综合网| 久久久久久久大尺度免费视频| 欧美97在线视频| 肉色欧美久久久久久久蜜桃 | 男女边吃奶边做爰视频| 久久精品夜色国产| 中国美白少妇内射xxxbb| 国产免费一级a男人的天堂| 久久精品夜夜夜夜夜久久蜜豆| 六月丁香七月| 欧美高清性xxxxhd video| 日本熟妇午夜| 精品一区二区免费观看| 春色校园在线视频观看| 精品亚洲乱码少妇综合久久| 成人毛片60女人毛片免费| 亚洲av福利一区| 18禁在线无遮挡免费观看视频| 成年免费大片在线观看| 国产一区有黄有色的免费视频 | 亚洲欧美成人综合另类久久久| 搡老妇女老女人老熟妇| 狠狠精品人妻久久久久久综合| 国内精品一区二区在线观看| 婷婷六月久久综合丁香| 乱人视频在线观看| 日韩欧美国产在线观看| 日韩欧美 国产精品| 亚洲国产欧美在线一区| 男女视频在线观看网站免费| 97超碰精品成人国产| 日产精品乱码卡一卡2卡三| 精品一区在线观看国产| 色视频www国产| 国产免费又黄又爽又色| 亚洲精品,欧美精品| 超碰av人人做人人爽久久| 非洲黑人性xxxx精品又粗又长| 国产伦精品一区二区三区四那| 国产永久视频网站| 少妇的逼好多水| 婷婷色av中文字幕| 网址你懂的国产日韩在线| 日本午夜av视频| 午夜亚洲福利在线播放| 亚洲久久久久久中文字幕| 午夜视频国产福利| av在线蜜桃| 欧美另类一区| 99热这里只有是精品在线观看| 熟妇人妻久久中文字幕3abv| 五月伊人婷婷丁香| 在线观看美女被高潮喷水网站| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 亚洲在线观看片| 天天躁夜夜躁狠狠久久av| 男女国产视频网站| 亚洲av电影不卡..在线观看| 欧美精品一区二区大全| 国产黄a三级三级三级人| 午夜爱爱视频在线播放| 久久久久免费精品人妻一区二区| 男人舔女人下体高潮全视频| 亚洲欧美一区二区三区国产| 内射极品少妇av片p| 少妇丰满av| 国产黄色小视频在线观看| 国产精品熟女久久久久浪| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 99久久精品国产国产毛片| 亚洲在线观看片| 不卡视频在线观看欧美| 91久久精品国产一区二区三区| 老司机影院成人| 极品少妇高潮喷水抽搐| 91久久精品电影网| 97在线视频观看| 国产精品一二三区在线看| 人妻一区二区av| 久久久久久伊人网av| 欧美丝袜亚洲另类| av天堂中文字幕网| 国产黄片视频在线免费观看| 中文字幕av在线有码专区| av免费在线看不卡| 欧美日韩亚洲高清精品| 国产黄色免费在线视频| av在线播放精品| 免费看av在线观看网站| 欧美bdsm另类| 波多野结衣巨乳人妻| 午夜福利高清视频| 日本爱情动作片www.在线观看| 青春草国产在线视频| 中国美白少妇内射xxxbb| av网站免费在线观看视频 | 欧美xxxx性猛交bbbb| 亚洲精品一二三| av国产免费在线观看| 欧美成人一区二区免费高清观看| 国产高清三级在线| 3wmmmm亚洲av在线观看| 国产老妇伦熟女老妇高清| 最近最新中文字幕大全电影3| 久久久久久久久久久免费av| 亚洲国产精品专区欧美| 在线观看美女被高潮喷水网站| 狠狠精品人妻久久久久久综合| 亚洲电影在线观看av| 日韩av免费高清视频| 免费看av在线观看网站| 亚洲国产高清在线一区二区三| 大话2 男鬼变身卡| 99久国产av精品| 国产精品.久久久| 日韩欧美精品v在线| 国产淫片久久久久久久久| av线在线观看网站| 日产精品乱码卡一卡2卡三| 久久6这里有精品| 亚洲精品一区蜜桃| 肉色欧美久久久久久久蜜桃 | 51国产日韩欧美| 国产av码专区亚洲av| 搡老妇女老女人老熟妇| 精品少妇黑人巨大在线播放| 乱人视频在线观看| 亚洲av中文字字幕乱码综合| 国产亚洲午夜精品一区二区久久 | 身体一侧抽搐| 国语对白做爰xxxⅹ性视频网站| 日韩国内少妇激情av| 国产亚洲午夜精品一区二区久久 | 中文字幕制服av| 成年版毛片免费区| av免费观看日本| 免费观看性生交大片5| 亚洲国产精品国产精品| 亚洲精品aⅴ在线观看| 超碰av人人做人人爽久久| 亚洲精品久久午夜乱码| 蜜桃久久精品国产亚洲av| 精品一区二区三区人妻视频| 久久精品久久精品一区二区三区| 亚洲三级黄色毛片| 日本-黄色视频高清免费观看| 深爱激情五月婷婷| 伊人久久国产一区二区| 日韩一本色道免费dvd| 搡老妇女老女人老熟妇| 韩国av在线不卡| h日本视频在线播放| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 日日啪夜夜爽| 精华霜和精华液先用哪个| 免费观看的影片在线观看| 中文资源天堂在线| 黄片wwwwww| 麻豆av噜噜一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 男人和女人高潮做爰伦理| 尾随美女入室| 久久这里只有精品中国| 久久久久九九精品影院| 一个人看视频在线观看www免费| 亚洲av.av天堂| 在现免费观看毛片| 国产有黄有色有爽视频| 免费高清在线观看视频在线观看| 免费观看无遮挡的男女| 美女被艹到高潮喷水动态| 免费观看在线日韩| 久久韩国三级中文字幕| 亚洲精品一二三| 免费电影在线观看免费观看| 国产精品国产三级国产专区5o| av天堂中文字幕网| 一级毛片久久久久久久久女| 99热这里只有精品一区| 伦理电影大哥的女人| 2022亚洲国产成人精品| 色播亚洲综合网| 赤兔流量卡办理| 最新中文字幕久久久久| 五月天丁香电影| 久久99热6这里只有精品| 日日啪夜夜撸| 最近手机中文字幕大全| 欧美高清成人免费视频www| 久久久午夜欧美精品| 美女内射精品一级片tv| 免费av观看视频| 高清日韩中文字幕在线| 久热久热在线精品观看| 肉色欧美久久久久久久蜜桃 | 国产高清有码在线观看视频| 少妇熟女欧美另类| 成人欧美大片| 欧美激情国产日韩精品一区| 亚洲国产精品成人综合色| 欧美极品一区二区三区四区| 亚洲无线观看免费| 成年版毛片免费区| 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久v下载方式| www.av在线官网国产| 插逼视频在线观看| 亚洲精品久久久久久婷婷小说| 最近最新中文字幕免费大全7| 99九九线精品视频在线观看视频| 成人av在线播放网站| 精品一区二区三卡| 亚洲电影在线观看av| 精品一区二区免费观看| 精品久久久久久电影网| 最近中文字幕高清免费大全6| 美女cb高潮喷水在线观看| 免费看av在线观看网站| 成年av动漫网址| 人妻一区二区av| 国产精品国产三级专区第一集| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 青春草国产在线视频| 欧美xxⅹ黑人| 国产高清有码在线观看视频| 白带黄色成豆腐渣| 亚洲av电影在线观看一区二区三区 | 久久久久久国产a免费观看| 国产大屁股一区二区在线视频| 国产男女超爽视频在线观看| 99热网站在线观看| 精品久久久久久成人av| 只有这里有精品99| 亚洲国产色片| 看黄色毛片网站| 欧美bdsm另类| 日本黄色片子视频| ponron亚洲| 一区二区三区乱码不卡18| 看黄色毛片网站| 成人国产麻豆网| 亚洲精品一区蜜桃| 黄片wwwwww| 51国产日韩欧美| 亚洲成色77777| 国模一区二区三区四区视频| 国内精品一区二区在线观看| 欧美一区二区亚洲| 日韩电影二区| 又大又黄又爽视频免费| 午夜久久久久精精品| 日韩精品有码人妻一区| 久久99热6这里只有精品| av天堂中文字幕网| 日日干狠狠操夜夜爽| 国产成人午夜福利电影在线观看| 亚洲av福利一区| 欧美日韩亚洲高清精品| 日韩av在线免费看完整版不卡| 五月天丁香电影| 国产男女超爽视频在线观看| 国产男人的电影天堂91| 男人狂女人下面高潮的视频| 国产黄频视频在线观看| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 舔av片在线| 亚洲av中文字字幕乱码综合| av又黄又爽大尺度在线免费看| 干丝袜人妻中文字幕| 在线观看一区二区三区| 国产v大片淫在线免费观看| 久久久午夜欧美精品| 成年人午夜在线观看视频 | 美女国产视频在线观看| 中文字幕制服av| 亚洲乱码一区二区免费版| 成年人午夜在线观看视频 | 五月伊人婷婷丁香| 免费观看性生交大片5| 91精品一卡2卡3卡4卡| 日本猛色少妇xxxxx猛交久久| 日本三级黄在线观看| 街头女战士在线观看网站| 精品久久久久久久末码| 永久免费av网站大全| 最近2019中文字幕mv第一页| 亚洲自拍偷在线| 亚洲不卡免费看| 久久久久久久久大av| 中文字幕久久专区| 搡老乐熟女国产| 久久精品国产亚洲av涩爱| 亚洲精品影视一区二区三区av| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 2022亚洲国产成人精品| 国精品久久久久久国模美| 国产精品不卡视频一区二区| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆| 日本黄大片高清| 久久久久久久午夜电影| 国产精品麻豆人妻色哟哟久久 | 嘟嘟电影网在线观看| 神马国产精品三级电影在线观看| 国产精品不卡视频一区二区| 亚洲成人av在线免费| 国内精品宾馆在线| 插逼视频在线观看| 亚洲成人久久爱视频| 欧美日韩在线观看h| 国产乱人偷精品视频| 国产黄a三级三级三级人| 国产精品99久久久久久久久| 亚洲高清免费不卡视频| 22中文网久久字幕| 久热久热在线精品观看| 亚洲精品色激情综合| 亚洲精品一二三| 一个人观看的视频www高清免费观看| 亚洲欧美精品自产自拍| 美女高潮的动态| 日韩不卡一区二区三区视频在线| 啦啦啦韩国在线观看视频| 少妇裸体淫交视频免费看高清| 久久久久久久久久人人人人人人| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩无卡精品| 十八禁国产超污无遮挡网站| 亚洲av电影在线观看一区二区三区 | 色5月婷婷丁香| 久久韩国三级中文字幕| 久久久久久九九精品二区国产| 久热久热在线精品观看| av又黄又爽大尺度在线免费看| 日韩成人伦理影院| 大片免费播放器 马上看| 国产精品综合久久久久久久免费| 久久久久久久大尺度免费视频| 国产精品人妻久久久影院| 亚洲无线观看免费| 亚洲精品国产av成人精品| 蜜臀久久99精品久久宅男| 熟妇人妻不卡中文字幕| 国产国拍精品亚洲av在线观看| 亚洲无线观看免费| 午夜免费激情av| 亚洲精品自拍成人| 精品久久久久久成人av| 色哟哟·www| 亚洲最大成人中文| 亚洲av免费高清在线观看| 中文欧美无线码| 天堂影院成人在线观看| 综合色av麻豆| 欧美日韩综合久久久久久| 国产中年淑女户外野战色| 久久久久精品久久久久真实原创| 99久久精品一区二区三区| 69人妻影院| av在线蜜桃| 天堂俺去俺来也www色官网 | 亚洲久久久久久中文字幕| 国产黄a三级三级三级人| 国产成人精品福利久久| 国产片特级美女逼逼视频| 五月伊人婷婷丁香| 国产淫语在线视频| 国产av国产精品国产| 国产高潮美女av| 久久久久久久国产电影| 国国产精品蜜臀av免费| 免费av不卡在线播放| 秋霞在线观看毛片| 免费电影在线观看免费观看| 国产伦在线观看视频一区| 最近手机中文字幕大全| 欧美日韩在线观看h| 国产日韩欧美在线精品| 免费观看性生交大片5| 永久免费av网站大全| 午夜精品国产一区二区电影 | 亚洲欧美日韩卡通动漫| 在线观看av片永久免费下载| 色尼玛亚洲综合影院| 激情 狠狠 欧美| 久久久久久久久久人人人人人人| a级毛色黄片| 97超视频在线观看视频| 国产伦一二天堂av在线观看| 国内精品宾馆在线| 看非洲黑人一级黄片| 国产伦精品一区二区三区四那| 六月丁香七月| 亚洲精品,欧美精品| 国产日韩欧美在线精品| 女人十人毛片免费观看3o分钟| 纵有疾风起免费观看全集完整版 | 国产一区二区三区av在线| 精品久久久噜噜| 亚洲aⅴ乱码一区二区在线播放| 精品人妻偷拍中文字幕| 精品久久久久久久久久久久久| 18禁在线播放成人免费| 亚洲av日韩在线播放| 亚洲人成网站在线播| 国产一区二区亚洲精品在线观看| 99久久九九国产精品国产免费| 狂野欧美白嫩少妇大欣赏| 日日撸夜夜添| 伊人久久精品亚洲午夜| 久久久久精品性色| 菩萨蛮人人尽说江南好唐韦庄| 久久久久免费精品人妻一区二区| 日本欧美国产在线视频| 美女主播在线视频| 亚洲人成网站在线观看播放| 久久久久九九精品影院| 日本午夜av视频| 成人午夜精彩视频在线观看| 精品久久久精品久久久| 亚洲欧美精品专区久久| 国产探花在线观看一区二区| 六月丁香七月| 国产午夜精品久久久久久一区二区三区| 十八禁国产超污无遮挡网站| 亚洲欧美日韩卡通动漫| 国产黄色免费在线视频| 国产精品麻豆人妻色哟哟久久 | 国产精品一区二区性色av| 欧美极品一区二区三区四区| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久电影| 18禁动态无遮挡网站| 精品一区在线观看国产| 成年女人看的毛片在线观看| 久久久久久久国产电影| 听说在线观看完整版免费高清| 精品熟女少妇av免费看| 熟妇人妻不卡中文字幕| 久久综合国产亚洲精品| 国产精品综合久久久久久久免费| av黄色大香蕉| 久久97久久精品| 亚洲精品久久久久久婷婷小说| 国产免费又黄又爽又色| 亚洲婷婷狠狠爱综合网|