異質(zhì)生境下山東銀蓮花的轉(zhuǎn)錄組分析及EST-SSR分子標(biāo)記開發(fā)

摘 要： "為探究稀有物種山東銀蓮花在全光照的山頂灌叢和陰暗的針闊混交林下兩種不同生境中的生態(tài)適應(yīng)機制，并開發(fā)其EST-SSR分子標(biāo)記，該研究利用Illumina高通量測序技術(shù)對開花期的山東銀蓮花葉片進行轉(zhuǎn)錄組測序，獲取其功能注釋和差異表達基因（DEGs）。結(jié)果表明：（1）轉(zhuǎn)錄組測序共得到53 536條Unigenes序列，其中27 448條成功獲得注釋。（2）差異表達基因5 635個，1 600個在山頂灌叢的山東銀蓮花中上調(diào)表達，其余4 035個下調(diào)表達。有2 460個差異表達基因注釋到GO數(shù)據(jù)庫2 533個三級條目中，1 051個差異表達基因注釋到KEGG數(shù)據(jù)庫的113條代謝通路中。（3）山東銀蓮花適應(yīng)于異質(zhì)生境的代謝通路主要涉及光合作用-天線蛋白通路和類黃酮生物合成通路，光合作用-天線蛋白通路中l(wèi)hca5基因上調(diào)表達，lhcb1、lhcb2和lhcb3基因下調(diào)表達，類黃酮生物合成通路中chs、c4h、f3′h、f3h、fls、ans、chi、ccoaomt和hct基因均上調(diào)表達。（4）從山東銀蓮花轉(zhuǎn)錄組數(shù)據(jù)中共獲得7 146個SSR位點分布于6 006條Unigenes序列中，共計106種重復(fù)基序，優(yōu)勢重復(fù)基序為單核苷酸重復(fù)。設(shè)計合成100對EST-SSR引物中共有68對引物具有有效性，其中11對具有多態(tài)性，共擴增24個多態(tài)性片段。該研究首次開發(fā)了山東銀蓮花EST-SSR分子標(biāo)記，為山東銀蓮花的保護和利用提供了重要的分子標(biāo)記資源。

關(guān)鍵詞： 異質(zhì)生境， 山東銀蓮花， 轉(zhuǎn)錄組， EST-SSR， 生態(tài)適應(yīng)機制

中圖分類號： "Q943

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

文章編號： "1000-3142（2025）01-0095-13

基金項目： "國家自然科學(xué)基金面上項目（31971546）。

第一作者： 單筱涵（1998—），碩士研究生，研究方向為植物生物學(xué)，（E-mail）shanxiaohan0920@163.com。

*通信作者： "卞?；ǎ┦?，教授，研究方向為植物生物學(xué)，（E-mail）fh_bian@163.com。

Transcriptome analysis and development of EST-SSR

molecular markers in Anemone shikokiana under

heterogeneous habitats

SHAN Xiaohan， AN Kang， ZHOU Chunxia， ZHANG Xin，

PANG Yujuan， LI Lixia， BIAN Fuhua

（ College of Life Sciences， Yantai University， Yantai 264005， Shandong， China ）

Abstract： "In order to explore， the ecological adaptation mechanisms of Anemone shikokiana in two distinct habitats， namely full-light hilltop scrub and shady mixed broadleaved-coniferous forest， and to develop its molecular markers， Illumina high-throughput sequencing technology for leaves of A. shikokiana collected during the flowering stage to obtain functional annotation and differentially expressed genes （DEGs）. The results were as follows： （1） A total of 53 536 Unigenes sequences were obtained， of which 27 448 were successfully annotated. （2） A total of 5 635 DEGs were obtained， 1 600 up-regulated and 4 035 down-regulated genes comparing A. shikokiana in full-light hilltop scrub and in shady mixed coniferous forest. A total of "2 460 DEGs were annotated to 2 533 tertiary entries in GO database. In addition， 1 051 DEGs were involved in 113 KEGG pathways. （3） The metabolic pathways adapted to heterogeneous habits in A. shikokiana mainly involved photosytheis-antenna protein pathway and flavonoid biosynthesis pathway. In the "photosynthesis-antenna protein pathway， the expression of lhca5 was up-regulated， while the

expressions of lhcb lhcb2 and lhcb3 were "down-regulated. Meanwhile， in the flavonoid biosynthesis pathway， the expressions of chs， c4h， f3′h， f3h， fls， ans， chi， ccoaomt and hct were all up-regulated. （4） A total of 6 006 Unigenes containing 7 146 SSRs were obtained from the transcriptome data of A.shikokiana. In the identified SSRs， the dominant repeat motifs were single nucleotide repeats in 106 repetitive motif types. Among the 100 pairs of EST-SSR primers， a total of 68 pairs were effective and 11 pairs with polymorphism， and 24 polymorphic fragments were amplified. Overall， in this paper， "for the first time EST-SSR molecular markers were developed， which would provide important molecular marker resource for the conservation and utilization of A. shikokiana.

Key words： heterogeneous habitat， Anemone shikokiana， transcriptome， EST-SSR， ecological adaptation mechanisms

山東銀蓮花（Anemone shikokiana）是毛茛科銀蓮花屬（Anemone）的多年生草本植物，間斷分布于中國膠東半島和日本四國島（王鷙等，2014）。在中國被列入《世界自然保護聯(lián)盟瀕危物種紅色名錄》（IUCN）易危（VU）等級，為典型的稀有物種（陳春利等，2018）。其根可入藥，有清熱解毒、止血除濕的功效（劉瓊，2014）；因其株型優(yōu)美、花大、花期長，在園藝上具有潛在的開發(fā)價值（侯元同和劉冰，2010）。

