Length variation, heteroplasmy and sequence divergence in the mitochondrial DNA of four species of sturgeon (Acipenser)

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

First photographic evidence of Indian Pangolin Manis crassicaudata Geoffroy, 1803 (Mammalia: Pholidota: Manidae), in Colonel Sher Jung National Park, Himachal Pradesh, India

The Indian Pangolin Manis crassicaudata is ‘Endangered’ on the IUCN Red List. We report the first photographic evidence of its presence in Colonel Sher Jung National Park of Paonta Valley in Himachal Pradesh. Camera-traps (n = 64) were used to assess the diversity and abundance of terrestrial mammals. Over 1,912 trap nights, we documented 16 mammals, of which the Indian Pangolin was photo-captured at a single camera-trap station with two independent records. Although the presence of species has been indicated in a few studies conducted in the study area, no photographic evidence of the Indian Pangolin from the National Park has previously been reported. We strongly recommend additional fieldwork in and around the protected area to learn more about the Indian Pangolin’s distribution, habitat use, and ecology.

Comparative analyses of the complete mitochondrial genomes of two southern African endemic guitarfish, Acroteriobatus annulatus and A. blochii

Shark-like rays (order Rhinopristiformes) are among the most threatened cartilaginous fish globally. Despite this, unresolved taxonomic issues still exist within the group. To date, no studies have used complete mitochondrial genomes to assess the phylogenetic placement of Acroteriobatus within the non-monophyletic family Rhinobatidae. The current study reports the first complete mitochondrial genomes for Acroteriobatus annulatus and A. blochii. Similar to other rhinopristiforms, the complete sequences of A. annulatus (16,773 bp) and A. blochii (16,771 bp) were circular molecules with gene organisations identical to that of the typical vertebrate mitogenome. The A + T content was higher than the G + C content, with a bias towards A and C nucleotides observed in all complete mitogenomes. The stem-and-loop secondary structures of the putative origin of light-strand replication were found to have highly conserved synthesis and stem regions, with all substitutions and indels restricted to the loop structure. The ratios of non-synonymous to synonymous substitution rates indicated that purifying selection has been the dominant driver of evolution in rhinopristiform mitogenomes. Phylogenetic reconstructions placed Acroteriobatus as a sister-group to Rhinobatos, confirming its affiliation with the family Rhinobatidae. However, based on its apparent polyphyly with the aforementioned genera, the familial assignment of Pseudobatos is not fully resolved and requires further investigation.

A New Mitochondrial Genome of Sogatella furcifera (Horváth) (Hemiptera: Delphacidae) and Mitogenome-Wide Investigation on Polymorphisms

Simple Summary

We completed one mitogenome of white-backed planthopper (WBPH), Sogatella furcifera (Horváth), with finding heteroplasmy phenomenon confirmed by PCR reaction and Sanger sequencing method. This heteroplasmy was not observed in WBPHs (n = 24) collected from the fields, suggesting that it may be uncommon in fields. We also analyzed single nucleotide polymorphisms, insertion and deletions, and simple sequence repeats among three currently available WBPH mitogenomes of Korea and China, suggesting that identified intraspecific variations could be potential candidates for developing markers to distinguish geographical populations of WBPH including Korean and Chinese. Phylogenetic analysis of 32 mitogenomes of Delphacidae including the three WBPH mitogenomes suggested that Delphacinae seems to be monophyletic and Sogatella species including WBPH are clearly formed as one clade.

Abstract

White-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the major sap-sucking rice pests in East Asia. We have determined a new complete mitochondrial genome of WBPH collected in the Korean peninsula using NGS technology. Its length and GC percentages are 16,613 bp and 23.8%, respectively. We observed one polymorphic site, a non-synonymous change, in the COX3 gene with confirmation heteroplasmy phenomenon within individuals of WBPH by PCR amplification and Sanger sequencing, the first report in this species. In addition, this heteroplasmy was not observed in wild WBPH populations, suggesting that it may be uncommon in fields. We analyzed single nucleotide polymorphisms, insertion, and deletions, and simple sequence repeats among the three WBPH mitogenomes from Korea and China and found diverse intraspecific variations, which could be potential candidates for developing markers to distinguish geographical populations. Phylogenetic analysis of 32 mitogenomes of Delphacidae including the three WBPH mitogenomes suggested that Delphacinae seems to be monophyletic and Sogatella species including WBPH are clearly formed as one clade. In the future, it is expected that complete mitogenomes of individuals of geographically dispersed WBPH populations will be used for further population genetic studies to understand the migration pathway of WBPH.

