云南高原異龍湖噬藻體psbA基因遺傳多樣性研究

李春筱，劉婷婷，劉玉珊，劉 麗

(昆明理工大學(xué)生命科學(xué)與技術(shù)學(xué)院，云南 昆明 650500)

噬藻體(Cyanophage)是一類(lèi)能夠特異性感染藍(lán)藻的雙鏈DNA病毒類(lèi)群，廣泛存在于海洋、淡水等自然環(huán)境中，是水生微生物群落的重要組成部分，能在不同時(shí)間和空間范圍內(nèi)感染藍(lán)藻，具有豐富的遺傳多樣性[1]。作為水生態(tài)系統(tǒng)中的重要活躍因子，噬藻體與宿主藍(lán)藻的新陳代謝和生命循環(huán)密切相關(guān)，具有重要的生態(tài)地位；在調(diào)控藍(lán)藻種群結(jié)構(gòu)及多樣性、水體初級(jí)生產(chǎn)力、生物地球化學(xué)循環(huán)以及介導(dǎo)微生物之間水平基因轉(zhuǎn)移等方面起到了重要的作用[2-4]。

噬藻體是高度差異性的病毒類(lèi)群，尚未發(fā)現(xiàn)其基因中存在類(lèi)似原核生物16SrDNA或真核生物18SrDNA的共有保守序列，因而缺乏廣泛通用的遺傳標(biāo)記基因。基于噬藻體的分離純化及全基因組測(cè)序研究發(fā)現(xiàn)部分種的噬藻體中存在個(gè)別共有保守基因，如g20基因[5-6]、psbA光合基因[7]和pol基因[8]等。g20基因編碼肌尾噬藻體衣殼組裝蛋白主要存在于肌尾科噬藻體中。pol基因編碼DNA聚合酶主要存在于短尾科噬藻體。psbA基因是編碼光合系統(tǒng)Ⅱ反應(yīng)中心D1蛋白的基因，在肌尾科和短尾科噬藻體均有發(fā)現(xiàn)[9]。許多噬藻體中都包含編碼光合系統(tǒng)Ⅱ中D1蛋白的psbA基因，88%的感染海洋聚球藻和原綠球藻的噬藻體都含有psbA基因[10]，全球海洋表層水中60%的psbA基因來(lái)源于病毒[11]。SULLIVAN等[12]根據(jù)光合基因的保守序列設(shè)計(jì)引物，對(duì)噬藻體的遺傳多樣性進(jìn)行了研究，發(fā)現(xiàn)噬藻體的核酸序列具有相當(dāng)豐富的多樣性。作為噬藻體的一種新型靶標(biāo)基因，光合作用基因psbA不僅可用于揭示噬藻體的遺傳多樣性，還可利用它探索噬藻體與其宿主藍(lán)藻之間的關(guān)系[2]。

目前，國(guó)內(nèi)對(duì)噬藻體遺傳多樣性的研究主要集中在中國(guó)近海海域[13]、東北稻田水體[14]、武漢東湖[15]等，有關(guān)淡水湖泊噬藻體多樣性的相關(guān)報(bào)道仍十分有限。該研究主要是運(yùn)用以PCR技術(shù)為核心的分子生物學(xué)方法對(duì)云南高原富營(yíng)養(yǎng)化湖泊異龍湖春季和秋季水樣噬藻體psbA基因的多樣性進(jìn)行研究，分析異龍湖水樣中噬藻體的遺傳多樣性，為噬藻體基因多樣性研究提供依據(jù)，并豐富富營(yíng)養(yǎng)化湖泊微型生態(tài)系統(tǒng)的理論體系。

1 材料與方法

1.1 水樣的采集及預(yù)處理

