菲律宾蛤仔(Ruditapes philippinarum)隶属于软体动物门(Mollusca)，双壳纲(Bivalvia)，异齿亚纲(Heterodonta)，帘蛤目(Veneroida)，帘蛤科(Veneridae)，蛤仔属(Ruditapes)，是我国传统四大养殖贝类之一，广泛分布在我国南北沿海。据统计，2010年菲律宾蛤仔世界产量达到360多万t^[1]，我国年产量在300 万t左右，占世界菲律宾蛤仔产量的90%以上^[2]。由于其营养丰富、味道鲜美，近年来市场需求量也日益增大。然而由于人工苗种的累代养殖，菲律宾蛤仔养殖群体的种质逐渐下降，养殖产业迫切需要菲律宾蛤仔的优良品种。近年来，菲律宾蛤仔的遗传改良技术取得了一定的成绩，如选择育种和杂交育种等^{[2, 3, 4, 5, 6]}。传统育种方法对目标性状缺乏准确可靠的早期预选手段，造成育种盲目性大、效率低。而紧密连锁的DNA标记为实现水产动物目标性状的早期选择提供了有效途径。因此，研究与水产动物重要经济性状紧密连锁的分子标记，特别是对缺乏遗传连锁图谱的水产动物DNA标记，就显得十分重要。

微卫星(microsatellites)又称简单序列重复(Simple Sequence Repeat,SSR)，一般由1—6 bp的简单序列重复排列组成。由于微卫星具有多态性丰富、共显性遗传、遵循孟德尔分离定律、易于(Polymerase Chain Reaction) PCR扩增等优点，是分子遗传学研究的理想分子标记，已广泛应用于水产动物的重要经济性状相关分子标记筛选、(Quantitative Trait Locus) QTL定位等方面^{[7, 8, 9, 10, 11, 12, 13, 14, 15]}。在海洋贝类，虾夷扇贝^[16]和马氏珠母贝^[17]等已报道了生长性状与微卫星标记的相关性分析。目前，对菲律宾蛤仔的研究，主要集中在生物学、生态学和养殖技术等方面^[2]。然而，有关菲律宾蛤仔重要经济性状的分子基础研究还刚刚起步，标记与生长性状的相关分析迄今还未见报道。

本实验利用20 对微卫星标记对菲律宾蛤仔斑马蛤家系进行遗传多样性分析，并对标记位点与生长性状相关性进行分析，筛选出对斑马蛤生长有利的基因型，实现分子标记和生长性状的连锁分析。为菲律宾蛤仔优良品系的进一步选育，以及遗传图谱的构建和QTL定位等提供参考。

实验所用斑马蛤材料为2011年大连石河菲律宾蛤仔野生群体中挑选的斑马蛤为亲本，交配产生F1代群体，上选F1代群体中大的个体，采用巢式设计建立了30 个F1家系，本实验所用家系为第8 组家系，随机选取其中107 个个体，对所有个体的壳长、壳宽、壳高和体重进行测量。取部分足尖组织，固定于95%乙醇中用于DNA的提取。

DNA的提取采用苯酚-氯仿-异戊醇法抽提DNA，并用1%琼脂糖凝胶电泳进行检测。DNA保存于-20 ℃备用。PCR反应体系：反应总体积为10 μL，其中DNA模板1 μL、上下游引物(10 μmol/L) 各0.4 μL、easytaq MIX 5 μL、ddH₂O 3.2 μL。实验中所用微卫星引物的核心序列及退火温度见表 1。PCR反应程序为：94 ℃ 5 min，94 ℃ 40 s，退火温度40 s，72 ℃ 40 s，30 个循环；72 ℃ 延伸10 min，4 ℃保存。PCR扩增产物经12%聚丙烯酰胺凝胶电泳检测，银染显色。凝胶在扫描仪上成像，拍照保存后对电泳谱带进行分析。等位基因大小用100 bp marker为参照标准对电泳条带进行判读。

