Genetic variation and population structure of marbled sole Pleuronectes yokohamae and cresthead flounder P. schrenki in Japan inferred from mitochondrial DNA analysis

Nucleotide sequence analysis in the first half of the mitochondrial control region was used to estimate genetic variation and population structure of marbled sole Pleuronectes yokohamae and cresthead flounder Pleuronectes schrenki, obtained from 10 localities in Japan. Overall estimate of pairwise population F_ST and NJ tree suggested genetic distinction between P. yokohamae and P. schrenki. SAMOVA also demonstrated two distinct genetic groups of the samples, MS (marbled sole, Honshu and southernmost Hokkaido) and CF(cresthead flounder, Hokkaido), which was pertinent to the natural range of P. yokohamae and P. schrenki, respectively. However, the haplotypes common to P. yokohamae and P. schrenki were found in all four samples of the CF group and two of six samples of the MS group. Haplotype and nucleotide diversities were high within MS and CF groups, and the mismatch distribution analysis showed a bimodal profile in both groups, suggesting secondary contact between allopatrically evolved MS and CF groups. Thus, the observed haplotype sharing between two species might have been the consequence of interspecific hybridization after their divergence. In addition, F_ST estimates and Mantel test suggested that the population structure of P. yokohamae might have been influenced by low or restricted gene flow with isolation by distance, whereas no such population structuring was implied for P. schrenki. The observed difference in population structuring between two species was surmised by the difference of migration character.     

Molecular characteristics of Tibetan antelope (Pantholops hodgsonii) mitochondrial DNA control region and phylogenetic inferences with related species

Although Tibetan antelope (Pantholops hodgsonii) is a distinctive wild species inhabiting the Tibet-Qinghai Plateau, its taxonomic classification within the Bovidae is still unclear and little molecular information has been reported to date. In this study of Tibetan antelope, the complete control regions of mtDNA were sequenced and compared to those of Tibetan sheep (Ovis aries) and goat (Capra hircus). The length of the control region in Tibetan antelope, sheep and goat is 1067, 1181/1106 and 1121 bp, respectively. A 75-bp repeat sequence was found near the 5′ end of the control region of Tibetan antelope and sheep, the repeat numbers of which were two in Tibetan antelope and three or four in sheep. Three major domain regions, including HVI, HVII and central domain, in Tibetan antelope, sheep and goat were outlined, as well as other less conserved blocks, such as CSB-1, CSB-2, ETAS-1 and ETAS-2. NJ cluster analysis of the three species revealed that Tibetan antelope was more closely related to Tibetan sheep than Tibetan goat. These results were further confirmed by phylogenetic analysis using the partial control region sequences of these and 13 other antelope species. Tibetan antelope is better assigned to the Caprinae rather than the Antilopinae subfamily of the Bovidae.     

Structure of mitochondrial DNA control region of Fenneropenaeus chinensis and phylogenetic relationship among different populations

This paper deals with the structure of mitochondrial DNA control region of Fenneropenaeus chinensis. The termination-associated sequence (TAS), cTAS, CSB-D-CSB-F, and CSB-1 are detected in the species. The results indicate that the structures of these parts are similar to those of most marine organisms. Two conserved regions and many stable conserved boxes are found in the extended TAS area, central sequences blocks, and conserved sequences blocks (CSBs). This is the special character of F. chinensis. All the mtDNA control region sequences do not have CSB2 and CSB3 blocks, which is quite different from most vertebrates. In addition, the complete mtDNA control region sequences are used to analyze the phylogenetic relationships of F. chinensis. The phylogenetic trees show a lack of genetic structure among populations, which is similar to many previous studies.     

Complete sequence of the mitochondrial genome of Odontamblyopus rubicundus (Perciformes: Gobiidae): genome characterization and phylogenetic analysis

