Skip to main content

Characterization of polymorphic microsatellite markers for the Japanese endangered land snail Mandarina

Abstract

Objective

Mandarina is an endangered land snail genus of the oceanic Ogasawara archipelago. On Chichijima Island, the largest inhabited island in Ogasawara, this genus is almost extinct in the wild due to predation by invasive species. Although ex situ conservation programs started in 2010, genetic diversity and population structure remain unclear due to a lack of genetic markers with sufficient genetic variation. In this study, we designed polymorphic microsatellite markers of Mandarina to enable genetic analysis and to develop appropriate conservation plans.

Results

Twenty-three polymorphic microsatellite markers were identified from the genomic DNA of wild samples of Mandarina mandarina. We assessed the genetic diversity of each marker. In 16 markers, neither linkage disequilibrium nor deviation from Hardy–Weinberg equilibrium was detected. These 16 markers were tested for multiplex PCR using low-density DNA extracted non-lethally from captive samples of M. mandarina, M. chichijimana and M. suenoae. Of the 16 markers, 15, 12 and 9 were usable for multiplex PCR, respectively. Genetic analysis using these microsatellite loci will be an important resource for the conservation of Mandarina.

Introduction

Mandarina is a Japanese endemic genus of land snails that is highly diversified in the oceanic Ogasawara archipelago [1, 2]. This genus symbolizes the Ogasawara Islands’ value as a World Natural Heritage site, but its population has declined seriously due to the impact of invasive species [3]. In particular, Mandarina species on Chichijima Island, the largest inhabited island in Ogasawara, are near extinction in the wild due to predation by the invasive malacophilous flatworm Platydemus manokwari [2]. Ex situ conservation programs for the Mandarina of Chichijima Island were initiated in 2010 [4].

Genetic analysis is an essential part of ex situ conservation. Bottlenecks at the start of captive breeding and inbreeding in captivity can alter the genetic structure and genetic diversity of populations [5,6,7,8]. Therefore, genetic analysis should be conducted in parallel with captive breeding, but genetic markers with sufficient intraspecific variation have not been developed for Mandarina on Chichijima Island. Although microsatellite markers were used for population genetic analysis of Mandarina species on Hahajima Island [9, 10], these markers did not work for the species on Chichijima Island. In the ex situ conservation of the Hawaiian Achatinella lila, also an oceanic land snail, the population began to decline, likely due to bottlenecks, approximately 10 years after the start of reproduction [7]. In Mandarina, 12 years after the start of the artificial breeding program, the lack of genetic analysis is an issue that must be addressed urgently.

Microsatellite markers are suitable for genetic profiling in the conservation of non-model organisms. Their advantages include a fast evolutionary rate, high intraspecific variation, low-cost experimental methods using multiplex PCR, and the possibility of diverting markers among related species [11, 12]. In addition, the ability to analyze tiny amounts of DNA collected in a non-lethal manner allows for genetic analysis in parallel with captive breeding [7]. Genetic analysis using microsatellite markers is expected to make a significant contribution to the development of breeding plans that avoid risks such as genetic degradation and changes in population structure. In this study, we developed polymorphic microsatellite markers for Mandarina mandarina and evaluated their potential use in the related species, M. chichijimana and M. suenoae. In addition, we examined the feasibility of using these markers in multiplex PCR using tiny amounts of DNA as a non-lethal, low-cost experimental method that is expected in conservation settings.

Main text

Materials and methods

