Tandemly repeated sequences in the mitochondrial DNA control region and phylogeography of the Pike-Perches Stizostedion

Comparative mitochondrial genome brings insights to slight variation in gene proportion and large intergenic spacer and phylogenetic relationship of mudskipper species

Fish mitochondrial genome have been largely studied worldwide for evolutionary and other genetic purposes and the structure and gene organization are commonly conservative. However, several studies have demonstrated that this scenario may present variations in some taxa, showing differentiation on the gene rearrangement. In this study, the complete mitogenome of terrestrial fish Boleophthalmus dussumieri was generated and compared with other species of the Exudercidae fishes. The newly complete mitogenome generated is circular and 16,685 bp of length, and it contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA genes (tRNAs), and one control region (CR), with high conservative structure, like other Mudskippers. Most of the PCG showed similar codon usage bias. The gene length was found to be different specially for the CR, 12S rRNA gene and ND5 gene in some taxon. All the Boleophthalmus species showed a gene duplication in the CR, except for B. dussumieri, and they presented a long intergenic spacer specially on the tRNA-Pro/ OH Tandem duplication/random loss (TDRL) and dimer-mitogenome and nonrandom loss (DMNL) are suitable to explain the mitogenome rearrangement observed in this study. The phylogenetic analysis well supported the monophyly of all mudskipper species and the analysis positioned the Periophthalmus clade as the most basal of the terrestrial fishes. This finding provides basis and brings insights for gene variation, gene rearrangements and replications showing evidence for variety of mitochondrial structure diversity within mudskippers.

Variants in the mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) infected with Oryctes rhinoceros nudivirus in oil palm and coconut plantations

The CRB (coconut rhinoceros beetle) haplotype was classified into CRB-S and CRB-G, based on the presence of single nucleotide polymorphisms (SNPs) in the mitochondrial cox1 gene. Mitochondrial genomes (mitogenomes) are the most widely used genetic resources for molecular evolution, phylogenetics, and population genetics in relation to insects. This study presents the mitogenome CRB-G and CRB-S which were collected in Johor, Malaysia. The mitogenome of CRB-G collected from oil palm plantations in 2020 and 2021, and wild coconut palms in 2021 was 15,315 bp, 15,475 bp, and 17,275 bp, respectively. The CRB-S was discovered in coconut and oil palms in 2021, and its mitogenome was 15,484 bp and 17,142 bp, respectively. All the mitogenomes have 37 genes with more than 99% nucleotide sequence homology, except the CRB-G haplotype collected from oil palm in 2021 with 89.24% nucleotide sequence homology. The mitogenome of Johor CRBs was variable in the natural population due to its elevated mutation rate. Substitutions and indels in cox1, cox2, nad2 and atp6 genes were able to distinguish the Johor CRBs into two haplotypes. The mitogenome data generated in the present study may provide baseline information to study the infection and relationship between the two haplotypes of Johor CRB and OrNV in the field. This study is the first report on the mitogenomes of mixed haplotypes of CRB in the field.

Chromosome-Level Genome Assembly of the Bioluminescent Cardinalfish Siphamia tubifer: An Emerging Model for Symbiosis Research

The bioluminescent symbiosis involving the sea urchin cardinalfish Siphamia tubifer and the luminous bacterium Photobacterium mandapamensis is an emerging vertebrate model for the study of microbial symbiosis. However, little genetic data are available for the host, limiting the scope of research that can be implemented with this association. We present a chromosome-level genome assembly for S. tubifer using a combination of PacBio HiFi sequencing and Hi-C technologies. The final assembly was 1.2 Gb distributed on 23 chromosomes and contained 32,365 protein coding genes with a BUSCO score of 99%. A comparison of the S. tubifer genome to that of another nonluminous species of cardinalfish revealed a high degree of synteny, whereas a comparison to a more distant relative in the sister order Gobiiformes revealed the fusion of two chromosomes in the cardinalfish genomes. The complete mitogenome of S. tubifer was also assembled, and an inversion in the vertebrate WANCY tRNA genes as well as heteroplasmy in the length of the control region were discovered. A phylogenetic analysis based on whole the mitochondrial genome indicated that S. tubifer is divergent from the rest of the cardinalfish family, highlighting the potential role of the bioluminescent symbiosis in the initial divergence of Siphamia. This high-quality reference genome will provide novel opportunities for the bioluminescent S. tubifer-P. mandapamensis association to be used as a model for symbiosis research.

Relative species abundance estimation in artificial mixtures of insects using mito-metagenomics and a correction factor for the mitochondrial DNA copy number

