Variation in salmonid mitochondrial DNA: evolutionary constraints and mechanisms of substitution

Waves Out of the Korean Peninsula and Inter- and Intra-Species Replacements in Freshwater Fishes in Japan

The Japanese archipelago is located at the periphery of the continent of Asia. Rivers in the Japanese archipelago, separated from the continent of Asia by about 17 Ma, have experienced an intermittent exchange of freshwater fish taxa through a narrow land bridge generated by lowered sea level. As the Korean Peninsula and Japanese archipelago were not covered by an ice sheet during glacial periods, phylogeographical analyses in this region can trace the history of biota that were, for a long time, beyond the last glacial maximum. In this study, we analyzed the phylogeography of four freshwater fish taxa, Hemibarbus longirostris, dark chub Nipponocypris temminckii, Tanakia ssp. and Carassius ssp., whose distributions include both the Korean Peninsula and Western Japan. We found for each taxon that a small component of diverse Korean clades of freshwater fishes migrated in waves into the Japanese archipelago to form the current phylogeographic structure of biota. The replacements of indigenous populations by succeeding migrants may have also influenced the phylogeography.

GLOBAL SURVEY OF MITOCHONDRIAL DNA SEQUENCES IN THE THREESPINE STICKLEBACK: EVIDENCE FOR RECENT MIGRATIONS

Phylogenetic analyses of mitochondrial DNA (mtDNA) sequences were used to assess the matriarchal genetic structure of the threespine stickleback, Gasterosteus aculeatus. A 747 base-pair (bp) fragment of the cytochrome b was sequenced from 36 individuals collected from 25 localities in Europe, North America, and Japan. Two major divergent clades were revealed: one widespread in Japan but with representatives in some Alaskan and British Columbian lakes and the other common in Europe and North America. A simple diagnostic test using the polymerase chain reaction (PCR) and a restriction enzyme was used to assay additional individuals, confirming the absence of the Japanese clade in the Atlantic basin. Geographic distribution of mtDNA variation suggests (1) a recent origin of the Atlantic populations, and (2) support for previous hypotheses about the existence of Pleistocene refugia for freshwater fishes in Alaska and British Columbia. Silent substitution rates were used to date the colonization of the Atlantic at 90,000 to 260,000 yr before present, which conflicts with earlier dates implied by the fossil record. The recent replacement of Atlantic mitochondrial lineages suggested by our data may be explained by severe reduction or extinction of northern Atlantic populations during the Pleistocene, followed by a recent reinvasion from the Pacific. With a global perspective of the distribution of genetic variation as a framework, meaningful comparisons at a smaller geographical scale will now be possible.

INFERRING PHYLOGENIES FROM mtDNA VARIATION: MITOCHONDRIAL-GENE TREES VERSUS NUCLEAR-GENE TREES

An accurately resolved gene tree may not be congruent with the species tree because of lineage sorting of ancestral polymorphisms. DNA sequences from the mitochondrially encoded genes (mtDNA) are attractive sources of characters for estimating the phylogenies of recently evolved taxa because mtDNA evolves rapidly, but its utility is limited because the mitochondrial genes are inherited as a single linkage group (haplotype) and provide only one independent estimate of the species tree. In contrast, a set of nuclear genes can be selected from distinct chromosomes, such that each gene tree provides an independent estimate of the species tree. Another aspect of the gene-tree versus species-tree problem, however, favors the use of mtDNA for inferring species trees. For a three-species segment of a phylogeny, the branching order of a gene tree will correspond to that of the species tree if coalescence of the alleles or haplotypes occurred in the internode between the first and second bifurcation. From neutral theory, it is apparent that the probability of coalescence increases as effective population size decreases. Because the mitochondrial genome is maternally inherited and effectively haploid, its effective population size is one-fourth that of a nuclear-autosomal gene. Thus, the mitochondrial-haplotype tree has a substantially higher probability of accurately tracking a short internode than does a nuclear-autosomal-gene tree. When an internode is sufficiently long that the probability that the mitochondrial-haplotype tree will be congruent with the species tree is 0.95, the probability that a nuclear-autosomalgene tree will be congruent is only 0.62. If each of k independently sampled nuclear-gene trees has a probability of congruence with the species tree of 0.62, then a sample of 16 such trees would be required to be as confident of the inference based on the mitochondrial-haplotype tree. A survey of mtDNA-haplotype diversity in 34 species of birds indicates that coalescence is generally very recent, which suggests that coalescence times are typically much shorter than internodal branch lengths of the species tree, and that sorting of mtDNA lineages is not likely to confound the species tree. Hybridization resulting in transfer of mtDNA haplotypes among branches could also result in a haplotype tree that is incongruent with the species tree; if undetected, this could confound the species tree. However, hybridization is usually easy to detect and should be incorporated in the historical narrative of the group, because reticulation, as well as cladistic events, contributed to the evolution of the group.

ALLOPATRIC ORIGIN OF SYMPATRIC POPULATIONS OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS) AS REVEALED BY MITOCHONDRIAL-DNA RESTRICTION ANALYSIS