DNA was extracted from the muscle tissues from the foot of M. mandarina using the NucleoSpin Tissue kit (TaKaRa). Microsatellite sequences were isolated by Ecogenics GmbH (Switzerland). Size-selected fragments from the genomic DNA of M. mandarina were enriched for simple sequence repeat (SSR) content using magnetic streptavidin beads and biotin-labelled GATA and GTAT repeat oligonucleotides (Ecogenics GmbH). The SSR-enriched library was analyzed on the Illumina MiSeq platform using the Nano 2 × 250 v2 format (Ecogenics GmbH). Microsatellite loci were selected if the number of repeat motifs was 2–4 and if the number of repeats was greater than 10. In total, 54 primer pairs were developed (Additional file 1) and tested for amplification and polymorphism using 24 wild samples of M. mandarina collected from Chichijima island (Additional file 2). Fluorescent-labelled universal primers were added to the forward primers for fluorescent dye labelling [13]. PIG-tails were attached to the reverse primers to reduce adenylation [14]. Three primer sequences, namely FAM-tail (5ʹ-FAM-GCCTCCCTCGCGCCA-3ʹ), VIC-tail (5ʹ-VIC-GCCTTGCCAGCCCGC-3ʹ), NED-tail (5ʹ-NED-CAGGACCAGGCTACCGTG-3ʹ) [15], were employed as universal primers. PCR amplifications were performed for single primer sets. The PCR solutions were 2 µL and contained the following components: 0.1–5 ng of genomic DNA dried at the bottom of the reaction tube in advance, 1.2 µL of Type-it Multiplex PCR Master Mix (QIAGEN), 0.03 µM tailed forward primer, 0.13 µM fluorescent-labelled universal primer and 0.13 µM reverse primer. The reactions had an initial denaturation step at 95 °C for 15 min, followed by 40 cycles at 94 °C for 30 s, 60 °C for 90 s, 72 °C for 60 s and finally 60 °C for 30 min. Product sizes were determined using the ABI 3130xl Genetic Analyzer and Peak Scanner software (Applied Biosystems) with GeneScan 500 LIZ dye Size Standard v2.0 (Applied Biosystems).

Genetic variability was calculated using GenALEx 6.5 [16]. Linkage disequilibrium and deviations from Hardy–Weinberg equilibrium (HWE) were calculated using Genepop Web version 4.2 [17]. Probabilities of individual identity (PID) and probabilities of sibling individual identity (PID [sib]) were estimated using CERVUS v3.0.7 [18]. For each microsatellite loci and cumulative PID and PID (sib) were calculated by multiplying PID or PID (sib) value across all loci. The multi-locus PID was calculated by multiplying PID for each locus.

Polymorphic microsatellite loci that showed neither linkage disequilibrium nor deviation from HWE were tested for multiplex PCR using captive samples of M. mandarina (N = 34), M. chichijimana (N = 69), and M. suenoae (N = 19) (Additional file 2). M. chichijimana is a sister species of M. mandarina, while M. suenoae is a more distantly related species [2]. DNA was extracted from the muscle tissues approximately 2 mm from the tip of the foot. This method can be used to collect DNA from captive individuals non-lethally, but the quantity of extracted DNA is low. Multiplex PCR was performed by simultaneously amplifying 2–4 sets of markers labelled with different fluorescent colours. (Additional file 3). The PCR solutions were 2 µL and contained the following components: 0.01–5 ng of genomic DNA dried at the bottom of the reaction tube in advance, 1.0 µL of Type-it Multiplex PCR Master Mix (QIAGEN), 0.2 µM fluorescent-labelled forward primer and 0.2 µM reverse primer for each primer sets. The reactions had an initial denaturation step at 95 °C for 15 min, followed by 35 cycles at 94 °C for 30 s, 60 °C for 90 s, 72 °C for 60 s and finally 60 °C for 30 min. Product sizes were determined using an ABI 3130xl Genetic Analyzer and Peak Scanner software (Applied Biosystems) with GeneScan 500 LIZ dye Size Standard v2.0 (Applied Biosystems).

Results

Twenty-three polymorphic microsatellite loci were identified from wild sample of M. mandarina (Table 1). Six to twenty alleles (average = 13.9) were detected; observed heterozygosity were 0.33–0.78 (average = 0.61), while expected heterozygosity showed higher value range of 0.50–0.92 (average = 0.87). PID ranged from 0.011 to 0.060, and PID(sib) ranged from 0.29 to 0.57. Multi-locus PID values, when using eight or more were sufficiently low to distinguish individuals or siblings (P < 0.0001 [19]) (Table 1). Linkage disequilibrium was not detected in any marker, while seven loci showed deviation from HWE (P < 0.05). We tested the remaining 16 loci for multiplex PCR using captive samples.

Table 1 Characterization of the polymorphic microsatellite markers for Mandarina mandarina