山東銀蓮花生長在海拔600 m以上的山頂灌叢和針闊混交林下兩種異質(zhì)生境中，山東灌叢常年全光照，土壤是以碎石、沙粒為主的砂質(zhì)性土壤；而針闊混交林下陰暗潮濕，土壤為營養(yǎng)豐富且疏松的腐殖土（Pang et al.， 2020）。岳喜元等（2023）研究表明，黃河三角洲的蘆葦通過調(diào)整個體大小與葉性狀適應(yīng)異質(zhì)生境下水分條件變化，王妍方等（2023）通過木棉表型可塑性探究其適應(yīng)干熱河谷和熱帶雨林兩種極端環(huán)境的機制。山東銀蓮花則通過改變?nèi)~片、根系的形態(tài)結(jié)構(gòu)及生理等策略適應(yīng)完全不同的生長環(huán)境（于文英等，2019；Pang er al.， 2020；逄玉娟等，2022）。然而，至今對該種適應(yīng)異質(zhì)生境的分子機制尚不清楚。植物適應(yīng)環(huán)境的本質(zhì)是基因特異性表達的結(jié)果，基因的特異性表達在一定程度上由轉(zhuǎn)錄水平進行調(diào)控（張椿雨等，2007）。轉(zhuǎn)錄組作為特定組織或細胞在某一發(fā)育階段或功能狀態(tài)下所有mRNA的集合，是研究基因功能與結(jié)構(gòu)的有效手段（魯艷輝等，2021）。運用Illumina高通量測序技術(shù)對植物轉(zhuǎn)錄組進行測序分析，可高效、精準(zhǔn)地獲得大量特定組織或細胞中表達的基因序列，實現(xiàn)在無參考基因組的基礎(chǔ)上進行分析（郭連安等，2021）。同時，轉(zhuǎn)錄組分析應(yīng)用于兩色金雞菊（孫浩男等，2024）、直立型扁蓿豆（烏日娜等，2022）等多種植物，并且已廣泛應(yīng)用于探究不同生境下植物的適應(yīng)性（周春苗等，2022；韋陳彬等，2022）。

分子標(biāo)記可用于植物的系統(tǒng)進化分析、功能基因標(biāo)記和分子輔助育種等，隨著高通量測序技術(shù)的發(fā)展，由轉(zhuǎn)錄組測序得到的EST-SSR分子標(biāo)記被作為最常用的分子標(biāo)記技術(shù)之一。其操作簡單，具有共顯性遺傳、穩(wěn)定性好和準(zhǔn)確度高等優(yōu)于其他分子標(biāo)記的特點（岳遠灝等，2022）。EST-SSR標(biāo)記開發(fā)成本低，適用于無參考基因組的物種，可有效地推動分子標(biāo)記技術(shù)在植物遺傳多樣性分析、指紋圖譜構(gòu)建、種質(zhì)資源保護、品種鑒定等方面的應(yīng)用（楊雄等，2021），目前已在大麻（邊境等，2023）、辣椒（田懷志等，2023）、楓香（李輝等，2023）等植物中開發(fā)與應(yīng)用。然而，關(guān)于山東銀蓮花分子標(biāo)記的報道尚不多見，未見公開報道的可利用 SSR 標(biāo)記，使得山東銀蓮花在分子標(biāo)記方面的研究難以深入。

本研究以異質(zhì)生境下山東銀蓮花為研究對象，利用Illumina高通量測序技術(shù)獲取其轉(zhuǎn)錄組數(shù)據(jù)，并進行組裝和功能注釋，得到其差異表達基因和大量豐富的SSR位點，隨后進行EST-SSR引物設(shè)計并加以驗證，以期探討山東銀蓮花在異質(zhì)生境中的生態(tài)適應(yīng)機制，并開發(fā)適用的EST-SSR分子標(biāo)記，填補山東銀蓮花在以上兩方面的研究空白，為今后其遺傳多樣性分析、種質(zhì)資源評價、功能基因標(biāo)記和分子輔助育種等研究奠定基礎(chǔ)，進一步實現(xiàn)物種的保護與利用。

1 材料與方法

1.1 實驗材料

本研究所用實驗材料均為開花期的山東銀蓮花葉片，其中轉(zhuǎn)錄組測序材料于2021年6月取自山東省煙臺市昆崳山國家級自然保護區(qū)的泰礴頂和寒風(fēng)嶺（表1），引物篩選材料于2021年9月取自青島嶗山和煙臺昆崳山（表2）。采集后用蒸餾水清洗干凈，立即放入液氮中冷凍。轉(zhuǎn)錄組測序由上海歐易生物醫(yī)學(xué)科技有限公司完成，引物篩選材料帶回實驗室后采用改良的CTAB法（閆苗苗等，2008）提取DNA，并利用1%瓊脂糖凝膠電泳對所提取的DNA進行檢測，將檢測合格的DNA統(tǒng)一稀釋至20 ng·μL^-1，放入-20 ℃冰箱中保存?zhèn)溆谩?/p>

1.2 轉(zhuǎn)錄組分析

De novo拼接：使用Trinity（Manfred et al.， 2011）軟件paired-end的拼接方法得到Transcript序列，根據(jù)序列相似性以及長度，挑選出最長的一條作為Unigene，之后再利用CD-HIT軟件聚類去冗余得到一套最終的Unigene，以此作為后續(xù)分析的參考序列。

功能注釋：利用DIAMOND（Buchfink et al.， 2014）軟件將Unigenes比對到NR、KOG、GO、Swiss-Prot、eggNOG、KEGG數(shù)據(jù)庫以及利用HMMER軟件比對Pfam數(shù)據(jù)庫來進行Unigenes的功能分析。

差異表達基因篩選：利用DESeq（Anders amp; Huber， 2012）軟件對各個樣本基因的Counts數(shù)目進行標(biāo)準(zhǔn)化處理，采用BaseMean值來估算表達量，計算差異倍數(shù)，并采用NB（負二項分布檢驗的方式）對Reads數(shù)進行差異顯著性檢驗，以差異倍數(shù)|log2FoldChange|gt;1及差異顯著性檢驗結(jié)果Plt;0.05為標(biāo)準(zhǔn)篩選差異表達基因。

1.3 引物設(shè)計與篩選

