紫花风铃木(Handroanthus impetiginosus)是紫葳科(Bignoniaceae)风铃木属落叶乔木，原产于南美和中美地区，可作为用材树种和观赏树种。盛花期时观赏效果极佳，花量大，花色多样^[1–3]。自20世纪70年代引入我国，紫花风铃木常被用作行道树、庭院树、园林观花树种且表现良好^[4–5]。但目前, 我国紫花风铃木园林应用具有如下问题：一是引种来源较为复杂，前人研究已发现我国园林栽培种中存在紫花风铃木和红花风铃木的杂交种^[6]以及种质资源和命名混乱的情况^[7]。二是紫花风铃木良种选育工作比较滞后，园林应用中可观测到不同个体在观赏效果上的差异，有的个体花量小、花色不显著, 群植时花期不统一。因此亟待开展紫花风铃木的遗传多样性研究，有助于推动繁育和保存紫花风铃木优质种质资源工作的开展，同时为优质种质创新提供支撑，提升景观观赏效果。

目前关于紫花风铃木的研究主要集中在栽种、观赏性、抗逆性等研究^[8–14]，关于分子水平遗传多样性在风铃木类中开展较多，单以紫花风铃木作为研究对象开展研究较少。对风铃木类植物的DNA分子水平遗传多样性的探究，已采用多种分子标记技术对风铃木类植物进行分类与鉴别^[6–7]和遗传多样性差异分析，如单核苷酸多态性(single nucleotide polymorphism, SNP)位点结合测序基因分型技术(genotyping-by-sequencing, GBS)、叶绿体trnL-trnF、gapC序列分子遗传分析、简单重复序列(simple sequence repeats, SSR)分子标记^[15]和区间简单重复序列(inter simple sequence repeats, ISSR)分子标记^[16]等。

在分子标记技术中，SSR分子标记技术因具有高灵活性、成本低、易推广、多态性丰富、重复性好的优势而得到广泛应用^[18–19]。SSR又叫微卫星DNA，SSR的产生是在DNA复制或修复过程中DNA滑动和错配或者有丝分裂、减数分裂姐妹染色单体不均等交换的结果，具有较高的物种特异性^[20–21]。SSR分子标记技术目前已广泛应用于植物的种群遗传学、遗传图谱构建、指纹分析等研究领域^[19–21]。当前国内还没有针对紫花风铃木进行基于SSR标记的群体遗传多样性研究, 因此本文利用SSR分子标记技术对广东省内6市12个群体的紫花风铃木进行遗传多样性分析研究, 探究其遗传多样性丰富度、遗传结构以及类群的划分，为紫花风铃木优良种质资源的保存和优质繁育提供理论依据。

紫花风铃木(Handroanthus impetiginosus)材料采自广东省6市12个群体，共72份，来源于同一个采样地的样本为1个群体(表 1)。于2020—2021年采集样本，采集新鲜叶片后，放入-80 ℃冰箱保存。

表 1(Table 1)