Complete Mitochondrial DNA Genome of Nine Species of Sharks and Rays and Their Phylogenetic Placement among Modern Elasmobranchs

Chondrichthyes occupy a key position in the phylogeny of vertebrates. The complete sequence of the mitochondrial genome (mitogenome) of four species of sharks and five species of rays was obtained by whole genome sequencing (DNA-seq) in the Illumina HiSeq2500 platform. The arrangement and features of the genes in the assembled mitogenomes were identical to those found in vertebrates. Both Maximum Likelihood (ML) and Bayesian Inference (BI) analyses were used to reconstruct the phylogenetic relationships among 172 species (including 163 mitogenomes retrieved from GenBank) based on the concatenated dataset of 13 individual protein coding genes. Both ML and BI analyses did not support the “Hypnosqualea” hypothesis and confirmed the monophyly of sharks and rays. The broad notion in shark phylogeny, namely the division of sharks into Galeomorphii and Squalomorphii and the monophyly of the eight shark orders, was also supported. The phylogenetic placement of all nine species sequenced in this study produced high statistical support values. The present study expands our knowledge on the systematics, genetic differentiation, and conservation genetics of the species studied, and contributes to our understanding of the evolutionary history of Chondrichthyes.

Population genetic structure and phylogeography of sterlet (Acipenser ruthenus, Acipenseridae) in the Ob and Yenisei river basins

The sterlet (Acipenser ruthenus Linnaeus, 1758) is a relatively small sturgeon widely distributed in Eurasian rivers from the Danube to the Yenisei. During the twentieth century, all wild sterlet populations have declined due to anthropogenic factors including: overfishing, poaching, construction of dams, and pollution. Despite the necessity of characterization both wild and captive stocks, few studies of population genetics have been performed thus far. Here we studied the genetic diversity and geographic structure of sterlet populations across the eastern range - Ob-Irtysh and Yenisei basins - by sequencing a 628-bp fragment of mitochondrial DNA control region. We identified 98 new haplotypes, delineated 12 haplogroups and estimated the time of basal haplogroup divergence within the species as over 8 million years ago. Our data suggest that Ob-Irtysh and Yenisei populations are isolated from each other and much lower genetic diversity is present in the Yenisei population than in the Ob-Irtysh population. Our data imply that sterlet populations in Siberian rivers underwent bottleneck or fragmentation, followed by subsequent population expansion. The data obtained here are important for sterlet population monitoring and restocking management.

PARALLEL EVOLUTION OF LAKE-STREAM PAIRS OF THREESPINE STICKLEBACKS (GASTEROSTEUS) INFERRED FROM MITOCHONDRIAL DNA VARIATION

Three drainage systems in British Columbia, Canada, contain divergent parapatric lake-stream pairs of threespine sticklebacks (Gasterosteus aculeatus): Drizzle and Mayer Lakes on Graham Island, Queen Charlotte Islands, and Misty Lake on northeastern Vancouver Island. Ecological and morphological differences between members of all three lake-stream pairs are strikingly similar; lake fish are melanistic and slim bodied with smaller mouths and more gill rakers than the mottled-brown and robust-bodied stream sticklebacks. We estimated the level of genetic divergence between lake and stream fish in Misty Lake and tested hypotheses of single versus multiple origins of the pairs by assaying mitochondrial DNA (mtDNA) restriction site variation in samples from the three lake systems. MtDNA analysis revealed the existence of two highly divergent lineages differing by 2.7% in sequence. One lineage predominated in Misty stream fish (73%), whereas the other lineage predominated in Misty Lake samples (96%). Comparable forms (lake or stream) in the different lakes did not cluster together in terms of mtDNA nucleotide divergence, suggesting that the pairs have had independent origins. We concluded that: (1) divergent mtDNA lineages in North Pacific sticklebacks stem from historical isolation in the two major glacial refugia proposed for the North Pacific (Beringia and Cascadia); (2) the stream and lake pair in Misty Lake are distinct gene pools; (3) the divergence between parapatric lake and stream Gasterosteus represents parallel evolution having occurred at least twice in the North Pacific; and (4) different scales of evolutionary divergence exist in North Pacific Gasterosteus, that is, a relatively ancient divergence of mtDNA clades as well as recent (i.e., postglacial) divergence of ecotypes within major clades.

