鹞落坪保护区典型农林交错带传粉昆虫多样性及其与环境因子间的关系.pdf

doi:10.3969/j.issn.2095-1736.2023.04.070鹞落坪保护区典型农林交错带传粉昆虫多样性及其与环境因子间的关系姚晨晨1,2,侯银续3,郭东旭1,2,万 霞1,2(1.安徽大学 资源与环境工程学院,合肥 230601;2.湿地生态保护与修复安徽省重点实验室,合肥 230601;3.安徽省疾病预防控制中心,合肥 230601)摘 要2019 年 10 月至 2021 年 10 月,在鹞落坪保护区典型农林交错带采集传粉昆虫 2 551 只,通过形态学及DNA Barcoding 互证,鉴定出传粉昆虫 92 种,隶属于 4 目 34 科 73 属。其中,中华蜜蜂(Apis cerana)、三条熊蜂(Bom-bus trifasciatus)、灰带管食蚜蝇(Eristalis cerealis)和东方菜粉蝶(Pieris canidia)是优势种。物种多样性分析结果表明,以农田为主的样线比以森林为主的样线和的多样性更高,多样化的农作物能为传粉昆虫提供更好的栖息地和食物来源;样线则因临近旅游区域而受到较多人为干扰,多样性显著较低。此外,从 WorldClim 数据库中下载 19 个环境因子(Bio1Bio19),分析物种及系统发育多样性与研究区域内的环境因子间的相关性关系:其中,有14 个环境因子(Bio1,Bio413,Bio16,Bio1819)与系统发育多样性指数以及物种丰富度指数之间存在线性相关关系,冗余分析结果则支持温度因子对系统发育多样性及丰富度影响更显著,温度可能是影响鹞落坪保护区传粉昆虫多样性的重要生态因子。关键词 传粉昆虫;物种丰富度;DNA barcoding;系统发育多样性;环境因子中图分类号 Q968文献标识码 A文章编号 2095-1736(2023)04-0070-06Diversity of pollinators and their relationships with environmental factors in typical forest-agriculture ecotone of Yaoluoping Nature ReserveYAO Chenchen1 2 HOU Yinxu3 GUO Dongxu1 2 WAN Xia1 2 1.School of Resources and Environmental Engineering Anhui University Hefei 230601 China 2.Anhui Province Key Laboratory of Wetland Ecosystem Protection and Restoration Hefei 230601 China 3.Anhui Center for Disease Control and Prevention Hefei 230601 China Abstract From October 2019 to October 2021 2 551 insect pollinators were collected in the typical forest-agriculture ecotone of Ya-oluoping Nature Reserve.92 species of pollinators were identified by morphology and DNA Barcoding belonging to 73 genera 34 fam-ilies and 4 orders among which Apis cerana Bombus trifasciatus Eristalis cerealis and Pieris canidia were the dominant species.The results of species diversity analysis showed that transect I dominated by farmland had higher species diversity than transects II and dominated by forests because a variety of crops can provide better habitat and food source for pollinators while transect which was near a tourist area severely disturbed by human had significantly lower diversity.In addition 19 Bio1-Bio19 environment factors were downloaded from WorldClim database to analyze the correlation between diversity including species and phylogenetic diversity and environmental factors.There was a linear correlation between 14 environmental factors Bio1 Bio4-13 Bio16 Bio18-19 and di-versity index including phylogenetic diversity and species richness.The results of redundancy analysis supported that temperature factors had a more significant effect on phylogenetic diversity and richness.Temperature could be the key ecological factors influencing the diversity of pollinators in the reserve.07第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023收稿日期:2022-03-04;最后修回日期:2022-09-15基金项目:科技部科技基础资源调查专项(2018FY100405)作者简介:姚晨晨,硕士,研究方向为环境生态与资源管理,E-mail:3315447844 通信作者:万霞,博士,教授,研究方向为生物多样性保护与湿地生态,E-mail:wanxia Keywords pollinator species richness DNA barcoding phylogenetic diversity environmental factor 地球上约 87.5%的被子植物依赖动物传粉,昆虫在所有传粉动物中占据主导地位1。传粉昆虫与开花植物在长期协同进化过程中形成了自然界中主要的能量传递,对维持全球生物多样性、生态系统的平衡和相对稳定发挥重要作用2-3。然而,随着人为干扰加剧、气候变暖、栖息地大量被破坏以及杀虫剂的不合理使用,全球范围内的传粉昆虫多样性正急剧减少和丧失,对生态系统传粉服务功能和人类福祉产生显著的负面影响4-5。我国是传粉昆虫资源最丰富的国家之一,但由于前期研究基础薄弱,在传粉昆虫的现状、作用及保护策略等方面的研究远落后于欧美国家。已有研究也多关注常见的、具有特定生态与经济价值的传粉蜂类及一些有指示作用的蝶类6-7,缺乏大范围的野生传粉昆虫多样性调查研究,存在现状不清、对未来发展预测缺乏等问题8。因此,在当今全球传粉昆虫多样性下降的大趋势下,许多未经调查区域的传粉昆虫也面临着即将丧失的危险,开展传粉昆虫多样性的研究、保护与利用刻不容缓。传统生物多样性研究中常用的测度指标包括物种丰富度指数、Shannon 指数、Simpson 指数、均匀度指数等,这些指标主要基于经典形态学的物种识别,包含了群落中物种个体数目和相对多度的信息。然而,随着各项研究方法和技术的发展,有研究认为这些指标已不再完全适用9。准确的物种鉴定是生物多样性估计、物种保护的必要前提。传粉昆虫种类繁多,同一物种内往往存在外部形态差异显著的雌雄二型、不同季节的色斑型等现象,给基于形态学分类的物种鉴定带来了前 所 未 有 的 挑 战10-11。近 年 来,DNA 条 形 码(DNA Barcording)技术及系统发育关系构建在群落生态学、生物多样性调查和保护中日益受到关注,已广泛应用于物种的识别或新种和隐种的发现12-14。基于DNA 条形码构建的系统发育树是研究系统发育多样性的重要基础,能够反映物种间的系统发育关系与进化信息,因此,Faith 提出将系统发育多样性(Phyloge-netic diversity,PD)作为度量指标,用来衡量群落中不同物种的系统进化程度,PD 指数的计算主要依赖所构建的系统发育树各分支的长度,是对基于经典形态学分类的多样性指数的有效补充,为更好评估区域生物多样性提供新的视角和方向15。鹞落坪国家级自然保护区位于安徽省岳西县境内(11602 116 11 E,30 57 31 06 N),总 面 积123 km2,地理上具有“南北过渡,襟带东西”的显著特征,是次生植被保存较完好的北亚热带森林生态系统,其复杂的生态环境、丰富的食物资源为各种昆虫的生存、繁衍等提供了良好的条件16。但目前有关保护区昆虫多样性(尤其是传粉昆虫多样性)的研究仅有一些零星的报道17-18,未将形态学分类方法和 DNA Bar-coding 技术联合用于传粉昆虫的分类鉴定中,并且保护区内传粉昆虫的物种多样性、系统发育多样性及其与环境因子间的关系也缺少分析。因此,鉴于保护区农林交错带蕴含的丰富生物多样及其特有的地理区位,非常有必要开展交错带的传粉昆虫多样性研究,以期为大别山资源昆虫的保护与利用、生态环境监测以及北亚热带生物多样性保护等提供科学依据。1 材料与方法1.1 调查方法选取保护区内两个具有典型农林交错带特征的区域:中畈村(海拔高度 500800 m)和鹞落坪村(海拔高度 800 1 200 m),共设置 4 条样线,每条样线约2 km,从 2019 年 10 月2021 年 10 月采用网捕法沿样线对生境区域内的传粉昆虫进行了调查,采样时间选择在植物花期内及昆虫活动较频繁的 9:0017:00,以15 min 为一个时间段,观察、记录在规定时间内遇到的每种开花植物植株上的传粉昆虫种类及数量,拍照并摄像后,将其捕捉存于 75%的酒精离心管中;利用 GPS 记录样线的海拔,并在离心管上标明日期、地点和采集者,带回实验室保存、整理及鉴定。各样线详细信息见表 1。1.2 DNA 提取、PCR 扩增及测序取昆虫样本胸部或腿部肌肉组织,使用血液/细胞/组织试剂盒提取样品的总 DNA。COI 基因片段的扩增引物选择通用引物 LCO-1490 和 HCO-2198 19。PCR 反应总体系为 25 L,包含 10buffer 2.5 L,dNTPs 2 L,上 下 游 引 物 各 1.0 L,Taq 聚 合 酶0.15 L,DNA 模板 3.0 L,ddH2O 补齐体积。PCR 反应条件:94 预变性 2 min;94 变性 20 s,47 退火50 s,72 延伸 1 min,35 个循环;最后 72 再延伸9 min。扩增产物由通用生物(安徽)股份有限公司进行 COI 基因片段纯化及双向测序。1.3 物种鉴定依据形态学分类特征以及将所扩增的 COI 基因序列分别在 BOLD 系统(http:www.barcodinglife.org)中的“IDENTIFICATION”功能进行比对、或在 GenBank 数据库(http:www.ncbi.nlm.nih.gov)中进行 Blast 相似性搜索对物种进行鉴定。17第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023表 1 样线生境特征及植物群落组成Table 1 Habitat characteristic and constitution of plant communities in transects样线Transect海拔Altitude/m生境特征Habitat characteristic主要植被组成Main constitution of plant738791以农田为主乔木层:马尾松(Pinus massoniana);灌木层:茅栗(Castanea seguinii)、短柄袍栎(Quer-cus glandulifera);草本层:以十字花科 Brassicaceae、豆科 Fabaceae、菊科 Asteraceae 为主;农作物以小麦(Triticum aestivum)、玉米(Zea mays)、大豆(Glycine max)为主9971 068以森林为主乔木层:江南恺木(Alnus trabeculosa),香果树林(Emmenopterys henryi);灌木层:茅栗(Castanea seguinii)、胡颓子(Elaeagnus pungens);草本层:以罂粟科 Papaveraceae、菊科Asteraceae 为主;农作物以小麦(Triticum aestivum)、玉米(Zea mays)为主1 0611 149以森林为主,临近旅游区乔木层:黄山松(Pinus taiwanensis),厚朴(Houpoea officinalis);灌木层:多枝杜鹃(Rho-dodendron polycladum);草本层:以菊科 Asteraceae、十字花科 Brassicaceae 为主;农作物种植以玉米(Zea mays)为主1 0621 184以森林为主乔木层:黄山松(Pinus taiwanensis),栓皮栎(Quercus variabilis);灌木层:茅栗(Casta-nea seguinii)、短柄袍栎(Quercus glandulifera);草本层:以罂粟科 Papaveraceae、菊科Asteraceae、十字花科 Brassicaceae 为主1.4 数据分析1.4.1 条形码序列分析及系统发育树构建使用 Chromas 软件查看所测序列的峰图质量并进行人工校对,确保 DNA 序列的可靠性。经校对的序列用 MEGA 6.0 进行比对分析,共得到无碱基缺失并能翻译通过的序列 288 条,序列长度为 654 bp。利用DnaSP V6.1 软件计算数据集中的碱基组成频率(A/T/C/G)、变异位点、简约信息位点、核苷酸多态性及单倍型组成等信息20。从每个物种中随机挑选 1 条序列,共 92 条序列,利用 IQ-TREE(http:iqtree.cibiv.univie.ac.at)基于最大似然法(Maximum likelihood,ML)构建系统发育树21,选择半翅目同蝽科黑刺同蝽(Acanthosoma nigrospina)作为外群。1.4.2 物种多样性指数和系统发育多样性指数物种多样性分析选择 Marglef 丰富度指数(AMg)、Shannon 指数(H)、Simpson 指数(R)及 Invsimpson 指数(R)4 个指标,计算公式如下:Marglef 丰富度指数(AMg)AMg=S-1ln N(1)其中,S 代表物种数,N 代表总体个体数;Shannon 指数(H)H=-Piln Pi(2)其中,Pi=Ni/N,Pi为样本中昆虫 i 的个体数比例,Ni为群落中昆虫 i 的个体数,N 为昆虫个体总数;Simpson 指数(R)R=1-P2i(3)Invsimpson 指数(R)R=1P2i(4)其中,Pi为样本中昆虫 i 的个体数比例。利用 R 语言中的 vegan 包、picante 包分别计算不同样线中传粉昆虫的物种及系统发育多样性指数。1.4.3 物种多样性、系统发育多样性、环境因子之间的关系 从世界气候数据库 WorldClim v1.4(http:www.worldclim.org)中选取 19 个环境因子(Bio1Bio19),包括年平均气温(Bio1)、年平均日较差(Bio2)、等温性(Bio3)、温 度 季 节 变 化(Bio4)、最 热 月 份 最 高 温(Bio5)、最冷月份最低温(Bio6)、年较差(Bio7)、最湿润季度平均气温(Bio8)、最干旱季度平均气温(Bio9)、最温暖季度的平均温度(Bio10)、最冷季度的平均温度(Bio11)、年降水量(Bio12)、最湿润月份的降水量(Bio13)、最干旱季度的降水量(Bio14)、季节性降水量(Bio15)、最湿润季度的降水量(Bio16)、最干旱季度的降水量(Bio17)、最温暖季度的降水量(Bio18)、最冷季度的降水量(Bio19)。利用函数 cor.test(),lm()对各样线传粉昆虫物种多样性、系统发育多样性与该区域环境因子间进行 Pearson 相关性分析并构建回归模型。同时,为避免各环境因子间存在共线性,本研究基于相关性和冗余分析(RDA)对 19 个环境因子进行筛选。2 结果与分析2.1 DNA 条形码序列遗传特征研究共得到 COI 序列 228 条,长度为 654 bp。利用 DnaSP V6.1 对数据进行分析,序列中 G+C 含量为0.30,多态性位点(Number of polymorphic/segregating sites)为 428,总突变数目(Total number of mutations)为819,核苷酸多态性(Nucleotide diversity)为 0.218 89,单倍型数目(Number of Haplotypes)为 149,单倍型多态性(Haplotype diversity)为 0.994 0。2.2 传粉昆虫群落组成及优势种依据 BOLD 数据库和 GenBank 数据库对序列进行比对,并结合形态学分类特征对物种进行鉴定,结果显27第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023示 4 个农林交错带共采集到传粉昆虫 2 551 只,隶属于4 目(膜翅目、双翅目、鞘翅目、鳞翅目)的 34 科 73 属92 种(表 2)。其中,中华蜜蜂 Apis cerana(362 只)、三条熊蜂 Bombus trifasciatus(243 只)、灰带管食蚜蝇 Eri-stalis cerealis(733 只)和东方菜粉蝶 Pieris canidia(165只)均超过总个体数的 6%,为区内的优势种(表 3),所有优势种在各样线中均有分布。表 2 不同样线传粉昆虫群落组成Table 2 Community composition of pollinators in different transects 样线Transect科数Numberof family属数Number ofgenus种数Number ofspecies个体数Number of individuals315469710253953832202833282234454727合计3473922 551表 3 鹞落坪自然保护区传粉昆虫优势种(个体数30 只)Table 3 The dominant species(individual 30)ofpollinators in Yaoluoping Nature Reserve科Family优势种Dominant species个体数Number ofindividuals蜜蜂科Apidae三条熊蜂 Bombus trifasciatus243中国红光熊蜂 Bombus ignitus64中华蜜蜂 Apis cerana362胡蜂科Vespidae墨胸胡蜂 Vespa velutina30黑腹虎头蜂 Vespa basalis33隧蜂科Halictidae日本淡脉隧蜂 Lasioglossum occidens51食蚜蝇科Syrphidae印度细腹食蚜蝇 Sphaerophoria indiana58细扁食蚜蝇 Episyrphus balteatus59灰带管食蚜蝇 Eristalis cerealis733大灰优食蚜蝇 Eupeodes corollae55丽蝇科Calliphoridae不显口鼻蝇 Stomorhina obsoleta80伪尼等彩蝇 Isomyia pseudonepalana104大蚊科Tipulidae花蚊 Tipula oleracea34粉蝶科Pieridae东方菜粉蝶 Pieris canidia165宽边黄粉蝶 Eurema hecabe452.3 物种多样性、系统发育多样性及相关性分析物种多样性指数计算结果表明(表 4),样线(以农田为主)具有比样线和(以森林为主)明显更高的物种多样性,并且样线的物种多样性较低;Simp-son 指数与 Invsimpson 指数之间呈现显著正相关(r0.8,P0.05)。采用基于 ML 法构建的系统发育树拓扑结构对 PD 指数进行了计算(表 4),并对物种多样性指数和 PD 指数进行相关性检验,表明 PD 指数与物种多样性之间不存在显著相关性(P0.05),与物种丰富度指数存在显著正相关(P0.05)。2.4 多样性指数与环境因子间的相关性关系为分析鹞落坪自然保护区传粉昆虫物种多样性指数、系统发育多样性指数(PD)及其与环境因子之间的关系,利用 R 语言中 Vegan 包对 19 个环境变量与 4 条样线传粉昆虫 PD 指数之间进行 Pearson 相关性分析,结果表明,所有环境因子中有 14 个环境因子(Bio1,Bio413,Bio16,Bio18 19)与 PD 指数存在线性相关关系(P0.05)。通过冗余分析(Redundan-cy analysis,RDA)对环境因子进行筛选,结果表明环境因子 Bio1(年平均气温)和 Bio4(季节性温度变化方差)对 PD 指数具有显著影响(图 1)。对 PD 指数与环境因子之间进行相关性分析和线性回归,结果表明 2 个环境因子与系统发育多样性指数之间具有显著正相关性,分别是年平均气温 Bio1(R=0.992 3,P=0.007 7)和 季 节 性 温 度 变 化 方 差 Bio4(R=0.991 2,P=0.008 8)(图 2)。物种多样性指数与环境因子之间的 Pearson 相关性分析则表明,传粉昆虫物种多样性指数(Shannon 指数和 Simpson 指数)与环境因子之间不存在相关性,但物种丰富度指数(物种丰富度和 Marglef 丰富度指数)与温度因子和湿度因子均存在相关性,冗余分析支持温度因子对物种丰富度影响更显著。表 4 物种多样性指数和系统发育多样性指数Table 4 Species diversity index and phylogenetic diversity index 样线Transect丰富度Richness系统发育多样性指数 PDMarglef 指数Marglef indexShannon 指数Shannon indexSimpson 指数Simpson indexInvsimpson 指数Invsimpson index6912.7810.363.040.909.725310.537.732.570.846.09337.035.672.540.867.17549.848.042.810.899.163 讨论与结论本文较系统地研究了鹞落坪国家级自然保护区农林交错带传粉昆虫多样性,并分析物种多样性、系统发育多样性及其与环境因子间的关系。4 个农林交错带共采集到传粉昆虫 2 551 只,经鉴定隶属于 4 目 34 科73 属 92 种,其中,4 个物种即中华蜜蜂、三条熊蜂、灰37第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023带管食蚜蝇和东方菜粉蝶是该区开花植物上的优势种。本次调查的传粉昆虫种类数量足以说明昆虫在显花植物授粉、维护农林生态系统平衡和稳定等方面的重要性。图 1 系统发育多样性与环境因子间 RDA 分析Figure 1 RDA analysis of phylogenetic diversity andenvironmental factors图 2 系统发育多样性与不同环境因子之间相关性分析Figure 2 Correlation analysis between phylogenetic diversity anddifferent environmental factors传统生物多样性测量常用物种丰富度指数等作为测度指标,这些指标将物种作为等价的研究单元,容易忽略不同物种对生物多样性的贡献大小10。系统发育多样性则从物种的亲缘关系层面来衡量一个区域内物种间的分歧程度,考虑不同物种间的系统发育关系,能更准确和全面地预测总体生物多样性格局16。本研究分别对物种多样性及系统发育多样性指数进行了计算,结果表明样线比样线、样线的物种多样性更高,这主要是样线为保护区内的主要农业区,种植有小麦、玉米、大豆等多种农作物,作物的多样化为传粉昆虫提供了良好的食物源和栖息地,因而该区域内传粉昆虫多样性明显较高;同时,多种传粉昆虫也为农作物提供了更好的授粉服务,保障了作物的品质和产量,对维持农田生态系统的平衡具有重要作用4。而样线临近保护区内的旅游区,受到较为剧烈的人为干扰的影响,导致传粉昆虫多样性明显较低。不同样线的传粉昆虫系统发育多样性与物种丰富度指数之间存在显著相关性,而与物种多样性指数之间不存在相关性。这一结论也在动植物等多个类群的研究中都得到了证实22-24,表明对分类地位不同的物种或区域开展评估和制定保护策略时,应综合考虑物种多样性和系统发育多样性。因气候改变、人类干扰加剧、土地利用集约化以及杀虫剂的大量使用导致传粉昆虫正面临灭绝的威胁5,深入分析环境因素对传粉昆虫多样性的影响日益受到重视。本文的分析结果表明所选择的 19 个环境因子中的 14 个与系统发育多样性指数存在线性相关关系;冗余分析(RDA)则进一步表明年平均气温、温度季节变化方差与系统发育多样性间存在显著正相关性(P0.05)。物种丰富度指数与温度和湿度因子也均存在相关性,而温度因子对丰富度影响更显著,支持温度是影响传粉昆虫多样性的重要生态因子25。温度升高和气候变暖则通过影响传粉昆虫的生存、发育、栖息环境以及开花植物的物候期,进而导致其丰富度、分布格局和种间关系等发生变化26-27。在全球气候变暖的大背景下,鹞落坪国家级自然保护区也显现了北亚热带森林生态系统的脆弱性,其蕴含的丰富生物多样性受到极大的威胁,进而也会对该区传粉昆虫及开花植物群落造成严重影响。因此,在后续有关传粉昆虫的保护中,不仅要选取传粉昆虫作为指示物种,也要选取响应的开花或蜜源植物作为指示物种,通过它们的协同或不协同的响应来监测和评估保护区的生态环境变化,为更好地制定资源昆虫保护策略提供科学依据。47第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023参考文献 1 OLLERTON J.Pollinator diversity distribution ecological func-tion and conservation J.Annual Review of Ecology Evolution and Systematics 2017 48 353-376.2 POTTS S G IMPERATRIZ-FONSECA V NGO H T et al.Safe-guarding pollinators and their values to human well-being J.Na-ture 2016 540 7632 220-229.3 杨大荣.传粉昆虫的研究进展 J.云南农业科技 2018 S1 3-7.4 RHODES C J.Pollinator decline-an ecological calamity in the mak-ing J.Science Progress 2018 101 2 121-160.5 DICKS L V BREEZE T D NGO H T et al.A global-scale expert assessment of drivers and risks associated with pollinator decline J.Nature Ecology&Evolution 2021 1453-1461.6 杨培 彭艳琼 赵荣华 等.大蜜蜂的生物学特性、面临威胁与保护策略 J.生物多样性 2018 26 5 476-485.7 FANG L ZHANG Y GAO K et al.Butterfly communities along the Heihe River Basin in Shaanxi Province a biodiversity conserva-tion priority area in China J .Journal of Insect Conservation 2019 23 5 873-883.8 TEICHROEW J L XU J AHRENDS A et al.Is China s unparal-leled and understudied bee diversity at risk J.Biological Conser-vation 2017 210 19-28.9 MIROCHNITCHENKO N A STUBER E F FONTAINE J J.Biodi-versity scale-dependence and opposing multi-level correlations under-lie differences among taxonomic phylogenetic and functional diversi-ty J.Journal of Biogeography 2021 48 12 2989-3003.10 SILVA M NOLL F B CASTRO A C M C.Phylogeographic analy-sis reveals high genetic structure with uniform phenotypes in the pa-per wasp Protonectarina sylveirae Hymenoptera Vespidae J.PLoS One 2018 13 3 e0194424.11 FRANCUSKI L LUDOKI J LUKAC M et al.Fine scale popula-tion structure of hoverfly pollinator Eristalis arbustorum an inte-grative study J.Journal of Insect Conservation 2020 24 1 49-63.12 ZHANG A B HAO M D YANG C Q et al.Barcoding R an in-tegrated R package for species identification using DNA barcodes J.Methods in Ecology and Evolution 2017 8 5 627-634.13 VU H T LE L.Bioinformatics analysis on DNA barcode sequences for species identification a review J.Annual Research&Review in Biology 2019 34 1 1-12.14 LIANG H ZHAO Y H RAFFERTY N E et al.Evolutionary and ecological factors structure a plant-bumblebee network in a biodiver-sity hotspot the Himalaya-Hengduan Mountains J.Functional E-cology 2021 35 11 2523-2535.15 FAITH D P.Systematics and conservation on predicting the fea-ture diversity of subsets of taxa J.Cladistics 1992 8 4 361-373.16 刘乃一 刘子豪 虞磊 等.鹞落坪国家级自然保护区蝶类的物种多样性及区系组成 J.基因组学与应用生物学 2018 37 5 1941-1946.17 王松 鲍方印 梅百茂 等.安徽鹞落坪国家级自然保护区蝶类的垂直分布及其群落多样性 J.应用生态学报 2009 20 9 2262-2270.18 LI W B LIU Y WU Y H et al.Community composition and diversi-ty of ground beetles Coleoptera Carabidae in Yaoluoping National Nature Reserve J.Journal of Insect Science 2017 17 6 1-8.19 SIMON C FRATI F BECKENBACH A et al.Evolution weighting and phylogenetic utility of mitochondrial gene sequences and a compila-tion of conserved polymerase chain reaction primers J.Annals of the Entomological Society of America 1994 87 6 651-701.20 ROZAS J FERRER-MATA A SNCHEZ-DELBARRIO J C et al.DnaSP 6 DNA sequence polymorphism analysis of large data sets J.Molecular Biology and Evolution 2017 34 12 3299-3302.21 NGUYEN L T SCHMIDT H A VON HAESELER A et al.IQ-TREE a fast and effective stochastic algorithm for estimating maxi-mum-likelihood phylogenies J.Molecular Biology and Evolution 2015 32 1 268-274.22 STEUDEL B HALLMANN C LORENZ M et al.Contrasting biodiversity-ecosystem functioning relationships in phylogenetic and functional diversity J.New Phytologist 2016 212 2 409-420.23 MATOS F A R MAGNAGO L F S GASTAUER M et al.Effects of landscape configuration and composition on phylogenetic dive