Mito-metagenomics (MMG) is becoming an alternative to amplicon metabarcoding for the assessment of biodiversity in complex biological samples using high-throughput sequencing. Whereas MMG overcomes the biases introduced by the PCR step in the generation of amplicons, it is not yet a technique free of shortcomings. First, as the reads are obtained from shotgun sequencing, a very low proportion of reads map into the mitogenomes, so a high sequencing effort is needed. Second, as the number of mitogenomes per cell can vary among species, the relative species abundance (RSA) in a mixture could be wrongly estimated. Here, we challenge the MMG method to estimate the RSA using artificial libraries of 17 insect species whose complete genomes are available on public repositories. With fresh specimens of these species, we created single-species libraries to calibrate the bioinformatic pipeline and mixed-species libraries to estimate the RSA. Our results showed that the MMG approach confidently recovers the species list of the mixtures, even when they contain congeneric species. The method was also able to estimate the abundance of a species across different samples (within-species estimation) but failed to estimate the RSA within a single sample (across-species estimation) unless a correction factor accounting for the variable number of mitogenomes per cell was used. To estimate this correction factor, we used the proportion of reads mapping into mitogenomes in the single-species libraries and the lengths of the whole genomes and mitogenomes.

Novel gene rearrangement in the mitochondrial genome of Siliqua minima (Bivalvia, Adapedonta) and phylogenetic implications for Imparidentia

Siliqua minima (Gmelin, 1791) is an important economic shellfish species belonging to the family Pharidae. To date, the complete mitochondrial genome of only one species in this family (Sinonovacula constricta) has been sequenced. Research on the Pharidae family is very limited; to improve the evolution of this bivalve family, we sequenced the complete mitochondrial genome of S. minima by next-generation sequencing. The genome is 17,064 bp in length, consisting of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), and two ribosomal RNA genes (rRNA). From the rearrangement analysis of bivalves, we found that the gene sequences of bivalves greatly variable among species, and with closer genetic relationship, the more consistent of the gene arrangement is higher among the species. Moreover, according to the gene arrangement of seven species from Adapedonta, we found that gene rearrangement among families is particularly obvious, while the gene order within families is relatively conservative. The phylogenetic analysis between species of the superorder Imparidentia using 12 conserved PCGs. The S. minima mitogenome was provided and will improve the phylogenetic resolution of Pharidae species.

Population genetics and comparative mitogenomic analyses reveal cryptic diversity of Amphioctopus neglectus (Cephalopoda: Octopodidae)

This study presented 96 cox1 and 76 cox3 genes of Amphioctopus neglectus populations. Three distinct lineages were formed from phylogenetic trees and networks constructed using haplotypes. Mitogenomes of A. neglectus-a and A. neglectus-b as the representatives of two lineages separated from population genetics were sequenced to compare with A. neglectus at the genome-level. Amphioctopus neglectus-a showed significant differences with A. neglectus, mainly reflected in gene length, intergenic regions and the secondary structure of tandem repeat motifs. Notably, two sequence deletions in mitogenomes of the two representative species were detected in different positions of major non-coding regions, which were the most distinct differences with A. neglectus. Pairwise genetic distances and the phylogenetic analysis supported the relationship of (A. neglectus-a + (A. neglectus + A. neglectus-b)). This study suggested that A. neglectus-a should be considered as a potential cryptic species of this complex, while A. neglectus-b needed further verification to be defined.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

Mitochondrial genome annotation and phylogenetic placement of Oreochromis andersonii and O. macrochir among the cichlids of southern Africa

Genetic characterization of southern African cichlids has not received much attention. Here, we describe the mitogenome sequences and phylogenetic positioning of Oreochromis andersonii and O. macrochir among the African cichlids. The complete mitochondrial DNA sequences were determined for O. andersonii and O. macrochir, two important aquaculture and fisheries species endemic to southern Africa. The complete mitogenome sequence lengths were 16642 bp and 16644 bp for O. andersonii and O. macrochir respectively. The general structural organization follows that of other teleost species with 13 protein-coding genes, 2 rRNAs, 22 tRNAs and a non-coding control region. Phylogenetic placement of the two species among other African cichlids was performed using Maximum Likelihood (ML) and Bayesian Markov-Chain-Monte-Carlo (MCMC). The consensus trees confirmed the relative positions of the two cichlid species with O. andersonii being very closely related to O. mossambicus and O. macrochir showing a close relation to both species. Among the 13 mitochondrial DNA protein coding genes ND6 may have evolved more rapidly and COIII was the most conserved. There are signs that ND6 may have been subjected to positive selection in order for these cichlid lineages to diversify and adapt to new environments. More work is needed to characterize the southern Africa cichlids as they are important species for capture fisheries, aquaculture development and understanding biogeographic history of African cichlids. Bio-conservation of some endangered cichlids is also essential due to the threat by invasive species.

Structure, evolution and phylogenetic informativeness of eelpouts (Cottoidei: Zoarcales) mitochondrial control region sequences

Control region (CR) is a major non-coding domain of mitochondrial DNA in vertebrates which contains the promoters for replication and transcription of mitochondrial genome along with the binding sites for metabolic machinery and, hence, is a vital element for the integrity of mitochondrial genome as a biological replicator. The origin and diversity of structural elements within CR have been intensively studied in recent years with the involvement of new diverse taxa. In this paper, we provide new data on the nucleotide and structural patterns of CR evolution and phylogenetic suitability among eelpouts (Cottoidei: Zoarcales). To achieve this, we carried out a comparative phylogenetic and structural analysis of 29 CR sequences belonging to the long shanny Stichaeus grigorjewi together with nine sequences of other eelpouts taxa representing four families in contrast to mitochondrial protein-coding fragments. The CR organization within S. grigorjewi, as well as in all other eelpouts, is consistent with the common three-domain structure known from most vertebrates. We found a hidden CR variation constrains on the landscape level and a lack of nucleotide saturation. Finally, our results demonstrate the advantage of the length variation in CR sequences for phylogenetic reconstructions among eelpouts.

