Comparative Phylogenetic Analyses of the Adaptive Radiation of Lake Tanganyika Cichlid Fish: Nuclear Sequences Are Less Homoplasious But Also Less Informative Than Mitochondrial DNA

Over 200 described endemic species make up the adaptive radiation of cichlids in Lake Tanga-nyika. This species assemblage has been viewed as both an evolutionary reservoir of old cichlid lineages and an evolutionary hotspot from which the modern cichlid lineages arose, seeding the adaptive radiations in Lakes Victoria and Malawi. Here we report on a phylogenetic analysis of Lake Tanganyika cichlids combining the previously determined sequences of the mitochondrial ND2 gene (1047 bp) with newly derived sequences of the nuclear RAG1 gene (∼700 bp of intron 2 and ∼1100 bp of exon 3). The nuclear data—in agreement with mitochondrial DNA—suggest that Lake Tanganyika harbors several ancient lineages that did not undergo rampant speciation (e.g., Bathybatini, Trematocarini). We find strong support for the monophyly of the most species-rich Tanganyikan group, the Lamprologini, and we propose a new taxonomic group that we term the C-lineage. The Haplochromini and Tropheini both have an 11-bp deletion in the intron of RAG1, strongly supporting the monophyly of this clade and its derived position. Mapping the phylogenetically informative positions revealed that, for certain branches, there are six times fewer apomorphies in RAG1. However, the consistency index of these positions is higher compared to the mitochondrial ND2 gene. Nuclear data therefore provide, on a per–base pair basis, less but more reliable phylogenetic information. Even if in our case RAG1 has not provided as much phylogenetic information as we expected, we suggest that this marker might be useful in the resolution of the phylogeny of older groups. Reviewing Editor: Dr. Rafael Zardoya     

基于ND4和ND5基因序列分析的鳅超科鱼类系统发育关系

ND4和ND5是线粒体基因组中编码NADH脱氢酶亚基4和亚基5的两个蛋白质编码基因.该研究以鳅超科鱼类为研究对象,新测定了10个物种的ND4和ND5基因全序列以及中间的3个tRNA基因共212 bp的序列,结合从GenBank 下载的15个物种的15条序列进行序列比较和系统发育关系分析.结果显示:鳅超科鱼类ND4基因全长1380～1387 bp,以ATG为起始密码子,终止密码子为不完全终止信号;ND5基因全长1821～1839bp,同样起始密码子为ATG,终止密码子为TAA或TAG;ND4和ND5基因之间插入了3个tRNA基因,分别编码携带组氨酸、丝氨酸、亮氨酸的tRNA.ND4和ND5基因(包含3个tRNA基因)中A、T、G、C的平均含量分别为30.4%、27.3%、14.2%、28.1%,A+T(57.7%)的含量高于G+C(42.3%)的含量.转换与颠换比(Ti/Tv)平均值为1.586.选取斑马鱼和鲤鱼作为外类群,采用最大简约法(MP)、最大似然法(ML)和贝叶斯推断法(BI)进行系统发育树的重建.三种方法的系统发育分析结果都显示:花鳅亚科、条鳅亚科、沙鳅亚科、平鳍鳅科及Vaillantellidae分别构成单系;它们的系统发育关系为:(Vaillantellidae+(沙鳅亚科+(花鳅亚科+(条鳅亚科+平鳍鳅科).这与线粒体全基因组和某些核基因(如RAG1基因)的研究结果类似,且支持率较高,表明ND4和ND5基因用于鳅超科鱼类的系统发育分析是可行的;但是该研究的结果有别于其他线粒体基因的分析结果,如基于cytb和D-loop基因进行的系统发育分析表明,条鳅亚科和花鳅亚科聚为姐妹群,再和平鳍鳅科聚在一起.这种差异可能是由于使用的基因长度差异造成的,长度越长,信息量越大,所反映的系统发育结果可能更加接近真实情况.     

华北地区中华蒙潮虫(甲壳动物亚门:等足目)种群遗传多样性研究

中华蒙潮虫Mongoloniscus sinensis(Dollfus,1901)隶属于甲壳动物亚门Crustacea等足目Isopoda潮虫亚目Oniscidea,中国特有种。为了探究中华蒙潮虫的种群遗传分化和系统进化关系,采用PCR对采自华北地区10个地理种群89只个体线粒体2个基因COⅠ和ND5进行联合分析。结果表明:1)中华蒙潮虫COⅠ部分基因长604 bp,ND5部分基因长615 bp,拼接序列长1 219 bp,T、C、A和G含量分别为41.0%、11.2%、30.8%和17.0%,具有显著的A+T偏倚;变异位点503个(占总核苷酸序列的41.3%),序列间的转换/颠换比值为2.8。2)89只个体共45种单倍型,单倍型多样性0.964,核苷酸多样性0.005 6,整体遗传多样性水平中等;单倍型H1、H15、H16、H21、H41为2～3个种群共享单倍型。3)联合基因(COⅠ+ND5)系统发育树表明,最早出现的是华北以北地区(山西大同、河北石家庄),最晚分化出的是华北以南地区(山西临汾、陕西西安未央区、河南新乡),演化路线为从北向南,个别种群单倍型未按地理来源形成明显的簇群。4)平均遗传分化指数为0.513,基因流为0.24;分子变异分析结果表明,种群的变异与分化主要来自种群内部,错配分布呈多峰,结合中性检验(Tajima's D=-1.429;Fu's F_s=6.499),发现中华蒙潮虫近期未经历扩张,但种群内部分化显著,增长平稳。本研究首次基于线粒体多基因联合分析了中华蒙潮虫种群遗传多样性。     

