Mitochondrial DNA rates and biogeography in European newts (genus Euproctus)

Sequence divergence for segments of three mitochondrial DNA (mtDNA) genes encoding the 12S and 16S ribosomal RNA and cytochrome b was examined in newts belonging to the genus Euproctus (E. asper, E. montanus, E. platycephalus) and in three other species belonging to the same family (Salamandridae), Triturus carnifex, T. vulgaris, and Pleurodeles waltl. The three Euproctus species occur (one species each) in Corsica, Sardinia, and the Pyrenees. This vicariant distribution is believed to have been determined by the disjunction and rotation of the Sardinia-Corsica microplate from the Pyreneean region. Because time estimates are available for the tectonic events that led to the separation of the three landmasses, we used sequence data to estimate rates of evolution for the three gene fragments and investigated whether they conform to the rate-constancy hypothesis. By the Tajima (1993, Genetics 135:599-607) test, we could not detect rate heterogeneities for the ribosomal genes and for transversions in the cytochrome b gene. Assuming that these sites are evolving linearly over time and that cessation of gene flow occurred simultaneously with vicariant events, we compared the time of divergence estimated by molecular distances with the divergence times based on the geological estimates. Because we have two estimates of divergence time from the geological record, the split of Corsica/Sardinia from the Pyrenees and the split of Corsica from Sardinia, we could compare ratios of molecular divergence with the ratio of geological time divergence. The ratios are very similar, indicating that the molecular clock hypothesis cannot be rejected. These geological events also allowed us to calculate absolute rates of evolution for ribosomal and cytochrome b genes and compare them to rates for the same regions in other salamandrids and other vertebrates. Ribosomal mtDNA rates are comparable to those reported for other vertebrates, but cytochrome b rates are 3-7 times lower in salamanders than in other ectotherms. From a phylogenetic perspective, our data suggest that the cladogenic events leading to species formation in Euproctus and Triturus occurred very closely in time, indicating that the two genera may not be monophyletic. A duplication of the cytochrome b gene in T. carnifex was found, and the implications of this finding for mtDNA phylogenetic studies are discussed.     

Variability of nucleotide sequences of the mitochondrial DNA cytochrome c gene in dolly varden and taranetz char

Nucleotide sequence of the 307-bp fragment of the mitochondrial DNA cytochrome b gene was determined in representatives of the three species of the Salvelinus genus, specifically, dolly varden char (S. malma), taranetz char (S. taranetzi), and white-spotted char (S. leucomaenis). These results pointed to a high level of mitochondrial DNA (mtDNA) divergence between white-spotted char and dolly varden char, on the one hand, and taranetz char, on the other (the mean d value was 5.45%). However, the divergence between the dolly varden char and taranetz char was only 0.81%, which is comparable with the level of intraspecific divergence in the dolly varden char (d = 0.87%). It was shown that the dolly varden char mitochondrial gene pool contained DNA lineages differing from the main mtDNA pool at least in the taranetz char-specific mitochondrial lineages. One of these dolly varden char mtDNA lineages was characterized by the presence of the restriction endonuclease MspI-D variant of the cytochrome b gene. This lineage was widely distributed in the Chukotka populations but it was not detected in the Yana River (Okhotsk sea) populations. These findings suggest that dolly varden char has a more ancient evolutionary lineage, diverging from the common ancestor earlier than did taranetz char.     

Analysis of genetic diversity of domestic water buffaloes and anoas based on variations in the mitochondrial gene for cytochrome b

