# Isolation and Characterization of Microsatellite Loci for the Isopod Crustacean Armadillidium vulgare and Transferability in Terrestrial Isopods.

Isabelle Giraud, Victorien Valette, Nicolas Bech, Frédéric Grandjean, Richard Cordaux

<figure><img src="image_1.png" /></figure>

## To cite this version:

Isabelle Giraud, Victorien Valette, Nicolas Bech, Frédéric Grandjean, Richard Cordaux. Isolation and Characterization of Microsatellite Loci for the Isopod Crustacean Armadillidium vulgare and Transferability in Terrestrial Isopods. PLoS ONE, 2013, 8 (10), pp.e76639. 10.1371/journal.pone.0076639. hal-00872658

HAL Id: hal-00872658

https://hal.science/hal-00872658

Submitted on 6 Jul 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Isolation and Characterization of Microsatellite Loci for the Isopod Crustacean *Armadillidium vulgare* and Transferability in Terrestrial Isopods

Isabelle Giraud*, Victorien Valette*, Nicolas Bech, Frédéric Grandjean, Richard Cordaux*

Université de Poitiers, UMR CNRS 7267 Ecologie et Biologie des Interactions, Equipe Ecologie Evolution Symbiose, Poitiers, France

## Abstract

*Armadillidium vulgare* is a terrestrial isopod (Crustacea, Oniscidea) which harbors *Wolbachia* bacterial endosymbionts. *A. vulgare* is the major model for the study of *Wolbachia*-mediated feminization of genetic males in crustaceans. As a consequence of their impact on host sex determination mechanisms, *Wolbachia* endosymbionts are thought to significantly influence *A. vulgare* evolution on various grounds, including population genetic structure, diversity and reproduction strategies. To provide molecular tools for examining these questions, we isolated microsatellite loci through 454 pyrosequencing of a repeat-enriched *A. vulgare* genomic library. We selected 14 markers and developed three polymorphic microsatellite multiplex kits. We tested the kits on two *A. vulgare* natural populations and found high genetic variation, thereby making it possible to investigate the impact of *Wolbachia* endosymbionts on *A. vulgare* nuclear variation at unprecedented resolution. In addition, we tested the transferability of these kits by cross-species amplification in five other terrestrial isopod species harboring *Wolbachia* endosymbionts. The microsatellite loci showed good transferability in particular in *Armadillidium nasatum* and *Chaetophiloscia elongata*, for which these markers represent promising tools for future genetic studies.

**Citation:** Giraud I, Valette V, Bech N, Grandjean F, Cordaux R (2013) Isolation and Characterization of Microsatellite Loci for the Isopod Crustacean *Armadillidium vulgare* and Transferability in Terrestrial Isopods. PLoS ONE 8(10): e76639. doi:10.1371/journal.pone.0076639

**Editor:** Mark A. Batzer, Louisiana State University, United States of America

**Received** July 26, 2013; **Accepted** August 30, 2013; **Published** October 3, 2013

**Copyright:** © 2013 Giraud et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was funded by a European Research Council Starting Grant (FP7/2007–2013, grant 260729 EndoSexDet) to RC. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have read the journal's policy and have the following conflict: RC is a member of the editorial board of PLOS ONE. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

## Introduction

*Armadillidium vulgare* is a terrestrial isopod (Crustacea, Oniscidea) which exhibits a worldwide distribution. *A. vulgare* harbors alpha-proteobacterial endosymbionts of the genus *Wolbachia* [1,2]. These maternally-inherited, intracytoplasmic bacteria are known to manipulate host reproduction to enhance their own transmission through four different mechanisms: cytoplasmic incompatibility, thelytokous parthenogenesis, male killing and feminization [3,4]. *Wolbachia* endosymbionts are prevalent in terrestrial isopods [5], and *A. vulgare* has emerged as a major model for studying *Wolbachia*-mediated feminization [3,6–8]. In *A. vulgare*, zygotes carrying *Wolbachia* develop a female phenotype, whatever their sex chromosome composition. In particular, genetic males harboring *Wolbachia* are converted into functional females. As a consequence, *A. vulgare* populations in which *Wolbachia* are present show sex ratio distortions towards females, thereby enhancing *Wolbachia* spread in infected populations. In addition, some *A. vulgare* individuals carry another feminizing factor, known as the f element, which may be a fragment of the *Wolbachia* genome carrying feminization information and transferred into the host nuclear genome [3,9]. Furthermore, the occurrence of multiple feminizing factors has generated genetic conflicts in this system, which resulted in the selection of *A. vulgare* nuclear genes resisting feminization [3,10,11]. Thus, sex determination mechanisms are very dynamic

