齐鲁工业大学学报   2021, Vol. 35 Issue (5): 31-35
0
喜马拉雅旱獭线粒体基因组与低氧适应研究进展[PDF全文]
南新营, 李优, 李耀东     
青海大学, 西宁 810016
摘要:喜马拉雅旱獭(Marmota himalayana)作为世居高原的一个特有物种,在长期的进化过程中,已获得了一整套完善的、严密的对于低氧环境的适应机制。线粒体作为动物的"能量工厂",通过氧化磷酸化产生能量,在氧的利用和能量代谢中发挥重要作用,成为学者研究生物适应高原低氧环境机制的首选。本文就喜马拉雅旱獭线粒体基因组与低氧适应进行综述,以增进对该物种应对低氧机理的了解,为进一步研究提供参考借鉴。
关键词喜马拉雅旱獭    啮齿目    低氧适应    线粒体基因组    
Advances in Mitochondrial Genome and Hypoxic Adaptation of Marmota Himalayana
NAN Xin-ying, LI You, LI Yao-dong     
Qinghai University, Xining 810016, China
Abstract: Marmot Himalayan, as a unique species living on the plateau, has acquired a complete set of strict adaptation mechanism to the hypoxic environment in the long course of evolution.As the "energy factory" of animals, mitochondria produce energy through oxidative phosphorylation and play an important role in oxygen utilization and energy metabolism, which has become the first choice for researchers to study the mechanism of organisms adapting to the low oxygen environment at high altitude.In this paper, the mitochondrial genome and hypoxia adaptation of Marmota himalayana were reviewed to enhance the understanding of the mechanism of hypoxia in this species and provide reference for further research.
Key words: marmota himalayana    rodentia    hypoxia adaptation    mitochondrial genome    

高原低氧适应是指世居高原的人或者动物为能在高原环境中生存, 经长期自然选择所产生的一种非可逆但能稳定遗传的形态结构及生理生化等方面的改变, 是机体为适应高原低氧环境而形成的一种保护机制。青藏高原是世界上最高的高原, 有“世界屋脊”之称, 平均海拔在4 000米以上, 总面积约240万平方公里, 生活着大约210种野生哺乳动物, 占全国总种数的50%左右, 为人类提供了高原自然界的原始“本底”, 为开展有关青藏高原生物学、高原医学、动植物遗传多样性及保护等学科的研究, 提供了理想的基地和天然实验室。目前学者对于高原低氧适应的研究多集中于藏羊[1]、藏猪[2]、藏鸡[3]、牦牛[4]、高原鼠兔[5]、蜥蜴[6]等物种。

1 旱獭的分类与分布

旱獭, 又名土拨鼠、草地獭、哈拉、雪猪、曲娃(藏语)等, 属哺乳纲(Mammalia), 啮齿目(Rodentia), 松鼠科(Sciuridae), 旱獭属(Marmota), 其体形肥大, 是松鼠科中体型最大的一种, 成年旱獭体重可达3.9~9 kg, 颈部粗短, 耳壳短小, 四肢短粗, 尾短而扁平, 体背棕黄色, 是陆生和穴居的草食性、冬眠性野生动物, 主要栖息于山地的阳坡、阶地、谷地、山麓平原等环境, 过家族生活, 个体接触密切[7]

目前世界上已知的旱獭属动物有14种, 我国有4种, 分别是蒙古旱獭、长尾旱獭、喜马拉雅旱獭、阿尔泰旱獭[8]。其中喜马拉雅旱獭是分布范围最广、数量最多的一种, 主要栖息于海拔2 500~5 000 m范围的高山草地和草原, 广泛分布于我国的青藏高原及毗邻的印度、尼泊尔及巴基斯坦等地的各类高寒草甸草原[9], 寄居洞穴可深达十数米, 在这样极低的氧分压下, 依然有很强的生理代谢能力和生殖能力。

2 线粒体基因组与低氧适应研究