Mitochondrial heteroplasmy in vertebrates using ChIP-sequencing data

BACKGROUND

Mitochondrial heteroplasmy, the presence of more than one mitochondrial DNA (mtDNA) variant in a cell or individual, is not as uncommon as previously thought. It is mostly due to the high mutation rate of the mtDNA and limited repair mechanisms present in the mitochondrion. Motivated by mitochondrial diseases, much focus has been placed into studying this phenomenon in human samples and in medical contexts. To place these results in an evolutionary context and to explore general principles of heteroplasmy, we describe an integrated cross-species evaluation of heteroplasmy in mammals that exploits previously reported NGS data. Focusing on ChIP-seq experiments, we developed a novel approach to detect heteroplasmy from the concomitant mitochondrial DNA fraction sequenced in these experiments.

RESULTS

We first demonstrate that the sequencing coverage of mtDNA in ChIP-seq experiments is sufficient for heteroplasmy detection. We then describe a novel detection method for accurate detection of heteroplasmies, which also accounts for the error rate of NGS technology. Applying this method to 79 individuals from 16 species resulted in 107 heteroplasmic positions present in a total of 45 individuals. Further analysis revealed that the majority of detected heteroplasmies occur in intergenic regions.

CONCLUSION

In addition to documenting the prevalence of mtDNA in ChIP-seq data, the results of our mitochondrial heteroplasmy detection method suggest that mitochondrial heteroplasmies identified across vertebrates share similar characteristics as found for human heteroplasmies. Although largely consistent with previous studies in individual vertebrates, our integrated cross-species analysis provides valuable insights into the evolutionary dynamics of mitochondrial heteroplasmy.

Species identification through mitochondrial rRNA genetic analysis

Inter-species and intraspecific variations in mitochondrial DNA (mtDNA) were observed in a bioinformatics analysis of the mitochondrial genomic sequences of 11 animal species. Some highly conserved regions were identified in the mitochondrial 12S and 16S ribosomal RNA (rRNA) genes of these species. To test whether these sequences are universally conserved, primers were designed to target the conserved regions of these two genes and were used to amplify DNA from 21 animal tissues, including two of unknown origin. By sequencing these PCR amplicons and aligning the sequences to a database of non-redundant nucleotide sequences, it was confirmed that these amplicons aligned specifically to mtDNA sequences from the expected species of origin. This molecular technique, when combined with bioinformatics, provides a reliable method for the taxonomic classification of animal tissues.

Lack of Structural Variation but Extensive Length Polymorphisms and Heteroplasmic Length Variations in the Mitochondrial DNA Control Region of Highly Inbred Crested Ibis, Nipponia nippon

The animal mitochondrial DNA (mtDNA) length polymorphism and heteroplasmy are accepted to be universal. Here we report the lack of structural variation but the presence of length polymorphism as well as heteroplasmy in mtDNA control region of an endangered avian species - the Crested Ibis (Nipponia nippon). The complete control region was directly sequenced while the distribution pattern and inheritance of the length variations were examined using both direct sequencing and genotyping of the PCR fragments from captive birds with pedigrees, wild birds and a historical specimen. Our results demonstrated that there was no structural variation in the control region, however, different numbers of short tandem repeats with an identical motif of CA3CA2CA3 at the 3'-end of the control region determined the length polymorphisms among and heteroplasmy within individual birds. There were one to three predominant fragments in every bird; nevertheless multiple minor fragments coexist in all birds. These extremely high polymorphisms were suggested to have derived from the 'replication slippage' of a perfect microsatellite evolution following the step-wise mutational model. The patterns of heteroplasmy were found to be shifted between generations and among siblings but rather stable between blood and feather samples. This study provides the first evidence of a very extensive mtDNA length polymorphism and heteroplasmy in the highly inbred Crested Ibis which carries an mtDNA genome lack of structural genetic diversity. The analysis of pedigreed samples also sheds light on the transmission of mtDNA length heteroplasmy in birds following the genetic bottleneck theory. Further research focusing on the generation and transmission of particular mtDNA heteroplasmy patterns in single germ line of Crested Ibis is encouraged by this study.

Site heteroplasmy in the mitochondrial cytochrome b gene of the sterlet sturgeon Acipenser ruthenus

