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

    家養(yǎng)雙峰駝SRY基因的序列分析與拷貝數(shù)研究

    2017-11-29 03:34:34陳慧玲任戰(zhàn)軍葉翔楊趙紀(jì)萍張成東
    西北農(nóng)業(yè)學(xué)報 2017年11期
    關(guān)鍵詞:雙峰駝家養(yǎng)拷貝

    陳慧玲,任戰(zhàn)軍,葉翔楊,趙紀(jì)萍,張成東

    (西北農(nóng)林科技大學(xué) 動物科技學(xué)院,陜西楊陵 712100)

    家養(yǎng)雙峰駝SRY基因的序列分析與拷貝數(shù)研究

    陳慧玲,任戰(zhàn)軍,葉翔楊,趙紀(jì)萍,張成東

    (西北農(nóng)林科技大學(xué) 動物科技學(xué)院,陜西楊陵 712100)

    為探究中國家養(yǎng)雙峰駝SRY基因的序列以及群體間拷貝數(shù)變異,對7個中國家養(yǎng)雙峰駝群體共94個個體的SRY基因部分序列進(jìn)行測序,分析SRY基因序列與其他物種的相似性;并且以 CYP2A基因?yàn)閮?nèi)參基因,運(yùn)用實(shí)時熒光定量PCR技術(shù)檢測SRY基因的拷貝數(shù)。結(jié)果表明,中國家養(yǎng)雙峰駝與人、馬、豬、羊及普通牛的SRY基因相似度分別為67.7%~68.1%、68.5%~69.4%、73.8%~74.4%、73.8%~74.8%、74.1%~74.8%。在家養(yǎng)雙峰駝SRY基因的多個拷貝中,除1個拷貝外,其余拷貝的HMG序列與羊駝、人、馬、豬、綿羊及普通牛的HMG序列相比,相似度分別為97.3%、84.2%、85.6%、86.3%、87.0%、85.6%。SRY基因在中國家養(yǎng)雙峰駝Y染色體上以多拷貝形式存在,變幅為 1~9 個拷貝。表明中國家養(yǎng)雙峰駝SRY基因的HMG序列是高度保守的,而且該基因?qū)儆诙嗫截惢颉?/p>

    家養(yǎng)雙峰駝;SRY基因;拷貝數(shù)

    性別決定基因(Sex-determining-region on Y chromosome,SRY)存在于哺乳動物Y染色體的雄性特異區(qū),具有性別決定功能,參與睪丸組織的發(fā)育過程。該基因編碼的蛋白包含C端、N端和高泳動類非組蛋白(High-mobility group, HMG)區(qū)域,C端與N端在物種間保守性較弱,而HMG區(qū)域在物種間具有高度保守性[1-2]。早期對SRY基因的研究發(fā)現(xiàn),SRY保守區(qū)的突變會破壞其雄性性別決定功能,引起性別反轉(zhuǎn)[3-4],SRY基因也能夠誘導(dǎo)基因型為XX的雌性小鼠向雄性發(fā)育[5]。SRY蛋白能夠識別并結(jié)合 SOX9(Sry-related HMG box-9)基因上游的增強(qiáng)子[6],使睪丸支持細(xì)胞的前體細(xì)胞中 SOX9基因的轉(zhuǎn)錄水平顯著上調(diào),進(jìn)而促進(jìn)下游基因的轉(zhuǎn)錄,誘導(dǎo)睪丸支持細(xì)胞的分化。隨后,分化后的睪丸支持細(xì)胞組裝成睪丸索,并且刺激生殖細(xì)胞、睪丸間質(zhì)細(xì)胞、睪丸血管細(xì)胞以及其他間質(zhì)細(xì)胞的分化,最終形成睪丸組織。

    在多數(shù)哺乳動物中,SRY基因是以單拷貝的形式存在[7],但是少數(shù)物種中該基因存在多拷貝現(xiàn)象。貓、小鼠、家兔和水牛中SRY為多拷貝基因[8-12]。接觸輻射或性染色體異常的人類SRY基因會出現(xiàn)拷貝數(shù)變異現(xiàn)象[13]。迄今為止,未見中國家養(yǎng)雙峰駝SRY基因研究的報道。本試驗(yàn)以中國4個地區(qū)7個家養(yǎng)雙峰駝群體為研究對象,對其SRY基因部分序列進(jìn)行測序,并且運(yùn)用實(shí)時熒光定量PCR技術(shù)檢測其拷貝數(shù),以期對家養(yǎng)雙峰駝SRY基因有初步了解。

    1 材料與方法

    1.1 試驗(yàn)材料

    按照隨機(jī)整群抽樣的方法,采集內(nèi)蒙古(阿拉善、蘇尼特)、甘肅(河西)、青海(海西州)、新疆(南疆、北疆和東疆)4個地區(qū)94峰雄性家養(yǎng)雙峰駝的血樣(表1)。另外,采集5峰雌性家養(yǎng)雙峰駝血樣作為陰性對照。頸靜脈采集血液于已滅菌的離心管中,加入ACD抗凝劑后在-80 ℃保存,備用。

    表1 家養(yǎng)雙峰駝樣本來源與樣本量Table 1 Source and number of domestic bactrian camel

    1.2 SRY引物設(shè)計(jì)與合成

    依據(jù)GenBank數(shù)據(jù)庫中雄性野生雙峰駝的相應(yīng)序列(GenBank登錄號:NW_006220067),運(yùn)用Primer 5.0軟件設(shè)計(jì)預(yù)期擴(kuò)增片段長度為708 bp 的引物用于SRY基因的測序。其中上游引物:5′-ATGCTTCTGCTATGTTTGCG-3′,下游引物:5′-ACCAAAAGTAACGGTGAGAATG-3′。引物由上海生工生物工程股份有限公司合成。

    1.3 PCR擴(kuò)增與測序

    利用酚—氯仿法提取基因組DNA[14],采用10 g/L瓊脂糖凝膠電泳和Nanodrop 2000檢測其純度和質(zhì)量濃度。將提取效果良好的基因組DNA樣品統(tǒng)一稀釋至20 ng/μL,置于-20 ℃保存,備用。

    挑選6~10個雄性個體的基因組DNA混合成DNA池,共構(gòu)建13個DNA池,-4 ℃保存。DNA池中每個樣品的質(zhì)量濃度盡可能保持一致。

    首先用SRY引物對雄性DNA池進(jìn)行PCR擴(kuò)增(以雌性雙峰駝DNA為陰性對照),采用10 g/L瓊脂糖凝膠電泳檢測引物的雄性特異性,并且確定其最適退火溫度為65.3 ℃。反應(yīng)體系為10 μL:EsTaqMaster Mix 5 μL,上下游引物各0.2 μL,DNA模板0.5 μL,ddH2O 4.1 μL(北京康為世紀(jì)有限公司)。 PCR反應(yīng)程序:95 ℃變性5 min;94 ℃ 30 s,65 ℃ 40 s, 72 ℃ 50 s,34個循環(huán);72 ℃ 10 min。雄性特異的擴(kuò)增產(chǎn)物由上海生工生物工程股份有限公司進(jìn)行測序。

    如果測序結(jié)果出現(xiàn)雙峰,則對該DNA池中的每個個體分別進(jìn)行PCR擴(kuò)增測序,以再次確定SNP位點(diǎn)。如果單個個體的序列峰圖仍在相同位點(diǎn)處出現(xiàn)雙峰,則需要進(jìn)行TA克隆測序予以驗(yàn)證。在每個群體中分別挑選1個個體進(jìn)行PCR擴(kuò)增,用瓊脂糖凝膠DNA回收試劑盒(Takara公司)純化回收PCR產(chǎn)物,產(chǎn)物經(jīng)過純化與回收之后,與pGEM-T載體連接,然后轉(zhuǎn)化至大腸桿菌感受態(tài)細(xì)胞,在37 ℃平板培養(yǎng)。每個個體隨機(jī)挑選10個陽性菌落在具有Amp+的LB培養(yǎng)基中震蕩培養(yǎng)過夜,吸取部分菌液送至上海生工生物工程股份有限公司進(jìn)行克隆測序。

    1.4 SRY序列相似性分析

    用DNASTAR軟件對所有個體的序列進(jìn)行比對,分析個體間序列的差異。運(yùn)用MegAlign軟件分析人(GenBank登錄號:JQ811934.1)、馬(GenBank登錄號:NM_001081810.1)、豬(GenBank登錄號:U49860.3)、綿羊(GenBank登錄號:Z30265.1)及普通牛(GenBank登錄號:Z30327.1)的SRY基因序列與測序所得的SRY基因部分序列的相似度以及物種間HMG區(qū)DNA序列的相似度。

    1.5 熒光定量引物的設(shè)計(jì)與合成

    TA克隆測序完成后,運(yùn)用MegAlign軟件比對所有序列,選擇在SRY基因的保守區(qū)設(shè)計(jì)SRYa引物用于實(shí)時熒光定量PCR擴(kuò)增。選擇常染色體上的單拷貝基因細(xì)胞色素P4502A( CYP2A)為內(nèi)參基因[14],計(jì)算SRY基因的拷貝數(shù)。引物信息見表2。

    1.6 熒光定量PCR擴(kuò)增

    采用實(shí)時熒光定量PCR技術(shù)檢測基因的拷貝數(shù)。首先挑選質(zhì)量較好的DNA模板用于標(biāo)準(zhǔn)曲線的制作,質(zhì)量濃度梯度依次為40、20、10、5和2.5 ng/μL,每個質(zhì)量濃度設(shè)置3個重復(fù)。熒光定量PCR程序完成之后,獲得各質(zhì)量濃度相對應(yīng)的Ct值,分別制作CYP2A和SRYa的標(biāo)準(zhǔn)曲線。將待測樣本均稀釋至5 ng/μL,每個樣本在同一板中同時進(jìn)行2對引物的擴(kuò)增,而且每對引物對應(yīng)的每個DNA樣本均設(shè)置3個重復(fù),每板都設(shè)有相同的校正樣本以減小板間的系統(tǒng)誤差。

    表2 SRYa和CYP2A引物信息Table 2 Information of SRYa and CYP2A

    反應(yīng)體系:DNA模板0.8 μL,SYBR Premix EXTaqⅡ(2×)(Takara) 5 μL,上下游引物(10 μmol/L)各0.4 μL,RNase-Free Water 3.4 μL。熒光定量反應(yīng)程序?yàn)椋?5 ℃ 預(yù)變性30 s;95 ℃變性5 s,退火30 s,35個循環(huán);95 ℃ 10 s,65 ℃ 5 s,65 ℃~95 ℃(0.5 ℃/s)。

    1.7 基因拷貝數(shù)的計(jì)算

    參照HAMILTON等[15]的方法計(jì)算SRY基因的拷貝數(shù)。計(jì)算公式如下:

    依據(jù)獲取的各質(zhì)量濃度梯度所對應(yīng)的Ct值,繪制標(biāo)準(zhǔn)曲線,可以獲得標(biāo)準(zhǔn)曲線的斜率k,進(jìn)而計(jì)算出引物的擴(kuò)增效率:E=10-1/k。

    依據(jù)所有96孔板中校正樣本的Ct值,計(jì)算校正樣本Ct值的平均值。再根據(jù)平均值計(jì)算校正樣本的拷貝數(shù):CN校正樣本=(Er)Ct r/(Et)Ct t。

    每個96孔板中校正樣本都設(shè)置3個重復(fù),分別計(jì)算不同板中校正樣本3個Ct值的平均值,再將不同板的待測樣本校正到同一水平,待測樣本與校正樣本之間校正系數(shù)的計(jì)算:

    R=(Et)△Ct t/(Er)△Ct r。

    待測樣本拷貝數(shù)計(jì)算:CN=(CN校正樣本) ×R。

    k代表標(biāo)準(zhǔn)曲線的斜率,t代表目的基因SRY,r代表內(nèi)參基因 CYP2A,E代表擴(kuò)增效應(yīng),Et代表SRYa引物的擴(kuò)增效率,Er代表CYP2A 引物的擴(kuò)增效率,Ctr代表CYP2A引物的Ct值,Ctt代表SRYa引物的Ct值,△Ctt=校正樣本SRYa的Ct值-待測樣本SRYa 的Ct值,△Ctr=校正樣本CYP2A 的Ct值-待測樣本CYP2A 的Ct值。

    運(yùn)用SPSS 18.0軟件Kolmogorov-Smirnov標(biāo)準(zhǔn)對所研究個體的基因拷貝數(shù)進(jìn)行正態(tài)分布檢驗(yàn)。

    2 結(jié)果與分析

    2.1 SRY基因的PCR擴(kuò)增

    Touchdown PCR擴(kuò)增后,擴(kuò)增產(chǎn)物經(jīng)10 g/L的瓊脂糖凝膠電泳檢測,結(jié)果如圖1所示,陰性對照與空白對照無擴(kuò)增產(chǎn)物,而所有雄性個體的條帶明亮清晰,片段長度與預(yù)期長度一致,并且無非特異性擴(kuò)增,表明SRY引物具有雄性特異性。

    圖1 SRY引物的雄性特異性檢測Fig.1 Detection of SRY specificity in male camel

    2.2 SRY基因序列分析

    運(yùn)用DNASTAR軟件對SRY引物擴(kuò)增所得序列進(jìn)行比對,結(jié)果如圖2所示,13個DNA池中SRY基因擴(kuò)增片段的測序峰圖均出現(xiàn)雙峰,本圖只展示其中3個DNA池的比對結(jié)果,其余DNA池測序結(jié)果與該圖一致,因此初步推斷SRY基因存在SNP位點(diǎn)以及堿基的插入。

    以DNA池中各單樣為模板分別對SRY引物進(jìn)行PCR擴(kuò)增并對擴(kuò)增產(chǎn)物測序,經(jīng)過序列比對,對應(yīng)DNA池中SNP位點(diǎn)處尋找各個單樣之間堿基的差異。結(jié)果如圖2所示,每個個體的測序峰圖在與池DNA序列的相同位點(diǎn)處存在雙峰,而且也具有堿基插入位點(diǎn),其他個體測序結(jié)果也與該結(jié)果一致,由此推測在家養(yǎng)雙峰駝中SRY基因可能是多拷貝基因。

    TA克隆測序結(jié)果如圖3所示,同一個體10個克隆之間的序列存在變異位點(diǎn)和1個AT二堿基插入,本圖只展示3個克隆序列的比對結(jié)果,其他群體的測序結(jié)果與青海駝一致。由于Y染色體MSY區(qū)不與X染色體重組,因而證明在家養(yǎng)雙峰駝中單個個體出現(xiàn)SNP位點(diǎn)的SRY基因是多拷貝基因。

    圖2 3個DNA池的SRY基因測序Fig.2 Sequencing result of the SRY gene of three gene pools

    圖3 1峰青海駝的SRY基因克隆測序Fig.3 SRY gene of a Qinghai sample by PCR cloning and sequencing

    2.3 不同物種間序列相似性分析

    對中國家養(yǎng)雙峰駝的SRY基因部分序列進(jìn)行克隆測序,共獲得708 bp的序列。運(yùn)用MegAlign軟件對家養(yǎng)雙峰駝SRY基因不同拷貝的序列與人、馬、豬及普通牛的SRY基因序列進(jìn)行比對分析,結(jié)果發(fā)現(xiàn),中國家養(yǎng)雙峰駝SRY基因與人、馬、豬、綿羊及普通牛的SRY基因相似度約為67.7%~68.1%、68.5%~69.4%、73.8%~74.4%、73.8%~74.8%、74.1%~74.8%。

    不同物種的SRY蛋白都含有保守的HMG結(jié)構(gòu)域,該結(jié)構(gòu)域是SRY蛋白重要的功能區(qū)。因此本試驗(yàn)比對羊駝HMG序列和家養(yǎng)雙峰駝708 bp 的序列,結(jié)果顯示在家養(yǎng)雙峰駝SRY基因的多個拷貝中,除1個拷貝外,其余拷貝中有146 bp的序列與羊駝HMG序列相似度高達(dá)97.3%,與人、馬、豬、綿羊及普通牛的HMG序列相比,相似度分別為84.2%、85.6%、86.3%、87.0%、85.6%。因此,可以初步確定中國家養(yǎng)雙峰駝SRY基因的HMG區(qū)域。由于拷貝之間的HMG序列有5個突變位點(diǎn),因而認(rèn)為家養(yǎng)雙峰駝的HMG序列為:

    TCATTGTATGGGCTCGTGATCAAAGGCGAAAGGTGGCTCTASAGAATMCCAAA- ATGCAGAACKCAGAGATCAGCAAGCGGC-TGGGATACCAGTGGAAATTACTTACAGA-AGSTGAAAAGCGGCCGTTCTTCGAGGAG-GCRCAGAGA。

    結(jié)合野生雙峰駝SRY基因的預(yù)測序列、氨基酸序列(GenBank登錄號:XP_006195462)以及羊駝HMG的氨基酸序列(GenBank登錄號:ABI31754.1),推測家養(yǎng)雙峰駝HMG的氨基酸序列為:

    IVWARDQRRKVAL(E/Q)N(T/P)KMQN(A/S)EISKRLGYQWKLLTE(A/G)EKRPFFEE AQR。其與野生雙峰駝和羊駝序列相似度約為95.8%。

    2.4 SRYa和CYP2A引物的PCR擴(kuò)增

    Touchdown PCR擴(kuò)增后,用10 g/L的瓊脂糖凝膠電泳檢測(圖4),結(jié)果顯示SRYa引物只在雄性個體中有目的條帶,雌性個體無條帶。CYP2A引物的擴(kuò)增產(chǎn)物為單一條帶,且無引物二聚體。

    圖4 SRYa(a)和CYP2A(b)引物的特異性檢測Fig.4 Detection of SRYa(a) and CYP2A(b) specificity

    2.5 溶解曲線

    引物的熒光定量溶解曲線如圖5所示,2對引物的擴(kuò)增產(chǎn)物單一,無引物二聚體及非特異性條帶擴(kuò)增。由此可見,該引物均可用于后續(xù)的熒光定量PCR擴(kuò)增。

    2.6 標(biāo)準(zhǔn)曲線

    CYP2A和SRYa引物的標(biāo)準(zhǔn)曲線如圖6所示,CYP2A和SRYa引物的標(biāo)準(zhǔn)曲線斜率分別為-3.207 1和-3.223 3。由此計(jì)算出CYP2A和SRYa引物的擴(kuò)增效率分別為2.05和2.04,R2均大于0.99,表明2對引物的擴(kuò)增效率較好,可用于后續(xù)的試驗(yàn)。

    2.7 拷貝數(shù)計(jì)算

    根據(jù)公式計(jì)算出校正樣本SRY基因的拷貝數(shù)是2,隨后計(jì)算94個樣本的拷貝數(shù)。結(jié)果顯示,在94個個體中,SRY基因拷貝數(shù)中值數(shù)為4,變幅為1~9。各群體拷貝數(shù)見表3。

    圖5 CYP2A和SRYa的 qPCR 溶解曲線Fig.5 Dissociation curves of CYP2A and SRYa

    圖6 CYP2A 和SRY擴(kuò)增引物的標(biāo)準(zhǔn)曲線Fig.6 Standard curves of CYP2A and SRY primers

    表3 7個家養(yǎng)雙峰駝群體SRY基因的拷貝數(shù)中值Table 3 Copy median number (CN) of SRY in seven populations

    94個個體基因拷貝數(shù)的正態(tài)分布分析結(jié)果顯示所有個體的拷貝數(shù)均不服從正態(tài)分布(圖7-a)。箱線圖結(jié)果發(fā)現(xiàn)少數(shù)雙峰駝個體的SRY基因拷貝數(shù)偏高或偏低,大約占群體的1.06%。其中,蘇尼特駱駝有1個個體與其他個體離散,拷貝數(shù)為9(圖7-b)。

    3 討 論

    平均值=4.11,標(biāo)準(zhǔn)差 = 1.44 Mean = 4.11,standard deviation = 1.44

    SRY基因是哺乳動物性別決定的開關(guān)基因,它的表達(dá)能夠啟動雄性發(fā)育的級聯(lián)信號通路,誘導(dǎo)性別特異的性腺發(fā)育,因而在物種的性別分化過程中具有至關(guān)重要的作用。本試驗(yàn)通過對中國家養(yǎng)雙峰駝的SRY基因部分序列進(jìn)行測序,而后將所得序列與其他物種SRY基因序列進(jìn)行比對,結(jié)果顯示其與牛和豬的相似度很高。Jirimutu等[14]通過對動物基因組中的 2 345 個單拷貝直系同源構(gòu)建系統(tǒng)進(jìn)化樹,結(jié)果發(fā)現(xiàn)駱駝與牛和豬的關(guān)系最近。雖然本試驗(yàn)所研究的SRY基因是多拷貝基因,但是物種間的序列相似度也能在一定程度上反映駱駝與牛、豬的進(jìn)化關(guān)系,進(jìn)化關(guān)系越近,序列變異程度越小。因此,該結(jié)果與Jirimutu等[14]的研究結(jié)果一致。

    本試驗(yàn)運(yùn)用實(shí)時熒光定量PCR技術(shù)檢測中國家養(yǎng)雙峰駝SRY基因的拷貝數(shù),結(jié)果發(fā)現(xiàn)拷貝數(shù)變幅為1~9。和其他物種相比,拷貝數(shù)也有差異。目前研究發(fā)現(xiàn),SRY基因在人類、馬、牛、黑猩猩、豬中均為單拷貝基因[7,15-17],而在家貓和水牛中則為多拷貝基因[10,12]。大鼠的SRY基因有11個拷貝,拷貝序列之間有一定的差異[11],11個拷貝中至少有6個拷貝能夠正常表達(dá)[18]。不同物種拷貝數(shù)的差異進(jìn)一步表明Y染色體在不同的物種間表現(xiàn)出異質(zhì)性。

    HMG區(qū)域在物種間相對保守,是SRY基因發(fā)揮性別決定功能的區(qū)域。因此,本試驗(yàn)對中國家養(yǎng)雙峰駝SRY基因的HMG區(qū)域進(jìn)行預(yù)測,結(jié)果顯示,HMG的堿基序列位于本試驗(yàn)所擴(kuò)增片段的193~338 bp,而二堿基插入在該區(qū)域的下游,因而不會引起該區(qū)域序列的移碼突變。由于不同拷貝的序列之間有變異位點(diǎn)存在,部分氨基酸會發(fā)生改變,氨基酸的改變是否會影響SRY基因的功能還是未知數(shù)。目前,對于大鼠SRY基因序列的分析結(jié)果顯示,11個拷貝的SRY序列之間的突變位點(diǎn)導(dǎo)致氨基酸序列改變,進(jìn)而引起核定位以及轉(zhuǎn)錄活性改變[11]。前人[19-20]研究發(fā)現(xiàn)SRY蛋白的功能并不局限于性別決定,其在成年雄性大腦、腎臟和腎上腺中也有表達(dá)。本試驗(yàn)所得708 bp的序列在拷貝序列間也存在多個突變位點(diǎn)及1個二堿基的插入位點(diǎn),推測這些位點(diǎn)的改變可能會影響睪丸組織的發(fā)育。這些變異也可能與SRY基因的其他功能相關(guān)聯(lián),然而這些問題還有待于進(jìn)一步探討。

    Reference:

    [1] LUNDRIGAN B L,TUCKER P K.Tracing paternal ancestry in mice,using the Y-linked,sex-determining locus,Sry[J].MolecularBiologyandEvolution,1994,11(3):483-492.

    [2] NACHMAN M W,AQUADRO C F.Polymorphism and divergence at the 5’ flanking region of the sex-determining locus,Sry,in mice[J].MolecularBiologyandEvolution,1994,11(3):539-547.

    [3] HARLEY V R,JACKSON D T,HEXTALL P J,etal.DNA binding activity of recombinant SRY from normal males and XY females[J].Science,1992,255(5043):453-456.

    [4] HAWKINS J R,TAILOR A,BERTA P,etal.Mutational analysis of the SRY:nonsense and missense mutations in XY sex reversal[J].HumanGenetics,1992,88(4):471-474.

    [5] KOOPMA P,GUBBAY J,VIVIAN N,etal.Male development of chromosomally female mice transgenic for SRY[J].Nature,1991,351(6322):117-121.

    [6] SEKIDO R,LOVELL-BADGE R.Sex determination involves synergistic action of SRY and SF1 on a specific SOX9 enhancer[J].Nature,2008,453(7197):930-934.

    [7] PARIA N,RAUDSEPP T,WILKERSON A J P,etal.A gene catalogue of the euchromatic male-specific region of the horse Y chromosome:comparison with human and other mammals[J].PLoSOne,2011,6(7):e21374.

    [8] TOURE A,SZOT M,MAHADEVAIAH S K,etal.A new deletion of the mouse Y chromosome long arm associated with the loss of Ssty expression,abnormal sperm development and sterility[J].Genetics,2004,166 (2):901-912.

    [9] GERALDES A,FERRAND N.A 7-bp insertion in the 3′ untranslated region suggests the duplication and concerted evolution of the rabbitSRYgene[J].GeneticsSelectionEvolution,2006,38:313-320.

    [10] WILKERSON A J P,RAUDSEPP T,GRAVES T,etal.Gene discovery and comparative analysis of X-degenerate genes from the domestic cat Y chromosome[J].Genomics,2008,92(5):329-338.

    [11] PROKOP J W,UNDERWOOD A C,TURNER M E,etal.Analysis of SRY duplications on the Rattus norvegicus Y-chromosome[J].BMCGenomics,2013,14(1):117-124.

    [12] 徐舒遠(yuǎn).水牛Y-SNPs篩選及多拷貝基因鑒定[D].陜西楊陵:西北農(nóng)林科技大學(xué),2014.

    XU SH Y.Y-SNPs scanning and multi-copy genes identification in water buffalo[D].Yangling Shaanxi:Northwest Aamp;F University,2014 (in Chinese with English abstract).

    [13] PREMI S,SRIVASTAVA J,CHANDY S P,etal.Tandem duplication and copy number polymorphism of theSRYgene in patients with sex chromosome anomalies and males exposed to natural background radiation[J].MolecularHumanReproduction,2006,12(2):113-121.

    [14] JIRIMUTU,WANG Z,DING G,etal.Genome sequences of wild and domestic bactrian camels[J].NatureCommunications,2012,3:1202-1202.

    [15] HAMILTONA C K,FAVETTA L A,DI MEO G P,etal.Copy number variation of testis-specific protein,Y-encoded (TSPY) in 14 different breeds of cattle (Bostaurus)[J].SexualDevelopment,2009,3(4):205-213.

    [16] SKALETSKY H,KURODA-KAWAGUCHI T,MINX P J,etal.The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes[J].Nature,2003,423(6942):825-837.

    [17] CHANG T C,YANG Y,RETZEL E F,etal.Male-specific region of the bovine Y chromosome is gene rich with a high transcriptomic activity in testis development[J].ProceedingsoftheNationalAcademyofSciences,2013,110(30):12373-12378.

    [18] TURNER M E,MARTIN C,MARTINS A S,etal.Genomic and expression analysis of multipleSRYloci from a single Rattus norvegicus Y chromosome[J].BMCGenetics,2007,8(1):1.

    [19] CZECH D P,LEE J,SIM H,etal.The human testis-determining factor SRY localizes in midbrain dopamine neurons and regulates multiple components of catecholamine synthesis and metabolism[J].JournalofNeurochemistry,2012,122(2):260-271.

    [20] DEWINGE P,CHIANG C W K,SINCHAK K,etal.Direct regulation of adult brain function by the male-specific factor SRY[J].CurrentBiology,2006,16(4):415-420.

    CorrespondingauthorREN Zhanjun, male, associate professor. Research area:economic zoology research.E-mail:renzhanjun@nwsuaf.edu.cn

    (責(zé)任編輯:顧玉蘭Responsibleeditor:GUYulan)

    CopyNumberandSequenceAnalysisofSRYGeneinChineseDomesticBactrianCamel

    CHEN Huiling, REN Zhanjun, YE Xiangyang, ZHAO Jiping and ZHANG Chengdong

    (College of Animal Science and Technology, Northwest Aamp;F University, Yangling Shaanxi 712100,China)

    The study aims to exploreSRYgene sequence of Chinese domestic bactrian camel and the copy number variation of SRY among populations. Partial sequence ofSRYgene of 94 samples from seven populations were sequenced to analyze the similarity among species. Additionally, CYP2A gene as a reference gene,SRYgene copy number was detected by real-time PCR. The similarity ofSRYgene between Chinese bactrian camel and human, horse, pig, sheep as well as bovine was 67.7%-68.1%, 68.5%-69.4%, 73.8%-74.4%, 73.8%-74.8% and 74.1%-74.8%, respectively. Except for one copy, the similarity between HMG of the rest copy of Chinese bactrian camel and HMG of alpaca, human, horse, pig, sheep as well as bovine was 97.3%, 84.2%, 85.6%, 86.3%, 87.0% and 85.6%, respectively. TheSRYgene was a multicopy gene in Chinese domestic bactrian camel and copy number variation ranged from 1 to 9. HMG sequences of Chinese domestic bactrian camel were highly conserved. TheSRYgene was a multicopy gene in Chinese domestic bactrian camel.

    Domestic bactrian camel;SRYgene; Copy number

    2016-09-20

    2016-10-25

    The National Natural Science Foundation of China (No. 31172178).

    CHEN Huiling, female, master student. Research area:evaluation, protection and utilization of animal genetic resources. E-mail:chlachj@nwsuaf.edu.cn

    日期:2017-11-17

    網(wǎng)絡(luò)出版地址:http://kns.cnki.net/kcms/detail/61.1220.S.20171117.1101.004.html

    2016-09-20

    2016-10-25

    國家自然科學(xué)基金(31172178)。

    陳慧玲,女,碩士研究生,從事動物遺傳資源評價、保護(hù)和利用研究。E-mail:chlachj@nwsuaf.edu.cn

    任戰(zhàn)軍,男,副教授,主要從事經(jīng)濟(jì)動物研究。E-mail:renzhanjun@nwsuaf.edu.cn

    Q951

    A

    1004-1389(2017)11-1569-08

    猜你喜歡
    雙峰駝家養(yǎng)拷貝
    極危物種
    ——野雙峰駝
    更正
    中國生殖健康(2018年1期)2018-11-06 07:14:38
    看!這些家養(yǎng)植物,賞心悅目但要小心
    百科知識(2018年9期)2018-05-28 19:11:44
    雞吃石子的秘密
    蘇尼特家養(yǎng)雙峰駝mtDNA Cytb基因和D-loop序列的遺傳多樣性
    我家的小烏龜
    朝三暮四
    文件拷貝誰最“給力”
    羅布泊無人區(qū)
    知識窗(2001年7期)2001-03-27 15:27:58
    日韩大尺度精品在线看网址| 久久久久精品国产欧美久久久| 免费看日本二区| 真人做人爱边吃奶动态| 国产熟女欧美一区二区| 国产精品国产三级国产av玫瑰| 狠狠狠狠99中文字幕| 久久久久性生活片| 国产伦在线观看视频一区| 老熟妇仑乱视频hdxx| 婷婷六月久久综合丁香| 99riav亚洲国产免费| 精品99又大又爽又粗少妇毛片| 久久久久久久久中文| 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| 男插女下体视频免费在线播放| 国产精品久久久久久久久免| 老熟妇乱子伦视频在线观看| 97人妻精品一区二区三区麻豆| 精品久久久久久久久av| 美女 人体艺术 gogo| 国产精品电影一区二区三区| 免费观看在线日韩| 床上黄色一级片| 精品欧美国产一区二区三| 老女人水多毛片| 日韩成人伦理影院| 国产亚洲精品综合一区在线观看| 丝袜美腿在线中文| 超碰av人人做人人爽久久| 一级毛片久久久久久久久女| 男女视频在线观看网站免费| 亚洲18禁久久av| 97超视频在线观看视频| 亚洲国产欧美人成| 午夜爱爱视频在线播放| 国产欧美日韩精品一区二区| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 国产亚洲av嫩草精品影院| 国产aⅴ精品一区二区三区波| av天堂中文字幕网| 亚洲熟妇中文字幕五十中出| 热99re8久久精品国产| a级一级毛片免费在线观看| 91久久精品电影网| 露出奶头的视频| 成人av在线播放网站| 黄色欧美视频在线观看| 日韩av不卡免费在线播放| 男女边吃奶边做爰视频| 在线免费观看不下载黄p国产| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线| 一级毛片电影观看 | 欧美三级亚洲精品| 99热只有精品国产| 人妻少妇偷人精品九色| 日韩欧美在线乱码| 日韩中字成人| 中国美女看黄片| 精品人妻一区二区三区麻豆 | 日韩成人av中文字幕在线观看 | 国内揄拍国产精品人妻在线| 99riav亚洲国产免费| 级片在线观看| 午夜福利高清视频| 最近视频中文字幕2019在线8| 国产日本99.免费观看| 中国美白少妇内射xxxbb| 国产av在哪里看| 久久久精品欧美日韩精品| 麻豆乱淫一区二区| 99久久九九国产精品国产免费| 热99在线观看视频| 精品欧美国产一区二区三| 免费人成在线观看视频色| 成人av一区二区三区在线看| 亚洲国产欧美人成| 国产精品伦人一区二区| 一级黄色大片毛片| 国产成人福利小说| 桃色一区二区三区在线观看| 亚洲国产精品合色在线| 久久久久久伊人网av| 午夜福利视频1000在线观看| 午夜视频国产福利| 亚洲成人中文字幕在线播放| 日韩欧美精品免费久久| 久久久精品大字幕| 91午夜精品亚洲一区二区三区| 九九热线精品视视频播放| 一级黄片播放器| 中国美白少妇内射xxxbb| 日本色播在线视频| 亚洲婷婷狠狠爱综合网| 欧美最黄视频在线播放免费| 免费看a级黄色片| 日韩成人av中文字幕在线观看 | 天堂√8在线中文| 成人性生交大片免费视频hd| 亚洲,欧美,日韩| 国产真实伦视频高清在线观看| 国产精品一区二区三区四区久久| 俄罗斯特黄特色一大片| 天天躁日日操中文字幕| 日韩欧美精品免费久久| 中文在线观看免费www的网站| 国产熟女欧美一区二区| 激情 狠狠 欧美| 一级av片app| 国产黄片美女视频| 欧美精品国产亚洲| 韩国av在线不卡| 国产伦在线观看视频一区| 国产一区二区在线观看日韩| 亚洲婷婷狠狠爱综合网| 日韩大尺度精品在线看网址| 美女大奶头视频| 亚洲国产精品成人综合色| av黄色大香蕉| 波多野结衣高清无吗| 欧美一区二区亚洲| 国产aⅴ精品一区二区三区波| 12—13女人毛片做爰片一| a级毛片a级免费在线| 秋霞在线观看毛片| h日本视频在线播放| 日韩欧美三级三区| a级毛色黄片| 国产精品乱码一区二三区的特点| 精品乱码久久久久久99久播| 简卡轻食公司| 天天一区二区日本电影三级| 久久中文看片网| 别揉我奶头~嗯~啊~动态视频| 欧美高清成人免费视频www| 亚洲自拍偷在线| 日韩av在线大香蕉| 日韩一区二区视频免费看| 亚洲美女黄片视频| 亚洲四区av| 亚洲av中文字字幕乱码综合| 美女免费视频网站| 久久久色成人| 国内精品一区二区在线观看| 国产av不卡久久| 成年免费大片在线观看| 插逼视频在线观看| 国产精品精品国产色婷婷| 成熟少妇高潮喷水视频| 91久久精品电影网| 嫩草影院入口| 日韩欧美在线乱码| 少妇丰满av| 久久午夜亚洲精品久久| 国产免费男女视频| 亚洲熟妇熟女久久| 亚洲美女搞黄在线观看 | 亚洲成a人片在线一区二区| 黄色配什么色好看| 色综合站精品国产| 欧美成人精品欧美一级黄| 欧美日韩一区二区视频在线观看视频在线 | 白带黄色成豆腐渣| 国产国拍精品亚洲av在线观看| 国产精华一区二区三区| 欧美绝顶高潮抽搐喷水| 亚洲内射少妇av| 国产单亲对白刺激| 麻豆久久精品国产亚洲av| 99九九线精品视频在线观看视频| 精品久久久久久久久久久久久| 99久久久亚洲精品蜜臀av| 日本免费a在线| 精品人妻一区二区三区麻豆 | 久久午夜亚洲精品久久| 99久久九九国产精品国产免费| 深夜精品福利| 亚洲av电影不卡..在线观看| 国产乱人视频| 毛片一级片免费看久久久久| 亚洲成av人片在线播放无| 热99在线观看视频| 色视频www国产| 在线a可以看的网站| 欧美一区二区国产精品久久精品| 亚洲精品乱码久久久v下载方式| 国产免费一级a男人的天堂| 我要搜黄色片| 又黄又爽又刺激的免费视频.| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 国产黄片美女视频| 少妇裸体淫交视频免费看高清| 老师上课跳d突然被开到最大视频| 欧美+亚洲+日韩+国产| 男女视频在线观看网站免费| 久久久久免费精品人妻一区二区| 天堂动漫精品| 国产精品永久免费网站| 国产午夜精品久久久久久一区二区三区 | 午夜福利成人在线免费观看| 午夜福利高清视频| 男插女下体视频免费在线播放| 日日摸夜夜添夜夜添av毛片| 国产高清视频在线观看网站| 亚洲国产精品久久男人天堂| 国产精品福利在线免费观看| 国产真实乱freesex| 国内精品一区二区在线观看| 乱人视频在线观看| 中文字幕久久专区| 色综合亚洲欧美另类图片| 99热全是精品| 一级黄色大片毛片| av在线播放精品| 可以在线观看毛片的网站| eeuss影院久久| 国产人妻一区二区三区在| 伦理电影大哥的女人| 国产精品日韩av在线免费观看| 色av中文字幕| 日产精品乱码卡一卡2卡三| 午夜福利在线观看吧| 午夜免费男女啪啪视频观看 | 国产在线男女| 日韩三级伦理在线观看| 欧美丝袜亚洲另类| 国产女主播在线喷水免费视频网站 | 麻豆一二三区av精品| 一区福利在线观看| 性欧美人与动物交配| 成人特级黄色片久久久久久久| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| 欧美日韩在线观看h| 成人无遮挡网站| 日韩成人伦理影院| 丰满乱子伦码专区| 国产在视频线在精品| 免费看a级黄色片| 亚洲av免费高清在线观看| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 久久久久久久亚洲中文字幕| 国产高清激情床上av| 久久人人爽人人片av| 插逼视频在线观看| 一区二区三区四区激情视频 | 男人舔奶头视频| 一级毛片久久久久久久久女| 男人狂女人下面高潮的视频| 老师上课跳d突然被开到最大视频| 免费看光身美女| 国产一区亚洲一区在线观看| 久久99热6这里只有精品| h日本视频在线播放| 国产男靠女视频免费网站| 中国美白少妇内射xxxbb| 久久精品夜夜夜夜夜久久蜜豆| 久久久久国内视频| 一本精品99久久精品77| 97碰自拍视频| 亚洲性夜色夜夜综合| 晚上一个人看的免费电影| 久久久久免费精品人妻一区二区| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 在线a可以看的网站| 免费大片18禁| 人人妻人人看人人澡| 日韩三级伦理在线观看| 久久久久久久久久久丰满| 无遮挡黄片免费观看| 麻豆国产av国片精品| 成人亚洲欧美一区二区av| 一夜夜www| 国产激情偷乱视频一区二区| 国产黄a三级三级三级人| 国产精品永久免费网站| 在线观看一区二区三区| 精品日产1卡2卡| 日产精品乱码卡一卡2卡三| 午夜精品在线福利| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 亚洲欧美日韩东京热| a级毛片a级免费在线| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 国产单亲对白刺激| 深夜a级毛片| 一级黄片播放器| 国产高清视频在线播放一区| 1024手机看黄色片| 日本黄大片高清| 97热精品久久久久久| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 欧美潮喷喷水| 精品久久久久久久久av| 日韩大尺度精品在线看网址| 五月玫瑰六月丁香| 亚洲在线观看片| 99久久精品热视频| 国产乱人偷精品视频| 欧美bdsm另类| 亚洲av不卡在线观看| 最好的美女福利视频网| 日本黄色片子视频| 亚洲欧美日韩无卡精品| 国产v大片淫在线免费观看| 免费看a级黄色片| 最近的中文字幕免费完整| 不卡一级毛片| 午夜精品在线福利| 国产免费一级a男人的天堂| 久久久久九九精品影院| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看| 国产视频一区二区在线看| 日日摸夜夜添夜夜爱| 欧美+日韩+精品| 亚洲精品国产成人久久av| 中文字幕久久专区| 国产一区二区在线av高清观看| 国产高清三级在线| 天堂动漫精品| 一级毛片我不卡| 99久久精品一区二区三区| av国产免费在线观看| 97超级碰碰碰精品色视频在线观看| 22中文网久久字幕| 最近的中文字幕免费完整| 国产女主播在线喷水免费视频网站 | av在线老鸭窝| 久久久久久久久久黄片| 又黄又爽又刺激的免费视频.| 99热网站在线观看| 尤物成人国产欧美一区二区三区| 99久久中文字幕三级久久日本| 亚洲图色成人| 久久99热这里只有精品18| 亚洲欧美精品综合久久99| 亚洲国产精品sss在线观看| 麻豆精品久久久久久蜜桃| 亚洲激情五月婷婷啪啪| 中文字幕精品亚洲无线码一区| 久久热精品热| 久久久久性生活片| 欧美日韩在线观看h| 国产69精品久久久久777片| 人人妻,人人澡人人爽秒播| 在线免费观看不下载黄p国产| 日本免费一区二区三区高清不卡| 欧美xxxx性猛交bbbb| 两个人的视频大全免费| 国产精品美女特级片免费视频播放器| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 国产三级中文精品| 最新在线观看一区二区三区| 美女大奶头视频| 国产v大片淫在线免费观看| 美女免费视频网站| 亚洲精品乱码久久久v下载方式| 日韩精品有码人妻一区| 直男gayav资源| 久久99热这里只有精品18| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 色尼玛亚洲综合影院| 九九久久精品国产亚洲av麻豆| 欧美国产日韩亚洲一区| 一进一出好大好爽视频| 99热这里只有是精品50| 男女边吃奶边做爰视频| 中文资源天堂在线| 一个人免费在线观看电影| 国产精品亚洲美女久久久| 午夜亚洲福利在线播放| 特级一级黄色大片| 一个人观看的视频www高清免费观看| 老熟妇乱子伦视频在线观看| 久久久久国内视频| 我的老师免费观看完整版| 高清毛片免费观看视频网站| 国产伦在线观看视频一区| av国产免费在线观看| 亚洲国产精品sss在线观看| 日日干狠狠操夜夜爽| 日本在线视频免费播放| 人妻夜夜爽99麻豆av| 亚洲高清免费不卡视频| 久久久欧美国产精品| 久久精品国产自在天天线| 午夜爱爱视频在线播放| 久久久久免费精品人妻一区二区| 久久亚洲精品不卡| 中出人妻视频一区二区| 中文字幕av成人在线电影| 18禁裸乳无遮挡免费网站照片| 国产黄片美女视频| 最近中文字幕高清免费大全6| 淫妇啪啪啪对白视频| 久久国产乱子免费精品| 午夜精品在线福利| 九九在线视频观看精品| 欧美性猛交黑人性爽| 校园春色视频在线观看| 国产一区二区三区av在线 | 欧美一区二区亚洲| 精品日产1卡2卡| 一级a爱片免费观看的视频| а√天堂www在线а√下载| 美女被艹到高潮喷水动态| 国产精品av视频在线免费观看| 女的被弄到高潮叫床怎么办| 国产午夜精品久久久久久一区二区三区 | 欧美成人a在线观看| 国产伦在线观看视频一区| avwww免费| 麻豆精品久久久久久蜜桃| 色5月婷婷丁香| 亚洲一级一片aⅴ在线观看| aaaaa片日本免费| 久久鲁丝午夜福利片| 91久久精品国产一区二区三区| 观看免费一级毛片| 搡老岳熟女国产| 在线看三级毛片| 国产精品综合久久久久久久免费| 一级毛片我不卡| 精品福利观看| 偷拍熟女少妇极品色| 赤兔流量卡办理| 哪里可以看免费的av片| 直男gayav资源| 精品人妻视频免费看| 韩国av在线不卡| 日日啪夜夜撸| 国产真实乱freesex| 一级黄色大片毛片| 精品久久久久久久末码| 搞女人的毛片| 欧美一区二区国产精品久久精品| av在线老鸭窝| 一级a爱片免费观看的视频| 国产精品嫩草影院av在线观看| 麻豆乱淫一区二区| 久久热精品热| 亚洲婷婷狠狠爱综合网| 在线免费观看的www视频| 国产亚洲91精品色在线| 啦啦啦啦在线视频资源| 1000部很黄的大片| 国产黄a三级三级三级人| 观看免费一级毛片| 精品国产三级普通话版| 一级av片app| 午夜福利在线在线| 成人av在线播放网站| 国产极品精品免费视频能看的| 久久久久精品国产欧美久久久| 91av网一区二区| 久久99热6这里只有精品| 97碰自拍视频| 一进一出抽搐动态| 日韩欧美 国产精品| 春色校园在线视频观看| 国产视频一区二区在线看| 毛片一级片免费看久久久久| 亚洲欧美日韩东京热| 最近中文字幕高清免费大全6| 中国美白少妇内射xxxbb| 久久久久久久午夜电影| 看片在线看免费视频| 综合色丁香网| 你懂的网址亚洲精品在线观看 | 露出奶头的视频| 少妇高潮的动态图| 国产淫片久久久久久久久| av中文乱码字幕在线| 亚洲不卡免费看| 变态另类丝袜制服| 不卡视频在线观看欧美| 亚洲四区av| 久久人人精品亚洲av| 亚洲图色成人| 色综合站精品国产| 国产白丝娇喘喷水9色精品| 麻豆国产av国片精品| 日韩欧美精品v在线| 中文字幕人妻熟人妻熟丝袜美| 色哟哟·www| 午夜日韩欧美国产| 18禁在线播放成人免费| 搞女人的毛片| 亚洲av熟女| 人人妻人人看人人澡| 久久精品91蜜桃| 欧美日韩在线观看h| 欧美中文日本在线观看视频| 最新中文字幕久久久久| 狠狠狠狠99中文字幕| 又黄又爽又免费观看的视频| 久久韩国三级中文字幕| 女同久久另类99精品国产91| 春色校园在线视频观看| 欧美三级亚洲精品| 在线免费十八禁| 欧美+亚洲+日韩+国产| 欧美xxxx性猛交bbbb| 日韩av不卡免费在线播放| 午夜激情福利司机影院| 人人妻,人人澡人人爽秒播| 国产亚洲欧美98| 村上凉子中文字幕在线| 蜜桃亚洲精品一区二区三区| 亚洲成人中文字幕在线播放| 非洲黑人性xxxx精品又粗又长| 欧美又色又爽又黄视频| 如何舔出高潮| 高清日韩中文字幕在线| 久久久精品大字幕| 天堂动漫精品| 大又大粗又爽又黄少妇毛片口| 床上黄色一级片| 午夜视频国产福利| 淫秽高清视频在线观看| 国产亚洲精品综合一区在线观看| a级毛片a级免费在线| 亚洲性久久影院| 一个人免费在线观看电影| 国产精品乱码一区二三区的特点| 69av精品久久久久久| 波野结衣二区三区在线| 性色avwww在线观看| 永久网站在线| 精品一区二区三区视频在线| 久久精品人妻少妇| 色视频www国产| 看非洲黑人一级黄片| 别揉我奶头 嗯啊视频| 久久久久久伊人网av| 午夜精品一区二区三区免费看| 亚洲内射少妇av| 小蜜桃在线观看免费完整版高清| 狂野欧美激情性xxxx在线观看| 亚洲成av人片在线播放无| 欧美一区二区亚洲| 伊人久久精品亚洲午夜| 18禁在线无遮挡免费观看视频 | 国产精品久久视频播放| 蜜臀久久99精品久久宅男| 99久久成人亚洲精品观看| 国产精品精品国产色婷婷| 精品一区二区三区视频在线观看免费| 欧美成人精品欧美一级黄| 麻豆国产97在线/欧美| 99久久九九国产精品国产免费| 亚洲人与动物交配视频| 亚洲av.av天堂| 国产精品99久久久久久久久| 1024手机看黄色片| 久久精品夜夜夜夜夜久久蜜豆| 在线国产一区二区在线| 久久久久九九精品影院| 18禁黄网站禁片免费观看直播| 99久久精品一区二区三区| videossex国产| 色在线成人网| 亚洲三级黄色毛片| 亚洲高清免费不卡视频| 2021天堂中文幕一二区在线观| 国产在视频线在精品| 天美传媒精品一区二区| 午夜福利高清视频| 亚洲天堂国产精品一区在线| 观看美女的网站| av在线蜜桃| 日本在线视频免费播放| 一级毛片久久久久久久久女| 国产高清视频在线播放一区| 直男gayav资源| 日韩欧美免费精品| 一区二区三区高清视频在线| 99久久成人亚洲精品观看| 一个人免费在线观看电影| 日本免费一区二区三区高清不卡| 日韩欧美国产在线观看| 69人妻影院| 久久热精品热| 老女人水多毛片| 亚洲欧美日韩卡通动漫| 亚洲精品在线观看二区| 国产精品一区二区免费欧美| 少妇猛男粗大的猛烈进出视频 | 男女之事视频高清在线观看| 悠悠久久av| 久久久久久久久久久丰满| 午夜亚洲福利在线播放| 在线观看av片永久免费下载| 国产综合懂色| 99国产精品一区二区蜜桃av| 国产精品久久久久久亚洲av鲁大| 淫妇啪啪啪对白视频| 嫩草影院入口| 精品人妻视频免费看| 波多野结衣巨乳人妻| 久久九九热精品免费| 亚洲美女视频黄频| 国产伦一二天堂av在线观看|