Eucommia bark/leaf extract improves HFD-induced lipid metabolism disorders via targeting gut microbiota to activate the Fiaf-LPL gut-liver axis and SCFAs-GPR43 gut-fat axis
本研究探讨了杜仲皮（EBE）、叶（ELE）调节肠道菌群以改善脂质代谢疾病的潜在机制。建立了高脂小鼠模型，并使用粪便微生物群移植（FMT）实验来靶向高脂血症的潜在关键微生物群。确定了肠道菌群改变与脂质代谢之间的因果关系，并且可以提供一种可行的高脂血症治疗方法。此外，获得的结果阐明了EBE和ELE的抗高脂作用的强度，以及杜仲叶是否可以代替杜仲树皮发挥抗高脂血症作用，从而可持续利用杜仲。
实验设计：
1.杜仲皮叶分别制成提取液，和绿原酸分别喂给第一部分不同饲养方式8组的小鼠，4周后每天收集新鲜粪便制成FMT喂给第二部分小鼠。
2.第二部分小鼠喂养FMT8周后，取材并分析小鼠脂代谢数据、脂代谢相关基因以及肠道微生物群变化。
实验流程：
1. 杜仲树皮和杜仲叶分别粉碎后以 1：10的比例在水中煮沸 4 h，将汤剂过滤、收集并干冻干燥，分别得到EBE和ELE。
2. 第一部分：C57BL/6 J小鼠分为8组NC、HFD、NC+CGA、NC+EBE、NC+ELE、HFD+CGA、HFD+EBE 和 HFD+ELE。饲喂第4周后，每天收集小鼠新鲜粪便加入PBS，将混合物离心获得粪便细菌悬浮液，其收集作为实验第 2 部分的供体粪便细菌溶液。
第二部分：C57BL/6 J小鼠高脂喂养并接受8周的第一部分的FMT，形成对应的8组：（HFD-（NC）FMT、HFD-（HFD）FMT、HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT、HFD-（NC+ELE）FMT、HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT 和 HFD-（HFD+ELE）FMT。每周测试一次小鼠体重和饲喂量，每4周测量一次小鼠体长。在第4周和第8周分别进行尾血采集，用于小鼠空腹血糖试验和口服葡萄糖耐量试验。第8周后，采集血液、肝脏、结肠、附睾脂肪、盲肠和粪便样本。
图1：实验设计与时间表
3. 进行组织学分析、生化分析、SCFA 分析、RNA 分离和实时定量 PCR （qRT-PCR）、肠道菌群分析。
4. 用ThermoQ-Autonomous 组合四极杆 Orbitrap 质谱仪与加热电喷雾电离（Thermo Fisher Scientific）确定 EBE 和 ELE 的组成。
实验结果：
图2：
图 2.EBE 和 ELE 的总离子通量模式。（A） EBE 的总离子通量图;（B） ELE 的总离子通量图。
EBE包含18个单体，包括7 种环烯醚类化合物、4 种酚酸、3 种属于木脂素和 4 种类黄酮。ELE包含 11 种单体化合物包括一种木脂素;三种环烯醚类化合物、三种酚酸类和四种类黄酮。EBE和ELE含8种相同的单体化合物，即杜仲苷、京单苷酸、曲霉酸、绿原酸、己醇C、槲皮素3-O-β-D-三丁基二糖苷、异槲皮素和山奈酚-3-O-（6′'-乙酰基）-β-D-吡喃葡萄糖苷。
图 3.粪便细菌悬浮液减轻了 HFD 诱导的肥胖。（A）C57BL/6 小鼠的代表性图像;（B）食物摄入量;（C）小鼠体重在 8 周内的变化;（D）体重和体长的交互密度图;（E）肝脏的代表性图像;（F）用苏木精和伊红染色的肝组织切片（放大倍数：×200，比例尺：100 μm）;（G）用油红O染色的肝组织切片（放大倍数：×200，比例尺：100 μm）;（H）肝脏重量;（I）油红 O 染色的脂滴结果;（J）附睾脂肪的代表性图像;（K）用苏木精和伊红染色的附睾脂肪组织切片（放大倍数：×200，比例尺：100 μm）;（l）附睾脂肪重量;（M）附睾脂肪细胞区。
在为期 8 周的 FMT 实验中，小鼠全程喂食 HFD，每组小鼠的每日喂食量表明小鼠每天摄入 2 至 3 克食物（图 3C）。所有小鼠的体重和身长逐渐增加（图 3B 和 3D）。然而，HFD-（NC）FMT 组和 HFD-（HFD）FMT 组小鼠的体重增加速度比其他组的小鼠快，第 8 周后，这些组小鼠的平均体重分别为 36.33 g 和 38.50 g，对应于与初始体重（第 0 周）相比分别增加了 10.28 g 和 11.87 g（图 3B）第8周后，HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT和HFD-（NC+ELE）FMT组小鼠体重显著低于HFD-（NC）FMT和HFD-（HFD）FMT组小鼠，HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT和HFD-（HFD+ELE）FMT组小鼠体重显著降低。HFD-（HFD）FMT 组肝组织平均重量为 1.07 g（图 3E 和 3H）。
HFD-（HFD）FMT 组肝组织平均重量为 1.07 g（图 3E 和 3H）。在健康小鼠中含有粪便细菌悬浮液的干预措施将高脂肪小鼠的肝组织重量降低到 1.01 g，但与 HFD-（HFD）FMT 组相比，这种减少并不显着（图 3、E 和 3H）。与 HFD-（HFD）FMT 组相比，含有 CGA 、 EBE 和 ELE 粪细菌悬浮液的干预措施显着减轻了高脂肪小鼠的肝组织重量。肝组织切片染色显示，HFD-（NC）FMT 和 HFD-（HFD）FMT 组分别在高脂小鼠肝组织中显示 15.90% 和 17.10% 的脂滴，这与它们的肝脏重量一致（图 3E-I）。含有 CGA 、 EBE 和 ELE 的粪便细菌悬浮液显着减少了脂滴的积累，从而减轻了肝组织的重量并预防了肝脂肪变性 （图 3E-I）。
附睾脂肪的结果表明，HFD-（NC）FMT 和 HFD-（HFD）FMT 组其重量显著增加，而 HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组的重量显著下降（图 3J 和 3L）。H&E 染色显示 HFD-（HFD）FMT 组高脂小鼠附睾脂肪细胞面积高达 9545.29 μm2含有 CGA 、 EBE 和 ELE 的粪便细菌悬浮液处理显着降低了附睾脂肪的细胞面积 （图 3K 和 3M）。此外，含有 CGA 、 EBE 和 ELE 的健康小鼠粪便细菌悬浮液的改善效果最强。
图 4.粪便细菌悬浮液对高脂血症小鼠血糖水平和脂质代谢的影响。（A）第4周空腹血糖水平;（B）第4周空腹血糖密度图;（C）第 4 周的OGTT;（D）第4周AUC值和30 min血糖值的交互密度图;（E）第8周空腹血糖水平;（F）第8周空腹血糖密度图;（G）第8周的OGTT;（H）第8周AUC值和30 min血糖值的交互密度图;（I）血浆 TC;（J）血浆 TG;（K）血浆 LDL-c;（l）血浆 HDL-c;（M）肝脏中的 TC;（N）肝脏中的 TG;（O）肝脏中的 ALT;（P）肝脏中的AST。数据是 SEM ±平均值（n = 5），根据单因素方差分析，不同的字母表示不同组之间的显着差异，然后是 Tukey 检验 （p < 0.05）。
通过在第 4 周和第 8 周测量高脂肪小鼠的空腹血糖和 OGT，发现在 HFD 下长期胃给药含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液显着影响了高脂肪小鼠血糖水平的变化（图 4A-H）。与 HFD-（HFD）FMT 组相比，HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组小鼠的空腹血糖值和曲线下面积 （AUC） 值显著降低 （图 4A-H）。第 4 周时，HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT 和 HFD-（HFD+ELE）FMT 组小鼠的空腹血糖值和 AUC 值与 HFD-（HFD）FMT 组相比显著降低，但只有空腹血糖值在第 8 周显著降低（图 4A-H）。然而，HFD-（NC）FMT 和 HFD-（HFD）FMT 组之间的空腹血糖和 AUC 值均无显著差异（图 4AH）。此外，高脂肪小鼠瘦素和 ADPN 的血清水平 显示，含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液抑制高脂肪小鼠瘦素和 ADPN 的合成，防止胰岛素抵抗，并降低血糖水平（图 S6C-D）。
通过测量血脂代谢和肝脏脂质代谢水平，我们发现给予含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液可显著改善 HFD 诱导的脂质代谢紊乱（图 4I-P 和图 S6A-B）。HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组血清 TC、TG 和 ALT 水平显著低于 HFD-（NC）FMT 和 HFD-（HFD）FMT 组，而 HDL-c 水平显著升高（图 4I-L 和图 S6A-B）。HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组肝脏中 TC、TG、AST 和 ALT 水平显著低于 HFD-（NC）FMT 和 HFD-（HFD）FMT 组（图 4M-P）。HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT 和 HFD-（HFD+ELE）FMT 组的血清 TG 水平显著低于 HFD-（NC）FMT 和 HFD-（HFD）FMT 组，而 HDL-c 水平显著升高（图 4I-L）。HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT 和 HFD-（HFD+ELE）FMT 组高脂小鼠肝脏中的 TC、TG、AST 和 ALT 水平显著低于 HFD-（NC）FMT 和 HFD-（HFD）FMT 组（图 4M-P）。身体的 TG 水平由 TBA 水平精确反映，TBA 水平与脂质代谢直接相关（图 S7）。这进一步表明，施用含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液可恢复高脂肪小鼠的正常脂质代谢水平，并且这些粪便细菌悬浮液在健康小鼠中比在高脂肪小鼠中更有效。
图 5.粪便细菌悬浮液对高脂小鼠结肠组织结构和 SCFA 产生的影响。（A）结肠的代表性图像;（B）用苏木精和伊红染色的结肠组织切片（放大倍数：×200，比例尺：100 μm，黑色箭头：中性粒细胞浸润，绿色箭头：细胞排列不紧密，红色箭头：粘膜下层破坏）;（C）结肠长度;（D）乙酸的浓度;（E）丙酸的浓度;（F）正丁酸的浓度;（G）异丁酸的浓度;（H）正戊酸的浓度;（I）i-戊酸的浓度;（J）总SCFA;（K）丁酸与血浆中TC和TG的相关性;（l）丁酸与肝脏中TC和 TG的相关性。数据是SEM ±平均值（n = 5），根据单因素方差分析后跟Tukey检验（p < 0.05），不同的字母代表不同组之间的显着差异。
粪便细菌悬浮液对 HFD 诱导的肠粘膜损伤的保护作用。结果表明，含有 CGA 、 EBE 或 ELE 的健康小鼠粪便细菌悬浮液可以改善小鼠的结肠组织，其结构相对完整，没有破坏粘膜下层（图 5A-C）。此外，HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组高脂小鼠结肠组织中闭合小带-1 （ZO-1） 和密蛋白-1 基因的表达水平与 HFD-（HFD）FMT 组小鼠相比显著上调，而干扰素-γ （INF-γ） 和白细胞介素 （IL）-1β 基因的表达水平显著下调（图 6A-C 和图 S8）。在高脂肪小鼠中观察到含有 CGA 、 EBE 和 ELE 的粪便细菌悬浮液的类似结果，所有这些都逆转了 HFD 诱导的结肠组织结构损伤（图 5A-C）。与 HFD-（HFD）FMT 组相比，在健康小鼠中含有粪便细菌悬浮液的干预措施显着上调了 Claudin-1 和 IL-6 的表达，但显着下调了 INF-γ 的表达（图 6A-C 和图 S8）。在健康小鼠中含有粪便细菌悬浮液的干预措施也未能显着生长结肠组织并保持结肠组织的结构完整性（图 5A-C）。因此，在提高 HFD 诱导的肠粘膜屏障损伤效果方面，粪便细菌悬浮液在健康小鼠中表现出最弱的效果，健康含有 CGA 、 EBE 和 ELE 的粪便细菌悬浮液的效果优于高脂小鼠。
HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT组和HFD-（NC+ELE）FMT组小鼠粪便丁酸和总酸含量较HFD-（HFD）FMT组小鼠显著升高，HFD-（NC+ELE）FMT组小鼠的乙酸含量较HFD-（HFD）FMT组显著升高（图5D-J）。结果还表明，含有 CGA 、 EBE 或 ELE 的粪细菌悬浮液促进了健康小鼠中 SCFA 的产生，这反过来又可能调节产生 SCFA 的肠道微生物群的相对丰度。然而，健康小鼠的粪便细菌悬液和含有 CGA 、 EBE 和 ELE 的高脂肪小鼠的粪便细菌悬液未观察到这种影响，只有 HFD-（HFD+CGA）FMT 组的总酸含量显著高于 HFD-（HFD）FMT 组（图 5DJ）。在血清和肝组织中正丁酸含量与 TC 和 TG 水平的相关性分析中，正丁酸含量与血清和肝脏组织中的 TC 和 TG 水平呈显著负相关 （p < 0.001） （图 5K-L）。因此，粪便细菌悬浮液除了调节血清和肝脏中的脂质代谢外，还可以通过调节肠道菌群促进丁酸的产生。
图 6.粪便细菌悬浮液对（A）结肠ZO-1表达的影响;（B）结肠 Claudin-1;（C）结肠 Occludin;（D）肝脏 LPL;（E）结肠Fiaf;（F）结肠GPR41;（G）结肠GPR43;（H）附睾脂肪GPR41;（I）mRNA 水平的附睾脂肪GPR43。
研究发现，与 HFD-（HFD）FMT 组相比，含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液可显著上调结肠组织中空腹诱导的脂肪因子 （Fiaf） 基因表达，显著下调肝组织中脂蛋白脂肪酶 （LPL） 基因表达，并通过 Fiaf-LPL 肠-肝轴循环改善 HFD 诱导的脂质代谢紊乱（图 6D-E）。SCFA 含量的增加上调了结肠组织中 G 蛋白偶联受体 （GPR） 基因 GPR41 和 GPR43 的表达（图 6FG）。当在不同组中测量高脂肪小鼠结肠组织中 GPR41 和 GPR43 的表达水平时，发现含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液显着上调 GPR43 的表达，而不是 GPR41 的表达 （图 6FG）。通过测定 HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT、HFD-（NC+ELE）FMT、HFD-（HFD+CGA）FMT、HFD-（HFD+EBE）FMT 和 HFD-（HFD+ELE）FMT 组附睾脂肪组织中 GPR41 和 GPR43 的表达观察到类似的结果（图 6H-I）。与 HFD-（HFD）FMT 组相比，FMT 组高脂小鼠附睾脂肪组织中 GPR43 的表达显著上调，但 GPR41 的表达未显著上调（图 6H-I）。因此，含有 CGA 、 EBE 或 ELE 的粪便细菌悬浮液可能会通过 SCFAs-GPR43 肠脂肪轴抑制脂质在脂肪细胞中的积累，从而抑制脂肪合成。粪便细菌悬浮液给药建立的 Fiaf-LPL 肠肝轴和 SCFAs-GPR43 肠脂肪轴取决于肠道微生物群结构的变化。因此，我们通过 FMT 将肠道微生物群和 CGA 、 EBE 和 ELE 的次生代谢物组合成一批新的小鼠，以基于 16S rRNA 测序探索微生物作用的潜在靶点。
图7.FMT下家庭水平实验组肠道菌群组成。（A）科水平物种组成分析;（B）科级分类等级树形图;（C）科水平物种组成热图。
基于粪菌悬液改善 HFD 诱导的脂质代谢紊乱的分析结果，含有 CGA 、 EBE 和 ELE 的粪菌悬液对健康小鼠的干预效果显著优于含有 CGA 、 EBE 和 ELE 的粪菌悬液对高脂小鼠的干预效果。结果，选择 HFD-（NC）FMT、HFD-（HFD）FMT、HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组的粪便进行 16S rRNA 测序。
粪便样品的稀疏曲线和物种积累曲线表现出先增加后趋于平缓的趋势。序列总数为 30,109,857，平均序列为 423，这表明本实验获得的数据具有很高的置信度（图 S9A-C）。ASV-Venn 图和 α 多样性分析结果显示，HFD-（HFD）FMT 组高脂小鼠的物种多样性和群落丰富度显著降低，而通过胃冲洗提供给高脂小鼠的含有 CGA 、 EBE 或 ELE 的健康小鼠粪便细菌悬浮液显着增加了肠道菌群的多样性和群落丰富度（图 S9D-E）。基于 Bray-Curtis 距离算法的非度量多维缩放 （NMDS） 和主成分分析 （PCA） 结果显示，HFD-（HFD）FMT 组小鼠与 HFD-（NC）FMT、HFD-（NC+CGA） FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组小鼠的距离显著，其中 PC1 和 PC2 解释了 72.6% 和 12.9% 的变异， 分别占总变化的 85.5%（图 S9F、图 S9H）。此外，我们基于 Bray-Curtis 距离算法和相似性分析方法分析了组间差异，结果显示 HFD-（NC）FMT、HFD-（HFD）FMT、HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组与 HFD-（HFD）FMT 组显著不同（图 S9G）。这与 NMDS 分析中观察到的结果一致。因此，用于 FMT 实验的高脂小鼠分组良好，并基于这些组进行微生物群分析。
分类组成、物种组成热图、分层聚类和门系统发育 krona 物种组成分析揭示了基于门、科和属水平的差异（图 7-8 和图 S10-11）。厚壁菌门、拟杆菌门和变形菌门是最重要的三个门。厚壁菌门/拟杆菌门值是肥胖抑制的重要指标，HFD-（NC）FMT、HFD-（HFD）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组显著低于 HFD-（HFD）FMT 组。变形菌门的增加通常与炎症反应有关，但通过灌胃提供给高脂肪小鼠的含有 CGA、EBE 或 ELE 的健康小鼠粪便细菌悬浮液显着降低了高脂肪小鼠粪便中变形菌门的相对丰度，而健康小鼠粪便细菌悬浮液的干预显着增加了变形菌门的相对丰度。
图 8.分析门、科和属水平不同物种的相对丰度。
在科水平上，HFD-（NC）FMT、HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组高脂小鼠粪便中丹毒丝菌科的相对丰度与 HFD-（HFD）FMT 组小鼠相比显著降低。HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT和HFD-（NC+ELE）FMT组与HFD-（HFD）FMT组相比，毛螺菌科和瘤胃球菌科的相对丰度显著增加，但脱硫弧菌科的相对丰度无显著变化。尽管每天提供相同体积的粪便细菌悬浮液，但粪便细菌悬浮液的差异可导致某些特定菌株丰度的显着变化。例如，与 HFD-（HFD）FMT 组相比，HFD-（NC+CGA）FMT 组毛螺菌科的相对丰度显著降低，Porphyromonadaceae的相对丰度显著增加。在属水平上，与 HFD-（HFD）FMT 组相比， HFD-（NC+CGA）FMT、HFD-（NC+EBE）FMT 和 HFD-（NC+ELE）FMT 组别的相对丰度显著降低，颤螺菌属体的相对丰度显著增加。然而，Desulfovibrio 的相对丰度没有显着变化。此外，含有 EBE 或 ELE 的健康小鼠的粪便细菌悬浮液显着增加了高脂肪小鼠粪便中 Erysipelotrichaceae_Clostridium 和瘤胃球菌的相对丰度;含有 ELE 的健康小鼠的粪便细菌悬浮液也显着增加了高脂肪小鼠粪便中 Lachnospiraceae_Clostridium 的相对丰度。
图 9.正丁酸和脂质代谢水平与微生物群的相关性分析。（A） 正丁酸与肠道菌群的相关性分析;（B） 科水平系统发育树图;（C） 不同物种 ASV 值与脂质代谢的相关性分析。
根据高脂小鼠粪便中正丁酸的含量以及毛螺菌科、瘤胃球菌科、颤螺菌属和普拉梭菌的相对丰度，正丁酸水平与瘤胃球菌科的相对丰度呈正相关，但与 Lachonospiraceae 的相对丰度没有显着相关性。 在属水平上，颤螺菌属和普拉梭菌的相对丰度与正丁酸含量不显著相关（图 9A）。
基于科水平的进化树分析，利用 Pearson 的相关性来检查各种物种的 ASV 与脂质代谢指标之间的关系（图 9B-C）。研究结果表明，Erysipelotrichaceae 、 Allobaculum 、 Lactobacillaceae 和 Lactobacillus 的变化与脂质合成呈正相关。瘤胃球菌科、 Oscillospira、瘤胃球菌属 和 Pophyromonadaceae 的变化与脂质合成呈负相关 （图 9C）。根据 Pearson 的相关性，我们进一步确定了 Erysipelotrichaceae、Allobaculum、Ruminococcaceae 和 Oscillospira 受 FMT 的显著影响（图 9C）。通过下调丹毒丝菌科和异杆菌属的相对丰度，上调瘤胃球菌科和颤螺菌属的相对丰度，显著抑制高脂小鼠血液和肝脏中脂质合成，从而改善长期 HFD 诱导的脂质代谢紊乱。
Erysipelotrichaceae 和 Ruminococcaceae 的 ASV 与 Fiaf 、 LPL 、 GPR41 和 GPR43 的表达水平相关 （图 10A）。丹毒科物种的相对丰度与结肠组织中Fiaf、GPR41、GPR43的表达水平、附睾脂肪组织中GPR43的表达水平呈显著负相关，与肝组织中LPL的表达水平呈显著正相关。瘤胃球菌科与结肠组织中 Fiaf 、 GPR41 和 GPR43 的表达水平以及附睾脂肪组织中 GPR43 的表达水平呈显著正相关，与肝组织中 LPL 表达水平呈显著负相关。因此，对含有 CGA 、 EBE 或 ELE 的健康小鼠的粪便细菌悬浮液获得的结果证明，对丹毒丝科和瘤胃球菌科相对丰度的调节可以成为通过 Fiaf-LPL 肠-肝轴和 SCFAs-GPR43 肠-脂肪轴改善脂质代谢紊乱的潜在靶点（图 10B）。
文章结论：
含有 EBE 或 ELE 的健康小鼠粪便混悬液可通过调节丹毒科和瘤胃球菌科的相对丰度来建立 Fiaf-LPL 肠-肝轴和 SCFAs-GPR43 肠-脂肪轴循环，从而发挥抗高脂血症作用。还发现肠道菌群变化与高脂血症之间的因果关系为治疗高脂血症提供了新的作用途径和靶点。还证明了通过杜仲饮食疗法预防高脂血症的理论基础。此外，我们发现 ELE 具有与 EBE 相似的抗高脂血症作用。因此，杜仲叶可以替代中药中的杜仲皮，从而实现植物资源的可持续利用。
文章来源：柏冬实验室
咨询电话：辉哥18911968507（微信同号）