There are two major groups of domestic water buffaloes in East and Southeast Asia: swamp buffaloes and river buffaloes. Genetic diversity among swamp and river buffaloes was studied by DNA sequence analysis of the mitochondrial gene for cytochrome b. The results showed that each of the two groups has mitochondrial DNA (mtDNA) with a specific cytochrome b haplotype. The pairwise nucleotide sequence divergence was calculated to be 2·67% between swamp and river buffaloes, suggesting that they might have diverged from the ancestral populations of Asiatic domestic water buffaloes, approximately 1 million years ago. In addition, the sequences of the same gene from three subspecies of anoa (lowland, mountain and quarles anoa) were determined and compared with that of a domestic water buffalo. The sequence divergence was 1·2% for mountain anoa vs quarles anoa, 3·6% for mountain anoa vs lowland anoa and 3·3% for quarles anoa vs lowland anoa. Moreover, the sequence divergence between water buffaloes and anoas was found to be approximately 3·33%. Our results provide molecular evidence to support the taxonomic classification, namely, that Asiatic buffaloes may be classified into four lineages, swamp buffalo, river buffalo, lowland anoa and mountain plus quarles anoa. However, the sequence divergence values among these four groups were lower than the sequence divergence values found in the genus and subgenus levels within the subfamily Bovinae. In particular, in contrast to some proposed taxonomic classifications, our results indicated that mtDNA in the water buffaloes and anoas did not diverge at the genus level.     

Deep genetic divergences among morphologically similar and parapatric Skistodiaptomus (Copepoda: Calanoida: Diaptomidae) challenge the hypothesis of Pleistocene speciation