基于山東銀蓮花轉(zhuǎn)錄組測序結(jié)果，使用軟件MISA（Beier et al.， 2017）進行SSR預(yù)測，設(shè)置重復(fù)基序的長度為1、2、3、4、5和6，最小重復(fù)次數(shù)為10次、6次、5次、5次、5次和5次；兩個微衛(wèi)星之間的序列長度不超過100 bp即被認定為復(fù)合微衛(wèi)星。SSR作為一種用途廣泛的分子標(biāo)記，其多態(tài)性水平越高則應(yīng)用價值越大（鄭燕等，2012）。當(dāng)SSR長度≥20 bp時可顯示出較高的多態(tài)性，長度在12～20 bp之間相對較低，在12 bp以下則很低（Temnykh et al.， 2001）。根據(jù)SSR長度≥20 bp、去除復(fù)合型SSR位點和容易與多腺苷化作用相混雜的A/T重復(fù)基序的原則，將所有SSR位點進行初步篩選后隨機選擇100對SSR位點利用軟件Primer3（Rozen amp; Skaletsky， 2000）進行引物設(shè)計，命名為ANS001-ANS100，送至北京六合華大基因科技有限公司合成。

利用青島嶗山滑溜口所采集樣品通過1%瓊脂糖凝膠電泳對100對引物進行篩選，以檢測引物有效性。PCR擴增體系為20 μL：2×Master Mix 4 μL、上下游引物（10 μmol·L^-1）各2 μL、DNA模板2 μL和ddH2O10 μL。反應(yīng)程序：94 ℃預(yù)變性3 min；94 ℃變性30 s，55 ℃退火35 s，72 ℃延伸60 s，35個循環(huán)；72 ℃終延伸10 min，4 ℃保存。隨后將無法正確擴增或擴增產(chǎn)物大小與預(yù)期不符的引物篩除，可正確擴增的引物利用樣品YD1（煙臺昆崳山泰礴頂）、YL1（煙臺昆崳山老鐵山）、QD1（青島嶗山嶗頂景區(qū)丹爐峰附近）、QL1（青島嶗山滑溜口）通過8%非變性聚丙烯酰胺凝膠電泳和銀染顯色的方法進行目的條帶的分離和檢測，篩選出具有多態(tài)性的EST-SSR引物。

2 結(jié)果與分析

2.1 轉(zhuǎn)錄組數(shù)據(jù)分析

2.1.1 總體特征 對山東銀蓮花進行無參轉(zhuǎn)錄組分析，共獲得 42.34 G的Clean Data，各樣本的有效數(shù)據(jù)量分布在6.98～7.10 G，Q30堿基分布在93.68%～94.20%之間，證明本次測序數(shù)據(jù)準(zhǔn)確性良好，可用于后續(xù)分析。利用De novo拼接出Unigenes 53 536條，總長度為57 909 453 bp，最長Unigene為15 502 bp，最短Unigene為301 bp，平均長度為1 081.69 bp。用皮爾遜相關(guān)系數(shù)表示D和L樣品間基因表達水平相關(guān)性（低于0.8表示相關(guān)性較差），D組間生物學(xué)重復(fù)的相關(guān)性均大于0.868 9，L組間生物學(xué)重復(fù)之間相關(guān)性均大于0.831 8，表明樣品組間相關(guān)性較高，可以此作為后續(xù)分析的參考序列。

2.1.2 功能注釋 對拼接組裝的53 536條Unigenes進行常見功能數(shù)據(jù)庫注釋，在7個數(shù)據(jù)庫中均注釋成功的Unigenes有3 505條，占總數(shù)的6.55%，至少在1個數(shù)據(jù)庫中注釋成功的Unigenes有27 448條，占總數(shù)的51.27%（圖1）。在NR數(shù)據(jù)庫中成功注釋的條數(shù)最多，共有27 044條，占總數(shù)的50.52%，其余分別為eggNOG 庫（23 934，44.71%）、Swiss-Prot 庫（19 637，36.68%）、Pfam 庫（17 859，33.36%）、GO 庫（17 540，32.76%）、KOG 庫（5 481，29.19%）和KEGG 庫（5 48 10.23%）。按照NR數(shù)據(jù)庫中山東銀蓮花的Unigenes與其他物種匹配可得，相似程度最高的物種為耬斗菜（Aquilegia coerulea，56.72%），其他按匹配度由高到低依次為博落回（Macleaya cordata，6.08%）、鴉片罌粟（Papaver somniferum，4.06%）、蓮（Nelumbo nucifera，3.29%）、葡萄（Vitis vinifera，1.92%）、歐洲栓皮櫟（Quercus suber，0.94%）、月季花（Rosa chinensis，0.9%）、Sphaerulina musiva SO2202（0.73%）、中華獼猴桃（Actinidia chinensis var. chinensis，0.58%）、向日葵（Helianthus annuus，0.55%）和其他物種（24.24%）。

2.2 差異表達Unigene分析

2.2.1 差異表達Unigene篩選 對山東銀蓮花的轉(zhuǎn)錄組數(shù)據(jù)進行篩選，共挖掘差異表達Unigenes 5 635個，數(shù)量較多，其中與林下（L）的山東銀蓮花相比，山頂（D）的山東銀蓮花中有上調(diào)基因1 600個，下調(diào)基因4 035個（圖2），說明不同生長環(huán)境對山東銀蓮花基因表達影響較大。

2.2.2 差異表達Unigene富集分析 將山東銀蓮花的差異表達Unigenes進行GO（圖3）和KEGG富集分析（圖4），結(jié)果顯示共有2 460個差異表達基因被注釋到GO數(shù)據(jù)庫2 533個三級條目中，共有1 051個差異表達基因被注釋到KEGG數(shù)據(jù)庫的113條代謝通路中。

