Effects of two breeding methods on growth and bacterial structure of Nitzschia closterium f. minutissima
-
摘要:
文章以小新月菱形藻 (Nitzschia closterium f. minutissima)为研究对象，分析比较了小新月菱形藻在负压光生物反应器与开放式桶培养下，藻密度、pH、溶解氧及菌落结构的变化情况。结果表明，在负压光生物反应器培养下的藻密度可达到1.33×107个·mL–1，明显高于开放式培养的藻密度 (8.36×106个·mL–1)。藻液中pH随藻密度增加而升高，两者呈显著正相关 (P<0.01)，在负压光生物反应器及开放式培养环境中pH最高值分别为10.3和9.3。溶解氧与pH变化趋势相反，在负压光生物反应器内溶解氧随藻密度增加而降低，最后稳定在6.5 mg·L–1，溶解氧的下降可能与玫瑰杆菌 (Roseobacter)成为优势细菌有关。利用16S rDNA基因的高通量测序技术，分析在培养过程中藻际菌群的结构变化，发现菌落的多样性显著下降 (P<0.05)，培养前期主要以变形杆菌 (Proteobacteria)和拟杆菌 (Bacteroidetes)为优势细菌，在负压光生物反应器内培养后期主要以蓝细菌 (Cyanobacteria)与玫瑰杆菌为优势细菌，其菌落结构与开放式桶存在明显差异。
Abstract:We compared the density, pH, dissolved oxygen and colony structure of Nitzschia closterium f. minutissima bred in negative pressure photobioreactor and open bucket, respectively. It is shown that the algae grew fast in negative pressure photobioreactor with the highest density of 1.33×107 cells·mL–1, significantly higher than that in open bucket (8.36×106 cells·mL–1). The pH in algal solution increased with increase of algal density, and there was a significant positive correlation between them (P<0.01). The highest pH values in negative pressure photobioreactor and open bucket were 10.3 and 9.3, respectively. Oppositely, in negative pressure photobioreactor, the dissolve oxygen of algea decreased with increase of algal density, which finally stabilized at about 6.5 mg·L–1. The decline of dissolved oxygen might be related to the fact that Roseobacter had become the dominant bacteria. Results of 16S rDNA gene high-throughput sequencing show that the bacterial diversity decreased significantly (P<0.05). Proteobacteria and Bacteroidetes were dominant bacteria at early stage of breeding. However, the dominant bacteria were Cyanobacteria and Roseobacter at late stage in negative pressure photobioreactor, which was obviously different from that in open bucket.
-
疣吻沙蚕 (Tylorrhynchus heterochaetus) 又名疵吻沙蚕，属环节动物门、多毛纲、叶须虫目、沙蚕科，俗称禾虫、流蜞，在中国东南沿海河口地区的泥沙质浅滩或稻田中广泛分布[1-2]。作为一种经济价值极高的多毛类，其味道鲜美且营养丰富，是中国广东、广西、福建、香港、澳门及东南亚各地独具特色的水产品，素有“水中冬虫夏草”之美誉[2-3]。更重要的是，疣吻沙蚕与水稻具有天然的共生关系，发展疣吻沙蚕稻田综合种养模式可显著增加水稻种植的经济效益[4]。“水稻+疣吻沙蚕”生态综合种养技术正日趋成熟，自2021年起连续3年入选广东省农业主推技术，在助力乡村振兴、保障国家粮食安全等方面应用潜力巨大。然而，疣吻沙蚕具有生境狭窄、种群地理隔离、群体恢复力低等物种特性，受栖息地破坏、过度捕捞、环境污染等因素的影响其自然种群呈逐年衰退的趋势，在有些地区甚至已经绝迹[4-5]。在突破疣吻沙蚕全人工繁殖技术的基础上，近年已实现其苗种的规模化培育，增养殖业的发展以及技术的革新培育了大量的人工繁育群体，但因遗传育种研究的滞后，疣吻沙蚕种质开始出现生长性能及抗逆性下降等退化现象[4]。为建立疣吻沙蚕资源管理策略并促进其可持续开发利用，特别是避免常年的产业化增养殖引起遗传结构单一化，防止人工繁育群体污染自然种群，有必要采用分子手段科学地分析及评估其种质现状，揭示种群遗传多样性水平与遗传结构，并以此指导新品种的培育与遗传改良研究。
目前，国内外针对疣吻沙蚕的研究主要集中在形态学、生活史、繁殖生物学、增养殖技术、营养与活性成分分析等方面[3,5]。Chen等[5-6]完成了疣吻沙蚕线粒体基因组测序，并基于线粒体COI序列分析了7个地理群体的遗传结构，除此之外，遗传学相关的研究鲜有报道；分子遗传信息的匮乏制约了疣吻沙蚕资源保护利用研究的深入开展。微卫星标记具有多态性丰富、共显性遗传、基因组覆盖广泛等优点[7]，该技术操作简便、检测周期短、稳定性高、重复性好，已广泛应用于谱系鉴定、亲缘关系检测、遗传多样性分析、种质资源评价、遗传图谱构建及基因定位等领域[8]。基因组survey是利用高通量测序技术实施小片段、低深度测序，然后基于Lander-Waterman模型进行k-mer分析，并根据k-mer频率和深度的统计结果评估物种基因组的大小与复杂程度等关键特征[9]。基因组survey不仅能为全基因组测序及高质量组装提供科学依据，而且对无参考基因组序列的物种而言，利用基因组survey数据大规模开发微卫星标记是当前最高效的策略之一，较传统方法具有效率高、周期短、成本低的优势[9]。基因组微卫星鉴定研究在水产经济动物中已广泛开展，筛选的标记获得了良好的遗传分析效果。例如，张永德等[9]在卵形鲳鲹 (Trachinotus ovatus) 基因组survey数据中检测到190 121个微卫星位点并成功开发了多态分子标记，为全基因组测序与组装、渔业资源保护利用及良种选育奠定了基础；上官清等[10]分析了斑鳢 (Channa maculata) 基因组微卫星特征，并筛选到20个多态性位点用于群体遗传多样性和遗传结构分析，为该物种的遗传监测、亲缘关系鉴定、种质资源养护及管理提供了技术支持。
本研究对疣吻沙蚕基因组进行低深度高通量测序，通过k-mer分析预测基因组大小、杂合度和重复比例等信息，对测序数据初步组装后搜索组装序列中的微卫星位点，分析其特征与分布规律，初步验证标记的有效性和多态性，以期指导基因组精细图谱的绘制，并为群体遗传学研究提供可靠的标记资源。
1. 材料与方法
1.1 实验材料
疣吻沙蚕采集自阳江市广东阳海农业技术发展有限公司疣吻沙蚕增养殖试验基地(111°55'23"E、 21°49'16"N) 保种的越南海防群体。用于基因组测序的疣吻沙蚕幼体用适量无菌水洗涤3次，解剖后剪取体壁肌肉组织装入2 mL冻存管，液氮速冻后置于−80 ℃保存。微卫星标记多态性验证群体的30尾个体经无菌水洗涤后，解剖剪取体壁肌肉组织，置于体积分数为95%的乙醇中在−20 ℃下保存。
1.2 基因组测序及生物信息学分析
运用苯酚-氯仿法提取疣吻沙蚕基因组DNA，经1.0% (w) 琼脂糖凝胶电泳检测完整度后，用NanoDrop 2000超微量分光光度计 (ThermoFisher，美国) 检测浓度和纯度。利用Covaris超声波破碎仪将基因组DNA片段化，筛选合适长度的DNA片段，经末端修复、加A尾、添加测序接头、纯化及PCR扩增等步骤建立350 bp小片段文库。测序文库经Qubit 2.0 (ThermoFisher，美国) 和Agilent 2100 Bioanalyzer (Agilent，美国) 检测浓度和插入片段大小后，用Illumina HiseqTM X Ten平台进行双末端测序。基因组测序工作委托北京诺禾致源科技股份有限公司完成。原始数据经质控和过滤后，用GCE 1.0.0[11]软件对有效数据进行k-mer分析。运用SOAPdenovo 2.01[12]软件，选择k-mer=41将有效数据组装至contig和scaffold级别，并统计GC含量和覆盖深度等信息。
1.3 微卫星位点搜索与引物设计
使用MISA (Microsatellite identification tool) 软件(http://pgrc.ipk-gatersleben.de/misa/misa.html)在长度大于500 bp的序列中搜索微卫星位点。运行参数为：重复基序长度1~6 bp，单碱基重复次数≥12次，二碱基重复次数≥6次，三碱基、四碱基重复次数≥5次，五碱基、六碱基重复次数≥4次；将间隔区域长度小于100 bp的相邻微卫星归为1个复合型位点。基于微卫星侧翼序列，用Primer 5[13]批量设计引物，主要参数为：引物长度18~27 bp，扩增产物100~300 bp，退火温度(Tm) 55~65 ℃，GC含量40%~60%，正、反向引物退火温差≤5 ℃；尽量避免出现发卡结构、二聚体、错配和引物二聚体；每个微卫星位点生成3~5对候选引物。随机选取50对微卫星引物委托生工生物工程(上海)股份有限公司合成。
1.4 多态微卫星标记的验证与筛选
使用天根生化科技(北京)有限公司的海洋动物组织基因组DNA提取试剂盒从疣吻沙蚕体壁肌肉组织中提取基因组DNA。PCR反应体系为20 μL，正、反向引物 (10 μmol·L−1) 各0.5 μL，DNA模板 (15 ng·μL−1) 1.0 μL，2×PCR Mix 10 μL，用超纯水补至20 μL。扩增反应由Bio-Rad My Cycler Thermal Cycler (Bio-Rad，美国) 完成，运行程序为：95 ℃预变性5 min；95 ℃变性30 s，59~60 ℃退火30 s，72 ℃延伸30 s，25个循环；72 ℃延伸5 min。PCR产物先进行1.5% (w) 琼脂糖凝胶电泳分析，参照预期片段大小筛选出有效扩增引物。有效扩增引物加荧光接头后 (正向引物5'端添加FAM荧光素)，在30尾疣吻沙蚕个体的基因组DNA中扩增以验证其多态性。PCR产物送至上海翼禾应用生物技术有限公司用3730XL测序分析仪 (Applied Biosystems，美国) 进行毛细管电泳分析，使用GeneMapper 3.2软件 (Applied Biosystems，美国) 进行基因分型。
1.5 数据分析
使用Excel 2016软件完成微卫星位点分布特征信息的统计分析和图表绘制。微卫星发生频率=含微卫星的序列总数/序列总数×100%；微卫星出现频率=微卫星总数/序列总数×100%；微卫星丰度 (个·Mb−1)=微卫星总数/序列总长度[14]。用GenAlEx 6.5软件[15]计算等位基因数 (Na)、有效等位基因数 (Ne)、观测杂合度 (Ho)、期望杂合度 (He)、多态信息含量 (PIC)。用Genepop在线软件 (https://genepop.curtin.edu.au/) 检验位点间的连锁不平衡及群体的哈迪-温伯格平衡 (Hardy-Weinberg Equilibrium, HWE)，并用Bonferroni法对显著性阈值进行校正。
2. 结果
2.1 疣吻沙蚕的基因组survey测序
低深度高通量测序产生的原始数据经质控后共获得57.48 Gb有效数据，碱基错误率为0.05%，Q20和Q30分别为95.56%和89.70%，表明基因组测序质量较高 (表1)。对有效数据进行k-mer分析 (k=17)，结果显示在深度为57时出现主峰值 (图1)，总k-mer为44 257 233 158，排除错误k-mer的误差影响后得到修正的基因组大小为759.53 Mb，杂合率为1.41%，重复序列比例为45.92%。选择k-mer=41将有效读段初步组装至contig和scaffold水平，最终获得的contig总长度为821 637 022 bp，最大长度为84 713 bp，N50为548 bp；scaffold总长度为840 375 821 bp，最大长度为89 326 bp，N50为662 bp。
表 1 疣吻沙蚕基因组 survey 测序数据统计Table 1. Statistics of genomic survey sequencing data of T. heterochaetus测序文库
Sequencing library原始数据量
Raw base/Gb有效数据比
Effective rate/%有效数据量
Clean base/Gb碱基错误率
Error rate/%Q20/% Q30/% GC 含量
GC content/%L1 33.12 99.68 33.01 0.04 96.23 90.77 39.13 L2 24.53 99.75 24.47 0.05 94.88 88.63 39.00 总计 Total 57.65 — 57.48 — — — — 均值 Mean — 99.72 — 0.05 95.56 89.70 39.07 ![]()
图 1 疣吻沙蚕基因组 k-mer 种类频率分布Figure 1. Frequency distribution of k-mer species in genome of T. heterochaetus2.2 疣吻沙蚕基因组微卫星位点特征
用MISA软件在109 881条组装序列中检测到 130 216个微卫星位点，长度共计2 341 179 bp；微卫星发生频率5.04%，出现频率5.97%，分布丰度为154.9个·Mb−1。疣吻沙蚕微卫星位点以单碱基和二碱基重复最为丰富，分别有45 582和42 298条，各占35.00%和32.48%；其次是三碱基重复 (18 782条)，占14.42%；六碱基数量最少，仅占2.44% (图2)。微卫星序列的重复数范围为4~56拷贝，主要集中在4~18拷贝，重复19次及以上的有3 703条，仅占3.50%；单碱基重复以12~16次最常见 (40 706条，89.30%)，二碱基重复集中在6~18次 (42 035条，99.38%)，三碱基重复以5~12次为主 (18 526条，98.64%)，四碱基重复以5~10次为主 (12 173条，98.91%)，而五碱基和六碱基重复主要为4~8次 (共11 176条，99.37%) (图2—图3)。
![]()
图 2 疣吻沙蚕基因组 6 种类型微卫星的数量与比例Figure 2. Number and proportion of six motif types of microsatellite loci in genome of T. heterochaetus![]()
图 3 疣吻沙蚕基因组微卫星重复数分布特征Figure 3. Distribution pattern of microsatellite repeat number in genome of T. heterochaetus疣吻沙蚕基因组微卫星共包含320种重复基序类型：单碱基2种、二碱基4种、三碱基10种、四碱基31种、五碱基91种、六碱基182种。单碱基重复基序拷贝数集中在12~15次，以C/G为主，占比58.02%；二碱基重复拷贝数多为6~10次，以AT/AT最为丰富，占比62.38%，其次是AC/GT (27.37%)、AG/CT (10.22%)，而CG/CG数量稀少 (0.02%)；在三碱基重复中，拷贝数多为5~12次，占比最高的是AAT/ATT，有6253条 (33.29%)，其次是ATC/GAT (32.15%)，而CCG/CGG数量最少，仅占0.06%；四碱基拷贝数多为5~8次，优势重复基序是AAAT/ATTT，有3567条 (28.98%)，其次是AATC/GATT (14.08%)、ACTC/GAGT (9.60%)、ATCC/GGAT (7.38%)；五碱基、六碱基拷贝数集中在4~7次，优势重复基序分别是AAAAT/ATTTT (17.07%) 和AACCCT/AGGGTT (11.73%) (图4)。
![]()
图 4 疣吻沙蚕基因组微卫星重复基序类型分布特征Figure 4. Distribution pattern of microsatellite motif types in genome of T. heterochaetus2.3 疣吻沙蚕基因组微卫星位点长度特征
根据重复序列长度可将微卫星位点分为两类：一类是长度达20 bp及以上的高度多态I型，另一类则是长度介于12~19 bp的中度多态II型[16]。疣吻沙蚕基因组微卫星位点的长度分布区间为12~336 bp，I型微卫星位点有43 240条，占33.21%，剩余的II型位点占66.79%；绝大部分I型微卫星位点的长度为20~39 bp，超过50 bp的位点仅占0.51% (图5-a)。进一步分析发现，在I型微卫星位点中，二、三、四、五碱基是最主要的重复类型，占85.92%，在后续多态标记的筛选中有较高的开发价值(图5-b)。
![]()
图 5 疣吻沙蚕基因组微卫星长度分布特征注：a. 不同长度区间微卫星数量及比例；b. 不同类型微卫星长度分布特征。Figure 5. Distribution pattern of length of microsatellite loci genome of T. heterochaetusNote: a. Number and percentage of microsatellite loci at different length intervals; b. Length distribution of the six motif types of microsatellite loci.2.4 疣吻沙蚕基因组微卫星多态标记验证筛选
使用Primer 5成功对37 370个微卫星位点设计了引物。随机选取50个位点合成引物并进行PCR验证，共获得41对 (82%) 有效扩增引物，表明MISA软件鉴定的微卫星位点具有较高有效性。多态性筛选结果显示，15个 (30%) 微卫星位点的引物表现出稳定且可重复的多态性 (表2，图6)。在30尾疣吻沙蚕中，15对引物共检测到87个等位基因，平均等位基因数5.800，等位基因频率为0.060~0.400，其中ThGM021位点的等位基因数最少 (Na=2.000)，ThGM004位点的等位基因数最多 (Na=12.000)。Ne为1.164~6.713，平均值为3.328；Ho为0.050~0.879，平均值为0.487；He为0.141~0.789，平均值为0.561；PIC为0.136~0.776，平均值为0.511 (表3)。在15个位点中，8个属高度多态性位点 (PIC>0.5)，5个属中度多态性位点 (0.25<PIC≤0.5)，2个属低度多态性位点 (PIC≤0.25)，表明该群体的遗传多样性较丰富。经Bonferroni校正后，有3个位点 (ThGM006、ThGM011、ThGM040) 偏离HWE，其中ThGM040属于高度多态性位点。各位点间无连锁不平衡现象。
表 2 疣吻沙蚕 15 对多态微卫星引物信息Table 2. Information of 15 polymorphic microsatellite loci in genome of T. heterochaetus位点
Locus引物序列 (5'—3')
Primer sequence (5'–3')重复单元
Repeat unit产物大小
Size/bp退火温度
Annealing temperature/℃ThGM004 F: TGCTGCTACTGCTACAGCTACTATG (TAC)18 289 60.0 R: CTGACAAAGTTTGGTGGCTG ThGM006 F: TGAAAATTAGTGTGATTTTGTCCC (CA)11 260 59.0 R: AGCCAACCAGAACATGAACA ThGM011 F: AACTTGGACTAAGGCTATCAAAAA (AG)17 220 59.0 R: CTTGGGGTTCATGCATCATT ThGM015 F: TTGGTTGTTATCCATGCACC (TAT)12 279 59.5 R: AGACAGCAGTGAAATAGCACCA ThGM017 F: ATTCGATAAGCATTCCACCG (ATGG)8 215 60.0 R: CTTGGTAGCTGGCCTGTCTC ThGM021 F: TGCGAAATGAGAAGTGAGCA (TA)10 277 60.0 R: TGCCTGTGTGGAATACCAAG ThGM024 F: ACCTGTCCACCCGTCATTTA (TAT)14 294 59.5 R: CCTTTAGGGGATGGCTACAA ThGM029 F: GAGCAAAATATTCAAGTTGGCA (ATT)12 243 59.0 R: TTGTTTGTCATATCTTCTAAAGAGCA ThGM033 F: GGAGTGGGGAGGATTTTAGC (TG)18 277 60.0 R: CCATGTACAGCATTCAGCCA ThGM035 F: GTAAGGGCAAGGGTTGTGAA (AG)13 226 60.0 R: ACCGTTACCCTAACCCCAAC ThGM038 F: TTACCCTGCCATCCTACCAG (TG)20 157 60.0 R: CTATTCTGCCAGTGGTCGCT ThGM040 F: GGATCCAGAAGGGGTAAAGC (TTA)11 239 59.5 R: GTTGGTCATGTTCCTGTTGC ThGM041 F: ACCAGCTGCTAGAGGCAGAC (ATG)7 260 60.0 R: TTAGGTCCTCACCCAGGGAT ThGM043 F: AAAAGCAAGTGGTAACACAAAATG (TCAT)11 272 59.5 R: CATTGGGCTCTGGGAATAAA ThGM047 F: CGACCTGCGGATTTAATTTG (TGG)12 148 60.0 R: ATATCTTGGCGGCGGATAG 注：F. 正向引物；R. 反向引物。 Note: F. Forward primer; R. Reverse primer. ![]()
图 6 部分多态微卫星标记的毛细管电泳分型结果Figure 6. A set of polymorphic microsatellite loci visualized by high-resolution capillary electrophoresis表 3 15 个多态微卫星位点在疣吻沙蚕群体中的遗传特征Table 3. Genetic characteristics of 15 polymorphic microsatellite loci in a T. heterochaetus population位点 Locus 等位基因数 Na 有效等位基因数 Ne 观测杂合度 Ho 期望杂合度 He 多态信息含量 PIC 哈迪-温伯格平衡的P 值 PHWE ThGM004 12 6.672 0.697 0.775 0.726 0.275 ThGM006 4 1.642 0.367 0.388 0.372 0.001* ThGM011 5 1.608 0.257 0.349 0.377 0.026* ThGM015 8 4.933 0.438 0.789 0.776 0.225 ThGM017 3 1.521 0.066 0.271 0.245 0.148 ThGM021 2 1.593 0.367 0.508 0.375 1.000 ThGM024 11 6.713 0.879 0.742 0.682 0.541 ThGM029 10 5.647 0.697 0.658 0.599 0.140 ThGM033 5 3.102 0.167 0.772 0.720 0.069 ThGM035 3 2.164 0.050 0.141 0.136 0.086 ThGM038 8 4.878 0.576 0.545 0.489 0.221 ThGM040 4 2.441 0.733 0.718 0.652 0.008* ThGM041 3 2.155 0.417 0.431 0.336 0.503 ThGM043 5 2.727 0.724 0.682 0.617 0.148 ThGM047 4 2.224 0.867 0.642 0.569 0.267 均值 Mean 5.800 3.328 0.487 0.561 0.511 — 注：*. Bonferroni法校正后显著偏离哈迪-温伯格平衡(P<0.05)；n=30。 Note: *. Significant departure from Hardy-Weinberg equilibrium after Bonferroni's correction (P<0.05); n=30. 3. 讨论
3.1 疣吻沙蚕基因组基本特征
不同物种间基因组的大小和复杂程度有明显差异，会直接影响到测序策略的选择及基因组的组装效果。因此，进行全基因组测序前须先评估物种基因组的基本特征。疣吻沙蚕基因组survey测序及k-mer分析估计其基因组为759.53 Mb，远大于水蛭 (Helobdella robusta)[17]和宽体金线蛭 (Whitmania Pigra)[18]等所有已报道的蛭纲物种，也较多毛纲的海蠕虫 (Capitella teleta)[17]、欧文虫 (Owenia fusiformis)[19]、巨型管虫 (Riftia pachyptila)[20]、Lamellibrachia luymesi[21]的大；与寡毛纲通俗腔蚓 (Metaphire vulgaris)[22]及多毛纲搓稚虫 (Streblospio benedicti)[23]的大小 (约0.7 Gb) 相当；但明显小于安德爱胜蚓 (Eisenia andrei)[24]、赤子爱胜蚓 (E. fetida)[25]，以及多毛纲的旋鳃虫(Spirobranchus lamarcki)[26]和深海管虫 (Paraescarpia echinospica)[27] (1.0~1.3 Gb)。物种间基因组大小的差异性在一定程度上与重复序列比例的高低有关。疣吻沙蚕基因组的重复序列比例 (45.92%) 低于旋鳃虫[26]和深海管虫[27]，但显著高于海蠕虫[17]、巨型管虫[20]，以及所有已公布的蛭纲动物 (重复比例均低于34%)[15-18]；其基因组杂合率 (1.41%) 与安德爱胜蚓[24]和赤子爱胜蚓[25]相当，远远高于L. luymesi (0.60%)[21]、搓稚虫 (0.29%)[23]与深海管虫 (0.63%)[27]，是目前已报道的杂合率最高的多毛类。此外，组装的contig和scalffold总长分别为821 637 022和840 375 821 bp，N50分别为548和662 bp，提示过高的杂合率会造成组装序列长度大于预估基因组大小，并导致显著偏低的N50指标[28]。综上所述，疣吻沙蚕基因组属复杂基因组类型，绘制染色体级别高质量全基因组图谱应优先考虑“PacBio+Illumina+Hi-C”策略。该研究数据为后续全基因组测序与组装提供了基础资料。
3.2 疣吻沙蚕基因组微卫星分布特征
水生动物基因组中单碱基重复微卫星占优势的现象较少，已报道的有中华绒螯蟹 (Eriocheir sinensis)[29]、鲤 (Cyprinus carpio)[30]、脊尾白虾 (Exopalaemon carinicauda)[31]和胡鲶 (Clarias batrachus)[32]。在疣吻沙蚕基因组中检测到 130 216个微卫星位点，位点数量随重复基序碱基数的增加而迅速减少，单碱基重复微卫星最丰富 (35.00%)，其次为二碱基重复 (32.48%)、三碱基重复 (14.42%)，这与以二碱基或三碱基为优势类型的鱼[9,33-35]、虾[36-37]、贝类[38]明显不同，特别是与亲缘关系较近的蛭类、寡毛类相比也表现出较大差异。宽体金线蛭[39]、天锡杜拉蚓 (Drawida gisti)[40]基因组中三碱基重复占绝对优势，主要基序类型分别是AAT、ATA和ATT、AAT。王斌等[41]对4种蛭类的转录组数据进行分析，同样发现三碱基重复微卫星占优势，在宽体金线蛭中甚至高达68.00%。作为最特殊的微卫星重复类型，三碱基重复微卫星可形成复杂的环-折叠构型来稳定DNA结构，从而更有利于转录过程中的解旋和蛋白质识别[33]。此外，疣吻沙蚕转录组中二碱基AT/TA重复基序频率 (32.19%) 最高，单碱基A/T重复类型(21.17%)其次，之后是三碱基AGC/GCT (29.06%) 基序 (未发表数据)，表明微卫星类型特征在转录组和基因组水平上存在差异，类似的差异性在圆鳍鱼 (Cyclopterus lumpus)[42]、凡纳滨对虾[36,43]、宽体金线蛭[39,41]中也被证实。不同物种表现出不同的微卫星优势类型分布规律，可能与其进化水平有关。
疣吻沙蚕基因组微卫星核心区集中在4~18拷贝，随着重复次数的增加，微卫星数量呈显著降低趋势，这可能与重复单元长度的不断增加使得微卫星稳定性降低或基序高频次重复导致更高的突变率有关[44]。进一步分析发现，二碱基至六碱基微卫星重复单元中A/T含量显著高于G/C含量，即微卫星序列表现出明显的A/T碱基优势，这与许多水生动物如脊尾白虾[31]、凡纳滨对虾[36]、虾夷扇贝 (Mizuhopecten yessoensis)[38]、宽体金线蛭[39]及旧金山湾卤虫 (Artemia franciscana)[45]中的研究结果一致。早期观点认为，微卫星富含A/T的原因可能是由于CpG甲基化后的C易脱氨基转变为T，导致G/C比例不断缩小，而突变的A/T碱基类型相应增多[33]；后来却发现这也可能与微卫星位点的产生方式，即与DNA的复制滑动存在一定关系[29,46]。微卫星富含A/T的原因可能是A/T含量高则Tm值降低，其序列容易发生DNA解链，并通过复制滑动机制和重组机制产生高A/T含量的重复类型的概率更高[38]。
3.3 疣吻沙蚕微卫星标记开发及其遗传分析效果
微卫星是重要的分子遗传学研究工具，开发多态标记是其广泛应用的关键。与传统方法相比，基于高通量测序技术开发微卫星标记更具优势，产生的基因组或转录组数据是标记开发的重要资源。然而与转录组微卫星相比，基因组微卫星的多态性往往更高且分布更广泛，可获得更优的基因组覆盖率[42]。Temnykh等[16]认为长度≥20 bp的微卫星位点多态性较高，长度10~20 bp的为中等多态性，小于10 bp的多态性低。在疣吻沙蚕基因组中共筛选到43 240条序列长度≥20 bp的I型微卫星位点，占比超过33%，远高于其转录组序列中20%的I型位点比例 (未发表数据)。挑选50个微卫星位点进行有效性及多态性验证，有41个 (82%) 位点的引物可扩增出特异性条带，其中15个 (30%) 表现出稳定且可重复的多态性。疣吻沙蚕基因组多态微卫星标记的筛选成功率与卵形鲳鲹[9]、凡纳滨对虾[36]、宽体金线蛭[39]中的结果相当。由此可见，疣吻沙蚕微卫星标记筛选成功率较高，获得的微卫星位点具有良好的开发潜力，是丰富且可靠的标记资源。
用微卫星标记开展群体遗传分析时，一般认为有效等位基因数越接近观测等位基因数，群体等位基因分布就越均匀；然而在实际分析中，通常会将全部条带视为有效等位基因，无效等位基因过剩便会造成等位基因分布不均[47]。从疣吻沙蚕基因组筛选的15个多态微卫星位点中，仅3个位点 (ThGM021、ThGM035、ThGM041) 的等位基因数接近有效等位基因数，提示大部分位的等位基因分布不均，这可能是因为验证群体的样本容量偏小引起了主效等位基因的缺失[9]。有研究表明，群体的最小样本容量一般与遗传分析选用的参数有关。利用He、PIC以及香农指数等评估遗传多样性时，群体样本容量达到27便能使分析结果接近总体水平的95%；但选择Na时，样本容量则需达到52以上[48]。疣吻沙蚕的群体样本量为30，故而选择He、PIC可更好地评价微卫星位点的多态性。15个位点的平均Ho (0.487)、平均He (0.561)和平均PIC (0.511)共同表明，疣吻沙蚕验证群体具有较丰富的遗传多样性，这与Chen等[5]利用COI标记的分析结果并不完全一致，可能是因为微卫星属于核标记，多态性更高且受选择性作用更小，其揭示遗传变异的灵敏度高于线粒体标记。通常认为，PIC越接近1，则群体杂合个体的比例越大、多态性越高：PIC<0.25为低度多态性，0.25≤PIC<0.5为中度多态性，PIC≥0.5为高度多态性[47]。在15个多态微卫星位点中，仅2个为低度多态性，其余13个 (86.7%) 均为高度或中度多态性。因此，基于疣吻沙蚕基因组序列筛选的微卫星标记多态性高且遗传信息丰富，可为群体遗传分析、种质资源评价、分子育种研究提供优质工具。
4. 结论
疣吻沙蚕基因组属于高杂合、高重复的复杂基因组，其微卫星位点的类型丰富且具备较好的多态性潜能，可作为有效资源用于微卫星标记的大规模开发，对种质资源评价与保护利用、种群遗传学以及分子育种研究具有实际价值。
-
![]()
图 2 小新月菱形藻密度与比生长率变化
Figure 2. Variation of density and specific growth rate of Nitzschia closterium f. minutissima
![]()
图 3 不同培养环境pH与溶解氧变化
Figure 3. Change of pH and dissolved oxygen in different culture environments
![]()
图 8 藻液菌群PCoA分析
Figure 8. Principal coordinates analysis of microflora in algae solution
表 1 负压反应器内部参数
Table 1 Internal parameters of negative pressure photobioreactor
直流电压/V
DC voltage负载电流/A
load current功率/W
power流量/L·min−1 flow 相对真空度/kPa
relative vacuum峰值流量 平均流量 24 <0.8 <19.2 20 13 ≈ −60 
下载: 导出CSV
表 2 不同培养环境藻液菌群α多样性指数比较
Table 2 α diversity index of microflora in algae solution in different breeding environments
组别
treatment多度涵盖估计量
Ace丰富度
Chao1香农指数
Shannon辛普森指数
Simpson开放式培养
open culture第1天 AT 71.9±2.21 71.6±1.18 4.23±0.08a 0.91±0.01a 第7天 BT 75.5±3.57 76.1±5.44 3.59±0.23ab 0.82±0.03ab 第10天 CT 73.1±5.78 75.2±8.50 3.03±0.04b 0.71±0.02b 第14天 DT 75.9±4.05 74.7±3.62 3.30±1.02ab 0.75±0.17ab 反应器培养
reactor culture第1天 AC 77.6±2.60 80.9±10.5 4.25±0.22a 0.91±0.01a 第7天 BC 73.4±1.96 72.8±2.00 3.87±0.20a 0.86±0.02a 第10天 CC 77.2±2.94 81.0±3.46 2.98±0.02b 0.69±0.01b 第14天 DC 71.1±4.36 73.2±5.01 2.78±0.29b 0.74±0.07b 注：同列数据上标不同字母表示平均值差异显著 (P<0.05)
Note: The data in the same column with different superscripts are significantly different (P<0.05).
下载: 导出CSV
-
[1] 李炳乾, 刘颖芬, 刘洪岩, 等. 小新月菱形藻生长条件及半连续培养条件研究[J]. 水产科技情报, 2012(2): 55-58. doi: 10.3969/j.issn.1001-1994.2012.02.001 [2] 刘娟妮, 胡萍, 姚领, 等. 微藻培养中光生物反应器的研究进展[J]. 食品科学, 2006, 27(12): 772-777. doi: 10.3321/j.issn:1002-6630.2006.12.203 [3] 张芬芬, 马晓建, 常春, 等. 气升式微藻光生物反应器的设计研究进展[J]. 现代化工, 2016, 36(10): 46-49. [4] 刘玉环, 黄磊, 王允圃, 等. 大规模微藻光生物反应器的研究进展[J]. 生物加工过程, 2016, 14(1): 65-73. doi: 10.3969/j.issn.1672-3678.2016.01.013 [5] 周进, 林光辉, 蔡中华. 微生物在藻际环境中的物质循环作用[J]. 应用生态学报, 2016, 27(8): 2708-2716. [6] 王少沛, 曹煜成, 李卓佳, 等. 水生环境中细菌与微藻的相互关系及其实际应用[J]. 南方水产, 2008, 4(1): 76-80. doi: 10.3969/j.issn.2095-0780.2008.01.014 [7] 熊青山, 潘秋玲, 王琴, 等. 溶藻微生物菌群的富集及其溶藻因素[J]. 应用与环境生物学报, 2016, 22(6): 1140-1144. [8] SHIN H, LEE E, SHIN J, et al. Elucidation of the bacterial communities associated with the harmful microalgae Alexandrium tamarense and Cochlodinium polykrikoides using nanopore sequencing[J]. Sci Rep, 2018, 8(1): 5323. doi: 10.1038/s41598-018-23634-6
[9] WANG H, HILL R T, ZHENG T, et al. Effects of bacterial communities on biofuel-producing microalgae: stimulation, inhibition and harvesting[J]. Crit Rev Biotechnol, 2014, 36(2): 1-12.
[10] JUAN F, INÉS G, MARÍA C, et al. Impact of microalgae-bacteria interactions on the production of algal biomass and associated compounds[J]. Mar Drug, 2016, 14(5): 100. doi: 10.3390/md14050100
[11] 贾雨川, 阎华, 关昌峰, 等. 内置组合转子新型管式光生物反应器及微藻光暗周期研究[J]. 现代化工, 2017, 37(9): 168-170. [12] De los DANIEL C V, GARCíA-CRUZ E L, FRANCOMORGADO M, et al. Short-term evaluation of the photosynthetic activity of an alkaliphilic microalgae consortium in a novel tubular closed photobioreactor[J]. J Appl Phycol, 2016, 28(2): 795-802. doi: 10.1007/s10811-015-0612-7
[13] 张成会. 利用管式光生物反应器培养湛江等鞭金藻的研究[D]. 湛江: 广东海洋大学, 2015: 32-40. [14] 肖玉朋, 戴玉杰, 申世刚, 等. 户外管式光生物反应器培养发状念珠藻细胞[J]. 食品科技, 2014, 39(3): 2-6. [15] 朱军保. 沙漠微藻管式光生物反应器研制及其相关研究[D]. 石河子: 石河子大学, 2015: 41-48. [16] ZHU J Y, RONG J F, ZONG B N. Factors in mass cultivation of microalgae for biodiesel[J]. Chin J Catal, 2013, 34(1): 80-100. doi: 10.1016/S1872-2067(11)60497-X
[17] 万晓安, 杨正健, 杨林. 光生物反应器中微藻生长影响因子研究进展[J]. 应用化工, 2016, 45(6): 1140-1145, 1154. [18] 石娟, 潘克厚. 不同光照条件对小新月菱形藻和等鞭金藻8701生长及生化成分的影响[J]. 中国水产科学, 2004, 11(2): 121-128. doi: 10.3321/j.issn:1005-8737.2004.02.007 [19] 梁英, 刁永芳, 陈书秀, 等. 温度对小新月菱形藻叶绿素荧光特性及生长的影响[J]. 水产科学, 2011, 30(8): 435-440. doi: 10.3969/j.issn.1003-1111.2011.08.001 [20] 张海阳. 基于射流流场的微藻混凝共聚气浮采收基础研究[D]. 徐州: 中国矿业大学, 2013: 44. [21] 黄新雪. 螺旋转子管式光生物反应器流体动力学特性研究[D]. 北京: 北京化工大学, 2017: 39. [22] 游亮, 崔莉凤, 刘载文, 等. 藻类生长过程中DO、pH与叶绿素相关性分析[J]. 环境科学与技术, 2007, 30(9): 42-44. doi: 10.3969/j.issn.1003-6504.2007.09.015 [23] VALDES F J, HEMANDEZ M R, CATALA L, et al. Estimation of CO2 stripping/CO2 microalgae consumption ratios in a bubble column photobioreactor using the analysis of the pH profiles.Application to Nannochloropsis oculata microalgae culture[J]. Bioresour Technol, 2012, 119: 1-6. doi: 10.1016/j.biortech.2012.05.120
[24] VADLAMANI A, VIAMAJALA S, PENDYALA B, et al. Cultivation of microalgae at extreme alkaline pH conditions: a novel approach for biofuel production[J]. ACS Sustainable Chem Eng, 2017, 5(8): 7284-7294. doi: 10.1021/acssuschemeng.7b01534
[25] LAKANIEMI A M, HULATT C J, WAKEMAN K D, et al. Eukaryotic and prokaryotic microbial communities during microalgal biomass production[J]. Bioresour Technol, 2012, 124: 387-393. doi: 10.1016/j.biortech.2012.08.048
[26] LURIA C M, AMARAL-ZETTLER L A, DUCKLOW H W, et al. Seasonal shifts in bacterial community responses to phytoplankton-derived dissolved organic matter in the western antarctic peninsula[J]. Front Microbiol, 2017, 8: 2117. doi: 10.3389/fmicb.2017.02117
[27] 王剑. 典型赤潮藻类可培养藻际细菌的分子鉴定与藻菌关系研究[D]. 广州: 暨南大学, 2014: 41. [28] 苗祯. 北极微藻藻际细菌类群结构及其相互作用机制研究[D]. 济南: 山东大学, 2013: 42-44. [29] 苗祯, 杜宗军, 李会荣, 等. 5株北极微藻藻际环境的细菌多样性[J]. 生态学报, 2015, 35(5): 1587-1600. [30] FANG F, GAO Y, GAN L, et al. Effects of different initial pH and irradiance levels on cyanobacterial colonies from Lake Taihu, China[J]. J Appl Phycol, 2018, 30(3): 1777-1793. doi: 10.1007/s10811-018-1394-5
[31] ZHANG Z M, YU Z D, ZHU L, et al. Gradient reduced aeration in an enhanced aerobic granular sludge process optimizes the dominant microbial community and its function[J]. Environ Sci, 2018, 4(5): 680-688.
[32] PRABAGARAN S R, MANORAMA R, DELILLE D, et al. Predominance of Roseobacter, Sulfitobacter, Glaciecola and Psychrobacter in seawater collected off Ushuaia, Argentina, Sub-Antarctica[J]. FEMS Microbiol Ecol, 2010, 59(2): 342-355.
[33] IVANOVA E P, GORSHKOVA N M, SAWABE T, et al. Sulfitobacter delicatus sp. nov. and Sulfitobacter dubius sp. nov. respectively from a starfish (Stellaster equestris) and sea grass (Zostera marina)[J]. Int J System Evolution Microbiol, 2004, 54(2): 475-480. doi: 10.1099/ijs.0.02654-0
[34] PARK J R, BAE J W, NAM Y D, et al. Sulfitobacter litoralis sp. nov., a marine bacterium isolated from the East Sea, Korea[J]. Int J Syst Evol Microbiol, 2007, 57(4): 692-695. doi: 10.1099/ijs.0.64267-0
[35] PALACIOS L, ARAHAL D R, REGUERA B, et al. Hoeflea alexandrii sp. nov., isolated from the toxic dinoflagellate Alexandrium minutum AL1V[J]. Int J Sys Evolution Microbiol, 2006, 56(8): 1991-1995. doi: 10.1099/ijs.0.64238-0
[36] JUNG M Y, SHIN K S, KIM S, et al. Hoeflea halophila sp. nov., a novel bacterium isolated from marine sediment of the East Sea, Korea[J]. Antonie Van Leeuwenhoek Int J Gen Mol Microbiol, 2013, 103(5): 971-978. doi: 10.1007/s10482-013-9876-6
[37] YANG Q, JIANG Z W, HUANG C H, et al. Hoeflea prorocentri sp. nov. isolated from a culture of the marine dinoflagellate Prorocentrum mexicanum PM01[J]. Antonie Van Leeuwenhoek Int J Gen Mol Microbiol, 2018, 111(10): 1-9.
[38] 陈正浩, 张永雨, 杨素萍. 海洋玫瑰杆菌类群研究进展[J]. 生态学报, 2015, 35(5): 1620-1629. [39] GENG H F, BELAS R. Molecular mechanisms underlying roseobacter-phytoplankton symbioses[J]. Curr Opin Biotechnol, 2010, 21(3): 332-338. doi: 10.1016/j.copbio.2010.03.013
[40] CROFT M T, LAWRENCE A D, RAUX-DEERY E, et al. Algae acquire vitamin B-12 through a symbiotic relationship with bacteria[J]. Nature, 2005, 438(764): 90-93.
[41] 杨小茹, 苏建强, 郑小伟, 等. 基于分子技术的1株产毒藻藻际细菌多样性分析[J]. 环境科学, 2009, 30(1): 271-279. doi: 10.3321/j.issn:0250-3301.2009.01.046 [42] 王剑, 王朝晖, 熊毅俊. 不同时期条纹环沟藻可培养藻际细菌研究[J]. 中国环境科学, 2014, 34(6): 1540-1547. -
期刊类型引用(1)
1. 赵玉薇，黄成，邢佳敏，陈兴汉，张虹. 疣吻沙蚕(Tylorrhynchus heterochaetus)成熟精子及卵子的超微结构观察. 福建农林大学学报(自然科学版). 2024(06): 792-796 . 
百度学术
其他类型引用(1)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 2
粤公网安备 44010502001741号