中华雏蝗（Chorthippus chinensis Tarb）线粒体基因组分析

采用Lon-PCR扩增线粒体全基因组和保守引物步移法结合克隆方法测定并拼接获得了中华雏蝗（Chorthippus chinensis Tarb）线粒体基因组全序列．序列的注释和分析结果表 明,中华雏蝗线粒体基因组序列全长15 599 bp,共有13个编码蛋白质基因、22个tRNA基因、2个rRNA基因和1个A+T富集区．基因顺序与非洲飞蝗（Locusta migratoria）相同,也发生了2个 tRNA Asp(D)和tRNALys(K)的倒置．13个编码蛋白质基因都使用了ATN作为起始密码子．除ND1以TAG和ND5的终止密码子为不完全的T外,其余11个编码蛋白质基因的终止密码子都为完整的TAA．6种直翅类昆虫13个蛋白质的氨基酸序列的联合数据集构建的系统树与形态分类系统一致,中华雏蝗与非洲飞蝗为姐妹群,并与东方蝼蛄构成一单系群．     

Molecular phylogeny of hipposiderid bats (Chiroptera: Hipposideridae) from Solomon Islands and Cape York Peninsula,Australia

Throughout the Old World Tropics, examples of cryptic diversity have been demonstrated in leaf‐nosed bats from the Family Hipposideridae. In addition to cryptic diversity, the widespread distributions and graded morphologies of some species render the taxonomy of this family poorly resolved. We used sequences of the mitochondrial gene ND2 and the nuclear gene RAG1 to reconstruct molecular phylogenetic relationships for hipposiderid species. Specimens from Solomon Island and Cape York Peninsula, Australia, were sequenced and combined with published sequences from Southeast Asia. Our results suggest the presence of undescribed diversity in northern Australia. Sequences of ND2 and RAG1 for samples of Cape York Hipposideros ater were distinct from those of Southeast Asia, with corrected sequence divergence estimates for the mitochondrial gene exceeding 23%. In comparison, there was no cryptic diversity detected in Solomon Island specimens when compared to conspecifics from Southeast Asia. On the contrary, three species members of the diadema morphological group from the region exhibit variability in the mitochondrial gene ND2, far below that of other species included in the analyses.     