Odontamblyopus rubicundus is a species of gobiid fishes, inhabits muddy-bottomed coastal waters. In this paper, the first complete mitochondrial genome sequence of O. rubicundus is reported. The complete mitochondrial genome sequence is 17119 bp in length and contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a control region and an L-strand origin as in other teleosts. Most mitochondrial genes are encoded on H-strand except for ND6 and seven tRNA genes. Some overlaps occur in protein-coding genes and tRNAs ranging from 1 to 7 bp. The possibly nonfunctional L-strand origin folded into a typical stem-loop secondary structure and a conserved motif (5^′-GCCGG-3^′) was found at the base of the stem within the tRNA ^Cys gene. The TAS, CSB-2 and CSB-3 could be detected in the control region. However, in contrast to most of other fishes, the central conserved sequence block domain and the CSB-1 could not be recognized in O. rubicundus, which is consistent with Acanthogobius hasta (Gobiidae). In addition, phylogenetic analyses based on different sequences of species of Gobiidae and different methods showed that the classification of O. rubicundus into Odontamblyopus due to morphology is debatable.     

Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.     

Assessment of the State of Flounders from Sivuch'ya Bight (Peter the Great Bay, Sea of Japan) Based on Histopathological Parameters

The liver and gills of five flounder species (Pleuronectes obscurus, P. punctatissimus, P. yokahamae, P. herzensteini, and Platichthys stellatus) from Sivuch'ya Bight were histologically studied. Specimens of all species were characterized by a wide spectrum of histomorphological changes. The proportions of most specific pathologies in P. obscurusfrom Sivuch'ya Bight were estimated. The pathologies comprised two major groups: those of the gills (edema and fusion of the respiratory lamellae; lifting and hyperplasia of the respiratory epithelium; vacuolization of the epithelium cells; and the presence of clavate lamellae, aneurysmal areas, and parasites) and those of the liver (vacuolization of the hepatocytes, aneurysms, necrosis of the parenchyma tissue, pathologies of the hepatocyte nuclei, pigment accumulation, and an increase in the number and size of melanomacrophage centers). The estimated proportions of pathological changes were compared with those observed in P. obscurusspecimens from the most contaminated area of Amurskii Bay. In both Sivuch'ya Bight and Amurskii Bay, the proportion of fish with various pathologies was high. The comparison of pathological changes in the gills and liver of P. obscurusand P. stellatusfrom Sivuch'ya Bight showed that, in terms of many parameters, the former species is more sensitive to pollution of the marine environment than the latter. Pathomorphological changes in the organs of the flounders studied are apparently caused by chronic pollutant action and suggest that the fish of both areas are in a damaged state.     

The mitochondrial DNA control region comparison studies of four hinged turtles and its phylogentic significance of the genusCuora sensu lato (Testudinata: Geoemydidae)

The complete sequences of the mitochondrial DNA (mtDNA) control region (CR) ofCistoclemmys flavomarginata, Cistoclemmys galbinifrons, Cuora aurocapitata andCyclemys atripons were amplified by long-polymerase chain reaction (Long-PCR). The lengths were 1207 bp, 1722 bp, 1379 bp and 980 bp, respectively. Combining with the CR sequence ofPyxidea mouhotii (DQ659152), we compared the CR structure, and identified three functional domains (TAS, CD and CSB) in which the conservation sequences (TAS, CSB-F, CSB-1, CSB-2 and CSB-3) were also successfully identified according to their homology to those of other turtles. These 5 turtles have the identical CSB-2 and CSB-3 sequences, and 4 of them have the same CSB-1 sequence while there is one base transversion (T → A) inCy. atripons. We analyzed the variable number of tandem repeat (VNTR) sequences or microsatellites at the 3′ end of CR. The motifs of tandem repeats (7 types) are made up of 2–8 nucleotides, and the copy numbers are from 4 to 48. All of the 5 turtles exceptCy. atripons have the “TATTATAT” repeats and are ended by TA. The results of CR structure analysis displayed that theCuora, Cistoclemmys, andPyxidea have many similarities, but differ fromCyclemys. WithIndotestudo elongate (DQ080043) andIndotestudo forstenii (DQ080044) as outgroups, using the CR sequences (1123bp) excluded the microsatellites at the 3′ end of CR, we constructed the molecular phylogenetic trees using the MP, ML and BI methods. The results showed that there was a strong support to the monophyly of theCuora group consisting ofCuora,Cistoclemmys andPyxidea, which has a close relationship withMauremys andChinemys but far fromCyclemys, which are consistent with the analysis of the CR structure of the 5 turtles.     

