Reconstructing the evolutionary history of China: A caveat about inferences drawn from ancient DNA

The decipherment of the meager information provided by short fragments of ancient mitochondrial DNA (mtDNA) is notoriously difficult but is regarded as a most promising way toward reconstructing the past from the genetic perspective. By haplogroup-specific hypervariable segment (HVS) motif search and matching or near-matching with available modem data sets, most of the ancient mtDNAs can be tentatively assigned to haplogroups, which are often subcontinent specific. Further typing for mtDNA haplogroup-diagnostic coding region polymorphisms, however, is indispensable for establishing the geographic/genetic affinities of ancient samples with less ambiguity. In the present study, we sequenced a fragment (similar to 982 bp) of the mtDNA control region in 76 Han individuals from Taian, Shandong, China, and we combined these data with previously reported samples from Zibo and Qingdao, Shandong. The reanalysis of two previously published ancient mtDNA population data sets from Linzi (same province) then indicates that the ancient populations had features in common with the modem populations from south China rather than any specific affinity to the European mtDNA pool. Our results highlight that ancient mtDNA data obtained under different sampling schemes and subject to potential contamination can easily create the impression of drastic spatiotemporal changes in the genetic structure of a regional population during the past few thousand years if inappropriate methods of data analysis are employed.

鱇(鱼良)白鱼线粒体DNA控制区结构和种群遗传多样性分析

用特异性引物对鱇(鱼良)白鱼(Anabarilius grahami)DNA进行PCR扩增,获得了鱇(鱼良)白鱼线粒体DNA控制区基冈全序列(930bp).控制区T、C、A和G碱基组成为29.8%、22.5%、33.0和14.7%.对照其他已报道的鱼类控制区结构,对鱇(鱼良)白鱼控制区结构进行了分析,识别了其终止序列区、中央保守区和保守序列区,找到了终止相关的序列TAS以及保守序列(CSB-F、CSB-D、CSB-1、CSB-2、CSB-3).同时运用DNA分析软件对鱇(鱼良)白鱼一个驯养种群(中国科学院昆明动物研究所珍稀鱼类繁育中心)及两个自然地理种群(江川县明星鱼洞、江川县牛摩村)进行了遗传多样性分析.结果显示:两个自然种群存在较强基因交流,未出现遗传分化;人工驯养种群遗传多样性最高,种群复壮程度较好.

Identification, characteristic and phylogenetic analysis of type II DNA topoisomerase gene in Giardia lamblia

The genes encoding type II DNA topoisomerases were investigated in Giardia lamblia genome, and a type IIA gene, GlTop 2 was identified. It is a single copy gene with a 4476 by long ORF without intron. The deduced amino acid sequence shows strong homology to eukaryotic DNA Top 2. However, some distortions were found, such as six insertions in the ATPase domain and the central domain, a similar to 100 as longer central domain; a similar to 200 as shorter C-terminal domain containing rich charged residues. These features revealed by comparing with Top 2 of the host, human, might be helpful in exploiting drug selectivity for antigiardial therapy. Phylogenetic analysis of eukaryotic enzymes showed that kinetoplastids, plants, fungi, and animals were monophyletic groups, and the animal and fungi lineages shared a more recent common ancestor than either did with the plant lineage; microsporidia grouped with fungi. However, unlike many previous phylogenetic analyses, the "amitochondriate" G. lamblia was not the earliest branch but diverged after mitochondriate kinetoplastids in our trees. Both the finding of typical eukaryotic type IIA topoisomerase and the phylogenetic analysis suggest G. lamblia is not possibly as primitive as was regarded before and might diverge after the acquisition of mitochondria. This is consistent with the recent discovery of mitochondrial remnant organelles in G. lamblia.

High Altitude Adaptation and Phylogenetic Analysis of Tibetan Horse Based on the Mitochondrial Genome

研究测定了西藏那曲(4,500 m)、云南中甸(3,300 m)、云南德钦(3,300 m)地区3匹藏马线粒体全基因组序列.3个地区的藏马线粒体基因组全长以及结构均与韩国济州岛的马类似,但比瑞典马线粒体基因组短.藏马基因组在DNA序列上的两两相似性达99.3%.通过对线粒体蛋白编码区的分析发现,NADH6基因的蛋白序列在三匹藏马中均表现快速进化的现象.这表明NADH6基因在藏马高原适应进化过程中扮演着重要角色.此外,利用7匹藏马的D-loop区域序列以及与其亲缘关系较近的马的序列首次构建的藏马的系统发育树显示,那曲藏马与中甸、德钦藏马属于不同的分支,且存在较大的遗传多样性,表明藏马可能为多地区起源.

