Sun skink diversification across the Indian-Southeast Asian biogeographical interface

Aim Widespread, transcontinental vertebrate groups represent ideal systems for biogeographical studies, because they can shed light on a wide range of questions relating to species diversification across the geographical template. We combined extensive geographical and genetic sampling from across multiple biogeographical realms to examine the timing and location of diversification in Asian sun skinks, a clade characterized by problematic species boundaries and a particularly enigmatic evolutionary history. Location Indian subcontinent, the Philippines, Southeast Asia and Sundaland. Methods We sequenced one mitochondrial and nine nuclear genes for most species in the genus Eutropis, and estimated phylogenetic relationships and divergence times using coalescent methods. To investigate the location of diversification events, we also estimated ancestral geographical ranges using several methods. Finally, we explored patterns of genetic diversity within several poorly understood, but widely distributed species. Results Divergence-time estimates indicate that Eutropis began to diversify during the Eocene. Biogeographical reconstructions show that species diversification was associated with dispersal into three biogeographical realms: India, Sundaland and the Philippines. Main conclusions The results of this study clarify several questions related to the evolutionary history of Eutropis, and place them in the context of classic Southeast Asian biogeography. Our study represents one of the first to compile a heavily sampled multilocus dataset ranging across international boundaries in southern Asia that have historically prevented a unified understanding of biogeographical and evolutionary processes involving the Indian subcontinent, mainland southern Asia and the island archipelagos of Southeast Asia.

分子系统学分析平台的建立和应用

分子系统学建立在实验和计算的基础之上。DNA快速测序技术的普及为分子系统学家提供了大量数据，而序列分析技术则是探索数据发现知识的重要工具。在基因组时代，随着大量模式生物完整基因组序列的获得，分子系统学正面临着前所未有的机遇和挑战。一方面，生命之树计划有助于确定新的模式生物和开展相应的基因组计划;另一方面，模式生物的基因组计划有助于阐明它们之间的进化关系和基因组的进化模式。更为重要的是，分子系统学序列分析技术已经发展成为探索与整合基因组数据的强有力工具，从而在生命科学中发挥重要作用。事实上，分子系统学和基因组学的相互渗透正在形成一门崭新的交叉学科——系统发育基因组学。 为了奠定分子系统学研究中信息管理和数据分析工作的坚实基础，我们建立了分子系统发育分析平台。该平台为研究人员提供专业数据库服务和数据分析技术支持，以及相关的网络资源。 分子系统发育分析平台包括了3个专业数据库。第一个是DNA凭证标本数据库。该数据库中的记录包括了7项字段：英文科名、中文科名、物种拉丁名、采集人、采集号、采集地和采集时间。用户可以通过设定单个或多个字段的取值进行检索。截止2004年6月1日，该数据库共包括3491条标本记录。第二个是引物数据库。PCR引物是分子系统学实验的重要条件之一。该数据库中的记录包括3项字段：引物名称、序列内容和退火温度。用户可以通过设定单个或多个字段的取值进行检索。截止2004年6月1日，该数据库共包括170条用于扩增植物细胞核、叶绿体和线粒体基因组DNA序列的引物记录。第三个是生物计算数据库。该数据库为研究人员提供传输和保存序列分析数据和结果文件的服务。 为了确保数据库的安全性和使用性，我们开发了数据库的接口和检索工具，以及系统管理员和用户资格认证程序。通过前者，使用者可以进行数据的上传、下载、管理和检索等操作。而后者则是对不同使用者身份和权限进行设定。管理员的权限高于用户，主要负责本系统的日常维护和管理工作，以及对新增管理员和用户进行资格认证。 分析技术支持旨在帮助用户快速掌握常用的系统发育分析方法，进行有效的数据分析，从复杂的统计学算法和计算机程序中解放出来，将精力集中于计算结果的生物学解释。在该部分中，我们首先简要介绍了常用的分析方法，并且针对分子系统学中的不同问题提供了相应的解决方案。这些问题包括：系统发育重建、替代速率和分歧时间的估计、祖先分布区的重建、性状进化假说的检验、以及密码子水平适应性进化的检测。我们特别强调了似然比检验和贝叶斯推测作为方法论上的重要进展在分子系统学中所发挥的关键作用。本部分还包括大量常用的分子系统学程序或软件包及其快速使用说明和命令模块。下载安装之后，用户即可按照说明使用命令模块进行数据分析。 此外，该平台还提供了一些常用的网络资源地址，如生物信息中心、分子进化和系统发育实验室、专业期刊和相关数据库等。 最后还给出了4个应用实例，即针对特定分子系统学问题的解决方案和初步的分析结果。 第一个例子说明系统发育重建方法的应用。为了确定杨梅科的系统学位置，对6种DNA序列和叶绿体trnL-F区内的间隔性状进行了分析。单个分析表明这6种序列之间在系统学信息上存在显著差异。叶绿体基因组序列的合并分析强烈支持杨梅科和（木麻黄科，（桦木科，核果桦科））的姐妹群关系，而间隔性状的存在能够充分提高其分辨率和支持率。 第二个例子说明如何推测历史生物地理学过程。我们对壳斗目8科25属植物叶绿体基因组的trnL-F、matK、rbcL和atpB的合并序列进行了最大简约分析，得到唯一的最大简约树。基于该系统树和25属植物的地理分布数据，采用扩散－替代分析方法重建了系统树每个节点上的祖先分布区，推测了壳斗目的分布历史。结果表明，壳斗目的历史生物地理学过程由3次替代事件和20次扩散事件组成。其中最重要的替代事件是由于冈瓦纳大陆和劳亚大陆分离所导致的南青冈科及其姐妹群之间的分化。另外，在壳斗科和核心高等金缕梅类中多次发生从欧亚大陆到北美洲、甚至南美洲的平行扩散事件。 第三个例子说明如何估计分歧时间。我们仍然使用扩散－替代分析中所用的最大简约树作为分析的依据，并根据等级制似然比检验确定的最优替代模型对该系统树的支长进行了最大似然优化。似然比检验表明，该系统树不服从分子钟假说。我们以冈瓦纳大陆和劳亚大陆分离的地质事件和5个属的最早化石记录作为标定点，采用罚分似然法在没有分子钟的条件下估计了壳斗目的科间分歧时间。结果表明，绝大多数科间分歧事件都发生在白垩纪。 第四个例子说明如何检测密码子水平的适应性进化。分支间可变选择压力模型的似然比检验表明SARS冠状病毒的S基因在跨种传播过程中发生了正选择。