in *A. vulgare*, outlining the prime influence of *Wolbachia*. These endosymbionts are thought to impact *A. vulgare* evolution on various additional grounds, including population genetic structure, diversity and reproduction strategies [12].

Mitochondrial DNA markers have been used in several studies, indicating relatively high variability in *A. vulgare* [1,13–15] as compared to other terrestrial isopod species such as *Porcellionides pruinosus* [16]. A more elaborate understanding of *A. vulgare/Wolbachia* interactions would benefit from information on *A. vulgare* nuclear variation. Recently, five microsatellite markers [17] were used to investigate nuclear variation in *A. vulgare* populations from western France, suggesting a genetic structure compatible with isolation by distance [14]. Although polymorphic, these markers may not be in sufficient number to offer the desired resolution for detecting a possibly subtle impact of *Wolbachia* on *A. vulgare* nuclear variation, population dynamics and evolution. To provide tools for examining these questions, we isolated microsatellite loci through 454 pyrosequencing of a repeat-enriched *A. vulgare* genomic library. We selected markers yielding clear amplification signals and showing appropriate polymorphism levels. Next, we used the candidate loci to develop three polymorphic microsatellite multiplex kits. The transferability of these kits was tested by cross-species amplification in other terrestrial isopod species harboring *Wolbachia* endosymbionts [5,18].

<footnote>* E-mail: richard.cordaux@univ-poitiers.fr</footnote>

<footnote>● These authors contributed equally to this work.</footnote>
# Materials and Methods

## Ethics Statement

No ethics statement was required for the described study. No specific permission was required for sampling the two *A. vulgare* field populations (La Crèche and Beauvoir-sur-Niort, France) because they were located in public areas. Field populations of *Armadillidium nasatum* (Poitiers, France) and *P. pruinosus* (Buxerolles, France) were sampled on private lands after the land owners gave permission to conduct sampling on the sites. None of these species is an endangered or protected species.

## Microsatellite isolation

For genomic library construction, we maximized genomic diversity by using eight *A. vulgare* female individuals selected from the following laboratory lines: BF (bac 377), BH (bac 366), CY (bac 291), CW (bac 49), POA (bac 42), WS (bac 45), WX (matricule 1288) and ZM (bac 47). Total genomic DNA was obtained for each individual by standard phenol-chloroform extraction [19] followed by RNase (10 mg/ml) treatment. DNA concentration was measured using a picogreen assay and equimolar amounts of the eight samples were pooled. The pooled sample was used by GenoScreen (Lille, France) to construct a microsatellite-enriched genomic library, as previously described [20]. The library was sequenced by GenoScreen in a partial 454 GS FLX sequencer run with Titanium chemistry, as previously described [20]. The resulting reads were analyzed with the QDD software [21] to identify reads containing microsatellite motifs and design primers for PCR amplification.

## Locus validation and polymorphism tests

All microsatellite loci with PCR primers designed using QDD were initially tested using two *A. vulgare* female individuals from our laboratory line BF (matricule 2756). Total genomic DNA from the two samples was extracted as above and subjected to whole genome amplification using the illustra GenomiPhi V2 DNA Amplification Kit (GE Healthcare) to generate large enough amounts of template DNA for microsatellite testing. To reduce genotyping costs, each locus was amplified and fluorescently labeled using the M13(-21) primer genotyping protocol [22]. This PCR method uses three primers: a locus-specific forward primer with M13(-21) tail at its 5' end, a locus-specific reverse primer and a universal 6_FAM-labeled M13(-21) primer. PCR amplification was performed in 10 µL reactions, using 0.5 µM of both 6_FAM-M13(-21) and reverse primers, 0.125 µM forward primer, 0.25 U GoTaq DNA Polymerase (Promega), 1X PCR reaction buffer (Promega), 0.2 mM dNTPs (Promega) and 1 µL DNA template. PCR thermal conditions were as previously described [22]. Subsequently, 0.5 µL PCR products were added to 9 µL formamide and 0.35 µL ROX standard (Life Technologies), and resolved by electrophoresis on an ABI PRISM 3130 Genetic Analyzer. Product sizes were determined using the GeneMapper software (Applied Biosystems), followed by eye verification.

Microsatellite loci amplifying in at least one of the two tested individuals and yielding unambiguous amplification signals were further evaluated for their informativeness. Amplification success rates and number of different alleles at each locus were assessed by genotyping a panel of 24 *A. vulgare* individuals (12 males and 12 females) from five laboratory lines: BF (n = 5), BFog (n = 4), WXa (n = 5), BG (n = 5) and ZM (n = 5). First, the 24 samples were subjected to whole genome amplification as described above. Next, microsatellite loci were genotyped using the M13 (-21) primer protocol described above.