线粒体基因组具有独特的结构特点, 存在于几乎所有动物的各组织细胞内, 由于其基因组较小、基因结构及排列简单、多拷贝、易于分离、严格的母系遗传、mtDNA序列在同一个体内为单倍型、没有种内重组、无内含子及非编码序列少(D-loop环除外), 已成为生物进化研究的非常有用的分子标识和最适遗传标志[10]。在1981年Anderson等测定出人类的线粒体基因组全序列以来[11], 已有十余种哺乳动物的线粒体基因组全序列被测定出来, 用以研究这些动物的遗传结构和起源分化[12]。Chao等[13]通过对喜马拉雅旱獭线粒体基因组测序, 已获得喜马拉雅旱獭线粒体基因组全序列(序列号: JX069958), 该报告显示喜马拉雅旱獭线粒体基因组长度为16 443bp, 由13个蛋白编码基因、2个核糖体RNA(rRNA)基因、22个转移RNA(tRNA)基因和一个典型的控制区组成, 13个编码蛋白质的功能基因包括7个NADH脱氢酶复合体基因(ND1-ND6, ND4L), 1个细胞色素b基因(cytB), 3个细胞色素氧化酶复合体基因(COX1-COX3), 2个ATP合成酶复合体基因(ATPase6、ATPase8)。这是松鼠科第一个完整线粒体基因组的报告, 为该物种的遗传功能及进化分析提供了理论依据。

过去学者关于喜马拉雅旱獭对低氧适应的研究多集中在生理、生化、形态方面, 贾荣莉[14]认为喜马拉雅旱獭的红细胞数目及心肌细胞的形态结构、生理功能与生化酶均未受到高原低氧的不良影响, 表明其对高原低氧有很好的适应性。马嵋等[15]研究显示, 新疆灰旱獭为适应高原低氧环境, 与SD大鼠相比较, 其心肌细胞明显肥大, 体重指数、右心室肥厚指数及纤维化程度也明显增高, 右心室发生重构。马英等[16]对喜马拉雅旱獭DNA控制区研究表明地理隔离是影响其种群动态变化的主要因素。何建等[17]实验表明喜马拉雅旱獭的平均红细胞体积(MCV)随海拔增高而降低, 血小板(PLT)随海拔增高而增加, 且不同海拔地区的心肺组织切片均有不同程度的损伤。

线粒体能量代谢在生物进化和物种形成过程中起着重要作用。早在1992年Wallace[18]提出人类mtDNA编码的13个氧化磷酸化必需基因的变异可能是人类适应缺氧的遗传因素之一。为了验证这一说法, 1994年Torroni[19]采用两种方法来确定高海拔地区人群mtDNA变异的性质是否与他们对高海拔的适应有关。第一种是比较生活在不同海拔(4 500 m, 3 600 m, 3 000 m)的藏族人mtDNA变异, 在此比较中, 单倍型和单倍群在海拔高度上均无显著差异, 此外在多个受试者单倍型(AS62, AS65, AS131, AS143, ASA53)在不同海拔高度的藏族个体中进行观察, 结果mtDNA单倍型与不同海拔之间没有联系, 当然这并不排除特定的mtDNA基因型有助于适应高海拔的可能性, 因为所有的研究样本都来自于相对较高的海拔(3 000 m以上);第二种是通过比较藏族人和其他亚洲人mtDNA单倍型和单倍群的频数分析表明在高海拔地区任何单倍群均未出现大的选择压力。这两种方法证实mtDNA突变不太可能在藏族人适应高海拔的过程中起主要作用。