Phylogenetic relationships among the genera of the Penaeidae (Crustacea : Decapoda) revealed by mitochondrial 16S rRNA gene sequences

The phylogenetic relationships within the family Penaeidae are examined based on mitochondrial 16S rRNA gene sequence analysis of 30 species from 20 genera. The analysis generally supports the three- tribe scheme proposed by Burkenroad ( 1983) but it is not consistent with the five- group classification of Kubo ( 1949). Three clades are resolved: ( Penaeus sensu stricto + Fenneropenaeus + Litopenaeus + Farfantepenaeus + Marsupenaeus + Melicertus + Funchalia + Heteropenaeus), ( Metapenaeus + Parapenaeopsis + Xiphopenaeus + Rimapenaeus + Megokris + Trachysalambria) and ( Metapenaeopsis + Penaeopsis + Parapenaeus), corresponding to the Penaeini, Trachypenaeini and Parapenaeini respectively, while the affinities of Atypopenaeus and Trachypenaeopsis are obscure. The molecular data support that Miyadiella represents the juvenile stage of Atypopenaeus. Within the Trachypenaeini, Trachypenaeus sensu lato is clearly paraphyletic, while the monophyly of Penaeus sensu lato in the Penaeini is questionable.

The complete mitochondrial genome of the helmet catfish Cranoglanis bouderius (Silurifonnes : Cranoglanididae) and the phylogeny of otophysan fishes

The complete sequence of the 16,539 nucleotide mitochondrial genome from the single species of the catfish family Cranoglanididae, the helmet catfish Cranoglanis bouderius, was determined using the long and accurate polymerase chain reaction (LA PCR) method. The nucleotide sequences of C. bouderius mitochondrial DNA have been compared with those of three other catfish species in the same order. The contents of the C. bouderius mitochondrial genome are 13 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, and a non-coding control region, the gene order of which is identical to that observed in most other vertebrates. Phylogenetic analyses for 13 otophysan fishes were performed using Bayesian method based on the concatenated mtDNA protein-coding gene sequence and the individual protein-coding gene sequence data set. The competing otophysan topologies were then tested by using the approximately unbiased test, the Kishino-Hasegawa test, and the Shimodaira-Hasegawa test. The results show that the grouping ((((Characifonnes, Gymnotiformes), Siluriformes), Cyprinifionnes), outgroup) is the most likely but there is no significant difference between this one and the other alternative hypotheses. In addition, the phylogenetic placement of the family Cranoglanididae among siluriform families was also discussed. (c) 2006 Elsevier B.V. All rights reserved.

Mitochondrial d-loop sequence variation and phylogeny of gobiobotine fishes

The mitochondrial DNA control region is amplified and sequenced from 8 genera and 10 species of gobiobotine fishes. The phylogenetic tree of Gobiobotinae and some representative species of other Cyprinid subfamilies obtained by the method of neighborhood joining, maximum likelihood and maximum parsimony with Danio rerio as an outgroup indicates that Gobiobotinae fishes are a monophyletic group which is close to Gobioninae subfamily. Gobiobotinae should be included into subfamily Gobioninae in terms of phylogenetic analysis. The research result supports that Gobiobotinae can be divided into genus Xenophysogobio and Gobiobotia. Xenophysogabio is the most primitive genera in the subfamily.

Mitochondrial cytochrome b sequence variations and phylogeny of the East Asian bagrid catfishes

The mitochondrial DNA cytochrome b gene was sequenced from 8 bagrid catfishes in China. Aligned with cytochrome b sequences from 9 bagrid catfishes in Japan, Korea and Russia retrieved from GenBank, and selected Silurus meridionalis, Liobagrus anguillicauda, Liobagrus reini and Phenacogrammus interruptus as outgroups, we constructed a matrix of 21 DNA sequences. The Kimura's two-parameter distances were calculated and molecule phylogenetic trees were constructed by using the maximum parsimony (MP) and neighbor-joining (NJ) methods. The results show that (i) there exist 3-bp deletions of mitochondrial cytochrome b gene compared with cypriniforms and characiforms; (ii) the molecular phylogenetic tree suggests that bagrid catfishes form a monophyletic group, and the genus Mystus is the earliest divergent in the East Asian bagrid catfishes, as well as the genus Pseudobagrus is a monophyletic group but the genus Pelteobagrus and Leiocassis are complicated; and 60 the evolution rate of the East Asian bagrids mitochondrial cytochrome b gene is about 0.18%-0.30% sequence divergence per million years.