異龍湖水樣采集于2015年4月(春季)與2015年9月(秋季)。參照HJ 495—2009《水質(zhì) 采樣方案設(shè)計(jì)技術(shù)規(guī)定》，使用采水器采集表層水樣(距水體表面約0.5 m處)和底層水樣(距水體表面約2 m處)各1.5 L，現(xiàn)場(chǎng)測(cè)定水溫、pH值和透明度。采集水樣后立即帶回實(shí)驗(yàn)室，將水樣混合后經(jīng)0.75 μm孔徑濾紙、0.45和0.22 μm孔徑微孔濾膜〔上海市新亞凈化器件廠(chǎng)微孔濾膜(滬Q/YY8-1-88)〕過(guò)濾，將過(guò)濾后的水樣用6萬(wàn)分子量的中空纖維膜組件進(jìn)行超濃縮，濃縮比例為200∶1，濃縮液分裝后于-20 ℃ 冰箱避光保存待用。

1.2 水樣中噬藻體基因組提取及PCR擴(kuò)增

1.2.1水樣中游離藍(lán)藻DNA的去除

為了去除濃縮水樣中游離的藍(lán)藻DNA片段，采用DNaseI處理濃縮水樣，并進(jìn)一步用藍(lán)藻16S-23SrDNA ITS基因[16-18]進(jìn)行PCR擴(kuò)增，若為陰性結(jié)果，則去除完全，陽(yáng)性則要使用DNaseI進(jìn)行再處理。

ITS基因由引物[19]ITS-F(5′-TGTACACACCGCCCGTCACACCA-3′)和ITS-R(5′-GTGGATAC-CTAGGCACACAGAG-3′)進(jìn)行PCR擴(kuò)增。PCR擴(kuò)增體系：2×Premix Taq (TAKARA)10 μL，F(xiàn)-Primer(10 μM) 1 μL，R-Primer (10 μM) 1 μL，模板DNA 2 μL，ddH2O 6 μL，總體積20 μL。陰性對(duì)照無(wú)菌ddH2O為模板。PCR擴(kuò)增條件：94 ℃預(yù)變性5 min；94 ℃變性30 s；55 ℃退火30 s；72 ℃延伸90 s；35個(gè)循環(huán)；72 ℃延伸7 min。

1.2.2水樣中病毒基因組提取及psbA基因的PCR擴(kuò)增

取經(jīng)處理的水樣，使用病毒基因組提取試劑盒Viral DNA Kit(OMEGA)提取噬藻體DNA。從GenBank中下載噬藻體psbA基因序列進(jìn)行BALST比對(duì)，采用Primer 5.0和Oligo 6分析評(píng)價(jià)，針對(duì)psbA基因設(shè)計(jì)巢式引物：psbA-F1(5′-GARTGGYTNTAYAAYGGNG-3′)和psbA-R1(5′-CCRTTNAGGTTRAANGCCA-3′)，psbA-F2(5′-CAYTTCTAYCCNATYTGG-3′)和psbA-R2(5′-TCYN GACTGGTTGAAGTTGA-3′)進(jìn)行PCR擴(kuò)增。PCR擴(kuò)增體系：2×Premix Taq (TaKaRa) 10 μL，F(xiàn)-Primer(10 μM) 1 μL，R-Primer (10 μM) 1 μL，模板DNA 2 μL，ddH2O 6 μL，總體積20 μL。陰性對(duì)照無(wú)菌ddH2O為模板。PCR擴(kuò)增條件：94 ℃預(yù)變性5 min；94 ℃變性30 s；退火30 s(psbA147 ℃，psbA250 ℃)；72 ℃ 延伸90 s；35個(gè)循環(huán)；72 ℃延伸7 min。

1.3 膠回收、TA克隆及序列測(cè)定

從w=2.0%的凝膠上切取580 bp左右的PCR產(chǎn)物，用普通瓊脂糖凝膠DNA回收試劑盒(TIANGEN)進(jìn)行膠純化，目的DNA片段連接至pEASY-T1 Cloning Vector，轉(zhuǎn)化大腸桿菌JM109感受態(tài)細(xì)胞。無(wú)菌條件下挑取Amp+的LB平板上50個(gè)樣品克隆，再用通用引物M13F、M13R對(duì)菌液進(jìn)行陽(yáng)性重組子的菌落PCR快速鑒定，退火溫度55 ℃，循環(huán)數(shù)30。選陽(yáng)性克隆的培養(yǎng)過(guò)夜菌液，送至測(cè)序公司(昆明碩擎生物科技)進(jìn)行一代測(cè)序。