NADH脱氢酶又称复合体Ⅰ或NADH-泛醌还原酶, 是线粒体呼吸链中第一个也是最大的酶复合物[20], 它接收从NADH氧化得到的电子并转移给泛醌, 发生在该复合物中的残基取代可干扰质子泵传送效率从而影响代谢表现[21]。Shi等[22]研究表明mt-ND1和mt-ND2基因的遗传多样性与高原适应相关, 可作为与高海拔适应相关的候选基因。Zhou等[23]研究发现ND2、ND4、ATPase6在鸟类适应高海拔过程中有积极选择的重要证据。Bao等[24]在藏鸡品种中发现了NADH脱氢酶mt-ND5亚基的错义突变, 通过测定不同基因型藏鸡的线粒体呼吸控制率(RCR)和NADH脱氢酶活性, 发现缺氧条件下藏鸡ND5基因突变与线粒体呼吸功能密切相关, 认为ND5是藏鸡适应低氧环境的候选基因。LiYu等[25]对2只中国高海拔地区(3 500 m)和其他8只越南东部低海拔地区(750 m)塌鼻猴线粒体基因组进行测序, 在中国塌鼻猴ND2和ND6中分别发现了2个重要的适应性突变位点(533和3 307), 认为可能是中国塌鼻猴适应低氧环境的突变位点。Xu等[26]对不同海拔(4 500 m, 3 300 m, 3 100 m)3匹藏马线粒体基因组进行测序分析, 发现ND6基因在3匹马中都具有较高的非同义突变率, 表明ND6基因可能在西藏马的高海拔适应过程中起重要作用。Dehean等[27]在人神经胶质瘤细胞中发现, 编码线粒体ND6基因的T14634C发生突变, 预测这种突变会破坏ND6蛋白的构象, 从而使细胞对低氧敏感丧失应答能力。Luo等[28]在高原鼠兔的低氧适应研究中发现在3只鼠兔的mtDNA编码蛋白中, 有186个基因内氨基酸变异, 鉴定出15种鼠兔专有氨基酸变异, 且这些氨基酸在11种非鼠兔哺乳动物和人类的系统发育分析中都是100%保守的, 受影响的鼠兔基因包括cytB、COX1、COX2、ND2、ND3、ND5和ND6, 推测这些氨基酸位点的变异可能是鼠兔适应低温低氧环境的遗传基础。

细胞色素c氧化酶(Cytochrome c oxidase, COX)即复合物Ⅳ, 是呼吸链中的末端复合物, 是氧化磷酸化过程中氧获取和利用的唯一部位。真核细胞中COX由9~13个亚基多肽组成, 三个线粒体编码的亚基, 即COX1-3构成催化核心, 它包含4个氧化还原中心, 即COX1中的CuB、Hema、Hema3和COX2中的CuA, 并催化电子从细胞色素c转化为分子氧, 其他亚基是核编码的, 在组装和稳定全酶和调节COX活性方面具有重要的作用[29]。COX催化氧的末端还原, 可能是高海拔适应过程中一个特别重要的选择目标。Zhang[30]等对青藏高原和华北地区迁徙蝗虫种群的线粒体做缺氧处理, 发现只有在西藏蝗虫中COX表现出明显较高的活性, 且COX活性的增强不是由蛋白质引起的, 而是由COX的铁细胞色素亲和力的增加和通过催化氧化还原中心的电子传输效率的增加而导致催化效率提高, 表明COX的功能调节是西藏蝗虫适应高原缺氧的基础。Scott等[31]为避免后天环境因素的干扰, 将斑头雁和2个亲缘关系较远的品种在低海拔地区统一饲养, 并对3个品种代谢途径中关键酶的活性进行测定, 发现斑头雁的COX功能有很大差异, 斑头雁的COX最大酶活性显著低于低海拔品种, 继而对COX1-4亚基蛋白和mRNA进行定量比较, 发现3个品种之间的4个基因表达量没有明显差异, 排除了酶表达量变化对COX酶活性的影响因素, 又对COX的4个亚基基因进行测序, 结果一个非同义突变位点在斑头雁的COX3序列上被发现, 导致COX3的第116位氨基酸残基由Trp突变成了Arg, 这种突变在脊椎动物中属于保守位点, 通过结构模型预测该突变可能改变了COX1和COX3亚基之间的相互作用, 降低了氧化磷酸化对细胞内缺氧的敏感性, 从而适应极度缺氧环境。