根據(jù)GO數(shù)據(jù)庫一級條目生物學(xué)過程、細胞成分和分子功能的不同，篩選3種分類中對應(yīng)差異基因數(shù)目大于2的GO二級條目，按照每個條目對應(yīng)的P值由小到大排序可得前10個富集分類，山東銀蓮花的差異表達基因主要顯著富集在脫落酸代謝過程、薄荷醇生物合成過程、硫脂生物合成過程等生物學(xué)過程中，胞質(zhì)大核糖體亞基、胞質(zhì)小核糖體亞基、核糖體等細胞成分中，芳醇脫氫酶（NADP+）活性、異戊烯醇脫氫酶活性、卡維醇脫氫酶活性等分子功能中。

將山東銀蓮花的差異表達基因比對到KEGG數(shù)據(jù)庫的參考代謝通路中，可劃分到6個類別中，其中與代謝相關(guān)通路的Unigenes數(shù)量最多，有480個，占比45.67%；其次遺傳信息加工相關(guān)通路的Unigenes有425個，占比40.44%；與環(huán)境信息處理相關(guān)通路的Unigenes有58個，占比5.52%；與生物系統(tǒng)相關(guān)通路的Unigenes有50個，占比4.76%；與細胞過程相關(guān)通路的Unigenes有37個，占比3.52%。共有10條KEGG通路達到極顯著水平（Plt;0.01），分別為核糖體、光合作用-天線蛋白、類黃酮生物合成、倍半萜和三萜生物合成、酮體的合成與降解、谷胱甘肽代謝、植物-病原互作、MAPK信號通路-植物、黃酮和黃酮醇的生物合成、二苯乙烯類、二芳基庚烷類和姜辣素的生物合成。

2.2.3 關(guān)鍵差異表達Unigene信息 從眾多通路中篩選出與山東銀蓮花適應(yīng)于異質(zhì)生境的主要通路（表3），其中15個差異表達Unigene屬于光合作用-天線蛋白通路，分別為1個上調(diào)基因lhca5（light-harvesting complex I chlorophyll a/b binding protein 5）、7個下調(diào)基因lhcb1（light-harvesting complex Ⅱ chlorophyll a/b binding protein 1）、6個下調(diào)基因lhcb2（light-harvesting complex Ⅱ chlorophyll a/b binding protein 2）、1個下調(diào)基因lhcb3（light-harvesting complex Ⅱ chlorophyll a/b binding protein 3）；11個差異表達Unigene屬于類黃酮生物合成通路，分別為1個chs（chalcone synthase）、1個c4h（trans-cinnamate 4-monooxygenase）、2個f3′h（flavonoid 3′-monooxygenase）、1個f3h（naringenin 3-dioxygenase）、1個fls（flavonol synthase）、1個ans（anthocyanidin synthase）、1個chi（chalcone isomerase）、1個ccoaomt（caffeoyl-CoA O-methyltransferase）、2個hct（shikimate O-hydroxycinnamoyltransferase），以上基因皆為上調(diào)表達。

2.3 SSR位點分析與引物開發(fā)

2.3.1 SSR位點總體特點 對山東銀蓮花轉(zhuǎn)錄組的53 536條Unigenes序列進行搜索，共檢出6 006條Unigenes序列中含有7 146個SSR位點，平均分布距離為8.1 kb，SSR長度在10～225 bp之間，SSR出現(xiàn)頻率（=SSR個數(shù)/總Unigene數(shù)）為13.35%，SSR發(fā)生頻率（=含SSR的Unigene數(shù)/總Unigene數(shù)）為11.22%。山東銀蓮花含有2個及以上的SSR位點的序列有907條，占全部序列的1.68％，含有復(fù)雜重復(fù)類型的SSR位點的序列有420條，占全部序列的0.78%。山東銀蓮花SSR類型較為豐富，共計106種重復(fù)基序（表4），優(yōu)勢重復(fù)基序為單核苷酸重復(fù)（3 248，45.59%），其次分別為三核苷酸重復(fù)（2 041，28.56%）、二核苷酸重復(fù)（1 645，23.02%）、六核苷酸重復(fù)（101，1.41%）、四核苷酸重復(fù)（81，1.13%）、五核苷酸重復(fù)（20，0.28%）（圖5）。

2.3.2 EST-SSR引物的開發(fā)與篩選 利用1%瓊脂糖凝膠電泳對從山東銀蓮花轉(zhuǎn)錄組數(shù)據(jù)中挑選的100對EST-SSR引物進行檢驗，可得68對有效擴增引物。為探究68對EST-SSR引物的多態(tài)性，采用8%非變性聚丙烯酰胺凝膠電泳檢驗來源于煙臺昆崳山泰礴頂（YD1）、煙臺昆崳山老鐵山（YL1）、青島嶗山嶗頂景區(qū)丹爐峰附近（QD1）、青島嶗山滑溜口（Q L1）的山東銀蓮花DNA擴增結(jié)果，68對EST-SSR引物均能在上述材料中擴增得到清晰條帶，其中有11對EST-SSR引物具有多態(tài)性（表5），共擴增24個多態(tài)性片段。以上結(jié)果表明，基于山東銀蓮花轉(zhuǎn)錄組數(shù)據(jù)開發(fā)EST-SSR引物是一種經(jīng)濟、可行的方法，可通過簡單的PCR擴增效果驗證得出EST-SSR引物的穩(wěn)定性、可重復(fù)性和多態(tài)性。

3 討論與結(jié)論

3.1 異質(zhì)生境中山東銀蓮花的適應(yīng)機制分析