Sturgeons are fish species with a complex biology. They are also characterized by complex aspects including polyploidization and easiness of hybridization. As with most of the Ponto-Caspian sturgeons, the populations of Acipenser ruthenus from the Danube have declined drastically during the last decades. This is the first report on mitochondrial point heteroplasmy in the cytochrome b gene of this species. The 1141 bp sequence of the cytb gene in wild sterlet sturgeon individuals from the Lower Danube was determined, and site heteroplasmy evidenced in three of the 30 specimens collected. Two nucleotide sequences were identified in these heteroplasmic individuals. The majority of the heteroplasmic sites are synonymous and do not modify the sequence of amino acids in cytochrome B protein. To date, several cases of point heteroplasmy have been reported in animals, mostly due to paternal leakage of mtDNA. The presence of specific point heteroplasmic sites might be interesting for a possible correlation with genetically distinct groups in the Danube River.

Population genetic structure of Earth's largest fish, the whale shark (Rhincodon typus)

Large pelagic vertebrates pose special conservation challenges because their movements generally exceed the boundaries of any single jurisdiction. To assess the population structure of whale sharks (Rhincodon typus), we sequenced complete mitochondrial DNA control regions from individuals collected across a global distribution. We observed 51 single site polymorphisms and 8 regions with indels comprising 44 haplotypes in 70 individuals, with high haplotype (h = 0.974 +/- 0.008) and nucleotide diversity (pi = 0.011 +/- 0.006). The control region has the largest length variation yet reported for an elasmobranch (1143-1332 bp). Phylogenetic analyses reveal no geographical clustering of lineages and the most common haplotype was distributed globally. The absence of population structure across the Indian and Pacific basins indicates that oceanic expanses and land barriers in Southeast Asia are not impediments to whale shark dispersal. We did, however, find significant haplotype frequency differences (AMOVA, Phi(ST) = 0.107, P < 0.001) principally between the Atlantic and Indo-Pacific populations. In contrast to other recent surveys of globally distributed sharks, we find much less population subdivision and no evidence for cryptic evolutionary partitions. Discovery of the mating and pupping areas of whale sharks is key to further population genetic studies. The global pattern of shared haplotypes in whale sharks provides a compelling argument for development of broad international approaches for management and conservation of Earth's largest fish.

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.

Inheritance of microsatellite loci in the white sturgeon (Acipenser transmontanus)

Nine tetramer motif (GATA)n microsatellite systems were developed for use in the white sturgeon, Acipenser transmontanus. We report inheritance patterns for these nine systems, which range from one possible disomic system to tetrasomy and octosomy, with some systems containing null alleles. Because of the complex modes of inheritance underlying these systems and the highly duplicated nature of the genome, we propose each allele be scored as its own dominant marker, similar to AFLPs or RAPDs. The utility of this method is validated by the observation that individual alleles within a microsatellite system generally fit the expectation for independent transmission and fit the expected transmission frequency for single copy nuclear markers.

Population genetics of shortnose sturgeon Acipenser brevirostrum based on mitochondrial DNA control region sequences

Shortnose sturgeon is an anadromous North American acipenserid that since 1973 has been designated as federally endangered in US waters. Historically, shortnose sturgeon occurred in as many as 19 rivers from the St. John River, NB, to the St. Johns River, FL, and these populations ranged in census size from 10(1) to 10(4), but little is known of their population structure or levels of gene flow. We used the polymerase chain reaction (PCR) and direct sequence analysis of a 440 bp portion of the mitochondrial DNA (mtDNA) control region to address these issues and to compare haplotype diversity with population size. Twenty-nine mtDNA nucleotide-substitution haplotypes were revealed among 275 specimens from 11 rivers and estuaries. Additionally, mtDNA length variation (6 haplotypes) and heteroplasmy (2-5 haplotypes for some individuals) were found. Significant genetic differentiation (P < 0.05) of mtDNA nucleotide-substitution haplotypes and length-variant haplotypes was observed among populations from all rivers and estuaries surveyed with the exception of the Delaware River and Chesapeake Bay collections. Significant haplotype differentiation was even observed between samples from two rivers (Kennebec and Androscoggin) within the Kennebec River drainage. The absence of haplotype frequency differences between samples from the Delaware River and Chesapeake Bay reflects a probable current absence of spawning within the Chesapeake Bay system and immigration of fish from the adjoining Delaware River. Haplotypic diversity indices ranged between 0.817 and 0.641; no relationship (P > 0.05) was found between haplotype diversity and census size. Gene flow estimates among populations were often low (< 2.0), but were generally higher at the latitudinal extremes of their distribution. A moderate level of haplotype diversity and a high percentage (37.9%) of haplotypes unique to the northern, once-glaciated region suggests that northern populations survived the Pleistocene in a northern refugium. Analysis of molecular variance best supported a five-region hierarchical grouping of populations, but our results indicate that in almost all cases populations of shortnose sturgeon should be managed as separate units.