线粒体 DNA 相关疾病的系统发育学研究

线粒体是细胞内提供能量的细胞器，并负责调节细胞的程序化死亡。因遗传 缺陷引起的线粒体功能障碍会导致ATP 合成障碍、能量产生不足而出现一系列 病症。线粒体DNA 相关疾病目前日益受到广泛的关注，然而在线粒体DNA 相 关疾病领域常用的病例对照法容易受到遗传背景、群体分层、数据质量等方面因 素的影响以致经常得到假阳性结果。系统发育分析方法有助于解决这些方面的问 题，因此我们以此分析方法开展了如下工作： 首先我们对线粒体DNA C1494T 突变及其所属单倍型类群和氨基糖苷类药 物性耳聋之间的关联进行了研究。之前有研究报道了两个中国氨基糖苷类药物性 耳聋的家系，经过对先征者的线粒体DNA 全序列测定在这两个家系中都发现了 C1494T 突变。我们采用系统发育方法分析这两个家系先征者的全序列后，发现 这两个个体都属于线粒体单倍型A。巧合的是，在我们之前的研究中，在一个来 自武汉的汉族样本WH6980 中也有C1494T 突变，而且该个体同样也属于单倍型 类群A。这不由得使人想到：C1494T 突变可能是单倍型类群A 中一个分支的界 定位点，或者，受母系遗传背景的影响，该突变偏好于在A 类群中发生。那么 很有可能单倍型类群A对能够引发氨基糖苷类药物性耳聋的C1494T 突变的发生 有促进作用。 为了验证这个假设，我们从三个省份随机选取了553 个正常个体来检测 C1494T 突变，以调查该突变在普通人群中的发生频率。另外，我们从1823 个中 国人样本中筛选出属于单倍型类群A 的111 个体，在这111 个个体中检测C1494T 突变，以调查该突变是否为单倍型类群A 特有的变异位点。我们的结果表明： 在553 个随机样本群中没有检测到C1494T 突变，这说明该突变是一个稀有的突 变，在正常人群中发生的频率极低。另外，在111 个属于单倍型类群A 的样本 群中，我们也没有检测到C1494T 突变的存在，这说明该突变并非单倍型类群A 特有的变异位点。经过对带有C1494T 突变的全序列进行综合的系统发育学分析 表明，这三条序列的C1494T 突变应是来源于同一次突变事件。同时，根据序列 之间共享变异位点的状况推断，两个耳聋家系具有很近的母系亲缘关系，我们推 测这两个家系的C1494T 突变是来源于一位共同的近期母系祖先。综合这些结果 表明，C1494T 突变的发生与单倍型类群即母系遗传背景没有相关性，是属于在人群中频率极低的散发性突变。 我们又将系统发育分析的方法应用于肿瘤相关的线粒体DNA 突变的研究 中。线粒体在细胞的自由基产生以及细胞凋亡中扮演重要角色；有研究报道线粒 体功能的缺陷能导致癌症的发生，同时又有很多报道指出癌变组织的线粒体基因 组存在异常。为探讨癌组织中线粒体基因组的变异情况及其在癌症的发生和发展 中所扮演的角色，我们选取乳腺癌早期患者为研究对象。 近年来，有大量的研究都报导了在癌组织中存在高频率高密度的线粒体 DNA 体细胞突变，并认为这些突变在肿瘤发生过程中可能具有功能相关性。但 这些研究存在着不可回避的问题，总结来讲归为三个方面：①数据质量差：很多 前面报道的突变数据，经系统发育分析，发现存在不少因样本交叉污染所造成的 假阳性突变。②所测片断太短。大多数研究只是检测了D-loop 区，只占了整个 基因组的约十六分之一，很难全面反映问题。③对照设置有问题，许多研究进行 简单的case-control 分析，忽略了不同个体间母系遗传背景的差异，导致大量假 阳性突变的产生。 针对以上问题的存在，我们基于系统发育分析的方法在中国的乳腺癌病人中 开展了一项研究。分别取得10 例乳腺癌早期患者的癌组织、癌旁组织、以及远 端正常组织；提取了总DNA，对每一份组织样品进行线粒体全基因组测序。这 样做有两个研究目的：第一是调查在严格的质量控制手段下，排除交叉污染所造 成的假阳性突变后，在癌症病人中是否仍然能观察到高频率高密度的线粒体 DNA 体细胞突变。第二个目的是调查中国乳腺癌患者线粒体DNA 体细胞突变情 况，为乳腺癌早期诊断提供有效的信息。 为了避免过去部分研究中出现的问题，我们采取了相应的措施。首先是在严 格的质量控制下，我们对同一个病人的三种不同组织，即癌组织，癌旁组织，正 常组织中的线粒体基因组进行全序列测定的方式来检测体细胞突变。同时将所测 得线粒体全基因组序列进行系统发育的分析；结果表明属于同一个病人的不同组 织都能忠实的聚到一支，而且每一条全序列都完整的带有所属单倍型特有的界定 位点，不存在任何样本交叉或者污染的情况；这样进一步确证了我们数据的可靠 性。同时经过这样的比较，可以非常清晰的筛选出在某一种组织类型中所发生的 体细胞突变。来自10 例病人的癌组织，癌旁组织和正常组织的29 条（7#患者远端正常组 织没有得到）线粒体基因组全序列中，我们只检测到了两个体细胞突变（T2275C 和A8601G）。这两个突变经PCR 克隆试验的验证，在远端正常组织中均不存在， 是真实的体细胞突变。突变A8601G 同时出现在3#患者的癌组织和癌旁组织中； 提示线粒体DNA 的体细胞突变可能先于细胞的癌变，这可能对乳腺癌的早期临 床诊断具有一定的意义。突变T2275C 只出现在6#患者的癌组织中，该突变位于 16S rRNA 基因非常保守的位点上，我们对该突变位点潜在的生物学作用作了进 一步的讨论。 在我们研究的病例中体细胞突变的比率要远远低于先前的报道；造成差异的 主要原因可能是因为我们采取了严格的质量控制手段。综合我们的研究结果表 明，先前所报道的乳腺癌组织中的线粒体基因体细胞突变的频率存在被高估的现 象。目前线粒体突变与肿瘤发生的关系还有待于更深入的研究，同时我们呼吁本 领域的研究要加强数据质量的可靠性。我们的研究表明，系统发育的分析方法在 线粒体DNA 相关疾病的研究中是有效的，因此我们推荐在线粒体DNA 相关疾 病的研究中进行推广。

