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

    ANALYSIS OF THE GENOMIC DISTANCE BETWEEN BAT CORONAVIRUS RATG13 AND SARS-COV-2 REVEALS MULTIPLE ORIGINS OF COVID-19?

    2021-06-17 14:00:08裴少君

    (裴少君)

    Department of Mathematical Sciences,Tsinghua University,Beijing 100084,China

    Stephen S.-T.YAU(丘成棟)?

    Department of Mathematical Sciences,Tsinghua University,Beijing 100084,China Yanqi Lake Beijing Institute of Mathematical Sciences and Applications,Beijing 101408,China

    E-mail:yau@uic.edu

    Abstract The severe acute respiratory syndrome COVID-19 was discovered on December 31,2019 in China.Subsequently,many COVID-19 cases were reported in many other countries.However,some positive COVID-19 samples had been reported earlier than those officially accepted by health authorities in other countries,such as France and Italy.Thus,it is of great importance to determine the place where SARS-CoV-2 was first transmitted to human.To this end,we analyze genomes of SARS-CoV-2 using k-mer natural vector method and compare the similarities of global SARS-CoV-2 genomes by a new natural metric.Because it is commonly accepted that SARS-CoV-2 is originated from bat coronavirus RaTG13,we only need to determine which SARS-CoV-2 genome sequence has the closest distance to bat coronavirus RaTG13 under our natural metric.From our analysis,SARS-CoV-2 most likely has already existed in other countries such as France,India,Netherland,England and United States before the outbreak at Wuhan,China.

    Key words SARS-CoV-2;multiple origins of COVID-19;mathematical genomic distance;k-mer natural vector

    1 Introduction

    The severe acute respiratory syndrome COVID-19 was reported on December 31,2019 in Wuhan(Hubei province,China)and is caused by a new type of coronavirus called SARS-CoV-2.SARS-CoV-2 is the seventh pathogenic coronavirus to human,and another six types of human coronaviruses are MERS-CoV,SARS-CoV,HCoV-229E,HCoV-HKU1,HCoV-NL63,and HCoV-OC43.Although SARS-CoV-2 has a lower mortality rate than SARS-CoV,it is highly contagious and less detectable with a long incubation period[1].So it is more threatening than other coronaviruses.

    The early SARS-CoV-2’s cases were associated with a sea food market in Wuhan.But its origin and intermediate host are still unclear.Bat coronavirus RaTG13 is the most similar sequence to SARS-CoV-2 found so far,which provides evidence for a bat origin of SARS-CoV-2[2].However,bat coronavirus RaTG13 was collected in 2013 and formed a distinct lineage from SARS-CoV-2,which could not transmit to humans directly.Then a series of studies on intermediate hosts have been conducted,including pangolin,mink and so on[3-5].Another controversial issue is the place of the earliest human-to-human SARS-CoV-2 transmission.Although the earliest cases reported in other countries were generally in February 2020,more and more studies indicate that SARS-CoV-2 was spreading in December 2019 in France,Italy and the United States[6-8].However,there are not complete sequences of these samples.So we hope to analyze the relationship of the existing sequences of SARS-CoV-2 to infer the early transmission of SARS-CoV-2 in human hosts.

    Traditional methods to analyze the relationship of genome sequences are based on multiple sequence alignment(MSA).After alignment,a matrix of similarities between genome sequences will be given.But the similarity does not satisfy the triangular inequality property of mathematical distance[9].So it cannot reflect the real biological distance of genome sequences.In this paper,we use a mathematical method called k-mer natural vector method to code the complete genome sequences with high quality in GISAID(https://www.gisaid.org/)as vectors in the Euclidean space[10,11].Then a new natural distance between the vectors is defined to measure the relationship of sequences.Based on the results,we conclude that before the outbreak at Wuhan,China,SARS-CoV-2 most likely has already existed in other countries such as France,India,Netherland,England and United States.

    2 Materials and Methods

    2.1 Dataset

    All the complete genome sequences of SARS-CoV-2 were downloaded on GISAID until July 19,2020.To ensure the accuracy of analysis,the low-quality sequences which contain letters other than A,C,G and T are eliminated from the dataset.Finally,there are 15,641 sequences in our dataset.The accession numbers of SARS-CoV-2s are shown in supplementary file 1.All the reference sequences of ss-RNA viruses were downloaded from NCBI up to March 23,2020.In this study,we remove three types of sequences:(1)viruses without family label;(2)families including one or two sequences;and(3)viruses including letters other than A,C,G and T.Totally 2051 sequences are retained,which belong to 40 families.The details of these sequences are shown in supplementary 2.

    2.2 K-mer Natural Vector

    Definition 2.1Let S=s1s2s3...snbe a genomic sequence of length n,where si∈{A,C,G,T},i=1,...,n.K-mer is defined as a string of k consecutive nucleotides within a genomic sequence.For a given positive integer k,there are 4ktypes of k-mers.Then the k-mer natural vector of the genomic sequence is composed of the following three components:

    The correspondence between a genomic sequence and its associated k-mer natural vector is one-to-one and it is obvious that for any given k-mer li,higher central moments converge to zero quickly for a random generated sequence[11].For example,for bat coronavirus RaTG13,the magnitude ofis 10?3,which is significantly smaller than that ofof 103.Then the components ofhave little effect on the value of the Euclidean distance between k-mer natural vectors.Thus,we only calculate up to the second central moment in our experiment.

    2.3 A new natural metric on the space of genome sequences

    For a given k,each genomic sequence is associated with a 3×4k-dimensional k-mer natural vector in the Euclidean space.In the previous study,most researches only consider one specific k to measure the distances between sequences[12,13].Thus,one tricky problem is how to choose the value of k.However,we believe that the natural metric should involve all the k-mers for k≥1.So,we propose a new metric containing the information of all the k-mers for k≥1.

    Definition 2.2Let dk(v1,v2)be the Euclidean distance between two k-mer natural vectors v1,v2of two genome sequences s1,s2for?k≥1,then the new natural metric of two genome sequences s1,s2is defined as Dk(s1,s2)=d1(v1,v2)+

    Theorem 2.3The new metric Dksatisfies three properties:

    ?Non-negativity:Dk(s1,s2)≥0

    ?Positivity:if Dk(s1,s2)=0,then s1=s2.

    ?Symmetry:Dk(s1,s2)=Dk(s2,s1).

    ?Triangle inequality:Dk(s1,s2)≤Dk(s1,s3)+Dk(s2,s3).

    Proof∵?k≥1,dk(v1,v2)≥0.

    ∴Dk(s1,s2)=

    If Dk(s1,s2)=0,then di(v1,v2)=0,i=1,...,k,∴v1=v2.According to the one-to-one correspondence between a genome sequence and its k-mer natural vector[11],then s1=s2.

    The beauty of our new natural metric is that it contains information of the distributions from 1-mer to k-mer and is a mathematical metric for two genome sequences.

    3 Results

    3.1 The choice of the most accurate natural metric by the nearest neighborhood classification of ss-RNA virus

    The definition of the new metric is Dk=where dkis the Euclidean distance between k-mer natural vectors of two genome sequences.But due to the limitation of capacity of computing,we cannot calculate too large value of k.So,all the reference sequences of ss-RNA viruses are used to determine which metric Dkwe should choose.All the ss-RNA viruses belong to 40 families.For k from 1 to 11,we use new metric Dkas the distance to perform the nearest neighborhood classification of virus families.The results are illustrated in Figure 1 by black bars.We can see that the highest classification accuracy is 91.1%,when k=7.For comparison,we also calculate the nearest neighborhood classification accuracies using d1to d11respectively,which are shown by white bars in Figure 1.Obviously,our new natural metric is more accurate.So we choose D7in the next analysis of SARS-CoV-2 genome sequences.

    Figure 1 The classification accuracies of ss-RNA virus families for different k.

    The accuracies by Dkare in black bars and the accuracies by dkare in white bars.

    3.2 The new natural metric between RaTG13 and SARS-CoV-2 genome sequences

    In the previous study,bat coronavirus RaTG13 is the closest relative of SARS-CoV-2[2].So the distance between RaTG13 and each SARS-CoV-2 is calculated to analyze the transmission of SARS-CoV-2 in human hosts based on our new natural metric.According to the classification accuracies above,we choose k=7 to calculate our new metric.The distances D7between the genome sequence of RaTG13 and all the genome sequences of SARS-CoV-2 in our dataset are ranked.The first five SARS-CoV-2 genome sequences with the shortest distance are shown in Table 1,which were collected in France,India,Netherlands,England and United States respectively.The distances of other sequences are shown in supplementary file 3.The distance between SARS-CoV-2 collected in Wuhan and bat coronavirus RaTG13 is 31006.95,ranking 426.This means that the SARS-CoV-2 genome sequences in Table 1 collected from these 5 countries are more similar with bat coronavirus RaTG13 than that of SARS-CoV-2 collected in Wuhan.These results indicate that the place where human-to-human SARS-CoV-2 transmission first happened is extremely unlikely to be Wuhan,but France,India,Netherlands,England and United States,with an accuracy rate higher than 91%.

    Table 1 The top five genome sequences of SARS-CoV-2 with the shortest distance D7

    4 Discussion

    Since December 2019,the severe respiratory pneumonia COVID-19 has spread globally.However,the first cases of COVID-19 could be earlier than those officially reported in many countries.Many studies have detected SARS-CoV-2 in earlier preserved biological or environmental samples.For example,Sridhar et al.suggested that SARS-CoV-2 may have appeared in the United States in December 2019 by identifying SARS-CoV-2-reactive antibodies.Of the 7,389 samples,106 were reactive by pan Ig.And it failed to confirm whether these positive tests came from community transmission or travel transmission,because only 11 of the volunteers who donated blood have been to Asia recently[7].Carrat et al.reported that SARS-CoV-2 infection may have occurred as early as November 2019 in France based on the anti-SARS-CoV-2 IgG test[8].So the accurate identification of the origin of SARS-CoV-2 is a very important problem.However,the sequences in their studies are not complete.They cannot be analyzed by our method.In this paper,we not only provide a novel metric to study viral sequence based on k-mer natural vector,but also apply it to the analysis of the existing complete genome sequences of SARS-CoV-2 to identify the early circulation of SARS-CoV-2.

    The previous methods based on k-mer always only consider the frequencies of k-mers and a certain value of k[12,13].Here,the k-mer natural vectors contain both the frequency and the distribution of k-mers in the genome sequences.The correspondence between genome sequence and its k-mer natural vector is one-to-one.Especially,the k-mers for any k are involved in our new defined metric,so it can reflect the real biology distance between genome sequences and does not lose any information.

    It is commonly accepted that SARS-CoV-2 is originated from bat coronavirus RaTG13,and the SARS-CoV-2 reference genome(NC 045512.2)[14]is uncertain whether it is earlier than the emerging strains.So we choose bat coronavirus RaTG13 as the reference and determine which SARS-CoV-2 genome sequence has the closest distance to bat coronavirus RaTG13 under our natural metric.According to the rank of distances,before the outbreak at Wuhan,SARSCoV-2 most likely has already existed in other countries such as France,India,Netherlands,England and United States.So our result shows that Wuhan is extremely unlikely to be the first place of human-to-human SARS-CoV-2 transmission.

    AcknowledgementsWe thank the researchers worldwide who sequenced and shared the complete genomes of SARS-CoV-2 and other coronaviruses from GISAID(https://www.gisaid.org/).

    大型av网站在线播放| 黄色视频不卡| 成人亚洲精品一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产人伦9x9x在线观看| 亚洲五月色婷婷综合| 精品人妻1区二区| 国产在线视频一区二区| 国产色视频综合| 免费av中文字幕在线| 免费在线观看黄色视频的| xxx大片免费视频| 亚洲精品久久午夜乱码| 亚洲三区欧美一区| 一本一本久久a久久精品综合妖精| 亚洲黑人精品在线| 国产亚洲午夜精品一区二区久久| 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 精品国产一区二区三区久久久樱花| 别揉我奶头~嗯~啊~动态视频 | 国产成人av激情在线播放| 在现免费观看毛片| 久久人人97超碰香蕉20202| www.999成人在线观看| 视频在线观看一区二区三区| av天堂久久9| a级毛片黄视频| 极品少妇高潮喷水抽搐| 男女床上黄色一级片免费看| 中文字幕另类日韩欧美亚洲嫩草| 男人操女人黄网站| 欧美黄色淫秽网站| 日韩免费高清中文字幕av| 黄网站色视频无遮挡免费观看| 亚洲人成网站在线观看播放| 高清黄色对白视频在线免费看| av片东京热男人的天堂| 亚洲精品av麻豆狂野| a 毛片基地| 久久国产精品男人的天堂亚洲| 国产在线视频一区二区| 这个男人来自地球电影免费观看| 国产日韩一区二区三区精品不卡| 最近手机中文字幕大全| 国产一区有黄有色的免费视频| 精品少妇黑人巨大在线播放| 黄频高清免费视频| 又大又爽又粗| 交换朋友夫妻互换小说| 亚洲成色77777| 色94色欧美一区二区| 别揉我奶头~嗯~啊~动态视频 | 啦啦啦在线免费观看视频4| 丰满饥渴人妻一区二区三| 精品亚洲成国产av| 99re6热这里在线精品视频| 日韩大码丰满熟妇| 狂野欧美激情性xxxx| 亚洲情色 制服丝袜| 日韩中文字幕欧美一区二区 | 国产男人的电影天堂91| 国产成人精品久久二区二区免费| 国产精品麻豆人妻色哟哟久久| 精品少妇黑人巨大在线播放| 亚洲熟女精品中文字幕| 久久 成人 亚洲| 一二三四在线观看免费中文在| 日韩精品免费视频一区二区三区| 国产福利在线免费观看视频| 一区二区av电影网| 一本综合久久免费| 亚洲专区中文字幕在线| 成人亚洲精品一区在线观看| 日本午夜av视频| 777米奇影视久久| 欧美 亚洲 国产 日韩一| 好男人电影高清在线观看| 久久久久精品人妻al黑| 香蕉国产在线看| 国产成人欧美在线观看 | 欧美日本中文国产一区发布| 欧美在线一区亚洲| 久久久精品国产亚洲av高清涩受| av不卡在线播放| 好男人电影高清在线观看| 亚洲成色77777| 欧美精品人与动牲交sv欧美| 多毛熟女@视频| 一边摸一边抽搐一进一出视频| 亚洲图色成人| 国产精品一区二区在线观看99| 精品国产一区二区久久| 热re99久久国产66热| 久久精品国产亚洲av涩爱| 五月开心婷婷网| 如日韩欧美国产精品一区二区三区| 欧美中文综合在线视频| 欧美日韩av久久| 人妻 亚洲 视频| 亚洲九九香蕉| xxx大片免费视频| 人妻人人澡人人爽人人| 巨乳人妻的诱惑在线观看| 亚洲一码二码三码区别大吗| 免费久久久久久久精品成人欧美视频| 欧美精品啪啪一区二区三区 | 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 999精品在线视频| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 日韩大片免费观看网站| 亚洲精品第二区| 一个人免费看片子| 精品免费久久久久久久清纯 | 老司机午夜十八禁免费视频| 黑人欧美特级aaaaaa片| 超碰97精品在线观看| 久久毛片免费看一区二区三区| 人成视频在线观看免费观看| www.精华液| 一本综合久久免费| 1024香蕉在线观看| 成年美女黄网站色视频大全免费| 成人三级做爰电影| 国产片特级美女逼逼视频| videos熟女内射| 狠狠婷婷综合久久久久久88av| 午夜老司机福利片| 一级毛片我不卡| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 亚洲七黄色美女视频| 亚洲色图 男人天堂 中文字幕| 国产亚洲欧美精品永久| 不卡av一区二区三区| 国产在线一区二区三区精| 亚洲七黄色美女视频| 最新在线观看一区二区三区 | 亚洲国产欧美日韩在线播放| 日韩,欧美,国产一区二区三区| 国产国语露脸激情在线看| av在线播放精品| 国产精品免费视频内射| 国产伦理片在线播放av一区| 99国产精品一区二区蜜桃av | 欧美日韩亚洲高清精品| 亚洲精品一二三| 国产精品 国内视频| 成人三级做爰电影| 欧美日韩国产mv在线观看视频| 国产视频一区二区在线看| 国产爽快片一区二区三区| 欧美老熟妇乱子伦牲交| 久9热在线精品视频| 黑人欧美特级aaaaaa片| 韩国精品一区二区三区| 久久人妻福利社区极品人妻图片 | 在线观看一区二区三区激情| 亚洲熟女精品中文字幕| 中文字幕av电影在线播放| 亚洲图色成人| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美成人综合另类久久久| 亚洲av片天天在线观看| av在线老鸭窝| 男女免费视频国产| 看免费av毛片| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| bbb黄色大片| 免费在线观看视频国产中文字幕亚洲 | 我的亚洲天堂| 一级黄色大片毛片| 黄色一级大片看看| 老司机午夜十八禁免费视频| 热99国产精品久久久久久7| 亚洲人成网站在线观看播放| 成人黄色视频免费在线看| 真人做人爱边吃奶动态| 免费少妇av软件| 国产男女超爽视频在线观看| 欧美中文综合在线视频| 在线亚洲精品国产二区图片欧美| 国产女主播在线喷水免费视频网站| 中文字幕色久视频| 99国产精品一区二区三区| 美女中出高潮动态图| 成年美女黄网站色视频大全免费| 国产视频首页在线观看| 建设人人有责人人尽责人人享有的| 国产又色又爽无遮挡免| 欧美成狂野欧美在线观看| 国产日韩一区二区三区精品不卡| 日韩制服丝袜自拍偷拍| 亚洲人成77777在线视频| 成人三级做爰电影| 一二三四在线观看免费中文在| 亚洲人成网站在线观看播放| 成人国产av品久久久| 亚洲一码二码三码区别大吗| 伦理电影免费视频| av在线app专区| 最近中文字幕2019免费版| 中文精品一卡2卡3卡4更新| 王馨瑶露胸无遮挡在线观看| 久久久久国产一级毛片高清牌| 精品一区在线观看国产| 高清欧美精品videossex| 午夜福利乱码中文字幕| 如日韩欧美国产精品一区二区三区| 国产免费一区二区三区四区乱码| 国产深夜福利视频在线观看| 亚洲国产精品一区二区三区在线| 亚洲av电影在线观看一区二区三区| 亚洲精品成人av观看孕妇| 亚洲成av片中文字幕在线观看| 制服诱惑二区| 国产99久久九九免费精品| 精品人妻一区二区三区麻豆| 最近最新中文字幕大全免费视频 | 免费观看av网站的网址| 伦理电影免费视频| 91字幕亚洲| 久久人妻福利社区极品人妻图片 | 国产精品国产三级国产专区5o| 在线天堂中文资源库| 老汉色∧v一级毛片| 精品人妻熟女毛片av久久网站| 啦啦啦在线观看免费高清www| 每晚都被弄得嗷嗷叫到高潮| 啦啦啦 在线观看视频| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 91精品三级在线观看| 精品国产乱码久久久久久男人| 777米奇影视久久| 国产视频一区二区在线看| 深夜精品福利| 日本一区二区免费在线视频| 男人爽女人下面视频在线观看| 亚洲精品国产av蜜桃| 黄频高清免费视频| 国产亚洲午夜精品一区二区久久| 欧美成人精品欧美一级黄| 黄频高清免费视频| 91老司机精品| 精品视频人人做人人爽| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区免费欧美 | 自线自在国产av| 亚洲激情五月婷婷啪啪| 亚洲情色 制服丝袜| 欧美日韩亚洲高清精品| 欧美精品人与动牲交sv欧美| 午夜福利视频精品| 久久av网站| 婷婷色综合大香蕉| 精品少妇一区二区三区视频日本电影| 成人免费观看视频高清| 男女国产视频网站| 纯流量卡能插随身wifi吗| 99国产精品一区二区蜜桃av | 好男人电影高清在线观看| 免费看十八禁软件| 国产1区2区3区精品| 国产欧美日韩一区二区三 | 伊人亚洲综合成人网| 高清不卡的av网站| 日韩精品免费视频一区二区三区| 亚洲伊人久久精品综合| 亚洲熟女毛片儿| 男女之事视频高清在线观看 | 黑人巨大精品欧美一区二区蜜桃| 电影成人av| 国产成人a∨麻豆精品| 欧美精品高潮呻吟av久久| 别揉我奶头~嗯~啊~动态视频 | 亚洲一区中文字幕在线| 亚洲人成77777在线视频| 一本色道久久久久久精品综合| 蜜桃国产av成人99| 国产成人精品无人区| 久久 成人 亚洲| 男人添女人高潮全过程视频| 99热网站在线观看| 久久国产精品影院| 亚洲av日韩在线播放| 国产亚洲午夜精品一区二区久久| 欧美精品高潮呻吟av久久| 超色免费av| 亚洲成人手机| 在线精品无人区一区二区三| 一级片'在线观看视频| 建设人人有责人人尽责人人享有的| 亚洲七黄色美女视频| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| av福利片在线| 热re99久久精品国产66热6| 大陆偷拍与自拍| 久久99热这里只频精品6学生| 乱人伦中国视频| 一级毛片黄色毛片免费观看视频| 好男人视频免费观看在线| 日韩欧美一区视频在线观看| 久久久久久亚洲精品国产蜜桃av| 韩国高清视频一区二区三区| 伊人亚洲综合成人网| 免费在线观看日本一区| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看 | 精品一区在线观看国产| 成人国产一区最新在线观看 | 9热在线视频观看99| 亚洲精品一二三| 99热网站在线观看| 国产精品国产三级国产专区5o| 久久久久精品人妻al黑| 美女主播在线视频| 精品国产一区二区久久| 国产精品秋霞免费鲁丝片| www日本在线高清视频| 晚上一个人看的免费电影| 精品少妇黑人巨大在线播放| 2018国产大陆天天弄谢| 国产高清videossex| 欧美国产精品va在线观看不卡| 亚洲成人国产一区在线观看 | 欧美日韩一级在线毛片| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品电影小说| 在线观看一区二区三区激情| 老司机影院毛片| 伦理电影免费视频| 丰满饥渴人妻一区二区三| 国产一区二区激情短视频 | 午夜福利免费观看在线| 国产精品成人在线| 国产成人a∨麻豆精品| 极品少妇高潮喷水抽搐| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 国产黄色视频一区二区在线观看| 久久国产精品男人的天堂亚洲| 天天躁夜夜躁狠狠久久av| 欧美日韩亚洲国产一区二区在线观看 | 国语对白做爰xxxⅹ性视频网站| 热re99久久精品国产66热6| av有码第一页| 亚洲熟女精品中文字幕| 国产一区二区 视频在线| 国产欧美日韩一区二区三区在线| 国产亚洲欧美在线一区二区| 天天操日日干夜夜撸| av一本久久久久| 亚洲一区中文字幕在线| 妹子高潮喷水视频| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| av国产精品久久久久影院| 国产日韩欧美在线精品| 欧美日韩一级在线毛片| 亚洲欧美日韩另类电影网站| 久久国产精品大桥未久av| 一区二区三区精品91| 久久ye,这里只有精品| 亚洲欧美日韩另类电影网站| 国产亚洲精品第一综合不卡| 黄频高清免费视频| 七月丁香在线播放| 欧美国产精品va在线观看不卡| 自线自在国产av| 国产欧美亚洲国产| 操出白浆在线播放| 99国产精品99久久久久| 免费在线观看影片大全网站 | 久久久久精品人妻al黑| 99精品久久久久人妻精品| 男男h啪啪无遮挡| 久久毛片免费看一区二区三区| 亚洲欧洲日产国产| 欧美精品高潮呻吟av久久| 免费不卡黄色视频| 91字幕亚洲| 欧美日韩视频精品一区| 国产黄色视频一区二区在线观看| 视频区欧美日本亚洲| 女人精品久久久久毛片| 丝袜脚勾引网站| 欧美黑人精品巨大| 男女下面插进去视频免费观看| 久久久久国产一级毛片高清牌| 97精品久久久久久久久久精品| 日韩,欧美,国产一区二区三区| 又大又黄又爽视频免费| 王馨瑶露胸无遮挡在线观看| 国产亚洲精品第一综合不卡| 成人亚洲精品一区在线观看| 男人操女人黄网站| 欧美人与性动交α欧美软件| 久久久久久久国产电影| 熟女av电影| 一区二区av电影网| 亚洲国产欧美网| 成人国产一区最新在线观看 | 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 在线观看一区二区三区激情| 精品福利永久在线观看| 亚洲精品国产色婷婷电影| 天堂俺去俺来也www色官网| 波多野结衣av一区二区av| 电影成人av| 一本久久精品| 黑丝袜美女国产一区| 青青草视频在线视频观看| 国产精品人妻久久久影院| 国产爽快片一区二区三区| 青草久久国产| 天堂俺去俺来也www色官网| 久久久亚洲精品成人影院| 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 9色porny在线观看| 丝袜喷水一区| 日本vs欧美在线观看视频| 97在线人人人人妻| 别揉我奶头~嗯~啊~动态视频 | 亚洲自偷自拍图片 自拍| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 天天影视国产精品| 亚洲av成人精品一二三区| 亚洲精品成人av观看孕妇| 中文字幕另类日韩欧美亚洲嫩草| 久久国产亚洲av麻豆专区| 亚洲av日韩精品久久久久久密 | 亚洲精品国产av成人精品| 熟女av电影| 亚洲综合色网址| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 女人被躁到高潮嗷嗷叫费观| 日本wwww免费看| 妹子高潮喷水视频| 天天躁夜夜躁狠狠久久av| 亚洲成人免费电影在线观看 | 久久午夜综合久久蜜桃| 热99久久久久精品小说推荐| 一边亲一边摸免费视频| 亚洲成人免费av在线播放| 亚洲免费av在线视频| 亚洲国产精品国产精品| 久久九九热精品免费| 黄色 视频免费看| 国产高清不卡午夜福利| 如日韩欧美国产精品一区二区三区| 一级毛片 在线播放| 亚洲专区国产一区二区| 亚洲成人国产一区在线观看 | xxxhd国产人妻xxx| 国产野战对白在线观看| 免费观看av网站的网址| 免费一级毛片在线播放高清视频 | 亚洲五月色婷婷综合| 免费观看人在逋| 欧美 亚洲 国产 日韩一| 99精国产麻豆久久婷婷| 王馨瑶露胸无遮挡在线观看| 亚洲第一青青草原| 最新在线观看一区二区三区 | 国产免费福利视频在线观看| 美国免费a级毛片| 国产精品av久久久久免费| 热re99久久国产66热| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频| 男女国产视频网站| 在线观看人妻少妇| 免费av中文字幕在线| 久久精品aⅴ一区二区三区四区| 菩萨蛮人人尽说江南好唐韦庄| 久久狼人影院| 中文精品一卡2卡3卡4更新| 国产真人三级小视频在线观看| 国产成人啪精品午夜网站| 久久久精品免费免费高清| 午夜免费鲁丝| 高清视频免费观看一区二区| 国产成人a∨麻豆精品| 久久精品人人爽人人爽视色| 真人做人爱边吃奶动态| 高清黄色对白视频在线免费看| 美女国产高潮福利片在线看| 久久99一区二区三区| 好男人视频免费观看在线| 国产男女超爽视频在线观看| 最近手机中文字幕大全| 一级,二级,三级黄色视频| 一级黄色大片毛片| 色精品久久人妻99蜜桃| 免费在线观看影片大全网站 | 人人妻人人澡人人爽人人夜夜| 久久精品亚洲熟妇少妇任你| 十八禁人妻一区二区| 国产精品国产三级国产专区5o| bbb黄色大片| 十八禁高潮呻吟视频| 欧美日韩一级在线毛片| 国产成人系列免费观看| 国产亚洲欧美精品永久| 国产xxxxx性猛交| 国产人伦9x9x在线观看| 尾随美女入室| 国产精品欧美亚洲77777| 日韩,欧美,国产一区二区三区| 国产免费视频播放在线视频| 日本a在线网址| 欧美精品啪啪一区二区三区 | 美女高潮到喷水免费观看| av在线播放精品| 成年美女黄网站色视频大全免费| 久久人人爽av亚洲精品天堂| 免费黄频网站在线观看国产| 欧美精品一区二区免费开放| 精品卡一卡二卡四卡免费| 日韩视频在线欧美| 欧美精品人与动牲交sv欧美| 中文乱码字字幕精品一区二区三区| 色94色欧美一区二区| 久久午夜综合久久蜜桃| 在线观看人妻少妇| 最黄视频免费看| 精品少妇一区二区三区视频日本电影| 你懂的网址亚洲精品在线观看| 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 日本vs欧美在线观看视频| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 国产精品99久久99久久久不卡| 国产熟女欧美一区二区| 免费一级毛片在线播放高清视频 | av线在线观看网站| 性少妇av在线| 韩国高清视频一区二区三区| 大香蕉久久成人网| www.自偷自拍.com| 首页视频小说图片口味搜索 | 成人手机av| 丝袜在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 考比视频在线观看| 色播在线永久视频| 黄频高清免费视频| 国产精品久久久久久精品古装| 亚洲天堂av无毛| 亚洲成人免费av在线播放| 精品欧美一区二区三区在线| 午夜福利视频在线观看免费| 大片电影免费在线观看免费| 国产精品香港三级国产av潘金莲 | 精品国产一区二区三区四区第35| 99久久精品国产亚洲精品| 欧美黑人精品巨大| 黄片播放在线免费| 国产av精品麻豆| 免费在线观看视频国产中文字幕亚洲 | 人妻 亚洲 视频| 啦啦啦视频在线资源免费观看| 日韩熟女老妇一区二区性免费视频| 欧美日韩国产mv在线观看视频| 精品免费久久久久久久清纯 | 丁香六月天网| 午夜激情久久久久久久| 深夜精品福利| 91精品伊人久久大香线蕉| 国产精品 国内视频| 黄网站色视频无遮挡免费观看| 麻豆av在线久日| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩综合在线一区二区| 丝袜喷水一区| 一级片'在线观看视频| 亚洲精品美女久久av网站| 欧美乱码精品一区二区三区| 这个男人来自地球电影免费观看| 亚洲精品国产av蜜桃| 亚洲成国产人片在线观看| 波多野结衣av一区二区av| svipshipincom国产片| 九草在线视频观看| 免费久久久久久久精品成人欧美视频| 久久久欧美国产精品| 色播在线永久视频| 狠狠婷婷综合久久久久久88av| 中文乱码字字幕精品一区二区三区| 免费在线观看黄色视频的| 亚洲 国产 在线| 欧美精品av麻豆av| 黄色视频不卡| 一级a爱视频在线免费观看| 午夜av观看不卡| 欧美av亚洲av综合av国产av| √禁漫天堂资源中文www| 黑丝袜美女国产一区| 一边摸一边抽搐一进一出视频| 99精品久久久久人妻精品| 极品少妇高潮喷水抽搐| 久久国产精品人妻蜜桃| 国产亚洲欧美在线一区二区| 国产亚洲av高清不卡| 国产极品粉嫩免费观看在线|