## Multiplexing and cross-species amplification

Based on the genotyping results of the 24-individual panel, we selected 14 microsatellite markers for which locus-specific forward primers (without M13(-21) tail) were ordered with labeled dyes (6_FAM, HEX or NED) (Table 1). First, we verified amplification of the 14 markers in simplex PCR conditions on three individuals from each of two *A. vulgare* field populations (La Crèche and Beauvoir-sur-Niort, France, Table 2). Total genomic DNA was extracted as above. All PCR reactions were carried out using the QIAGEN multiplex PCR kit according to the manufacturer's standard microsatellite amplification protocol in a final volume of 10 µL, with an annealing temperature of 58°C and a final concentration of 0.2 µM for each primer. DNA concentrations were adjusted for all individuals between 20 and 60 ng/µL. Next, 1 µL PCR product was added to 18 µL formamide and 0.5 µL ROX standard (Life Technologies). PCR products were resolved by electrophoresis and their size determined as described above.

After simplex PCR verification, we pooled the 14 markers in three multiplex kits (Table 1) according to amplified fragment sizes and dyes to maximize efficiency and minimize costs. The multiplex kits were tested with the same three individuals used for simplex PCR reactions, using the QIAGEN multiplex PCR kit as described above. Identical results were obtained for both simplex and multiplex PCR conditions, thereby validating the use of the three multiplex kits in subsequent analyses. Polymorphism of the 14 microsatellite loci in *A. vulgare* field populations was evaluated by genotyping 20 individuals from each of two populations (La Crèche and Beauvoir-sur-Niort, France) using the multiplex kits.

To investigate transferability of the 14 microsatellite markers, cross-species amplifications were performed in five terrestrial isopod species related to *A. vulgare* and known to harbor *Wolbachia* endosymbionts [5,18]: *A. nasatum* (n = 8) and *P. pruinosus* (n = 8), which were sampled in the field in 2012, and *Chaetophiloscia elongata* (n = 8), *Porcellio scaber* (n = 8) and *Oniscus asellus* (n = 8) from laboratory populations (Table 3). Total genomic DNA was extracted as above. Genotyping was performed using the three multiplex kits.

## Data analyses

To assess genetic variability and transferability of our microsatellite markers, we calculated number of alleles (Na), unbiased expected heterozygosity (He) [23], observed heterozygosity (Ho) and Fis [24] using Genetix version 4.05.2. We computed these genetic indices from two *A. vulgare* populations and from individuals of the five other aforementioned species. Departure from Hardy-Weinberg expectations was assessed for each microsatellite marker using exact tests (5000 permutations), as implemented in GENEPOP version 3.4 [25]. Linkage disequilibrium was assessed for each microsatellite marker using FSTAT version 2.9.3.2 [26] and with 1000 permutations. The level of significance was adjusted for multiple testing using a sequential Bonferroni correction technique [27].

# Results and Discussion

## Locus identification, validation and polymorphism