用PopGene3.2软件进行数据处理、遗传多样性分析,计算等位基因数(Na)、有效等位基因数(Ne)、观测杂合度(Ho)及期望杂合度(He)。根据Botstein等^[18]的方法可计算多态信息含量(PIC)。利用PIC-CALC计算软件，根据等位基因频率来计算每个位点的多态信息含量，其计算公式如下：

式中，n为某一位点上的等位基因数；p_i、p_j分别为群体中第i和第j个位点的等位基因频率，j=i+1。

用SPSS 20.0软件中的一般线性模型(GLM) 对菲律宾蛤仔斑马蛤家系的生长相关性状与20个微卫星标记相关性进行分析。其线性模型如下：

y_ij=μ+g_i+ε_ij

式中,y_ij为斑马蛤某性状第i个标记第j个个体表型值,μ为群体平均值，g_i为第i个标记的效应值，ε_ij为随机误差。

经测量得到的壳长、壳宽、壳高和体重4 个生长性状，均符合正态分布，可直接进行相关性分析。正态分布信息见表 2。

利用父母本和8 个子代进行引物筛选，经PCR 扩增、聚丙烯酰胺凝胶电泳，从32 对微卫星引物中筛选得到20 对扩增条带清晰、稳定，具有多态性的位点。检测了父母本和107 个F2子代个体在20 个微卫星位点的基因型(表 3)。SSR372和SSR373位点在菲律宾蛤仔斑马蛤家系部分个体中的扩增结果见图 1。

图 1 引物SSR372和SSR373在斑马蛤F2 家系部分个体的聚丙烯电泳图谱 Fig. 1 Electrophoresis pattern of PCR product of partial individuals in the F2 family of Ruditapes philippinarum amplified with microsatellite marker SSR372 and SSR373 图中M为100 bp DNA Marker，A、B为不同等位基因

图选项

20 个微卫星位点共检测到41 个等位基因，片段大小为109—430 bp，各位点的等位基因数2—3 个，平均等位基因数2.05 个，平均有效等位基因数为1.71 个。各位点的期望杂合度、观测杂合度、基因型分布和多态信息含量等值见表 3。20 个位点的期望杂合度平均值为0.431，观测杂合度平均值为0.504，多态信息含量平均值为0.324。在20 个位点中，有2 个位点SSR213和SSR261的观测杂合度低于期望杂合度，表现出纯合子过剩、杂合子不足现象。另外18 个位点的观测杂合度则高于期望杂合度，其中SSR11和SSR164的观测杂合度与期望杂合度的差异较大。经卡方检验，有17 个位点的基因型比例符合孟德尔分离定律，另外3 个位点SSR11，SSR164和SSR213的基因型分布显著偏离了孟德尔定律(P<0.01)(表 3)。

用SPSS 20.0软件中的一般线性模型(General Linear Model,GLM)，对微卫星位点基因型与斑马蛤主要生长性状(壳长、壳宽、壳高和体重)之间的相关性进行分析，经方差分析检验呈显著性差异的位点。在20 对微卫星位点中，共找出4 个与斑马蛤生长相关的位点(表 4)。SSR9位点与壳高呈显著关联，SSR135和SSR164与壳宽显著相关，SSR142位点与体重呈显著性相关。由于相关分析得到的微卫星位点除SSR9号引物外其他3 个位点基因型个数都小于3，故使用SPSS中compare means过程进行不同位点表型值的比较。在SSR9位点上，AA基因型的壳长，壳宽，壳高和体重4 项生长指标平均值均高于其他基因型个体，说明AA基因型对斑马蛤生长起正效应。另外，SSR135位点上的BB基因型、SSR142位点上的AA基因型和SSR164位点上的AB基因型，在4 项生长指标中的平均值均高于各自位点中的其他基因型，说明4个基因型在各自位点上对斑马蛤的生长具有正效应。

