Molecular Biology and Evolution

MBE Advance Access originally published online on June 29, 2005
Molecular Biology and Evolution 2005 22(10):2104-2112; doi:10.1093/molbev/msi204

This Article Abstract FREE Full Text (PDF) All Versions of this Article: 22/10/2104    most recent msi204v1 Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (31) Request Permissions Google Scholar Articles by Mueller, R. L. Articles by Boore, J. L. Search for Related Content PubMed PubMed Citation Articles by Mueller, R. L. Articles by Boore, J. L. Social Bookmarking          What's this?

© The Author 2005. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oupjournals.org

Research Article

Molecular Mechanisms of Extensive Mitochondrial Gene Rearrangement in Plethodontid Salamanders

Rachel Lockridge Mueller^*,,1 and Jeffrey L. Boore^,

^* Museum of Vertebrate Zoology, University of California, Berkeley; Evolutionary Genomics Department, Department of Energy Joint Genome Institute and University of California Lawrence Berkeley National Laboratory, Walnut Creek, California; and Department of Integrative Biology, University of California, Berkeley

E-mail: rmueller{at}uchicago.edu.

	Abstract

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 
Extensive gene rearrangement is reported in the mitochondrial genomes of lungless salamanders (Plethodontidae). In each genome with a novel gene order, there is evidence that the rearrangement was mediated by duplication of part of the mitochondrial genome, including the presence of both pseudogenes and additional, presumably functional, copies of duplicated genes. All rearrangement-mediating duplications include either the origin of light-strand replication and the nearby tRNA genes or the regions flanking the origin of heavy-strand replication. The latter regions comprise nad6, trnE, cob, trnT, an intergenic spacer between trnT and trnP and, in some genomes, trnP, the control region, trnF, rrnS, trnV, rrnL, trnL1, and nad1. In some cases, two copies of duplicated genes, presumptive regulatory regions, and/or sequences with no assignable function have been retained in the genome following the initial duplication; in other genomes, only one of the duplicated copies has been retained. Both tandem and nontandem duplications are present in these genomes, suggesting different duplication mechanisms. In some of these mitochondrial DNAs, up to 25% of the total length is composed of tandem duplications of noncoding sequence that includes putative regulatory regions and/or pseudogenes of tRNAs and protein-coding genes along with the otherwise unassignable sequences. These data indicate that imprecise initiation and termination of replication, slipped-strand mispairing, and intramolecular recombination may all have played a role in generating repeats during the evolutionary history of plethodontid mitochondrial genomes.

Key Words: complete mitochondrial genomes • gene rearrangement • tandem duplication • pseudogenes • plethodontids

	Introduction

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 
In contrast to the conclusions drawn from early and limited sampling, animal mitochondrial genomes possess unexpected diversity both in gene order and in the presence, extent, and distribution of noncoding sequence (Inoue et al. 2003; Boore, Medina, and Rosenberg 2004; Miller et al. 2004; Yokobori et al. 2004). Groups of organisms with highly differing gene orders are appropriate model systems for studying the factors that effect mitochondrial gene rearrangement (Dowton and Campbell 2001). Such groups have been identified among invertebrates (Yamazaki et al. 1997; Dowton and Austin 1999; Shao et al. 2001), and testable hypotheses for the causes of frequent gene rearrangement have been generated; for example, highly rearranged parasitic wasp mitochondrial genomes may result from oxidative stress imposed by the host immune response (Dowton and Campbell 2001). Within the largest salamander family, Plethodontidae, 6 of the 22 mitochondrial genomes examined have rearrangements of independent phylogenetic origin, and 12 of these 22, plus 1 from the family Rhyacotritonidae, have tandem repeats in noncoding sequences. This high instance of gene rearrangement makes plethodontid salamanders an excellent model system for examining the mechanisms of such vertebrate mitochondrial genome instability.

In the most commonly invoked model of mitochondrial gene rearrangement, a region of the genome is duplicated and the supernumerary genes, no longer maintained by selection, are eliminated; the original gene arrangement is either restored or altered, depending on the pattern of gene loss (Moritz, Dowling, and Brown 1987; Boore 2000). Genomic evidence for this model of rearrangement includes (1) repeated motifs, which can cause duplication via slipped-strand mispairing; (2) stem-loop structures, which can cause duplication via intramolecular recombination (Stanton et al. 1994; Moore, Gudikote, and Van Tuyle 1998); and (3) pseudogenes, which may be ancestral duplications in the process of being eliminated (Arndt and Smith 1998; Macey et al. 1998). This mechanism cannot explain gene inversions, which have also been reported (Smith et al. 1989; Boore 1999; Dowton et al. 2003; Miller et al. 2004). In each of the six rearranged genomes reported here, evidence of at least one historical duplication event has been retained. Plethodontid rearrangements involve regions of the genome that are independently rearranged or duplicated in other vertebrate taxa (Moritz and Brown 1986; Desjardins and Morais 1990; Moritz 1991; Pääbo et al. 1991; Quinn and Mindell 1996). The molecular mechanisms of gene rearrangement at work in plethodontids may therefore be the same mechanisms acting across vertebrates.