In the paper, restriction-fragment length polymorphisms in mitochondrial DNA (mtDNA) were studied to test the hypothesis that sympatric populations of lake whitefish in the Allegash basin have recently diverged through sympatric speciation. Thirteen restriction enzymes were used to analyze mtDNA of 156 specimens representing 13 populations from eastern Canada and northern Maine where normal and dwarf phenotypes of whitefish exist in sympatry and allopatry. Two monophyletic assemblages of populations that exhibit different geographic distributions were identified. One showed an eastern distribution that expands from Cape Breton to the Allegash basin and the other exhibits a more western distribution. The Allegash basin was the only area of overlap. The western assemblage exhibited the normal size phenotype in all cases, whereas the eastern assemblage exhibited the normal size phenotype in allopatric conditions and the dwarf size phenotype in sympatry. The existence of sympatric pairs in the Allegash basin result from the secondary contact of two monophyletic groups of whitefish that evolved allopatrically in separate refugia during the last glaciation events. The weak mtDNA difference of sympatric pairs suggests that speciation of lake whitefish in eastern North America was accompanied by only minor alterations of the ancestral gene pool.

Thermal habit, metabolic rate and the evolution of mitochondrial DNA

The hallmarks of animal mitochondrial DNA (mtDNA) are a rapid rate of sequence evolution, a small genome carrying the same set of homologous genes, maternal inheritance and lack of recombination. Over the past few years, a variety of different observations has challenged these accepted notions of mitochondrial biology. Notable examples include evidence for variable rates of mtDNA sequence evolution among taxa, evidence for large and variable mitochondrial genome sizes in certain groups, and a growing number of cases in metazoans of 'paternal leakage' in the inheritance of mtDNA. Several recent studies have uncovered different lines of evidence suggesting that an organism's thermal habit, or metabolic rate, can influence the evolution of mtDNA.

Basal jawed vertebrate phylogenomics using transcriptomic data from Solexa sequencing

The traditionally accepted relationships among basal jawed vertebrates have been challenged by some molecular phylogenetic analyses based on mitochondrial sequences. Those studies split extant gnathostomes into two monophyletic groups: tetrapods and piscine branch, including Chondrichthyes, Actinopterygii and sarcopterygian fishes. Lungfish and bichir are found in a basal position on the piscine branch. Based on transcriptomes of an armored bichir (Polypterus delhezi) and an African lungfish (Protopterus sp.) we generated, expressed sequences and whole genome sequences available from public databases, we obtained 111 genes to reconstruct the phylogenetic tree of basal jawed vertebrates and estimated their times of divergence. Our phylogenomic study supports the traditional relationship. We found that gnathostomes are divided into Chondrichthyes and the Osteichthyes, both with 100% support values (posterior probabilities and bootstrap values). Chimaeras were found to have a basal position among cartilaginous fishes with a 100% support value. Osteichthyes were divided into Actinopterygii and Sarcopterygii with 100% support value. Lungfish and tetrapods form a monophyletic group with 100% posterior probability. Bichir and two teleost species form a monophyletic group with 100% support value. The previous tree, based on mitochondrial data, was significantly rejected by an approximately unbiased test (AU test, p = 0). The time of divergence between lungfish and tetrapods was estimated to be 391.8 Ma and the divergence of bichir from pufferfish and medaka was estimated to be 330.6 Ma. These estimates closely match the fossil record. In conclusion, our phylogenomic study successfully resolved the relationship of basal jawed vertebrates based on transtriptomes, EST and whole genome sequences.

Monoclonal antibody based immunodot for specific detection of proteins of the shrimp Penaeus species

Frozen shrimp continued to be the single largest item of export from India in terms of value accounting for about 44% of the total marine export earnings. Headless, peeled frozen shrimp is a common and dominant item in the market and there is need for differentiating peeled Penaeus sp from Metapenaeus, Parapenopsis and Macrobrachium sp as consumer preference and price vary. Furthermore, there is need to find out original species used in value addition of shrimp products. Hence, it is essential for development of simple and consumer friendly technique for the identification of shrimp and their products in the market. Two monoclonal antibodies (MAbs) C-15 (IgG3) and C-52 (IgG2a) reacting with 65 and 47 kD proteins of Penaeus monodon respectively in the Western blot were selected. In epitope analysis by immunodot, the two MAbs reacted and recognized specific proteins of P. monodon, Fenneropenaeus indicus and Littopenaeus vannamei and not that of Metapenaeus, Parapenopsis, Macrobrachium rosenbergii, crabs and fishes. The immunodot required 120 min for completion. The sensitivity of the immunodot to detect proteins of P. monodon was 0.225 mg with MAb C-15 and 0.028 mg with MAb C-52. The MAb based immunodot developed, could be used for identifying and differentiating meat of P. monodon, F. indicus, and L. vannamei from that of Metapenaeus, Parapenopsis, M. rosenbergii, crabs and fishes.

Cytochrome-b sequence variation among parrots