表 1 紫花风铃木群体信息 Table 1 Information of sample collection from 12 populations of Handroanthus impetiginosus 代码 Code 形态特征 Morphological character 数量 Number 地点 Location sf 叶尖为渐尖或长渐尖，叶基楔形，树皮为灰褐色，树皮为深纵裂或纵裂，花色为深粉红色 5 佛山市顺德区 Shunde Dist., Foshan fs 叶形为椭圆或长椭圆，花色为深粉红色或紫红色 4 佛山市三水区Sanshui Dist., Foshan zc 叶尖为渐尖，树皮光滑，花色为深粉红色或白色 4 广州市增城区Zengcheng Dist., Guangzhou ns 叶形为椭圆形，叶基楔形，花色以深粉红色为主 5 广州市南沙区Nansha Dist., Guangzhou th 叶缘细锯齿或锯齿，叶基楔形为主，树皮颜色灰白色为主，花色以深粉红色和粉红色为主 5 广州市天河区天河公园 Tianhe Park, Tianhe Dist., Guangzhou ld 叶尖尾尖，叶缘具锯齿，叶基楔形，树皮颜色为灰褐色，树皮以纵裂为主，花色深粉红色 8 广州市天河区龙洞 Longdong, Tianhe Dist., Guangzhou hn 叶尖为渐尖，叶缘具锯齿。叶基楔形，树皮光滑，花色粉红色或紫红色 3 广州市天河区华南农业大学 SCAU, Tianhe Dist., Guangzhou jl 叶尖尾尖或渐尖为主，叶缘具锯齿，叶形椭圆，树皮灰棕色，花色粉红色和深粉红色为主 10 深圳市龙岗区金龙路 Jinlong Road, Longgang Dist., Shenzhen xm 叶缘具锯齿，叶形为卵圆或椭圆，树皮灰褐色或灰棕色，花色以紫红色为主 8 深圳市福田区香蜜湖路 Xiangmihu Road, Futian Dist., Shenzhen jm 叶形以椭圆和长椭圆为主，叶基楔形为主，树皮灰棕色和灰白色为主，树皮光滑为主，花色以深粉红色为主 8 江门市 Jiangmen hd 叶缘锯齿，叶形为长椭圆，树皮颜色为灰白色，树皮纵裂，花色粉红色 2 惠州市惠东县Huidong County, Huizhou zs 叶缘具锯齿，叶基以楔形为主，树皮以灰褐色为主，树皮以深纵裂为主，花色为粉红色和深粉红色 10 中山市翠沙路 Cuishau Road, Zhongshan 总计Total 72

表 1 紫花风铃木群体信息 Table 1 Information of sample collection from 12 populations of Handroanthus impetiginosus

风铃木叶片中含有较多多糖和多酚物质^[22]，为减少DNA提取时多糖和多酚的干扰，参考周仙莉等^[23]的改良CTAB法提取DNA, 然后用1.5%琼脂糖凝胶电泳和核酸检测仪检测DNA质量和浓度, 后保存于4 ℃冰箱中待用。

选取同属植物H. billbergii的30对引物(mTb002~mTb052)^[24]和紫葳科印叶藤属植物Stizophyllum riparium的28对引物(Stiz2~Stiz45)^[25]，共计58对引物进行聚合酶链反应(polymerase chain reaction, PCR)扩增。所用引物由生工生物工程(上海)股份有限公司提供。随机选取8个样本进行引物的特异性和多态性筛选，最终筛选出9对杂带较少、主带清晰、多态性较好、条带差异明显的SSR引物进行后续实验(表 2)。

表 2(Table 2)

表 2 紫花风铃木的SSR引物 Table 2 SSR primers of Handroanthus impetiginosus 引物Primer 正向引物Forward primer (5′~3′) 反向引物Reverse primer (5′~3′) 荧光标记Fluorescent labelling mTb002 F: GGAACGTGCTAGTGTGTGTG R: AGAGAGTGAGTTGCAACAAAAGT 羧基荧光素Carboxy fluorescein mTb013 F: CTTCTCATTCATTTTGGTGG R: GCTTCAACACTTTCACACATC 羧基荧光素Carboxy fluorescein mTb015 F: AGCAACACAAGGAGCATTT R: GTCAGACCCAATAACTTACCTTC 羧基荧光素Carboxy fluorescein mTb018 F: GTGGTGCAGCGACTTCT R: ACATCATCGTCATCCTCATC 羧基荧光素Carboxy fluorescein mTb020 F: CATTGACTCGTTGTCCC R: GATCCTACAGTCTCACATAGAAG 羧基荧光素Carboxy fluorescein mTb024 F: CCATTTGCTTGCCTTACCCA R: AAGCAAACAACCACTCTGCA 六氯荧光素Hexachloro fluorescein mTb032 F: CGTCGAATATCTAGTGTGGG R: ACAGATGAAGAGAAAACCAAAG 六氯荧光素Hexachloro fluorescein mTb041 F: CGACATTCTTGCTCCCAATCA R: AAACAGCGGCAAGAAAGGTT 羧基荧光素Carboxy fluorescein mTb043 F: GGTCTAGCACGTGACTAACC R: CCCAATACGAGGCATATGTGA 羧基荧光素Carboxy fluorescein

