摘 要 黑線姬鼠（Apodemus agrarius）是農(nóng)林業(yè)典型害鼠，為探究黑線姬鼠對(duì)食物資源利用的變化特征，選取張廣才嶺東部地區(qū)三道關(guān)針闊葉混交林為研究地點(diǎn)，分別于夏季（2020年6—8月）和秋季（2020年9—11月），采用鋏日法和樣線法采集黑線姬鼠樣本，利用DNA宏條形碼技術(shù)進(jìn)行黑線姬鼠食性分析。結(jié)果表明：黑線姬鼠的夏季食物主要為草本植物，胃容物包含14科19屬21種，其中以稗（Echinochloa crus-galli，87. 0%）和堪察加費(fèi)菜（Phedimus kamtschaticus，7. 0%）為主；秋季食物主要為喬木植物果實(shí)，胃容物包含15科40屬44種，其中以蒙古櫟（Quercus mongolica，48. 3%）、濕地風(fēng)毛菊（Saussurea umbrosa，17. 0%）、樟子松（Pinus sylvestris var. mongolica，7. 9%）和葎草（Humulus scandens，6. 2%）為主；夏秋兩季均攝食的有14種，以濕地風(fēng)毛菊為主。通過(guò)Shannon-Wiener指數(shù)（夏季，H′ = 0. 50；秋季，H′ = 0. 83）、Simpson指數(shù)（夏季，D = 0. 77；秋季，D = 0. 62）和Chao1 指數(shù)（夏季，Chao1 = 12. 29；秋季，Chao1 =13. 50）等多樣性指數(shù)分析，發(fā)現(xiàn)黑線姬鼠的秋季食物組成較夏季更豐富，食物物種多樣性更高。研究結(jié)果探明了張廣才嶺三道關(guān)林區(qū)黑線姬鼠夏秋兩季的食物資源利用特征，為科學(xué)開(kāi)展鼠害防治提供基礎(chǔ)資料。

關(guān)鍵詞：黑線姬鼠；食性；季節(jié)變化；DNA宏條形碼；張廣才嶺

中圖分類號(hào)：Q958. 1

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：2310 - 1490（2024）- 04 - 0727 - 07

DOI：10.12375/ysdwxb.20240405

食物是物種生存和群落動(dòng)態(tài)變化的物質(zhì)與能量基礎(chǔ)，是生存的資源［1］。對(duì)于嚙齒動(dòng)物來(lái)說(shuō)，食物資源是否充沛不僅影響個(gè)體的行為特征，也會(huì)影響種群數(shù)量的動(dòng)態(tài)變化。研究發(fā)現(xiàn)，在食物資源相對(duì)匱乏時(shí)，營(yíng)養(yǎng)不良的鼠類比營(yíng)養(yǎng)正常的鼠類更傾向于打斗，以獲得更多的食物資源［2］。通過(guò)附加食物的實(shí)驗(yàn)研究表明，相比無(wú)附加食物的種群，附加食物能有效降低鼠類種群內(nèi)部的競(jìng)爭(zhēng)行為，并使種群數(shù)量維持在較高水平［3］。上述研究表明食物資源的可利用性對(duì)鼠類種群數(shù)量變化具有潛在影響。