Multiple independent structural dynamic events in the evolution of snake mitochondrial genomes

BACKGROUND

Mitochondrial DNA sequences have long been used in phylogenetic studies. However, little attention has been paid to the changes in gene arrangement patterns in the snake's mitogenome. Here, we analyzed the complete mitogenome sequences and structures of 65 snake species from 14 families and examined their structural patterns, organization and evolution. Our purpose was to further investigate the evolutionary implications and possible rearrangement mechanisms of the mitogenome within snakes.

RESULTS

In total, eleven types of mitochondrial gene arrangement patterns were detected (Type I, II, III, III-A, III-B, III-B1, III-C, III-D, III-E, III-F, III-G), with mitochondrial genome rearrangements being a major trend in snakes, especially in Alethinophidia. In snake mitogenomes, the rearrangements mainly involved three processes, gene loss, translocation and duplication. Within Scolecophidia, the O_L was lost several times in Typhlopidae and Leptotyphlopidae, but persisted as a plesiomorphy in the Alethinophidia. Duplication of the control region and translocation of the tRNA^Leu gene are two visible features in Alethinophidian mitochondrial genomes. Independently and stochastically, the duplication of pseudo-Pro (P*) emerged in seven different lineages of unequal size in three families, indicating that the presence of P* was a polytopic event in the mitogenome.

CONCLUSIONS

The WANCY tRNA gene cluster and the control regions and their adjacent segments were hotspots for mitogenome rearrangement. Maintenance of duplicate control regions may be the source for snake mitogenome structural diversity.

Complete mitochondrial genome sequences of the northern spotted owl (Strix occidentalis caurina) and the barred owl (Strix varia; Aves: Strigiformes: Strigidae) confirm the presence of a duplicated control region

We report here the successful assembly of the complete mitochondrial genomes of the northern spotted owl (Strix occidentalis caurina) and the barred owl (S. varia). We utilized sequence data from two sequencing methodologies, Illumina paired-end sequence data with insert lengths ranging from approximately 250 nucleotides (nt) to 9,600 nt and read lengths from 100–375 nt and Sanger-derived sequences. We employed multiple assemblers and alignment methods to generate the final assemblies. The circular genomes of S. o. caurina and S. varia are comprised of 19,948 nt and 18,975 nt, respectively. Both code for two rRNAs, twenty-two tRNAs, and thirteen polypeptides. They both have duplicated control region sequences with complex repeat structures. We were not able to assemble the control regions solely using Illumina paired-end sequence data. By fully spanning the control regions, Sanger-derived sequences enabled accurate and complete assembly of these mitochondrial genomes. These are the first complete mitochondrial genome sequences of owls (Aves: Strigiformes) possessing duplicated control regions. We searched the nuclear genome of S. o. caurina for copies of mitochondrial genes and found at least nine separate stretches of nuclear copies of gene sequences originating in the mitochondrial genome (Numts). The Numts ranged from 226–19,522 nt in length and included copies of all mitochondrial genes except tRNA^Pro, ND6, and tRNA^Glu. Strix occidentalis caurina and S. varia exhibited an average of 10.74% (8.68% uncorrected p-distance) divergence across the non-tRNA mitochondrial genes.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

Population structure and genetic diversity of Sinibrama macrops from Ou River and Ling River based on mtDNA D-loop region analysis, China

In order to understand the influence of human activities such as habitat fragmentation on freshwater fish population evolution, we investigated and compared the genetic diversity and phylogeography of Sinibrama macrops populations in the Oujiang River and Ling River. Mitochondrial control region sequences (D-loop region) of 131 specimens from six populations were obtained and analyzed. The diversity of main stream in the Ou River was lower than that in Ling River. Changtan population showed the lowest diversity (H = 0.646 ± 0.077; π = 0.00060 ± 0.00820). Pairwise F_ST, gene flow (Nm), and genetic distance (Da) indicated that Longquan and Changtan significantly differentiate from other populations. Nested clade phylogeographical analysis (NCPA) showed some clades and total cladogram experienced isolation by distance. In conclusion, the populations from severely fragmented Ou River have the lower diversity and more intense differentiation than that from the mainstream of Ling River, Changtan population present the lowest diversity and were isolated by the dam construction.

Selfish drive can trump function when animal mitochondrial genomes compete