The nucleotide sequence of a 307 bp fragment of the mitochondrial cytochrome-b gene was determined for 12 species of parrot, using the polymerase chain reaction and direct sequencing. Sequence divergence ranged from 26-54 differences in pairwise comparisons, with the majority of base substitutions occurring at third positions of codons. The transition:transversion ratio was determined to be higher (approximately 24.3:1) in recently divergent parrot lineages than has generally been observed in other groups. Strongly biased base composition, particularly at the third position of codons, is evident among the sequences. Phylogenetic relationships among more divergent taxa were estimated, using only transversion substitutions, while all the substitutions were useful for closely related taxa. The African genera Psittacus and Poicephalus are closely related, in contrast to the Australian genera Nymphicus, Purpureicephalus and Melopsittacus, which represent more divergent lineages. The cockatoos appear to represent an ancient lineage within the parrots.

Molecular systematics of Salmonidae: combined nuclear data yields a robust phylogeny

The phylogeny of salmonid fishes has been the focus of intensive study for many years, but some of the most important relationships within this group remain unclear. We used 269 Genbank sequences of mitochondrial DNA (from 16 genes) and nuclear DNA (from nine genes) to infer phylogenies for 30 species of salmonids. We used maximum parsimony and maximum likelihood to analyze each gene separately, the mtDNA data combined, the nuclear data combined, and all of the data together. The phylogeny with the best overall resolution and support from bootstrapping and Bayesian analyses was inferred from the combined nuclear DNA data set, for which the different genes reinforced and complemented one another to a considerable degree. Addition of the mitochondrial DNA degraded the phylogenetic signal, apparently as a result of saturation, hybridization, selection, or some combination of these processes. By the nuclear-DNA phylogeny: (1) (Hucho hucho, Brachymystax lenok) form the sister group to (Salmo, Salvelinus, Oncorhynchus, H. perryi); (2) Salmo is the sister-group to (Oncorhynchus, Salvelinus); (3) Salvelinus is the sister-group to Oncorhynchus; and (4) Oncorhynchus masou forms a monophyletic group with O. mykiss and O. clarki, with these three taxa constituting the sister-group to the five other Oncorhynchus species. Species-level relationships within Oncorhynchus and Salvelinus were well supported by bootstrap levels and Bayesian analyses. These findings have important implications for understanding the evolution of behavior, ecology and life-history in Salmonidae.

Cytochrome b nucleotide sequence variation among the Atlantic Alcidae

Analysis of cytochrome b nucleotide sequences of the six extant species of Atlantic alcids and a gull revealed an excess of adenines and cytosines and a deficit of guanines at silent sites on the coding strand. Phylogenetic analyses grouped the sequences of the common (Uria aalge) and Brünnich's (U. lomvia) guillemots, followed by the razorbill (Alca torda) and little auk (Alle alle). The black guillemot (Cepphus grylle) sequence formed a sister taxon, and the puffin (Fratercula arctica) fell outside the other alcids. Phylogenetic comparisons of substitutions indicated that mutabilities of bases did not differ, but that C was much more likely to be incorporated than was G. Imbalances in base composition appear to result from a strand bias in replication errors, which may result from selection on secondary RNA structure and/or the energetics of codon-anticodon interactions.

Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes

Background

Phylogenetic analyses of jawed vertebrates based on mitochondrial sequences often result in confusing inferences which are obviously inconsistent with generally accepted trees. In particular, in a hypothesis by Rasmussen and Arnason based on mitochondrial trees, cartilaginous fishes have a terminal position in a paraphyletic cluster of bony fishes. No previous analysis based on nuclear DNA-coded genes could significantly reject the mitochondrial trees of jawed vertebrates.

Results

We have cloned and sequenced seven nuclear DNA-coded genes from 13 vertebrate species. These sequences, together with sequences available from databases including 13 jawed vertebrates from eight major groups (cartilaginous fishes, bichir, chondrosteans, gar, bowfin, teleost fishes, lungfishes and tetrapods) and an outgroup (a cyclostome and a lancelet), have been subjected to phylogenetic analyses based on the maximum likelihood method.

Conclusion

Cartilaginous fishes have been inferred to be basal to other jawed vertebrates, which is consistent with the generally accepted view. The minimum log-likelihood difference between the maximum likelihood tree and trees not supporting the basal position of cartilaginous fishes is 18.3 ± 13.1. The hypothesis by Rasmussen and Arnason has been significantly rejected with the minimum log-likelihood difference of 123 ± 23.3. Our tree has also shown that living holosteans, comprising bowfin and gar, form a monophyletic group which is the sister group to teleost fishes. This is consistent with a formerly prevalent view of vertebrate classification, although inconsistent with both of the current morphology-based and mitochondrial sequence-based trees. Furthermore, the bichir has been shown to be the basal ray-finned fish. Tetrapods and lungfish have formed a monophyletic cluster in the tree inferred from the concatenated alignment, being consistent with the currently prevalent view. It also remains possible that tetrapods are more closely related to ray-finned fishes than to lungfishes.

Erosion of interspecific reproductive barriers resulting from hatchery supplementation of rainbow trout sympatric with cutthroat trout

The frequency of hybridization between cutthroat (Onchorhynchus clarki clarki) and rainbow (O. mykiss irideus) trout from coastal habitats in British Columbia, Canada, was examined in seven populations where the two species are sympatric with no history of rainbow trout stocking and compared with areas where native rainbow trout populations have been supplemented with hatchery fish (three populations). Four nuclear markers were used to identify each species and interspecific hybrids and one mitochondrial marker showed the direction of gene exchange between species. The frequency of hybrids was significantly higher (Fisher exact test, P < 0.001) in river systems where hatchery rainbow trout have been introduced (50.6% hybrids) than in populations where the two species naturally co-occur without supplementation (9.9% hybrids).