黑線姬鼠（Apodemus agrarius）是一種小型農(nóng)林業(yè)害鼠［4］，遍布于我國(guó)東部各地及新疆西北部［5］。黑線姬鼠屬雜食性偏植物食性，具有較寬泛的食物資源利用類型［6?7］。根據(jù)本課題組2014—2018年對(duì)張廣才嶺東部地區(qū)的持續(xù)監(jiān)測(cè)可知，初期黑線姬鼠主要分布于農(nóng)田生境［8］，而后期針闊葉混交林中的黑線姬鼠數(shù)量急劇增加，并成為三道關(guān)林區(qū)的優(yōu)勢(shì)物種［9］。在黑線姬鼠分布區(qū)域轉(zhuǎn)變以及種群數(shù)量增加這一過(guò)程中，黑線姬鼠的食物資源利用特征如何變化尚未可知。食性分析是掌握鼠類對(duì)食物資源利用和環(huán)境適應(yīng)能力情況的重要手段［10?11］。傳統(tǒng)食性分析方法一般包括糞便顯微組織分析法［12］、胃容物分析法、剖胃目視法［13?15］、野外觀察法［16］和幾何形態(tài)學(xué)法［17］，其優(yōu)點(diǎn)是成本低、易操作［18］，但具有工作量大、食物組分難以準(zhǔn)確區(qū)分的缺點(diǎn)［19］。隨著分子生物技術(shù)的發(fā)展，目前多數(shù)研究采用基于高通量測(cè)序的DNA宏條形碼技術(shù)進(jìn)行食性分析，這大大提高了食物分析的準(zhǔn)確性［20?23］。本研究應(yīng)用DNA宏條形碼技術(shù)分析張廣才嶺東部地區(qū)黑線姬鼠夏秋兩季的植物性食性特征，在野生動(dòng)物食性研究方面突破了傳統(tǒng)的顯微分析方法，為野生動(dòng)物營(yíng)養(yǎng)生態(tài)學(xué)的研究提供了新的方法和思路，以期為科學(xué)制定鼠害防治策略提供基礎(chǔ)資料。

1 材料與方法

1. 1 研究區(qū)概況

三道關(guān)林場(chǎng)（44°41′—44°51′ N，129°17′—129°35′ E）位于牡丹江市西北部，屬長(zhǎng)白山系張廣才嶺東坡。該區(qū)域?qū)贉貛Т箨懶约撅L(fēng)氣候，年平均降水量為537 mm，海拔380 ～ 550 m，屬于流水侵蝕山地［24］。研究區(qū)為針闊葉混交林生境，主要有紅松（Pinus koraiensis）、蒙古櫟（Quercus mongolica）、胡桃楸（Juglans mandshurica）、樟子松（Pinus sylvestrisvar. mongolica）、紫椴（Tilia amurensis）和毛榛（Cory?lus mandshurica）等［25?26］。

1. 2 樣品采集

分別于夏季（2020年6—8月）和秋季（2020年9—11月）采用鋏日法（15 cm × 8 cm鐵板捕鼠鋏）采集黑線姬鼠樣本。在三道關(guān)林場(chǎng)設(shè)置2塊生境環(huán)境相同的調(diào)查樣地，每塊樣地按3條樣線布鋏，樣線間隔20 m，鋏距5 m，每條樣線放置100個(gè)捕鼠鋏。以炒熟的白瓜子為誘餌，24 h后查看捕獲情況并補(bǔ)充誘餌。排除捕獲的其他鼠類物種后，共采集黑線姬鼠21只，其中夏季 7只（樣本名稱：XH1、XH2、XH3、XH4、XH5、XH6和XH7）、秋季14只（樣本名稱：QH1、QH2、QH3、QH4、QH5、QH6、QH7、QH8、QH9、QH10、QH11、QH12、QH13和QH14）。本研究已通過(guò)牡丹江師范學(xué)院實(shí)驗(yàn)動(dòng)物福利倫理委員會(huì)的審查（2023002）。

1. 3 樣品制備與DNA 宏條形碼分析

在無(wú)菌環(huán)境下解剖黑線姬鼠獲取胃容物，取0. 3 ～ 0. 5 g 胃容物置于2 mL 離心管中，放入冰箱-80 ℃保存，以待后續(xù)分析。在無(wú)菌條件下向胃容物中加入TE緩沖液進(jìn)行研磨，加入10 μL蛋白酶K，于56 ℃恒溫條件下消化12 h 后，使用FAST 土壤DNA提取試劑盒（MP Biomedicals， Santa Ana，CA）提取DNA。所得產(chǎn)物溶于TE緩沖液中，經(jīng)2%瓊脂糖凝膠電泳檢測(cè)（110 V，20 min）后于冰箱-20 ℃保存。

