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

    葉綠體序列和核微衛(wèi)星顯示同域分布潤楠屬物種雜交嚴(yán)重

    2025-03-29 00:00:00蔣凱王正偉陳小勇
    廣西植物 2025年3期

    摘 要:" 樟科潤楠屬(Machilus Nees)植物為常綠喬木,是亞熱帶常綠闊葉林常見的優(yōu)勢種,具有重要的經(jīng)濟價值,但屬內(nèi)分類和種間系統(tǒng)進化關(guān)系存在較大爭議。為了解其可能的原因,該文選擇12個同域分布的潤楠屬物種,采用72個葉綠體DNA序列片段和10個核微衛(wèi)星位點進行分析。結(jié)果表明:(1) 72個葉綠體DNA序列片段共44 025 bp的序列上僅有48個變異位點,多態(tài)性極低。(2) 10個核微衛(wèi)星標(biāo)記在12個潤楠屬物種中通用性好、多態(tài)性高,不同物種間共享大量等位基因,獨有等位基因少。(3) 基于葉綠體DNA不能得到支持率高的一致系統(tǒng)進化樹,微衛(wèi)星標(biāo)記支持建潤楠(M. oreophila)、狹葉潤楠(M. rehderi)和紅楠(M. thunbergii)的分類地位,但其他物種間關(guān)系混亂且核-質(zhì)系統(tǒng)進化關(guān)系沖突,其原因可能是種間頻繁的雜交漸滲模糊了種間的界限。后續(xù)需綜合形態(tài)學(xué)和多種分子標(biāo)記,開發(fā)有效的分析工具,以理解潤楠屬內(nèi)復(fù)雜的系統(tǒng)進化關(guān)系。

    關(guān)鍵詞: 潤楠屬, 葉綠體序列, 核微衛(wèi)星, 系統(tǒng)進化, 漸滲雜交

    中圖分類號:" Q949

    文獻標(biāo)識碼:" A

    文章編號:" 1000-3142(2025)03-0555-12

    Chloroplast sequences and nuclear microsatellites reveal extensive hybridization among sympatric species of the genus Machilus

    JIANG Kai1,2, WANG Zhengwei1, CHEN Xiaoyong2*

    ( 1. Eastern China Conservation Centre for Wild Endangered Plant Resources, Shanghai Chenshan Botanical Garden, Shanghai 201602,China;

    2. School of Ecological and Environmental Sciences, East China Normal University, Shanghai 200062, China )

    Abstract:" Plants of the genus Machilus Nees (Lauraceae) are evergreen trees and are common dominant species in subtropical evergreen broad-leaved forests, with significant economic value. However, there is considerable debate regarding the classification within the genus and the phylogenetic relationships among species. To explore the potential reasons behind this, this study selected 12 sympatric species of Machilus and analyzed 72 chloroplast fragments and 10 nuclear microsatellite loci. The results were as follows: (1) A total of only 48 variable sites were found across the 44 025 base pairs of the 72 chloroplast fragments, indicating extremely low polymorphism. (2) The 10 nuclear microsatellite markers showed good universality across the 12 species of Machilus, with high polymorphism, a large number of shared alleles among different species, and very few private alleles. (3) A high-support consensus phylogenetic tree could not be obtained based on chloroplast DNA, while microsatellite markers supported the taxonomic status of M. oreophila, M. rehderi, and M. thunbergii. However, those relationships among other species were confused and there was incongruence between nuclear and plastid phylogenies. This confusion is likely attributed to frequent interspecific hybridization" that have blurred the boundaries between species. Further research is needed to integrate morphological data and multiple molecular markers, and to develop effective analytical tools to understand the complex phylogenetic relationships within the genus Machilus.

    Key words: Machilus, chloroplast sequence, nuclear microsatellites, phylogenetics, introgressive hybridization

    潤楠屬(Machilus Nees)隸屬于樟科(Lauraceae),全球已知100余種,廣泛分布于中國的大陸和臺灣島,以及朝鮮半島、日本群島、中南半島、印度次大陸、馬來西亞、印度尼西亞和菲律賓等地。我國有82種,其中63種為我國特有(Wei et al., 2008)。大部分潤楠屬植物為東亞熱帶、亞熱帶常綠闊葉森林的建群種,是研究生物地理的重要類群,對于理解亞熱帶常綠闊葉林的起源和演化以及生物多樣性保護具有重要意義(林夏珍,2007a, b; Zhu et al., 2017; Jiang et al., 2021; Fan et al., 2022; Hu et al., 2022; Xiao et al., 2022)。此外,潤楠屬植物還具有重要的經(jīng)濟價值和觀賞價值,為江南四大名木——楠木的主要來源。屬內(nèi)不少物種為優(yōu)良園林觀賞植物,其樹皮和葉可制作褐色染料以及各種熏香的調(diào)和劑或飲用水的純化劑,全株可提取芳香油,種子油可作潤滑劑(Cheng et al., 2009; Yang SS et al., 2022; Yang SY et al., 2022; Yang et al., 2023; To et al., 2023; Huong et al., 2024)。

    潤楠屬植物的準(zhǔn)確分類和鑒定是保護和可持續(xù)利用潤楠屬資源的基礎(chǔ)。形態(tài)學(xué)和分子標(biāo)記證據(jù)已確立潤楠屬的屬級地位(李樹剛,1982; 王玉國和韋發(fā)南,2003; 韋發(fā)南和唐賽春,2006; Rohwer et al., 2009)。經(jīng)典形態(tài)分類系統(tǒng)依據(jù)花被片外面有無毛以及果實形狀和大小將該屬分為6個大組(李樹剛,1982)?;诖耍顦鋭偤晚f發(fā)南(1984)對我國潤楠屬進行了修訂,后韋發(fā)南等(2001)對一些潤楠屬物種進行了合并,如嘉道理楠(Persea kadooriei)=信宜潤楠(Machilus wangchiana)、長序潤楠(M. longipedunculata) = 浙江潤楠(M. chekiangensis)、密脈潤楠(M. polyneura)=刨花潤楠(M. pauhoi)等。然而,潤楠屬內(nèi)雖有一些形態(tài)特征十分相似但不易區(qū)分的物種,如華潤楠(M. chinensis)與紅楠(M. thunbergii)、絨毛潤楠(M. velutina)與黃絨潤楠(M. grijsii)等種類的葉形比較相近,常發(fā)生誤定現(xiàn)象(林夏珍,2007a,b)。因此,莊雪影(1997)、莊雪影(2002)、林夏珍(2007a,b)以及鐘義和夏念和(2010)對該屬部分物種進行孢粉學(xué)、胚胎學(xué)、葉表皮和木材結(jié)構(gòu)的比較,嘗試為潤楠屬屬下分類尋找輔助特征。然而,組內(nèi)不同物種間花粉特征存在交叉,屬內(nèi)一些組的劃分也是武斷的(湯庚國和向其柏,1995)。此外,這些研究所選材料不一致,很難將其結(jié)果進行統(tǒng)一分析。

    近幾十年來,分子標(biāo)記廣泛應(yīng)用于物種界定和系統(tǒng)進化分析,并在潤楠屬中進行了探討。例如,陳俊秋等(2009)運用ITS序列、基于貝葉斯進化樹的分析結(jié)果表明,潤楠屬組成鱷梨屬群中一獨立分支是一個自然的單系類群,但潤楠屬內(nèi)物種分支混亂,與傳統(tǒng)植物學(xué)分類中對潤楠屬內(nèi)各組的劃分相矛盾。因此,不支持前人依據(jù)花被片外面被毛情況及果實大小所建立的潤楠屬屬下系統(tǒng)。Li等(2011)采用ITS和LEAFY序列分析鱷梨屬群的系統(tǒng)進化關(guān)系,得到類似結(jié)論。Song等(2017)基于樟科27個代表類群共47個葉綠體基因組序列構(gòu)建系統(tǒng)樹,確定了樟科不同屬間的進化關(guān)系,進一步支持潤楠屬為單系類群。Xiao等(2022)使用鱷梨族(Perseeae)53個物種(包括25個潤楠屬物種)的葉綠體基因組構(gòu)建了屬級水平的系統(tǒng)進化樹,并確認潤楠屬為單系,但屬內(nèi)不同物種間分支在系統(tǒng)進化樹上的支持率仍不高。這些研究結(jié)果顯示,潤楠屬種間關(guān)系不明晰,其屬下物種系統(tǒng)進化關(guān)系復(fù)雜,還需要深入探討這一現(xiàn)象產(chǎn)生的原因。

    造成潤楠屬系統(tǒng)進化關(guān)系復(fù)雜、物種難以界定的原因:可能是該屬起源較晚,物種間尚未形成明顯的系統(tǒng)進化關(guān)系;也可能是因這個屬種間雜交非常嚴(yán)重而導(dǎo)致種間形態(tài)特征相似、遺傳關(guān)系不穩(wěn)(Rohwer et al., 2009)。因此,本研究使用72個葉綠體DNA序列片段和10對基于紅楠基因組開發(fā)出的核微衛(wèi)星位點分析同域分布的潤楠屬物種,以期了解這些物種間的進化關(guān)系,初步分析不同標(biāo)記在潤楠屬分類中的適用性,并探討造成潤楠屬種間系統(tǒng)分類困難的可能原因。

    1 材料與方法

    1.1 研究地點和樣品采集

    我們主要以同域分布的潤楠屬物種作為研究對象,從分布于廣西壯族自治區(qū)、廣東省和湖南省的潤楠屬植物中選取了12個分布區(qū)重疊的物種用于本研究。不同物種之間的地理位置距離5~450 km不等(圖1)。這12個物種分別是華潤楠、基脈潤楠(M. decursinervis)、建潤楠(M. oreophila)、狹葉潤楠(M. rehderi)、瓊桂潤楠(M. foonchewii)、短序潤楠(M. breviflora)、粗壯潤楠(M. robusta)、薄葉潤楠(M. leptophylla)、黃絨潤楠、黔桂潤楠(M. chienkweiensis)、紅楠和刨花潤楠。它們分屬于潤楠屬5個組,即毛花組6個物種、光花組3個物種、大果組1個物種、絨毛組和滇黔桂組各1個物種。每個物種至少采集1份植物標(biāo)本,樣品用硅膠干燥處理,并保存于上海辰山植物標(biāo)本館中。

    1.2 同域分布潤楠屬物種的潛在分布區(qū)模擬

    為預(yù)測12個潤楠屬植物的潛在分布區(qū),使用MaxEnt v.3.3.1軟件(Phillips amp; Dudík, 2008)中的最大熵方法來模擬當(dāng)前的分布。從科學(xué)引文索引數(shù)據(jù)庫、中國國家知識基礎(chǔ)設(shè)施、全球生物多樣性信息設(shè)施和中國虛擬植物標(biāo)本館進行搜索,獲得了12個潤楠屬物種的分布記錄。排除明顯超出自然分布范圍的記錄,并移除同一個物種中相近(48 km以內(nèi))的記錄點。從WorldClim數(shù)據(jù)庫(http://www.worldclim.org/)中提取當(dāng)前19個生物氣候因子數(shù)據(jù)。使用ENMtools v 1.3軟件 (Warren et al., 2010)排除高度相關(guān)的生物氣候變量(rgt;0.7),最終保留了4個氣候變量(Bio1、Bio3、Bio16和Bio19)。我們進行20次獨立的模型交叉驗證,其中75%的數(shù)據(jù)用于訓(xùn)練,剩余25%用于測試。采用特征曲線下面積(area under the curve,AUC)評估模型性能(Elith et al., 2011)。

    1.3 DNA提取和PCR擴增

    稱取硅膠干燥的葉片樣品約35 mg,使用天根植物基因組提取試劑盒[DP305,天根生化科技(北京)有限公司]提取基因組總DNA。潤楠屬植物多糖多酚含量較高,采用氯仿抽提2次以去除雜質(zhì)和提高DNA質(zhì)量。于4 ℃冰箱靜置1 d降解RNA,吸取3 μL DNA原液,用0.8%瓊脂糖凝膠電泳進行質(zhì)量檢測。根據(jù)紅楠葉綠體基因組的長單拷貝區(qū)序列,設(shè)計了72對引物,用于潤楠屬12個物種的葉綠體DNA擴增和測序。使用Mastercycler nexus PCR擴增儀(Eppendorf, Hamburg, Germany)進行PCR擴增。將擴增條帶清晰且長度與設(shè)計的目標(biāo)片段大小相同的PCR擴增產(chǎn)物,送生工生物工程(上海)股份有限公司進行雙向測序。

    根據(jù)發(fā)表的紅楠微衛(wèi)星引物(Kaneko et al., 2012),隨機選取18對微衛(wèi)星引物序列,交由生工生物工程(上海)股份有限公司合成。每對引物使用2個不同紅楠種群個體的基因組DNA進行溫度梯度擴增,以優(yōu)化最適退火溫度。PCR產(chǎn)物經(jīng)1.5%瓊脂糖凝膠電泳,在Gel Doc2000TM凝膠成像系統(tǒng)(Bio-RAD, California, USA)上拍照記錄。通過溫度梯度的方法,我們從18對引物中共篩選出10對條帶清晰且穩(wěn)定性好的引物用于后續(xù)研究(Jiang et al., 2021)。使用這10對引物對潤楠屬12個物種的種群樣品進行擴增,每個種群擴增11~28個樣品。擴增產(chǎn)物送生工生物工程(上海)股份有限公司用Applied BiosystemsTM 3730XL儀器(Thermo Fisher, the USA)進行片段掃描,獲得樣品的原始片段長度數(shù)據(jù)。隨后,我們使用GeneMarker軟件(Holland amp; Parson, 2011)對片段長度數(shù)據(jù)進行讀取和判定,確定個體的基因型。

    1.4 葉綠體基因組序列數(shù)據(jù)分析

    1.4.1 序列拼接及校正 使用Geneious R11軟件(https://www.geneious.com)讀取測序獲得的原始序列峰圖,并對序列進行拼接、 手工校正序列存在的沖突位點,最終獲得一致序列。使用MAFFT v.6.240軟件(Katoh et al., 2019),對12個物種的72個葉綠體DNA序列片段分別進行對齊。其中,59個片段完全一致,13個片段存在變異。因此,后續(xù)主要使用這13個葉綠體DNA序列片段進行分析。本研究將12個物種的13個多態(tài)葉綠體DNA序列片段首尾拼接后,獲得序列矩陣。同時,從NCBI數(shù)據(jù)庫中下載2個潤楠屬植物 [滇潤楠(Machilus yunnanensis)]和銹毛潤楠(M. balansae)和5個楠屬植物 [閩楠(Phoebe bournei)、楨楠(P. zhennan)、浙江楠(P. chekiangensis)、竹葉楠(P. omeiensis)和紫楠(P. sheareri)]的葉綠體基因組序列,并從這些基因組上提取上述13個葉綠體DNA序列片段所在區(qū)域。最終將這些序列進行首尾對齊,相連后組成一個包含19個物種10 756 bp長度的葉綠體DNA序列片段組成的序列矩陣,用于后續(xù)系統(tǒng)進化分析。

    1.4.2 葉綠體DNA序列突變網(wǎng)絡(luò)分析 使用Network v.5.0軟件 (Bandelt et al., 1999),分析葉綠體DNA序列的突變網(wǎng)絡(luò)。對于每個潤楠屬物種,將13條具多態(tài)性的葉綠體DNA片段首尾拼接成1條序列,每條序列表示1個單倍型。使用Median-Joining網(wǎng)絡(luò)算法(Weight = 10, Epsilon = 0)對這14個單倍型的突變狀態(tài)進行計算,改變Epsilon參數(shù)為10、20、30分別計算單倍型間的關(guān)系,最終獲得單倍型間的突變網(wǎng)絡(luò)圖。

    1.4.3 系統(tǒng)進化關(guān)系分析 使用鄰接法(neighbor joining method, NJ)和非加權(quán)配對算術(shù)平均法(unweighted paired group method using arithmetic average, UPGMA)構(gòu)建物種的系統(tǒng)進化樹,自展重復(fù)(bootstrap)1 000次。相關(guān)分析在MEGA v.10.0.5軟件中完成(Kumar et al., 2018)。

    使用MEGA軟件進行最大簡約法(maximum parsimony,MP)分析。基于啟發(fā)式搜索,隨機逐步添加物種序列,并將這些序列逐個添加到能獲得最優(yōu)樹的位置,重復(fù)這一步驟直到所有序列全部添加進來建立初始樹。在初始樹的基礎(chǔ)上使用二分樹再接法(tree bisection reconnection, TBR)中的分支交換法尋找最優(yōu)樹。使用1 000次自展重復(fù)抽樣來評估系統(tǒng)進化樹每支的支持率。

    根據(jù)前面模型選擇結(jié)果確定適合本研究數(shù)據(jù)的最優(yōu)進化模型,以最大簡約樹為初始樹采用最大似然法(maximum likelihood,ML)估計模型中的參數(shù)值,并尋找似然值最大時的系統(tǒng)進化樹,使用RAxML v.8.0軟件 (Stamatakis, 2014)基于最適模型GTR+G+I進行系統(tǒng)進化樹構(gòu)建,使用1 000次自展重復(fù)抽樣評估系統(tǒng)進化樹上每支的支持率。

    使用MrBayes v.3.2軟件(Ronquist et al., 2012)基于貝葉斯方法構(gòu)建系統(tǒng)進化樹?;隈R爾科夫鏈蒙特卡洛法(Markov chain Monte Carlo approach, MCMC)進行107代運算,每次運行均可得到1棵系統(tǒng)樹,每隔1 000代取樣1次。本研究將運行的前25%馬爾科夫鏈作為老化值(burn-in)丟棄,最終使用剩下的樹進行一致樹分析。使用Tracer v.1.7軟件 (Rambaut et al., 2018)對似然值分布和有效樣本大小進行分析,以判斷各參數(shù)值是否收斂。另外,通過判斷二分樹頻率的平均標(biāo)準(zhǔn)差是否小于或等于0.01來確定運行是否達到穩(wěn)態(tài)。最終以50%多數(shù)規(guī)則合一樹作為最優(yōu)樹。

    1.4.4 物種分化時間推斷 選取楠屬5個物種(閩楠、楨楠、浙江楠、竹葉楠和紫楠)作為外類群。使用BEAST v2.5軟件(Bouckaert et al., 2019)對物種分化時間進行估計。使用BEAUti軟件創(chuàng)建BEAST的輸入文件。使用不相關(guān)對數(shù)正態(tài)分布的松散分子鐘(uncorrelated lognormal relaxed clock)模型,核苷酸的替代模型選擇GTR,選用Yule模型。到目前為止,潤楠屬中尚未發(fā)現(xiàn)可信的化石。因此,根據(jù)樟科系統(tǒng)進化樹定年結(jié)果,將潤楠屬和楠屬的分化時間設(shè)置為23 Ma,使用正態(tài)分布,其他參數(shù)采用默認設(shè)置。運行1×107次馬爾科夫鏈,每隔1 000次進行1次取樣,并將最初的1×106次作為老化值舍棄。使用Tracer v.1.7軟件 (Rambaut et al., 2018)對運行結(jié)果進行分析,并判斷各參數(shù)計算結(jié)果的收斂性。對獲得的后綴為.tree的文件使用TreeAnnotator v.2.4.7軟件進行注釋,使用Figtree軟件對獲得注釋的系統(tǒng)樹文件進行顯示。

    1.5 微衛(wèi)星數(shù)據(jù)分析

    1.5.1 異常位點、連鎖不平衡檢驗和啞等位基因檢測 使用Lositan軟件 (Antao et al., 2008)進行中性檢驗。使用FSTAT v2.9.3.2軟件 (http://www.unil.ch/popgen/softwares/fstat.htm)對篩選出的10個位點進行連鎖不平衡的顯著性檢驗。

    由于本研究所使用的微衛(wèi)星引物是針對紅楠開發(fā),采用這些引物對同屬其他物種進行研究時,可能存在啞等位基因的情況。因此, 在進行數(shù)據(jù)分析之前我們使用Microchecker軟件 (Van Oosterhout et al., 2004)對每個種群的啞等位基因進行檢測。

    1.5.2 種群遺傳多樣性分析 計算平均等位基因數(shù)(mean number of alleles, A)和等位基因豐富度(allele richness, AR)來衡量種群遺傳多樣性。使用TFPGA v1.3軟件(http://www.marksgenetic-software.net/)計算期望雜合度(unbiased heterozyousity, HE)和觀察雜合度(observed heterozyousity, HO)。在FSTAT v2.9.3.2軟件中經(jīng)過1 000次模擬來檢驗總的近交系數(shù)(inbreeding coefficient, FIS)是否顯著偏離0。

    1.5.3 遺傳聚類分析 采用STRUCTURE v2.3 (Pritchard et al., 2000)對12個物種進行聚類分析,推斷種群遺傳結(jié)構(gòu)及每個個體在種群中分布的概率。在軟件中設(shè)置K=1到K=15,基于以上設(shè)置的K值對微衛(wèi)星數(shù)據(jù)分別進行迭代周期(length of period)為1×105、MCMC為1×106次的運算,計算時選擇混合模型(admixture model)。使用Evanno等(2005)的方法得出最優(yōu)物種數(shù)K,使用最優(yōu)K值對分析的12個物種進行分類。本文使用基于DA遺傳距離的鄰接樹(NJ)法獲得每個物種間遺傳距離之間的關(guān)系。相關(guān)分析使用Population 1.2.30軟件 (https://bioinformatics.org/populations/)計算。每個節(jié)點使用基于1 000次重復(fù)的自展概率值進行估算。

    2 結(jié)果與分析

    2.1 葉綠體基因組序列分析

    2.1.1 葉綠體基因組基本信息 使用72對引物對潤楠屬12個物種的長單拷貝區(qū)域進行擴增測序,最終獲得長44 025 bp的序列。12個物種共44 025 bp長度的序列上只有48個變異位點(占總序列長度的0.11%),其中單變異位點為26個(占總序列長度的0.06%)、簡約信息位點為22個(占0.05%)。在72個葉綠體DNA序列片段中,只有13個片段具多態(tài)性。因此,后續(xù)使用這13個葉綠體DNA序列片段組成長度為10 756 bp的序列矩陣進行系統(tǒng)進化分析。對12個潤楠屬物種13個序列的異質(zhì)性檢驗表明,序列間替換型式異質(zhì)性在多數(shù)物種對間不顯著,不會影響系統(tǒng)進化分析。對已測序的紅楠、刨花潤楠、銹毛潤楠和滇潤楠的葉綠體基因組比較發(fā)現(xiàn),葉綠體基因組之間序列的相似度為99.9%,遠高于其他屬內(nèi)物種間葉綠體基因組的相似度。

    2.1.2 潤楠屬系統(tǒng)進化關(guān)系 由于物種間很小的序列差異,因此基于多態(tài)葉綠體DNA片段采用不同的構(gòu)建進化樹方法未得到具一致拓撲結(jié)構(gòu)的系統(tǒng)進化樹(圖2)。UPGMA樹顯示紅楠、建潤楠和狹葉潤楠聚為一支(支持率為0.99),其他的潤楠屬物種聚為一支(支持率為0.97),除這兩大支以外,其他支系的自展支持率都很低。ML樹表明建潤楠、紅楠和狹葉潤楠聚為一支,剩下的物種聚為一支,但自展支持率都很低。NJ樹表明狹葉潤楠、建潤楠和紅楠單獨成為一支(支持率為1.00),其他的9個物種聚為一支且無明顯的進化關(guān)系,支持率很低。這兩種基于距離矩陣方法給出的結(jié)果基本一致,紅楠、建潤楠和狹葉潤楠聚為一支的支持率都很高,其他物種沒有明顯進化關(guān)系(圖2)。MP樹顯示建潤楠、紅楠和狹葉潤楠聚為一支,剩下的其他物種聚為一支。

    2.2 微衛(wèi)星位點分析

    2.2.1 中性位點及啞等位基因分析 使用10個微衛(wèi)星位點對234個個體進行PCR擴增,均獲得清晰的產(chǎn)物。對異常位點的檢測發(fā)現(xiàn),有4個位點分布于中性區(qū)間(95%置信區(qū)間)之外。這些位點的遺傳分化系數(shù)FST大于中性標(biāo)記模擬的FST值的0.975分位數(shù),表明這些位點可能受到歧化選擇的作用。

    10個位點兩兩之間均不存在顯著的連鎖關(guān)系,表明它們?yōu)楠毩⑦z傳的。對啞等位基因的檢測表明,10個微衛(wèi)星位點中有3個可能存在啞等位基因,對這3個可能存在啞等位基因的位點的等位基因頻率校正后用于后續(xù)遺傳多樣性的計算。

    2.2.2 遺傳多樣性 10個微衛(wèi)星位點在12個物種中共檢測到179個等位基因,平均每位點等位基因數(shù)為7.8,建潤楠種群最少(平均為5.0),紅楠種群最多(10.4)。物種的獨有等位基因數(shù)都很少,其中紅楠種群的獨有等位基因最多(平均為每位點0.7個),基脈潤楠和建潤楠中沒有獨有等位基因。

    12個物種中觀察雜合度和期望雜合度中等,期望雜合度為0.501(粗壯潤楠)~0.838 (瓊桂潤楠),觀察雜合度為0.442 (粗壯潤楠)~0.718 (紅楠) (表1)。近交系數(shù)都為正值,最大的為0.277 (刨花潤楠),最小的為0.103(薄葉潤楠),表明這些物種存在較高的近交(表1)。

    所有物種總的遺傳分化系數(shù)FST為0.156。兩兩間的遺傳分化系數(shù)范圍為0.006~0.304,其中華潤楠和瓊桂潤楠間最小,粗壯潤楠和短序潤楠間最大。

    2.3 不同分子標(biāo)記的物種進化關(guān)系及其比較

    通過葉綠體基因組序列分析發(fā)現(xiàn),在這14個潤楠屬物種組成的突變網(wǎng)絡(luò)圖中,共有5個推斷的共同祖先節(jié)點,在圖中分別為V1、V2、V3、V4和V5 (圖3:A)。銹毛潤楠和短序潤楠具有共同祖先V1,黃絨潤楠和基脈潤楠具有共同祖先V2,滇潤楠、黔桂潤楠、薄葉潤楠、粗壯潤楠和瓊桂潤楠具有共同祖先V3,刨花潤楠和華潤楠具有共同祖先V4,建潤楠、紅楠和狹葉潤楠具有共同祖先V5。從V1到V2經(jīng)過2步突變,從V2到V3以及從V3到V4都只經(jīng)歷了1步突變,而從V2到V5經(jīng)歷了14步突變。

    貝葉斯系統(tǒng)樹顯示短序潤楠單獨成為一支,紅楠、建潤楠和狹葉潤楠聚為一支,黃絨潤楠和基脈潤楠聚為一支,華潤楠和刨花潤楠聚為一支,剩下的5個物種聚為一支,但這些系統(tǒng)樹節(jié)點的后驗概率值都不高?;谡量葡到y(tǒng)進化樹確定的潤楠屬和楠屬分化時間為2 300萬年前(23 Ma),潤楠屬的樹冠時間(crown age)為10.75 Ma(圖3:B)。根據(jù)這2個估計時間對潤楠屬系統(tǒng)樹定年,結(jié)果表明V1形成時間大約為7.69 Ma,V2形成的時間大約為4.65 Ma, V3形成的時間大約為3.99 Ma, V4形成的時間大約為1.61 Ma,V5的形成時間大約為2.58 Ma。

    使用10個中性微衛(wèi)星位點對12個潤楠屬物種的個體進行聚類, 最佳K值為2、5和8 (ΔK值最大時的K值)。在這3種情景下,聚類結(jié)果與基于葉綠體序列得到的系統(tǒng)進化關(guān)系都存在明顯的不一致,即核-質(zhì)沖突。當(dāng)K值等于2時,華潤楠、瓊桂潤楠、短序潤楠、狹葉潤楠和紅楠聚為一支,基脈潤楠、黃絨潤楠、粗壯潤楠和薄葉潤楠聚為一支,而黔桂潤楠和刨花潤楠表現(xiàn)出兩支的混雜;當(dāng)K=5時,華潤楠、瓊桂潤楠、短序潤楠和刨花潤楠聚為一支,基脈潤楠、黃絨潤楠、粗壯潤楠和薄葉潤楠聚為一支,而狹葉潤楠、建潤楠和紅楠各成為一支;當(dāng)K=8時,狹葉潤楠、建潤楠和紅楠能夠很好地各成一支,而其他幾個物種混雜比較嚴(yán)重(圖3:C)。

    基于遺傳距離的最近鄰接法構(gòu)建的無根樹顯示,華潤楠、瓊桂潤楠和狹葉潤楠聚成一個支系,基脈潤楠、黃絨潤楠、粗壯潤楠和薄葉潤楠聚為一個支系,而其他的物種介于這兩個支系之間(圖3:D)。最近鄰體法計算出的各個物種之間的關(guān)系與Structure聚類的結(jié)果比較相近,華潤楠、瓊桂潤楠和狹葉潤楠親緣關(guān)系比較近,基脈潤楠、黃絨潤楠、粗壯潤楠和薄葉潤楠親緣關(guān)系比較近。

    3 討論與結(jié)論

    本研究發(fā)現(xiàn)潤楠屬下5組12個物種的葉綠體基因序列突變率較低,據(jù)此得到的系統(tǒng)進化關(guān)系不明晰,與前人的研究結(jié)果類似(Rohwer et al., 2009; Li et al., 2011; Liu et al., 2022)。采用突變率較高的10個潤楠屬內(nèi)通用的核微衛(wèi)星位點,可以清晰地將12個潤楠屬物種分為5組,并且建潤楠、狹葉潤楠和紅楠可單獨分出來,其他物種分為2組?;诤宋⑿l(wèi)星數(shù)據(jù)發(fā)現(xiàn)12個潤楠屬物種共享等位基因占比高,不同物種獨有等位基因占比很低。此外,根據(jù)核微衛(wèi)星聚類結(jié)果與基于cpDNA序列的系統(tǒng)進化關(guān)系存在明顯的不一致。這些信息表明在本研究中同域分布的潤楠屬物種間存在較強的雜交漸滲。

    潤楠屬葉綠體基因組多態(tài)性極低。本研究分析的潤楠屬葉綠體基因組長單拷貝區(qū)約4.4萬個bp,但這12個物種包含的變異位點僅48個,僅占整個葉綠體基因組序列的0.1%左右,與前人的研究結(jié)果基本一致(Song et al., 2017; Liu et al., 2021, 2022)。在其他被子植物中,也有學(xué)者提出以往研究可能高估了葉綠體基因組的多態(tài)性(Smith, 2015)。Liu等 (2022)通過對樟科25個屬131個物種的191個葉綠體基因組序列進行分析,指出盡管葉綠體基因組序列可提高樟科系統(tǒng)進化關(guān)系的分辨率,并可為一些物種分子鑒定提供信息,但增加葉綠體基因組序列只部分提高了物種鑒別能力,仍需要大規(guī)模的核基因數(shù)據(jù)用于確定種間進化關(guān)系。

    使用不同系統(tǒng)進化樹構(gòu)建方法,并基于葉綠體基因得到的潤楠屬種間關(guān)系存在一定差異。上述方法都能將建潤楠、狹葉潤楠和紅楠這3個物種與其他物種分開。從單倍型網(wǎng)絡(luò)圖可以看出,這3個物種之間的差異相對于它們的共同祖先而言有2步(建潤楠)、4步(狹葉潤楠)和5步(紅楠)突變。但是,其他9個物種間的關(guān)系非常混亂,不同數(shù)據(jù)及不同方法得到的結(jié)果差異很大。這些不一致的結(jié)果主要是因葉綠體分辨率較低而造成的。樟科植物葉綠體片段的分辨率都較低(Liu et al., 2022),即使采用更多的葉綠體基因組片段和物種,最終也難以得到支持率較高的系統(tǒng)進化關(guān)系(Liu et al., 2021)。因此,僅根據(jù)葉綠體基因組序列很難構(gòu)建可信的潤楠屬種間系統(tǒng)進化樹。

    采用突變率較高的核微衛(wèi)星位點,可以提升潤楠屬物種分辨率。首先,基于核微衛(wèi)星分析的結(jié)果表明建潤楠、短序潤楠和紅楠可單獨分出來,這與基于葉綠體序列構(gòu)建系統(tǒng)進化樹得到的結(jié)果基本一致,表明這3個物種之間的進化關(guān)系在葉綠體與核基因中都有較高的支持率。然后,與葉綠體序列得到的結(jié)果相比,核微衛(wèi)星可進一步將華潤楠、瓊桂潤楠、狹葉潤楠和刨花潤楠分為1組,將黔桂潤楠、基脈潤楠、黃絨潤楠、粗壯潤楠和薄葉潤楠分為另外1組。最后,當(dāng)進一步劃分時,上述2組物種內(nèi)則出現(xiàn)明顯的個體混合現(xiàn)象,提示這些物種間可能存在較強的雜交漸滲。

    核微衛(wèi)星數(shù)據(jù)顯示不同物種間共享大量等位基因且獨有等位基因很少,同時也存在明顯的核-質(zhì)沖突。上述結(jié)果可能是因物種間頻繁的漸滲雜交而導(dǎo)致,也可能是共同祖先等位基因保留(Duan et al., 2023; Feng et al., 2023)。由于核微衛(wèi)星突變速率較快(10-6~10-3每位點每代)且為中性標(biāo)記,大部分位點不受選擇影響(Vieira et al., 2016),因此共同祖先保留的微衛(wèi)星位點在不同物種中各自發(fā)生突變的可能性較高。在隨機漂變且缺少種間漸滲的情況下,不同物種共享相同微衛(wèi)星等位基因的可能性較低。因此,潤楠屬物種共享大量的等位基因不太可能因共同祖先等位基因的保留而引起。這表明潤楠屬物種間廣泛的漸滲雜交是不同物種共享大量等位基因且獨有等位基因少的合理解釋。

    由于種間漸滲雜交的廣泛存在,需要綜合利用形態(tài)特征、分子標(biāo)記對潤楠屬進行分類修訂。漸滲雜交現(xiàn)象在很多植物中普遍存在(Suarez-Gonzalez et al., 2018),尤其是在多樣性較高的屬內(nèi)。例如,F(xiàn)ogelqvist等(2015)使用STRUCTURE軟件采用384個SNPs和4個葉綠體SSR位點研究3種同域分布的柳屬(Salix)植物,發(fā)現(xiàn)一個柳屬雜交區(qū)具有較多的雜交漸滲個體;基于核SSR數(shù)據(jù),López-Caamal等 (2018)分析了Tithonia tubaeformis和圓葉腫柄菊(T. rotundifolia)的雜交漸滲現(xiàn)象,發(fā)現(xiàn)很多雜交漸滲個體具有不同親本的次生代謝產(chǎn)物。多數(shù)潤楠屬物種都分布于我國華南和西南地區(qū),本研究特意選取12個地理分布重疊的潤楠屬物種進行研究。此外,潤楠屬物種多為蟲媒傳粉(Watanabe et al., 2015; Tong et al., 2023),重疊的地理分布為潤楠屬物種共享傳粉昆蟲提供了條件。因此,潤楠屬內(nèi)同域分布物種間廣泛的漸滲雜交對分類和系統(tǒng)進化關(guān)系界定產(chǎn)生了嚴(yán)重影響。種間漸滲雜交通過基因流動導(dǎo)致不同物種基因庫的混合,從而模糊物種間的界線,導(dǎo)致形態(tài)特征重疊、相似。種間漸滲雜交形成復(fù)雜的系統(tǒng)進化關(guān)系網(wǎng)絡(luò),而非簡單的樹狀結(jié)構(gòu)(Stull et al., 2023),這種網(wǎng)絡(luò)結(jié)構(gòu)反映了不同物種間多次基因交換事件。多次漸滲雜交事件使得傳統(tǒng)的系統(tǒng)進化樹不再適用,需要使用網(wǎng)絡(luò)分析方法來更準(zhǔn)確地描述這些物種之間的演化關(guān)系。復(fù)雜系統(tǒng)進化關(guān)系的界定不僅有助于識別出新的雜交物種,而且還可揭示出傳統(tǒng)方法無法檢測的基因流路徑。此外,廣泛的漸滲雜交對潤楠屬的分類也提出了更高的要求,需要結(jié)合多種數(shù)據(jù)源,包括形態(tài)學(xué)特征、分子標(biāo)記、生態(tài)位信息等,更需要開發(fā)新的分析工具,綜合分析以得到更為可信的分類結(jié)果和種間漸滲雜交關(guān)系。

    綜上所述,增加葉綠體基因組序列并不能提高同域分布的潤楠屬種間系統(tǒng)進化關(guān)系的分辨率,微衛(wèi)星數(shù)據(jù)可以獲得分辨率更高的種間關(guān)系,但物種間的漸滲雜交對利用核基因數(shù)據(jù)進行系統(tǒng)進化關(guān)系分析的影響較大;不同潤楠屬物種間存在大量的共享等位基因及很少的獨有等位基因,表明潤楠屬物種間發(fā)生了廣泛的漸滲雜交,模糊了種間的界限。后續(xù)仍需使用具更高分辨率的核基因組序列,并結(jié)合傳統(tǒng)的形態(tài)分類方法,開發(fā)新的有效的分析工具,以理解潤楠屬內(nèi)復(fù)雜的系統(tǒng)進化關(guān)系,為該類群的保護和可持續(xù)利用提供基礎(chǔ)。

    致謝 感謝上海辰山植物園的胡永紅教授級高工、黃衛(wèi)昌教授級高工和上海植物園肖月娥教授級高工在論文準(zhǔn)備和修改中提供的幫助;感謝廣西植物研究所胡興華研究員對本研究提出的寶貴建議;感謝上海植物園于鳳揚在分子實驗中給予的幫助。

    參考文獻:

    ANTAO T, LOPES A, LOPES RJ, et al., 2008. LOSITAN: A workbench to detect molecular adaptation based on a Fst-outlier method" [J]. BMC Bioinformatics, 9(1): 323.

    BANDELT HJ, FORSTER P, RHL A, 1999. Median-joining networks for inferring intraspecific phylogenies" [J]. Molecular Biology and Evolution, 16(1): 37-48.

    BOUCKAERT R, VAUGHAN TG, BARIDO-SOTTANI J, et al., 2019. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis [J]. PLOS Computational Biology, 15(4): e1006650.

    CHEN JQ, LI L, LI J, et al., 2009. Bayesian inference of nrDNA ITS sequences from Machilus (Lauraceae) and its systematic significance" [J]. Plant Diversity, 31(2): 117-126." [陳俊秋, 李朗, 李捷, 等, 2009. 樟科潤楠屬植物ITS序列貝葉斯分析及其系統(tǒng)學(xué)意義 [J]. 云南植物研究, 31(2): 117-126.]

    CHENG W, ZHU C, XU W, et al., 2009. Chemical constituents of the bark of Machilus wangchiana and their biological activities" [J]. Journal of Natural Products, 72(12): 2145-2152.

    DUAN L, FU L, CHEN HF, 2023. Phylogenomic cytonuclear discordance and evolutionary histories of plants and animals" [J]. Science China Life Sciences, 66(12): 2946-2948.

    ELITH J, PHILLIPS SJ, HASTIE T, et al., 2011. A statistical explanation of MaxEnt for ecologists [J]. Diversity and Distributions, 17(1): 43-57.

    EVANNO GS, REGNAUT SJ, GOUDET J, 2005. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study [J]. Molecular Ecology, 14(8): 2611-2620.

    FAN D, LEI S, LIANG H, et al., 2022. More opportunities more species: Pleistocene differentiation and northward expansion of an evergreen broad-leaved tree species Machilus thunbergii (Lauraceae) in Southeast China" [J]. BMC Plant Biology, 22(1): 35.

    FENG C, WANG J, LISTON A, et al., 2023. Recombination variation shapes phylogeny and introgression in wild diploid strawberries [J]. Molecular Biology and Evolution, 40(3): msad049.

    FOGELQVIST J, VERKHOZINA AV, KATYSHEV AI, et al., 2015. Genetic and morphological evidence for introgression between three species of willows" [J]. BMC Evolutionary Biology, 15(1): 193.

    HOLLAND MM, PARSON W, 2011. GeneMarker HID: A reliable software tool for the analysis of forensic STR data" [J]. Journal of Forensic Sciences, 56(1): 29-35.

    HU JL, CI XQ, LIU ZF, et al., 2022. Assessing candidate DNA barcodes for Chinese and internationally traded timber species [J]. Molecular Ecology Resources, 22(4): 1478-1492.

    HUONG LT, CHAU DTM, AN NTG, et al., 2024. Essential oils of Lauraceae: Antimicrobial activity and constituents of essential oil from two Machilus species from Vietnam" [J]. Journal of Essential Oil Bearing Plants, 27(1): 177-187.

    JIANG K, TONG X, DING YQ, et al., 2021. Shifting roles of the East China Sea in the phylogeography of red nanmu in East Asia" [J]. Journal of Biogeography, 48(10): 2486-2501.

    KANEKO Y, LIAN C, WATANABE S, et al., 2012. Development of microsatellites in Machilus thunbergii (Lauraceae), a warm-temperate coastal tree species in Japan" [J]. American Journal of Botany, 99(7): E265-E267.

    KATOH K, ROZEWICKI J, YAMADA KD, 2019. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization" [J]. Briefings in Bioinformatics, 20(4): 1160-1166.

    KUMAR S, STECHER G, LI M, et al., 2018. MEGA X: Molecular evolutionary genetics analysis across computing platforms [J]. Molecular Biology and Evolution, 35(6): 1547-1549.

    LI L, LI J, ROHWER JG, et al., 2011. Molecular phylogenetic analysis of the Persea group (Lauraceae) and its biogeographic implications on the evolution of tropical and subtropical Amphi-Pacific disjunctions" [J]. American Journal of Botany, 98(9): 1520-1536.

    LI SG, 1982. Machilus" [M]// LI XW. Flora Reipublicae Popularis Sinicae. Beijing: Science Press: 201-224." [李樹剛, 1982. 潤楠屬 [M]//李錫文. 中國植物志. 北京: 科學(xué)出版社: 201-224.]

    LI SG, WEI FN, 1984. The research on Machilus Nees plants in South China and its neighboring areas" [J]. Guihaia, 4(2): 93-106. [李樹剛, 韋發(fā)南, 1984. 華南及其鄰近地區(qū)潤楠屬(Machilus Nees)植物的研究 [J]. 廣西植物, 4(2): 93-106.]

    LIN XZ, 2007a. Quantitative classification of plants of Machilus in Zhejiang Province" [J]. Scientia Silvae Sinicae, 43(11): 151-156." [林夏珍, 2007a. 浙江潤楠屬植物的數(shù)量分類" [J]. 林業(yè)科學(xué), 43(11): 151-156.]

    LIN XZ, 2007b. The systematic study of Machilus plants in Zhejiang Province [M]. Beijing: China Forestry Press." [林夏珍, 2007b. 浙江省潤楠屬植物系統(tǒng)研究 [M]. 北京: 中國林業(yè)出版社.]

    LIU ZF, MA H, CI XQ, et al., 2021.Can plastid genome sequencing be used for species identification in Lauraceae?" [J]. Botanical Journal of the Linnean Society, 197(1): 1-14.

    LIU ZF, MA H, ZHANG XY, et al., 2022. Do taxon-specific DNA barcodes improve species discrimination relative to universal barcodes in Lauraceae? [J]. Botanical Journal of the Linnean Society, 199(4): 741-753.

    LPEZ-CAAMAL A, REYES-CHILPA R, TOVAR-SNCHEZ E, 2018. Hybridization between Tithonia tubaeformis and T. rotundifolia (Asteraceae) evidenced by nSSR and secondary metabolites" [J]. Plant Systematics and Evolution, 304(3): 313-326.

    PHILLIPS SJ, DUDK M, 2008. Modeling of species distributions with MaxEnt: New extensions and a comprehensive evaluation" [J]. Ecography, 31(2): 161-175.

    PRITCHARD JK, STEPHENS MJ, DONNELLY PJ, 2000. Inference of population structure using multilocus genotype data" [J]. Genetics, 155(2): 945-959.

    RAMBAUT A, DRUMMOND AJ, XIE D, et al., 2018. Posterior summarization in Bayesian phylogenetics using tracer 1.7 [J]. Systematic Biology, 67(5): 901-904.

    ROHWER JG, LI J, RUDOLPH B, et al., 2009. Is Persea (Lauraceae) monophyletic? Evidence from nuclear ribosomal ITS sequences" [J]. Taxon, 58(4): 1153-1167.

    RONQUIST F, TESLENKO M, VAN DER MARK P, et al., 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space" [J]. Systematic Biology, 61(3): 539-542.

    SMITH DR, 2015. Mutation rates in plastid genomes: They are lower than you might think" [J]. Genome Biology and Evolution, 7(5): 1227-1234.

    SONG Y, YU WB, TAN Y, et al., 2017. Evolutionary comparisons of the chloroplast genome in Lauraceae and insights into loss events in the Magnoliids" [J]. Genome Biology and Evolution, 9(9): 2354-2364.

    STAMATAKIS A, 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies" [J]. Bioinformatics, 30(9): 1312-1313.

    STULL GW, PHAM KK, SOLTIS PS, et al., 2023. Deep reticulation: The long legacy of hybridization in vascular plant evolution" [J]. The Plant Journal, 114(4): 743-766.

    SUAREZ-GONZALEZ A, LEXER C, CRONK QCB, 2018. Adaptive introgression: a plant perspective" [J]. Biology Letters, 14(3): 20170688.

    TANG GG, XIANG QB, 1995. Pollen morphology of the family Lauraceae in China" [J]. Acta Phytotaxonomica Sinica, 33(2): 161-171." [湯庚國, 向其柏, 1995. 樟科植物花粉形態(tài)研究 [J]. 植物分類學(xué)報, 33(2): 161-170.]

    TO DC, HOANG LM, NGUYEN HT, et al., 2023.Dataset on the compounds from the leaves of Vietnamese Machilus thunbergii and their anti-inflammatory activity" [J]. Data in Brief, 51: 109713.

    TONG ZY, WU LY, FENG HH, et al., 2023. New calculations indicate that 90% of flowering plant species are animal-pollinated [J]. National Science Review, 10(10): nwad219.

    VAN OOSTERHOUT C, HUTCHINSON WF, WILLS DPM, et al., 2004. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data" [J]. Molecular Ecology Notes, 4(3): 535-538.

    VIEIRA MLC, SANTINI L, DINIZ AL, et al., 2016.Microsatellite markers: What they mean and why they are so useful" [J]. Genetics and Molecular Biology, 39(3): 312-328.

    WANG YG, WEI FN, 2003. A comparative study of pollen morphology between the closely related genera Machilus and Phoebe of the Lauraceae family" [J]. Guihaia, 23(1): 27-30." [王玉國, 韋發(fā)南, 2003. 樟科兩個近緣屬潤楠屬和楠木屬花粉形態(tài)的比較研究 [J]. 廣西植物, 23(1): 27-30.]

    WARREN DL, GLOR RE, TURELLI M, 2010. ENMTools: A toolbox for comparative studies of environmental niche models" [J]. Ecography, 33(3): 607-611.

    WATANABE S, NOMA N, NISHIDA T, 2015. Flowering phenology and mating success of the heterodichogamous tree Machilus thunbergii Sieb. et Zucc (Lauraceae)" [J]. Plant Species Biology, 31(1): 29-37.

    WEI FN, VAN DER WERFF H, WU Z, et al., 2008. Machilus [M]//WU ZY, RAVEN PH, HONG DY. Flora of China. Beijing: Science Press; St. Louis: Missouri Botanical Garden Press, 7: 102-254.

    WEI FN, TANG SC, 2006. On the circumscription of Machilus and Persea" [J]. Acta Phytotaxonomica Sinica, 44(4): 437-442." [韋發(fā)南, 唐賽春, 2006. 關(guān)于樟科潤楠屬和鱷梨屬的分類界線問題 [J]. 植物分類學(xué)報, 44(4): 437-442.]

    WEI FN, WANG YG, HE SQ, 2001.A taxonomic revision on some species of Machilus Nees (Lauraceae) from China" [J]. Guihaia, 21(3):191-194." [韋發(fā)南, 王玉國, 何順清, 2001. 中國樟科潤楠屬植物一些種類修訂 [J]. 廣西植物, 21(3): 191-194. ]

    XIAO TW, YAN HF, GE XJ, 2022. Plastid phylogenomics of tribe Perseeae (Lauraceae) yields insights into the evolution of East Asian subtropical evergreen broad-leaved forests" [J]. BMC Plant Biology, 22(1): 32.

    YANG SS, WU HC, HWANG TL, et al., 2022.Anti-inflammatory butanolides and lignanoids from the root of Machilus zuihoensis var. mushaensis [J]. Bioorganic Chemistry, 129: 106166.

    YANG SY, JIANG GH, SUN QH, et al., 2022.Chemical constituents from Machilus salicina Hance" [J]. Biochemical Systematics and Ecology, 103: 104432.

    YANG SY, JIANG GH, WANG CX, et al., 2023. Compounds from the leaves and stems of Machilus tenuipilis [J]. Chemistry of Natural Compounds, 59: 406-409.

    ZHUANG XY, 1997. The taxonomic study of Machilus species in Hong Kong [J]. Guihaia, 17(4): 4-7." [莊雪影, 1997. 香港潤楠屬植物的分類研究 [J]. 廣西植物, 17(4): 4-7.]

    ZHUANG XY, ZHANG Y, SUN TX, 2002. Leaf epidermis and their taxonomic significance in Machilus of Hong Kong [J]. Journal of South China Agricultural University, 23(1): 52-54." [莊雪影, 張粵, 孫同興, 2002. 香港潤楠屬植物葉表皮形態(tài)及分類學(xué)意義 [J]. 華南農(nóng)業(yè)大學(xué)學(xué)報, 23(1): 52-54.]

    ZHONGY, XIA NH, 2010. Leaf epidermal features of Machilus (Lauraceae) from China and their systematic significance [J]. Journal of Tropical and Subtropical Botany, 18(2): 109-121." [鐘義, 夏念和, 2010. 國產(chǎn)潤楠屬植物的葉表皮特征及其系統(tǒng)學(xué)意義 [J]. 熱帶亞熱帶植物學(xué)報, 18(2): 109-121.]

    ZHU Q, LIAO BY, LI P, et al., 2017. Phylogeographic pattern suggests a general northeastward dispersal in the distribution of Machilus pauhoi in South China" [J]. PLoS ONE, 12: e0184456.

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

    基金項目:" 上海市綠化和市容管理局項目(G232408,F(xiàn)132415)。

    第一作者: 蔣凱 (1987—),博士,工程師,主要從事植物保育遺傳學(xué)研究,(E-mail)jiangkai@csnbgsh.cn。

    *通信作者:" 陳小勇,博士,教授,研究方向為動植物互作及進化,(E-mail)xychen@des.ecnu.edu.cn。

    精品久久久久久久人妻蜜臀av| 男人舔奶头视频| 国产精品av视频在线免费观看| 欧美性感艳星| 下体分泌物呈黄色| 有码 亚洲区| 在线观看美女被高潮喷水网站| 国产人妻一区二区三区在| 啦啦啦中文免费视频观看日本| 久久国内精品自在自线图片| 日韩视频在线欧美| 边亲边吃奶的免费视频| 国产成人精品婷婷| 乱系列少妇在线播放| 国产大屁股一区二区在线视频| 国产白丝娇喘喷水9色精品| 久久99蜜桃精品久久| 国产一区二区亚洲精品在线观看| 少妇人妻精品综合一区二区| 97超碰精品成人国产| 男的添女的下面高潮视频| 亚洲精品乱久久久久久| 日韩亚洲欧美综合| av线在线观看网站| 成人免费观看视频高清| 在线观看一区二区三区| 欧美xxxx黑人xx丫x性爽| 中文字幕制服av| 亚洲性久久影院| 久久久亚洲精品成人影院| 自拍欧美九色日韩亚洲蝌蚪91 | 菩萨蛮人人尽说江南好唐韦庄| 国产综合精华液| 韩国av在线不卡| 少妇的逼水好多| 91久久精品电影网| 亚洲最大成人av| 黄色视频在线播放观看不卡| 成人综合一区亚洲| 亚洲精品自拍成人| 免费大片18禁| 欧美极品一区二区三区四区| 在线a可以看的网站| 日本wwww免费看| av在线播放精品| 国产精品久久久久久精品古装| 久久亚洲国产成人精品v| 男人舔奶头视频| 精品一区在线观看国产| av网站免费在线观看视频| 亚洲自偷自拍三级| 熟女人妻精品中文字幕| 欧美日本视频| 亚洲精品,欧美精品| 久久久久久久午夜电影| 久久久久性生活片| www.av在线官网国产| 日韩伦理黄色片| 在线观看免费高清a一片| 国产爽快片一区二区三区| 岛国毛片在线播放| 亚洲综合色惰| 一区二区三区乱码不卡18| 又爽又黄a免费视频| av天堂中文字幕网| 亚洲精品乱码久久久久久按摩| av免费观看日本| 国产欧美另类精品又又久久亚洲欧美| 久久久精品94久久精品| 黄色怎么调成土黄色| av国产久精品久网站免费入址| 成人漫画全彩无遮挡| 午夜免费鲁丝| 91精品伊人久久大香线蕉| 最近中文字幕高清免费大全6| 一二三四中文在线观看免费高清| 91久久精品国产一区二区成人| 亚洲av电影在线观看一区二区三区 | a级毛片免费高清观看在线播放| 国产成人a∨麻豆精品| 五月天丁香电影| 亚洲精品影视一区二区三区av| 午夜福利视频精品| 一本久久精品| 人妻少妇偷人精品九色| 乱系列少妇在线播放| 成年免费大片在线观看| 亚洲成人一二三区av| 女人被狂操c到高潮| 国产爽快片一区二区三区| 国产欧美日韩精品一区二区| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 亚洲国产精品国产精品| 边亲边吃奶的免费视频| 丰满乱子伦码专区| 成人亚洲精品av一区二区| 黄色一级大片看看| 亚洲av免费在线观看| 国产精品嫩草影院av在线观看| 免费人成在线观看视频色| 男的添女的下面高潮视频| a级毛片免费高清观看在线播放| 久久久久性生活片| 亚洲成色77777| 少妇人妻久久综合中文| 日韩免费高清中文字幕av| 国产色爽女视频免费观看| 免费观看在线日韩| 国产视频首页在线观看| 精品99又大又爽又粗少妇毛片| 欧美最新免费一区二区三区| 七月丁香在线播放| 国产精品嫩草影院av在线观看| a级一级毛片免费在线观看| 欧美亚洲 丝袜 人妻 在线| 一级毛片久久久久久久久女| 日本三级黄在线观看| 精品国产一区二区三区久久久樱花 | 午夜福利在线在线| 久久久久久久精品精品| 亚洲在线观看片| 日韩欧美 国产精品| 久久久精品欧美日韩精品| 女人被狂操c到高潮| 亚洲色图av天堂| 你懂的网址亚洲精品在线观看| eeuss影院久久| 久久韩国三级中文字幕| 人妻 亚洲 视频| 中文欧美无线码| 国产成人免费无遮挡视频| 极品少妇高潮喷水抽搐| 亚洲激情五月婷婷啪啪| 男男h啪啪无遮挡| 在线观看人妻少妇| 亚洲精品乱码久久久久久按摩| 午夜亚洲福利在线播放| 涩涩av久久男人的天堂| av网站免费在线观看视频| 亚洲av电影在线观看一区二区三区 | 国内少妇人妻偷人精品xxx网站| 久久久精品欧美日韩精品| 中文字幕亚洲精品专区| 国产熟女欧美一区二区| 亚州av有码| 日韩 亚洲 欧美在线| 亚洲精品一二三| 精华霜和精华液先用哪个| 大片免费播放器 马上看| 日本一本二区三区精品| 只有这里有精品99| 亚洲av免费高清在线观看| 国产精品人妻久久久久久| 亚洲熟女精品中文字幕| 禁无遮挡网站| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 成年免费大片在线观看| 欧美xxⅹ黑人| 国产黄片视频在线免费观看| 在线观看一区二区三区激情| 国产精品无大码| 嫩草影院入口| 人人妻人人看人人澡| 日韩伦理黄色片| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 最近手机中文字幕大全| 亚洲,欧美,日韩| 亚洲综合色惰| 国产精品久久久久久精品古装| 成人亚洲欧美一区二区av| 日韩av免费高清视频| 啦啦啦啦在线视频资源| 1000部很黄的大片| 国产精品久久久久久av不卡| 国产精品不卡视频一区二区| www.色视频.com| 乱系列少妇在线播放| 国产欧美亚洲国产| 亚洲av男天堂| 99热全是精品| 欧美日韩一区二区视频在线观看视频在线 | 国产淫片久久久久久久久| 日韩成人伦理影院| 国产一区亚洲一区在线观看| 蜜桃亚洲精品一区二区三区| 午夜福利网站1000一区二区三区| 欧美一区二区亚洲| 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 乱系列少妇在线播放| 九草在线视频观看| 国产亚洲最大av| 国产91av在线免费观看| 只有这里有精品99| 成年人午夜在线观看视频| 高清视频免费观看一区二区| 亚洲国产欧美在线一区| 久久久久久久国产电影| 免费看a级黄色片| 亚洲婷婷狠狠爱综合网| 在线精品无人区一区二区三 | 禁无遮挡网站| 大香蕉久久网| 一级二级三级毛片免费看| 高清在线视频一区二区三区| 91aial.com中文字幕在线观看| 性色avwww在线观看| 免费av观看视频| 亚洲av成人精品一二三区| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看| 午夜激情久久久久久久| 久久久久久久大尺度免费视频| 1000部很黄的大片| 一区二区三区免费毛片| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 亚洲在线观看片| 亚洲精品日韩av片在线观看| 51国产日韩欧美| 成人免费观看视频高清| 午夜亚洲福利在线播放| 久久久欧美国产精品| 亚洲综合色惰| 亚洲人成网站在线播| 精品人妻偷拍中文字幕| 婷婷色综合www| 一个人看视频在线观看www免费| 亚洲av福利一区| 成年免费大片在线观看| 日韩人妻高清精品专区| videossex国产| 中国美白少妇内射xxxbb| 99久久中文字幕三级久久日本| 最近的中文字幕免费完整| 亚洲经典国产精华液单| 观看美女的网站| 亚洲精品一区蜜桃| 在线看a的网站| 又粗又硬又长又爽又黄的视频| 免费黄色在线免费观看| 国产探花极品一区二区| 国产精品秋霞免费鲁丝片| 插阴视频在线观看视频| 尤物成人国产欧美一区二区三区| 大香蕉久久网| 久久国产乱子免费精品| 91久久精品国产一区二区三区| 水蜜桃什么品种好| av在线亚洲专区| 午夜日本视频在线| 一级a做视频免费观看| 禁无遮挡网站| 男女无遮挡免费网站观看| 一级毛片黄色毛片免费观看视频| 亚洲av男天堂| 一级二级三级毛片免费看| 女人被狂操c到高潮| 欧美最新免费一区二区三区| 听说在线观看完整版免费高清| 日本一本二区三区精品| 高清日韩中文字幕在线| 亚洲图色成人| 中文欧美无线码| 人妻少妇偷人精品九色| 我的老师免费观看完整版| 一二三四中文在线观看免费高清| 搞女人的毛片| 我要看日韩黄色一级片| 成年女人在线观看亚洲视频 | 国产淫语在线视频| 国产一区二区三区av在线| 偷拍熟女少妇极品色| 日韩免费高清中文字幕av| 又大又黄又爽视频免费| 看免费成人av毛片| 日韩不卡一区二区三区视频在线| 亚洲婷婷狠狠爱综合网| 亚洲精品一区蜜桃| 成人亚洲欧美一区二区av| 亚洲精品乱久久久久久| 午夜免费鲁丝| 欧美高清性xxxxhd video| 日韩国内少妇激情av| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 国产成人a区在线观看| 欧美xxⅹ黑人| 久久久久久伊人网av| 最后的刺客免费高清国语| 国产 一区 欧美 日韩| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av天美| 99久久九九国产精品国产免费| 国产成人91sexporn| 男插女下体视频免费在线播放| 国产一区有黄有色的免费视频| 身体一侧抽搐| 国产精品国产三级国产av玫瑰| 久久久色成人| 成人欧美大片| av在线观看视频网站免费| 毛片一级片免费看久久久久| 免费看光身美女| 国产亚洲最大av| 男女无遮挡免费网站观看| 99久久中文字幕三级久久日本| 久久久久久久国产电影| 男女下面进入的视频免费午夜| 免费看av在线观看网站| 欧美日韩视频精品一区| 欧美3d第一页| 国产精品无大码| 久久精品国产亚洲网站| 久久久午夜欧美精品| 国产 精品1| 制服丝袜香蕉在线| 国产精品嫩草影院av在线观看| 精品国产露脸久久av麻豆| 卡戴珊不雅视频在线播放| 日本熟妇午夜| 亚洲欧美日韩卡通动漫| 女人久久www免费人成看片| 精品少妇久久久久久888优播| 欧美日本视频| 久久综合国产亚洲精品| 在线精品无人区一区二区三 | 王馨瑶露胸无遮挡在线观看| 国产黄片美女视频| 六月丁香七月| 欧美区成人在线视频| 午夜福利视频精品| .国产精品久久| 亚洲色图综合在线观看| 国产精品一区二区三区四区免费观看| 欧美zozozo另类| 特级一级黄色大片| 少妇裸体淫交视频免费看高清| 熟女av电影| 精品一区在线观看国产| 尾随美女入室| 校园人妻丝袜中文字幕| 亚洲欧美日韩东京热| 国产综合精华液| 亚洲性久久影院| av在线天堂中文字幕| 在线播放无遮挡| 日韩强制内射视频| 人人妻人人爽人人添夜夜欢视频 | 日本午夜av视频| 色视频www国产| av在线观看视频网站免费| 免费在线观看成人毛片| 国产黄频视频在线观看| 91精品国产九色| 日本与韩国留学比较| a级毛色黄片| freevideosex欧美| 麻豆精品久久久久久蜜桃| 汤姆久久久久久久影院中文字幕| 韩国高清视频一区二区三区| 国产午夜精品久久久久久一区二区三区| 青春草国产在线视频| freevideosex欧美| 成年免费大片在线观看| 美女主播在线视频| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 少妇人妻一区二区三区视频| 精品99又大又爽又粗少妇毛片| 男女下面进入的视频免费午夜| 99热这里只有精品一区| 亚洲av.av天堂| 在线观看人妻少妇| 午夜激情久久久久久久| 国产精品国产av在线观看| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美在线精品| 精品久久久久久久久亚洲| 久久久久精品性色| 免费观看av网站的网址| 丰满少妇做爰视频| 尾随美女入室| 久久人人爽人人爽人人片va| 欧美极品一区二区三区四区| 亚洲精品,欧美精品| 最近2019中文字幕mv第一页| 国产老妇女一区| 毛片一级片免费看久久久久| 内射极品少妇av片p| 亚洲精品国产成人久久av| 国产久久久一区二区三区| 2021天堂中文幕一二区在线观| 性插视频无遮挡在线免费观看| 不卡视频在线观看欧美| 久久影院123| 精品视频人人做人人爽| 免费av毛片视频| 婷婷色av中文字幕| 国产免费福利视频在线观看| 日日摸夜夜添夜夜爱| 国产在线男女| 成人国产av品久久久| 亚洲伊人久久精品综合| 99热这里只有是精品在线观看| 国产成人精品久久久久久| 一个人观看的视频www高清免费观看| 99久久精品一区二区三区| 少妇人妻久久综合中文| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av片在线观看秒播厂| 国产男人的电影天堂91| 久久久精品免费免费高清| 亚洲最大成人av| 欧美激情国产日韩精品一区| 天天躁日日操中文字幕| 国产免费福利视频在线观看| 秋霞在线观看毛片| 看非洲黑人一级黄片| 亚洲国产精品国产精品| 又粗又硬又长又爽又黄的视频| 赤兔流量卡办理| 综合色丁香网| 夫妻性生交免费视频一级片| 国产老妇伦熟女老妇高清| 男女边摸边吃奶| 啦啦啦啦在线视频资源| 九九久久精品国产亚洲av麻豆| 真实男女啪啪啪动态图| 成人特级av手机在线观看| 人人妻人人爽人人添夜夜欢视频 | 99热国产这里只有精品6| 日日啪夜夜爽| 国产av不卡久久| 欧美变态另类bdsm刘玥| 国产精品女同一区二区软件| 国产精品一区二区在线观看99| 精品久久久精品久久久| 国产永久视频网站| 在线看a的网站| 久久精品夜色国产| 国产一区二区在线观看日韩| 99久国产av精品国产电影| 交换朋友夫妻互换小说| 亚洲成人一二三区av| 少妇丰满av| 在线观看人妻少妇| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 最近中文字幕高清免费大全6| 三级男女做爰猛烈吃奶摸视频| 大香蕉97超碰在线| 汤姆久久久久久久影院中文字幕| 狠狠精品人妻久久久久久综合| 成人国产av品久久久| 在线观看国产h片| 成人高潮视频无遮挡免费网站| 一区二区三区免费毛片| 22中文网久久字幕| 嘟嘟电影网在线观看| 久久久久久伊人网av| 成人欧美大片| 一区二区三区四区激情视频| 在线天堂最新版资源| 狂野欧美白嫩少妇大欣赏| 哪个播放器可以免费观看大片| 国产爽快片一区二区三区| 精品午夜福利在线看| 在线播放无遮挡| 一个人观看的视频www高清免费观看| 免费观看av网站的网址| 国产精品国产三级国产专区5o| 最新中文字幕久久久久| 青春草视频在线免费观看| 欧美一区二区亚洲| 美女脱内裤让男人舔精品视频| 国产 精品1| 精华霜和精华液先用哪个| 丝袜喷水一区| 听说在线观看完整版免费高清| 亚洲色图av天堂| 亚洲精品日本国产第一区| 只有这里有精品99| av黄色大香蕉| 搡老乐熟女国产| 三级经典国产精品| 熟女人妻精品中文字幕| eeuss影院久久| 男女那种视频在线观看| 久久影院123| 老司机影院毛片| 国产美女午夜福利| 精品国产三级普通话版| 国产欧美日韩精品一区二区| 女人十人毛片免费观看3o分钟| 亚洲欧美精品自产自拍| 又爽又黄无遮挡网站| 日韩大片免费观看网站| 中国美白少妇内射xxxbb| 久久午夜福利片| 男女啪啪激烈高潮av片| 免费黄网站久久成人精品| 高清av免费在线| 国产探花极品一区二区| 国产精品嫩草影院av在线观看| 亚洲自拍偷在线| 最近最新中文字幕大全电影3| 激情 狠狠 欧美| 国产男女超爽视频在线观看| 亚洲精品第二区| 99久久精品热视频| 久久久欧美国产精品| 亚洲婷婷狠狠爱综合网| 亚洲av国产av综合av卡| 尤物成人国产欧美一区二区三区| 久久久久久久午夜电影| 久久女婷五月综合色啪小说 | 成年女人看的毛片在线观看| 国产69精品久久久久777片| 美女国产视频在线观看| 熟女av电影| 一本久久精品| 91aial.com中文字幕在线观看| 久久久a久久爽久久v久久| 一级二级三级毛片免费看| 插逼视频在线观看| 日本黄色片子视频| 色吧在线观看| 在线精品无人区一区二区三 | 男女无遮挡免费网站观看| 国产黄片视频在线免费观看| 亚洲av日韩在线播放| 如何舔出高潮| 精品久久久久久久久av| 白带黄色成豆腐渣| 最近2019中文字幕mv第一页| 九九久久精品国产亚洲av麻豆| 国产午夜精品一二区理论片| 18禁在线无遮挡免费观看视频| 亚洲一区二区三区欧美精品 | 综合色av麻豆| 九九在线视频观看精品| 天堂俺去俺来也www色官网| 国产成人免费观看mmmm| 人人妻人人爽人人添夜夜欢视频 | 国产男女内射视频| 亚洲色图av天堂| 97精品久久久久久久久久精品| www.色视频.com| 中国美白少妇内射xxxbb| 免费观看的影片在线观看| 日本与韩国留学比较| 国产成人a区在线观看| 91在线精品国自产拍蜜月| 午夜福利网站1000一区二区三区| 狠狠精品人妻久久久久久综合| 夫妻午夜视频| 人人妻人人澡人人爽人人夜夜| 最近2019中文字幕mv第一页| 免费av毛片视频| 亚洲,欧美,日韩| 五月天丁香电影| 99热国产这里只有精品6| 午夜视频国产福利| 久久午夜福利片| 亚洲欧美精品专区久久| 精品一区二区三卡| 少妇猛男粗大的猛烈进出视频 | 亚洲精品国产av蜜桃| 一级毛片电影观看| 国产有黄有色有爽视频| 亚洲精品中文字幕在线视频 | 国产黄色免费在线视频| 老女人水多毛片| 国产黄片视频在线免费观看| 在线观看av片永久免费下载| 尤物成人国产欧美一区二区三区| 亚洲av一区综合| 成人毛片a级毛片在线播放| 国产久久久一区二区三区| 听说在线观看完整版免费高清| 亚洲欧美一区二区三区黑人 | 91精品国产九色| 国产一区二区亚洲精品在线观看| 色网站视频免费| 婷婷色av中文字幕| 国产精品三级大全| av福利片在线观看| 成年免费大片在线观看| 啦啦啦在线观看免费高清www| 国产探花极品一区二区| 亚洲人与动物交配视频| 亚洲一级一片aⅴ在线观看| 日韩,欧美,国产一区二区三区| 最近2019中文字幕mv第一页| 最近最新中文字幕免费大全7| 久久鲁丝午夜福利片| 最近的中文字幕免费完整| 免费不卡的大黄色大毛片视频在线观看| 久久亚洲国产成人精品v| 久久久久网色| 国产精品久久久久久精品电影| 不卡视频在线观看欧美| 少妇人妻精品综合一区二区| 日韩三级伦理在线观看| 网址你懂的国产日韩在线| 成人鲁丝片一二三区免费|