稻族的分子系统学研究

稻族（Oryzeae）是禾本科中包含多种经济植物的重要类群，现有大约12个属，广布于全球的热带和温带地区。虽然有证据表明稻族是一个单系类群，但稻族的分类处理和属间系统发育关系以及稻族的生物地理学等方面仍存在许多悬而未决的问题。本研究利用了7个DNA片段，包括2个叶绿体基因片段（trnL和matK）、1个线粒体基因片段（nad1）和4个分布在不同染色体上的核基因片段（Adh1、Adh2、GPA1和Waxy）以及形态性状（87个）重建了稻族的系统发育关系，并在此基础上利用6个DNA片段（trnL、matK、nad1、Adh1、Adh2和GPA1）探讨了稻属基部类群的系统发育关系。同时在系统发育研究的基础上，进一步探讨了稻族物种间的分歧时间及生物地理学问题。主要结果如下： 1.稻族的系统学 多基因和形态性状的分析表明：1）稻族是单系类群可分为两个主要分支，相当于传统的两个亚族，第一个亚族（Zizaniinae）包括稻属﹑假稻属和Porteresia，而第二个亚族（Oryzinae）包括其余8个属。两性花是稻族的原始状态，而单性花多次起源，共起源了3次，因此单性花和两性花的区别不宜作为划分亚族的依据;2）Zizania 与分布在南美的单型属Rhynchoryza关系最近;3）一些单型属（Hydrochloa、Porteresia和Prosphytochloa）的确立得不到分子证据的支持，特别是Hydrochloa，形态证据分析的结果明确不支持建立该单型属;4）形态分析结果支持Zizania 4个种聚为一支，而东亚的一个种（Zizania latifolia）为该属基部类群。 2.稻属基部类群的系统发育关系 多基因系统发育分析表明：1）稻属是一个单系类群;2）所有数据都支持稻属10种染色体组类型各自为单系，包括稻属中最后确定的染色体组类型HK（Oryza schlechteri和Porteresia coarctata）;3）F染色体组与HK和HJ中的H染色体组关系最近，类似于E染色体组与异源四倍体的CD染色体组中的D染色体组关系最近; 4）双亲遗传的核基因和母系遗传的线粒体和叶绿体基因的对比分析表明，HJ染色体组两个物种的母系来源是H，而HK染色体组两个物种母系来源是K。三个染色体组相比较，H分化较早，J次之，K分化最晚;5）在系统发育研究基础上，本研究认为稻属4个复合体的划分是可信的，原来未划分到任何一个复合体中的两个种O. brachyantha和O. schlechteri以及新近归入稻属的Oryza coarctata都应归入O. ridleyi复合体。 3.稻族各谱系的分化时间及其生物地理学意义 利用分子钟及其改进方法对稻族各谱系的起源和分化时间进行了估测，并在此基础上探讨了形成各类群地理分布式样的可能原因。结果表明，稻族在始新世（Eocene）晚期（3640万年前）起源于东南亚随后分为两大支。稻属基部类群G染色体组物种与稻属其他物种在1200万年前分开;稻属中A/B/C/E染色体组类群在700万年前开始分化。 在稻族的第一大支Oryzinae亚族中，稻属和假稻属1400万年前分开后，通过远程扩散及随后的漫长的历史过程形成了目前的分布格局。稻属从东南亚起源并逐渐扩散到非洲、澳洲及美洲;而假稻属则从非洲出发扩散到全球的热带和亚热带地区。这两个属的演化历史非常相似，都是全球广泛分布且起源时间较晚。 在稻族第二大支Zizaniinae亚族中，Zizania是稻族中唯一欧亚-北美间断分布属，与其分布在南美的姊妹群Rhynchoryza在2554万年前分开，这一支可能是先从南美扩散到东亚，然后再从东亚扩散到北美;Zizaniopsis和Luziola两属在南美洲和北美洲都有分布，其分歧时间为2180万年前，这两属具有类似的进化历史，即通过上新世（Pliocene）末期隆起的巴拿马陆桥形成现在的分布格局;Chikusichola与Potamophila/Prosphytochloa这一支从中新世（Miocene）早期（2270万年前）开始分化，由于Chikusichola（分布在东南亚）、Potamophila（澳大利亚）和Prosphytochloa（非洲）相距很远，这3个属的扩散必然与跨洋远距离扩散有关。