Mitochondrial genomes compete for transmission from mother to progeny. We explored this competition by introducing a second genome into Drosophila melanogaster to follow transmission. Competitions between closely related genomes favored those functional in electron transport, resulting in a host-beneficial purifying selection. In contrast, matchups between distantly related genomes often favored those with negligible, negative or lethal consequences, indicating selfish selection. Exhibiting powerful selfish selection, a genome carrying a detrimental mutation displaced a complementing genome, leading to population death after several generations. In a different pairing, opposing selfish and purifying selection counterbalanced to give stable transmission of two genomes. Sequencing of recombinant mitochondrial genomes showed that the noncoding region, containing origins of replication, governs selfish transmission. Uniparental inheritance prevents encounters between distantly related genomes. Nonetheless, in each maternal lineage, constant competition among sibling genomes selects for super-replicators. We suggest that this relentless competition drives positive selection, promoting change in the sequences influencing transmission.

Assembly and variation analyses of Clarias batrachus mitogenome retrieved from WGS data and its phylogenetic relationship with other catfishes

Whole genome sequencing (WGS) using next generation sequencing technologies paves the way to sequence the mitochondrial genomes with greater ease and lesser time. Here, we used the WGS data of Clarias batrachus, generated from Roche 454 and Ion Torrent sequencing platforms, to assemble the complete mitogenome using both de novo and reference based approaches. Both the methods yielded almost similar results and the best assembled mitogenome was of 16,510 bp size (GenBank Acc. No. KM259918). The mitogenome annotation resulted in 13 coding genes, 22 tRNA genes, 2 rRNA genes and one control region, and the gene order was found to be identical with other catfishes. Variation analyses between assembled and the reference (GenBank Acc. No. NC_023923) mitogenome revealed 51 variations. The phylogenetic analysis of coding DNA sequences and tRNA supports the monophyly of catfishes. Two SSRs were identified in C. batrachus mitogenome, out of which one was unique to this species. Based on the relative rate of gene evolution, protein coding mitochondrial genes were found to evolve at a much faster pace than the d-loop, which in turn are followed by the rRNAs; the tRNAs showed wide variability in the rate of sequence evolution, and on average evolve the slowest. Among the coding genes, ND2 evolves most rapidly. The variations present in the coding regions of the mitogenome and their comparative analyses with other catfish species may be useful in species conservation and management programs.

The complete mitochondrial genome of pike-perch, Sander lucioperca (Perciformes: Percidae)

The complete mitochondrial genome of Sander lucioperca has been sequenced and analyzed in this study. It was a circular double-stranded DNA molecule of 16,541 base pairs (bp) in length and exhibited 37 typical animal mitochondrial genes. The gene order and base composition were similar to those of other percid species. All protein-coding genes were initiated with ATG except for COX 1, which began with GTG instead. However, the termination codons of 13 protein-coding genes varied with TAG, TAA, TA or T. Within CR, we detected five copies of 10 bp tandemly repeated sequences domain, which immediately followed the tRNA(Pro). These mitogenome sequence data would contribute to better understanding phylogenetic relationships and population genetics of the family Percidae.

Genetic multipartitions based on D-Loop sequences and chromosomal patterns in brown chromis, Chromis multilineata (Pomacentridae), in the Western Atlantic

Connectivity levels among Brazilian reef fish fauna populations have attracted growing interest, mainly between mainland shores and oceanic islands. The Pomacentridae, whose phylogeographic patterns are largely unknown in the Atlantic, are a family of dominant fish in reef regions. We present data on the variability and population structure of damselfish Chromis multilineata in different areas along the northeast coast of Brazil and in the waters around the oceanic islands of Fernando de Noronha (FNA) and Saint Peter and Saint Paul Archipelago (SPSPA) through analysis of the HVR1 mtDNA sequence of the control region. The remote SPSPA exhibits the highest level of genetic divergence among populations. Conventional and molecular cytogenetic analysis showed similar karyotype patterns (2n = 48 acrocentrics) between these insular areas. Our estimates reveal three genetically different population groups of C. multilineata on the Brazilian coast. The level of genetic structure is higher than previous data suggested, indicating complex panel of interactions between the oceanic island and coastal populations of Brazil.

Control region length dynamics potentially drives amino acid evolution in tarsier mitochondrial genomes

Patterns and processes of molecular evolution critically influence inferences in phylogeny and phylogeography. Within primates, a shift in evolutionary rates has been identified as the rationale for contrasting findings from mitochondrial and nuclear DNA studies as to the position of Tarsius. While the latter now seems settled, we sequenced complete mitochondrial genomes of three Sulawesi tarsiers (Tarsius dentatus, T. lariang, and T. wallacei) and analyzed substitution rates among tarsiers and other primates to infer driving processes of molecular evolution. We found substantial length polymorphism of the D-loop within tarsier individuals, but little variation of predominant lengths among them, regardless of species. Length variation was due to repetitive elements in the CSB domain-minisatellite motifs of 35 bp length and microsatellite motifs of 6 bp length. Amino acid evolutionary rates were second highest among major primate taxa relative to nucleotide substitution rates. We observed many radical possibly function-altering amino acid changes that were rarely driven by positive selection and thus potentially slightly deleterious or neutral. We hypothesize that the observed pattern of an increased amino acid evolutionary rate in tarsier mitochondrial genomes may be caused by hitchhiking of slightly deleterious mutations with favored D-loop length variants selected for maximizing replication success within the cell or the mitochondrion.

Geographic variation within a tandemly repeated mitochondrial DNA D-loop region of a North American freshwater fish, Pylodictis olivaris