對(duì)合格樣本的檢測(cè)區(qū)域設(shè)置3個(gè)重復(fù)實(shí)驗(yàn)，進(jìn)行高保真PCR擴(kuò)增，同時(shí)以標(biāo)準(zhǔn)基因組DNA Mix作為陽(yáng)性對(duì)照，擴(kuò)增引物選用葉綠體基因trnL（UAA）內(nèi)含子序列中的P6環(huán)（10 143 bp）。以胃容物DNA為模板，16S UAA為目的區(qū)域，用引物trnL（UAA）-F（5′-GGGCAATCCTGAGCCAA-3′）和trnL（UAA）-R（5′-CCATTGAGTCTCTGCACCTATC-3′）［27］對(duì)trnL 葉綠體基因片段進(jìn)行擴(kuò)增。PCR 擴(kuò)增體系總體積為25. 00 μL，包括5×反應(yīng)緩沖液5. 00 μL，5 × GC緩沖液5. 00 μL，dNTP（2. 5 mmol/L）2. 00 μL，正向引物（10 μmol/L）1. 00 μL，反向引物（10 μmol/L）1. 00 μL，模板DNA（DNA Template）2. 00 μL，雙蒸餾水8. 75 μL，Q5 DNA聚合酶0. 25 μL。PCR擴(kuò)增程序：98 ℃預(yù)變性2 min；98 ℃變性15 s，55 ℃退火30 s，72 ℃延伸30 s，共20個(gè)循環(huán)；72 ℃終延伸5 min，4 ℃保存。將PCR產(chǎn)物用磁珠純化回收，將上清液移到1. 5 mL離心管中保存，采用Illumina MiSeq/NovaSeq高通量測(cè)序平臺(tái)制備文庫(kù)以及對(duì)群落DNA片段進(jìn)行雙端（paired-end）測(cè)序。測(cè)序由上海天昊生物科技有限公司完成。

1. 4 數(shù)據(jù)分析

原始數(shù)據(jù)經(jīng)過(guò)QIIME2軟件的cutadapt插件去除可能包含的adapter序列和引物后，統(tǒng)計(jì)原始數(shù)據(jù)質(zhì)量。為得到高質(zhì)量的測(cè)序數(shù)據(jù)，以提高后續(xù)生物信息分析的準(zhǔn)確性，使用QIIME2軟件的DADA2插件對(duì)數(shù)據(jù)進(jìn)行質(zhì)量過(guò)濾、降噪、拼接及去嵌合體。采用擴(kuò)增序列變體（amplicon sequence variant，ASV）僅去掉重復(fù)序列，不進(jìn)行相似度聚類，使用QIIME2 的DADA2插件進(jìn)行質(zhì)控后直接生成分類單元，并對(duì)每個(gè)樣品進(jìn)行ASV抽平，然后選取有效的OTU代表性序列，與數(shù)據(jù)庫(kù)（16S默認(rèn)Greengenes Database 13_8版本）進(jìn)行比對(duì)，采用樸素貝葉斯分類器（na?veBayes classifier）注釋被挑選出來(lái)的有效特征序列，將每個(gè)樣本的分類信息按照不同的分類水平整合到一起，導(dǎo)入QIIME2 軟件，再運(yùn)用QIIME2 軟件的DADA2插件進(jìn)行質(zhì)控后直接生成物種注釋信息，在物種注釋信息的基礎(chǔ)上，去掉注釋信息為葉綠體和不能注釋到界級(jí)別的OTU及其包含的序列，并刪除豐度小于0. 005% 的特征，基于OTU 的絕對(duì)豐度及注釋信息，對(duì)每個(gè)樣品在7個(gè)分類水平（界、門(mén)、綱、目、科、屬、種）上的所有序列數(shù)目所占總序列數(shù)的比例進(jìn)行統(tǒng)計(jì)，根據(jù)ASV/OTU表，統(tǒng)計(jì)研究樣本在界、門(mén)、綱、目、科、屬、種7 個(gè)分類水平上的物種豐度。用R軟件包中的ggplot2和VennDiagram分析物種相對(duì)豐度柱狀圖，進(jìn)而分析黑線姬鼠夏秋季植物性食物組成，由于在物種注釋時(shí)可能存在錯(cuò)配現(xiàn)象，導(dǎo)致出現(xiàn)的物種不是本地物種，在這種情況下，將非本地物種所在屬作為注釋結(jié)果。采用觀測(cè)物種數(shù)、Shannon-Wiener指數(shù)（H'）、Simpson指數(shù)（D）和Chao1指數(shù)進(jìn)行兩季食物多樣性比較，公式為