高通量測序技術(shù)能全面快速地獲得植物的轉(zhuǎn)錄本序列信息，為植物基因發(fā)掘和次生代謝物調(diào)控等研究提供了巨大便利（Lu et al.， 2010）。本研究利用高通量測序技術(shù)對異質(zhì)生境中的山東銀蓮花進行轉(zhuǎn)錄組測序，通過對其差異表達基因進行富集分析，發(fā)現(xiàn)山東銀蓮花適應(yīng)兩種異質(zhì)生境的機制主要表現(xiàn)在光合作用和類黃酮生物合成兩方面。光合作用作為植物生長和一切代謝活動的生理基礎(chǔ)，受光照、水分、溫度等不同環(huán)境因子的影響（Fang et al.， 2022），在其調(diào)控過程中捕光復(fù)合物Ⅱ葉綠素a/b結(jié)合蛋白（light-harvesting complex Ⅱ chlorophyll a/b binding protein， LHC）發(fā)揮重要作用。LHC由核基因（lhc）編碼，具有迅速將光能傳到PSⅠ和PSⅡ反應(yīng)中心使光能轉(zhuǎn)化為化學(xué)能的作用，可以促進光合作用進行（Sapeta et al.， 2023）。lhc基因家族主要包含lhca和lhcb兩個進化類群，前者編碼植物光系統(tǒng)I中的LHCⅠ膜蛋白，捕光截面較大且可同時捕獲光譜中不同波長的光（Mozzo et al.， 2010），后者編碼植物光系統(tǒng)Ⅱ中的LHCⅡ膜蛋白，可結(jié)合類囊體膜上50%的色素形成葉綠體類囊體中含量最豐富的捕光天線蛋白（王云鶴，2020）。本研究發(fā)現(xiàn)，相比針闊混交林下生長的山東銀蓮花，山頂灌叢中的山東銀蓮花中l(wèi)hca5顯著上調(diào)，與Ganeteg等（2004）研究表明lhca5在強光條件下表達水平顯著提高的研究結(jié)果一致。此外，本研究還發(fā)現(xiàn)山頂?shù)纳綎|銀蓮花中l(wèi)hcb1、lhcb2和lhcb3顯著下調(diào)，蔣倩（2019）研究發(fā)現(xiàn)在干旱脅迫下水芹的lhcb基因表達下調(diào)可阻止葉綠體接收過多光能，使光合作用減弱，Rahele等（2022）的研究也表明lhcb基因表達量的減少導(dǎo)致番茄葉綠素含量降低以響應(yīng)干旱脅迫。因此，可推測lhcb1、lhcb2和lhcb3的下調(diào)可減少捕光天線蛋白的合成，減弱光合作用，減少水分的消耗，從而降低干旱對山東銀蓮花自身產(chǎn)生的影響。山頂灌叢與針闊混交林下比較，光照強、土壤濕度小是最明顯的不利于山東銀蓮花生存的逆境因素。本研究表明，lhca和lhcb兩類基因的上調(diào)和下調(diào)很巧妙地解決了這一難題，適應(yīng)了山頂灌叢不利的環(huán)境，實現(xiàn)非生物脅迫響應(yīng)。

類黃酮是植物體內(nèi)廣泛存在的次生代謝物質(zhì)，可吸收280～315 nm波長的紫外光，保護植物體器官，尤其是光合組織免受或少受輻射傷害（鄒鳳蓮等，2004），在植物生長發(fā)育和熱脅迫等非生物脅迫響應(yīng)過程中發(fā)揮多種生理功能（Chen et al.， 2023；Guo et al.， 2023；Zhuang et al.， 2023）。類黃酮生物合成的前體物質(zhì)來自于苯丙烷代謝，并受一系列關(guān)鍵酶的調(diào)控，如查爾酮合成酶（CHS）、查爾酮異構(gòu)酶（CHI）、黃酮醇合成酶（FLS）、黃酮3-羥基化酶（F3H）、黃酮3′-羥化酶（F3′H）、花色素合成酶 （ANS）等（Xiong et al.， 2016；Shen et al.， 2022；戴明潔，2022）。本研究發(fā)現(xiàn)山頂灌叢的山東銀蓮花中編碼相關(guān)酶的基因均較林下顯著上調(diào)，肖韻錚等（2020）研究表明chs和chi基因表達水平與類黃酮含量呈顯著正相關(guān)，表明山頂灌叢中的山東銀蓮花在生長過程中可以合成較多的類黃酮。山頂灌叢中的強光照，除了可見光，還包括紫外光，同時伴隨著午間的高溫，長時間的紫外光照射會造成植物DNA及其他結(jié)構(gòu)的破壞（Batschauer， 1993），加上高溫的影響，對植物將造成重大傷害甚至死亡。本研究中，山頂灌叢山東銀蓮花大量類黃酮的合成，可吸收多余的紫外光，減少紫外光對植株的破壞，同時緩解短時間的高溫，以適應(yīng)不良的生活環(huán)境。

3.2 山東銀蓮花的SSR位點分析與引物開發(fā)

山東銀蓮花的SSR類型相對豐富，從單核苷酸重復(fù)序列到六核苷酸重復(fù)序列不等。研究表明三核苷酸基序中AAG/CTT是主要的基序（Xiang et al.， 2023）， 與本研究中山東銀蓮花SSR位點分析結(jié)果一致，說明EST-SSR在其發(fā)生和進化過程中具有高度保守性。EST-SSR分子標(biāo)記源于功能基因表達序列，具有保守程度高、種間及屬間通用性好等優(yōu)點（婁永峰等，2022）。本研究基于拼接質(zhì)量較高的轉(zhuǎn)錄組數(shù)據(jù)，利用SSR位點設(shè)計并合成100對引物，其中有68對引物可以實現(xiàn)有效擴增，有11對引物在不同地域的山東銀蓮花材料中表現(xiàn)出多態(tài)性，具有多態(tài)性的EST-SSR分子標(biāo)記可用于遺傳多樣性分析與分子標(biāo)記輔助育種。因為國內(nèi)山東銀蓮花材料稀缺，獲取日本四國島的山東銀蓮花材料受限，所以在本研究中未進行相關(guān)聚類分析以及系統(tǒng)進化樹的構(gòu)建，后續(xù)將通過銀蓮花屬植物對本研究所得的山東銀蓮花EST-SSR多態(tài)性引物的通用性進行驗證，同時利用挖掘的差異表達基因設(shè)計EST-SSR引物用于區(qū)分不同生境的山東銀蓮花居群，為其進化研究和種質(zhì)資源保護奠定基礎(chǔ)。