群体的遗传多样性主要表现在杂合度、等位基因数和多态信息含量三个方面。杂合度大小可以反映出遗传结构变异程度的高低，杂合度期望越高，群体的遗传结构越复杂^[19]。本实验中，观测杂合度和期望杂合度的平均值分别为0.504和0.431，说明实验家系的遗传变异和遗传多样性较丰富。从等位基因来看，有效等位基因越接近所检测到的等位基因的绝对值，说明等位基因在群体中的分布越均匀^[20]。本研究在20 个微卫星位点共检测到41 个等位基因，其中有效等位基因数为34 个，接近于实际等位基因数，表明所检测的微卫星位点的等位基因在群体中分布较均匀。多态信息含量可作为一个衡量遗传标记所包含的或能提供的遗传信息容量的指标。一般认为，在某一群体中，当PIC > 0.5时该位点表现为高度多态，0.25<PIC<0.5时该位点表现为中度多态，当PIC<0.25时该位点表现为低度多态^[21]。本实验中，19 个位点均表现为中度多态位点，1 个位点表现为高度多态位点，表明所选微卫星标记多态性相对丰富，可作为有效的遗传标记进行遗传多样性分析。

在许多海洋贝类中，经常会出现偏分离的等位基因，这似乎是海洋贝类中普遍存在的现象。在本实验中，发现3 个微卫星位点的基因型分离比显著偏离孟德尔遗传定律(P<0.01)，部分位点偏分离的原因可能与生活环境和人工选择有关，再有可能是由隐性致死基因、取样或基因型判读错误等因素引起^[22]。大部分微卫星位点的期望杂合度与观测杂合度基本接近，说明群体内的基因型频率未发生大的改变。但有2 个位点SSR11和SSR164的观测杂合度明显大于期望杂合度，在SSR164发现了一种纯合子完全缺失的现象(基因型AA为0)，说明该位点附近可能存在着与之紧密连锁的隐性致死基因^[23]。

菲律宾蛤仔的壳长、壳宽、壳高、体重等性状属于数量性状，数量性状受多个基因的控制，其遗传基础比较复杂且容易受环境的影响。标记与性状之间的连锁分析，是根据标记位点的基因型以及数量性状的表型对个体进行显著性检验，差异显著则说明标记与数量性状存在关联^[24]。因此，如果一个群体的性状差异显著，或两个群体的差异很大，就可以通过标记与性状的相关分析，找出性状与一个或多个标记的遗传相关，一旦发现显著相关，即可认为存在一个数量性状位点，从而实现从表型到基因型选择育种的转变^[25]。借助微卫星标记进行性状连锁分析或 QTL 定位，在水产动物育种中已也有较多报道^{[10, 11, 12, 13, 14, 15, 16, 17]}。

目前，菲律宾蛤仔遗传标记的开发较少,不足以满足遗传作图和QTL 定位的需要。所以，有必要应用有限的微卫星标记对菲律宾蛤仔生长相关性状进行连锁分析。在本研究中，共找到了4 个与生长性状相关的位点，SSR9位点与壳高存在显著的相关关系(P<0.05)，SSR135和SSR164与壳宽显著相关(P<0.05)，SSR142与体重显著相关(P<0.05)。在这些关联中，出现了几个标记同一个性状相关，说明这些位点存在多因一效的现象，并符合数量性状的相关理论。对上述位点不同基因型的多重比较表明，SSR9位点的AA基因型、SSR135位点的BB基因型、SSR142位点的AA基因型和SSR164位点的AB基因型都与壳长、壳宽、壳高和体重性状呈正相关。目前，分子标记在菲律宾蛤仔中的应用主要集中于遗传多样性研究^{[26, 27]}，有关菲律宾蛤仔性状相关分子标记的筛选迄今未见报道。本文筛选出菲律宾蛤仔生长性状相关的微卫星标记，为菲律宾蛤仔分子标记辅助选育提供参考。