式中：Pi為第i 種植物的相對(duì)豐度；N 為所有序列數(shù)，ni指第i 個(gè)ASV所含的序列數(shù)。

2 結(jié)果

2. 1 黑線姬鼠夏季植物性食性特征分析

黑線姬鼠夏季胃容物中共鑒定出14科19屬21種植物，植物相對(duì)豐度（relative abundance，RD）大于1. 0%的植物有5種。黑線姬鼠夏季優(yōu)勢(shì)食物為禾本科（Poaceae）稗屬（Echinochloa）稗（E. crus-galli），達(dá)植物總豐度的87. 0%，其后依次為景天科（Crassu?laceae）費(fèi)菜屬（Phedimus）堪察加費(fèi)菜（Ph. kamts?chaticus）、菊科（Asteraceae）風(fēng)毛菊屬（Saussurea）濕地風(fēng)毛菊（S. umbrosa）、菊科蒿屬（Artemisia）蔞蒿（A. selengensis）和松科（Pinaceae）松屬（Pinus）樟子松，分別占植物總豐度的7. 0%、1. 4%、1. 4% 和1. 1%（圖1）。

2. 2 黑線姬鼠秋季植物性食性特征分析

黑線姬鼠秋季胃容物中共鑒定出15科40屬44種植物，其中RD gt; 1. 0%的植物有10種。黑線姬鼠秋季優(yōu)勢(shì)食物為殼斗科（Fagaceae）櫟屬（Quercus）蒙古櫟，RD = 48. 3%。此外，RD大于5. 0%的植物包括濕地風(fēng)毛菊（RD = 17. 0%）、樟子松（RD = 7. 9%）、大麻科（Cannabaceae）葎草屬（Humulus）葎草（H. scan?dens，RD = 6. 2%）（圖2）。

2. 3 黑線姬鼠夏秋兩季植物性食性的差異

經(jīng)Illumina高通量測(cè)序后，在21份樣本中共得到1 517 392條reads，經(jīng)降噪等處理后得到1 418 287條有效序列用于數(shù)據(jù)分析。經(jīng)聚類獲得77個(gè)ASV，其中黑線姬鼠夏季ASV有33個(gè)（占比42. 9%）、秋季ASV有58個(gè)（占比75. 3%，圖3）。夏季特有ASV共19個(gè)，秋季特有ASV共44個(gè)，夏秋兩季共有14個(gè)ASV，均少于兩季特有ASV數(shù)，說(shuō)明黑線姬鼠夏秋季植物性食物組成有差異。黑線姬鼠夏季主要以草本植物為食，取食種類主要有5種，而秋季主要以喬木植物果實(shí)為食，取食種類主要有10種。兩季相同的食物種類為濕地風(fēng)毛菊和樟子松，但夏季取食量較少（圖4）。

從黑線姬鼠的植物性食物多樣性特征來(lái)看，秋季的觀測(cè)物種數(shù)、Chao1指數(shù)和Shannon-Wiener多樣性指數(shù)均高于夏季，而夏季Simpson 指數(shù)高于秋季（表1），表明黑線姬鼠秋季食物組成多樣性明顯高于夏季，取食植物種類和數(shù)量更為豐富。

3 討論