Molecular phylogeny of Japanese Catocala moths based on nucleotide sequences of the mitochondrial ND5 gene

Phylogenetic relationships of 31 Japanese Catocala species were analyzed based on the partial nucleotide sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene (762 bp). When several non-Catocala Noctuidae moths were designated as the outgroup, these Catocala species formed a monophyletic group. However, divergences between these Catocala species were very deep, and no close phylogenetic relationships were recognized among them except for that between the two recently separated species, C. xarippe and C. fulminea. The remote relationships implied for several pairs of species suggest that the color of the hindwings is a changeable characteristic, and does not reflect phylogenetic lineage. Continental specimens were analyzed in 20 of 31 Catocala species, and all of them showed a close relationship with their Japanese counterpart. However, the closeness of the nucleotide sequences between the Japanese and continental individuals of the same species varied from species to species, indicating that isolation between the Japanese and continental populations of these species occurred at many different times. The two analyzed species endemic to North America showed a close relationship with their morphologically inferred Japanese counterparts, indicating that the geographic separation and following speciation between these Eurasian and American species occurred much more recently compared with the speciation events among the Catocala species now found in Japan.     

Phylogenetic relationships of the Cobitoidea (Teleostei: Cypriniformes) inferred from mitochondrial and nuclear genes with analyses of gene evolution

The superfamily Cobitoidea of the order Cypriniformes is a diverse group of fishes, inhabiting freshwater ecosystems across Eurasia and North Africa. The phylogenetic relationships of this well-corroborated natural group and diverse clade are critical to not only informing scientific communities of the phylogeny of the order Cypriniformes, the world's largest freshwater fish order, but are key to every area of comparative biology examining the evolution of traits, functional structures, and breeding behaviors to their biogeographic histories, speciation, anagenetic divergence, and divergence time estimates. In the present study, two mitochondrial gene sequences (COI, ND4+5) and four single-copy nuclear gene segments (RH1, RAG1, EGR2B, IRBP) were used to infer the phylogenetic relationships of the Cobitoidea as reconstructed from maximum likelihood (ML) and partitioned Bayesian Analysis (BA). Analyses of the combined mitochondrial/nuclear gene datasets revealed five strongly supported monophyletic Cobitoidea families and their sister-group relationships: Botiidae+(Vaillantellidae+(Cobitidae+(Nemacheilidae+Balitoridae))). These recovered relationships are in agreement with previous systematic studies on the order Cypriniformes and/or those focusing on the superfamily Cobitoidea. Using these relationships, our analyses revealed pattern lineage- or ecological-group-specific evolution of these genes for the Cobitoidea. These observations and results corroborate the hypothesis that these group-specific-ancestral ecological characters have contributed in the diversification and/or adaptations within these groups. Positive selections were detected in RH1 of nemacheilids and in RAG1 of nemacheilids and genus Vaillantella, which indicated that evolution of RH1 (related to eye's optic sense) and RAG1 (related to immunity) genes appeared to be important for the diversification of these groups. The balitorid lineage (those species inhabiting fast-flowing riverine habitats) had, as compared with other cobitoid lineages, significantly different dN/dS, dN and dS values for ND4 and IRBP genes. These significant differences are usually indicative of weaker selection pressure, and lineage-specific evolution on genes along the balitorid lineage. Furthermore, within Cobitoidea, excluding balitorids, species living in subtropics had significantly higher dN/dS values in RAG1 and IRBP genes than those living in temperate and tropical zones. Among tropical cobitoids, genes COI, ND5, EGR2B, IRBP and RH1, had a significantly higher mean dS value than those species in subtropical and temperate groups. These findings suggest that the evolution of these genes could also be ecological-group-specific and may have played an important role in the adaptive evolution and diversification of these groups. Thus, we hypothesize that the genes included in the present study were actively involved in lineage- and/or ecological-group-specific evolutionary processes of the highly diverse Cobitoidea. These two evolutionary patterns, both subject to further testing, are hypothesized as integral in the diversification with this major clade of the world's most diverse group of freshwater fishes.     

