Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Quantifying the Elevation of Mitochondrial DNA Evolutionary Substitution Rates Over Nuclear Rates in the Intertidal Copepod Tigriopus californicus

Mitochondrial DNA (mtDNA) genomes generally evolve rapidly in animals, but considerable variation in the rates of evolution of mtDNA occurs among taxa. Higher levels of mutation will tend to increase the amount of polymorphism, which should also scale with population size, but there are mixed signals from previous studies on the evolutionary outcomes of the interactions of these processes. The copepod Tigriopus californicus provides an interesting model in which to study the evolution of mtDNA because it has high levels of divergence among populations and there is the suggestion that this divergence could be involved in reproductive isolation. This species also appears to have an elevated mtDNA substitution rate, but previous studies did not provide an accurate measurement. This article examines the rate of mtDNA substitution versus nuclear substitution in T. californicus and finds that the mtDNA rate for synonymous sites averages 55-fold higher, a level that exceeds the rates found in most other invertebrates. Levels of polymorphism are also examined in both mtDNA and nuclear genes, and it is shown that the effective population size of mtDNA genes is much lower than that of nuclear genes. In addition, no correlation between polymorphism in mtDNA and nuclear genes is found across populations, which suggest factors other than demography may shape polymorphism in this species. The results from this study suggest that mtDNA is evolving at a very rapid rate in this copepod species, and this could increase the likelihood that mtDNA evolution is involved in the generation of reproductive isolation.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

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

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

Mining the mammalian genome for artiodactyl systematics

A total of 7,806 nucleotide positions derived from one mitochondrial and eight nuclear DNA segments were used to provide a robust phylogeny for members of the order Artiodactyla. Twenty-four artiodactyl and two cetacean species were included, and the horse (order Perissodactyla) was used as the outgroup. Limited rate heterogeneity was observed among the nuclear genes. The partition homogeneity tests indicated no conflicting signal among the nuclear genes fragments, so the sequence data were analyzed together and as separate loci. Analyses based on the individual nuclear DNA fragments and on 34 unique indels all produced phylogenies largely congruent with the topology from the combined data set. In sharp contrast to the nuclear DNA data, the mtDNA cytochrome b sequence data showed high levels of homoplasy, failed to produce a robust phylogeny, and were remarkably sensitive to taxon sampling. The nuclear DNA data clearly support the paraphyletic nature of the Artiodactyla. Additionally, the family Suidae is diphyletic, and the nonruminating pigs and peccaries (Suiformes) were the most basal cetartiodactyl group. The morphologically derived Ruminantia was always monophyletic; within this group, all taxa with paired bony structures on their skulls clustered together. The nuclear DNA data suggest that the Antilocaprinae account for a unique evolutionary lineage, the Cervidae and Bovidae are sister taxa, and the Giraffidae are more primitive.     

Stochasticity overrules the "three-times rule": genetic drift,genetic draft,and coalescence times for nuclear loci versus mitochondrial DNA

Abstract.— Palumbi et al. (2001) proposed a "three-times rule" that uses mitochondrial DNA (mtDNA) sequences to predict probabilities of monophyly for nuclear loci (i.e., whether the alleles within a taxon coalesce with one another before they coalesce with alleles from a sister taxon). They use neutral coalescent theory to infer these probabilities from the ratio of interspecific divergence to intraspecific variation of mtDNA. We show that the estimated probabilities have very wide confidence intervals because of the inherent stochasticity of the mtDNA coalescent process. Under neutrality, the true probability of monophyly can be much higher, or much lower, than predicted by the three-times rule. We also review recent empirical and theoretical studies that refute neutrality-based predictions concerning mtDNA variation and divergence. We conclude that the three-times rule is neither a useful test for neutral molecular evolution nor a reliable guide to genealogical species.     

Rates of nucleotide substitution in Cornaceae (Cornales)-Pattern of variation and underlying causal factors