Sequencing of the microsatellite-enriched library yielded 18,511 reads. The sequence dataset is available in the Dryad database at http://doi.org/10.5061/dryad.md545. Of these, 5073 (27%) reads contained microsatellite motifs according to QDD analysis. Primer pairs were designed for all loci fulfilling our criteria for primer design [20]. The 146 resulting loci comprised 93 di-, 43 tri-, 5 tetra-, 1 penta- and 4 hexanucleotide repeat microsatellites with
<table><caption>**Table 1.** Microsatellite multiplex kits developed for the terrestrial isopod *Armadillidium vulgare*.</caption><thead><tr><th>Locus name</th><th>Repeat motif</th><th>Forward primer (5'-3')</th><th>Reverse primer (5'-3')</th><th>Dye</th></tr></thead><tbody><tr><td colspan="5"><strong>Multiplex kit #1</strong></td></tr><tr><td>AV0023</td><td>AG</td><td>TGGAATTTATGTTGGAGAGGG</td><td>GAGGTTAAGTCTGGGGTCGG</td><td>HEX</td></tr><tr><td>AV0056</td><td>GTT</td><td>TTCAAAGGAGCGTTTGACCT</td><td>AACCACAGCAACAACAGCAG</td><td>6_FAM</td></tr><tr><td>AV0085</td><td>GTA</td><td>CATGCCGTAAGTCCTCTAGACA</td><td>TGTGTTATGGTAATTACATTGAAGTTT</td><td>NED</td></tr><tr><td>AV0086</td><td>TTC</td><td>CCCTTGGCTTCCGATACTT</td><td>TGTCCACAAAGCCAAAAATGA</td><td>HEX</td></tr><tr><td>AV0096</td><td>AAC</td><td>TGGCATAAACCAGCTATAAACC</td><td>TAGTTGCTTTTCCCCTACTTTTG</td><td>6_FAM</td></tr><tr><td colspan="5"><strong>Multiplex kit #2</strong></td></tr><tr><td>AV0002</td><td>ACTCCG</td><td>CGACTCCGACTCCGAATG</td><td>TTCCGACATGTACGATTTTATCA</td><td>6_FAM</td></tr><tr><td>AV0016</td><td>TC</td><td>GCTCATTTATGATCTCGTCGC</td><td>CTCCCACGTGGTTGATCTTC</td><td>6_FAM</td></tr><tr><td>AV0018</td><td>CAA</td><td>GAAGAAATTCAAACTTCACCATCA</td><td>CTTTGAACAGACTTACGAATAACATC</td><td>HEX</td></tr><tr><td>AV0032</td><td>TC</td><td>TTTCAACCTTCCTAACCAAACC</td><td>TTGTTTTATATCCACGACCATCC</td><td>NED</td></tr><tr><td>AV0099</td><td>TG</td><td>CCCCATTTGTGCATGTAGTG</td><td>ACCCCTCGCTTACATTACCC</td><td>HEX</td></tr><tr><td colspan="5"><strong>Multiplex kit #3</strong></td></tr><tr><td>AV0061</td><td>CT</td><td>GTTTGTATGCATTTACCCCTCTTC</td><td>GTATGGAACGAAGGGACCG</td><td>HEX</td></tr><tr><td>AV0063</td><td>TACA</td><td>CAAAACATCTGTACGGATTCCC</td><td>GCCAAACATAAATGCTCGCT</td><td>NED</td></tr><tr><td>AV0089</td><td>CTA</td><td>TTGTTACTTCTACCACACTATTGC</td><td>TGGCTCTATAATGATCAATGGAA</td><td>HEX</td></tr><tr><td>AV0128</td><td>GAT</td><td>TGTCGTTGTGAACAGGCTAAA</td><td>CGTCCGTCGAATGATATTTGT</td><td>6_FAM</td></tr></tbody></table>

doi:10.1371/journal.pone.0076639.t001

5 to 22 repeat units (Table S1). Two *A. vulgare* individuals were genotyped for the 146 microsatellite loci and 41 loci were validated under our amplification conditions and criteria (Table S1). Out of the 41 loci, a first polymorphism analysis based on a 24-individual

panel allowed us to identify 33 polymorphic loci (i.e. 80%) (Table S1). Among these 33 polymorphic loci, we selected 14 loci for inclusion in multiplex kits, according to the following criteria: (i) repeat type (larger motifs favored), (ii) number of different alleles