Ancient sex-specific extension of the cytochrome c oxidase II gene in bivalves and the fidelity of doubly-uniparental inheritance

Bivalves of the families Mytilidae and Unionidae show a unique mode of mitochondrial DNA inheritance called doubly uniparental inheritance. In addition to receiving the maternally transmitted mtDNA lineage, males receive a separate mtDNA genome from their fathers. This system is sometimes compromised, in that female genomes are occasionally recruited into the male cycle of inheritance. These masculinization events are common in the Mytilidae but have not been reported in the Unionidae. In order to estimate the age of the male and the female lineages in the Unionidae and to look for evidence of masculinization, we sequenced the junction between the cytochrome c oxidase II gene and the cytochrome c oxidase I gene. The unionid male and female lineages diverged approximately 450 MYA. There is no evidence for masculinization during this period, suggesting that there are taxon-specific differences in the rate of masculinization. Coincidentally, a 200-codon extension of the COII gene is present in the male genome of the Unionidae and may be responsible for the absence of masculinization.

Sturgeon genetics and cytogenetics: recent advancements and perspectives

The aim of this review is to introduce current knowledge in the field of sturgeon genetics. The first section deals with sturgeon cytogenetics, reviewing karyotype organization and polyploidization events during evolution of Acipenseriformes. The second section concerns the results of applications of molecular biology to studies of phylogenetic relationships between extant species, intraspecific analysis of wild populations and stocks for conservation purposes, together with characterization of molecular markers for species identification, relevant to forensic and conservation issues.

Evidence for a slowed rate of molecular evolution in the order acipenseriformes

A test of the hypothesis that the members of the order Acipenseriformes (sturgeons and paddlefishes) possess a slowed rate of molecular evolution was carried out by conducting relative-rate comparisons with representatives of four groups of teleost fishes (Cypriniformes, Elopomorpha, Salmonidae, and Percomorpha) using 21 nuclear or mitochondrial protein loci and the nuclear and mitochondrial small subunit rRNA genes, obtained from the literature or our own research. In 70 out of 81 comparisons between individual taxa (86%), acipenseriform sequences showed slower rates of change than the homologous teleost loci examined. When teleost sequences are considered together, 21 of the 23 loci show slower rates of substitution in the acipenseriform lineage. Teleost proteins show 1.85 times as many unique amino acid differences as acipenseriform proteins, when both are compared with outlier sequences. These results support a hypothesis of slowed molecular evolutionary rate in the Acipenseriformes.

Genome duplication events and functional reduction of ploidy levels in sturgeon (Acipenser, Huso and Scaphirhynchus)

Sturgeon (order Acipenserformes) provide an ideal taxonomic context for examination of genome duplication events. Multiple levels of ploidy exist among these fish. In a novel microsatellite approach, data from 962 fish from 20 sturgeon species were used for analysis of ploidy in sturgeon. Allele numbers in a sample of individuals were assessed at six microsatellite loci. Species with approximately 120 chromosomes are classified as functional diploid species, species with approximately 250 chromosomes as functional tetraploid species, and with approximately 500 chromosomes as functional octaploids. A molecular phylogeny of the sturgeon was determined on the basis of sequences of the entire mitochondrial cytochrome b gene. By mapping the estimated levels of ploidy on this proposed phylogeny we demonstrate that (I) polyploidization events independently occurred in the acipenseriform radiation; (II) the process of functional genome reduction is nearly finished in species with approximately 120 chromosomes and more active in species with approximately 250 chromosomes and approximately 500 chromosomes; and (III) species with approximately 250 and approximately 500 chromosomes arose more recently than those with approximately 120 chromosomes. These results suggest that gene silencing, chromosomal rearrangements, and transposition events played an important role in the acipenseriform genome formation. Furthermore, this phylogeny is broadly consistent with previous hypotheses but reveals a highly supported oceanic (Atlantic-Pacific) subdivision within the Acipenser/Huso complex.

Genetic evidence for limited dispersal in the coastal California killifish, Fundulus parvipinnis

