曾 帆,吳田振,于振朋,孫林霞,郭天慧

(南京師范大學(xué)生命科學(xué)學(xué)院,江蘇省生物多樣性與生物技術(shù)重點實驗室,江蘇 南京 210023)

趨同進化(convergent evolution)是進化生物學(xué)的研究熱點,是指親緣關(guān)系較遠的物種在自然選擇的作用下,為了適應(yīng)相似的生境而產(chǎn)生的相同或相似的表型、生理以及行為特征[1]. 盡管自然界中趨同進化的例子很多,但是趨同進化的分子機制研究還很少. 大熊貓(Ailuropodamelanoleuca)和小熊貓(Ailurusfulgensstyani),親緣關(guān)系較遠,通過基因組的測定,鑒定了多個與肢體發(fā)育和蛋白消化相關(guān)的基因出現(xiàn)趨同信號,為大、小熊貓演化出趨同的食竹食性以及具有抓握竹子的偽拇指提供了分子證據(jù)[2-3]. 蝙蝠和海豚在系統(tǒng)發(fā)生樹上相距甚遠,但均演化出回聲定位的能力,是趨同進化的又一典型,近些年鑒定了與聽覺相關(guān)的基因在二者之間具有平行/趨同進化位點,且在Perestin基因的系統(tǒng)發(fā)育樹中聚在一起,為二者趨同進化提供了很強的分子證據(jù)[4-5].

海洋哺乳動物(marine mammals),也稱為“海獸”,是哺乳動物中進化十分獨特的類群,它們由不同的陸生祖先分別從陸地重返海洋,親緣關(guān)系較遠,但為了適應(yīng)海洋環(huán)境,進化出一系列趨同的表型和生理特征,是研究趨同進化的理想模型. 它們都擁有與魚類類似的流線型體型以減少在水中游動時的阻力[6];皮下脂肪層增厚用以保持體溫,減少溫度散失;滲透調(diào)節(jié)能力增強用以適應(yīng)高滲的海洋環(huán)境;儲氧能力提高用以增加在水中的屏息時長[7];骨密度的變化用以適應(yīng)不同的潛水深度[8]. 除此以外,海洋哺乳動物普遍發(fā)生了味覺和嗅覺的喪失[9];鯨類和海牛發(fā)生了毛發(fā)的缺失[10],為趨同進化的研究提供了理想的模型.

滲透調(diào)節(jié)是生境轉(zhuǎn)換后海洋哺乳動物共同面臨的一個最主要的挑戰(zhàn). 形態(tài)學(xué)研究表明,許多海洋哺乳動物的腎臟演化出復(fù)腎結(jié)構(gòu),如鯨類、鰭足類、北極熊和水獺擁有數(shù)十甚至上千個小腎,且每一個小腎在形態(tài)結(jié)構(gòu)和生理作用上都相當于陸生哺乳動物的獨立腎,以此增加了腎臟的表面積,提高了腎臟的排泄能力,為海洋哺乳動物形成高濃度的濃縮尿提供了前提. 其中,完全營水生生活的鯨類的腎最為特殊,呈葡萄狀,小腎數(shù)目在種間差別很大,體型較大的須鯨可達6 000～7 000個[11-12],而體型較小的淡水豚類小腎一般僅有65～170個[13]. 前人研究發(fā)現(xiàn)鯨類尿液的滲透濃度以及Na+濃度均顯著高于大小與之相近或比之大的陸生哺乳動物,如鯨類尿液中尿素濃度比牛高4倍[14],提示鯨類通過排高度濃縮的尿液以調(diào)節(jié)水鹽平衡. 然而,海洋哺乳動物滲透調(diào)節(jié)的分子機制的研究才剛剛開始. 研究發(fā)現(xiàn)腎臟中與滲透調(diào)節(jié)相關(guān)的基因在海洋哺乳動物中出現(xiàn)了高表達的現(xiàn)象[15],進一步為海洋哺乳動物具有高效的水鹽調(diào)節(jié)提供了分子證據(jù). Xu等發(fā)現(xiàn)許多與滲透調(diào)節(jié)相關(guān)的基因在鯨類譜系中顯著受正選擇作用,如水通道蛋白2(aquaporin2,AQP2)和尿素轉(zhuǎn)運蛋白(urea transporters,UTA)編碼基因SLC14A2顯著受到正選擇,提示鯨類可能通過增強對水和尿素的轉(zhuǎn)運以產(chǎn)生高濃縮的尿液,以此適應(yīng)高滲的海洋環(huán)境,進一步證實鯨類產(chǎn)生濃縮尿是維持體內(nèi)水平衡的一個有效機制. 水鹽調(diào)節(jié)系統(tǒng)——腎素-血管緊張素-醛固酮系統(tǒng)(RAAS)中兩個關(guān)鍵的成員管緊張素轉(zhuǎn)換酶(angiotensin converting enzyme,ACE)和管緊張素原(angiotensinogen,AGT)鑒定到多個正選擇位點,特別是長江江豚中多個基因受正選擇,提示鯨類水鹽調(diào)節(jié)系統(tǒng)是維持滲透調(diào)節(jié)的又一重要機制[16-17].

