<title>Introduction</title> <p>Metazoan mitochondrial genomes evolve rapidly, thus creating extensive sequence variation among individuals within species. While mitochondrial DNA (mtDNA) mutations are widely used in evolutionary studies, they also play a central role in apoptosis, disease and aging ##REF##10101183##[1]##, ##REF##11932010##[2]##. Due to essentially maternal inheritance and an mtDNA bottleneck during early oogenesis, most individuals are homoplasmic for a single mtDNA genome (i.e. they carry a single mtDNA genotype) ##REF##11932010##[2]##–##REF##18227871##[5]##. This indicates that newly arising mtDNA mutations segregate rapidly in the female germline. Consistently, evidence from cattle, mice and humans indicates that heteroplasmy (i.e. the mixture of mtDNA genotypes within an organism) is generally resolved to homoplasmy within a few generations ##REF##11932010##[2]##–##REF##9399894##[6]##. Therefore, heteroplasmy is considered as a transitional and short-lived state of metazoan mitochondrial evolution. In sharp contrast, we report an unusual transpecific heteroplasmy that has been stably inherited in oniscids (terrestrial isopod crustaceans) for at least 30 million years.</p>
<title>Materials and Methods</title> <p>Total DNA was extracted from gonads, fat body and nerve tissues of each individual as previously described ##REF##2762322##[19]##. PCR reactions were performed using GoTaq polymerase (Promega), primer combinations shown in ##TAB##0##Table 1## and the following program: 95°C for 3 min, 35 cycles of 95°C for 30 sec, 51∼54°C for 50 sec, 72°C for 1 min 30, and 72°C for 5 min. Purified PCR products were directly sequenced with the BigDye Terminator kit (Applied Biosystems) and analyzed on an ABI Prism 3130 Genetic Analyzer.</p> <p>Endonuclease digestion of PCR products were performed with HpyCH4V (New England BioLabs). This enzyme recognizes a T<underline>G</underline>CA restriction motif (the site corresponding to the 11973 mitochondrial polymorphism is underlined). Digestions were performed with 3.5 units of enzyme (i.e. in excess to avoid partial digestions) at 37°C for 1 h. Digested products were separated on 2% agarose gels in TAE buffer for 1 h 30 at 35 V. Gels were stained with ethidium bromide and examined under UV light. HpyCH4V <italic>in silico</italic> digestion of <italic>A. vulgare</italic> mitochondrial genome ##REF##17906827##[7]## predicted a single undigested band (399 bp) in 11973A individuals and two digested bands (250 bp and 149 bp) in 11973G individuals. 11973G/A heteroplasmic individuals exhibited the three bands.</p>
<title>Results and Discussion</title> <p>We have previously identified by direct DNA sequencing a putative G/A heteroplasmic point mutation at nucleotide position 11973 relative to the mitochondrial genome sequence of the oniscid <italic>Armadillidium vulgare</italic>; this polymorphism occurs in a tRNA anticodon and predicts an alanine/valine alloacceptor tRNA change ##REF##17906827##[7]##. Interestingly, the 11973G/A heteroplasmy was found in all six individuals tested (five females and one male from various geographic origins) ##REF##17906827##[7]##. To independently confirm these results, we designed a PCR-RFLP assay (##FIG##0##Fig. 1##) and verified the 11973G/A heteroplasmy in these samples. PCR-RFLP analysis of three additional <italic>A. vulgare</italic> individuals (two females and one male) indicated that they all shared the 11973G/A heteroplasmy (data not shown).</p> <p>To test whether the observed polymorphism is indeed of mitochondrial origin rather than derived from a nuclear mitochondrial pseudogene, we isolated mtDNA of four <italic>A. vulgare</italic> female samples using an alkaline lysis protocol ##UREF##0##[8]##. Specific primers F1 (<named-content content-type="gene">5′-ACCATCTATGAAGGTTACTGTC</named-content>) and R1 (<named-content content-type="gene">5′-ACCAACAGATTCCATACCCAAG</named-content>) targeting the nuclear actin gene were used in a PCR assay to confirm the absence of nuclear DNA in the samples. We found the 11973G/A heteroplasmy to be present in all four individuals, thus confirming its mitochondrial origin (##FIG##0##Fig. 1A##).</p> <p>Presence of the 11973G/A heteroplasmy in all <italic>A. vulgare</italic> individuals analyzed, regardless of their geographic origin and sex, raises the question of its evolutionary origin. To address this issue, we screened 20 isopod crustacean species (from 13 different families) for the presence of the 11973G/A heteroplasmy by direct DNA sequencing and PCR-RFLP. We detected the 11973G/A heteroplasmy in 16 species comprising 10 divergent families of the crustacean suborder Oniscidea ##UREF##1##[9]## (##TAB##0##Table 1##). These results indicated that the heteroplasmy originally identified in <italic>A. vulgare</italic> actually is a transpecific polymorphism with a very deep evolutionary origin. The fossil record indeed suggests that oniscids arose at least ∼100 million years ago ##UREF##2##[10]## and extent oniscid families are definitely attested by at least ∼34 million years ago ##UREF##3##[11]##. Thus, fossil evidence conservatively suggests that the 11973G/A heteroplasmy may have been present in oniscids for at least ∼30 million years.</p> <p>Sharing of a heteroplasmy in such a wide range of species suggests that it may be stably inherited by oniscids. To gain insight into this question, we investigated transmission of the 11973G/A heteroplasmy across three generations of an <italic>A. vulgare</italic> maternal lineage (##FIG##1##Fig. 2##). The heteroplasmy was found in all individuals of the pedigree (8 females and 6 males), thus demonstrating 100% transmission to both males and females across three generations. Constitutive heteroplasmy has previously been reported in male bivalves ##REF##17681397##[12]##, ##REF##11280005##[13]##. In these species, male and female germlines are homoplasmic for different mtDNA genotypes, each of which is uniparentally transmitted ##REF##17681397##[12]##, ##REF##11280005##[13]##. As a result, females inherit mtDNA following usual strict maternal inheritance and males inherit mtDNA from both parents and, therefore, are heteroplasmic (except in the germline). In sharp contrast with the doubly uniparental mode of mtDNA inheritance of bivalves, both male and female <italic>A. vulgare</italic> gonadic tissues are heteroplasmic (data not shown). Thus, unlike the bivalve system, paternal mtDNA transmission does not need to be invoked in <italic>A. vulgare</italic>. This view is further supported by the fact that we found the 11973G/A heteroplasmy in the parthenogenetic oniscid <italic>Trichoniscus pusillus pusillus</italic> (i.e. females produce eggs that develop without male fertilization) ##UREF##4##[14]##. Therefore, classical strict maternal inheritance of mtDNA in principle can account for the presence of the 11973G/A heteroplasmy in both <italic>A. vulgare</italic> males and females.</p> <p>How, then, does the 11973G/A heteroplasmy persist in oniscids on the long term? Because of mtDNA bottlenecks during oogenesis, persistence of two distinct populations of mtDNA genotypes within all individuals is very unlikely, unless the homoplasmic state at nucleotide position 11973 is lethal in oniscid species. Yet, mtDNA bottlenecks would predict lethal homoplasmic eggs to be frequently produced, which would seem to induce unbearable reproduction costs to females. Alternatively, stable transmission across generations can be achieved if oniscid mitochondria are constitutively heteroplasmic.</p> <p>Consistently, at least two terrestrial isopod species (<italic>A. vulgare</italic> and <italic>Porcellionides pruinosus</italic>) possess an atypical mitochondrial genome formed by a circular ∼28-kb dimer consisting of two ∼14-kb monomers fused in opposite polarities and a linear ∼14-kb monomer ##REF##9872960##[15]##, ##REF##10447705##[16]## (##FIG##2##Fig 3##). The circular and linear molecules may not be at equilibrium and the circular state might be generated from the linear state ##REF##9872960##[15]##. Such trimeric structure could provide an adequate substrate for the emergence of a constitutive intra-mitochondrial heteroplasmy. Strikingly, the three monomer units of the <italic>A. vulgare</italic> mitochondrial genome are identical at the nucleotide level, except for the 11973G/A heteroplasmy ##REF##17906827##[7]##. It is thought that the extreme sequence homogeneity among monomers is maintained by concerted evolution or it may be linked to the replication mechanism of this atypical mitochondrial genome ##REF##17906827##[7]##. If so, maintaining a heteroplasmic state for ∼30 million years despite evolutionary forces otherwise tending to homogenize monomer sequences necessarily requires a selective process.</p> <p>Interestingly, the 11973G/A heteroplasmy occurs in a tRNA anticodon specifying tRNA<sup>Ala</sup> and tRNA<sup>Val</sup>
##REF##17906827##[7]##. This 65 bp-long tRNA exhibits a typical cloverleaf secondary structure (##FIG##3##Fig. 4##), as predicted by the tRNAScan-SE program ##REF##9023104##[17]##, suggestive of biological functionality. In addition, an analysis of the <italic>A. vulgare</italic> mitochondrial genetic code using the GenDecoder program ##REF##16845034##[18]## reveals that <italic>A. vulgare</italic> possesses the standard invertebrate mitochondrial genetic code, confirming the use of codons GCA for alanine (used 118 times) and GTA for valine (used 51 times). Thus, tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> are both necessary for <italic>A. vulgare</italic> mitochondria to perform proper mRNA translation. Because the tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> are apparently encoded nowhere else in the <italic>A. vulgare</italic> mitochondrial genome ##REF##17906827##[7]##, the tRNA locus exhibiting the 11973G/A heteroplasmy is the sole endogenous source of both tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> for <italic>A. vulgare</italic> mitochondria. Since lack of either tRNA could have dramatic consequences on mitochondrial metabolism, the 11973G/A heteroplasmy may have been maintained by a process of balancing selection for the past ∼30 million years of oniscid evolution.</p>
<title>Results and Discussion</title> <p>We have previously identified by direct DNA sequencing a putative G/A heteroplasmic point mutation at nucleotide position 11973 relative to the mitochondrial genome sequence of the oniscid <italic>Armadillidium vulgare</italic>; this polymorphism occurs in a tRNA anticodon and predicts an alanine/valine alloacceptor tRNA change ##REF##17906827##[7]##. Interestingly, the 11973G/A heteroplasmy was found in all six individuals tested (five females and one male from various geographic origins) ##REF##17906827##[7]##. To independently confirm these results, we designed a PCR-RFLP assay (##FIG##0##Fig. 1##) and verified the 11973G/A heteroplasmy in these samples. PCR-RFLP analysis of three additional <italic>A. vulgare</italic> individuals (two females and one male) indicated that they all shared the 11973G/A heteroplasmy (data not shown).</p> <p>To test whether the observed polymorphism is indeed of mitochondrial origin rather than derived from a nuclear mitochondrial pseudogene, we isolated mtDNA of four <italic>A. vulgare</italic> female samples using an alkaline lysis protocol ##UREF##0##[8]##. Specific primers F1 (<named-content content-type="gene">5′-ACCATCTATGAAGGTTACTGTC</named-content>) and R1 (<named-content content-type="gene">5′-ACCAACAGATTCCATACCCAAG</named-content>) targeting the nuclear actin gene were used in a PCR assay to confirm the absence of nuclear DNA in the samples. We found the 11973G/A heteroplasmy to be present in all four individuals, thus confirming its mitochondrial origin (##FIG##0##Fig. 1A##).</p> <p>Presence of the 11973G/A heteroplasmy in all <italic>A. vulgare</italic> individuals analyzed, regardless of their geographic origin and sex, raises the question of its evolutionary origin. To address this issue, we screened 20 isopod crustacean species (from 13 different families) for the presence of the 11973G/A heteroplasmy by direct DNA sequencing and PCR-RFLP. We detected the 11973G/A heteroplasmy in 16 species comprising 10 divergent families of the crustacean suborder Oniscidea ##UREF##1##[9]## (##TAB##0##Table 1##). These results indicated that the heteroplasmy originally identified in <italic>A. vulgare</italic> actually is a transpecific polymorphism with a very deep evolutionary origin. The fossil record indeed suggests that oniscids arose at least ∼100 million years ago ##UREF##2##[10]## and extent oniscid families are definitely attested by at least ∼34 million years ago ##UREF##3##[11]##. Thus, fossil evidence conservatively suggests that the 11973G/A heteroplasmy may have been present in oniscids for at least ∼30 million years.</p> <p>Sharing of a heteroplasmy in such a wide range of species suggests that it may be stably inherited by oniscids. To gain insight into this question, we investigated transmission of the 11973G/A heteroplasmy across three generations of an <italic>A. vulgare</italic> maternal lineage (##FIG##1##Fig. 2##). The heteroplasmy was found in all individuals of the pedigree (8 females and 6 males), thus demonstrating 100% transmission to both males and females across three generations. Constitutive heteroplasmy has previously been reported in male bivalves ##REF##17681397##[12]##, ##REF##11280005##[13]##. In these species, male and female germlines are homoplasmic for different mtDNA genotypes, each of which is uniparentally transmitted ##REF##17681397##[12]##, ##REF##11280005##[13]##. As a result, females inherit mtDNA following usual strict maternal inheritance and males inherit mtDNA from both parents and, therefore, are heteroplasmic (except in the germline). In sharp contrast with the doubly uniparental mode of mtDNA inheritance of bivalves, both male and female <italic>A. vulgare</italic> gonadic tissues are heteroplasmic (data not shown). Thus, unlike the bivalve system, paternal mtDNA transmission does not need to be invoked in <italic>A. vulgare</italic>. This view is further supported by the fact that we found the 11973G/A heteroplasmy in the parthenogenetic oniscid <italic>Trichoniscus pusillus pusillus</italic> (i.e. females produce eggs that develop without male fertilization) ##UREF##4##[14]##. Therefore, classical strict maternal inheritance of mtDNA in principle can account for the presence of the 11973G/A heteroplasmy in both <italic>A. vulgare</italic> males and females.</p> <p>How, then, does the 11973G/A heteroplasmy persist in oniscids on the long term? Because of mtDNA bottlenecks during oogenesis, persistence of two distinct populations of mtDNA genotypes within all individuals is very unlikely, unless the homoplasmic state at nucleotide position 11973 is lethal in oniscid species. Yet, mtDNA bottlenecks would predict lethal homoplasmic eggs to be frequently produced, which would seem to induce unbearable reproduction costs to females. Alternatively, stable transmission across generations can be achieved if oniscid mitochondria are constitutively heteroplasmic.</p> <p>Consistently, at least two terrestrial isopod species (<italic>A. vulgare</italic> and <italic>Porcellionides pruinosus</italic>) possess an atypical mitochondrial genome formed by a circular ∼28-kb dimer consisting of two ∼14-kb monomers fused in opposite polarities and a linear ∼14-kb monomer ##REF##9872960##[15]##, ##REF##10447705##[16]## (##FIG##2##Fig 3##). The circular and linear molecules may not be at equilibrium and the circular state might be generated from the linear state ##REF##9872960##[15]##. Such trimeric structure could provide an adequate substrate for the emergence of a constitutive intra-mitochondrial heteroplasmy. Strikingly, the three monomer units of the <italic>A. vulgare</italic> mitochondrial genome are identical at the nucleotide level, except for the 11973G/A heteroplasmy ##REF##17906827##[7]##. It is thought that the extreme sequence homogeneity among monomers is maintained by concerted evolution or it may be linked to the replication mechanism of this atypical mitochondrial genome ##REF##17906827##[7]##. If so, maintaining a heteroplasmic state for ∼30 million years despite evolutionary forces otherwise tending to homogenize monomer sequences necessarily requires a selective process.</p> <p>Interestingly, the 11973G/A heteroplasmy occurs in a tRNA anticodon specifying tRNA<sup>Ala</sup> and tRNA<sup>Val</sup>
##REF##17906827##[7]##. This 65 bp-long tRNA exhibits a typical cloverleaf secondary structure (##FIG##3##Fig. 4##), as predicted by the tRNAScan-SE program ##REF##9023104##[17]##, suggestive of biological functionality. In addition, an analysis of the <italic>A. vulgare</italic> mitochondrial genetic code using the GenDecoder program ##REF##16845034##[18]## reveals that <italic>A. vulgare</italic> possesses the standard invertebrate mitochondrial genetic code, confirming the use of codons GCA for alanine (used 118 times) and GTA for valine (used 51 times). Thus, tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> are both necessary for <italic>A. vulgare</italic> mitochondria to perform proper mRNA translation. Because the tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> are apparently encoded nowhere else in the <italic>A. vulgare</italic> mitochondrial genome ##REF##17906827##[7]##, the tRNA locus exhibiting the 11973G/A heteroplasmy is the sole endogenous source of both tRNA<sup>Ala</sup> and tRNA<sup>Val</sup> for <italic>A. vulgare</italic> mitochondria. Since lack of either tRNA could have dramatic consequences on mitochondrial metabolism, the 11973G/A heteroplasmy may have been maintained by a process of balancing selection for the past ∼30 million years of oniscid evolution.</p>
<p>Conceived and designed the experiments: VD RC IM. Performed the experiments: VD IM. Analyzed the data: VD CSG DB RC IM. Contributed reagents/materials/analysis tools: CSG DB RC. Wrote the paper: VD RC.</p> <p>Due to essentially maternal inheritance and a bottleneck effect during early oogenesis, newly arising mitochondrial DNA (mtDNA) mutations segregate rapidly in metazoan female germlines. Consequently, heteroplasmy (i.e. the mixture of mtDNA genotypes within an organism) is generally resolved to homoplasmy within a few generations. Here, we report an exceptional transpecific heteroplasmy (predicting an alanine/valine alloacceptor tRNA change) that has been stably inherited in oniscid crustaceans for at least thirty million years. Our results suggest that this heteroplasmy is stably transmitted across generations because it occurs within mitochondria and therefore escapes the mtDNA bottleneck that usually erases heteroplasmy. Consistently, at least two oniscid species possess an atypical trimeric mitochondrial genome, which provides an adequate substrate for the emergence of a constitutive intra-mitochondrial heteroplasmy. Persistence of a mitochondrial polymorphism on such a deep evolutionary timescale suggests that balancing selection may be shaping mitochondrial sequence evolution in oniscid crustaceans.</p>
<p>We thank Houda Ben Afia, Christine Braquart-Varnier, Frédéric Grandjean, Maureen Labarussias, Marion Lachat, Roland Raimond and Sébastien Verne for providing samples. We thank Daniel Guyonnet and Frédéric Chevalier for technical assistance and Laurence Maréchal-Drouard (IBMP Strasbourg, France) for discussions.</p>
<fig id="pone-0002938-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002938.g001</object-id><label>Figure 1</label><caption><title>Assessment of character state at nucleotide position 11973 relative to the <italic>Armadillidium vulgare</italic> mitochondrial genome sequence.</title><p>(A) Presence of the 11973G/A heteroplasmy in mitochondrial DNA extracts of <italic>A. vulgare</italic> individuals. Direct sequencing of a heteroplasmic female sample (left) with forward (top) and reverse (bottom) primers. Partial trace files are shown, with the 11973G/A heteroplasmy highlighted with black triangles. HpyCH4V digestion profiles (right). Lanes: (1) 100 bp molecular ladder (BenchTop, Promega), (2) undigested control amplicon, (3) and (4) digested amplicons of two heteroplasmic female samples. (B) Presence of the 11973G homoplasmic state in <italic>Chaetophiloscia elongata</italic> individuals. Direct sequencing of a homoplasmic individual (left) with forward (top) and reverse (bottom) primers. Partial trace files are shown, with the 11973G homoplasmy highlighted with black triangles. HpyCH4V digestion profiles showing complete digestion of the PCR products (right). Lanes: (1) 200 bp molecular ladder (SmartLadder, Eurogentec), (2) undigested control amplicon, (3) and (4) digested amplicons of two homoplasmic samples.</p></caption></fig> <fig id="pone-0002938-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002938.g002</object-id><label>Figure 2</label><caption><title>Stable inheritance of the 11973G/A heteroplasmy across three generations (I–III) of an <italic>Armadillidium vulgare</italic> maternal lineage (top).</title><p>HpyCH4V digestion profiles showing that all individuals carry the 119973G/A heteroplasmy (bottom). M: 100 bp molecular ladder (BenchTop, Promega).</p></caption></fig> <fig id="pone-0002938-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002938.g003</object-id><label>Figure 3</label><caption><title>Schematic representation of the atypical structure of <italic>Armadillidium vulgare</italic> mitochondrial genome, highlighting the relative position of the 11973G/A heteroplasmy.</title><p>The actual replication mechanism of this mitochondrial genome currently is unknown.</p></caption></fig> <fig id="pone-0002938-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002938.g004</object-id><label>Figure 4</label><caption><title>Secondary structure of the dual tRNA<sup>Ala/Val</sup> from the <italic>Armadillidium vulgare</italic> mitochondrial genome.</title><p>The 11973G/A heteroplasmy is highlighted in the red rectangle, in the anticodon shown in yellow. Orange circles denote conserved nucleotides for mitochondrial tRNAs ##REF##7508516##[20]##.</p></caption></fig>
<table-wrap id="pone-0002938-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002938.t001</object-id><label>Table 1</label><caption><title>Distribution of the 11973G/A mitochondrial heteroplasmy in 20 isopod crustacean species.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">SUBORDER, Family</td><td align="left" rowspan="1" colspan="1">Species</td><td align="left" rowspan="1" colspan="1">Origin</td><td align="left" rowspan="1" colspan="1">N<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">11973<xref ref-type="table-fn" rid="nt102">b</xref>
</td><td align="left" rowspan="1" colspan="1">Assay<xref ref-type="table-fn" rid="nt103">c</xref>
</td><td align="left" rowspan="1" colspan="1">PCR primers<xref ref-type="table-fn" rid="nt104">d</xref>
</td></tr></thead><tbody><tr><td colspan="7" align="left" rowspan="1">FLABELLIFERA</td></tr><tr><td align="left" rowspan="1" colspan="1">Sphaeromatidae</td><td align="left" rowspan="1" colspan="1">
<italic>Dynamene bidentata</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">F, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td colspan="7" align="left" rowspan="1">ONISCIDEA</td></tr><tr><td align="left" rowspan="1" colspan="1">Ligiidae</td><td align="left" rowspan="1" colspan="1">
<italic>Ligia oceanica</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF1/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1">Philosciidae</td><td align="left" rowspan="1" colspan="1">
<italic>Chaetophiloscia elongata</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF1/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Philoscia muscorum</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF1/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1">Agnaridae</td><td align="left" rowspan="1" colspan="1">
<italic>Hemilepistus reaumuri</italic>
</td><td align="left" rowspan="1" colspan="1">Tunisia</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Armadillidae</td><td align="left" rowspan="1" colspan="1">
<italic>Armadillo officinalis</italic>
</td><td align="left" rowspan="1" colspan="1">Turkey</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Cubaris murina</italic>
</td><td align="left" rowspan="1" colspan="1">Guadeloupe</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Armadillidiidae</td><td align="left" rowspan="1" colspan="1">
<italic>Armadillidium assimile</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">2138ND/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Armadillidium depressum</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">2138ND/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Armadillidium nasatum</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF1/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Armadillidium vulgare</italic>
</td><td align="left" rowspan="1" colspan="1">Brazil, France, Greece, Tunisia</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">2138ND/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1">Balloniscidae</td><td align="left" rowspan="1" colspan="1">
<italic>Balloniscus sellowii</italic>
</td><td align="left" rowspan="1" colspan="1">Brazil</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Cylisticidae</td><td align="left" rowspan="1" colspan="1">
<italic>Cylisticus convexus</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Platyarthridae</td><td align="left" rowspan="1" colspan="1">
<italic>Platyarthrus hoffmannseggii</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Platyarthrus caudatus</italic>
</td><td align="left" rowspan="1" colspan="1">Italy</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Porcellionidae</td><td align="left" rowspan="1" colspan="1">
<italic>Porcellio gallicus</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Porcellio spinicornis</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Trachelipodidae</td><td align="left" rowspan="1" colspan="1">
<italic>Trachelipus rathkii</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R, D</td><td align="left" rowspan="1" colspan="1">VF1/1737ND</td></tr><tr><td align="left" rowspan="1" colspan="1">Trichoniscidae</td><td align="left" rowspan="1" colspan="1">
<italic>Trichoniscus pusillus pusillus</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">F, R</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr><tr><td align="left" rowspan="1" colspan="1">Tylidae</td><td align="left" rowspan="1" colspan="1">
<italic>Helleria breviconis</italic>
</td><td align="left" rowspan="1" colspan="1">France</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">G/A</td><td align="left" rowspan="1" colspan="1">R, D</td><td align="left" rowspan="1" colspan="1">VF2/VR2</td></tr></tbody></table></alternatives></table-wrap>
<table-wrap-foot><fn id="nt101"><label>a</label><p>Number of individuals tested.</p></fn><fn id="nt102"><label>b</label><p>Indicates whether all individuals within species were homoplasmic for guanine (G) or heteroplasmic for guanine and adenine (G/A) at mitochondrial genome position 11973.</p></fn><fn id="nt103"><label>c</label><p>Assay used to ascertain 11973 heteroplasmic status: direct sequencing with forward (F) or reverse (R) primer, and HpyCH4V digestion (D).</p></fn><fn id="nt104"><label>d</label><p>Primer pairs used to amplify the 11973 mitochondrial genome region by PCR: VF1 (<named-content content-type="gene">5′-CCCGTTTGAGTGTGGGTTTGA</named-content>), VF2 (<named-content content-type="gene">5′- TGGTTTTTGATGTTGAGATT</named-content>), 2138ND (<named-content content-type="gene">5′-TCCTAGGGATTGGCCATTTA</named-content>), VR2 (<named-content content-type="gene">5′- CGCTTACGTTACGATAAACT</named-content>) and 1737ND (<named-content content-type="gene">5′- TATTTGGGTGCGAGGAACTC</named-content>).</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This research was funded by the Centre National de la Recherche Scientifique (CNRS) and the French Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche. VD was supported by a Ph.D. fellowship from Région Poitou-Charentes. RC was supported by a CNRS Young Investigator ATIP award. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0002938.g001"/> <graphic id="pone-0002938-t001-1" xlink:href="pone.0002938.t001"/> <graphic xlink:href="pone.0002938.g002"/> <graphic xlink:href="pone.0002938.g003"/> <graphic xlink:href="pone.0002938.g004"/>
[{"label": ["8"], "element-citation": ["\n"], "surname": ["Souty-Grosset", "Raimond", "Tourte"], "given-names": ["C", "R", "M"], "year": ["1992"], "article-title": ["D\u00e9terminisme \u00e9pig\u00e9n\u00e9tique du sexe et divergence g\u00e9n\u00e9tique de l'ADN mitochondrial chez Armadillidium vulgare Latr. (Crustac\u00e9 Onisco\u00efde): variabilit\u00e9 inter et intrapopulations."], "source": ["C R Acad Sc Paris"], "volume": ["314"], "fpage": ["119"], "lpage": ["125"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["Schmalfuss"], "given-names": ["H"], "year": ["2003"], "article-title": ["World catalog of terrestrial isopods (Isopoda: Oniscidea)."], "source": ["Stuttgarter Beitr\u00e4ge zur Naturkunde A"], "volume": ["654"], "fpage": ["1"], "lpage": ["341"]}, {"label": ["10"], "element-citation": ["\n"], "surname": ["Alonso", "Arillo", "Barron", "Corral", "Grimalt"], "given-names": ["J", "A", "E", "JC", "JO"], "year": ["2000"], "article-title": ["A new fossil resin with biological inclusions in Lower Cretaceous deposits from Alava (northern Spain, Basque-Cantabrian Basin)."], "source": ["Journal of Paleontology"], "volume": ["74"], "fpage": ["158"], "lpage": ["178"]}, {"label": ["11"], "element-citation": ["\n"], "surname": ["Benton"], "given-names": ["MJ"], "year": ["1993"], "source": ["The fossil record 2"], "publisher-loc": ["London"], "publisher-name": ["Chapman and Hall"], "fpage": ["845"]}, {"label": ["14"], "element-citation": ["\n"], "surname": ["Legrand", "Juchault", "Forest", "von Vaupel Klein"], "given-names": ["JJ", "P", "J", "JC"], "year": ["2006"], "article-title": ["The ontogeny of sex and sexual physiology."], "source": ["The Crustacea"], "publisher-loc": ["Leiden"], "publisher-name": ["Brill Academic Publishers"], "fpage": ["353"], "lpage": ["460"]}]
PLoS One. 2008 Aug 13; 3(8):e2938