The California killifish, Fundulus parvipinnis, is a marine species that lives in salt marshes, estuaries and wetlands along the California and Baja California coasts. In order to estimate levels of dispersal between different coastal habitats over its range, we have studied six populations using morphological and genetic markers. Lateral line scale and vertebrae counts showed significant differences between individuals collected north of Punta Eugenia and south of Punta Eugenia. Morphological differences across Punta Eugenia were accompanied by large genetic differences at the mitochondrial control region (5.8%). Gene flow was in general very reduced over the range of the species (pairwise average F(st)=0.70, Nm=0.30), with a strong break at Punta Eugenia (F(st)=0.95, Nm=0.03). Such limited interchanges between coastal habitats have important theoretical and conservation implications.

Heteroplasmy in the mtDNA control region of sturgeon (Acipenser, Huso and Scaphirhynchus)

Data from 1238 fishes from 19 sturgeon species and 1 paddlefish were used to analyze heteroplasmy in sturgeon. Lengths of central repeat units ranged from 74 to 83 bp among sturgeon species. No repeat sequence was found in the paddlefish, Polyodon spathula. A general feature of the repeat units was the presence of termination associated sequence (TAS) motifs. About 50% of 138 interspecific mutations observed among the D-loop sequences are located 10 bp down- and upstream from these TAS motifs. Interestingly, most homoplasmic species showed deletions upstream to the TAS motifs, whereas deletions downstream to the TAS motifs observed in two species do not seem to preclude heteroplasmy. Calculations of secondary structures and thermal stabilities of repeat units showed DeltaG values for all heteroplasmic species to be <-8 and for most homoplasmic species DeltaG value to be >-8. Most heteroplasmic fishes had two and/or three repeat units. No homoplasmic sturgeon with >2 repeat units were observed. Molecular phylogeny based on the entire cytochrome b showed that heteroplasmy probably resulted from a single evolutionary event. Our data demonstrate that heteroplasmy is present in most sturgeon species and suggest that the thermal stability of the secondary structure of the repeat unit in combination with mutations downstream of the TAS sequences influences heteroplasmy.

Mitochondrial genome organization and vertebrate phylogenetics

With the advent of DNA sequencing techniques the organization of the vertebrate mitochondrial genome shows variation between higher taxonomic levels. The most conserved gene order is found in placental mammals, turtles, fishes, some lizards and Xenopus. Birds, other species of lizards, crocodilians, marsupial mammals, snakes, tuatara, lamprey, and some other amphibians and one species of fish have gene orders that are less conserved. The most probable mechanism for new gene rearrangements seems to be tandem duplication and multiple deletion events, always associated with tRNA sequences. Some new rearrangements seem to be typical of monophyletic groups and the use of data from these groups may be useful for answering phylogenetic questions involving vertebrate higher taxonomic levels. Other features such as the secondary structure of tRNA, and the start and stop codons of protein-coding genes may also be useful in comparisons of vertebrate mitochondrial genomes.

Identification of Nile perch (Lates niloticus), grouper (Epinephelus guaza), and wreck fish (Polyprion americanus) by polymerase chain reaction-restriction fragment length polymorphism of a 12S rRNA gene fragment

Restriction site analysis of polymerase chain reaction (PCR) products from a conserved region of the 12S rRNA gene has been used for the specific identification of Nile perch (Lates niloticus), grouper (Epinephelus guaza), and wreck fish (Polyprion americanus). Amplification of DNA isolated from muscle samples was carried out using a set of primers flanking a region of 436 bp from the mitochondrial 12S rRNA gene. Digestions of the PCR products with RsaI and Sau96I endonucleases, followed by agarose gel electrophoresis of the digested PCR products, yielded specific profiles that enabled direct identification of each species analyzed.

Phylogenetic relationships of the North American sturgeons (order Acipenseriformes) based on mitochondrial DNA sequences

The evolutionary relationships of the extant species within the order Acipenseriformes are not well understood. Nucleotide sequences of four mitochondrial genes (12S rRNA, COII, tRNA(Phe), and tRNA(Asp) genes) in North American sturgeon and paddlefish were examined to reconstruct a phylogeny. Analysis of the combined gene sequences suggests a basal placement of the paddlefish with regard to the sturgeons. Nucleotide sequences of all four genes for the three Scaphirhynchus species were identical. The position of Scaphirhynchus based on our data was uncertain. Within the genus Acipenser, the two Acipenser oxyrinchus subspecies were very similar in sequence and found to be basal to the remaining Acipenser species examined. Based on our data, Acipenser transmontanus and Acipenser medirostris were sister taxa, as were Acipenser fulvescens and Acipenser brevirostrum. Comparison of our results with hypotheses of sturgeon relationships proposed by previous authors is presented. The sequence data presented here are phylogenetically useful and provide a solid foundation of genetic information for the North American Acipenseriformes that can be expanded to include Eurasian species to provide a global picture of sturgeon evolution.