Macroevolutionary relationships of species of Drosophila melanogaster group based on mtDNA sequences

The phylogenetic relationships among the Drosophila melanogaster group species were analyzed using approximately 1700 nucleotide-long sequences of the mitochondrial DNA. Phylogenetic analysis was performed using this region consisting of a part of the cytochrome b (cytb) coding gene, the entire coding sequences of tRNA-Leu, tRNA-Ser and the first subunit of NADH dehydrogenase (NADH1), and a part of the 16S-rRNA gene. The study of these sequences showed that this region of mtDNA is very invariable, as regards with the type of the genes that it contains, as well as the order that they are located on it. The resulting phylogenetic trees reveal a topology that separates the species into three main ancestral lines, leading to the following subgroups: (a) ananassae subgroup, (b) montium subgroup, and (c) melanogaster and Oriental subgroups. The inferred topology complements and generally agrees with previously proposed classifications based on morphological and molecular data.

Sequence polymorphisms within the human mitochondrial genes MTATP6, MTATP8 and MTND4

By sequencing the control region of mitochondrial DNA, the majority of human DNA samples can be differentiated. A further increase in differentiation probability may be possible, e.g. by extending the sequenced region to coding regions of the mitochondrial genome. Restriction to those positions that do not result in a change of the amino acids guarantees that the information thus obtained does not refer to phenotypically relevant information. In the present study the sequence data of the mitochondrial genes MTATP6, MTATP8 and MTND4 were collected from 109 subjects and analyzed in order to define variable positions suitable for identification purposes. There were 32 variable base positions among 850 bases studied from MTATPase genes and 1,200 bases of the MTND4 gene showed 28 variable positions. "Hot spots" for base exchanges were found in both regions and one position (position 11,719 in the MTND4 gene) seems to be suitable for SNP investigation for forensic purposes.

On combining protein sequences and nucleic acid sequences in phylogenetic analysis: the homeobox protein case

Amino acid encoding genes contain character state information that may be useful for phylogenetic analysis on at least two levels. The nucleotide sequence and the translated amino acid sequences have both been employed separately as character states for cladistic studies of various taxa, including studies of the genealogy of genes in multigene families. In essence, amino acid sequences and nucleic acid sequences are two different ways of character coding the information in a gene. Silent positions in the nucleotide sequence (first or third positions in codons that can accrue change without changing the identity of the amino acid that the triplet codes for) may accrue change relatively rapidly and become saturated, losing the pattern of historical divergence. On the other hand, non-silent nucleotide alterations and their accompanying amino acid changes may evolve too slowly to reveal relationships among closely related taxa. In general, the dynamics of sequence change in silent and non-silent positions in protein coding genes result in homoplasy and lack of resolution, respectively. We suggest that the combination of nucleic acid and the translated amino acid coded character states into the same data matrix for phylogenetic analysis addresses some of the problems caused by the rapid change of silent nucleotide positions and overall slow rate of change of non-silent nucleotide positions and slowly changing amino acid positions. One major theoretical problem with this approach is the apparent non-independence of the two sources of characters. However, there are at least three possible outcomes when comparing protein coding nucleic acid sequences with their translated amino acids in a phylogenetic context on a codon by codon basis. First, the two character sets for a codon may be entirely congruent with respect to the information they convey about the relationships of a certain set of taxa. Second, one character set may display no information concerning a phylogenetic hypothesis while the other character set may impact information to a hypothesis. These two possibilities are cases of non-independence, however, we argue that congruence in such cases can be thought of as increasing the weight of the particular phylogenetic hypothesis that is supported by those characters. In the third case, the two sources of character information for a particular codon may be entirely incongruent with respect to phylogenetic hypotheses concerning the taxa examined. In this last case the two character sets are independent in that information from neither can predict the character states of the other. Examples of these possibilities are discussed and the general applicability of combining these two sources of information for protein coding genes is presented using sequences from the homeobox region of 46 homeobox genes from Drosophila melanogaster to develop a hypothesis of genealogical relationship of these genes in this large multigene family.

The mitochondrial genomes of the iguana (Iguana iguana) and the caiman (Caiman crocodylus): implications for amniote phylogeny

The complete mitochondrial genomes of two reptiles, the common iguana (Iguana iguana) and the caiman (Caiman crocodylus), were sequenced in order to investigate phylogenetic questions of tetrapod evolution. The addition of the two species allows analysis of reptilian relationships using data sets other than those including only fast-evolving species. The crocodilian mitochondrial genomes seem to have evolved generally at a higher rate than those of other vertebrates. Phylogenetic analyses of 2889 amino-acid sites from 35 mitochondrial genomes supported the bird-crocodile relationship, lending no support to the Haematotherma hypothesis (with birds and mammals representing sister groups). The analyses corroborated the view that turtles are at the base of the bird-crocodile branch. This position of the turtles makes Diapsida paraphyletic. The origin of the squamates was estimated at 294 million years (Myr) ago and that of the turtles at 278 Myr ago. Phylogenetic analysis of mammalian relationships using the additional outgroups corroborated the Marsupionta hypothesis, which joins the monotremes and the marsupials to the exclusion of the eutherians.

