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

    基于無損檢測技術(shù)的濕地松生長及材性性狀遺傳變異分析*

    2017-07-18 12:10:14張帥楠欒啟福姜景民
    林業(yè)科學(xué) 2017年6期
    關(guān)鍵詞:濕地松遺傳變異遺傳力

    張帥楠 欒啟福 姜景民

    (中國林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所 浙江省林木育種技術(shù)研究重點實驗室 杭州 311400)

    ?

    基于無損檢測技術(shù)的濕地松生長及材性性狀遺傳變異分析*

    張帥楠 欒啟福 姜景民

    (中國林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所 浙江省林木育種技術(shù)研究重點實驗室 杭州 311400)

    【目的】 基于無損檢測技術(shù)研究濕地松活立木基本密度(ρ)、彈性模量(MOE)和生長量的遺傳變異規(guī)律及其遺傳和表型相關(guān)關(guān)系,在保證生長量遺傳改良的同時加強(qiáng)木材物理力學(xué)性質(zhì)的選育力度,為速生優(yōu)質(zhì)濕地松人工林建設(shè)提供種質(zhì)基礎(chǔ)?!痉椒ā?以浙江杭州長樂林場22年生濕地松半同胞家系測定林為研究對象,采用木材阻力儀Pilodyn和應(yīng)力波測量儀Hitman等無損檢測儀器測定2 000多株濕地松活立木ρ、MOE等物理力學(xué)性狀,采用測桿等標(biāo)準(zhǔn)測具測定對應(yīng)單株的枝下高、冠幅、分枝大小、分枝角以及樹高和胸徑等形質(zhì)和生長性狀,并通過ASReml-R軟件混合線性模型的限制性極大似然估計法(REML)對各性狀進(jìn)行方差分析、遺傳力估算和相關(guān)性分析?!窘Y(jié)果】 1)彈性模量相對值MOEP與應(yīng)力波速v之間呈極顯著(P<0.001)遺傳相關(guān)(R=0.947 0±0.020 1)和表型相關(guān)(R=0.948 0±0.002 7),與基本密度相對值ρP之間呈顯著(P<0.01)遺傳(R=-0.447 0±0.154 0)和表型(R=-0.538 0±0.019 1)相關(guān),以v對MOEP進(jìn)行間接選擇,其有效性達(dá)到了97.42%。2)材性性狀的遺傳力在0.292 0~0.305 0之間,生長性狀遺傳力在0.062 5~0.216 0之間; 以1%的選擇強(qiáng)度計算,材積和MOE能獲得30%的遺傳增益。3)樹高、胸徑、材積、冠幅和分枝大小兩兩之間均呈顯著的(P<0.05)遺傳正相關(guān)關(guān)系,枝下高與樹高正相關(guān),分枝角與分枝大小負(fù)相關(guān),但相關(guān)不顯著,而胸徑、樹高、枝下高、材積、冠幅和分枝大小兩兩之間均具有極顯著的表型正相關(guān)性; 材性和生長性狀之間普遍存在著顯著的表型相關(guān)關(guān)系(P<0.05),但遺傳相關(guān)不顯著。【結(jié)論】 僅利用應(yīng)力波速v即能快速評估群體內(nèi)每木彈性模量(MOE)相對大小,比傳統(tǒng)上僅用木材基本密度(ρ)評估MOE效果更好;ρ和MOE受中等程度的遺傳控制,其遺傳力明顯高于生長性狀; 材性與生長性狀之間遺傳相關(guān)不顯著,能夠?qū)崿F(xiàn)兩類性狀的同步遺傳改良。

    濕地松; 生長性狀; 木材密度; 彈性模量; 無損評估; 遺傳分析

    濕地松(Pinuselliottii)是我國南方紅壤低丘地區(qū)引種成功的速生用材、采脂樹種,其成林早、生長快、樹干通直、產(chǎn)脂量高(朱志淞等, 1993),當(dāng)前人工林面積達(dá)到300萬hm2。由于我國濕地松人工林主要位于亞熱帶季風(fēng)氣候區(qū),容易受到起伏不平的地形和夏季臺風(fēng)、冬季低溫冰凍影響,在極端氣候條件下彎折、倒伏風(fēng)險較高(欒啟福等, 2008),因此針對我國造林地區(qū)特點,加強(qiáng)濕地松木材強(qiáng)度性狀選育對于我國濕地松人工林健康持續(xù)發(fā)展具有積極意義。

    木材基本密度(basic wood density,ρ)、彈性模量(modulus of elasticity,MOE)等性狀是木材強(qiáng)度的重要測量指標(biāo),也是建筑材、紙漿材重要選擇性狀(Zobeletal., 1989)。然而受限于傳統(tǒng)測量方法,在濕地松材性性狀改良進(jìn)程中,尤其是涉及到破壞性取樣的彈性模量等力學(xué)性狀每次測量的樣品量較少(吳際友等, 2000; 徐有明等, 2001; 駱秀琴等, 2003)。從統(tǒng)計角度而言其結(jié)果誤差較大,這增大了遺傳選擇錯誤或遺漏的概率。因此,開發(fā)利用操作簡單方便、破壞性小且結(jié)果可靠的無損評估技術(shù),以大樣本試驗材料為基礎(chǔ),研究彈性模量等材性性狀的遺傳變異規(guī)律,對于提高遺傳測定的選擇精度與效率,促進(jìn)遺傳改良進(jìn)程具有重大意義。近年來基于力學(xué)、聲學(xué)的無損檢測技術(shù)為大樣本量估測活立木基本密度和彈性模量奠定了基礎(chǔ),逐漸成為國際上相關(guān)性狀測量、選擇的常規(guī)手段(Whiteetal., 2007)。其中已成功應(yīng)用于活立木彈性模量無損評估的是基于應(yīng)力波速法的AV技術(shù),它通過測量應(yīng)力波在開始探頭和結(jié)束探頭之間的傳播時間來測得應(yīng)力波速v,該數(shù)值與立木強(qiáng)度具有顯著的相關(guān)性,可以作為活立木強(qiáng)度的檢測工具; Pilodyn是應(yīng)用最為普遍的一種基于阻力的活立木密度無損快速檢測儀器(Lenzetal., 2013; 易敏等, 2014; Chenetal., 2015)。

    基于這2種簡潔快速測量方法,本研究著重加強(qiáng)濕地松彈性模量遺傳變異研究,在國內(nèi)首次實現(xiàn)了2 000多株22年生濕地松半同胞家系全林每木基本密度、彈性模量的測量,同時按照傳統(tǒng)方法測量了樹高、胸徑,并估測了其冠幅、枝下高和分枝大小等性狀。通過無偏估計濕地松遺傳群體主要物理力學(xué)性狀、形態(tài)性狀和生長量性狀之間的遺傳及其相關(guān)特性,為濕地松大群體活立木多性狀綜合評價及聯(lián)合改良提供依據(jù),這對于我國濕地松人工林健康持續(xù)發(fā)展具有積極意義。

    1 材料與方法

    1.1 試驗林概況

    試驗地位于浙江省杭州市余杭區(qū)長樂林場,30°27′ N,119°48′E,亞熱帶濕潤季風(fēng)性氣候,年均溫16.1 ℃,年均降水量1 399 mm,屬低山緩坡,土壤以紅壤為主。試驗林營建于1994年春,為濕地松子代測定林,含有158個半同胞家系,試驗布置采取完全隨機(jī)區(qū)組設(shè)計,6株單行小區(qū),6個重復(fù),株行距2 m×3 m。根據(jù)培育規(guī)程,在林分10年生時對其進(jìn)行隔株間伐,現(xiàn)保留單株共計2 500余株。

    1.2 數(shù)據(jù)測定方法

    1)生長性狀測定: 2015年12月,在林分停止生長期,采用測桿等標(biāo)準(zhǔn)測具對胸徑、樹高、枝下高、冠幅、分枝性狀等主要生長性狀進(jìn)行全林每木調(diào)查,其中分枝大小和分枝角均按3級評測標(biāo)準(zhǔn)(1表示分枝較小或分枝角為0°~30°,3表示分枝較大或分枝角為60°~90°),單株材積按V=0.375H×DBH2估算(欒啟福等, 2011)。

    2)材性無損測定: 利用Pilodyn(6J,PROCEQ,Switzerland)測定獲得樣木Pilodyn值(P)(朱景樂等, 2008),P與木材基本密度(ρ)呈顯著線性負(fù)相關(guān)關(guān)系,因此本文中將P作為基本密度相對值ρP(欒啟福等, 2011); 利用活立木應(yīng)力波速測量專用工具Hitman(Fibre-gen,Christchurch,New Zealand)對選擇的樣木進(jìn)行測量獲得波速v(Fibre-gen, 2004; Carteretal., 2005),并利用彈性模量(MOE)、ρ和v之間的關(guān)系式MOE=ρ×v2(Rossetal., 1988)計算出彈性模量相對值MOEP(MOEP=1/P*v2)(Chenetal., 2015)。

    1.3 統(tǒng)計分析方法

    1)描述性數(shù)據(jù)分析利用R軟件(v 3.1.3,CRAN,http://cran.r-project.org)中pastecs程序包的stat.desc()函數(shù)進(jìn)行(林元震等, 2014)。

    2)通過ASReml-R軟件,利用其混合線性模型的限制性極大似然估計法(REML)對各性狀的數(shù)據(jù)進(jìn)行方差分析、遺傳力估算和相關(guān)性分析(Butleretal., 2009; 林元震等, 2014)。該線性模型如下:

    (1)

    式中:Yijk表示第ijk株樹的性狀觀測值,μ表示該性狀觀察值的總體平均值,Ri表示重復(fù)效應(yīng),F(xiàn)j表示家系效應(yīng),RFij表示重復(fù)和家系交互效應(yīng),eijk表示誤差效應(yīng),其中Ri作為模型固定效應(yīng),F(xiàn)j和RFij作為模型隨機(jī)效應(yīng)。

    性狀狹義遺傳力估算公式為:

    (2)

    性狀遺傳相關(guān)分析計算公式為:

    (3)

    性狀表型相關(guān)分析計算公式為:

    (4)

    性狀指標(biāo)遺傳增益計算公式(Namkoong, 1981)為:

    (5)

    相關(guān)性狀間間接選擇有效性計算公式(Whiteetal., 2007)為:

    (6)

    式中:rg表示性狀x、y之間的遺傳相關(guān)關(guān)系,hx、hy分別表示性狀x、y的遺傳力的平方根,其中y是目標(biāo)性狀。

    2 結(jié)果與分析

    2.1 生長和材性性狀各指標(biāo)的表型變異分析

    2 000余株濕地松的生長和材性性狀的基本概況見表1。22年生濕地松的胸徑、樹高、枝下高、材積和冠幅分別為20.45 cm,16.35 m,9.55 m,0.277 m3和3.26 m,枝下高約占整個樹高的2/3左右,即濕地松冠層較小,單株樹木建筑材質(zhì)量和產(chǎn)量較高。其中胸徑、樹高、枝下高和冠幅變異系數(shù)(14.15%~26.99%)較小,材積變異系數(shù)達(dá)到了50%以上,說明材積觀測值變異幅度大,具有較好的選擇效果。彈性模量相對值MOEP的變異系數(shù)(31.60%)遠(yuǎn)大于基本密度相對值ρP(10.06%),這表明了與ρ相比,MOE可能具有更豐富的遺傳變異性。

    2.2 生長和材性性狀的遺傳力和選擇增益

    表2顯示了通過方差組分估算出的各性狀狹義遺傳力大小。材性性狀指標(biāo)的遺傳力為0.292 0~0.309 0,明顯高于生長性狀的遺傳力(0.062 5~0.216 0),即ρ和MOE受相對較高的遺傳控制,而生長性狀更易受環(huán)境影響,適宜的生長環(huán)境或栽培管理技術(shù)可以獲得一定的生長增益。表2還顯示了以1%入選率對各性狀進(jìn)行直接選擇時,各性狀所能獲得的遺傳增益大小(其中,ρP的遺傳增益計算時以ρP值最小的前1%個個體作為入選樣本)??梢钥闯?,材積(31.47%)和MOEP(29.11%)能獲得較大增益,具有良好的選擇效應(yīng);ρP的遺傳增益(14.57%)較小,但是ρ的微小變化對最終木材產(chǎn)品的變化和效益仍是值得重視的(Zobeletal., 1972),因此濕地松ρ也具有較大的選擇效益。

    表1 濕地松半同胞家系生長和材性性狀各指標(biāo)基本概況Tab. 1 The basic situations of growth and wood properties of half-sib families of slash pine

    表2 濕地松半同胞家系生長和材性性狀的遺傳力及其 選擇增益①Tab. 2 Narrow-sense heritabilities and genetic gains of growth and wood properties of half-sib families of slash pine

    ① 性狀全稱見表1; 遺傳增益Ga是以1%入選率計算的。Character details see Tab.1; The genetic gainGais calculated under the selection rate 1%.

    2.3 生長和材性性狀間的相關(guān)性分析

    對濕地松生長和材性性狀各指標(biāo)進(jìn)行遺傳和表型相關(guān)分析(表3)。就生長性狀而言,胸徑、材積、冠幅和分枝大小兩兩之間均呈顯著的遺傳正相關(guān)關(guān)系(R=0.440 0~0.995 0,P<0.05); 樹高與胸徑(R=0.753 0 ± 0.134 0)、材積(R=0.865 0 ± 0.085 7)之間呈極顯著(P<0.001)遺傳正相關(guān)關(guān)系; 枝下高僅與樹高呈顯著的遺傳正相關(guān)性(R=0.436 0 ± 0.236 0,P<0.05)。分枝角除與分枝大小呈現(xiàn)極顯著表型負(fù)相關(guān)關(guān)系(R=-0.147 0 ± 0.025 1,P<0.001)外,與其他性狀間均無顯著表型或遺傳相關(guān)關(guān)系; 而胸徑、樹高、枝下高、材積、冠幅和分枝大小兩兩之間均具有極顯著(P<0.001)的表型正相關(guān)性。

    同時,相關(guān)分析結(jié)果還顯示材性和生長性狀之間遺傳相關(guān)不顯著,但普遍存在著顯著的表型正相關(guān),顯示材性和生長性狀的遺傳改良可同步進(jìn)行,培育措施與生長和材性表現(xiàn)有關(guān)。對于材性性狀無損測定指標(biāo),ρP與v之間遺傳相關(guān)不顯著,僅表型相關(guān)(R=-0.264 0 ± 0.025 1,P<0.001)達(dá)到了顯著性水平; MOEP和v之間呈極顯著(P<0.001)的遺傳和表型正相關(guān)關(guān)系,相關(guān)系數(shù)分別達(dá)到0.947 0和0.948 0,以v對MOEP進(jìn)行間接選擇,其有效性達(dá)到了97.42%。MOEP和ρP之間呈顯著(P<0.01)遺傳負(fù)相關(guān)(R=-0.447 0±0.154 0)和極顯著(P<0.001)表型負(fù)相關(guān)(R=-0.538 0±0.019 1),即MOE和ρ之間顯著正相關(guān),但其與ρ的相關(guān)系數(shù)遠(yuǎn)小于與v的相關(guān)系數(shù)。

    表3 濕地松半同胞家系各性狀間遺傳相關(guān)和表型相關(guān)性①Tab. 3 Genetic and phenotypic correlations between traits of growth and wood properties of half-sib families of slash pine

    ①性狀全稱見表1。上三角為遺傳相關(guān),下三角為表型相關(guān),括號中為該相關(guān)系數(shù)的標(biāo)準(zhǔn)誤;*,**,***分別表示相關(guān)顯著性水平為P<0.05,P<0.01,P<0.001。Character details see Tab. 1. Upper triangle is genetic correlations, and lower triangle is phenotypic correlations; Standard errors are given in parenthesis.*,**and***representP<0.05,P<0.01 andP<0.001, respectively.

    3 討論

    彈性模量(MOE)是林木建筑用材最重要的力學(xué)性質(zhì)(Raymondetal., 2004; Lasserreetal., 2007),然而受傳統(tǒng)評估方法等限制,不同地區(qū)對濕地松MOE的改良均處于非優(yōu)先地位。美洲濕地松人工林林地地勢平坦,集約經(jīng)營強(qiáng)度高,濕地松遺傳改良主要集中在材積和抗銹病選擇育種上(Whiteetal., 2004),兼顧密度(ρ)等材性性狀。國內(nèi)對濕地松木材材性方面的相關(guān)研究也集中在生長性狀,兼顧密度性狀,對MOE等力學(xué)性狀的遺傳分析研究鮮有報道。單方面對生長性狀的追求則被認(rèn)為會嚴(yán)重影響木材的質(zhì)量(Kumar, 2004),不利于木材質(zhì)量的提高。特別是我國南方濕地松人工林林地主要是低山丘陵,且集約化經(jīng)營程度不高,容易受到起伏不平的地形和夏季臺風(fēng)、冬季低溫冰凍影響,在極端氣候條件下彎折、倒伏風(fēng)險較高(欒啟福等, 2008)。對ρ、MOE等影響木材強(qiáng)度的性狀的改良對于濕地松木材質(zhì)量提高及其樹種的健康持續(xù)發(fā)展具有積極意義。

    本試驗基于全林單株測定,采用無損評估方法,在國內(nèi)首次實現(xiàn)了大規(guī)模濕地松活立木生長和材性的遺傳變異分析。遺傳力分析結(jié)果表明,濕地松ρ和MOE狹義遺傳力約為0.30,受中等程度的遺傳控制,與Chen等(2015)基于無損評估方法對5 000余株歐洲云杉(Piceaabies)材性遺傳分析結(jié)果較一致; 而明顯高于駱秀琴等(2003)基于常規(guī)檢測方法得出的研究結(jié)果(濕地松彈性模量廣義遺傳力為0.108 2),這可能是由于在其研究中所取樣木少(18個種源,共計90個單株)而導(dǎo)致的統(tǒng)計結(jié)果誤差大造成的,而且本文中估算MOE的狹義遺傳力標(biāo)準(zhǔn)誤較低(SE=0.078 3),從側(cè)面反映了大樣本量下估算遺傳參數(shù)的精確度和可靠性更高。其中,研究結(jié)果得出對MOE進(jìn)行直接選擇所能獲得的遺傳增益較大,這可能是由于此前樹種改良工作中受研究條件限制沒有對MOE進(jìn)行直接選育,造成試驗林內(nèi)各家系或單株間MOE差異較大引起的。這表明急需加強(qiáng)濕地松MOE的遺傳改良研究,加快濕地松材性選育工作進(jìn)程,摒棄MOE表現(xiàn)不良的單株或家系,進(jìn)而提高濕地松MOE或材性性狀的整體水平,為濕地松建筑用材的選育奠定基礎(chǔ)。

    同時,了解生長與材性性狀的遺傳相關(guān)性有助于速生優(yōu)質(zhì)林木的選育,然而不同樹種或研究方法對生長與材性性狀的相關(guān)性研究結(jié)果不一致(Fries, 1986; 姜景民等, 1996; Ukrainetzetal., 2008; 劉青華等, 2009),因此一直以來這都是育種工作者研究的重點。本文中相關(guān)性分析結(jié)果表明濕地松材性和生長性狀之間遺傳相關(guān)不顯著,這與Fries(1986)和Wang等(1999)對小干松(Pinuscontortassp.latifolia)以及吳際友等(2010)對濕地松生長和材性性狀的相關(guān)性研究結(jié)果相同,可以實現(xiàn)二者的同步改良。同時,濕地松材性和生長性狀之間還普遍存在著顯著的表型正相關(guān)性,說明在特定的環(huán)境條件影響下,材性和生長性狀之間能夠呈現(xiàn)一定的正相關(guān)關(guān)系,這有利于二者的同步選擇,即在實際生產(chǎn)上能夠?qū)崿F(xiàn)生長量大且材質(zhì)優(yōu)良濕地松家系或單株的選育。

    相關(guān)性分析還發(fā)現(xiàn)應(yīng)力波速(v)和MOE相對值(MOEP)之間呈高度的極顯著遺傳和表型正相關(guān),以v對MOEP進(jìn)行間接選擇,其有效性達(dá)到了97.42%。這表明在比較同一測定群體MOE值相對大小時直接比較v值即可,即僅利用AV技術(shù)即能快速評估MOE相對值,這大大減少了同一群體內(nèi)MOE性狀評估與選擇的工作量。然而,在此遺傳變異分析研究過程中,所取樣本量越大,研究所需的人力、物力就相對越多,因此,比較分析生長、材性等性狀在不同取樣規(guī)模下遺傳參數(shù)估算的差異性,論證取樣規(guī)模對其性狀遺傳參數(shù)的影響,進(jìn)而確定各性狀遺傳參數(shù)評估所需的最適樣本規(guī)模便成為了本研究進(jìn)一步探討的關(guān)鍵問題。

    4 結(jié)論

    本文通過Pilodyn和應(yīng)力波速技術(shù)等無損方法,以全林單株測定為基礎(chǔ),大規(guī)模評估了濕地松半同胞家系活立木的基本密度(ρ)和彈性模量(MOE)等材性性狀,并通過遺傳和相關(guān)分析,揭示了濕地松生長和材性的遺傳變異及其相關(guān)關(guān)系。

    1) 22年生濕地松的胸徑、樹高、枝下高、材積和冠幅等主要生長性狀分別為20.45 cm、16.35 m、9.55 m、0.277 m3和3.26 m,枝下高約占整個樹高的2/3左右,即濕地松冠層較小,有利于建筑材產(chǎn)量提高。材積觀測值變異幅度大,具有較好的選擇效果。與ρ相比,MOE具有更豐富的遺傳變異性。

    2)ρ和MOE受中等程度的遺傳控制,其遺傳力明顯高于生長性狀,這說明相較于材性性狀,生長性狀更易受環(huán)境影響。通過一定的選擇強(qiáng)度,材積、ρ和MOE能獲得較大增益,具有良好的選擇效果。同時,本研究結(jié)果還顯示材性和生長性狀之間僅存在著顯著的表型相關(guān),能夠?qū)崿F(xiàn)兩類性狀同步遺傳改良。

    3) 應(yīng)力波速(v)和MOE相對值(MOEP)之間呈高度的遺傳和表型正相關(guān)關(guān)系,相關(guān)系數(shù)均達(dá)到了極顯著水平,僅利用v即能快速評估群體內(nèi)每木MOE相對大小,該方法比傳統(tǒng)上僅用ρ評估木材MOE準(zhǔn)確率高。

    姜景民,孫海菁,劉昭息. 1996. 火炬松紙漿材優(yōu)良家系多性狀選擇. 林業(yè)科學(xué)研究,9(5): 455-460.

    (Jiang J M, Sun H Q, Liu Z X. 1996. Multitraits selection of loblolly pine families for pulpwood. Forest Research, 9(5): 455-460.[in Chinese])

    林元震,陳曉陽. 2014. R與ASReml-R統(tǒng)計分析教程. 北京: 中國林業(yè)出版社, 36-251.

    (Lin Y Z, Chen X Y.2014. R & ASReml-R statistical analysis tutorial. Beijing: China Forestry Publishing House, 36-251. [in Chinese])

    劉青華,金國慶,張 蕊,等. 2009.24年生馬尾松生長、形質(zhì)和木材基本密度的種源變異與種源區(qū)劃. 林業(yè)科學(xué),45(10): 55-61.

    (Liu Q H, Jin G Q, Zhang R,etal. 2009. Provenance variation in growth, stem-form and wood density of Masson pine at 24-year-old and the provenance division. Scientia Silvae Sinicae, 45(10): 55-61. [in Chinese])

    欒啟福,姜景民,張建忠,等. 2011. 火炬松×加勒比松F1代生長、樹干通直度和基本密度遺傳和配合力分析. 林業(yè)科學(xué),47(3): 178-183.

    (Luan Q F, Jiang J M, Zhang J Z,etal. 2011. Estimation of heritability and combining ability for growth,stem-straightness and wood density of the F1generation ofPinustaeda×P.caribaea. Scientia Silvae Sinicae, 47(3): 178-183. [in Chinese])

    欒啟福,盧 萍,肖復(fù)明,等. 2008. 雨雪冰凍重災(zāi)區(qū)濕地松受害情況調(diào)查及其原因初步分析. 林業(yè)科學(xué),44(11): 50-54.

    (Luan Q F, Lu P, Xiao F M,etal. 2008. Investigation on the damage ofPinuselliottiiin the freezing rain and snow area and the analysis on the reason. Scientia Silvae Sinicae, 44(11): 50-54. [in Chinese])

    駱秀琴,姜笑梅,殷亞方,等. 2003. 濕地松15個家系木材材性遺傳變異及優(yōu)良家系評估. 林業(yè)科學(xué)研究,16(6): 694-699.

    (Luo X Q, Jiang X M, Yin Y F,etal. 2003. Genetic variation and comprehensive assessment in wood properties of 15 families ofPinuselliottii. Forest Research, 16(6): 694-699. [in Chinese])

    吳際友,李志輝,龍應(yīng)忠,等. 2010. 濕地松全同胞家系主要經(jīng)濟(jì)性狀的遺傳變異與選擇研究. 中南林業(yè)科技大學(xué)學(xué)報,30(8): 1-4.

    (Wu J Y, Li Z H, Long Y Z,etal. 2010. Genetic variation and selection of main economic characteristics of slash pine full-sib family. Journal of Central South University of Forestry & Technology, 30(8): 1-4. [in Chinese])

    吳際友,龍應(yīng)忠,余格非,等. 2000. 濕地松半同胞家系主要經(jīng)濟(jì)性狀的遺傳分析及聯(lián)合選擇. 林業(yè)科學(xué),36(1): 56-61.

    (Wu J Y, Long Y Z, Yu G F,etal. 2000. Genetics analysis and combined selection of main economic characters of half-sib families for slash pine. Scientia Silvae Sinicae, 36(1): 56-61. [in Chinese])

    徐有明,鮑春紅,周志翔,等. 2001. 濕地松種源生長量、材性的變異與優(yōu)良種源綜合選擇. 東北林業(yè)大學(xué)學(xué)報,29(5): 18-21.

    (Xu Y M, Bao C H, Zhou Z X,etal. 2001. Variation analyses of tree growth and wood properties among slash pine provenances and comprehensive selection of superior provenance. Journal of Northeast Forestry University, 29(5): 18-21. [in Chinese])

    易 敏,賴 猛,孫曉梅,等. 2014. 無損檢測在日本落葉松材性育種中的應(yīng)用前景探討. 林業(yè)科學(xué),50(11): 96-103.

    (Yi M, Lai M, Sun X M,etal. 2014. Application prospect on applying nondestructive testing technology in wood property breeding ofLarixkaempferi. Scientia Silvae Sinicae, 50(11): 96-103. [in Chinese])

    朱景樂,王軍輝,張守攻,等. 2008. 毛白楊材性指標(biāo)預(yù)測及選擇. 林業(yè)科學(xué),44(7): 23-28.

    (Zhu J L, Wang J H, Zhang S G,etal. 2008. Wood property estimation and selection ofPopulustomentosa. Scientia Silvae Sinicae, 44(7): 23-28. [in Chinese])

    朱志淞,丁衍疇. 1993. 濕地松. 廣州: 廣東科技出版社,62-298.

    (Zhu Z S, Ding Y C. 1993. Slash pine. Guangzhou: Guangdong Science & Technology Press, 62-298.[in Chinese])

    Butler D G, Cullis B R, Gilmour A R,etal. 2009. ASReml-R reference manual. Queensland Department of Primary Industries.

    Carter P, Briggs D, Ross R J,etal. 2005. Acoustic testing to enhance Western forest values and meet customer wood quality needs//Harrington C A, Schoenholtz S H. Productivity of Western forests: a forest products focus. General technical report PNW-GTR-642. Portland: U.S. Department of Agriculture, Forest Service, Pacific Northwest Research,121-129.

    Chen Z Q, Karlsson B, Lundqvist S O,etal. 2015. Estimating solid wood properties using Pilodyn and acoustic velocity on standing trees of Norway spruce. Annals of Forest Science, 72(4):1-10.

    Fibre-gen. 2004. Director ST300-operating manual. http://www.fibre-gen.com/documents/ST300OperatingManual.pdf.

    Fries A. 1986. Volume growth and wood density of plus tree progenies ofPinuscontortain two Swedish field trials. Scand J For Res,1(1/4): 403-419.

    Kumar S. 2004. Genetic parameter estimates for wood stiffness, strength, internal checking, and resin bleeding for radiata pine. Canadian Journal of Forest Research, 34(12): 2601-2610.

    Lasserre J P, Mason E G, Watt M S. 2007. Assessing corewood acoustic velocity and modulus of elasticity with two impact based instruments in 11-year-old trees from a clonal-spacing experiment ofPinusradiataD. Don. Forest Ecology and Management, 239: 217-221.

    Lenz P, Auty D, Achim A,etal. 2013. Genetic improvement of white spruce mechanical wood traits—early screening by means of acoustic velocity. Forests, 4(3): 575-594.

    Namkoong G. 1981. Introduction to quantitative genetics in forestry. London: Castle House Publications.

    Raymond C A, Kube P D, Pinkard L,etal. 2004. Evaluation of non-destructive methods of measuring growth stress inEucalyptusglobulus: relationships between strain, wood properties and stress. Forest Ecology and Management, 190(2): 187-200.

    Ross R J, Pellerin R F. 1988. NDE of wood-based composites with longitudinal stress wave. Forest Products Journal, 38(5): 39-45.

    Ukrainetz N K, Kang K Y, Aitken S N,etal. 2008. Heritability and phenotypic and genetic correlations of coastal Douglas-fir (Pseudotsugamenziesii) wood quality traits. Canadian Journal of Forest Research, 38(6): 1536-1546.

    Wang T, Aitken S N, Rozenberg P,etal. 1999. Selection for height growth and Pilodyn pin penetration in lodgepole pine: effects on growth traits, wood properties and their relationships. Can J For Res, 29: 434-445.

    White T L, Adams W T, Neale D B.2007. Forest genetics. CABI, Wallingford.

    White T L, Byram T D. 2004. Slash pine tree improvement//Dickens E D, Barnett J P, Hubbard W G,etal. Proceedings of the slash pine symposium. Slash pine: still growing and growing! SRS-76. Asheville, NC: U.S. Department of Agriculture, Forest Service, Southern Research Station, 7-19.

    Zobel B J, Kellison R C, Matthias M F,etal. 1972. Wood density of the southern pines. North Carolina Agricultural Experiment Station.Tech Bul 208:56.

    Zobel B J, Van Buijtenen J P. 1989. Wood variation: its causes and control. Springer Science & Business Media.

    (責(zé)任編輯 徐 紅)

    Genetic Variation Analysis for Growth and Wood Properties of Slash Pine Based onThe Non-Destructive Testing Technologies

    Zhang Shuainan Luan Qifu Jiang Jingmin

    (ZhejiangProvincialKeyLaboratoryofTreeBreedingResearchInstituteofSubtropicalForestry,ChineseAcademyofForestryHangzhou311400)

    【Objective】 The genetic variability and genetic and phenotypic correlations among Basic Wood Density (ρ), Modulus of Elasticity (MOE) and growth were studied to strengthen the selection and breeding intensity of the physical and mechanical wood properties in order to provide the better germplasms for growing slash pine(Pinuselliottii)plantation.【Method】 Theρa(bǔ)nd MOE were determined by the non-destructive evaluation technique, Pilodyn and acoustic velocity (AV), studied on the half-sib families of slash pine at 22 years old in Changle State Forest Farm in northern Zhejiang province. The growth traits such as diameter at breast height (DBH), height (H), under-branch height (Hb), crown width (CW), branch size (B) and branch angle (Ba) were measured by conventional method. The heritabilities and genetic- and phenotypic-correlations for growth and wood properties were estimated using Residual Maximum Likelihood (REML) in the flexible mixed modeling program ASReml-R.【Result】 1) The results of nondestructive testing of wood properties showed that relative value of MOE (MOEP)had highly positive and significant genetic (R=0.947 0±0.020 1,P<0.001) and phenotypic (R=0.948 0±0.002 7,P<0.001) correlations with wave velocity (v), and negative genetic (R=-0.447 0±0.154 0,P<0.01) and phenotypic (R=-0.538 0±0.019 1,P<0.001) correlations with relative value ofρ(ρP) (It indicated that MOE andρwere closely related. But the correlation coefficient was much less than the correlation withv). 2) The heritabilities of wood properties was between 0.292 0~0.305 0, and the heritabilities of growth traits was between 0.062 5~0.216 0. The genetic gains of volume (V) and MOEPcould be 30% (the selection rate=1%). 3) For growth traits, the genetic correlations among DBH,V, CW andBwere significantly positive. And the genetic correlations betweenHand DBH,Vwere highly positive significantly.Hbjust had moderate positive genetic correlation withH.Bashowed weak phenotypic correlation withBnegatively only.Other traits showed no significant correlations. It is worth mentioning that the phenotypic correlations among DBH,H,Hb,VandBwere extremely significantly positive. The genetic correlations between growth and wood properties were not significant. However, the phenotypic correlations between them were significant but the correlation coefficients were low.【Conclusion】 1) Only usevcould quickly assess the relative value of MOE of each tree within the population. The selection effect is better than that only usingρto evaluate MOE. 2)ρa(bǔ)nd MOE were controlled by a moderate degree of genetic control, and the heritabilities were higher than growth traits obviously. 3) Wood properties just had significant phenotypic correlations with growth traits. Improvement of these traits simultaneously was possible.

    Pinuselliottii; growth traits; wood density; modulus of elasticity; non-destructive evaluation; genetic analysis

    10.11707/j.1001-7488.20170604

    2016-06-16;

    2016-07-15。

    國家自然科學(xué)基金項目(31570668); “十二五”國家林業(yè)科技支撐計劃專題(2012BAD01B0203)。

    S722.3+3

    A

    1001-7488(2017)06-0030-07

    *欒啟福為通訊作者。

    猜你喜歡
    濕地松遺傳變異遺傳力
    某國家生豬核心育種場加系大白母豬繁殖性狀遺傳參數(shù)估計
    美系大白種豬生長性狀和繁殖性狀遺傳力估計
    養(yǎng)豬(2022年4期)2022-08-17 07:07:02
    濕地松豐產(chǎn)造林技術(shù)探討
    濕地松和火炬松栽培技術(shù)差異
    先導(dǎo)編輯技術(shù)可編輯近90%的人類遺傳變異
    基于改進(jìn)遺傳變異算子的海島算法
    電子制作(2019年24期)2019-02-23 13:22:18
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    濕地松人工林生長規(guī)律
    通過胎盤效率改良母豬繁殖性能的研究進(jìn)展
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    毛片女人毛片| 国产午夜精品论理片| 亚洲在久久综合| 免费播放大片免费观看视频在线观看| 亚洲国产精品成人综合色| 日韩电影二区| 亚洲欧美成人综合另类久久久| 亚洲一区高清亚洲精品| 1000部很黄的大片| 国产高清国产精品国产三级 | 国产av码专区亚洲av| 亚洲美女视频黄频| 国产69精品久久久久777片| 欧美xxⅹ黑人| 国产淫片久久久久久久久| av在线蜜桃| 国产人妻一区二区三区在| freevideosex欧美| 国产一级毛片在线| 国产真实伦视频高清在线观看| 97热精品久久久久久| 久久久久久久久久成人| 国产av国产精品国产| 亚洲无线观看免费| 亚洲国产精品专区欧美| 久久人人爽人人片av| 欧美激情国产日韩精品一区| 听说在线观看完整版免费高清| 伊人久久国产一区二区| 最近最新中文字幕免费大全7| 国产亚洲一区二区精品| 久久久久久九九精品二区国产| 欧美成人一区二区免费高清观看| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 免费看光身美女| 男女边吃奶边做爰视频| 午夜免费观看性视频| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕| 能在线免费看毛片的网站| 国产在线男女| 亚洲精品视频女| 成年av动漫网址| 寂寞人妻少妇视频99o| 一区二区三区免费毛片| 伊人久久精品亚洲午夜| 国产黄a三级三级三级人| 日本黄色片子视频| 精品人妻视频免费看| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 一级毛片我不卡| 免费av观看视频| 我要看日韩黄色一级片| 日韩欧美国产在线观看| 亚洲av成人av| 国产91av在线免费观看| 国产亚洲最大av| 麻豆av噜噜一区二区三区| 成人毛片60女人毛片免费| 亚洲精品久久午夜乱码| 麻豆av噜噜一区二区三区| 亚洲国产日韩欧美精品在线观看| 99久久精品一区二区三区| 91精品一卡2卡3卡4卡| 内射极品少妇av片p| 国产av码专区亚洲av| 国产免费一级a男人的天堂| 汤姆久久久久久久影院中文字幕 | 亚洲国产精品成人久久小说| 又黄又爽又刺激的免费视频.| 91狼人影院| 波野结衣二区三区在线| 美女主播在线视频| 亚洲熟妇中文字幕五十中出| 免费播放大片免费观看视频在线观看| 一级a做视频免费观看| 少妇裸体淫交视频免费看高清| 午夜福利视频精品| 97精品久久久久久久久久精品| 亚洲精品国产av蜜桃| 不卡视频在线观看欧美| 亚洲国产高清在线一区二区三| 少妇人妻一区二区三区视频| a级毛片免费高清观看在线播放| 黄片无遮挡物在线观看| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 少妇丰满av| 毛片一级片免费看久久久久| 精品一区二区免费观看| 男的添女的下面高潮视频| 国产探花在线观看一区二区| 麻豆av噜噜一区二区三区| 精品欧美国产一区二区三| 91久久精品国产一区二区三区| 国产亚洲5aaaaa淫片| 亚洲乱码一区二区免费版| 国产午夜精品久久久久久一区二区三区| av天堂中文字幕网| 秋霞伦理黄片| 简卡轻食公司| 内地一区二区视频在线| 亚洲欧美日韩东京热| 人体艺术视频欧美日本| 久久精品国产自在天天线| 久久久久久国产a免费观看| 午夜福利在线观看吧| 成人二区视频| 亚洲在线观看片| 亚洲经典国产精华液单| 国内揄拍国产精品人妻在线| 免费观看无遮挡的男女| 亚洲国产成人一精品久久久| 中文资源天堂在线| 国产极品天堂在线| 亚洲精品国产成人久久av| 99久久精品国产国产毛片| 熟妇人妻不卡中文字幕| 亚洲最大成人中文| 国产真实伦视频高清在线观看| 色哟哟·www| 国语对白做爰xxxⅹ性视频网站| 国产v大片淫在线免费观看| 国产成人午夜福利电影在线观看| 亚洲精品亚洲一区二区| 久久久久久久午夜电影| 国产黄色视频一区二区在线观看| 日韩电影二区| 日韩av不卡免费在线播放| 99久久精品热视频| 亚洲高清免费不卡视频| av黄色大香蕉| 成人亚洲欧美一区二区av| 国产 一区精品| 男女边吃奶边做爰视频| 久99久视频精品免费| 亚洲人成网站高清观看| 亚洲第一区二区三区不卡| 在线观看av片永久免费下载| 97在线视频观看| videossex国产| 99久国产av精品| 欧美激情国产日韩精品一区| 国产乱人视频| 精品一区在线观看国产| 国产精品国产三级国产专区5o| 成年女人在线观看亚洲视频 | 国产又色又爽无遮挡免| 2022亚洲国产成人精品| 欧美性猛交╳xxx乱大交人| 男人爽女人下面视频在线观看| 天堂影院成人在线观看| 最近中文字幕2019免费版| 美女主播在线视频| 18禁在线无遮挡免费观看视频| 欧美 日韩 精品 国产| 日韩在线高清观看一区二区三区| 久久99精品国语久久久| 深爱激情五月婷婷| 成人二区视频| 亚洲av二区三区四区| av国产免费在线观看| 乱人视频在线观看| 麻豆乱淫一区二区| 色5月婷婷丁香| 男的添女的下面高潮视频| 国产欧美另类精品又又久久亚洲欧美| 国内精品宾馆在线| 成年免费大片在线观看| 日韩成人伦理影院| av国产免费在线观看| 久久久久久久午夜电影| 久久精品夜色国产| 国产成人一区二区在线| 少妇人妻一区二区三区视频| 亚洲精品乱久久久久久| 国产精品美女特级片免费视频播放器| 人人妻人人看人人澡| ponron亚洲| 亚洲av成人av| 午夜爱爱视频在线播放| 国产亚洲av片在线观看秒播厂 | 91aial.com中文字幕在线观看| 一级a做视频免费观看| 精品人妻熟女av久视频| 秋霞伦理黄片| 亚洲欧美日韩无卡精品| 一个人观看的视频www高清免费观看| 欧美成人精品欧美一级黄| 十八禁国产超污无遮挡网站| 97热精品久久久久久| 国产黄片美女视频| 成年女人在线观看亚洲视频 | 国产成人精品婷婷| 内射极品少妇av片p| 中文字幕久久专区| 国产v大片淫在线免费观看| 国模一区二区三区四区视频| 亚洲国产欧美人成| 国产欧美另类精品又又久久亚洲欧美| 久久久久免费精品人妻一区二区| 亚洲在线观看片| videos熟女内射| 一级爰片在线观看| 免费人成在线观看视频色| av在线亚洲专区| 黄色日韩在线| av在线蜜桃| 久久99精品国语久久久| 久久久久久伊人网av| 日韩欧美精品免费久久| 精品久久久久久成人av| 三级男女做爰猛烈吃奶摸视频| 最新中文字幕久久久久| 一区二区三区乱码不卡18| 天堂影院成人在线观看| 免费av不卡在线播放| 亚洲国产精品专区欧美| 午夜激情欧美在线| www.av在线官网国产| 国内揄拍国产精品人妻在线| 2021少妇久久久久久久久久久| 成人欧美大片| 国产真实伦视频高清在线观看| 久久精品久久精品一区二区三区| 国产麻豆成人av免费视频| 99久国产av精品| 一二三四中文在线观看免费高清| 色哟哟·www| 男人和女人高潮做爰伦理| 久久人人爽人人片av| 观看免费一级毛片| 日韩视频在线欧美| 最近中文字幕2019免费版| 久久精品国产自在天天线| 亚洲内射少妇av| 欧美一区二区亚洲| 99热这里只有是精品在线观看| 国产在线男女| 九九在线视频观看精品| 婷婷色综合大香蕉| 国产精品一及| 国产视频首页在线观看| 啦啦啦中文免费视频观看日本| 亚洲国产精品成人综合色| 日本一二三区视频观看| 久久99热6这里只有精品| 精品99又大又爽又粗少妇毛片| 国产伦精品一区二区三区四那| 久久鲁丝午夜福利片| 三级国产精品欧美在线观看| 欧美bdsm另类| 99热网站在线观看| 国产精品女同一区二区软件| 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 欧美潮喷喷水| 美女脱内裤让男人舔精品视频| av国产免费在线观看| 午夜免费激情av| 国产精品麻豆人妻色哟哟久久 | 日本爱情动作片www.在线观看| 特大巨黑吊av在线直播| .国产精品久久| 九色成人免费人妻av| 青春草国产在线视频| 白带黄色成豆腐渣| 色吧在线观看| 人妻系列 视频| 精品酒店卫生间| 一级毛片 在线播放| 男人爽女人下面视频在线观看| 国产精品不卡视频一区二区| 亚洲怡红院男人天堂| 国内精品宾馆在线| 国产探花极品一区二区| 亚州av有码| 99热6这里只有精品| 天天一区二区日本电影三级| 亚洲美女视频黄频| 99热6这里只有精品| 国产精品一区二区性色av| 一个人观看的视频www高清免费观看| 国产伦理片在线播放av一区| 舔av片在线| 精品久久久噜噜| 成人国产麻豆网| 国产精品国产三级国产专区5o| 麻豆久久精品国产亚洲av| 丰满少妇做爰视频| 成人亚洲欧美一区二区av| 91在线精品国自产拍蜜月| 国产又色又爽无遮挡免| 美女国产视频在线观看| 亚洲成人久久爱视频| 欧美潮喷喷水| 亚洲av日韩在线播放| 亚洲av.av天堂| 欧美日韩视频高清一区二区三区二| 韩国高清视频一区二区三区| 97超碰精品成人国产| 久久久午夜欧美精品| 97在线视频观看| 三级国产精品片| 国产精品人妻久久久久久| 亚洲激情五月婷婷啪啪| 日韩不卡一区二区三区视频在线| 一区二区三区四区激情视频| 春色校园在线视频观看| 久99久视频精品免费| 亚洲国产精品国产精品| 欧美97在线视频| 黄色一级大片看看| 国产精品无大码| 中文字幕av在线有码专区| 色网站视频免费| 国产午夜精品久久久久久一区二区三区| 精品久久久久久电影网| 亚洲精品国产成人久久av| 男插女下体视频免费在线播放| 国产一级毛片在线| 不卡视频在线观看欧美| 少妇人妻一区二区三区视频| 成人漫画全彩无遮挡| 麻豆乱淫一区二区| 伊人久久精品亚洲午夜| 黄片wwwwww| 欧美xxxx黑人xx丫x性爽| av线在线观看网站| 成人高潮视频无遮挡免费网站| 男女边摸边吃奶| 精品久久国产蜜桃| 亚洲aⅴ乱码一区二区在线播放| 国产人妻一区二区三区在| 国产色婷婷99| 免费av不卡在线播放| 久久久久精品久久久久真实原创| 亚洲欧美日韩卡通动漫| 亚洲精品456在线播放app| 日日摸夜夜添夜夜添av毛片| 午夜福利在线在线| av一本久久久久| 99热全是精品| kizo精华| 自拍偷自拍亚洲精品老妇| 高清av免费在线| 非洲黑人性xxxx精品又粗又长| 久久久久久久久中文| 国内少妇人妻偷人精品xxx网站| 一级a做视频免费观看| 人妻少妇偷人精品九色| av专区在线播放| 日韩精品青青久久久久久| 一级a做视频免费观看| 99热6这里只有精品| 欧美人与善性xxx| 免费高清在线观看视频在线观看| 99久久中文字幕三级久久日本| av国产久精品久网站免费入址| 亚洲美女搞黄在线观看| 亚洲av二区三区四区| 亚洲怡红院男人天堂| 亚洲国产最新在线播放| 极品教师在线视频| 三级男女做爰猛烈吃奶摸视频| 久久精品国产自在天天线| 一二三四中文在线观看免费高清| 亚洲国产成人一精品久久久| 免费av观看视频| 久久久亚洲精品成人影院| 久久人人爽人人爽人人片va| 午夜福利视频1000在线观看| 日韩一区二区视频免费看| 欧美xxⅹ黑人| 精品一区在线观看国产| 韩国高清视频一区二区三区| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄| 能在线免费观看的黄片| 小蜜桃在线观看免费完整版高清| 亚洲内射少妇av| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 久久久久网色| 九九在线视频观看精品| 18禁在线播放成人免费| 一级毛片aaaaaa免费看小| 婷婷六月久久综合丁香| 极品教师在线视频| 日日啪夜夜撸| 国产麻豆成人av免费视频| 日本免费在线观看一区| 一二三四中文在线观看免费高清| 日韩在线高清观看一区二区三区| 亚洲精品色激情综合| 欧美 日韩 精品 国产| 中文字幕av成人在线电影| 99久久精品国产国产毛片| 精品酒店卫生间| 亚洲av电影在线观看一区二区三区 | 黄色配什么色好看| 成人漫画全彩无遮挡| 日韩视频在线欧美| 精品一区二区免费观看| 又爽又黄无遮挡网站| 日韩av在线免费看完整版不卡| videossex国产| 一个人看视频在线观看www免费| 男女下面进入的视频免费午夜| 韩国av在线不卡| 国产一级毛片七仙女欲春2| 日日摸夜夜添夜夜爱| 日日摸夜夜添夜夜添av毛片| 嫩草影院精品99| 国产一区有黄有色的免费视频 | 精华霜和精华液先用哪个| 最近中文字幕高清免费大全6| 国产在视频线在精品| 人体艺术视频欧美日本| 亚洲精品国产成人久久av| 国产精品国产三级专区第一集| 欧美一区二区亚洲| 国产伦一二天堂av在线观看| 毛片一级片免费看久久久久| 日本黄大片高清| 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 禁无遮挡网站| 女人久久www免费人成看片| 亚洲av电影在线观看一区二区三区 | 少妇裸体淫交视频免费看高清| 日韩欧美国产在线观看| 国产精品福利在线免费观看| 午夜福利视频1000在线观看| 大香蕉久久网| 日本猛色少妇xxxxx猛交久久| 中文字幕免费在线视频6| 中文字幕av成人在线电影| 一区二区三区乱码不卡18| 欧美人与善性xxx| 淫秽高清视频在线观看| 中文在线观看免费www的网站| 亚洲欧美日韩无卡精品| 亚洲成人久久爱视频| 午夜福利在线观看免费完整高清在| 中国国产av一级| 中文字幕久久专区| 精品一区二区免费观看| 亚洲国产成人一精品久久久| 国产单亲对白刺激| 少妇熟女aⅴ在线视频| av在线天堂中文字幕| 国产在线男女| 国产欧美另类精品又又久久亚洲欧美| 最近2019中文字幕mv第一页| 国产精品一区二区三区四区免费观看| 亚洲最大成人av| 亚洲久久久久久中文字幕| 欧美精品一区二区大全| 久久久久九九精品影院| 岛国毛片在线播放| 亚洲四区av| 神马国产精品三级电影在线观看| 中文字幕制服av| 国产 亚洲一区二区三区 | 又爽又黄a免费视频| 久久人人爽人人爽人人片va| 国产伦一二天堂av在线观看| 一级毛片我不卡| 国内少妇人妻偷人精品xxx网站| 干丝袜人妻中文字幕| 草草在线视频免费看| 中文字幕久久专区| 日本免费a在线| 色网站视频免费| 国产精品一区www在线观看| 亚洲欧美精品专区久久| 天堂影院成人在线观看| 国产亚洲最大av| 九草在线视频观看| 日韩亚洲欧美综合| 国产三级在线视频| 亚洲国产精品国产精品| 一级毛片电影观看| 国产成人精品福利久久| 精品一区二区三区视频在线| 亚洲精品日韩av片在线观看| 国产 一区 欧美 日韩| 淫秽高清视频在线观看| 人人妻人人澡人人爽人人夜夜 | 老司机影院毛片| 国产人妻一区二区三区在| 国产 一区精品| 午夜久久久久精精品| 激情五月婷婷亚洲| 国产黄片美女视频| 久99久视频精品免费| 五月伊人婷婷丁香| 国产免费又黄又爽又色| 搡老妇女老女人老熟妇| h日本视频在线播放| 国产精品久久久久久av不卡| 成年版毛片免费区| 欧美成人a在线观看| 国产高清不卡午夜福利| 国产亚洲av嫩草精品影院| 亚洲精品自拍成人| 亚洲精品乱码久久久v下载方式| 少妇的逼好多水| 国产一区二区三区综合在线观看 | 一夜夜www| 免费观看在线日韩| 又爽又黄无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品专区欧美| 久久久欧美国产精品| 午夜福利视频精品| 精品国产一区二区三区久久久樱花 | 免费观看a级毛片全部| videos熟女内射| 美女国产视频在线观看| 国产探花在线观看一区二区| 亚洲精品久久午夜乱码| 亚洲在久久综合| 欧美日韩国产mv在线观看视频 | 国产成人精品一,二区| 久久精品综合一区二区三区| 日韩国内少妇激情av| 丰满少妇做爰视频| 日韩av不卡免费在线播放| 国产单亲对白刺激| 七月丁香在线播放| 亚洲精品一二三| 校园人妻丝袜中文字幕| 欧美+日韩+精品| 青青草视频在线视频观看| 亚洲久久久久久中文字幕| 国产伦在线观看视频一区| 国产精品精品国产色婷婷| 成人高潮视频无遮挡免费网站| 91精品伊人久久大香线蕉| av专区在线播放| 国产av不卡久久| 国产免费又黄又爽又色| 国产成人福利小说| 免费看不卡的av| 乱人视频在线观看| 午夜视频国产福利| 日本猛色少妇xxxxx猛交久久| 国产伦精品一区二区三区四那| 国产男人的电影天堂91| 天堂网av新在线| 一边亲一边摸免费视频| 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| av在线亚洲专区| 激情 狠狠 欧美| 国产精品99久久久久久久久| 国产高清国产精品国产三级 | 成人毛片a级毛片在线播放| av国产久精品久网站免费入址| 欧美激情久久久久久爽电影| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 色播亚洲综合网| 欧美日韩一区二区视频在线观看视频在线 | 国产三级在线视频| 精品99又大又爽又粗少妇毛片| 婷婷六月久久综合丁香| 亚洲精品一区蜜桃| 亚洲一级一片aⅴ在线观看| 成人欧美大片| 超碰av人人做人人爽久久| www.av在线官网国产| 97热精品久久久久久| 午夜福利网站1000一区二区三区| 亚洲精品第二区| 欧美日韩在线观看h| 亚洲久久久久久中文字幕| 久久久久久久久久黄片| 中文字幕av在线有码专区| 国产黄片视频在线免费观看| 夫妻性生交免费视频一级片| 久久国内精品自在自线图片| 亚洲国产精品sss在线观看| 亚洲丝袜综合中文字幕| 成人无遮挡网站| 亚洲欧美精品专区久久| 国产淫语在线视频| 六月丁香七月| 国产黄频视频在线观看| 天堂网av新在线| 欧美成人a在线观看| 亚洲乱码一区二区免费版| 免费人成在线观看视频色| av国产久精品久网站免费入址| 欧美日韩在线观看h| av专区在线播放| 18禁裸乳无遮挡免费网站照片| 亚洲成人中文字幕在线播放| 精品一区二区三卡| 18禁在线播放成人免费| 久久久久久久久久人人人人人人| 国产精品嫩草影院av在线观看| 久久99精品国语久久久| 极品教师在线视频| 国产成人免费观看mmmm| 亚洲最大成人中文| 国产中年淑女户外野战色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 九草在线视频观看| 成人性生交大片免费视频hd|