ATP合酶(ATPsynthase)即复合物Ⅴ, 由F0和F1两个功能域组成, 与复合物Ⅰ-Ⅳ一起成为氧化磷酸化过程中ATP产生的基础[32]。李耀东等[33]通过比较喜马拉雅旱獭与平原地区同种动物松鼠线粒体ATPase6基因表达的差异, 认为旱獭组ATPase6基因的表达高于松鼠组, 且旱獭组大都有ATPase6基因片段点突变。Zhao等[34]在藏鸡的ATPase6基因中检测到6个SNPs, 发现1个错义突变(9441G>A), 导致了氨基酸的替换, 推测该突变可能在藏鸡能量转换和代谢中起重要作用。李娜等[35]研究发现不同物种间高海拔地区西藏沙蜥线粒体基因组中蛋白质编码基因进化最快, 不同基因间ATPase8进化速率最快, 且在ATPase8基因上存在两个正选择位点, 因此认为ATPase8可能在西藏沙蜥高海拔适应中发挥重要作用。

细胞色素b(cytochrome b) 是复合物Ⅲ的一个催化活性亚基, 两端有2个CoQ结合位点, 其中一个是位于细胞质侧的QH2氧化位点Q0, 与低电位的细胞色素bL有关, 另一个是基质侧的Q还原位点Qi, 与高电位细胞色素bH有关[36]。目前未见有关于细胞色素b与低氧适应相关的报道。

3 展望

近年来, 高原世居动物的低氧适应分子机制已成为高原医学研究的热点, 目前国内外关于动物线粒体基因组对高原低氧适应的研究已取得初步成果, 但仍未全面解析高原动物线粒体基因组适应低氧环境的分子机制。喜马拉雅旱獭能良好地适应高原低氧、低温等气候, 是研究高原低氧适应机制的理想动物模型。高通量测序技术加速了线粒体基因组的测序, 揭示了真核生物组织结构、基因内容、复制和转录模式的多样性, 为学者从分子层面研究该机制提供了更为便捷有效的途径。学者应结合目前研究成果, 从DNA-RNA-蛋白质水平深入了解线粒体基因组在低氧适应过程中的重要作用。明确喜马拉雅旱獭低氧适应的分子机制不仅有助于该物种种群多样性的保护, 更为高原疾病的预防和治疗提供了新的方向和思路。