稻属A基因组的分子系统发育--兼论亚洲栽培稻的起源

稻属Oryza隶属禾本科Poaceae，包括20多个野生种和2个栽培种(亚洲栽培稻O. sativa L和非洲栽培稻O. glaberrima Steud) ，广泛分布于全球热带和亚热带。稻属物种可划分为10个基因组(又称染色体组)类型：A, B, C, BC, CD, E, F, G, HJ 和 HK。栽培稻所属的A基因组是稻属中物种数目最多、地理分布最广的基因组类型，由8个种组成。由于栽培稻属于A基因组，故A基因组物种是栽培稻遗传改良的巨大基因源。数十年来，国际上许多学者对A基因组类群开展了大量涉及形态、细胞、同工酶和分子标记方面的研究，但由于A基因组物种间遗传关系十分接近，形态上差异小且地理分布重叠，使得A基因组物种的系统发育、物种起源和生物地理学等方面存在诸多悬而未决的问题，是稻属中分类和鉴定困难较多的类群。本文利用核基因内含子序列，结合转座子插入分析，重建了A基因组的系统发育，估测了各类群的分化时间;与此同时，基于多克隆测序和基因谱系分析，探讨了O. rufipogon和O. nivara遗传关系以及亚洲栽培稻起源。主要研究结果如下： 1. A基因组的系统发育 在水稻全基因组数据库搜索的基础上，测定了4个单拷贝核基因（Adh1 及3个未注释基因）的内含子序列，构建了稻属A基因组8个种的系统发育关系。基于最大简约法和贝叶斯法的系统发育分析表明：1）澳大利亚的O. meridionalis为A基因组的基部类群;2）亚洲栽培稻两个亚种O. sativa ssp. japonica 和 O. sativa ssp. indica分别和不同的野生类群聚为独立的两个分支，支持japonica 和 indica为多次起源;3）O. rufipogon和O. nivara在系统发育树上完全混在一起，显示出二者间不存在遗传分化;4）非洲一年生野生种O. barthii是非洲栽培稻O. glaberrima的祖先，而非洲多年生野生种O. longistaminata与O. glaberrima/O. barthii.亲缘关系较远;5）分子钟方法估测A基因组类群约在2百万年前（2.0MYA）开始分化，亚洲栽培稻和非洲栽培稻，以及亚洲栽培稻的两个亚种则分别在0.7和 0.4 MYA左右开始分化。此外，通过核基因内含子序列与其它常用片段如ITS，matK等对比分析表明，进化速率相对较快的核基因内含子序列可以有效地用于近缘类群的系统发育研究。 2. Oryza rufipogon 和O. nivara群体遗传研究及亚洲栽培稻起源 对于亚洲野生类群O. rufipogon和O. nivara是合并为一个种还是处理为两个独立的种一直存在争议。在系统发育研究基础上，我们选取4个核基因内含子或5’-UTR区（Waxy, LHS，CatA和1个未注释基因），对采自整个分布区的群体样品进行了多克隆测序，结果表明：1）检测到O. rufipogon和O. nivara均有较高的核苷酸多态性，4个位点上π值和θw值平均分别为0.011和0.014;2）且二者在遗传上没有明显分化，两个类群在4个核基因位点上均检测到大量共享多态（shared polymorphism），未发现固有差异(fixed difference)，表明它们历史上可能属于一个大群体，支持将二者作为种内不同生态型或亚种处理;3）基因谱系树表明亚洲栽培稻的两个亚种indica和japonica分别和不同的O. rufipogon (包括O. nivara)群体聚在一起，进一步从基因谱系角度支持亚洲栽培稻多次起源假说。 3．转座子在群体遗传与系统发育研究中的应用 鉴于目前植物谱系地理学研究中缺乏具有足够信息量的分子标记用于检测种内遗传变异，我们选取3个核基因中的转座子，通过对取自O. rufipogon和O. nivara整个分布区的37份样品的克隆测序，探讨了进化速率快、信息含量丰富的转座子序列在群体遗传上的应用。结果表明：1）无论在物种水平还是群体水平，转座子能检测到比包括内含子在内的其它DNA区域高得多的遗传变异;2）在物种水平上，异交多年生的O. rufipogon和自交一年生的O. nivara多样性均较高，且2个种间相差很小，二者在3个位点上平均核苷酸多样性π值均为0.013，差别主要表现在O. rufipogon杂合位点比例（46.1%）明显高于O. nivara（9.1%），说明交配系统不同并不一定和物种多样性水平相关;3）是否发生转座子序列插入是有价值的系统发育信息，发生在不同染色体上3个基因中的转座子插入进一步证实A基因组基部类群是O. meridionalis;通过叶绿体中3个转座子的插入现象推断了稻族一些四倍体物种，如稻属BC基因组的一些类群的母本来源。

