基于线粒体基因组的沙蜥高海拔适应研究

高海拔地区的物种容易受到低温、低氧、强紫外辐射等极端环境因素的影响。研究这些物种对特殊环境的应对反应,能够为进一步理解适应进化的机制提供重要线索。线粒体是细胞的能量代谢中心,因此线粒体基因组很可能在动物高原适应中起着重要作用。鬣蜥科Agamidae沙蜥属Phrynocephalus物种广泛分布于海拔1000~5300 m范围内,是研究高原适应的良好材料。本研究对2种高海拔沙蜥和6种低海拔沙蜥的线粒体基因组进行了比较研究,检测了可能经历过正选择的蛋白编码基因,探讨了线粒体基因在沙蜥高海拔适应中的作用。结果发现,在不同物种间,高海拔西藏沙蜥Phrynocephalus theobaldi的线粒体基因组中蛋白编码基因的进化速率最快;在不同基因间,ATP8具有最快的进化速率。使用分支-位点模型进行正选择检测,发现ATP8基因在西藏沙蜥中存在明显的正选择信号(P<0.05,ω>1)。通过贝叶斯方法进一步计算每个位点的后验概率,发现在ATP8基因上存在2个正选择位点。这些结果说明ATP8基因可能在西藏沙蜥高海拔适应中起到了重要的作用。但在同为高海拔的青海沙蜥Phrynocephalus vlangalii中,却没有发现类似的正选择信号,这揭示不同物种高海拔适应的分子机制可能不同。     

M、F型ATP8基因在雌雄三角帆蚌不同发育时期的表达差异研究

双壳贝类中有两种线粒体遗传方式,一种通过精子传递线粒体基因组,一种通过卵子传递线粒体基因组。这种现象被称为线粒体基因组的双单亲遗传(doubly uniparental inheritance, DUI)。本实验首先对三角帆蚌(Hyriopsis cumingii) F、M-type ATP8基因进行了克隆,得到了它们的cDNA全序列。F-type ATP8基因全长275 bp,编码49个氨基酸;M-type ATP8基因全长222 bp,编码58个氨基酸。在此基础上研究了三角帆蚌从胚胎时期到三龄成熟时期时F、M-type ATP8基因的表达差异。胚胎期到8月龄的三角帆蚌性腺部位组织团中,M、F-type ATP8基因在各时期均有表达,同时期对比发现,F-type ATP8基因的表达量明显高于M-type ATP8基因。12月龄雌雄三角帆蚌各组织中,F-type ATP8表达量高于M-type ATP8,处于主导地位。24月龄雌雄三角帆蚌各组织中,M、F-type ATP8均能检测到表达,M-type ATP8仅在性腺中高表达。36月龄雌雄三角帆蚌各组织中,F-type ATP8均能检测到,M-type ATP8表达量都非常低;而在雄性中,仅在性腺中高表达。研究表明,随着三角帆蚌的发育,M-type ATP8仅在雄性性腺中存在,在其他组织中选择性降解。     

抑制AX-3表达改变本氏烟氧化和光合磷酸化ATP积累

基因组注释预测了AX-3普遍存在于不同植物种类,但其功能鲜有报道。本研究通过StAX-3基因反义cDNA在本氏烟(Nicotiana benthamiana)体内表达,实现了其AX-3基因表达下调;蛋白质组表达谱测定分析显示,伴随AX-3积累减少,168个DAPs (表达倍绝对值大于1.5)的积累随之改变;其中15个密切关联DAPs网络分析显示,AX-3下调表达,主要降低了9种线粒体内代谢酶类的积累,包括氧化磷酸化ATP合成酶和TCA循环核心脱氢酶,同时,增加了光合作用代谢酶类合成,包括光合磷酸化ATP合成酶和叶绿素合成底物萜类酶等。叶绿素和ATP含量等表型测定结果,证实了萜类底物和叶绿素合成积累增加,促进光合磷酸化ATP合成积累。结果表明,StAX-3调控氧化和光合磷酸化ATP合成积累模式。该结果不仅明确了St AX-3具有调控线粒体和叶绿体能量代谢的功能,而且为进一步研究植物ataxin-3功能提供理论指导。     

