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

    基于ITS和trnl-trnF序列的新疆扁桃種質(zhì)資源親緣關(guān)系研究

    2017-03-30 09:41:13徐崇志
    關(guān)鍵詞:扁桃親緣葉綠體

    朱 玲,徐崇志,張 銳,周 燕,高 山

    (1.塔里木大學(xué) 植物科學(xué)學(xué)院,新疆阿拉爾 843300;2.塔里木盆地生物資源保護(hù)利用兵團(tuán)重點(diǎn)實(shí)驗(yàn)室,新疆阿拉爾 843300;3.塔里木大學(xué) 生命科學(xué)學(xué)院,新疆阿拉爾 843300)

    基于ITS和trnl-trnF序列的新疆扁桃種質(zhì)資源親緣關(guān)系研究

    朱 玲1,2,徐崇志1,2,張 銳2,3,周 燕3,高 山1

    (1.塔里木大學(xué) 植物科學(xué)學(xué)院,新疆阿拉爾 843300;2.塔里木盆地生物資源保護(hù)利用兵團(tuán)重點(diǎn)實(shí)驗(yàn)室,新疆阿拉爾 843300;3.塔里木大學(xué) 生命科學(xué)學(xué)院,新疆阿拉爾 843300)

    為探究新疆扁桃種質(zhì)的親緣關(guān)系,以新疆‘紙皮’與‘阿曼尼沙’等24個(gè)扁桃種質(zhì)為材料,對(duì)其葉片及果實(shí)的植物學(xué)性狀進(jìn)行對(duì)比分析。提取總DNA后分別對(duì)ITS與trnl-trnF序列進(jìn)行PCR擴(kuò)增、純化及測序,將所測序列進(jìn)行拼接及同源性比對(duì)。利用MEGA5軟件構(gòu)建基于ITS和trnl-trnF序列的不同扁桃種質(zhì)間的系統(tǒng)發(fā)育樹。通過植物學(xué)性狀對(duì)比發(fā)現(xiàn),‘矮豐’的葉形指數(shù)最高為4.12,‘葉爾羌’的出仁率為72.9%,在供試24個(gè)扁桃種質(zhì)中最高。經(jīng)相關(guān)生物信息學(xué)軟件分析表明,ITS序列長度為608~610 bp,包括42個(gè)變異位點(diǎn)和9個(gè)簡約信息位點(diǎn),G+C含量為61.16%~61.90%。trnl-trnF序列長度為915~933 bp,其變異位點(diǎn)和簡約信息位點(diǎn)豐富,分別為457個(gè)和332個(gè),分別占總長的48.8%和35.7%,G+C含量較低,為32.46%~32.94%?;贗TS和trnl-trnF序列的系統(tǒng)發(fā)育樹表明,ITS序列間的遺傳距離為0.000~0.008,trnl-trnF序列間的遺傳距離為0.000~0.400,‘雙薄’與‘雙軟’具有相同的遺傳背景。

    扁桃;ITS;trnl-trnF;植物學(xué)性狀;系統(tǒng)發(fā)育

    扁桃(AmygdaluscommunisL.)屬薔薇科(Rosaceae)、李亞科(Prunoideae)、桃屬(AmygdalusL.)、扁桃亞屬(Amygdalussubgen.)樹種,又名巴旦杏、美國大杏仁、巴旦姆(新疆)。其作為栽培植物,種仁營養(yǎng)價(jià)值很高,并具有藥理特性,是世界著名的干果和木本油料樹種[1]。扁桃原產(chǎn)于中亞、西亞和非洲北部山區(qū),世界上共約有52~53個(gè)種[2]。中國新疆也是扁桃原產(chǎn)地之一,在新疆天山山區(qū)有野生扁桃林分布,主栽區(qū)在新疆南部喀什地區(qū)的莎車縣、英吉沙縣、疏附縣、疏勤縣、葉城縣、澤普縣和喀什市[3-4]。扁桃的栽培歷史悠久,至今約有6 000 a,中國的扁桃栽培歷史已逾1 300 a,由于長期的自然演化和自播繁衍,形成豐富的資源[5]。中國扁桃共有6個(gè)種,即普通扁桃(A.communisL.)、西康扁桃(A.tanguticaKorsh.)、長柄扁桃(A.pedunculataPall.)、野扁桃(A.ledebourianaSchleche.)、蒙古扁桃(A.mongolicaRicker)和榆葉梅[A.triloba(Lindl.) Ricker]及10個(gè)變種[2-3],本研究供試24份扁桃種質(zhì)均屬于普通扁桃(A.communisL.)。新疆喀什地區(qū)的莎車和英吉沙2個(gè)縣,于20世紀(jì)70年代建立扁桃品種類型匯集圃,大概有34個(gè)扁桃品種類型。近年來,從美國引進(jìn)‘濃帕爾’、‘布特’和‘卡買爾’等[6]11個(gè)扁桃品種。目前新疆本地扁桃種質(zhì)資源混亂。因此,如何有效地進(jìn)行種質(zhì)資源的分子鑒定就成為現(xiàn)階段扁桃種質(zhì)資源研究的重要目標(biāo)。

    核糖體rDNA ITS(Internal transcribed spacer)序列在基因內(nèi)保守但因基因間存在一定變異,進(jìn)化速率較編碼區(qū)快且便于設(shè)計(jì)通用引物,所以常被用來分析植物屬間及屬下水平的系統(tǒng)發(fā)育關(guān)系[7-8]。葉綠體基因組(Choroplast DNA cpDNA)其遺傳方式為單親遺傳,具有較高的序列保守性且進(jìn)化速率恒定,葉綠體rDNA trnl-trnF序列片段大小適中,不受功能限制,其進(jìn)化速率大于功能編碼區(qū),能提供較多系統(tǒng)學(xué)意義的信息位點(diǎn),因此也常被用于植物屬間及屬下水平的系統(tǒng)發(fā)育分析[9-10]。楊波等[11]對(duì)18個(gè)扁桃種質(zhì)進(jìn)行遺傳多樣性和親緣關(guān)系RAPD分析,結(jié)果表明,供試的18個(gè)扁桃種質(zhì)分為2大類,聚類結(jié)果與種質(zhì)分布結(jié)果基本一致。呂志江等[12]采用ISSR分子標(biāo)記技術(shù),對(duì)中國新疆分布的野扁桃5個(gè)居群共120個(gè)個(gè)體的遺傳多樣性進(jìn)行研究。目前,有關(guān)結(jié)合植物學(xué)性狀和核糖體基因及葉綠體基因來分析扁桃親緣的研究較少。邱蓉等[13]以新疆塔城的野扁桃為試驗(yàn)材料,通過擴(kuò)增其細(xì)胞核ITS序列和葉綠體psbA-trnH序列與GeneBank中下載的矮扁桃的ITS和psbAtrnH 序列進(jìn)行比對(duì),綜合分析認(rèn)為野扁桃和矮扁桃為同一個(gè)種。本研究采用植物學(xué)性狀和核糖體ITS及葉綠體trnl-trnlF序列比對(duì)的方法,對(duì)新疆喀什地區(qū)莎車縣的24個(gè)扁桃種質(zhì)構(gòu)建系統(tǒng)進(jìn)化樹,分析其種質(zhì)間的親緣關(guān)系,以期為今后新疆扁桃種質(zhì)資源的開發(fā)、保護(hù)及利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 植物材料

    本試驗(yàn)收集并分析了‘紙皮’、‘阿曼尼沙’及‘雙薄’等24個(gè)扁桃樣本,均取自新疆喀什地區(qū)莎車縣二林場(扁桃種質(zhì)資源圃)。8月初采樣,隨機(jī)選取當(dāng)年新梢頂端健康成熟葉片及果實(shí),部分葉片在液氮中冷凍保存,備用。

    1.2 植物學(xué)性狀測定

    清洗收集的葉片及果實(shí),將其晾干。利用數(shù)顯游標(biāo)卡尺(精度為0.01 cm)、直尺和電子天平(精度為0.01 g)對(duì)扁桃葉片及果實(shí)的主要植物學(xué)性狀進(jìn)行測定。利用DPS 6.55軟件對(duì)所得數(shù)據(jù)進(jìn)行顯著性分析,對(duì)比不同扁桃種質(zhì)間的植物學(xué)性狀差異。

    1.3 DNA提取

    采用CTAB法[14]并加以改良,提取扁桃葉片總DNA,提取的質(zhì)量用8 g/L瓊脂糖凝膠電泳及NanoDrop2000超微量分光光度計(jì)進(jìn)行檢測,置于-20 ℃的冰箱存放,備用。

    1.4 PCR 擴(kuò)增及基因測序

    PCR擴(kuò)增引物及反應(yīng)程序如表1所示。PCR的反應(yīng)體系為20 μL:MgCl21.5 μL,10×buffer 2.0 μL,10 mmol/L dNTP 0.4 μL,10 μmol/L正反向引物各0.4 μL,滅菌ddH2O 13.1 μL,Taq酶0.2 μL,DNA 模板2.0 μL。PCR擴(kuò)增產(chǎn)物先在10 g/L瓊脂糖凝膠上檢測,用溴化乙錠染色,在紫外燈下觀察并拍照查看有無擴(kuò)增產(chǎn)物,然后將PCR產(chǎn)物送到上海生工生物工程技術(shù)服務(wù)有限公司進(jìn)行測序。

    表1 PCR擴(kuò)增所用引物及反應(yīng)程序Table 1 Primers and protocols of PCR

    1.5 序列比對(duì)分析及聚類

    DNA序列的拼接應(yīng)用SeqMan ( DNASTAR Inc.,Madison,Wisconsin,USA) 軟件,用EditSeq( DNASTAR Inc.,Madison,Wisconsin,USA)軟件對(duì)拼接好的序列進(jìn)行分析,利用MEGA 5軟件進(jìn)行多重序列的比對(duì)并對(duì)ITS和trnl-trnlF序列分別構(gòu)建UPGMA系統(tǒng)進(jìn)化樹。

    2 結(jié)果與分析

    2.1 主要植物學(xué)性狀比較

    對(duì)扁桃葉片及果實(shí)的植物性狀進(jìn)行比較(表2和表3),便于對(duì)不同扁桃種質(zhì)的異同點(diǎn)進(jìn)行分析。由表2可以看出,大部分扁桃種質(zhì)的葉形為披針形;‘紙皮’為闊橢圓形;‘葉爾羌’和‘寒豐’為長橢圓形。同時(shí),‘寒豐’和‘巴旦王’的葉長、葉寬及葉柄長度與‘紙皮’和‘阿曼尼沙’等多數(shù)扁桃種質(zhì)相比具有顯著差異。‘矮豐’的葉形指數(shù)最高,‘雙果’和‘苦巴旦’次之,‘鷹嘴’的葉形指數(shù)最低。由表3可以看出,‘白薄’與‘尖軟’的果實(shí)長度較其他種質(zhì)略長,‘巴旦桃’內(nèi)果皮厚度最大,‘白薄’次之,‘葉爾羌’的內(nèi)果皮厚度最小。在本試驗(yàn)24個(gè)扁桃種質(zhì)中,多數(shù)種質(zhì)內(nèi)果皮呈淡褐色,‘鷹嘴’和‘晚豐’等少數(shù)種質(zhì)則呈現(xiàn)褐色,而‘白薄’呈黃白色,其內(nèi)果皮顏色較其他種質(zhì)偏黃,‘小石子’則呈灰褐色?!~爾羌’出仁率最高為72.9%,‘紙皮’出仁率為68.5%排列第2位,‘巴旦桃’的出仁率僅有19.3%,為出仁率最低的種質(zhì)。

    表2 不同扁桃種質(zhì)葉片比較分析Table 2 Comparative analysis of germplasm leaves of different species in Amygdalus communis

    注:同列中不同字母表示樣本之間該指標(biāo)存在極顯著差異(P<0.01)。下表同。

    Note:Differences,different letters in same column indicate significant difference among samples(P<0.01).The same as below.

    表3 不同扁桃種質(zhì)果實(shí)比較分析Table 3 Comparative analysis of different germplasm fruits in Amygdalus communis

    2.2 核糖體ITS序列比對(duì)及系統(tǒng)發(fā)育樹的構(gòu)建

    將試驗(yàn)測得‘紙皮’、‘阿曼尼沙’和‘雙薄’等24個(gè)扁桃樣本的ITS序列用EditSeq軟件進(jìn)行分析。結(jié)果表明,ITS序列長度為608~610 bp,G+C含量為61.16%~61.9%。用MEGA5進(jìn)行多重序列比對(duì),對(duì)位后長度為611 bp,包括42個(gè)變異位點(diǎn)和9個(gè)簡約信息位點(diǎn)。重復(fù)次數(shù)為1 000 次,構(gòu)建UPGMA樹(圖1),結(jié)果顯示,在遺傳距離為0.004處,‘巴旦王’、‘巴旦桃’和‘尖軟’親緣關(guān)系最近,同時(shí)這3個(gè)種質(zhì)又可與‘鷹嘴’、‘白雙’和‘矮豐’等共14個(gè)種質(zhì)聚為一枝,表明其親緣關(guān)系較近?!妆 汀庾臁蹫橐恢Γ⒂凇喟偷瘶?gòu)成姐妹群?!p薄’和‘晚豐’聚為一枝,并與‘小石子’構(gòu)成姐妹群。而‘小雙’、‘皮斯特’‘石子’和‘紙皮’各自為一枝,表明這4個(gè)扁桃種質(zhì)與其他種質(zhì)間的親緣關(guān)系較遠(yuǎn)。

    圖1 MEGA5軟件構(gòu)建基于ITS序列的扁桃UPGMA樹Fig.1 UPGMA tree of Amygdalus communis based on ITS sequences by MEGA5

    2.3 葉綠體trnl-F序列同源性分析及系統(tǒng)發(fā)育樹的構(gòu)建

    將試驗(yàn)測得‘紙皮’ 、‘阿曼尼沙’和‘雙薄’等24個(gè)扁桃樣本的trnl-trnF序列拼接后用EditSeq軟件分析表明,trnl-trnF序列長度為915~933 bp,G+C含量為32.46%~32.94%。用MEGA5進(jìn)行多重序列比對(duì),對(duì)位后長度為937 bp,其中包含豐富的變異位點(diǎn)和簡約信息位點(diǎn),分別為457個(gè)和332個(gè),分別占總長的48.8%和35.7%。并構(gòu)建UPGMA樹(圖2),重復(fù)1 000次。聚類結(jié)果表明,24個(gè)扁桃種質(zhì)最終被聚為兩枝,‘雙薄’與‘雙軟’聚為一枝,表明其親緣關(guān)系最近,而其他22個(gè)則聚為一大枝,表明‘矮豐’、‘巴旦桃’和‘鷹嘴’等22個(gè)扁桃種質(zhì)之間的親緣關(guān)系較近,且與‘雙薄’與‘雙軟’的親緣關(guān)系較遠(yuǎn)。

    3 討 論

    3.1 主要植物學(xué)性狀差異的分析

    從植物學(xué)性狀對(duì)比結(jié)果來看,供試24個(gè)扁桃種質(zhì)在葉片外形指標(biāo)、內(nèi)果皮外形指標(biāo)及出仁率方面存在一定差異?!埰ぁ娜~形為闊橢圓形,‘葉爾羌’和‘寒豐’為長橢圓形,這3 個(gè)扁桃種質(zhì)在葉形上與其他種質(zhì)間存在明顯物種差異性?!⒙嵘场绕溆?1個(gè)扁桃種質(zhì)的葉形均為披針形,其中‘鷹嘴’和‘矮豐’為長披針形,‘雙果’和‘白薄’為長橢披針形,‘公巴旦’則為卵披針形。在24份扁桃種質(zhì)中,‘矮豐’的葉形指數(shù)最高,達(dá)到4.12,‘雙果’和‘苦巴旦’次之,分別為4.06和4.02,‘鷹嘴’的葉形指數(shù)最低,為2.27,說明‘矮豐’、‘雙果’和‘苦巴旦’的葉形是供試種質(zhì)中最為細(xì)長的3個(gè)種質(zhì),而‘鷹嘴’的葉形則既寬又扁。

    圖2 MEGA5軟件構(gòu)建基于trnl-trnF序列的扁桃UPGMA樹Fig.2 UPGMA tree of Amygdalus communis based on trnl-trnF sequences by MEGA5

    根據(jù)果實(shí)外部性狀比較分析可知,在供試24個(gè)扁桃種質(zhì)中,除‘鷹嘴’‘晚豐’‘苦巴旦’‘皮斯特’‘巴旦桃’和‘尖軟’的內(nèi)果皮顏色為褐色,多數(shù)種質(zhì)內(nèi)果皮呈淡褐色。‘白薄’呈黃白色,其內(nèi)果皮顏色較其他種質(zhì)偏黃,‘小石子’呈灰褐色。‘白薄’與‘尖軟’的果實(shí)長度與其他扁桃種質(zhì)間存在顯著差異,較其他種質(zhì)略長?!~爾羌’出仁率最高為72.9%,‘紙皮’出仁率為68.5%,排列第2位,其內(nèi)果皮厚度分別為1.03 mm與1.36 mm,是供試扁桃種質(zhì)中內(nèi)果皮厚度最小的2個(gè)種質(zhì),這表明‘葉爾羌’和‘紙皮’的商品性最高?!偷┨摇某鋈事蕛H有19.3%,為供試扁桃種質(zhì)中出仁率最低的種質(zhì),并且其內(nèi)果皮厚度也為供試扁桃種質(zhì)中最厚的,說明‘巴旦桃’較難取食,其商品性最低。

    3.2 ITS序列及trnl-trnF的分子進(jìn)化比較分析

    在目前對(duì)扁桃種質(zhì)資源及遺傳多樣性的研究中,主要以傳統(tǒng)形態(tài)學(xué)特征及RAPD與ISSR等分子標(biāo)記為研究手段。通過核糖體rDNA ITS序列與葉綠體cpDNA trnl-trnF序列基因測序、序列比對(duì)及構(gòu)建系統(tǒng)發(fā)育樹來探討果樹植物親緣關(guān)系的研究還較少。王化坤等[15]通過測定6個(gè)核果類果樹ITS序列,用最大簡約法構(gòu)建桃、李、梅、杏、櫻桃的系統(tǒng)發(fā)育樹。結(jié)果發(fā)現(xiàn),核果類果樹ITS1和ITS2的分子進(jìn)化速率不同和信息量存在差異,桃演化順序?yàn)榘偷┬?山桃-普通桃、新疆桃,并支持將核果類果樹分成4個(gè)亞屬。本研究依據(jù)核基因ITS序列和葉綠體trnl-trnF序列構(gòu)建24個(gè)新疆扁桃不同種質(zhì)的系統(tǒng)發(fā)育樹。供試的24個(gè)扁桃種質(zhì)的ITS及trnl-trnF序列經(jīng)相關(guān)生物信息學(xué)軟件分析,結(jié)果顯示ITS序列長度為608~610 bp,對(duì)位后長度為611 bp,包括42個(gè)變異位點(diǎn)和9個(gè)簡約信息位點(diǎn),G+C含量為61.16%~61.90%。trnl-trnF序列長度為915~933 bp,對(duì)位后長度為937 bp,其中包含豐富的變異位點(diǎn)和簡約信息位點(diǎn),分別為457個(gè)和332個(gè),分別占總長的48.8%和35.7%,G+C含量較低,為32.46%~32.94%。ITS序列間的遺傳距離為0.000~0.008,trnl-trnF序列間的遺傳距離為0.000~0.400,ITS序列總體遺傳距離較trnl-trnF序列小。在本研究中扁桃ITS序列遺傳距離較小,出現(xiàn)的遺傳變異位點(diǎn)較少,信息含量偏低,說明24種扁桃起源較為相似。相比較而言,trnl-trnF序列則具有豐富的變異位點(diǎn)和簡約信息位點(diǎn),說明trnl-trnF序列在扁桃中存在豐富的遺傳變異。而趙海光等[16]及陳仁芳等[17]研究表明,在種屬進(jìn)化關(guān)系中ITS序列變異位點(diǎn)高于trnL-trnF,與本研究結(jié)論不一致,原因可能在于本研究是同一種內(nèi)不同地域間扁桃種質(zhì)的親緣關(guān)系,體現(xiàn)了核糖體ITS序列和葉綠體trnL-trnF序列在不同分類水平的鑒別能力。

    核糖體上的基因通常為多拷貝數(shù)、中度重復(fù)序列,1個(gè)重復(fù)單位由5.8S、18S、26S編碼區(qū)以及一些間隔區(qū)組成。真核生物核糖體內(nèi)轉(zhuǎn)錄間隔區(qū)ITS (Internal transcribed spacer) 位于18S和26S基因之間,中部被5.8S分為ITS1區(qū)與ITS2區(qū)。ITS序列既有核苷酸序列的高度變異性又有長度上的保守性,說明這些序列很容易在近緣類群間排序,并且因其豐富的變異性可在種間、屬間及屬下水平解決植物系統(tǒng)發(fā)育問題。研究表明,在被子植物ITS1區(qū)與ITS2區(qū)的序列長度變異很小(均小于300 bp)[8]。在探討一些被子植物種間、屬間及屬下水平的系統(tǒng)發(fā)育時(shí),ITS區(qū)提供的信息對(duì)于被子植物系統(tǒng)發(fā)育研究是很有價(jià)值的,目前已被廣泛地應(yīng)用于被子植物種間及近緣屬間的系統(tǒng)發(fā)育和分類研究中[18-20]。本研究對(duì)于扁桃屬不同扁桃種質(zhì)的ITS序列聚類分析得到,‘巴旦王’、‘巴旦桃’和‘尖軟’聚為一枝,表明其親緣關(guān)系最近。同時(shí)這3個(gè)種質(zhì)又可與‘鷹嘴’、‘白雙’和‘矮豐’等共14個(gè)種質(zhì)聚為一枝姐妹群,表明其親緣關(guān)系較近?!妆 汀庾臁蹫橐恢?,并與‘苦巴旦’構(gòu)成姐妹群?!p薄’和‘晚豐’聚為一枝,并與‘小石子’構(gòu)成姐妹群。而‘小雙’‘皮斯特’‘石子’和‘紙皮’各自為一枝,表明這4個(gè)扁桃種質(zhì)與其他種質(zhì)間的親緣關(guān)系較遠(yuǎn)。

    葉綠體DNA(chloroplast DNA,cpDNA)最早是在藻類中發(fā)現(xiàn)的,葉綠體基因序列相當(dāng)保守,其非編碼區(qū)DNA進(jìn)化速率普遍高于功能編碼區(qū)。cpDNA遺傳方式為母系遺傳,遺傳性狀穩(wěn)定,有獨(dú)立的進(jìn)化路線,不依賴于其他任何數(shù)據(jù)即可構(gòu)建分子系統(tǒng)進(jìn)化樹。葉綠體基因trnL-trnF(包括trnL內(nèi)含子和trnL-F基因間隔區(qū)),在進(jìn)化上具有選擇壓力小和進(jìn)化速率較快的特點(diǎn),較為適合種間的系統(tǒng)發(fā)育和遺傳分析[21],目前cpDNA trnL-trnF已被廣泛用于植物屬間、種間關(guān)系的鑒別。如宮霞等[22]利用葡萄trnL-trnF的序列,對(duì)于葡萄科各屬間的分子系統(tǒng)學(xué)進(jìn)行了研究。本試驗(yàn)利用葉綠體基因通用引物對(duì)24種扁桃樣品進(jìn)行擴(kuò)增,測序結(jié)果采用相關(guān)生物信息學(xué)軟件進(jìn)行序列拼接、比對(duì)及系統(tǒng)發(fā)育樹構(gòu)建。根據(jù)葉綠體基因母系遺傳的特點(diǎn),可以推測‘矮豐’、‘巴旦王’和‘鷹嘴’等22個(gè)扁桃種質(zhì)具有相同的遺傳背景,而‘雙薄’與‘雙軟’則可能為同一母本。

    綜上所述,從植物學(xué)性狀對(duì)比結(jié)果來看,供試24個(gè)扁桃種質(zhì)在葉片外形指標(biāo)、內(nèi)果皮外形指標(biāo)及出仁率方面存在差異。利用ITS序列及trnL-trnF序列對(duì)供試24個(gè)扁桃種質(zhì)構(gòu)建系統(tǒng)發(fā)育樹,聚類結(jié)果明顯不同。利用ITS序列所構(gòu)建的系統(tǒng)發(fā)育樹可聚為多枝,并且‘雙薄’與‘雙軟’種質(zhì)親緣關(guān)系較遠(yuǎn)。而基于葉綠體trnL-trnF序列構(gòu)建的系統(tǒng)發(fā)育樹中‘雙薄’與‘雙軟’能夠聚為一枝,說明它們具有相同的遺傳背景。有研究表明被子植物中葉綠體DNA的進(jìn)化速率比核糖體DNA低[23],這可能是導(dǎo)致本研究中ITS序列和trnL-trnF序列構(gòu)建的系統(tǒng)發(fā)育樹有差異的原因。

    Reference:

    [1] 成健紅,譚敦炎,艾爾肯,等.巴旦杏花物候?qū)W與形態(tài)學(xué)研究[J].西北植物學(xué)報(bào),2001,21(2):339-344.

    CHENG J H,TAN D Y,Aierken,etal.Flowering phenology and morphology on almond[J].ActaBotanyBoreali-OccidentalisSinica,2001,21(2):339-344(in Chinese with English abstract).

    [2] 曾 斌,李 疆,羅淑萍,等.扁桃屬植物種質(zhì)資源鑒定的SSR分析研究[J].新疆農(nóng)業(yè)科學(xué),2009,46 (1):18-22.

    ZENG B,LI J,LUO SH P,etal.Identification of genetic relationship of amygdalus plants by SSR[J].XinjiangAgriculturalSciences,2009,46 (1):18-22(in Chinese with English abstract).

    [3] 汪祖華,莊恩及.中國果樹志·桃卷[M].北京:中國林業(yè)出版社,2001:81-85.

    WANG Z H,ZHUANG E J.Chinese Fruit Trees·Peach Roll [M].Beijing:China Forestry Publishing House,2001:81-85(in Chinese).

    [4] 蘭彥平,吐拉克孜,郭文英,等.巴旦杏的研究現(xiàn)狀及開發(fā)利用前景[J].林業(yè)科學(xué)研究,2004,17(5):674-679.

    LAN Y P,Tulakezi,GUO W Y,etal.The present situation and prospect research of development and utilization on almond[J].ForestResearch,2004,17(5):674-679(in Chinese).

    [5] 王 琳,姜 喜,李志軍.中國扁桃種質(zhì)資源研究進(jìn)展[J].北方園藝,2013(20):178-181.

    WANG L,JIANG X,LI ZH J.Research progress of almond germplasm resource in China [J].NorthernHorticulture,2013(20):178-181(in Chinese with English abstract).

    [6] 梁艷霞,王占和,張志玲,等.我國扁桃發(fā)展現(xiàn)狀、問題及發(fā)展前景[J].現(xiàn)代園藝,2015(6):9.

    LIANG Y X,WANG ZH H,ZHANG ZH L,etal.The development present situation,problems and development prospects of almond in our country[J].XiandaiHorticulture,2015(6):9(in Chinese).

    [7] BRUCE G,BALDWIN.Phylogenetic utility of the internal transcribed spacers of nuclear ribosomal DNA in plants:An example from the compositae[J].MolecularPhylogenetisandEvolution,1992,1(1):3-16.

    [8] 趙志禮,徐珞珊,董 輝,等.核糖體DNA ITS區(qū)序列在植物分子系統(tǒng)學(xué)研究中的價(jià)值[J].植物資源與環(huán)境學(xué)報(bào),2000,9(2):50-54.

    ZHAO ZH L,XU L SH,DONG H,etal.Evaluation of ITS sequence of nrDNA in plant molecular systematics[J].JournalofPlantResourcesandEnvironment,2000,9(2):50-54(in Chinese with English abstract).

    [9] BAYER B,STARR J R.Tribal phylogeny of the Asteraceae based on two non-coding chloroplast sequences,the trnL intron and trnL/trnF intergenic spacer[J].AnnalsoftheMissouriBotanicalGarden,1998,85:242-256.

    [10] NISHIKAWA T,SALOMON B,KOMATSUDA T,etal.Molecular phylogeny of the genus Hordeum using three chloroplast DNA sequences[J].Genome,2002,45:1157-1166.

    [11] 楊 波,龔 鵬,車玉紅,等.新疆引進(jìn)扁桃品種與當(dāng)?shù)仄贩N間種質(zhì)資源的RAPD分析[J].新疆農(nóng)業(yè)科學(xué),2009,46(1):91-94.

    YANG B,GONG P,CHE Y H,etal.RAPD analysis of local and introduct-ed almond varieties in Xinjiang[J].XinjiangAgriculturalSciences,2009,46(1):91-94(in Chinese with English abstract).

    [12] 呂志江,李 疆,吾買爾夏提塔漢,等.新疆野扁桃種質(zhì)資源遺傳多樣性的ISSR分析[J].果樹學(xué)報(bào),2010,27(6):918-923.

    Lü ZH J,LI J,Wumaier Xiatitahan,etal.ISSR analysis for genetic diversity ofAmygdalusledebourianagermplasm from Xingjiang,China[J].JournalofFruitScience,2010,27(6):918-923(in Chinese with English abstract).

    [13] 邱 蓉,程中平,王章利,等.利用核ITS序列和葉綠體psbA-trnH序列探討野扁桃和矮扁桃分類學(xué)關(guān)系[J].果樹學(xué)報(bào),2012,29(3):387-392.

    QIU R,CHENG ZH P,WANG ZH L,etal.Study on taxonomic relation ofAmygdalusledebourianaandA.nanabased on ITS and psbA-trnH sequences[J].JournalofFruitScience,2012,29(3):387-392(in Chinese with English abstract).

    [14] WEISING K H,WOLFF K,MEYER W.DNA Fingerprinting in Plants and Fungi[M].Boca Raton,Florida,US:CRC Press,1995:50-54.

    [15] 王化坤,陶建敏,渠慎春,等.核果類果樹ITS序列分子進(jìn)化及系統(tǒng)發(fā)育關(guān)系研究[J].園藝學(xué)報(bào),2010,37(3):363-374.

    WANG H K,TAO J M,QU SH CH,etal.Molecular evolution and phylogeny of stone fruit trees based on sequences of the internal transcribed spacers(ITS)of nuclear ribosomal DNA [J].ActaHorticulturaeSinica,2010,37(3):363-374(in Chinese with English abstract).

    [16] 趙海光,周建建,曹珊珊,等.基于ITS和trnL-F序列堿基差異的繁縷及其近緣種的親緣關(guān)系分析[J].植物資源與環(huán)境學(xué)報(bào),2009,18(1):1-5.

    ZHAO H G,ZHOU J J,CAO SH SH,etal.Analysis of interspecific relationship amongStellariamediaand its related species based on ITS and trnL-F sequence differences[J].JournalofPlantResourcesandEnvironment,2009,18(1):1-5(in Chinese with English abstract).

    [17] 陳仁芳,張 澤,唐 洲,等.桑屬ITS、trnL-F、rps16序列與進(jìn)化分析[J].中國農(nóng)業(yè)科學(xué),2011,44(8):1553-1561.

    CHEN R F,ZHANG Z,TANG ZH,etal.Morus ITS,trnL-F,rps16 sequence and phylogenetic analysis of mulberry resources [J].ScientiaAgriculturaSinica,2011,44(8):1553-1561.(in Chinese with English abstract).

    [18] 李云霞,蘭芙蓉,常朝陽,等.基于ITS和trnL-F序列的青藏高原棘豆屬植物分子系統(tǒng)學(xué)研究[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,39(11):187-192.

    LI Y X,LAN F R,CHANG ZH Y,etal.Molecular phylogeny ofOxytroposDC.of Qinghai-Tibetan Plateau by ITS and trnL-F sequences[J].JournalofNorthwestA&FUniversity(NaturalScienceEdition),2011,39(11):187-192(in Chinese with English abstract).

    [19] 李曉芹,馬 婧,李 敏,等.基于ITS和trnL-F以及rps4基因片段的粗疣蘚屬分類學(xué)地位研究[J].西北植物學(xué)報(bào),2014,34(1):72-76.

    LI X Q,MA J,LI M,etal.Taxonomic status ofFauriellabesch.inferred from its,trnl-f and rps4 gene regions[J].ActaBotanyBoreali-OccidentalisSinica,2014,34(1):72-76(in Chinese with English abstract).

    [20] 徐玲玲,李同建,張美云,等.基于核ITS與葉綠體trnL-F序列分析12種紫金牛屬植物的種間關(guān)系與變異[J].園藝學(xué)報(bào),2009,36(10):1531-1537.

    XU L L,LI T J,ZHANG M Y,etal.Interspecific relationships and variation of 12 species in Ardisia Sw.(Myrsinaceae) based on ITS and trnL-F Data Sets [J].ActaHorticulturaeSinica,2009,36(10):1531-1537(in Chinese with English abstract).

    [21] 曹慶芹,甄 貞,姜 潔,等.葉綠體DNA分析技術(shù)及其在栗屬植物中的應(yīng)用[J].果樹學(xué)報(bào),2008,25(3):396-399.

    CAO Q Q,ZHEN ZH,JIANG J,etal.Chloroplast DNA analysis technology and its application inCastanea[J].JournalofFruitScience,2008,25(3):396-399(in Chinese with English abstract).

    [22] 宮 霞,王煥沖,楊 勇,等.基于TrnL-F序列的葡萄科各屬分子特征及系統(tǒng)發(fā)育學(xué)研究[J].果樹學(xué)報(bào),2013,30(1):28-36.

    GONG X,WANG H CH,YANG Y,etal.Molecular characteristics and phylogenitic analysis ofVitaceaebased on TrnL-F spacer sequences.[J].JournalofFruitScience,2013,30(1):28-36(in Chinese with English abstract).

    [23] 唐 萍,彭 程.葉綠體基因組進(jìn)化的速率和方式[J].生物學(xué)通報(bào),2010,45(6):8-10.

    TANG P,PENG CH.The rate and patterns of chloroplast DNA evolution[J].BulletinofBiology,2010,45(6):8-10(in Chinese ).

    (責(zé)任編輯:潘學(xué)燕 Responsible editor:PAN Xueyan)

    Study on Genetic Relationship ofAmygdaluscommunisPopulations Based on ITS and trnl-trnF Sequences

    ZHU Ling1,2,XU Chongzhi1,2,ZHANG Rui2,3,ZHOU Yan3and GAO Shan1

    (1.College of Plant Sciences,Tarim University,Alar Xinjiang 843300,China;2.Key Laboratory of Protection and Utilization of Bological Resources in Tarim Basin of Xinjiang Production and Construction Corps,Alar Xinjiang 843300,China; 3.College of Life Sciences,Tarim University,Alar Xinjiang 843300,China)

    The purpose of the paper is to explore the genetic relationship and germplasm resources ofAmygdaluscommunisin Xinjiang for providing a scientific basis for protection and utilization in the future.With ‘Zhipi’ and ‘Amannisha’ and the other 24Amygdaluscommunisvarieties in Xinjiang as materials,and then we compared and analysed its botanical characters in leaves and fruits.After the total DNA was extracted,the ITS and trnl-trnF sequences were respectively amplified and sequenced.These sequences were joined together,then the homology of the sequences were compared.Using MEGA5 software to build the phylogenetic tree of differentAmygdaluscommunisvarieties based on ITS trnl-trnF sequences.Compared with botanical properties,the results found that ‘Aifeng’ had the highest leaf shape index of 4.12,the highest kernel percentage was 72.9% in ‘Yeerqiang’,it was the highest in the selected 24Amygdaluscommunisgermplasm resources.The analysis by related bioinformatics software showed that ITS was in 608-610 bp sequence length,including 42 mutation loci and 9 simple information sites,G+C content was at 61.16%-61.90%.Sequence length of trnl-trnF was in 915-933 bp and it had abundant variation and simple information sites which was 457 and 332 respectively,accounting for 48.8% and 35.7% of the total length respectively,its G+C content was low at 32.46%-32.94%.ITS and trnl-trnF sequences and,phylogenetic tree showed that the genetic distance between ITS sequence was between 0.000-0.008,genetic distance between trnl-trnF sequences was between 0.000-0.400,‘Shuangbao’ and ‘Shuangruan’ had the same genetic background.

    Amygdaluscommunis; ITS; trnl-trnF; Botanical characters; Phylogenetic analysis

    ZHU Ling,female,master student.Research area:fruit genetics and breeding.E-mail:bettyzl2015@163.com

    GAO Shan,male,master,associate professor.Research area:cultivation techniques of agroforestry.E-mail:1227081916@qq.com

    日期:2017-03-03

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

    2016-02-20

    2016-03-20

    新疆生產(chǎn)建設(shè)兵團(tuán)塔里木盆地生物資源保護(hù)利用重點(diǎn)實(shí)驗(yàn)室開放基金課題(BRYB1007)。

    朱 玲,女,在讀碩士,研究方向?yàn)楣麡溥z傳育種。E-mail:bettyzl2015@163.com

    高 山,男,碩士,副教授,主要從事林農(nóng)間作高效栽培技術(shù)研究。E-mail:1227081916@qq.com

    S662.9

    A

    1004-1389(2017)03-0412-08

    Received 2016-02-20 Returned 2016-03-20

    Foundation item Open Fund Project of Xinjiang Production and Construction Corps Key Laboratory of Protection and Utilization of Bio-resources in Tarim Basin(No.BRYB1007).

    猜你喜歡
    扁桃親緣葉綠體
    谷子近緣野生種的親緣關(guān)系及其利用研究
    中國醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    扁桃斑鳩菊揮發(fā)性成分和無機(jī)元素分布
    中成藥(2018年6期)2018-07-11 03:01:20
    小白菜種質(zhì)遺傳多樣性與親緣關(guān)系的SRAP 和SSR分析
    大漠瀚海中的珍稀瀕危野生植物之蒙古扁桃
    脂肪酶Novozyme435手性拆分(R,S)-扁桃酸
    南方紅豆杉葉綠體非編碼序列PCR體系優(yōu)化及引物篩選
    茶樹葉綠體DNA的PCR-RFLP反應(yīng)體系優(yōu)化
    煙草葉綠體密碼子的偏好性及聚類分析
    亚洲第一区二区三区不卡| 国产爽快片一区二区三区| 婷婷成人精品国产| 老鸭窝网址在线观看| 久久人人爽人人片av| 国产精品 国内视频| kizo精华| 咕卡用的链子| 欧美激情高清一区二区三区 | 欧美日本中文国产一区发布| 黄色 视频免费看| 亚洲成色77777| av在线app专区| 在线 av 中文字幕| 午夜免费鲁丝| 9色porny在线观看| 日韩av在线免费看完整版不卡| 成人午夜精彩视频在线观看| 丝袜美足系列| 卡戴珊不雅视频在线播放| 熟女少妇亚洲综合色aaa.| 又大又爽又粗| 亚洲精品美女久久久久99蜜臀 | 下体分泌物呈黄色| 国产高清国产精品国产三级| 成年动漫av网址| 在线观看www视频免费| av国产精品久久久久影院| 久久精品久久久久久噜噜老黄| 国产熟女欧美一区二区| 精品人妻在线不人妻| 悠悠久久av| 夜夜骑夜夜射夜夜干| 日本猛色少妇xxxxx猛交久久| 国产精品成人在线| 国产成人精品在线电影| 日韩 欧美 亚洲 中文字幕| 亚洲三区欧美一区| 午夜日韩欧美国产| 午夜久久久在线观看| 久久99精品国语久久久| 乱人伦中国视频| 高清不卡的av网站| 视频在线观看一区二区三区| 99国产精品免费福利视频| 亚洲精品一区蜜桃| 久久青草综合色| 国产一区亚洲一区在线观看| 超色免费av| 国产成人精品在线电影| 久久精品熟女亚洲av麻豆精品| 精品福利永久在线观看| 久久久精品免费免费高清| 黄色 视频免费看| 日本欧美国产在线视频| 欧美久久黑人一区二区| 七月丁香在线播放| av在线观看视频网站免费| 免费高清在线观看日韩| 人妻 亚洲 视频| 自线自在国产av| 午夜激情av网站| 少妇猛男粗大的猛烈进出视频| 一边亲一边摸免费视频| 中文字幕人妻熟女乱码| 51午夜福利影视在线观看| 高清视频免费观看一区二区| 日本91视频免费播放| 啦啦啦啦在线视频资源| 一区二区三区乱码不卡18| 制服人妻中文乱码| 桃花免费在线播放| 国产黄频视频在线观看| 国产有黄有色有爽视频| 人人妻人人爽人人添夜夜欢视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成人一二三区av| 成年人午夜在线观看视频| av国产精品久久久久影院| 午夜久久久在线观看| 午夜老司机福利片| 丝袜人妻中文字幕| 在线天堂最新版资源| 99国产精品免费福利视频| 精品视频人人做人人爽| 国产深夜福利视频在线观看| 九色亚洲精品在线播放| 亚洲成人手机| 在线观看国产h片| 免费不卡黄色视频| 免费高清在线观看视频在线观看| 久久人妻熟女aⅴ| 不卡视频在线观看欧美| 少妇精品久久久久久久| 欧美日韩亚洲高清精品| 久热爱精品视频在线9| 国产精品久久久久久精品古装| 日韩免费高清中文字幕av| 国产一区二区 视频在线| 亚洲国产av新网站| 亚洲av中文av极速乱| 国产精品久久久久久人妻精品电影 | 亚洲精品aⅴ在线观看| 九草在线视频观看| 一区福利在线观看| 国产深夜福利视频在线观看| 一边摸一边做爽爽视频免费| av电影中文网址| 午夜免费鲁丝| √禁漫天堂资源中文www| 久久鲁丝午夜福利片| 国产成人系列免费观看| 亚洲国产成人一精品久久久| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区久久| 黑人欧美特级aaaaaa片| 又粗又硬又长又爽又黄的视频| 亚洲欧洲国产日韩| 日韩一本色道免费dvd| 精品国产乱码久久久久久男人| 亚洲精品美女久久久久99蜜臀 | 国产免费又黄又爽又色| 久久久久久久国产电影| 侵犯人妻中文字幕一二三四区| 一级毛片我不卡| 日本av手机在线免费观看| 多毛熟女@视频| 色婷婷av一区二区三区视频| 国产精品一国产av| 精品福利永久在线观看| 国产精品.久久久| 老汉色∧v一级毛片| 日本av手机在线免费观看| 亚洲av男天堂| 日本午夜av视频| 永久免费av网站大全| 成人亚洲欧美一区二区av| 国产在线视频一区二区| 最近最新中文字幕免费大全7| 街头女战士在线观看网站| 久久精品国产亚洲av高清一级| www.av在线官网国产| 少妇被粗大的猛进出69影院| 国产亚洲av片在线观看秒播厂| 欧美成人精品欧美一级黄| 亚洲欧美精品综合一区二区三区| 最近手机中文字幕大全| 久久 成人 亚洲| 青青草视频在线视频观看| 夫妻午夜视频| 成人国产av品久久久| 波野结衣二区三区在线| 男女床上黄色一级片免费看| 久热这里只有精品99| 亚洲伊人色综图| 久久久久久久大尺度免费视频| 黄网站色视频无遮挡免费观看| 97精品久久久久久久久久精品| 女性被躁到高潮视频| 久久精品亚洲熟妇少妇任你| 高清不卡的av网站| 日本av手机在线免费观看| 欧美久久黑人一区二区| 国产高清国产精品国产三级| √禁漫天堂资源中文www| 大片电影免费在线观看免费| 国产精品亚洲av一区麻豆 | 性色av一级| 亚洲欧美激情在线| 一级a爱视频在线免费观看| 色婷婷久久久亚洲欧美| 亚洲自偷自拍图片 自拍| 99久久99久久久精品蜜桃| 一级毛片我不卡| 亚洲图色成人| 在线精品无人区一区二区三| 亚洲精品美女久久av网站| 久久久亚洲精品成人影院| 午夜久久久在线观看| 一区二区三区乱码不卡18| 老熟女久久久| 在线精品无人区一区二区三| 欧美亚洲 丝袜 人妻 在线| 少妇人妻久久综合中文| 国产成人av激情在线播放| 久久久久网色| 亚洲少妇的诱惑av| 国产精品香港三级国产av潘金莲 | 国产在线一区二区三区精| 考比视频在线观看| 黄频高清免费视频| 国产黄色免费在线视频| 亚洲精品乱久久久久久| 亚洲精品国产一区二区精华液| 日本爱情动作片www.在线观看| 美国免费a级毛片| 曰老女人黄片| 国产精品蜜桃在线观看| 午夜福利视频在线观看免费| 在线天堂中文资源库| 肉色欧美久久久久久久蜜桃| 欧美日韩视频精品一区| 午夜激情久久久久久久| 咕卡用的链子| 丝瓜视频免费看黄片| 热99久久久久精品小说推荐| 99国产综合亚洲精品| 欧美日韩福利视频一区二区| 免费av中文字幕在线| 精品一区二区三区av网在线观看 | 亚洲综合精品二区| 免费看av在线观看网站| av在线app专区| 尾随美女入室| 天堂中文最新版在线下载| 久久99精品国语久久久| 又粗又硬又长又爽又黄的视频| 丁香六月欧美| 免费黄频网站在线观看国产| 午夜av观看不卡| 亚洲精品国产av成人精品| 亚洲人成网站在线观看播放| 别揉我奶头~嗯~啊~动态视频 | 久久 成人 亚洲| 国产成人系列免费观看| 精品国产乱码久久久久久小说| 最近手机中文字幕大全| 国产在线视频一区二区| 国产精品欧美亚洲77777| 美女视频免费永久观看网站| 天天操日日干夜夜撸| 精品少妇久久久久久888优播| 国产亚洲av片在线观看秒播厂| 男女之事视频高清在线观看 | 亚洲中文av在线| 欧美精品人与动牲交sv欧美| 成人国产av品久久久| 国产精品免费大片| 美国免费a级毛片| 精品亚洲成国产av| 国产成人精品福利久久| 成人毛片60女人毛片免费| 国产97色在线日韩免费| 黄色毛片三级朝国网站| 高清在线视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 视频在线观看一区二区三区| 国产成人av激情在线播放| 丝袜在线中文字幕| 亚洲男人天堂网一区| 色精品久久人妻99蜜桃| 中文精品一卡2卡3卡4更新| 一本大道久久a久久精品| 丁香六月天网| 99国产综合亚洲精品| 国产精品无大码| 亚洲欧洲日产国产| 久久韩国三级中文字幕| 在现免费观看毛片| 精品亚洲乱码少妇综合久久| 亚洲在久久综合| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 亚洲国产av影院在线观看| 国产欧美日韩综合在线一区二区| 国产av国产精品国产| 久久久久国产精品人妻一区二区| tube8黄色片| 精品人妻在线不人妻| 男的添女的下面高潮视频| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 免费黄色在线免费观看| 亚洲少妇的诱惑av| 亚洲中文av在线| 国产成人91sexporn| 热99国产精品久久久久久7| 国产成人系列免费观看| 亚洲欧美一区二区三区黑人| 99热网站在线观看| 91老司机精品| 青春草国产在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 丝袜美腿诱惑在线| av福利片在线| 亚洲欧洲精品一区二区精品久久久 | av又黄又爽大尺度在线免费看| 一本大道久久a久久精品| 久久精品熟女亚洲av麻豆精品| 中文天堂在线官网| 黄片播放在线免费| 九九爱精品视频在线观看| 在线看a的网站| 免费久久久久久久精品成人欧美视频| 天堂俺去俺来也www色官网| 色94色欧美一区二区| 日日撸夜夜添| 国产日韩欧美视频二区| 夫妻性生交免费视频一级片| 中文字幕最新亚洲高清| 中文字幕亚洲精品专区| 丝袜美腿诱惑在线| 国产高清国产精品国产三级| 如何舔出高潮| 亚洲国产精品一区二区三区在线| 久久99热这里只频精品6学生| 午夜影院在线不卡| 中文欧美无线码| 9色porny在线观看| 精品国产超薄肉色丝袜足j| 99re6热这里在线精品视频| 亚洲国产av新网站| 男人爽女人下面视频在线观看| tube8黄色片| 久久久久久久久久久久大奶| 九色国产91popny在线| 最好的美女福利视频网| 日本a在线网址| 亚洲第一电影网av| 午夜两性在线视频| 亚洲自拍偷在线| 亚洲欧美一区二区三区黑人| 色综合欧美亚洲国产小说| 一本大道久久a久久精品| 国产av在哪里看| 亚洲精品中文字幕一二三四区| 少妇的丰满在线观看| 亚洲av成人一区二区三| 国产亚洲精品第一综合不卡| 日本免费a在线| 精品久久蜜臀av无| 久久精品国产综合久久久| 亚洲欧美日韩无卡精品| 搡老妇女老女人老熟妇| www.精华液| 妹子高潮喷水视频| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久久久电影| 国产精品亚洲美女久久久| 亚洲少妇的诱惑av| 中文字幕另类日韩欧美亚洲嫩草| 成人三级做爰电影| 一级a爱视频在线免费观看| 欧美日本中文国产一区发布| 人人妻人人澡人人看| 窝窝影院91人妻| 亚洲精品中文字幕在线视频| 久久精品aⅴ一区二区三区四区| 一级片免费观看大全| 国产精品 欧美亚洲| 久久久久久免费高清国产稀缺| 亚洲熟妇熟女久久| 日本撒尿小便嘘嘘汇集6| 午夜激情av网站| 亚洲片人在线观看| 大香蕉久久成人网| 久久久久久久久中文| 午夜福利18| 最近最新中文字幕大全免费视频| 成人国语在线视频| 美女国产高潮福利片在线看| av视频免费观看在线观看| 18美女黄网站色大片免费观看| 亚洲第一电影网av| 欧美一区二区精品小视频在线| 麻豆av在线久日| 日韩精品中文字幕看吧| 满18在线观看网站| 9热在线视频观看99| 国产高清有码在线观看视频 | 黄色片一级片一级黄色片| 麻豆久久精品国产亚洲av| 一级,二级,三级黄色视频| 午夜激情av网站| 欧美成人性av电影在线观看| 韩国av一区二区三区四区| 99在线人妻在线中文字幕| 日本在线视频免费播放| 一区二区日韩欧美中文字幕| 国产精品 欧美亚洲| 中文字幕人成人乱码亚洲影| 亚洲av电影不卡..在线观看| 精品一品国产午夜福利视频| 国产av一区在线观看免费| 国产欧美日韩一区二区三| www.熟女人妻精品国产| 国产91精品成人一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 久久婷婷人人爽人人干人人爱 | 老汉色av国产亚洲站长工具| 看免费av毛片| 国产黄a三级三级三级人| 怎么达到女性高潮| 变态另类丝袜制服| 色综合婷婷激情| 女生性感内裤真人,穿戴方法视频| 怎么达到女性高潮| av有码第一页| 色播在线永久视频| 黄色视频,在线免费观看| 亚洲国产日韩欧美精品在线观看 | 99久久精品国产亚洲精品| 国产精品久久久久久亚洲av鲁大| 十分钟在线观看高清视频www| 亚洲伊人色综图| 久久这里只有精品19| 亚洲 国产 在线| 可以在线观看的亚洲视频| 一区二区三区高清视频在线| 女警被强在线播放| 神马国产精品三级电影在线观看 | 成人手机av| 亚洲精品国产区一区二| 可以在线观看的亚洲视频| 日韩精品青青久久久久久| 国产成人av教育| 午夜精品国产一区二区电影| 非洲黑人性xxxx精品又粗又长| 大型黄色视频在线免费观看| 很黄的视频免费| 国产三级在线视频| 亚洲成国产人片在线观看| 巨乳人妻的诱惑在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品影院久久| 亚洲一码二码三码区别大吗| 夜夜躁狠狠躁天天躁| 亚洲国产精品sss在线观看| 免费看十八禁软件| 亚洲专区字幕在线| 19禁男女啪啪无遮挡网站| 人妻久久中文字幕网| 少妇 在线观看| 国产一区二区三区在线臀色熟女| 国产aⅴ精品一区二区三区波| 国产精品影院久久| 如日韩欧美国产精品一区二区三区| 亚洲av五月六月丁香网| 亚洲无线在线观看| 人人澡人人妻人| 男女下面进入的视频免费午夜 | 母亲3免费完整高清在线观看| 欧美乱码精品一区二区三区| 国产一区二区在线av高清观看| 久久午夜综合久久蜜桃| 国产真人三级小视频在线观看| 欧美日韩福利视频一区二区| 中文亚洲av片在线观看爽| 亚洲人成77777在线视频| 日韩欧美国产在线观看| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 午夜激情av网站| 亚洲七黄色美女视频| 又大又爽又粗| 亚洲最大成人中文| 视频在线观看一区二区三区| 老司机靠b影院| 国产麻豆69| 亚洲国产精品sss在线观看| 亚洲男人的天堂狠狠| 国产一卡二卡三卡精品| 日日摸夜夜添夜夜添小说| 欧美激情高清一区二区三区| 欧美精品亚洲一区二区| tocl精华| 女人高潮潮喷娇喘18禁视频| 亚洲国产中文字幕在线视频| 中文字幕人妻丝袜一区二区| 黑人操中国人逼视频| 亚洲欧美一区二区三区黑人| 日本黄色视频三级网站网址| 国产成人av教育| 成人三级黄色视频| 黄色丝袜av网址大全| 99国产精品免费福利视频| 丝袜美足系列| 久久亚洲真实| 身体一侧抽搐| 亚洲免费av在线视频| 一区福利在线观看| 黄色a级毛片大全视频| or卡值多少钱| 亚洲av五月六月丁香网| 午夜精品在线福利| 人人妻人人澡欧美一区二区 | 人人妻人人爽人人添夜夜欢视频| 亚洲中文日韩欧美视频| av网站免费在线观看视频| 美女国产高潮福利片在线看| 午夜精品在线福利| 最好的美女福利视频网| 成人18禁在线播放| 国产一区二区在线av高清观看| 香蕉久久夜色| 一级片免费观看大全| 老汉色∧v一级毛片| 欧美丝袜亚洲另类 | 国产激情欧美一区二区| 国产伦人伦偷精品视频| 动漫黄色视频在线观看| 日本vs欧美在线观看视频| 精品一品国产午夜福利视频| 91大片在线观看| 久久精品国产清高在天天线| 午夜久久久在线观看| 成人三级黄色视频| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美日韩在线播放| 午夜精品在线福利| 12—13女人毛片做爰片一| 黄片大片在线免费观看| 老司机深夜福利视频在线观看| 精品久久久久久久毛片微露脸| 一级毛片女人18水好多| 欧美激情 高清一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 亚洲 国产 在线| av福利片在线| 丝袜在线中文字幕| tocl精华| 色哟哟哟哟哟哟| 欧美成狂野欧美在线观看| 在线观看www视频免费| 禁无遮挡网站| 动漫黄色视频在线观看| a在线观看视频网站| 亚洲专区中文字幕在线| 黄色片一级片一级黄色片| ponron亚洲| 99在线视频只有这里精品首页| 成年版毛片免费区| 一夜夜www| 欧美日韩瑟瑟在线播放| 狠狠狠狠99中文字幕| 亚洲自偷自拍图片 自拍| av在线天堂中文字幕| 丝袜美足系列| 国产精品爽爽va在线观看网站 | 人妻久久中文字幕网| av视频免费观看在线观看| 制服人妻中文乱码| 老汉色av国产亚洲站长工具| 黄色毛片三级朝国网站| 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 69av精品久久久久久| av免费在线观看网站| 国产亚洲av嫩草精品影院| 久久人人97超碰香蕉20202| 性欧美人与动物交配| 国产精品香港三级国产av潘金莲| 中文字幕最新亚洲高清| 亚洲九九香蕉| 色哟哟哟哟哟哟| 成熟少妇高潮喷水视频| 女性生殖器流出的白浆| 亚洲欧美日韩另类电影网站| av片东京热男人的天堂| 免费在线观看完整版高清| 咕卡用的链子| 高清毛片免费观看视频网站| 欧美av亚洲av综合av国产av| 久久青草综合色| 一个人观看的视频www高清免费观看 | 99国产精品免费福利视频| 少妇熟女aⅴ在线视频| 精品熟女少妇八av免费久了| 日本 av在线| 亚洲av片天天在线观看| 桃红色精品国产亚洲av| 一级a爱片免费观看的视频| 窝窝影院91人妻| 久久天躁狠狠躁夜夜2o2o| 国产三级在线视频| 人妻丰满熟妇av一区二区三区| 日韩av在线大香蕉| 亚洲伊人色综图| 黄片播放在线免费| 欧美在线一区亚洲| 波多野结衣高清无吗| 天堂影院成人在线观看| 国产伦人伦偷精品视频| 午夜福利高清视频| 天堂影院成人在线观看| 性色av乱码一区二区三区2| ponron亚洲| 国产人伦9x9x在线观看| 亚洲av第一区精品v没综合| 18禁观看日本| 成人18禁高潮啪啪吃奶动态图| 亚洲七黄色美女视频| 欧美+亚洲+日韩+国产| xxx96com| 免费搜索国产男女视频| 婷婷精品国产亚洲av在线| 国产成人av激情在线播放| 免费搜索国产男女视频| 国产三级在线视频| 啦啦啦韩国在线观看视频| 视频在线观看一区二区三区| 成人三级黄色视频| 成人精品一区二区免费| 久久久久九九精品影院| 亚洲国产看品久久| netflix在线观看网站| 好男人在线观看高清免费视频 | 色综合亚洲欧美另类图片| 精品一区二区三区av网在线观看| 19禁男女啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品成人综合色| 在线观看www视频免费|