• <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則
    日本欧美国产在线视频| 观看美女的网站| 男女啪啪激烈高潮av片| 精品国产露脸久久av麻豆| 人成视频在线观看免费观看| 91aial.com中文字幕在线观看| 一区在线观看完整版| 少妇的丰满在线观看| 亚洲国产欧美在线一区| 欧美另类一区| 国产一区二区 视频在线| 国产1区2区3区精品| 高清欧美精品videossex| 欧美精品亚洲一区二区| 亚洲国产欧美日韩在线播放| 精品亚洲成a人片在线观看| 精品一品国产午夜福利视频| 熟女av电影| 色94色欧美一区二区| 国产成人精品久久二区二区91 | 日韩av在线免费看完整版不卡| 你懂的网址亚洲精品在线观看| 满18在线观看网站| 超色免费av| 视频在线观看一区二区三区| 久久免费观看电影| freevideosex欧美| 91国产中文字幕| 99香蕉大伊视频| 少妇被粗大猛烈的视频| 亚洲欧美成人综合另类久久久| 性少妇av在线| 亚洲国产av新网站| 久久久久久伊人网av| 精品一区二区免费观看| 有码 亚洲区| 在线观看免费高清a一片| 丝袜在线中文字幕| videosex国产| 成人国产麻豆网| 高清不卡的av网站| 女人高潮潮喷娇喘18禁视频| 肉色欧美久久久久久久蜜桃| √禁漫天堂资源中文www| 五月天丁香电影| 欧美av亚洲av综合av国产av | 日本vs欧美在线观看视频| 综合色丁香网| av电影中文网址| 久久国内精品自在自线图片| 国产免费一区二区三区四区乱码| 国产成人aa在线观看| 天堂俺去俺来也www色官网| www.自偷自拍.com| 亚洲av综合色区一区| 我要看黄色一级片免费的| 丰满乱子伦码专区| 熟女少妇亚洲综合色aaa.| 国产精品一区二区在线不卡| 久久青草综合色| 嫩草影院入口| 少妇人妻精品综合一区二区| 女人被躁到高潮嗷嗷叫费观| 热99久久久久精品小说推荐| 天堂中文最新版在线下载| 欧美日韩一级在线毛片| 欧美精品av麻豆av| 热re99久久国产66热| 侵犯人妻中文字幕一二三四区| 男女午夜视频在线观看| 中文字幕人妻丝袜一区二区 | 在线观看免费日韩欧美大片| 啦啦啦在线观看免费高清www| av国产精品久久久久影院| 777久久人妻少妇嫩草av网站| av线在线观看网站| 欧美xxⅹ黑人| 久久热在线av| 亚洲欧美精品综合一区二区三区 | 十八禁网站网址无遮挡| 人人妻人人澡人人爽人人夜夜| 母亲3免费完整高清在线观看 | 亚洲精品国产av蜜桃| 伦精品一区二区三区| 深夜精品福利| 亚洲成色77777| 国产极品天堂在线| 精品一区二区免费观看| 99久久综合免费| 乱人伦中国视频| 99国产综合亚洲精品| 日韩大片免费观看网站| 国产老妇伦熟女老妇高清| 天天躁狠狠躁夜夜躁狠狠躁| 人成视频在线观看免费观看| 日本爱情动作片www.在线观看| 午夜精品国产一区二区电影| 欧美xxⅹ黑人| 美女xxoo啪啪120秒动态图| 热re99久久精品国产66热6| 国产精品久久久av美女十八| 亚洲视频免费观看视频| 亚洲人成网站在线观看播放| 大香蕉久久网| 亚洲激情五月婷婷啪啪| 亚洲成国产人片在线观看| 国产精品欧美亚洲77777| 国产高清国产精品国产三级| 久久韩国三级中文字幕| 日韩免费高清中文字幕av| 曰老女人黄片| 国产精品久久久久久久久免| 国产成人精品一,二区| 亚洲视频免费观看视频| 国产精品一区二区在线观看99| 伊人亚洲综合成人网| 少妇人妻 视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品一二三| 大香蕉久久成人网| 中文字幕人妻丝袜一区二区 | 亚洲av电影在线观看一区二区三区| 久久久精品区二区三区| 亚洲经典国产精华液单| 各种免费的搞黄视频| 人妻少妇偷人精品九色| 捣出白浆h1v1| 国产精品人妻久久久影院| 9热在线视频观看99| 亚洲欧美成人精品一区二区| 精品一区二区三区四区五区乱码 | 久久久国产一区二区| 亚洲欧美日韩另类电影网站| 伦理电影大哥的女人| 99热网站在线观看| 九色亚洲精品在线播放| 色94色欧美一区二区| 亚洲国产欧美在线一区| 18禁动态无遮挡网站| 黑丝袜美女国产一区| 亚洲综合色网址| 精品少妇一区二区三区视频日本电影 | 久久久国产欧美日韩av| 只有这里有精品99| 最近手机中文字幕大全| 欧美人与善性xxx| av在线app专区| 叶爱在线成人免费视频播放| 久久这里只有精品19| 成人18禁高潮啪啪吃奶动态图| 午夜老司机福利剧场| 美女中出高潮动态图| 在线观看免费日韩欧美大片| 日韩免费高清中文字幕av| 两个人免费观看高清视频| 波多野结衣一区麻豆| 欧美日韩国产mv在线观看视频| 久久精品熟女亚洲av麻豆精品| 久久99精品国语久久久| 欧美人与性动交α欧美软件| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 高清欧美精品videossex| 少妇的逼水好多| 久久精品国产自在天天线| 一本久久精品| 深夜精品福利| 国产精品亚洲av一区麻豆 | 亚洲国产毛片av蜜桃av| 国产成人一区二区在线| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| 久久综合国产亚洲精品| av女优亚洲男人天堂| av在线老鸭窝| 国产白丝娇喘喷水9色精品| 中文欧美无线码| 搡老乐熟女国产| 亚洲欧美清纯卡通| 在线精品无人区一区二区三| 波多野结衣一区麻豆| xxxhd国产人妻xxx| 在线观看免费视频网站a站| 在线看a的网站| 欧美少妇被猛烈插入视频| 欧美日韩亚洲高清精品| 精品人妻偷拍中文字幕| 久久 成人 亚洲| 国产xxxxx性猛交| 两性夫妻黄色片| 美国免费a级毛片| 国产又色又爽无遮挡免| 久久久久视频综合| 自线自在国产av| 亚洲精品aⅴ在线观看| 我的亚洲天堂| 久久久久久久亚洲中文字幕| 黄频高清免费视频| 亚洲综合色网址| 欧美日韩精品网址| 亚洲精品美女久久av网站| 精品久久久久久电影网| 丝袜在线中文字幕| 国产黄色免费在线视频| av不卡在线播放| 久久av网站| 七月丁香在线播放| 美女国产高潮福利片在线看| 午夜免费观看性视频| 免费播放大片免费观看视频在线观看| 日韩免费高清中文字幕av| 国产精品蜜桃在线观看| 亚洲在久久综合| 满18在线观看网站| 18禁裸乳无遮挡动漫免费视频| 少妇猛男粗大的猛烈进出视频| 超碰97精品在线观看| 久久亚洲国产成人精品v| 麻豆av在线久日| 久久久国产一区二区| 日本欧美国产在线视频| 丝袜喷水一区| 亚洲欧美精品自产自拍| 2018国产大陆天天弄谢| 国产1区2区3区精品| 三级国产精品片| 精品亚洲乱码少妇综合久久| 午夜免费观看性视频| 肉色欧美久久久久久久蜜桃| 十分钟在线观看高清视频www| 9191精品国产免费久久| 免费黄网站久久成人精品| 亚洲中文av在线| 久久久精品国产亚洲av高清涩受| 日日爽夜夜爽网站| 国产精品久久久久久久久免| 亚洲精品美女久久久久99蜜臀 | 亚洲国产看品久久| 久久免费观看电影| 男女国产视频网站| 亚洲av综合色区一区| 岛国毛片在线播放| 亚洲国产av影院在线观看| 99久久精品国产国产毛片| 男女无遮挡免费网站观看| 国产黄色免费在线视频| 欧美少妇被猛烈插入视频| 久久久久久免费高清国产稀缺| 亚洲精品美女久久av网站| 国产激情久久老熟女| 欧美成人午夜免费资源| 五月伊人婷婷丁香| 免费看av在线观看网站| 国产极品天堂在线| 一二三四在线观看免费中文在| 亚洲av电影在线进入| 久久久久人妻精品一区果冻| 日本欧美视频一区| 国产精品亚洲av一区麻豆 | 捣出白浆h1v1| 成人影院久久| h视频一区二区三区| 日本wwww免费看| 久久久久久免费高清国产稀缺| 亚洲av欧美aⅴ国产| 久久午夜福利片| 深夜精品福利| 亚洲国产精品一区三区| 考比视频在线观看| 国产男女内射视频| 一级爰片在线观看| 七月丁香在线播放| 久久精品国产鲁丝片午夜精品| 久久99精品国语久久久| 久久久久国产精品人妻一区二区| 国产免费福利视频在线观看| 丰满少妇做爰视频| 夫妻性生交免费视频一级片| 亚洲色图综合在线观看| 91精品三级在线观看| 欧美另类一区| 久久久久久久亚洲中文字幕| 色吧在线观看| 男女午夜视频在线观看| 日本vs欧美在线观看视频| 飞空精品影院首页| 男女啪啪激烈高潮av片| 亚洲一区中文字幕在线| 最近最新中文字幕大全免费视频 | 人人妻人人添人人爽欧美一区卜| 蜜桃在线观看..| 国产男人的电影天堂91| 国产免费又黄又爽又色| 国产一区二区在线观看av| 国产精品久久久久久精品古装| 美女大奶头黄色视频| 日韩免费高清中文字幕av| 女人久久www免费人成看片| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久精品电影小说| 欧美日韩综合久久久久久| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区| 久久精品久久精品一区二区三区| 深夜精品福利| 天堂8中文在线网| 久久久欧美国产精品| 高清av免费在线| 老女人水多毛片| 久久久久久伊人网av| 妹子高潮喷水视频| 如日韩欧美国产精品一区二区三区| 黄色视频在线播放观看不卡| 美女大奶头黄色视频| 不卡av一区二区三区| 欧美最新免费一区二区三区| 9191精品国产免费久久| 色网站视频免费| 午夜福利视频精品| 女性生殖器流出的白浆| av卡一久久| 十八禁网站网址无遮挡| av在线播放精品| 亚洲成av片中文字幕在线观看 | 亚洲第一av免费看| xxx大片免费视频| 亚洲少妇的诱惑av| 日韩大片免费观看网站| 日日啪夜夜爽| 男女边摸边吃奶| av国产久精品久网站免费入址| 黄色配什么色好看| 欧美日韩视频高清一区二区三区二| 另类亚洲欧美激情| 中国国产av一级| 人妻 亚洲 视频| 午夜激情久久久久久久| h视频一区二区三区| 啦啦啦在线观看免费高清www| 波多野结衣av一区二区av| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 汤姆久久久久久久影院中文字幕| 一本大道久久a久久精品| 国产在线视频一区二区| 亚洲 欧美一区二区三区| 中文字幕人妻丝袜一区二区 | 国产一级毛片在线| av免费观看日本| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 免费看av在线观看网站| 美女主播在线视频| 中文字幕亚洲精品专区| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三区在线| 一区二区av电影网| 晚上一个人看的免费电影| 午夜福利网站1000一区二区三区| 国产成人免费无遮挡视频| 欧美变态另类bdsm刘玥| 色哟哟·www| 亚洲国产精品国产精品| 亚洲人成电影观看| 我要看黄色一级片免费的| 久久久久精品人妻al黑| 美女xxoo啪啪120秒动态图| 在线观看一区二区三区激情| 搡女人真爽免费视频火全软件| www日本在线高清视频| 国产成人精品福利久久| 欧美精品av麻豆av| 国产精品麻豆人妻色哟哟久久| 免费观看a级毛片全部| 国产精品嫩草影院av在线观看| 亚洲经典国产精华液单| 黄片播放在线免费| 成年动漫av网址| 一本—道久久a久久精品蜜桃钙片| 丰满饥渴人妻一区二区三| 久久精品国产亚洲av高清一级| 看免费av毛片| 考比视频在线观看| 99久久中文字幕三级久久日本| 国产精品久久久久久精品古装| 亚洲精品自拍成人| 一本大道久久a久久精品| av有码第一页| 久久久精品94久久精品| 精品亚洲成a人片在线观看| 欧美日韩精品成人综合77777| 久久精品国产综合久久久| 欧美黄色片欧美黄色片| 777米奇影视久久| 狠狠精品人妻久久久久久综合| 天天躁夜夜躁狠狠久久av| 久久午夜福利片| 亚洲,欧美精品.| 中文字幕色久视频| 午夜精品国产一区二区电影| 亚洲av中文av极速乱| 亚洲精品国产av成人精品| 精品人妻熟女毛片av久久网站| 日本猛色少妇xxxxx猛交久久| 亚洲美女搞黄在线观看| 一区二区三区精品91| 99久久精品国产国产毛片| 国产在线免费精品| 搡老乐熟女国产| 日韩视频在线欧美| 日韩一卡2卡3卡4卡2021年| 国产精品蜜桃在线观看| 超色免费av| 国产精品免费视频内射| 熟女电影av网| 丝袜美足系列| 捣出白浆h1v1| 成人影院久久| 秋霞伦理黄片| 成人漫画全彩无遮挡| 中文字幕人妻丝袜一区二区 | 亚洲精品国产av蜜桃| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av | 午夜激情av网站| 亚洲一区二区三区欧美精品| 飞空精品影院首页| 一级毛片 在线播放| 欧美在线黄色| 丰满迷人的少妇在线观看| 久久鲁丝午夜福利片| 好男人视频免费观看在线| 黄频高清免费视频| 精品少妇一区二区三区视频日本电影 | 欧美日韩一区二区视频在线观看视频在线| 国产免费视频播放在线视频| 欧美激情 高清一区二区三区| 久久青草综合色| 夫妻性生交免费视频一级片| 国产老妇伦熟女老妇高清| 久久久亚洲精品成人影院| 国产精品国产三级国产专区5o| 制服诱惑二区| 涩涩av久久男人的天堂| 97在线视频观看| 国产成人一区二区在线| 9191精品国产免费久久| 男的添女的下面高潮视频| 天堂俺去俺来也www色官网| 26uuu在线亚洲综合色| 精品亚洲成国产av| av一本久久久久| 亚洲在久久综合| 色播在线永久视频| 亚洲av日韩在线播放| 日韩一区二区三区影片| 曰老女人黄片| 最新的欧美精品一区二区| 欧美亚洲日本最大视频资源| 女人精品久久久久毛片| 赤兔流量卡办理| 校园人妻丝袜中文字幕| 久久久久久久久久人人人人人人| 9热在线视频观看99| 久久精品亚洲av国产电影网| 亚洲av在线观看美女高潮| 亚洲综合色网址| 欧美激情 高清一区二区三区| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜制服| 韩国av在线不卡| 美女主播在线视频| 国产日韩欧美在线精品| 亚洲情色 制服丝袜| 久久女婷五月综合色啪小说| 女性被躁到高潮视频| 在线观看免费日韩欧美大片| 亚洲国产精品国产精品| 欧美另类一区| 女人久久www免费人成看片| av不卡在线播放| 国产极品天堂在线| 亚洲国产精品国产精品| www.熟女人妻精品国产| 黄色配什么色好看| 中文乱码字字幕精品一区二区三区| 免费高清在线观看视频在线观看| 99久久人妻综合| 高清av免费在线| 精品少妇黑人巨大在线播放| 高清欧美精品videossex| 这个男人来自地球电影免费观看 | 一二三四在线观看免费中文在| 成年美女黄网站色视频大全免费| 国产乱来视频区| 精品国产一区二区久久| 国产精品香港三级国产av潘金莲 | 夫妻午夜视频| 五月开心婷婷网| 中文乱码字字幕精品一区二区三区| 黄片无遮挡物在线观看| 久久久a久久爽久久v久久| 成人影院久久| 大陆偷拍与自拍| 一二三四在线观看免费中文在| 欧美日韩一级在线毛片| 日韩电影二区| 亚洲精品av麻豆狂野| 久久精品国产亚洲av高清一级| 亚洲精品在线美女| 中国三级夫妇交换| 18禁国产床啪视频网站| 一级毛片我不卡| 久久久亚洲精品成人影院| 狠狠精品人妻久久久久久综合| av网站在线播放免费| 在线精品无人区一区二区三| 欧美 亚洲 国产 日韩一| 国产极品天堂在线| videos熟女内射| 啦啦啦在线观看免费高清www| 在线精品无人区一区二区三| 99久国产av精品国产电影| 国产精品偷伦视频观看了| 精品亚洲成a人片在线观看| 国产精品一区二区在线不卡| 国产精品av久久久久免费| 狂野欧美激情性bbbbbb| 国产一区二区在线观看av| 国产精品国产av在线观看| 大片免费播放器 马上看| 男女啪啪激烈高潮av片| 久久精品国产a三级三级三级| 久久久精品94久久精品| 免费高清在线观看日韩| 人妻系列 视频| av免费在线看不卡| a级片在线免费高清观看视频| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 一本久久精品| 亚洲精品aⅴ在线观看| 如何舔出高潮| 国产男女内射视频| 看十八女毛片水多多多| 少妇 在线观看| 美女国产视频在线观看| 亚洲精华国产精华液的使用体验| 一级毛片我不卡| 国产黄频视频在线观看| 国产精品久久久av美女十八| 久久99热这里只频精品6学生| 亚洲欧美精品自产自拍| 欧美日韩亚洲高清精品| 亚洲欧洲日产国产| 国产高清不卡午夜福利| 一边摸一边做爽爽视频免费| av片东京热男人的天堂| 黄片小视频在线播放| 99久久精品国产国产毛片| 欧美日韩视频精品一区| 巨乳人妻的诱惑在线观看| 国产av国产精品国产| 男女下面插进去视频免费观看| 精品99又大又爽又粗少妇毛片| 免费av中文字幕在线| 久久国产精品大桥未久av| 建设人人有责人人尽责人人享有的| 色吧在线观看| 亚洲精品美女久久av网站| 人妻少妇偷人精品九色| 极品人妻少妇av视频| 国产成人一区二区在线| 老司机影院成人| 日韩中字成人| 少妇的逼水好多| 亚洲综合精品二区| 69精品国产乱码久久久| 久久久久久久亚洲中文字幕| 亚洲精品av麻豆狂野| 亚洲精品视频女| 一个人免费看片子| 免费观看性生交大片5| 电影成人av| 精品亚洲成a人片在线观看| 校园人妻丝袜中文字幕| 电影成人av| 国产野战对白在线观看| 最近中文字幕高清免费大全6| 成人毛片60女人毛片免费| 国产深夜福利视频在线观看| 精品第一国产精品| 久久99一区二区三区| 韩国精品一区二区三区| 精品第一国产精品| 国产成人精品久久二区二区91 | 午夜福利乱码中文字幕| 亚洲国产最新在线播放| 制服诱惑二区| 久久久精品免费免费高清| 午夜福利影视在线免费观看| 免费观看a级毛片全部| 少妇被粗大的猛进出69影院| 国语对白做爰xxxⅹ性视频网站| 国产熟女午夜一区二区三区| 少妇的逼水好多| 2021少妇久久久久久久久久久| 免费播放大片免费观看视频在线观看| 熟女av电影| 国产精品久久久久久精品电影小说| 少妇人妻久久综合中文| 捣出白浆h1v1| 少妇熟女欧美另类| 又粗又硬又长又爽又黄的视频| 国产精品一区二区在线观看99|