DIFFERENT RATES OF MITOCHONDRIAL-DNA SEQUENCE EVOLUTION IN KIRK DIK-DIK (MADOQUA-KIRKII) POPULATIONS

We have investigated evolutionary rates of the mitochondrial genome among individuals of Madoqua kirkii using the relative rate test. Our results demonstrate that individuals of two chromosome races, East African cytotype A and Southwest African cytotype D, evolve about 2.3 times faster than East African cytotype B. Cytogenetic changes, DNA repair efficiency, mutagens, and more likely, hitherto unrecognized factors will account for the rate difference we have observed. Our results suggest additional caution when using molecular clocks in the estimation of divergence time, even within lineages of closely related taxa. Rate heterogeneity in microevolutionary timescales represents a potentially important aspect of basic evolutionary processes and may provide additional insights into factors which affect genome evolution. (C) 1995 Academic Press, Inc.

Phylogenetic relationships of five pika species (genus Ochotona) based on mitochondrial DNA restriction maps

Restriction site mapping of mitochondrial DNA (mtDNA) with 16 restriction endonucleases was used to examine the phylogenetic relationships of Ochotona cansus, O. huangensis, O. thibetana, O. curzoniae and O. erythrotis. A 1-kb length variation between O. erythrotis of subgenus Pika and other four species of subgenus Ochotona was observed, which may be a useful genetic marker for identifying the two subgenera. The phylogenetic tree constructed using PAUP based on 61 phylogenetically informative sites suggests that O. erythrotis diverged first, followed by O. cansus, while O. curzoniae and O. huangensis are sister taxa related to O. thibetana, The results indicate that both O. cansus and O. huangensis should be treated as independent species. If the base substitution rate of pikas mtDNA was 2% per million years, then the divergence time of the two subgenera, Pika and Ochotana, is about 8.8 Ma ago of late Miocence, middle Bao-dian of Chinese mammalian age, and the divergence of the four species in subgenus Ochotona would have occurred about 2.5 - 4.2 Ma ago, Yushean of Chinese mammalian age. This calculation appears to be substantiated by the fossil record.