尿素轉(zhuǎn)運蛋白(Urea transporters)是特異性轉(zhuǎn)運尿素的膜轉(zhuǎn)運蛋白,對機體的滲透調(diào)節(jié)起著關(guān)鍵性的作用[18]. 其中,UT-A是目前已克隆出的尿素轉(zhuǎn)運蛋白亞家族之一,共6個成員(UT-A1至UT-A6),由同一基因SLC14A2編碼并經(jīng)過不同的啟動子啟動轉(zhuǎn)錄后剪接形成. 在UT-A亞家族的成員中,UT-A1、UT-A2、UT-A3和UT-A4均在腎臟中表達,且SLC14A2基因敲除后尿素轉(zhuǎn)運蛋白表達缺陷的小鼠表現(xiàn)出嚴重的尿素濃縮功能缺陷,說明此基因在機體滲透調(diào)節(jié)過程中發(fā)揮了及其重要的作用[19-20]. 鯨類中該基因顯著受正選擇,提示鯨類通過提高尿素轉(zhuǎn)運能力,從而增強了機體的水鹽調(diào)節(jié)功能. 然而,該基因在其他海洋哺乳動物中的進化以及是否與鯨類一樣具有趨同的排高濃縮尿液的滲透調(diào)節(jié)機制尚不清楚. 據(jù)此,本研究下載了73個哺乳動物的SLC14A2基因,運用比較基因組學(xué)方法探討了以下科學(xué)問題:(1)滲透調(diào)節(jié)基因SLC14A2在海洋哺乳動物不同類群中是否受到相似的選擇壓力?(2)海洋哺乳動物間是否具有尿素轉(zhuǎn)運的趨同分子機制?

1 材料與方法

1.1 SLC14A2基因序列數(shù)據(jù)