Identifying causes of genetic divergence is a central goal in evolutionary biology. Although rates of nucleotide substitution vary among taxa and among genes, the causes of this variation tend to be poorly understood. In the present study, we examined the rate and pattern of molecular evolution for five DNA regions over a phylogeny of Cornus, the single genus of Cornaceae. To identify evolutionary mechanisms underlying the molecular variation, we employed Bayesian methods to estimate divergence times and to infer how absolute rates of synonymous and nonsynonymous substitutions and their ratios change over time. We found that the rates vary among genes, lineages, and through time, and differences in mutation rates, selection type and intensity, and possibly genetic drift all contributed to the variation of substitution rates observed among the major lineages of Cornus. We applied independent contrast analysis to explore whether speciation rates are linked to rates of molecular evolution. The results showed no relationships for individual genes, but suggested a possible localized link between species richness and rate of nonsynonymous nucleotide substitution for the combined cpDNA regions. Furthermore, we detected a positive correlation between rates of molecular evolution and morphological change in Cornus. This was particularly pronounced in the dwarf dogwood lineage, in which genome-wide acceleration in both molecular and morphological evolution has likely occurred.     

Phylogenetic relationships of five pika species (genus Ochotona) based on mitochondrial DNA restriction maps

Restriction site mapping of mitochondrial DNA (mtDNA) with 16 restriction endonucleases was used to examine the phylogenetic relationships of Ochotona cansus, O. huangensis, O. thibetana, O. curzoniae and O. erythrotis. A 1-kb length variation between 0. erythrotis of subgenus Pika and other four species of subgenus Ochotona was observed, which may be a useful genetic marker for identifying the two subgenera. The phylogenetic tree constructed using PAUP based on 61 phylogenetically informative sites suggests that O. aythrotis diverged first, followed by O. cansus, while O. atrzoniae and O. huangensis are sister taxa related to O. thibetana. The results indicate that both O. cansus and O. huangensis should be treated as independent species. If the base substitution rate of pikas mtDNA was 2% per million years, then the divergence time of the two subgenera, Pika and Ochotona, is about 8.8 Ma ago of late Miocence, middle Bao-dian of Chinese mammalian age, and the divergence of the four species in subgenus     

Test for simultaneous divergence using approximate Bayesian computation

Comparative phylogeographic studies often reveal disparate levels of sequence divergence between lineages spanning a common geographic barrier, leading to the conclusion that isolation was nonsynchronous. However, only rarely do researchers account for the expected variance associated with ancestral coalescence and among-taxon variation in demographic history. We introduce a flexible approximate Bayesian computational (ABC) framework that can test for simultaneous divergence (TSD) using a hierarchical model that incorporates idiosyncratic differences in demographic history across taxon pairs. The method is tested across a range of conditions and is shown to be accurate even with single-locus mitochondrial DNA (mtDNA) data. We apply this method to a landmark dataset of putative simultaneous vicariance, eight geminate echinoid taxon pairs thought to have been split by the Isthmus of Panama 3.1 million years ago. The ABC posterior estimates are not consistent with a history of simultaneous vicariance given these data. Subsequent ABC estimates under a constrained model that assumes two divergence times across the eight taxon pairs suggests simultaneous divergence 3.1 million years ago in seven of the taxon pairs and a more recent divergence in the remaining taxon pair. These ABC estimates on the simultaneous divergence of the seven taxon pairs correspond to a DNA substitution rate of approximately 1.59% per lineage per million years at the mtDNA cytochrome oxidase I gene. This ABC framework can easily be modified to analyze single taxon-pair datasets and/or be expanded to include multiple loci, migration, recombination, and other idiosyncratic demographic histories. The flexible aspect of ABC and its built-in evaluation of estimator bias and statistical power has the potential to greatly enhance statistical rigor in phylogeographic studies.     

Decreased diversity but increased substitution rate in host mtDNA as a consequence of Wolbachia endosymbiont infection

A substantial fraction of insects and other terrestrial arthropods are infected with parasitic, maternally transmitted endosymbiotic bacteria that manipulate host reproduction. In addition to imposing direct selection on the host to resist these effects, endosymbionts may also have indirect effects on the evolution of the mtDNA with which they are cotransmitted. Patterns of mtDNA diversity and evolution were examined in Drosophila recens, which is infected with the endosymbiont Wolbachia, and its uninfected sister species D. subquinaria. The level of mitochondrial, but not nuclear, DNA diversity is much lower in D. recens than in D. subquinaria, consistent with the hypothesized diversity-purging effects of an evolutionarily recent Wolbachia sweep. The d(N)/d(S) ratio in mtDNA is significantly greater in D. recens, suggesting that Muller's ratchet has brought about an increased rate of substitution of slightly deleterious mutations. The data also reveal elevated rates of synonymous substitutions in D. recens, suggesting that these sites may experience weak selection. These findings show that maternally transmitted endosymbionts can severely depress levels of mtDNA diversity within an infected host species, while accelerating the rate of divergence among mtDNA lineages in different species.     