小鼠2-细胞胚胎ATP合成酶6基因特异表达分析及鉴定

合子基因组活化是小鼠胚胎早期发育由细胞质调控向核调控转变的关键 .小鼠合子基因组活化发生在 2 细胞胚胎阶段 ,通过对 2 细胞胚胎阶段特异性表达基因的分析 ,可以从分子水平上揭示早期小鼠胚胎的发育机理 .用DD RTPCR技术 ,从单个小鼠 2 细胞胚胎与成熟卵母细胞 (MII细胞 )中分离了 2个差异片段 ,片段 2同小鼠睾丸中表达的一个未知片段具有高度同源性 .经过cDNA文库构建、筛选 ,分离到其全长cDNA .序列分析结果表明 ,该基因为小鼠ATP合成酶亚单位 6基因 .ATP合成酶亚单位 6基因由线粒体DNA编码 ,与细胞内ATP的合成相关 .小鼠 2 细胞胚胎特异表达的ATP合成酶亚单位 6基因可能与胚胎正常发育相关     

孟加拉笛鲷线粒体基因组序列结构及其进化

采用Long-PCR扩增线粒体全基因组方法得到了孟加拉笛鲷线粒体基因组全序列.序列分析结果表明,孟加拉笛鲷线粒体基因组序列全长16 511 bp,共有13个编码蛋白质基因、22个tRNA基因、2个rRNA基因和1个D-loop区.在编码蛋白质基因中,除COⅠ是以GTG作为起始密码子外,其它均是以ATG起始,NDⅠ、COⅡ、ND3以TAG作为终止密码子,而ND4、Cyt b则以不完全的T为终止密码子,其余8个蛋白质基因的终止密码子均为TAA.孟加拉笛鲷线粒体基因组各基因长度、位置与典型的脊椎动物相似,其编码蛋白质基因和rRNA基因与其它硬骨鱼类具有很高的同源性.基于14种笛鲷线粒体区段COⅠ、COⅡ和Cyt b基因的全序列合并成的一个组合数据集构建系统进化树,显示孟加拉笛鲷与四带笛鲷关系最为密切.     

太湖新银鱼线粒体D-loop和Cytb片段序列结构与进化速率比较

共获得49个太湖新银鱼(Neosalanx taihuensis)个体的线粒体细胞色素b(Cyt b)全序列和控制区(D-loop)部分序列。所测线粒体D-loop部分序列长度变化范围为648～680bp,识别到位于前端的一个串联重复序列、一个终止相关序列(ETAS),3个中央保守区保守序列(CSB-F、CSB-E、CSB-D)及一个保守序列区保守序列(CSB-1),结构与其他鱼类的研究结果类似。太湖新银鱼线粒体Cyt b和D-loop片段的相对进化速率的比较研究结果表明,太湖新银鱼D-loop总的序列多态性位点的比例为0.83%,低于线粒体Cyt b部分总的序列多态性位点的比例(1.31%)。假设太湖新银鱼Cyt b基因平均进化速率相对值为1,贝叶斯(Bayes)MCMC模拟给出Cyt b基因的相对速率区间估计为1.000±0.131,而D-loop基因的相对速率为0.859±0.261,表明太湖新银鱼D-loop基因的进化速率低于Cyt b基因,同时,后验概率分布的变异方差也比较大。说明Cyt b基因比D-loop基因具有相对较高的进化速率,也相对更接近分子钟假设。因此,可以认为Cyt b基因比D-loop基因更适于太湖新银鱼种内及近缘种间相关分子生态及系统地理格局的研究。     

麦穗鱼线粒体基因组序列测定及分析