動(dòng)物食性分析方法是科學(xué)掌握鼠類食物資源利用特征和環(huán)境適應(yīng)方式的重要手段。傳統(tǒng)食性分析方法簡(jiǎn)單、成本低，但在確定食物組成方面準(zhǔn)確性不高，很多植物碎屑難于鑒定［13］，而近年來(lái)應(yīng)用較多的分子分析方法能較好地解決這一問(wèn)題［28］，因此當(dāng)前嚙齒動(dòng)物的食性分析研究中多采用DNA宏條形碼技術(shù)［29?30］，但該方法也存在一定的局限性，如在實(shí)驗(yàn)操作時(shí)，樣本里混入了污染物或測(cè)序時(shí)發(fā)生序列錯(cuò)誤等，都會(huì)影響后續(xù)測(cè)序結(jié)果的準(zhǔn)確性［19］。Brandon-Mong et al. ［31］研究發(fā)現(xiàn)，DNA宏條形碼分析結(jié)果中出現(xiàn)高水平OTUs的原因一方面是文庫(kù)測(cè)序時(shí)引物雜交引起的錯(cuò)誤，另一方面是在檢測(cè)中發(fā)現(xiàn)存在殘留組織，如卵、微小標(biāo)本和內(nèi)共生細(xì)菌（如沃爾巴克氏體）等。

不同生境類型中的鼠類可利用的食物資源類型存在明顯差異。張旭新［32］研究發(fā)現(xiàn)，鼢鼠（Myospa?lacinae）在不同生境類型樣地中的食物種類組成明顯不同，說(shuō)明鼢鼠的取食在很大程度上受生境中植被類型的影響。本研究結(jié)果發(fā)現(xiàn)，在三道關(guān)林場(chǎng)針闊混交林生境中，黑線姬鼠夏季以草本植物為主要食物資源，秋季以喬木植物果實(shí)為主要食物資源。與以往研究相比，張廣才嶺農(nóng)田生境中黑線姬鼠在夏秋季主要利用人工種植的大豆（Glycine max）、馬鈴薯（Solanum tuberosum）等作為食物資源，在冬春季則以動(dòng)物性食物和綠色植物為食［23］。這一結(jié)果表明，一方面黑線姬鼠對(duì)果實(shí)等食物資源存在偏好，會(huì)優(yōu)先利用果實(shí)作為主要能量補(bǔ)充［33］，另一方面，不同生境中的食物資源類型會(huì)驅(qū)動(dòng)黑線姬鼠取食發(fā)生變化。

本研究表明，季節(jié)變化也會(huì)引起黑線姬鼠可利用的植物性食物資源類型發(fā)生轉(zhuǎn)變。三道關(guān)林場(chǎng)屬于溫帶大陸性季風(fēng)氣候，本研究所選樣地氣候濕潤(rùn)，低矮植被茂盛，喬木主要以蒙古櫟和樟子松為主。根據(jù)研究結(jié)果可知，黑線姬鼠夏秋兩季食物種類組成差異較大，依據(jù)植物科學(xué)數(shù)據(jù)中心（https：//www.plantplus. cn/cn）查詢得知，對(duì)于本研究所鑒定出的植物，夏季大部分植物的花期在6—7月，果期在8—10月，此時(shí)食物資源豐富，黑線姬鼠食物組成豐富，多為草本植物；秋季喬木的果期在9—10月，草本植物的果期在8—9月，食物種類較少，且多為果實(shí)，該查詢結(jié)果與夏秋兩季生長(zhǎng)的植物花果期相符，說(shuō)明了該數(shù)據(jù)的準(zhǔn)確性。有研究表明，黑線姬鼠主要食用植物的根、莖和葉，較少食用種子［13］，而本研究夏季食物中發(fā)現(xiàn)樟子松，可能是秋冬季節(jié)果實(shí)成熟自然脫落和貯食遺留。

黑線姬鼠為雜食性鼠類，除了以植物為主要食物外，還以蝗蟲(chóng)、蚯蚓等無(wú)脊椎動(dòng)物為食［13，23］。本研究未進(jìn)行動(dòng)物性食性檢測(cè)，尚不能獲取黑線姬鼠動(dòng)物性食物構(gòu)成的有效信息，后續(xù)需開(kāi)展相關(guān)研究。

參考文獻(xiàn)：

［1］ LI D W， ZHANG C Z， CAO Y W， et al. Food preference strategyof four sympatric rodents in a temperate forest in northeast China［J］. ZooKeys， 2023， 1158： 163-177.

［2］ WHATSON T S， SMART J L， DOBBING J. Social interactionsamong adult male rats after early undernutrition ［J］. British Jour?nal of Nutrition， 1974， 32（2）： 413-419.