表 2 紫花风铃木的SSR引物 Table 2 SSR primers of Handroanthus impetiginosus

采用荧光标记法进行基因分型，在每对引物的5′端标记荧光染料。PCR反应体系为30 μL，包含2×SanTaq PCR Master Mix 15 μL，正、反向引物(10 μmol/L)各1 μL，模板DNA (50 ng/μL) 1 μL, ddH₂O 12 μL。PCR扩增程序为95 ℃预变性5 min，然后95 ℃变性30 s，57 ℃复性90 s，72 ℃延伸30 s, 共28个循环，最后72 ℃延伸30 min，扩增产物于4 ℃保存备用。

利用Genemaker v2.2.0软件对基因分型结果进行判读，按引物预期扩增长度范围读出符合条件的数值并进行整理。利用GenAlEx 6.5进行遗传多样性统计、分子变异方差检验分析(analysis of molecular variance, AMOVA)和主坐标分析(principal coordinates analysis, PCoA)。利用Structure 2.3.4软件对群体进行遗传结构分析，并用Structure Harvester计算△K和LnP(D)，从而计算出最佳群体数目(K值)。利用PowerMaker 3.25计算多态性信息量(polymorphism information content, PIC)和计算Nei’s遗传距离和遗传一致度，并使用MEGA X 10.1进行邻接聚类分析(neighboring-joining clustering analysis, NJ聚类)。

从表 3可见，9对引物在72份材料中检测出123个等位基因，观测等位基因数(N_a)和有效等位基因数(N_e)分别为6.000~22.000和2.317~12.314，均值分别为13.667和6.072。Shannon指数(I)为1.133~2.519，平均为1.964，遗传多样性较为丰富。平均观测杂合度(H_o)和期望杂合度(H_e)分别为0.310和0.775。基因流(N_m)为0.519~1.611，平均为0.941。

表 3(Table 3)

表 3 9对SSR引物的遗传多样性参数 Table 3 Genetic diversity parameters of 9 pairs of SSR primer 引物 Primer 观测等位基因数 Number of observed allele gene (Na) 有效等位基因数 Number of effective allele gene (Ne) Shannon指数 Shannon index (I) 观测杂合度 Observed heterozygosity (Ho) 期望杂合度 Expected heterozygosity (He) 基因流 Gene flow (Nm) 多态信息量 Polymorphic information content (PIC) mTb002 8.000 3.624 1.511 0.014 0.724 0.519 0.682 mTb013 9.000 3.851 1.675 0.403 0.740 0.741 0.713 mTb015 18.000 3.708 1.966 0.347 0.730 1.611 0.717 mTb018 11.000 2.817 1.515 0.194 0.645 0.719 0.620 mTb020 15.000 7.058 2.199 0.236 0.858 1.078 0.843 mTb024 15.000 9.341 2.414 0.681 0.893 1.241 0.884 mTb032 19.000 9.618 2.519 0.222 0.896 0.902 0.888 mTb041 6.000 2.317 1.133 0.264 0.568 0.855 0.529 mTb043 22.000 12.314 2.743 0.431 0.919 0.802 0.913 平均Mean 13.667 6.072 1.964 0.310 0.775 0.941 0.754

表 3 9对SSR引物的遗传多样性参数 Table 3 Genetic diversity parameters of 9 pairs of SSR primer

从表 4可见，群体间的平均N_e为3.272个, 其中jm群体的最大，hd群体最小；I为0.448~1.630，平均为1.159，最大、最小的群体分别为jm和hd，与N_e的一致；H_e为0.292~0.734，H_o为0.222~0.389。

表 4(Table 4)