Complete mitochondrial genome of the boreal digging frog Kaloula borealis (Anura, Microhylidae)

We sequenced the complete mitochondrial genome from the boreal digging frog Kaloula borealis. The genome sequence was 17,173 bp in size, and the gene order and contents were identical to those of previously reported amphibian mitochondrial genomes. Of 13 protein-coding genes (PCGs), 5 genes (CO2, ATPase 6, CO3, ND3, and ND4) had incomplete stop codons. Also ND1 gene used GTG as a start codon, while CO1 and ND5 genes used AGG as a stop codon. The base composition of K. borealis mitogenome showed a strong anti-G bias (6.11%) on the 3rd position of PCGs.     

长蠹科几种检疫性害虫的ND4基因序列及系统进化

长蠹科昆虫严重危害林木和仓贮物品。应用非损伤性DNA测序技术测定了来自不同国家的长蠹科害虫的线粒体DNA ND4 基因的部分序列。在获得的204 bp的序列中,5种昆虫的序列变异丰富,多数变异发生在密码子的第3位点上。用PAUP3.1.1数据分析软件构建了5个种的合意简约树。并将实验结果与形态学特征比较分析,探讨5个种及所在属的系统进化。结果表明：双钩异翅长蠹所在的异翅长蠹属分化最早,其次是竹大长蠹所在的大长蠹属、双棘长蠹和黑双棘长蠹所在的双棘长蠹属及红艳长蠹所在的钻木长蠹属。双棘长蠹和黑双棘长蠹隶属同一个属,遗传关系最近,分化最晚,与形态学研究结果相吻合。     

基于核基因RAG2部分序列探索菊头蝠总科、蝙蝠科和鞘尾蝠科物种的系统发育关系

采用PCR技术获得了翼手目(菊头蝠总科、蝙蝠科以及鞘尾蝠科)16种19个个体的核基因重组激活基因2(Recombination activating gene,RAG2)部分序列,长度为730～760 bp.结合从Genbank中提取的翼手目3科14个体的RAG2序列,通过构建贝叶斯(Bayesian inference,BI)和最大似然树对翼手目种属间的系统进化关系进行研究.研究表明:在菊头蝠总科中,菊头蝠科和蹄蝠科是两个独立的科,且在蹄蝠科中,大蹄蝠、中蹄蝠、普氏蹄蝠3种之间的亲缘关系非常近;蝙蝠科中的南蝠属与棕蝠属是姊妹群;长翼蝠亚科应提升为长翼蝠科.此外,鞘尾蝠科与犬吻蝠科形成姊妹群关系.     

Mechanical barriers to introgressive hybridization revealed by mitochondrial introgression patterns in Ohomopterus ground beetle assemblages

To reveal the role of diverged body size and genital morphology in reproductive isolation among closely related species, we examined patterns of, and factors limiting, introgressive hybridization between sympatric Ohomopterus ground beetles in central Japan using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. We sampled 17 local assemblages that consisted of two to five species and estimated levels of interspecific gene flow using the genetic distance, D(A), and maximum-likelihood estimates of gene flow. Sharing of haplotypes or haplotype lineages was detected between six of seven species that occurred in the study areas, indicating mitochondrial introgression. The intensity and direction of mitochondrial gene flow were variable among species pairs. To determine the factors affecting introgression patterns, we tested the relationships between interspecific D(A) and five independent variables: difference in body size, difference in genital size, phylogenetic relatedness (nuclear gene sequence divergence), habitat difference, and species richness of the assemblage. Body and genital size differences contributed significantly to preventing gene flow. Thus, mechanical isolation mechanisms reduce the chance of introgressive hybridization between closely related species. Our results highlight the role of morphological divergence in speciation and assemblage formation processes through mechanical isolation.     

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.     