［3］ 聶海燕， 劉季科， 蘇建平. 小型嚙齒動(dòng)物種群系統(tǒng)調(diào)節(jié)復(fù)合因子理論的野外實(shí)驗(yàn)研究： 食物可利用性和捕食對(duì)根田鼠種群空間行為的作用模式及其對(duì)種群調(diào)節(jié)的探討［J］. 獸類學(xué)報(bào)，1995， 15（1）： 41-52.

NIE H Y， LIU J K， SU J P. Field experimental studies on themultifactorial hypothesis of population system regulation for smallrodents： the effect pattern of food availability and predation onspacing behaviour of root voles and the function of spacing behav?iour in population regulation［J］. Acta Theriologica Sinica， 1995，15（1）： 41-52.

［4］ 李殿偉. 張廣才嶺森林嚙齒動(dòng)物分散貯食行為策略研究［D］.哈爾濱： 東北林業(yè)大學(xué)， 2019.

LI D W. Scatter-hoarding strategy of rodents to the forest treeseeds in the temperate forest of Zhangguangcai Mountains in north?eastern China［D］. Harbin： Northeast Forestry University， 2019.

［5］ 楊再學(xué)， 金星， 郭永旺， 等. 貴州省不同地區(qū)黑線姬鼠種群數(shù)量動(dòng)態(tài)分析［J］. 山地農(nóng)業(yè)生物學(xué)報(bào)， 2015， 34（1）： 13-17；27； 2.

YANG Z X， JIN X， GUO Y W， et al. Analysis of populationquantity dynamics of Apodemus agrarius in different regions of Gui?zhou Province［J］. Journal of Mountain Agriculture and Biology，2015， 34（1）： 13-17； 27； 2.

［6］ 牛翠娟， 婁安如， 孫儒泳， 等. 基礎(chǔ)生態(tài)學(xué)［M］. 3版. 北京：高等教育出版社， 2015.

NIU C J， LOU A R， SUN R Y， et al. Foundations in ecology［M］. 3rd ed. Beijing： Higher Education Press， 2015.

［7］ 楊再學(xué)， 鄭元利， 金星， 等. 黑線姬鼠胎仔數(shù)變動(dòng)規(guī)律及對(duì)種群數(shù)量的影響［J］. 中國(guó)農(nóng)學(xué)通報(bào)， 2013， 29（33）： 363-367.

YANG Z X， ZHENG Y L， JIN X， et al. Fluctuation law of littersize and its influence on population quantity of Apodemus agrarius［J］. Chinese Agricultural Science Bulletin， 2013， 29（33）：363-367.

［8］ 金志民. 黑龍江省張廣才嶺林區(qū)嚙齒動(dòng)物群落格局動(dòng)態(tài)研究［D］. 哈爾濱： 東北林業(yè)大學(xué)， 2019.

JIN Z M. A study on dynamics of rodent community patterns in for?ests of Zhangguangcai Mountain， Heilongjiang Province， China［D］. Harbin： Northeast Forestry University， 2019.

［9］ 李殿偉， 姚旭， 金志民， 等. 張廣才嶺東部山地嚙齒動(dòng)物群落多樣性調(diào)查［J］. 中國(guó)森林病蟲(chóng)， 2019， 38（5）： 28-33； 37.

LI D W， YAO X， JIN Z M， et al. Investigation of rodent commu?nity diversity in the eastern forest of Zhangguangcai Mountains［J］. Forest Pest and Disease， 2019， 38（5）： 28-33； 37.

［10］ 李天保. 濟(jì)源市王屋山地區(qū)林棲鼠類食性分析［J］. 現(xiàn)代農(nóng)業(yè)科技， 2020（20）： 83-85； 88.

LI T B. Dietary analysis of forest rodents in Wangwu Mountainarea of Jiyuan City［J］. Modern Agricultural Science and Tech?nology， 2020（20）： 83-85； 88.

［11］ 雍仲禹， 郭聰， 張美文， 等. 嚙齒動(dòng)物食性研究的意義及方法評(píng)述［J］. 生態(tài)學(xué)雜志， 2011， 30（11）： 2637-2645.