中国都江堰汉族人群线粒体DNA 控制区遗传变异与长寿的相关性研究

随着社会的进步和医疗卫生水平的不断提高，人类获得了更高的平均寿 命，很多国家都步入了老龄化社会的行列。由于长寿具有遗传的倾向，所以 科学家们致力于人类长寿及衰老性疾病发生机理的研究，目的是为了使人类 在获得更长寿命的同时能够抵御或减缓老年性疾病的侵袭，远离衰老带来的 困扰，享受高质量的生活。 线粒体是真核生物的重要细胞器，具有长度为约为16569bp 的环状 DNA 分子。在人类群体，特别是欧洲群体的相关性研究中，线粒体DNA （mtDNA）编码区和控制区的一些多态性位点显示出与长寿及衰老性疾病的 相关。特别是mtDNA 控制区的C150T 变异除了在多个长寿人群中富集之外， 更是具有改变mtDNA 重链复制起始位点的功能。 为了探讨mtDNA 控制区多态性位点与中国汉族长寿人群是否存在相关 性，本研究在中国四川省都江堰地区采集了556 名年龄90 至108 岁的互无 关系的长寿老人血液样本，其中男性202 名，女性354 名。同时还采集了 214 名长寿老人的亲属和312 名无关对照的血液样本，年龄分别在10 至69 岁之间和22 到73 岁之间。我们对这些样本的mtDNA 通过测序和RFLP 等 手段进行了扫描，采集并记录了mtDNA 单倍型类群信息和控制区位点多态 信息。 在该人群中，本研究发现mtDNA 的主要单倍型类群与长寿没有显著的 相关性，总体单倍型类群频率分布在三个组别中基本一致（p=0.318）。对 mtDNA 控制区C150T 变异的频率在三个组别中做了包括总体频率差异，区 分mtDNA 单倍型类群的频率差异，区分样本性别的频率差异以及mtDNA 单倍型类群与性别信息联合的频率差异的分析。虽然在个别的比较中得到了 显著差异，但经过多重检验校正后，结果均变得不显著。此外，对146、152189 和195 等四个同样处于mtDNA 控制区复制起始区域的变异位点的初步 分析，同样没有获得显著的差异。不支持此前在欧洲长寿人群和日本长寿人 群得到的结论。 综上所述，本研究第一次在中国汉族人群中对mtDNA 控制区多态性与 长寿的相关性进行了研究。mtDNA 与人类长寿的关系还有待于更深层次的 机理性研究和功能性研究来揭示。

