Fermentation optimization of Penicillium citrinum PA-33 strain by response surface method
-
摘要:目的
优化桔青霉Penicillium citrinum PA-33发酵培养基和发酵条件，以提高桔青霉PA-33的抗菌活性。
方法采用单因素试验确定桔青霉PA-33发酵所需最适基础培养基、碳氮源和无机盐，并利用响应面法设计确定最适发酵培养基配方；在发酵条件单因素试验基础上，采用三元二次通用旋转组合设计和频率分析法优化其最适发酵条件组合。
结果经优化后，最佳发酵培养基配方为马铃薯汁液219.91 g·L−1、甘露醇34.11 g·L−1、黄豆粉6.25 g·L−1；最适发酵条件为装液量50 mL、接种量3.5% (φ)、发酵温度28 ℃，摇床转速150 r·min−1、发酵时间12 d。优化后发酵液对大肠埃希菌Escherichia coli的抑菌圈直径为28.99 mm，较优化前抑菌圈直径 (18.73 mm) 增加了10.26 mm。
结论采用响应面法、三元二次通用旋转组合设计和频率分析法优化发酵工艺，显著提高了桔青霉PA-33发酵液的抗菌活性，为该菌株的抗菌活性物质的分离以及工业化生产提供依据。
Abstract:ObjectiveTo optimize the medium composition and fermentation conditions of Penicillium citrinum PA-33 to improve its antibacterial activity.
MethodSingle factor experiments were used to determine the optimum basic medium, carbon sources, nitrogen sources and inorganic salts. The optimum formulation of fermentation medium was determined by response surface methodology. On the basis of single factor tests of fermentation conditions, the optimal combination of fermentation conditions was optimized by ternary quadratic rotation unitized design and frequency analysis method.
ResultThe optimal medium composition were: Potato juice 219.91 g·L−1, mannitol 34.11 g·L−1 and soybean powder 6.25 g·L−1. The optimal fermentation conditions were: Liquid volume 50 mL, inoculation concentration 3.5% (φ), culture temperature 28 ℃, rotation speed 150 r·min−1 and fermentation for 12 d. The inhibition zone diameter of fermentation broth after optimization on Escherichia coli reached 28.99 mm, increased by 10.26 mm compared with 18.73 mm under the original conditions.
ConclusionResponse surface methodology, ternary quadratic rotation unitized design and frequency analysis method significantly enhance the antibacterial activity of fermentation broth of Penicillium citrinum PA-33, and this study provides a basis for isolation of antibacterial active substances and industrial production of this strain.
-
种群是物种存在、进化及种间关系的基本单位，是连接群落和个体的纽带，也是生态系统的基本组成部分[1]。扩散是几乎所有生物生活史的一个基本特征，生物的扩散过程大大影响生物组织的形成，如：个体的存活、生长和繁殖，种群和群落的组成、结构和动态，物种的维持、进化和地理分布等[2-3]。目前，种群扩散研究日益受到广泛关注，已成为生态学研究中的一个热点领域[4-7]。在自然界中，大部分植物通过花粉和种子扩散来完成种群的扩张和迁移，花粉和种子扩散是关乎植物群落动态的一个重要生态过程[8]，它与生境、种源特征和传播者有着密切的关系，因植物生活型和所处生境的不同而存在明显的差别[9-10]；同时，它能够直接影响幼苗的命运，进而影响种群的更新动态和分布格局[11]。因此，种群扩散机制的研究，对种群的维持、物种的进化、物种的成种和预测其种群演化趋势等方面具有重要意义。
黄牛木Cratoxylum cochinchinense是亚热带地区的乡土树种[12]，其抗逆性好、适生性强，在酸性、砂质等贫瘠土壤中生长良好，是我国华南地区的优良先锋树种[13-14]。黄牛木还是极具潜力的城市景观树种，其黄褐色树干形成的特色林内景观以及春夏花期形成的彩色景观具有较高观赏价值[15]。目前，对黄牛木的研究主要集中于成分分析及药用价值方面[16-17]，其种群生态学的研究较少，鲜见有种群扩散方面的研究报道。纵观我国华南珠三角等发达地区，尽管城市内的绿地景观已发展得较为成熟，但城市周边仍有一些生态脆弱的山林地，这些困难立地林相不稳定，景观效果有待提升。黄牛木具备优良的品种特性，可作为该地区植被修复和生态改造工程的优选树种。因此，黄牛木种群扩散动态的研究，不仅可为其种群的天然更新和保护管理提供理论依据，为华南地区乡土树种的种群扩散研究提供参考，同时也可促进这一优良乡土树种的开发和利用，从而为森林群落的生态修复及生态系统的可持续发展作出贡献。
1. 材料与方法
1.1 研究地概况
茶山公园(简称茶山)位于广州天河区的五山街(原五山镇，茶山是五山中的其中一山)，北与沙河镇长湴村(现天河区长兴街)相邻，南至广九铁路。地理位置为113°21′~113°22′E和23°9′~23°8′N，低丘地形，海拔93.6 m。本文调研样地选自鸡笼顶，样地土壤为花岗岩、变质石英砂岩、流纹岩和红色砾岩发育的赤红壤。
白云山风景名胜区(简称白云山)座落于广州市东北部，位于北纬23°09′~23°13′，东经113°16′~113°19′，为九连山脉向东南延伸的支脉，属丘陵地形，素有“南粤第一山”之称。白云山景区面积28 km2，最高峰为摩星岭，海拔382 m，属南亚热带季风海洋性气候，年均气温21.4~21.8 ℃，年降水量1 689.3~1 876.5 mm，雨季集中在4—9月份，且降水量占全年降水量的85%以上，年均湿度79%。本文调研样地选自竹园景区，样地土壤为流纹状花岗岩、砂页岩、石英岩及混合岩母质发育而成的酸性赤红壤[18]。
石花山森林公园(简称石花山)位于广东省台山市的东北部，地理位置为22°26′~22°81′N和112°81′~112°9′E，石花山虽高300余m，却融奇、秀、幽于一山。台山属亚热带海洋性季风气候，年平均气温21.8 ℃，年平均日照2 006 h，年均降雨量1 936 mm。公园内随处可见未被人工化的自然群落，本文调研样地选自西严寺东面山顶，样地土壤类型为花岗岩发育的酸性赤红壤。
1.2 研究方法
1.2.1 外业调查
在林分踏查的基础上，采用典型取样法，选择具有代表性的黄牛木天然群落地段作为样地布设区域。于2015年12月至2016年2月期间，在广州茶山公园的鸡笼顶(东南坡)、广州白云山风景名胜区的竹园景区(西北坡)分别选取3 900 m2的样地，并在台山石花山森林公园的西严寺东面选取3 200 m2的样地(其中东南坡2 000 m2，西北坡1 200 m2)，总面积11 000 m2。对样地内所有高度≥1.3 m的黄牛木个体的胸径、树高、冠幅、枝下高等生长因子进行测定，高度<1.3 m的黄牛木个体测定除胸径之外的其他生长因子，并记录黄牛木个体之间的相对坐标。
1.2.2 种群径阶划分
参照有关种群大小级结构代替年龄结构的分析方法[19]，以“空间替代时间”，将林木的胸径大小作为度量其年龄大小的指标。把树高<1.3 m和胸径<1 cm的黄牛木幼苗记为1径阶；树高≥1.3 m且胸径≥1 cm时，每隔2 cm增加一个径阶，将黄牛木依次划分为若干径阶。
1.2.3 种子传播与种群扩散
运用ArcGIS分别绘制样地黄牛木种群个体和冠幅空间分布专题图，分析黄牛木个体的聚集情况和种子的传播方向及种群的扩散方式，然后采用缓冲区分析与空间计算功能测定种群各径阶的累计冠层投影面积，分析其累计面积的增长规律，结合SPSS 17.0软件进行回归模拟分析，构建种群扩散速率相关模型，进而预测黄牛木种群扩散的动态规律。
1.2.4 扩散速率
借鉴胡梦宇[20]对藜蒴Castanopsis fissa种群扩散速率的研究方法，通过计算种群达到某一冠幅时投影增长量所需要的时间，测算黄牛木种群的扩散速率，即单位时间内种群冠层投影的增长面积。本文以胸径代替时间变化，种群扩散速率=种群冠层增长面积/胸径。
2. 结果与分析
2.1 种群扩散方式分析
绘制各样地黄牛木种群扩散分布图(图1)。从图1可以看出，黄牛木个体总是在母树某一方向上呈现明显的聚集分布并逐渐向外扩散，但这一方向不全为下坡方向，白云山的黄牛木种群甚至沿斜上坡方向扩散。可见，黄牛木扩散不只是靠种子重力传播。从图1还可以看出，3个样地的黄牛木种群均表现出不同程度的聚集分布，而且在常年风向的下风向的聚集现象更为明显，说明黄牛木种子的传播受风向影响的可能性较大。风传播是常见的非生物传播方式，许多植物的果实和种子以此方式传播[21]，它们一般具有下述特征：1)细小质轻，能悬浮在空气中被风力吹送至远处，如兰科植物稗Echinochloa crusgalli、马唐Digitaria sanguinalis等；2)果实或种子的表面常生有絮毛、果翅，或其他有助于承受风力飞翔的特殊构造, 如马兜铃Aristolochia debilis、云杉属Picea、滇油杉Keteleeria evelyniana等。黄牛木种子为蒴果，一侧具翅，成熟期为夏季，成熟后小而轻[22]。初步推断，风力是黄牛木种群扩散传播的主要动力。
![]()
图 1 黄牛木种群扩散分布情况图中每一景观栅格为10 m×10 m，着色区域代表黄牛木个体的空间分布情况，红色区域为其种群斑块内的胸径最大个体，即种群母树Figure 1. Diffusion and distribution of Cratoxylum cochinchinense populationEach lattice is 10 m × 10 m in the graph，the coloring part indicates spatial distribution of Cratoxylum cochinchinense, and the red part indicates the parent tree with the largest DBH in population文献显示，广州的风向常年为偏北或东北方向，夏季盛行风向为东风和东南风[23]。从图1a可以看出，在夏季盛行风向(偏东和东南方向)和常年风向(偏北或东北方向)上，茶山样地的植株呈现明显的聚集现象和扩散趋势，在偏东方向上尤为明显，因此认为种子的扩散受风力的影响极大。当然，由于该样地夏季盛行风向也是该坡的斜下坡方向，故不能排除种子由于受重力作用沿斜坡滚落散布在母树周边区域，使得幼苗在母株周围呈现密集的簇状分布。但从图1b可以看出，白云山种群个体扩散方向为该坡的斜上坡方向，正是该地种子成熟期盛行风向(东南风)的下风向。因此可基本断定，黄牛木种群靠风力扩散。
台山的风向常年为东北风、偏北风，夏季盛行南风[24]。从图1c可以看出，石花山样地的植株在南坡沿正南方向扩散，即种子成熟期盛行风向；在北坡沿东北偏北方向扩散，为该地常年风向。故进一步确定，黄牛木种群靠风力扩散。
2.2 种群扩散动态分析
上述分析表明，种群在不同的地域均呈现出聚集分布，这种分布特性直接关系到林木冠层空间分布的形成。为更直观地呈现黄牛木种群在不同发育阶段的扩散动态，结合样地中黄牛木个体的相对位置和各径阶平均冠幅，运用地理信息系统(GIS)数字化和输出专题图功能绘制黄牛木种群的冠层扩散分布图(图2)。
![]()
图 2 黄牛木种群冠层扩散情况图中每一景观栅格为10 m×10 m，着色区域代表黄牛木冠层投影的空间分布Figure 2. Canopy diffusion of Cratoxylum cochinchinense populationEach lattice is 10 m × 10 m in the graph，the coloring part indicates projection of Cratoxylum cochinchinense canopy结合茶山样地的黄牛木种群冠层扩散分布图(图2a)，运用GIS软件中的合并计算功能，将各径阶分布依次合并，通过图层属性表的几何计算功能得到茶山黄牛木各个径阶冠层投影面积，茶山黄牛木种群第1~11径阶冠层投影增长面积依次为69.245、16.871、182.252、262.519、232.929、258.751、136.663、102.866、64.192、125.410和49.288 m2(表1)。由此可见，黄牛木种群冠层投影面积的增长量随径阶的增加而先增大后减小，说明种群扩散程度随径阶增大先上升后下降，在第4径阶时冠层面积增长量达到最大值262.519 m2，此时种群的扩散程度最高。
表 1 黄牛木种群生长因子统计与扩散速率测定结果Table 1. Growth factor statistics and diffusion rate of Cratoxylum cochinchinense population调查地点
Sample
plot径阶1)
Diameter
class平均胸径/cm
Average
DBH平均冠幅/m
Average crown width投影面积/m2
Canopy projection area累计投影面积/m2
Cumulative projection area投影增长面积/m2
Growth of projected area扩散速率2)/
(m2·cm-1)
Diffusion rate茶山
Chashan1 0.7 0.64 57.970 1 500.986 69.245 54.73 2 2.0 1.54 278.237 1 431.741 16.871 9.24 3 3.8 2.25 452.280 1 414.870 182.252 90.15 4 5.8 3.15 506.477 1 232.618 262.519 132.91 5 7.8 3.54 386.610 970.099 232.929 128.33 6 9.6 4.23 317.200 737.170 258.751 122.28 7 11.8 4.36 161.358 478.419 136.663 82.22 8 13.4 5.99 146.789 341.756 102.866 42.27 9 15.9 6.38 64.192 238.890 64.192 35.66 10 17.7 6.58 125.410 174.698 125.410 67.79 11 19.5 7.90 49.288 49.288 49.288 − 白云山
Baiyun mountain1 0.8 0.79 74.644 1 374.181 17.406 15.27 2 1.9 1.24 301.900 1 356.775 102.391 54.52 3 3.8 1.94 458.132 1 254.384 232.904 112.39 4 5.9 2.61 517.476 1 021.480 366.714 189.88 5 7.8 3.28 340.751 654.766 262.731 137.16 6 9.7 3.84 189.505 392.035 170.537 98.44 7 11.5 3.99 148.960 221.498 136.265 49.55 8 14.2 4.55 16.249 85.233 16.249 8.78 9 16.1 4.69 68.984 68.984 68.984 − 石花山
Shihua mountain1 0.8 0.38 0.775 709.426 0.115 0.09 2 2.1 1.11 16.533 709.311 6.207 3.14 3 4.1 2.10 91.076 703.104 42.808 24.96 4 5.8 2.66 177.765 660.296 100.895 45.04 5 8.0 3.49 244.906 559.401 152.533 91.28 6 9.7 4.23 220.326 406.868 186.753 79.85 7 12.1 5.25 194.700 220.115 194.700 76.35 8 14.6 1.45 1.673 25.415 1.672 1.11 9 16.1 5.50 23.743 23.743 23.743 − 1)第1径阶为树高 <1.3 m 和胸径 <1 cm 的黄牛木个体，其余径阶为树高 ≥1.3 m 且胸径≥1 cm 的黄牛木个体，胸径每隔 2 cm 增加一个径阶；2) “−” 代表需测得下一径阶数据方可进行计算
1) The first diameter class includes Cratoxylum cochinchinense trees with height <1.3 m and DBH <1 cm, the other diameter classes include trees with height≥1.3 m and DBH≥1 cm, and the diameter class increases one class with every 2 cm increase in DBH; 2) “−” indicates the value can not be calculated without data of the next diameter class结合白云山样地的黄牛木种群冠层扩散分布图(图2b)，计算得到白云山黄牛木种群第1~9径阶冠层投影增长面积依次为17.406、102.39、232.904、366.714、262.731、170.537、136.265、16.249和68.984 m2(表1)。可见，黄牛木种群的冠层投影面积的增长量随径阶的增加而先增大后减小，和茶山的冠层投影面积增长趋势相似，在第4径阶种群扩散程度最高。
结合石花山样地的黄牛木种群冠层扩散分布图(图2c)，计算得到黄牛木种群从第1~9径阶冠层投影增长面积依次为0.115、6.207、42.808、100.895、152.533、186.753、194.700、1.672和23.743 m2(表1)。可见，该地黄牛木种群的冠层扩散趋势与另外2个样地相似，冠层投影面积的增长量随径阶的增加而先增大后减小，但种群的扩散峰值有所延后，在第7径阶达到最大值，种群扩散程度最高。
2.3 种群扩散速率分析
黄牛木种群始终处于动态变化中，不同时期的种群斑块有不同的形状及大小。经统计相关数据，结合GIS空间分析与几何计算功能测定黄牛木种群冠层投影面积，进而计算种群的扩散速率，测定结果见表1。
运用SPSS 17.0软件进行曲线回归分析，得到黄牛木种群扩散速率与平均胸径的函数拟合曲线模型，茶山的相关参数估计值见表2(其他省略)。R2值越大、P值越小，两者的相关性越显著，据此选取黄牛木种群扩散速率的最适函数模型：茶山y1=−1.018 1x12+ 18.071x1+ 27.365(R2=0.443，P=0.129)；白云山y2=−3.086 2x22+ 43.308x2 − 8.6201( R2=0.821，P=0.000)；石花山y3=−1.5698x32 + 26.826x3 −39.059( R2=0.790，P=0.048)。其中x是各个径阶对应的平均胸径，cm；y是各个径阶对应的扩散速率，m2/cm。
表 2 茶山样地黄牛木种群扩散速率与平均胸径的函数模型方程及参数估计Table 2. Function model and parameter estimation of diffusion rate - average DBH for Cratoxylum cochinchinense population in Chashan sample plot方程
Equation参数估计值 Estimated parameter R2 F ν1 ν2 P a b1 b2 线性 Linear 0.003 0.023 1 8 0.882 80.041 −0.393 对数 Logarithm 0.037 0.310 1 8 0.593 61.672 8.042 倒数 Reciprocal 0.082 0.717 1 8 0.422 85.001 −30.106 二次 Quadratic 0.443 2.786 2 7 0.129 27.365 18.071 −1.018 复合 Composite 0.018 0.145 1 8 0.713 52.329 1.019 幂 Power 0.082 0.719 1 8 0.421 40.396 0.229 S 0.065 0.560 1 8 0.476 4.266 −0.514 增长 Growth 0.018 0.145 1 8 0.713 3.958 0.019 指数 Index 0.018 0.145 1 8 0.713 52.329 0.019 Logistic 0.018 0.145 1 8 0.713 0.019 0.982 模型构建结果表明，各样地种群扩散速率的函数模型均为一元二次方程。可以计算出，当种群的扩散速率无限接近0时，黄牛木种群冠幅的扩散受到阻碍，此时必须对黄牛木进行间伐，以促进种群的生长与更新。将y=0分别代入各样地的函数模型，得到x1=19.2 cm，x2=13.8 cm，x3=15.5 cm；即黄牛木的茶山种群处于第10径阶，白云山种群处于第8径阶，石花山种群处于第9径阶时，应对黄牛木种群进行间伐。因此，在目前的生长状态下，3个样地均应采取适当的间伐措施，以促进黄牛木天然林的更新生长和保护管理。
从模型构建结果还可以看出，各样地黄牛木种群扩散趋势相似，可以整理获得黄牛木种群扩散的一般规律。计算黄牛木种群胸径和扩散速率的平均值，并运用SPSS软件进行曲线拟合回归分析，得到黄牛木种群的扩散速率模型：y =−0.969 2x2 + 17.336x + 20.531(R2=0.943，P=0.020)。
2.4 单株黄牛木扩散速率分析
整理外业调查数据，得到单株黄牛木各径阶的胸径、冠幅、冠幅面积以及冠幅增长面积，结果见表3。通过对冠幅与胸径进行回归分析，得到二次函数的R2最大，为0.998，P为0.000，说明冠幅(y)与胸径(x)呈显著正相关，故单株黄牛木的冠幅与胸径的最佳拟合模型是一元二次方程，函数公式为y4=−0.009 6x42 + 0.481 9x4 + 0.303 9(R2=0.999，P=0.000)。当黄牛木的胸径为3.03 cm时，黄牛木的模拟冠幅为1.18 m，与实际调查数据1.67 m相差不大，拟合情况较好。
表 3 黄牛木单株生长因子调查与模拟结果Table 3. Survey and simulation results of growth factors of individual Cratoxylum cochinchinense tree径阶1)
Diameter
class数量/株
Number胸径/cm
DBH实际调查结果 Survey result 模拟分析结果 Simulation analysis result 冠幅/m
crown width冠幅面积/m2
crown area增长面积/m2
growth area冠幅/m
crown width冠幅面积/m2
crown area增长面积/m2
growth area扩散速率/
(m2·cm−1) Diffusion rate1 422 0.8 0.60 0.286 0.286 0.68 0.366 0.366 0.4 2 607 2.0 1.29 1.316 1.030 1.24 1.209 0.843 0.7 3 389 3.9 2.10 3.451 2.135 2.02 3.209 2.000 1.1 4 250 5.8 2.81 6.189 2.738 2.79 6.114 2.905 1.5 5 123 7.9 3.43 9.259 3.070 3.51 9.681 3.567 1.7 6 59 9.7 4.10 13.190 3.930 4.08 13.060 3.379 1.9 7 34 11.7 4.53 15.658 2.468 4.63 16.828 3.769 1.9 8 8 14.1 5.27 21.809 6.151 5.19 21.126 4.298 1.8 9 7 16.0 5.52 23.926 2.118 5.56 24.238 3.112 1.6 1) 第 1 径阶为树高 <1.3 m 和胸径 <1 cm 的黄牛木个体，树高≥1.3 m 且胸径≥1 cm 的黄牛木个体，其胸径每隔 2 cm 增加一个径阶
1) The first diameter class includes Cratoxylum cochinchinense trees with height <1.3 m and DBH <1 cm, the other diameter classes include trees with height≥1.3 m and DBH≥1 cm, and the diameter class increases one class with every 2 cm increase in DBH由冠幅−胸径函数得到单株黄牛木冠幅模拟数据，计算出冠幅扩散速率的数值，模拟结果见表3。运用扩散速率与胸径作回归分析，其最佳拟合模型亦是一元二次方程，函数公式为：y0=−0.013 5x02 + 0.310 6x0 + 0.111 3(R2=0.999，P=0.000)。
从表3中数据可以看出，单株黄牛木的扩散速率随径阶增大先增大后减小，说明黄牛木在生长发育的初期，营养、空间等资源较充足，种群扩散比较顺利；随着黄牛木年龄的增长，种群扩散速率趋势由上升趋势变为下降，直至无限接近于0，其冠幅扩散速率同样无限接近于0，此时黄牛木冠幅达到最大值，若不采取适当措施，黄牛木的冠幅增长将受到限制。
3. 讨论与结论
本文结合种群生态学和景观生态学的研究方法，通过GIS技术，对黄牛木种群的扩散规律进行系列分析，从种子的传播方式到种群的扩散动态，然后建立种群扩散速率的函数模型，进而模拟预测种群未来的发展趋势。本研究为植物种群扩散和空间分布研究提供了一种新的思路和方法。
植物种群扩散通常与种子扩散关系密切，而种子扩散往往需要借助动物、风、水、机械及人等媒介。在热带雨林中，满身是刺的黑尾狐棕靠刺鼠偷食其种子得以扩散传播[25]；胡梦宇[20]则通过对不同样地藜蒴种群的景观物种流研究，确定了藜蒴种子靠重力传播。黄牛木种子质轻具翅，其种子扩散方向多为当地夏季(黄牛木种子成熟期)盛行风向或常年风向的下风向，这与风传播的特征吻合[21, 26]，因此认为黄牛木种群靠风力扩散。但是，不排除其他方式传播的可能性，例如，在种子和根茎的扩散过程中遇到障碍使大量繁殖体集中在一处而其他地方很少或完全没有，动物及人为活动对聚集分布也有加强作用。今后可以采用标记黄牛木种子的方法来精确、直接地研究种子的传播，进一步确定种子在不同途径下的传播机制，丰富黄牛木的生物学理论，为种群空间分布动态和扩散规律提供更深厚的研究基础。
黄牛木的冠层扩散情况表明，冠层扩散程度随径阶增大先上升后下降，原因是群落环境水分、光照、养分等自然条件从充足到紧缺、种内和种间的竞争加大，最终导致扩散程度的变化，这种变化为物种生长扩散过程中的正常现象，该结果支持了前人的研究[20, 27-28]。但石花山的扩散程度峰值有所延后，原因可能是在生长发育初期，种群基数小，产生的种子也较少，以至于扩散程度较低，随着种群数量的增加，种群基数也逐年增加，扩散程度亦逐渐上升；也有可能是种群在扩散过程中受到环境、人为等方面的影响，具体原因有待进一步研究。研究黄牛木种群的冠层变化规律，有利于深入了解植物上层空间随时间变化的特点，为种群空间分布研究提供理论基础。
运用冠层面积增长动态推算种群扩散规律，结果表明，种群扩散速率随种群年龄增加先增大后减小，各样地种群扩散速率与胸径的函数模型为一元二次方程，这与胡梦宇[20]对珠三角地区不同样地藜蒴景观林的研究结果基本一致，说明不同种群的扩散趋势有相似之处。但各个模拟函数的系数均不相同，说明不同的生长环境其生长情况可能不同，模拟方程呈现的是特定区域种群生长扩散的规律，不能盲目套用，一概而论。赵强民[28]通过对广州白云山木荷种群的扩散速率进行研究，得到木荷种群的冠层累计投影面积与径阶的函数模型为一元一次方程，说明种群扩散速率的模型构建有不同的拟合选择，应充分根据不同研究对象进行合理选择。未来，还要进一步考虑边缘效应、林窗效应、自然稀疏的随机性以及种子传播过程中风速产生的影响，提升黄牛木种群扩散速率模拟预测的精确度，为种群的合理间伐及抚育更新提供理论依据。
通过对单株黄牛木的胸径增长和冠幅扩散规律进行研究，利用扩散速率与胸径的函数模型，可进一步确定人工经营下黄牛木种植密度以及间伐时间，由扩散速率-胸径函数模型计算得到，当y达到最大值1.90 m2/cm时，x =11.5 cm，即黄牛木的胸径为11.5 cm时，种群处于第6径阶，其冠幅扩散速率达到最大值1.90 m2/cm，此时应对种群进行适当间伐，以保证其最大效益；将y=0代入函数，计算得到x=23.4，即黄牛木胸径达到23.4 cm时，种群处于第12径阶，此时必须对黄牛木进行间伐，否则将会阻碍种群的生长。当然，树种的生长状况处于不断变化中，以上推论为理想状态下黄牛木的间伐时间，但可以肯定的是，适时进行合理的人工疏伐，有利于实现人工经营下黄牛木种群的最大景观价值和森林生产效益。
黄牛木作为华南地区的优良乡土树种，目前仍处于天然群落生长状态，对黄牛木扩散机制的研究还不够深入，对其种子传播后种子的命运、种子适应各种传播途径的机理以及扩散速率与风速之间的关系等方面的研究都有待加强。因此，只有了解了黄牛木的扩散机制，才能采取有针对性的措施，实现黄牛木天然林的合理抚育和有效管理。从人工经营管理的角度而言，分析黄牛木的生长规律和扩散动态，可加快黄牛木在景观绿化及植被修复中的推广应用，从而促进困难立地的生态平衡与发展。
-
![]()
图 1 不同培养基成分对发酵液抗菌活性的影响
B图中横坐标的数字1、2、3、4、5、6、7、8、9和10分别代表葡萄糖、蔗糖、淀粉、玉米粉、乳糖、果糖、糊精、甘油、甘露醇和麦芽糖；C图中横坐标的数字1、2、3、4、5、6、7、8、9、10和11分别代表对照、蛋白胨、酵母浸粉、牛肉膏、胰蛋白胨、黄豆粉、尿素、KNO3、NH4NO3、NH4Cl和陈皮；各图中柱子上方的不同小写字母表示差异显著(P<0.05, Duncan’s法)
Figure 1. Effects of different medium ingredients on the antibacterial activities of fermentation broth
In figure B，number 1，2，3，4，5，6，7，8，9 and 10 represent glucose, sucrose, starch, corn flour, lactose, fructose, dextrin, glycerin, mannitol and maltose, respectively; In figure B，number 1，2，3，4，5，6，7，8，9, 10 and 11 represent control, peptone, yeast extract, beef extract, tryptone, soybean powder, urea, KNO3, NH4NO3, NH4Cl and bran, respectively; In each figure, different lowercase letters on the bars indicate significant difference (P<0.05, Duncan’s test)
![]()
图 2 菌株PA-33抗菌活性的响应曲面图
A1、A2和A3为响应面3D图；B1、B2和B3为等高线图
Figure 2. Response surface plots of antibacterial activity of strain PA-33
A1, A2, and A3: 3D maps of response surfaces; B1, B2, and B3: Contour maps
![]()
图 3 不同发酵条件对菌株PA-33抑菌活性的影响
各图中, 柱子上方的不同小写字母表示差异显著(P<0.05, Duncan’s法)
Figure 3. Effect of different dermentation condition on antimicrobial activity of strain PA-33
In each figure, different lowercase letters on the bars indicate significant difference (P<0.05, Duncan’s test)
表 1 不同质量浓度的培养基成分对菌株PA-33抗菌活性的影响1)
Table 1 Effects of different concentrations of medium components on antibacterial activities of fermentation broth
ρ(马铃薯)/(g·L−1) Patato content 抑菌圈直径/mm Inhibition zone diameter ρ(甘露醇)/(g·L−1) Mannitol content 抑菌圈直径/mm Inhibition zone diameter ρ(黄豆粉)/((g·L−1) Soybean powder content 抑菌圈直径/mm Inhibition zone diameter 100 22.93±0.52ab 20 22.03±0.12c 5.0 23.83±0.69c 150 23.13±0.25ab 25 23.33±0.09b 7.5 27.93±0.19a 200 23.47±0.25a 30 26.60±0.83a 10.0 25.77±0.31b 250 22.53±0.19b 35 25.50±0.08a 12.5 25.13±0.09b 300 21.80±0.16c 40 23.43±0.87b 15.0 22.43±0.17d 1) 同列数据后的不同小写字母表示差异显著 (P<0.05，Duncan’s 法 )
1) Different lowercase letters in the same column indicate significant difference (P<0.05，Duncan’s test)
下载: 导出CSV
表 2 Box-Behnken试验方差分析
Table 2 Analysis of variance for Box-Behnken experiment
来源1) Source SS DF MS F P2) 模型 Model 76.25 9 8.47 22.24 0.000 2** A 1.07 1 1.07 2.80 0.138 3 B 15.62 1 15.62 41.01 0.000 4** C 14.36 1 14.36 37.70 0.000 5** AB 0.07 1 0.07 0.17 0.691 9 AC 9.52 1 9.52 24.98 0.001 6** BC 0.01 1 0.01 0.02 0.894 4 A2 14.56 1 14.56 38.22 0.000 5** B2 12.31 1 12.31 32.30 0.000 7** C2 5.19 1 5.19 13.61 0.007 8** 残差
Residual2.67 7 0.38 失拟
Lack of fit1.96 3 0.65 3.71 0.118 9 误差
Error0.71 4 0.18 总变异
Total variation78.92 16 1) A、B 和 C 分别表示马铃薯、甘露醇和黄豆粉；2) “**” 表示模型具有极显著影响 (P<0.01)
1) A, B and C represent potato, mannitol and soybean powder, respectively; 2) “**” indicates the model is highly significant(P<0.01)
下载: 导出CSV
表 3 三元二次通用旋转组合试验设计及结果1)
Table 3 Design and results for quadratic general rotary composite
试验号
Test No.θ/℃ 接种量(φ)/ %
Inoculation concentration装液量/mL
Liquid volume抑菌圈直径/mm
Inhibition zone diameter1 30(1) 5(1) 60(1) 28.30 2 30(1) 5(1) 40(−1) 28.37 3 30(1) 2(−1) 60(1) 27.80 4 30(1) 2(−1) 40(−1) 27.93 5 26(−1) 5(1) 60(1) 28.73 6 26(−1) 5(1) 40(−1) 28.87 7 26(−1) 2(−1) 60(1) 28.40 8 26(−1) 2(−1) 40(−1) 28.57 9 24.6(−1.682) 3.5(0) 50(0) 26.88 10 31.4(1.682) 3.5(0) 50(0) 27.00 11 28(0) 0.977(−1.682) 50(0) 28.80 12 28(0) 6(1.682) 50(0) 29.23 13 28(0) 3.5(0) 33.2(−1.682) 28.53 14 28(0) 3.5(0) 66.8(1.682) 28.87 15 28(0) 3.5(0) 50(0) 28.80 16 28(0) 3.5(0) 50(0) 28.93 17 28(0) 3.5(0) 50(0) 29.20 18 28(0) 3.5(0) 50(0) 29.13 17 28(0) 3.5(0) 50(0) 29.00 18 28(0) 3.5(0) 50(0) 28.90 1) 括号中的数值为编码值
1) Numbers in the brackets are coded values
下载: 导出CSV
表 4 三元二次通用旋转组合试验结果方差分析
Table 4 Analysis of variance for results from quadratic general rotary united design
来源 Source SS DF MS PCC F P X1 0.283 6 1 0.283 6 −0.515 4 3.617 2 0.086 4 X2 0.385 1 1 0.385 1 0.573 9 4.910 4 0.051 0 X3 0.000 3 1 0.000 3 0.018 9 0.003 6 0.953 5 X12 6.406 3 1 6.406 3 −0.943 9 81.696 5 0.000 1 X22 0.064 5 1 0.064 5 0.275 6 0.822 2 0.385 9 X32 0.028 5 1 0.028 5 −0.187 3 0.363 6 0.559 9 X1X2 0.012 0 1 0.012 0 0.122 8 0.153 2 0.703 7 X1X3 0.001 5 1 0.001 5 0.043 9 0.019 3 0.892 3 X2X3 0.001 0 1 0.001 0 0.035 9 0.012 9 0.911 8 回归 Regression 7.349 1 9 0.816 6 10.413 2 0.000 8 残差 Residual 0.784 2 10 0.078 4 失拟 Lack of fit 0.541 6 5 0.108 3 2.233 2 0.130 9 误差 Error 0.242 5 5 0.048 5 总和 Total 8.133 2 19 1) X1、X2 和 X3 分别表示温度、接种量和装液量；PCC 为偏相关系数
1) X1，X2 andX3 represent temperature，inoculation concentration and liquid volume, respectively; PCC is patial correlation coefficient
下载: 导出CSV
表 5 各变量取值频率分布
Table 5 Frequency distribution of variable value
水平
Levelθ 接种量 Inoculation concentration 装液量 Liquid volume 次数 Times 频率 Frequency 次数 Times 频率 Frequency 次数 Times 频率 Frequency −1.682 0 0.000 5 0.143 7 0.200 −1.000 10 0.286 5 0.143 7 0.200 0.000 25 0.714 5 0.143 7 0.200 1.000 0 0.000 10 0.286 7 0.200 1.682 0 0.000 10 0.286 7 0.200
下载: 导出CSV
-
[1] 罗清, 彭程, 叶波平. 青霉属真菌研究新进展[J]. 药物生物技术, 2016, 5(23): 452-456. [2] 郑婕施, 江北, 张琪, 等. 青霉属真菌活性代谢产物研究进展[J]. 工业微生物, 2017(06): 50-56. doi: 10.3969/j.issn.1001-6678.2017.06.009 [3] LU Z Y, LIN Z J, WANG W L, et al. Citrinin dimers from the halotolerant fungus Penicillium citrinum B-57[J]. J Nat Prod, 2008, 71(4): 543-546. doi: 10.1021/np0704708
[4] DU L, LI D, ZHANG G, et al. Novel carbon-bridged citrinin dimers from a volcano ash-derived fungus Penicillium citrinum and their cytotoxic and cell cycle arrest activities[J]. Tetrahedron, 2010, 66(47): 9286-9290. doi: 10.1016/j.tet.2010.09.036
[5] CLARK B R, CAPON R J, LACEY E, et al. Citrinin revisited: From monomers to dimers and beyond[J]. Org. Biomol Chem, 2006, 4(8): 1520-1528. doi: 10.1039/b600960c
[6] LIU H C, DU L, ZHU T J, et al. Two new citrinin dimers from a volcano ash-derived fungus, Penicillium citrinum HGY1-5[J]. Helv Chim Acta, 2010, 93(11): 2224-2230. doi: 10.1002/hlca.v93.11
[7] NI M, LIN W L, YANG P, et al. A novel citrinin derivative from the marine-source fungus Penicillium citrinum[J]. Acta Pharm, 2015, 50(2): 203-206.
[8] LIU Y, LI X M, MENG L H, et al. Bisthiodiketopiperazines and acorane sesquiterpenes produced by the marine-derived fungus Penicillium adametzioides AS-53 on different culture media[J]. J Nat Prod, 2015, 78(6): 1294-1299. doi: 10.1021/acs.jnatprod.5b00102
[9] 申光辉, 郑丽君, 张志清, 等. 解淀粉芽胞杆菌PC2产抑菌物质培养基及发酵条件优化[J]. 微生物学通报, 2017, 44(6): 1358-1369. [10] 王昊鹏, 吴黎明, 赵柳微, 等. 烟曲霉素发酵培养基的优化研究[J]. 食品工业科技, 2017: 1-14. [11] RUQAYYAH T I D, JAMAL P, ALAM M Z, et al. Application of response surface methodology for protein enrichment of cassava peel as animal feed by the white-rot fungus Panus tigrinus M609RQY[J]. Food Hydrocolloid, 2014, 42: 298-303. doi: 10.1016/j.foodhyd.2014.04.027
[12] MUTHUKUMAR V, RAJESH N, VENKATASAMY R, et al. Mathematical modeling for radial overcut on electrical discharge machining of incoloy 800 by response surface methodology[J]. Procedia Materials Science, 2014, 6: 1674-1682. doi: 10.1016/j.mspro.2014.07.153
[13] YANG P, FANG M, LIU Y W. Optimization of a phase adjuster in a thermo-acoustic stirling engine using response surface methodology[J]. Energy Procedia, 2014, 61: 1772-1775. doi: 10.1016/j.egypro.2014.12.209
[14] 李莉, 张赛, 何强, 等. 响应面法在试验设计与优化中的应用[J]. 实验室研究与探索, 2015, 34(8): 41-45. doi: 10.3969/j.issn.1006-7167.2015.08.011 [15] 梁新乐, 黄莹莹, 张虹, 等. 响应面法优化桔青霉F-5-5核酸酶P1发酵培养基碳氮源[J]. 核农学报, 2011, 25(1): 57-61. [16] 谢祥聪, 刘琴英, 蒋冬花, 等. 淡色生赤壳菌(Bionectriao chroleuca)Bo-1菌株产生抗菌物质的发酵条件优化[J]. 生物技术通报, 2014(1): 166-170. [17] 郭雷, 朱文成, 刘玮炜, 等. 抗菌活性海洋真菌HN4-13的鉴定及其发酵优化[J]. 微生物学通报, 2013, 40(6): 951-958. [18] 胡一峰, 曹一岚, 黄俊龙. 美伐他汀发酵菌种筛选和发酵条件优化[J]. 浙江化工, 2012, 43(2): 4-6. doi: 10.3969/j.issn.1006-4184.2012.02.002 [19] 喻晨, 张亚雄, 赵劼, 等. 响应面法优化桔青霉产核酸酶P1培养基[J]. 食品科学, 2011, 32(17): 283-286. [20] MENG L H, LIU Y, LI X M, et al. Citrifelins A and B, citrinin adducts with a tetracyclic framework from cocultures of marine-derived isolates of Penicillium citrinum and Beauveria felina[J]. J Nat Prod, 2015, 78(9): 2301-2305. doi: 10.1021/acs.jnatprod.5b00450
-
期刊类型引用(3)
1. 杨洲，区湘颖，张志宏，罗宇康，何俊宇，欧俊辰，段洁利. 早期香蕉枯萎病Foc4双探针核酸纸基检测传感器研制. 农业工程学报. 2024(05): 38-46 . 
百度学术
2. 杨亚娥，卢加举，王忠泽，刘凡值，罗立娜. 香蕉主要病虫害及防治措施研究进展. 现代农业科技. 2024(17): 86-93 . 百度学术
3. 梁基校，陈东，利齐欣，林鸿生，张景欣，杨祁云，王飞钊，蒲小明. 香蕉种苗中枯萎病菌和细菌性软腐病菌实时荧光PCR检测. 农业研究与应用. 2023(05): 8-14 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1586
- HTML全文浏览量: 1
- PDF下载量: 1277
- 被引次数: 3