綜上所述，本研究對山頂灌叢和針闊混交林下山東銀蓮花的轉(zhuǎn)錄組進行分析，表明其通過調(diào)控光合作用途徑和類黃酮生物合成途徑涉及的基因以適應(yīng)異質(zhì)生境， 并首次開發(fā)多態(tài)性較好的EST-SSR分子標(biāo)記，填補山東銀蓮花在此方面的空白，對其保護與利用具有重要意義。

參考文獻：

ANDERS S， HUBER W， 2012. BioC 2012： Differential expression of RNA-Seq data at the gene level "［J］. EMBL. https：//www.bioconductor.org/help/course-materials/2012/BioC2012/DESeq.pdf.

BATSCHAUER A， 1993. A plant gene for photolyase： An enzyme catalyzing the repair of UV-light-induced DNA damage ［J］. The Plant Journal， 4（4）： 705-709.

BEIER S， THIEL T， MNCH T， et al.， 2017. MISA-web： a web server for microsatellite prediction ［J］. Bioinformatics， 33（16）： 2583-2585.

BIAN J， WANG XN， CAO K， et al.， 2023. Genetic structure analysis and fingerprint construction of hemp by EST-SSR analysis ［J］. Journal of "Plant Genetic Resources， 24（6）： 1794-1808. "［邊境， 王曉楠， 曹焜， 等， 2023. 大麻EST-SSR遺傳結(jié)構(gòu)分析及指紋圖譜構(gòu)建 ［J］. 植物遺傳資源學(xué)報， 24（6）： 1794-1808.］

BUCHFINK B， XIE C， HUSON DH， 2014. Fast and sensitive protein alignment using DIAMOND ［J］. Nature Methods， 12（1）： 59-60.

CHEN LC， GU DZ， LIU Y， et al.， 2018. Establishment of tissue culture regeneration system of Anemone shikokiana ［J］. Jiangsu Agricultural Sciences， 46（24）： 47-50. "［陳春利， 顧德政， 劉毓， 等， 2018. 山東銀蓮花組織培養(yǎng)植株再生體系的建立 ［J］. 江蘇農(nóng)業(yè)科學(xué)， 46（24）： 47-50.］

CHEN YY， LU HQ， JIANG KX， et al.， 2023. The flavonoid biosynthesis and regulation in Brassica napus： A review ［J］. International Journal of "Molecular Sciences， 24： 357.

DAI MJ， 2022. Molecular cloning and functional analysis of CHIL and F3H in flavonoid biosynthesis pathway in mulberry ［D］. Zhenjiang： Jiangsu University of Science and Technology. "［戴明潔， 2022. 桑樹類黃酮合成途徑關(guān)鍵基因CHIL和F3H的分子克隆及功能分析 ［D］. 鎮(zhèn)江： 江蘇科技大學(xué).］

FANG L， MARTRE P， JIN K， et al.， 2022. Neglecting acclimation of photosynthesis under drought can cause significant errors in predicting leaf photosynthesis in wheat ［J］. Global Change Biology， 29： 505-521.

GANBARI MSR， TOHIDFAR M， AHMADIKHAH A， 2022. Gene profile analysis and molecular-physiological evaluation of tomato genotypes under drought stress ［J］. Journal of Applied Biotechnology Reports， 9（3）： 74.

GANETEG U， KIMMEK F， JANSSON S， 2004. Lhca5-an LHC-type protein associated with photosystem I ［J］. Plant Molecular Biology， 54（5）： 641-651.

GHANBARI MSR，TOHIDFAR M， AHMADIKHAH A， 2022. Gene profile analysis and molecular-physiological evaluation of tomato genotypes under drought stress ［J］. Journal of "Applied Biotechnology Reports， 9（3）： 747-62.

GUO LA， MO RY， TAN J， et al.， 2021. Analysis on differential expressed genes in the transcriptome of Atractylodes macrocephalafrom different producing areas ［J］. Chinese Medicinal Materials， 44（12）： 2787-2792. "［郭連安， 莫讓瑜， 譚均， 等， 2021. 不同產(chǎn)地白術(shù)轉(zhuǎn)錄組差異表達基因分析 ［J］. 中藥材， 44（12）： 2787-2792.］

GUO LA， TAN J， DENG XS， et al.， 2023. Integrated analysis of metabolome and transcriptome reveals key candidate genes involved in flavonoid biosynthesis in Pinellia ternata under heat stress ［J］. Journal of Plant Research， 136（3）： 359-369.

HOU YT， LIU B， 2010. Anemone shikokiana（Makino） Makino ［J］. Encyclopedia of "Knowledge （12）： 40. "［侯元同， 劉冰， 2010. 山東銀蓮花 ［J］. 百科知識 （12）： 40.］

JIANG Q， 2019. Genomic analysis and response factors under abiotic stress in Oenanthe javanica ［D］. Nanjing： Nanjing Agricultral University. "［蔣倩， 2019. 水芹基因組研究及非生物脅迫下的響應(yīng)因子分析 ［D］. 南京： 南京農(nóng)業(yè)大學(xué).］

LI H， FENG YH， TANG SS， et al.， 2023. Development and validity evaluation of Liquidambar formosana EST-SSR primers based on transcriptome sequencing ［J］. Guihaia， 43（2）： 327-335. "［李輝， 馮源恒， 唐生森， 等， 2023. 基于轉(zhuǎn)錄組測序的楓香EST-SSR引物開發(fā)及有效性評價 ［J］. 廣西植物， 43（2）： 327-335.］

LIU Q， 2014. Research on genetic diversity of Anemone shikokiana ［D］. Jinan： Shandong Normal University. "［劉瓊， 2014. 山東銀蓮花的遺傳多樣性研究 ［D］. 濟南： 山東師范大學(xué).］

LOU YF， GAN R， ZHU CL， et al.， 2022. Development and application of EST-SSR markers based on transcriptome of Phyllostachys vivax McClure f. aureocaulis N. X. Ma ［J］. Plant Science Journal， 40（3）： 355-364. "［婁永峰， 甘然， 朱成磊， 等， 2022. 黃稈烏哺雞竹轉(zhuǎn)錄組EST-SSR分子標(biāo)記開發(fā)與應(yīng)用 ［J］. 植物科學(xué)學(xué)報， 40（3）： 355-364.］