利用麦穗鱼Pseudorasbora parva和相关鱼类的部分线粒体基因序列,设计出2对长批引物和30对短批引物,采用基于长PCR的2次PCR扩增法测定并注释麦穗鱼线粒体基因组全序列。结果表明,麦穗鱼线粒体基因组长16600bp,A+T含量为58.9%,37个基因位置及组成与其它硬骨鱼一致,均由13个蛋白编码基因、22个tRNA、2个rRNA基因和1个控制区(D-loop)组成。其中L链仅含8个tRNA(Pro、T yr、Ser、Ala、Asn、Cys、Glu、Gln)及ND6基因,其余基因皆由H链编码。基因排列紧密,间隔序列共计13处64bp,长度从1～32bp不等;基因重叠区7处23bp,重叠碱基数在1～7bp之间。13个蛋白编码基因中,除COI起始密码子为GTG外,其余均以ATG为起始密码子;有8个基因(ND1、ND2、COI、ATP6、ATP8、ND4L、ND5、ND6)3’端有完全的TAA或TAG终止密码子,其它5个基因终止密码子为不完整的TA(ND3和ND4)或T(COⅡ,COⅢ,Cyt b)。除tRNASer(AGY)外,其余21个tRNA基因的二级结构均为典型的三叶草结构。预测的lrRNA二级结构共有6个结构域,53个茎环结构,srRNA二级结构包含43个茎环结构。控制区(D-loop)存在3个结构区:终止序列区(TAS)、中央保守区(CSB-F、CSB-D)和保守序列区(CSB-1、CSB-2、CSB-3),其中TAS与DNA复制终止相关,出现茎环结构。     

SEQUENCE AND ANALYSIS OF COMPLETE MITOCHONDRIAL GENOME OF PSEUDORASBORA PARVA

利用麦穗鱼Pseudorasbora parva和相关鱼类的部分线粒体基因序列,设计出2对长批引物和30对短批引物,采用基于长PCR的2次PCR扩增法测定并注释麦穗鱼线粒体基因组全序列.结果表明,麦穗鱼线粒体基因组长16600 bp,A+T含量为58.9％,37个基因位置及组成与其它硬骨鱼一致,均由13个蛋白编码基因、22个tRNA、2个rRNA基因和1个控制区(D-loop)组成.其中L链仅含8个tRNA(Pro、Tyr、Ser、Ala、Asn、Cys、Glu、Gln)及ND6基因,其余基因皆由H链编码.基因排列紧密,间隔序列共计13处64 bp,长度从1～32 bp不等;基因重叠区7处23 bp,重叠碱基数在1～7bp之间.13个蛋白编码基因中,除COI起始密码子为GTG外,其余均以ATG为起始密码子;有8个基因(ND1、ND2、COI、ATP6、ATP8、ND4L、ND5、ND6)3端有完全的TAA或TAG终止密码子,其它5个基因终止密码子为不完整的TA (ND3和ND4)或T(COⅡ,COⅢ,Cyt b).除tRNAser(AGY)外,其余21个tRNA基因的二级结构均为典型的三叶草结构.预测的lrRNA二级结构共有6个结构域,53个茎环结构,srRNA二级结构包含43个茎环结构.控制区(D-loop)存在3个结构区:终止序列区(TAS)、中央保守区( CSB-F、CSB-D)和保守序列区(CSB-1、CSB-2、CSB-3),其中TAS与DNA复制终止相关,出现茎环结构.     

小鼠单个GV期卵母细胞中ATP8基因的表达分析

目的:探讨小鼠GV期卵母细胞线粒体中ATP8(ATP合酶亚基8)基因的表达情况。方法:应用挤压法从卵巢中分离获得生发泡期(germinal vesicle,GV)卵母细胞;用RT-PCR检测GV期单个卵母细胞中ATP8基因的表达:其中cDNA的合成分两种方法进行:一是将GV期单个卵母细胞直接进行RT合成cDNA,二是先用DNA酶加EcoRⅠ酶祛除mtDNA和核DNA后再进行RT;回收产物构建克隆质粒并测序。结果:1.5%琼脂糖电泳显示、测序结果均表明ATP8基因在GV期卵母细胞中有表达。结论:小鼠GV期卵母细胞特异表达的ATP8基因可能与卵母细胞的正常发育成熟相关。     