Organismal-, cellular-, and genomic-level properties may be linked to mitochondrial genome instability, both in these salamanders and more generally. In this paper, we (1) describe mitochondrial genome rearrangements and expansions of noncoding DNA in plethodontid and related salamanders, (2) infer possible molecular mechanisms by which these rearrangements and expansions originated, and (3) outline possible explanations for the extensive mitochondrial genome instability in this clade.

	Materials and Methods

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 
Twenty-four complete mitochondrial genomes of plethodontid and related salamanders were sequenced, assembled, and annotated as described elsewhere (Mueller et al. 2004) (GenBank accession numbers AY728212–AY728235). Annotated genomes were examined for novel gene orders and the presence of unassignable sequences. Putative pseudogene sequences in each rearranged genome were identified by position and aligned with the corresponding functional gene sequences using ClustalW (gap parameters set to default: open = 10, extend = 5) and manual adjustment. Two other noncoding regions of the genome were also examined: the control region (CR hereafter), and an intergenic spacer between trnT and trnP (IGS hereafter) present in diverse salamander clades that may be a remnant of at least one older duplication (McKnight and Shaffer 1997; Zardoya et al. 2003; Zhang et al. 2003a,b). Because these regions are highly variable with only a few, or no, short conserved sequences (Shadel and Clayton 1997), we defined these two regions by their genome position. All noncoding regions were tested for the presence of tandem repeat elements and for shared repeats among the different regions by the construction of Pustell DNA-DNA matrices using MacVector (window size = 30, minimum percent score = 80, hash = 6) (Accelrys, San Diego, Calif.). The genome sequence of Hydromantes italicus is incomplete (trnT, the IGS, trnP, the CR, and trnF are missing) and therefore was not examined for repeats. Finally, each genome was tested for possible heteroplasmy in the numbers of repeated regions by examining the sequences of individual clones for their ability to produce assemblies with identical end sequences but different numbers of repeats from the primary assembly.

	Results and Discussion

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 
Duplications Have Mediated Most Rearrangements
Six of the 22 plethodontid genomes have novel gene orders and arrangements of noncoding sequences that are consistent with the duplication-random loss model of gene rearrangement (Boore 2000). In this model, a portion of the genome containing at least two genes is duplicated. One of the two copies of each duplicated gene eventually loses function, becomes a pseudogene, and is excised from the genome by mutational processes because it is not maintained by selection. Which gene copy is ultimately lost is determined by the first loss-of-function mutation, with some patterns restoring the original order and others leading to rearrangement. In this study, two such rearrangements include the origin of light-strand replication (O_L), two include the IGS, and two include both the origin of heavy-strand replication (O_H, contained within the CR) and the IGS. In some cases, there are evident pseudogenes that signal historical duplications. In others, these pseudogenes have decayed beyond 60% identity to the functional copy and are inferred based solely on their genome position. In all cases, the exact boundaries of the duplicated fragments have likely been obscured by deletion, and the extent of the original duplication may therefore be underestimated. The phylogenetic positions of the taxa with these rearranged genomes indicate separate rearrangement events, with the possible exception of one rearrangement shared by Aneides flavipunctatus and Aneides hardii, although the genome of A. hardii also appears to have undergone a second duplication-mediated rearrangement (fig. 1). For each rearranged mitochondrial genome, we infer the minimum set of contiguous genes that, when duplicated, creates an intermediate sequence from which gene losses could have established the observed gene arrangement.

View larger version (45K): [in this window] [in a new window]   FIG. 1.— Phylogenetic relationships among plethodontid salamanders and two out-groups from other families as indicated (from Mueller et al. 2004). Asterisks indicate lineages that have experienced a duplication-mediated rearrangement. Aneides hardii and Aneides flavipunctatus may share one synapomorphic rearrangement; however, A. hardii underwent a second duplication-mediated rearrangement, not present in A. flavipunctatus. Numbers to the right of species names are mitochondrial genome sizes. The sequences of two species' genomes, Nototriton abscondens and Hydromantes italicus, are incomplete.

 
Rearrangements Including the O_L
Batrachoseps attenuatus
The hypothesized duplicated region in this species' mitochondrial DNA (mtDNA) includes a small fragment of trnW and all of trnA, trnN, O_L, trnC, and trnY (300 bp in total) (fig. 2). For the purposes of these analyses, O_L refers to the entire noncoding region normally found between trnN and trnC that has the potential to form a large stem-loop structure known to function as an origin of replication in some systems. We refer to a copy of O_L by pseudogene annotation if it is greatly deficient in potential for forming this structure relative to the other copy, although no experimental data are available to infer which, if either, actually functions as an origin of replication. Deletions in the pseudogenes have reduced their sizes to the following lengths and percentages of original length: partial trnW + trnA, 77 bp (the percentage of original length cannot be determined because the amount of trnW included in the initial duplication is unknown); trnN, 65 bp (93%); O_L, 8 bp (27%); trnC, 56 bp (86%); and trnY, 28 bp (amount of trnY duplicated is unknown).

View larger version (32K): [in this window] [in a new window]   FIG. 2.— Mitochondrial gene rearrangements in Batrachoseps attenuatus. OL = origin of light-strand replication. Single letters refer to tRNA genes for the corresponding amino acid. Shaded boxes represent recognizable pseudogenes and question marks indicate stretches of sequence that cannot be assigned based on sequence similarity. Lengths of genes and pseudogenes are not to scale. (A) Current gene order. (B) Hypothesized duplication-random loss model for deriving this gene order. trnN retains 60% identity to the functional copy overall, not including two deletions (1 and 3 bp) in the pseudogene; the last 39 bp retain 74% to the functional copy. The remaining pseudogenes have decayed beyond recognition.

 
Hydromantes brunus
The hypothesized duplicated region in this species' mtDNA is slightly larger than in B. attenuatus, encompassing the end of nad2 and all of trnW in addition to trnA, trnN, O_L, trnC, and trnY (525 bp) (fig. 3). Deletions in the pseudogenes have reduced their sizes to the following lengths and percentages of original length: partial nad2, 127 bp (amount of nad2 duplicated is unknown); trnW, 60 bp (88%); trnA, 55 bp (81%); trnN, 65 bp (94%); O_L, 27 bp (73%); trnC, 57 bp (86%); and trnY, 39 bp (58%). Although the initial duplications of B. attenuatus and H. brunus involved similar regions of the genome, different copies of the redundant genes were subsequently lost; this is consistent with random loss, in which the copy of the gene that sustains the first disabling mutation continues to decay.

View larger version (37K): [in this window] [in a new window]   FIG. 3.— Mitochondrial gene rearrangements in Hydromantes brunus. Designations are as in figure 2. (A) Current gene order. (B) Hypothesized duplication-random loss model for deriving this gene order. trnW retains 69% identity to the functional copy, not including two deletions (3 and 5 bp) in the pseudogene. trnA retains 56% identity to the functional copy overall, not including two deletions (3 and 5 bp) in the pseudogene; the last 44 bp retain 66% identity. trnN retains 82% identity to the functional copy, not including one 5-bp deletion in the pseudogene. OL retains 70% identity to the functional copy overall, not including 10 bp missing from its end; the first 18 bp retain 94% identity. trnC retains 87% identity to the functional copy, not including three deletions (2, 1, and 2 bp) in the pseudogene. trnY and nad2 have decayed beyond recognition.

 
Rearrangements Including the CR and/or the IGS
Stereochilus marginatus
The hypothesized duplication in this species' mtDNA includes nad6, trnE, cob, trnT, and the IGS for a total of 1,790 bp plus the unknown length of the IGS (fig. 4). Deletions in the pseudogenes have reduced their sizes to the following lengths and percentages of original lengths: nad6, 186 bp (36%); trnE, 66 bp (98.5%); and cob, 42 bp (3.7%). trnT, if it still exists, cannot be identified. Two different sets of repeats exist in the two regions of the genome inferred to be copies of the IGS.