Molecular phylogenetics and biogeography of Lepus in Eastern Asia based on mitochondrial DNA sequences

In spite of several classification attempts among taxa of the genus Lepus, phylogenetic relationships still remain poorly understood. Here, we present molecular genetic evidence that may resolve some of the current incongruities in the phylogeny of the leporids. The complete mitochondrial cytb, 12S genes, and parts of ND4 and control region fragments were sequenced to examine phylogenetic relationships among Chinese hare taxa and other leporids throughout the World using maximum parsimony, maximum likelihood, and Bayesian phylogenetic reconstruction approaches. Using reconstructed phylogenies, we observed that the Chinese hare is not a single monophyletic group as originally thought. Instead, the data infers that the genus Lepus is monophyletic with three unique species groups: North American, Eurasian, and African. Ancestral area analysis indicated that ancestral Lepus arose in North America and then dispersed into Eurasia via the Bering Land Bridge eventually extending to Africa. Brooks Parsimony analysis showed that dispersal events followed by subsequent speciation have occurred in other geographic areas as well and resulted in the rapid radiation and speciation of Lepus. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches estimated the divergence time between the three major groups within Lepus. The genus appears to have arisen approximately 10.76 MYA (+/- 0.86 MYA), with most speciation events occurring during the Pliocene epoch (5.65 +/- 1.15 MYA similar to 1.12 +/- 10.47 MYA). (c) 2005 Elsevier Inc. All rights reserved.

Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China

Leptobrachium ailaonicum is a vulnerable anuran restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes at mid-elevations (approximately 2000-2600 m) in the subtropical Mount Wuliang and Mount Ailao ranges in southwest China (Yunnan Province) and northern Vietnam. Given high habitat specificity and lack of suitable habitat in lower elevations between these ranges, we hypothesized limited gene flow between populations throughout its range. We used two mitochondrial genes to construct a phylogeographic pattern within this species in order to test our hypothesis. We also examined whether this phylogeographic pattern is a response to past geological events and/or climatic oscillations. A total of 1989 base pairs were obtained from 81 individuals of nine populations yielding 51 unique haplotypes. Both Bayesian and maximum parsimony phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographical regions separated by deep river valleys. These results suggest a long history of allopatric separation by vicariance. The distinct geographic distributions of four major clades and the estimated divergence time suggest spatial and temporal separations that coincide with climatic and paleogeographic changes following the orogeny and uplift of Mount Ailao during the late Miocene to mid Pliocene in southwest China. At the southern distribution, the presence of two sympatric yet differentiated clades in two areas are interpreted as a result of secondary contact between previously allopatric populations during cooler Pleistocene glacial cycles. Analysis of molecular variance indicates that most of the observed genetic variation occurs among the four regions implying long-term interruption of maternal gene flow, suggesting that L ailaonicum may represent more than one distinct species and should at least be separated into four management units corresponding to these four geographic lineages for conservation. (C) 2009 Elsevier Inc. All rights reserved.

Variation in mitochondrial DNA and phylogeny of six species of pikas (Ochotona)

Phylogenetic relationships of six species of Ochotona were investigated using mitochondrial DNA (mtDNA) restriction-site analysis. The phylogenetic tree constructed using PAUP was based on 62 phylogenetically informative sites with O. erythrotis designated as an outgroup. Two clades were evident. One contained O. curzoniae, O. huangensis, and O. thibetana. in the second, O. daurica was a sister taxon of O. cansus. The five species appear to come from different maternal lineages. The branching structure of the tree and sequence divergence confirm that huangensis and cansus are distinct species, and that mol-osa is a synonym of O. cansus rather than O. thibetana. Divergence time, estimated from genetic distances, indicates that the ancestors of the two maternal lineages diverged ca. 6.5 x 10(6) years ago. O. curzoniae diverged from O. huangensis, and O. daurica diverged from O. cansus, at about the same time (ca. 3.4 x 10(6) years ago). These data suggest a period of rapid radiation of the genus Ochotona in China, perhaps during the late Pliocene. These calculations correspond roughly to tectonic events and environmental changes in China throughout this period, and appear to be substantiated by the fossil record.