构建钙荧光探针细胞用于探究胞浆和线粒体钙对ATP刺激的响应

目的:构建表达基因编辑钙探针(GECIs)的细胞系HeLa-GECIs,探究细胞应答外界ATP刺激中钙离子在细胞内的响应和变化。方法:分别用能够直接通过荧光强度反映细胞胞浆内和线粒体内钙离子相对浓度的2种钙探针cyto-GCaMP6和4mt-GCaMP6感染HeLa细胞,获得2种表达钙离子探针的HeLa细胞系;在感染了2种腺病毒探针24 h后,用共聚焦荧光显微镜检测荧光探针在HeLa细胞内的表达情况;在表达2种钙探针的细胞的培养基中加入外源ATP,用Time-lapse成像动态观测技术观察HeLa细胞内钙离子对外环境中ATP的响应。结果:共聚焦荧光显微镜观察,确定95%以上的细胞表达了对应的钙离子指示荧光探针;Time-lapse成像动态观测技术观察发现,在细胞培养基中加入ATP后,细胞胞浆钙探针荧光强度瞬时(3～6 s)升至10倍,200 s后逐渐降低到基础水平;线粒体钙到达峰值(4倍)的时间稍滞后(5～8 s),并且回落更慢,300 s时至1.5倍。在ATP受体P2X7抑制剂A438079预处理的实验组,上述胞浆钙和线粒体钙浓度上升不明显。结论:构建了能在活体细胞内通过荧光探针实时监测钙离子响应胞外ATP刺激的细胞实验体系,为进一步深入探究ATP等危险信号导致细胞的炎性损伤机制奠定了基础。     

Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude

We performed a mitochondrial whole-genome comparison study in 40 Tibetan and 50 Han Chinese. All subjects could be classified into 13 haplogroups pertained to the Macrohaplogroup M and N that pitched different quadrants by principal component analysis. We observed a difference in the M9 haplogroup and identified 18 significant variants by comparing whole sequences between Tibetan and Han populations. Variants in ND2, COX2, tRNA alanine and 12S rRNA were predicted to confer increased protein stability in Tibetans. We compared the base substitutions of nonsynonymous (NS) versus synonymous (S) of 13 protein-encoding genes and found the NS/S values of the ATP6, ATP8, and Cyt b genes were larger (>1) in Tibetans than that in Han population. Our findings provide clues for the existence of adaptive selection for the ATP6, ATP8, Cyt b, ND2, COX2, tRNA alanine and 12S rRNA genes in Tibetans which likely contributed to adaptation to their specific geographic environment, such as high altitude.     

Complete nucleotide sequences and gene organization of mitochondrial genome of Bufo gargarizans

The complete mitochondrial genome of Bufo gargarizans was sequenced using overlapping polymerase chain reaction (PCR) amplicons (GenBank Accession No. DQ275350). The genome is 17,277 base pairs in length, containing 13 protein-coding genes (ATP6, ATP8, COI-III, ND1-6, ND4L, Cyt b), 2 ribosomal RNAs (12S rRNA and 16S rRNA), 22 transfer RNAs and a putative control region. We analyzed the sequence using bioinformatics methods comparing the obtained mtDNA sequence with other toads and frogs. Based on the concatenated nucleotide sequences of protein-coding genes, we constructed a phylogenetic tree with maximum likelihood (ML) and maximum parsimony (MP) methods and discussed the phylogenetic relationships among 11 species of Anura.     

Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Amplification of multiple copies of mitochondrial Cytochrome b gene fragments in the Australian freshwater crayfish,Cherax destructor Clark (Parastacidae: Decapoda)