YONG Z Y， GUO C， ZHANG M W， et al. Significance andmethodology of rodent’s food habit research： a review［J］. Chi?nese Journal of Ecology， 2011， 30（11）： 2637-2645.

［12］ ELLIOTT T F， VERNES K. Notes on the diets of four rodent spe?cies from Goodenough Island［J］. Australian Mammalogy，2020， 43（2）： 256-259.

［13］ 雍仲禹， 張美文， 郭聰， 等. 洞庭湖區(qū)黑線姬鼠食性調(diào)查［J］. 動(dòng)物學(xué)雜志， 2012， 47（3）： 115-121.

YONG Z Y， ZHANG M W， GUO C， et al. Dietary habit ofApodemus agrarius in Dongting Lake area［J］. Chinese Journal ofZoology， 2012， 47（3）： 115-121.

［14］ MULUNGU L S， MASSAWE A W， KENNIS J， et al. Differ?ences in diet between two rodent species， Mastomys natalensisand Gerbilliscus vicinus， in fallow land habitats in central Tanza?nia［J］. African Zoology， 2011， 46（2）： 387-392.

［15］ 金志民， 劉鑄， 李殿偉， 等. 黑龍江省東南部林區(qū)6種嚙齒動(dòng)物食性分析［J］. 畜牧與飼料科學(xué)， 2011， 32（5）： 11-13.

JIN Z M， LIU Z， LI D W， et al. Analysis on the feeding habitsof six kinds of rodents in the southeast forest region of Heilongji?ang Province［J］. Animal Husbandry and Feed Science， 2011，32（5）： 11-13.

［16］ YIP S J S， DICKMAN C R. Foraging and food selection in a des?ert rodent： diet shifts of the sandy inland mouse between popula?tion booms and busts［J］. Animals， 2023， 13（10）： 1702.

［17］ CANO A R G， FERNáNDEZ M H， áLVAREZ-SIERRA M á.Dietary ecology of Murinae （Muridae， Rodentia）： a geometricmorphometric approach［J］. PLoS One， 2013， 8（11）： e79080.

［18］ 湛振杰， 張超， 陳敏豪， 等. 基于DNA宏條形碼技術(shù)的大興安嶺北部歐亞水獺冬季食性分析［J］. 生物多樣性， 2023， 31（6）： 85-96.

ZHAN Z J， ZHANG C， CHEN M H， et al. DNA metabarcodingbasedwinter diet analysis of Eurasian otter （Lutra lutra） in thenorthern Greater Khingan Mountains［J］. Biodiversity Science，2023， 31（6）： 85-96.

［19］ 楊舒晗， 蔡延芳， 王野影， 等. 貴州茂蘭喀斯特地區(qū)隱紋花松鼠秋季食性分析［J］. 野生動(dòng)物學(xué)報(bào)， 2022， 43（2）：337-344.

YANG S H， CAI Y F， WANG Y Y， et al. Autumn diet analysison Swinhoe’s striped squirrel（Tamiops swinhoei） in MaolanKarst area， Guizhou Province， China［J］. Chinese Journal ofWildlife， 2022， 43（2）： 337-344.

［20］ 嚴(yán)麗君， 王普， 施啟龍， 等. 動(dòng)物食性分析在生態(tài)學(xué)中的應(yīng)用研究進(jìn)展：基于DNA宏條形碼技術(shù)［J］. 生態(tài)學(xué)報(bào)， 2023，43（8）： 3007-3019.

YAN L J， WANG P， SHI Q L， et al. Applications of animal dietanalysis based on DNA metabarcoding in ecological research［J］.Acta Ecologica Sinica， 2023， 43（8）： 3007-3019.

［21］ ZHANG X X， ZOU Y， ZOU X， et al. DNA metabarcoding un?covers the diet of subterranean rodents in China［J］. PLoS One，2022， 17（4）： e0258078.

［22］ LOPES C M， DE BARBA M， BOYER F， et al. DNA metabar?coding diet analysis for species with parapatric vs sympatric dis?tribution： a case study on subterranean rodents［J］. Heredity，2015， 114（5）： 525-536.