The present study reports the distribution of a 35-bp mitochondrial DNA (mtDNA) D-loop tandemly repeated sequence in the populations of a North American freshwater catfish, Pylodictis olivaris, and the important role of a past geological event in the phylogeographic pattern of this species. A total of 330 individuals of flathead catfish, representing 34 drainages throughout the species' native range in the United States, were collected. While more than 70% of individuals sampled from the Southeastern Gulf Coast drainages were characterized by the presence of a 35-bp mtDNA D-loop tandem repeat proximal to the 5' end, more than 95% of samples from the Mississippi River and its tributaries, as well as from the drainages of the Southwest Gulf Coast region, lack this tandem repeat. Concomitantly, phylogenetic analyses revealed the existence of two distinct matrilineal lineages (lineage I and II) of P. olivaris, which were estimated to have diverged from a common ancestor sometime between 0.70 and 2.05myr ago. While one lineage is comprised of samples from the Mississippi River and its tributaries and rivers draining to the Southwest Gulf Coast, the other lineage is comprised of samples from the Southeastern Gulf Coast drainages. Each lineage also has two sub-lineages, which also showed geographic specificity.

Live sharksucker Echeneis naucrates (Linnaeus 1758) mitochondrial genome: the first report of Echeneidae complete mitochondrial genome

The complete mitochondrial genome of Echeneis naucrates is 16,611 bp in length. It comprises a control region, 13 protein-coding genes, 2 ribosomal RNAs (rRNAs) and 22 transfer RNAs (tRNAs), with an arrangement typical of vertebrate mitochondrial genomes. Base composition on the heavy strand is 30.24% A, 25.45% C, 15.02% G and 29.29% T. The control region is 940 bp in length, containing putative termination associated sequences (TASs) and conserved sequence blocks (CSBs). Two copies of a tandem repeat (AATATTAT) were found in all six individuals investigated. The hypothesis of selection for an optimal number of repeats as well as the evolutionary dynamics of tandem repeats in E. naucrates control regions await further investigations.

A comprehensive description and evolutionary analysis of 22 grouper (perciformes, epinephelidae) mitochondrial genomes with emphasis on two novel genome organizations

Groupers of the family Epinephelidae are a diverse and economically valuable group of reef fishes. To investigate the evolution of their mitochondrial genomes we characterized and compared these genomes among 22 species, 17 newly sequenced. Among these fishes we identified three distinct genome organizations, two of them never previously reported in vertebrates. In 19 of these species, mitochondrial genomes followed the typical vertebrate canonical organization with 13 protein-coding genes, 22 tRNAs, two rRNAs, and a non-coding control region. Differing from this, members of genus Variola have an extra tRNA-Ile between tRNA-Val and 16S rRNA. Evidence suggests that this evolved from tRNA-Val via a duplication event due to slipped strand mispairing during replication. Additionally, Cephalopholisargus has an extra tRNA-Asp in the midst of the control region, likely resulting from long-range duplication of the canonical tRNA-Asp through illicit priming of mitochondrial replication by tRNAs. Along with their gene contents, we characterized the regulatory elements of these mitochondrial genomes' control regions, including putative termination-associated sequences and conserved sequence blocks. Looking at the mitochondrial genomic constituents, rRNA and tRNA are the most conserved, followed by protein-coding genes, and non-coding regions are the most divergent. Divergence rates vary among the protein-coding genes, and the three cytochrome oxidase subunits (COI, II, III) are the most conserved, while NADH dehydrogenase subunit 6 (ND6) and the ATP synthase subunit 8 (ATP8) are the most divergent. We then tested the phylogenetic utility of this new mt genome data using 12 protein-coding genes of 48 species from the suborder Percoidei. From this, we provide further support for the elevation of the subfamily Epinephelinae to family Epinephelidae, the resurrection of the genus Hyporthodus, and the combination of the monotypic genera Anyperodon and Cromileptes to genus Epinephelus, and Aethaloperca to genus Cephalopholis.

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: an ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8-10 length variants. The motifs of the HTR arrays are about 353-362 bp and the number of repeats ranges from 1 to 11.

Limited phylogenetic distribution of a long tandem-repeat cluster in the mitochondrial control region in Bubo (Aves, Strigidae) and cluster variation in Blakiston's fish owl (Bubo blakistoni)