Non-coding copies of fragments of the mitochondrial genome translocated to the nucleus or pseudogenes are being found with increasing frequency in a diversity of organisms. As part of a study to evaluate the utility of a range of mitochondrial gene regions for population genetic and systematic studies of the Australian freshwater crayfish, Cherax destructor (the yabby), we report the first detection of Cytochrome b (Cyt b) pseudogenes in crustaceans. We amplified and sequenced fragments of the mitochondrial Cyt b gene from 14 individuals of C. destructor using polymerase chain reaction (PCR) with primers designed from conserved regions of Penaeus monodon and Drosophila melanogaster mitochondrial genomes. The phylogenetic tree produced from the amplified fragments using these primers showed a very different topology to the trees obtained from sequences from three other mitochondrial genes, suggesting one or more nuclear pseudogenes have been amplified. Supporting this conclusion, two highly divergent sequences were isolated from each of two single individuals, and a 2 base pair (bp) deletion in one sequence was observed. There was no evidence to support inadvertent amplification of parasite DNA or contamination of samples from other sources. These results add to other recent observations of pseudogenes suggesting the frequent transfer of mitochondrial DNA (mtDNA) genes to the nucleus and reinforces the necessity of great care in interpreting PCR-generated Cyt b sequences used in population or evolutionary studies in freshwater crayfish and crustaceans more generally.     

Differential gene expression during apoptosis induced by a serum factor: Role of mitochondrial F<Subscript>0</Subscript>-F<Subscript>1</Subscript> ATP synthase complex

The number of genes that are up regulated or down regulated during apoptosis is large and still increasing. In an attempt to characterize differential gene expression during serum factor induced apoptosis in AK-5 cells (a rat histiocytoma), we found subunit 6 and subunit 8 of the transmembrane proton channel and subunit alpha of the catalytic core of the mitochondrial F₀-F₁ ATP synthase complex to be up regulated during apoptosis. The increase in the expression levels of these subunits was concomitant with a transient increase in the intracellular ATP levels, suggesting that the increase in cellular ATP content is a result of the increase in the expression of ATP synthase subunits' gene and de novo protein synthesis. Depleting the cellular ATP levels with oligomycin inhibited apoptosis significantly, pointing to the requirement of ATP during apoptosis. Caspase 1 and caspase 3 activity and the loss of mitochondrial membrane potential were also inhibited by oligomycin during apoptosis in these cells, suggesting that the oligomycin induced inhibition of apoptosis could be due to inhibition of caspase activity and inhibition of mitochondrial depolarization. However, cytochrome C release during apoptosis was found to be completely independent of intracellular ATP content. Besides the ATP synthase complex genes, other mitochondrial genes like cytochrome C oxidase subunit II and III also showed elevated levels of expression during apoptosis. This kind of a mitochondrial gene expression profile suggests that in AK-5 cells, these genes are upregulated in a time-linked manner to ensure sufficient intracellular ATP levels and an efficient functioning of the mitochondrial respiratory chain for successful completion of the apoptotic pathway.     

Identification of camelid specific residues in mitochondrial ATP synthase subunits

ATP synthase is an enzyme involved in oxidative phosphorylation from prokaryotic to eukaryotic cells. In mammals it comprises at least 16 subunits from which the mitochondrial encoded ATP6 and ATP8 are essential. Mitochondrial genes variations have been suggested to allow rapid human and animal adaptation to new climates and dietary conditions (Mishmar et al. 2003). Camelidae taxa are uniquely adapted to extremely hot and dry climates of African-Asian territories and to cold and hypoxic environments of the South American Andean region. We sequenced and analyzed ATP6 and ATP8 genes in all camelid species. Based on the available structural data and evolutionary conservation of the deduced proteins we identified features proper of the group. In Old World camels the ATP8, important in the assembly of the F0 complex, showed a number of positively charged residues higher than in the other aligned species. In ATP6 we found the camelid specific substitutions Q47H and I106V that occur in sites highly conserved in other species. We speculate that these changes may have functional importance.     

Complete mitochondrial genome of the bullhead torrent catfish, Liobagrus obesus (Siluriformes, Amblycipididae): Genome description and phylogenetic considerations inferred from the Cyt b and 16S rRNA genes