1.4 psbA基因的系統(tǒng)進(jìn)化分析

待測(cè)序完成后，對(duì)異龍湖水域水樣中psbA環(huán)境序列的單克隆序列分別用Chromas、BioEdit(Ver 7.0.0)進(jìn)行序列校正、編輯。使用BLAST將得到的序列與數(shù)據(jù)庫(kù)中收錄的噬藻體psbA基因進(jìn)行初步的比對(duì)分型，用MEGA 6.0軟件繪制系統(tǒng)進(jìn)化樹(shù)，根據(jù)系統(tǒng)進(jìn)化樹(shù)來(lái)分析異龍湖中噬藻體的遺傳多樣性。噬藻體序列標(biāo)注為“YL年月S+第幾條序列”。該研究獲得的psbA序列已提交至NCBI網(wǎng)站，其序列號(hào)為MH822090-MH822104。

2 結(jié)果與分析

2.1 異龍湖水樣噬藻體psbA基因的PCR擴(kuò)增及克隆測(cè)序

將采集的異龍湖水樣經(jīng)0.75 μm孔徑濾紙，0.45和0.22 μm微孔濾膜過(guò)濾和超濃縮后，DNase處理去除濃縮水樣中藍(lán)藻DNA，并用藍(lán)藻16S-23SrDNA ITS基因進(jìn)行PCR擴(kuò)增，從結(jié)果為陰性的樣品中提取噬藻體DNA；PCR擴(kuò)增psbA基因于580 bp位置有明顯亮帶(圖1)，經(jīng)膠回收獲取目的基因，樣品經(jīng)克隆篩選，將陽(yáng)性克隆送測(cè)序公司測(cè)序，共獲得15條不同的psbA基因，其中4月(春季)水樣中獲得3條，9月(秋季)水樣中獲得12條psbA環(huán)境序列。

圖1 異龍湖水樣中噬藻體psbA基因的PCR擴(kuò)增及陽(yáng)性克隆檢測(cè)Fig.1 PCR amplification and positive clone of psbA gene from cyanophage in Yilong Lake

2.2 異龍湖水樣中psbA基因親緣序列分析

異龍湖水樣中擴(kuò)增出的psbA基因長(zhǎng)度為580 bp(引物除外)，可翻譯成186個(gè)氨基酸殘基。對(duì)獲得的psbA基因序列進(jìn)行BLAST比對(duì)后，剔除部分非噬藻體源的psbA基因序列。該研究從異龍湖水樣中獲得15條psbA基因環(huán)境序列，經(jīng)NCBI網(wǎng)站上基于氨基酸水平進(jìn)行BLAST比對(duì)發(fā)現(xiàn)(表1)：有4條psbA基因序列(YL04S-1～3和YL09S-6)與哈薩克斯坦淡水中噬藻體psbA序列(KZ5-52073/3)相似性最高，相似率為99%；有10條psbA基因序列(YL09S-1～5，YL09S-7～8，YL09S-10～12)與中國(guó)東北稻田水體克隆(PW-DA-II-21 ，PW-MDJ-II-4，PW-DA-II-35)具有最高的相似性，相似率在98%～99%之間；僅有1條psbA基因序列(YL09S-9)與日本稻田克隆(NoF6)具有最高的相似性，相似率為98%。

表1 異龍湖水體psbA基因與NCBI數(shù)據(jù)庫(kù)在氨基酸水平的同源性比較

Table 1 Comparison of homology at amino acid level betweenpsbAgene and NCBI database in Yilong Lake

