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

    野生與養(yǎng)殖厚殼貽貝體框特征和線粒體16S rRNA 基因的比較分析

    2023-10-24 03:18:58夏飛宇張秀梅王一航
    南方水產科學 2023年5期
    關鍵詞:特征生長差異

    夏飛宇,張秀梅,徐 翩,徐 穎,王一航

    浙江海洋大學 水產學院,浙江 舟山 316022

    厚殼貽貝 (Mytiluscoruscus) 是廣泛分布于黃海和東海沿岸海域的重要經濟貝類。浙江省嵊泗縣是我國厚殼貽貝的原產地和主要養(yǎng)殖區(qū)域[1],20 世紀70 年代前,嵊泗列島產出的貽貝一直為野生厚殼貽貝[2],它們通常生活在低潮帶至水深20 m 處,附著于海浪擊打的巖礁上和溝縫中。隨著厚殼貽貝人工育苗技術研究的開展和推廣[1,3],其養(yǎng)殖業(yè)在嵊泗縣逐漸發(fā)展起來,使用筏或延繩的懸浮養(yǎng)殖技術迅速普及。較其他養(yǎng)殖方式(如木樁養(yǎng)殖和底播養(yǎng)殖),懸浮養(yǎng)殖下的貽貝生長速度更快,軟組織飽滿程度更高[4]。

    潮間帶是陸地與海洋的過渡地帶,嵊泗縣周邊潮間帶以巖相為主,與養(yǎng)殖海域相比,這是種截然不同的生境。對巖相潮間帶的大量研究表明,海浪拍擊、巖縫擠壓、空氣暴露時長以及生長密度會對軟體動物的外殼產生深遠影響[5-6];這些復雜的因素使得無法直接預測其形態(tài)變異。因此,了解厚殼貽貝的生活環(huán)境是如何影響其生長的,對貽貝養(yǎng)殖業(yè)的高質量發(fā)展極為重要。

    線粒體DNA (mtDNA)由于具有突變率高、母系遺傳幾乎不發(fā)生重組等特點,是研究物種遺傳多樣性和分子進化的有效工具。對于雙殼貝類形態(tài)學的研究,國內外已有較多報道,多與遺傳多樣性結合(通過微衛(wèi)星標記和DNA 條形碼)[7-9]。如Yuan 等[10]利用形態(tài)學、微衛(wèi)星和mtDNA 標記對日本野生和中國養(yǎng)殖群體櫛孔扇貝 (Chlamys farreri) 進行研究,發(fā)現它們之間無顯著的形態(tài)差異,但日本野生扇貝可分為南北兩個種群且和中國養(yǎng)殖群體存在遺傳差異;He 等[11]通過全基因組測序和殼形分析證實太平洋牡蠣 (Crassostreagigas) 的外殼形狀會在環(huán)境應激反應中發(fā)生變化;徐義平等[12]結合形態(tài)學、SRAP 標記兩種手段,探究了中國沿海4 個省份的11 個縊蟶 (Sinonovaculaconstricta) 群體,指出浙江樂清灣野生群體間具有較高的形態(tài)和遺傳相似性。目前關于野生和養(yǎng)殖厚殼貽貝形態(tài)差異的研究尚未見報道,尚不清楚這種差異是否由遺傳造成。基于此,本研究采集了嵊泗縣嵊山、枸杞島巖相潮間帶 (簡稱野生) 和筏架養(yǎng)殖區(qū) (簡稱養(yǎng)殖) 不同生境下的厚殼貽貝樣本,測量分析其殼長、殼質量,構建體框特征并進行形態(tài)特征分析,同時對比其線粒體16SrRNA基因序列,以期為厚殼貽貝的種質資源保護及其養(yǎng)殖產業(yè)發(fā)展提供參考。

    1 材料與方法

    1.1 調查區(qū)域及樣品采集

    厚殼貽貝于2021 年8 月至2022 年4 月采自浙江省嵊泗縣嵊山鎮(zhèn)周邊海域 (122.73°E—122.79°E,30.69°N—30.75°N)。其中養(yǎng)殖貽貝取自枸杞島干斜村金盟厚殼貽貝筏架養(yǎng)殖區(qū);野生樣品為枸杞島巖相潮間帶自然生長的貽貝。根據殼長確定厚殼貽貝年齡[13]:1 齡殼長介于24.0~56.0 mm,2 齡殼長介于56.0~75.0 mm,3 齡殼長介于75.0~85.0 mm,4 齡貝殼長介于85.0~100.0 mm,5 齡殼長介于100.0~110.0 mm。為避免異速生長帶來的影響,體框特征分析選取的樣品為2+ 齡。

    1.2 樣品處理

    采集貝體完整無損的厚殼貽貝,刮干凈殼表面附著的生物及足絲,帶回實驗室解剖。烘干后用電子天平稱取殼干質量 (精度0.01 g)。從中隨機挑選157 只厚殼貽貝 (養(yǎng)殖78 只,野生79 只),參考白曉倩等[14]的方法并略做改進,構建體框特征模型。體框特征點如圖1 所示,包括殼寬(W)、殼高 (H)、殼長頂 (A)、殼背頂 (B)、鉸合部上端(C)、殼頂 (D)、寬頂 (O)、足絲孔 (E) 及殼頂至腹緣的垂足(F)。通過電子數顯游標卡尺 (精度0.01 mm) 測量它們的長度和距離,以此獲得W (L1)、H (L2)、AB (L3)、BC (L4)、CD (L5)、DE (L6)、AE (L7)、OA (L8)、OB (L9)、OC (L10)、OD (L11)、OE (L12)、OF (L13) 13 個貝體框特征指標。

    圖1 體框特征測量點Fig.1 Measurement points of shell frame characteristics

    1.3 線粒體DNA 提取、PCR 擴增及測序

    隨機額外選擇野生和養(yǎng)殖厚殼貽貝各20 只,解剖并提取其閉殼肌mtDNA。步驟如下:取約100 mg 閉殼肌組織,采用標準的酚-氯仿法提取基因組DNA,將乙醇沉淀后的基因組DNA 加入100 μL 胰蛋白酶/EDTA (TE) 緩沖液溶解。擴增16SrRNA基因片段的引物序列如下:16S ar:5'-CGCCTGTTTATCAAAAACAT-3';16S br:5'-CCGGTCTGAACTCAGATCACGT-3'。擴增產物純化及測序(Sanger) 由上海生工生物公司完成。

    1.4 數據處理與分析

    1.4.1 形態(tài)差異分析

    使用SPSS 26.0 軟件分別計算野生和養(yǎng)殖厚殼貽貝的體框特征指標平均值、標準差以及變異系數,不同生長區(qū)域的貽貝通過獨立樣本t檢驗判斷顯著性 (P<0.05 表示差異顯著)。結果由GraphPad Prism 9.0 軟件繪圖。此外,使用R-4.1.3 軟件對13 個體框特征參數進行Pearson 相關性分析;通過主成分分析 (Principal component analysis,PCA) 對數據進行降維,保留特征根大于1 的主成分 (Principal component,PC) 用于后續(xù)分析;運用逐步判別法,構建野生和養(yǎng)殖厚殼貽貝的判別函數并計算判別準確率,計算公式參照李思發(fā)等[15]。

    1.4.2 遺傳多樣性分析

    將測序結果用MEGA 7.0 剪切、比對并輔以人工校對分析核苷酸變異,以Kimura 雙參數法 (Kimura-2-parameter)為替代模型并基于鄰接距離法 (Neighbor-Joining) 構建野生和養(yǎng)殖厚殼貽貝的系統發(fā)育樹[16](重抽樣3 000 次)。使用DNASP 5.10 計算單倍型數量 (N)、單倍型多樣性 (Hd) 和核苷酸多樣性 (π) 來估計野生和養(yǎng)殖種群的遺傳多樣性。由POPART 繪制median joining network (MJN),構建養(yǎng)殖和野生厚殼貽貝的單倍型網絡圖。

    2 結果

    2.1 野生與養(yǎng)殖厚殼貽貝殼長與殼質量的關系

    對已收集到的所有厚殼貽貝殼長 (L)、殼干質量 (G) 數據 (養(yǎng)殖316 只、野生393 只) 進行回歸分析,得到兩條回歸關系式:G野= 0.000 3L2.5937(R2=0.930),G養(yǎng)=0.000 2L2.5537(R2=0.877)。從回歸式和回歸曲線 (圖2) 可以看出,厚殼貽貝的殼長-殼質量關系式基本符合冪函數增長規(guī)律。野生貽貝殼質量與殼長的冪指數b(2.593 7) 要略大于養(yǎng)殖貽貝(b=2.553 7),意味著在相同體長下野生群體的殼質量更大,且這種趨勢隨著貽貝個體的增長更為明顯。

    圖2 野生與養(yǎng)殖厚殼貽貝殼長和殼質量的關系Fig.2 Relationship between shell length and shell mass of wild and cultured M.coruscus

    2.2 野生與養(yǎng)殖厚殼貽貝體框特征參數比較

    不同生長區(qū)域厚殼貽貝的貝殼體框特征參數見表1。從圖3 可以看出,除L4、L11和L13外,野生群體和養(yǎng)殖群體的各項體框特征差異顯著 (P<0.01)。與自然生長在潮間帶的野生貽貝相比,養(yǎng)殖群體各項體框特征長度均相對較短,僅L7大于野生群體。統計結果顯示,野生群體L4的變異系數最大 (21.22%),遠高于養(yǎng)殖群體及其他特征指標,說明野生群體貝殼在L4這一指標上的選擇潛力很大。13 項體框指標中,野生群體有9 項指標 (L1、L3、L4、L5、L6、L7、L9、L11、L13) 的變異系數大于養(yǎng)殖群體,說明野生群體的貝殼在生長過程中的可塑性要略高于養(yǎng)殖群體。上述結果表明,不同生長區(qū)域的厚殼貽貝外觀上的區(qū)分度較大。

    表1 野生與養(yǎng)殖厚殼貽貝體框特征參數Table 1 Shell frame characteristics parameters of wild and cultured M.coruscus

    圖3 野生與養(yǎng)殖厚殼貽貝體框特征差異Fig.3 Differences in shell frame characteristics between wild and cultured M.coruscus

    2.3 體框特征之間的相關性分析和主成分分析

    由圖4 可知,野生與養(yǎng)殖貽貝各體框特征間存在相關性 (P<0.05),但所有體框特征間的相關系數絕對值均未大于0.7。L1和L2作為衡量厚殼貽貝外部形態(tài)的重要指標,與其他變量的正相關最多。而L7與其他變量均為負相關,且與L6的負相關系數最大 (-0.65)。各體框特征中,L11與其他特征的相關程度最低,僅與L8、L9顯著相關。L9與L2、L8存在最高的正相關性 (相關系數均為0.63)。

    圖4 體框特征間的相關性系數Fig.4 Correlation coefficient between shell frame characteristics

    通過對體框特征進行KMO 適合度和Bartlett 球形檢驗,發(fā)現各特征值的相關系數矩陣和零矩陣存在顯著性差異 (Bartlett 球形檢驗P<0.01),且主成分分析的適合度為適中 (KMO=0.742>0.7),因此對13 個體框特征進一步做主成分分析。提取的4 個特征根大于1 的主成分中,PC1 解釋的方差貢獻率為32.44% (表2),遠高于另外3 個主成分 (分別為14.02%、11.33%和9.70%)。PC1 由L1、L2、L3、L6、L7、L8、L9、L10、L12決定 (載荷絕對值大于0.5),第二個主成分主要由L5(載荷絕對值等于0.633) 決定 (表2)。L11主要與PC3 相關,而PC4 主要受L4影響。由PC1 和PC2 的得分圖 (圖5) 可見,養(yǎng)殖群體的數據點較為集中,說明其形態(tài)一致性較高;而野生群體的數據點分散,說明其形態(tài)更為多樣。

    表2 野生與養(yǎng)殖厚殼貽貝體框特征之間的主成分分析Table 2 Principal component analysis of shell frame characteristics between wild and cultured M.coruscus

    圖5 野生與養(yǎng)殖厚殼貽貝主成分分析得分圖Fig.5 Principal component scores of wild and cultured M.coruscus

    2.4 體框特征之間的判別分析

    采用逐步判別分析法,從13 個體框特征中選擇對判別結果有顯著影響的。F檢驗結果顯示L1、L2、L10、L12這4 個體框特征對判別結果有顯著性影響 (P<0.01)。以此得到一個典型判別函數LD1=30.60L1+15.12L2+14.34L10+6.99L12-285.02,典型判別函數相關性為0.848,Wilks'lambda檢驗顯著性<0.01,說明該判別函數具有顯著性。

    使用Fisher's 判別法,得到各生長區(qū)域厚殼貽貝的分類函數如下:

    養(yǎng)殖群體:y養(yǎng)=573.15L1+448.99L2+437.26L10+132.74L12-285.02

    厚殼群體:y野=670.28L1+497.00L2+482.79L10+154.93L12-364.06

    依靠判別函數,對養(yǎng)殖貽貝的判別準確率達到了98.7%,而對潮間帶野生貽貝的判別準確率也有91.1%,綜合判別準確率為94.9%??梢娕袆e公式能有效區(qū)分養(yǎng)殖和野生群體。

    2.5 野生與養(yǎng)殖群體16S rRNA 基因差異

    對擴增和測序成功的24 只厚殼貽貝序列進行比對后,獲得507 bp 的16SrRNA基因片段。樣品的平均核苷酸組成為A=30.8%,C =15.9%,G=23.3%,T=30.0%,G+C 平均含量39.2%。共檢測出12 個單倍型,其中養(yǎng)殖群體獨有單倍型4 個,野生群體獨有單倍型6 個,兩群體共享單倍型2 個。構建的單倍型網絡圖 (圖6) 顯示,厚殼貽貝的單倍型明顯分為了兩簇 (分別以Hap-1 和Hap-2 為中心),各簇之內的單倍型堿基變異步數均較低,但兩簇之間的堿基變異步數達到了101。本次樣品未檢測到Hap-7 和Hap-8間的連接節(jié)點單倍型。由系統發(fā)育樹 (圖7) 可知,每個分支均包含來自野生和養(yǎng)殖的厚殼貽貝,表明野生和養(yǎng)殖厚殼貽貝基本無分子水平上的遺傳分化。

    圖6 野生與養(yǎng)殖厚殼貽貝16S rRNA基因的單倍型網絡圖Fig.6 Haploid network diagram of 16S rRNA gene in wild and cultured M.coruscus

    圖7 鄰接距離法構建的厚殼貽貝16S rRNA 系統發(fā)育樹注:YS 代表野生厚殼貽貝,YZ 代表養(yǎng)殖厚殼貽貝。Fig.7 Neighbour-joining tree based on 16S rRNA gene of M.coruscusNote: YS represents wild M.coruscus,and YZ represents cultured M.coruscus.

    由DNASP 分析可見,養(yǎng)殖貽貝的Hd為0.682,明顯低于野生貽貝 (0.894)。在兩個群體中一共檢測到112 個多態(tài)位點,其中簡約信息位點在野生和養(yǎng)殖群體中的數量均達到105。兩個群體的π較為接近 (養(yǎng)殖群體為0.087,野生群體0.089),展現出較高的遺傳多樣性。

    3 討論

    3.1 野生和養(yǎng)殖厚殼貽貝的形態(tài)結構比較

    殼形既是雙殼貝類選育的重要指標,也是實際銷售過程中消費者推斷其經濟價值的重要判斷依據。本研究中,L1、L2代表的變量分別是殼寬、殼高,其大小直接影響殼的整體外形。相同體長下野生厚殼貽貝的殼更寬、更厚(P<0.01),意味著野生群體形態(tài)可能更圓鈍,而養(yǎng)殖群體相對更扁平。這與Grefsrud 和Oivind[17]對巨海扇蛤 (Pectenmaximus) 的研究一致,他們發(fā)現與同齡野生巨海扇蛤相比,養(yǎng)殖海扇蛤的殼更脆、更小、更輕,在應對潮汐和海浪時,較大的殼質量及尺寸可以起到更有效的保護作用。L10、L12是通過逐步判別篩選出的另兩個重要變量,分別表示殼背頂距離足絲孔和鉸合部上端的距離。結合L1、L7、L8與L11的值,可以看出野生厚殼貽貝殼寬頂O 的位置更靠近前端,它將殼面分為上 (四邊形AEOB)、下 (四邊形EOCD) 兩部分,養(yǎng)殖群體有寬大扁平的上部外殼,而野生群體則有較為膨脹的下部外殼。通常來說,扁平的殼形能減小水流阻力,而厚重圓潤的殼形使野生厚殼貽貝更能抵抗海浪的拍擊。顧志峰等[18]在對馬氏珠母貝 (Pinctada martensi) 的研究中也發(fā)現大亞灣野生種群的貝體外形較凸,而印度養(yǎng)殖種群的貝體外形顯得扁平。L4的大小代表了貽貝背角彎曲的程度,L4越大,貽貝的背角越不明顯,背緣呈弧形三角形。L4與其他指標的相關性相對較低且有著極高的變異系數,說明在后天的生長發(fā)育過程中,背角可能更易受外界環(huán)境干擾。已有研究也報道了厚殼貽貝幼年個體背角很明顯,但隨著生長會產生不同的變化[19]。

    3.2 野生和養(yǎng)殖厚殼貽貝的遺傳多樣分析

    測序結果表明厚殼貽貝出現明顯的堿基A/T 傾向性(A/T 占比達到60%),這與軟體動物線粒體基因堿基組成特點相一致[20-21]。遺傳多樣性是生物多樣性的重要組成部分,其高低決定了物種的生存進化能力。本研究中野生與養(yǎng)殖厚殼貽貝呈現出高核苷酸多樣性和單倍型多樣性,其16SrRNA核苷酸多樣性均高于0.08,大于四角蛤蜊 (Mactraveneriformis)、西施蛤 (Coelomactraantiquate)[22]、菲律賓蛤仔 (Ruditapesphilippinarum)[23],說明嵊泗厚殼貽貝種質資源保護較好,且野生群體的遺傳多樣性要高于養(yǎng)殖群體。系統發(fā)育結果顯示所有厚殼貽貝形成了兩個分支,但野生和養(yǎng)殖厚殼貽貝在每個分支均有出現,可能是因為它們之間存在基因流,說明單從遺傳差異無法完全解釋野生和養(yǎng)殖厚殼貽貝的形態(tài)差異。管成林等[8]、Shen 等[24]在對線粒體COI基因的研究中同樣發(fā)現野生和養(yǎng)殖厚殼貽貝無明顯的遺傳分化。嵊泗養(yǎng)殖海域厚殼貽貝苗種部分引自福建,葉瑩瑩等[25]對福州、舟山、溫州、寧德4 個地理群體厚殼貽貝的16SrRNA序列分析發(fā)現,福州群體與其他群體間遺傳距離較大,這可能是造成本研究中單倍型兩簇間出現較大堿基差異的原因。

    3.3 野生與養(yǎng)殖厚殼貽貝形態(tài)差異的解析

    雙殼貝類的殼表現出一種極強的表型可塑性[26],這種可塑性可能是種內遺傳變異所致[27-28],也可以歸因于某些環(huán)境因素,這是對棲息環(huán)境長期適應的結果[29-31]。本研究中,即使是距潮間帶僅幾百米的養(yǎng)殖厚殼貽貝,殼的形態(tài)也與野生群體有顯著性差異,這可能是由波浪暴露影響所致。固著生物在生命周期會產生位移風險,這種風險隨著流速的增加而上升,隨著貽貝韌性或附著強度的增加而降低[32]。很多學者就海洋波浪作用對貽貝殼形態(tài)的影響做了研究,發(fā)現波浪暴露區(qū)域的貽貝有著和非波浪暴露區(qū)貽貝不同的殼分泌模式 (形狀和厚度)。例如Steffani 和Branch[33]發(fā)現紫貽貝 (M.galloprovincialis) 在受到波浪中度暴露時的生長速度及殼厚度,大于無暴露和高強度波浪暴露時。Akester和Martel[34]發(fā)現暴露在潮間帶中的貽貝 (M.trossulus) 的殼高/寬比較低,還有更堅固的鉸合齒和更厚的韌帶。因此,波浪暴露可能是影響貽貝殼形態(tài)的重要環(huán)境因素。

    造成野生與養(yǎng)殖貽貝形態(tài)差異的另一個原因可能是生長密度導致的營養(yǎng)條件差異。雖然已有研究表明過高的生長密度不會使貽貝受到物理上的直接擠壓、磨損[35],但大面積的養(yǎng)殖活動直接影響了海區(qū)的水動力結構[36]。溫英[37]通過對流場的模擬證實了枸杞島貽貝養(yǎng)殖區(qū)的水層流速會受到貽貝串阻力的衰減作用。由于貽貝是濾食性生物,水流交換強度會直接影響其棲息地的餌料質量和數量。相較于流速平緩的養(yǎng)殖區(qū),潮間帶的營養(yǎng)物質擴散更快,沿岸的波浪可以快速彌補貽貝附著地的顆粒有機物消耗[38]。這種營養(yǎng)條件差異會影響貝類的外殼形態(tài),如Watanabe 和Katayama[39]發(fā)現棲息于營養(yǎng)條件較好環(huán)境中的菲律賓蛤仔個體比生長于貧營養(yǎng)環(huán)境中的個體殼寬更小,殼形橢圓化程度更高。由于雙殼貝類的殼形與其出肉率關系密切[40-41],在波浪作用頻繁的潮間帶,食物供應的增加可能使厚殼貽貝生長更快,進而可以分配更多的能量用于軟組織發(fā)育。但與此同時,海水動力也可能增加厚殼貽貝對足絲附著強度和貝殼厚度的需求,且貽貝的攝食率在食物濃度達到一定閾值時會保持相對穩(wěn)定,新陳代謝加快,消耗的能量也越來越多[42]。因此推測在一定閾值內,波浪暴露和營養(yǎng)條件會對貽貝生長產生積極影響,但過度的波浪作用和過快的水流速度則可能產生負面影響。

    4 結論

    綜上所述,枸杞島海域厚殼貽貝種質資源保護較好,但不同的生長環(huán)境使浮筏養(yǎng)殖厚殼貽貝較同齡野生群體的貝殼更為輕薄,16SrRNA基因測序結果證明了這種形態(tài)差異不是由遺傳造成。本研究所選取的13 個體框特征指標中有11 個存在顯著性差異,關鍵判別變量為L1(殼寬)、L2(殼高)、L10(殼頂至鉸合部上端距離) 和L12(殼頂至足絲孔距離),這4 個體框特征可作為評價厚殼貽貝整體形狀、生長潛力和營養(yǎng)狀況的重要指標。生長密度和水動力差異可能是造成厚殼貽貝野生和養(yǎng)殖群體形態(tài)差異的主要原因,今后可從這兩個角度出發(fā),合理規(guī)劃枸杞島貽貝筏架區(qū)的布局。

    猜你喜歡
    特征生長差異
    相似與差異
    音樂探索(2022年2期)2022-05-30 21:01:37
    碗蓮生長記
    小讀者(2021年2期)2021-03-29 05:03:48
    如何表達“特征”
    生長在哪里的啟示
    華人時刊(2019年13期)2019-11-17 14:59:54
    找句子差異
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    生長
    文苑(2018年22期)2018-11-19 02:54:14
    生物為什么會有差異?
    抓住特征巧觀察
    《生長在春天》
    亚洲美女黄片视频| av天堂久久9| 在线观看免费日韩欧美大片| 亚洲精品自拍成人| 欧美乱码精品一区二区三区| 窝窝影院91人妻| 精品少妇久久久久久888优播| 最新美女视频免费是黄的| √禁漫天堂资源中文www| 午夜老司机福利片| 久久这里只有精品19| 国产欧美日韩一区二区精品| 亚洲精品中文字幕一二三四区 | 俄罗斯特黄特色一大片| 日韩三级视频一区二区三区| 午夜福利免费观看在线| 日韩欧美一区二区三区在线观看 | 最近最新中文字幕大全电影3 | 亚洲精品av麻豆狂野| 国产高清国产精品国产三级| 50天的宝宝边吃奶边哭怎么回事| 国产在线精品亚洲第一网站| 丰满迷人的少妇在线观看| 欧美精品av麻豆av| 正在播放国产对白刺激| 午夜久久久在线观看| bbb黄色大片| 丝袜人妻中文字幕| 日韩视频一区二区在线观看| a级片在线免费高清观看视频| 国产日韩欧美在线精品| www.自偷自拍.com| 一级a爱视频在线免费观看| 下体分泌物呈黄色| 精品国内亚洲2022精品成人 | cao死你这个sao货| 十八禁网站网址无遮挡| 无遮挡黄片免费观看| 日韩中文字幕欧美一区二区| 91精品国产国语对白视频| 老司机福利观看| 久久香蕉激情| 少妇裸体淫交视频免费看高清 | 女性被躁到高潮视频| 黑人操中国人逼视频| 一本久久精品| 无限看片的www在线观看| 下体分泌物呈黄色| 9色porny在线观看| 777米奇影视久久| 日韩大片免费观看网站| 色综合欧美亚洲国产小说| 欧美国产精品一级二级三级| 国产男女内射视频| 黄色丝袜av网址大全| 成人18禁在线播放| 人人妻人人添人人爽欧美一区卜| 91麻豆av在线| 久久久久精品国产欧美久久久| 我的亚洲天堂| 热re99久久国产66热| 国产在视频线精品| 在线 av 中文字幕| 五月天丁香电影| 国产在视频线精品| 久久久久精品国产欧美久久久| 成人黄色视频免费在线看| 欧美亚洲 丝袜 人妻 在线| 亚洲熟女精品中文字幕| aaaaa片日本免费| 欧美一级毛片孕妇| 一本—道久久a久久精品蜜桃钙片| 久9热在线精品视频| 久久精品91无色码中文字幕| 色精品久久人妻99蜜桃| 久久久水蜜桃国产精品网| 人人妻人人澡人人爽人人夜夜| 国产亚洲精品久久久久5区| 黄频高清免费视频| 狠狠婷婷综合久久久久久88av| 亚洲国产精品一区二区三区在线| 丝袜美腿诱惑在线| 亚洲精品美女久久久久99蜜臀| 久久久久精品人妻al黑| 天天添夜夜摸| 人成视频在线观看免费观看| 香蕉国产在线看| 免费在线观看黄色视频的| 黄频高清免费视频| 天堂动漫精品| 国产在线免费精品| 最近最新中文字幕大全免费视频| 国产黄色免费在线视频| 中文字幕人妻熟女乱码| 中亚洲国语对白在线视频| 黄色视频,在线免费观看| 亚洲精品美女久久久久99蜜臀| 91麻豆精品激情在线观看国产 | 久久久久国内视频| 在线观看免费日韩欧美大片| 99热网站在线观看| 一本色道久久久久久精品综合| 欧美亚洲 丝袜 人妻 在线| 满18在线观看网站| 99精国产麻豆久久婷婷| 国产欧美日韩精品亚洲av| 考比视频在线观看| 国产成人系列免费观看| 嫩草影视91久久| 精品国产超薄肉色丝袜足j| 老司机亚洲免费影院| 搡老熟女国产l中国老女人| 国精品久久久久久国模美| 啦啦啦在线免费观看视频4| 老司机福利观看| 午夜激情久久久久久久| 国产精品电影一区二区三区 | 不卡av一区二区三区| 午夜日韩欧美国产| 亚洲欧美日韩高清在线视频 | 亚洲少妇的诱惑av| 两个人免费观看高清视频| 亚洲av日韩精品久久久久久密| 久久久久久久久久久久大奶| 别揉我奶头~嗯~啊~动态视频| 日韩有码中文字幕| 国产麻豆69| 自拍欧美九色日韩亚洲蝌蚪91| 黄色怎么调成土黄色| 色婷婷久久久亚洲欧美| 中文字幕制服av| 精品人妻在线不人妻| 国产主播在线观看一区二区| 老司机靠b影院| 精品福利观看| 国产精品98久久久久久宅男小说| 成人国产一区最新在线观看| 人妻 亚洲 视频| 他把我摸到了高潮在线观看 | 亚洲天堂av无毛| 俄罗斯特黄特色一大片| 亚洲伊人色综图| 日韩成人在线观看一区二区三区| 国产高清激情床上av| 日韩欧美一区视频在线观看| 国产在线免费精品| 一区二区三区精品91| 日韩大片免费观看网站| av一本久久久久| 亚洲国产欧美网| 亚洲熟女毛片儿| 人人妻,人人澡人人爽秒播| 男人舔女人的私密视频| 在线观看一区二区三区激情| 一二三四社区在线视频社区8| 亚洲av日韩精品久久久久久密| 久久精品人人爽人人爽视色| 天天操日日干夜夜撸| 69精品国产乱码久久久| 老司机福利观看| 午夜日韩欧美国产| 丰满饥渴人妻一区二区三| 自线自在国产av| 黑人巨大精品欧美一区二区mp4| 国产精品1区2区在线观看. | 久久人妻av系列| 热99国产精品久久久久久7| 精品高清国产在线一区| 色综合婷婷激情| 高清欧美精品videossex| 曰老女人黄片| 黑人巨大精品欧美一区二区mp4| 免费在线观看影片大全网站| 日韩制服丝袜自拍偷拍| 少妇精品久久久久久久| 久久精品熟女亚洲av麻豆精品| 人妻一区二区av| a级毛片在线看网站| 亚洲一卡2卡3卡4卡5卡精品中文| av超薄肉色丝袜交足视频| 大型av网站在线播放| a级毛片在线看网站| 欧美亚洲 丝袜 人妻 在线| 黄网站色视频无遮挡免费观看| 黄色成人免费大全| 一边摸一边做爽爽视频免费| 成人特级黄色片久久久久久久 | 久久精品国产亚洲av香蕉五月 | 久热这里只有精品99| 叶爱在线成人免费视频播放| 欧美精品一区二区免费开放| 欧美国产精品一级二级三级| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品二区激情视频| 欧美精品人与动牲交sv欧美| 色婷婷av一区二区三区视频| 一级毛片女人18水好多| 色播在线永久视频| 国产精品一区二区精品视频观看| 精品亚洲乱码少妇综合久久| 视频区欧美日本亚洲| 免费观看a级毛片全部| 两个人看的免费小视频| 麻豆成人av在线观看| 欧美日韩亚洲高清精品| 国产精品久久久久久精品古装| 久久婷婷成人综合色麻豆| 亚洲av成人不卡在线观看播放网| 蜜桃在线观看..| 狠狠狠狠99中文字幕| 精品久久久久久电影网| av国产精品久久久久影院| 免费在线观看影片大全网站| 久久精品国产亚洲av香蕉五月 | 12—13女人毛片做爰片一| 精品亚洲成a人片在线观看| 国产精品九九99| 久久久久国内视频| 国产97色在线日韩免费| 国产1区2区3区精品| 久久久久久人人人人人| kizo精华| 亚洲一区二区三区欧美精品| 亚洲精品乱久久久久久| 啦啦啦中文免费视频观看日本| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 日韩精品免费视频一区二区三区| 老汉色av国产亚洲站长工具| 黄频高清免费视频| 在线观看免费午夜福利视频| 黄色怎么调成土黄色| 在线观看66精品国产| 精品一品国产午夜福利视频| 大陆偷拍与自拍| 久久国产精品大桥未久av| 亚洲av电影在线进入| 制服诱惑二区| 久久午夜综合久久蜜桃| 欧美精品亚洲一区二区| 国产日韩欧美视频二区| 女人高潮潮喷娇喘18禁视频| 热99久久久久精品小说推荐| 成在线人永久免费视频| 欧美日韩成人在线一区二区| 午夜免费鲁丝| 亚洲国产欧美日韩在线播放| 欧美日本中文国产一区发布| 国内毛片毛片毛片毛片毛片| 在线观看免费日韩欧美大片| 欧美黄色淫秽网站| 我的亚洲天堂| 亚洲三区欧美一区| 他把我摸到了高潮在线观看 | 久久 成人 亚洲| 久久久久网色| 美女视频免费永久观看网站| 免费黄频网站在线观看国产| a级片在线免费高清观看视频| videosex国产| 最近最新免费中文字幕在线| 人人妻人人澡人人看| 国产成人啪精品午夜网站| 一级毛片精品| 亚洲欧美精品综合一区二区三区| 国产深夜福利视频在线观看| 国产高清激情床上av| 亚洲av日韩精品久久久久久密| a级毛片在线看网站| 亚洲国产av影院在线观看| 成人影院久久| 国产在视频线精品| 欧美亚洲日本最大视频资源| 国产精品一区二区在线观看99| 狂野欧美激情性xxxx| 少妇猛男粗大的猛烈进出视频| 久久久国产成人免费| 亚洲av片天天在线观看| 99国产极品粉嫩在线观看| 欧美激情高清一区二区三区| 久久人妻av系列| 三上悠亚av全集在线观看| 亚洲久久久国产精品| av福利片在线| 一本综合久久免费| 99re在线观看精品视频| 高清毛片免费观看视频网站 | 自拍欧美九色日韩亚洲蝌蚪91| 国产激情久久老熟女| 水蜜桃什么品种好| 亚洲精品一二三| bbb黄色大片| 欧美国产精品va在线观看不卡| 国产日韩欧美亚洲二区| 国产日韩欧美视频二区| 亚洲三区欧美一区| 大香蕉久久网| 2018国产大陆天天弄谢| 丝袜在线中文字幕| 叶爱在线成人免费视频播放| 日本wwww免费看| 99国产精品一区二区三区| 成人国产一区最新在线观看| 精品人妻在线不人妻| 久久精品国产亚洲av高清一级| 婷婷丁香在线五月| 亚洲精品久久午夜乱码| 亚洲成人免费av在线播放| 亚洲七黄色美女视频| 首页视频小说图片口味搜索| 在线 av 中文字幕| 久久精品国产a三级三级三级| 18禁美女被吸乳视频| 久久人妻av系列| 乱人伦中国视频| 亚洲一码二码三码区别大吗| av网站在线播放免费| 国产亚洲一区二区精品| 黄色视频不卡| 99精品久久久久人妻精品| 我的亚洲天堂| 亚洲专区中文字幕在线| 日韩视频一区二区在线观看| 亚洲人成伊人成综合网2020| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 麻豆av在线久日| 亚洲精品久久午夜乱码| 国产精品久久久久久精品古装| 俄罗斯特黄特色一大片| 日韩精品免费视频一区二区三区| 日韩人妻精品一区2区三区| av国产精品久久久久影院| tube8黄色片| 亚洲人成电影观看| 男男h啪啪无遮挡| 男女无遮挡免费网站观看| 99精品久久久久人妻精品| 在线十欧美十亚洲十日本专区| 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 国产主播在线观看一区二区| 丰满人妻熟妇乱又伦精品不卡| 一区福利在线观看| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费视频网站a站| 亚洲综合色网址| 搡老岳熟女国产| 成人国产av品久久久| 国产在线免费精品| 美女福利国产在线| 国产三级黄色录像| 大片免费播放器 马上看| 亚洲精品粉嫩美女一区| 美女福利国产在线| 99精国产麻豆久久婷婷| 欧美人与性动交α欧美软件| 在线看a的网站| 国产精品熟女久久久久浪| 欧美激情高清一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲av第一区精品v没综合| 久久精品亚洲精品国产色婷小说| 国产亚洲精品久久久久5区| 成人永久免费在线观看视频 | 久热爱精品视频在线9| av超薄肉色丝袜交足视频| 午夜激情av网站| av网站免费在线观看视频| 欧美日韩av久久| av网站免费在线观看视频| 水蜜桃什么品种好| 高清毛片免费观看视频网站 | 怎么达到女性高潮| 国产三级黄色录像| 丰满少妇做爰视频| 亚洲七黄色美女视频| 欧美精品啪啪一区二区三区| cao死你这个sao货| 亚洲精品成人av观看孕妇| 免费看十八禁软件| 不卡一级毛片| 国产国语露脸激情在线看| 韩国精品一区二区三区| 日韩大片免费观看网站| 国产淫语在线视频| 成人影院久久| 精品人妻熟女毛片av久久网站| 青青草视频在线视频观看| 十八禁网站免费在线| av国产精品久久久久影院| 亚洲精品粉嫩美女一区| 99九九在线精品视频| 叶爱在线成人免费视频播放| 久久久久网色| 日韩大码丰满熟妇| 丝袜在线中文字幕| 老司机午夜十八禁免费视频| 免费少妇av软件| 免费久久久久久久精品成人欧美视频| 91麻豆精品激情在线观看国产 | 夫妻午夜视频| 91麻豆精品激情在线观看国产 | svipshipincom国产片| 国产精品久久久人人做人人爽| 制服人妻中文乱码| 精品人妻熟女毛片av久久网站| 啦啦啦在线免费观看视频4| 国产免费现黄频在线看| 夜夜骑夜夜射夜夜干| 高清黄色对白视频在线免费看| 精品福利观看| 久久九九热精品免费| 欧美日韩亚洲高清精品| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 午夜激情久久久久久久| 久久中文看片网| 精品少妇久久久久久888优播| 美女高潮到喷水免费观看| 99九九在线精品视频| 天堂中文最新版在线下载| 下体分泌物呈黄色| 欧美精品一区二区免费开放| www.自偷自拍.com| 91麻豆av在线| 18禁国产床啪视频网站| 国产1区2区3区精品| 久久国产精品大桥未久av| 欧美 亚洲 国产 日韩一| 满18在线观看网站| 日韩中文字幕欧美一区二区| 女警被强在线播放| 国产在线观看jvid| 搡老乐熟女国产| 国产在线视频一区二区| 香蕉久久夜色| 夜夜爽天天搞| 老汉色∧v一级毛片| 国产熟女午夜一区二区三区| 建设人人有责人人尽责人人享有的| 国产成人精品久久二区二区91| 精品国产乱码久久久久久小说| 丝袜美腿诱惑在线| 俄罗斯特黄特色一大片| 黄色怎么调成土黄色| 国产又色又爽无遮挡免费看| 极品人妻少妇av视频| 蜜桃国产av成人99| 男女边摸边吃奶| av不卡在线播放| 成人亚洲精品一区在线观看| 久久人妻av系列| 亚洲第一av免费看| aaaaa片日本免费| 欧美 日韩 精品 国产| tocl精华| 日韩欧美国产一区二区入口| 熟女少妇亚洲综合色aaa.| 亚洲视频免费观看视频| 午夜福利视频精品| 99国产精品一区二区蜜桃av | 少妇精品久久久久久久| 免费av中文字幕在线| 十八禁网站网址无遮挡| 日日摸夜夜添夜夜添小说| 黑丝袜美女国产一区| 亚洲va日本ⅴa欧美va伊人久久| 日韩人妻精品一区2区三区| 亚洲精品国产区一区二| 欧美黄色淫秽网站| 亚洲视频免费观看视频| 亚洲成人手机| 亚洲一码二码三码区别大吗| 国产精品久久久久久精品古装| 日韩一卡2卡3卡4卡2021年| 日韩精品免费视频一区二区三区| xxxhd国产人妻xxx| 一夜夜www| 亚洲成人免费电影在线观看| 亚洲人成伊人成综合网2020| 人妻 亚洲 视频| 免费在线观看日本一区| 51午夜福利影视在线观看| 日韩一区二区三区影片| 一区二区av电影网| av在线播放免费不卡| 久久人人爽av亚洲精品天堂| 99国产综合亚洲精品| 水蜜桃什么品种好| 99热国产这里只有精品6| 免费少妇av软件| 黑人欧美特级aaaaaa片| 国产精品av久久久久免费| 国产精品免费视频内射| 黄色片一级片一级黄色片| 精品国产一区二区三区四区第35| 亚洲欧美日韩高清在线视频 | 变态另类成人亚洲欧美熟女 | 国产主播在线观看一区二区| 国精品久久久久久国模美| 精品免费久久久久久久清纯 | 亚洲精品美女久久av网站| 成人永久免费在线观看视频 | av天堂在线播放| 国产精品久久电影中文字幕 | 国产成人精品无人区| 超碰97精品在线观看| 欧美黄色淫秽网站| 午夜福利,免费看| 国产亚洲精品一区二区www | 精品少妇久久久久久888优播| 精品国产乱码久久久久久男人| 国产亚洲av高清不卡| 精品一区二区三区四区五区乱码| 国产主播在线观看一区二区| 91精品国产国语对白视频| 纵有疾风起免费观看全集完整版| 老司机靠b影院| 大码成人一级视频| 欧美老熟妇乱子伦牲交| netflix在线观看网站| 视频在线观看一区二区三区| 国产精品影院久久| 欧美大码av| 天天添夜夜摸| 中文亚洲av片在线观看爽 | 免费在线观看视频国产中文字幕亚洲| 成人手机av| 国产一区有黄有色的免费视频| 黑人欧美特级aaaaaa片| 黄频高清免费视频| 美女视频免费永久观看网站| 大片电影免费在线观看免费| 大香蕉久久网| 黄色成人免费大全| av片东京热男人的天堂| 欧美黄色片欧美黄色片| 高清欧美精品videossex| 国产精品偷伦视频观看了| 精品亚洲成国产av| 淫妇啪啪啪对白视频| 国产伦人伦偷精品视频| 丝袜在线中文字幕| 亚洲avbb在线观看| 大香蕉久久网| 一级a爱视频在线免费观看| 三级毛片av免费| 在线永久观看黄色视频| 欧美乱码精品一区二区三区| 青草久久国产| 精品国产一区二区三区四区第35| 日日摸夜夜添夜夜添小说| 老鸭窝网址在线观看| 亚洲专区中文字幕在线| 一级毛片电影观看| 精品国产一区二区三区久久久樱花| 成人特级黄色片久久久久久久 | 日韩欧美国产一区二区入口| 久久ye,这里只有精品| 欧美日韩亚洲综合一区二区三区_| 欧美成人午夜精品| 久久精品国产99精品国产亚洲性色 | 狠狠婷婷综合久久久久久88av| 超碰97精品在线观看| 老司机福利观看| 成人免费观看视频高清| 欧美成人午夜精品| 日韩中文字幕视频在线看片| 岛国毛片在线播放| 久热爱精品视频在线9| 国产成人免费观看mmmm| 十八禁网站免费在线| 下体分泌物呈黄色| 性少妇av在线| 国产伦理片在线播放av一区| 757午夜福利合集在线观看| 一进一出抽搐动态| 一进一出好大好爽视频| 一本大道久久a久久精品| 国产欧美日韩一区二区精品| 成年动漫av网址| 黑人操中国人逼视频| 欧美成狂野欧美在线观看| av在线播放免费不卡| 精品卡一卡二卡四卡免费| 免费人妻精品一区二区三区视频| 丝瓜视频免费看黄片| 国产亚洲欧美精品永久| 宅男免费午夜| 亚洲中文日韩欧美视频| 国产精品免费一区二区三区在线 | 成人免费观看视频高清| 女性被躁到高潮视频| 欧美日韩成人在线一区二区| 国产深夜福利视频在线观看| 午夜成年电影在线免费观看| 黄色视频不卡| 女性被躁到高潮视频| 大香蕉久久成人网| 亚洲专区中文字幕在线| 国产成人啪精品午夜网站| 在线天堂中文资源库| 欧美 亚洲 国产 日韩一| 高清毛片免费观看视频网站 | 在线观看免费午夜福利视频| 国产真人三级小视频在线观看| 久久人妻av系列| 精品国产乱码久久久久久男人| 天堂俺去俺来也www色官网| 亚洲成a人片在线一区二区| 国产不卡一卡二| 制服人妻中文乱码| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲一级av第二区|