View larger version (38K): [in this window] [in a new window]   FIG. 4.— Mitochondrial gene rearrangements in Stereochilus marginatus. Designations are as in figure 2 except that IGS = intergenic spacer between trnT and trnP and CR = control region. (A) Current gene order. (B) Hypothesized duplication-random loss model for deriving this gene order. nad6 retains 52% identity to the inferred corresponding 186-bp portion of the functional copy, not including two insertions (2 and 1 bp) in the pseudogene; base pairs 55–132 retain 74% identity. trnE retains 85% identity to the functional copy, not including one 1-bp deletion in the pseudogene. The 42-bp cob retains 79% identity to the first 44 bp of the functional copy, not including one 2-bp deletion in the pseudogene. The first putative copy of the IGS, between trnT and nad6, is 1,982 bp in length and contains three complete, and one partial, >95% identical tandem repeats of 437 bp each. The sequence between cob and trnP contains five >99% identical copies of a different 110-bp tandem repeat sequence and a sixth 87% identical copy. Notably, following an intervening 102-bp stretch of nonrepetitive sequence, there is a complete copy of the 437-bp repeat unit found in the putative IGS between trnT and nad6. This repeat unit is 85% identical to copies in the other putative IGS, despite being separated from them by 1,630 bp. No recognizable trnT remains.

 
Plethodon elongatus
This species' mitochondrial genome has a gene order and pattern of noncoding sequence consistent with two separate duplications. The initial hypothesized duplication spans the region of the genome including nad6, trnE, cob, trnT, the IGS, trnP, and the CR for a total of 2,811 bp plus the unknown preduplication length of the IGS (fig. 5). Different copies of the redundant genes subsequently decayed to pseudogenes in P. elongatus than in S. marginatus. Deletions in the inferred redundant genes have reduced their sizes to the following lengths and percentages of original lengths: nad6 + trnE + cob, 47 bp (2.7%); trnP, length unknown because the boundary between the IGS and trnP is unassignable; and trnT + IGS, 107 bp (original length of the IGS is unknown). Notably, two 959-bp 97% identical copies of the putative CR are retained in the genome, indicating that the two may be undergoing concerted evolution as has been reported for several other taxa (Arndt and Smith 1998; Kumazawa et al. 1998; Lee et al. 2001; Inoue et al. 2003).

View larger version (38K): [in this window] [in a new window]   FIG. 5.— Mitochondrial gene rearrangements in Plethodon elongatus. Designations are as in figure 4. (A) Current gene order. (B) Two hypothesized duplication events mediating rearrangement in the P. elongatus mitochondrial genome. The first duplication resulted in tandem repeats of the region spanning from nad6 to the CR. trnP retains 61% identity to the functional copy overall, not including two insertions (3 bp each) in the pseudogene; the last 44 bp retain 70% identity. nad6, trnE, cob, and trnT have all decayed beyond recognition. In contrast, the two copies of the CR are 97% identical. Later, a duplicate copy of the fragment comprising the end of nad5, nad6, trnE, cob, trnT, the IGS, trnP, and a portion of the CR was inserted between the second CR and trnF, likely by intramolecular recombination. The two copies of this fragment are 99% identical.

 
Following this initial rearrangement by duplication-random loss, we hypothesize that a second duplication gave rise to two copies of the region including the last 57 bp of nad5, nad6 + trnE + cob, trnT, the IGS, trnP, and the first 761 bp of one CR (1,150 bp in total). These two copies are not adjacent to one another; rather, they are separated by 3,054 bp. The two copies are >99% identical in sequence, implying either (1) very recent duplication or (2) concerted evolution, which is also operating to maintain the sequence identity of the two full-length CRs.

Aneides flavipunctatus
The hypothesized duplication in this species' mtDNA includes nad6, trnE, cob, trnT, the IGS, and trnP for a total of 1,860 bp plus the unknown length of the IGS (fig. 6). The same copies of the redundant genes subsequently decayed to pseudogenes in A. flavipunctatus as in S. marginatus. Deletions in nad6 and trnE have completely excised them from the genome. cob and trnT, if they still exist, cannot be identified. As seen in the S. marginatus genome, two different sets of repeats exist in the two regions of the genome inferred to contain copies of the IGS.

View larger version (33K): [in this window] [in a new window]   FIG. 6.— Mitochondrial gene rearrangements in Aneides flavipunctatus. Designations are as in figure 4. (A) Current gene order. (B) Hypothesized duplication-random loss model for deriving this gene order. nad6 and trnE have been excised from the genome. The region between trnT and nad6, which comprises a putative IGS and trnP, is 3,050 bp in length and contains seven complete, and one partial, >96% identical 388-bp tandem repeats. Each tandem repeat contains a 72-bp trnP with 54% overall sequence identity to trnP, not including two insertions in the pseudogene (2 and 1 bp); the last 51 bp are 63% identical to the functional trnP. The region between trnE and trnP contains ten >92% identical copies of a 66-bp tandem repeat, a 357-bp stretch of variable numbers of short repeats (5–10 bp), and 71 bp of nonrepetitive sequence. No recognizable cob or trnT exists. The two repetitive regions of the genome resemble neither one another nor the CR.

 
Aneides hardii
This mitochondrial genome has a gene order and pattern of noncoding sequence consistent with two separate tandem duplications (fig. 7). The initial hypothesized duplication includes a near-identical region of the genome duplicated in A. flavipunctatus (from nad6 through the IGS), suggesting that this duplication may be a synapomorphy of the Aneides clade; however, unlike in A. flavipunctatus, no evidence remains that trnP was duplicated in A. hardii. The absence of an identifiable trnP in this duplicated region of the A. hardii genome suggests that (1) the trnP identified in A. flavipunctatus is false, (2) trnP has been excised or completely degraded from A. hardii, or (3) a different duplication-mediated rearrangement occurred independently in each lineage. The second hypothesized duplication in A. hardii includes nad6—the IGS, plus a second IGS, trnP, the CR, trnF, rrnS, trnV, rrnL, trnL1, and nad1 for a total of 6,227 bp and the unknown lengths of the two IGS regions. The copies of nad6 and trnE involved in the second round of duplication must have been functional at the time of reduplication, based on the current position of functional copies in the genome. The nad6 and trnE between nad5and cob, which were not involved in the second round of duplication, have been reduced to 11 bp (2%).