<table><caption>**Table 2.** Characterization of the 14 microsatellite loci used in multiplex kits in two *Armadillidium vulgare* populations.</caption><thead><tr><th rowspan="2">Populations</th><th colspan="3">La Crèche (N = 20)</th><th colspan="3">Beauvoir-Sur-Niort (N = 20)</th><th colspan="4">Overall population (N = 40)</th></tr><tr><th colspan="3">GPS coordinates</th><th colspan="3">46°21'40.08011"N, 00°18'21.95247"W</th><th colspan="3">46°10'35.91493"N, 00°28'30.45661"W</th><th colspan="4">Parameters</th></tr><tr><th>Parameters</th><th>Na</th><th>He/Ho</th><th>Fis</th><th>Na</th><th>He/Ho</th><th>Fis</th><th>Na</th><th>Size range (bp)</th><th>He/Ho</th><th>Fis</th><th>Multiplex kit #1</th><th>Multiplex kit #2</th><th>Multiplex kit #3</th></tr></thead><tbody><tr><td>AV0023</td><td>1</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>1</td><td>193</td><td>-</td><td>-</td><td>AV0056</td><td>4</td><td>0.57/0.70</td><td>-0.23</td></tr><tr><td>AV0056</td><td>4</td><td>0.57/0.70</td><td>-0.23</td><td>5</td><td>0.46/0.50</td><td>-0.09</td><td>6</td><td>198-219</td><td>0.52/0.60</td><td>-0.15</td><td>AV0085</td><td>2</td><td>0.10/0.10</td><td>-0.03</td></tr><tr><td>AV0085</td><td>2</td><td>0.10/0.10</td><td>-0.03</td><td>4</td><td>0.28/0.30</td><td>-0.09</td><td>4</td><td>175-190</td><td>0.19/0.20</td><td>-0.06</td><td>AV0086</td><td>3</td><td>0.10/0.10</td><td>-0.01</td></tr><tr><td>AV0086</td><td>3</td><td>0.10/0.10</td><td>-0.01</td><td>3</td><td>0.23/0.25</td><td>-0.09</td><td>3</td><td>113-122</td><td>0.17/0.18</td><td>-0.06</td><td>AV0096</td><td>2</td><td>0.51/1.00</td><td><strong>-1.00</strong></td></tr><tr><td>AV0096</td><td>2</td><td>0.51/1.00</td><td><strong>-1.00</strong></td><td>3</td><td>0.55/1.00</td><td><strong>-0.87</strong></td><td>3</td><td>83-107</td><td>0.52/1.00</td><td><strong>-0.94</strong></td><td>AV0002</td><td>5</td><td>0.51/0.45</td><td>0.12</td></tr><tr><td>AV0016</td><td>3</td><td>0.34/0.30</td><td>0.11</td><td>5</td><td>0.54/0.61</td><td>-0.15</td><td>5</td><td>260-308</td><td>0.52/0.53</td><td>-0.02</td><td>AV0018</td><td>4</td><td>0.71/0.70</td><td>0.02</td></tr><tr><td>AV0018</td><td>4</td><td>0.71/0.70</td><td>0.02</td><td>2</td><td>0.14/0.15</td><td>-0.06</td><td>3</td><td>114-126</td><td>0.24/0.23</td><td>0.08</td><td>AV0032</td><td>4</td><td>0.43/0.45</td><td>-0.04</td></tr><tr><td>AV0032</td><td>4</td><td>0.43/0.45</td><td>-0.04</td><td>5</td><td>0.77/0.85</td><td>-0.10</td><td>5</td><td>97-136</td><td>0.74/0.78</td><td>-0.04</td><td>AV0099</td><td>5</td><td>0.49/0.20</td><td><strong>0.60</strong></td></tr><tr><td>AV0099</td><td>5</td><td>0.49/0.20</td><td><strong>0.60</strong></td><td>3</td><td>0.38/0.35</td><td>0.07</td><td>4</td><td>89-105</td><td>0.40/0.40</td><td>0.00</td><td>AV0061</td><td>1</td><td>-</td><td>-</td></tr><tr><td>AV0061</td><td>1</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>1</td><td>138</td><td>-</td><td>-</td><td>AV0063</td><td>3</td><td>0.45/0.50</td><td>-0.12</td></tr><tr><td>AV0063</td><td>3</td><td>0.45/0.50</td><td>-0.12</td><td>3</td><td>0.50/0.40</td><td>0.20</td><td>3</td><td>129-137</td><td>0.48/0.45</td><td>0.05</td><td>AV0089</td><td>2</td><td>0.10/0.10</td><td>-0.03</td></tr><tr><td>AV0089</td><td>2</td><td>0.10/0.10</td><td>-0.03</td><td>1</td><td>-</td><td>-</td><td>2</td><td>85-88</td><td>0.05/0.05</td><td>-0.01</td><td>AV0128</td><td>2</td><td>0.51/0.95</td><td><strong>-0.90</strong></td></tr><tr><td>AV0128</td><td>2</td><td>0.51/0.95</td><td><strong>-0.90</strong></td><td>3</td><td>0.45/0.61</td><td>-0.37</td><td>3</td><td>120-129</td><td>0.49/0.79</td><td><strong>-0.61</strong></td><td></td><td></td><td></td></tr></tbody></table>

