Genetic diversity of selected GIFT strains of Oreochromis niloticus
-
摘要:
用微卫星和扩增片段长度多态性(amplified fragment length polymorphism,AFLP)分子标记对海南、广州和青岛吉富品系尼罗罗非鱼的遗传多样性和遗传分化进行了研究。9对微卫星引物和4对AFLP引物的分析结果一致。微卫星分析表明,青岛吉富鱼的平均等位基因数(4.8)、平均观测杂合度(0.528)、平均多态信息含量(0.605)最高,海南吉富鱼的平均等位基因数(4.4)、平均观测杂合度(0.479)、平均多态信息含量(0.549)最低。AFLP分析显示,青岛吉富鱼的多态位点比例(48.4%)和基因多样性(0.245)最高,海南吉富鱼的多态位点比例(36.3%)和基因多样性(0.147)最低。这些表明,青岛吉富鱼的遗传多样性最高,海南吉富鱼最低。遗传分化分析表明,两两群体间遗传分化显著(微卫星FST为0.07~0.11,P < 0.01;AFLP FST为0.24~0.29,P < 0.01)。AMOVA分析显示,大部分遗传变异(微卫星的分析结果为91.26%;AFLP的为67.6%)来源于群体内个体间,表明吉富鱼选育品系尚具有进一步选育的潜力。
-
关键词:
- 吉富罗非鱼 /
- 遗传多样性 /
- 遗传分化 /
- 微卫星DNA /
- 扩增片段长度多态性(AFLP)
Abstract:Genetic diversity and differentiation among three stocks of selected GIFT strains of Oreochromis niloticus was analyzed using microsatellite DNA and AFLP markers. Nine microsatellite loci and four primer pairs of AFLP revealed consistent results. Microsatellite data showed that Qingdao GIFT stock had the highest mean number of alleles (4.8), mean observed heterozygosity (0.528) and mean PIC (0.605), whereas Hainan GIFT the lowest. AFLP analysis demonstrated that Qingdao GIFT had the highest proportion of polymorphic loci (48.4%) and gene diversity within populations (0.245), whereas Hainan GIFT the lowest. These observations indicated that Qingdao GIFT has the highest genetic diversity, whereas Hainan GIFT the lowest. Genetic differentiation among the three GIFT strains was significant as revealed by microsatellite (FST=0.07~0.11, P < 0.01) and AFLP (FST =0.24~0.29, P < 0.01).AMOVA indicated that the majority of genetic variation (91.26% by microsatellite analysis and 67.6% by AFLP) was attributed to difference of individuals within stocks, suggesting that GIFT strains of O.niloticus has some potential for further selective breeding.
-
阳宗海是云南九大高原湖泊之一,属珠江流域南盘江水系,位于中国云南省昆明市,承担着所在流域及下游地区饮用、农业灌溉和工业用水三大功能[1]。随着流域开发强度的增加,阳宗海水体营养盐浓度逐渐升高,20世纪90年代曾多次出现水华现象。2008年砷 (As) 污染发生后,管理部门采取了一系列措施,治理及改善阳宗海水质[2]。截至2020年,阳宗海水质已从2008年的劣Ⅴ类恢复到III类。阳宗海水质虽通过治理取得了显著成效,但仍然面临富营养化风险[3]。多年来,管理部门通过增殖放流向阳宗海投放大量鱼苗来净化水质。据不完全统计,截至2022年,已累计放流以鲢 (Hypophthalmichthys molitrix)、鳙 (Aristichthys nobilis) 为主的1 000多万尾鱼苗 (阳宗海管理部门数据),并实施了严格的禁渔措施。以鱼净水是常用的生物调控方法之一,具有环保、投入低和持续时间长的优点,但前提是须明确鱼类群落结构及其资源现状[4]。安莉等[5]于2011—2012年对阳宗海鱼类资源进行了调查,然而至今未见新的报道,目前湖区渔业资源状况不明。
鱼类资源评估是河流渔业管理和资源保护的基本工作[6],目前常用的方法有数学模型评估[7-9]及水声学评估法[10-11]。数学模型评估主要利用鱼类体长、年龄等较容易获取的捕捞数据进行估算。水声学评估通过对研究区域的监测,可获得鱼类的尺寸规格、时空动态分布、行为特征及资源量等。目前声学技术已成为观测和掌握渔业资源变动的重要手段[12],能够同时在多个营养级 (大型食肉鱼、小型觅食鱼和浮游动物) 进行采样,结合网具捕捞数据,便可进行生物量评估[13-14]。
相关研究表明,鱼类存在昼夜迁移行为,在日落前开始扩散并向上移动,日出后开始聚集并下降[15]。例如,清河水库鲢、鳙昼间聚集于底层,夜间向上层迁移,且分散分布[16]。鱼类的昼夜行为变化导致资源评估结果在昼夜间出现差异。因此,本研究首次利用水声学技术分别于昼间和夜间对阳宗海鱼类资源进行评估,旨在掌握阳宗海鱼类资源现状,分析鱼类昼夜行为差异,评估阳宗海鱼类资源调查的适宜时间,以期为高原湖泊及封闭型水库渔业资源的保护及管理提供参考。
1. 材料与方法
1.1 研究区域
阳宗海 (102°58'E—103°01'E, 24°51'N—24°58'N) 位于昆明市东南,地跨宜良、呈贡和澄江三县,距阳宗镇约4 km,距昆明市区约36 km。阳宗海南北长12 km,东西平均宽2.65 km,流域面积约192 km2,平均水深约20 m,最大水深约30 m,是由地面断裂凹陷发育形成的地堑式断陷湖泊。其最大库容为6.17×109 m3,调查期间库容为5.96×109 m3。
1.2 鱼类调查
结合阳宗海形态特征,在主要入湖口、出湖口、湖湾和洞穴等不同生境处设置了10个采样点 (图1),于2022年3、5、8和10月进行了4次采样,使用网具为单层刺网 (网目0.65、1.2 、1.8 cm)、三层刺网 (网目3.5、6.0、15.0 cm) 和地笼网 (0.4 m×0.4 m×长20 m、1.2 m×1.2 m×长30 m),每个采样点放置不同规格网具各1张。于18:30—19:30下网,次日05:00—07:00收网,采集水深范围约为0~10 m,刺网采集水层为中上层,地笼为底层。渔获物现场鉴定、测量全长、体长、体质量等信息,鉴定主要依据文献 [17-20]。
1.3 声学调查
水声学设备为分裂波束鱼探仪 EY60 (Simrad,挪威),换能器中心频率为120 kHz,发射功率为200 W,−3 dB波束宽为7.0°,脉冲宽度为256 μs。通过不锈钢支架将换能器固定在探测船只右舷,距离船首1/3船体长度,入水深度0.5 m,波束垂直向下发射。通过Garmin 60CSx来实时确定探测位置,鱼探仪与笔记本电脑相连,实时显示和存储数据。探测之前采用碳化钨标准球 (Φ 23 mm) 对换能器进行校正[21]。
夜间和昼间的探测时间分别为2022年6月20日20:00—24:00和21日8:30—12:30。调查船为阳宗海管理部门执法船只,全长约10 m,调查期间湖面较平静,探测时以10 km·h−1的航速匀速行进。由于库区整体呈狭长型,且岸线整齐,采用非固定间距平行航线,在实际走航中略有偏差 (图1),总航程约38.8 km,调查覆盖率为6.47,达到了Aglen[22]提出的覆盖率在6.0 以上的标准。
1.4 数据分析
用Pinkas的相对重要性指数 (Index of relative importance, IRI) 研究鱼类群落优势种, 计算公式如下[23]:
$$ \mathrm{I}\mathrm{R}\mathrm{I}=(N+W)\times F $$ (1) 式中:N为某一种鱼类的尾数占总尾数的百分比 (%);W为某一种鱼类的生物量占总生物量的百分比 (%);F为某一种鱼类出现的站位数占总站位数的百分比 (%)。本文选取IRI≥500的鱼类为优势种。
声学原始数据 (.raw) 使用Echoview 5.4 软件处理。从声学映像中可看出调查区域鱼类主要为分散分布,因此采用回波计数法处理数据,目标提取使用单体目标检测及鱼体轨迹追踪实现[24]。通过定义表面线和底部线来排除表层1 m以上、底部0.3 m以下的所有声学信号。目标强度 (Target Strength, TS) 最小阈值设置为−58 dB,以确保大部分噪声信号被排除在外,其他区域噪声手动剔除。单体检测、鱼体追踪的参数设置参照Rudstam等[25]的研究。数据分析过程中,将150 m作为一个基本调查单元 (EDSU),单独评估鱼类密度、个体分布和行为特征等。每个 EDSU 的鱼类密度(尾·m−3) 的计算均采用回波计数 (EchoCounting) 方法进行[26],公式如下:
$$ \varphi = \frac{N}{{PV}} $$ (2) 式中:N为探测到的鱼类数目;
$ \varphi $ 为单位体积水体的鱼类数量,即鱼类体积密度;V为每一个ping探测的水体体积;P为分析数据的ping数量。体长与TS的关系使用基于基尔霍夫近似模型的鲢的目标强度经验公式计算[27]:
$$ S_{\mathrm{t}}=20 \log L-67.35 $$ (3) 式中:St为目标强度 (dB);L为体长 (cm)。
为分析不同规格鱼类昼夜行为差异,人为将鱼类回波信号分为3类:小个体 (−60~−53 dB)、中等个体 (−53~−47 dB) 和大个体 (>−47 dB),对应鱼类长度分别为小个体 (2~5 cm)、中等个体 (5~10 cm)和大个体 (>10 cm)。
采用 ArcGIS 10.1 软件进行鱼类资源空间分布建模,将鱼类密度及对应的经纬度数据导入ArcGIS,采用反距离加权插值方法进行栅格插值运算,导出每个栅格的鱼类密度,乘以相应的栅格体积,最后将所有栅格值求和得到阳宗海鱼类资源量。
通过TS、相对游泳速度、运动方向 (垂直和水平方向变化) 4个因子来分析调查期间阳宗海鱼类昼夜行为特征。TS可以表征昼夜鱼类体长差异,相对游泳速度可以体现鱼类昼间游泳速度的差异,鱼类在水平和垂直方向的运动情况可以通过鱼体与换能器之间的角度 (水平方向角度变化范围 0°~360°,垂直方向为−90°~90°) 来判定。
2. 结果
2.1 鱼类种类组成
本年度共采集到鱼类24种,32976尾,总质量约792 kg,隶属于7目14科23属 (表1)。其中泥鳅 (Misgurrnus anguillicaudatus)、黄鳝 (Monopterus albus)、鲫 (Carasius auratus)、花䱻 (Hemibarbus maculatus) 和云南光唇鱼 (Acrossocheilus yunnanensis) 为土著种,其他种类均为外来引入种。根据鱼类IRI值的大小排序,IRI≥500的鱼类为优势种。结果显示,太湖新银鱼 (Neosalanx taihuensis)、鲢、鳙、子陵吻鰕虎 (Rhinogobius giurinus)、䱗 (Hemiculter leucisculus)、间下鱵 (Hyporhamphus intermedius) 为阳宗海鱼类优势种。其中,太湖新银鱼、子陵吻鰕虎、间下鱵等小个体鱼类 (平均体长<6 cm) 在数量上占比超90%;从体质量组成来看,鲢、鳙等大个体鱼类 (平均体质量>1000 g) 占比较大。
表 1 阳宗海鱼类种类组成Table 1. Composition of fish species in Yangzonghai Lake种类
Species数量百分比
Quantity
proportion/%体质量量百分比
Body mass
ratio/%体长范围
Range of body
length/mm体质量范围
Range of body
mass/g相对重要
性指数
IRI栖息水层
Habitat太湖新银鱼 Neosalanx taihuensis 65.18 1.19 25~61 0.16~0.99 6 637 中上层 子陵吻鰕虎 Rhinogobius giurinus 15.43 0.37 20~62 0.10~3.57 1 580 底层 间下鱵 Hyporhamphus intermedius 4.85 0.84 85~170 1.72~11.28 569 上层 䱗 Hemiculter leucisculus 2.88 4.39 45~195 1.20~145.34 727 上层 黄颡鱼 Pelteobagrus fulvidraco 2.43 2.36 33~195 0.71~90.76 479 底层 黄䱂 Hypseleotris swinhonis 2.22 0.06 23~53 0.13~2.35 228 底层 高体鳑鲏 Rhodeus ocellatus 2.06 0.18 19~74 0.16~9.45 224 中下层 瓦氏黄颡鱼 Pelteobagrus vachelli 1.26 1.10 33~165 0.85~72.1 236 底层 鲢 Hypophthalmichthys molitrix 1.23 64.99 77~600 2.34~3 930.00 6 622 中上层 麦穗鱼 Pseudorasbora parva 0.72 0.14 35~94 0.60~15.22 86 底层 鲫 Carasius auratus 0.47 2.09 47~245 3.73~496.18 256 中下层 鳙 Aristichthys nobilis 0.31 17.84 77~620 10.00~4 950.00 1 815 中上层 尼罗罗非鱼 Oreochromis nilotica 0.29 1.55 25~270 0.55~680.00 147 中下层 粘皮鲻鰕虎 Mugilogobius myxodermus 0.17 0.00 20~35 0.20~0.60 5 底层 花䱻 Hemibarbus maculatus 0.14 0.40 103~220 16.61~172.67 43 中下层 大鳞副泥鳅 Paramisgurnus dabryanus 0.09 0.06 87~150 3.47~28.70 12 底层 云南光唇鱼 Acrossocheilus yunnanensis 0.08 0.09 65~141 4.34~130.00 17 中下层 泥鳅 Misgurrnus anguillicaudatus 0.05 0.02 85~130 3.61~21.40 4 底层 青鱼 Mylopharyngodon piceus 0.04 0.64 220~395 138.07~430.70 34 中下层 鲤 Cyprinus carpio 0.04 0.43 173~270 137.94~605.00 24 中下层 草鱼 Ctenopharyngodon idellus 0.04 0.70 137~490 45.90~2 660.00 22 中上层 鲇 Silurus asotus 0.04 0.54 220~500 120.00~2 000.00 12 底层 食蚊鱼 Gambusia affinis 0.01 <0.01 3.4 0.70~0.76 0.2 上层 黄鳝 Monopterus albus <0.01 <0.01 210 9.00 0.01 底层 2.2 鱼类水平分布
从声学映像图可以看出,在深水区,夜间大量弱声信号 (TS −70~−58 dB) 聚集在5~10 m水层,形成密集的声散射层 (图2-a),而昼间鱼类回波信号主要呈分散分布 (图2-b)。对探测江段鱼类资源水平分布进行GIS插值,得到探测江段的鱼类水平分布模型 (图3),昼夜鱼类水平分布主要表现为南北2个浅水区为高密度区,在湖区中心鱼类密度较低。独立样本t检验结果显示,阳宗海昼夜鱼类密度均值无显著性差异 (P>0.05)。对不同调查时间的鱼类密度进行方差齐性检验 (P<0.05),表明在显著性水平0.05下鱼类密度不服从正态分布。使用非参数检验进行分析,Kurskal-Wallis结果显示,阳宗海昼夜鱼类密度存在显著性差异 (P<0.05)。
2.3 鱼类垂直分布
整体上,鱼类主要分布在10 m以浅水层,10 m以下水层的密度极低,大部分区域密度为0。鱼类密度昼夜间分布存在显著性差异 (F=265.36, P<0.001),表现为鱼类在昼间相对于夜间整体向下游动。其中夜间鱼类平均水深为7.52 m,分布在5~10 m水层;昼间鱼类平均水深为9.99 m,分布在6~14 m水层。相关分析结果可知,夜间鱼类密度与水深呈显著性正相关 (P=0.01, r=0.22),而昼间鱼类密度与水深呈显著性负相关 (P<0.01, r=−0.43,图4)。
从个体来看,大个体鱼类在夜间分布范围较大 (四分位距为6.31 m,中位数为9.13 m),在昼间分布范围较小且向上游动 (四分位距为3.00 m,中位数为8.26 m),表明大个体鱼类在垂直方向上的迁移更为活跃。小个体和中等个体鱼类在昼间均向下游动了约2.00 m (图5)。
2.4 鱼类行为特征
如图6所示,鱼类TS、相对游泳速度、运动方向 (垂直和水平方向变化) 4个行为因子在昼夜间均存在显著性差异 (P<0.001)。昼间鱼类平均TS大于夜间,均值分别为−54.18和−55.52 dB,四分位数分别为−54.28和−55.93 dB。水平方向上,昼夜间50%鱼类游动方向均与调查船航向相反,昼间四分位距位于120°~180°,夜间位于180°~190°。垂直方向上,昼间第三分位数为负值 (−3.72),表明昼间至少有75%的鱼类向下游动;夜间第三分位数为正值 (0.02),表明夜间至少有75%的鱼类向上游动。相对游泳速度昼间大于夜间,平均值分别为2.80和2.57 m·s−1。
2.5 鱼类资源量估算
昼夜鱼类平均密度分别为 (0.033±0.063) 和 (0.038±0.152) 尾·m−3,鱼类密度统计信息见表2。图4为鱼类密度频率分布图,从图中可以看出,鱼类密度右偏,表明阳宗海鱼类低密度单元格较多。基于GIS反距离权重插值结果,对调查区域资源量进行估算。对阳宗海底图进行投影,计算得出整个调查区域面积约为37 km2,调查期间库容为5.96×109 m3,初步估算,阳宗海资源量为1.96×107~2.25×107 尾。
表 2 不同时间鱼类密度分布Table 2. Fish density distribution in different periods变量
Variable鱼类平均密度
Average fish density/(尾·m−3)95%置信区间
95% confidence interval/(尾·m−3)偏度
Skewness峰度
Kurtosis资源量
Abundance/(107 尾)夜间 Nighttime 0.0378±0.153 0.020~0.045 7.964 71.824 2.25 昼间 Daytime 0.0329±0.064 0.001~0.066 2.683 6.849 1.96 3. 讨论
3.1 阳宗海鱼类种群变化
云南高原湖泊和水库经历了从自然渔业生产到现代渔业生产的转变,鱼类区系改变巨大[9]。历史资料记载的阳宗海有鱼类18种[19-20,28],由于过度捕捞、环境污染、外来鱼种的引入,阳宗海鱼类群落结构现已发生巨大变化[4,29]。2022年不同月份的4次采样,共采集到鱼类24种,其中土著种仅5种,分别为泥鳅、黄鳝、鲫、花䱻和云南光唇鱼,均为广布种,其中花䱻和云南光唇鱼为近年增殖放流引入。从种类组成来看,土著种占20.8%;从渔获量看,土著鱼类的数量仅占总量的0.88%,除鲫有一定的资源量外,其他种类数量极少。阳宗海特有种中仅阳宗白鱼 (Anabarilius yangzongensis) 于2010年采集到1尾,其他种类如阳宗金线鲃 (Sinocyclocheilus yangzongensis)、短尾鱊 (Acheilognathus brevicaudatus)、阳宗海云南鳅 (Yunnanilus yangzonghaiensis) 和纺锤云南鳅 (Y. elakatis) 已多年未见[4]。阳宗海土著鱼类资源严重衰退,目前已形成了以外来鱼类为优势种、土著鱼类种类和数量不断减少的新群落结构,这也与云南其他高原湖泊鱼类种群的变迁基本一致[30-32]。
从鱼类生态习性角度来看,主要优势种子陵吻鰕虎为底栖鱼类,由于渔探仪存在水底探测盲区,加之在数据处理中去除了底部30 cm以上的水体,因此该生态类型鱼类未包括在声学结果分析中;另外,太湖新银鱼为弱声学散射体,且在夜间聚集分布,因此声学数据处理时需单独分析。阳宗海其他种类均适用于水声学方法进行资源评估。
3.2 鱼类空间分布的影响因素
一般认为鱼类群落的动态变化受内源性和外源性因素影响[33]。内源性因素指鱼类自身繁殖、觅食、洄游及种内、种间竞争等行为导致的种群在数量、时空上的变化;外源性因素指栖息地、水文、环境等改变引起的变化[34]。阳宗海南部有阳宗大河、七星河、摆衣河等入湖河流,且湖泊南部人口密集,人类活动造成的污染较大,北部有火电厂的冷却水排入,导致湖泊南北两侧水体营养盐浓度高于其他区域,可为鱼类提供更为丰富的饵料生物[3]。因此鱼类高密度区主要集中在其南部和北部。
阳宗海鱼类垂直分布的主要影响因子是水深和水温。水深是鱼类生活空间的反映,鱼的体型大小、生活习性等决定了鱼类对水深的需求[35]。Miranda等[36]研究表明,鱼类分布与水深呈现出非随机关系,并根据季节和种类而变化,约75%的物种分布在12 m以浅。阳宗海鱼类主要分布在水体中上层,深水区密度较低,这与云南其他高原湖库的研究结果类似[5,14]。
鱼类缺乏调节体温的能力,其体温随着环境水温的变化而改变,并对温度变化产生有限的适应,水温过高或过低均会对鱼类的生长、活动产生影响[5]。阳宗海是暖单次混合型湖泊[37],存在明显的季节性分层与混合现象。在夏季存在明显的温跃层,表层相对恒温,介于21~24 ℃,温跃层位于9~13 m,与表层温度相差高达7~8 ℃[3];随着水温出现分层,溶解氧也呈现出与水温分层相似的垂直分层结构[38]。本次声学调查时间为夏季,此时温跃层的水温和溶解氧急剧下降,因此阳宗海鱼类大多分布在温跃层以上,即10 m以内水层。
3.3 鱼类行为昼夜差异的影响
鱼类行为在昼夜间存在显著性差异,主要表现为中小型鱼类在夜间上浮、昼间下沉的特点,而夜间大个体鱼类垂直迁移更为活跃,符合鲢、鳙等鱼类的生态习性[16]。鱼类昼夜行为节律差异的成因主要是昼夜光照和温度差异,促使浮游动植物等饵料生物呈现周期性昼夜垂直迁移行为[39],间接影响到鱼类的活动。另外,小型鱼类为逃避敌害白天会躲避于水底,夜间上浮觅食[40]。但这些差异也并非绝对,对突尼斯深水型水库的研究表明,不同季节生物量在昼夜间存在差异,表现为在冬春季,白天高于夜间,夏季则夜间高于白天,在多种鱼类混栖的水体中,昼夜差异虽然存在,但规律不同[41-42]。基于阳宗海鱼类昼夜分布情况,未来阳宗海鱼类声学调查可基于不同目的选择不同的调查策略。如:以太湖新银鱼等小型鱼类为主要调查对象,可在夜间进行调查,声学数据处理方法可选用回波积分法;以鲢、鳙等大中型鱼类为调查对象,可选择在白天调查,声学数据处理方法可选用回波计数法。
夜间大量弱声信号 (TS<−58 dB) 聚集在5~10 m水层,形成密集的声散射层 (图2-a)。本研究主要针对鲢、鳙等大中型鱼类资源进行估算,因此在数据分析时,将TS最小阈值设置为−58 dB,在资源量估算中,亦未将这部分弱声信号包含在内。结合渔获调查结果及鱼类生态习性,推测是太湖新银鱼等小型鱼类。
3.4 鱼类目标强度的影响因素
种类鉴别是渔业资源声学评估中的一个难点。由于淡水鱼类分布特征为多种类混栖,且多为分散分布,通常研究人员采用声学调查结合渔获物采样来分配积分值,以评估资源量[43]。此次调查过程中虽未同步进行渔获物采样,但通过不同季节的采样,基本上掌握了阳宗海鱼类的种类组成及群落结构。除鲢、鳙外,目前缺乏阳宗海其他主要优势种类 TS 值与鱼体大小关系的系统研究,无法准确根据鱼类TS值计算出相应的鱼类大小,因此本研究选择主要优势种鲢的体长与TS的经验公式。鱼类TS值受脉冲发射频率、有无鱼鳔、鱼类在波束中的位置、鱼类在波束中的倾角和河流环境等因素影响[44-45],因此文中估算出的鱼类尺寸会与实际有所偏差。本研究中鱼类TS值主要用于探讨不同规格的鱼类在昼夜行为上的差异,鱼类资源评估主要是通过鱼类密度估算得出。因此TS经验公式的选取对于讨论鱼类资源变化的影响较小,研究结果能够较为准确地反映湖泊鱼类的昼夜行为差异和资源空间分布结构。
3.5 阳宗海渔业资源保护与利用建议
随着渔业方式的转型升级,近年来大水面生态渔业逐步推广。其主要目的是通过自然修复和人工群落调控的手段维护湖泊生态系统健康、实现营养盐移出和净水渔业功能[46]。结果显示,小型鱼类在阳宗海鱼类群落中的占比很高,仅太湖新银鱼、子陵吻鰕虎在渔获物中的数量占比就超过了80% (表1);在声学调查结果中,小于−51 dB (6 cm)的个体占总数的94.5%,这些信号大部分来源于高体鳑鲏 (Rhodeus ocellatus)、黄䱂 (Hypseleotris swinhonis)、麦穗鱼 (Pseudorasbora parva) 等小型鱼类,这些鱼类的主要摄食对象为浮游动物,导致浮游植物生物量增加[47]。小型鱼类过多会对湖泊浮游动物种群造成巨大压力,导致湖泊浮游动物种类、丰度出现不同程度的下降[1],不利于浮游植物的调控[48]。研究表明,在洱海中通过对西太公鱼 (Hypomesus nipponensis) 的定量清除来加强浮游动物的牧食,进而控制浮游植物的密度和生物量,对于防控蓝藻水华爆发具有重要意义[49]。据了解,除太湖新银鱼的专项捕捞外,阳宗海实施了严格的禁渔制度。另外,资源评估结果显示,鲢、鳙的生物量为1.83~2.11 g·m−3,该生物密度远未达到可以有效遏制水华的生物水平[48]。本研究建议:1) 在每年开展太湖新银鱼捕捞的基础上,同时对其他小型鱼类进行生态捕捞;可投放阳宗海所属水系的中上层土著肉食性鱼类,将小型鱼类转化为优质水产品;2) 持续进行鲢、鳙等鱼类的增殖放流,以达到以鱼净水的目的;3) 对阳宗海生态容量进行评估,以制定更加科学合理的湖泊保护利用管理策略。
-
表 1 筛选的罗非鱼微卫星引物的序列
Table 1 Sequences of screened SSR primers of tilapia
位点 locus 引物序列5′-3′
primer sequence 5′-3′重复单元
repeat unit注册号
GeneBank accession no.UNH735 F:GTGCACACAGCAGAGGCTAA
R:CTTTGCTGCTGTCGGAGTTTgt G63984 UNH846 F:TGGAGCAGCTTCTTCTACATCA
R:CACATGATGGAAGCCGTGTAca G68185 UNH849 F:CCCAGGTGCATATTTCCCTA
R:CTTGCTCTGCTTTGCTGAAAgt G68187 UNH855 F:ACTCCCGCTGTTGCTGTTAG
R:GAGGGGAGCCTACAACGTAAac G68191 UNH871 F:CGTGACAGCTGTGTATGCAA
R:TGGCATTCAGATTGAAATGGTac G68201 UNH954 F:GGAAAACGTTTGGAGAGACG
R:AAACGGAGCTCCTGTCTGAAtg G68250 UNH985 F:GCGTCTTGATGCAGGATACA
R:TCCCGACGAGCAACTGTTATca G68266 UNH995 F:CCAGCCCTCTGCATAAAGAC
R:GCAGCACAACCACAGTGCTAgt G68274 UNH1007 F:TTCTCTACCTGTGAACATTTGC
R:AAGGCAGTCCTGCTCTATGCca G68283 表 2 3个吉富鱼群体9个微卫星位点的等位基因频率
Table 2 Allele frequencies at nine microsatellite loci among three GIFT strains of O.niloticus
位点 locus 等位基因 allele 海南吉富鱼 Hainan GIFT 广州吉富鱼 Guangzhou GIFT 青岛吉富鱼 Qingdao GIFT UNH985 A 0.071 0.094 - B 0.095 - - C 0.179 0.234 0.114 D 0.524 0.438 0.443 E 0.131 - - F - 0.234 0.443 UNH871 A 0.571 0.516 0.671 B 0.262 0.047 0.100 C 0.155 0.375 0.200 D 0.012 0.063 0.029 UNH995 A 0.036 0.016 0.243 B 0.167 0.297 0.114 C 0.119 0.031 0.143 D 0.679 0.578 0.343 E - 0.078 0.157 UNH855 A 0.167 0.234 - B 0.786 0.766 0.600 C 0.048 - 0.400 UNH1007 A 0.048 - - B 0.036 0.016 0.143 C 0.571 0.469 0.214 D 0.262 0.516 0.600 E 0.083 - 0.043 UNH846 A - 0.094 0.114 B - 0.031 0.114 C - 0.125 0.057 D 0.441 0.234 0.229 E 0.083 0.391 0.243 F 0.476 0.125 0.243 UNH735 A 0.036 0.031 0.515 B 0.167 0.297 0.162 C 0.274 0.031 0.059 D 0.191 0.063 0.088 E 0.155 0.188 0.059 F 0.095 - 0.074 G 0.083 0.391 0.044 UNH849 A 0.071 0.094 0.031 B 0.107 0.063 0.016 C 0.238 0.094 0.172 D 0.191 0.406 0.219 E 0.143 0.297 0.375 F 0.250 0.047 0.188 UNH954 A 0.012 0.032 0.100 B 0.024 - 0.057 C 0.571 0.452 0.357 D 0.381 0.387 0.257 E - 0.129 0.129 F 0.012 - 0.100 表 3 3个吉富鱼选育群体的遗传多样性
Table 3 Genetic diversity of three selected GIFT strains of O.niloticus
位点 locus 参数 parameters 海南吉富鱼 Hainan GIFT 广州吉富鱼 Guangzhou GIFT 青岛吉富鱼 Qingdao GIFT UNH985 A 5 4 3 Ho 0.595 0.844 0.343 He 0.670 0.701 0.603 PIC 0.625 0.637 0.508 UNH871 A 4 4 4 Ho 0.571 0.656 0.486 He 0.588 0.597 0.506 PIC 0.517 0.507 0.451 UNH995 A 4 5 5 Ho 0.214 0.250 0.457 He 0.502 0.579 0.776 PIC 0.455 0.505 0.728 UNH855 A 3 2 2 Ho 0.429 0.406 0.800 He 0.357 0.365 0.487 PIC 0.316 0.295 0.365 UNH1007 A 5 3 4 Ho 0.571 0.719 0.486 He 0.601 0.522 0.580 PIC 0.541 0.397 0.521 UNH846 A 3 6 6 Ho 0.452 0.656 0.886 He 0.579 0.763 0.812 PIC 0.478 0.717 0.770 UNH735 A 7 6 7 Ho 0.381 0.594 0.500 He 0.830 0.730 0.697 PIC 0.796 0.671 0.660 UNH849 A 6 6 6 Ho 0.595 0.281 0.313 He 0.817 0.735 0.757 PIC 0.780 0.681 0.705 UNH954 A 5 4 6 Ho 0.500 0.065 0.486 He 0.534 0.639 0.778 PIC 0.433 0.555 0.734 平均average A 4.7 4.4 4.8 Ho 0.479 0.497 0.528 He 0.609 0.626 0.666 PIC 0.549 0.552 0.605 注:A. 等位基因数;Ho. 观测杂合度;He. 期望杂合度;PIC.多态信息含量
Note: A. number of alleles;Ho. obsevered heterozygosity;He. expected heterozygosity;PIC. polymorphism information content表 4 3个吉富鱼养殖群体的Hardy-Weinberg平衡的χ2检验
Table 4 Chi-square test for Hardy-Weinberg equilibrium of three selected GIFT strains of O.niloticus
位点
locus海南吉富鱼
Hainan GIFT广州吉富鱼
Guangzhou GIFT青岛吉富鱼
Qingdao GIFTUNH985 18.063 16.012* 15.028* UNH871 8.803 9.509 11.218 UNH995 39.337* 95.618* 41.604* UNH855 2.924 0.443 14.927* UNH1007 31.248* 6.009 16.412* UNH846 5.719 9.992 24.058 UNH735 99.374* 18.235 63.103* UNH849 29.923* 93.452* 78.821* UNH954 5.707 67.313* 47.030* 总位点
overall loci39.58 64.71* 68.49* 注:*. P < 0.05
Note:*. P < 0.05表 5 吉富鱼群体间的遗传距离(对角线下方)与遗传分化指数FST(对角线上方)
Table 5 Pairwise Nei′s genetic distance (below diagonal) and FST (above diagonal) among three GIFT stocks
群体 stock 1 2 3 1 海南吉富鱼 HainanGIFT - 0.070* 0.108* 2 广州吉富鱼 GuangzhouGIFT 0.138 - 0.080* 3 青岛吉富鱼 ShanghaiGIFT 0.245 0.183 - 注:*. P < 0.01
Note:*. P < 0.01表 6 吉富鱼群体的AMOVA分析
Table 6 AMOVA of GIFT strain of O.niloticus based on microsatellite marker
变异来源
source of variation自由度
df方差总和
sum of squares变异组分
variance components所占比例/%
percentage of variation群体间 between stocks 2 44.448 0.269 8.74 群体内 within stocks 215 603.909 2.809 91.26 总计 total 217 648.358 3.078 100 表 7 4对引物在吉富鱼群体内检出的扩增片段总数、多态性片段比例和基因多样性
Table 7 Number of bands, percentage of polymorphic bands, and genetic diversity within stocks observed in GIFT strains of O.niloticus
群体
stock位点数
No. of bands多态位点比例
(P0.05)hS VarI% VarL% Hainan GIFT 248 36.3 0.179±0.011 21.0 79.0 Guangzhou GIFT 243 39.9 0.210±0.012 18.6 81.4 Qingdao GIFT 244 48.4 0.245±0.011 26.3 73.7 注:P0.05. 在5%水平上的多态位点比例(%);hS. 群体内基因多样性;S.E.. 标准误差;VarI. 样本方差;VarL. 位点方差
Note:P0.05. proportion of polymorphic loci at 5% level (%);hS. gene diversity within populations;S.E. standard error;VarI. variance due to sampling of individuals;VarL. variance due to sampling of loci表 8 吉富鱼群体遗传变异的分布
Table 8 Partitioning of genetic variation within and between three selected GIFT strains of O.niloticus based on AFLP marker
变异来源
source of variation自由度
df方差总和
sum of squares变异组分
variance components所占比例/%
percentage of variation群体间 between stocks 2 567.171 7.428 32.4 群体内 within stocks 106 1642.994 15.500 67.6 总计 total 108 2210.165 22.928 100 表 9 吉富鱼群体间遗传分化指数FST(对角线下)和遗传距离(对角线上)
Table 9 Pairwise FST(below diagonal)and Reynolds et al genetic distance between stocks (above diagonal)of GIFT strains of O.niloticus
群体 stock 1 2 3 1 海南吉富鱼 HainanGIFT - 0.274 0.329 2 广州吉富鱼 GuangzhouGIFT 0.240* - 0.348 3 青岛吉富鱼 QingdaoGIFT 0.280* 0.294* - 注:*. P < 0.01
Note:*. P < 0.01 -
[1] EKNATH A E, TAYAMEN M M, PALADA-VERA M S, et al. Genetic improvement of farmed tilapias: the growth performance of eight strains of Oreochromis niloticus tested in different farm environments[J]. Aquac, 1993, 111(1/4): 171-188. doi: 10.1016/0044-8486(93)90035-W
[2] 赵金良, 李思发, 李晨虹, 等. 吉富等五品系尼罗罗非鱼生化遗传标志[J]. 上海水产大学学报, 1997, 6(3): 166-171. [3] 杜诚, 卢迈新, 叶星, 等. 五个罗非鱼群体的遗传分析与分子标记[J]. 大连水产学院学报, 2005, 20(1): 26-28. [4] 张传涛, 李思发, 颉晓勇. 吉富品系尼罗罗非鱼选育系F6~F9遗传变异的RAPD分析[J]. 上海水产大学学报, 2007, 16(1): 7-10. doi: 10.3969/j.issn.1004-7271.2007.01.002 [5] RUTTEN M J M, KOMEN H, DEERENBERGR M, et al. Genetic characterization of four strains of Nile tilapia using microsatellite markers[J]. Animal Genetics, 2004, 35(2): 93-97. doi: 10.1111/j.1365-2052.2004.01090.x
[6] ROMANA-EGUIA M R R, IKEDA M, BASIAO Z U, et al. Genetic diversity in farmed Asian Nile and red hybrid tilapia stocks evaluated from microsatellite and mitochondrial DNA analysis [J]. Aquac, 2004, 236(1/4): 131-150. doi: 10.1016/j.aquaculture.2004.01.026
[7] MARIETTE S, LE-CORRE V, AUSTERLITZ F, et al. Sampling within the genome for measuring within-population diversity: trade-offs between markers [J]. Mol Ecol, 2002, 11(7): 1 145-1 156. doi: 10.1046/j.1365-294x.2002.01519.x
[8] 王志勇, 柯才焕, 王艺磊, 等. 从AFLP指纹和标记基因序列看我国养殖的四种鲍的亲缘关系[J]. 高技术通讯, 2005, 15(12): 93-98. doi: 10.3321/j.issn:1002-0470.2004.12.020 [9] 喻达辉, 朱嘉濠. 利用ITS和AFLP标记探讨中国、日本和澳大利亚3种珠母贝的亲缘关系[J]. 南方水产, 2006, 2(5): 36-44. doi: 10.3969/j.issn.2095-0780.2006.05.007 [10] 岳志芹, 王伟继, 孔杰, 等. 用AFLP方法分析中国对虾抗病选育群体的遗传变异[J]. 水产学报, 2005, 29(1): 13-19. doi: 10.3321/j.issn:1000-0615.2005.01.003 [11] 郭奕惠, 黄桂菊, 喻达辉. 合浦珠母贝DNA的抽提和RAPD反应体系的优化[J]. 南方水产, 2006, 2(4): 59-64. doi: 10.3969/j.issn.2095-0780.2006.04.010 [12] 王小玉, 喻达辉, 郭奕惠, 等. 七种珍珠贝RAPD鉴别标记的初步研究[J]. 南方水产, 2006, 2(1): 18-22. https://www.schinafish.cn/cn/article/id/3aa2630e-b18b-49bd-9950-c0eaf166fc74 [13] CARLETON K L, STREELMAN J T, LEE B Y, et al. Rapid isolation of CA microsatellites from the tilapia genome[J]. Anim Genetics, 2002, 33(2): 140-144. doi: 10.1046/j.1365-2052.2002.00817.x
[14] VOS P, HOGERS R, BLEEKER M, et al. AFLP: new technique for DNA fingerpringting[J]. Nucleic Acids Res, 1995, 23(21): 4 407-4 414. doi: 10.1093/nar/23.21.4407
[15] 王小玉, 喻达辉, 黄桂菊, 等. 合浦珠母贝印度家系F1代的AFLP分析[J]. 中国水产科学, 2007, 14(1): 52-58. doi: 10.3321/j.issn:1005-8737.2007.01.008 [16] YEH F C, BOYLE T J B. Population genetic analysis of co-dominant and dominant markers and quantitative traits[J]. Belg J Botany, 1997, 129(2): 157. https://www.semanticscholar.org/paper/Population-genetic-analysis-of-codominant-and-and-Yeh/013886d718112ec35e2585b9de38a35444cf3a67
[17] EXCOFFIER L, LAVAL G, SCHNEIDER S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis[J]. Evolutionary Bioinformatics Online, 2005, 1(3): 47-50. doi: 10.1177/117693430500100003
[18] EXCOFFIER L, SMOUSE P, QUATTRO J. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data[J]. Genetics, 1992, 131(2): 479-491. doi: 10.1093/genetics/131.2.479
[19] VEKEMANS X, BEAUWENS T, LEMAIRE M, et al. Data from amplified fragment length polymorphism (AFLP) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size[J]. Mol Ecol, 2002, 11(1): 139-151. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsRU5HTmV3UzIwMjQwOTEwMTY1MjU1EiA2NGQ5ZjM1NmVjNjYxMmQyZmIxODQxNzZhY2FhZGY3MRoIaHNxbnoya3k%3D
[20] REYNOLES J, WEIR B S, COCKERHAMC C. Estimation of the coancestry coefficient: basis for a short-term genetic distance[J]. Genetics, 1983, 105(3): 767-779. doi: 10.1093/genetics/105.3.767
[21] 韩志强, 高天翔, 王志勇, 等. 黄姑鱼群体遗传多样性的AFLP分析[J]. 水产学报, 2006, 30(5): 640-646. doi: 10.3321/j.issn:1000-0615.2006.05.010 [22] 王志勇, 王艺磊, 林利民, 等. 利用AFLP指纹技术研究中国沿海真鲷群体的遗传变异和趋异[J]. 水产学报, 2001, 25(4): 289-293. doi: 10.3321/j.issn:1000-0615.2001.04.001 [23] MICKETT K, MORTON C, FENG J, et al. Assessing genetic diversity of domestic populations of channel catfish (Ictalurus punctatus) in Alabama using AFLP markers[J]. Aquac, 2003, 228(1): 91-105. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsRU5HTmV3UzIwMjQwOTEwMTY1MjU1EiA2OTNlYjU3MTg5NDc0NzkwZjdiYzk3N2ZkY2M5NGYzNxoIbXd3cXF0ZnU%3D
[24] YUE G H, LI Y, LIM L C, et al. Monitoring the genetic diversity of three Asian arowana (Scleropages formosus) captive stocks using AFLP and microsatellites[J]. Aquac, 2004, 37(1/4): 89-102. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsRU5HTmV3UzIwMjQwOTEwMTY1MjU1EiBiMDg0NjZhNTU4MDU3YWJiODEwYjYwMzA1OWIyYjhjMxoIcXQ4d296aW4%3D
[25] 李斌, 顾万春, 陈晓阳. AFLP标记在生物种质资源研究中的应用[J]. 世界林业研究, 2001, 14(4): 11-17. doi: 10.3969/j.issn.1001-4241.2001.04.002 [26] 夏德全, 曹萤, 吴婷婷, 等. 用RAPD分析对罗非鱼遗传变异的研究及其对杂种优势的应用[J]. 水产学报, 1999, 23(1): 27-32. https://www.cqvip.com/doc/journal/969518668 -
期刊类型引用(3)
1. 吕硕,童剑锋,朱振宏,薛铭华,邱悦,李柏萱,刘必林. 基于声学的2023年夏季西北太平洋小型中上层鱼类时空分布. 上海海洋大学学报. 2025(02): 394-402 . 百度学术
2. 程高,陈国宝,陈丕茂,佟飞,牛麓连,陈钰祥. 基于声学技术定点监测海洋牧场鱼类资源昼夜变化研究. 南方水产科学. 2024(05): 63-70 . 本站查看
3. 王文卓,张纯,薄萍,王浩然,贾烁,王娜娜. 基于基尔霍夫近似模型的南海大黄鱼声学目标强度研究. 南方水产科学. 2024(06): 95-103 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 5454
- HTML全文浏览量: 146
- PDF下载量: 1916
- 被引次数: 3