Molecular phylogeny and biogeography of woolly flying squirrel (Rodentia : Sciuridae), inferred from mitochondrial cytochrome b gene sequences

To investigate the genetic diversity between the populations of woolly flying squirrels (Eupetaurus) from the eastern and western extremes of the Himalayas, partial mitochondrial cytochrome b gene sequences (390-810bp) that were determined from the museum specimens were analyzed using maximum parsimony (MP) and maximum likelihood (ML) methods. The molecular data reveal that the two specimens that were collected in northwestern Yunnan (China) are members of the genus Eupetaurus. Reconstructed phylogenetic relationships show that the populations of Eupetaurus in the eastern and western extremes of the Himalayas are two distinct species with significant genetic differences (12%) and diverged about 10.8 million years ago. Eupetaurus is significantly different from Petaurista and Pteromys. The level of estimated pairwise-sequence divergence observed between Eupetaurus and Petaurista or Pteromys is greater than that observed between Eupetaurus and Trogopterus, Belomys, Glaucomys, or Hylopetes. Considering the divergence time of the two Eupetaurus groups, the glaciations and the uplift of the Himalayas and Qinghai-Tibet plateau during the Pliocene-Pleistocene period might be the major factors affecting the present distribution of Eupetaurus along the Himalayas. (C) 2004 Elsevier Inc. All rights reserved.

Speciation and Phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: Analysis of the Cytochrome b Gene of mtDNA from Diverse Populations

Gao-Yan Li, Xu-Zhen Wang, Ya-Hui Zhao, Jie Zhang, Chun-Guang Zhang, and Shun-Ping He (2009) Speciation and phylogeography of Opsariichthys bidens (Pisces: Cypriniformes: Cyprinidae) in China: analysis of the cytochrome b gene of mtDNA from diverse populations. Zoological Studies 48(4): 569-583. The cyprinid fish Opsariichthys bidens Gunther is distributed in all major river systems of continental East Asia, and represents an attractive model for phylogeographic studies among cyprinid species or within a given species. In this study, we investigated the phylogeographic and demographic history of this species, using partial sequences of the cytochrome (cyt) b gene in mitochondrial (mt)DNA. Fish samples were collected from almost all major river systems where O. bidens is distributed in China. Sequence analysis showed remarkably high polymorphism, with 125 haplotypes in the 234 specimens examined, and with 89.8% of haplotypes occurring in only 1 specimen. A neutrality test indicated that some groups were not at mutation-drift equilibrium, suggesting a past population expansion. These results were supported by a mismatch distribution analysis. Based on our analysis, O. bidens consists of 4 groups belonging to 2 clades. The divergence time of the 2 clades was estimated to be 11.06-8.04 my. This value corresponds to the time of the 2nd uplift of the Qinghai-Tibet Plateau, the emergence of the East Asian monsoon, and the Epoch-6 Event. A two species scheme is proposed. http://zoolstud.sinica.edu.tw/Journals/48.4/569.pdf

The complete mitochondrial DNA sequence and the phylogenetic position of Achalinus meiguensis (Reptilia: Squamata)

The mitochondrial genome complete sequence of Achalinus meiguensis was reported for the first time in the present study. The complete mitochondrial genome of A. meiguensis is 17239 bp in length and contains 13 protein-coding genes, 22 tRNA, 2 rRNA, and 2 non-coding regions (Control regions). On the basis of comparison with the other complete mitochondrial sequences reported, we explored the characteristic of structure and evolution. For example, duplication control regions independently occurred in the evolutionary history of reptiles; the pseudo-tRNA of snakes occurred in the Caenophidia; snake is shorter than other vertebrates in the length of tRNA because of the truncations of T psi C arm (less than 5 bp) and "DHU" arm. The phylogenic analysis by MP and BI analysis showed that the phylogenetic position of A. meiguensis was placed in Caenophidia as a sister group to other advanced snakes with the exclusion of Acrochordus granulatus which was rooted in the Caenophidia. Therefore we suggested that the subfamily Xenodermatinae, which contains A. meiguensis, should be raised to a family rank or higher rank. At the same time, based on the phylogenic statistic test, the tree of Bayesian was used for estimating the divergence time. The results showed that the divergence time between Henophidia and Caenophidia was 109.50 Mya; 106.18 Mya for divergence between Acrochordus granulatus and the other snakes of the Caenophidia; the divergence time of A. meiguensis was 103 Mya, and Viperidae diverged from the unilateral of Elapidae and Colubridae was 96.06 Mya.

