食性分析是动物生态学中理解物种生存适应机制的重要手段,有助于理解捕食者、食物资源与环境条件之间的相互关系
[1-3]。传统食性研究通过摄食观察或对胃容物、粪便分析来直接了解食物资源信息,但受观察周期、食物分类经验与消化程度等影响,在确定食物资源的准确性和全面性上有所欠缺
[4-5]。近年来,DNA宏条形码技术为食性研究提供了新视角,该技术通过对食物残存DNA进行高通量测序与序列比对,能够实现对食物组成定性与相对定量分析,具有高灵敏度、高分辨率和广谱检测能力,可精准识别微量、高度降解或形态相近的食物成分,极大降低了人工视觉分类造成的误差
[6-8]。目前,该技术已广泛应用于啮齿类动物食性研究中
[9-10]。
大林姬鼠(
Apodemus peninsulae)是我国北方地区一种常见的杂食性偏植食性的小型啮齿类动物
[11]。多数研究发现大林姬鼠主要分布于林区生境,偏好郁闭度和落叶覆盖度高的乔木生境
[12],其植物性食物组成多与所处生境中的植物类型与资源量有关,多以植物果实和种子为食
[13-14]。然而,动物在不同生境中的摄食策略可能因其独特的食物资源格局而发生适应性调整。本课题组2020—2022年持续监测发现,熔岩台地生境有大林姬鼠分布,其种群数量逐年增加,并已发展为优势鼠种。熔岩台地由火山熔岩溢流冷凝形成,整体地势平坦,伴有凸凹起伏的地表构造,由于土壤层较薄,其植被类型以灌木和草本为主
[15]。熔岩台地生境类型独特,植物性食物资源有别于林区生境,大林姬鼠在此生境中的食物组成特征尚未可知。因此,本研究运用DNA宏条形码技术对镜泊湖熔岩台地生境中大林姬鼠春、秋两季的植物性食物组成进行了分析,旨在揭示该物种在熔岩台地生境中不同季节的植物性食物组成特征及差异,为啮齿动物的摄食生态学研究提供基础资料,并为理解动物对特殊生境的适应机制提供新的案例。
1 研究区概况
镜泊湖熔岩台地位于黑龙江省牡丹江市西南部的地质公园自然保护区内(44°03'46″—44°03
′49″N,128°58
′40″—129°00
′20″E),平均海拔320 m。该区属温带大陆性气候,年平均气温3.6 ℃,降水多在夏季,年降水量为506.4 mm
[16]。受风化程度及成土过程差异的影响,研究区植被群落自火山口中心台地向外围逐步分化,依次为地衣-苔藓群落、矮草本群落、高草本群落、矮灌木群落、高灌木群落和小乔木-灌木群落
[17]。
2 材料与方法
2.1 样本采集
于秋季(2021年9—11月)和春季(2022年4—6月)采用铗捕法捕捉大林姬鼠。选取包含灌木丛、草丛等生境的3块样地,每块样地设计3条长200 m的样线,样线间隔20 m,每条样线等间距设置40个捕鼠铗。于每月中上旬投放熟瓜子为诱饵,24 h后检查并记录捕获数量并补充诱饵,持续捕捉1周。累计捕获27只躯体完整(未受降水腐蚀、微生物降解或动物啃食)的大林姬鼠样本,包括春季13只(SC1 ~ SC13)、秋季14只(QC1~QC14)。样本带回实验室后,在无菌条件下解剖,取0.3~0.5 g胃容物置于2 mL离心管中,放入-80 ℃冰箱保存,待下一步分析。本研究已通过牡丹江师范学院实验动物福利伦理委员会审批(批准号:20210066)。
2.2 DNA提取
使用OMEGA M5635-02试剂盒提取胃容物样本总DNA。采用荧光分光光度计(Quantifluor ST fluorometer,Promega,E6090)测定荧光值以检测DNA浓度,并使用1.2%琼脂糖凝胶电泳检测DNA质量。
2.3 PCR扩增与测序
对于合格的样本检测区域进行高保真PCR扩增,设置3个重复,以标准基因组DNA Mix作为阳性对照。以胃容物DNA为模板,用引物
trnL(UAA)-F(5'-GGGCAATCCTGAGCCAA-3')和
trnL(UAA)-R(5'-CCATTGAGTCTCTGCACCTATC-3')
[18]扩增叶绿体
trnL(UAA)基因片段。PCR扩增总体积为25.00 μL,包含5 × 反应缓冲液5.00 μL,5 × GC缓冲液5.00 μL,dNTP(2.5 mmol/L)2.00 μL,正、反向引物(10 μmol/L)各1.00 μL,DNA模板2.00 μL,双蒸水8.75 μL,Q5DNA聚合酶0.25 μL。PCR反应程序为98 ℃预变性2 min;98 ℃变性15 s,55 ℃退火30 s,72 ℃延伸30 s,共20个循环;72 ℃终延伸5 min,4 ℃保存。PCR产物经1.5%琼脂糖凝胶电泳检测后,采用磁珠法纯化回收,所得上清液转移至1.5 mL离心管保存。对PCR产物进行定量后,根据所需测序量按相应比例混合,送至上海派诺森生物科技股份有限公司测序。
2.4 测序数据处理
测序所得原始数据使用cutadapt软件切除序列的引物片段,并弃去未匹配引物的序列;使用VSEARCH软件对序列进行拼接、质控、去重和去除嵌合体,获得高质量序列;使用QIIME2软件在97%相似度水平对高质量序列聚类,生成OTU表,并将该表分别与Greengenes数据库(
https://ngdc.cncb.ac.cn/databasecommons/database/id/3120)和SILVA数据库(
http://www.arb-silva.de)进行比对与物种注释。
2.5 数据分析
为反映大林姬鼠对不同植物性食物的摄食偏好,计算各食物类群的相对丰度(relative abundance,RA)和出现频率(frequency of occurrence,FO)。使用Primer 5软件计算摄取植物性食物的观测物种指数及Shannon、Simpson等α多样性指数,用以比较春秋两季植物性食物多样性的差异,采用Mann-Whitney U检验判断其差异显著性。为分析不同样本间植物性食物组成的β多样性,使用R软件ape包进行基于Bray-Curtis距离的主坐标分析(principal co-ordinates analysis,PCoA),结合Adonis非参数检验分析春秋两季植物性食物组成差异的显著性。
3 结果
3.1 大林姬鼠植物性食物组成
大林姬鼠春季植物性食物共有82属,隶属于30目46科。在科水平上,相对丰度大于1.00%的植物有11种,其中主要优势类群7种(
RA ≥ 5.00%),分别为榆科(Ulmaceae,
RA = 31.48%,)、百合科(Liliaceae,
RA = 13.70%)、菊科(Asteraceae,
RA = 10.16%)、豆科(Fabaceae,
RA = 8.34%)、薯蓣科(Dioscoreaceae,
RA = 7.26%)、牻牛儿苗科(Geraniaceae,
RA = 7.05%)和十字花科(Brassicaceae,
RA = 5.00%)(
图1A)。从出现频率来看,榆科与菊科出现频率最高,均为100.00%,其次是百合科、豆科、无患子科(Sapindaceae)和蔷薇科(Rosaceae),均为92.31%(
图1B)。
在属水平上,相对丰度大于1.00%的植物有12种,其中主要优势类群7种(
RA ≥ 5.00%),分别为榆属(
Ulmus,
RA = 31.48%)、洼瓣花属(
Lloydia,
RA = 13.67%)、小滨菊属(
Leucanthemella,
RA = 9.95%)、薯蓣属(
Dioscorea,
RA = 7.26%)、老鹳草属(
Geranium,
RA = 7.05%)、大豆属(
Glycine,
RA = 5.94%)和葶苈属(
Draba,
RA = 5.00%)(
图1C)。在出现频率上,榆属和小滨菊属出现频率最高,为100.00%,其次是洼瓣花属、大豆属和槭属(
Acer),为92.31%(
图1D)。
大林姬鼠秋季植物性食物共有116属,隶属于35目58科。在科水平上,相对丰度大于1.00%的植物有16种,主要优势类群7种(
RA ≥ 5.00%),分别为豆科(
RA = 15.56%)、桦木科(Betulaceae,
RA = 13.80%)、葡萄科(Vitaceae,
RA = 11.85%)、禾本科(Poaceae,
RA = 9.90%)、壳斗科(Fagaceae,
RA = 9.09%)、百合科(
RA = 7.01%)和菊科(
RA = 6.06%)(
图2A)。出现频率最高的为豆科与菊科,出现频率为100.00%,其次是葡萄科、禾本科、无患子科(Sapindaceae)和润肺木科(Strombosiaceae),出现频率为92.86%(
图2B)。
在属水平上,相对丰度大于1.00%的有18种,其中主要优势类群7种(
RA ≥ 5.00%),分别为榛属(
Corylus,
RA = 13.80%)、蛇葡萄属(
Ampelopsis,
RA = 11.85%)、大豆属(
RA = 9.26%)、栎属(
Quercus,
RA = 9.09%)、野黍属(
Eriochloa,
RA = 7.63%)、洼瓣花属(
RA = 7.00%)和小滨菊属(
RA = 5.61%)(
图2C)。出现频率最高的为大豆属,出现频率为100.00%,其次为蛇葡萄属、小滨菊属和槭属,出现频率为92.86%(
图2D)。
3.2 植物性食物组成的季节性差异
数据分析结果表明,大林姬鼠在春秋两季植物性食物共有130属,隶属于36目66科,其中两季共有植物有68属,春季特有14属,秋季特有48属,分别占食物总属数的52.31%、10.77%和36.92%(
图3)。共有类群包括大豆属、洼瓣花属和小滨菊属等,春季特有类群主要包括绢藓属(
Entodon)、松属(
Pinus)和苦参属(
Sophora)等,秋季特有类群主要包括堇菜属(
Viola)、蒿属(
Artemisia)和柳属(
Salix)等。
结合实地调查与测序结果,将植物性食物种类按植物生活型进一步划分。结果发现,春季,大林姬鼠的植物性食物以草本植物为主(
RA = 46.99%),木本植物次之(
RA = 42.97%),苔藓(
RA = 2.16%)和藻类(
RA = 0.25%)的取食占比较低;秋季,木本植物成为其主要食物来源(
RA = 50.38%),草本植物次之(
RA = 49.22%),苔藓与藻类的取食仅为微量(
RA < 0.01%)(
图4)。
3.3 植物性食物多样性的季节性差异
从大林姬鼠植物性食物多样性特征来看(
图5A),秋季的观测物种指数、Shannon指数、Simpson指数和Chao1指数均高于春季,表明秋季植物性食物多样性、丰富度高于春季,但差异不显著(
P > 0.05)。基于主坐标分析(PCoA)的空间分布特征显示(
图5B),植物性食物组成在不同季节间呈现明显的组间分异趋势,说明大林姬鼠植物性食物组成随季节更替发生变化。Adonis检验结果表明,季节变化对大林姬鼠的植物性食物组成具有极显著影响(
R² = 0.09,
P = 0.005)。
4 讨论
食物资源利用方式是动物在不同环境条件下形成的一种适应性策略。Soininen
et al.
[19]对加拿大极地地区的研究发现,在食物资源匮乏时,啮齿动物会取食苔藓以维持能量代谢。本研究利用DNA宏条形码技术分析发现,镜泊湖熔岩台地大林姬鼠的植物性食物以草本和木本植物为主,部分时段也取食苔藓和藻类,其中苔藓和藻类仅在春季植物性食物组成中有较高相对丰度和出现频率,这可能与当地食物资源的季节可利用性有关。我国东北地区春季气温波动大,可利用植物资源相对有限
[20],而苔藓和藻类有较强的低温适应能力
[21-22],因而可能成为春季大林姬鼠的部分替代性植物性食物,以减缓春季植物性食物资源短缺所带来的生存压力。
本研究结果表明,大林姬鼠的植物性食物组成在不同季节呈现明显的变化特征,春季时其植物性食物以草本植物为主。春季植被相对稀疏且种类较少,草本植物因生长周期短且易获取
[23],对大林姬鼠来说是资源量丰富的食物来源。尽管草本植物整体相对丰度较高,但从单一食物组分来看,木本植物中的榆属是大林姬鼠春季的主要植物性食物,相对丰度为31.48%。实地调查发现,研究区域分布有较多数量的榆树,每年4—5月是其果实成熟期,这为此时的大林姬鼠提供了丰富的食物来源。秋季大林姬鼠则主要取食木本植物,食物丰富度和多样性均高于春季。根据植物科学数据中心(
https://www.plantplus.cn/cn)资料,该区域大部分木本植物果期为9—10月,与本研究秋季采样时间相符。此时食物资源充足,大林姬鼠可能倾向于根据食物资源可利用性贮食,以应对冬季低温和食物资源匮乏的影响
[13,24]。
食物资源的质量和适口性是影响动物食性的关键因素
[9]。大豆属、小滨菊属和洼瓣花属等植物在春秋两季均具有较高的相对丰度与出现频率,这与3种植物在研究区域分布广泛有关。其中,大豆属植物的蛋白质和油脂含量高、种皮薄而适口性好,小滨菊属和洼瓣花属植物分布范围较广且易于采食,因此在大林姬鼠春秋两季植物性食物资源中占据一定比例。
已有研究表明,啮齿动物会结合所在生境类型展示出多样化的摄食策略
[25]。本研究在比较大林姬鼠春秋两季植物性食物多样性指数时发现,部分样本(春季SC1、SC9及秋季QC2、QC11)在多样性指数分布上明显偏离整体趋势。从食性组成来看,春季异常样本食物种类丰富,表现出较宽的食谱广度,而秋季异常样本则呈现出摄食对象高度集中、食物结构单一的特征。这一现象可能与样本采集的微生境食物资源分布格局有关。调查显示,春季异常样本采集于草丛生境,秋季异常样本采集于小乔木与灌木交错生境。我们推测,在不同微生境中,植物性食物资源的分布格局可能驱动了大林姬鼠摄食行为的差异,在资源稀疏或分散的草丛生境中,大林姬鼠倾向于拓宽取食生态位以满足营养需求,而在资源相对丰富且集中分布的灌丛生境中,则偏好取食少数优势类群食物以提高摄食效率。值得注意的是,个体摄食偏好、种间竞争压力等因素也会对啮齿动物食性组成产生影响
[26-27],这一推测仍有待进一步研究验证。
本研究区域地形地貌复杂,植被群落类型和微生境的多样性为大林姬鼠提供了丰富的植物性食物资源和庇护场所,使其能根据季节变化灵活调整摄食策略以利于种群规模逐步扩大。本研究仅限于大林姬鼠植物性食物成分分析,尚不能反映其完整的食性结构,未来将开展动物性食物检测,进一步解析熔岩台地生境中啮齿动物食物组成特征。
黑龙江省省属高等学校基本科研业务费科研项目(1453ZD026)
牡丹江师范学院教育教学改革项目(SJGY20220619)