克隆名稱(chēng)序列比對(duì)最近親緣序列序列編號(hào)同源性/%序列來(lái)源序列提交者YL04S-1185/186KZ5-52073/3AQT1911399哈薩克斯坦淡水ALEXYUK等,2016YL04S-2185/186KZ5-52073/3AQT1911399哈薩克斯坦淡水ALEXYUK等,2016YL04S-3185/186KZ5-52073/3AQT1911399哈薩克斯坦淡水ALEXYUK等,2016YL09S-1184/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-2185/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-3185/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-4185/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-5185/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-6185/186KZ5-52073/3AQT1911399哈薩克斯坦淡水ALEXYUK等,2016YL09S-7185/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-8184/186PW-DA-II-21AJO7061299中國(guó)稻田水WANG等,2015YL09S-9182/186NoF6BAH2318698日本稻田水WANG等,2009YL09S-10183/186PW-MDJ-II-4AJO7065198中國(guó)稻田水WANG等,2015YL09S-11183/186PW-DA-II-35AJO7062698中國(guó)稻田水WANG等,2015YL09S-12184/186PW-DA-II-21AJO7061298中國(guó)稻田水WANG等,2015

2.3 異龍湖水樣psbA基因系統(tǒng)進(jìn)化樹(shù)分析

psbA基因存在光合作用生物中，因此區(qū)分psbA基因是否為病毒來(lái)源至關(guān)重要。SHARON等[20]通過(guò)統(tǒng)計(jì)分析發(fā)現(xiàn)，psbA基因編碼D1蛋白的氨基酸序列模型R/KETTXXXSQ/H中可變的三肽元件XXX的變化可以反映出不同噬藻體所處的分類(lèi)地位，其中EQE、ENE、EEV、EEE、EDV、EVE、EQV和EDE型等幾種三肽元件更易來(lái)自于病毒D1亞基，而在宿主藍(lán)藻中極少發(fā)現(xiàn)。SULLIVAN等[12]發(fā)現(xiàn)，聚球藻噬藻體的基因組序列中鳥(niǎo)嘌呤和胞嘧啶所占的比例(GC含量)比聚球藻的GC含量少10%左右，即噬藻體GC含量低于宿主藻的GC含量。

從異龍湖濃縮水樣中獲得的15條psbA基因環(huán)境系列，構(gòu)建進(jìn)化樹(shù)(圖2)，可將psbA基因序列劃分成3簇，分別命名為ClusterⅠ～Cluster Ⅲ，對(duì)psbA基因序列的GC含量和編碼D1蛋白的氨基酸序列中三肽元件R/KETTXXXSQ/H的變化進(jìn)行分析發(fā)現(xiàn)，ClusterⅠ中GC含量為47.9%～50.0%，序列中三肽元件全部為EEV型，ClusterⅠ中包括的10條序列為病毒來(lái)源，即來(lái)源于噬藻體；Cluster Ⅱ中只含1條序列，GC含量為48.8%，序列中三肽元件為ETE型；Cluster Ⅲ中GC含量為44.5%～44.6%，4條序列中三肽元件全部為ETE型。研究發(fā)現(xiàn)，在聚球藻菌株及感染聚球藻的肌尾和短尾科噬藻體中均存在三肽元件ETE型[8]，結(jié)合水樣處理方法、GC含量及序列的系統(tǒng)進(jìn)化位置及氨基酸同源性比較，推測(cè)ClusterⅡ和Cluster Ⅲ的psbA基因來(lái)自于噬藻體。因此，該研究獲得的15條psbA基因序列均來(lái)自于噬藻體。

圖2 異龍湖噬藻體psbA序列進(jìn)化關(guān)系、GC含量及D1蛋白氨基酸三肽元件模型Fig.2 Phylogenetic relationship, GC content and the amino acid tripletsequences from D1 protein motifs of cyanophage psbA gene in Yilong Lake