Phylogeography of the tailed frog (Ascaphus truei): implications for the biogeography of the Pacific Northwest

Tailed frogs are distributed in high-gradient streams within the disjunct mesic forests of the Pacific Northwest and represent the basal lineage of the anurans. We sequenced 1,530 nucleotides of the mitochondrial cytochrome b and NADH dehydrogenase subunit two genes from 23 populations and used parsimony, maximum-likelihood, and nested-clade analyses to estimate relationships among populations and infer evolutionary processes. We found two divergent haplotype clades corresponding with inland Rocky Mountain populations and coastal populations and separated by up to 0.133 substitutions per site. Within the coastal assemblage, haplotypes formed clades by mountain range with 0.010-0.024 substitutions per site divergence among populations. Inland haplotypes exhibited minimal genetic structure, with the exception of 0.021 substitutions per site distance between populations from the East Fork of the South Fork of the Salmon River and all other inland haplotypes. The magnitude of divergence between inland and coastal populations, as well as the paleobotanical record, suggest isolation of these lineages occurred during the late Miocene to early Pliocene, probably in response to the rise of the Cascade Mountains. Genetic structure within coastal and inland populations is consistent with isolation in refugia during the late Pliocene and early Pleistocene. Closely related inland haplotypes reflect range expansion following glaciation. The depth of divergence between inland and coastal populations supports the persistence of mesic forests within the inland Pacific Northwest throughout the Pleistocene and is congruent with patterns found in several other mesic forest species. Based on mitochondrial divergence and previous allozyme and morphological data, we recommend recognition of inland populations as a distinct species, Ascaphus montanus.

Molecular phylogeny of osteoglossoids: a new model for Gondwanian origin and plate tectonic transportation of the Asian arowana

One of the traditional enigmas in freshwater zoogeography has been the evolutionary origin of Scleropages formosus inhabiting Southeast Asia (the Asian arowana), which is a species threatened with extinction among the highly freshwater-adapted fishes from the order Osteoglossiformes. Dispersalists have hypothesized that it originated from the recent (the Miocene or later) transmarine dispersal of morphologically quite similar Australasian arowanas across Wallace's Line, but this hypothesis has been questioned due to their remarkable adaptation to freshwater. We determined the complete nucleotide sequences of two mitochondrial protein genes from 12 osteoglossiform species, including all members of the suborder Osteoglossoidei, with which robust molecular phylogeny was constructed and divergence times were estimated. In agreement with previous morphology-based phylogenetic studies, our molecular phylogeny suggested that the osteoglossiforms diverged from a basal position of the teleostean lineage, that heterotidines (the Nile arowana and the pirarucu) form a sister group of osteoglossines (arowanas in South America, Australasia, and Southeast Asia), and that the Asian arowana is more closely related to Australasian arowanas than to South American ones. However, molecular distances between the Asian and Australasian arowanas were much larger than expected from the fact that they are classified within the same genus. By using the molecular clock of bony fishes, tested for its good performance for rather deep divergences and calibrated using some reasonable assumptions, the divergence between the Asian and Australasian arowanas was estimated to date back to the early Cretaceous. Based on the molecular and geological evidence, we propose a new model whereby the Asian arowana vicariantly diverged from the Australasian arowanas in the eastern margin of Gondwanaland and migrated into Eurasia on the Indian subcontinent or smaller continental blocks. This study also implicates the relatively long absence of osteoglossiform fossil records from the Mesozoic.

Association of a novel human mtDNA ATPase6 mutation with immature sperm cells

This study reports the first clearly defined heteroplasmic mutation in immature human sperm cells. The human sperm mitochondrial genome from residue 8186-9341 was analysed with the aim of identifying point mutations which may be associated with human male infertility. The semen samples analysed were obtained from 88 fertile men, 19 with oligozoospermia, and 12 with severe oligozoospermia. Using single strand conformation polymorphism analysis a heteroplasmic T to C transition was detected in the ATPase6 gene, at nucleotide position 8821, in semen samples from one out of 12 (8%) severely oligozoospermic men, but not in oligozoospermic men or normospermic men. This mutation changed the amino acid serine to proline at residue 99 of the mitochondrial ATPase6 in a region which is highly conserved in other vertebrates including rat, bovine, chicken, salmonids and Xenopus. The mutation was detected in semen samples collected from the same man 9 months apart and in peripheral blood lymphocytes. Single sperm cell analyses did not find this mutation in the mature sperm, but the mutation was detected in 7% of immature spermatids. Our finding suggests that immature spermatids with this mutation fail to develop fully.

Sequence polymorphism in the mtDNA HV1 region in Japanese and Chinese