参考文献
[1]
MA L, SHAO X, WANG Y, et al. Molecular cloning, characterization and expression of myoglobin in tibetan antelope (pantholops hodgsonii), a species with hypoxic tolerance[J]. Gene, 2014, 533(2): 532-537. DOI:10.1016/j.gene.2013.09.030
[2]
ZHANG B, BAN D, GOU X, et al. Genome-wide DNA methylation profiles in tibetan and yorkshire pigs under high-altitude hypoxia[J]. J Anim Sci Biotechnol, 2019, 10: 25. DOI:10.1186/s40104-019-0316-y
[3]
ZHANG Y, ZHANG X, ZHANG Y, et al. Comparative transcriptomic and proteomic analyses provide insights into functional genes for hypoxic adaptation in embryos of tibetan chickens[J]. Sci Rep, 2020, 10(1): 11213. DOI:10.1038/s41598-020-68178-w
[4]
DING Y, LIU J, XU Y, et al. Evolutionary adaptation of aquaporin-4 in yak (Bos grunniens) brain to high-altitude hypoxia of qinghai-tibetan plateau[J]. High Alt Med Biol, 2020, 21(2): 167-175. DOI:10.1089/ham.2019.0076
[5]
LI J, YANG Q, BAI Z, et al. Chronic cold exposure results in subcutaneous adipose tissue browning and altered global metabolism in qinghai-tibetan plateau pika (Ochotona curzoniae)[J]. Biochem Biophys Res Commun, 2018, 500(2): 117-123. DOI:10.1016/j.bbrc.2018.03.147
[6]
CORDERO GA, ANDERSSON BA, SOUCHET J, et al. Physiological plasticity in lizard embryos exposed to high-altitude hypoxia[J]. Journal of Experimental Zoology, 2017, 327(7): 423-432. DOI:10.1002/jez.2115
[7]
何清. 生态与环境[M]. 石家庄: 河北科学技术出版社, 2013.
[8]
孟德荣. 我国旱獭的种类及国内分布研究[J]. 毛皮动物饲养, 1996(3): 18-19.
[9]
刘海青, 范微, 张静宵, 等. 喜马拉雅旱獭实验动物化的研究进展[J]. 中国比较医学杂志, 2015, 25(11): 64-68. DOI:10.3969/j.issn.1671-7856.2015.11.014
[10]
SCHAACK S, HO EKH, MACRAE F. Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome[J]. Philos Trans R Soc Lond B Biol Sci, 2020, 375(1790): 20190173. DOI:10.1098/rstb.2019.0173
[11]
ANDERSON S, A T BANKIER, B G BARRELL, et al. Sequence and organization of the human mitochondrial genome[J]. Nature, 1981, 290: 457-465. DOI:10.1038/290457a0
[12]
KAMALAKKANNAN R, BHAVANA K, PRABHU VR, et al. The complete mitochondrial genome of indian gaur, bos gaurus and its phylogenetic implications[J]. Sci Rep, 2020, 10(1): 11936. DOI:10.1038/s41598-020-68724-6
[13]
CHAO Q J, LI Y D, GENG X X, et al. Complete mitochondrial genome sequence of marmota himalayana (rodentia: sciuridae) and phylogenetic analysis within rodentia[J]. Genet Mol Res, 2014, 13(2): 2739-2751. DOI:10.4238/2014.April.14.3
[14]
贾荣莉. 低氧适应动物喜马拉雅旱獭[J]. 解剖科学进展, 1997(1): 75-77.
[15]
马嵋, 姜涛, 张春, 等. 高原低氧适应动物新疆灰旱獭右心室重构的组织学改变[J]. 中国实验动物学报, 2014, 22(2): 45-47. DOI:10.3969/j.issn.1005-4847.2014.02.010
[16]
马英, 李海龙. 喜马拉雅旱獭线粒体DNA控制区遗传多样性及系统发育[J]. 兽类学报, 2019, 39(3): 285-294.
[17]
何建, 陈芳, 李海龙, 等. 喜马拉雅旱獭对高原低氧环境的生理适应性研究[J]. 野生动物学报, 2020, 41(1): 63-69. DOI:10.3969/j.issn.1000-0127.2020.01.009
[18]
WALLACE D C. Mitochondrial DNA complex Ⅰ and Ⅲ mutations associated with Leber's hereditary optic neuropathy[J]. Genetics, 1992, 730: 163-173.
[19]
TORRONI A, MILLER J A, MOORE L G, et al. Mitochondrial DNA analysis in Tibet: implications for the origin of the tibetan population and its adaptation to high altitute[J]. Am J Phys Anthropol, 1994, 93(2): 189-199. DOI:10.1002/ajpa.1330930204
[20]
DA FONSECA RR, JOHNSON WE, O'BRIEN SJ, et al. The adaptive evolution of the mammalian mitochondrial genome[J]. BMC Genomics, 2008, 9: 119. DOI:10.1186/1471-2164-9-119
[21]
HASSANIN A, ROPIQUET A, COULOUX A, et al. Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe caprini (bovidae, antilopinae)[J]. J Mol Evol, 2009, 68(4): 293-310. DOI:10.1007/s00239-009-9208-7
[22]
SHI Y, HU Y S, WANG J, et al. Genetic diversities of MT-ND1 and MT-ND2 genes are associated with high-altitude adaptation in yak[J]. Mitochondrial DNA Part A, 2018, 29(3): 485-494. DOI:10.1080/24701394.2017.1307976
[23]
ZHOU T, SHEN X, IRWIN DM, et al. Mitogenomic analyses propose positive selection in mitochondrial genes for high-altitude adaptation in galliform birds[J]. Mitochondrion, 2014, 18: 70-75. DOI:10.1016/j.mito.2014.07.012
[24]
BAO H G, ZHAO C J, WU C X, et al. Associationof MT-ND5gene variation with mitochondrial respiratory control ratio and NADH dehydrogenase activityin Tibet chicken embryos[J]. Anim Genet, 2007, 38: 514-516. DOI:10.1111/j.1365-2052.2007.01622.x
[25]
YU L, WANG X P, ZHANG Y P, et al. Mitogenomic analysis of chinese snubnosed monkeys: evidence of positive selection in NADH dehydrogenasegenes in high-altitude adaptation[J]. Mitochondrion, 2011, 11(3): 497-503. DOI:10.1016/j.mito.2011.01.004
[26]
XU S Q, HUA S, HE J, et al. High altitude adaptation and phylogeny analysis of tibetan horse based on the mitochondrial genome[J]. J Genomics, 2007, 34(8): 720-729.
[27]
DEHAAN C, HABIBI-NAZHAD B, YAN E, et al. Mutation in mitochondrial complex Ⅰ ND6 subunit isassociated with defective response to hypoxia in human glioma cells[J]. Mol Cancer, 2004, 3: 19. DOI:10.1186/1476-4598-3-19
[28]
LUO Y J, GAO W X, GAO Y Q, et al. Mitochondrial genome analysis of ochotona curzoniae and implication of cytochrome c oxidase in hypoxic adaptation[J]. Mitochondrion, 2008, 8(5-6): 352-357. DOI:10.1016/j.mito.2008.07.005
[29]
CAPALDI, R A. Structure and function of cytochrome c oxidase[J]. Annu Rev Biochem, 1990, 59: 569-596. DOI:10.1146/annurev.bi.59.070190.003033
[30]
ZHANG ZY, CHEN B, ZHAO DJ, et al. Functional modulation of mitochondrial cytochrome c oxidase underlies adaptation to high-altitude hypoxia in a Tibetan migratory locust[J]. Proceedings: Biological Sciences, 2013, 280(1756): 20122758.
[31]
SCOTT G R, SCHULTE P M, EGGINTON S, et al. Molecular evolution of cytochrome c oxidase underlies high-altitude adaptation in the bar-headedgoose[J]. Mol Biol Evol, 2011, 28(1): 351-363. DOI:10.1093/molbev/msq205
[32]
ZHAO RZ, JIANG S, ZHANG L, et al. Mitochondrial electron transport chain, ROS generation and uncoupling (review)[J]. Int J Mol Med, 2019, 44(1): 3-15.
[33]
李耀东, 顾存林, 龙启福. 低氧对喜马拉雅旱獭线粒体ATPase6基因表达的影响[J]. 山东医药, 2012, 52(24): 30-31. DOI:10.3969/j.issn.1002-266X.2012.24.011
[34]
XIAOLING ZHAO. High-altitude adaptation of Tibetan chicken from MT-COI and ATP-6 perspective[J]. Mitochondrial DNA Part A, 2015, 27(5): 3280-3288.
[35]
李娜, 杨伟钊, 傅金钟. 基于线粒体基因组的沙蜥高海拔适应研究[J]. 四川动物, 2015, 34(6): 810-816.
[36]
GAO X, WEN X, ESSER L, et al. Structural basis for the quinone reduction in the bc1 complex: a comparative analysis of crystal structures of mitochondrial cytochrome bc1 with bound substrate and inhibitors at the qi site[J]. Biochemistry, 2003, 42(30): 9067-9080. DOI:10.1021/bi0341814