將從異龍湖水樣中獲得的噬藻體psbA序列與具有代表性的藍(lán)藻、綠藻、噬藻體分離株及噬藻體psbA環(huán)境序列等參考序列在核苷酸水平構(gòu)建系統(tǒng)進(jìn)化樹(shù)。如圖3所示，所有的序列在進(jìn)化樹(shù)結(jié)構(gòu)上可將其劃分為10個(gè)進(jìn)化簇(Cluster A～Cluster J)。該研究獲得的15條異龍湖噬藻體psbA序列全部分布在Cluster A和C中。其中，Cluster A包括來(lái)自中國(guó)東北稻田[21]和日本稻田水體psbA序列[22]和該研究得到的11條異龍湖序列；Cluster B主要來(lái)自海洋[23]的噬藻體分離株；Cluster C單獨(dú)成簇，包括該研究得到的4條異龍湖序列；Cluster D主要來(lái)自于北冰洋[24]和海洋原綠球藻(Prochlorococcus)及其噬藻體類(lèi)群；Cluster E和F主要來(lái)自于日本稻田水體及奧杜藻屬(Odontellarhomboides)的噬藻體psbA序列類(lèi)群[25]；Cluster G主要來(lái)自于淡水伊利湖的噬藻體類(lèi)群[26]；Cluster H、I和J包括來(lái)自藍(lán)藻[27]、淡水湖[28]、原綠球藻、海洋聚球藻(Synechococcus)、日本稻田水體及其噬藻體分離株的psbA序列。

圖3 異龍湖與稻田生態(tài)系統(tǒng)及其他環(huán)境psbA序列構(gòu)建的系統(tǒng)進(jìn)化樹(shù)Fig.3 Neighbor-joining phylogenetic tree showing the relationship of psbA amino acid sequence from Yilong Lake, paddy floodwater ecosystems and other environments

Cluster A的克隆來(lái)自稻田生態(tài)系統(tǒng)和該研究的11條序列，形成6個(gè)進(jìn)化亞簇，分別為來(lái)自中國(guó)東北稻田水體中的psbA基因序列形成的4個(gè)進(jìn)化亞簇A2、A3、A4和A6；1個(gè)來(lái)自日本稻田水體的psbA基因序列形成的進(jìn)化亞簇A5；其中該研究獲得的9條異龍湖psbA基因序列形成一個(gè)進(jìn)化小亞簇(YLW-1)，與中國(guó)東北稻田水體亞簇的psbA序列進(jìn)化距離相近，和日本稻田水體共同形成進(jìn)化亞簇A1；1條異龍湖序列YL09S-9與中國(guó)東北稻田水體psbA基因序列同處于A4亞簇；還有1條異龍湖序列YL09S-10與中國(guó)東北稻田水體psbA基因序列形成的A6進(jìn)化亞簇處于同一進(jìn)化亞簇。Cluster C包括該研究得到的4條異龍湖序列(YL04S-1、YL04S-2、YL04S-3和YL09S-6)，形成一個(gè)相對(duì)獨(dú)立的進(jìn)化簇YLW-2，與Cluster D～J在進(jìn)化距離上相隔較遠(yuǎn)，推測(cè)為新的噬藻體psbA基因類(lèi)群。