Phylogeny and speciation of the eastern Asian cyprinid genus Sarcocheilichthys

The genus Sarcocheilichthys is a group of small cyprinid fishes comprising 10 species/sub-species widely distributed in East Asia, which represents a valuable model for understanding the speciation of freshwater fishes in East Asia. In the present study, the molecular phylogenetic relationship of the genus Sarcocheilichthys was investigated using a 1140 bp section of the mitochondrial cytochrome b gene. Two different tree-building methods, maximum parsimony (MP) and Bayesian methods, yielded trees with almost the same topology, yielding high bootstrap values or posterior probabilities. The results showed that the genus Sarcocheilichthys consists of two large clades, clades I and II. Clade I contains Sarcocheilichthys lacustris, Sarcocheilichthys sinensis and Sarcocheilichthys parvus, with S. parvus at a basal position. In clade II, Sarcocheilichthys variegatus microoculus is at a basal position; samples of the widespread species, Sarcocheilichthys nigripinnis, form a large subclade containing another valid species Sarcocheilichthys czerskii. Sarcocheilichthys kiangsiensis is retained at an intermediate position. Since S. czerskii is a valid species in the S. nigripinnis clade, remaining samples of S. nigripinnis form a paraphyly. This speciation process is attributed to geographical isolation and special environmental conditions experienced by S. czerskii and stable environments experienced by the other S. nigripinnis populations. This type of speciation process was suggested to be very common. Samples of Sarcocheilichthys sinensis sinensis and Sarcocheilichthys sinensis fukiensis that did not form their own monophyletic groups suggest an early stage of speciation and support their sub-species status. Molecular clock analysis indicates that the two major lineages of the genus Sarcocheilichthys, clades I and II diverged c. 8.89 million years ago (mya). Sarcocheilichthys v. microoculus from Japan probably diverged 4.78 mya from the Chinese group. The northern-southern clades of S. nigripinnis began to diverge c. 2.12 mya, while one lineage of S. nigripinnis evolved into a new species, S. czerski, c. 0.34 mya. (C) 2008 The Authors Journal compilation (C) 2008 The Fisheries Society of the British Isles.

Molecular phylogeny of three subspecies of common carp Cyprinus carpio, based on sequence analysis of cytochrome b and control region of mtDNA

The complete cytochrome b and the control region of mtDNA (about 2070 bp in total) of 10 strains belonging to three subspecies of the common carp, including three wild subspecies (the Yangtze River wild common carp - Cyprinus carpio haematopterus, Yuanjiang River wild common carp Cyprinus carpio rubrofuscus and Volga River wild common carp - Cyprinus carpio carpio) and seven domestic strains (Xingguo red carp, Russian scattered scaled mirror carp, Qingtian carp, Japanese Koi carp, purse red carp, Big-belly carp, German mirror carp) were sequenced. Phylogenetic analysis indicated that the 10 strains form three distinct clades, corresponding to C. c. haematopterus, C. c. rubrofuscus and C. c. carpio respectively. Purse red carp, an endemic domestic strain in Jiangxi province of China, showed a higher evolution rate in comparison with the other strains of C. c. haematopterus, most probably because of intensive selection and a long history of domestication. Base variation ratios among the three subspecies varied from 0.78% (between C. c. haematopterus and C. c. rubrofuscus) to 1.47%(between C. c. carpio and C. c. rubrofuscus). The topography of the phylogenetic tree and the geographic distribution of three subspecies closely resemble each other. The divergence time between C. c. carpio and the other two subspecies was estimated to be about 0.9 Myr and about 0.5 Myr between C. c. haematopterus and C. c. rubrofuscus. Based on phylogenetic analysis, C. c. rubrofuscus might have diverged from C. c. haematopterus.