Sampled locations, GPS coordinates (longitude and latitude in World Geodetic System 1984) and number of sampled individuals (N) are shown. Number of alleles (Na), size range of alleles (bp), unbiaised expected heterozygosity (He), observed heterozygosity (Ho), and Fis are shown. Significant values (P≤0.01) are shown in bold.
doi:10.1371/journal.pone.0076639.t002
<table><caption>**Table 3.** Transferability of the 14 microsatellite loci used in multiplex kits in five terrestrial isopod species.</caption><thead><tr><th>Species</th><th>Armadillidium nasatum</th><th>Porcellionides pruinosus</th><th>Chaetophiloscia elongata</th><th>Oniscus asellus</th><th>Porcellio scaber</th></tr><tr><th>Populations</th><th>Poitiers</th><th>Buxerolles</th><th>Laboratory line</th><th>Laboratory line</th><th>Laboratory line</th></tr><tr><th>GPS coordinates</th><th colspan="5">46°35'3.77006"N, 00°22'16.07919"E</th></tr><tr><th>Parameters</th><th>Na</th><th>Size range</th><th>He/Ho</th><th>Fis</th><th>Na</th><th>Size range</th><th>He/Ho</th><th>Fis</th><th>Na</th><th>Size range</th><th>He/Ho</th><th>Fis</th></tr></thead><tbody><tr><td colspan="14">Multiplex kit #1</td></tr><tr><td>AV0023</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>1</td><td>193</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td></tr><tr><td>AV0056</td><td>3</td><td>198-209</td><td>0.24/0.13</td><td>0.5</td><td>IN</td><td>-</td><td>2</td><td>209-219</td><td>0.13/0.13</td><td>IN</td><td>-</td><td>6</td><td>231-279</td><td>0.85/0.43</td><td><strong>0.51</strong></td></tr><tr><td>AV0085</td><td>4</td><td>169-178</td><td>0.74/0.38</td><td><strong>0.51</strong></td><td>IN</td><td>-</td><td>2</td><td>175-181</td><td>0.13/0.13</td><td>IN</td><td>-</td><td>-</td><td>1</td><td>175</td><td>-</td><td>-</td></tr><tr><td>AV0086</td><td>1</td><td>110</td><td>-</td><td>-</td><td>3</td><td>107-113</td><td>0.62/0.29</td><td>0.56</td><td>2</td><td>110-113</td><td>0.48/0.67</td><td>-0.43</td><td>2</td><td>110-113</td><td>0.53/0.17</td><td>0.71</td><td>2</td><td>110-113</td><td>0.13/0.13</td><td>-</td></tr><tr><td>AV0096</td><td>4</td><td>77-89</td><td>0.70/0.71</td><td>-0.02</td><td>3</td><td>74-89</td><td>0.44/0.50</td><td>-0.15</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>2</td><td>74-83</td><td>0.53/0.50</td><td>0.06</td><td>3</td><td>74-89</td><td>0.59/0.17</td><td><strong>0.74</strong></td></tr><tr><td colspan="14">Multiplex kit #2</td></tr><tr><td>AV0002</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>AV0016</td><td>2</td><td>94-116</td><td>0.36/0.14</td><td>0.63</td><td>IN</td><td>-</td><td>-</td><td>2</td><td>98-116</td><td>0.53/0.83</td><td>-0.67</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>AV0018</td><td>3</td><td>100-106</td><td>0.54/0.38</td><td>0.32</td><td>IN</td><td>-</td><td>-</td><td>3</td><td>97-135</td><td>0.59/0.67</td><td>-0.14</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>AV0032</td><td>3</td><td>89-99</td><td>0.73/0.60</td><td><strong>0.2</strong></td><td>IN</td><td>-</td><td>-</td><td>2</td><td>89-99</td><td>0.47/0.60</td><td>-0.33</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td><td>95-99</td><td>0.32/0.33</td><td>-0.05</td></tr><tr><td>AV0099</td><td>3</td><td>160-174</td><td>0.69/0.63</td><td>0.1</td><td>IN</td><td>-</td><td>-</td><td>1</td><td>160</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="14">Multiplex kit #3</td></tr><tr><td>AV0061</td><td>2</td><td>144-146</td><td>0.50/0.25</td><td>0.52</td><td>IN</td><td>-</td><td>-</td><td>2</td><td>138-146</td><td>0.13/0.13</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td></tr><tr><td>AV0063</td><td>2</td><td>133-137</td><td>0.20/0.20</td><td>-</td><td>3</td><td>129-137</td><td>0.54/0.43</td><td>0.22</td><td>3</td><td>129-137</td><td>0.67/0.67</td><td>0</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td></tr><tr><td>AV0089</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>2</td><td>85-88</td><td>0.53/0.50</td><td>0.06</td><td>2</td><td>85-88</td><td>0.33/0.38</td><td>-0.17</td><td>2</td><td>85-88</td><td>0.36/0.40</td><td>-0.14</td><td>2</td><td>85-88</td><td>0.17/0.17</td><td>-</td></tr><tr><td>AV0128</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>2</td><td>126-129</td><td>0.53/0.80</td><td>-0.6</td><td>IN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>IN</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Parameters are as described in the legend to Table 2. IN indicates inconsistent amplification of the locus in one species.
doi:10.1371/journal.pone.0076639.t003
scored in 24-individual panel (higher number favored), and (iii) amplification success rate in 24-individual panel (higher number of successfully genotyped individuals favored). The 14 loci were combined in three multiplex kits according to ranges of amplification sizes (Table 1).