View larger version (26K): [in this window] [in a new window]   FIG. 7.— Mitochondrial gene rearrangements in Aneides hardii. Designations are as in figure 4. (A) Current gene order. (B) Two hypothesized duplication events in the history of the A. hardii mitochondrial genome. The first duplication comprised the region spanning from nad6 to the IGS; this duplication may be a synapomorphy shared by A. hardii and Aneides flavipunctatus. The second duplication comprised the region spanning from IGS to nad1. The putative IGS (IGSa) that bounds the two copies of this duplication fragment (Copy I and Copy II) contains varying numbers of a 343-bp tandem repeat. Copy I has one complete and one partial tandem repeat units, and Copy II has four complete >95% identical tandem repeat units and one partial tandem repeat unit. Adjacent to these IGSa sequences, Copy II contains a functional nad6, whereas Copy I possesses neither a functional nad6 nor a recognizable nad6. The two copies of trnE, adjacent to nad6 (Copy II) and IGSa (Copy I), are 97% identical to one another excluding a 3-bp deletion in the D-stem of Copy II. cob, trnT, and an additional IGS (IGSb), if retained in the genome, should be found in each fragment between trnE and trnP. There is no recognizable cob or trnT in this region in either fragment, nor is there any sequence similarity with the 343-bp repeat unit found in IGSa. Rather, this region of each fragment contains 67-bp tandem repeats, followed by 300 bp of nonrepetitive sequence. Copy I has six complete tandem repeats; Copy II has four complete repeats and a partial fifth repeat. Not including this indel, the two copies of this region are 96% identical. Both Copy I and Copy II contain trnP, CR, and trnF; the two copies of this region are 96% identical overall, and both copies of both tRNAs have viable secondary structures. Adjacent to trnF, Copy I has a complete rrnS, trnV, rrnL, and trnL1. In contrast, Copy II contains 111 bp that are 86% identical to the first 111 bp of the functional rrnS, followed by a trnL1 that is 97% identical to the trnL1 in Copy I; both trnL1s have viable secondary structure. Adjacent to trnL1, Copy I has a 672-bp nad1 that is 92% identical to the first 672 bp of the functional copy but contains multiple stop codons; Copy II contains the functional nad1.

 
The two large, duplicated fragments resulting from the second hypothesized round of duplication (Copy I and Copy II; see fig. 7) share high levels of sequence similarity over much of their length, although nad6, much of rrnS, and all of trnV and rrnL have been excised from the genome and nad1 has accumulated multiple stop codons. In contrast, two very similar, presumably functional, copies of trnP, trnF, trnL1, the CR, and possibly trnE have been retained in this genome. There are several possible explanations for this pattern. Consistent with the duplication-random loss model, this duplication event may be sufficiently recent that mutations have not yet disrupted trnP, trnF, trnL1, the CR, or trnE. Alternatively, these copies may have been maintained by selection or may have had their mutations corrected by gene conversion (Eberhard, Wright, and Bermingham 2001). Selection to retain duplicate copies of these genes or the CR may act either on the production of their gene products or their ability to form secondary structures necessary for message processing (Kumazawa et al. 1998) or on regulatory function, respectively. However, the presence of tandem repeats highly conserved between Copy I and Copy II is difficult to explain by selection in the absence of a clear function assignable to these sequences. Rather, it suggests either a recent duplication event or a purely mechanistic cause of concerted evolution between portions of these two fragments (Moritz and Brown 1987).

Possible Mechanisms Yielding Tandem Duplications
Several mechanisms may produce duplications in a circular molecule: recombination, slipped-strand mispairing, errors in synchronizing the points of initiation and termination of replication, or some combination of these processes (Macey et al. 1997a,b; Boore and Brown 1998; Boore 2000; Dowton and Campbell 2001). At least two such mechanisms are plausible explanations for the patterns observed in the order of genes and noncoding DNA in most plethodontid rearrangements. In the genomes of H. brunus, B. attenuatus, and A. hardii (second duplication), the replication origin is in the center of the duplicated region, indicating that both initiation and termination of replication occurred at alternative sites; alternate initiation or termination alone would duplicate only the region of the genome upstream or downstream, respectively, of the replication origin. A presumably functional O_L with viable potential secondary structure (Pääbo et al. 1991; Macey, Schulte, and Larson 2000) exists in both H. brunus and B. attenuatus. This suggests that light-strand replication initiation has not been completely transferred to alternate structures in these genomes and that alternate initiation may be causally linked to imprecise termination. In the genomes of S. marginatus, A. flavipunctatus, and A. hardii (first duplication), no evidence remains that a replication origin was involved in the duplication. The patterns in these genomes are still consistent with initiation and termination of replication at alternate sites, but both sites are located downstream of the O_H. The IGS, which is the putative alternative initiation site for S. marginatus and A. hardii, may be an ancient duplicated CR based on its position in the genome (McKnight and Shaffer 1997) and thus may retain limited initiation capability. Imprecise termination alone is also a plausible duplication mechanism for S. marginatus, A. flavipunctatus, and A. hardii based on the proximity of these duplicated regions to the CR and the possible erosion of the ends of the duplicated region. The first duplication in the P. elongatus genome is bounded by the CR; depending on the location of the O_H within the CR, imprecise termination alone, or imprecise termination with initiation at an alternate site, is a possible duplication mechanism. We cannot eliminate the alternative possibility that either intramolecular recombination or slipped-strand mispairing caused any of these duplications, nor can we rule them out as mechanisms for excision of redundant genes following these duplications (Holt, Dunbar, and Jacobs 1997; Lunt and Hyman 1997; Kajander et al. 2000; Tang et al. 2000a; Miller et al. 2004).