地中海实蝇幼虫分子鉴定

通过对从秘鲁进口的葡萄中截获实蝇类幼虫进行ITS区和线粒体COⅠ,COⅡ,COⅢ、ND5基因序列的扩增和测序,并与GenBank中对应的序列进行比对,结果表明,截获样品ITS区序列和地中海实蝇Ceratitis capitata(Wiedemann)同源性为95·16%(其中ITS1为99·52%,ITS2为86·2%),线粒体COⅠ,COⅡ,COⅢ,ND5基因序列和地中海实蝇C.capitata同源性为100%,99·9%,99·5%,99·8%;基于COⅠ序列构建的系统发育树中,幼虫样品和地中海实蝇最为接近。根据序列分析和系统发育关系分析的结果,将截获的实蝇类幼虫鉴定为地中海实蝇C.capitata。     

Geographic mode of speciation in a mountain specialist Avian family endemic to the Palearctic

Mountains host greater avian diversity than lowlands at the same latitude due to their greater diversity of habitats stratified along an elevation gradient. Here we test whether this greater ecological heterogeneity promotes sympatric speciation. We selected accentors (Prunellidae), an avian family associated with mountains of the Palearctic, as a model system. Accentors differ in their habitat/elevation preferences and south‐central Siberia and Himalayan regions each host 6 of the 13 species in the family. We used sequences of the mtDNA ND2 gene and the intron 9 of the Z chromosome specific ACO1 gene to reconstruct a complete species‐level phylogeny of Prunellidae. The tree based on joint analysis of both loci was used to reconstruct the family's biogeographic history and to date the diversification events. We also analyzed the relationship between the node age and sympatry, to determine the geographic mode of speciation in Prunellidae. Our data suggest a Miocene origin of Prunellidae in the Himalayan region. The major division between alpine species (subgenus Laiscopus) and species associated with shrubs (subgenus Prunella) and initial diversification events within the latter happened within the Himalayan region in the Miocene and Pliocene. Accentors colonized other parts of the Palearctic during the Pliocene‐Pleistocene transition. This spread across the Palearctic resulted in rapid diversification of accentors. With only a single exception dating to 0.91 Ma, lineages younger than 1.5 Ma are allopatric. In contrast, sympatry values for older nodes are >0. There was no relationship between node age and range symmetry. Allopatric speciation (not to include peripatric) is the predominant geographic mode of speciation in Prunellidae despite the favorable conditions for ecological diversification in the mountains and range overlaps among species.     

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.     

Rejecting strictly allopatric speciation on a continental island: prolonged postdivergence gene flow between Taiwan (Leucodioptron taewanus,Passeriformes Timaliidae) and Chinese (L. canorum canorum) hwameis

Allopatry is conventionally considered the geographical mode of speciation for continental island organisms. However, strictly allopatric speciation models that assume the lack of postdivergence gene flow seem oversimplified given the recurrence of land bridges during glacial periods since the late Pliocene. Here, to evaluate whether a continental island endemic, the Taiwan hwamei (Leucodioptron taewanus, Passeriformes Timaliidae) speciated in strict allopatry, we used weighted‐regression‐based approximate Bayesian computation (ABC) to analyse the genetic polymorphism of 18 neutral nuclear loci (total length: 8500 bp) in Taiwan hwamei and its continental sister species, the Chinese hwamei (L. canorum canorum). The nonallopatry model was found to fit better with observed genetic polymorphism of the two hwamei species (posterior possibility = 0.82). We also recovered unambiguous signals of nontrivial bidirectional postdivergence gene flow (N_em » 1) between Chinese hwamei and Taiwan hwamei until 0.5 Ma. Divergence time was estimated to be 3.5 to 2 million years earlier than that estimated from mitochondrial cytochrome b sequences. Finally, using the inferred nonallopatry model to simulate genetic variation at 24 nuclear genes examined showed that the adiponectin receptor 1 gene may be under divergent adaptation. Our findings imply that the role of geographical barrier may be less prominent for the speciation of continental island endemics, and suggest a shift in speciation studies from simply correlating geographical barrier and genetic divergence to examining factors that facilitate and maintain divergence, e.g. differential selection and sexual selection, especially in the face of interpopulation gene flow.     