本研究在哺乳動物中選取了73個物種,涵蓋了15個目,其中海洋哺乳動物包括鯨類25個代表性物種,鰭足類6個物種以及海牛類的1個物種(海牛).SLC14A2基因序列是在NCBI數(shù)據(jù)庫和OrthoMaM(Orthologous Mammalian Markers,http://orthomam2.mbb.univ-montp2.fr/OrthoMaM)數(shù)據(jù)庫中下載得到. 對于基因組尚未注釋的物種,通過本地BLAST,最終獲得了73個哺乳動物SLC14A2基因的序列. 由于本地BLAST可能會造成相鄰?fù)怙@子拼接處存在非3倍數(shù)堿基的插入或缺失,所以采用了一系列針對于存在移碼突變的序列比對方法,運用Macse軟件中的AlignSequence、Exportalignment程序、Hmmcleaner程序、Prank程序以及Gblocks程序,確定獲得序列為同源序列[21-23].

1.2 選擇壓力分析

氨基酸的非同義替代(non-synonymous substitution rate,dN)與同義替代(synonymous substitution rate,dS)的比值ω(dN/dS)常用作選擇壓力(selective pressure)的檢測指標(Kimura,1979). 如果ω=1為中性選擇(neutral selection);ω>1為正選擇(positive selection);ω<1為純化選擇(negative selection). 由于影響功能基因的選擇作用往往都只是發(fā)生在一些特定的支系的特定位點上,因此,為了檢測海洋哺乳動物進化過程中經(jīng)歷了怎樣的選擇壓力,本研究采用了基于密碼子的最大似然法(maximum likelihood,ML)的PAML4.7軟件包中的codeML程序,檢測特定支系或者位點上所受的選擇壓力[24].

首先,僅針對33種海洋哺乳動物SLC14A2基因序列的數(shù)據(jù)集,運用位點模型(site modal)中的三對模型(M7vs M8、M8 vs M8a和M1a vs M2a)檢測正選擇位點. 運用似然率檢測法(likelihoodratio tests,LRTs)評估兩個模型之間是否存在顯著差異,計算2ΔL值與自由度之間的卡方分布情況,若P值<0.05,則認為兩個模型之間存在顯著差異. 利用經(jīng)驗貝葉斯(bayes empirical bayes,BEB)方法計算M8和M2a模型中的正選擇位點的后驗概率(Posterior probability),如果后驗概率≥0.8,即認為該位點是潛在的正選擇位點[25]. 另外,使用在線網(wǎng)站Datamonkey(http://www.datamonkey.org/)中的固定效應(yīng)似然法(fixed effects likelihood,FEL)進一步檢測正選擇位點,設(shè)置FEL顯著性水平小于0.2[26]. 同時被兩種ML方法檢測的位點,為強烈的正選擇位點.

其次,在所有哺乳動物73個物種的SLC14A2基因序列的數(shù)據(jù)集中,選擇分支模型(branch model)和支位點模型(branch-site model)進一步探究正選擇作用是否局限于某些特定的進化譜系. 分支模型主要包括free ratio模型和two ratio模型. Free ratio模型允許系統(tǒng)發(fā)育樹上每一個進化譜系具有不同的ω值,而Two ratio模型允許前景支(即所感興趣的支系)與背景支之間存在不同的ω值. 運用LRT分別檢測Free ratio和Two ratio模型與零假設(shè)one ratio模型(即假設(shè)系統(tǒng)發(fā)生樹上所有支系具有相同的ω值)是否差異顯著. 如果模型顯著,則通過Free ratio檢測ω>1的正選擇支系,運用Two ratio模型檢測前景支海洋哺乳動物是否出現(xiàn)加速進化. 此外,進一步運用更加嚴格的分支位點模型(Ma)分析,該模型允許在特定支上的每一個位點都有獨立的ω值,將進化譜系分為前景枝(foreground,即感興趣的枝系)和背景枝(background),使用LRT對備擇假設(shè)Ma(正選擇模型:0<ω0<1,ω1>1和ω2≥1)和零假設(shè)Ma0(中性模型:0<ω0<1,ω1=1和ω2=1)進行比較. 當成對模型(Ma vs Ma0)差異顯著且ω>1時,接受正選擇模型Ma,且所檢測的位點P>0.8時,被認為是強烈的正選擇位點. 衡量枝位點模型顯著性的P值則通過多重檢驗錯誤發(fā)現(xiàn)率(false discovery rate,FDR)予以矯正,校正后的P<0.05時則為模型顯著.

1.3 海洋哺乳動物SLC14A2基因平行/趨同位點檢測

為了驗證海洋哺乳動物在二次入水過程中,是否發(fā)生了趨同進化,本研究使用PAMLv4.7中Codeml程序?qū)γ總€基因進行祖先序列重建. 通過對海洋哺乳動物各枝系中尿素轉(zhuǎn)運蛋白氨基酸的序列分析,尋找趨同氨基酸替代,得到海洋哺乳動物支系間從祖先到現(xiàn)生氨基酸變化. 根據(jù)預(yù)測的祖先位點的差異,編碼區(qū)基因序列的氨基酸替換事件分為平行替換和趨同替換兩種模式. 如果兩個進化枝具有相同的祖先氨基酸狀態(tài),并衍生出相同的現(xiàn)生氨基酸,便認為該位點為平行替換位點;反之,如果兩個進化枝具有不同的祖先氨基酸狀態(tài),并衍生出相同的現(xiàn)生氨基酸,則為趨同替換位點,這兩種模式都被認為是海洋哺乳動物趨同進化的證據(jù). 考慮到偶然氨基酸替代引起的干擾,本研究在JTT-fgen氨基酸替換模型下進行統(tǒng)計檢驗,比較實際觀察到的趨同位點數(shù)目與隨機期望出現(xiàn)的趨同位點數(shù). 若觀測值大于期望值,則認為存在氨基酸趨同或平行替換;反之,則不存在[27].

1.4 蛋白質(zhì)三維結(jié)構(gòu)

位于關(guān)鍵區(qū)域的氨基酸改變可能會通過影響構(gòu)型,進而對蛋白質(zhì)的功能產(chǎn)生影響. 為深入探究這些正選擇位點是否位于重要的功能域,本研究使用印太瓶鼻海豚(Tursiopsaduncus)SLC14A2基因的核苷酸序列在I-TASSER網(wǎng)站上(http://zhanglab.ccmb.med.umich.edu/I-TASSER/)進行了蛋白三維結(jié)構(gòu)的預(yù)測,并將檢測到的正選擇位點用PyMOL軟件(http://pymol.sourceforge.net/)定位到了其空間結(jié)構(gòu)中. 使用uniprot網(wǎng)站(http://www.uniprot.org/)查閱獲取更多功能域的信息.

2 結(jié)果與討論

2.1 選擇壓力分析

首先,運用位點模型對33個海洋哺乳動物數(shù)據(jù)集進行選擇壓力分析,LRTs結(jié)果顯示在SLC14A2基因中正選擇模型M8顯著優(yōu)于中性模型M7和M8a(P<0.001),且ω=1.955,提示這個基因顯著受正選擇作用. 然后,在M8模型中,使用BEB方法在SLC14A2基因中鑒定到了9個顯著受到正選擇的位點(35、462、463、467、468、469、470、512、741)且后驗概率均大于0.8. 位點模型中的M1a vs. M2a較M8 vs. M8a模型所得到的結(jié)果更為穩(wěn)定,因此又使用前者進一步檢測,結(jié)果表明M2a(ω=6.517)模型明顯優(yōu)于M1a模型(P<0.001),并鑒定了2個正選擇的位點(35、463),這兩個位點同時也被M8模型鑒定到(表1). 另外,進一步使用DataMonkey中的FEL模型鑒定了4個正選擇位點(35、90、276、735). 綜合以上3種ML方法,共有 2個位點(35、463)至少被兩種以上ML方法鑒定到,認為是強烈的正選擇位點(表1).

其次,為了檢測SLC14A2基因受到的正選擇是否局限于特定的進化譜系,在所有哺乳動物譜系73個物種中,LRT檢測發(fā)現(xiàn)free ratio模型顯著優(yōu)于one ratio模型(P<0.001),且在鯨類的多個支系中檢測到顯著的正選擇信號,即須鯨類的共同祖先支、瓶鼻海豚和印太瓶鼻海豚最近共同祖先以及鼠海豚和小頭鼠海豚共同祖先(ω>1,P<0.001),而在其他海洋哺乳動物類群以及陸生哺乳動物中均沒有鑒定出顯著的正選擇. 另外,以32種海洋哺乳動物作為前景支(ω=0.17211),其他陸生哺乳動物為背景支(ω=0.13012),運用two ratio模型檢測發(fā)現(xiàn)SLC14A2基因在海洋哺乳動物支系具有顯著的加速進化(P<0.001)(圖1).

由于分支模型和位點模型并不能同時考慮位點和支系發(fā)生的正選擇作用,因此進一步采用了更為嚴格的支位點模型檢測選擇壓力. 結(jié)果在小須鯨鑒定到顯著的正選擇,且鑒定了2個顯著的正選擇位點(592、598,表1). 綜合M8、M2a、FEL以及支位點模型的方法,共有4個正選擇位點(35、463、592、598)被認為是強烈的正選擇位點.

表1 在海洋哺乳動物中用兩種ML方法檢測到的正選擇位點Table 1 Positively selected sites detected using two maximum likelihood(ML)methods across cetacean phylogeny

2.2 海洋哺乳動物SLC14A2基因特異氨基酸突變位點和平行/趨同進化分析

為了探究海洋哺乳動物適應(yīng)高滲的海洋環(huán)境是否進化出了特異位點,對SLC14A2基因進行掃描后發(fā)現(xiàn)海洋哺乳動物具有11個特異的氨基酸突變位點. 其中,鯨類具有5個特異突變位點(Y14H、T398A、L455H、Q485R、A707S);鰭足類具有6個特異突變位點(G64T、S99G、T111S、K518M、A527T、S529H,圖1),其中第398、707位氨基酸位于尿素轉(zhuǎn)運蛋白的跨膜結(jié)構(gòu)域中,提示海洋哺乳動物可能通過尿素轉(zhuǎn)運蛋白特異的氨基酸突變,提高機體在高滲環(huán)境下的適應(yīng)能力.

進一步檢測SLC14A2基因在海洋哺乳動物類群間是否存在趨同的分子機制,重建了該基因的祖先序列,并進行非同義平行/趨同氨基酸位點檢測,鑒定了6個平行替換位點在海洋哺乳動物類群間表現(xiàn)出顯著的趨同進化信號. 如圖1所示,在鯨類譜系中鑒定了2個平行替換位點(L20S、R532S),如在20位氨基酸位點處具有一個平行氨基酸替代,中華白海豚的最近共同祖先支(n)與白鯨和一角鯨最近共同祖先支(v)均為絲氨酸(S),而其他哺乳動物支系均為亮氨酸(L). 在鯨類和鰭足類兩個海洋哺乳動物類群中共鑒定了4個平行替換位點(H1P、A219V、L508T、H814I),例如在第1位氨基酸位點處具有一個平行氨基酸替代,本研究中使用的鰭足類所有物種和抹香鯨(Physetercatodon)支系均為P,而其他哺乳動物物種此位點為H;在219位點處,小須鯨(Balaenopteraacutorostrata)和夏威夷僧海豹(Neomonachusschauinslandi)鑒定了1個平行氨基酸替代(A219V). 這些位點經(jīng)過JTT-fgen氨基酸替換模型統(tǒng)計檢驗,觀測值大于期望值,表明檢測到的位點是可靠的平行/趨同氨基酸替代位點.

不同海洋哺乳動物類群之間鑒定到的平行/趨同氨基酸提示為了適應(yīng)共同的高滲環(huán)境海洋哺乳動物可能進化出了趨同的尿素轉(zhuǎn)運機制,其中第219和814位氨基酸位點位于尿素轉(zhuǎn)運蛋白重要的跨膜結(jié)構(gòu)域,因此海洋哺乳動物可能通過跨膜結(jié)構(gòu)域的改變,增強了對尿素的轉(zhuǎn)運能力,從而對機體的水鹽平衡進行調(diào)控.

Free-ratio和支位點模型鑒定的正選擇分別用紅線和紅五星標注. 右側(cè)為SLC14A2基因在海洋哺乳動物中檢測到的特異氨基酸位點(黃色)和平行/趨同位點(藍色). 字母表示祖先枝. 作為外群的陸生哺乳動物包括鯨偶蹄目、食肉目、嚙齒目等.圖1 Free-ratio和支位點模型鑒定的海洋哺乳動物SLC14A2基因的正選擇Fig.1 Evidence of positive selection across the phylogeny of marine mammals identified by free-ratio and branch-site models in SLC14A2 gene

強烈正選擇位點和位于跨膜結(jié)構(gòu)域的特異氨基酸位點分別用紫色和藍色圓球表示.圖2 SLC14A2蛋白質(zhì)三維結(jié)構(gòu)Fig.2 The three-dimensional structure of the SLC14A2 protein

2.3 SLC14A2基因三維結(jié)構(gòu)分析

印太瓶鼻海豚(T.aduncus)的SLC14A2基因核苷酸序列在I-TASSER網(wǎng)站上(http://zhanglab.ccmb.med.umich.edu/I-TASSER/)進行了蛋白三維結(jié)構(gòu)的預(yù)測,并運用PyMOL軟件(http://pymol.sourceforge.net/)將強烈正選擇位點和重要的特異氨基酸突變定位到3D結(jié)構(gòu)中(圖2). 這些位點大多位于重要的結(jié)構(gòu)域或結(jié)構(gòu)域附近,如第463位正選擇位點位于蛋白激酶C的磷酸化位點S490附近;第35位正選擇位點位于蛋白的極性氨基酸(polar amino acid)區(qū)域,此區(qū)域使得此段蛋白穩(wěn)定位于膜外,同時此位點也與蛋白激酶C的磷酸化位點S26鄰近[28],提示正選擇位點可能與蛋白尿素轉(zhuǎn)運功能相關(guān). 另外,第398、707位特異氨基酸位點存在于蛋白重要的跨膜結(jié)構(gòu)域. 跨膜結(jié)構(gòu)域的改變通常與通道蛋白開放時間和轉(zhuǎn)運效率密切相聯(lián)系[29-30]. 因此海洋哺乳動物可能通過跨膜結(jié)構(gòu)域的改變,增強了對尿素的轉(zhuǎn)運能力,從而對機體的水鹽平衡進行調(diào)控.

3 結(jié)論

趨同進化是進化生物學(xué)中的核心科學(xué)問題. 海洋哺乳動物不同類群是獨立起源的,由不同的陸生祖先演化而來,親緣關(guān)系相距甚遠. 但是,為了適應(yīng)高滲的海洋環(huán)境,海洋哺乳動物產(chǎn)生了趨同的滲透調(diào)節(jié)適應(yīng)機制,具有很強的排鹽保水能力,且能夠產(chǎn)生高度濃縮的尿液. 先前的解剖學(xué)證據(jù)表明海洋哺乳動物具有復(fù)腎結(jié)構(gòu),通過增加腎臟表面積以提高腎臟的排泄能力[11-12],為海洋哺乳動物形成高濃度的濃縮尿提供了前提. 本研究發(fā)現(xiàn)尿素轉(zhuǎn)運基因——SLC14A2在海洋哺乳動物發(fā)生了顯著的加速進化(P<0.001),且有2個正選擇位點至少被兩種ML方法檢測到,這兩個氨基酸位點都位于蛋白激酶C的磷酸化位點附近,推測與蛋白尿素轉(zhuǎn)運功能相關(guān). 因為高滲透的內(nèi)環(huán)境主要通過蛋白激酶C介導(dǎo)增加 UT-A1 和UT-A1A3的磷酸化,從而增強尿素的轉(zhuǎn)運,且已有研究證實敲除SLC14A2基因則出現(xiàn)不能濃縮尿液的生理缺陷[19]. 因此本研究在海洋哺乳動物中鑒定的加速進化和正選擇,提示可能增強了尿素轉(zhuǎn)運能力從而維持了極高的滲透壓以濃縮尿液,為海洋哺乳動物產(chǎn)生高濃縮尿液提供了分子證據(jù). 另外,本研究通過對73個物種的SLC14A2基因進行比較,在鯨類和鰭足類中共鑒定出11個特異的氨基酸位點,其中2個氨基酸位于蛋白質(zhì)重要的功能結(jié)構(gòu)域,提示海洋哺乳動物可能通過這些氨基酸突變使尿素轉(zhuǎn)運蛋白功能發(fā)生適應(yīng)性的改變,從而增強尿素轉(zhuǎn)運能力以適應(yīng)高滲的海洋環(huán)境. 特別重要的是,在海洋哺乳動物三個類群間鑒定了6個平行/趨同進化位點,且其中兩個位點位于蛋白質(zhì)的跨膜結(jié)構(gòu)域中,有研究表明跨膜結(jié)構(gòu)域中單個氨基酸的改變即可改變通道蛋白的動力學(xué)和滲透特性[29-30]. 因此,推測海洋哺乳動物可能通過改變跨膜結(jié)構(gòu)域的三維結(jié)構(gòu)進而影響通道蛋白打開的時間和尿素轉(zhuǎn)運的效率,進一步為海洋哺乳動物尿素轉(zhuǎn)運蛋白在滲透調(diào)節(jié)中發(fā)揮重要功能以產(chǎn)生高度濃縮的尿液提供趨同的分子證據(jù).

特別需要指出的是,用free-ratio模型和更嚴格的支位點模型鑒定的正選擇均發(fā)生在鯨類譜系,可能與鯨類完全營水生生活面臨較強的滲透調(diào)節(jié)壓力有關(guān). 之前生理學(xué)研究表明,鯨類尿液中尿素濃度比其近緣的陸生哺乳動物——牛高4倍[14],說明鯨類通過排高度濃縮的尿液以調(diào)節(jié)水鹽平衡. 鯨類SLC14A2顯著的正選擇進一步證實鯨類通過增強尿素轉(zhuǎn)運以產(chǎn)生高度濃縮的尿液. 另外,比較基因組和轉(zhuǎn)錄組學(xué)也在長江淡水江豚和海生江豚中分別鑒定了與尿素轉(zhuǎn)運蛋白SLC14A2、ACE相關(guān)基因的特異的SNP,且發(fā)現(xiàn)兩類江豚存在SLC14A2基因的表達差異,進一步為海生江豚產(chǎn)生高濃縮尿液(東亞江豚尿液:(74.96±20.06)mmol/L;長江江豚(40.7±15.78)mmol/L)提供了分子證據(jù)[17]. 先前研究已經(jīng)證實水通道蛋白AQP2在鯨類中顯著受正選擇,AQP2基因的表達定位于腎臟的集合管,通過主動運輸對水分進行重吸收,該基因顯著的正選擇提示鯨類可能通過增強水通道蛋白的功能而增強集合管對水的重吸收的能力,進一步為產(chǎn)生高度濃縮尿液提供了分子證據(jù). 滲透調(diào)節(jié)是一個綜合調(diào)控的過程,腎素—血管緊張素—醛固酮系統(tǒng)(RAAS)在哺乳動物水鹽代謝調(diào)節(jié)中起重要作用. 該系統(tǒng)中關(guān)鍵基因血管緊張素轉(zhuǎn)化酶ACE和血管緊張素原AGT在多個研究中均證實鯨類譜系顯著受正選擇[16-17],表明鯨類可能還通過RAAS調(diào)節(jié)系統(tǒng)更高效地維持鯨類的水鹽平衡.

綜上,本研究通過對尿素轉(zhuǎn)運蛋白編碼基因SLC14A2在海洋哺乳動物中的進化分析,揭示海洋哺乳動物SLC14A2基因的進化與滲透調(diào)節(jié)機制的關(guān)系. 分析結(jié)果表明海洋哺乳動物類群間SLC14A2基因的加速進化、正選擇位點以及平行/趨同位點的鑒定,均為海洋哺乳動物趨同的產(chǎn)生高度濃縮尿液以保水提供了分子證據(jù). 另外,正選擇主要集中在鯨類進化譜系,可能與鯨類完全營水生生活面臨更大的滲透調(diào)節(jié)的壓力有關(guān). 當然,SLC14A2基因在海洋哺乳動物中滲透調(diào)節(jié)的具體功能,將來還需要進一步的功能驗證,以揭示該基因與滲透調(diào)節(jié)的關(guān)系.