We investigated the nucleotide substitution and insertion/deletion polymorphism of the HV1 region in mtDNA by sequencing blood samples from 150 unrelated Japanese and 120 unrelated Chinese and revealed 108 sequence types from the Japanese group and 87 sequence types from the Chinese. Some substitutions were characteristic of East Asian populations as compared with data reported on Caucasian populations, and some were area-specific among East Asians. The level of genetic diversity and genetic identity revealed by this system was superior to that obtained by VNTR systems for nuclear DNA. These results show the usefulness of mtDNA sequencing in forensic examination for individual identification. We also found some sequence variations in the homopolymeric tract of cytosine (np16180-16194 in the Anderson's reference sequence) that might suggest some hints regarding the mechanisms for and the development of heteroplasmic length variations in this tract.

Zoogeographical implications of variation in mitochondrial DNA of Arctic grayling (Thymallus arcticus)

Most of the northern half of North America's freshwater fish survived the last glacial period in unfrozen refugia peripheral to the ice sheets. In our study, the question of which refugia Arctic grayling (Thymallus arcticus) inhabited during the Wisconsinan Ice Age, and how they subsequently dispersed to their present geographical range, was examined using mitochondrial DNA (mtDNA) analysis. mtDNA from 12 T. arcticus populations was analysed by direct sequencing and restriction fragment length polymorphism analysis (RFLP). Our data support the hypotheses that T. arcticus had a large refugial population in the Bering Refuge (shown by high mtDNA diversity in extant Alaskan populations) and that British Columbia was colonized from the Bering Refuge (shown by mtDNA haplotype similarities). Our data also show that a disjunct southern set of populations in Montana is significantly different from the northern grayling, in terms of restriction haplotype frequency and distinguishing sequence characteristics. Sequence results yielded an estimated divergence time of 370,000 years between the northern and Montana grayling haplotypes. We conclude that T. arcticus survived the Wisconsinan glaciation in at least two refugia: (i) the Bering Refuge north of the ice sheets; and (ii) either the Upper Missouri or the southwest Alberta Refuge, south of the ice sheet.

The complete nucleotide sequence of the mitochondrial DNA of the dogfish, Scyliorhinus canicula

We have determined the complete nucleotide sequence of the mitochondrial DNA (mtDNA) of the dogfish, Scyliorhinus canicula. The 16,697-bp-long mtDNA possesses a gene organization identical to that of the Osteichthyes, but different from that of the sea lamprey Petromyzon marinus. The main features of the mtDNA of osteichthyans were thus established in the common ancestor to chondrichthyans and osteichthyans. The phylogenetic analysis confirms that the Chondrichthyes are the sister group of the Osteichthyes.

The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions

The 17,191-bp mitochondrial DNA (mtDNA) of a Japanese colubrid snake, akamata (Dinodon semicarinatus), was cloned and sequenced. The snake mtDNA has some peculiar features that were found in our previous study using polymerase chain reaction: duplicate control regions that have completely identical sequences over 1 kbp, translocation of tRNALeu(UUR) gene, shortened TpsiC arm for most tRNA genes, and a pseudogene for tRNAPro. Phylogenetic analysis of amino acid sequences of protein genes suggested an unusually high rate of molecular evolution in the snake compared to other vertebrates. Southern hybridization experiments using mtDNAs purified from multiple akamata individuals showed that the duplicate state of the control region is not a transient or unstable feature found in a particular individual, but that it stably occurs in mitochondrial genomes of the species. This may, therefore, be regarded as an unprecedented example of stable functional redundancy in animal mtDNA. However, some of the examined individuals contain a rather scanty proportion of heteroplasmic mtDNAs with an organization of genes distinct from that of the major mtDNA. The gene organization of the minor mtDNA is in agreement with one of models that we present to account for the concerted evolution of duplicate control regions.

The rate of mitochondrial 12S rRNA gene evolution is similar in freshwater turtles and marsupials

Assertions that the "conventional" rate of mitochondrial DNA (mtDNA) evolution is reduced in poikilotherms in general and turtles in particular were tested for side-necked turtles (Pleurodira: Chelidae). Homologous data sets of mitochondrial 12S rRNA gene sequences were used to compare the average divergence between the Australian and South American species for two Gondwanan groups: the chelid turtles and the marsupials. The mean nucleotide divergences between continental groups for both the turtles and the marsupials are remarkably similar. These data suggest that the rate of evolution of mitochondrial 12S rRNA gene is not substantially slower in turtles than in the homeothermic marsupials.

Phylogenetic relationships within genus Leuciscus (Pisces, Cyprinidae) in Portuguese fresh waters, based on mitochondrial DNA cytochrome b sequences

To investigate phylogenetic relationships among Leuciscus species occurring in Portuguese inland waters, the cytochrome b gene was sequenced from representatives of the main rivers. This study supports the recognition of the species level for L. pyrenaicus, including populations from the southern Portuguese drainages (Tejo, Sado, and Guadiana drainages), and for L. carolitertii, including populations from the northern Portuguese drainages. The existence of two new species occurring in the extreme southwestern drainages of Mira and Arade is also suggested. The present results support the monophyly of the Mira and the Arade populations, as well as an early divergence of these two lineages. The present-day distribution of Leuciscus species is seen as a consequence of Pliocene and Pleistocene events, such as river disjunctions and posterior confluence in epicontinental seas and river captures. A mixture of haplotypes was observed in the Mondego and the Tejo drainages, which could be a consequence of ancient river captures, with a possible mitochondrial DNA introgression in the Tejo drainage and a recent introduction by man in the Mondego drainage. The pattern of differentiation among mtDNA haplotypes and their geographic distribution is discussed in terms of evolutionary aspects.

Phylogenetic relationship of the genus Oncorhynchus species inferred from nuclear and mitochondrial markers

The phylogenetic relationship among the salmonid fishes of the genus Oncorhynchus has been analyzed using various kinds of markers for a long time. However, there are three major disagreements among those studies; (1) the authenticity of the Pacific salmon group as a monophyletic cluster, (2) the phylogenetic relationship among three Pacific salmons (pink salmon, sockeye salmon, and chum salmon), and (3) the phylogenetic position of masu salmon. We used allozyme electrophoresis to clarify the phylogenetic relationship between the Pacific salmon group and the Pacific trout group. Furthermore, we reanalysed published mitochondrial DNA D-loop sequences (Shedlock et al., 1992). Allozymic data and mtDNA data indicated the following consistent results; (1) all Pacific salmons formed a monophyletic cluster, (2) chum salmon and pink salmon were clustered within those Pacific salmons, (3) masu salmon formed a cluster with other Pacific salmons and diverged first in this group.

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II B1 and B2 peptides, and of the intervening intron (548-773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for the B1 sequence than for B2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the B1 sequence, the level of nonsynonymous nucleotide substitution (dN) exceeded the level of synonymous substitution (dS) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of dN indicates that balancing selection maintains diversity at this salmonid Mhc class II locus, as is common for Mhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelic B1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100000 years in the Pleistocene.

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (micros Ps)/(microv Pv), where micros and microv denote mutation rate of transitions and transversions, respectively, and Ps and Pv denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions, Ps/Pv approximately 1 at fourfold degenerate sites, so that the s/v ratio at fourfold degenerate sites is approximately micros/microv, which is a measure of mutational contribution to transition bias. Similarly, the s/v ratio at nondegenerate sites is also an estimate of micros/microv if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochrome b (cyt-b) in pocket gophers, the s/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, the s/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the micros/microv ratio is about 2, the Ps/Pv ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content on s/v ratio of protein-coding genes.

Genetic relationship among three subspecies of Oncorhynchus masou determined by mitochondrial DNA sequence analysis

It is generally accepted that there are 3 subspecies of Oncorhynchus masou in Japan, namely, Masu salmon (Oncorhynchus masou masou (Brevoort)), Amago salmon (O. masou ishikawae Jordan & McGregor), and Biwa salmon (O. masou rhodurus Jordan & McGregor or O. masou subsp. Kimura). Since the genetic relationship of these three taxa is not well known, there has been considerable confusion over their nomenclature. We have clarified the genetic relationship among these three taxa by partially sequencing their mitochondrial DNA. Sequences of 948 base pairs from the 3' region of the ATPase subunit 6 gene to the 5' region of the cytochrome oxidase subunit 3 gene were obtained for 20 individuals including wild Biwa salmon, wild and farmed Amago and Masu salmon. Furthermore, 2,162 base pairs from the 3' region of ATPase subunit 6 gene to the 5' region of NADH dehydrogenase subunit 4L gene were determined in 4 individuals. In total, there were 26 sites of base substitutions. The haplotypes of Masu salmon and Amago salmon were similar. On the other hand, 17 of the 26 sites had substitutions characteristic of Biwa salmon. A matrix of genetic distances and maximum parsimony analysis among the haplotypes indicated that Biwa salmon is genetically more distant from Masu and Amago salmon, than Masu salmon is from Amago salmon. This means that Biwa salmon diverged from the common ancestor of the Oncorhynchus masou complex before the divergence between Masu salmon and Amago salmon.

Mitochondrial DNA sequences reveal close relationships between social parasitic ants and their host species

In the tribe Leptothoracini, the phylogenetic relationship of socially parasitic ants (Doronomyrmex kutteri, D. goesswaldi and Harpagoxenus sublaevis) and their host species Leptothorax acervorum has been controversial. Even more controversial is the relationship between the socially parasitic ant Chalepoxenus muellerianus and its host species Leptothorax unifasciatus, L. nigriceps, L. interruptus and L. recedens. On the basis of morphological, ecological and ethological criteria it has been argued that socially parasitic ants and their respective hosts always evolved from common ancestors, and hence it has been postulated that these species should be included in common taxonomical groups. This would require the division of the tribe Leptothoracini into two subgroups, one comprising the subgenus Leptothorax (s. str.) and the other the subgenus Myrafant, together with their respective parasitic genera. We have used the polymerase chain reaction (PCR) to compare a 360-bp sequence of the mitochondrial cytochrome b gene of 14 species belonging to the tribe Leptothoracini and an outgroup species Tetramorium impurum (Tetramoriini). The results generally agree with the morphological studies which suggest that a common ancestral species differentiated into host and parasite species. This relationship is very obvious within the Leptothorax (s. str.) group but less pronounced in the species belonging to the Myrafant group. Leptothorax (Temnothorax) recedens shows a greater sequence divergence than the outgroup species T. impurum.

Cytochrome b sequence variation and a molecular phylogeny of the live-bearing fish genus Gambusia (Cyprinodontiformes: Poeciliidae)

Nucleotide sequences for a 402-base segment of the mitochondrial cytochrome b gene were determined from 25 species of live-bearing fishes. A total of 34 sequences representing 24 species of the genus Gambusia and 1 species of Belonesox were generated via the polymerase chain reaction. The levels of overall variation were consistent with those from other genera of fishes. In total, 137 of 402 (34.1%) nucleotides exhibited variation within or among the species. Observed differences at 24 (17.9%) of the 134 cytochrome b codons would result in amino acid replacements. Phylogenetic analyses employing various weighting schemes resulted in several clades representing traditionally recognized taxonomic groups. However, precise relationships among species-groups remained uncertain. Randomization tests indicated that these topologies contained significant nonrandom phylogenetic information. As with other fishes, the overall rate of divergence appeared to be slower than that of other vertebrates and the overall replacement/substitution pattern was suggestive of nonrandom evolutionary input.

Evolutionary analysis of cytochrome b sequences in some Perciformes: evidence for a slower rate of evolution than in mammals

To obtain information relative to the phylogenesis and microevolutionary rate of fish mitochondrial DNA, the nucleotide sequence of cytochrome b gene in seven fish species belonging to the order of Perciformes was determined. Sequence analysis showed that fish mitochondrial DNA has a nucleotide compositional bias similar to that of sharks but lower compared to mammals and birds. Quantitative evolutionary analysis, carried out by using a markovian stochastic model, clarifies some phylogenetic relationships within the Perciformes order, particularly in the Scombridae family, and between Perciformes, Gadiformes, Cypriniformes, and Acipenseriformes. The molecular clock of mitochondrial DNA was calibrated with the nucleotide substitution rate of cytochrome b gene in five shark species having divergence times inferred from paleontological estimates. The results of such analysis showed that Acipenseriformes diverged from Perciformes by about 200 MY, that the Perciformes common ancestor dates back to 150 MY, and that fish mitochondrial DNA has a nucleotide substitution rate three to five times lower than that of mammals.

Phylogeny of the Drosophila obscura species group deduced from mitochondrial DNA sequences

Approximately 2 kb corresponding to different regions of the mtDNA of 14 different species of the obscura group of Drosophila have been sequenced. In spite of the uncertainties arising in the phylogenetic reconstruction due to a restrictive selection toward a high mtDNA A+T content, all the phylogenetic analysis carried out clearly indicate that the obscura group is formed by, at least, four well-defined lineages that would have appeared as the consequence of a rapid phyletic radiation. Two of the lineages correspond to monophyletic subgroups (i.e., affinis and pseudoobscura), whereas the obscura subgroup remains heterogeneous assemblage that could be reasonably subdivided into at least two complexes (i.e., subobscura and obscura).

Cloning and characterization of the European seabass, Dicentrarchus labrax, mitochondrial genome

Mitochondrial DNA (mtDNA) from the European seabass, Dicentrarchus labrax, has been cloned and characterized. Its gene organization was deduced by a comparison of the sequenced termini of different subclones obtained from European seabass mtDNA to the completely-sequenced mtDNAs from carp and freshwater loach. The difference in genome size between the European seabass mtDNA (approximately 18 kb) and most of the other characterized fish mtDNAs (approximately 16.5 kb) is accounted for by the displacement-loop (D-loop). Comparisons have been performed between the derived amino-acid sequences of three sequenced genes, cytochrome c oxidase subunit 2 (COII), NADH dehydrogenase subunit 4L (ND4L) and ATP synthase subunit 8 (ATPase8), from D. labrax, and their counterparts in other fishes and Xenopus laevis.

Variation in mitochondrial DNA and population structure of the Taipei treefrog Rhacophorus taipeianus in Taiwan

Mitochondrial DNA (mtDNA) from 40 samples of the Taipei treefrog Rhacophorus taipeianus collected from seven populations in Taiwan were sequenced to document the DNA sequence variation in anuran mtDNA and to elucidate the phylogeographic population structure in the Taipei treefrog. Sequences of 722-764 bases in length, including a 108-bp segment of the cytochrome b gene and a 614-656-bp D-loop segment, were obtained by direct sequencing using polymerase chain reaction (PCR). The variation in length was due to a 40-bp region that tandemly repeated four to five times in the D-loop region. The first repeat is the most conserved one among the five repeats because there are no variable sites in this repeat. Besides the 40-bp length variation, 28 positions in the 764-bp sequences are variable and distributed evenly in the cytochrome b gene fragment and D-loop region. Variation in the D-loop of the Taipei treefrog is comparable to those of other vertebrates. Two well-differentiated lineages (northern and central) differing by mean sequence divergence of 1.7% are identified and concordant with their geographic distributions. The two lineages are inferred to have split from a common ancestral population in the early Pleistocene. However, the interpopulation divergence of the northern lineage (< 0.33%) is apparently lower than that of the central lineage (1.11%), implying that the two lineages evolved independently and had different demographic histories after divergence. This study reveals that anuran D-loop has potential as a genetic marker in phylogenetic and population genetic analyses of anurans.