To investigate the phylogenetic position of Blakiston's fish owl (Bubo blakistoni), we sequenced the mitochondrial (mt) DNA control region and cytochrome b (cyt b) for nine Bubo species. Maximum-likelihood analyses of combined control region and cyt b sequences, and cyt b sequences alone, showed that species formerly placed in genus Ketupa comprise a monophyletic group. Unexpectedly, we discovered a long cluster of 20-25 tandem repeat units 77 or 78bp long in the third control region domain in four of the nine Bubo species for which the control region was sequenced (B. blakistoni, B. flavipes, and B. ketupu in the Ketupa clade; B. lacteus), leading to overall control region lengths of 3.0-3.8kpb estimated from agarose gel electrophoresis. The control region in B. lacteus is the longest (3.8kbp) reported to date in vertebrates. Sequencing of eight repeat units at each end of the cluster in 20 B. blakistoni individuals detected several types of repeat units 77 or 78bp long, and six patterns in the order of unit types. The occurrence of a repeat cluster in all three species examined in the Ketupa clade suggests their common ancestor also had a cluster, whereas a maximum parsimony tree showed repeat-unit types grouping by species, rather than by paralog groups, suggesting independent origins of the clusters. We reconcile these results with a turnover model, in which the range in cluster-length variation and unit types at the 5' end are hypothetically functionally constrained by the protein-binding function of the control region, but otherwise there is a continual turnover of units in evolutionary time, with new unit types arising through mutations, proliferating by duplication of single and double repeat blocks, and being lost through deletion. Estimated free energies for reconstructed secondary structures of single and especially pairs of repeat units were higher than for homologous single-unit blocks in species lacking a repeat cluster, supporting slipped-strand mispairing as the mechanism of cluster turnover.

Mitochondrial DNA STR analysis as a tool for studying the green sea turtle (Chelonia mydas) populations: the Mediterranean Sea case study

The Mediterranean population of the green sea turtle Chelonia mydas is critically endangered. Genetic analysis of this population using the ordinary haplotyping system, based on sequence analysis of a segment of the mitochondrial DNA (mtDNA) D-loop (control region), revealed very little variation. The most common haplotype, CM-A13, was observed in all but three individuals in hundreds of samples in previous studies. In search for a more informative marker we sequenced the 3' of the mitochondrial control region which contains an AT-rich microsatellite. We found a unique pattern that consists of four AT short tandem repeats (STRs) with varying copy numbers. This allowed us to construct a new haplotyping system composed of four different STR sizes for each mtDNA sequence. Our new mitochondrial STR (mtSTR) haplotyping approach revealed 33 different haplotypes within the nesting and stranded sea turtles along the Mediterranean Israeli seashore. The Israeli coast nesting females had 10 different haplotypes that can be used for monitoring and conservation purposes. The mtSTR haplotyping system can clearly assist in fingerprinting of individual turtles. Moreover, it can be used for estimating phylogenetic distances within populations. This case study shows that the mtSTR haplotyping is applicable for the study of global green sea turtle populations and could also be considered as markers of genetic variability in other species.

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

BACKGROUND

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

PRINCIPAL FINDINGS

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

CONCLUSIONS/SIGNIFICANCE

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.

Complete mitochondrial DNA sequence of Crassostrea nippona: comparative and phylogenomic studies on seven commercial Crassostrea species

The complete mitochondrial genome of Crassostrea nippona was determined and compared with six other Crassostrea mitogenomes from GenBank in an attempt to shed light on the evolutionary relatedness within Crassostrea. The total length of the mitogenome was 20,030 bp for C. nippona, which was the largest among seven Crassostrea mitogenomes. Among six Asian oysters, the gene order of mitochondrial DNA was identical except for C. nippona with a transposition of trnG. While the American oyster C. virginica and Asian oysters showed broad differences in gene order with relocation of most tRNA genes and indels of duplicated tRNAs and rrnS, indicating the relatively distant relationships between the American oyster and Asian oysters. Different from other six Crassostrea oysters, C. nippona had two repeats of 66 bp in non-coding regions. Pairwise divergence among the seven Crassostrea oysters based on DNA sequences of 12 protein-coding genes ranged from 3.1 to 44.4% (Kimura two-parameter distance, K2P). The close relationship between C. nippona and C. hongkongensis was revealed by K2P of 18.9%. Phylogenetic analyses robustly revealed Crassostrea monophyly, with C. virginica at the basal position. The results of phylogenetic analyses strongly supported C. gigas and C. angulata had the closest relationship, with C. sikamea being the sister taxon. These findings presented here provide a better insight into the relationships within Crassostrea and will be useful for further evolution studies of oysters.

[Analysis of complete mitochondrial genome of Cristaria plicata]

The complete mitochondrial genome of Cristaria plicata was obtained using long amplification polymerase chain reaction (LA-PCR). Analysis showed that the total length of sequence was 15 712 bp, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 26 non-coding regions ranged from 2 bp to 328 bp in size. The nucleotide composition of A, T, C, G were 36.54% 27.22% 23.22% and 13.02%, respectively. Most genes were encoded on the L strand, While ND3-ND5, ND4L, COI-COIII, ATP6, ATP8, tRNAAsp, and tRNAHis were encoded on the H strand. The arrangement of genes in C. plicata was identical to that of Lampsilis ornata, but was different from that of Hyriopsis cum-ingii between COII and 12S rRNA in the same family. Thirteen protein genes contained 3 initiation codons, i.e., I (AUU, AUC), V (GUG), and M (AUA, AUG), all of which had complete stop codons (UAA or UAG), except for the stop codon in ND2 that had an incomplete T. Fifteen tRNAs had typical cloverleaf structure, except for tRNAThr, tRNALys, tRNASer(UCN), tRNAAsp, tRNAArg, tRNATyr, and tRNAMet. Like other freshwater bivalvia, C. plicata contained ATP8 gene, which might be related to the balance of osmotic pressure in the cytoplasm.

