薛文凱,朱攀,德吉,郭小芳
西藏大學理學院,西藏 拉薩 850000
青藏高原是世界上海拔最高、面積最大的高原,其高海拔、低氧、復雜的地貌和脆弱的生態(tài)環(huán)境形成了獨特的生態(tài)系統(tǒng)(Tang et al.,2017)。青藏高原的湖泊系統(tǒng)由1500多個大小不同的湖泊組成,其面積占中國湖泊總面積的49.5%(Wan et al.,2014)。青藏高原具有山脈眾多,太陽輻射強、降水少、蒸發(fā)量大等地理和氣候特征,從而形成了以咸水湖和鹽湖為主的湖泊系統(tǒng)(Zhang et al.,2019)。納木措作為青藏高原代表性的咸水湖,海拔高達4720 m,水面積超過500 km2,是地球上最高的湖泊之一(Spencer et al.,2014)。水體中的微生物在有機物和污染物的生物降解或轉(zhuǎn)化中起著關鍵作用,維持著水體的生態(tài)平衡。某些微生物的存在與特定類型的污染密切相關,是作為評估水質(zhì)和氣候的重要生物指標(Glasl et al.,2017)。
絲狀真菌作為組成水體微生物的重要類群,它們適應力強,廣泛分布于各種生境中,即便在高鹽度的死海中仍能存活較長時間(Kinsey et al.,2003;Kis-Papo et al.,2003)。目前關于青藏高原水體微生物的研究主要包括細菌和真菌的多樣性及影響因子(王艷紅等,2021;鄭艷艷等;2022);自然生態(tài)系統(tǒng)條件下細菌群落的多樣性(Wang et al.,2018;Antunes et al.,2022);生態(tài)位方面的研究以原生動物為主(安瑞志等,2021;張鵬等,2022)。但用水體絲狀真菌的群落演變特征來反映水文環(huán)境變化的研究較少。冰川和高山湖泊生態(tài)系統(tǒng)是對氣候變化和人為影響極為敏感的環(huán)境。隨著旅游業(yè)的發(fā)展,人類活動對高原環(huán)境的影響日益凸顯,已有研究表明近年來青藏高原周圍地區(qū)的大氣污染和黑煙越來越多,導致高原上污染物的沉積增加(Xu et al.,2009;Cong et al.,2013;Spencer et al.,2014)。受全球變暖的影響,降水、冰川和永久凍土的消融增強,導致湖泊儲水量增加和水文條件發(fā)生改變(Zhang et al.,2020),高原湖泊正面臨著污染和氣候變化帶來的威脅。
優(yōu)勢種是群落中占優(yōu)勢地位的種群,研究優(yōu)勢種及其生態(tài)位對理解群落中各種群間在時空上的相互關系以及對資源的競爭程度具有重要意義(Grinnell,1917;Helaou?t et al.,2009)。同時低溫的湖水環(huán)境也有較大幾率分離到耐冷的菌株。因此本研究以納木措春、夏、秋 3個季節(jié)的沿岸水體為樣本,分離純化出其中的可培養(yǎng)絲狀真菌,探究絲狀真菌優(yōu)勢種的時空生態(tài)位和季節(jié)變化特征以及其群落中物種組成的變化規(guī)律,以期為利用微生物物種組成特征作為水質(zhì)和氣候變化的生物指標提供理論基礎,為低溫環(huán)境污水的生物修復提供菌株資源。
在2020年8月以及2021年5月和10月,沿湖區(qū)周邊選擇 20 個樣點(91°02′06″—90°15′31″E,30°54′32″—30°30′14″N),各樣點相距 15 km 左右(納木措禁止船只航行;為了在納木措湖周邊均勻設置樣點,在 Google地圖上按納木措湖周長與樣點數(shù)量的比值定位標記各采樣點)。用消毒過的水樣采集器,采集離湖岸15 m處的水樣,各樣點采5 L水,盛于2個無菌桶中,用于環(huán)境因子的測定和絲狀真菌的分離。在樣點信息記錄中,除 19號樣點為觀景臺,樣點存在地表徑流外,他樣點周邊植被稀疏,明顯特征。
測定納木措水體的 8種環(huán)境因子,水體溫度(t)、pH、總?cè)芙夤腆w量(TDS)、電導率(EC)和鹽濃度(Salt),以多功能參數(shù)儀(Multi-Parameter PCSTestrTM35)進行原位測定,3次重復。銨態(tài)氮(NH4+-N)、總氮(TN)以及總磷(TP)的測定,由西藏博源檢測公司完成,其中,銨態(tài)氮(NH4+-N)測定采用納氏試劑分光光度法;總氮(TN)采用堿性過硫酸鉀消解紫外分光光度法;總磷(TP)采用鉬酸銨分光光度法(邢殿樓等,2016)。
用孔徑0.45 μm、直徑50 mm的無菌醋酸纖維素濾膜,抽濾50 mL和100 mL的湖水,每梯度各抽濾3個重復,抽濾后將濾膜貼于固體RBCH、酸性PDA、酸性YM平板培養(yǎng)基上進行20 ℃的恒溫培養(yǎng),培養(yǎng)基配比參照亞萍等(2013)。3—5 d后統(tǒng)計平板上的菌落數(shù),挑取單個菌落接種在RBCH斜面培養(yǎng)基上純化,繼續(xù)恒溫培養(yǎng)相同時間后觀察并記錄其形態(tài)特征(若有雜菌污染則再次純化),將純化好的菌株保存在體積分數(shù)為 30%的甘油管中(?20 ℃),用于高原微生物種質(zhì)資源庫的建設。
分別用點植法和插片培養(yǎng)法培養(yǎng)純化菌種,觀察絲狀真菌的菌落形態(tài)和顯微特征,根據(jù)真菌的分類鑒定指標對所屬類群做初步鑒定(Ainsworth,2008)。分子鑒定由上海生工完成,其中DNA的提取、PCR的擴增采用CTAB法(孫立夫等,2009),引物為 ITS1(5′-TCCGTAGGTGAACCTGCGG-3′)和 ITS4(5′-TCCTCCGCTTATTGATATGC-3′),用BioEdit軟件校對 DNA序列并去除引物和雙峰序列,然后在NCBI核酸序列數(shù)據(jù)庫中搜索同源性菌株(選相似度98%以上的模式菌株)(Landeweert et al.,2003),最終確定絲狀真菌的分類地位。
利用R 4.0.2和Excel 2018軟件計算絲狀真菌的多樣性指數(shù)、優(yōu)勢種的更替率(r)、優(yōu)勢度(Y)、出現(xiàn)頻率(fi)以及生態(tài)位寬度(Bi)、生態(tài)位重疊(Oik),并對水環(huán)境因子的差異性、絲狀真菌的多樣性、優(yōu)勢種的時空異質(zhì)性以及優(yōu)勢種與水環(huán)境因子的相關性進行分析。
1.4.1 多樣性指數(shù)
通過計算豐富度指數(shù)(R)、Shannon-Wiener多樣性指數(shù)(H')、Pielou均勻度指數(shù)(J')、Simpson多樣性指數(shù)(D)來判斷納木措水體絲狀真菌的多樣性。
式中:
S——絲狀真菌物種數(shù);
Pi——i種的個體數(shù)與總個體數(shù)的比值;
N——某樣點物種的總數(shù)。
1.4.2 優(yōu)勢種的劃分及更替率
確定優(yōu)勢種采用McNaughton優(yōu)勢度計算方法(McNaughton,1967;Habib et al.,1997),優(yōu)勢種的季節(jié)更替率計算參考文獻(張皓等,2016;梁淼等,2018)。
式中:
fi——i種的出現(xiàn)頻率;
N——樣點內(nèi)的物種總數(shù);
ni——i物種的總數(shù);
a、b——相鄰兩季優(yōu)勢種的種類數(shù);
c——相鄰兩季共有優(yōu)勢種的種類數(shù),(Y>0.02為優(yōu)勢種,Y>0.1為絕對優(yōu)勢種)。
1.4.3 生態(tài)位
生態(tài)位寬度值的計算采用Colwell et al.(1971)加權(quán)修正的Levins指數(shù)。
生態(tài)位重疊值的計算采用Pianka指數(shù)(1973)。
式中:
Bi——物種i的生態(tài)位寬度值;
Pij——空間生態(tài)位中代表i種在j樣點的個體數(shù)占該樣點總個體數(shù)的比例;
r——總樣點數(shù);
Pij——時間生態(tài)位中代表i種在j季節(jié)的個體數(shù)占i種在所有季節(jié)總個體數(shù)的比值;
r——總季節(jié)數(shù);
Oik——物種i與物種k的生態(tài)位重疊指數(shù),時空生態(tài)位寬度=時間生態(tài)位寬度×空間生態(tài)位寬度。
對納木措3個季節(jié)的8種水環(huán)境因子做差異性分析(圖1)。pH在春季和夏季、春季和秋季間存在極顯著差異(P<0.01),其中春季高于秋季、夏季;溫度在3個季節(jié)間的差異均為極顯著(P<0.01),且夏季>秋季>春季;NH4+-N在春季和秋季、夏季和秋季間呈顯著差異(P<0.05);TP在春季和秋季間呈顯著差異(P<0.05),在夏季和秋季間呈極顯著差異(P<0.01);TN僅在春季和夏季間存在顯著差異(P<0.05),以上3種環(huán)境因子均為春季最高,秋季最低,此外,TDS、EC和Salt在3個季節(jié)中均無顯著差異。總體而言,湖水中的部分環(huán)境因子在季節(jié)更替中出現(xiàn)了顯著變化。
圖1 納木措不同季節(jié)水體環(huán)境因子差異Figure 1 Differences of water environmental factors in different seasons in Nam Co Lake
在納木措3個季節(jié)中共鑒定出3359株絲狀真菌,如圖2a所示,絲狀真菌總豐度在夏季與春季、夏季與秋季間呈極顯著差異(P<0.01),且夏季最多,為1412株,其次為秋季(1026株),春季最少(921株)。絲狀真菌物種數(shù)在秋季與夏季、秋季與春季間呈極顯著差異(P<0.005),其中春季包括 20屬62種,優(yōu)勢種9個,夏季包括22屬47種,優(yōu)勢種13個,秋季包括13屬47種,優(yōu)勢種7個。絲狀真菌的多樣性指數(shù)(圖2b),除Shannon-Wiener多樣性指數(shù)、Simpson多樣性指數(shù)在春季與夏季間無顯著差異外,其他季節(jié)間的各指數(shù)均存在顯著(P<0.05)和極顯著差異(P<0.01),其中秋季的Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)均為最高,春季的Simpson多樣性指數(shù)最高。
圖2 納木措不同季節(jié)水體絲狀真菌多樣性指數(shù)的比較Figure 2 Comparison of diversity indices of filamentous fungi in different seasons in Nam Co Lake
為進一步分析納木措各季節(jié)水體絲狀真菌優(yōu)勢種的時空異質(zhì)性(圖 3),采用 Bray-Curtis距離的非度量多維尺度(NMDS)分析,按季節(jié)分為 3組,發(fā)現(xiàn)秋季的2號和15號以及夏季的14號樣點中優(yōu)勢種群落結(jié)構(gòu)與本季節(jié)其他樣點差異較大,其余各樣點的優(yōu)勢種群落結(jié)構(gòu)在同一季節(jié)內(nèi)的差異較小。ANOSIM 結(jié)果顯示,R=0.939,P<0.01,說明各季節(jié)間的差異顯著大于季節(jié)內(nèi)樣點間的差異。
圖3 納木措不同季節(jié)水體絲狀真菌優(yōu)勢種的空間變異分析Figure 3 Spatial variation analysis of the dominant species of filamentous fungi in different seasons in Nam Co Lake
對納木措3個季節(jié)的水體絲狀真菌優(yōu)勢種進行劃分,共得到優(yōu)勢種23種,如表1所示,普通青霉Penicillium commune為3個季節(jié)的共有優(yōu)勢種,酒色青霉Penicillium vinaceum、凍土毛霉Mucor hiemalis為春季和夏季的共有優(yōu)勢種,殼青霉Penicillium crustosum、灰玫瑰青霉Penicillium griseoroseum為春季和秋季的共有優(yōu)勢種。春季和秋季有獨有優(yōu)勢種4個,夏季有獨有優(yōu)勢種10個。由春季進入夏季時,優(yōu)勢屬Penicillium的占比下降,優(yōu)勢種的數(shù)量增加;由夏季進入秋季時,優(yōu)勢屬Penicillium的占比上升,優(yōu)勢種的數(shù)量下降。經(jīng)計算,夏季與秋季優(yōu)勢種的更替率最高,為95%,春季與夏季優(yōu)勢種的更替率次之,為84%,春季與秋季優(yōu)勢種的更替率最低,為77%,說明納木措水體絲狀真菌優(yōu)勢種的季節(jié)性變化較大。
表1 納木措水體絲狀真菌優(yōu)勢種在不同季節(jié)的優(yōu)勢度(Y)與出現(xiàn)頻率(fi)Table 1 Dominance (Y) and frequency (fi) of dominant species of filamentous fungi in different seasons in Nam Co Lake
如表2所示,納木措水體絲狀真菌優(yōu)勢種的空間生態(tài)位寬度值區(qū)間在 0.017—0.942,其中殼霉菌Stagonosporopsissp.的空間生態(tài)位寬度值最高,比阿婁維扎青霉Penicillium bialowiezense的空間生態(tài)位寬度值最低。時間生態(tài)位寬度值區(qū)間在 0.333—0.979,其中P.commune的時間生態(tài)位寬度值最高,波蘭青霉Penicillium polonicum等19個菌種的時間生態(tài)位寬度值均較低。絲狀真菌優(yōu)勢種的時空生態(tài)位寬度值均較小,區(qū)間在 0.006—0.314,其中P.bialowiezense的時空生態(tài)位寬度值最低,Stagonosporopsissp.的時空生態(tài)位寬度值最高,且處于窄生態(tài)位(0≤Bi<0.15)的菌種占比較大,這表明納木措水體絲狀真菌優(yōu)勢種易受時空變化的影響。
表2 納木措水體絲狀真菌優(yōu)勢種在不同季節(jié)的生態(tài)位寬度(Bi)Table 2 Niche width of dominant species of filamentous fungi in different seasons in Nam Co Lake
在空間維度上(表3),納木措各樣點絲狀真菌優(yōu)勢種所組成的群落間的生態(tài)位重疊度較高,其中絲狀真菌群落在160對樣點間高度重疊(Oik>0.8),占總對數(shù)的 84%;在 24對樣點間中度重疊(0.6 表3 納木措各樣點水體絲狀真菌群落的生態(tài)位重疊(Oik)Table 3 Niche overlap of filamentous fungal communities at different sites in Nam Co Lake 表4 納木錯各季節(jié)水體絲狀真菌群落的生態(tài)位重疊(Oik)Table 4 Niche overlap of filamentous fungal communities in different seasons in Nam Co Lake 分析納木措不同季節(jié)的水環(huán)境因子,發(fā)現(xiàn)總?cè)芙夤腆w量、電導率和鹽濃度在3個季節(jié)中均無顯著差異,其他 5種環(huán)境因子在各季間存在顯著(P<0.05)或極顯著差異(P<0.01)。水環(huán)境因子與絲狀真菌優(yōu)勢種的Spearman等級相關性分析(圖4)顯示,在春季,優(yōu)勢種P.vinaceum、M.hiemalis分別與全磷、銨態(tài)氮呈極顯著負相關和極顯著正相關(P<0.01)。在夏季,優(yōu)勢種Stagonosporopsissp.、高山被孢霉Mortierella alpina、M.hiemalis與電導率、鹽濃度、總?cè)芙夤腆w量呈顯著(P<0.05)或極顯著負相關(P<0.01);優(yōu)勢種兩型蠟蚧菌Lecanicillium dimorphum與溫度、pH分別呈顯著正相關和顯著負相關(P<0.05),溫度又與小不整球殼菌Plectosphaerella plurivora呈顯著正相關(P<0.05),pH與微結(jié)節(jié)霉菌Microdochium bolleyi呈顯著負相關(P<0.05);優(yōu)勢種微糙枝孢Cladosporium asperulatum與全氮、銨態(tài)氮呈顯著負相關(P<0.05),同樣全氮與鏈格孢Alternaria chlamydosporigena呈顯著負相關(P<0.05);優(yōu)勢種M.hiemalis與全磷呈顯著正相關(P<0.05),與銨態(tài)氮呈極顯著正相關(P<0.01)。秋季,優(yōu)勢種枝狀枝孢Cladosporium cladosporioides與鹽濃度、電導率、總?cè)芙夤腆w量呈顯著正相關(P<0.05),P.bialowiezense與銨態(tài)氮、pH以及P.griseoroseum與全磷均呈顯著負相關(P<0.05)。 圖4 納木措不同季節(jié)水體環(huán)境因子與絲狀真菌優(yōu)勢種豐度的相關性Figure 4 Correlation between water environmental factors and the abundance of dominant species of filamentous fungi in different seasons in Nam Co Lake 納木措3個季節(jié)的水體中共鑒定出3359株絲狀真菌,各季節(jié)的多樣性均較高,且α多樣性指數(shù)間大多存在顯著差異(P<0.05)。時空異質(zhì)性的結(jié)果表明,絲狀真菌優(yōu)勢種在各季節(jié)間的差異顯著大于季節(jié)內(nèi)各樣點的差異。有相似研究表明,生境差異能使湖泊中細菌群落的多樣性和組成產(chǎn)生異質(zhì)性(Jiao et al.,2020;Zhang et al.,2021);納木措不同季節(jié)間絲狀真菌的多樣性和群落結(jié)構(gòu)存在顯著的時間異質(zhì)性,且時間異質(zhì)性大于空間異質(zhì)性,各季節(jié)間環(huán)境因子的差異形成生境的差異,可能是造成這一結(jié)果的重要原因。絲狀真菌優(yōu)勢種的數(shù)量和物種組成隨著季節(jié)的改變發(fā)生顯著變化,各季節(jié)優(yōu)勢種的更替率也較高,如其中的優(yōu)勢屬Penicillium,由春季進入夏季,夏季進入秋季,其占比先下降后上升,這可能與Penicillium在低溫環(huán)境中的優(yōu)勢更強有關(Cantrell et al.,2011)。結(jié)合不同季節(jié)間環(huán)境因子的變化幅度來看,環(huán)境的變化越大可能是群落越不穩(wěn)定,更替率越大的重要原因。 生態(tài)位寬度是反映種群對環(huán)境適應能力和資源利用能力的重要指標(李德志等,2006;彭松耀等,2015)。在空間維度上,Stagonosporopsissp.的空間生態(tài)位寬度值最高,受空間資源的限制較小;P.bialowiezense的空間生態(tài)位寬度值最低,受空間資源的限制較大(何雄波等,2018)。在時間維度上,P.commune的時間生態(tài)位寬度值最高,受季節(jié)變化的限制較小;P.polonicum等19個菌種的時間生態(tài)位寬度值均較低,受季節(jié)變化的限制較大(黃保宏等,2005;馬玲等,2012)。絲狀真菌優(yōu)勢種的時空生態(tài)位寬度值均較小,且窄生態(tài)位的菌種占比較大,總體上納木措水體絲狀真菌優(yōu)勢種易受時空變化的影響。生態(tài)位重疊值反映著種群間對資源的競爭強弱(求錦津等,2018),納木措各樣點絲狀真菌優(yōu)勢種所組成的群落間的生態(tài)位重疊度較高,在空間維度上,各樣點絲狀真菌優(yōu)勢種組成的群落高度重疊的占84%;在時間維度上,3個季節(jié)的絲狀真菌優(yōu)勢種所組成的群落間的生態(tài)位重疊度同樣較高(77.2%—86.4%)。表明各絲狀真菌群落的競爭較為激烈,且群落結(jié)構(gòu)不穩(wěn)定正處于發(fā)展的狀態(tài)。雖然不同季節(jié)的絲狀真菌群落間不會直接產(chǎn)生競爭關系,但它們之間較高的生態(tài)位重疊正是構(gòu)建各個季節(jié)不同群落結(jié)構(gòu)的重要動力。總體而言,不同水域或不同季節(jié)塑造了水體中絲狀真菌不同的群落結(jié)構(gòu),不同微生物種群對各種環(huán)境因子有著不同的耐受度和適應力(Williams et al.,2008),競爭優(yōu)勢隨著環(huán)境變化而改變,群落發(fā)生演替,表現(xiàn)為各季節(jié)絲狀真菌的多樣性和物種組成產(chǎn)生異質(zhì)性以及優(yōu)勢種在各季節(jié)間較高的更替率。 微生物群落的多樣性和物種組成受到CO2、有機物、氮、磷和溫度等各種環(huán)境因素的影響(Zhang et al.,2014;Li et al.,2019;Xia et al.,2020)。時空變化對絲狀真菌群落結(jié)構(gòu)的影響本質(zhì)上是不同水域和不同季節(jié)的水環(huán)境因子與絲狀真菌各種群相互作用的體現(xiàn),組成微生物群落的各物種彼此之間對資源存在廣泛的競爭(Esser et al.,2015;Hassani et al.,2018),環(huán)境的變化能直接影響它們的競爭力,因此,環(huán)境條件的變化造成了微生物群落的演替。相關研究表明水體真菌的組成和分布等受環(huán)境因子的影響(王艷紅等,2021),本研究中,不同季節(jié)的優(yōu)勢種與水環(huán)境因子的相關性差別較大,如夏季的兩個優(yōu)勢種與溫度具有顯著的正相關(P<0.05),且夏季湖水的溫度顯著高于春季和秋季,說明溫度的升高可能是造成二者逐漸成為優(yōu)勢種的原因之一。 對環(huán)境響應敏感,受水體環(huán)境限制較強的某些原生動物類群作為環(huán)境指示種被認為有著可靠的生態(tài)學意義(汪志聰?shù)龋?010)。通過與其他學者相似的研究比較發(fā)現(xiàn)(表5),本研究中的絲狀真菌優(yōu)勢類群對水環(huán)境的變化有著更敏感的響應,在季節(jié)變化中,群落的更替率和構(gòu)成群落的物種變化更大。目前已有研究表明,青藏高原的變暖率遠高于全球平均水平(Liu et al.,2000),全球變暖造成冰川消融增加,冰川融水的大量輸入改變了高原湖泊的鹽度、養(yǎng)分和生物多樣性(Milner et al.,2009),這些水文條件的變化可能對湖泊生態(tài)系統(tǒng)的結(jié)構(gòu)和功能產(chǎn)生嚴重影響(Lin et al.,2017)。且微生物對環(huán)境變化的高度敏感性,使得用生物方法評判水質(zhì)和氣候的細微變化成為可能,雖然微生物檢測法需要在統(tǒng)一的標準下(如采樣方法、培養(yǎng)基的種類、稀釋倍數(shù)、培養(yǎng)環(huán)境等)才有可比性,但隨著技術(shù)的發(fā)展和分離方法的完善,利用微生物的多樣性、群落結(jié)構(gòu)、生態(tài)位等信息作為環(huán)境變化的指標在未來具有巨大的應用潛力。 表5 納木措微生物與其他區(qū)域微生物的生態(tài)位研究對比Table 5 Comparison of niche studies of microorganisms in Nam Co Lake with those in other regions (1)納木措不同季節(jié)間絲狀真菌的多樣性和群落結(jié)構(gòu)存在顯著的異質(zhì)性,時間資源比空間資源的影響更顯著,各季節(jié)環(huán)境因子的差異以及不同時期環(huán)境因子與絲狀真菌之間相關性的差異是驅(qū)動其群落構(gòu)建的主要動力。 (2)在時空生態(tài)位寬度以及生態(tài)位重疊上,大多優(yōu)勢種處于窄生態(tài)位,且群落間重疊度較高,絲狀真菌的群落結(jié)構(gòu)不穩(wěn)定,在不同季節(jié)的環(huán)境條件下,一些菌種的競爭優(yōu)勢隨之改變,群落發(fā)生演替,塑造了各季節(jié)絲狀真菌的多樣性和組成差異,表現(xiàn)為優(yōu)勢種在各季節(jié)間較高的更替率。 (3)相比于其他類型的水體微生物生態(tài)位,絲狀真菌優(yōu)勢類群對水環(huán)境的變化有著更敏感的響應,這種高度敏感性,可以用來評估水質(zhì)和氣候的微小變化,具有較大的應用潛力。2.5 納木措不同季節(jié)水體環(huán)境因子與絲狀真菌優(yōu)勢種的相關性分析
3 討論
4 結(jié)論