Slipped-strand mispairing at noncontiguous repeats in Poecilia reticulata: a model for minisatellite birth

The standard slipped-strand mispairing (SSM) model for the formation of variable number tandem repeats (VNTRs) proposes that a few tandem repeats, produced by chance mutations, provide the "raw material" for VNTR expansion. However, this model is unlikely to explain the formation of VNTRs with long motifs (e.g., minisatellites), because the likelihood of a tandem repeat forming by chance decreases rapidly as the length of the repeat motif increases. Phylogenetic reconstruction of the birth of a mitochondrial (mt) DNA minisatellite in guppies suggests that VNTRs with long motifs can form as a consequence of SSM at noncontiguous repeats. VNTRs formed in this manner have motifs longer than the noncontiguous repeat originally formed by chance and are flanked by one unit of the original, noncontiguous repeat. SSM at noncontiguous repeats can therefore explain the birth of VNTRs with long motifs and the "imperfect" or "short direct" repeats frequently observed adjacent to both mtDNA and nuclear VNTRs.

Restricted gene flow and incipient speciation in disjunct Pacific Ocean and Sea of Cortez populations of a reef fish species, Girella nigricans

Population disjunctions have been proposed to play an important role in speciation processes. In this study, we have examined the possible role of the Pacific Ocean-Sea of Cortez disjunction as a contributing factor to cryptic speciation in a reef fish, the opaleye, Girella nigricans. Mitochondrial control region (D-loop) sequences (380 bp) of 117 individuals completely separated opaleye populations from the Pacific Ocean and the Sea of Cortez. Although opaleye exhibit pelagic larval stages that remain in the water column for several months, gene flow between the Pacific Ocean and the Sea of Cortez was found to be extremely limited (FST = 0.84, Nm = 0.10). Whereas limited gene flow and reciprocal monophyly suggest that the observed physical and genetic disjunction are potentially contributing to the incipient speciation of Pacific and Sea of Cortez opaleye, moderate levels of D-loop sequence divergence (3.3%) and the absence of fixed allozyme markers challenge this idea. Pacific Coast populations also exhibited restricted gene flow levels (FST = 0.25, Nm = 1.49) across Punta Eugenia, a recognized oceanographic boundary along the Baja California coast. Thus, opaleye individuals grouped into three clades: one clade in the Sea of Cortez, one Pacific clade south of Punta Eugenia, and one Pacific clade north of Punta Eugenia. Future work in this region will determine if our results can be generalized to other disjunct populations.

PCR-RFLP analysis of mitochondrial DNA: a reliable method for species identification

A method for identification of game species has been developed on the basis of the amplification of a specific part of the mitochondrial genome (tRNA(Glu)/cytochrome b) using the polymerase chain reaction (PCR). To distinguish between several game species, the obtained 464-bp-long PCR products were cut with different restriction endonucleases (RE) resulting in species-specific restriction fragment length polymorphism (RFLP). Even closely related deer species could be distinguished by application of one or two RE. Natural polymorphisms of the target sequence within one species were examined for red deer (Cervus elaphus), and base pair substitutions were identified affecting the RFLP pattern.

Organization of the large mitochondrial genome in the isopod Armadillidium vulgare

The mitochondrial DNA (mtDNA) in animals is generally a circular molecule of approximately 15 kb, but there are many exceptions such as linear molecules and larger ones. RFLP studies indicated that the mtDNA in the terrestrial isopod Armadillidium vulgare varied from 20 to 42 kb. This variation depended on the restriction enzyme used, and on the restriction profile generated by a given enzyme. The DNA fragments had characteristic electrophoretic behaviors. Digestions with two endonucleases always generated fewer fragments than expected; denaturation of restriction profiles reduced the size of two bands by half; densitometry indicated that a number of small fragments were present in stoichiometry, which has approximately twice the expected concentration. Finally, hybridization to a 550-bp 16S rDNA probe often revealed two copies of this gene. These results cannot be due to the genetic rearrangements generally invoked to explain large mtDNA. We propose that the large A. vulgare mtDNA is produced by the tripling of a 14-kb monomer with a singular rearrangement: one monomer is linear and the other two form a circular dimer. Densitometry suggested that these two molecular structures were present in different proportions within a single individual. The absence of mutations within the dimers also suggests that replication occurs during the monomer phase.

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