表 4 紫花风铃木的群体遗传多样性 Table 4 Genetic diversity of Handroanthus impetiginosus populations 群体 Population 数量 Number Na Ne I Ho He sf 5 4.222 3.209 1.246 0.222 0.658 fs 4 3.222 2.865 1.012 0.306 0.583 zc 4 3.333 2.857 1.057 0.222 0.604 ns 5 4.778 3.821 1.334 0.356 0.669 th 5 3.889 3.067 1.143 0.378 0.611 ld 8 4.444 3.178 1.213 0.319 0.632 hn 3 2.889 2.422 0.851 0.370 0.481 jl 10 6.333 4.426 1.525 0.356 0.718 xm 8 5.444 3.900 1.391 0.222 0.672 jm 8 6.778 4.924 1.630 0.389 0.734 hd 2 1.778 1.659 0.448 0.278 0.292 zs 10 3.889 2.941 1.062 0.278 0.559 平均Mean 6 4.250 3.272 1.159 0.308 0.601 sf、fs、zc、ns、th、ld、hn、jl、xm、jm、hd、zs见表 1; Na、Ne、I、Ho、He见表 3; 下同 sf, fs, zc, ns, th, ld, hn, jl, xm, jm, hd, zs see Table 1. Na, Ne, I, Ho, He see Table 3. The same below

表 4 紫花风铃木的群体遗传多样性 Table 4 Genetic diversity of Handroanthus impetiginosus populations

AMOVA分析表明(表 5)，紫花风铃木群体间的遗传变异仅占8%，而大部分遗传变异发生在群体内(92%)，说明紫花风铃木群体间结构弱，群体间遗传分化程度相对较低，群体内遗传分化程度相对较高。从表6可见，紫花风铃木群体间的分化系数(F_st)为0.077，居于0.050~0.150，表明所有样本的遗传分化有7.7%来自群体间，有92.3%来自群体内部，其结果与反映的结果一致。来自群体内近交系数(F_is)是0.574，群体间近交系数(F_it)为0.607。

表 5(Table 5)

表 5 紫花风铃木群体的AMOVA分析 Table 5 AMOVA analysis of Handroanthus impetiginosus populations 来源 Source 自由度 df 均方和 SS 均方MS 方差分量 Variance component % 群体间 Among populations 11 92.289 8.390 0.274 8 群体内 Within populations 132 490.850 6.552 3.276 92 总计Total 143 502.139 14.942 3.550 100

表 5 紫花风铃木群体的AMOVA分析 Table 5 AMOVA analysis of Handroanthus impetiginosus populations

用Structure软件中基于贝叶斯聚类方法对紫花风铃木群体的遗传结构进行初步分析(图 1)，当K=2时最佳。此时可将12个紫花风铃木群体分为2个不同类群(图 2)。其中类群1包括sf、ld、zc、zs共4个群体，彼此之间遗传背景较为相似，尤其是ld和zs结合前文简要收集的形态学信息也观察到这2个群体表型性状较接近；类群2包括ns、th、hn、jl、xm、jm、hd和fs共8个群体。

图 1(Fig. 1)

图 1 基于△K=mean(|L*(K)|)/sd(L(K))准则选择K值 Fig. 1 Choose K value based on criterion of △K=mean(|L*(K)|)/sd(L(K))

图 2(Fig. 2)

图 2 紫花风铃木群体的遗传结构图 Fig. 2 Population genetic structure maps of Handroanthus impetiginosus

主坐标分析大致将紫花风铃木群体分为2个类群，类群A有sf、ld、zs、zc群体，类群B有ns、th、hn、hd、fs群体，而jm、jl、xm群体的个体分散在A、B类群中，但大多数个体分布在类群B中(图 3)。主坐标分析中个体位置越近说明亲缘关系较近，紫花风铃木不同群体的个体间存在相互交叉, 说明来自不同群体的个体间可能存在较近的亲缘关系。主坐标分析结果与Structure结构分析的结果基本一致，在Structure结构分析中jm、jl、xm则与ns、th、hn、hd、fs群体归属一类。但结合形态学特征来看，类群与形态特征之间并没有明显的关系，不同形态特征的群体被分为同一类。

