Spatio-temporal changes of bacterioplankton communities in Litopenaeus vannamei desalinated ponds and their responses to physicochemical factors
-
摘要: 浮游细菌调控是对虾养殖水体环境控制策略的核心内容，探究浮游细菌群落构建的一般规律，可进一步推动对虾养殖水体水质调控技术的研究。运用16S rRNA高通量测序技术，对凡纳滨对虾 (Litopenaeus vannamei) 淡化养殖池塘进行了12次周际调查。结果显示，48个样品共获得2 854个操作分类单元 (Operational taxonomic units, OTU，97%相似性)，序列比对发现古细菌2门1纲1科1属，细菌30门59纲98目199科433属，其中优势菌群25属。优势菌群在组成上有较高的相似性，但各池优势菌群在分布和相对丰度变动上有较大差异。各池系统发育多样性指数总平均值为77.57，变幅为24.39~111.65；香农多样性指数总平均值为3.96，变幅为2.64~5.06；物种丰富度指数总平均值为716，变幅为229~1 054。非度量多维标度分析 (Non-metric multidimensional scaling, NMDS) 表明各池塘浮游细菌群落在养殖初期差异较大，中、后期差异减小，冗余分析 (Redundancy analysis, RDA) 显示活性磷、碱度、溶解氧和硫化物可显著影响浮游细菌的群落结构。Abstract: The regulation of bacterioplankton is the core content of the environmental control strategy in shrimp aquaculture. Exploring the general rules of the construction of bacterioplankton community can further promote the research on water quality regulation for shrimp ponds. Using 16s rRNA high-throughput sequencing technology, we conducted 12 weekly surveys in Litopenaeus vannamei desalinated ponds.The results indicate that a total of 2 854 OTUs (97% similarity) were obtained from 48 samples (Archaea belonged to 2 phyla, 1 class, 1 family and 1 genus; and bacteria belonged to 30 phyla, 59 classes, 98 orders, 199 families and 433 genera, among which 25 genera were dominant flora). The dominant flora had high similarity in the composition, but with great differences in the distribution and relative abundance in each pond. The total average phylogenetic diversity index was 77.57, ranging from 24.39 to 111.65; the total average Shannon diversity index was 3.96, ranging from 2.64 to 5.06; the total average species richness index was 716, ranging from 229 to 1 054. NMDS analysis shows that the community structure of bacterioplankton in each pond varied greatly at the early stage of aquaculture, but not so greatly at the middle and late stages. The results of redundancy analysis shows that the labile phosphorus, total alkalinity, dissolved oxygen and sulfide were the main environmental factors affecting the distribution characteristics of community structure of bacterioplankton.
-
近年来，虾类养殖已成为水产养殖发展较为迅速的产业之一[1]。脊尾白虾 (Exopalaemon carinicauda) 是黄、渤海沿岸主要经济虾类之一，具有广盐、广温性，抗病力和繁殖力较强等优点，已逐渐成为沿海滩涂混合养殖的重要品种；脊尾白虾对低盐耐受性非常强，甚至可以在淡水中生活，是目前低盐混养的重要虾类品种[2-3]。水体是虾类健康养殖的关键因素，水体生态环境的变化会影响虾类的生长和存活。急性低盐胁迫对凡纳滨对虾 (Litopenaeus vannamei)、克氏原螯虾 (Procambarus clarkii)、斑节对虾 (Penaeus monodon)、罗氏沼虾 (Macrobrachium rosenbergii) 和东方新糠虾 (Neomysis orientalis) 等的生长和成活率均有显著影响[4-7]。低盐环境下，凡纳滨对虾和克氏原螯虾的肝胰腺均受到损伤，抗氧化能力和免疫酶活力下降，并使凡纳滨对虾卵巢发育相关基因受到抑制[8-10]。免疫相关基因C型凝集素在斑节对虾肝胰腺和肠中的表达量显著上升，酚氧化酶原基因在脊尾白虾血淋巴和鳃组织中表达量显著升高，参与低盐胁迫的应激反应[11-13]。
转录组测序可在各种环境条件下对物种进行高通量测序，通过测序结果进行基因结构分析和基因功能注释，分析特定条件下相关基因的表达情况，以揭示其代谢网络及调控响应机理[14-15]。转录组测序已在临床诊断、农业基因组学和法医学科学等多个领域，发挥着越来越重要的作用，其测序成本也大幅降低[16]。目前高通量转录组测序技术在鱼类、贝类和虾蟹类中均有应用，主要集中在各种胁迫对水产动物生理机制的影响、环境胁迫相关基因筛选、生长发育相关基因筛选、饲料成分对水产动物转录组的影响、微卫星育种标记开发等[17-21]。王日芳[22]进行脊尾白虾转录组测序，开发了33个微卫星标记进行遗传多样性分析，发现3个近交家系的近交程度较高。孙政[23]利用新一代测序技术 (ROCHE454) 对脊尾白虾胚胎发育不同时期以及成虾应对各种病原刺激时的转录组进行测序，发现成虾在应对各种病原刺激时免疫相关基因表达量增加，在胚胎阶段参与外皮发生的基因十分活跃。冯宁宁[24]对脊尾白虾WSSV潜伏组与急性感染组的转录组进行测序，筛选出濒死与对照群体中的差异表达基因。为探索脊尾白虾低盐度人工养殖的可能性，筛选脊尾白虾在自然海水和淡水环境下的差异表达基因，本研究运用第二代Illumina测序技术对脊尾白虾进行了转录组测序，为进一步研究其在低盐环境下的调控机理提供理论依据。
1. 材料与方法
1.1 实验设计
实验在江苏海洋大学海水养殖动物病害与生态免疫学实验室进行，实验用脊尾白虾来自连云港南极路水产市场，规格整齐、体质健壮，养殖1周以适应实验室环境，暂养水体温度24~26 ℃、盐度31、pH 8.25、溶解氧质量浓度7.6~8.1 mg·L−1。
实验1开始时随机挑选脊尾白虾120尾，平均湿质量 (4.12±0.56) g，分为自然海水 (盐度31) 和低盐胁迫 (盐度0.2) 2个处理组，每组设3个重复，共6个水族箱，每个水族箱 (水体积约30 L) 放养20尾脊尾白虾，人工充氧。胁迫约5.5 h后虾体开始从尾部变白，约6 h后虾开始侧卧游动，估计此时脊尾白虾处于昏厥状态。为确保样品质量，低盐胁迫组脊尾白虾侧卧约5 min开始取样 (低盐胁迫约6 h)，海水组与低盐胁迫组每个水族箱分别选取5尾，整虾速冻于液氮中，备用于转录组测序。实验2设计同上，随机挑选脊尾白虾120尾，分别在第0、第0.5、第1、第2、第4和第6小时，从每个水族箱取出2尾，解剖出肌肉、鳃和肝胰脏于液氮中速冻保存，备用于实时定量验证实验。
1.2 方法
1.2.1 RNA提取与文库构建
将冻存的实验1低盐胁迫组和海水组样品分别于液氮中全虾研磨混合均匀后进行脊尾白虾总RNA提取，检测RNA样品纯度、浓度和完整性。随后用带有Oligo (dT) 的磁珠富集脊尾白虾mRNA，加入破碎液将mRNA进行随机打断，以mRNA为模板，反转录合成双链cDNA链，利用AMPure XP beads纯化cDNA，纯化的双链cDNA再进行末端修复、加尾并连接测序接头，然后用AMPure XP beads进行片段大小选择，最后通过PCR富集构建cDNA文库，检测文库的浓度和插入片段大小，并通过Q-PCR准确定量有效浓度，以保证文库质量。常规方法提取实验2样品RNA，用实时定量反转录试剂盒制成cDNA备用。
1.2.2 转录组测序与分析
实验1获得的检测合格的cDNA文库送样到北京百迈客生物科技有限公司进行HiSeq2500高通量测序，测序读长为PE125。对原始数据进行过滤，滤掉原始数据中的接头序列和低质量读数，获得高质量的数据后，通过Trinity软件进行序列组装得到脊尾白虾的Unigene库。使用BLAST软件将Unigene序列与NR、Swiss-Prot、GO、COG、KOG和KEGG等多个数据库进行比对得到Unigene注释信息。将低盐胁迫组和海水组高质量的测序数据与组装得到的Unigene库进行序列比对，使用EBSeq进行差异表达分析，将错误发现率FDR (False discovery rate)<0.01且差异倍数FC (Fold change)≥2作为筛选标准，最终得到自然海水组和低盐胁迫组样品间的差异表达基因[25-27]，进行差异表达基因的功能注释和富集分析。利用topGO软件对海水组和低盐胁迫组样品间差异表达基因进行GO富集分析，以有向无环图展示各富集节点在GO体系中的关联，挖掘显著富集节点中的差异表达基因，分析差异表达基因可能行使的分子功能。以差异表达基因钠/唾液酸共转体 (Sodium/sialic acid cotransporter, Sialin) 为例，对所在的溶酶体信号通路进行展示并分析。
1.2.3 差异表达基因的筛选与验证
从转录组分析得到的自然海水组和低盐胁迫组样品间差异表达基因中，挑选生殖蛋白A (Reproductive-related protein A, RPA)、ETS转录因子4 (ETS-related transcription factor 4, Elf 4)、钠/唾液酸共转体、胰蛋白酶 (Tripsin)、Smads蛋白4 (Mothers against decapentaplegic protein 4, SMAD 4)、钙结合蛋白 (Calcium-binding protein, CBP) 和α-葡糖苷酶 (Alpha-glucosidase, AGC) 等7条差异表达基因进行引物设计，以实验2中反转录得到的cDNA为模板，通过RT-PCR方法进行验证，表1为验证引物序列，其中Actin为内参基因。利用SPSS 11.0 软件进行统计分析，P<0.05表示差异显著。
表 1 RT-PCR验证引物Table 1. Primers for real-time quantitative verification引物名称
Primer name引物序列 (5'–3')
Primer sequenceActin F CCGAGACATCAAGGAGAAGC Actin R ATACCGCAAGATTCCATACCC RPA RT F GGAAATGCCAACACCTACGC RPA RT R TCCACGTAGCTTTCCTGAGC Elf 4 RT F ACATTGAGGAGCTGGATCGC Elf 4 RT R AGGCGAACTCTGAATGACGG Sialin RT F TTAGTGCATGGAGGCTGGTATT Sialin RT R ATTGACGGCTACGTGTTGAGGT Tripsin RT F CGTCCTCACTTGTTCTCCTTGA Tripsin RT R TGCCGTTACCTAATGTCTGGTT SMAD 4 RT F CTGTCACGGTTGATGGTTACGT SMAD 4 RT R CTTGCTCGGTCACTCTGCTCT CBP RT F CTCAAGAAGGGCACAGGATTT CBP RT R TACCCTGAGTCATGGAACTGC AGC RT F TCACCTTCTGCCATACTTATACTCT AGC RT R TCAGCACCTCGTAACTTCCTC 注：RPA. 生殖蛋白A；Elf 4. ETS转录因子4；Sialin. 钠/唾液酸共转体；SMAD 4. Smads蛋白4；CBP. 钙结合蛋白；AGC. α-葡糖苷酶 Note: RPA. Reproductive-related protein A; Elf 4. ETS-related transcription factor; Sialin. Sodium/sialic acid cotransporter; SMAD 4. Mothers against decapentaplegic protein 4; CBP. Calcium-binding protein; AGC. Alpha-glucosidase 2. 结果
2.1 测序结果及产出质量
运用Illumina HiSeq2500高通量测序共得到13.92 Gb高质量数据，海水组与低盐胁迫组样品Q30碱基百分比不小于87.93% (表2)，质量较高。
表 2 低盐胁迫组与海水组数据产出质量Table 2. Data output quality of freshwater and seawater groups样品
Sample总读数
Read Number总碱基数
Base NumberGC含量
GC Content/%Q30碱基百分比
Q30 base percentage/%低盐胁迫组 (盐度0.2) Low salinity stress group (Salinity 0.2) 25 971 965 6 540 875 181 44.07 88.79 海水组 (盐度31) Seawater group (Salinity 31) 29 309 634 7 381 520 261 44.14 87.93 注：Q30碱基百分比指质量值大于或等于30的碱基所占的百分比，而Q30对应的碱基识别出错的概率为1/1 000 Note: Q30 base percentage refers to the percentage of base with quality score greater than or equal to 30, and the probability of incorrect base call corresponding to Q30 is 1/1 000. 测序数据组装共得到111 618条转录本和72 734条Unigene，转录本与Unigene的N50分别为2 341和1 502。Unigene平均长度为951.16，其中长度在300~500 bp的占比最高 (47.56%)，长度在1 kb以上的占24.93% (表3)。
表 3 组装结果统计Table 3. Assembly results长度范围
Length range/bp转录本
Transcript单基因簇
Unigene300~500 43 156 (38.66%) 34 590 (47.56%) 500~1 000 29 350 (26.30%) 20 011 (27.51%) 1 000~2 000 18 777 (16.82%) 10 181 (14.00%) >2 000 20 335 (18.22%) 7 952 (10.93%) 总数 Total number 111 618 72 734 总长度 Total length 142 324 737 69 181 596 N50长度 N50 length 2 341 1 502 平均长度 Mean length 1 275 11 951.16 2.2 Unigenes功能注释
将Unigene序列与各数据库进行比对，获得Unigene注释信息。脊尾白虾2个样品共有72 734个Unigene，通过筛选最终获得22 879条有注释信息的Unigene，约31.45%，其中有21 931条比对到Nr数据库 (表4)。
表 4 Unigene注释统计表Table 4. Unigene annotation statistics数据库
DatabaseUnigene数量
Number of annotated unigene长度≥300 nt的Unigene数量
Number of Unigene whose length≥300 nt长度≥1 000 nt的Unigene数量
Number of Unigene whose length≥1 000 ntCOG 7 035 7 035 4 381 GO 9 100 9 100 5 046 KEGG 6 645 6 645 4 137 KOG 14 515 14 515 8 732 Pfam 15 403 15 403 9 720 Swiss-Prot 14 940 14 940 9 093 Nr 21 931 21 931 11 913 All 22 879 22 879 12 058 2.2.1 Nr数据库注释
蚤状溞 (Daphnia pulex) 的同源序列最多，2 143条，占注释序列总数的9.77%；其次为赤拟谷盗 (Tribolium castaneum) 1 157条，文昌鱼 (Branchiostoma floridae) 987条，紫色球海胆 (Strongylocentrotus purpuratus) 927条，虱子 (Pediculus humanus) 755条，小头虫 (Capitella teleta) 693条，囊舌虫 (Saccoglossus kowalevskii) 618条，牡蛎 (Crassostrea gigas) 569条，丽蝇蛹集金小蜂 (Nasonia vitripennis) 562条，海兔 (Aplysia californica) 550条和其他物种12 813条。
2.2.2 GO功能分类
GO功能注释结果见图1，共有9 100条Unigenes被注释分类，从细胞组分可细分为15类，占比最多的是细胞组成 (31.50%)，其次为细胞 (31.12%)。从分子功能细分为16类，催化活性类包含Unigenes最多 (49.45%)，结合次之 (45.13%)。依据其参与的生物学过程可细分为22类，代谢过程类包含Unigenes最多 (55.59%)，细胞过程类次之 (49.26%)。
![]()
2.2.3 COG功能分类
通过与COG数据库进行比对，22 879条Unigenes中有7 035条 (30.75%) 比对到了同源序列，共包含9 379条功能注释信息，分布于25个大类 (图2)。其中功能预测类共有2 430条注释信息，占比最大 (25.91%)，其次为复制、重组和修复，808条注释信息 (8.61%)。
2.2.4 Unigene的KEGG代谢途径分析
KEGG分析显示，共有6 645条Unigenes被注释，分布于216个已知途径中，其中前12个代谢途径见表5，注释基因数占总量的31.12%。前5个途径分别是核糖体 (ko03010) 312条、溶酶体 (ko04142) 225条、内质网中的蛋白质加工 (ko04141) 223条、RNA转运 (ko03013) 221条和剪接体 (ko03040) 171条。
表 5 前12个代谢途径Table 5. First 12 metabolic pathways代谢通路
Pathway代谢通路
名称
Pathway ID基因数量
Gene number核糖体 Ribosome ko03010 312 溶酶体 Lysosome ko04142 225 内质网中的蛋白质加工
Protein processing in endoplasmic reticulumko04141 223 RNA转运 RNA transport ko03013 221 剪接 Spliceosome ko03040 171 嘌呤代谢 Purine metabolism ko00230 149 泛素-蛋白酶体
Ubiquitin mediated proteolysisko04120 147 吞噬体 Phagosome ko04145 141 内噬作用 Endocytosis ko04144 140 氧化磷酸化 Oxidative phosphorylation ko00190 125 mRNA监测途径
mRNA surveillance pathwayko03015 108 过氧物酶体 Peroxisome ko04146 106 2.3 差异表达基因功能注释及富集分析
低盐胁迫产生1 492条差异表达基因，有829条显著上调，663条显著下调。共810条差异表达基因得到注释，其中252条注释到GO数据库，165条注释到KEGG数据库，222条注释到COG数据库。
2.3.1 差异表达基因GO功能富集
差异表达基因GO分类注释图，其横坐标为三大分类下的二级节点，纵坐标为该节点中被注释基因的数量和百分比 (图3)。分析得出，差异表达基因主要注释到细胞组成、催化活性和代谢过程中。而百分比明显大于所有基因占比的节点为转运活性和色素沉着，富集显著性较高，但色素沉着的注释百分比较低。
利用topGO对海水组和低盐胁迫组的差异表达基因进行功能富集，得到的有向无环图 (分子功能)，可以清晰地显示出各富集节点的层级关系，分子功能越具体层级越低 (图4)。最显著的10个节点在图中用方框表示。各节点的富集显著性以颜色区分，由高到低为红色>橘色>黄色。
![]()
图 4 差异表达基因GO富集有向无环图 (分子功能)Figure 4. Directed acyclic map of GO enrichment of differentially expressed genes (Molecular function)将最显著的8个节点按富集显著性排序，各节点的详细信息见表6。GO:0003824富集程度最高，该节点中表达量上调幅度最大的Unigene为α-己二酸氨基转移酶，log2FC数值为3.73；下调幅度最大的为磷脂酶，log2FC数值为−2.73 (图4)。
表 6 差异表达基因topGO富集结果示意表 (分子功能)Table 6. Schematic diagram of topGO enrichment results of differential expression gene (Molecular function)节点编号
GO ID节点名称
GO term总基因数
Total gene numberDEG基因数
DEG gene number期望值
Expected显著性
KSGO:0003824 催化活性 Catalytic activity 4 507 125 120.43 3.1×10−9 GO:0042302 表皮结构成分 Structural constituent of cuticle 29 5 0.77 0.000 22 GO:0070011 肽酶活性 Peptidase activity 224 5 5.99 0.002 85 GO:0016787 水解酶活性 Hydrolase activity 1 754 48 46.87 0.003 26 GO:0004879 配体激活序列特异性 DNALigand-activated sequence-specific DNA 12 0 0.32 0.004 06 GO:0005230 细胞外配体门控离子通道 Extracellular ligand-gated ion channel 37 0 0.99 0.004 17 GO:0005319 脂质转运蛋白活性 Lipid transporter activity 23 8 0.61 0.004 38 GO:0022857 跨膜转运蛋白活性 Transmembrane transporter activity 453 18 12.1 0.005 09 注：DEG指差异表达基因 Note: DEG is differentially expressed gene. 2.3.2 差异表达基因KEGG注释
根据KEGG数据库注释结果，分析代谢通路中差异基因表达情况，通过通路探索差异的根源。为了识别差异表达基因参与的代谢途径，对差异表达基因KEGG的注释结果按照通路类型进行分类，发现注释到溶酶体信号通路的基因最多 (13条)，其次为赖氨酸降解信号通路 (7条)，光转导和色氨酸代谢信号通路分别有5条 (图5)。
![]()
图 5 差异表达基因KEGG分类图1. 遗传信息处理；2—6. 有机系统；7—9. 细胞过程；10—15. 环境信息处理；16. 人类疾病；17—44. 代谢Figure 5. KEGG classification of differential expression gene1. Genetic information processing; 2–6. Organismal systems; 7–9. Cellular processes; 10–15. Environmental information processing; 16. Human diseases; 17–44. Metabolism2.3.3 KEGG通路富集分析
通路的富集程度由富集因子体现，值越小富集水平越显著；而Q值 (多重假设检验校正后的P值) 越大富集显著性越可靠。筛选出20条差异表达基因显著富集的代谢通路，其中富集程度最显著的8个代谢通路是光转导、赖氨酸降解、赖氨酸生物合成、溶酶体、血管内皮生长因子信号途径、叶酸生物合成、色氨酸代谢和血管平滑肌收缩 (表7)。
表 7 KEGG富集部分结果Table 7. Partial results of KEGG enrichment代谢通路
Metabolic pathwayKEGG通路编号
KO pathway No.富集因子
Enrichment factorQ 光转导 Phototransduction-fly ko04745 0.13 3.412 2×10−2 赖氨酸降解 Lysine degradation ko00310 0.24 1.095 9×10−1 赖氨酸生物合成 Lysine biosynthesis ko00300 0.10 2.608 2×10−1 溶酶体 Lysosome ko04142 0.43 2.947 3×10−1 血管内皮生长因子信号途径 VEGF signaling pathway ko04370 0.05 3.363 9×10−1 叶酸生物合成 Folate biosynthesis ko00790 0.18 5.565 8×10−1 色氨酸代谢 Tryptophan metabolism ko00380 0.24 5.799 3×10−1 血管平滑肌收缩 Vascular smooth muscle contraction ko04270 0.07 8.143 4×10−1 对溶酶体信号通路进行展示，可以看出溶酶体信号通路中的溶酶体酸性水解酶包括蛋白酶、糖苷酶和鞘磷脂酶等 (图6，表8)。低盐胁迫条件下，组织蛋白酶L显著上调，log2FC数值为2.48，而组织蛋白酶A和B显著下调，log2FC数值分别为−2.46、−1.28和−1.38；相对而言，葡糖脑苷脂酶、鞘磷脂磷酸二酯酶和鞘磷脂磷酸二酯酶2均显著下调，log2FC数值分别为−1.26、−1.89和−1.31。溶酶体信号通路主要溶酶体膜蛋白 (Major lysosomal membrane proteins) 中，溶酶体膜蛋白2在低盐胁迫条件下显著下调，log2FC数值为−1.64；次要溶酶体膜蛋白 (Minor lysosomal membrane proteins) 中，NPC细胞内胆固醇转运体1和2在低盐胁迫条件下均下调，log2FC数值分别为−1.16和−2.11，而钠/唾液酸共转体为混合调节，log2FC数值分别为3.02、1.35和−2.01。
表 8 溶酶体信号通路包含Unigenes序列信息Table 8. Lysosome signaling pathway containing Unigenes sequence information基因描述
Gene description海水组FPKM值
FPKM in seawater group低盐胁迫组FPKM值
FPKM in low salinity stress group差异倍数对数值
log2FC上调/下调
Up/Down钠/唾液酸共转体
Sodium/sialic acid cotransporter, Sialin0.72 6.87 3.02 上调 钠/唾液酸共转体
Sodium/sialic acid cotransporter, Sialin3.70 10.83 1.35 上调 钠/唾液酸共转体
Sodium/sialic acid cotransporter, Sialin4.10 1.14 −2.01 下调 组织蛋白酶A Cathepsin A 0.23 1.63 −2.46 下调 组织蛋白酶B Cathepsin B 0.92 0.16 −1.28 下调 组织蛋白酶B Cathepsin B 66.66 31.23 −1.38 下调 组织蛋白酶L Cathepsin L 60.39 26.51 2.48 上调 鞘磷脂磷酸二酯酶
Sphingomyelin phosphodiesterase, SMPD15.50 1.67 −1.89 下调 鞘磷脂磷酸二酯酶2
Sphingomyelin phosphodiesterase 2, SMPD116.44 7.57 −1.31 下调 NPC细胞内胆固醇转运体1
NPC intracellular cholesterol transporter 1, NPC147.00 38.91 −1.16 下调 NPC细胞内胆固醇转运体2
NPC intracellular cholesterol transporter 2, NPC26.45 13.48 −2.11 下调 溶酶体膜蛋白2
Lysosome membrane protein 2, LIMP1.70 0.60 −1.64 下调 葡糖脑苷脂酶
Glucocerebrosidase, GBA14.69 6.97 −1.26 下调 2.4 差异表达基因的验证
运用RT-PCR测得差异表达基因在脊尾白虾不同部位 (鳃、肝胰腺和肌肉) 中的表达量，并得到差异倍数对数值 (log2FC) 随时间的变化。图7参考值为转录组测序数据中该差异表达基因，根据log2转换后的差异倍数对数值 (log2FC)，详细信息见表9。生殖蛋白A基因在肝胰腺中表达量显著升高，第6小时log2FC为4.23，与转录组分析结果较为一致；但在鳃和肌肉中表达下调 (图7-A)。ETS转录因子4基因表达量在鳃中呈先下降再上升的趋势，第4、第6小时log2FC约为2.58，与转录组分析结果较为一致，而在肝胰腺和肌肉中呈波动变化 (图7-B)。钠/唾液酸共转体基因在鳃、肝胰腺和肌肉中差异表达均不显著 (图7-C)。胰蛋白酶基因在鳃和肌肉中表达下调，与转录组分析结果较为一致；但在肝胰腺中表达上调，呈先升后降的趋势 (图7-D)。Smads蛋白4基因在鳃和肌肉中表达下调，与转录组分析结果较为一致，而在肝胰腺中呈波动变化 (图7-E)。钙结合蛋白基因表达量在鳃、肝胰腺和肌肉3个检测组织下调幅度均较转录组分析结果大 (图7-F)。α-葡糖苷酶基因在鳃中表达上调，而在肝胰腺和肌肉中表达下调 (图7-G)。
![]()
图 7 低盐胁迫对脊尾白虾log2FC的影响log2FC. 差异倍数对数值Figure 7. Effects of low salinity stress on log2FC of E. carinicaudalog2FC. The pair value of fold change表 9 验证的Unigenes序列信息Table 9. Sequence information of verified unigenes基因描述
Gene description海水组FPKM值
FPKM in seawater group低盐胁迫组FPKM值
FPKM in low salinity stress group差异倍数对数值
log2FC上调/下调
Up/Down生殖蛋白A Reproductive-related protein A, RPA 8.01 190.85 4.37 上调 ETS转录因子4 ETS-related transcription factor 4, Elf 4 0.45 4.34 2.97 上调 钠/唾液酸共转体 Sodium/sialic acid cotransporter, Sialin 0.72 6.87 3.02 上调 胰蛋白酶 Tripsin 8.07 0.88 −3.20 下调 Smads蛋白4 Mothers against decapentaplegic protein 4, SMAD 4 1.589 0.56 −1.36 下调 钙结合蛋白 Calcium-binding protein, CBP 0.502 0.451 −0.13 下调 α-葡糖苷酶 Alpha-glucosidase, AGC 3.538 5.228 0.55 上调 3. 讨论
高通量转录组测序技术在水产养殖中已得到广泛应用，刘九美[28]研究开发出脊尾白虾多态性微卫星位点及SNP位点，以揭示其回交家系的优良性状。通过转录组测序及后续分析，利用RACE技术可克隆得出脊尾白虾全长基因[29-30]。王传聪等[31]进行罗氏沼虾肝胰腺组织转录组测序，获得33 450条Unigenes。孙健[32]用Illumina Hiseq 4000测序平台，对氨氮胁迫后日本沼虾 (M. nipponense) 的肝胰腺进行转录组测序，共获得75 742条转录本，63 453条Unigenes，注释率100%。陈雪峰等[33]采用Illumina HiSeqTM4000高通量测序平台，对罗氏沼虾卵巢发育4个时期的卵巢组织进行转录组测序，共获得95 379条Unigenes，注释率54.55%。董丽君等[34]进行凡纳滨对虾肝胰腺组织Illumina HiSeq 2500测序，筛选低温胁迫下的差异表达基因，共获得50 921条Unigenes，平均长度为828 bp，N50为1 589 bp，其中37.28%得到注释。本研究运用Illumina HiSeq2500测序平台，对脊尾白虾进行转录组测序，共获得13.92 Gb高质量数据，海水组与低盐胁迫组样品共获得72 734条Unigenes，平均长度为951.16 bp，N50为1 502 bp，其中22 879条得到注释，占31.45%。注释率较低，可能是由于相近物种的基因注释信息在数据库中收录较少，且脊尾白虾自身转录组测序研究较少。
陈科[35]对凡纳滨对虾低盐度胁迫转录组jin'x测序，发现最显著变化的GO通路为催化活性、组蛋白H4乙酰化和表皮结构成分等；显著变化的KEGG通路有鞘糖脂生物合成、赖氨酸降解和糖胺聚糖生物合成信号通路等。本研究通过差异表达基因GO功能富集分析，发现低盐胁迫对脊尾白虾催化活性、表皮结构成分、肽酶活性、水解酶活性、配体激活序列特异性、脂质转运蛋白活性、跨膜转运蛋白活性等影响显著。KEGG通路富集分析结果显示，主要富集到的8个代谢通路，分别为光转导、赖氨酸降解、赖氨酸生物合成、溶酶体、血管内皮生长因子信号通路、叶酸生物合成、色氨酸代谢、血管平滑肌收缩，其中注释到溶酶体信号通路的差异表达基因最多。
溶酶体是细胞内具有降解、胞外分泌和信号传递功能的一类细胞器，包括溶酶体外膜和溶酶体内腔[36]。细胞受到胁迫时激活自噬功能，可通过溶酶体降解来达到代谢平衡[37]。溶酶体中含有很多蛋白，半胱氨酸蛋白酶是其中含量较多的一种组织蛋白酶，可降解并消化溶酶体内腔中的生物大分子[36]。神经酰胺是细胞凋亡的第二信使，而葡萄糖神经酰胺合成酶催化其糖基化，使之无细胞毒性，促进鞘糖脂合成[38]。对溶酶体信号通路挖掘发现，低盐度胁迫条件下，注释到该信号通路的13个差异表达基因中，钠/唾液酸共转体为混合调节，组织蛋白酶L显著上调，而其他基因均显著下调，推测低盐度胁迫影响了溶酶体对胆固醇的清除和转运，可能通过调控神经酰胺相关酶的活性，影响细胞的生长和信号传导，并在一定程度上影响了溶酶体的降解功能。
ETS转录调控因子家族参与调节造血细胞、神经细胞和血管内皮细胞等的分化，ETS转录因子4是ETS家族重要成员，参与细胞多种生理或病理过程[39-41]。Smads蛋白4是SMADs家族中公共介体性SMAD (Common mediator SMAD) 的一种，具有保持SMAD复合物转录活性的作用，参与转化生长因子β的细胞内信号转导，与组织纤维化有一定关系[42-43]。钙结合蛋白能与钙离子形成复合物，维护细胞内的钙稳态；钙离子通过钙结合蛋白在生物体内发挥作用，是肌肉纤维运动的主要调控信号，影响肌肉纤维的收缩性能[44]。胰蛋白酶是一种丝氨酸蛋白酶，它不仅起到消化酶的作用，还起到活化作用，激活其他蛋白酶类，对胶原纤维结构有一定破坏作用，可能参与了海参自溶[45-46]。正常情况下胰蛋白酶活性被适当抑制，少量胰蛋白酶原转化为活性胰蛋白酶，被胰蛋白酶抑制剂灭活，防止细胞受损。胰蛋白酶的激活超过胰蛋白酶抑制剂的能力，会导致损伤细胞的各种蛋白酶的激活，这是胰腺炎发病的主要原因[47]。α-葡糖苷酶参与淀粉等糖原的代谢，具有水解和转糖苷的双重作用，抑制其活性可减缓葡萄糖的生成，同时调节体内的脂肪水平[48]。本研究转录组数据分析发现，低盐胁迫显著提高α-葡糖苷酶和ETS转录因子4基因的表达，而Smads蛋白4、钙结合蛋白和胰蛋白酶表达显著降低，推测脊尾白虾可能通过促进ETS转录因子4产生，提高α-葡糖苷酶活性以维持机体的正常生理功能，抑制胰蛋白酶活性以维持脊尾白虾低盐应激反应的稳态；而钙结合蛋白的表达下调表明低盐胁迫影响了脊尾白虾肌肉的纤维收缩能力。实时定量PCR结果显示，基因在鳃、肝胰腺和肌肉中的表达情况不尽相同，并且随时间变化，基因的表达量也有所波动，可能与转录组测序的时间点只有一个临界时间点，而且是全虾整体混合测序，而实时定量验证所用模板为脊尾白虾的各个组织，且覆盖了低盐胁迫过程中的各个阶段有关。
然而，图7验证试验中生殖蛋白A、ETS转录因子4、钠/唾液酸共转体、胰蛋白酶、Smads蛋白4、钙结合蛋白和α-葡糖苷酶等7条差异表达基因偏差较大，可能有2个原因：1) 脊尾白虾低盐胁迫实验结果可能受采样操作过程和自身盐度耐受差异等因素的影响，造成一定偏差；2) 可能与使用的分析软件有关，本文所用的分析软件优点是通过输入RT-PCR测定数值，可以直接算出log2转换后的差异倍数对数值 (log2FC)，基因上调和下调趋势一目了然，弥补了先前使用2−△△CT方法对下调基因展示的不足，但也带来了标准偏差较大的问题。
本研究挖掘了脊尾白虾低盐环境下的相关信号通路及差异表达基因数据，为脊尾白虾分子生物学研究提供了丰富的基因资源，也可为今后进行脊尾白虾低盐胁迫生理机制的探讨和关键差异表达基因的深入挖掘提供技术支撑。
-
![]()
图 1 凡纳滨对虾养殖池塘门水平的浮游细菌组成
Figure 1. Composition of bacteriaoplankon at phylum level in L. vannamei ponds
![]()
图 3 3口池塘的α多样性系数分析
注：箱体上中下线分别为75、50 (中位数) 和25分位数，轴须线最长不超过1.5倍箱体范围，黑色空心圆表示平均数；差异显著性用* (P<0.05)、** (P<0.01)以及*** (P<0.001) 表示；图中的样本量：A：n=22、B：n=12、C：n=12。
Figure 3. α diversity index analysis of bacterioplankton in three ponds
Note: The upper, middle and lower lines of the box are 75, 50 (Median) and 25 quantiles, respectively. The maximum length of whiskers shall not exceed 1.5 times of the box range. The black hollow circles represent the average values. The significant differences were represented by * (P<0.05), ** (P<0.01) and *** (P<0.001). The numbers of replicated samples in this figure are: A: n=22; B: n=12; C: n=12.
![]()
图 4 凡纳滨对虾养殖池塘浮游细菌群落NMDS分析
Figure 4. NMDS analysis of bacterioplankton community in L. vannamei pond
![]()
图 5 各池浮游细菌属水平上共有和特有属的数量
Figure 5. Number of common and endemic genera of bacterioplankton at genus level in each pond
![]()
图 6 各池浮游细菌群落的线性判别分析 (菌群LDA>4)
Figure 6. Linear discriminant analysis of bacterioplankton community in each pond (Bacterial flora with LDA>4)
![]()
图 7 优势菌群分布和群落多样性 (H') 与主要理化因子之间的关系 (II型标尺)
注：Aci. 不动杆菌属；Aer. 气单胞菌属；Fla. 黄杆菌属；Gem. 芽殖杆菌属；GpI. GpIIa；Ilu. 微酸菌属；Lim. 湖栖菌属；Pol. 多核杆菌；Rho. 红杆菌属；Sed. 沉积物杆状菌属；Spa. 发光细菌属；Sph. 鞘脂菌属。
Figure 7. Species associations of dominan flora and diversity (H') with environmental factors (Scaling II)
Note: Aci. Acinetobacter sp.; Aer. Aeromonas sp.; Fla. Flavobacterium sp.; Gem. Gemmobacter sp.; GpI. GpIIa; Ilu. Ilumatobacter sp.; Lim. Limnohabitans sp.; Pol. Polynucleobacter sp.; Rho. Rhodobacter sp.; Sed. Sediminibacterium sp.; Spa. Spartobacteria genera incertae sedis; Sph. Sphingomonas sp..
表 1 凡纳滨对虾养殖池塘水体主要理化因子
Table 1 Environmental factors in L. vannamei ponds
环境因子
Environmental factor池塘A
Pond A池塘B
Pond B池塘C
Pond C水温 Temperature/℃ 26.19±2.01 26.15±1.57 26.37±2.00 pH 8.10±0.50 8.29±0.37 8.30±0.38 溶解氧质量浓度 DO/(mg·L−1) 8.00±0.68 8.45±1.17 8.11±1.15 铵态氮质量浓度 NH4-N/(mg·L−1) 0.45±0.73 0.44±0.39 0.33±0.27 亚硝酸氮质量浓度 NO2-N/(mg·L−1) 0.01±0.01 0.01±0.01 0.01±0.01 硝酸态氮质量浓度 NO3-N/(mg·L−1) 0.34±0.44 0.41±0.42 0.48±0.53 活性磷质量浓度 PO4-P/(mg·L−1) 0.12±0.14 0.37±0.24 0.22±0.28 活性硅酸盐质量浓度 SiO3-Si/(mg·L−1) 11.2±3.55 12.38±7.38 4.03±2.98 高锰酸盐指数 CODMn/(mg·L−1) 8.34±7.19 6.12±3.18 5.34±2.91 叶绿素 a 质量浓度 Chl-a/(mg·L−1) 75.78±131.01 56.13±42.81 58.36±61.48 硫化物质量浓度 Sul/(mg·L−1) 0.03±0.04 0.02±0.01 0.02±0.01 矿化度质量浓度 Mineralization degree/(mg·L−1) 1209.44±2445.58 813.39±171.69 814.35±41.91 碱度 ALK/(mg·L−1) 138.18±12.06 115.32±33.82 84.46±14.73 总硬度 Total hardness/(mg·L−1) 362.76±395.28 476.5±227.34 349.36±60.18 
下载: 导出CSV
表 2 各样本有效序列数据统计
Table 2 Valid sequences of each sample
样品
Sample条形码
Barcode有效序列
Valid sequence/条碱基数
Base number平均长度
Mean length/bp最短序列长度
Min. length/bp最长序列长度
Max. length/bpS1 GTAACA 82 736 34 604 966 418.26 363 469 S2 CCAGAC 81 905 34 619 470 422.68 361 476 S3 GGTGAA 56 859 23 683 959 416.54 367 475 S4 TGCATC 85 212 35 667 198 418.57 353 469 S5 TCGACC 83 246 34 547 936 415.01 365 470 S6 GTCGCG 73 829 30 547 465 413.76 362 453 S7 CGGATG 83 662 34 668 763 414.39 356 465 S8 GTGAAA 84 384 34 811 366 412.54 350 466 S9 ATCTTG 100 945 41 590 813 412.01 350 452 S10 TATGCA 73 801 30 659 241 415.43 352 459 S11 GTAACA 86 471 35 924 249 415.45 354 471 S12 GCGAGG 90 432 37 722 526 417.14 351 465 S13 CACGAT 54 439 22 735 194 417.63 373 465 S14 GCGGTA 44 227 18 478 585 417.81 352 471 S15 TATCGA 61 431 25 638 627 417.36 359 473 S16 ATCACG 52 384 21 792 829 416.02 376 471 S17 CGGATG 99 211 41 215 038 415.43 360 473 S18 CGCATA 100 001 41 427 908 414.27 372 476 S19 TGCATC 65 712 27 246 907 414.64 359 476 S20 TCAGTA 76 029 31 690 841 416.83 362 471 S21 CGGCAC 75 202 31 405 455 417.61 364 469 S22 ATCACG 71 605 30 328 221 423.55 352 470 S23 CGGATG 63 358 26 635 233 420.39 367 470 S24 GTGAAA 56 429 23 698 925 419.98 350 471 S25 TCAGTA 92 400 38 487 742 416.53 355 435 S26 GAAGTG 87 305 37 043 958 424.31 368 448 S27 TCGACC 94 075 39 472 386 419.58 360 464 S28 CTTGTA 52 273 22 054 548 421.91 373 472 S29 GTTTCG 44 697 18 684 447 418.02 360 462 S30 ATCTTG 59 737 25 072 871 419.72 365 465 S32 GCCATC 78 462 32 993 239 420.5 357 474 S33 TGTGTT 69 613 29 167 734 419.00 351 474 S34 CTTGTA 56 083 23 551 548 419.94 356 471 S35 GTTTCG 45 701 18 993 334 415.60 359 464 S36 TTCGTA 46 337 19 212 573 414.63 356 468 S37 CCAGAC 51 978 21 598 326 415.53 372 436 S38 AGCAGT 77 233 31 835 991 412.21 350 470 S39 GAGGAA 75 226 31 060 291 412.89 370 468 S40 AAGGTA 46 920 19 633 227 418.44 352 450 S41 ATCACG 43 983 18 242 559 414.76 352 469 S42 TAGGAC 66 639 27 702 230 415.71 356 470 S43 TGGACG 49 206 20 369 621 413.97 357 472 S44 AGAACA 50 322 20 765 128 412.65 356 469 S45 GGTGTG 41 169 17 011 770 413.22 350 471 S46 AACTAT 67 691 28 070 176 414.68 357 468 S47 ACTGCG 60 220 25 559 117 424.43 359 474 S48 TGTGTT 94 261 39 467 506 418.7 355 474 S49 TAGGAC 86 904 36 530 365 420.35 353 467
下载: 导出CSV
-
[1] CABELLO F C. Heavy use of prophylactic antibiotics in aquaculture: a growing problem for human and animal healthy and for the environment[J]. Environ Microbiol, 2006, 8: 1137-1144. doi: 10.1111/j.1462-2920.2006.01054.x
[2] MA Z, SONG X, WAN R, et al. A modified water quality index for intensive shrimp ponds of Litopenaeus vannamei[J]. Ecol Indic, 2013, 24: 287-293. doi: 10.1016/j.ecolind.2012.06.024
[3] GONZALEZ J M, SIMO R, MASSANA R, et al. Bacterial community structure associated with a dimethylsulfo niopropionate-producing North Atlantic algal bloom[J]. Appl Environ Microb, 2000, 66: 4237-4246. doi: 10.1128/AEM.66.10.4237-4246.2000
[4] WEI C L, ZENG Y H, TANG K X, et al. Comparison of bacterioplankton communities in three mariculture ponds farming different commercial animals in subtropical Chinese coast[J]. Hydrobiologia, 2009, 632: 107-126. doi: 10.1007/s10750-009-9831-8
[5] CARON D A. Inorganic nutrients, bacteria, and the microbial loop[J]. Microb Ecol, 1994, 28: 295-298. doi: 10.1007/BF00166820
[6] COTNER J B, BIDDANDA B. Small players, large role: microbial influence on biogeochemical processes in pelagic aquatic ecosystems[J]. Ecosystems, 2002, 5: 105-121. doi: 10.1007/s10021-001-0059-3
[7] LEMONNIER H, COURTIES C, MUGNIER C, et al. Nutrient and microbial dynamics in eutrophying shrimp ponds affected or unaffected by vibriosis[J]. Mar Pollut Bull 2010, 60: 402-411.
[8] AMANN R I. Phylogenetic identification and in situ detection of individual microbial cells without cultivation[J]. Microbiol Rev, 1995, 59: 143-169. doi: 10.1128/mr.59.1.143-169.1995
[9] LI C, LIU J, CHEN X, et al. Determinism governs the succession of disturbed bacterioplankton communities in a coastal maricultural ecosystem[J]. Sci Total Environ, 2022, 828: 154457. doi: 10.1016/j.scitotenv.2022.154457
[10] LIU C H, CHANG F H, YANG J W, et al. Free-living marine bacterioplankton composition and diversity along the Kuroshio region[J]. Deep-Sea Res I, 2022, 183: 103741. doi: 10.1016/j.dsr.2022.103741
[11] LU Z B, REN Z M, LIN W H, et al. Succession, sources, and assembly of bacterial community in the developing crab larval microbiome[J]. Aquaculture, 2021, 548: 737600.
[12] LIU Y P, PAN B Z, ZHU X Z, et al. Patterns of microbial communities and their relationships with water quality in a large-scale water transfer system[J]. J Environ Manage, 2022, 319: 115678. doi: 10.1016/j.jenvman.2022.115678
[13] 靳燕, 邱莹, 董志, 等. 北运河浮游细菌集合群落空间变化的环境解释[J]. 中国环境科学, 2021, 41(3): 1378-1386. doi: 10.3969/j.issn.1000-6923.2021.03.041 [14] 杨淑芳. 凡纳滨对虾养殖池塘水质动态变化及细菌多样性研究[D]. 厦门: 集美大学, 2017: 10-19. [15] 张皓. 养殖环境微生物群落的动态变化及与水环境的互作[D]. 苏州: 苏州大学, 2015: 16-29. [16] YANG W, ZHENG C, ZHENG Z M, et al. Nutrient enrichment during shrimp cultivation alters bacterioplankton assemblies and destroys community stability[J]. Ecotox Environ Safe, 2018, 156: 366-374. doi: 10.1016/j.ecoenv.2018.03.043
[17] YANG W Z, ZHENG J Y, CHEN Q, et al. Succession of bacterioplankton community in intensive shrimp (Litopenaeus vannamei) aquaculture systems[J]. Aquaculture, 2018, 497: 200-213. doi: 10.1016/j.aquaculture.2018.07.053
[18] 国家环境保护总局水和废水检测分析方法编委会. 水和废水检测分析方法(第四版)[M]. 北京: 中国环境科学出版社, 2002: 30-152. [19] EDGAR R C. Search and clustering orders of magnitude faster than BLAST[J]. Bioinformatics, 2010, 26(19): 2460-2461. doi: 10.1093/bioinformatics/btq461
[20] AKSNES D L, WASSMANN P. Modeling the significance of zooplankton grazing for export production[J]. Limnol Oceanogr, 1993, 38(5): 978-985. doi: 10.4319/lo.1993.38.5.0978
[21] FAITH D P. Conservation evaluation and phylogenetic diversity[J]. Biol Conserv, 1992, 61: 1-10. doi: 10.1016/0006-3207(92)91201-3
[22] SHANNON C E, WEAVER W. The mathematical theory of communication[M]. Urbana, IL: University of Illinois Press, 1949: 55-187.
[23] 杨家新. 微生物生态学[M]. 北京: 化学工业出版社, 2004: 45-134. [24] HAN J H, BAEK K, LEE M H. Lacibacter nakdongensis sp. nov., isolated from river sediment[J]. Int J Syst Evol Microbiol, 2017, 67: 352-356. doi: 10.1099/ijsem.0.001631
[25] WITCZAK P, CYDZIK-KWIATKOWSKA A, ZIELIŃSKA M. Treatment of the liquid phase of digestate from a biogas plant for water reuse[J]. Bio Tech, 2019, 276: 226-235. doi: 10.1016/j.biortech.2018.12.077
[26] XUAN L, SHENG Z, LU J, et al. Bacterioplankton community responses and the potential ecological thresholds along disturbance gradients[J]. Sci Total Environ, 2019, 696: 134015. doi: 10.1016/j.scitotenv.2019.134015
[27] YANG Y, GAO Y, HUANG X, et al. Adaptive shifts of bacterioplankton communities in response to nitrogen enrichment in a highly polluted river[J]. Environ Pollut, 2019, 245: 290-299. doi: 10.1016/j.envpol.2018.11.002
[28] COMTE J, del GIORGIO P A. Composition influences the pathway but not the outcome of the metabolic response of bacterioplankton to resource shifts[J]. PLoS One, 2011, 6: e25266. doi: 10.1371/journal.pone.0025266
[29] ADAMOVICH B V, MEDVINSKY A B, NIKITINA L V, et al. Relations between variations in the lake bacterioplankton abundance and the lake trophic state: evidence from the 20-year monitoring[J]. Ecol Indic, 2019, 7: 120-129.
[30] 吕锡斌, 吴云成, 陈良强, 等. 赤水河流域浮游细菌群落特征及其与水质的关系[J]. 环境科学学报, 2021, 41(11): 4596-4605. doi: 10.13671/j.hjkxxb.2021.0125 [31] JING X Y, GOU H, GONG Y, et al. Seasonal dynamics of the coastal bacterioplankton at intensive fish-farming areas of the Yellow Sea, China revealed by high-throughput sequencing[J]. Mar Pollut Bull, 2019, 139: 366-375. doi: 10.1016/j.marpolbul.2018.12.052
[32] SHADE A, HANNES P, STEVEN D A, et al. Fundamentals of microbial community resistance and resilience[J]. Front Microbiol, 2012, 19(3): 417.
[33] 郭远涛. 吸附固定化对藻菌共生体在畜禽养殖废水处理中的强化效果研究[D]. 南昌: 南昌大学, 2021: 12-36. [34] CAVENDER-BARES J, KOZAK K H, FINE P V A, et al. The merging of community ecology and phylogenetic biology[J]. Ecol Lett, 2009, 12: 693-715. doi: 10.1111/j.1461-0248.2009.01314.x
[35] CHESSON P. Mechanisms of maintenance of species diversity[J]. Annu Rev Ecol Syst, 2000, 31: 343-366. doi: 10.1146/annurev.ecolsys.31.1.343
[36] FAN L M, HU G D, QIU L P, et al. Variations in bacterioplankton communities in aquaculture ponds and the influencing factors during the peak period of culture[J]. Environ Pollut, 2020, 258: 113656. doi: 10.1016/j.envpol.2019.113656
[37] IVES A R, CARPENTER S R. Stability and diversity of ecosystems[J]. Science, 2007, 317: 58-62. doi: 10.1126/science.1133258
[38] YACHI S, LOREAU M. Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis[J]. Proc Natl Acad Sci USA, 1999, 96: 1463-1468. doi: 10.1073/pnas.96.4.1463
[39] 胡晓娟, 李卓佳, 曹煜成, 等. 强天气干扰条件下粤西凡纳滨对虾养殖池塘细菌群落动态特征[J]. 南方水产科学, 2012, 8(5): 52-59. doi: 10.3969/j.issn.2095-0780.2012.05.008 [40] CONNELL J H. Diversity in tropical rain forests and coral reefs: high diversity of trees and corals is maintained only in a non-equilibrium state[J]. Science, 1978, 199: 1302-1310. doi: 10.1126/science.199.4335.1302
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
粤公网安备 44010502001741号