重离子束辐照诱变啤酒酵母线粒体DNA研究

摘要 II 5. 在不同注量离子束辐照后筛选出来的呼吸缺陷型酵母菌株中扩增获得位 于第 12 染色体上的 SOF1 基因，而在同样的扩增体系中没有得到野生型 菌株的该基因。 6. 选取离子束辐照后筛选出来的呼吸缺陷型酵母菌株再次进行辐照，发现 其在低剂量范围（＜0.93Gy）辐照下非常敏感，而在高剂量范围（＞ 0.93Gy）又表现出一定程度的辐射抗性。 结论： 1. 离子束辐照酵母细胞，直接或间接作用于酵母线粒体DNA，导致线粒体 DNA损伤，形成呼吸缺陷的酵母菌株。 2. I 类内含子和 II 类内含子对于离子束辐照的敏感性不同： II 类内含子比较 稳定，II 类内含子可能利用自身编码的反转录酶通过目的DNA引导的反 转录机制对受到辐照损伤的II 类内含子进行修复。 3. 离子束辐照后 SOF1 基因可能发生了突变，影响酵母细胞的生长。 4. 呼吸缺陷型酵母菌株因其线粒体 DNA发生变化及线粒体功能的改变， 使 呼吸缺陷型酵母菌株在不同剂量区的离子束辐照下表现不同辐射敏感 性。目的： 研究啤酒酵母的线粒体 DNA 在重离子辐照作用下的突变效应及其突变机 理。 材料与方法： 利用兰州重离子研究装置(HIRFL)加速的氖、碳离子辐照酵母细胞，用 TTC 显色培养基筛选呼吸缺陷型酵母菌株，并用 mtDNA 限制性酶切手段分析其突变 规律。采用 PCR扩增并对目的产物测序的方法对辐照后线粒体DNA上的 I 类内 含子和 II类内含子进行研究。 结果： 1. TTC 显色实验表明：离子束辐照导致酵母线粒体上的电子传递链发生改 变，产生的还原氢减少，造成呼吸缺陷。 2. 利用限制性酶切实验对线粒体 DNA进行研究，结果表明：离子束辐照诱 变筛选出来的呼吸缺陷型酵母菌株其线粒体DNA变化明显： 主要表现为 酶切条带缺失严重。即使在同一注量下筛选出来的呼吸缺陷型酵母菌株， 其酶切图谱也不相同。 3. 通过 PCR 手段对辐照后酵母线粒体 DNA 碱基序列进一步进行分析，发 现经不同注量离子束辐照后筛选出来的呼吸缺陷型酵母菌株，其I 类内含 子（ai4 and ai5）经设计不同引物进行扩增，没有获得目的条带，说明此 序列发生了突变，可能对离子束辐照比较敏感。 4. 经不同注量离子束辐照后筛选出来的呼吸缺陷型酵母菌株，其 II 类内含 子（ai2）的碱基序列与野生型相比没有变化，表现出在离子束辐照作用 下比较稳定的特性。

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): Unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L crocea individuals and two L polyactis individuals. Interestingly, all seven CR sequences from L crocea were identical while the three sequences from L polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys. (c) 2008 Elsevier B.V. All rights reserved.

Tandem duplication-random loss is not a real feature of oyster mitochondrial genomes

Duplications and rearrangements of coding genes are major themes in the evolution of mitochondrial genomes, bearing important consequences in the function of mitochondria and the fitness of organisms. Yu et al. (BMC Genomics 2008, 9: 477) reported the complete mt genome sequence of the oyster Crassostrea hongkongensis (16,475 bp) and found that a DNA segment containing four tRNA genes (trnK(1), trnC, trnQ(1) and trnN), a duplicated (rrnS) and a split rRNA gene (rrnL5') was absent compared with that of two other Crassostrea species. It was suggested that the absence was a novel case of "tandem duplication-random loss" with evolutionary significance. We independently sequenced the complete mt genome of three C. hongkongensis individuals, all of which were 18,622 bp and contained the segment that was missing in Yu et al.'s sequence. Further, we designed primers, verified sequences and demonstrated that the sequence loss in Yu et al.'s study was an artifact caused by placing primers in a duplicated region. The duplication and split of ribosomal RNA genes are unique for Crassostrea oysters and not lost in C. hongkongensis. Our study highlights the need for caution when amplifying and sequencing through duplicated regions of the genome.