Complete mitochondrial DNA sequence and phylogenetic analysis of Zhikong scallop Chlamys farreri (Bivalvia: Pectinidae)

The complete mitochondrial genome of Zhikong scallop Chlamys farreri is 21,695 bp in length and contains 12 protein-coding genes (the atp8 gene is absent, as in most bivalves), 2 ribosomal RNA genes, and 22 transfer RNA genes. The heavy strand has an overall A+T content of 58.7%. GC and AT skews for the mt genome of C. farreri are 0.337 and -0.184, respectively, indicating the nucleotide bias against C and A. The mitochondrial gene order of C. farreri differs drastically from the scallops Argopecten irradians, Mimachlamys nobilis and Placopecten magellanicus, which belong to the same family Pectinidae. 6623 bp non-coding nucleotides exist intergenically in the mitogenome of C. farreri, with a large continuous sequence (4763 bp) between tRNA ( Val ) and tRNA ( Asn ). Two repeat families are found in the large continuous sequence, which seems to be a common feature of scallops. Phylogenetic analysis based on 12 concatenated amino acid sequences of protein-coding genes supports the monophyly of Pectinidae and paraphyletic Pteriomorphia with respect to Heteroconchia.

Structure and conservation of tandem repeats in the mitochondrial DNA control region of the Least Brook lamprey (Lampetra aepyptera)

Tandemly repeated sequences are a common feature of vertebrate mitochondrial DNA control regions. However, questions still remain about their mode of evolution, function, and phylogenetic distribution. We report phylogenetic and geographic patterns of variation of control region repeat sequence and number in a nonparasitic lamprey, Lampetra aepyptera. A survey of populations from throughout the species' range revealed remarkably low repeat sequence polymorphism but some interpopulation variation in repeat number. The high sequence similarity extended to repeats observed in other species in the genus Lampetra and other lamprey genera. The very low levels of variation suggest a high copy turnover. Our data are consistent with the illegitimate elongation model of repeat gain and loss and further suggest that repeat change occurs at internal copies. However, the limited variation across some species of lamprey suggests that functional constraints may further limit variation.

Comparative analysis of complete mitochondrial DNA control region of four species of Strigiformes

The sequence of the whole mitochondrial (mt) DNA control region (CR) of four species of Strigiformes was obtained. Length of the CR was 3,290 bp, 2,848 bp, 2,444 bp, and 1,771 bp for Asio flammeus, Asio otus, Athene noctua, and Strix aluco, respectively. Interestingly, the length of the control region was maximum in Asio flammeus among all the avian mtDNA control regions sequenced thus far. In addition, the base composition and organization of mtDNA CR of Asio flammeus were identical to those reported for other birds. On the basis of the differential frequencies of base substitutions, the CR may be divided two variable domains, I and III, and a central conserved domain, II. The 3' end of the CR contained many tandem repeats of varying lengths and repeat numbers. In Asio flammeus, the repeated sequences consisted of a 126 bp sequence that was repeated seven times and a 78 bp sequence that was repeated 14 times. In Asio otus, there were also two repeated sequences, namely a 127 bp sequence that was repeated eight times and a 78 bp sequence that was repeated six times. The control region of Athene noctua contained three sets of repeats: a 89 bp sequence that was repeated three times, a 77 bp sequence that was repeated four times, and a 71 bp sequence that was repeated six times. Strix aluco, however, had only one repeated sequence, a 78 bp sequence that was repeated five times. The results of this study seem to indicate that these tandem repeats may have resulted from slipped-strand mispairing during mtDNA replication. Moreover, there are many conserved motifs within the repeated units. These sequences could form stable stem-loop secondary structures, which suggests that these repeated sequences play an important role in regulating transcription and replication of the mitochondrial genome.

Clear genetic structure of Pinus kwangtungensis (Pinaceae) revealed by a plastid DNA fragment with a novel minisatellite

BACKGROUND AND AIMS

Pinus kwangtungensis is a five-needled pine, inhabiting isolated mountain tops, cliffs or slopes in the montane areas of southern China and northern Vietnam. Global warming and long-term deforestation in southern China threaten its existence and genetic integrity, and this species is listed as vulnerable in the China Species Red List. However, the level and distribution of genetic diversity in this vulnerable species are completely unknown. In this paper, the genetic diversity and structure are examined using paternally inherited plastid markers to shed light on its evolutionary history and to provide a genetic perspective for its conservation.

METHODS

By means of direct sequencing, a new polymorphic fragment containing a minisatellite site was identified within the plastid genome of P. kwangtungensis. Using the minisatellite site along with five SNPs (one indel and four substitutions) within the same fragment, the population genetic structure and pollen flow were analysed in 17 populations of P. kwangtungensis in southern China.

KEY RESULTS