图 3(Fig. 3)

图 3 紫花风铃木群体的主坐标分析。黄色和蓝色圆圈分别为类群A和类群B。 Fig. 3 PCoA of Handroanthus impetiginosus populations. Yellow and blue circles are group A and group B, respectively.

应用NJ聚类法对12个紫花风铃木群体基于Nei’s遗传距离进行群体聚类分析，可划分为2个类群(图 4)，第一类包括sf、ld、zs、zc、xm、jl、jm群体，来自佛山、广州、中山、深圳和江门；第二类包括ns、th、hn、hd、fs，来自广州、惠州和佛山。这与Structure结构分析和主坐标分析的结果相近, 均可被分为2个群体。xm、jl、jm这3个群体在不同的分析方法中归在不同的类群里，说明紫花风铃木遗传距离与地理距离之间可能无显著相关性。

图 4(Fig. 4)

图 4 基于Nei’s遗传距离的紫花风铃木群体的NJ聚类分析 Fig. 4 NJ cluster analysis of Handroanthus impetiginosus populations based on Nei's genetic distance

生物的遗传多样性是物种多样性和生态系统多样性的基础，也是生物多样性的基础。一般情况下，狭义遗传多样性侧重于描述物种之间的遗传变异。丰富的遗传多样性是物种适应环境变化的物质基础, 物种为了适应环境变化而进化, 遗传变异与进化相伴相生。变异的丰富度越大意味着物种的遗传多样性越高, 表明物种适应能力越强。评价物种遗传多样性水平，可以从形态学或DNA分子等角度切入研究^[26–27]，而随着分子标记技术的发展, 其在物种鉴定、种群遗传多样性研究及物种起源等方面得到了广泛的应用^[28]。相比于使用形态学评价遗传多样性，DNA分子水平分析受环境影响较小，稳定性高，结果更为可靠^[29]。而群体遗传多样性主要分析和把握群体内和群体间的遗传多样性分布情况，可以用来衡量种群的数量变化、遗传变异、生态习性等情况^{[26–27, 30–31]}。

对紫花风铃木遗传多样性开展DNA分子水平研究，能为科学利用和培育优良紫花风铃木品种提供更多的参考。目前已经有多种分子标记技术在风铃木类植物中应用于遗传多样性研究，叶绿体trnL-trnF和gapC序列标记和SNP分子标记用于对风铃木类树种进行分子水平的鉴种和分类，纠正采集样本中风铃木类树种同名异物或异物同名的问题，同时也表明紫花风铃木遗传多样性较为丰富，种下变异较多^[6–7]。Collevatti等^[15]对分布在巴西Altamiro de Moura Pacheco生态公园的4种风铃木植物进行SSR分子标记研究，采集了75个紫花风铃木的样本, 平均H_e和H_o分别为0.857和0.703，结合种群结构分析表明紫花风铃木具有较高的遗传多样性, 但多态性水平较低。Alves等^[16]采用8对ISSR引物对巴西Açu国家森林公园内30个紫花风铃木样本进行研究，共检测出62个等位基因位点，PIC为0.34~0.49，平均为0.43，表明ISSR标记的多态性为中等，能有效反映紫花风铃木的遗传多样性。同时，当Nei’s遗传距离为0.70时可将30份样本划分成12个群体，但Shannon指数仅为0.52，说明多样性指数较低。

一般认为PIC > 0.5时该标记是高度多态的；PIC为0.25~0.5时具有中度多态；PIC < 0.25时为低多态性的^[32]。本研究中9对引物的PIC为0.529~0.913，说明具有高度多态性，比Collevatti等^[15]基于SSR分子标记和Alves等^[16]基于ISSR分子标记的紫花风铃木样本遗传多样性分析的多态性较高，说明筛选的9对引物可以用于对紫花风铃木群体的遗传多样性研究。