LU YH， YU FB， GUO JH， et al.， 2021. Analysis of transcriptomic differences between winged and wingless bird cherry-oat aphid Rhopalosiphum padi on wheat ［J］. Journal of Plant Protection， 48（6）： 1303-1309. "［魯艷輝， 虞飛博， 郭嘉雯， 等， 2021. 有翅型和無翅型小麥禾谷縊管蚜轉(zhuǎn)錄組差異分析 ［J］. 植物保護學(xué)報， 48（6）： 1303-1309.］

LUCINSK R， JACKOWSKI G， 2006. The structure， functions and degradation of pigment-binding proteins of Photosystem II ［J］. Acta Biochimica Polonica， 53 （4）： 693-708.

LU TT， LU GJ， FAN DL， et al.， 2010. Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-Seq ［J］. Genome Research， 20（9）： 1238-1249.

MANFRED G， BRIAN J， MORAN Y， et al.， 2011. Trinity： reconstructing a full-length transcriptome without a genome from RNA-Seq data ［J］. Nature Biotechnology， 29（7）： 644.

MOZZO M， MANTELLI M， PASSARINI F， et al.， 2010. Functional analysis of Photosystem I light-harvesting complexes （Lhca） gene products of Chlamydomonas reinhardtii ［J］. Biochimic Biophysica Acta， 1797（2）： 212-221.

PANG YJ， LI LX， BIAN FH， 2020. Photosynthetic and ultrastructural adaptability of Anemone shikokiana leaves to heterogeneous habitats ［J］. Brazilian Journal of "Botany， 43： 979-988.

PANG YJ， SUN W， WANG XS， et al.， 2022. Influence of soil heterogeneity on root structure of Anemone shikokiana ［J］. Journal of "Yantai University， 35（1）： 28-34. "［逄玉娟， 孫維， 王行松， 等， 2022. 土壤異質(zhì)性對山東銀蓮花根結(jié)構(gòu)的影響 ［J］. 煙臺大學(xué)學(xué)報， 35（1）： 28-34.］

ROZEN S， SKALETSKY H， 2000. Primer3 on the WWW for general users and for biologist programmers ［J］. Bioinformatics Methods and Protocols， 132： 365-86.

SAPETA H， YOKONO M， TAKABAYASHI A， et al.， 2023. Reversible down-regulation of photosystems I and II leads to fast photosynthesis recovery after long-term drought in Jatropha curcas ［J］. Journal of "Experimental Botony， 74（1）： 336-351.

SHEN N， WANG TF， GAN Q， et al.， 2022. Plant flavonoids： Classification， distribution， biosynthesis， and antioxidant activity ［J］. Food Chemistry， 383： 132531.

SUN HN， LIU DY， LI MY， et al.， 2024. Transcriptome data analysis of Coreopsis tinctoria based on high-throughput sequencing ［J］. Molecular Plant Breeding， 22（20）： 6683-6697. "［孫浩男， 劉冬云， 李明陽， 等， 2024. 基于高通量測序的兩色金雞菊轉(zhuǎn)錄組數(shù)據(jù)分析 ［J］. 分子植物育種， 22（20）： 6683-6697.］

TEMNYKH S， DECLERCK G， LUKASHOVA A， 2001. Computational and experimental analysis of microsatellites in rice （Oryza sativa L.）： Frequency， length variation， transposon associations， and genetic marker potential ［J］. Genome Research， 11（8）： 1441-1452.

TIAN HZ， GUO H， TIAN H， et al.， 2023. Transcriptome characteristics analysis and EST-SSR marker development of capsicum annuum under waterlogging stress ［J］. Guihaia， 43（11）： 2052-2064. ［田懷志， 郭豪， 田浩， 等， 2023. 水澇脅迫下辣椒轉(zhuǎn)錄組特征分析及 EST-SSR 標(biāo)記開發(fā) ［J］. 廣西植物， 43（11）： 2052-2064.］

WANG YF， ZHU RJ， CHENG XP， 2023.Research on phenotypic plasticity of different diameter class Bombax ceiba in heterogeneous habitats ［J］. Ecological Science， 42（3）： 61-66. "［王妍方， 朱潤軍， 程希平， 2023. 異質(zhì)生境下不同徑級木棉的表型可塑性研究 ［J］. 生態(tài)科學(xué)， 42（3）： 61-66.］

WANG YH， 2020. Identification and expression profiling of Lhc genes in Ramie ［D］. Wuhan： Huazhong Agricultural University. "［王云鶴， 2020. 苧麻Lhc基因的鑒定和表達模式研究 ［D］. 武漢： 華中農(nóng)業(yè)大學(xué).］

WANG Z， PANG YJ， LIU CL， et al.， 2014. Distribution pattern and analysis of influencing factors of rare plant Anemone shikokiana （Makino） Makino ［J］. Bulletin of Botanical Research， 34（4）： 440-445. "［王鷙， 逄玉娟， 劉傳林， 等， 2014. 稀有植物山東銀蓮花（Anemone shikokiana （Makino） Makino）的分布格局及影響因子分析 ［J］. 植物研究， 34（4）： 440-445.］

WEI CB， GU W， TIAN R， et al.， 2022. Transcriptomic data analysis of Alisma orientale tubers grew in different ecological environments by high-throughput sequencing ［J］. Modern Chinese Medicine， 24（12）： 2342-2349. "［韋陳彬， 谷巍， 田榮， 等， 2022. 基于高通量測序的不同生態(tài)環(huán)境東方澤瀉塊莖轉(zhuǎn)錄組分析 ［J］. 中國現(xiàn)代中藥， 24（12）： 2342-2349.］