西北太平洋海岸带大弹涂鱼复合体的隐存种与进化历史

大弹涂鱼Boleophthalmus pectinirostris间断分布于西太平洋海岸带东亚与马来西亚马六甲海峡,但马来西亚种群的分类地位尚存争议。研究使用线粒体ND5基因序列（718 bp）与核位点Rag1基因序列（1395 bp）对西北太平洋海岸带11个地点的45尾大弹涂鱼属鱼类进行系统发育关系重建,结果表明大弹涂鱼包括东亚与马来西亚两个单系群,两者形成姊妹群关系。GMYC分析、*BEAST物种树支持大弹涂鱼东亚谱系和马来西亚谱系是不同种。分子测定年龄分析表明大弹涂鱼东亚谱系与马来西亚谱系之间的分化时间为2.73百万年。因此,西北太平洋海岸带大弹涂鱼是复合体,包括两个物种:东亚种群是大弹涂鱼Boleophthalmus pectinirostris sensu stricto,而马来西亚种群是隐存种Boleophthalmus sp.。大弹涂鱼与隐存种之间的物种分化可能是晚上新世冰期海平面下降产生的地理隔离以及间冰期洋流对基因交流的阻碍两方面相互作用的结果。       

The youngest split in sympatric schizothoracine fish (Cyprinidae) is shaped by ecological adaptations in a Tibetan Plateau glacier lake

Although new empirical evidence shows that sympatric speciation has occurred in some species, there are few indisputable model organisms for this process of speciation. The two subspecies ( Gymnocypris eckloni eckloni and G . e . scoliostomus ) of the schizothoracine Gymnocypris fish species complex from a small glacier lake in the Tibetan Plateau, Lake Sunmcuo, fit several of the key characteristics of the sympatric speciation model. We used combined mitochondrial control region sequences and the cytochrome b gene (1894 bp) to address the phylogenetics and population genetics of 232 specimens of G . e . eckloni and G . e . scoliostomus , as well as all of its closely related sister species. We found that: (i) a total of four old lineages were uncovered in the widespread G . e . eckloni , of which only one was shown to be shared with all G . e . scoliostomus individuals and (ii) the new subspecies ( G . e . scoliostomus ) evolved in Lake Sunmcuo from the ancestral G . e . eckloni population within approximately 0.057 Ma. These two taxa of the species complex are morphologically distinct, and reproductive isolation is further suggested. Ecological disruptive selection based on morphological traits (e.g. mouth cleft characters) and food utilization may be a mechanism of incipient speciation of two sympatric populations within Lake Sunmcuo. This study provides the first genetic evidence for sympatric speciation in the schizothoracine fish.     

Molecular phylogenetic relationships among species of the Malagasy-Comoran gecko genus Paroedura (Squamata: Gekkonidae)

We use approximately 3100bp of mitochondrial (ND2, ND4) and nuclear (RAG1, phosducin) DNA sequence data to recover phylogenetic relationships among 14 of the 16 recognized taxa of the lizard genus Paroedura as well as two undescribed forms. These geckos are endemic to Madagascar and the Comores and are popularly kept and bred by herpetoculturalists. The closest relative of Paroedura is another Indian Ocean leaf-toed gecko, Ebenavia. Both Bayesian inference and maximum parsimony strongly support the monophyly of two major clades within Paroedura that conflict with existing species group assignments based on scale characteristics. Our well-resolved tree elucidates a biogeographic pattern in which eastern Paroedura are most basal and western and south-western species form a monophyletic group. Our data demonstrate the phylogenetic utility of phosducin, a novel marker in squamate phylogenetics, at the intrageneric level.