In all of the species, multiple loci were identified for which multiplex PCR with low-density DNA was possible. Among captive samples of M. mandarina, one of the 16 loci, Manman_4034, did not show distinct amplification, probably due to interference between primers in multiplex PCR or a low density of DNA. In the other 15 loci, the average number of alleles was 12.2, and the average observed and expected heterozygosity values were 0.63 and 0.83, respectively (Table 2). In the case of M. chichijimana, a sister species, 12 loci showed distinct amplification and polymorphism. In these 12 loci, the average number of alleles was 17.8, and the average observed and expected heterozygosity values were 0.68 and 0.88, respectively (Table 2). Among the samples of M. suenoae, a more distant species from M. mandarina, nine loci showed distinct amplification and polymorphism. Of these nine loci, the average number of alleles was 6.7, and the observed and expected heterozygosity values were both 0.75 (Table 2).

Table 2 Genetic diversity of captive samples of M. manadrina, M. chichijimana and M. suenoae

Discussion

The large number of alleles in the 23 loci (an average of 13.3 for 24 wild samples) was consistent with a previous study on Mandarina on Hahajima Island [10]. Mandarina is estimated to have evolved rapidly in Ogasawara, and high allele diversity may be related to this history. The number of markers available in common corresponded to the phylogenetic relationships among the three species, with more loci available for M. chichijimana, which is closely related to M. mandarina. In addition to the three species in this study, M. hirasei and M. tomiyamai—distantly related species of M. mandarina—are begging benefit from breeding programs. Further testing is needed to determine whether the respective markers can be used in these related species.

We generated polymorphic microsatellite markers that will be useful for the conservation of Mandarina. These markers will contribute to the clarification of the unknown population genetic structure of Mandarina and will help in the establishment of conservation units. Moreover, these markers were available for multiplex PCR using non-lethally obtained DNA. This method allows low-cost analysis of changes in genetic diversity in parallel with breeding, and it is expected to contribute to the development of appropriate breeding plans. Furthermore, the gene flow and population dynamics of Mandarina revealed by microsatellite analysis will be an important resource for elucidating the adaptive radiation process of these land snails on oceanic islands.

Limitation

We did not test the optimal primer set combination for multiplex PCR. Some combinations may increase or decrease the amplification quality of microsatellite loci.

It is possible that deviations from HWE were detected at unproblematic markers because of the small number of wild samples compared to allelic diversity. Finally, these microsatellite markers have not been tested for species from islands other than Chichijima.

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly due to their use in ongoing publication but are available from the corresponding author on reasonable request.

Abbreviations

bp:

Base pair

Na:

Number of alleles

Ho:

Observed heterozygosity

He:

Expected heterozygosity

Fis:

Inbreeding coefficient

HWE:

Hardy–Weinberg equilibrium

PID:

Probabilities of individual identity

PID (sib):

Probabilities of sibling individual identity

PCR:

Polymerase chain reaction

