孫金魁,許厚強,石鵬飛,阮 涌
(1.貴州大學(xué) 高原山地動物遺傳育種與繁殖教育部重點實驗室/貴州省動物遺傳育種與繁殖重點實驗室,貴陽 550025; 2.貴州大學(xué)動物科學(xué)學(xué)院,貴陽550025)
關(guān)嶺牛具有體質(zhì)強壯、抗病力強等優(yōu)良特性,其自身的獨特基因是生物多樣性的重要組成部分[1-2]。與其它黃牛品種相比,關(guān)嶺牛產(chǎn)肉率低、生長遲緩,限制了高原山地畜牧業(yè)的發(fā)展。肌原性細胞的增殖分化是影響肌肉生長發(fā)育的關(guān)鍵,肌細胞增強因子2(myocyte enhancer factor2, MEF2)在肌肉細胞中廣泛存在,可將該基因視作肌源性基因表達的主要調(diào)節(jié)物[3-4]。
研究表明,MEF2通過調(diào)節(jié)肌肉特異基因的表達影響骨骼肌發(fā)育分化進程,MEF2屬于轉(zhuǎn)錄調(diào)節(jié)因子MADS-Box家族,該家族包括微小染色體維持蛋白(MCM1)、AGAMOUS、DEFICIENS和血清反應(yīng)因子(SRF)等成員[5]。MEF2蛋白具備DNA結(jié)合活性,擁有結(jié)合肌肉肌酸激酶(MCK)基因啟動子中富含A/T的DNA序列的能力[6]。前期研究中,MEF2基因家族在脊椎動物中包括MEF2A、MEF2B、MEF2C和MEF2D四個成員[7],該家族成員結(jié)構(gòu)相似,皆含有中心MEF2結(jié)構(gòu)域、N端結(jié)合DNA的MADS區(qū)域和C端轉(zhuǎn)錄激活區(qū)域;C末端區(qū)域主要由選擇性拼接而成,MEF2因子可發(fā)生同源二聚化和多種異源二聚化,以具體的方式參與細胞生命周期中許多重要的轉(zhuǎn)錄過程[8]。MEF2主要通過MADS-box結(jié)構(gòu)域和MEF2結(jié)構(gòu)域與其他轉(zhuǎn)錄輔因子的互作來調(diào)控其下游靶基因的轉(zhuǎn)錄活性,具備調(diào)控細胞凋亡(包括神經(jīng)元、肌肉、血管、淋巴細胞)和骨骼在內(nèi)的多種組織的分化、形態(tài)發(fā)生和維持的作用[9-10]。
MEF2A基因是具有bHLH(basic helix-loop-helix)結(jié)構(gòu)域的轉(zhuǎn)錄因子,包含MADS-box和HJURP-C兩個保守結(jié)構(gòu)域,兩個結(jié)構(gòu)域可識別肌肉調(diào)控基因啟動區(qū)中E-box元件(CANNTG),進而激活下游基因的特異性表達[11-12]。Edmondson等[13]研究表明,在肌再生過程中,MEF2A是最早可檢測到的基因之一,在骨骼肌中高度表達,對肌纖維聚合形成、骨骼肌干細胞和成肌細胞增殖分化具有不可或缺的作用[14-15]。Wu等[16]研究發(fā)現(xiàn),MEF2A具備承接結(jié)構(gòu)蛋白與調(diào)控因子的功能,該基因異位表達可結(jié)合Capn3啟動子區(qū)域的預(yù)測位點,二者互作可影響L6成肌細胞的分化模型。MEF2A具備獨特的轉(zhuǎn)錄活性,能夠維持出生后心肌細胞結(jié)構(gòu)的完整性,Clark和Naya[17]研究證實,MEF2A可通過調(diào)節(jié)Gtl2/Dio3位點,進而刺激轉(zhuǎn)錄輔因子CITED2來促進心肌細胞的增殖;此外,MEF2A可通過激活SIRT1的表達來保護血管內(nèi)皮細胞[18]。Blixt等[19]研究報道,MEF2A基因在小鼠破骨細胞模型中具有性別特異性調(diào)節(jié)作用。因此,MEF2A對于肌肉組織生長發(fā)育的調(diào)控作用因組織類型而異,對骨骼肌成肌細胞增殖分裂具有重要作用。
目前,有關(guān)MEF2A基因在關(guān)嶺牛(Bostaurus)成肌細胞上的研究鮮有報道,影響成肌細胞的機理尚不明確。因此,本試驗以關(guān)嶺牛成肌細胞作為研究對象,構(gòu)建篩選MEF2A基因的最佳干擾載體,轉(zhuǎn)染載體至細胞后,qRT-PCR法檢測成肌細胞中肌生成相關(guān)基因的表達;檢測細胞周期與活力特性。為進一步探索MEF2A基因?qū)﹃P(guān)嶺牛肌肉發(fā)育的分子效應(yīng)及調(diào)控機制提供數(shù)據(jù)參考。
Lipofectamine TM3000CD Reagent脂質(zhì)體轉(zhuǎn)染試劑盒購自Thermo Fisher Scientific(美國);DMEM/F12培養(yǎng)基、Gibco澳洲胎牛血清均購自重慶奧怡生物技術(shù)有限公司;PBS溶液、青鏈霉素、Opti-MEMTM均購自廣西卓一生物技術(shù)有限公司;Trizol試劑購自英駿生物技術(shù)有限公司(美國);SYBR Green qPCR Master Mix、逆轉(zhuǎn)錄試劑盒均購自Genstar公司(北京)。超微量紫外分光光度計購自美國Thermo Fisher公司;酶標儀購自美國Bio Tek公司;實時熒光定量PCR儀(型號為CFX96 Touch)、多功能梯度擴增PCR儀(型號為ABIVeriti-TM)、凝膠成像系統(tǒng)(型號為Uni-versal Hood Ⅱ)均購自美國BIO-RAD有限公司。
試驗動物來自貴州省安順市西秀區(qū)幺鋪鎮(zhèn)屠宰場,選擇健康的3日齡犢牛,采集犢牛背最長肌,生理鹽水(0.9% NaCl)沖洗后放入含雙抗的PBS溶液中保存帶回實驗室。成肌細胞分離培養(yǎng)與鑒定參照宋林錦等[20]的操作方法。率先清理筋膜與雜質(zhì),剪碎肌肉組織后加入Ⅱ型膠原酶,放置37 ℃恒溫搖床消化1 h,加入培養(yǎng)基混勻,依次通過滅菌紗布、200目和400目細胞篩,濾液以8 000 r·min-1離心8 min,棄上清液后加入10% DMEM/F-12培養(yǎng)基重懸細胞,放置于37 ℃,5% CO2的細胞培養(yǎng)箱中。利用免疫熒光法鑒定成肌細胞純度,將細胞懸液接種在置有玻璃片的24孔板中,細胞培養(yǎng)箱中培養(yǎng)2 h;細胞爬片后,加入4%多聚甲醛于4 ℃條件下固定35 min,PBS洗滌3次后置于培養(yǎng)皿支撐物上。吸取60 μL破膜封閉液滴于防水膜上,2 h破膜封閉,取60 μL一抗(α-actin)于防水膜上,蓋上玻片于4 ℃放置;避光孵育二抗(goat anti-rabbit lgG)2 h后,PBS洗滌3次,每次4 min,DAPI染色5 min,PBS洗滌3次,每次4 min。加1滴熒光封片劑于玻片上,將有細胞的一面蓋上,用于后續(xù)拍照儲存。
參考NCBI數(shù)據(jù)庫中的MEF2A基因編碼區(qū)序列(GenBank登錄號:NM_001083638.2)與短發(fā)夾RNA(short hairpin RNA, shRNA)設(shè)計原則,設(shè)計MEF2A基因的shRNA片段引物,共設(shè)置4個干擾組靶序列和1個NC組序列,分別命名為shRNA-MEF2A-1、shRNA-MEF2A-2、shRNA-MEF2A-3、shRNA-MEF2A-4與shRNA-NC。于序列兩端添加EcoR Ⅰ與SalⅠ酶切位點,序列中段插入Loop (TCAAGAG)結(jié)構(gòu),序列最后加上6個T,可作為RNA聚合酶Ⅲ的轉(zhuǎn)錄終止子發(fā)揮作用。由上海吉瑪制藥技術(shù)有限公司合成片段引物,引物序列信息見表1。
表1 shRNA信息
復(fù)蘇空載體質(zhì)粒,利用EcoRⅠ與SalⅠ同時進行酶切,取1.5 μL pGPU6-GFP-Neo、3 μL10×Green Buffer、1.5 μLEcoRⅠ、1.5 μLSalⅠ,加ddH2O補足30 μL,37 ℃條件下酶切2 h,酶切產(chǎn)物用1.2%的瓊脂糖凝膠電泳檢測并切膠回收;將shRNA片段連接到pGPU6-GFP-Neo載體,16 ℃金屬浴過夜連接。連接產(chǎn)物轉(zhuǎn)化到DH5α感受態(tài)細胞中,挑選單克隆菌落擴大培養(yǎng),37 ℃,200 r·min-1振蕩培養(yǎng)16 h,將菌液送至上海生工生物工程股份有限公司進行測序驗證。
成肌細胞接種于6孔細胞培養(yǎng)板中,待細胞生長至80%~90%的融合度時,利用Lipofectamine 3000進行轉(zhuǎn)染,分別將重組載體pGPU6-GFP-Neo-NC、pGPU6-GFP-Neo-MEF2A-1、pGPU6-GFP-Neo-MEF2A-2、pGPU6-GFP-Neo-MEF2A-3、pGPU6-GFP-Neo-MEF2A-4轉(zhuǎn)染至成肌細胞內(nèi),各孔轉(zhuǎn)染質(zhì)粒2 500 ng,空白組不轉(zhuǎn)染質(zhì)粒,各轉(zhuǎn)染組設(shè)置3個生物學(xué)重復(fù),放置于37 ℃,5% CO2的培養(yǎng)箱中。48 h后將細胞置于熒光倒置顯微鏡下觀察生長狀態(tài)并檢測MEF2A基因shRNA干擾效率。
Trizol法提取細胞總RNA,反轉(zhuǎn)錄得到cDNA作為模板。依據(jù)NCBI中公布的基因序列,利用Primer 5.0設(shè)計牛肌細胞增強因子家族基因以及細胞生長相關(guān)基因CDK2、BCL2與CCNA2的熒光引物(表2)。qRT-PCR檢測肌細胞增強因子家族基因(MEF2A、MEF2B、MEF2C、MEF2D)以及BCL2、CDK2與CCNA2基因在最佳干擾組和對照組中的表達水平,GAPDH為內(nèi)參基因。PCR反應(yīng)體系10 μL:SYBR Green Mix 5 μL,上、下游引物各0.5 μL,cDNA 0.5 μL,ddH2O 3.5 μL。每個樣品做3個生物學(xué)重復(fù),反應(yīng)條件為95 ℃預(yù)變性30 s;95 ℃變性5 s,退火5 s,95 ℃延伸5 s,共39個循環(huán)。
表2 引物信息
將最佳干擾組與對照組載體分別轉(zhuǎn)染至成肌細胞,用4 ℃預(yù)冷的PBS洗滌,加入1 mL 0.25%不含EDTA的胰酶進行消化,1 000 r·min-1離心5 min,棄上清液,加入預(yù)冷的1.5 mL的70%乙醇吹打細胞,4 ℃儲存過夜。離心乙醇固定過的細胞,棄上清液;PBS洗滌細胞3次;加入100 μL RNase A重懸細胞,37 ℃水浴鍋中放置30 min,加入500 μL PI染色液,4 ℃避光反應(yīng)30 min;利用流式細胞儀進行檢測。
最佳干擾組與對照組載體分別轉(zhuǎn)染至成肌細胞后,消化計數(shù)并接種于96孔細胞板中,各組別設(shè)置3個重復(fù),在0、6、12、24、48和72 h時間段內(nèi)進行檢測,各孔添加10 μL的CCK8試劑,放置于37 ℃,5% CO2培養(yǎng)箱中培養(yǎng)2 h,利用酶標儀(SynerGyH1, BioTek)檢測各組別在450 nm波長處的吸光度。
利用在線軟件ProtParam(https://web.expasy.org/protparam/)、PSORT II Preadict (https://psort.hgc.jp/form2.html)與STRING(https://string-db.org/)分別對關(guān)嶺牛MEF2A的蛋白理化特性、亞細胞結(jié)構(gòu)定位與蛋白網(wǎng)絡(luò)圖譜進行預(yù)測分析。
采用2-ΔΔCt法計算基因表達量,各組別設(shè)置3個生物學(xué)重復(fù),數(shù)據(jù)均用“平均值±標準誤”表示,采用獨立樣本t檢驗和單因素方差分析進行數(shù)據(jù)處理(SPSS 25.0),多重比較采用Duncan’s多重極差檢驗法,運用Graphpad 6.0作圖。
如圖1所示,固定細胞后,α-actin經(jīng)異硫氰酸熒光素(fluorescein isothiocyanate, FITC)標記后在細胞質(zhì)內(nèi)陽性表達為綠色發(fā)光狀態(tài)。細胞核為橢圓形狀態(tài),經(jīng)4′, 6-二脒基-2-苯基吲哚(DAPI)染色后發(fā)藍色光。細胞核與細胞質(zhì)染色重合良好,表明成功培養(yǎng)鑒定關(guān)嶺牛成肌細胞,進一步鑒定細胞純度達到90%左右,可開展后續(xù)試驗。
圖1 成肌細胞免疫熒光鑒定(200×)
將酶切純化回收的pGPU6-GFP-Neo載體與退火后的shRNA置于16 ℃恒溫金屬浴中過夜連接,連接產(chǎn)物轉(zhuǎn)化到DH5α感受態(tài)細胞中,挑選單克隆菌落進行搖菌培養(yǎng),篩選陽性菌液送至上海生工生物工程股份有限公司測序驗證。測序結(jié)果如圖2所示,插入的重組質(zhì)粒片段與設(shè)計合成的序列一致,說明成功構(gòu)建了重組質(zhì)粒。
圖2 重組質(zhì)粒測序分析結(jié)果
將成肌細胞接種到6孔板中,細胞融合度達到90%左右時,將各干擾組與陰性對照組轉(zhuǎn)染至成肌細胞,放置于細胞培養(yǎng)箱中48 h后用熒光倒置顯微鏡觀察各組別綠色熒光蛋白的發(fā)光情況。結(jié)果如圖3所示,在熒光倒置顯微鏡下的各干擾組、陰性對照組均有清晰明亮的綠色熒光,其中空白組不發(fā)光,說明轉(zhuǎn)染成功。
A. shRNA-MEF2A-1;B. shRNA-MEF2A-2;C. shRNA-MEF2A-3;D. shRNA-MEF2A-4;E. shRNA-NC;F. 空白對照
將各重組質(zhì)粒分別轉(zhuǎn)染至牛成肌細胞48 h,提取細胞RNA并逆轉(zhuǎn)錄為cDNA,qRT-PCR法檢測MEF2A基因在不同組別中的表達量特性,篩選干擾效果最佳的試驗組別。如圖4所示,shRNA-MEF2A-1、shRNA-MEF2A-2與shRNA-MEF2A-3和shRNA-MEF2A-4干擾效率依次為74.45%、58.67%、94.66%與62.98%。結(jié)果表明,4個干擾載體均能極顯著下調(diào)MEF2A基因表達量(P<0.01),選用干擾效率最高的shRNA-MEF2A-3干擾載體開展后續(xù)試驗。
*. 表示差異顯著(P<0.05);**. 表示差異極顯著(P<0.01),下同
分別提取轉(zhuǎn)染shRNA-MEF2A-3最佳干擾組和shRNA-NC對照組的細胞總RNA,對肌生成基因及周期與凋亡相關(guān)基因表達量進行qRT-PCR檢測。結(jié)果如圖5所示,干擾MEF2A基因表達后,肌細胞增強因子MEF2B、MEF2C與MEF2D基因表達量均極顯著上調(diào)(P<0.01);細胞周期因子CDK2表達量顯著下調(diào)(P<0.05),CCNA2表達量極顯著下調(diào)(P<0.01);抗凋亡因子BCL2表達量顯著下調(diào)(P<0.05)。推斷MEF2A基因在關(guān)嶺牛成肌細胞中作用顯著,與家族其他成員具有相互調(diào)控關(guān)系;干擾MEF2A可延長關(guān)嶺牛成肌細胞周期,對其增殖具有抑制作用。
圖5 干擾MEF2A對成肌細胞基因表達的影響
利用流式細胞儀檢測轉(zhuǎn)染shRNA-MEF2A-3干擾組和shRNA-NC對照組的細胞周期特性。結(jié)果如圖6所示,與shRNA-NC相比,轉(zhuǎn)染shRNA-MEF2A-3的成肌細胞數(shù)量在G1期顯著增多(P<0.05),S期與G2/M期細胞數(shù)量極顯著減少(P<0.01)。研究結(jié)果表明,干擾MEF2A后,細胞周期G1至S期的進程被阻滯,細胞增殖受到抑制。
A. 干擾組細胞周期;B. NC組細胞周期;C. 細胞周期比例
分別在0、6、12、24、48和72 h檢測shRNA-MEF2A-3干擾組和shRNA-NC對照組的細胞生長增殖特性。CCK8檢測結(jié)果如圖7所示,自6 h后,shRNA-MEF2A-3干擾組和shRNA-NC對照組在5個不同時間段內(nèi)的細胞生長情況均達到極顯著差異水平(P<0.01)。
圖7 干擾MEF2A對成肌細胞增殖的影響
在線軟件預(yù)測MEF2A(No.:NM_001083638.2)基因的mRNA序列為5 568 bp,編碼492個氨基酸。MEF2A蛋白分子式為C2282H3649N659O737S21,相對分子質(zhì)量約為52 782.31 u;絲氨酸(13.8%)在氨基酸中組成含量最高,色氨酸(0.4%)組成含量最低;帶正電荷的氨基酸殘基數(shù)量(精氨酸+賴氨酸)為45個;帶負電荷的氨基酸殘基數(shù)量(天冬氨酸+谷氨酸)為39個。亞細胞定位顯示,MEF2A主要位于細胞核(60.9%)。STRING預(yù)測牛MEF2A蛋白互作網(wǎng)絡(luò)圖譜,結(jié)果如圖8所示,包含組蛋白去乙酰化酶(HDAC)家族蛋白成員HDAC4、HDAC5與HDAC7,參與核心組蛋白(H2A、H2B、H3和H4)N末端部分賴氨酸殘基的去乙酰化[21];絲裂原蛋白活化激酶(MAPK)家族蛋白成員MAPK11、MAPK14,可通過激活環(huán)中TGY基序的雙重磷酸化發(fā)揮功效[22],調(diào)節(jié)細胞分化與生長凋亡,在細胞反應(yīng)級聯(lián)中起重要作用,以及影響心肌發(fā)育和成肌細胞分化的MEF2C、EP300、MEF2D、MYOD1、NKX2-5等相關(guān)蛋白。
圖8 牛MEF2A蛋白互作網(wǎng)絡(luò)預(yù)測分析
MEF2本身不具備肌源性活性,當與bHLH(basic helix-loop-helix)蛋白家族互作時,可驅(qū)動并放大肌源性分化程序[23]。MEF2蛋白N端具備結(jié)合DNA功能的MADS區(qū)域和中心MEF2結(jié)構(gòu)域,C端為轉(zhuǎn)錄激活區(qū)域,該區(qū)域包含大量的核定位序列和磷酸化位點,在MEF2基因家族成員及不同物種間的保守性較低且序列變化多樣,與基因轉(zhuǎn)錄活性關(guān)系密切,可進一步調(diào)控MEF2蛋白活性[24-25]。MEF2主要通過CaMK-HDACs(Ca2+-鈣調(diào)素依賴的蛋白激酶-組蛋白去乙?;?、Calcineurin(鈣調(diào)磷酸酶)和MAPK(絲裂原活化蛋白激酶)3種信號通路發(fā)揮作用[26]。MEF2基因大多通過MADS-box和MEF2結(jié)構(gòu)域與其他轉(zhuǎn)錄輔因子協(xié)同作用,家族成員與許多參與肌生成的基因相結(jié)合,以啟動或加強肌源性基因的表達[27-28],通過控制下游靶基因的轉(zhuǎn)錄活性來調(diào)控細胞的分裂、分化和凋亡過程。目前已證實,牛MEF2A基因被定位于21號染色體[29],是一種廣譜表達基因,在肌肉發(fā)育過程中,該基因表達具有發(fā)育時空性與組織特異性。MEF2A基因的多態(tài)性對畜禽屠宰性能存在一定的影響,且該基因為影響畜禽屠宰性狀的候選基因之一[30]。Liu等[31]研究發(fā)現(xiàn),同時敲除MEF2A、MEF2C和MEF2D后,小鼠肌損傷修復(fù)過程受到顯著抑制,成肌細胞增殖受到嚴重阻礙。而本試驗在抑制MEF2A基因表達后,關(guān)嶺牛成肌細胞生長活性下降,細胞增殖同樣受到嚴重阻礙。同時,有研究表明MEF2A、MEF2C表達水平在蛙類動物冬眠期具有重要作用[32]。
MEF2A作為一種較為保守的轉(zhuǎn)錄因子,廣泛分布于各種組織或器官中,在包括決定細胞命運、遷移和形狀在內(nèi)的多種生物學(xué)過程中發(fā)揮關(guān)鍵作用[33]。MEF2A是肌生成的關(guān)鍵調(diào)節(jié)因子和肌再生的重要激活劑,也可參與脂質(zhì)代謝途徑[34]。Wang等[35]研究發(fā)現(xiàn),MEF2A下游的MEG3-DIO3 miRNA在牛成肌細胞可調(diào)節(jié)PP2A信號;MEF2A-MEG3/DIO3-PP2A軸在肌再生過程中發(fā)揮重要作用,為骨骼肌的發(fā)育機制研究提供了新的線索。Hennebry等[36]研究報道,MSTN可通過調(diào)控MEF2與MyoD基因家族表達進而調(diào)節(jié)骨骼肌纖維類型組成;在轉(zhuǎn)基因小鼠中表達過活性形式的MEF2可促進慢肌纖維形成,但MEF2的失活會導(dǎo)致慢肌纖維形成的喪失[37]。Schiaffino與Serrano[38]研究發(fā)現(xiàn),細胞內(nèi)Ca2+濃度較低時,HDAC、HDAC5、HDAC7與MEF2蛋白N末端的MADS/MEF2結(jié)構(gòu)域相互作用而抑制MEF2的活性;當Ca2+濃度升高時,Ca2+與鈣調(diào)蛋白結(jié)合并激活CaMK,其中CaMKIV發(fā)揮磷酸化作用使HDACs失去活性,MEF2從HDACs上脫離,與本研究中預(yù)測的MEF2A蛋白網(wǎng)絡(luò)圖譜結(jié)果相吻合。此外,有研究表明MEF2A在抗炎癥和治療癲癇等方面有著重要作用[39-41]。
本試驗初步分析了關(guān)嶺牛MEF2A基因的蛋白理化性質(zhì),利用qRT-PCR法檢測發(fā)現(xiàn)抑制MEF2A基因表達后,肌細胞增強因子MEF2B、MEF2C與MEF2D的表達量均極顯著高于shRNA-NC對照組,進一步推測該家族成員具有協(xié)同作用的特性,且MEF2A基因在成肌細胞中的作用比家族其它成員更加顯著。CDK2作為細胞周期蛋白依賴性激酶家族中的一員,與細胞周期蛋白CCNA2相互作用可調(diào)控細胞周期進程[42]。MEF2A基因被抑制后,相較于shRNA-NC對照組,成肌細胞中CDK2與CCNA2基因表達量皆降低;流式細胞結(jié)果顯示,G1期細胞數(shù)量顯著增多,S期細胞數(shù)量極顯著減少,表明細胞周期進程受到阻滯。同時發(fā)現(xiàn)干擾MEF2A基因可導(dǎo)致細胞抗凋亡因子BCL2基因表達量下降,進一步研究發(fā)現(xiàn)轉(zhuǎn)染干擾組的成肌細胞自6 h后的細胞增殖能力相較于對照組下降明顯。由此推斷,干擾MEF2A基因?qū)е鲁杉〖毎芷谘娱L,細胞增殖受到抑制,與Wang等[43]的研究結(jié)果相似。然而,MEF2A影響成肌細胞的具體分子機制還需要更進一步的探究。
本研究成功構(gòu)建了關(guān)嶺牛MEF2A干擾載體,干擾MEF2A基因明顯抑制關(guān)嶺牛成肌細胞的增殖。同時推測,MEF2A通過調(diào)控成肌細胞中肌細胞增強因子基因和相關(guān)生長凋亡基因的表達來影響肌肉發(fā)育性能,為進一步完善MEF2A基因?qū)﹃P(guān)嶺牛生長發(fā)育的影響和探究地方優(yōu)質(zhì)種質(zhì)資源提供了理論參考。