［23］ SOININEN E M， RAVOLAINEN V T， BR?THEN K A， et al.Arctic small rodents have diverse diets and flexible food selection［J］. PLoS One， 2013， 8（6）： e68128.

［24］ 李殿偉， 孟繁星， 劉佳慧， 等. 張廣才嶺東部?jī)煞N姬鼠兩性異形比較研究［J］. 湖北農(nóng)業(yè)科學(xué)， 2022， 61（20）： 105-110； 187.

LI D W， MENG F X， LIU J H， et al. Comparative study of twosexual dimorphisms of two species of Apodemus in the east ofZhangguangcailing［J］. Hubei Agricultural Sciences， 2022， 61（20）： 105-110； 187.

［25］ 王悅山， 王麗麗， 范春楠， 等. 張廣才嶺森林植物區(qū)系研究［J］. 北華大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2012， 13（5）： 573-577.

WANG Y S， WANG L L， FAN C N， et al. Study of flora inZhangguangcai Mountains［J］. Journal of Beihua University（Natural Science）， 2012， 13（5）： 573-577.

［26］ 孫越， 夏富才， 趙秀海， 等. 張廣才嶺溫帶次生針闊混交林物種組成和群落結(jié)構(gòu)特征［J］. 生態(tài)學(xué)報(bào)， 2017， 37（10）：3425-3436.

SUN Y， XIA F C， ZHAO X H， et al. The composition and com?munity structure of coniferous and broadleaved mixed forest spe?cies in the Zhangguangcai Mountains， China［J］. Acta EcologicaSinica， 2017， 37（10）： 3425-3436.

［27］ TABERLET P， COISSAC E， POMPANON F， et al. Power andlimitations of the chloroplast trnL （UAA） intron for plant DNAbarcoding［J］. Nucleic Acids Research， 2007， 35（3）： e14.

［28］ KHANAM S， HOWITT R， MUSHTAQ M， et al. Diet analysisof small mammal pests： a comparison of molecular and microhis?tological methods［J］. Integrative Zoology， 2016， 11（2）：98-110.

［29］ KLURE D M， GREENHALGH R， DEARING M D. Addressingnontarget amplification in DNA metabarcoding studies ofarthropod-feeding rodents［J］. Mammal Research， 2022， 67（4）： 499-509.

［30］ NEBY M， KAMENOVA S， DEVINEAU O， et al. Issues ofunder-representation in quantitative DNA metabarcoding weakenthe inference about diet of the tundra vole Microtus oeconomus［J］. PeerJ， 2021， 9： e11936.

［31］ BRANDON-MONG G J， GAN H M， SING K W， et al. DNA me?tabarcoding of insects and allies： an evaluation of primers andpipelines［J］. Bulletin of Entomological Research， 2015， 105（6）： 717-727.

［32］ 張旭新. 寧夏六盤(pán)山地區(qū)鼢鼠食性與胃內(nèi)容物微生物相關(guān)性研究［D］. 楊凌： 西北農(nóng)林科技大學(xué)， 2022.

ZHANG X X. Study on the correlation between diet and stomachcontents microbes of zokor in Liu-pan Mountains， Ningxia［D］.Yangling： Northwest Aamp;F University， 2022.

［33］ 李楠， 黎道洪， 袁萌. 亞熱帶森林嚙齒類動(dòng)物對(duì)林區(qū)兩種植物果實(shí)的捕食和擴(kuò)散研究： 以寬闊水自然保護(hù)區(qū)為例［J］. 四川動(dòng)物， 2011， 30（3）： 362-365； 371.

LI N， LI D H， YUAN M. Diffusion effects of plant seeds by ro?dents collecting fruits of two species plants in subtropical forest：a case of Kuankuoshui Nature Reserve， Guizhou of China［J］. Si?chuan Journal of Zoology， 2011， 30（3）： 362-365； 371.

基金項(xiàng)目：研究生科技創(chuàng)新項(xiàng)目（kjcx2023-118mdjnu）；黑龍江省自然科學(xué)基金項(xiàng)目（SS2021C006）；牡丹江師范學(xué)院橫向課題（2021H001）