References

  1. Chiba S. Ecological and morphological patterns in communities of land snails of the genus Mandarina from the Bonin Islands. J Evol Biol. 2004;17:131–43.

    Article  Google Scholar 

  2. Chiba S, Cowie RH. Evolution and extinction of land snails on Oceanic Islands. Annu Rev Ecol Syst. 2016;47:123–41.

    Article  Google Scholar 

  3. Ohbayashi T, Okochi I, Sato H, Ono T, Chiba S. Rapid decline of endemic snails in the Ogasawara Islands. Western Pacific Ocean Appl Entomol Zool. 2007;42(3):479–85.

    Article  Google Scholar 

  4. Mori H, Inada M, Chiba S. Conservation programmes for endemic land snails in the Ogasawara Islands: captive breeding and control of invasive species. Tentacle. 2020;28:23–7.

    Google Scholar 

  5. Bouzat JL. Conservation genetics of population bottlenecks: the role of chance, selection, and history. Conserv Genet. 2010;11:463–78.

    Article  Google Scholar 

  6. Campos CI, Martinez MA, Acosta D, Diaz-luque JA, Berkunsky I, Lamberski NL, Cruz-nieto J, Russello MA, Wright TF. Genetic diversity and population structure of two endangered neotropical parrots inform in situ and Ex Situ conservation strategies. Diversity. 2021;13:386.

    Article  Google Scholar 

  7. Price MR, Hadfield MG. Population genetics and the effects of a severe bottleneck in an Ex Situ population of critically endangered hawaiian tree snails. PLoS ONE. 2014;9(12):1–16.

    CAS  Google Scholar 

  8. Tong Y, Durka W, Zhou W, Zhou L, Yu D, Dai L. Forest ecology and management Ex Situ conservation of Pinus koraiensis can preserve genetic diversity but homogenizes population structure. For Ecol Manage. 2020;465:117820.

    Article  Google Scholar 

  9. Davison A, Chiba S. Contrasting response to Pleistocene climate change by ground living and arboreal Mandarina snails from the oceanic Hahajima archipelago. Philos Trans R Soc Lond. 2008;363:3391–400.

    Article  Google Scholar 

  10. Davison A, Chiba S, Kawata M. Characterization of 17 microsatellite loci in the Japanese land snail genera Mandarina, Ainohelix, and Euhadra (Mollusca, Gastropoda, Pulmonata). Mol Ecol Notes. 2004;4:423–5.

    CAS  Article  Google Scholar 

  11. Hodel RGJ, et al. The report of my death was an exaggeration: a review for researchers using microsatellites in the 21st century. Application in Plant Sciences. 2016;4:1600025.

    Article  Google Scholar 

  12. Witzenberger KA, Hochkirch A. Ex situ conservation genetics: a review of molecular studies on the genetic consequences of captive breeding programmes for endangered animal species. Biodivers Conserv. 2011;20(9):1843–61.

    Article  Google Scholar 

  13. Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000;18:233–4.

    CAS  Article  Google Scholar 

  14. Brownstein MJ, Carpten JD, Smith JR. Modulation of non-templated nucleotide addition by Taq DNA polymerase: primer modifications that facilitate genotyping. Biotechniques. 1996;20:1004–10.

    CAS  Article  Google Scholar 

  15. Blacket MJ, Robin C, Good RT, Lee SF, Miller AD. Universal primers for fluorescent labelling of PCR fragments: an efficient and cost-effective approach to genotyping by fluorescence. Mol Ecol Resour. 2012;12:456–63.

    CAS  Article  Google Scholar 

  16. Peakall ROD, Smouse PE. GenALEx 6: genetic analysis in excel. Popula- tion genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  17. Rousset F. genepop’007: a complete re-implementation of the genepop software for windows and linux. Mol Ecol Resour. 2008;8(1):103–6.

    Article  Google Scholar 

  18. Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16(5):1099–106.

    Article  Google Scholar 

  19. Waits LP, Luikart G, Taberlet P. Estimating the probability of identity among genotypes in natural populations: cautions and guidelines. Mol Ecol. 2001;10(1):249–56.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We are sincerely grateful to M. Inada, A. Wakui, K. Ashizawa, S. Wada for their help in providing the materials. We also thank T. Hirano and D. Yamazaki for their helpful supports of the experiments. This research was supported by the Kanto Regional Environment Office, the Japan Ministry of the Environment.

Funding

This work was supported by Ministry of the Environment, Japan.

Author information

Authors and Affiliations

Authors

Contributions

SN, HM and SC conceived of the study. HM and SC provided critical samples. SN and KM performed DNA extractions. KM designed the primers. SN undertook fragment analysis, genotyping and data analysis. SN wrote the initial manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shu Nishida.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

List of microsatellite markers tested in this study. The table includes locus name, repeat motif, amplicon size, forward and reverse primer sequences, and fluorescent label. A dagger (‡) indicates that loci were monomorphic. An asterisk (*) indicates that loci showed unclear or no amplification.

Additional file 2.

List of samples used in this study. The table includes sample name, species name, wild or captive, sampling year (or birth year in the case of captive samples), and sampling site (or founders’ sampling site in the case of captive samples).

Additional file 3.

List of combinations of markers in the multiplex PCR. The table includes combination name, loci in each combination, and fluorescent label. A dagger (‡) and an asterisk (*) indicate that loci were not available for M. chichijimana and M. suenoae, respectively.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nishida, S., Mimura, K., Mori, H. et al. Characterization of polymorphic microsatellite markers for the Japanese endangered land snail Mandarina. BMC Res Notes 15, 255 (2022). https://doi.org/10.1186/s13104-022-06147-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13104-022-06147-4

Keywords

  • Nuclear microsatellite
  • Mandarina
  • Mandarina mandarina
  • Endangered species
  • Oceanic land snails
  • Ogasawara islands