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

    中亞熱帶赤皮青岡天然種群表型多樣性分析

    2023-11-01 04:51:20秦之曠劉娜周霞李志輝曹福亮李何
    廣西植物 2023年9期
    關(guān)鍵詞:表型性狀種質(zhì)資源

    秦之曠 劉娜 周霞 李志輝 曹福亮 李何

    摘?要:赤皮青岡(Quercus gilva)是我國中亞熱帶地區(qū)極具經(jīng)濟價值的珍貴用材樹種。為揭示其表型多樣性與變異規(guī)律以及影響表型變異的關(guān)鍵地理氣候因子,該研究以14個赤皮青岡天然種群115個單株為材料,對其15個樹體和葉片性狀進行測量,并采用巢式方差分析、表型分化系數(shù)分析、多樣性指數(shù)分析、相關(guān)性分析、主成分分析以及聚類分析等方法,探究赤皮青岡的表型變異規(guī)律及其與地理氣候因子的相關(guān)性。結(jié)果表明:(1)赤皮青岡表型多樣性水平較高,15個表型性狀的變異系數(shù)和Shannon-Wiener指數(shù)平均值分別為35.070%和1.998,14個種群的變異系數(shù)在14.94%(洞口)~35.56%(龍山)之間、Shannon-Wiener指數(shù)在1.127(松陽)~1.980(常寧)之間。(2)15個表型性狀在種群間和種群內(nèi)均存在極顯著差異(P<0.01),種群平均表型分化系數(shù)為41.491%,表型變異主要來源于種群內(nèi)。(3)相關(guān)性分析發(fā)現(xiàn),赤皮青岡部分性狀間存在顯著或極顯著相關(guān)性,表型性狀與地理因子間相關(guān)性不明顯,降水是影響赤皮青岡表型性狀的主要氣候因子。(4)主成分分析發(fā)現(xiàn),基于11個表型可提煉出4個主成分,累計貢獻率為82.961%。(5)根據(jù)聚類分析可將14個種群分為三類,分別為小葉片小葉柄類、大葉片大葉柄類和中葉片小葉柄類,種群的聚類結(jié)果與其地理位置間無明顯相關(guān)性。該研究結(jié)果為赤皮青岡種質(zhì)資源的保護和利用提供了科學依據(jù),為其良種選育奠定了重要基礎(chǔ)。

    關(guān)鍵詞: 赤皮青岡, 珍貴樹種, 種質(zhì)資源, 表型性狀, 變異規(guī)律

    中圖分類號:Q948.12

    文獻標識碼:A

    文章編號:1000-3142(2023)09-1622-14

    收稿日期:2022-09-15

    基金項目:國家自然科學基金 (32201589); 湖南省自然科學基金 (2021JJ41069); 湖南省林業(yè)科技創(chuàng)新資金項目 (XLKY202218); 福建省公益類科研院所專項(2022R1010007)。

    第一作者: 秦之曠(1997-),碩士,研究方向為林木育種,(E-mail)1667830341@qq.com。

    *通信作者:李何,博士,講師,碩士研究生導師,研究方向為林木育種,(E-mail)lihecsuft@126.com。

    Phenotypic diversity of Quercus gilva natural

    populations in middle subtropical China

    QIN Zhikuang1, LIU Na1, ZHOU Xia1, LI Zhihui1, CAO Fuliang2, LI He1*

    ( 1. College of Forestry, Central South University of Forestry and Technology, Changsha 410004, China;

    2. College of Forestry, Nanjing Forestry University, Nanjing 210037, China )

    Abstract:Quercus gilva is a precious timber tree species in middle subtropical China. In order to investigate the phenotypic diversity, phenotypic variation pattern and the major geographic and climatic factors influencing phenotypic variation of Q. gilva populations, 15 growth and leaf traits of 115 individuals from 14 natural populations were measured, and then the phenotypic diversity, variation pattern and the correlation between phenotypic variation and geographic and climatic factors were studied using nested analysis of variance, phenotypic differentiation coefficient, diversity index, correlation analysis, principal component analysis and cluster analysis. The results were as follows: (1) The average variation coefficient and the average Shannon-Wiener index of 15 phenotypic traits were 35.070% and 1.998, respectively. Variation coefficient and Shannon-Wiener index of 14 populations ranged from 14.94% (Dongkou) to 35.56% (Longshan) and from 1.127 (Songyang) to 1.980 (Changning), respectively. These results indicated that Q. gilva showed a relatively high level of phenotypic diversity. (2) Significant differences in 15 phenotypic traits were found among and within populations (P<0.01), with an average phenotypic differentiation coefficient of 41.491%. This result elucidated that the phenotypic variation within populations was greater than that among populations. (3) Correlation analysis showed that there were significant or entremely significant correlations among parts of the traits, whereas no correlation was found between phenotypic traits and geographic factors. Precipitation was the major factor that influenced the phenotypic traits of Q. gilva. (4) Principle component analysis indicated that the first four principal components explained 82.961% of total phenotypic contribution rate of Q. gilva. (5) Cluster analysis revealed that 14 Q. gilva populations were divided into three groups according to their different leaf and petiole size, such as small-leaved micropetiole, large-leaved macropetiole and medium-leaved micropetiole, which was not related to the geographic distribution of populations. This study provides scientific basis for the protection and utilization of Q. gilva germplasm resources, and lays an important foundation for its improved variety breeding.

    Key words: Quercus gilva, precious tree species, germplasm resources, phenotypic traits, variation pattern

    赤皮青岡(Quercus gilva)是殼斗科青岡屬常綠喬木,其樹干通直,邊材呈黃褐色,心材呈暗紅褐色,紋理優(yōu)美具有光澤,強韌有彈性,有獨特香味,抗蟲蛀,不易腐蝕,是我國亞熱帶地區(qū)極具經(jīng)濟價值和發(fā)展前景的珍貴用材樹種之一(端木,1995;Zaynab et al., 2018)。其樹體高大,為深根性樹種,在水土保持和涵養(yǎng)水源方面有良好的表現(xiàn),從而具有很高的生態(tài)價值(夏蟬,2021)。赤皮青岡在我國中亞熱帶地區(qū)分布較廣,主要分布于湖南、福建、浙江、江西、貴州、臺灣等?。╖eng et al., 2019;歐陽澤怡等,2021)。然而,從20世紀開始,赤皮青岡天然群落受到明顯的人為干擾,絕大部分赤皮青岡天然林被過度采伐,采伐跡地用于城鄉(xiāng)發(fā)展、筑路修橋或人工林營造,導致大多赤皮青岡處于散生和小群落生存狀態(tài),其種質(zhì)資源大量流失(Deng et al., 2018)。

    種質(zhì)資源既是遺傳多樣性的載體,也是種質(zhì)創(chuàng)新和林木遺傳改良的物質(zhì)基礎(chǔ),闡明種質(zhì)資源的多樣性是對其合理保護和利用的前提(Bai et al., 2019)。植物表型是遺傳因子和環(huán)境在長期過程中共同作用的結(jié)果,表型多樣性是遺傳多樣性的直接體現(xiàn)(李洪果等,2019;郭燕等,2022)。利用植物的表型性狀檢測遺傳多樣性具有簡便、快捷的特點(尚帥斌等,2015),是長期以來最基本和最直接的方法。近年來,種質(zhì)資源的表型性狀分析已在閩楠(Phoebe bournei)(Li et al., 2018;歐漢彪等,2021)、黑楊(Populus nigra)(陳存等,2021)、板栗(Castanea mollissima)(郭燕等,2022)、茶條槭(Acer ginnala)(武艷虹等,2018)、唐古特白刺(Nitraria tanguticum)(吉德娟等,2021)等林木上開展,通過分析各表型性狀,揭示了表型多樣性和變異,挖掘了影響表型變異的關(guān)鍵因子,為相關(guān)樹種種質(zhì)資源的保護和利用奠定了重要基礎(chǔ)。

    目前,赤皮青岡天然種群仍處于持續(xù)減少的狀態(tài),天然種質(zhì)資源不斷流失,調(diào)查與闡明天然種質(zhì)的多樣性以合理制定其保護和利用策略是當前亟待開展的重要研究。但是,關(guān)于赤皮青岡種質(zhì)資源的相關(guān)研究非常有限,僅謝?。?011)和林文俊等(2017)對赤皮青岡天然群落進行了調(diào)查,并且調(diào)查群落僅分布在閩清和擎天巖村,很難具有代表性。在分子層面,雖然有基于ISSR分子標記技術(shù)(朱品紅,2014)和葉綠體基因組測序技術(shù)(Sugiura et al., 2015)的赤皮青岡遺傳多樣性和遺傳結(jié)構(gòu)的相關(guān)報道,但所調(diào)查的種群分布較集中,不足以揭示我國赤皮青岡種質(zhì)的遺傳多樣性;在表型層面,僅謝健(2021)對4個種源赤皮青岡幼林的生長性狀進行了分析。樹高、胸徑、材積和干型等是用材樹種的重要經(jīng)濟性狀,而葉片是植物進行光合和蒸騰作用的重要媒介,其形態(tài)直接影響生長發(fā)育和繁殖特性,并最終影響樹種的適應(yīng)和進化(郭燕等,2022)。

    本研究以具有代表性的、基本覆蓋赤皮青岡在大陸地區(qū)分布范圍的14個赤皮青岡天然種群115個單株為材料,對其15個樹體和葉片性狀進行測量,采用巢式方差分析、表型分化系數(shù)分析、多樣性指數(shù)分析、相關(guān)性分析、主成分分析以及聚類分析等方法,通過分析表型多樣性和變異及其地理變化,擬探討:(1)赤皮青岡天然種群表型多樣性水平;(2)赤皮青岡表型多樣性呈現(xiàn)的變異規(guī)律;(3)赤皮青岡表型變異與地理氣候因子之間的相關(guān)性。以期為赤皮青岡種質(zhì)資源的收集、保護和利用提供科學依據(jù),為其良種選育奠定重要基礎(chǔ)。

    1?材料與方法

    1.1 材料

    基于前期對赤皮青岡天然分布區(qū)內(nèi)種質(zhì)資源的調(diào)查,2021年5月至2022年3月對14個赤皮青岡天然種群共115個單株進行了調(diào)查和樣本采集,采樣種群基本覆蓋了我國中亞熱帶赤皮青岡天然分布區(qū),各種群的地理氣候信息和樣本數(shù)量如表1所示。對樣本數(shù)大于10個的種群,隨機選取10株樣本,并確保樣本間較大的距離;由于部分種群樣本較少或采樣困難,因此對樣本數(shù)小于10株的種群,則采集種群內(nèi)的所有樣本。用GPS記錄各樣本的經(jīng)度和緯度,利用測高儀、胸徑尺等工具測量且記錄樹體性狀。用高枝剪等工具在當?shù)?2:00—14:00時樹冠中下層的陽面,獲取健康枝條,隨機采集成熟葉片(李珊等,2016),并選取10片無病蟲害的完整葉片進行集中拍照、記錄。

    1.2 氣候數(shù)據(jù)的獲取

    從全球氣候數(shù)據(jù)網(wǎng)WorldClim(https://www.worldclim.org/)下載當前時期氣候數(shù)據(jù),共下載與溫度和降水相關(guān)的19個氣候因子,數(shù)據(jù)精度為2.5’,版本為2.1(Fick & Hijmans, 2017),使用R包Raster中extract函數(shù)提取種群所在位置的氣候數(shù)據(jù)。在19個氣候因子中,由于一些因子間可能存在高度相關(guān)性,因此用prcomp函數(shù)對種群所在地的氣候因子進行主成分分析,提取前三個主成分中載荷大于0.8的氣候因子,得到5個氣候因子(表1),用于相關(guān)性分析。

    1.3 表型性狀的測定

    對各種質(zhì)的樹高、胸徑、冠徑等樹體性狀進行測量,其中冠徑為東西冠徑和南北冠徑的平均值,每個性狀重復測量3次,取平均值作為實測值。

    計算各種質(zhì)的單株材積(v),公式為v=0.4π×(d1.3/2)×h+3×10(式中,d1.3為胸徑,h為樹高)(林昌庚,1964)。

    記錄樹干干型情況,對這一非數(shù)值性狀由3人分別進行1~10賦值統(tǒng)計。根據(jù)劉雨等(2011)的研究且結(jié)合赤皮青岡實際情況進行以下賦值:有明顯主干無分枝,樹干通直記為10;有明顯主干無分枝,但樹干輕微彎曲記為9;有明顯主干無分枝,但樹干彎曲明顯記為8;有明顯主干且出現(xiàn)1個分岔,但分岔和主干較為通直記為7;有明顯主干和1個分岔,樹干較為彎曲記為6;有明顯主干和多個分岔,主干分枝較直記為5;有主干但第1個分岔較早,枝下高小于1.3 m,主干分枝較直記為4;有主干但第1個分岔較早,枝下高小于1.3 m,主干分枝較彎曲記為3;無明顯主干,但主要枝條較直記為2;無明顯主干,枝條較彎曲記為1。

    使用ImageJ 圖像分析軟件分析各種質(zhì)的葉片圖像,測量獲取葉片的葉長、葉寬、葉面積、葉柄長、葉柄寬、葉柄面積、鋸齒數(shù)等數(shù)據(jù),重復3次,取平均值作為實測值。通過測量值計算葉長寬比、葉柄長寬比和葉脈間距。其中,葉脈間距計算公式為葉脈間距=葉長/(單邊葉脈數(shù)+1)(郭燕等,2020)。

    1.4 數(shù)據(jù)的處理和分析

    1.4.1 巢式方差分析和多重比較?在SPSS 26軟件中,采用巢式設(shè)計模型對赤皮青岡的15個表型性狀進行巢式方差分析,其線性模型為Yijk=m+ai+bj(i)+eijk。式中, Yijk為第 i個種群的第 j個單株的第 k個觀測值, m為總體平均值, ai為第 i個種群的效應(yīng)值,bj(i)為第 i個種群的第 j個單株的效應(yīng)值, eijk為第ijk個觀測值的試驗誤差。

    1.4.2 多樣性指數(shù)分析?在Excel軟件中,通過公式CV=δ/x× 100%和VST=δ2t/(δ2t+δ2s)計算各性狀的變異系數(shù)(CV)及表型分化系數(shù)(VST)(葛頌等,1988)。式中,δ2t表示種群間方差,δ2s表示種群內(nèi)方差。在Past 4.09中計算各種群表型性狀的Shannon-Wiener多樣性指數(shù),先根據(jù)平均值(mean,x)和標準差(SD,δ)將所有性狀數(shù)據(jù)分成 10級,轉(zhuǎn)換為質(zhì)量性狀,第一級為< x-2δ,中間每級增加0.5δ,直到第10級為≥ x+2δ,再計算Shannon-Wiener多樣性指數(shù)。

    1.4.3 相關(guān)性分析、主成分分析與聚類分析?使用SPSS 26軟件,通過Pearson相關(guān)系數(shù)分析法,研究赤皮青岡各表型性狀間以及性狀與地理位置和氣候因子間的相關(guān)性。在R軟件中,先通過prcomp函數(shù)表型數(shù)據(jù)進行標準化,再進行主成分分析?;谥鞒煞值梅窒禂?shù)矩陣,提取貢獻率最大的第一和第二主成分,采用平方歐式距離和組間連接法進行種群間聚類分析。赤皮青岡天然種群受破壞較嚴重,大多種群內(nèi)樣本數(shù)量較少,并且種群間樹齡結(jié)構(gòu)差異較大,靖州和松陽種群僅由3株古樹和少量次生幼樹組成,保靖、通道、正板、插柳、安化、龍山和周寧等種群內(nèi)僅存幾株古樹和中齡赤皮青岡,鳳凰、洞口、常寧、婺源和普陀種群為樹齡較大的同齡林。因此,在進行主成分分析、聚類分析以及表型性狀與地理和氣候因子相關(guān)性分析時,剔除樹高、胸徑、冠徑和材積4個受樹齡影響較大的表型性狀,使用其余11個表型性狀進行分析。

    2?結(jié)果與分析

    2.1 赤皮青岡表型性狀變異

    赤皮青岡15個表型性狀的變異系數(shù)平均值為35.070%(表2)。其中,材積的變異系數(shù)最大(為152.606%),變異幅度為0.05~58.72 m3。葉片長寬比變異幅度最小(為11.868%),變異幅度為2.540~4.534。15個表型性狀的變異系數(shù)由大到小依次為材積>胸徑>冠徑>樹高>葉柄面積>葉面積>葉柄長寬比>葉柄長>葉脈間距>鋸齒數(shù)>葉長>葉柄寬>葉寬>干型>葉長寬比。這說明赤皮青岡樹體性狀變異更加豐富。

    15個表型性狀的多樣性指數(shù)(H)在1.543~2.142之間,平均值為1.998(表2)。其中,材積的多樣性指數(shù)最低(為1.543),葉寬的多樣性指數(shù)最高(為2.142)。15個表型性狀的多樣性由高到低分別為葉寬>葉脈間距>葉長>葉長寬比>鋸齒數(shù)>葉柄寬>樹高>葉柄長>葉面積>葉柄長寬比>葉柄面積>冠徑>胸徑>干型>材積。這說明赤皮青岡在個體水平上存在較為豐富的表型變異。

    2.2 赤皮青岡表型變異來源及種群變異規(guī)律

    赤皮青岡15個表型性狀在種群間和種群內(nèi)的方差分析結(jié)果見表3。由表3的F檢驗結(jié)果可知,無論是在種群間還是種群內(nèi)15個表型性狀都存在極顯著差異(P<0.01),說明赤皮青岡的樹體和葉片性狀在種群間和種群內(nèi)的變異都較為豐富。

    由表4可知,保靖、插柳、安化、龍山、婺源、周寧、松陽種群的樹高、胸徑、冠徑和干型等樹體指標值都較大,其中松陽種群的樹高、材積、冠徑、干型均為最大,表現(xiàn)為樹形高大通直的特征;靖州種群的胸徑最大;鳳凰、洞口、常寧、普陀種群的平均材積相對較小,與各種群實際樹齡結(jié)構(gòu)相符。靖州、常寧、普陀種群的葉長寬比較小,葉片表現(xiàn)為短而寬的葉形;保靖、龍山、婺源種群的葉長寬比較大,葉片呈細長型;插柳種群的葉長和葉面積最大;正板赤皮青岡的葉片鋸齒數(shù)最多,葉面積較大,表現(xiàn)為葉片大且鋸齒多的特點。

    2.3 赤皮青岡種群內(nèi)表型變異特征

    變異系數(shù)(CV)表示種群內(nèi)表型性狀的變異情況,變異系數(shù)越大表明表型的離散程度越大。由表5可知,14個種群的變異系數(shù)在14.94%(洞口)~35.56%(龍山)之間,平均值為23.07%,說明在種群水平上赤皮青岡的表型性狀變異水平較高;14個種群中,僅洞口種群的15個表型性狀變異系數(shù)平均值未超過15%,說明該種群赤皮青岡表型變異水平較低;15個表型性狀中干型、葉長、葉寬、葉長寬比、葉柄寬和鋸齒數(shù)的變異系數(shù)平均值都未超過15%,說明這6個性狀較為穩(wěn)定。與表2結(jié)果一致,除干型外,樹高、胸徑、冠徑和材積等樹體性狀的變異系數(shù)平均值大于葉片性狀的變異系數(shù)平均值,說明赤皮青岡天然種質(zhì)的樹體性狀受樹齡影響過高而表現(xiàn)出較高變異水平。

    2.4 赤皮青岡種群內(nèi)表型多樣性

    Shannon-Wiener多樣性指數(shù)反映不同表型性狀的豐富和均勻程度,14個種群的15個表型性狀的多樣性指數(shù)見表6。由表6可知,Shannon-Wiener指數(shù)在1.127(松陽)~1.980(常寧)之間,14個種群的多樣性指數(shù)從高到低依次為常寧>普陀>洞口>插柳>保靖>鳳凰>周寧>通道>婺源>龍山>安化>正板>靖州>松陽。Shannon-Wiener多樣性指數(shù)與種群內(nèi)樣本數(shù)量大小有一定的關(guān)系,靖州、安化、松陽等樣本較少的種群表現(xiàn)出較低的Shannon-Wiener指數(shù),在種質(zhì)資源收集和保護工作中應(yīng)重視。

    2.5 赤皮青岡表型性狀變異來源及種群間表型分化

    表型分化系數(shù)(VST)是種群間表型變異占總表型的百分比,用于評估種群間的表型分化情況。由表7可知,15個赤皮青岡表型性狀的分化系數(shù)在23.420%(葉柄長)~57.359%(胸徑)之間。赤皮青岡種群間表型方差分量占總表型變異的31.682%,種群內(nèi)占44.677%,表型分化系數(shù)為41.491%,說明種群內(nèi)的變異在赤皮青岡表型變異中起主要作用。

    2.6 赤皮青岡表型性狀間的相關(guān)性及其與地理和氣候因子間的相關(guān)性

    15個赤皮青岡表型性狀間的相關(guān)性結(jié)果如表8所示。其中,部分性狀間存在顯著或極顯著相關(guān)性。樹高與胸徑、材積、干型、冠徑、葉長寬比呈極顯著正相關(guān)性,與葉長和葉脈間距呈顯著正相關(guān),說明較高的赤皮青岡往往具有較大的胸徑、較優(yōu)良的干型和冠形,葉片也較為狹長。干型與胸徑、材積、冠徑均呈極顯著正相關(guān),說明優(yōu)良的干型在一定程度上可以提高赤皮青岡的生長潛力。葉面積與葉長寬比無明顯相關(guān)性,說明不同種質(zhì)的葉長寬比表現(xiàn)出一定程度的穩(wěn)定性。葉寬與葉柄長、葉柄寬、葉柄面積均呈極顯著正相關(guān),葉長和葉面積與葉柄寬和葉柄面積呈極顯著正相關(guān),說明大的葉片通常伴隨葉柄粗長的特征。葉脈間距僅與鋸齒數(shù)呈極顯著負相關(guān),說明鋸齒越多葉脈間距越小。

    11個表型性狀與種群所在地的地理和氣候因子的相關(guān)性分析結(jié)果如表9所示。赤皮青岡表型性狀與地理氣候因子間的相關(guān)性不明顯。干型與最濕月份降水量、最濕季度降水量和最暖季度降水量呈顯著正相關(guān),葉寬與最暖季度降水量呈顯著負相關(guān),葉柄寬與最濕月份降水量和最暖季度降水量呈顯著負相關(guān);經(jīng)度、緯度、最冷月份最低溫度和最濕季度平均溫度與赤皮青岡表型性狀間無相關(guān)性。這表明降水對赤皮青岡表型性狀的影響大于溫度和地理因子。

    2.7 赤皮青岡表型性狀的主成分分析和聚類分析

    對赤皮青岡11個表型性狀進行主成分分析,得到4個特征根大于1的主成分,累計貢獻率為82.961% (表10)。由表10可知, PC1特征值為4.205,解釋總變異的38.231%, 主要包括葉長、 葉寬、葉面積(載荷的絕對值>0.8),主要表征葉片大?。籔C2的特征值為2.447,解釋了22.248%的總變異,主要包括葉柄長和葉柄長寬比,與葉柄性狀相關(guān);PC3特征值為1.344,主要表征干型,解釋了12.214%的總變異;PC4特征值為1.129,主要表征鋸齒數(shù),解釋了10.268%的總變異。由圖1可知,赤皮青岡表型性狀未呈現(xiàn)出明顯的地理變異規(guī)律,其中湖南鳳凰和福建周寧種群以其較低的PC1得分與其他種群區(qū)分開來。

    以平方歐式距離為14個赤皮青岡種群間的距離,采用組間連接法進行聚類(圖2)。當平方歐式距離為15時,根據(jù)11個表型性狀可將14個種群分為三類(與主成分分析結(jié)果相符):第一類為湖南鳳凰和福建周寧,主要特征為葉長、葉寬、葉面積、葉柄長、葉柄寬、葉柄面積和葉脈間距值均最?。ū?),表現(xiàn)為小葉片小葉柄類;第二類為湖南靖州、浙江松陽和湖南插柳,其葉寬、葉面積、葉柄長、葉柄寬、葉柄面積和葉脈間距的值最大,葉長寬比值最低,說明第二類呈現(xiàn)出葉片大、葉寬且葉柄粗長的特征,表現(xiàn)為大葉片大葉柄類;第三類包括湖南保靖、湖南通道、湖南正板、湖南洞口、湖南安化、湖南龍山、湖南常寧、江西婺源和浙江普陀9個種群,屬于第一類和第二類的銜接部分,葉片形態(tài)接近第二類,葉長、葉寬、葉面積的值普遍較大,但葉柄形態(tài)更接近第一類,普遍較小,總體表現(xiàn)為中葉片小葉柄類。通過比較種群的聚類結(jié)果和地理位置之間的關(guān)系,發(fā)現(xiàn)二者關(guān)聯(lián)并不顯著,說明赤皮青岡表型性狀與地理因子間的相關(guān)性不明顯。

    3?討論與結(jié)論

    3.1 赤皮青岡表型性狀的多樣性

    表型性狀的變異是環(huán)境和遺傳因子共同影響下的表達(曹德美等,2021),尤其是赤皮青岡這類珍稀植物,種質(zhì)資源流失嚴重且面臨復雜的生存環(huán)境和長期的地理隔離,極易發(fā)生樹體和葉片等表型性狀的變異。

    本研究結(jié)果表明,赤皮青岡15個表型性狀的Shannon-Wiener多樣性指數(shù)在1.543~2.142之間,平均值為1.998。高于瀕危植物云南黃連(1.634~1.839)(Coptis chinensis)(楊維澤等,2013),與瀕危植物格木(Erythrophleum fordii)(2.027 8)(李洪果等,2019)相近,低于四川珙桐(2.360)(Davidia involucrata)(張騰駒等,2019)。這表明赤皮青岡雖面臨種質(zhì)資源流失的問題,仍表現(xiàn)出較高的多樣性水平,特別是葉長、葉寬和葉脈間距等性狀。葉片的形態(tài)直接影響葉片生長發(fā)育特性,間接影響植物的適應(yīng)能力,赤皮青岡葉片表型性狀的多樣性說明其具有較強的環(huán)境適應(yīng)性。

    本研究中,赤皮青岡15個表型性狀的極大值是極小值的1.78(葉長寬比)~1 174.40(材積)倍,其中樹體性狀的極大值與極小值的比值普遍大于葉片性狀,表明樹體性狀較葉片性狀更易受到樹齡本身的影響。本研究中,赤皮青岡15個表型性狀的變異系數(shù)在11.868%~152.606%之間,天然種質(zhì)間表型性狀差異較大,尤其是樹高、胸徑、冠徑和材積等受樹齡影響大的樹體性狀;剔除以上4個性狀后,其余11個表型性狀的變異系數(shù)為11.868%~29.718%,與灰楸(Catalpa fargesii)(9.66%~30.39%) (于曉池等,2021) 和板栗(Castanea mollissima)(8.68%~32.73%)(郭燕等,2022)相近。這說明赤皮青岡干型和葉片性狀變異豐富,在不同環(huán)境因子和遺傳因子的影響下做出相應(yīng)的改變,這與Shannon-Wiener多樣性指數(shù)的結(jié)論一致,表明赤皮青岡在復雜生長環(huán)境中表現(xiàn)出較強的適應(yīng)能力。

    3.2 赤皮青岡表型性狀的變異來源和分化

    本研究中,赤皮青岡15個表型性狀在種群間和種群內(nèi)都表現(xiàn)出較豐富的變異,而種群內(nèi)的表型變異更為豐富(平均表型分化系數(shù)VST=41.491%)。這與大多數(shù)林木種質(zhì)資源的表型研究結(jié)果基本一致,如小果白刺(Nitraria sibirica)(40.71%)(董昕等,2017)、香椿(Toona sinensis)(40.209%)(吳軍等,2018)等樹種均存在較為豐富的表型變異。這表明赤皮青岡葉片等表型性狀在種群間的遺傳分化程度相對較低,可能是本研究所調(diào)查種群均位于我國亞熱帶地區(qū),分布區(qū)內(nèi)土壤理化性質(zhì)和氣候環(huán)境等條件差別不大,并且赤皮青岡在自然分布范圍內(nèi)呈碎片化分布的時間較短,暫未因生殖隔離而形成較大的差異。

    本研究中,赤皮青岡11個葉片及干型性狀無論是在種群間還是在種群內(nèi)都存在極顯著差異,而種群內(nèi)的表型變異為主要變異來源(表型分化系數(shù)在23.420%~56.757%之間)。同時,14個種群的11個葉片及干型性狀的變異系數(shù)平均值在7.35%(干型)~23.33%(葉柄面積)之間,并且同一表型性狀的Shannon-Wiener指數(shù)在種群間存在較大差異,提示該物種較易發(fā)生種群分化。本研究所調(diào)查的赤皮青岡天然種群大多零散分布于村落周圍,雖然個別天然種群保存較為完整,但大部分種群因受人為干擾較嚴重而造成天然生境的改變,加上生境不連續(xù)、樣本較少的零散分布,種群間的基因流受阻,從而引起種群分化和表型性狀變異。

    3.3 赤皮青岡表型變異與地理氣候因子間的相關(guān)性

    不同物種對環(huán)境因子的敏感度和適應(yīng)力有一定的差異,其表型性狀呈現(xiàn)出不同的地理變異規(guī)律,如板栗(郭燕等,2022)的葉片由南向北表現(xiàn)出長寬比變小的趨勢,唐古特白刺(吉德娟等,2021)的葉片厚度與海拔呈正相關(guān),青楊(Populus cathayana)(曹德美等,2021)的側(cè)脈數(shù)與經(jīng)緯度呈極顯著正相關(guān)。本研究中,在適宜赤皮青岡生長的環(huán)境下,經(jīng)緯度和溫度對赤皮青岡表型性狀的影響不顯著,而降水與赤皮青岡表型性狀間存在一定的相關(guān)性。這與姜小龍(2020)的研究結(jié)果類似,降水是影響南嶺青岡(Quercus championii)和福建青岡(Q. chungii)分布最大的氣候因子。本研究發(fā)現(xiàn),氣候因子中最濕月份降水量、最濕季度降水量以及最暖季度降水量僅與赤皮青岡干型呈正相關(guān),而最濕月份降水量和最暖季度降水量與葉寬和葉柄寬呈負相關(guān)。在適宜的降水范圍內(nèi),降水量的增加會促進維管形成層的活動,有助于樹干的均勻增粗,從而表現(xiàn)出優(yōu)良的干型。隨著降水量的增加,櫟類樹種的葉片和葉柄傾向于長條狀發(fā)展(李東勝等,2013),葉寬和葉柄寬表現(xiàn)出減小的趨勢,該現(xiàn)象與苦楝(Melia azedarach)(陳麗君等,2016)、板栗(郭燕等,2022)等的研究結(jié)果相似。因此,降水是影響赤皮青岡表型和分布最主要的因子,這為赤皮青岡的資源圃和人工林基地建設(shè)提供了科學依據(jù)。

    3.4 赤皮青岡的保護與開發(fā)利用策略

    赤皮青岡天然種群的生境不連續(xù)、跨度較大、碎片化較為嚴重,目前大多天然種群處于被孤立狀態(tài)且包含的種質(zhì)數(shù)量較少。赤皮青岡種子為頑拗性種子,外形較大,自然擴散距離十分有限,并且花粉的傳播受到我國南方復雜地形的阻礙,長此以往,高片段化分布的種群間基因流受阻(姜小龍,2020),加上環(huán)境異質(zhì)性的影響,從而引起種群分化和表型性狀變異。因此,僅開展就地保護既不能有效地為赤皮青岡種群間的基因交流提供廊道,也不利于種群的延續(xù)和遺傳多樣性的維持。為了更好地保護、管理和利用赤皮青岡種質(zhì)資源,應(yīng)加強對現(xiàn)有天然種群和單株的保護力度,特別是單株數(shù)量稀少的種群;加快開展種質(zhì)資源收集工作,鑒于赤皮青岡種群內(nèi)表型變異水平較高,要注重對種群內(nèi)每個單株的收集;應(yīng)將赤皮青岡作為鄉(xiāng)土樹種在其天然分布區(qū)內(nèi)推廣造林,以適應(yīng)能力較強的大葉片大葉柄類的靖州、松陽和插柳種群為種源進行采種育苗和造林,既可提高赤皮青岡的適應(yīng)性,又可改善種群片段化分布的現(xiàn)狀,為種群間的交流提供廊道,幫助該物種的延續(xù)。

    參考文獻:

    BAI QS, CAI YL, HE BX, et al., 2019. Core set construction and association analysis of Pinus massoniana from Guangdong Province in southern China using SLAF-seq[J]. Sci Rep, 9: 13157.

    CAO DM, ZHANG YH, CHENG XQ, et al., 2021. Genetic variation of leaf phenotypic traits in different populations of Populus cathayana [J]. Sci Silv Sin, 57(8): 56-67.[曹德美, 張亞紅, 成星奇, 等, 2021. 青楊不同種群葉片表型性狀的遺傳變異 [J]. 林業(yè)科學, 57(8): 56-67.]

    CHEN C, DING CJ, HUANG QJ, et al., 2021. Construction of phenotypic core collection of Populus deltoides [J]. For Res, 34(2): 1-11.[陳存, 丁昌俊, 黃秦軍, 等, 2021. 美洲黑楊表型核心種質(zhì)庫構(gòu)建 [J]. 林業(yè)科學研究, 34(2): 1-11.]

    CHEN LJ, LIU MQ, LIAO BY, et al., 2016. Variation of leaf morphological characters of Melia azedarach from different provenances [J]. J SW For Univ, 36(5): 16-20.[陳麗君, 劉明騫, 廖柏勇, 等, 2016. 苦楝不同種源葉片性狀變異研究 [J]. 西南林業(yè)大學學報, 36(5): 16-20.]

    DENG M, JIANG XL, HIPP AL, et al., 2018. Phylogeny and biogeography of East Asian evergreen oaks (Quercus section Quercus; Fagaceae): Insights into the Cenozoic history of evergreen broad-leaved forests in subtropical Asia [J]. Mol Phylogenet Evol, 119: 170-181.

    DONG X, WANG L, LU YZ, et al., 2017. Phenotypic variation of Nitraria sibirica natural populations in Shandong [J]. For Res, 30(2): 293-299.[董昕, 王磊, 魯儀增, 等, 2017. 山東稀有植物小果白刺天然群體表型變異研究 [J]. 林業(yè)科學研究, 30(2): 293-299.]

    DUAN M, 1995. Comprehensive utilization of Quercus resources in China [J]. J Beijing For Univ, 17(2): 109-110.[端木, 1995. 我國青岡屬資源的綜合利用 [J]. 北京林業(yè)大學學報, 17(2): 109-110.]

    FICK SE, HIJMANS RJ, 2017. WorldClim 2: new 1km spatial resolution climate surfaces for global land areas [J]. Int J Clim, 37(12): 4302-4315.

    GE S, WANG MX, CHEN YW, 1988. An analysis of population genetic structure of Masson pine by isozyme tech-nique [J]. Sci Silv Sin, 24(4): 399-409. [葛頌, 王明庥, 陳岳武, 1988. 用同工酶研究馬尾松群體的遺傳結(jié)構(gòu) [J]. 林業(yè)科學, 24(4): 399-409.]

    GUO Y, ZHANG SH, LI Y, et al., 2022. Diversity analysis of 36 leaf phenotypic traits of Chinese chestnut [J]. Sci Agric Sin, 55(5): 991-1009.[郭燕, 張樹航, 李穎, 等, 2022. 中國板栗36個葉片表型性狀的多樣性 [J]. 中國農(nóng)業(yè)科學, 55(5): 991-1009.]

    GUO Y, ZHANG SH, LI Y, et al., 2020. Studies on the leaf morphology, anatomical structure and drought resistance evaluation of 238 Chinese chestnut varieties (Strains) [J]. Acta Hortic Sin, 47(6): 1033-1046.[郭燕, 張樹航, 李穎, 等, 2020. 中國板栗238份品種(系)葉片形態(tài)、解剖結(jié)構(gòu)及其抗旱性評價 [J]. 園藝學報, 47(6): 1033-1046.]

    JIANG XL, 2020. Phylogenetic relationship and population genetic structure of Quercus chungii and Q. championii[D]. Changsha: Central South University of Forestry and Technology.[姜小龍, 2020. 福建青岡和嶺南青岡系統(tǒng)發(fā)育關(guān)系及居群遺傳結(jié)構(gòu) [D]. 長沙: 中南林業(yè)科技大學.]

    JI DJ, ZHANG DF, YU Q, 2021. Phenotypic diversity of Nitraria tangutorum Bobr. in Qaidam basin [J]. J Centr S Univ For Technol, 41(10): 57-66.[吉德娟, 張得芳, 于倩, 2021. 柴達木盆地唐古特白刺的表型多樣性 [J]. 中南林業(yè)科技大學學報, 41(10): 57-66.]

    LIN CG, 1964. Stem control in forest volume measurement technology [J]. Sci Silv Sin, (4): 87-97.[林昌庚, 1964. 林木蓄積量測算技術(shù)中的干形控制問題 [J]. 林業(yè)科學, (4): 87-97.]

    LI YG, LIU XH, MA JW, et al., 2018. Phenotypic variation in Phoebe bournei populations preserved in the primary distribution area [J]. J For Res, 29(1): 35-44.

    LI DS, SHI ZM, FENG QH, et al., 2013.Response of leaf morphometric traits of Quercus species to climate in the temperate zone of the North-South Transect of Eastern China [J]. Acta Phytoecol Sin, 37(9): 793-802.[李東勝, 史作民, 馮秋紅, 等, 2013. 中國東部南北樣帶暖溫帶區(qū)櫟屬樹種葉片形態(tài)性狀對氣候條件的響應(yīng) [J]. 植物生態(tài)學報, 37(9): 793-802.]

    LI S, GAN XH, HAN HY, et al., 2016. Leaf phenotypic traits of Tetracentron sinense, an endangered plant species [J]. For Res, 29(5): 687-697.[李珊, 甘小洪, 憨宏艷, 等, 2016. 瀕危植物水青樹葉的表型性狀變異 [J]. 林業(yè)科學研究, 29(5): 687-697.]

    LIU Y, LI DW, QIN YS, et al., 2011. Phenotypic diversity for natural populations of Juniperus rigida in the Helan Mountains of Ningxia [J]. J Zhejiang A & F Univ, 28(4): 619-627. [劉雨, 李登武, 秦廷松, 等, 2011. 寧夏賀蘭山杜松天然群體的表型多樣性 [J]. 浙江農(nóng)林大學學報, 28(4): 619-627.]

    LIN WJ, LIU QY, LUO X, et al., 2017. Research on species diversity and niche of natural small population of Quercus glauca [J]. Fujian For, (6): 38-41. [林文俊, 劉沁月, 羅蕭, 等, 赤皮青岡天然小居群物種多樣性與生態(tài)位研究 [J]. 福建林業(yè), (6): 38-41.]

    OU HB, DONG LJ, LI J, et al., 2021. Relationship between phenotypic diversity of Phoebe bournei deed traits and environmental geographical factors [J]. J NE For Univ, 49(8): 45-50.[歐漢彪, 董利軍, 李娟, 等, 2021. 閩楠種實性狀表型多樣性與環(huán)境地理因子的關(guān)系 [J]. 東北林業(yè)大學學報, 49(8): 45-50.]

    OUYANG ZY, OUYANG SL, WU JY, et al., 2021. Research progress of precious commercial tree species Quercus gilva [J]. Hunan For Sci Technol, 48(6): 74-79.[歐陽澤怡, 歐陽碩龍, 吳際友, 等, 2021. 珍貴用材樹種赤皮青岡研究進展 [J]. 湖南林業(yè)科技, 48(6): 74-79.]

    SHANG SB, 2015. Genetic diversity in natural populations of Vatica mangachapoi on Hainan Island, China [D]. Beijing: Chinese Academy of Forestry.[尚帥斌, 2015. 海南島青梅天然居群遺傳多樣性 [D]. 北京: 中國林業(yè)科學研究院.]

    SUGIURA N, TANG DQ, KUROKOCHI H, et al., 2015. Genetic structure of Quercus gilva Blume in Japan as revealed by chloroplast DNA sequences [J]. Botany, 93(12): 873-880.

    WU J, ZHONG ZZ, LOU J, et al., 2018. A study on phenotypic diversity of cone and seed in natural populations of Toona sinensis [J]. Acta Agric Univ Jiangxi, 40(2): 248-256.[吳軍, 鐘志真, 樓軍, 等, 2018. 香椿天然群體種實性狀表型多樣性研究 [J]. 江西農(nóng)業(yè)大學學報, 40(2):248-256.]

    WU YH, FAN ZL, LI J, et al., 2018. Phenotypic diversity of seeds and fruits in natural populations of Acer ginnala in China [J]. Guihaia, 38(6): 795-803.[武艷虹, 樊澤璐, 李佳, 等, 2018. 茶條槭自然種群種子和果實表型多樣性研究 [J]. 廣西植物, 38(6): 795-803.]

    XIA C, 2021. Effects of shading on the growth and physiological characteristics of Quercus gilva seedlings [D]. Changsha: Central South University of Forestry and Technology.[夏嬋, 2021. 遮陰對赤皮青岡幼苗生長及生理特性的影響 [D]. 長沙: 中南林業(yè)科技大學.]

    XIE J, 2021. Variations in growth traits and provenance selection of Quercus gilv [J]. J Fujian Agric For Univ (Nat Sci Ed), 50(5): 636-640.[謝健, 2021. 不同種源赤皮青岡幼林生長性狀變異及種源選擇 [J]. 福建農(nóng)林大學學報(自然科學版), 50(5): 636-640.]

    XIE J, 2011. Research on population ecology of Quercus gilvs [D]. Fuzhou: Fujian Agriculture and Forestry University.[謝健, 2011. 赤皮青岡種群生態(tài)學研究 [D]. 福州: 福建農(nóng)林大學.]

    YANG WZ, JIN H, LI WY, et al., 2013. Phenotypic diversity of different populations of endangered plant Coptis yunnanensis [J]. J Yunnan Univ (Nat Sci Ed), 35(5): 719-726.[楊維澤, 金航, 李晚誼, 等, 2013. 瀕危植物云南黃連不同居群表型多樣性研究 [J]. 云南大學學報(自然科學版), 35(5): 719-726.]

    YU XC, LI F, OU Y, et al., 2021. Construction of core germplasm of Sorbus mandshurica based on phenotype [J]. For Res, 34(6): 38-45.[于曉池, 李鳳, 歐陽, 等, 2021. 基于表型的灰楸核心種質(zhì)構(gòu)建 [J]. 林業(yè)科學研究, 34(6): 38-45.]

    ZAYNAB M, PAN D, NOMAN A, et al., 2018. Transcriptome approach to address low seed germination in Quercus gilva to save forest ecology [J]. Biochem Syst Ecol, 81: 62-69.

    ZENG QM, LIU B, LIU RQ, et al., 2019. The complete chloroplast genome sequence of Quercus gilva (Fagaceae)[J]. Mitchondrial DNA B, 4(2): 2493-2494.

    ZHANG TJ, CHEN XH, KANG XK, et al., 2019. Phenotypic diversity of leaf morphologic traits of Davidia involucrata natural populations in Sichuan Province [J]. Sin J Acta Ecol, 38(1): 35-43.[張騰駒, 陳小紅, 康喜坤, 等, 2019. 四川省珙桐天然種群葉表型多樣性 [J]. 生態(tài)學雜志, 38(1): 35-43.]

    ZHU PH, 2014. Analysis on genetic diversity and genetic structure of Quercus gilva population [D]. Changsha: Central South University of Forestry and Technology.[朱品紅, 2014. 赤皮青岡居群遺傳多樣性與遺傳結(jié)構(gòu)分析 [D]. 長沙: 中南林業(yè)科技大學.]

    (責任編輯?蔣巧媛?王登惠)

    猜你喜歡
    表型性狀種質(zhì)資源
    瓠瓜類型砧木重組近交系(RIL)遺傳作圖群體耐濕澇表型性狀分析
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    云南德宏羊奶果資源收集和評價簡報
    綠肥作物紫云英研究進展
    大白菜種質(zhì)資源抗根腫病基因CRa和CRb的分子標記鑒定與分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    玉米種質(zhì)資源抗旱性鑒定研究進展
    楊梅種質(zhì)資源遺傳多樣性研究進展
    茄子重要性狀表型多樣性及其在育種研究中的應(yīng)用
    日韩大片免费观看网站| 1024香蕉在线观看| 日本精品一区二区三区蜜桃| 亚洲精品中文字幕一二三四区 | 自拍欧美九色日韩亚洲蝌蚪91| 精品一区在线观看国产| 在线 av 中文字幕| 人人妻,人人澡人人爽秒播| 欧美日本中文国产一区发布| 亚洲国产欧美一区二区综合| av网站在线播放免费| 久久久久久久久免费视频了| 中文字幕av电影在线播放| 色视频在线一区二区三区| 免费av中文字幕在线| 国产成人欧美在线观看 | 国产日韩欧美在线精品| 一区福利在线观看| 亚洲成人国产一区在线观看| 精品少妇黑人巨大在线播放| 国产男女超爽视频在线观看| 美女大奶头黄色视频| 黄片大片在线免费观看| 国产极品粉嫩免费观看在线| 最新在线观看一区二区三区| 法律面前人人平等表现在哪些方面 | 亚洲成国产人片在线观看| 91麻豆av在线| 久久精品亚洲熟妇少妇任你| 十八禁高潮呻吟视频| videosex国产| 亚洲国产av新网站| 热re99久久国产66热| 久久久久久久大尺度免费视频| 午夜福利,免费看| 国产成人影院久久av| www.熟女人妻精品国产| 国产精品免费视频内射| 两个人免费观看高清视频| 大型av网站在线播放| 狂野欧美激情性xxxx| 国产成+人综合+亚洲专区| 午夜免费成人在线视频| 久久中文字幕一级| 成人国语在线视频| 一本久久精品| 国产一区有黄有色的免费视频| 日韩视频一区二区在线观看| 性少妇av在线| 精品少妇内射三级| 国产日韩欧美亚洲二区| 国产精品av久久久久免费| 国产一区二区在线观看av| 少妇猛男粗大的猛烈进出视频| 国产高清视频在线播放一区 | 国产亚洲精品一区二区www | 老熟女久久久| 国产男女超爽视频在线观看| 夜夜夜夜夜久久久久| 欧美日韩国产mv在线观看视频| 男人舔女人的私密视频| 色94色欧美一区二区| 日韩制服骚丝袜av| 91九色精品人成在线观看| 日本五十路高清| 最近中文字幕2019免费版| 少妇人妻久久综合中文| av电影中文网址| 日本91视频免费播放| 叶爱在线成人免费视频播放| 一本—道久久a久久精品蜜桃钙片| 午夜福利影视在线免费观看| 巨乳人妻的诱惑在线观看| 精品福利观看| 亚洲少妇的诱惑av| 精品国产乱码久久久久久男人| 两人在一起打扑克的视频| 一本综合久久免费| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 亚洲欧洲日产国产| 国产精品免费视频内射| 国产伦理片在线播放av一区| 欧美人与性动交α欧美精品济南到| 免费在线观看日本一区| 嫩草影视91久久| 叶爱在线成人免费视频播放| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡| av视频免费观看在线观看| 欧美亚洲 丝袜 人妻 在线| 青青草视频在线视频观看| 国产亚洲一区二区精品| 中文字幕人妻丝袜一区二区| 精品国产乱码久久久久久小说| 中文字幕高清在线视频| 乱人伦中国视频| 精品国产一区二区三区久久久樱花| 十分钟在线观看高清视频www| 亚洲欧美清纯卡通| 国产淫语在线视频| 亚洲人成电影观看| 一级毛片女人18水好多| 三上悠亚av全集在线观看| 色94色欧美一区二区| 日韩大片免费观看网站| 水蜜桃什么品种好| 高清视频免费观看一区二区| 久久精品久久久久久噜噜老黄| 欧美日韩国产mv在线观看视频| 国产精品久久久av美女十八| 交换朋友夫妻互换小说| 国产片内射在线| 老汉色∧v一级毛片| 女人久久www免费人成看片| 精品一区二区三区四区五区乱码| 无限看片的www在线观看| 日本91视频免费播放| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧洲日产国产| 国精品久久久久久国模美| 青青草视频在线视频观看| 精品一区二区三区av网在线观看 | 青草久久国产| 国产成人精品无人区| 国产精品久久久久久人妻精品电影 | 在线永久观看黄色视频| 久久久国产精品麻豆| 亚洲欧美一区二区三区黑人| 久久久久久久久久久久大奶| av福利片在线| 国产福利在线免费观看视频| 亚洲精品久久久久久婷婷小说| 国产精品偷伦视频观看了| 久久九九热精品免费| 午夜久久久在线观看| 男男h啪啪无遮挡| 女性被躁到高潮视频| 日韩一区二区三区影片| 亚洲 国产 在线| 美女中出高潮动态图| 亚洲,欧美精品.| 亚洲视频免费观看视频| 日本黄色日本黄色录像| 日韩视频一区二区在线观看| 亚洲精品自拍成人| 三级毛片av免费| 国产精品欧美亚洲77777| 黑人欧美特级aaaaaa片| 一二三四在线观看免费中文在| 国产淫语在线视频| 国产免费福利视频在线观看| 免费不卡黄色视频| 精品第一国产精品| 久热爱精品视频在线9| 亚洲国产看品久久| 亚洲人成77777在线视频| 免费高清在线观看视频在线观看| 免费观看a级毛片全部| 久久午夜综合久久蜜桃| 欧美另类亚洲清纯唯美| 一级片'在线观看视频| 一级片免费观看大全| 日本黄色日本黄色录像| 久久久久久久大尺度免费视频| 高潮久久久久久久久久久不卡| 久久精品亚洲av国产电影网| 欧美老熟妇乱子伦牲交| 国产亚洲一区二区精品| 淫妇啪啪啪对白视频 | 中国国产av一级| 丝袜美足系列| 岛国在线观看网站| 男人爽女人下面视频在线观看| av一本久久久久| 又大又爽又粗| 中文字幕人妻丝袜制服| 男女免费视频国产| 女警被强在线播放| 黄片播放在线免费| 岛国在线观看网站| 午夜免费成人在线视频| 99re6热这里在线精品视频| 女人久久www免费人成看片| 超色免费av| 欧美日韩一级在线毛片| 亚洲激情五月婷婷啪啪| 久久亚洲国产成人精品v| 亚洲精品国产区一区二| 国产精品香港三级国产av潘金莲| 在线永久观看黄色视频| 在线看a的网站| 极品人妻少妇av视频| 中文字幕高清在线视频| a 毛片基地| www.自偷自拍.com| 日韩电影二区| 狠狠婷婷综合久久久久久88av| 国产免费现黄频在线看| 男女床上黄色一级片免费看| 精品视频人人做人人爽| 性色av乱码一区二区三区2| 亚洲 国产 在线| 久久天堂一区二区三区四区| 精品福利永久在线观看| 国产成人欧美在线观看 | 1024香蕉在线观看| 免费看十八禁软件| 亚洲精品国产色婷婷电影| 精品少妇一区二区三区视频日本电影| 久久亚洲精品不卡| 大型av网站在线播放| www.精华液| av免费在线观看网站| 精品人妻在线不人妻| 国产精品久久久久久人妻精品电影 | 精品乱码久久久久久99久播| 国产伦人伦偷精品视频| 久久这里只有精品19| 免费在线观看日本一区| 免费在线观看影片大全网站| 亚洲精品美女久久久久99蜜臀| 免费人妻精品一区二区三区视频| 亚洲国产欧美在线一区| 老司机亚洲免费影院| 日本欧美视频一区| 国产成人精品久久二区二区91| 亚洲第一欧美日韩一区二区三区 | 欧美一级毛片孕妇| 午夜日韩欧美国产| 无遮挡黄片免费观看| 国产精品自产拍在线观看55亚洲 | 777米奇影视久久| 婷婷色av中文字幕| 啪啪无遮挡十八禁网站| 亚洲av成人一区二区三| 深夜精品福利| 欧美另类亚洲清纯唯美| videos熟女内射| 亚洲国产欧美一区二区综合| 国产激情久久老熟女| 午夜久久久在线观看| 欧美另类亚洲清纯唯美| 人妻人人澡人人爽人人| 天堂俺去俺来也www色官网| 亚洲精品美女久久久久99蜜臀| 老司机午夜十八禁免费视频| 久久精品亚洲熟妇少妇任你| 狠狠精品人妻久久久久久综合| 精品国产乱码久久久久久男人| 国产欧美日韩综合在线一区二区| 亚洲精品久久午夜乱码| 国产亚洲av高清不卡| 精品一区二区三卡| 91av网站免费观看| 亚洲国产av影院在线观看| 黄色怎么调成土黄色| 国产精品久久久久久精品电影小说| 狠狠精品人妻久久久久久综合| 19禁男女啪啪无遮挡网站| 波多野结衣一区麻豆| 中文字幕最新亚洲高清| 啦啦啦啦在线视频资源| 十八禁网站免费在线| 麻豆国产av国片精品| 91精品国产国语对白视频| 日韩欧美国产一区二区入口| 美女脱内裤让男人舔精品视频| 爱豆传媒免费全集在线观看| 9热在线视频观看99| 亚洲精品一二三| 一区在线观看完整版| 热re99久久国产66热| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 亚洲精品在线美女| 人成视频在线观看免费观看| 亚洲欧美清纯卡通| 国产在线视频一区二区| 久久热在线av| 久久精品aⅴ一区二区三区四区| 亚洲国产av影院在线观看| 妹子高潮喷水视频| 中国美女看黄片| 亚洲自偷自拍图片 自拍| 一二三四社区在线视频社区8| 国产精品免费大片| 巨乳人妻的诱惑在线观看| 菩萨蛮人人尽说江南好唐韦庄| 不卡一级毛片| a在线观看视频网站| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花| 久久99热这里只频精品6学生| 热99久久久久精品小说推荐| 大片电影免费在线观看免费| 亚洲欧美色中文字幕在线| 久久这里只有精品19| 飞空精品影院首页| 欧美成狂野欧美在线观看| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 亚洲国产欧美网| 亚洲精品美女久久av网站| 久久久久精品国产欧美久久久 | a在线观看视频网站| 9热在线视频观看99| 免费观看av网站的网址| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 在线观看免费高清a一片| 久久av网站| 国产片内射在线| 日韩 欧美 亚洲 中文字幕| 久久久久久亚洲精品国产蜜桃av| 狠狠精品人妻久久久久久综合| 亚洲美女黄色视频免费看| 久久国产精品影院| cao死你这个sao货| 国产精品熟女久久久久浪| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜添小说| 人妻一区二区av| 久久久久精品国产欧美久久久 | 日韩大片免费观看网站| 日韩欧美国产一区二区入口| 一区二区av电影网| 天天躁夜夜躁狠狠躁躁| 亚洲熟女精品中文字幕| 亚洲精品在线美女| 国产av又大| 一区二区三区激情视频| 亚洲精品久久久久久婷婷小说| 曰老女人黄片| 国产成人精品久久二区二区91| 国产精品久久久久成人av| 国产精品一区二区在线观看99| 久久青草综合色| 亚洲久久久国产精品| 精品第一国产精品| 国产精品二区激情视频| 黄片播放在线免费| 在线精品无人区一区二区三| 久久青草综合色| 在线精品无人区一区二区三| 久久99热这里只频精品6学生| 99精品欧美一区二区三区四区| av一本久久久久| 亚洲国产看品久久| 国产精品欧美亚洲77777| 在线观看免费视频网站a站| 精品久久蜜臀av无| 在线观看一区二区三区激情| 老司机福利观看| 久久精品国产a三级三级三级| 母亲3免费完整高清在线观看| 国产精品欧美亚洲77777| 97精品久久久久久久久久精品| 十八禁网站网址无遮挡| 99久久精品国产亚洲精品| 亚洲欧美精品综合一区二区三区| 岛国毛片在线播放| 色综合欧美亚洲国产小说| 久久久久久久大尺度免费视频| 纵有疾风起免费观看全集完整版| 黄色怎么调成土黄色| 日韩 欧美 亚洲 中文字幕| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看 | 精品久久久久久电影网| 久久国产精品影院| 精品国产一区二区三区久久久樱花| 亚洲国产精品一区三区| 在线观看一区二区三区激情| 在线观看免费午夜福利视频| 国产亚洲av高清不卡| 黄频高清免费视频| 9191精品国产免费久久| 精品久久久久久久毛片微露脸 | 中文字幕人妻丝袜一区二区| xxxhd国产人妻xxx| 国产三级黄色录像| 亚洲欧美一区二区三区久久| 天堂中文最新版在线下载| 久久精品亚洲熟妇少妇任你| 电影成人av| 欧美精品亚洲一区二区| 国产黄频视频在线观看| 中文字幕制服av| 三级毛片av免费| 午夜免费成人在线视频| 777米奇影视久久| 久久精品熟女亚洲av麻豆精品| 美女扒开内裤让男人捅视频| 国产91精品成人一区二区三区 | 欧美日韩黄片免| 久9热在线精品视频| 欧美日韩亚洲高清精品| 久久久精品国产亚洲av高清涩受| 少妇精品久久久久久久| 美女扒开内裤让男人捅视频| 黄片小视频在线播放| 人人妻人人添人人爽欧美一区卜| 久久久国产精品麻豆| 天天躁日日躁夜夜躁夜夜| 91av网站免费观看| 欧美精品亚洲一区二区| 欧美成狂野欧美在线观看| 午夜视频精品福利| 国产成人免费观看mmmm| 在线观看免费视频网站a站| 久久久久精品国产欧美久久久 | 国产精品 欧美亚洲| 国产精品久久久av美女十八| av一本久久久久| 久久av网站| 啦啦啦在线免费观看视频4| 女人久久www免费人成看片| 国产免费一区二区三区四区乱码| 成年美女黄网站色视频大全免费| a 毛片基地| 汤姆久久久久久久影院中文字幕| 午夜视频精品福利| 欧美黄色淫秽网站| 成年动漫av网址| 国产人伦9x9x在线观看| 久久久久精品人妻al黑| 天天躁日日躁夜夜躁夜夜| 亚洲av日韩在线播放| av免费在线观看网站| 超色免费av| 亚洲国产成人一精品久久久| 两人在一起打扑克的视频| 欧美97在线视频| 国产片内射在线| 99国产精品免费福利视频| 美女福利国产在线| 久久久久久久国产电影| 别揉我奶头~嗯~啊~动态视频 | 多毛熟女@视频| 91麻豆精品激情在线观看国产 | 黄频高清免费视频| 亚洲欧美日韩另类电影网站| 亚洲色图 男人天堂 中文字幕| 国产欧美日韩一区二区三区在线| 国产成人欧美在线观看 | 91成人精品电影| 国产成人欧美在线观看 | 免费少妇av软件| 一级毛片电影观看| 日韩三级视频一区二区三区| 亚洲专区字幕在线| 日韩人妻精品一区2区三区| 亚洲熟女毛片儿| 国产熟女午夜一区二区三区| 国产日韩欧美亚洲二区| 日本猛色少妇xxxxx猛交久久| 精品乱码久久久久久99久播| 男男h啪啪无遮挡| 中文字幕另类日韩欧美亚洲嫩草| 国产精品免费视频内射| 后天国语完整版免费观看| 国产高清视频在线播放一区 | 1024香蕉在线观看| 91精品伊人久久大香线蕉| 香蕉丝袜av| 狂野欧美激情性bbbbbb| 在线观看免费午夜福利视频| 亚洲第一av免费看| 黄片播放在线免费| 丝袜美腿诱惑在线| 国产精品免费视频内射| 亚洲人成电影免费在线| 国产又爽黄色视频| 两个人看的免费小视频| 国产成人精品在线电影| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| 精品人妻1区二区| 欧美中文综合在线视频| 精品人妻熟女毛片av久久网站| 天天躁日日躁夜夜躁夜夜| av超薄肉色丝袜交足视频| 亚洲精品第二区| 69精品国产乱码久久久| 欧美av亚洲av综合av国产av| 91大片在线观看| 大香蕉久久网| 亚洲午夜精品一区,二区,三区| 91国产中文字幕| 精品人妻1区二区| 欧美国产精品一级二级三级| 国产三级黄色录像| 曰老女人黄片| 在线观看免费午夜福利视频| 男人添女人高潮全过程视频| 久久青草综合色| 性色av乱码一区二区三区2| 日韩一区二区三区影片| 性少妇av在线| 亚洲视频免费观看视频| netflix在线观看网站| av网站免费在线观看视频| 国产成人欧美| 少妇的丰满在线观看| 国产伦理片在线播放av一区| av电影中文网址| 丝袜人妻中文字幕| 女人高潮潮喷娇喘18禁视频| 欧美另类一区| 国产成人精品在线电影| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 欧美午夜高清在线| 国产成人精品久久二区二区免费| a级毛片在线看网站| 欧美日韩亚洲高清精品| 2018国产大陆天天弄谢| 十八禁网站网址无遮挡| 女性被躁到高潮视频| 免费看十八禁软件| 日本五十路高清| 欧美精品av麻豆av| 欧美激情极品国产一区二区三区| 一边摸一边抽搐一进一出视频| 国产熟女午夜一区二区三区| 亚洲国产看品久久| 日日摸夜夜添夜夜添小说| 天天操日日干夜夜撸| 999精品在线视频| 亚洲第一欧美日韩一区二区三区 | 美女脱内裤让男人舔精品视频| 九色亚洲精品在线播放| 国产一区二区在线观看av| 国产精品国产av在线观看| 真人做人爱边吃奶动态| 欧美精品人与动牲交sv欧美| 成人国产一区最新在线观看| 国产精品 欧美亚洲| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 欧美精品一区二区免费开放| 蜜桃国产av成人99| 亚洲人成电影免费在线| 老汉色∧v一级毛片| 中文字幕人妻熟女乱码| 欧美+亚洲+日韩+国产| 极品人妻少妇av视频| 免费黄频网站在线观看国产| 色综合欧美亚洲国产小说| 在线观看免费视频网站a站| www.熟女人妻精品国产| 欧美激情 高清一区二区三区| 法律面前人人平等表现在哪些方面 | www.999成人在线观看| 一区二区三区四区激情视频| 国产亚洲欧美在线一区二区| 亚洲国产欧美网| 狠狠婷婷综合久久久久久88av| 在线十欧美十亚洲十日本专区| 亚洲一码二码三码区别大吗| 99国产精品一区二区三区| 纵有疾风起免费观看全集完整版| 日韩 亚洲 欧美在线| 一区二区三区精品91| 欧美日韩黄片免| 99久久国产精品久久久| av天堂在线播放| 精品欧美一区二区三区在线| 啦啦啦免费观看视频1| 亚洲国产中文字幕在线视频| 久久亚洲国产成人精品v| 窝窝影院91人妻| 久久女婷五月综合色啪小说| 亚洲精品成人av观看孕妇| 欧美黄色淫秽网站| 日韩,欧美,国产一区二区三区| 午夜激情av网站| 性少妇av在线| 99re6热这里在线精品视频| 黑人操中国人逼视频| 亚洲七黄色美女视频| 欧美成人午夜精品| 日本vs欧美在线观看视频| 久久中文看片网| 精品卡一卡二卡四卡免费| 天天操日日干夜夜撸| 热re99久久精品国产66热6| av又黄又爽大尺度在线免费看| 成年人免费黄色播放视频| 老司机在亚洲福利影院| 夜夜夜夜夜久久久久| 亚洲精品美女久久av网站| 一二三四在线观看免费中文在| 国产97色在线日韩免费| 老汉色av国产亚洲站长工具| 欧美日韩av久久| 桃花免费在线播放| 国产欧美日韩精品亚洲av| a 毛片基地| 十分钟在线观看高清视频www| 亚洲国产日韩一区二区| 久久国产精品影院| 可以免费在线观看a视频的电影网站| 国产在线免费精品| 秋霞在线观看毛片| 777米奇影视久久| 久久女婷五月综合色啪小说| 精品国产超薄肉色丝袜足j| 国产激情久久老熟女| 亚洲九九香蕉| 亚洲国产欧美日韩在线播放| 国产精品 欧美亚洲| 免费女性裸体啪啪无遮挡网站| 在线亚洲精品国产二区图片欧美|