Mitochondrial DNA (mtDNA) from the bullhead torrent catfish, Liobagrus obesus, was isolated by long-polymerase chain reaction (Long-PCR) with universal primers and was fully sequenced by primer working using flanking sequences. The complete mtDNA from L. obesus was 16,531 bp in length and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region, demonstrating a structure very similar to that of other bony fish. An analysis of the protein-coding genes revealed a statistically substantiated bias in (T+C): (A+G) content, supporting earlier findings regarding this peculiarity. As indicated by a chi-square test, the observed scores for pyrimidine and purine content were different from those expected assuming a 50:50 ratio: chi(2)=41.63, d.f.=5, p<0.000001 for three categories, including the 1st, 2nd, and 3rd codon positions. Further, there was a difference in nucleotide content between ND6 and the other 12 protein-coding genes in L. obesus. The values of p-distances, as summarized for different scales of evolutionary history at the Cyt b gene, revealed a clear pattern of increased nucleotide diversity at four levels: (1) intraspecies, (2) intragenus, (3) intrafamily, and (4) intraorder. Scores of average p-distances of the four categories in catfish were (1) 1.59+/-0.54%, (2) 5.28+/-1.72% (3) 16.37+/-1.26%, and (4) 19.81+/-0.14%, respectively. These data support the hypothesis that speciation in the order Siluriformes, in most cases, follows a geographic mode through the accumulation of a numerous small genetic changes over a long time period. A phylogenetic tree for the bullhead torrent catfish and several other fish species belonging to the order Siluriformes was developed on the basis of respective Cyt b sequences (1138 bp); the analysis revealed a monophyletic origin for the five examined families. A species-specific clustering of sequences from single species was obtained, supporting additionally basic phylogenetic information for the catfish and the barcoding suitability of Cyt b sequence data. Lastly, one of the well-supported properties of our phylogenetic tree (99% repetition level in our analysis) was the monophyletic placement of all catfish (order Siluriformes) among other ray-finned fish of the class Actinopterigii. Also discussed herein are the aspects of phylogeny based on the 16S rRNA gene.     

Complete nucleotide sequence and gene arrangement of the mitochondrial genome of the crab-eating frog Fejervarya cancrivora and evolutionary implications

The complete nucleotide sequence of the mitochondrial (mt) genome of the crab-eating frog, Fejervarya cancrivora Gravenhorst (Amphibia: Anura: Ranidae), was determined. The mt genome is 17,843 bp long and contains 13 protein-coding (ATP6, ATP8, COI-III, ND1-6 and 4L, and Cyt b) and two ribosomal RNA (12S and 16SrRNA) genes. Although metazoan mt genomes typically encode 22 transfer RNA genes (tRNAs), the F. cancrivora mtDNA contains 23 tRNAs due to the presence of an extra copy of tRNA(Met). A major noncoding region and a prominent intergenic spacer corresponding to the control region and light-strand replication origin were also found. To confirm the phylogenetic position of F. cancrivora, we compared the gene arrangement with that of other anurans and performed phylogenetic analyses based on mt genomic data. The genome organization of F. cancrivora mtDNA differs from that of typical vertebrates and neobatrachian frogs but is identical with that of F. limnocharis, suggesting that the unique gene arrangement occurred in the common ancestor of the genus. Phylogenetic analyses supported the monophyly of the Fejervarya species used here as well as the dicroglossini clade. Although the family Ranidae as previously recognized (= Ranidae, Discoglossidae, and some other natatanuran families; sensu Frost et al., 2006) is shown as a clade in the maximum parsimony analysis, the maximum likelihood and the Bayesian analyses suggest the paraphyly of the Ranidae with respect to the families, Mantellidae and Rhacophoridae. Three-tandem duplications of gene regions followed by subsequent deletions of supernumerary genes were proposed to explain the evolution of the extra tRNA(Met) and translocation of ND5 from the original neobatrachian gene order.