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

    濕加松無性系表型遺傳多樣性研究

    2019-03-19 08:54:14郭文麗李義良趙奮成廖仿炎鐘歲英林昌明葉威方
    植物研究 2019年2期
    關(guān)鍵詞:皮厚材積冠幅

    郭文麗 李義良 趙奮成 鐵 軍 廖仿炎鐘歲英 林昌明 葉威方

    (1.山西師范大學(xué)生命科學(xué)學(xué)院,臨汾 041004; 2.廣東省森林培育與保護利用重點實驗室,廣州 510520; 3.廣東省林業(yè)科學(xué)研究院,廣州 510520; 4.長治學(xué)院生物科學(xué)與技術(shù)系,長治 046011; 5.廣東省臺山紅嶺種子園,臺山 529200)

    濕加松(Pinuselliottii×P.caribaeaMorelet var.hondurensis)屬松科(Pinaceae)松屬(Pinus),為裸子植物。濕加松是經(jīng)控制授粉后得到的雜交子代,其母本為濕地松(PinuselliottiiEngelmann),父本為洪都拉斯加勒比松變種(PinuscaribaeaMorelet var.hondurensis)。濕加松不僅具有母本干型通直圓滿、分枝情況良好、抗風(fēng)抗病蟲害能力強、耐水漬、耐短期低溫等優(yōu)勢,還具有父本早期生長迅速、生長量大等優(yōu)勢[1]。從上世紀末開始,我國木材市場就出現(xiàn)了緊缺產(chǎn)脂力高和出材量大的松樹良種的現(xiàn)象,濕加松依靠獨有的優(yōu)勢已成為我國華南地區(qū)兼具有建材產(chǎn)脂和大徑材林培育潛力的重要樹種[2]。

    目前,我國濕加松無性系有兩個來源:其一為引種自澳大利亞昆士蘭州的濕加松,其二是上世紀九十年代以來廣東省林業(yè)科學(xué)研究院與臺山紅嶺種子園等單位聯(lián)合培育的濕加松[3]。經(jīng)過二十多年的大量研究,同等立地條件下我國廣東省培育的部分濕加松優(yōu)良家系在材積生長方面上優(yōu)于澳大利亞F1代家系。

    前人在研究不同植物表型性狀普遍采用了主成分分析和聚類分析的方法,例如豌豆[4]、茶樹[5~6]、谷子[7]、陸地棉[8]、甜高粱[9~10]、杜鵑紅山茶[11]等作物已經(jīng)進行了相關(guān)研究。本研究擬選取廣東省臺山市紅嶺種子園35個濕加松無性系作為試驗材料,測量其8個表型性狀,采用主成分分析和聚類分析相結(jié)合,綜合分析濕加松無性系表型遺傳多樣性,為更高效的選擇親本提供理論依據(jù),進一步為良種選育奠定基礎(chǔ)。

    1 試驗材料及方法

    1.1 試驗材料、方法

    濕加松無性系均來自廣東省臺山市紅嶺種子園,試驗于2017年8月25日~2017年9月14日在紅嶺種子園進行。該種子園營建于2000~2003年,總面積18.5 hm2,分為3個大區(qū)39個小區(qū),建園無性系共51個,每個小區(qū)16個無性系,無性系配置采用了三位順序錯位排列法,株行距為5 m×5 m。建園無性系均選自早期營建的濕加松遺傳測定林,從測定林中選擇生長量大、干形通直、分枝良好、結(jié)實正常、樹體健康的優(yōu)良單株,采集其穗條帶回苗圃中嫁接、繁殖,由此組成建園材料。該地屬南亞熱帶海洋性氣候,年均溫21.8℃,年降水量1 940 mm;平均海拔30 m,土壤為花崗巖發(fā)育的酸性磚紅壤,pH值5.0~5.5,土層深厚。自然條件下適合濕加松生長發(fā)育、開花結(jié)實。本試驗選定35個無性系,包括10個來自澳洲全同胞家系(EHA),2個來自澳洲半同胞家系F2(EHAF2),23個來自臺山自配全同胞家系,每個無性系從種子園隨機選取3~6分株不等,共160分株。試驗材料詳細信息見表1。

    1.2 表型性狀調(diào)查

    設(shè)計對35個濕加松無性系從生長性狀和形質(zhì)性狀兩個方面調(diào)查,生長性狀包括樹高、胸徑、材積;形質(zhì)性狀包括冠幅、分枝角、側(cè)枝直徑、輪枝距、皮厚。

    樹高:用測高儀(VertexⅣ超聲波樹木測高儀)測量選定濕加松的樹高(單位:m)。

    胸徑:用胸徑尺測量選定濕加松的胸徑(1.3 m處)(單位:cm)。

    材積:利用樹高、胸徑計算選定濕加松的材積(單位:m3,單株材積:V)的計算公式[12]如下:

    V=(f×π×H×D2)/40 000

    (1)

    式中,f為形數(shù),取f=0.5,π為圓周率,H為樹高,D為胸徑。

    冠幅:先測量選定材料的東西、南北向冠幅,再計算平均值(單位:m)。

    表135個濕加松無性系材料名稱表

    Table135clonesofP.elliottii×P.caribaeaMoreletvar.hondurensismaterialnametable

    選3個有代表性的側(cè)枝,距地面3~5 m,分別測量分枝角、側(cè)枝直徑、輪枝距。

    分枝角:測量選定的3個側(cè)枝的分枝角(單位:°)。

    側(cè)枝直徑:用胸徑尺測量選定的3個側(cè)枝直徑,再計算平均值(單位:cm)。

    輪枝距:分別測量選定的3個側(cè)枝之間距離,再計算平均值(單位:cm)。

    皮厚:位于樹高1.3 m處測量選定材料的樹皮厚度(單位:cm)。

    1.3 數(shù)據(jù)處理及統(tǒng)計分析

    本研究采用R軟件[13]對35個濕加松無性系的表型數(shù)據(jù)進行分析,使用summary()函數(shù)計算各表型性狀的最小值、最大值、平均值,使用stat.desc()函數(shù)計算了標準誤、方差、變異系數(shù),以及進行了多樣性指數(shù)、表型相關(guān)性分析、主成分和聚類分析。

    遺傳多樣性指數(shù)采用Shannon-Weaver多樣性指數(shù)(H),其計算公式為:

    H=-∑PilnPi

    (2)

    式中:H為任一組中某一性狀的多樣性程度;P為某一性狀第i個級別的種質(zhì)份數(shù)占總種質(zhì)材料的百分比;ln表示自然對數(shù)。

    其中表型性狀的相關(guān)性分析是間接選擇的重要參數(shù),記為rtt。計算公式如下:

    (3)

    式中:Covt1,t2類是兩性狀的協(xié)方差;Vt1、Vt2是兩性狀的遺傳方差。

    主成分分析采用降維把多維空間多個相關(guān)變量簡化為少量幾個相互獨立的綜合指標,可以更好地反映出各指標間的綜合情況[14]。本研究使用psych程序包對35個濕加松無性系進行了主成分分析。

    聚類分析采用了系統(tǒng)聚類法(層次聚類法),利用函數(shù)hclust()對35個濕加松無性系進行聚類分析。

    2 結(jié)果與分析

    2.1 35個濕加松無性系的表型性狀差異

    對35個濕加松無性系的8個表型性狀進行了描述性統(tǒng)計分析,分析結(jié)果見表2。

    表2可知,濕加松樹高變幅8.60~17.00 m,平均值為13.37 m;胸徑變幅12.70~36.50 cm,平均值為24.85 cm;材積變幅0.07~0.71 m3,平均值為0.34 m3;冠幅變幅3.00~10.40 m,平均值為5.45 m;分枝角變幅43.33~90.00°,平均值為75.92°;側(cè)枝直徑變幅2.33~7.74 cm,平均值為4.44 cm;輪枝距變幅25.00~90.00 cm,平均值為51.31 cm;皮厚變幅0.20~1.30 cm,平均值為0.67 cm。

    8個表型性狀變異范圍為12.31%~39.28%,均值為22.82%,其中生長性狀介于12.31~29.28%,均值為22.90%,形質(zhì)性狀介于13.34~47.25%,均值為26.86%。其變異系數(shù)由大到小依次為:材積>皮厚>側(cè)枝直徑>輪枝距>冠幅>胸徑> 分枝角>樹高。依據(jù)張巖等[15]可將變異系數(shù)分為三級,小(<25%)、中(25~50%)、大(>50%)。本研究只有小、中兩級。其中,材積的變異系數(shù)最大,呈現(xiàn)了豐富的變異性;樹高和分枝角的變異系數(shù)均小于15%,說明其遺傳穩(wěn)定性較好。

    表235個濕加松無性系8個表型性狀的描述性統(tǒng)計及遺傳多樣性

    Table2Descriptivestatisticsandgeneticdiversityofeightphenotypictraitsin35clonesofP.elliottii×P.caribaeaMoreletvar.hondurensis

    性狀Characteristic性狀Characteristic最小值Min最大值Max平均值±標準誤Mean±SE方差Variance變異系數(shù)CV(%)多樣性指數(shù)H生長性狀Growth traits樹高Tree height(m)8.6017.0013.37±0.132.7112.311.61胸徑Breast diameter(cm)12.7036.5024.85±0.3418.0517.101.67材積Volume(m3)0.070.710.34±0.010.0239.281.80平均值A(chǔ)verage value22.901.69形質(zhì)性狀Shape character冠幅Crown width(m)3.0010.405.45±0.101.4121.831.83分枝角Branching angle(°)43.3390.0075.92±0.81102.5913.341.50側(cè)枝直徑Lateral branch diameter(cm)2.337.474.44±0.091.1924.571.80輪枝距Wheel pitch(cm)25.0090.0051.31±1.01156.7424.401.78皮厚Skin thickness(cm)0.201.300.67±0.020.0429.761.75平均值A(chǔ)verage value22.781.73總均值Total mean22.821.72

    2.2 35個濕加松無性系表型性狀的遺傳多樣性

    對供試無性系的8個性狀進行了遺傳多樣性分析(表2)。從表2可知,8個性狀多樣性指數(shù)變化范圍在1.50~1.83,均值為1.72,生長性狀介于1.61~1.80,均值為1.69,生長性狀介于1.50~1.94,均值為1.77。其遺傳多樣性指數(shù)由高到低依次為:冠幅>材積>側(cè)枝直徑>輪枝距>皮厚>胸徑>樹高>分枝角,表明了表型性狀的遺傳多樣性均較豐富。

    2.3 35個濕加松無性系表型性狀的相關(guān)性分析

    對8對不同表型性狀進行相關(guān)性分析(表3)。從表3可知,各性狀之間存在著不同程度的相關(guān)性,且大部分呈顯著或極顯著相關(guān)。其中呈極顯著正相關(guān)的有11對,依次為材積與皮厚(0.98)、冠幅與側(cè)枝直徑(0.94)、樹高與皮厚(0.92)、樹高與材積(0.89)、材積與側(cè)枝直徑(0.82)、樹高與胸徑(0.81)、材積與冠幅(0.72)、胸徑與側(cè)枝直徑(0.71)、胸徑與冠幅(0.69)、樹高與輪枝距(0.62)、樹高與側(cè)枝直徑(0.46);呈顯著正相關(guān)的有5對,依次為輪枝距與皮厚(0.72)、胸徑與輪枝距(0.63)、側(cè)枝直徑與皮厚(0.55)、材積與輪枝距(0.51)、樹高與冠幅(0.36)。表明供試材料材積較高且皮較厚、冠幅與側(cè)枝直徑較大、樹體高大以及輪枝距較大;皮厚與材積、樹高呈極顯著相關(guān),與輪枝距、側(cè)枝直徑呈顯著正相關(guān),因此對供試材料的育種選擇中應(yīng)以材積為主要指標,其次為樹高。

    表3 8個表型性狀之間的相關(guān)性分析

    注:***,**和*為分別表示在0.001、0.01和0.05水平上存在顯著相關(guān)。

    Note:***,** and * indicate significant correlations at the levels of 0.001,0.01 and 0.05,respectively.

    表4濕加松無性系各類群表型性狀的特征

    Table4TheaverageandcoefficientvariationamountofphenotypictraitsindifferentgermplasmgroupsofP.elliottii×P.caribaeaMoreletvar.hondurensis

    性狀Traits資源數(shù)Number樹高Tree height(m)胸徑Breast diameter(cm)材積Volume(m3)冠幅Crown width(m)分枝角Branching angle(°)側(cè)枝直徑Lateral branch diameter(cm)輪枝距Wheel pitch(cm)皮厚Skin thickness(cm)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)均值Mean變異系數(shù)CV(%)種質(zhì)類群Germplasms groupsⅠ213.080.0424.090.050.320.114.730.0957.710.044.150.3346.690.050.690.41Ⅱ314.521.4526.811.280.423.115.820.6669.24138.894.540.0165.45500.790.005Ⅲ113.830.0226.471.280.380.0015.472.658034.723.510.0023412.50.670.005Ⅳ1313.1224.0610.580.320.0055.540.1573.512.54.424.9146.0912.50.640.02Ⅴ1213.720.7226.23320.380.025.660.00179.27868.064.852.4255.866.130.70.005Ⅵ412.811.2823.4280.650.270.054.851.8179.71312.53.80.6855.4100.620.005

    圖1 35個濕加松無性系的系統(tǒng)聚類圖Fig.1 Hierarchical Cluster Map of 35 Clones of P.elliottii×P.caribaea Morelet var. hondurensis

    2.4 聚類分析

    采用R軟件對35個濕加松無性系進行了系統(tǒng)聚類分析(圖1)。在歐氏距離為12時,可將35個無性系分為6大類群,每個類群的份數(shù)、表型性狀的平均值以及變異系數(shù)(表4)。

    第Ⅰ類:有2個無性系,均來源于澳洲全同胞家系,共8分株,占供試材料總數(shù)的5%。主要表現(xiàn)為冠幅較小,可為培育窄冠型良種提供種質(zhì)。冠幅和分枝角(<60°)的均值在六大類群里均最小,變異系數(shù)均較小。

    第Ⅱ類:有3個無性系,1個來源于臺山全同胞家系、1個來源于澳洲全同胞家系、1個來源于澳洲半同胞家系F2,共13分株,占總材料的8.125%。主要表現(xiàn)為樹體高大,樹冠寬,粗皮且材積大。樹高、胸徑、材積、冠幅、輪枝距和皮厚的均值在六大類群里均最大,且材積的變異系數(shù)最大。綜合分析可知,該類群長勢最優(yōu),可為培育速生良種提供種質(zhì)材料,有較高的用材經(jīng)濟價值。

    第Ⅲ類:只有1個無性系,來源于臺山全同胞家系,該無性系長勢良好,但分枝角的均值最大(80°),側(cè)枝直徑與輪枝距的均值最小,變異系數(shù)都較低,僅3分株,占總材料的1.875%。

    第Ⅳ類:共13個無性系,分別有3個來源于澳洲全同胞家系、10個均來源于臺山全同胞家系,共64分株,占總材料的40%。其整體長勢較差,且樹高和側(cè)枝直徑的變異系數(shù)均最大,需進一步培育從中選取優(yōu)良品種。

    第Ⅴ類:包括了12個無性系,有1個來源于澳洲半同胞家系F2、2個來源于澳洲全同胞家系、9個來源于臺山全同胞家系,共55分株,占總材料的34.375%。其側(cè)枝直徑的均值最大,分枝角及輪枝距的變異系數(shù)最大,說明選擇范圍較寬,整體上的長勢優(yōu)于第Ⅳ類,適合后期大量種植。

    第Ⅵ類:包括4個無性系,2個來源于澳洲全同胞家系、2個來源于臺山全同胞家系,共17分株占總材料的10.625%。其樹高、材積與皮厚的均值最小,胸徑的變異系數(shù)最大,說明該類群長勢最差,但皮薄,胸徑選擇范圍較寬,可為培育薄皮型良種提供種質(zhì)。

    2.5 主成分分析

    本研究針對8個性狀采用R軟件的程序包psych進行了主成分分析(表5)。結(jié)果表明,前3個主成分的特征值均大于1,固選為主成分,3個主成分的累積貢獻率達86.00%,基本上可表現(xiàn)6個性狀的綜合變異情況。第1主成分貢獻率為39.36%,樹高、胸徑和皮厚的系數(shù)明顯高于其它性狀,說明了濕加松表型遺傳多樣性的重要變異因子為樹高、胸徑和皮厚三個表型性狀;第2主成分貢獻率為30.32%,冠幅、側(cè)枝直徑的特征向量值均較高,說明第2主成分是冠幅、側(cè)枝直徑的綜合反映;第3主成分貢獻率為16.33%,系數(shù)最大的是分枝角(0.92),說明第3主成分主要反映了分枝角因子。

    表535個濕加松無性系基于8個表型性狀的5個主成分的特征向量和貢獻率

    Table5Eigenvectorandcontributionrateoffiveprincipalcomponentsbasedoneightphenotypictraitsin35clonesofP.elliottii×P.caribaeaMoreletvar.hondurensis

    表型性狀Phenotypic traits主成分Principal componentRC1RC2RC3樹高Tree height0.790.19-0.11胸徑Breast diameter0.70.50.05材積Volume0.690.58-0.11冠幅Crown width0.130.850.12分枝角Branching angle-0.12-0.020.92側(cè)枝直徑Lateral branch diameter0.160.85-0.15輪枝距Wheel pitch0.63-0.010.42皮厚Skin thickness0.80.11-0.19特征值Characteristic value2.72.081.12方差貢獻率Variance contribution rate(%)39.3630.3216.33累積方差貢獻率Cumulative variance contribution rate(%)39.3669.6786.00

    注:下劃線表示每個主成分較高的特征向量。

    Note:Underline represents the higher eigenvector for each principal component.

    3 討論

    3.1 濕加松無性系的遺傳多樣性分析

    遺傳多樣性作為種質(zhì)資源的重要評價因子,也是培育優(yōu)良品種的研究提供重要依據(jù)[16]。調(diào)查表型性狀的遺傳多樣性是衡量表型變異最直觀最簡易的途徑,但易受環(huán)境、個體差異或標記性狀較少等諸多外在因素影響。前人在表型性狀對松屬植物遺傳多樣性上的研究[17]較少,僅在賀蘭山杜松(Juniperusrigida)[18]、歐洲赤松(PinussylvestrisL.)[19]等少數(shù)松樹上有過相關(guān)報道,大多采用RFLP標記[20]、RAPD標記[21]、SSR標記[22~24]、ISSR[25]等分子標記法研究松樹遺傳變異規(guī)律。本文通過對35個濕加松無性系的8個表型性狀進行遺傳多樣性分析,結(jié)果發(fā)現(xiàn)8個表型性狀的平均遺傳多樣性指數(shù)較高,表明其有著較高的豐富度與均勻度,豐富了我國濕加松無性系的遺傳背景。另外,從計算結(jié)果可以看出大部分性狀的變異系數(shù)和遺傳多樣性基本統(tǒng)一,材積均較高,分枝角和樹高均較低,說明其離散程度相似。影響濕加松變異的有遺傳與非遺傳兩大因素,各占一半,但研究表明家系的不同對遺傳影響最強,其次是無性系的不同[26],所以該研究的遺傳多樣性分析結(jié)果可能受到了不同來源的無性系影響。

    3.2 濕加松無性系的相關(guān)性分析

    通過相關(guān)性分析,9對表型性狀之間存在著較高的相關(guān)性,其中有11對極顯著正相關(guān),5對顯著正相關(guān)。研究表明樹體越高,皮較厚,胸徑、材積、側(cè)枝直徑、輪枝距也較大,這與梁德洋等[27]對紅松無性系的研究結(jié)果相似,樹高、胸徑與材積之間均呈極顯著正相關(guān)。李義良等[12]研究表明濕地松、加勒比松的樹高、胸徑、材積、皮厚及冠幅兩兩間均呈極顯著正相關(guān),本研究中除胸徑與材積、胸徑與皮厚及冠幅與皮厚均呈不顯著相關(guān)其它性狀也呈顯著正相關(guān)。綜上所述,在選育濕加松時應(yīng)優(yōu)先參考材積為第一指標,同時應(yīng)參考樹高、皮厚、胸徑、冠幅及側(cè)枝直徑等綜合指標,為今后的選育工作給予理論指導(dǎo)。本研究中參試無性系數(shù)量不同可能對試驗結(jié)果產(chǎn)生一定的誤差,后期應(yīng)增加材料的數(shù)量,另外應(yīng)增加產(chǎn)脂力、木材密度、種子量等性狀的收集,進一步分析更具說服力[28]。此外,用R軟件進行相關(guān)性分析要比IBM SPSS Statisticsr類似軟件更為直觀,其中極顯著相關(guān)分成0~0.001和0.001~0.01兩個級別,更有利于最優(yōu)指標的提取。

    3.3 濕加松無性系的聚類分析

    利用濕加松無性系的8個表型性狀相關(guān)數(shù)據(jù)進行系統(tǒng)聚類分析,將表型差異較小的分為一個類群,把35個無性系分為6大類群。第Ⅰ類為選育窄冠型優(yōu)良基因的材料;第Ⅱ類為后期選育速生良種的材料,木材經(jīng)濟價值較高;第Ⅲ類雖長勢良好,但分枝角較大,應(yīng)縮短生長周期;第Ⅳ類長勢均較差,需進一步具體分析再選育良種;第Ⅴ類長勢優(yōu),為后期選優(yōu)的主要良種材料;第Ⅵ類為選育薄皮型良種的材料。通過對各類群的具體分析,部分來源相同的聚在一類,如來源于澳洲全同胞家系的EHA12與EHA08-2、EHA11與EHA14分別聚為一類,來源于臺山全同胞家系的父本為H01與H03、H05與H06的大多聚在一類,說明它們在表型上的變異較小。大部分來源相同的種質(zhì)聚在不同類群,說明無性系在表型上存在了較大的差異,與無性系的來源沒有必然的聯(lián)系,這與萬述偉等[4]和王麗俠等[29]的研究結(jié)果相同。在育種時可依據(jù)表型性狀的聚類結(jié)果,同時結(jié)合分子標記法,去合理分配種質(zhì),可有效避免近交,減少控制授粉的工作,大大縮短育種周期,提高濕加松無性系的改良效率。

    3.4 濕加松無性系的主成分分析

    本研究對濕加松的8個表型性狀進行主成分分析,提取了3個主成分,反映了6個性狀的大部分信息,累積貢獻率高達86.00%,與相關(guān)性分析結(jié)果相似。第1主成分主要代表了樹高、胸徑和皮厚三個生長性狀的綜合情況,體現(xiàn)出濕加松在樹高、胸徑、皮厚三方面明顯的雜種優(yōu)勢,這與李義良等[30]對濕加松雜種優(yōu)勢的研究結(jié)果相似。此外,第2主成分主要反映了冠幅、側(cè)枝直徑的綜合情況;第3主成分主要反映了分枝角的情況,所以這些性狀也可評價濕加松的優(yōu)劣。

    綜上所述,對不同來源濕加松無性系的表型性狀進行分析,初步了解了濕加松種質(zhì)的表型遺傳多樣性情況,為今后種質(zhì)利用及良種選育奠定了理論基礎(chǔ)。

    猜你喜歡
    皮厚材積冠幅
    登記簿里的信任
    黨員生活(2023年7期)2023-07-26 21:05:23
    無人機遙感影像提取的單木冠幅數(shù)據(jù)在桉樹林分蓄積量估測中的應(yīng)用1)
    城市綠地微環(huán)境對土壤動物群落多樣性的影響
    施肥對三江平原丘陵區(qū)長白落葉松人工林中齡林單木樹冠圓滿度影響
    基于無人機高分影像的冠幅提取與樹高反演
    5齡熱墾628材積量季節(jié)生長節(jié)律與氣象因子關(guān)聯(lián)度初步研究
    福建省闊葉樹二元材積方程修訂
    “臉皮厚”
    皮厚,不怕冷
    誤撞等3則
    一级毛片 在线播放| 亚洲一级一片aⅴ在线观看| xxx大片免费视频| 校园人妻丝袜中文字幕| 一级黄片播放器| 在线观看一区二区三区激情| 欧美最新免费一区二区三区| 哪个播放器可以免费观看大片| 国产精品欧美亚洲77777| 国产在视频线精品| 一本大道久久a久久精品| 中文字幕精品免费在线观看视频 | 免费观看a级毛片全部| 在线观看免费高清a一片| 久久97久久精品| 精品人妻熟女av久视频| 丝瓜视频免费看黄片| 另类亚洲欧美激情| 国产免费又黄又爽又色| 黄色欧美视频在线观看| 久久久久久久亚洲中文字幕| av在线老鸭窝| 亚洲第一区二区三区不卡| 极品人妻少妇av视频| 日韩精品免费视频一区二区三区 | 亚洲精品,欧美精品| 日日撸夜夜添| 免费观看的影片在线观看| 亚洲人成网站在线播| 交换朋友夫妻互换小说| 亚洲内射少妇av| 精品一区二区三区视频在线| 考比视频在线观看| 免费人妻精品一区二区三区视频| 黑人猛操日本美女一级片| a级毛片黄视频| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 高清欧美精品videossex| 日日啪夜夜爽| 午夜福利在线观看免费完整高清在| 亚洲综合色惰| 午夜91福利影院| 久久精品国产a三级三级三级| 大香蕉久久网| 亚洲欧美成人综合另类久久久| 国产高清有码在线观看视频| 欧美成人精品欧美一级黄| 欧美+日韩+精品| 亚洲精品,欧美精品| 国产日韩一区二区三区精品不卡 | 狠狠婷婷综合久久久久久88av| 搡女人真爽免费视频火全软件| 亚洲精品av麻豆狂野| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看 | 啦啦啦中文免费视频观看日本| 欧美一级a爱片免费观看看| 国产精品久久久久久久电影| 一边亲一边摸免费视频| 99热6这里只有精品| 国产av一区二区精品久久| 亚洲精品视频女| 亚洲精品国产色婷婷电影| 狂野欧美激情性bbbbbb| 国产精品国产av在线观看| 久久免费观看电影| 人人澡人人妻人| 久久久精品区二区三区| 国产成人精品福利久久| 精品久久久久久久久亚洲| xxxhd国产人妻xxx| 国产亚洲精品第一综合不卡 | 日韩三级伦理在线观看| 日韩亚洲欧美综合| 日韩欧美一区视频在线观看| 久久99蜜桃精品久久| 黄色怎么调成土黄色| 丝袜在线中文字幕| 国产高清不卡午夜福利| 午夜久久久在线观看| 人妻少妇偷人精品九色| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 精品少妇内射三级| 国产免费又黄又爽又色| 亚洲美女视频黄频| 日韩亚洲欧美综合| 国产精品国产av在线观看| 亚洲精品,欧美精品| 午夜免费鲁丝| 国产亚洲精品第一综合不卡 | 日韩av免费高清视频| 美女国产高潮福利片在线看| 伊人久久精品亚洲午夜| 亚洲欧美成人综合另类久久久| videossex国产| 精品亚洲成国产av| 国产精品 国内视频| 在线观看国产h片| 国产69精品久久久久777片| 日韩大片免费观看网站| 在线免费观看不下载黄p国产| 18+在线观看网站| 日产精品乱码卡一卡2卡三| 黑人欧美特级aaaaaa片| 中文字幕制服av| 成人亚洲精品一区在线观看| 久久人人爽人人爽人人片va| 国产乱来视频区| 不卡视频在线观看欧美| 七月丁香在线播放| 97精品久久久久久久久久精品| 亚洲欧美成人精品一区二区| 国产一区二区在线观看日韩| 在线 av 中文字幕| 午夜福利,免费看| 日本色播在线视频| 王馨瑶露胸无遮挡在线观看| 黄色欧美视频在线观看| 欧美人与性动交α欧美精品济南到 | 久久久久国产网址| 亚洲国产av影院在线观看| 精品久久久精品久久久| 精品久久久噜噜| 美女大奶头黄色视频| 国产片内射在线| 999精品在线视频| 五月开心婷婷网| 桃花免费在线播放| 97精品久久久久久久久久精品| 哪个播放器可以免费观看大片| 91精品国产国语对白视频| 亚洲久久久国产精品| 午夜日本视频在线| 嫩草影院入口| 国产熟女欧美一区二区| 男女国产视频网站| 久久精品国产a三级三级三级| 永久网站在线| 精品少妇久久久久久888优播| 成年av动漫网址| 全区人妻精品视频| 国产片内射在线| 18禁动态无遮挡网站| 91精品一卡2卡3卡4卡| 国产精品无大码| 久久韩国三级中文字幕| tube8黄色片| av网站免费在线观看视频| 青青草视频在线视频观看| 18禁动态无遮挡网站| 一级毛片黄色毛片免费观看视频| 成人影院久久| 狠狠精品人妻久久久久久综合| 99精国产麻豆久久婷婷| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品视频女| 中文字幕av电影在线播放| av在线app专区| 亚洲精品aⅴ在线观看| 美女xxoo啪啪120秒动态图| 三上悠亚av全集在线观看| 99视频精品全部免费 在线| 七月丁香在线播放| 麻豆精品久久久久久蜜桃| 黄色配什么色好看| 高清av免费在线| 韩国高清视频一区二区三区| 男女国产视频网站| 亚洲精品一区蜜桃| 久久久亚洲精品成人影院| 午夜免费鲁丝| 亚洲丝袜综合中文字幕| 成人黄色视频免费在线看| 午夜91福利影院| 成人午夜精彩视频在线观看| 男人操女人黄网站| 精品一区二区三区视频在线| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 一级毛片 在线播放| 在线免费观看不下载黄p国产| 大香蕉久久网| 国产亚洲最大av| 在线亚洲精品国产二区图片欧美 | 建设人人有责人人尽责人人享有的| 国产在线免费精品| 97在线视频观看| 晚上一个人看的免费电影| 在线天堂最新版资源| a级毛片免费高清观看在线播放| 婷婷成人精品国产| 制服丝袜香蕉在线| 少妇猛男粗大的猛烈进出视频| 亚洲av成人精品一区久久| 国产精品成人在线| 一本大道久久a久久精品| 尾随美女入室| 啦啦啦在线观看免费高清www| 韩国高清视频一区二区三区| 欧美三级亚洲精品| 午夜福利视频精品| 高清视频免费观看一区二区| 欧美精品国产亚洲| 久久久久久久大尺度免费视频| 三级国产精品欧美在线观看| 久久韩国三级中文字幕| 久久久欧美国产精品| 极品人妻少妇av视频| 99热6这里只有精品| 黑丝袜美女国产一区| av线在线观看网站| 91午夜精品亚洲一区二区三区| av在线播放精品| 天天影视国产精品| 亚洲av成人精品一二三区| 熟女人妻精品中文字幕| 国产av精品麻豆| 日韩强制内射视频| 国产色爽女视频免费观看| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品电影小说| 看免费成人av毛片| 欧美性感艳星| 日韩精品免费视频一区二区三区 | 九色亚洲精品在线播放| 99视频精品全部免费 在线| 国产精品久久久久久久久免| 久久热精品热| 日本欧美国产在线视频| 国产在线免费精品| 99九九线精品视频在线观看视频| 日韩电影二区| 日韩大片免费观看网站| 久久亚洲国产成人精品v| 国产精品秋霞免费鲁丝片| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区黑人 | 中文字幕最新亚洲高清| 国产精品麻豆人妻色哟哟久久| 国产黄频视频在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲成人av在线免费| 久久午夜综合久久蜜桃| 免费人成在线观看视频色| 26uuu在线亚洲综合色| 国产 精品1| 成人亚洲精品一区在线观看| 久久精品久久精品一区二区三区| 久久av网站| 欧美+日韩+精品| 大码成人一级视频| 乱人伦中国视频| 国产精品国产三级专区第一集| 午夜老司机福利剧场| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 人人妻人人澡人人爽人人夜夜| 一级毛片电影观看| 久久国产精品男人的天堂亚洲 | 草草在线视频免费看| 亚洲精品成人av观看孕妇| 夜夜看夜夜爽夜夜摸| 久久久久久久精品精品| 午夜免费观看性视频| 桃花免费在线播放| 亚洲精品一区蜜桃| 男人操女人黄网站| 妹子高潮喷水视频| 国产极品天堂在线| 欧美日韩av久久| 91久久精品电影网| videossex国产| 亚洲久久久国产精品| 少妇的逼好多水| av线在线观看网站| 少妇人妻 视频| 久久精品国产亚洲网站| 欧美精品人与动牲交sv欧美| 国模一区二区三区四区视频| www.av在线官网国产| a级毛片免费高清观看在线播放| 国产高清三级在线| 少妇人妻久久综合中文| 狠狠精品人妻久久久久久综合| av在线观看视频网站免费| 精品一区二区三卡| 精品久久蜜臀av无| 丝袜喷水一区| 哪个播放器可以免费观看大片| 亚洲人与动物交配视频| 纵有疾风起免费观看全集完整版| 一本大道久久a久久精品| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 午夜免费男女啪啪视频观看| 国产精品免费大片| 天天影视国产精品| 亚洲国产欧美在线一区| 国产高清三级在线| 亚州av有码| 久久久久久人妻| 国产黄频视频在线观看| 久久久精品区二区三区| 国产精品偷伦视频观看了| 亚洲精品国产av成人精品| 久久精品人人爽人人爽视色| 美女福利国产在线| 黄色一级大片看看| 在线观看美女被高潮喷水网站| 赤兔流量卡办理| 国产欧美另类精品又又久久亚洲欧美| 极品少妇高潮喷水抽搐| 九九久久精品国产亚洲av麻豆| 亚洲国产欧美日韩在线播放| 欧美人与性动交α欧美精品济南到 | 看非洲黑人一级黄片| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 日韩亚洲欧美综合| 国产女主播在线喷水免费视频网站| 精品久久久精品久久久| 国产成人精品福利久久| 久久青草综合色| 中国美白少妇内射xxxbb| 丰满乱子伦码专区| kizo精华| 久久综合国产亚洲精品| 国产精品99久久久久久久久| 天美传媒精品一区二区| 日韩,欧美,国产一区二区三区| 天堂俺去俺来也www色官网| 日本黄大片高清| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| av女优亚洲男人天堂| av.在线天堂| 人妻人人澡人人爽人人| 国产精品久久久久久久电影| av不卡在线播放| 99精国产麻豆久久婷婷| 婷婷色综合www| 精品国产一区二区三区久久久樱花| 91精品国产国语对白视频| 久久久久久伊人网av| 99热这里只有精品一区| 成人国产麻豆网| 丝袜在线中文字幕| 最新中文字幕久久久久| 亚洲av在线观看美女高潮| 99热国产这里只有精品6| 欧美丝袜亚洲另类| 亚洲熟女精品中文字幕| 久热久热在线精品观看| 美女主播在线视频| 亚洲av综合色区一区| 国产精品偷伦视频观看了| 国产视频内射| 国产熟女午夜一区二区三区 | 亚洲精品乱久久久久久| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 两个人免费观看高清视频| 又粗又硬又长又爽又黄的视频| 免费黄频网站在线观看国产| 国语对白做爰xxxⅹ性视频网站| 久久精品人人爽人人爽视色| 能在线免费看毛片的网站| 久久久久久久久久人人人人人人| 能在线免费看毛片的网站| 成人亚洲欧美一区二区av| 毛片一级片免费看久久久久| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 久久人妻熟女aⅴ| 日本欧美国产在线视频| 母亲3免费完整高清在线观看 | 国产亚洲精品第一综合不卡 | 国产精品久久久久成人av| 精品国产一区二区三区久久久樱花| 午夜老司机福利剧场| 国产熟女欧美一区二区| 国内精品宾馆在线| a级毛片在线看网站| 水蜜桃什么品种好| 大香蕉久久成人网| 特大巨黑吊av在线直播| 日韩不卡一区二区三区视频在线| 岛国毛片在线播放| 免费观看性生交大片5| 欧美日本中文国产一区发布| 丝袜美足系列| 2018国产大陆天天弄谢| 伊人亚洲综合成人网| 欧美 亚洲 国产 日韩一| 国产男女超爽视频在线观看| 国内精品宾馆在线| 亚洲av欧美aⅴ国产| 最近的中文字幕免费完整| 国产极品天堂在线| 亚洲精品中文字幕在线视频| 日本色播在线视频| 免费黄网站久久成人精品| 视频区图区小说| 亚洲不卡免费看| 亚洲伊人久久精品综合| 又黄又爽又刺激的免费视频.| 欧美精品国产亚洲| 考比视频在线观看| 天美传媒精品一区二区| 少妇高潮的动态图| 人妻人人澡人人爽人人| 成人毛片a级毛片在线播放| 国产精品久久久久久av不卡| 日韩在线高清观看一区二区三区| 天堂8中文在线网| 欧美人与善性xxx| 色哟哟·www| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看| 日韩精品免费视频一区二区三区 | 国产成人av激情在线播放 | 亚洲欧洲日产国产| 亚洲,欧美,日韩| 最近中文字幕高清免费大全6| 久久久久国产网址| 亚洲综合精品二区| 交换朋友夫妻互换小说| 丰满饥渴人妻一区二区三| 亚洲天堂av无毛| 中文字幕免费在线视频6| 一区二区av电影网| 精品久久久噜噜| a级毛片免费高清观看在线播放| 中国美白少妇内射xxxbb| 最近中文字幕高清免费大全6| videosex国产| 五月天丁香电影| 婷婷色麻豆天堂久久| 妹子高潮喷水视频| 日韩不卡一区二区三区视频在线| 一级爰片在线观看| av天堂久久9| 免费播放大片免费观看视频在线观看| 国产日韩欧美亚洲二区| 久久久国产一区二区| av视频免费观看在线观看| 青春草亚洲视频在线观看| 人妻一区二区av| 久久韩国三级中文字幕| 在线看a的网站| 2018国产大陆天天弄谢| 亚洲人成网站在线观看播放| av在线播放精品| 最近中文字幕2019免费版| 亚洲av免费高清在线观看| 中文天堂在线官网| 精品一品国产午夜福利视频| 99九九线精品视频在线观看视频| 国产在线免费精品| 热99国产精品久久久久久7| 亚洲av中文av极速乱| 美女视频免费永久观看网站| 极品人妻少妇av视频| 中文字幕久久专区| 黄片播放在线免费| 99国产综合亚洲精品| 人妻少妇偷人精品九色| 欧美人与性动交α欧美精品济南到 | 久热久热在线精品观看| 中国三级夫妇交换| 午夜免费男女啪啪视频观看| 卡戴珊不雅视频在线播放| 成年美女黄网站色视频大全免费 | 国产国语露脸激情在线看| tube8黄色片| 国产亚洲最大av| 精品酒店卫生间| av不卡在线播放| 亚洲精品aⅴ在线观看| 性色av一级| 日韩制服骚丝袜av| 两个人免费观看高清视频| 久久99热这里只频精品6学生| 免费高清在线观看日韩| 亚洲国产精品成人久久小说| 国产午夜精品久久久久久一区二区三区| 亚洲不卡免费看| 国产亚洲欧美精品永久| 国产探花极品一区二区| 国语对白做爰xxxⅹ性视频网站| 国产精品 国内视频| 特大巨黑吊av在线直播| xxx大片免费视频| 久久av网站| av国产久精品久网站免费入址| av网站免费在线观看视频| 美女视频免费永久观看网站| 久久久久视频综合| 97超碰精品成人国产| av卡一久久| 满18在线观看网站| 久热久热在线精品观看| 777米奇影视久久| 亚洲性久久影院| 日韩欧美精品免费久久| 汤姆久久久久久久影院中文字幕| 乱人伦中国视频| 夫妻午夜视频| 看十八女毛片水多多多| 黄色视频在线播放观看不卡| 少妇丰满av| 成人综合一区亚洲| 天美传媒精品一区二区| 天堂中文最新版在线下载| 丰满饥渴人妻一区二区三| 日本vs欧美在线观看视频| 另类亚洲欧美激情| 丝袜喷水一区| 99九九线精品视频在线观看视频| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 极品人妻少妇av视频| 赤兔流量卡办理| 精品少妇内射三级| 亚洲精品亚洲一区二区| 亚洲国产精品国产精品| 精品一区在线观看国产| 国产免费现黄频在线看| 少妇人妻精品综合一区二区| 精品一区二区免费观看| 97超碰精品成人国产| 欧美日韩国产mv在线观看视频| 男女啪啪激烈高潮av片| 免费黄网站久久成人精品| 国产亚洲午夜精品一区二区久久| 日韩熟女老妇一区二区性免费视频| 下体分泌物呈黄色| 国产欧美日韩综合在线一区二区| 大又大粗又爽又黄少妇毛片口| 久久久久久久精品精品| 两个人的视频大全免费| 狠狠精品人妻久久久久久综合| 亚洲内射少妇av| 女性生殖器流出的白浆| 九色亚洲精品在线播放| 精品久久久久久电影网| 国产精品国产av在线观看| 99九九在线精品视频| 国产精品国产三级专区第一集| 精品亚洲成a人片在线观看| 欧美性感艳星| 亚洲欧美一区二区三区黑人 | 少妇熟女欧美另类| 成人国产av品久久久| 精品视频人人做人人爽| 卡戴珊不雅视频在线播放| 色5月婷婷丁香| 亚洲欧洲精品一区二区精品久久久 | 各种免费的搞黄视频| 在线亚洲精品国产二区图片欧美 | 丰满饥渴人妻一区二区三| av播播在线观看一区| 香蕉精品网在线| 国产精品一国产av| 狠狠精品人妻久久久久久综合| 九九久久精品国产亚洲av麻豆| 精品人妻一区二区三区麻豆| 精品午夜福利在线看| 亚洲欧美色中文字幕在线| 一级片'在线观看视频| 久久久精品免费免费高清| 夫妻午夜视频| 国产伦理片在线播放av一区| 久久青草综合色| 色94色欧美一区二区| videossex国产| 高清黄色对白视频在线免费看| 一级黄片播放器| 欧美日韩视频精品一区| 精品99又大又爽又粗少妇毛片| 国产极品粉嫩免费观看在线 | www.av在线官网国产| 黄色毛片三级朝国网站| 成人免费观看视频高清| 国产av国产精品国产| 精品熟女少妇av免费看| 欧美精品亚洲一区二区| 午夜av观看不卡| 国产深夜福利视频在线观看| 黑人高潮一二区| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 国产av精品麻豆| 下体分泌物呈黄色| 蜜臀久久99精品久久宅男| 三级国产精品片| 哪个播放器可以免费观看大片| 久久久久久人妻| 成人国产av品久久久| 一级二级三级毛片免费看| 九草在线视频观看| 好男人视频免费观看在线| 寂寞人妻少妇视频99o| 狂野欧美激情性bbbbbb| 哪个播放器可以免费观看大片| 国产乱人偷精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 考比视频在线观看| 啦啦啦视频在线资源免费观看| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 成人黄色视频免费在线看| 精品亚洲乱码少妇综合久久| 在线观看三级黄色|