Possible Mechanisms Yielding Nontandem Duplications
The second duplication event in P. elongatus produced nontandem repeat fragments separated from one another by 3,054 bp, which nonetheless have 99.5% sequence identity. This pattern cannot be easily explained by slipped-strand mispairing, errors in initiation or termination of replication, or any combination of these processes. However, this pattern is consistent with the action of intramolecular recombination. One copy of the region comprising the end of nad5, nad6 + trnE + cob, trnT, the putative IGS, trnP, and part of one of the two CRs may have been excised from one genome, forming a separate minicircle that was then reintegrated into another genome between the second CR and trnF. Such intramolecular recombination has been reported in the mitochondrial genomes of both healthy and diseased tissue in several taxa (Holt, Dunbar, and Jacobs 1997; Lunt and Hyman 1997; Kajander et al. 2000; Tang et al. 2000a; Miller et al. 2004; Yokobori et al. 2004). In addition to this nontandem repeat in P. elongatus, the pattern of widely separated repeat units in S. marginatus is consistent with the action of intramolecular recombination. However, the pattern in S. marginatus is also consistent with retention of these widely separated repeat units since the original duplication event, coupled with slower accumulation of mutations in the repeat units relative to the remainder of the duplicated fragment.

Tandem Repeats in the CR and IGS That Do Not Effect Gene Rearrangement
Only 4 of the 24 genomes analyzed for this study contain a duplication-mediated rearrangement involving the CR or nearby regions of the genome. Thirteen of the remaining 20 genomes contain tandem repeats of noncoding sequence in the CR and/or the IGS, as reported in other taxa (Wallis 1987; Delarbre et al. 2001). These results are summarized in table 1. The length, number, and sequence of the repeat units vary within and among genomes, although in two cases, the same repeats are present in both the IGS and the CR. In the case of nontandem repeats, intramolecular recombination is a more plausible generation mechanism. In the case of tandem repeats, we cannot discriminate between slipped-strand mispairing and intramolecular recombination. The tandem repeats in the Rhyacotriton variegatus genome are unusual and are discussed in more detail below.

View this table: [in this window] [in a new window]   Table 1 Repeat Elements in the IGS between trnT and trnP and the CR Between trnP and trnF in 13 Species' Genomes

 
In the R. variegatus genome, the 5,527-bp IGS contains six copies of an 880-bp tandem repeat. The first 820 bp of this repeat are a cob, and the last 62 bp are an additional copy of trnT. cob is lacking the first 260 bp of the gene, but each of the six copies is >80% identical to the corresponding portion of the functional copy, not including 20 bp of deletions. All copies of the pseudogene contain numerous stop codons when translated in all three reading frames. The copies of trnT are 90% identical to the inferred original copy but, in contrast to cob, all but one may remain functional; their secondary structures appear viable. Following the sixth complete repeat unit, there are 247 bp which appear unrelated to the repeat sequence. In contrast, there are no repeats in the 755-bp CR, nor do the CR and IGS regions share any common sequence.

Heteroplasmy has been reported in some individuals whose mitochondrial genomes contain tandem repeats (Densmore, Wright, and Brown 1985; Wallis 1987; Townsend and Rand 2004). Similarly, patterns in the assembly of individual clone sequences into contigs may suggest heteroplasmy in seven plethodontid genome sequences: Eurycea bislineata, "Bolitoglossa sp. nov.," Batrachoseps wrightorum, Hemidactylium scutatum, Plethodon petraeus, Ensatina eschscholtzii, and Aneides flavipunctatus. However, we note that polymerase chain reaction (PCR)–based evidence for the presence or absence of heteroplasmy is imperfect. A nonheteroplasmic sequence may indicate that multiple mitochondrial genome haplotypes are absent or that, although present, the additional copy or copies were not amplified by PCR. Furthermore, products that appear heteroplasmic may be generated by strand switching during PCR, analogous to slipped-strand mispairing. Finally, heteroplasmy is predicted in genomes with high numbers of tandem repeats, and the assembly of sequences of individual clones into one contig may prove problematic for such genomes even in the absence of heteroplasmy.

Possible Explanations for Extensive Gene Rearrangement in Plethodontids
A high incidence of rearranged mitochondrial genomes may result from any or some combination of four factors: (1) a high rate of mitochondrial genome partial duplication; (2) a low rate of duplication excision; (3) a low mitochondrial genome effective population size, particularly a bottleneck in the number of mitochondrial genomes transmitted to offspring via maternal oocytes; and (4) selection for, or absence of selection against, rearranged/expanded mitochondrial genomes at either the cell or the organism level. No data addressing any of these four factors are currently available for plethodontids. However, studies of selection on mitochondrial genomes containing duplications have been carried out in other systems. Here, we briefly discuss their possible relevance to the high levels of plethodontid mitochondrial instability. We note, however, that inferences drawn from such studies should be considered as hypotheses directing further research and that studies explicitly measuring all four of these factors in plethodontids and other clades with high levels of mitochondrial genome rearrangement are required to draw any firm conclusions.

The fixation of structurally altered mitochondrial genomes in a population depends in part on whether they are selectively advantageous, neutral, or disadvantageous to the organism. Unlike deletions, large duplications of portions of the mitochondrial genome are generally not pathogenic in humans (Tang et al. 2000b; DiMauro and Schon 2003), suggesting that there may not be strong organism-level selection against the duplications in plethodontid mitochondrial genomes. However, a link between compact mitochondrial genomes and metabolic efficiency has been proposed (Selosse, Albert, and Godelle 2001), and high levels of mitochondrial duplications such as those seen in plethodontids lead to measurable reduction in respiratory chain efficiency in human cell lines (Holt, Dunbar, and Jacobs 1997). Salamanders have extremely low aerobic metabolic requirements compared to other ectotherms (Feder 1976). Organism-level selection against the reduction in aerobic efficiency that may result from a mitochondrial duplication is therefore unlikely to be as strong in plethodontids as in other, more highly aerobic tetrapods. Low energetic costs are also broadly correlated with the accumulation of noncoding DNA in the nuclear genome (Gregory 2003). If relatively weak organism-level selection against mitochondrial duplications is contributing to their high incidence in plethodontids, we would predict (1) high mitochondrial genome instability in other organisms with similarly low aerobic metabolic requirements and (2) small fitness differences between plethodontids with and without mitochondrial duplications relative to differences between other, more highly aerobic animals with and without mitochondrial duplications.