本文采用9对SSR引物对来自12个群体的72份紫花风铃木样本进行遗传多样性分析，Shannon多样性指数为1.964，H_e和H_o分别为0.775和0.310，表明紫花风铃木具有丰富的遗传多样性，采集自深圳和江门的样本群体遗传多样性相对较高，但均比Collevatti等^[15]研究中的H_e和H_o偏小。H_e明显大于H_o，说明供试紫花风铃木样本受到人为选择或近交等因素的影响较大。李蓉蓉^[6]的研究也表明，目前园林苗圃中存在紫花风铃木和红花风铃木的杂交个体，这值得进一步去探究园林应用中的紫花风铃木的杂交现象。

物种的群体遗传结构受地理位置、生境环境等多种因素的影响。Ayala^[33]对群体的遗传分化程度进行了划分，F_st为0~0.05时群体间遗传分化很小, 可不考虑；F_st为0.05~0.15时群体间存在中等程度的遗传分化；F_st为0.15~0.25时群体间遗传分化较大；F_st大于0.25表明群体间有很大的遗传分化。本研究中12个紫花风铃木群体的F_st为0.077，表明群体的遗传分化处于中等程度。AMOVA分析表明，紫花风铃木遗传变异仅有8%来自群体间，而大部分遗传变异发生在群体内。而群体的基因流为0.941，小于1。本研究中的紫花风铃木群体处于中等水平遗传分化和低基因流，这主要与紫花风铃木的引种与栽培繁殖方式有关。目前我国园林栽培应用的紫花风铃木引种来源不一，情况较为复杂。本研究选取的紫花风铃木样本多栽种于道路或城市绿地中，每个群体的个体构成受到人为因素影响较大，并不是以自然更新为主形成的自然群体，群体遗传分化规律性较弱。人们更侧重选择观赏效果优良的个体栽种在一起，因此同一样地采集的紫花风铃木可能有多个引种来源。再加上园林应用中为了加快优良紫花风铃木繁殖速度，常常采用嫁接、扦插等无性繁殖^[10]手段，更延缓了紫花风铃木遗传分化速度。此外，由于本研究在收集样本阶段没有记录样本对应植株是否为实生苗还是无性繁殖苗，因此若样本中有来自实生苗个体组成的群体，其在苗圃繁育阶段遗传分化也可能受到紫花风铃木传粉方式的影响。紫花风铃木授粉存在同株同花合子后自交不亲和的现象，但可以进行同株异花传粉自交^[34]，紫花风铃木簇生花相和花量大的特点也利于促进同株异花传粉自交。其次，紫花风铃木的花结构决定了其主要依靠具有长距离飞行能力且具有长喙的昆虫进行传粉^[35]，因此对异株传粉也有一定的限制。

依据Structure结构分析、主坐标分析以及聚类分析，均可将12个群体大致分为2大类群，相互补充，但与本研究中简要提及的各群体形态学上的特征分类并没有完全对应，因此日后对紫花风铃木开展遗传多样性分析可以将分子研究结合更多细致的形态学指标如雌雄蕊的细部结构特征，进行更加全面的遗传多样性分析。这3种方法中，来自深圳的xm、jl以及江门的jm这3个群体在不同方法中被归属于不同的类群，在主坐标分析中还呈现散布的特征，总体而言分布规律并不显著。推测这3个群体的个体存在种间杂交的可能，值得进一步探究。这也表明紫花风铃木群体结构与地理分布并没有完全统一，很可能还是因为不同地区的引种来源情况复杂的问题。本研究的不足之处在于选取的样本数量相对较少，且集中在广东。今后对紫花风铃木开展遗传多样性研究和选育工作时，要进一步扩充取样范围和取样数量，采集不同地区栽培种或有条件获取原产地的样本，为紫花风铃木的优质种质资源收集、科学育种等提供更多有利的遗传学数据。