Comparison research and phylogenetic implications of mitochondrial control regions in four soft-shelled turtles of Trionychia (Reptilia, Testudinata)

The mitochondrial control regions (CRs) and flanking sequences of Pelodiscus sinensis, Apalone ferox, Palea steindachneri and Carettochelys insculpta were obtained using Long-PCR with gene-specific primers. The CR lengths of the four species were 1843 bp, 1356 bp, 1725 bp, and 969 bp. The base composition percentages of A+T were 60.5%, 60.7%, 65.7%, 64.7%, respectively. Combined with CR sequences of other three soft-shelled turtles published in GenBank (Pelodiscus sinensis, Korea, AY962573; Dogania subplana, AF366350; Lissemys punctata, EF050073), we compared the CR structures and identified three functional domains (TAS, CD and CSB) in which conserved sequence blocks (TAS, CSB -F, CSB-1, CSB-2 and CSB-3) were also successfully identified according to their sequence similarities to those of other turtles. The variable numbers of tandem repeats (VNTRs 1) with 50–52 bp motif were identified at 5′-end of CR among the five soft-shelled turtles P. sinensis (China), P. sinensis (Korea), A. ferox, P. steindachneri, D. subplana. The copy number of the VNTRs varied from 5 to 15. VNTRs 2 with 2–11 bp motif were identified in the 3′- end of CR among all of the six soft-shelled turtles with variable number of motifs from 4 to 29. Moreover, VNTRs 3 with 6 bp motif were identified between CSB-1 and CSB-2 of CR both in P. sinensis (China) and P. sinensis (Korea), in which the number of motifs varied from 19 to 29. The types and distribution of VNTRs of the six soft-shelled turtles were also discussed. With Alligator mississippiensis as an outgroup, combined with the CR sequences (excluding VNTRs) of other five turtles which were published in GenBank, the molecular phylogenetic trees were constructed using PAUP 4.0b10 and MrBayes ver. 3.0. The results strongly supported the monophyly of Carretochelyidae and Carettochelyidae as sister group to an assemblage of Cryptodira. Our research suggested that the earliest phylogenetic tree splits into three separated basal branches; the Pelomedusidira (Pelomedusa subrufa), the Carettochelyidae (C. insculpta), and an assemblage of Cryptodira and the C. insculpta that might be a representation of distinctive suborder.     

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.     

Population genetic structure of the point-head flounder, <Emphasis Type="Italic">Cleisthenes herzensteini</Emphasis>, in the Northwestern Pacific

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. To determine the population genetic structure and the demographic and colonization history of Cleisthenes herzensteini in the Northwestern Pacific, one hundred and twenty-one individuals were sampled from six localities along the coastal regions of Japan and the Yellow Sea of China. Mitochondrial DNA variation was analyzed using DNA sequence data from the 5′ end of control region. High levels of haplotype diversity (>0.96) were found for all populations, indicating a high level of genetic diversity. No pattern of isolation by distance was detected among the population differentiation throughout the examined range. Analyses of molecular variance (AMOVA) and the conventional population statistic Fst revealed no significant population genetic structure among populations. According to the exact test of differentiation among populations, the null hypothesis that C. herzensteini within the examined range constituted a non-differential mtDNA gene pool was accepted. The demographic history of C. herzensteini was examined using neutrality test and mismatch distribution analyses and results indicated Pleistocene population expansion (about 94–376 kya) in the species, which was consistent with the inference result of nested clade phylogeographical analysis (NCPA) showing contiguous range expansion for C. herzensteini. The lack of phylogeographical structure for the species may reflect a recent range expansion after the glacial maximum and insufficient time to attain migration-drift equilibrium.     

Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences

