• 中国科学引文数据库(CSCD)来源期刊
  • 中文核心期刊
  • 中国农林核心期刊(A类)
  • 中国高校百佳科技期刊

基于AMPK/PGC-1α信号通路探讨菟丝子多糖改善大鼠运动能力的作用及机制

庞伊婷, 麻飞

庞伊婷, 麻飞. 基于AMPK/PGC-1α信号通路探讨菟丝子多糖改善大鼠运动能力的作用及机制[J]. 云南农业大学学报(自然科学), 2024, 39(2): 53−59. DOI: 10.12101/j.issn.1004-390X(n).202306031
引用本文: 庞伊婷, 麻飞. 基于AMPK/PGC-1α信号通路探讨菟丝子多糖改善大鼠运动能力的作用及机制[J]. 云南农业大学学报(自然科学), 2024, 39(2): 53−59. DOI: 10.12101/j.issn.1004-390X(n).202306031
PANG Yiting, MA Fei. Effects and Mechanism of Cuscuta chinensis Polysaccharides on Improving Motor Ability in Rats Based on AMPK/PGC-1α Signaling Pathway[J]. JOURNAL OF YUNNAN AGRICULTURAL UNIVERSITY(Natural Science), 2024, 39(2): 53-59. DOI: 10.12101/j.issn.1004-390X(n).202306031
Citation: PANG Yiting, MA Fei. Effects and Mechanism of Cuscuta chinensis Polysaccharides on Improving Motor Ability in Rats Based on AMPK/PGC-1α Signaling Pathway[J]. JOURNAL OF YUNNAN AGRICULTURAL UNIVERSITY(Natural Science), 2024, 39(2): 53-59. DOI: 10.12101/j.issn.1004-390X(n).202306031

基于AMPK/PGC-1α信号通路探讨菟丝子多糖改善大鼠运动能力的作用及机制

基金项目: 新世纪广西高等教育教学改革工程项目(2014JGB140)。
详细信息
    作者简介:

    庞伊婷(1982—),女,广西博白人,硕士,副教授,主要从事体育教学训练学、运动医学和体育人文社会学研究。E-mail:pangyiting2022@163.com

    通信作者:

    麻飞(1993—),男,安徽萧县人,硕士,讲师,主要从事体育教学训练学、运动医学和体育人文社会学研究。E-mail:jingwumen2023@163.com

  • 中图分类号: S567.219.099

摘要:
目的 

基于腺苷酸活化蛋白激酶(AMPK)/过氧化物酶体增殖活化受体γ共激活因子1α (PGC-1α)信号通路探讨菟丝子多糖改善大鼠运动能力的作用及潜在机制。

方法 

将大鼠50只分为正常对照组、运动模型组及菟丝子多糖低、中、高剂量组(10、200和400 mg/kg),除正常对照组外,其余各组大鼠进行游泳训练。游泳至力竭后,测定各组大鼠血清抗疲劳指标以及肝组织能量代谢和氧化应激指标,AMPKαPGC-1α mRNA的表达水平,磷酸化AMPKα (p-AMPKα)、AMPKα和PGC-1α的蛋白表达水平。

结果 

与运动模型组相比,菟丝子多糖低、中、高剂量组大鼠游泳力竭时间显著或极显著延长,血清尿素氮、乳酸、皮质酮以及肝组织丙二醛含量极显著减少,血清睾酮以及肝组织糖原、钠钾ATP酶、钙镁ATP酶、超氧化物歧化酶、谷胱甘肽过氧化物酶水平显著增加,肝组织AMPKα mRNA及其蛋白表达无明显变化,PGC-1α mRNA及p-AMPKα、PGC-1α蛋白表达极显著增加。

结论 

菟丝子多糖可提升高强度大运动量训练大鼠的运动能力,且该作用与抗疲劳、增加能量物质储备、抑制氧化应激及活化AMPK/PGC-1α信号通路等有关。

 

Effects and Mechanism of Cuscuta chinensis Polysaccharides on Improving Motor Ability in Rats Based on AMPK/PGC-1α Signaling Pathway

Abstract:
Purpose 

To study the effects and potential mechanism of Cuscuta chinensis polysaccharides (CCPs) on improving motor ability in rats, based on adenosine monophosphate activated protein kinase (AMPK)/peroxisome proliferator activator receptor γ coactivator 1α (PGC-1α) signaling pathway.

Methods 

Fifty rats were divided into normal control (NC) group, exercise model (EM) group and CCPs low-dose (CCPs-L), medium-dose (CCPs-M) and high-dose (CCPs-H) groups (10, 200 and 400 mg/kg). Except NC group, all rats in the other groups needed swimming training. After swimming to exhaustion, the anti-fatigue indexes in serum, the indexes of energy metabolism and oxidative stress in liver tissue, the mRNA expression levels of AMPKα and PGC-1α, and the protein expression levels of phosphorylated AMPKα (p-AMPKα), AMPKα and PGC-1α were measured.

Results 

Compared with EM group, swimming exhaustion time of CCPs groups was significantly or extremely significantly prolonged; the contents of blood urea nitrogen, lactic acid and corticosterone in serum and malondialdehyde in liver tissue were extremely significantly decreased; the contents of serum testosterone and liver glycogen, and the activities of Na+-K+-ATPase, Ca2+-Mg2+-ATPase, superoxide dismutase and glutathion peroxidase in liver tissue were significantly increased; the mRNA and protein levels of AMPKα in liver tissue had no significant changes, but the mRNA level of PGC-1α and protein levels of p-AMPKα and PGC-1α were extremely significantly increased.

Conclusion 