Phylogenetic Applications
The analyses presented in this study describe major genome-level mitochondrial mutations for one representative individual from each included lineage. Therefore, they cannot address the extent to which these rearrangements characterize populations, species, or more inclusive phylogenetic groups, with the possible exception of the two Aneides (Moritz and Brown 1987). Aneides flavipunctatus and A. hardii may share a synapomorphic rearrangement, although high levels of such rearrangements within plethodontids necessitate further sampling from this clade to eliminate the possibility of homoplasy. Additional sampling within all rearranged lineages and their close relatives will enable dating of the appearance and persistence of individual rearrangements, thereby addressing the rate at which pseudogenes decay and/or are excised from the genome. Similarly, it will allow dating of the appearance and persistence of tandem duplications within and around the CR and IGS. Finally, it will enable more accurate identification of the ends of duplicated fragments, allowing characterization of secondary structures, repeats, or other genomic sequences that may facilitate duplication both in plethodontids and more generally. In addition to providing further insight into the evolutionary history of vertebrate mitochondrial genomes, this understanding of molecular- and population-level dynamics of mitochondrial genome instability is critical for evaluating the strength of using mitochondrial genome rearrangements as a genome-level character for phylogenetic analysis (Sankoff et al. 1992; Boore and Brown 1998; Macey, Schulte, and Larson 2000).

	Acknowledgements

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 
We thank M. Phillips, R. Zardoya, M. Fujita, R. Gillespie, C. Moritz, D. B. Wake, and M. H. Wake for comments on the manuscript. Part of this work was performed by the University of California Lawrence Berkeley National Laboratory under the auspices of the U.S. Department of Energy, Office of Biological and Environmental Research, contract number DE-AC03-76SF00098. R.L.M. was supported by a National Science Foundation (NSF) predoctoral fellowship and National Institutes of Health training grant. Additional funds came from an NSF doctoral dissertation improvement grant to R.L.M. and David B. Wake (0105824) and the AmphibiaTree Project (NSF EF-0334939).

	Footnotes

 
¹ Present address: Department of Organismal Biology and Anatomy, University of Chicago.

Scott Edwards, Associate Editor

	References

TOP Abstract Introduction Materials and Methods Results and Discussion Acknowledgements References

 

Arndt, A., and M. J. Smith. 1998. Mitochondrial gene rearrangement in the sea cucumber genus Cucumaria. Mol. Biol. Evol. 15:1009–1016.[Abstract]

Boore, J. L. 1999. Animal mitochondrial genomes. Nucleic Acids Res. 27:1767–1780.[Abstract/Free Full Text]

———. 2000. The duplication/random loss model for gene rearrangement exemplified by mitochondrial genomes of deuterostome animals. Pp. 133–147 in D. Sankoff and J. Nadeau, eds. Computational biology series, Vol. 1. Kluwer Academic Publishers, Dordrecht, Netherlands.

Boore, J. L., and W. M. Brown. 1998. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr. Opin. Genet. Dev. 8:668–674.[CrossRef][Web of Science][Medline]

Boore, J. L., M. Medina, and L. A. Rosenberg. 2004. Complete sequences of the highly rearranged molluscan mitochondrial genomes of the scaphopod Graptacme eborea and the bivalve Mytilus edulis. Mol. Biol. Evol. 21:1492–1503.[Abstract/Free Full Text]

Delarbre, C., A.-S. Rasmussen, U. Arnason, and G. Gachelin. 2001. The complete mitochondrial genome of the hagfish Myxine glutinosa: unique features of the control region. J. Mol. Evol. 53:634–641.[CrossRef][Web of Science][Medline]

Densmore, L. D., J. W. Wright, and W. M. Brown. 1985. Length variation and heteroplasmy are frequent in mitochondrial DNA from parthenogenetic and bisexual lizards (genus Cnemidophorus). Genetics 110:689–707.[Abstract/Free Full Text]

Desjardins, P., and R. Morais. 1990. Sequence and gene organization of the chicken mitochondrial genome: a novel gene order in higher vertebrates. J. Mol. Biol. 212:599–634.[CrossRef][Web of Science][Medline]

DiMauro, S., and E. A. Schon. 2003. Mitochondrial respiratory-chain diseases. N. Engl. J. Med. 348:2656–2668.[Free Full Text]

Dowton, M., and A. D. Austin. 1999. Evolutionary dynamics of a mitochondrial rearrangement "hot spot" in the hymenoptera. Mol. Biol. Evol. 16:298–309.[Abstract]

Dowton, M., and N. J. H. Campbell. 2001. Intramitochondrial recombination: is it why some mitochondrial genes sleep around? Trends Ecol. Evol. 16:269–271.[CrossRef][Medline]

Dowton, M., L. R. Castro, S. L. Campbell, S. D. Bargon, and A. D. Austin. 2003. Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J. Mol. Evol. 56:517–526.[CrossRef][Web of Science][Medline]

Eberhard, J. R., T. F. Wright, and E. Bermingham. 2001. Duplication and concerted evolution of the mitochondrial control region in the parrot genus Amazona. Mol. Biol. Evol. 18:1330–1342.[Abstract/Free Full Text]

Feder, M. E. 1976. Lunglessness, body size, and metabolic rate in salamanders. Physiol. Zool. 49:398–406.

Gregory, T. R. 2003. Variation across amphibian species in the size of the nuclear genome supports a pluralistic, hierarchical approach to the C-value enigma. Biol. J. Linn. Soc. 79:329–339.[CrossRef]

Holt, I. J., D. R. Dunbar, and H. T. Jacobs. 1997. Behaviour of a population of partially duplicated mitochondrial DNA molecules in cell culture: segregation, maintenance and recombination dependent upon nuclear background. Hum. Mol. Genet. 6:1251–1260.[Abstract/Free Full Text]

Inoue, J. G., M. Miya, K. Tsukamoto, and M. Nishida. 2003. Evolution of the deep-sea gulper eel mitochondrial genomes: large-scale gene rearrangements originated within the eels. Mol. Biol. Evol. 20:1917–1924.[Abstract/Free Full Text]

Kajander, O. A., A. T. Rovio, K. Majamaa, J. Poulton, J. N. Spelbrink, I. J. Holt, P. J. Karhunen, and H. T. Jacobs. 2000. Human mtDNA sublimons resemble rearranged mitochondrial genoms found in pathological states. Hum. Mol. Genet. 9:2821–2835.[Abstract/Free Full Text]

Kumazawa, Y., H. Ota, M. Nishida, and T. Ozawa. 1998. The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions. Genetics 150:313–329.[Abstract/Free Full Text]