We used mitochondrial [cytochrome c oxidase subunit I (CO I ), cytochrome b , and 16S] and nuclear [internal transcribed spacer (ITS) phylogenies of Skistodiaptomus copepods to test hypotheses of Pleistocene divergence and speciation within the genus. Mitochondrial (mt)DNA sequence divergences do not support hypotheses for Pleistocene speciation and instead suggest much more ancient speciation events in the genus. Skistodiaptomus oregonensis and Skistodiaptomus pygmaeus (i.e. two morphologically similar and parapatric species) exhibited uncorrected mtDNA sequence divergences exceeding 20%. Similarly, we identified three divergent clades of Skistodiaptomus pallidus that exhibited mtDNA sequence divergences exceeding 15%, suggesting that even intraspecific divergence within this morphospecies predates the Pleistocene. We found clear evidence of CO I pseudogenes in S. pygmaeus , but their presence did not lead to significant overestimates of sequence divergences for this gene. Substitution saturation and strong purifying selection have most likely led to underestimates of sequence divergences and divergence times among Skistodiaptomus . The widespread phenomenon of morphological stasis among genetically divergent copepod groups indicates that speciation often occurs with little or no morphological change. Instead, morphological evolution may occur idiosyncratically after speciation and create discordant patterns of morphological similarity, shared ancestry and divergence time. Cryptic species complexes are therefore common in copepods, and morphological species concepts underestimate their true species diversity.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 150–165.     

SPECIATION BY POLYPLOIDY IN TREEFROGS: MULTIPLE ORIGINS OF THE TETRAPLOID,HYLA VERSICOLOR

Speciation by polyploidy is rare in animals, yet, in vertebrates, there is a disproportionate concentration of polyploid species in anuran amphibians. Sequences from the cytochrome b gene of the mitochondrial DNA (mtDNA) were used to determine phylogenetic relationships among 37 populations of the diploid-tetraploid species pair of gray treefrogs, Hyla chrysoscelis and Hyla versicolor. The diploid species, H. chrysoscelis, consists of an eastern and a western lineage that have 2.3% sequence divergence between them. The tetraploid species, H. versicolor, had at least three separate, independent origins. Two of the tetraploid lineages are more closely related to one or the other of the diploid lineages (0.18%–1.4% sequence divergence) than they are to each other (1.9%–3.4% sequence divergence). The maternal ancestor of the third tetraploid lineage is unknown. The phylogenetic relationships between the two species and among lineages within each species support the hypothesis of multiple origins of the tetraploid lineages.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the Hymenoptera

A higher AT content and rate of mtDNA sequence divergence was found in parasitic wasps (Apocrita) compared with nonparasitic wasps (Symphyta). The compositional bias was reflected in extreme codon bias for a cytochrome oxidase I protein coding gene fragment as well as in the types of amino acid substitutions that have occurred during the evolution of this gene fragment. In some instances, compositional bias influenced the definition of a conservative amino acid change. The increased rate of mtDNA sequence evolution probably arose during the early Jurassic, coincident with the first appearance of parasitic wasps in the fossil record. Our results suggest a causal link between the rate of sequence divergence and the parasitic lifestyle.Abbreviations AT adenosine-thymine - CO-1 cytochrome oxidase 1 - mtDNA mitochondrial DNA - Myr million years Correspondence to: M. Dowton     

The mtDNA sequence of the ostrich and the divergence between paleognathous and neognathous birds

The complete mitochondrial DNA (mtDNA) molecule of the ostrich, Struthio camelus, was sequenced. The size of the molecule is 16,591 nucleotides. Since the ostrich represents the paleognathous birds, comparison with the mtDNA of the neognathous chicken, the only avian species reported so far in databases, made it possible to identify common and, probably, general avian mtDNA characteristics. Relative to other vertebrates, the avian NADH6 and tRNA-Glu genes are positioned upstream of the control region rather than the cytochrome b gene. The NADH3 gene of the ostrich is terminated by a stop codon at position 207. Thus, the gene is about 140 nucleotides shorter than in other vertebrates. The sequence for L-strand origin of replication is missing in both birds, and four transfer RNA genes of the two avian mtDNAs deviate from common characteristics of tRNAs of vertebrate mtDNAs by having an adenine (and not a thymidine) at position 8. Due to the absence of suitable fossils, most paleontological datings of avian divergences are conjectural. Molecular dating of the divergence between the ostrich and the chicken indicates that these two avian lineages separated 80-90 MYA. Phylogenetic analysis of complete cytochrome b genes of six avian orders showed that Passeriformes represent the earliest divergence among recent birds, contradicting the commonly accepted notion of a basal position of the Palaeognathae among recent birds.      

Molecular analysis of an 11,700-year-old rodent midden from the Atacama Desert, Chile

DNA was extracted from an 11,700-year-old rodent midden from the Atacama Desert, Chile and the chloroplast and animal mitochondrial DNA (mtDNA) gene sequences were analysed to investigate the floral environment surrounding the midden, and the identity of the midden agent. The plant sequences, together with the macroscopic identifications, suggest the presence of 13 plant families and three orders that no longer exist today at the midden locality, and thus point to a much more diverse and humid climate 11,700 years ago. The mtDNA sequences suggest the presence of at least four different vertebrates, which have been putatively identified as a camelid (vicuna), two rodents (Phyllotis and Abrocoma), and a cardinal bird (Passeriformes). To identify the midden agent, DNA was extracted from pooled faecal pellets, three small overlapping fragments of the mitochondrial cytochrome b gene were amplified and multiple clones were sequenced. These results were analysed along with complete cytochrome b sequences for several modern Phyllotis species to place the midden sequence phylogenetically. The results identified the midden agent as belonging to an ancestral P. limatus. Today, P. limatus is not found at the midden locality but it can be found 100 km to the north, indicating at least a small range shift. The more extensive sampling of modern Phyllotis reinforces the suggestion that P. limatus is recently derived from a peripheral isolate.     

Phylogeography of the canyon treefrog, Hyla arenicolor (Cope) based on mitochondrial DNA sequence data

Patterns of phylogeography and gene flow were examined in the canyon treefrog, Hyla arenicolor. A total of 973 bp of mitochondrial cytochrome b sequence data were obtained for 65 individuals from 53 populations, yielding 50 unique haplotypes. Interpopulation sequence variation ranged from 0 to 13.7%. Phylogenetic analysis revealed three deeply divergent mtDNA lineages. These three Clades were mapped onto geography and found to represent completely concordant, nonoverlapping, geographical regions. Levels of sequence divergence between the three Clades were equal to or greater than published levels of divergence found in other vertebrate species and genera. Furthermore, one Clade of H. arenicolor was found to be more closely related to the outgroup H. eximia than to other H. arenicolor, suggesting that the taxonomy of this species may require revision.     

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.     

Phylogeography of the Alpine salamander, Salamandra atra (Salamandridae) and the influence of the Pleistocene climatic oscillations on population divergence

Fifty individuals of the endemic Alpine salamander, Salamandra atra, representing 13 populations throughout the range of the two currently recognized subspecies, atra and aurorae, were examined for sequence variation in a large portion (1050 bp) of the mitochondrial cytochrome b gene. We revealed a large number of mitochondrial DNA (mtDNA) haplotypes (10). Interpopulation sequence divergence was very low, ranging from 0 to 3.1%. The relationships among haplotypes were poorly resolved. The divergence time estimate between several mtDNA haplotypes suggested a pre-Pleistocene differentiation approximately 3 million years ago. Moreover, the impact of the Pleistocene glaciations on the phylogeographical patterns appears to have been secondary, although a somewhat reduced genetic variability was found in populations living in areas that were directly affected by the glaciation.     

Molecular evolution of mitochondrial 12S RNA and cytochrome b sequences in the pantherine lineage of Felidae

DNA sequence comparisons of two mitochondrial DNA genes were used to infer phylogenetic relationships among 17 Felidae species, notably 15 in the previously described pantherine lineage. The polymerase chain reaction (PCR) was used to generate sequences of 358 base pairs of the mitochondrial 12S RNA gene and 289 base pairs of the cytochrome b protein coding gene. DNA sequences were compared within and between 17 felid and five nonfelid carnivore species. Evolutionary trees were constructed using phenetic, cladistic, and maximum likelihood algorithms. The combined results suggested several phylogenetic relationships including (1) the recognition of a recently evolved monophyletic genus Panthera consisting of Panthera leo, P. pardus, P. onca, P. uncia, P. tigris, and Neofelis nebulosa; (2) the recent common ancestry of Acinonyx jubatus, the African cheetah, and Puma concolor, the American puma; and (3) two golden cat species, Profelis temmincki and Profelis aurata, are not sister species, and the latter is strongly associated with Caracal caracal. These data add to the growing database of vertebrate mtDNA sequences and, given the relatively recent divergence among the felids represented here (1-10 Myr), allow 12S and cytochrome b sequence evolution to be addressed over a time scale different from those addressed in most work on vertebrate mtDNA.      

Molecules and morphology: evidence for cryptic hybridization in African Hyalomma (Acari: Ixodidae)

The role of natural hybridization and introgression as part of the evolutionary process is of increasing interest to zoologists, particularly as more examples of gene exchange among species are identified. We present mitochondrial and nuclear sequence data for Hyalomma dromedarii, Hyalomma truncatum, and Hyalomma marginatum rufipes (Acari: Ixodidae) collected from one-humped camels in Ethiopia. These species are well differentiated morphologically and genetically; sequence data from the mitochondrial DNA (mtDNA) cytochrome oxidase I gene indicates 10-14% divergence between the species. However, incongruence between morphology and the mtDNA phylogeny was observed, with multiple individuals of H. dromedarii and H. truncatum present on the same mtDNA lineage as H. marginatum rufipes. Thus, individuals with morphology of H. dromedarii and H. truncatum are indistinguishable from H. marginatum rufipes on the basis of mtDNA. Multiple copies of ITS-2 were subsequently cloned and sequenced for a subset of individuals from the mtDNA phylogeny, representing both 'normal' and 'putative hybrid' individuals. Very low sequence divergence (0.3%) was observed within 'normal' individuals of both H. dromedarii and H. truncatum relative to the 'putative hybrid' individuals (6 and 2.7%, respectively). The pattern of intra-individual variation in ITS-2 within 'putative hybrid' individuals, particularly in H. dromedarii, strongly suggests that gene flow has occurred among these Hyalomma species, but no indication of this is given by the morphology of the individuals.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

Unusual mitochondrial DNA polymorphism in two local populations of blue tit Parus caeruleus

Mitochondrial DNA (mtDNA) from 25 blue tits Parus caeruleus sampled from two populations of the Grenoble region (France) was assayed for polymorphism with 17 restriction endonucleases. Nine genotypes were found. Several mtDNA genotypes were also analysed by amplification via the polymerase chain reaction (PCR) and direct sequencing of 903 bp of the cytochrome b gene. The mtDNA polymorphism is greater in P. caeruleus than in other comparable bird species and results from the presence of two clearly differentiated mitochondrial lineages. Using the data of restriction polymorphism, the mean sequence divergence between individuals of the two lineages is 1.23%. Therefore, P. caeruleus should fall into the category II of phylogeographic pattern sensu Avise et al. (1987): discontinuous mtDNA genotypes which co-occur in the same region. P. caeruleus, like humans and other mobile species with high gene flow, seems to have lost its geographic structure in terms of mtDNA phylogeny. This unusual mitochondrial polymorphism can be explained by the recent admixture of two long-term isolated populations. This could be accounted for by two different scenarios. One assumes a simultaneous post-glacial colonization of the Grenoble region by two isolated European populations of P. caeruleus. Alternatively, hybridization between P. caeruleus and P. cyanus could have caused the observed pattern of mtDNA variation.     

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.     

The largemouth bass cytochrome b gene

The cytochrome b gene is one of the protein-coding genes of the mitochondrial genome that has gained importance because of the ease with which molecular techniques can be applied to the analysis of its structure. The nucleotide sequence of the largemouth bass ( Microplerus salmoides ) cytochrome b gene was determined and the inferred amino acid structure is presented in the form of a structural model derived originally from rat cytochrome b . The inferred amino acid sequences from divergent animal species are aligned and compared in the context of this model. The data suggest that regions of the gene may be evolving at different rates due to different selection pressures associated with functional constraints. Conserved and variable regions of cytochrome b have been identified and can be targeted for species identification, the examination of intraspeciflc variation, and phylogenetic reconstructions in future research.     

Mitochondrial diversity of Opsariichthys bidens (Teleostei, Cyprinidae) in three Chinese drainages

We describe the phylogeographic structure of 28 Chinese populations of the cyprinid Opsariichthys bidens across three main Chinese river drainages. Our study is based on the phylogenetic analysis of the complete mitochondrial cytochrome b gene (1140 bp). We combined this analysis with population processes inferred from nested clade analysis (NCA) and mismatch distributions. Both analyses showed that Chinese O. bidens consists of five mtDNA lineages (Opsariichthys 1-5) with high genetic divergence among them. Molecular divergences (TrN+G) higher than 20% among the Opsariichthys 1-5 mtDNA lineages suggest a taxonomic underestimation at the species level. About 92% of the genetic variance among samples was explained by differences among Opsariichthys mtDNA lineages. Drainage-restricted haplotypes with high frequencies and moderate nucleotide diversity show that Opsariichthys populations have evolved independently. NCA results were congruent with the phylogeny, and unimodal mismatch distributions with negative Tajima's D values suggest population expansions in some Opsariichthys lineages. The phylogeographic structure of the Opsariichthys 1-5 mtDNA lineages appears to be related to their long-term interruption of gene flow (theta(ST)>0.97). Our results suggested that fragmentation of ancestral ranges might have caused Opsariichthys diversification in Chinese waters. However, current distribution of common haplotypes across the Yangtze and Pearl drainages suggests a recent river connection that could have favoured gene flow across drainages. Overall, the results indicated that the richness of current Asian widespread species might have been underestimated, and that the cyprinid populations of O. bidens in the Yangtze, Pearl and Hai He drainages may correspond to five species.