Haplotype and genetic diversity analysis of Spodoptera frugiperda invading in three provinces of central and southern China
-
摘要:目的
揭示入侵我国中南三省(区)的草地贪夜蛾Spodoptera frugiperda种群间可能存在的遗传变异，探讨其入侵来源、入侵路线和扩散方式，为制定有效控制和阻击草地贪夜蛾的方案提供理论依据和遗传信息。
方法利用分子标记技术从线粒体DNA和核DNA角度研究草地贪夜蛾种群的遗传多样性和种群结构。对采集自我国广东、广西和湖南3个省(区)的草地贪夜蛾，利用CO I和Tpi基因片段进行了单倍体型分析，利用9个微卫星标记分析种群间的遗传结构和遗传多样性。
结果供试样品中，除了湖南张家界群体有2个样品为CO I玉米型单倍体型、1个样品为杂合型，其他群体的个体均为CO I水稻型单倍体型。基于SSR标记的群体遗传分析结果表明，被检测的群体遗传多样性较巴西等国家的美洲群体低，各群体遗传分化不明显。广东东莞群体与湖南张家界群体的遗传距离最近，与广西南宁群体的遗传距离较远。
结论遗传距离与实际距离之间不一定呈正相关，通过遗传结构的分析有助于推测草地贪夜蛾扩散过程中，天气背景场、地理隔离和人为传播各自所起的作用。入侵湖南省张家界的草地贪夜蛾出现了低频率CO I玉米型和高频率的CO I水稻型单倍体型，说明虫源的非单一性；草地贪夜蛾在入侵地定殖时间较短，目前未形成明显的遗传分化。
Abstract:ObjectiveTo reveal possible genetic variation among Spodoptera frugiperda populations in three provinces of central and southern China, trace the invasion source and spread pathway of S. frugiperd, and provide the theoretical basis and genetic information for effectively controlling this insect.
MethodThe samples of S. frugiperda were collected from Guangdong, Guangxi and Hunan of China. Haplotype analysis was carried out using CO I and Tpigene markers. Nine microsatellite (SSR) markers were used to analyze the genetic diversity and genetic structure of S. frugiperda populations.
ResultsExpect that two samples of S. frugiperda populations from Zhangjiajie of Hunan were characterized to be CO I corn strain and one sample to be hybrid strain, other samples were CO I rice train. The results of SSR markers showed that lower genetic diversity was detected for each population compared with those in South America populations, and the genetic differentiation across populations was no significant. However, the genetic distance was the nearest between Guangdong Dongguan population and Hunan Zhangjiajie population, and the farthest between Guangdong Dongguan population and Guangxi Nanning population.
ConclusionThere is no significant relationship between genetic distance and geographic distance. The genetic structure analysis helps to speculate on the role of weather background field, geographical isolation and human transmission in the process of S. frugiperda spread. S. frugiperda invading in central-south China might be original from multiple locations. The short time of S. frugiperda colonization in the invasive lands leads to the low genetic diversity. No obvious genetic differentiation has formed at present.
-
Keywords:
- Spodoptera frugiperda /
- haplotype /
- SSR marker /
- genetic diversity
-
种群是物种存在、进化及种间关系的基本单位，是连接群落和个体的纽带，也是生态系统的基本组成部分[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 草地贪夜蛾群体间Nei’ s标准遗传距离的UPGMA法聚类图
Figure 1. UPGMA cluster analysis using the Nei’ s standard genetic distance indices for Spodoptera frugiperda
![]()
图 2 基于SSR标记的草地贪夜蛾群体的主坐标分析
1：湖南张家界群体为主；2：广东三地群体为主；3：广西南宁群体为主
Figure 2. Principal coordinate analysi based on SSR markers for Spodoptera frugiperda populations
1：The populations were mainly from Zhangjiajie of Hunan；2：The populations were mainly from three cities (district) of Guangdong；3：The populations were mainly from Nanning of Guangxi
表 1 供试幼虫样本信息表
Table 1 Information of larva samples used in this study
样品编号
Sample code采集地点
Sampling site采集时间
Sampling date虫态
Stage寄主植物
Host数量/头
NumberGDHD 广东省广州市花都花山镇儒林村 2019−05−10 幼虫 玉米 15 HNZJJ 湖南省张家界市慈利县通津钱铺 2019−05−20 幼虫 玉米 11 GXNN 广西壮族自治区南宁市宾阳县古济 2019−05−21 幼虫 玉米 7 GDDG 广东省东莞市 2019−05−01 幼虫 玉米 3 GDQY 广东省清远市连山壮族瑶族自治县 2019−06−16 幼虫 玉米 1 GDSG 广东省韶关市翁源县 2019−06−19 幼虫 甘蔗 5 合计 Total 42 
下载: 导出CSV
表 2 草地贪夜蛾2种单倍型CO I基因序列多态性位点比较
Table 2 Polymorphic site distribution of CO I gene sequence of haplotypes in Spodoptera frugiperda populations
单倍体型1)
Haplotype多态性位点 Polymorphic site 33 69 120 351 426 432 462 492 HM136592(CS) T T C T T C C T HNZJJ-30# T T C T T C C T HM136602(RS) C A T C C T T C GDHD-7# C A T C C T T C 1)HM136592(CS)：草地贪夜蛾玉米型；HNZJJ-30#：湖南张家界30#样本；HM136602(RS)：草地贪夜蛾水稻型；GDHD-7#：广州花都7#样本
1) HM136592(CS)：S. frugiperda corn strain；HNZJJ-30#：Sample 30# from Zhangjiajie of Hunan；HM136602(RS)：S. frugiperda rice strain；GDHD-7#：Sample 7# from Huadu of Guangdong
下载: 导出CSV
表 3 CO I和Tpi标记鉴定的草地贪夜蛾42个样品单倍型样品数量
Table 3 The number of haplotype of 42 Spodoptera frugiperda samples identified by COI and Tpi markers
采样地点
Sampling siteCO I 标记 CO I marker Tpi 标记 Tpi marker 水稻型
Rice strain玉米型
Corn strain杂合型 (S. f/M. s)
Hybrid strain玉米型 (AT)
Rice strain水稻型 (GA)
Corn strain杂合型 (AT/GA)
Hybrid strain广东花都 GDHD 15 0 0 9 3 3 湖南张家界 HNZJJ 8 2 1 6 0 5 广西南宁 GXNN 7 0 0 6 0 1 广东东莞 GDDG 3 0 0 3 0 0 广东清远 GDQY 1 0 0 1 0 0 广东韶关 GDSG 5 0 0 2 0 3 合计 Total 39 2 1 27 3 12
下载: 导出CSV
表 4 9个SSR位点的遗传变异参数、F-statistics统计值和基因流1)
Table 4 Indices of gene diversity，F-statistics values and gene flow deduced by nine SSR loci
SSR位点 Site NA NE I HO HE FIS FIT FST NM spf05 3.000 1.718 0.62 0.221 0.338 0.346 0.465 0.183 1.117 spf1502 3.500 1.832 0.647 0.323 0.334 0.035 0.271 0.244 0.773 spf1592 6.667 5.000 1.633 0.958 0.762 −0.258 −0.13 0.102 2.206 spf1706 2.167 1.283 0.337 0.217 0.186 −0.163 −0.048 0.099 2.269 spf343 4.667 3.524 1.291 0.92 0.686 −0.342 −0.188 0.115 1.929 spf670 5.500 3.887 1.376 0.824 0.681 −0.209 −0.015 0.160 1.312 spf789 5.500 4.340 1.437 0.93 0.707 −0.315 −0.111 0.156 1.357 spf918 4.167 2.737 1.125 0.824 0.608 −0.354 −0.227 0.094 2.407 spf997 4.500 2.856 1.135 0.536 0.612 0.125 0.303 0.204 0.976 平均值 Average 4.408 3.020 1.067 0.639 0.546 −0.126 0.036 0.151 1.594 1) NA: 观测等位基因数; NE: 有效等位基因数; I: Shannon信息指数; HO: 表观杂合度; HE: 预期杂合度; FIS: 种群内近交系数; FIT: 总近交系数; FST: 种群间分化系数；NM：基因流参数
1) NA : Allele number; NE : Effective allele number; I: Shannon index; HO : Observed heterozygosity; HE: Expected heterozygosity; FIS : Inbreeding coefficient within population; FIT : Total inbreeding coefficient; FST: Differentiation coefficient among populations; NM: Gene flow among populations
下载: 导出CSV
表 5 草地贪夜蛾群体遗传多样性指数1)
Table 5 The genetic diversity indexes of Spodoptera frugiperda populations
地理种群 Population NA NE I HO HE F 广东花都 GDHD 8.444±0.784 4.374±0.728 1.594±0.188 0.629±0.085 0.685±0.074 0.087±0.069 湖南张家界 HNZJJ 5.222±0.76 3.467±0.518 1.300±0.166 0.657±0.100 0.640±0.066 −0.028±0.099 广西南宁 GXNN 4.778±0.572 3.412±0.479 1.285±0.145 0.683±0.111 0.651±0.058 −0.033±0.118 广东东莞 GDDG 2.333±0.441 2.071±0.370 0.649±0.183 0.556±0.157 0.389±0.101 −0.426±0.151 广东韶关 GDSG 4.000±0.601 3.126±0.631 1.111±0.174 0.644±0.093 0.578±0.072 −0.113±0.097 平均值 Average 4.407±0.554 3.020±0.482 1.067±10.162 0.640±0.119 0.546±0.076 −0.252±0.089 1) NA: 观测等位基因数; NE: 有效等位基因数; I: Shannon信息指数; HO: 表观杂合度; HE: 预期杂合度; F: 固定指数
1) NA: Allele number; NE: Effective allele number; I: Shannon index; HO: Observed heterozygosity; HE: Expected heterozygosity; F: Fixation index
下载: 导出CSV
表 6 草地贪夜蛾不同地理种群的Nei’ s遗传距离(D)和遗传相似性(I)1)
Table 6 Genetic distance (D) and genetic identity (I) between different geo-populations of Spodoptera frugiperda
种群 Population 花都 Huadu 张家界 Zhangjiajie 南宁 Nanning 东莞 Dongguan 韶关 Shaoguan 花都 Huadu 0.895 0.863 0.831 0.866 张家界 Zhangjiajie 0.111 0.798 0.781 0.814 南宁 Nanning 0.147 0.225 0.670 0.771 东莞 Dongguan 0.185 0.247 0.400 0.740 韶关 Shaoguan 0.144 0.206 0.260 0.301 1)表中对角线以上为遗传相似性(I)，对角线以下为Nei’ s标准遗传距离(D)
1)The genetic identity data(I) are located above the diagonal line，and Nei’ s standard genetic distance(D) below the diagonal line
下载: 导出CSV
-
[1] TODD E L, POOLE R W. Keys and illustrations for the armyworm moths of the Noctuid genus Spodoptera Guenée from the Western Hemisphere[J]. Ann Entomol Soc Am, 1980, 73(6): 722-738. doi: 10.1093/aesa/73.6.722
[2] 联合国粮食及农业组织. 秋粘虫[Z/OL]. [2019-07-08]. http://www.fao.org/fall-armyworm/zh/. [3] UZAYISENGA B, WAWRU B, KAJUGA J, et al. First record of the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae), in Rwanda[J]. Afr Entomol, 2018, 26(1): 244-246. doi: 10.4001/003.026.0244
[4] 魏启文. 草地贪夜蛾入侵中国21个省份, 如何快速监测识别仍是防控痛点[Z/OL]. [2019-07-08]. http://finance.ifeng.com/c/7oNarkr1YVk. [5] HAJIBABAEI M, JANZEN D H, BURNS J M, et al. DNA barcodes distinguish species of tropical Lepidoptera[J]. P Natl Acac Sci USA, 2006, 103(4): 968-971. doi: 10.1073/pnas.0510466103
[6] NAGOSHI R N, GOERGEN G, TOUNOU K A, et al. Analysis of strain distribution, migratory potential, and invasion history of fall armyworm populations in northern Sub-Saharan Africa[J]. Sci Rep, 2018, 8(1): 3710. doi: 10.1038/s41598-018-21954-1
[7] 张磊, 靳明辉, 张丹丹, 等. 入侵云南草地贪夜蛾的分子鉴定[J]. 植物保护, 2019, 45(2): 19-24. [8] 陈冬平, 陈志廷, 徐汉虹, 等. 不同分子标记对入侵广州草地贪夜蛾的寄主型鉴别能力比较[J]. 华南农业大学学报, 2020, 41(1): 17-21. [9] 徐丽娜, 胡本进, 苏贤岩, 等. 入侵安徽省草地贪夜蛾的遗传分析[J/OL]. 植物保护, [2019-07-28]. https://doi.org/10.16688/j.zwbh.2019348. [10] ELLEGREN H. Microsatellites: Simple sequences with complex evolution[J]. Nat Genet, 2004, 5(6): 435-445. doi: 10.1038/nrg1348
[11] ARIAS R S, BLANCO C A, PORTILLA M, et al. First microsatellites from Spodoptera frugiperda (Lepidoptera: Noctuidae) and their potential use for population genetics[J]. Ann Entomol Soc Am, 2011, 104(3): 576-587. doi: 10.1603/AN10135
[12] PAVINATO V A C, MARTINELLI S, DE LIMA P F, et al. Methodology microsatellite markers for genetic studies of the fall armyworm, Spodoptera frugiperda[J]. Gent Mol Res, 2013, 12(1): 370-380. doi: 10.4238/2013.February.8.1
[13] ARIAS O C E C. Population genetic structure and demographic history of Spodoptera frugiperda (Lepidoptera: Noctuidae): Implications for insect resistance management programs[J]. Pest Manag Sci, 2019, 75(11): 2948-2957. doi: 10.1002/ps.5407
[14] NAGYLAKI T. Fixation indices in subdivided populations[J]. Genetics, 1998, 148(3): 1325-1332.
[15] 师沛琼, 杨茂发, 吕召云, 等. 贵州省烟青虫遗传多样性[J]. 中国农业科学, 2014, 47(9): 1836-1846. [16] 吕召云, 杨茂发, 师沛琼, 等. 贵州省烟蚜遗传多样性分析[J]. 中国农业科学, 2013, 46(13): 2685-2694. [17] SLATKIN M, BARTON N H. A comparison of three indirect methods for estimating average levels of gene flow[J]. Evolution, 1989, 43(7): 1349-1368. doi: 10.1111/j.1558-5646.1989.tb02587.x
[18] 张磊, 柳贝, 姜玉英, 等. 中国不同地区草地贪夜蛾种群生物型分子特征分析[J/OL]. 植物保护, [2019-07-08]. https://doi.org/10.16688/j.zwbh.2019296. -
期刊类型引用(3)
1. 杨洲，区湘颖，张志宏，罗宇康，何俊宇，欧俊辰，段洁利. 早期香蕉枯萎病Foc4双探针核酸纸基检测传感器研制. 农业工程学报. 2024(05): 38-46 . 
百度学术
2. 杨亚娥，卢加举，王忠泽，刘凡值，罗立娜. 香蕉主要病虫害及防治措施研究进展. 现代农业科技. 2024(17): 86-93 . 百度学术
3. 梁基校，陈东，利齐欣，林鸿生，张景欣，杨祁云，王飞钊，蒲小明. 香蕉种苗中枯萎病菌和细菌性软腐病菌实时荧光PCR检测. 农业研究与应用. 2023(05): 8-14 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1701
- HTML全文浏览量: 15
- PDF下载量: 1574
- 被引次数: 3