Lee, J.-S., M. Miya, Y.-S. Lee, C. G. Kim, E.-H. Park, Y. Aoki, and M. Nishida. 2001. The complete DNA sequence of the mitochondrial genome of the self-fertilizing fish Rivulus marmoratus (Cyprinodontiformes, Rivulidae) and the first description of duplication of a control region in fish. Gene 280:1–7.[CrossRef][Web of Science][Medline]

Lunt, D. H., and B. C. Hyman. 1997. Animal mitochondrial DNA recombination. Nature 387:247.[CrossRef][Medline]

Macey, J. R., A. Larson, N. B. Ananjeva, Z. Fang, and T. J. Papenfuss. 1997a. Two novel gene orders and the role of light-strand replication in rearrangement of the vertebrate mitochondrial genome. Mol. Biol. Evol. 14:91–104.[Abstract]

Macey, J. R., A. Larson, N. B. Ananjeva, and T. J. Papenfuss. 1997b. Replication slippage may cause parallel evolution in the secondary structures of mitochondrial transfer RNAs. Mol. Biol. Evol. 14:30–39.[Abstract]

Macey, J. R., J. A. Schulte II, and A. Larson. 2000. Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards. Syst. Biol. 49:257–277.[CrossRef][Web of Science][Medline]

Macey, J. R., J. A. I. Schulte, A. Larson, and T. J. Papenfuss. 1998. Tandem duplication via light-strand synthesis may provide a precursor for mitochondrial genomic rearrangement. Mol. Biol. Evol. 15:71–75.[Abstract]

McKnight, M. L., and H. B. Shaffer. 1997. Large, rapidly evolving intergenic spacers in the mitochondrial DNA of the salamander family Ambystomatidae (Amphibia: Caudata). Mol. Biol. Evol. 14:1167–1176.[Abstract]

Miller, A. D., T. T. T. Nguyen, C. P. Burridge, and C. M. Austin. 2004. Complete mitochondrial DNA sequence of the Australian freshwater crayfish, Cherax destructor (Crustacea: Decapoda: Parastacidae): a novel gene order revealed. Gene 331:65–72.[CrossRef][Web of Science][Medline]

Moore, C. A., J. Gudikote, and G. C. Van Tuyle. 1998. Mitochondrial DNA rearrangements, including partial duplications, occur in young and old rat tissues. Mutat. Res. 421:205–217.[Web of Science][Medline]

Moritz, C. 1991. Evolutionary dynamics of mitochondrial DNA duplications in parthenogenetic geckos, Heteronotia binoei. Genetics 129:221–230.[Abstract]

Moritz, C., T. E. Dowling, and W. M. Brown. 1987. Evolution of animal mitochondrial DNA: relevance for population biology and systematics. Ann. Rev. Ecol. Syst. 18:269–292.[CrossRef][Web of Science]

Moritz, C., and W. M. Brown. 1986. Tandem duplication of D-loop and ribosomal RNA sequences in lizard mitochondrial DNA. Science 233:1425–1427.[Abstract/Free Full Text]

———. 1987. Tandem duplications in animal mitochondrial DNAs: variation in incidence and gene content among lizards. Proc. Natl. Acad. Sci. USA 84:7183–7187.[Abstract/Free Full Text]

Mueller, R. L., J. R. Macey, M. Jaekel, D. B. Wake, and J. L. Boore. 2004. Morphological homoplasy, life history evolution, and historical biogeography of plethodontid salamanders inferred from complete mitochondrial genomes. Proc. Natl. Acad. Sci. USA 101:13820–13825.[Abstract/Free Full Text]

Pääbo, S., W. K. Thomas, K. M. Whitfield, Y. Kumazawa, and A. C. Wilson. 1991. Rearrangements of mitochondrial transfer RNA genes in marsupials. J. Mol. Evol. 33:426–430.[CrossRef][Web of Science][Medline]

Quinn, T. W., and D. P. Mindell. 1996. Mitochondrial gene order adjacent to the control region in crocodile, turtle, and tuatara. Mol. Phylogenet. Evol. 5:344–351.[CrossRef][Web of Science][Medline]

Sankoff, D., G. Leduc, N. Antoine, B. Paquin, B. F. Lang, and R. Cedergren. 1992. Gene order comparisons for phylogenetic inference evolution of the mitochondrial genome. Proc. Natl. Acad. Sci. USA 89:6575–6579.[Abstract/Free Full Text]

Selosse, M.-A., B. Albert, and B. Godelle. 2001. Reducing the genome size of organelles favours gene transfer to the nucleus. Trends Ecol. Evol. 16:135–141.[CrossRef][Medline]

Shadel, G. S., and D. A. Clayton. 1997. Mitochondrial DNA maintenance in vertebrates. Annu. Rev. Biochem. 66:409–435.[CrossRef][Web of Science][Medline]

Shao, R., N. J. H. Campbell, E. R. Schmidt, and S. C. Barker. 2001. Increased rate of gene rearrangement in the mitochondrial genomes of three orders of hemipteroid insects. Mol. Biol. Evol. 18:1828–1832.[Free Full Text]

Smith, M. J., D. K. Banfield, K. Doteval, S. Gorski, and D. J. Kowbel. 1989. Gene arrangement in sea star mitochondrial DNA demonstrates a major inversion event during echinoderm evolution. Gene 76:181–185.[CrossRef][Web of Science][Medline]

Stanton, D. J., L. L. Daehler, C. C. Moritz, and W. M. Brown. 1994. Sequences with the potential to form stem-and-loop structures are associated with coding-region duplications in animal mitochondrial DNA. Genetics 137:233–241.[Abstract]

Tang, Y., G. Manfredi, M. Hirano, and E. A. Schon. 2000a. Maintenance of human rearranged mitochondrial DNAs in long-term cultured transmitochondrial cell lines. Mol. Biol. Cell 11:2349–2358.[Abstract/Free Full Text]

Tang, Y., E. A. Schon, E. Wilichowski, M. E. Vazquez-Memije, E. Davidson, and M. P. King. 2000b. Rearrangements of human mitochondrial DNA (mtDNA): new insights into the regulation of mtDNA copy number and gene expression. Mol. Biol. Cell 11:1471–1485.[Abstract/Free Full Text]

Townsend, J. P., and D. M. Rand. 2004. Mitochondrial genome size variation in New World and Old World populations of Drosophila melanogaster. Heredity 93:98–103.[CrossRef][Web of Science][Medline]