We genotyped 40 *A. vulgare* individuals from two field populations with the three multiplex kits. No significant linkage disequilibrium between the different loci was observed after sequential Bonferroni correction. Two loci (AV0023 and AV0061) were monomorphic in the two tested populations (Table 2). However, we kept these two loci in our multiplex kits because they were polymorphic in our initial 24-individual panel, suggesting that they may be informative for other populations. The number of alleles for the 12 other loci ranged from 2 to 6 among the 40 *A. vulgare* individuals, with a mean of 3.9 alleles per locus. Microsatellite loci AV0002, AV0018, AV0056 and AV0099 were highly polymorphic in both populations and they may be particularly relevant markers for analyses requiring high discriminating power, e.g. to investigate paternity in *A. vulgare*. Observed heterozygosity levels varied from 0.1 to 1 in both populations. After sequential Bonferroni correction for multiple testing, loci AV0096 and AV0128 revealed significant excess of heterozygotes in both populations and locus AV0099 showed a significant deficit in La Crèche population. Thus, these loci departed from Hardy-Weinberg equilibrium in concerned populations, likely because of non-exhaustive population sampling. These loci may therefore turn out to be useful for future studies with a more important sampling. The other highlighted polymorphic microsatellites represent a useful set of markers to perform genetic studies on *A. vulgare* and to investigate the impact of *Wolbachia* endosymbionts on *A. vulgare* population genetic structure and evolution.

## Locus transferability in terrestrial isopod species

The 14 newly developed markers were tested in five isopod species. Results are summarized in Table 3. Among these species, *A. nasatum* and *C. elongata* revealed high cross-species transferability with amplification success of 71% (10/14 loci) and 86% (12/14 loci), respectively. Conversely, *O. asellus*, *P. pruinosus* and *P. scaber* revealed moderate amplification success with 21% (3/14 loci), 29% (4/14 loci) and 43% (6/14 loci), respectively. Depending on species, mean number of alleles ranged from 1.8 to 2.8 and mean

expected heterozygosity ranged from 0.40 to 0.53. Results from *A. nasatum* are not really surprising given the close phylogenetic relationship with *A. vulgare*. The important rate of successful cross-amplification found in *C. elongata* is more surprising but microsatellite loci show a reduced number of alleles and low heterozygosity indices in *C. elongata* relative to *A. vulgare*. We detected no linkage disequilibrium whereas departure from Hardy-Weinberg expectations was detected for loci AV0056 and AV0096 in *P. scaber* and for loci AV0032 and AV0085 in *A. nasatum*. The deficit in heterozygotes observed for these loci could be explained by the quite small sampling.

## Conclusions

In sum, our work highlights a large set of microsatellite markers useful for studies on *A. vulgare* and other terrestrial isopod species. The polymorphism of these markers now makes it possible to analyze genetic diversity, population structure and reproduction strategies of *A. vulgare* at unprecedented resolution. A study based on these markers is now underway to analyze the impact of *Wolbachia* bacterial endosymbionts on *A. vulgare* nuclear variation. Moreover, these microsatellite markers showed good transferability in five other terrestrial isopod species, in particular in *A. nasatum* and *C. elongata*, for which these microsatellite markers represent promising tools for future genetic studies.

## Supporting Information

<caption>**Table S1** Information on the 146 microsatellite loci from *Armadillidium vulgare* tested for inclusion in multiplex kits. For each locus the following information is provided: locus name, repeat motif, repeat number, reference sequence, forward primer, reverse primer, PCR product size in reference sequence, and results of experimental tests leading to final selection for three multiplex kits. (XLSX)</caption>

## Author Contributions

Conceived and designed the experiments: RC FG NB. Performed the experiments: IG VV. Analyzed the data: IG VV RC FG NB. Contributed reagents/materials/analysis tools: VV FG RC. Wrote the paper: RC. Coordinated the study: RC. Revised the manuscript: VV FG NB.

## References

1. Cordaux R, Michel-Salzat A, Frelon-Raimond M, Rigaud T, Bouchon D (2004) Evidence for a new feminizing Wolbachia strain in the isopod Armadillidium vulgare: evolutionary implications. Heredity 93: 78-84.

2. Verne S, Johnson M, Bouchon D, Grandjean F (2007) Evidence for recombination between feminizing Wolbachia in the isopod genus Armadillidium. Gene 397: 58-66.