結(jié)合圖2可知，Cluster A中的異龍湖psbA基因的氨基酸序列可變的三肽元件除序列YL09S-10之外均為EEV型；且YL09S-9與其他序列進(jìn)化距離較遠(yuǎn)，所以其余9條序列聚在一簇。在Cluster C中，4條異龍湖序列氨基酸序列可變的三肽元件均為ETE型，并且單獨(dú)成簇；與Cluster D～H中所包含的海洋、稻田等psbA序列的進(jìn)化距離相隔較遠(yuǎn)。云南高原湖泊異龍湖水樣中psbA基因序列和Cluster G中其他淡水湖泊(北美洲伊利湖)中的psbA基因序列在系統(tǒng)進(jìn)化樹(shù)上的進(jìn)化距離表明，雖都是淡水湖，但兩者進(jìn)化歷史相互獨(dú)立，僅存在有限的基因交流。此外，大部分異龍湖psbA克隆都在Cluster A中，且與中國(guó)東北稻田的psbA序列相隔較近，這表明異龍湖噬藻體psbA基因的分布狹窄，遠(yuǎn)離其他淡水湖泊和海洋的分布，與海洋、國(guó)外淡水湖和日本稻田水體不同，與中國(guó)東北稻田的生態(tài)類(lèi)群進(jìn)化距離較近，進(jìn)化相對(duì)獨(dú)立。從異龍湖不同季節(jié)春季和秋季采集水樣中分別獲得的3條和12條psbA基因序列，秋季水樣中psbA序列分布于Cluster A和Cluster C這2個(gè)進(jìn)化簇中，春季的3條psbA基因序列分布于Cluster C中，顯示秋季噬藻體psbA基因序列的遺傳多樣性更為豐富。

3 討論與結(jié)論

在該研究中，采用光合作用基因psbA作為靶標(biāo)基因，光合作用基因在噬藻體內(nèi)廣泛存在且十分保守。此前的研究中曾用到g20基因，但是相對(duì)于psbA基因來(lái)說(shuō)，g20基因只局限存在于肌尾病毒科噬藻體，對(duì)于短尾病毒科的噬藻體，g20無(wú)法檢測(cè)出來(lái)，因此，使用g20作為研究噬藻體遺傳多樣性的分子標(biāo)記有一定的局限性。CHENARD[24]在研究淡水和海洋中噬藻體的多樣性時(shí)使用psbA基因作為分子標(biāo)記的靶標(biāo)基因；該研究主要參照上述文獻(xiàn)的實(shí)驗(yàn)技術(shù)方法，選用psbA基因作為研究異龍湖噬藻體多樣性的分子標(biāo)記基因，結(jié)果證明psbA基因能很好地揭示噬藻體的遺傳多樣性。

為區(qū)分環(huán)境樣品中獲得的psbA基因是噬藻體來(lái)源還是宿主來(lái)源，目前國(guó)內(nèi)外相關(guān)研究中JING等[11]、SHARON等[20]均采用了CHENARD等[24]建議的方法：(1)將采集的水樣先經(jīng)濾紙過(guò)濾，再分別經(jīng)0.45和0.22 μm孔徑的濾膜過(guò)濾以除去藻類(lèi)宿主等細(xì)胞型生物，濾液通過(guò)超濾濃縮，對(duì)濃縮液進(jìn)行psbA基因的PCR擴(kuò)增前，先用DNase處理以除去濃縮水樣中可能含有的游離藍(lán)藻DNA，用藍(lán)藻ITS基因序列擴(kuò)增結(jié)果為陰性后，再提取濃縮水樣中的DNA；(2)基于噬藻體和宿主的psbA基因序列的GC含量分析表明，噬藻體的GC含量低于其宿主的GC含量；(3)基于編碼D1蛋白的氨基酸序列模型R/KETTXXXSQ/H中三肽元件XXX變化的區(qū)分方法表明,該方法從采集水樣的處理、已有噬藻體和宿主的psbA基因序列的GC含量以及基于編碼D1蛋白的氨基酸序列模型R/KETTXXXSQ/H中三肽元件XXX變化的統(tǒng)計(jì)分析后提出，可以區(qū)分出環(huán)境樣品中獲得的psbA基因來(lái)源。

在對(duì)高原湖泊異龍湖中噬藻體psbA基因多樣性進(jìn)行研究時(shí)發(fā)現(xiàn)，與淡水、海洋環(huán)境比較，異龍湖噬藻體序列與來(lái)自中國(guó)東北稻田水體中噬藻體psbA基因群集更為相似，與其他淡水湖泊、海洋噬藻體psbA基因進(jìn)化距離較遠(yuǎn)；不同環(huán)境水域中噬藻體含有其獨(dú)特的類(lèi)型，顯示不同自然環(huán)境中噬藻體psbA基因組成存在一定差異，psbA基因的分布與其存在的環(huán)境有很大的關(guān)系。ZHONG等[29]對(duì)法國(guó)的安納西湖和布爾熱湖的噬藻體多樣性進(jìn)行分析，發(fā)現(xiàn)這2個(gè)湖泊中的噬藻體psbA存在著一定的差異。