Wallis, G. P. 1987. Mitochondrial DNA insertion polymorphism and germ line heteroplasmy in the Triturus cristatus complex. Heredity 58:229–238.

Yamazaki, N., R. Ueshima, J. A. Terrett et al. (12 co-authors). 1997. Evolution of pulmonate gastropod mitochondrial genomes: comparisons of gene organizations of Euhadra, Cepaea and Albinaria and implications of unusual tRNA secondary structures. Genetics 145:749–758.[Abstract]

Yokobori, S.-I., N. Fukuda, M. Nakamura, T. Aoyama, and T. Oshima. 2004. Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes. Mol. Biol. Evol. 21:2034–2046.[Abstract/Free Full Text]

Zardoya, R., E. Malaga-Trillo, M. Veith, and A. Meyer. 2003. Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani. Gene 317:17–27.[CrossRef][Web of Science][Medline]

Zhang, P., Y.-Q. Chen, Y.-F. Liu, H. Zhou, and L.-H. Qu. 2003a. The complete mitochondrial genome of the Chinese giant salamander, Andrias davidianus (Amphibia: Caudata). Gene 311:93–98.[CrossRef][Web of Science][Medline]

Zhang, P., Y.-Q. Chen, H. Zhou, X.-L. Wang, and L.-H. Qu. 2003b. The complete mitochondrial genome of a relic salamander, Ranodon sibiricus (Amphibia: Caudata) and implications for amphibian phylogeny. Mol. Phylogenet. Evol. 28:620–626.[CrossRef][Web of Science][Medline]

Accepted for publication June 21, 2005.

CiteULike    Connotea    Del.icio.us    What's this?

This article has been cited by other articles:

				A. R. Jex, R. S. Hall, D. T. J. Littlewood, and R. B. Gasser An integrated pipeline for next-generation sequencing and annotation of mitochondrial genomes Nucleic Acids Res., November 5, 2009; (2009) gkp883v1. [Abstract] [Full Text] [PDF]

				M. Dowton, S. L. Cameron, J. I. Dowavic, A. D. Austin, and M. F. Whiting Characterization of 67 Mitochondrial tRNA Gene Rearrangements in the Hymenoptera Suggests That Mitochondrial tRNA Gene Position Is Selectively Neutral Mol. Biol. Evol., July 1, 2009; 26(7): 1607 - 1617. [Abstract] [Full Text] [PDF]

				N. C. Sheffield, H. Song, S. L. Cameron, and M. F. Whiting A Comparative Analysis of Mitochondrial Genomes in Coleoptera (Arthropoda: Insecta) and Genome Descriptions of Six New Beetles Mol. Biol. Evol., November 1, 2008; 25(11): 2499 - 2509. [Abstract] [Full Text] [PDF]

				A. J. Heidel and G. Glockner Mitochondrial Genome Evolution in the Social Amoebae Mol. Biol. Evol., July 1, 2008; 25(7): 1440 - 1450. [Abstract] [Full Text] [PDF]

				A. Kurabayashi, M. Sumida, H. Yonekawa, F. Glaw, M. Vences, and M. Hasegawa Phylogeny, Recombination, and Mechanisms of Stepwise Mitochondrial Genome Reorganization in Mantellid Frogs from Madagascar Mol. Biol. Evol., May 1, 2008; 25(5): 874 - 891. [Abstract] [Full Text] [PDF]

				M. K. Fujita, J. L. Boore, and C. Moritz Multiple Origins and Rapid Evolution of Duplicated Mitochondrial Genes in Parthenogenetic Geckos (Heteronotia binoei; Squamata, Gekkonidae) Mol. Biol. Evol., December 1, 2007; 24(12): 2775 - 2786. [Abstract] [Full Text] [PDF]

				S. Tang and B. C. Hyman Mitochondrial Genome Haplotype Hypervariation Within the Isopod Parasitic Nematode Thaumamermis cosgrovei Genetics, June 1, 2007; 176(2): 1139 - 1150. [Abstract] [Full Text] [PDF]

				S. A. M. Amer and Y. Kumazawa The Mitochondrial Genome of the Lizard Calotes versicolor and a Novel Gene Inversion in South Asian Draconine Agamids Mol. Biol. Evol., June 1, 2007; 24(6): 1330 - 1339. [Abstract] [Full Text] [PDF]

				G. C. Gibb, O. Kardailsky, R. T. Kimball, E. L. Braun, and D. Penny Mitochondrial Genomes and Avian Phylogeny: Complex Characters and Resolvability without Explosive Radiations Mol. Biol. Evol., January 1, 2007; 24(1): 269 - 280. [Abstract] [Full Text] [PDF]

				A. Burzynski, M. Zbawicka, D. O. F. Skibinski, and R. Wenne Doubly Uniparental Inheritance Is Associated With High Polymorphism for Rearranged and Recombinant Control Region Haplotypes in Baltic Mytilus trossulus Genetics, November 1, 2006; 174(3): 1081 - 1094. [Abstract] [Full Text] [PDF]

				D. San Mauro, D. J. Gower, R. Zardoya, and M. Wilkinson A Hotspot of Gene Order Rearrangement by Tandem Duplication and Random Loss in the Vertebrate Mitochondrial Genome Mol. Biol. Evol., January 1, 2006; 23(1): 227 - 234. [Abstract] [Full Text] [PDF]

This Article Abstract FREE Full Text (PDF) All Versions of this Article: 22/10/2104    most recent msi204v1 Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (31) Request Permissions Google Scholar Articles by Mueller, R. L. Articles by Boore, J. L. Search for Related Content PubMed PubMed Citation Articles by Mueller, R. L. Articles by Boore, J. L. Social Bookmarking          What's this?

Online ISSN 1537-1719 - Print ISSN 0737-4038

Copyright © 2009 Society for Molecular Biology and Evolution

Oxford Journals Oxford University Press

• Site Map
• Privacy Policy
• Frequently Asked Questions

Other Oxford University Press sites: Oxford University Press Oxford Journals China Oxford Journals Japan American National Biography Booksellers' Information Service Children's Fiction and Poetry Children's Reference Corporate & Special Sales Dictionaries Dictionary of National Biography Digital Reference English Language Teaching Higher Education Textbooks Humanities International Education Unit Law Medicine Music Online Products Oxford English Dictionary Reference Rights and Permissions Science School Books Social Sciences Very Short Introductions World's Classics