Analysis of 227 individuals from 17 populations revealed ten haplotypes at the minisatellite site. The haplotype diversity at species level was relatively high (0.629). Genetic diversity of each population ranged from 0 to 0.779, and the western populations harboured more genetic variation than the eastern and Hainan populations, although the former appeared to have experienced a bottleneck in recent history. Population subdivision based on this site was high (F(ST) = 0.540 under IAM; R(ST) = 0.677 under SMM). Three major clusters (eastern, western and Hainan) were identified based on a neighbor-joining dendrogram generated from genetic distances among the populations. The genetic structures inferred from all the polymorphic sites and the SNPs were in concordance with that from the minisatellite site.

CONCLUSIONS

The results suggest that there are at least three refugia for P. kwangtungensis and that populations in these refugia should be treated as separate evolutionarily significant units or conservation units. The high diversities in the western populations suggest that these were much larger in the past (e.g. glacial stages) and that the shrinking population size might have been caused by recent events (e.g. deforestation, global warming, etc.). The western populations should be given priority for conservation due to their higher genetic diversity and limited population sizes. It is concluded that the newly found minisatellite may serve as a novel and applicable molecular marker for unravelling evolutionary processes in P. kwangtungensis.

The complete mitochondrial genome of the Senegal sole, Solea senegalensis Kaup. Comparative analysis of tandem repeats in the control region among soles

The complete nucleotide sequence of the mitochondrial genome for the Senegal sole Solea senegalensis Kaup was determined. The mitochondrial DNA was 16,659 base pairs (bp) in length. Sequence features of the 13 protein-coding genes, two ribosomal RNAs and 22 transfer RNAs are described. The non-coding control region (1017 bp) was compared with those of the closely related soles Solea solea and Solea lascaris. The typical conservative blocks were identified. A cluster of 42 and 22 tandemly arrayed repeats was detected near the 3' end of control region in S. solea and S. lascaris, respectively. On the contrary, only two (93.8% of haplotypes) or three copies (6.2%) of an 8-bp repeated sequence motif was found in S. senegalensis. Phylogenetic analysis showed that 7 out of 9 of haplotypes bearing three copies grouped in a separate cluster. Possible mechanisms influencing the evolution of control region among soles are discussed.

Structure and evolution of the mitochondrial control region and flanking sequences in the European cave salamander Proteus anguinus

The European cave salamander Proteus anguinus Laurenti 1768 is one of the best-known subterranean animals, yet its evolutionary history and systematic relationships remain enigmatic. This is the first comprehensive study on molecular evolution within the taxon, using an mtDNA segment containing the control region (CR) and adjacent sequences. Two to seven tandem repeats of 24-32 bp were found in the intergenic spacer region (VNTR1), and three, four or six repeats, 59-77 bp each, in the 3' end of the CR (VNTR2). Different molecular mechanisms account for VNTR2 formation in different lineages of Proteus. The overall CR variation was lower than that of the spacer region, the 3' end of the cytb gene, or the tRNA genes. Individual genes and the concatenated non-repetitive sequences produced similar, well resolved maximum likelihood, Bayesian inference and parsimony trees. The numbers of repeat elements as well as the genealogy of the VNTR2 repeat units were mostly inconsistent with the groupings of the non-repetitive sequences. Different degrees of repeat array homogenization were detected in all major groups. Orthology was established for the first and the second VNTR2 elements of some populations. These two copies may therefore be used for analyses at the population level. The pattern of CR sequence variation points to strong genetic isolation of hydrographically separated populations. Genetic separation of the major groups of populations is incongruent with the current division into subspecies.

Characterization of a mosaic minisatellite locus in the mitochondrial DNA of Norway spruce [Picea abies (L.) Karst.]

A mosaic minisatellite region has been identified in the mitochondrial genome of Norway spruce (Picea abies). The array was composed of three tandem repeats PaTR1 (32 bp), PaTR2a (26 bp) and PaTR2b (26 bp). PaTR2a and PaTR2b differed by one base substitution. The analysis of 92 trees covering the whole natural distribution area of the species allowed detection of 11 length variants ranging from 131 bp to 447 bp. This high intra-specific polymorphism relies on variation in the number of the tandem repeats. Population genetic parameters estimated among 14 populations suggested high population differentiation (Gst=0.749). The phylogenetic analysis of the 11 sequenced length variants has been performed using a parsimony approach. The topology of the tree showed a good association of groups with geographical origin and a low level of size homoplasy. The phylogenetic reconstruction also suggests that this minisatellite locus has mainly evolved by an increase in the repeat copy number.

Species boundaries, populations and colour morphs in the coral reef three-spot damselfish (Dascyllus trimaculatus) species complex

Coloration patterns of tropical reef fishes is commonly used for taxonomic purposes, yet few studies have focused on the relationship between species boundaries and coloration types. The three-spot damselfish (Dascyllus trimaculatus) species complex comprises four species that vary both in geographical ranges and colour patterns making them an ideal model to study these relationships. We analysed the mitochondrial control region of 122 individuals from all four species collected from 13 localities. Individuals from two species (Dascyllus albisella and D. strasburgi) grouped into monophyletic clades, while the two other species (D. trimaculatus and D. auripinnis) were found to be paraphyletic. Coloration patterns were therefore not found to be good predictors of genetic isolation. In contrast, geographical origin was always consistent with the observed genetic pattern.