CHENARD等[24]對(duì)海洋和淡水中不同地點(diǎn)的噬藻體psbA基因進(jìn)行擴(kuò)增，并通過(guò)構(gòu)建噬藻體psbA基因的系統(tǒng)進(jìn)化樹(shù)發(fā)現(xiàn)，海洋中的克隆劃分為7個(gè)類(lèi)群，淡水中的克隆則為2個(gè)來(lái)自不同湖泊的獨(dú)立類(lèi)群，表明海洋和淡水中的噬藻體存在著較大的差異。WANG等[22]對(duì)日本稻田水中的噬藻體多樣性進(jìn)行了研究，利用PCR擴(kuò)增得到37條不同的噬藻體psbA序列，發(fā)現(xiàn)噬藻體的遺傳多樣性與其宿主有關(guān)，其群落結(jié)構(gòu)的豐富度與藍(lán)藻多樣性呈正比。荊瑞勇等[8]在對(duì)中國(guó)東北稻田水體的研究中獲得17條噬藻體psbA序列，并發(fā)現(xiàn)其中存在新的psbA基因類(lèi)群。該研究對(duì)富營(yíng)養(yǎng)化湖泊異龍湖水體的研究得到的結(jié)果與上述結(jié)果相似，并發(fā)現(xiàn)其中也存在新的psbA基因類(lèi)群；從2015年春季(4月)及秋季(9月)采集的水樣中分別獲得3條和12條psbA環(huán)境系列及其在進(jìn)化樹(shù)中的分布，顯示異龍湖秋季噬藻體psbA基因遺傳多樣性更加豐富，表明異龍湖噬藻體psbA基因序列在不同季節(jié)的水樣中存在差異；這可能與不同季節(jié)異龍湖的水溫、pH值等環(huán)境因素的變化有關(guān)。水溫是影響浮游藻類(lèi)生長(zhǎng)的重要因子，適宜的水溫有利于浮游藻類(lèi)的生長(zhǎng)[30]，即在一定的溫度范圍內(nèi)，水溫升高可促進(jìn)藻類(lèi)光合作用。異龍湖春季(4月)水溫(20.2 ℃)較低，秋季(9月)水溫為23.9 ℃，秋季水溫較高對(duì)浮游藻類(lèi)的生長(zhǎng)更加有利，同時(shí)異龍湖秋季水樣pH值(8.6)較春季水樣(7.8)高，弱堿性的水體也適宜藻類(lèi)生長(zhǎng)，而噬藻體群落結(jié)構(gòu)的豐富度與藍(lán)藻的多樣性呈正比，從而也影響了異龍湖噬藻體種類(lèi)的不同，顯示異龍湖噬藻體隨季節(jié)、溫度等環(huán)境因素的變化而導(dǎo)致的主要基因型的差異，表明異龍湖噬藻體基因遺傳多樣性與環(huán)境水質(zhì)參數(shù)有一定的關(guān)系。

綜上所述，以psbA基因?yàn)榘袠?biāo)基因，對(duì)高原湖泊異龍湖中噬藻體類(lèi)群的分析研究表明，異龍湖水體中噬藻體psbA基因類(lèi)群與海洋、其他淡水湖和日本稻田水體psbA基因類(lèi)群不同，與中國(guó)東北稻田的生態(tài)類(lèi)群進(jìn)化距離較近，進(jìn)化相對(duì)獨(dú)立，存在新的噬藻體psbA基因類(lèi)群，并且在不同季節(jié)的水樣中噬藻體psbA基因序列存在差異，秋季其遺傳多樣性更加豐富。