3. Cordaux R, Bouchon D, Greve P (2011) The impact of endosymbionts on the evolution of host sex-determination mechanisms. Trends Genet 27: 332-341.

4. Werren JH, Baldo L, Clark ME (2008) Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol 6: 741-751.

5. Cordaux R, Pichon S, Hatira HB, Doublet V, Greve P, et al. (2012) Widespread Wolbachia infection in terrestrial isopods and other crustaceans. Zookeys 176: 123-131.

6. Rigaud T (1997) Inherited microorganisms and sex determination of arthropod hosts. In: O'Neill SL, Hoffmann AA, Werren JH, Influential passengers: inherited microorganisms and arthropod reproduction. New York: Oxford University Press. pp. 81-101.

7. Rigaud T, Juchault P, Mocquard JP (1997) The evolution of sex determination in isopod crustaceans. Bioessays 19: 409-416.

8. Bouchon D, Cordaux R, Grève P (2008) Feminizing Wolbachia and the evolution of sex determination in isopods. In: Bourtzis K, Miller T, Insect Symbiosis, Volume 3. Boca Raton, FL: Taylor and Francis Group LLC. pp. 273-294.

9. Legrand JJ, Juchault P (1984) Nouvelles donnees sur le determinisme genetique et epigenetique de la monogenie chez le crustace isopode terrestre Armadillidium vulgare Latr. Genet Sel Evol 16: 57-84.

10. Rigaud T, Juchault P (1993) Conflict between feminizing sex ratio distorters and an autosomal masculinizing gene in the terrestrial isopod Armadillidium vulgare Latr. Genetics 133: 247-252.

11. Rigaud T, Juchault P (1992) Genetic control of the vertical transmission of a cytoplasmic sex factor in Armadillidium vulgare Latr. (Crustacea, Oniscidea). Heredity 68: 47-52.

12. Verne S, Moreau J, Caubet Y, Bouchon D, Johnson M, et al. (2007) Male mating success during parturial intermoult in the terrestrial isopod Armadillidium vulgare revealed by the use of a microsatellite locus. Journal of Crustacean Biology 27: 217-219.

13. Rigaud T, Bouchon D, Souty-Grosset C, Raimond R (1999) Mitochondrial DNA polymorphism, sex ratio distorters and population genetics in the isopod Armadillidium vulgare. Genetics 152: 1669-1677.

14. Verne S, Johnson M, Bouchon D, Grandjean F (2012) Effects of parasitic sex-ratio distorters on host genetic structure in the Armadillidium vulgare-Wolbachia association. J Evol Biol 25: 264-276.

15. Grandjean F, Rigaud T, Raimond R, Juchault P, Souty-Grosset C (1993) Mitochondrial DNA polymorphism and feminizing sex factors dynamics in a natural population of Armadillidium vulgare (Crustacea, Isopoda). Genetica 92: 55-60.

16. Michel-Salzat A, Cordaux R, Bouchon D (2001) Wolbachia diversity in the Porcellionides pruinosus complex of species (Crustacea: Oniscidea): evidence for host-dependent patterns of infection. Heredity 87: 428-434.

17. Verne S, Puillandre N, Brunet G, Gouin N, Samollow PB, et al. (2006) Characterization of polymorphic microsatellite loci in the terrestrial isopod Armadillidium vulgare. Molecular Ecology Notes 6: 328-330.
18. Cordaux R, Michel-Salzat A, Bouchon D (2001) Wolbachia infection in crustaceans: novel hosts and potential routes for horizontal transmission. *J Evol Biol* 14: 237-243.

19. Kocher TD, Thomas WK, Meyer A, Edwards SV, Paabo S, et al. (1989) Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc Natl Acad Sci U S A 86: 6196-6200.

20. Malausa T, Gilles A, Meglecz E, Blanquart H, Duthoy S, et al. (2011) High-throughput microsatellite isolation through 454 GS-FLX Titanium pyrosequencing of enriched DNA libraries. Mol Ecol Resour 11: 638-644.

21. Meglecz E, Costedoat C, Dubut V, Gilles A, Malausa T, et al. (2010) QDD: a user-friendly program to select microsatellite markers and design primers from large sequencing projects. Bioinformatics 26: 403-404.

22. Schuelke M (2000) An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol 18: 233-234.

23. Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583-590.

24. Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370.

25. Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Journal of Heredity 86: 248-249.

26. Goudet J (2001) FSTAT, a Program to Estimate and Test Gene Diversities and Fixation Indices. 2.9.3 ed.

27. Rice WR (1989) Analyzing tables of statistical tests. Evolution 43: 223-225.