GEOGRAPHIC DISTRIBUTION OF MOLECULAR VARIANCE WITHIN THE BLUE MARLIN (MAKAIRA NIGRICANS): A HIERARCHICAL ANALYSIS OF ALLOZYME,SINGLE-COPY NUCLEAR DNA,AND MITOCHONDRIAL DNA MARKERS

This study presents a comparative hierarchical analysis of variance applied to three classes of molecular markers within the blue marlin (Makaira nigricans). Results are reported from analyses of four polymorphic allozyme loci, four polymorphic anonymously chosen single-copy nuclear DNA (scnDNA) loci, and previously reported restriction fragment length polymorphisms (RFLPs) of mitochondrial DNA (mtDNA). Samples were collected within and among the Atlantic and Pacific Oceans over a period of several years. Although moderate levels of genetic variation were detected at both polymorphic allozyme (H = 0.30) and scnDNA loci (H = 0.37), mtDNA markers were much more diverse (h = 0.85). Allele frequencies were significantly different between Atlantic and Pacific Ocean samples at three of four allozyme loci and three of four scnDNA loci. Estimates of allozyme genetic differentiation (θ_O) ranged from 0.00 to 0.15, with a mean of 0.08. The θ_O values for scnDNA loci were similar to those of allozymes, ranging from 0.00 to 0.12 with a mean of 0.09. MtDNA RFLP divergence between oceans (θ_O = 0.39) was significantly greater than divergence detected at nuclear loci (95% nuclear confidence interval = 0.04–0.11). The fourfold smaller effective population size of mtDNA and male-mediated gene flow may account for the difference observed between nuclear and mitochondrial divergence estimates.     

A new timeframe for the diversification of Japan’s mammals

Aim The main Japanese islands are land‐bridge islands divided by the biogeographic division Blakiston’s Line and represent two natural laboratories for studying land‐bridge diversification. Colonization of the current mammal fauna has been dated to the middle to late Pleistocene using fossil evidence. The purpose of this paper is to apply a molecular clock to the genetic divergences between Japanese mammalian taxa and their sister mainland taxa to test the late Pleistocene land‐bridge colonization hypothesis. Location The main Japanese islands (Kyushu, Shikoku, Honshu and Hokkaido). Methods I used mitochondrial DNA (cytochrome b) and a species tree approach to estimate the divergence times of 24 Japanese non‐volant terrestrial mammal taxa and their mainland sister taxa using the program *beast . I then tested for evidence of non‐simultaneous divergence among these taxon‐pairs by controlling for expected coalescent stochasticity using the program Ms Bayes . Results Divergence events between taxa on Japan and their mainland sister taxa were significantly older than expected under the current paradigm, which is based on fossil data. Consistent with the land‐bridge colonization hypothesis, there was evidence of multiple divergence events. Main conclusions These results implicate a colonization timeframe that is older than posited by the current paradigm based on fossil evidence. However, these results are still consistent with the land‐bridge colonization hypothesis. Multiple periods of land‐bridge connectivity may account for the current mammalian fauna in Japan. In addition, half of the divergence time estimates in the Honshu–Shikoku–Kyushu region were clumped around 2.4 Ma, which might suggest a dramatic interchange period, concordant with a period of significant global cooling, when the first land bridge may have connected Japan to the mainland.     

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

ABSTRACT: BACKGROUND: The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs). RESULTS: To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are reciprocally monophyletic in the mtgenome phylogeny but a H. pileatus individual falls into the H. lar Y chromosome clade. CONCLUSIONS: Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.     

Correlated evolutionary divergence of egg size and a mitochondrial protein across the Isthmus of Panama

An explicit assumption of studies that employ a mitochondrial DNA (mtDNA) molecular clock is that mtDNA evolves independently of morphology. Here we report a very strong correlation between egg size divergence and cytochrome c oxidase-1 (CO1) amino acid sequence divergence among sister species of bivalve molluscs separated by the Central American Isthmus (i.e., "geminate" species). Analyses of the molecular data reveal that CO1 sequences likely did not diverge as a function of time or evolve in response to positive natural selection. Given that an excess of CO1 amino acid polymorphism exists within species (as expected if most mutations are only slightly deleterious), a third hypothesis is that reductions in effective population size could simultaneously increase the fixation rate of nearly neutral mtDNA polymorphisms and in some way also facilitate egg size evolution. The remarkable strength of the relationship between egg size and CO1 amino acid sequence demonstrates that, even in the absence of an obvious functional relationship or clock-like evolution, the amounts of molecular and morphological change can be tightly correlated, and therefore may reflect common processes. Accordingly, the assumption that the evolutionary divergence of molecules and morphology are independent must always be carefully examined.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Patterns of polymorphism and gene flow of gender-associated mitochondrial DNA lineages in European mussel populations