The molecular evidence of phylogenetic status regarding the Formosan serow (Capricornis swinhoei) is not robust and little is known about the genetic diversity of the Sumatran serow (Capricornis sumatraensis), which partly is due to the hardness in sample collection. Here we determined the sequences of the complete mitochondrial DNA control region (1,014 bp) of 19 Sumatran-serow individuals. Nine new haplotypes were defined based on 78 variable sites. Combined analysis with other 32 haplotypes downloaded from the public database, including 1 Sumatran-serow, 11 Formosan-serow and 20 Japanese-serow (Capricornis crispus) haplotypes, a relatively high level of nucleotide diversity was first observed in Sumatran serow (π = 0.0249). By comparative analysis with structural consensus sequences from other mammals, we have identified central, left and right domains and depicted the putative functional structure, including extend termination associated sequences and conserve sequence blocks, in mtDNA control region. The alignment of mtDNA control region revealed that both Sumatran and Japanese serow have two tandem repeats (TRs), but three TRs in Formosan serow. Phylogenetic analyses revealed that the Formosan serow is distinct species with the Japanese serow, but a sister group with the Sumatran serow. The divergence time estimated among three serow species revealed that Pleistocene climate changes and the uplift of Qinghai-Tibetan plateau might play an important role in the genetic differentiation of the serows. These results mainly provide the convinced evidence on the genetic relationship between three serow species.     

Recombination and Evolution of Duplicate Control Regions in the Mitochondrial Genome of the Asian Big-Headed Turtle,Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3′ end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.     

Assessing Nuclear and Mitochondrial DNA Sequence Variation Within Steinernema (Rhabditida: Steinernematidae)

DNA sequence analysis was used to characterize the nuclear ribosomal DNA ITS1 region and a portion of the COII and 16S rDNA genes of the mitochondrial genome from Steinernema entomopathogenic nematodes. Nuclear ITS1 nucleotide divergence among seven Steinernema spp. ranged from 6 to 22%, and mtDNA divergence among five species ranged from 12 to 20%. No intraspecific variation was observed among three S. feltiae strains. Phylogenetic analysis of both nuclear and mitochondrial DNA sequences confirms the existing morphological relationships of several Steinernema species. Both the rDNA ITS1 and mtDNA sequences were useful for resolving relationships among Steinernema taxa.     

Comparative Mitogenomics of the Assassin Bug Genus Peirates (Hemiptera: Reduviidae: Peiratinae) Reveal Conserved Mitochondrial Genome Organization of P. atromaculatus,P. fulvescens and P. turpis

In this study, we sequenced four new mitochondrial genomes and presented comparative mitogenomic analyses of five species in the genus Peirates (Hemiptera: Reduviidae). Mitochondrial genomes of these five assassin bugs had a typical set of 37 genes and retained the ancestral gene arrangement of insects. The A+T content, AT- and GC-skews were similar to the common base composition biases of insect mtDNA. Genomic size ranges from 15,702 bp to 16,314 bp and most of the size variation was due to length and copy number of the repeat unit in the putative control region. All of the control region sequences included large tandem repeats present in two or more copies. Our result revealed similarity in mitochondrial genomes of P. atromaculatus, P. fulvescens and P. turpis, as well as the highly conserved genomic-level characteristics of these three species, e.g., the same start and stop codons of protein-coding genes, conserved secondary structure of tRNAs, identical location and length of non-coding and overlapping regions, and conservation of structural elements and tandem repeat unit in control region. Phylogenetic analyses also supported a close relationship between P. atromaculatus, P. fulvescens and P. turpis, which might be recently diverged species. The present study indicates that mitochondrial genome has important implications on phylogenetics, population genetics and speciation in the genus Peirates.     

Organization of the mitochondrial genome of Atlantic cod, Gadus morhua.