CCPs can improve the exercise ability of rats trained with high intensity, which is related to anti-fatigue, increasing energy reserves, inhibiting oxidative stress and activating AMPK/PGC-1α signaling pathway.

 

  • 为推进山水林田湖草一体化保护和系统治理,国家层面出台了生态保护和修复领域综合型规划−《全国重要生态系统保护和修复重大工程总体规划(2021—2035年)》(简称“双重”规划),该规划立足国家生态安全战略格局,提出了“三区四带”为核心的全国重要生态系统保护和修复重大工程总体布局。生态系统保护和修复工程侧重项目落地性和实施成效[1],应具备特定空间范围属性[2],科学识别并确定生态保护与修复空间是“双重”规划落地实施的关键环节[2-3]。地方实施生态系统保护和修复工程在宏观尺度上支撑和服务国家生态安全格局重要骨架和关键节点,中观尺度上以流域为系统实施单元,微观尺度上以县域为具体实施单元。相关研究从流域[4-5]、省域[3]、区域[6-7]和县域[8]尺度开展了生态保护与修复分区、识别探索,研究方法多以生态功能重要性、敏感性评价[9-14]和最小阻力模型[6-8, 15-18]为主,为生态保护和修复提供了理论基础和科学借鉴;但仍存在一定不足,如生态功能重要性、敏感性评价方法和评价结果分级分类标准体系不统一,且与国土空间规划体系衔接不紧密,客观上限制了成果落地应用。桂黔滇喀斯特石漠化防治生态功能区被纳入“双重”规划布局,未来将实施一批重大生态保护和修复工程,而有关典型喀斯特地区生态保护与修复空间识别的相关研究尚未见报道。

    云南省昆明市石林彝族自治县(以下简称“石林县”)地处珠江上游,为中国西南喀斯特典型地区,是桂黔滇喀斯特石漠化防治国家重点生态功能区,受温暖潮湿季风气候和特殊地质背景制约[19],生态系统敏感脆弱,区内石漠化、水土流失、植被退化等生态环境问题突出,严重制约了区域可持续发展。按照《云南省国土空间规划(2021—2035年)》和《云南省重要生态系统保护和修复重大工程总体规划(2021—2035年)》布局分区,石林县属滇东滇东南石漠化带,是生态保护和修复重大工程落地实施的重点区域。自然资源环境本底条件客观上决定了生态保护与修复空间布局,石林县自然资源禀赋相对不足势必导致生态保护与修复难度大、成本高,因此,保护与修复应重点关注区域内受损严重、开展修复需求最迫切、恢复效果最明显的重点区域[20]。本研究基于ArcGIS 10.5平台,根据《生态保护红线划定指南》(环办生态〔2017〕48号)技术要求,采用NPP定量指标法,对石林县全域开展生态系统服务功能重要性和生态敏感性评价,结合碳酸盐岩分布、岩溶发育程度及土地利用类型,识别并确定生态保护与修复空间,以期为喀斯特地区生态保护与修复提供科学依据。

    石林县总面积1680.88 km2,地势自东向西倾斜,呈阶梯状逐级下降,区域内石灰岩广布,以石林、石芽、溶洞、洼地、落水洞、地下暗河、盲谷等喀斯特景观为主。全县土地资源匮乏,耕地以零星斑块状散布于峰丛洼地间,山地陡坡垦殖现象普遍[19, 21]。石漠化和水土流失是制约石林县社会经济可持续发展的重要因素,也是生态保护与修复需要克服的首要压力。

    本研究生态系统服务功能重要性和生态敏感性评价采用NPP定量指标评估法,涉及数据主要包括多年平均植被净初级生产力(net primary production,NPPm)、DEM、归一化植被指数(normalized difference vegetation index,NDVI)、土壤、气象、土地利用现状等。其中,多年平均NPPm (分辨率30 m)数据来源于昆明市自然资源局;DEM数据为GDEMV3 (分辨率30 m),来源于地理空间数据云(https://www.gscloud.cn);NDVI数据(分辨率30 m)由Landsat 8 OLI_TIRS (来源于地理空间数据云)计算得到;土壤数据来源于石林县农业农村局;多年平均气温、降水等气象数据来源于WorldClimate (http://www.worldclimate.com)和云南省124个气象站点逐日降雨数据(2000—2021年);土地利用现状数据来源于石林县自然资源局。所有数据投影坐标统一为WGS_1984_UTM_Zone_48N。

    生态系统服务功能重要性和生态敏感性评价方法直接影响生态保护与修复空间识别结果[22]。根据《云南省生态功能区划》,确定研究区生态系统服务功能重要性评价内容主要包括水源涵养、水土保持和生物多样性维护,生态敏感性评价内容主要包括水土流失和石漠化。《生态保护红线划定指南》是指导全国生态保护红线划定的技术文件,其生态系统服务功能重要性、敏感性评价方法体系被《资源环境承载能力和国土空间开发适宜性评价技术指南(试行)》直接引用,并作为指导地方编制国土空间规划和开展生态保护与修复的技术遵循。生态系统服务功能重要性评价一般采用模型评价法[23],但由于该方法所需模型参数多且对数据精度要求较高,难适用于监测体系不完善且监测能力较弱的区域。鉴于生态环境监测实际和数据可获取性,本研究以《生态保护红线划定指南》为技术参考,采用NPP定量指标法对研究区全域开展生态系统服务功能重要性评价,采用评估模型法开展生态敏感性评价。

    (1)水源涵养服务功能重要性评价

    水源涵养是生态系统通过其特有的结构拦截和滞蓄降水,并通过蒸散发实现水资源的调控,其主要功能表现为补充地下水、缓和地表径流、减缓河流流量季节波动等[3],其评价模型为:

    $$ {W}{R}={{\mathrm{NPP}}}_{{\mathrm{m}}}\times {F}_{{\mathrm{s}}}\times {F}_{{\mathrm{p}}}\times \left(1-{F}_{{\mathrm{slo}}}\right) 。 $$

    式中:WR为生态系统水源涵养服务能力指数;NPPm为多年平均净初级生产力;Fs为土壤渗流因子,按照土壤质地类型由重黏土到砂土分别在0~1之间均等赋值得到,如重黏土为1/13,黏土为2/13,粉砂质黏土为3/13,······,砂土为1;Fp为多年平均降水量因子,本文以研究区边界裁剪WordClimate气象数据中的年降水量图层,将其作为研究区多年平均降水量栅格;Fslo为坡度因子,通过DEM坡度分析得到。

    以上各因子均重采样为30 m×30 m的栅格,采用最大最小值法归一化,经栅格计算器利用Int [ (水源涵养服务能力栅格/水源涵养服务能力栅格数据最大值)×100]计算得到生态系统服务值栅格图,计算每个栅格像元的累加服务值,以累加服务值占生态系统服务总值比例的50%和80%为分级分界点,将水源涵养服务功能重要性分为极重要、重要和一般重要3个等级,将极重要和重要区作为生态保护空间候选区域。

    (2)水土保持服务功能重要性评价

    水土保持是生态系统通过其结构与过程减少降雨侵蚀、减少土壤流失、保持土壤的作用,是生态系统重要服务功能之一[3, 23],其评价模型为:

    $$ {S}_{ {\mathrm{pro}}}={{\mathrm{NPP}}}_{{\mathrm{m}}}\times \left(1-K\right)\times \left(1-{F}_{{\mathrm{slo}}}\right) 。 $$

    式中:Spro为土壤保持服务能力指数;K为土壤可蚀性因子,根据《生态保护红线划定指南》模型评价法土壤保持服务功能重要性评价K因子算法,利用石林县土壤质地数据经普通克里金空间插值得到。

    (3)生物多样性维护功能重要性评价

    生物多样性维护功能是生态系统在维持遗传多样性、物种多样性、生态系统多样性方面发挥的作用[3-4],其评价模型为:

    $$ {S}_{ {\mathrm{bio}}}={{\mathrm{NPP}}}_{{\mathrm{m}}}\times {F}_{{\mathrm{p}}}\times {F}_{{\mathrm{t}}}\times \left(1-{F}_{{\mathrm{a}}}\right) 。 $$

    式中:Sbio为生物多样性维护服务能力指数;Fp为多年平均降水量因子,Ft为多年平均气温;Fa为海拔因子。

    (4)生态系统服务功能综合评价

    开展生态系统服务功能重要性评价,有利于分析生态系统服务功能空间分布规律,明确生态服务功能重要区域,为国土空间开发和生态保护提供科学依据[24]。由于无法确定各生态系统服务功能在总生态系统服务功能中的贡献,因此不存在对各生态系统服务功能进行权重选择[23],为客观、中立、综合评价研究区生态系统服务功能,采用生态系统服务功能综合指数(comprehensive indexes of ecosystem services,CIES)作为评价指标,其计算公式为:

    $$ \mathrm{C}\mathrm{I}\mathrm{E}\mathrm{S}=\mathrm{W}\mathrm{R}+{S}_{ {\mathrm{pro}}}+{S}_{ {\mathrm{bio}}} 。 $$

    (1)水土流失敏感性评价

    水土流失按土壤侵蚀发生动力条件可分为水力侵蚀和风力侵蚀,水力侵蚀是造成西南喀斯特地区水土流失的主要原因[19],其评价模型为:

    $$ {S E}_{i}=\sqrt[4]{{R}_{i}\times {K}_{i}\times {LS}_{i}\times {C}_{i}} 。 $$

    式中:SEii空间单元水土流失敏感性指数; Ri为降雨侵蚀力因子,根据《生态保护红线划定指南》水土流失敏感性评估模型,利用云南省124个气象站点逐日降水数据(2000—2021年)分别计算各站点R值,再经协同克里金空间插值(辅助变量为高程)得到;Ki为土壤可蚀性因子;LSi为坡度坡长因子,Ci为植被覆盖度,均由Landsat 8 OLI_TIRS影像提取NDVI计算得到。水土流失敏感性评价因子分级标准见表1

    表  1  水土流失敏感性评价因子分级赋值
    Table  1.  Classification and assignment of soil erosion sensitivity
    指标
    indexes
    降雨侵蚀力
    rainfall erosivity
    土壤可蚀性
    soil erodibility
    地形起伏度
    relief amplitude
    植被覆盖度
    vegetation coverage
    分级赋值
    graded assignment
    一般敏感
    general sensitive
    <100石砾、沙、粗砂土、细砂土、黏土
    sand, loamy sand, sandy loam, sandy clay, clay
    <50≥0.61
    敏感
    sensitive
    100~600面砂土、壤土、沙壤土、粉黏土、壤黏土
    silt loam, loam, sandy clay loam, silt clay loam, clay loam
    50~3000.2~0.63
    极敏感
    extremely sensitive
    >600砂粉土、粉土
    silt loam, silt
    >300≤0.25
    下载: 导出CSV 
    | 显示表格

    (2)石漠化敏感性评价

    石漠化导致生态系统稳定性减弱、敏感性增强,影响区域生态系统服务和产品供给,开展石漠化敏感性评价可定量描述石漠化对人类活动的敏感程度,《生态保护红线划定指南》推荐评价模型为:

    $$ {S}_{ i}=\sqrt[3]{{D}_{i}\times {P}_{i}\times {C}_{i}} 。 $$

    式中:Sii空间单元石漠化敏感性指数;Di为碳酸岩出露面积百分比,通过遥感影像解译计算求得;Pi为地形坡度,利用DEM开展坡度计算得到。石漠化敏感性评价因子分级标准见表2

    表  2  石漠化敏感性评价因子分级赋值
    Table  2.  Classification and assignment of stony desertification sensitivity
    指标
    indexes
    碳酸岩出露
    面积百分比/%
    percentage of
    exposed area of carbonate rock
    地形
    坡度/(°)
    terrain slope
    植被覆盖度
    vegetation coverage
    分级赋值
    graded assignment
    一般敏感
    general sensitive
    ≤30≤8≥0.61
    敏感
    sensitive
    30~708~250.2~0.63
    极敏感
    extremely sensitive
    ≥70≥25≤0.25
    下载: 导出CSV 
    | 显示表格

    已有研究表明:建设用地、水体、居民地及道路在遥感影像解译提取研究区碳酸盐岩出露面积百分比时易产生误判[25]。石漠化是自然背景与人类活动共同作用下形成的一种特殊地质现象,岩性是发生石漠化的基础因素,而遥感影像解译只能提取出露地表的碳酸盐岩,忽视了掩藏于植被和土层之下的主体基岩,无法整体反映典型喀斯特地区岩性基底。为弥补以上不足,本研究在依据《生态保护红线划定指南》推荐方法开展石漠化敏感性评价的同时,基于《1∶20万宜良幅区域水文地质普查报告》,按照研究区碳酸盐岩的实际分布范围及其岩溶发育强弱程度(强、中、弱)开展石漠化敏感性评价分级,分为极敏感、敏感、一般敏感3个等级。碳酸盐岩岩溶发育强弱程度与石漠化敏感性对应关系及其赋值规则见表3

    表  3  碳酸盐岩岩溶发育程度与石漠化敏感性对应关系
    Table  3.  Corresponding relationship between the degree of carbonate rock development and the sensitivity of rocky desertification
    岩溶发育等级
    level of karst development
    主要碳酸盐岩地层时代
    main carbonate rock strata age
    敏感性
    sensitivity
    分级赋值
    graded assignment
    地质时代
    geological era
    地层
    strata

    weak
    泥盆系
    Devonian
    一打得组
    Yidade formation
    一般敏感
    general sensitive
    1
    石炭系下统
    the lower Carboniferous
    万寿山组
    Wanshoushan formation

    medium
    寒武系下统
    the lower Cambrian
    沧浪铺组
    Canglangpu formation
    敏感
    sensitive
    3
    筇竹寺组
    Qiongzhusi formation
    龙王庙组
    Longwangmiao formation
    寒武系上统
    the upper Cambrian
    陡坡寺组
    Douposi formation
    双龙潭组
    Shuanglongtan formation
    志留系上统
    the upper Silurian
    玉龙寺组
    Yulongsi formation
    泥盆系中统
    the middle Devonian
    曲靖组
    Qujing formation

    strong
    震旦系
    Sinian
    灯影组
    Dengying formation
    极敏感
    extremely sensitive
    5
    石炭系下统
    the lower Carboniferous
    上司段组
    Shangsiduan formation
    石炭系中统
    the middle Carboniferous
    威宁群
    Weining group
    石炭系上统
    the upper Carboniferous
    马平组
    Maping formation
    二叠系下统
    the lower Permian
    阳新组
    Yangxin formation
    下载: 导出CSV 
    | 显示表格

    (3)生态敏感性综合评价

    生态敏感性是指区域出现生态环境问题的可能性和可能程度[11, 24, 26],开展生态敏感性评价是对研究区各类生态敏感区进行识别并将潜在生态问题聚焦落地到具体空间区域的过程[24, 27-28],是区域开展生态修复的基础,其评价模型为:

    $$ \mathrm{S}\mathrm{S}\mathrm{I}={S E}_{i}+{S}'_{i} 。 $$

    式中:SSI为生态敏感性综合指数;$ {S}'_{i} $为石漠化敏感性指数,由碳酸盐岩岩溶发育强弱等级划分赋值得到。

    生态系统服务功能重要性和生态敏感性可直观表征区域生态系统质量和生态产品供给能力,利用ArcGIS软件的Union工具将研究区生态系统服务功能重要性、生态敏感性图层联合,可得出生态系统服务功能重要性和生态敏感性分区图。根据《自然生态空间用途管制办法(试行)》中自然生态空间的有关定义,提取土地利用数据中对应森林、草地、河流、湖泊、滩涂和荒地的地类作为研究区自然生态空间;以自然生态空间为边界,裁剪生态系统服务功能重要性、生态敏感性分区结果可得到自然生态空间管控分区。

    山水林田湖草生命共同体是特定条件下不同类型生态系统在空间范围内相互作用形成的空间组合[2],识别并确定生态保护与修复空间旨在保护和提升现有生态系统服务功能[3, 29-30]。区域生态系统服务功能重要性越高,生态敏感性越高,则其对国土空间开发的约束力越强[23],识别并确定生态保护空间和生态修复空间是区域生态系统服务功能“保优提质”的重要前提。生态系统服务功能极重要区和重要区是维系区域生态系统结构完整、生态服务功能持续稳定的关键区,是优质生态服务供给区,也是生态保护工程实施的重点区域;生态极敏感区和敏感区是保障区域生态系统结构完整、生态服务功能不降低的热点区,是生态服务提质区,也是生态修复工程落地的核心区域[27, 31-34]

    通过对研究区进行生态系统服务功能重要性和生态敏感性评价,将生态系统服务功能综合指数划分为极重要、重要和一般重要,将生态敏感性综合指数划分为极敏感、敏感和一般敏感,构建生态系统服务功能重要性和生态敏感性交叉判断矩阵(表4),根据各栅格单元生态系统服务功能重要性和生态敏感性组合,将自然生态空间管控分区划分为三类空间,即生态保护区、生态修复区、生态保护+生态修复区。

    表  4  生态系统服务功能重要性和生态敏感性交叉判断矩阵
    Table  4.  Crossed judgment matrix of ecological service function importance and ecological sensitivity
    敏感性
    sensitivity
    重要性 importance
    一般重要区
    general important area
    重要区
    important area
    极重要区
    extremely important area
    一般敏感区
    general sensitive area
    无需管控
    no need to control
    生态保护
    ecological protection
    生态保护
    ecological protection
    敏感区
    sensitive area
    生态修复
    ecological restoration
    生态保护+生态修复
    ecological protection and restoration
    生态保护+生态修复
    ecological protection and restoration
    极敏感区
    extremely sensitive area
    生态修复
    ecological restoration
    生态保护+生态修复
    ecological protection and restoration
    生态保护+生态修复
    ecological protection and restoration
    注:表中内容为生态系统服务功能重要性和生态敏感性分区两两交叉区域的管控策略。
    Note: The table contains are the control strategies of two intersecting areas of ecosystem service function importance and ecological sensitivity.
    下载: 导出CSV 
    | 显示表格

    研究区水源涵养功能重要性评价结果(图1a)总体为一般重要,其中,极重要区和重要区零星分布于西街口镇北部和鹿阜街道北部降雨量和植被覆盖度较高的区域,面积共计1.39 km2,仅占研究区面积的0.08%。水土保持功能重要性评价结果(图1b)总体为一般重要,其空间分布格局与水源涵养功能重要性非常相似,极重要区和重要区集中连片分布于鹿阜街道,零星分布于板桥、大可、长湖等6个乡镇(街道),面积共计7.81 km2,占比仅为0.46%,这主要是因为研究区为典型喀斯特地区,长期岩溶作用形成了特有的地表、地下立体双层结构,该结构渗透强烈,造成地表非地带性干旱[35],加之碳酸盐岩风化成土较慢[19],降雨侵蚀加剧土壤侵蚀,导致水源涵养和水土保持功能较弱。生物多样性维护功能重要性评价结果(图1c)显示:极重要区主要分布于圭山国家森林公园,重要区主要分布于西街口、圭山两镇交界植被覆盖度较高的杨梅山山区,极重要区和重要区面积共计24.35 km2,占研究区面积的1.45%。为系统表征研究区生态系统服务功能,将水源涵养、水土保持和生物多样性维护3种服务功能评价结果进行空间叠加计算并分级,结果(图1d)显示:研究区生态系统服务功能重要性总体一般,极重要区主要位于鹿阜街道北部区域和圭山国家森林公园一带,重要区分布范围与生物多样性维护功能重要区基本重叠,生态系统服务功能极重要区和重要区面积共计32.18 km2,占研究区总面积的1.91%,这主要是因为研究区特殊的喀斯特生态系统异常敏感脆弱,严重制约了其生态系统调节服务功能。

    图  1  生态系统服务功能重要性评价
    Figure  1.  Evaluation of the importance of ecological services

    水土流失敏感性评价结果(图2a)显示:研究区水土流失敏感性呈现明显的东西差异,即中部和东部水土流失敏感性明显高于西部,极敏感区和敏感区面积共计1414.08 km2,占研究区总面积的84.13%。石漠化敏感性评价结果(图2b)显示:石漠化极敏感区主要位于西街口镇,这与实际情况不符,主要是因为建设用地、裸地、耕地、道路等在提取碳酸盐岩时会产生误判,进而导致区域碳酸盐岩出露百分比可能高于碳酸盐岩出露面积较大的区域。因此,根据碳酸盐岩分布和岩溶发育强弱程度开展石漠化敏感性评价分级,结果(图2c)显示:极敏感区和敏感区广泛分布于研究区全境,面积共计1256.81 km2,占研究区总面积的 74.77%。为系统表征研究区生态敏感性,将水土流失敏感性和石漠化敏感性进行空间叠加计算并进行分级,结果(图2d)显示:研究区生态敏感性较高,极敏感区连片分布于中北部和西南部、斑块状分布于西部区域,敏感区广泛分布于大部分区域;极敏感区和敏感区面积共计1597.62 km2,占研究区总面积的95.05%。

    图  2  生态敏感性评价
    Figure  2.  Evaluation of ecological sensitivity

    生态系统服务功能重要性和生态敏感性分区见图3a,经土地利用数据提取的研究区自然生态空间见图3b。以自然生态空间为边界,裁剪生态系统服务功能重要性、生态敏感性分区结果得到自然生态空间管控分区(图3c)。

    图  3  自然生态空间管控分区
    Figure  3.  Natural ecological space control zones

    按照管控策略(表4),可进一步将自然生态空间(图3)划分为生态保护区、生态修复区、生态保护+生态修复区(图4)。生态保护区主要分布于石林街道中部、鹿阜街道北部、西街口镇中部及板桥街道东南部与大可乡接壤的山区,面积17.95 km2,占研究区自然生态空间管控分区的2.05%,占研究区总面积的1.07%。生态保护区水源涵养、水土保持、生物多样性维护等生态系统服务功能较其他区域更为显著,生态敏感性一般,是生态保护工程实施的重点区域。按照《自然生态空间用途管制办法(试行)》要求,该区应严格控制各类开发利用活动对生态空间的占用和扰动,确保生态系统结构稳定,生态功能不降低。生态修复区广泛分布于研究区大部区域,面积833.26 km2,占研究区自然生态空间管控分区的94.98%,占研究区总面积的49.57%。该区水源涵养、水土保持、生物多样性维护等生态系统服务功能一般,但生态敏感性较高,石漠化和水土流失并存,是生态修复工程布局的重点区域。生态保护+生态修复区主要分布于鹿阜街道北部、圭山国家森林公园及圭山镇与西街口北部接壤的杨梅山山区,面积26.08 km2,占研究区自然生态空间管控分区的2.97%,占研究区总面积的1.55%。该区生态系统服务功能显著且敏感脆弱,需同时实施生态保护工程和生态修复工程。

    图  4  生态保护与修复空间分区
    Figure  4.  Ecological protection and restoration spatial zoning

    以往生态修复相关研究多以单要素或单功能修复为主[3],缺乏从系统性和完整性角度识别生态保护与修复空间,注重单一工程典型示范但对生态系统认识不足,生态修复措施缺乏系统应对性[1]。当前,生态保护与修复的主体是山水林田湖草生命共同体[2],核心目标是通过人为干预维持生态系统结构稳定和提升服务功能。喀斯特生态系统是受特殊地质背景、水文条件及人为干扰多重因素制约的生态系统,其生态承载力低,对国土资源开发的干扰高度敏感,基于研究区全域生态系统服务功能重要性和生态敏感性评估识别生态保护与生态修复空间的方法,遵循山水林田湖草一体化保护和系统治理思路,综合考虑生态系统服务功能重要性和生态敏感性,将研究区生态保护与修复的重点聚焦于自然生态空间,较传统方法具有更强的系统性和针对性,为喀斯特地区生态环境保护和国土空间开发提供科学支撑。

    已有相关研究[15-18]在开展生态系统服务功能重要性和生态敏感性评价时,选取的评价方法各不相同,评估分级标准不统一。《生态保护红线划定指南》有关生态系统服务功能重要性、敏感性评价标准体系作为国土空间规划编制和开展生态保护与修复的技术遵循,在评价方法和评价分级方面作了具体要求。本研究以此作为参照,开展生态系统服务功能重要性和生态敏感性评价,发现圭山国家森林公园位于自然生态空间管控分区的极重要、极敏感区,充分说明基于生态系统服务功能重要性和生态敏感性评价识别生态保护与修复空间方法的科学性和适用性。本研究还发现《生态保护红线划定指南》推荐的石漠化敏感性评价模型有其局限性,采用碳酸盐岩出露百分比作为重要评价因子,关注了出露地表的碳酸盐岩,但对喀斯特地区岩性基底考虑不足,该模型基于遥感影像解译提取出露碳酸盐岩的方法易受遥感影像分辨率限制和地物光谱反射率不同而导致分类结果误判[25]。喀斯特地区石漠化差异程度主要受区域地质构造影响,本研究结合水文地质资料,综合考虑碳酸盐岩实际分布范围和岩溶发育程度,弥补了石漠化敏感性评价模型设计的不足。

    与国土空间规划体系衔接,是生态保护与修复空间识别研究成果落地实施的重要前提[1]。为确保生态保护与修复空间分区结果可用于地方生态保护和修复实际,本研究构建了生态系统服务功能重要性和生态敏感性交叉判断矩阵,将生态系统服务功能综合评价结果、生态敏感性综合评价结果与研究区自然生态空间进行空间叠加,将生态保护与修复空间重点聚焦到自然生态空间,对应市县国土空间规划分区的生态保护与保留区,该区核心管控目标为保护与修复,因此,本研究确定的生态保护与修复空间同国土空间规划实施保护与修复的空间范围一致,有利于指导生态保护与生态修复工程落地实施。

    研究区生态保护与修复三类空间中,生态修复区分布范围最广、面积最大,是生态保护和修复工程重点布局和实施的区域。岩性基底是影响岩溶地区生态环境的关键因素[36],由于生态修复区碳酸盐岩广泛发育,岩性较脆弱,易遭风化溶蚀,地表发育出宽大的溶蚀沟槽及溶蚀裂隙、落水洞、洼地等,地下发育出规模较大的溶洞、暗河管道系统,地表土壤随雨水运移进入地下溶洞暗河,加剧水土流失。因此,该区应严格禁止布局与生态修复目标不符的建设项目,特别是重污染涉水项目;同时,应严格管控农业种植,逐步减少化肥农药使用,避免造成地下水污染;大力实施退耕还林还草,提升森林质量,提高区域生态服务价值。生态修复工程实施应开展生态系统服务权衡分析,综合考虑用地类型、自然地理单元连续性和完整性,确定生态修复总体目标,确保主导生态服务功能及其他功能协同发挥[22],结合区域主要生态环境问题细致划分修复单元,进而确定生态修复单元目标和差异化修复措施,探索关键要素控制、全要素耦合的系统治理成效考核指标,制订基于山水林田湖草系统治理的修复方案,确保生态投资回报率在合理区间,提升生态修复成效。

    本研究为喀斯特地区生态保护与修复空间识别提供了科学参考,然而,研究成果能否应用于具体生态保护和修复决策,还需进一步科学论证。生态修复的对象本质上是生态系统,直观上则是作用于景观,因此在对自然生态空间进行人工修复和自然修复的同时,可借鉴景观生态学“斑块—廊道—基质”理论,以生态保护与修复空间为斑块,科学构建生态廊道,发挥非自然生态空间作为基质的景观动态主导作用,提升生态保护和修复成效。生态保护工程主要以生物多样性维护为主,科学建立和完善生态廊道,维持和提高生境斑块的景观连通性。生态修复工程以提升森林质量为基础,坚持保护优先,自然修复为主,开展人工促进、天然更新等措施。

    本研究基于碳酸盐岩分布和岩溶发育强弱程度设计的石漠化敏感性评价方法,弥补了《生态保护红线划定指南》石漠化敏感性评价模型在水文地质因素方面设计的不足;采用生态系统服务功能综合评价和生态敏感性综合评价相结合的方法,衔接国土空间规划自然生态空间,识别了生态保护区、生态修复区、生态保护+生态修复区三类空间,面积分别为17.95、833.26和26.08 km2,为研究区生态保护和修复工程布局和实施提供了落地空间参考;明确了研究区生态保护和修复工程的重点方向,针对研究区资源环境禀赋现状,从治理理念和修复策略探讨了生态修复基本思路,为确保和提升生态修复成效提供了思路。

  • 图  1   各组大鼠肝组织p-AMPKα、AMPKα及PGC-1α蛋白的表达

    注:NC. 正常组,EM. 运动模型组,CCPs-L. 菟丝子多糖低剂量组,CCPs-M. 菟丝子多糖中剂量组,CCPs-H. 菟丝子多糖高剂量组。

    Figure  1.   Expression levels of p-AMPKα, AMPKα and PGC-1α proteins in liver tissue of rats in each group

    Note: NC. normal control group, EM. exercise model group, CCPs-L. Cuscuta chinensis polysaccharides low-dose group, CCPs-M. C. chinensis polysaccharides medium-dose group, CCPs-H. C. chinensis polysaccharides high-dose group; p-AMPKα. phosphorylated AMPKα, AMPKα. adenosine monophosphate activated protein kinase α, PGC-1α. peroxisome proliferator activator receptor γ coactivator 1α, GAPDH. glyceraldehyde phosphate dehydrogenase.

    表  1   实时定量PCR引物信息

    Table  1   Primer information of real-time PCR

    基因名称
    gene name
    引物序列 (5′→3′)
    primer
    sequence
    扩增产物长度/bp
    length of amplification product
    AMPKαF:GGAGAAAGATGGCGACG124
    R:AACTCGTGCTTGCCCAC
    PGC-1αF:TGTTCCCGATCACCATATTCC219
    R:GTATTCATCCCTCTTGAGCCTTT
    GAPDHF:GGTCTCCTGCGACTTCA133
    R:TGGTCCAGGGTTTCTTACTCC
    下载: 导出CSV

    表  2   菟丝子多糖对运动模型大鼠体质量的影响

    Table  2   Effects of Cuscuta chinensis polysaccharides on the body weight of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    体质量/g body weight
    干预前 before intervention 干预后 after intervention
    NC 245.36±18.07 407.32±30.38
    EM 237.90±15.86 374.83±38.13*
    CCPs-L 100 240.85±22.75 372.57±29.76
    CCPs-M 200 242.31±20.12 366.52±32.24
    CCPs-H 400 238.43±17.30 369.74±43.18
    注:NC. 正常组,EM. 运动模型组,CCPs-L. 菟丝子多糖低剂量组,CCPs-M. 菟丝子多糖中剂量组,CCPs-H. 菟丝子多糖高剂量组;与正常对照组相比,*P<0.05,**P<0.01;与运动模型组相比,#P<0.05,##P<0.01;n=10;下同。
    Note: NC. normal control group, EM. exercise model group, CCPs-L. C. chinensis polysaccharides low-dose group, CCPs-M. C. chinensis polysaccharides medium-dose group, CCPs-H. C. chinensis polysaccharides high-dose group; compared with the normal control group, *P<0.05, **P<0.01; compared with the exercise model group, #P<0.05, ##P<0.01; the same as below.
    下载: 导出CSV

    表  3   菟丝子多糖对运动模型大鼠游泳力竭时间的影响

    Table  3   Effects of C. chinensis polysaccharides on the swimming exhaustion time of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    游泳力竭时间/min
    swimming exhaustion time
    NC
    EM18.75±3.22
    CCPs-L10021.93±2.36#
    CCPs-M20024.18±3.89##
    CCPs-H40024.64±2.07##
    下载: 导出CSV

    表  4   菟丝子多糖对运动模型大鼠血清抗疲劳指标的影响

    Table  4   Effects of C. chinensis polysaccharides on the serum anti-fatigue indexes of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    尿素氮含量/(ng·mL−1)
    content of BUN
    乳酸含量/(ng·mL−1)
    content of LA
    皮质酮含量/(ng·mL−1)
    content of CTC
    睾酮含量/(ng·mL−1)
    content of TTS
    NC162.35±11.970.46±0.0355.39±4.6112.44±1.37
    EM250.26±24.83**1.05±0.12**98.23±10.35**5.72±0.62**
    CCPs-L100219.61±18.23##0.81±0.09##70.43±8.26##6.59±0.80#
    CCPs-M200204.74±17.56##0.64±0.05##63.10±5.84##7.41±0.55##
    CCPs-H400167.90±15.31##0.62±0.07##56.06±7.53##7.57±1.08##
    注/Note: BUN. blood urea nitrogen; LA. lactic acid; CTC. corticosterone; TTS. testosterone.
    下载: 导出CSV

    表  5   菟丝子多糖对运动模型大鼠肝组织能量代谢指标的影响

    Table  5   Effects of C. chinensis polysaccharides on the indexes of energy metabolism in liver tissue of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    糖原含量/(mg·g−1)
    content of glycogen
    钠钾ATP酶活性/(μmol·g−1·h−1)
    activity of Na+-K+-ATPase
    钙镁ATP酶活性/(μmol·g−1·h−1)
    activity of Ca2+-Mg2+-ATPase
    NC 8.23±0.68 1.26±0.14 0.98±0.07
    EM 2.47±0.31** 0.80±0.09** 0.67±0.08**
    CCPs-L 100 4.86±0.53## 0.95±0.11## 0.74±0.05#
    CCPs-M 200 4.60±0.44## 1.13±0.15## 0.83±0.10##
    CCPs-H 400 6.53±0.75## 1.21±0.13## 0.89±0.09##
    下载: 导出CSV

    表  6   菟丝子多糖对运动模型大鼠肝组织氧化应激指标的影响

    Table  6   Effects of C. chinensis polysaccharides on the indexes of oxidative stress in liver tissue of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    丙二醛含量/(nmol·mg−1)
    content of MDA
    超氧化物歧化酶活性/(U·mg−1)
    activity of SOD
    谷胱甘肽过氧化物酶活性/(U·mg−1)
    activity of GSH-Px
    NC7.21±0.89319.51±28.30161.42±14.29
    EM14.70±1.67**204.08±19.72**95.78±8.62**
    CCPs-L10011.53±1.28##257.46±22.78##126.94±13.57##
    CCPs-M20010.04±1.32##276.33±30.24##145.26±15.93##
    CCPs-H4008.35±0.63##288.17±26.51##147.60±11.02##
    注/Note: MDA. malondialdehyde; SOD. superoxide dismutase; GSH-Px. glutathion peroxidase.
    下载: 导出CSV

    表  7   菟丝子多糖对运动模型大鼠肝组织AMPKαPGC-1α mRNA表达的影响

    Table  7   Effects of C. chinensis polysaccharides on the levels of AMPKα and PGC-1α mRNAs in liver tissue of exercise model rats

    分组
    groups
    剂量/
    (mg·kg−1)
    dosage
    腺苷酸活化
    蛋白激酶α
    AMPKα
    过氧化物酶体增殖活化
    受体γ共激活因子1α
    PGC-1α
    NC1.04±0.031.02±0.05
    EM0.98±0.110.33±0.04**
    CCPs-L1001.09±0.150.64±0.06##
    CCPs-M2001.03±0.120.87±0.10##
    CCPs-H4000.96±0.140.93±0.08##
    注/Note: AMPKα. adenosine monophosphate activated protein kinase α, PGC-1α. peroxisome proliferator activator receptor γ coactivator 1α; 下同/the same as below.
    下载: 导出CSV

    表  8   菟丝子多糖对运动模型大鼠肝组织p-AMPKα、AMPKα及PGC-1α蛋白表达的影响

    Table  8   Effects of C. chinensis polysaccharides on the levels of p-AMPKα, AMPKα and PGC-1α proteins in liver tissue of exercise model rats

    分组
    groups
    剂量/(mg·kg−1)
    dosage
    磷酸化腺苷酸活化蛋白激酶α/
    腺苷酸活化蛋白激酶α
    p-AMPKα/AMPKα
    腺苷酸活化蛋白激酶α/
    磷酸甘油醛脱氢酶
    AMPKα/GAPDH
    过氧化物酶体增殖活化受体γ共激活因子1α/
    磷酸甘油醛脱氢酶
    PGC-1α/GAPDH
    NC 0.75±0.06 0.82±0.05 0.91±0.10
    EM 0.14±0.03** 0.79±0.09 0.32±0.05**
    CCPs-L 100 0.31±0.04## 0.76±0.07 0.45±0.04##
    CCPs-M 200 0.35±0.06## 0.81±0.06 0.53±0.08##
    CCPs-H 400 0.54±0.05## 0.83±0.08 0.58±0.06##
    注/Note: p-AMPKα. phosphorylated AMPKα; GAPDH. glyceraldehyde phosphate dehydrogenase.
    下载: 导出CSV
  • [1]

    NILSSON M I, MIKHAIL A, LAN L, et al. A five-ingredient nutritional supplement and home-based resistance exercise improve lean mass and strength in free-living elderly[J]. Nutrients, 2020, 12(8): 2391. DOI: 10.3390/nu12082391.

    [2]

    WITHEE E D, TIPPENS K M, DEHEN R, et al. Effects of methylsulfonylmethane (MSM) on exercise-induced oxidative stress, muscle damage, and pain following a half-marathon: a double-blind, randomized, placebo-controlled trial[J]. Journal of the International Society of Sports Nutrition, 2017, 14(1): 24. DOI: 10.1186/s12970- 017-0181-z.

    [3] 李洪一, 蔡协清, 班允强. 鸡枞菌蛋白提取及对运动大鼠肌肉增长和运动能力的影响[J]. 中国食品添加剂, 2021, 32(12): 121. DOI: 10.19804/j.issn1006-2513.2021.12.016.
    [4]

    JIANG F, YANG R F, XUE D Y, et al. Effects of a natural nutritional supplement on immune cell infiltration and immune gene expression in exercise-induced injury[J]. Frontiers in Nutrition, 2022, 9: 987545. DOI: 10.3389/fnut.2022.987545.

    [5] 徐维茵, 董坤园, 李世昌, 等. 炒菟丝子工艺优选及质量标准研究[J]. 时珍国医国药, 2021, 32(5): 1130. DOI: 10. 3969/j.issn.1008-0805.2021.05.09.
    [6] 王尧尧, 王蕾, 戚莹雪, 等. 菟丝子药材化学成分研究进展[J]. 山东中医药大学学报, 2020, 44(6): 705. DOI: 10.16294/j.cnki.1007-659x.2020.06.022.
    [7] 郭爱民, 曹建民, 朱静, 等. 菟丝子对大鼠抗运动性疲劳能力及脑组织自由基的影响[J]. 中国实验方剂学杂志, 2013, 19(9): 274. DOI: 10.11653/syfj2013090274.
    [8] 燕宇真, 王慧春, 曾阳. 不同炮制方法对菟丝子总黄酮和多糖含量及抗氧化能力的影响[J]. 农产品加工, 2019, 2(3): 53. DOI: 10.16693/j.cnki.1671-9646(X).2019.02.0 14.
    [9] 徐先祥, 李道中, 彭代银, 等. 菟丝子多糖改善糖尿病大鼠糖脂代谢作用[J]. 中国实验方剂学杂志, 2011, 17(18): 232. DOI: 10.13422/j.cnki.syfjx.2011.18.071.
    [10] 蓝瑞高, 梁益军. 人参皂苷CK对力竭游泳大鼠抗疲劳作用及骨骼肌氧化应激的影响[J]. 云南农业大学学报(自然科学), 2022, 37(3): 491. DOI: 10.12101/j.issn.1004- 390X(n).202109040.
    [11] 肖红. 复方五味子提取物对力竭游泳大鼠的抗疲劳作用[J]. 基因组学与应用生物学, 2020, 39(2): 746. DOI: 10.13417/j.gab.039.000746.
    [12]

    CHU H B, ZHANG Z H, CHEN D D, et al. Content determination of phenylpropanoids and enhancing exercise ability of effective fractions in Pedicularis densispica[J]. Pharmacognosy Magazine, 2017, 13(50): 230. DOI: 10.4103/0973-1296.204567.

    [13]

    ZHANG J T, ZHANG C, LIU A F, et al. Synthesis of icariin-zinc and its protective effect on exercise fatigue and reproductive system related glands in male rats[J]. Frontiers in Pharmacology, 2021, 12: 611722. DOI: 10.3389/fphar.2021.611722.

    [14] 尹爱武, 田润, 王盼. 菟丝子粗提物抗疲劳作用研究[J]. 食品工业科技, 2012, 33(14): 164. DOI: 10.13386/j.issn10 02-0306.2012.14.014.
    [15]

    HARGREAVES M, SPRIET L L. Skeletal muscle energy metabolism during exercise[J]. Nature Metabolism, 2020, 2(9): 817. DOI: 10.1038/s42255-020-0251-4.

    [16] 王冬. 黑灵芝多糖对运动训练大鼠抗氧化水平与运动能力的影响[J]. 中国食用菌, 2020, 39(9): 65. DOI: 10.13629/j.cnki.53-1054.2020.09.016.
    [17]

    KRUK J, ABOUL-ENEIN B H, DUCHNIK E. Exercise-induced oxidative stress and melatonin supplementation: current evidence[J]. Journal of Physiological Sciences, 2021, 71(1): 27. DOI: 10.1186/s12576-021-00812-2.

    [18]

    WU F N, LI Z, CAI M X, et al. Aerobic exercise alleviates oxidative stress-induced apoptosis in kidneys of myocardial infarction mice by inhibiting ALCAT1 and activating FNDC5/irisin signaling pathway[J]. Free Radical Biology Medicine, 2020, 158(10): 171. DOI: 10.10 16/j.freeradbiomed.2020.06.038.

    [19] 祁小妮, 隋英, 李振亮, 等. 响应曲面优化菟丝子多糖提取工艺及抗氧化活性研究[J]. 中国酿造, 2015, 34(8): 35. DOI: 10.11882/j.issn.0254-5071.2015.08.007.
    [20] 叶春林, KHUDOYBERDIEV I, 陈颖, 等. 菟丝子多糖的抗氧化活性和抑制肿瘤细胞增殖的研究[J]. 河南工业大学学报(自然科学版), 2020, 41(5): 73. DOI: 10.16433/j.1673-2383.2020.05.012.
    [21]

    DRAKE J C, WILSON R J, LAKER R C, et al. Mitochondria-localized AMPK responds to local energetics and contributes to exercise and energetic stress-induced mitophagy[J]. Proceedings of the National Academy of Sciences of the United States of America, 2021, 118(37): e2025932118. DOI: 10.1073/pnas.2025932118.

    [22]

    THIRUPATHI A, SOUZA C T. Multi-regulatory network of ROS: the interconnection of ROS, PGC-1 alpha, and AMPK-SIRT1 during exercise[J]. Journal of Physiology and Biochemistry, 2017, 73(4): 487. DOI: 10.1007/s13105-017-0576-y.

    [23] 井宏颖, 吕克宁, 宋晓晨, 等. 地黄饮子对运动性疲劳小鼠运动能力的影响[J]. 湖南中医药大学学报, 2021, 41(1): 34. DOI: 10.3969/j.issn.1674-070X.2021.01.007.
    [24] 林致辉, 周庆莹, 王梦妮, 等. 小建中汤对运动性疲劳小鼠骨骼肌AMPK/PGC-1α信号通路的影响[J]. 中国实验方剂学杂志, 2020, 26(13): 73. DOI: 10.13422/j.cnki.syfjx.20201337.
  • 期刊类型引用(1)

    1. 周豹,陈国平,赵俊三,尹莹,于祖国. 滇东南喀斯特地区生态系统服务时空演变及权衡协同关系. 水利水电技术(中英文). 2024(11): 123-142 . 百度学术

    其他类型引用(0)

图(1)  /  表(8)
计量
  • 文章访问数:  813
  • PDF下载量:  9
  • 被引次数: 1
出版历程
  • 通信作者:  麻飞 jingwumen2023@163.com
  • 收稿日期:  2023-06-20
  • 修回日期:  2024-03-06
  • 网络首发日期:  2024-04-24

目录

/

返回文章
返回