Mussels of the genus Mytilus have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally and the F type is transmitted maternally. RFLP analysis is used to assess phylogenetic relationships and nucleotide diversity and divergence for both mtDNA genomes in European populations of M. edulis and Atlantic and Mediterranean forms of M. galloprovincialis. Ten restriction endonucleases were used to assay variation in regions of the ND2 and COIII genes for a total of 77 individuals. F and M genomes show a concordant phylogenetic split into two major divergent clades, one specific to Mediterranean M. galloprovincialis and the other containing haplotypes from the three taxa. For both genomes, the geographical distribution of mtDNA variation suggests: (i) extensive levels of mtDNA introgression; (ii) asymmetric mtDNA gene flow from Atlantic to Mediterranean populations; and (iii) recurrent historical hybridization events. Significantly higher mtDNA diversity and divergence are observed for the M than F genome in all three Mytilus taxa, although the evolutionary forces responsible for these differences cannot be resolved. The extensive mtDNA gene flow between European Mytilus taxa conflicts with the restricted mtDNA introgression observed in American mussels , implying geographical variation in the nature of nuclear/mtDNA interactions regulating biparental inheritance.     

DNA variability and recombination rates at X-linked loci in humans.

We sequenced 11,365 bp from introns of seven X-linked genes in 10 humans, one chimpanzee, and one orangutan to (i) provide an average estimate of nucleotide diversity (pi) in humans, (ii) investigate whether there is variation in pi among loci, (iii) compare ratios of polymorphism to divergence among loci, and (iv) provide a preliminary test of the hypothesis that heterozygosity is positively correlated with the local rate of recombination. The average value for pi was low 0.063%, SE = 0.036%, about one order of magnitude smaller than for Drosophila melanogaster, the species for which the best data are available. Among loci, pi varied by over one order of magnitude. Statistical tests of neutrality based on ratios of polymorphism to divergence or based on the frequency spectrum of variation within humans failed to reject a neutral, equilibrium model. However, there was a positive correlation between heterozygosity and rate of recombination, suggesting that the joint effects of selection and linkage are important in shaping patterns of nucleotide variation in humans.     

A comparative summary of genetic distances in the vertebrates from the mitochondrial cytochrome b gene

Mitochondrial cytochrome b (cytb) is among the most extensively sequenced genes to date across the vertebrates. Here, we employ nearly 2,000 cytb gene sequences from GenBank to calculate and compare levels of genetic distance between sister species, congeneric species, and confamilial genera within and across the major vertebrate taxonomic classes. The results of these analyses parallel and reinforce some of the principal trends in genetic distance estimates previously reported in a summary of the multilocus allozyme literature. In particular, surveyed avian taxa on average show significantly less genetic divergence than do same-rank taxa surveyed in other vertebrate groups, notably amphibians and reptiles. Various biological possibilities and taxonomic "artifacts" are considered that might account for this pattern. Regardless of the explanation, by the yardstick of genetic divergence in this mtDNA gene, as well as genetic distances in allozymes, there is rather poor equivalency of taxonomic rank across some of the vertebrates.     

No variation and low synonymous substitution rates in coral mtDNA despite high nuclear variation

Background  

The mitochondrial DNA (mtDNA) of most animals evolves more rapidly than nuclear DNA, and often shows higher levels of intraspecific polymorphism and population subdivision. The mtDNA of anthozoans (corals, sea fans, and their kin), by contrast, appears to evolve slowly. Slow mtDNA evolution has been reported for several anthozoans, however this slow pace has been difficult to put in phylogenetic context without parallel surveys of nuclear variation or calibrated rates of synonymous substitution that could permit quantitative rate comparisons across taxa. Here, I survey variation in the coding region of a mitochondrial gene from a coral species (Balanophyllia elegans) known to possess high levels of nuclear gene variation, and estimate synonymous rates of mtDNA substitution by comparison to another coral (Tubastrea coccinea).