The mitochondrial DNA (mtDNA) from the Atlantic cod, Gadus morhua, was mapped using 11 different restriction enzymes and cloned into plasmid vectors. Sequence data obtained from more than 10 kilobases of cod mtDNA show that the genome organization, genetic code, and the overall codon usage have been conserved throughout the evolution of vertebrates. Comparison of the derived amino acid sequences of proteins encoded by cod mtDNA to the ones encoded by Xenopus laevis mtDNA revealed that the amino acid identity range from 46% to 93% for the different proteins. ND4L is most divergent while COI is most conserved. GUG was found as the translation initiation codon of the COI gene, indicating a dual coding function for this codon. The sequences of the 997 base pair displacement-loop (D-loop)-containing region and the origin of L-strand replication (oriL), are presented. Only few of the primary and secondary structure features found to be conserved among mammalian mitochondrial D-loops, can be identified in cod. Presence of CSB-2 in the D-loop-containing region and the conserved hairpin structure at oriL, indicates that replication of bony fish mtDNA may follow the same general scheme as described for higher vertebrates.     

Mitochondrial genome of the eurasian otterLutra lutra (Mammalia, Carnivora, Mustelidae)

We determined the complete mitochondrial genome of the Eurasian otterLutra lutra, which is an endangered species in Korea. The circle genome (16,536 bp in size) consists of 13 protein-coding, 22 tRNA, and 2 rRNA genes, and a control region, as found in other metazoan animals. Out of the 37 genes, 28 are encoded on the H-strand, and the nine (ND6 and 8 tRNA genes) on the L-strand. Three overlaps among the 13 protein-coding genes were found: ATP8-ATP6, ND4L-ND4, and ND5-ND6. A control region (1090 bp) including the origin of H-strand replication (OH), TAS (a conserved motif TACAT-16bp-ATGTA) and CSB (CSB-1, CSB-2. and CSB-3) was observed between tRNA-Pro and tRNA-Phe genes, and O_L, with 36 highly conserved nucleotides between tRNA-Asn (N) and tRNA-Cys (C) within a cluster of five tRNA genes (WANCY), as typically found in vertebrates. The other important characteristics of theL. lutra mitochondrial genome were described in detail. In addition, a maximum likelihood and Bayesian trees of 9 mustelid species and 1 outgroup were reconstructed based on the nucleotide sequences of 11 protein-coding genes excluding ATP8 and ND6. It showed that Lutrinae formed a monophyletic group with Mustelinae that is not monophyletic. Within the subfamily Lutrinae,L. lutra andEnhydra lutris were grouped together and thenLontra canadentis placed as a sister of the clade. The present result is the first complete mitochondrial genome sequence reported from the genusLutra, and is applicable to molecular phylogenetic, phylogeographic, conservation biological studies for mustelid members. In particular, exploration of sequence variations of the control region may be helpful for analyzing inter-and intra-species variations in the genusLutra.     

Complete mitochondrial genomes of three sea basses Lateolabrax (Perciformes,Lateolabracidae) species: Genome description and phylogenetic considerations

Genus Lateolabrax consists of three species, Japanese sea bass Lateolabrax japonicus, spotted sea bass Lateolabrax maculatus and blackfin sea bass Lateolabrax latus. The complete mitochondrial DNA (mtDNA) of the three sea basses were amplified and sequenced to characterize and discuss their phylogenetic relationships. The length of mitogenomes was 16,593 bp, 16,479 bp and 16,600 bp, respectively, and all of them consisted of 13 protein-coding genes, 2 ribosomal RNA (rRNA), 22 transfer RNA (tRNA) and a control region, which are typical for mtDNA of vertebrate. Most genes were encoded on the H-strand, except for the ND6 and eight tRNA genes encoding on the L-strand. A significant variation among the three species was detected in length of the control region. Phylogenetic relationship among the three species was constructed based on the datasets, including the 12 protein-coding genes (except ND6 gene), 22 tRNA and 2 rRNA sequences. The results supported the sister taxon between L. japonicus and L. maculatus. The genetic resources reported here are useful for further studies in taxonomy and phylogeny of the three sea basses and related species.