《International Journal of Molecular Sciences》：Comparative Evolution of SLC12 Transporters in Tilapias and Salinity-Responsive Expression in Nile Tilapia
编辑推荐：
罗非鱼物种在盐度耐受性上差异显著，为研究渗透适应的分子基础提供了有用体系。溶质载体家族12（SLC12）编码阳离子-氯离子共转运蛋白，在离子稳态和渗透调节中发挥核心作用。本研究利用比较基因组、系统发育、蛋白质结构和染色体分析，比较了尼罗罗非鱼（Oreochro
罗非鱼物种在盐度耐受性上差异显著，为研究渗透适应的分子基础提供了有用体系。溶质载体家族12（SLC12）编码阳离子-氯离子共转运蛋白，在离子稳态和渗透调节中发挥核心作用。本研究利用比较基因组、系统发育、蛋白质结构和染色体分析，比较了尼罗罗非鱼（Oreochromis niloticus）、莫桑比克罗非鱼（O. mossambicus）和蓝罗非鱼（O. aureus）的SLC12基因家族。另外，研究人员重新分析了已发表的来自GIFT尼罗罗非鱼在0、17或27 ppt下维持3.5个月的长期盐度RNA-seq数据集，重点关注鳃、肠和皮肤中SLC12的表达。三种罗非鱼均含有17个SLC12基因座，亚家族拷贝数相同，染色体分布大致保守。在尼罗罗非鱼转录组数据集中，旁系同源的SLC12A2、SLC12A7和SLC12A10拷贝表现出不同的组织特异性和盐度依赖性表达模式。SLC12A8在尼罗罗非鱼中也表现出盐度响应性表达，促使进一步比较序列和结构分析。三个SLC12A8直系同源物在序列和预测结构上大体相似，但观察到若干物种特异性残基差异。这些发现揭示了罗非鱼中SLC12家族的整体保守性，识别了尼罗罗非鱼旁系同源成员间的盐度响应性表达差异以及SLC12A8的种间序列变异，并为进一步功能研究提供了候选特征。
论文解读：罗非鱼SLC12家族的比较进化与盐度响应
环境盐度直接决定鱼类体液与外界环境之间的渗透梯度，维持体液渗透压和无机离子稳态是鱼类适应不同水生环境的基本生理要求。淡水硬骨鱼面临持续的离子丢失和水渗入，而海水或高盐环境中的鱼类则需应对失水以及Na
+和Cl
?的被动内流。鳃、肠和肾通过协同调节离子吸收、离子排泄和水转运来维持体内稳态，其中鳃上皮的离子细胞是跨上皮Na
+和Cl
?转运的主要场所。溶质载体家族12（SLC12）编码电中性阳离子-氯离子共转运蛋白，是维持细胞内、外Na
+、K
+和Cl
?稳态、细胞体积以及跨上皮离子转运的重要膜转运系统。鱼类SLC12家族因硬骨鱼特异性基因组复制而更加复杂，复制拷贝的保留与丢失导致多个亚家族存在旁系同源基因，不同拷贝可能在组织表达、调控响应和生理功能上发生分化。罗非鱼属（Oreochromis）物种间盐度耐受性差异显著，莫桑比克罗非鱼（O. mossambicus）具有突出的广盐性和高盐耐受能力，蓝罗非鱼（O. aureus）也有较高耐盐性，而尼罗罗非鱼（O. niloticus）的高盐耐受性相对有限，具有较大遗传改良空间。以往研究多集中于单一物种、QTL候选区或经典离子调节通路，尚缺乏对SLC12家族中功能了解较少的成员及复制产生的旁系同源基因在蛋白质序列、膜拓扑和环境响应性方面分化的系统比较。因此，该研究以三种罗非鱼为比较体系，在基因组水平上系统比较SLC12家族，并结合已发表的尼罗罗非鱼长期盐度转录组数据，为理解盐度适应的分子机制提供候选特征。该论文发表在《International Journal of Molecular Sciences》。
研究人员首先利用比较基因组、系统发育、蛋白质结构和染色体分析，比较了尼罗罗非鱼、莫桑比克罗非鱼和蓝罗非鱼的SLC12基因家族；随后重新分析了已发表的GIFT尼罗罗非鱼长期盐度RNA-seq数据，表征SLC12成员在鳃、肠和皮肤中的组织及盐度依赖性表达；最后对SLC12A8进行了序列、跨膜拓扑和三维结构比较。结果表明，三种罗非鱼的SLC12家族组成和染色体分布总体保守，但旁系同源拷贝间存在明显的表达分化，且单拷贝SLC12A8存在种间序列变异。这些发现揭示了SLC12家族“整体保守、局部分化”的特征，为后续功能研究提供了候选位点。
关键技术方法方面，研究人员从NCBI、Ensembl和Figshare获取三种罗非鱼的染色体水平参考基因组及注释文件，采用TBLASTN和基于Pfam结构域谱的HMMER搜索鉴定候选SLC12基因，并以系统发育聚类最终确认。表达分析重新利用了已发表的GIFT尼罗罗非鱼长期盐度RNA-seq数据，来自DDBJ BioProject PRJDB9175和PRJDB10880，共36个文库，涵盖鳃、肠和皮肤，每种组织每个盐度4个生物学重复；数据经STAR比对、featureCounts定量和DESeq2进行差异表达分析。此外还使用MEME、Pfam、TMHMM、DeepTMHMM、Phobius、KaKs_Calculator和SWISS-MODEL等方法完成基序、结构域、跨膜拓扑、分子进化和三维结构比较。
2.1 三种罗非鱼SLC12家族成员的鉴定与蛋白质特征
三种罗非鱼各鉴定到17个SLC12家族成员，共51个SLC12基因。三个物种家族组成一致：SLC12A1、SLC12A3、SLC12A4、SLC12A6和SLC12A8各为单拷贝；SLC12A2、SLC12A5和SLC12A7各有两个拷贝；SLC12A9和SLC12A10各有三个拷贝。51个蛋白长度为757–1212个氨基酸，分子量82.55–133.66 kDa，理论等电点（pI）为5.62–8.16，不稳定指数为34.17–44.43，总平均亲水性（GRAVY）为?0.0342–0.4511。TMHMM预测所有蛋白均含多个跨膜螺旋，多数为10–12个。
2.2 SLC12家族的系统发育、保守基序、结构域和基因结构
系统发育分析将SLC12蛋白归入SLC12A1、A2、A3、A4、A5、A6、A7、A8、A9和A10亚家族，三种罗非鱼的直系同源蛋白分别聚在对应分支。同一亚家族或直系同源组内保守基序排列相似，而亚家族之间基序数量和组成存在差异。51个蛋白均含AA_permease结构域，其中48个还含SLC12相关结构域，21个含AA_permease_N结构域。三个SLC12A8蛋白仅检测到AA_permease结构域，未检测到SLC12相关和AA_permease_N结构域。外显子数量、内含子长度和基因组跨度在不同成员间存在差异。
2.3 SLC12基因在三种罗非鱼中的染色体分布
染色体定位显示，三种罗非鱼中17个SLC12基因座均分布在11条染色体或支架组上。尼罗罗非鱼和蓝罗非鱼的分布模式高度相似，莫桑比克罗非鱼对应直系同源基因座的相对分布也总体一致。尼罗罗非鱼的SLC12A7a和SLC12A9c位于LG18，该连锁群包含前期研究鉴定的盐度耐受相关QTL簇；SLC12A8位于LG16。
2.4 长期盐度处理下SLC12基因的组织表达模式
在GIFT尼罗罗非鱼长期盐度处理数据中，SLC12成员在鳃、肠和皮肤中呈现明显的组织和盐度依赖性表达。SLC12A10a在鳃中淡水条件下表达最高，高盐时显著降低；SLC12A10b整体表达水平低且随盐度变化较小。SLC12A2a在鳃和皮肤中随盐度升高总体下降，而鳃中SLC12A2b在17 ppt时较高、27 ppt时下降。KCC类成员SLC12A4、SLC12A6、SLC12A7a和SLC12A7b在不同组织和盐度中表现出不同变化。SLC12A8在肠中表达高于鳃和皮肤，鳃中随盐度升高逐渐下降，皮肤中从淡水到17 ppt下降、27 ppt略回升。SLC12A9a在鳃和皮肤中随盐度升高下降，肠中27 ppt时回升；SLC12A9b在鳃中表达较高。SLC12A9c因表达过低未在主要表达图中显示，其鳃中17 ppt相对淡水有升高趋势，但多重检验校正后不显著。
2.5 三种罗非鱼17个SLC12直系同源基因座的Ka、Ks和Ka/Ks比较
与5种外部参考鱼类的同源序列比较显示，Ka/Ks范围为0.023–0.279，均值为0.120，全部小于1，表明SLC12家族蛋白编码序列整体受较强的纯化选择。SLC12A4的平均Ka/Ks最低，为0.043；SLC12A5a、SLC12A5b和SLC12A6也较低。SLC12A8的平均Ka/Ks为0.169，SLC12A9a/b/c分别为0.152、0.174和0.155，SLC12A10a/b/c分别为0.163、0.128和0.151，虽相对较高但仍低于1。这些比较仅作描述性说明，未被用于推断进化速率的显著差异。
2.6 三种罗非鱼SLC12A8序列与三维结构比较
SLC12A8多序列比对显示三种蛋白整体相似。OnSLC12A8第443位为苏氨酸（Thr），而OmSLC12A8和OaSLC12A8对应位置分别为异亮氨酸Ile453和Ile449。SWISS-MODEL使用同一模板A0A665XEB9构建三个模型，GMQE约为0.66–0.67，与模板序列一致性约为79%。全局叠合得到RMSD约为0.37–0.52 ?，但三种模型来自同一模板，因此该值仅作为建模框架内的相似性描述。局部结构比较显示，OnSLC12A8的T443极性侧链靠近S439主链羰基，可能形成与含异亮氨酸位点不同的局部理化环境，但研究未据此推断特定氢键或原子相互作用。
2.7 三种罗非鱼SLC12A8跨膜拓扑比较
TMHMM预测OnSLC12A8含11个跨膜螺旋，而OaSLC12A8和OmSLC12A8含12个；DeepTMHMM和Phobius均预测三个蛋白含12个跨膜螺旋。因此11与12的差异依赖预测方法，不构成一致的种间跨膜拓扑差异。三种方法都将T443/I449/I453同源位点定位在对应跨膜螺旋C端边界附近，但精确边界随方法变化，说明该位点位于预测敏感的局部膜边界区域。
讨论部分指出，三种罗非鱼SLC12家族组成和染色体位置高度保守，说明盐度耐受差异并非源于近期SLC12成员的大规模获得或丢失，而可能涉及更细微的基因调控或蛋白功能差异。旁系同源拷贝间的表达差异提示复制后的调控亚功能化，即在长期盐度适应中，不同拷贝可在不同组织和盐度范围内被差异化募集。SLC12A8作为单拷贝基因整体受纯化选择，但存在种间氨基酸差异；其肠中较高表达和盐度响应性使其成为合适的候选功能研究对象，但现有证据不足以将其直接归类为典型盐离子转运蛋白，T443/I449/I453变异也不能直接与种间盐度耐受差异关联。局部异亮氨酸到苏氨酸的替换改变了跨膜边界附近的残基极性，可能影响局部结构环境，但需要定点突变和功能实验验证。
研究结论可翻译为：总之