DNA sequences from the mitochondrial DNA control region are used to test the phylogeographic relationships among the pike-perches, Stizostedion (Teleostei: Percidae) and to examine patterns of variation. Sequences reveal two types of variability: single nucleotide polymorphisms and 6 to 14 copies of 10- to 11-base-pair tandemly repeated sequences. Numbers of copies of the tandem repeats are found to evolve too rapidly to detect phylogenetic signal at any taxonomic level, even among populations. Sequence similarities of the tandem repeats among Stizostedion and other percids suggest concerted evolutionary processes. Predicted folding of the tandem repeats and their proximity to termination-associated sequences indicate that secondary structure mediates slipped-strand mispairing among the d-loop, heavy, and light strands. Neighbor-joining and maximum parsimony analyses of sequences indicate that the genus is divided into clades on the continents of North America and Eurasia. Calibrating genetic distances with divergence times supports the hypothesis that Stizostedion dispersed from Eurasia to North America across a North Pacific Beringial land bridge approximately 4 million years before present, near the beginning of the Pliocene Epoch. The North American S. vitreum and S. canadense appear separated by about 2.75 million years, and the Eurasian S. lucioperca and S. volgensis are diverged by about 1.8 million years, suggesting that speciation occurred during the late Pliocene Epoch.

Heteroplasmy, length and sequence variation in the mtDNA control regions of three percid fish species (Perca fluviatilis, Acerina cernua, Stizostedion lucioperca)

The nucleotide sequence of the control region and flanking tRNA genes of perch (Perca fluviatilis) mtDNA was determined. The organization of this region is similar to that of other vertebrates. A tandem array of 10-bp repeats, associated with length variation and heteroplasmy was observed in the 5' end. While the location of the array corresponds to that reported in other species, the length of the repeated unit is shorter than previously observed for tandem repeats in this region. The repeated sequence was highly similar to the Mt5 element which has been shown to specifically bind a putative D-loop DNA termination protein. Of 149 perch analyzed, 74% showed length variation heteroplasmy. Single-cell PCR on oocytes suggested that the high level of heteroplasmy is passively maintained by maternal transmission. The array was also observed in the two other percid species, ruffe (Acerina cernua) and zander (Stizostedion lucioperca). The array and the associated length variation heteroplasmy are therefore likely to be general features of percid mtDNAs. Among the perch repeats, the mutation pattern is consistent with unidirectional slippage, and statistical analyses supported the notion that the various haplotypes are associated with different levels of heteroplasmy. The variation in array length among and within species is ascribed to differences in predicted stability of secondary structures made between repeat units.

Molecular phylogeny of Acipenserinae

The family Acipenseridae consists of 25 extant sturgeon species (19 species of Acipenserinae and 6 species of Scaphirhynchinae). Together with two extant paddlefish species, Polyodon spathula and Psephurus gladius (Polyodontidae), it composes the order Acipenseriformes, the most numerous of all living "fossil" fishes. This paper presents results of sequencing of three regions of the cytochrome b gene (650 bp), and fragments of 12S (150 bp) and 16S (350 bp) rRNA genes, from all extant species of Acipenserinae (species of Acipenser and Huso) and Scaphirhynchus albus (Scaphirhynchinae). The phylogenetic tree obtained for combined data is the first comprehensive treatment of phylogeny within the Acipenserinae. Three general conclusions are inferred from the tree: (1) The pallid sturgeon, S. albus, is the sister-species of all species of Acipenser and Huso. (2) The two species of Huso are embedded within the genus Acipenser. It also appears that Huso is not a separate taxonomic unit. (3) There are at least three main clades within Acipenser: A. sturio-A. oxyrinchus, A. schrenckii-A. transmontanus, and all Ponto-Caspian species plus A. dabryanus and A. brevirostrum. There is congruence between ploidy and the branching patterns of the sturgeon species. A hypothetical evolutionary history of the Acipenseriformes based on the paleontological, geological, and molecular data is discussed.

Evolution of repeated sequence arrays in the D-loop region of bat mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, theta and pi, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.

Evolution of repeated sequence arrays in the D-loop region of bat mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, theta and pi, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.