WU RN， XU B， SHI FL， 2022. Analysis of SSR characteristics for Medicago ruthenica ‘Zhilixing’ based on transcriptome data ［J］. Journal of "Northwest A amp; F University （Natural Science Edition）， 50（5）： 1-8. "［烏日娜， 徐舶， 石鳳翎， 2022. 基于轉(zhuǎn)錄組數(shù)據(jù)的直立型扁蓿豆SSR序列特征分析 ［J］. 西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版）， 50（5）： 1-8.］

XIANG NY， LU BJ， YUAN T， et al.， 2023. De Novo transcriptome assembly and EST-SSR marker development and application in Chrysosplenium macrophyllum ［J］. Genesis， 14： 279.

XIAO YZ， HAN SM， QIN Z， et al.， 2020. Analysis of transcriptome sequencing and related genes of flavonoids biosynthesis from Polygonatum kingianum ［J］. Journal of "Henan Agricultural University， 54（6）： 931-940. "［肖韻錚， 韓世明， 秦昭， 等， 2020. 滇黃精轉(zhuǎn)錄組測序及類黃酮合成相關(guān)基因的分析 ［J］. 河南農(nóng)業(yè)大學(xué)學(xué)報， 54（6）： 931-940.］

XIONG S， TIAN N， LONG JH， et al.， 2016. Molecular cloning and characterization of a flavanone 3-hydroxylase gene from Artemisia annua L. ［J］. Plant Physiology and Biochemistry， 105（5）： 29-36.

YAN MM， WEI GC， PAN XH， et al.， 2008. A method suitable for extracting genomic DNA from animal and plant —— modifed CTAB Method agricultural ［J］. Journal of "Anhui Agricultural Science， 36（20）： 8488+8558. "［閆苗苗， 魏光成， 潘效紅， 等， 2008. 一種適用于動物與植物總DNA提取的方法——改良CTAB法 ［J］. 安徽農(nóng)業(yè)科學(xué)， 36（20）： 8488+8558.］

YANG X， YANG N， YUAN QH， et al.， 2021. Development and application of EST-SSR molecular markers in Pinus bungeana ［J］. Journal of "Beijing Forestry University， 43（7）： 1-11. "［楊雄， 楊寧， 袁啟華， 等， 2021. 白皮松 EST-SSR 分子標(biāo)記的開發(fā)及應(yīng)用 ［J］. 北京林業(yè)大學(xué)學(xué)報， 43（7）： 1-11.］

YU WY， GAO Y， PANG YJ， et al.， 2019. Response of leaf morphology and structure of Anemone shikokiana to heterogeneous habitats and altitude changes ［J］. Acta Ecologica Sinica， 39 （12）： 4413-4420. "［于文英， 高燕， 逄玉娟， 等， 2019. 山東銀蓮花葉片形態(tài)結(jié)構(gòu)對異質(zhì)生境和海拔變化的響應(yīng) ［J］. 生態(tài)學(xué)報， 39（12）： 4413-4420.］

YUE XY， YUE W， YAN K， 2023. Adaptation of functional traits of Phragmites australis to different water Habitats in the Yellow River Delta ［J］. Journal of Northeast Forestry University， 51（3）： 18-21. "［岳喜元， 岳偉， 嚴寬， 2023. 黃河三角洲蘆葦功能性狀對不同水分生境的適應(yīng)性 ［J］. 東北林業(yè)大學(xué)學(xué)報， 51（3）： 18-21.］

YUE YH， YAN LJ， HUANG X， et al.， 2022. Development EST-SSR markers based on transcriptome sequences of Tilia miqueliana Maxim ［J］. Molecular Plant Breeding： 1-18 ［2023-10-20］. http：//kns.cnki.net/kcms/detail/46.1068.S.20220715.1525.008.html. ［岳遠灝， 嚴靈君， 黃犀， 等， 2022. 基于全長轉(zhuǎn)錄組測序的南京椴EST-SSR標(biāo)記開發(fā) ［J］. 分子植物育種： 1-8 ［2023-10-20］. http：//kns.cnki.net/kcms/detail/46.1068.S.20220715.1525.008.html. ］

ZHANG CY， LONG Y， FENG J， et al.， 2007. Transcriptional regulation of plant genes and its significance in biology ［J］. Hereditas （7）： 793-799. ""［張椿雨， 龍艷， 馮吉， 等， 2007. 植物基因在轉(zhuǎn)錄水平上的調(diào)控及其生物學(xué)意義 ［J］. 遺傳 （7）： 793-799.］

ZHENG Y， JIN GL， WU WR， 2012. Relationship between sequence completeness and polymorphism of microsatellites ［J］. Genomics and Applied Biology， 31（6）： 587-591. "［鄭燕， 金谷雷， 吳為人， 2012. 微衛(wèi)星序列完整性與多態(tài)性的關(guān)系 ［J］. 基因組學(xué)與應(yīng)用生物學(xué)， 31（6）： 587-591.］

ZHOU CM， JIANG LP， DONG GQ， et al.， 2022. Analysis of differential expression genes of Chrysanthemum indicum var. aromaticum associated with different altitude habitats using RNA-seq technology ［J］. Modernization of Traditional Chinese Medicine and Materia Medica-World Science and Technology， 24（4）： 1356-1368. "［周春苗， 蔣莉萍， 董剛強， 等， 2022. 基于轉(zhuǎn)錄組測序研究不同海拔生境神農(nóng)香菊的基因表達差異 ［J］. 世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化， 24（4）： 1356-1368.］

ZHUANG WB， LI YH， SHU XC，et al.， 2023. The classification， molecular structure and biological biosynthesis of flavonoids， and their roles in biotic and abiotic stresses ［J］. Molecules， 28（8）： 3599.

ZOU FL， SHOU SY， YE WZ， et al.， 2004. Advances in the research on flavonoid biosynthesis and plant stress response ［J］. Chinese Journal of "Cell Biology， 26（1）： 39-44. "［鄒鳳蓮， 壽森炎， 葉紈芝， 等， 2004. 類黃酮化合物在植物脅迫反應(yīng)中作用的研究進展 ［J］. 細胞生物學(xué)學(xué)報， 26（1）： 39-44.］

（責(zé)任編輯 李 莉 王登惠）