• <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)用
    九色亚洲精品在线播放| 午夜日韩欧美国产| av福利片在线| 午夜福利免费观看在线| 亚洲成人免费av在线播放| 亚洲欧美一区二区三区久久| 亚洲精品乱久久久久久| 精品国产一区二区三区久久久樱花| 久久久久久人人人人人| 国产精品电影一区二区三区 | 国产亚洲欧美98| 黄频高清免费视频| 精品一品国产午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 在线看a的网站| 精品人妻在线不人妻| 亚洲熟女毛片儿| 久久国产亚洲av麻豆专区| aaaaa片日本免费| 欧美日韩亚洲综合一区二区三区_| 成年动漫av网址| 久热爱精品视频在线9| 亚洲综合色网址| 一级,二级,三级黄色视频| 美女高潮喷水抽搐中文字幕| 亚洲 国产 在线| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久国产一级毛片高清牌| 精品一区二区三区四区五区乱码| 黄频高清免费视频| 新久久久久国产一级毛片| 757午夜福利合集在线观看| 多毛熟女@视频| 国产熟女午夜一区二区三区| 久久中文字幕一级| 人妻 亚洲 视频| 国产深夜福利视频在线观看| 中文字幕高清在线视频| 国产区一区二久久| 99riav亚洲国产免费| 大香蕉久久网| avwww免费| 91字幕亚洲| 国产乱人伦免费视频| 黄色丝袜av网址大全| 中出人妻视频一区二区| 久久人妻熟女aⅴ| tube8黄色片| 午夜激情av网站| 国产1区2区3区精品| 亚洲成a人片在线一区二区| 交换朋友夫妻互换小说| 操美女的视频在线观看| 日本黄色日本黄色录像| 人人妻,人人澡人人爽秒播| 亚洲一区中文字幕在线| av天堂久久9| av片东京热男人的天堂| 成人亚洲精品一区在线观看| avwww免费| 伊人久久大香线蕉亚洲五| 亚洲 欧美一区二区三区| 日韩有码中文字幕| 大型av网站在线播放| 国产精品久久久人人做人人爽| 国产精品1区2区在线观看. | 国产成人免费无遮挡视频| 最新在线观看一区二区三区| 99国产精品一区二区蜜桃av | 精品熟女少妇八av免费久了| av有码第一页| ponron亚洲| 下体分泌物呈黄色| 黄片播放在线免费| 操出白浆在线播放| 亚洲av熟女| avwww免费| 欧美成人午夜精品| 黄色毛片三级朝国网站| 久久性视频一级片| 久久香蕉国产精品| 91字幕亚洲| 美女国产高潮福利片在线看| 国产精品一区二区免费欧美| 久久午夜亚洲精品久久| 亚洲精品国产精品久久久不卡| 丝瓜视频免费看黄片| 亚洲一区中文字幕在线| 国产在线精品亚洲第一网站| 久久久国产成人免费| 精品福利永久在线观看| 美女高潮喷水抽搐中文字幕| 国产高清激情床上av| 久久久精品区二区三区| 久久久久国产精品人妻aⅴ院 | 男人舔女人的私密视频| 欧美日韩av久久| 在线观看免费高清a一片| 丰满迷人的少妇在线观看| 视频区图区小说| 国产黄色免费在线视频| 精品国产乱码久久久久久男人| 日韩精品免费视频一区二区三区| 999久久久国产精品视频| 色老头精品视频在线观看| 久久午夜综合久久蜜桃| 精品久久蜜臀av无| 一区在线观看完整版| 大香蕉久久成人网| 久久国产精品大桥未久av| 亚洲av美国av| av天堂在线播放| 老司机在亚洲福利影院| 成人特级黄色片久久久久久久| 国产精品久久久久成人av| 丰满人妻熟妇乱又伦精品不卡| 最新在线观看一区二区三区| 精品熟女少妇八av免费久了| 高清黄色对白视频在线免费看| 中文字幕人妻丝袜一区二区| av在线播放免费不卡| 自线自在国产av| 亚洲中文字幕日韩| 午夜福利影视在线免费观看| 免费在线观看日本一区| 国产成人啪精品午夜网站| 女警被强在线播放| 男人操女人黄网站| cao死你这个sao货| 国产精品自产拍在线观看55亚洲 | 精品熟女少妇八av免费久了| 99riav亚洲国产免费| 18禁美女被吸乳视频| 国产三级黄色录像| 亚洲第一青青草原| 国产真人三级小视频在线观看| 国产深夜福利视频在线观看| 中文字幕制服av| 丝袜美足系列| 国产精品永久免费网站| 日本撒尿小便嘘嘘汇集6| 少妇的丰满在线观看| 久久这里只有精品19| 免费不卡黄色视频| 男人操女人黄网站| 久9热在线精品视频| 18禁裸乳无遮挡免费网站照片 | 国产精品久久久久成人av| 少妇猛男粗大的猛烈进出视频| 两个人看的免费小视频| 亚洲精品乱久久久久久| 亚洲人成电影免费在线| 欧美乱妇无乱码| 男人的好看免费观看在线视频 | 黑人欧美特级aaaaaa片| 可以免费在线观看a视频的电影网站| bbb黄色大片| 热99re8久久精品国产| 国产国语露脸激情在线看| 国产欧美日韩综合在线一区二区| 国产欧美日韩一区二区精品| 十八禁人妻一区二区| 国产xxxxx性猛交| 亚洲av熟女| 夫妻午夜视频| 在线观看免费视频日本深夜| 热99久久久久精品小说推荐| 成人国产一区最新在线观看| 99国产极品粉嫩在线观看| 黑人巨大精品欧美一区二区蜜桃| av中文乱码字幕在线| 免费看a级黄色片| videos熟女内射| 亚洲精品在线观看二区| 男人舔女人的私密视频| 美女视频免费永久观看网站| 国产高清视频在线播放一区| 纯流量卡能插随身wifi吗| 12—13女人毛片做爰片一| 精品电影一区二区在线| 婷婷成人精品国产| 精品国产亚洲在线| 女人被躁到高潮嗷嗷叫费观| 动漫黄色视频在线观看| 一区二区三区激情视频| 丝袜在线中文字幕| 日韩免费av在线播放| 中文字幕制服av| 一二三四在线观看免费中文在| 亚洲午夜理论影院| netflix在线观看网站| 午夜免费鲁丝| 国产日韩一区二区三区精品不卡| xxx96com| 怎么达到女性高潮| 国产精品一区二区免费欧美| 真人做人爱边吃奶动态| 在线国产一区二区在线| 久久精品国产综合久久久| 亚洲国产精品合色在线| 欧美日韩瑟瑟在线播放| 亚洲精品在线观看二区| 亚洲国产精品sss在线观看 | 看片在线看免费视频| 一区二区三区国产精品乱码| netflix在线观看网站| 又黄又粗又硬又大视频| 欧美乱色亚洲激情| a级毛片黄视频| 女性被躁到高潮视频| 亚洲欧美精品综合一区二区三区| 亚洲精品中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| 精品人妻1区二区| 精品视频人人做人人爽| 国产精品二区激情视频| 亚洲精品国产精品久久久不卡| 天堂中文最新版在线下载| 免费看十八禁软件| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜精品久久久久久毛片777| 国产亚洲欧美98| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕人妻熟女乱码| 亚洲成a人片在线一区二区| 午夜成年电影在线免费观看| 精品国产乱子伦一区二区三区| 精品卡一卡二卡四卡免费| 看黄色毛片网站| 在线观看免费高清a一片| 下体分泌物呈黄色| 欧美日韩黄片免| 新久久久久国产一级毛片| 久热这里只有精品99| www.999成人在线观看| 女性生殖器流出的白浆| 久久亚洲真实| 欧美日韩视频精品一区| 久久精品国产99精品国产亚洲性色 | 视频在线观看一区二区三区| 大香蕉久久成人网| 女人被狂操c到高潮| 人人澡人人妻人| 亚洲一区二区三区不卡视频| 国产亚洲av高清不卡| 欧美+亚洲+日韩+国产| 在线永久观看黄色视频| videosex国产| 电影成人av| 日本撒尿小便嘘嘘汇集6| 国产高清国产精品国产三级| 国产激情欧美一区二区| 欧美日韩瑟瑟在线播放| 免费在线观看影片大全网站| 日韩成人在线观看一区二区三区| 99热只有精品国产| 免费女性裸体啪啪无遮挡网站| 桃红色精品国产亚洲av| 久久久久视频综合| 亚洲精品av麻豆狂野| 午夜激情av网站| 亚洲成人免费av在线播放| 老汉色av国产亚洲站长工具| 91av网站免费观看| 久99久视频精品免费| 欧美久久黑人一区二区| 精品福利观看| 亚洲 国产 在线| 国产一卡二卡三卡精品| 欧美激情高清一区二区三区| 操出白浆在线播放| 大片电影免费在线观看免费| 国产激情欧美一区二区| 国产视频一区二区在线看| 亚洲美女黄片视频| 捣出白浆h1v1| 麻豆av在线久日| 大片电影免费在线观看免费| 好看av亚洲va欧美ⅴa在| 国产无遮挡羞羞视频在线观看| 日韩成人在线观看一区二区三区| 99国产精品免费福利视频| 18禁美女被吸乳视频| 国产麻豆69| 男女高潮啪啪啪动态图| 欧美 亚洲 国产 日韩一| 乱人伦中国视频| 啦啦啦 在线观看视频| av超薄肉色丝袜交足视频| 校园春色视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产精品免费视频内射| 女人高潮潮喷娇喘18禁视频| 日本精品一区二区三区蜜桃| 老鸭窝网址在线观看| 欧美丝袜亚洲另类 | 久久精品aⅴ一区二区三区四区| 欧美黄色片欧美黄色片| 久久99一区二区三区| x7x7x7水蜜桃| 91大片在线观看| 人妻丰满熟妇av一区二区三区 | 国产精品影院久久| 午夜福利,免费看| 麻豆av在线久日| 国产激情久久老熟女| 欧美乱码精品一区二区三区| 一进一出抽搐gif免费好疼 | 丝瓜视频免费看黄片| 日韩精品免费视频一区二区三区| 午夜福利,免费看| 黄频高清免费视频| 女人被狂操c到高潮| 少妇猛男粗大的猛烈进出视频| 亚洲三区欧美一区| 80岁老熟妇乱子伦牲交| 亚洲人成77777在线视频| 波多野结衣一区麻豆| 制服人妻中文乱码| 久久香蕉精品热| 国产xxxxx性猛交| 亚洲av电影在线进入| 久久久国产成人免费| 免费在线观看日本一区| 亚洲欧美激情在线| 在线观看免费视频网站a站| 亚洲性夜色夜夜综合| 视频区图区小说| 黄色毛片三级朝国网站| 美女视频免费永久观看网站| 宅男免费午夜| 麻豆av在线久日| 搡老乐熟女国产| 在线观看66精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 五月开心婷婷网| 亚洲一卡2卡3卡4卡5卡精品中文| 99精品久久久久人妻精品| 国产精品久久电影中文字幕 | 亚洲av成人不卡在线观看播放网| 99热只有精品国产| 啪啪无遮挡十八禁网站| 91麻豆av在线| 欧美 亚洲 国产 日韩一| 午夜影院日韩av| 狠狠婷婷综合久久久久久88av| 国产国语露脸激情在线看| 看黄色毛片网站| 亚洲美女黄片视频| 久久中文字幕人妻熟女| 女人被狂操c到高潮| 黄色毛片三级朝国网站| 亚洲精品中文字幕在线视频| 免费在线观看完整版高清| 两性夫妻黄色片| 法律面前人人平等表现在哪些方面| 热99国产精品久久久久久7| 国产日韩欧美亚洲二区| 女性生殖器流出的白浆| 中文字幕人妻熟女乱码| 看片在线看免费视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品二区激情视频| 欧美乱妇无乱码| xxxhd国产人妻xxx| 国产无遮挡羞羞视频在线观看| 热99久久久久精品小说推荐| 在线观看午夜福利视频| 两个人看的免费小视频| 中文亚洲av片在线观看爽 | 亚洲免费av在线视频| 黄色毛片三级朝国网站| 男人舔女人的私密视频| 女性生殖器流出的白浆| 99久久精品国产亚洲精品| 国产1区2区3区精品| 丝瓜视频免费看黄片| 国产精品一区二区在线观看99| 日本撒尿小便嘘嘘汇集6| 国产激情欧美一区二区| 久久人人97超碰香蕉20202| 午夜福利在线免费观看网站| 首页视频小说图片口味搜索| 国产精品久久久久成人av| 久久国产精品男人的天堂亚洲| 男女免费视频国产| 成熟少妇高潮喷水视频| 中文字幕精品免费在线观看视频| 欧美av亚洲av综合av国产av| 亚洲国产精品一区二区三区在线| 亚洲avbb在线观看| 无限看片的www在线观看| 国产精品久久久av美女十八| 久久人妻熟女aⅴ| 成人三级做爰电影| 欧美日韩av久久| 80岁老熟妇乱子伦牲交| 国产精华一区二区三区| 精品第一国产精品| 一进一出抽搐gif免费好疼 | 欧美黑人精品巨大| 亚洲五月天丁香| 久久影院123| 中文字幕制服av| 精品久久久久久久毛片微露脸| 亚洲第一欧美日韩一区二区三区| 老司机在亚洲福利影院| 国产欧美日韩一区二区三| 久久久久久久久久久久大奶| av天堂久久9| 不卡一级毛片| 91精品三级在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品高清国产在线一区| 国产精品久久久av美女十八| 最近最新免费中文字幕在线| 淫妇啪啪啪对白视频| 伊人久久大香线蕉亚洲五| 黄色视频,在线免费观看| 精品久久蜜臀av无| 欧美精品亚洲一区二区| 亚洲成av片中文字幕在线观看| 1024视频免费在线观看| 日本撒尿小便嘘嘘汇集6| 丰满人妻熟妇乱又伦精品不卡| 国产成人啪精品午夜网站| 国内久久婷婷六月综合欲色啪| 精品国内亚洲2022精品成人 | 不卡av一区二区三区| 一进一出抽搐gif免费好疼 | 日韩欧美一区视频在线观看| 夜夜爽天天搞| 久久影院123| 国产高清视频在线播放一区| 国产男女超爽视频在线观看| 欧美成人免费av一区二区三区 | 美女午夜性视频免费| 中文欧美无线码| 王馨瑶露胸无遮挡在线观看| 在线视频色国产色| 国产成人影院久久av| 亚洲午夜理论影院| 精品人妻在线不人妻| 久久香蕉激情| 免费日韩欧美在线观看| 丰满的人妻完整版| 久久久精品国产亚洲av高清涩受| 女性被躁到高潮视频| 欧美午夜高清在线| 欧美人与性动交α欧美精品济南到| 女人爽到高潮嗷嗷叫在线视频| av天堂在线播放| 人人妻人人澡人人爽人人夜夜| 免费看十八禁软件| 欧美久久黑人一区二区| 午夜精品国产一区二区电影| 国产成人精品无人区| 国产日韩一区二区三区精品不卡| 99国产极品粉嫩在线观看| 高清毛片免费观看视频网站 | 欧美亚洲 丝袜 人妻 在线| 男男h啪啪无遮挡| 黄色女人牲交| 热99re8久久精品国产| 国产精品99久久99久久久不卡| 成年人免费黄色播放视频| 99久久综合精品五月天人人| 午夜两性在线视频| 叶爱在线成人免费视频播放| av天堂久久9| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 在线国产一区二区在线| 欧美精品高潮呻吟av久久| 麻豆国产av国片精品| av电影中文网址| 亚洲视频免费观看视频| 国产精品美女特级片免费视频播放器 | 久久天躁狠狠躁夜夜2o2o| 男男h啪啪无遮挡| 亚洲片人在线观看| 美女 人体艺术 gogo| 美女国产高潮福利片在线看| 黄色视频不卡| 久久久久精品国产欧美久久久| www.自偷自拍.com| 国产成人影院久久av| 欧美乱色亚洲激情| 成年人免费黄色播放视频| 国产亚洲欧美精品永久| 亚洲精品在线观看二区| 亚洲五月婷婷丁香| 欧美日韩瑟瑟在线播放| 午夜福利在线免费观看网站| 欧美丝袜亚洲另类 | 人成视频在线观看免费观看| 欧美最黄视频在线播放免费 | 超碰成人久久| 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| 国产在线一区二区三区精| 在线观看午夜福利视频| 不卡av一区二区三区| 亚洲av日韩精品久久久久久密| 中文字幕最新亚洲高清| 日韩 欧美 亚洲 中文字幕| 久久久精品国产亚洲av高清涩受| 一级a爱视频在线免费观看| 香蕉丝袜av| 黄色毛片三级朝国网站| a级片在线免费高清观看视频| 国产97色在线日韩免费| 少妇粗大呻吟视频| 日韩大码丰满熟妇| 久久久久国产精品人妻aⅴ院 | 久久久久久久久免费视频了| 亚洲欧美色中文字幕在线| 久久久精品国产亚洲av高清涩受| 天堂俺去俺来也www色官网| 99国产精品一区二区蜜桃av | 黄色视频不卡| www.熟女人妻精品国产| 国产又爽黄色视频| 18禁黄网站禁片午夜丰满| 日本黄色视频三级网站网址 | 无限看片的www在线观看| av有码第一页| 50天的宝宝边吃奶边哭怎么回事| 黄色女人牲交| av网站在线播放免费| 国产野战对白在线观看| 操出白浆在线播放| 精品一区二区三区视频在线观看免费 | netflix在线观看网站| 国产一区二区三区综合在线观看| 丰满迷人的少妇在线观看| 久久狼人影院| 久久久久视频综合| 十八禁高潮呻吟视频| 国产欧美日韩综合在线一区二区| 日韩中文字幕欧美一区二区| 亚洲av熟女| 亚洲va日本ⅴa欧美va伊人久久| 国产男女超爽视频在线观看| 一级毛片高清免费大全| 国产一区二区三区综合在线观看| 免费av中文字幕在线| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 十八禁高潮呻吟视频| 国产国语露脸激情在线看| 老司机靠b影院| 女警被强在线播放| 好男人电影高清在线观看| 午夜精品久久久久久毛片777| 麻豆乱淫一区二区| 亚洲国产欧美一区二区综合| 久久久久国产精品人妻aⅴ院 | 50天的宝宝边吃奶边哭怎么回事| 国产不卡av网站在线观看| 多毛熟女@视频| 成年动漫av网址| 老司机午夜福利在线观看视频| 嫩草影视91久久| 国产欧美日韩综合在线一区二区| 在线观看www视频免费| 99riav亚洲国产免费| 欧美精品亚洲一区二区| av网站免费在线观看视频| 欧美激情久久久久久爽电影 | 久久精品亚洲精品国产色婷小说| 精品国产一区二区三区久久久樱花| 丝袜美足系列| 亚洲人成77777在线视频| 99国产精品99久久久久| 亚洲av日韩精品久久久久久密| 日日夜夜操网爽| 久久久水蜜桃国产精品网| 国产精品电影一区二区三区 | 亚洲欧洲精品一区二区精品久久久| 久久亚洲精品不卡| 夫妻午夜视频| 日韩欧美一区二区三区在线观看 | 成人永久免费在线观看视频| 国内毛片毛片毛片毛片毛片| 免费在线观看视频国产中文字幕亚洲| 91av网站免费观看| 欧美 亚洲 国产 日韩一| 99国产精品一区二区蜜桃av | 久久精品人人爽人人爽视色| 真人做人爱边吃奶动态| 国产精品自产拍在线观看55亚洲 | 女性生殖器流出的白浆| 看黄色毛片网站| 99精国产麻豆久久婷婷| 中文字幕精品免费在线观看视频| 久久中文字幕人妻熟女| av电影中文网址| 超色免费av| 欧美精品亚洲一区二区| 午夜福利在线观看吧| 成人特级黄色片久久久久久久| 一区二区三区精品91| 丝袜在线中文字幕| 日韩有码中文字幕| 日本五十路高清| 少妇的丰满在线观看| 国产不卡一卡二| 黄色丝袜av网址大全| 亚洲男人天堂网一区| 亚洲 欧美一区二区三区| 午夜福利在线免费观看网站| 无限看片的www在线观看| 大型av网站在线播放| 亚洲欧美日韩另类电影网站|