:: Korean Journal of Applied Entomology ::
Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1225-0171(Print)
ISSN : 2287-545X(Online)
Korean Journal of Applied Entomology Vol.56 No.3 pp.229-239
DOI : https://doi.org/10.5656/KSAE.2016.08.0.036

Identification of the South Korean Hermit Beetle (Coleoptera: Scarabaeidae: Cetoniinae)

Taeman Han, In Gyun Park, Ki-Kyoung Kim1, Sergey Ivanov2, Haechul Park*
Applied Entomology Division, Department of Agricultural Biology, National Institute of Agricultural Science, 166 Nongsaengmyeong-ro, Iseo-myeon, Wanju-gun, Jeollabuk-do 55365, Korea
1Biological Resources Research Department, National Institute of Biological Resources, Incheon 22689, Korea
2Far Eastern Branch, FSI “Russian Maritime Register of Shipping,” Str. Stanyukovicha, 29-A, Vladivostok 690950, Russia
Corresponding author:culent@korea.kr
20160628 20160801 20161102


The South Korean hermit beetle has previously been identified as Osmoderma opicum, which is distributed in Japan. Because of its rarity, this species is classified as an endangered species. To date, however, the identity of this species in South Korea has not been conclusively confirmed. To assess the taxonomic status of the hermit beetle occurring in South Korea, we performed a comparative study with the beetle’s Eurasian congeners, based on morphological examination and molecular analysis using COI gene sequences. The results clearly showed that the South Korean hermit beetle is identical to Osmoderma caeleste, which has been described from the Russian Far East. Therefore, we suggest that the taxonomic identity of the South Korean hermit beetle is O. caeleste.

남한산 큰자색호랑꽃무지의 종 동정 (딱정벌레목, 풍뎅이과, 꽃무지아과)

한 태만, 박 인균, 김 기경1, 세르게이 이바노브2, 박 해철*
국립농업과학원 농업생물부 곤충산업과
1국립생물자원관 동물자원과
2러시아 선박해양등록극동지소


남한산 큰자색호랑꽃무지는 그 동안 일본에 분포하는 Osmoderma opicum으로 동정되어 왔으며, 국내에서는 그 희귀성에 의해 멸종위기종으 로 다루어지고 있지만, 남한 개체군에 대한 면밀한 종의 실체에 대한 구명은 없었다. 이에 대해, 남한산 큰자색호랑꽃무지와 유라시안 근연종에 대 해 형태 및 COI유전자 염기서열을 비교를 실시하였다. 결과적으로, 남한산 개체군은 형태적으로 극동러시아에 분포하는 O. caeleste와 동일하였다. 따라서, 남한산 큰자색호랑꽃무지는 O. caeleste임을 제안한다.

    Rural Development Administration

    The genus Osmoderma LePeletier de Saint-Fargeau and Serville consists of 12 species, which occur throughout the Palearctic (9 species) and Nearctic (3 species) regions (Audisio et al., 2007). Species of Osmoderma are commonly referred to as “hermit beetles” because they live in hollows in old broadleaved trees and are rarely observed in nature (Audisio et al., 2009). European species of these beetles (mainly O. eremita) are regarded as one of the most important “flagship species” and “umbrella species” for the conservation of local saproxylic communities (Ranius, 2002a, 2002b; Audisio et al., 2009). Among these, five species are listed as near threatened to endangered in the IUCN Red List of Threatened Species by virtue of their rarity and association with a specific habitat type (Nardi and Micó, 2010; Nieto and Alexander, 2010).

    In the Korean Peninsula, there are a total of four species and one subspecies in the genus Osmoderma. O. opicumLewis, 1887 has been recorded in South Korea by Cho (1969), whereas three species, O. barnabitaMotschulsky, 1845 described by Stebnicka in 1980, O. davidis forma amurensisGusakov, 2002 described by Gusakov (2002), and O. caelesteGusakov, 2002 described by Audisio in 2007, are listed in North Korea, although without any reliable specimen identification. Furthermore, recently, a new subspecies, O. opicum coreanum, was also described from North Korea by Tauzin (2013).

    It is noteworthy that South Korean specimens of Osmoderma have been extremely rarely observed. Only four specimens, consisting of three females and a male, have been collected in South Korea since 1969 (Kim, 2011). The South Korean population, which is considered to be O. opicum, has accordingly been registered as a critically endangered species in the Korean Red Data Book (Lee, 2013). In order to determine accurate species traits of the South Korean hermit beetle, in 2012, we first obtained a Korean female specimen, which we assumed to be O. opicum (specimen no. 3257) and six authenticated Japanese specimens of O. opicum. However, these two geographically separated populations differed in several of their external features. It was accordingly suggested that the South Korean hermit beetle might have been incorrectly recorded as O. opicum. Therefore, over the past 4 years, we have attempted to obtain a larger number of specimens of the Korean hermit beetle as well as voucher specimens of the close relatives in neighboring regions. Recently, three appreciable studies on the European hermit beetles, carried out mainly using molecular and distributional data, have contributed to a reasonable and objectively supported taxonomic arrangement (Audisio et al., 2007, 2009; Landvik et al., 2013). It was therefore possible to compare our molecular data with the results obtained for European hermit beetles.

    The South Korean hermit beetle is an important insect species in the conservation management of threatened species comprising forest saproxylic communities. Nonetheless, as mentioned above, there is still a lack of accurate information on the exact species identity and ecological traits of the Korean population that we seek to conserve.

    Materials and Methods

    Specimen collection and morphospecies examination

    To assess the exact species status of the South Korean hermit beetle, we obtained a total of 16 dried specimens belonging to three Osmoderma species: six South Korean specimens, which have been recognized as O. opicum, two of which were previously recorded as voucher specimens (nos. 7299 and 9572) by Kim (2011); five specimens of true O. opicum from Honshu, Japan; three specimens of true O. caeleste from Far East Russia; and two specimens of O. barnabita from West Russia. We were, however, unable to obtain any specimens of O. opicum coreanum and O. davidis (Table 1).

    For simultaneous examination of both morphological and molecular characters, we first observed the general structure of the specimens under a stereoscopic microscope (MZ 16A and MZ 6; Leica, Solms, Germany) for morphospecies identification and diagnostic characteristic delimitation. All of the examined materials and genomic DNA stocks are preserved in the insect collection at the National Institute of Agricultural Science, Jeonju, Korea. The provincial abbreviations of the collecting sites within Korea are as follows: GW, Gangwondo; GB, Gyeongsangbuk-do.

    DNA extraction, PCR strategy, and DNA analyses

    To reconfirm the specific status of the specimens, we performed molecular analysis using sequences of the cytochrome c oxidase subunit I (COI) gene, which is useful in identifying species within a closely related species group (Han et al., 2012). For genomic DNA extraction, we removed the genitalia from the abdomens of both sexes and used them directly to extract genomic DNA using non-destructive methods (Han et al., 2012). The remaining genital structures were cleaned by heating with 10% KOH solution in a Wise Therm ®HB-48P heating block at 60℃ for 1-2 h and then preserved in microvials containing glycerine for morphological examination.

    For PCR based on dried specimens, we first attempted to amplify the DNA barcoding region (658 bp) of the COI gene using the universal primer set LCO1490/HCO2198 (Folmer et al., 1994); however, we were unable to obtain amplicons for most samples; this should be expected considering the likelihood of DNA degradation in specimen maintained under unknown dried conditions (Goldstein and Desalle, 2003; Hajibabaei et al., 2006; Wandeler et al., 2007). To retrieve the barcoding regions of the COI gene, we adopted the PCR methodology focused on old museum specimens, as described by Han et al. (2014), and designed 13 new specific primers for the COI gene of O. opicum and O. caeleste using Primer3 (Untergrasser et al., 2012) (Table 2). PCR was conducted using AccuPower ProFi Taq PCR PreMix (Bioneer, Daejeon, Korea). The amplification profile consisted of pre-denaturation for 5 min at 95℃, 35-45 cycles of 30 s at 95℃ for denaturation, 30 s at 48-53℃ for primer annealing, and then 45 s at 72℃ for extension, with a final extension for 5 min at 72℃. To confirm product size and purity, the PCR amplicons were assessed by 0.7% agarose gel electrophoresis using Safe-Pinky (Gen DEPOT, Barker, TX, USA). The amplicons successfully obtained were purified using the AccuPrep PCR Purification Kit (Bioneer, Daejeon, Korea) following the manufacturers protocol. Sequencing was performed commercially by Macrogen (Seoul, Korea). All sequences were generated In both directions. The quality and potential polymorphic sites of the analyzed sequences were checked using Chromas 2.33 (Technelysium Pty Ltd, Australia).

    COI sequence analysis

    In this study, sixteen COI sequences from 16 dried specimens were successfully generated. In addition, 32 COI sequences of five species published in two previous studies (Audisio et al., 2009; Landvik et al., 2013) were downloaded from the GenBank (http://www.ncbi.nlm.nih.gov/genbank) (Table 1). Data sets of the nucleotides of COI genes were aligned in MEGA 5.2 (Tamura et al., 2011) using ClustalW with the default settings (Gap Opening Penalty = 15, Gap Extension Penalty = 6.66 in both pairwise and multiple alignments). The anterior and posterior regions of uncertain alignment were eliminated from the data matrix. The COI alignment for reading frames was checked manually by translating sequences into amino acids to identify stop codons and potential shifts. Sequences were finally trimmed to 657-bp fragments. COI sequences generated in this study are available from GenBank under the accession numbers KX346561-76 (Table 1). For molecular species identification, a neighbor-joining (NJ) analysis (Saitou and Nei 1987) was performed using MEGA 5.2. Genetic distances within and between species of hermit beetle were calculated using Kimura’s two-parameter test (Kimura, 1980) in accordance with Nei’s empirical guidelines (Nei 1991, 1996).

    Results and Discussion

    Morphospecies identification of the Korean hermit beetle

    We examined diagnostic morphological features for each specimen and species. The South Korean hermit beetle (Fig. 1 A and B) was revealed to be identical to the true O. caeleste from Far East Russia (Fig. 1C and D). These taxa share the same characteristics, such as an arched frontal margin of the head, a reverse triangular depression on the pronotum (Fig. 2A-D) and shortly developed parameres of the aedeagus (Fig. 3A and B). In contrast, O. caeleste can easily be distinguished from the true O. opicum on the basis of the truncated frontal margin of the head, the longitudinal and narrow groove on the pronotum (Fig. 2E and F), and the distinctly elongated parameres of the aedeagus (Fig. 3C). O. barnabita (Fig. 1G and H) is also distinctly separated from the South Korean hermit beetle by the darker chestnut color of the body, the longitudinal and narrow groove, the rather spare punctures on the pronotum (Fig. 2G and H), and the shortly developed, but with the lobe likely expanded base of parameres of the aedeagus (Fig. 3D). In the Far East Asia, additional related species, O. davidis is distributed in Russia, E. China and Mongolia. The South Korean hermit beetle was compared with this species according to the literatures by Boucher (2002) and Tauzin (2013). O. davidis has the ginkgo leaf-liked pronotal groove and the elongate parameres of the aedeagus. It is indicated that the South Korean hermit beetle and O. davidis are a distinctly separated species, respectively.

    Tauzin (2013) described a subspecies, Osmoderma opicum coreanum, from Hamgyong province, North Korea based on a male holotype, which was provided Dr. Jingke Li (see Tauzin 2013, p. 2). However, we must point out to have doubt on the original locality of the specimen of the new subspecies, because Dr. Li has often provided to purchasers many insect specimens collected from adjacent countries of Asia such as China, Russia, and Japan as though they were from North or South Korea. The existence of this subspecies in Korea is needed to reconfirm in the further study.

    COI sequence analysis

    A total of 16 COI sequences representing three morphospeices, O. opicum, O. barnabita, and O. caeleste including the Korean population, were successfully generated. Before analyzing the COI sequences, we confirmed the absence of any putative Numts and heteroplasmies from the analyzed sequences. One sequence (specimen no. 7299), only partly sequenced in 226 bp fragment, was confirmed to be identical with the other sequences of the Korean hermit beetle, but not including in the calculation of genetic distance and in NJ analysis.

    The NJ tree topology (Fig. 4) and the pairwise distances of COI sequences (Table 3) showed that the Korean specimens were strongly clustered to those of the Far East Russian O. caeleste with low intraspecific (range: 0-0.3%). But, this clade was distinctly separated from other species clades by large interspecific distances (range: 13.7-19.4%). This result indicates that the South Korean hermit beetle and the true O. caeleste from the Far East Russia are the same species. The Japanese species, O. opicum, was also solely clustered with low intraspecific (range: 0-0.8%) and large interspecific distances (range: 12.1-16.0%). The interspecific distances between O. caeleste and O. opium were significantly large, ranging from 16.6% to 17.2%. This result insinuates that these two species have undergone an independent evolution, respectively, and more early separated from the extant European species. Two West Russian specimens of O. barnabita analyzed in this study were clustered with the previous sequences of O. barnabita (Audisio et al., 2009; Landvik et al., 2013). Therefore, the South Korean hermit beetle is identified to O. caeleste based on both of morphological and DNA barcoding analyses.

    Systematic Accounts

    Family Scarabaeidae Latreille

    Subfamily Cetoniinae Leach

    Tribe Trichiini Fleming

    Genus Osmoderma LePeletier and Audinet-Serville

    Osmoderma caeleste (Gusakov, 2002)

    Osmoderma barnabitum var. castaneumTauzin, 1994: 204; Tauzin, 2006: 43 and 44.

    Gymnodus caelestisGusakov, 2002: 13 and 39 (Primorski Krai, Far East Russia).

    Osmoderma sikhotenseBoucher, 2002: 426 (Sikhote-Aline, Far East Russia).

    [Korean records]

    Osmoderma opicum Lewis: Cho, 1969: 688 (Atlas with a female photo, misidentification); Kim, 1978: 379 (Distributional map, listed); Shin and Yoon, 1994: 153 (Check list); Kim, 1998: 358 (List); Kim, 2001: 132 (Examined five specimens, misidentification); Kim, 2011: 33 (Examined five specimens, misidentification).

    Material examined. South Korea: 1 male, Gangwon-do, Yanggu, 26. IX. 2011, collector unknown (voucher no. 3257); 1 female, Gangwon-do, Injae, Jogyeong-dong, Garigol, 26. VIII. 2000, H. N. Seok (voucher no. 7299); 1 female, Gangwondo, Injae, Buk-myeon, Yongdaeri, 31. VII. 2013, Lim et al. (voucher no. 9571); 1 male, Gangwon-do, Hwacheon, Sangseomyeon, Guunri, 9. VIII. 2000, J. R. Cha (voucher no. 9572); 1 female, Gangwon-do, Pyeongchang, Jinbu-myeon, Mt. Odae, 8. VII. 2010, collector unknown (voucher no. 7387); 1 female, Gyeongsangbuk-do, Yongyang, Ilweol-myeon, Mt. Ilweol, XI. 2008, Y. Y. Baik (voucher no. 9596).

    Distribution. South Korea, North Korea (?), E. China (?), Mongolia (?), and the Far East Russia.

    Remarks. We examined the female photo represented in Cho (1969) and two materials (nos. 7299 and 9572) used in Kim (2001, 2011). Based on our morphological and molecular examinations, three previous South Korean records were misidentifications of O. caeleste. Tauzin (2013) mentioned that O. caeleste is distributed in the Far East Russia, Beijing in China, Inner Mongolia and Korea. Among them, the distributions of China, Mongolia and Korea were listed having no the basis of examined specimen records. Therefore, it is needed that the exact distributional range of this species has to reexamine in the further study.


    We would like to thank Mr. Yoohyun Baek, Mr. Hyunkyu Jang, Mr. Minchul Kwon for providing valuable specimens for this study.

    This study was carried out with the support of the Cooperative Research Program for Agricultural Science & Technology Development (Project No. PJ00939502), Rural Development Administration, Republic of Korea.


    Osmoderma spp. A and B: O. caeleste from Korea (A: no. 9572, male; B: no. 7387, female); C and D: O. caeleste from the Far East Russia (C: no. 9412, male; D: no. 9411, female); E and F: O. opicum from Japan (E: no. 7295, male; F: no. 7298, female); G and H: O. barnabita from West Russia (G: no. 7362, male; H: no. 7372, female).


    The dorsal aspects of the head and the pronotum of Osmoderma spp. A and B: O. caeleste from Korea (A: no. 9572, male; B: no. 7387, female); C and D: O. caeleste from the Far East Russia (C: no. 9412, male; D: no. 9411, female); E and F: O. opicum from Japan (E: no. 7295, male; F: no. 7298, female); G and H: O. barnabita from West Russia (G: no. 7362, male; H: no. 7372, female).


    Aedeagus of Osmoderma spp. A: O. caeleste from Korea (A: no. 9572, male); B: O. caeleste from the Far East Russia (B: no. 9412); C: O. opicum from Japan (C: no. 7295, male); D: O. barnabita from West Russia (D: no. 7362, male).


    Neighbor-Joining tree of Osmoderma spp. inferred from cytochrome c oxidase subunit I (COI ). Terminal taxa in bold denote the examined specimens in this study.

    The examined specimens and information for the successful sequences of COI gene
    *Denotes outgroup taxon extracted from GenBank. GBAn is denoted the gene bank accession number. Numbers of 18 to 50 are extracted from GenBank.
    Primers used in the study
    Numbers refer to the position of the primer 5'-end mapped on DNA barcoding region in COI gene (Hebert et al., 2003).
    Pairwise distances of COI within and between Osmoderma spp.
    Numbers are indicated as mean (minimum-maximum) of the pairwise distance.


    1. Audisio P , Brustel H , Carpaneto G M , Coletti G , Mancini E , Piattella E , Trizzino M , Dutto M , Antonini G , Biase A (2007) Updating the taxonomy and distribution of the European Osmoderma, and strategies for their conservation , Fragmenta Entomol. Roma, Vol.39 (2) ; pp.273-290
    2. Audisio P , Brustel H , Capraneto GM , Goletti G , Mancini E , Trizzino M , Antonini G , De Biase A (2009) Data on molecular taxonomy and genetic diversification of the European Hermit beetles, a species complex of endangered insects (Coleoptera: Scarabaeidae, Cetoniinae, Osmoderma) , J. Zool. Sys. Evol. Res, Vol.47 (1) ; pp.88-95
    3. Boucher S (2002) Un nouvel Osmoderma endémique des monts Sikhote-Aline (Coleoptera, Cetoniidae) , Bull. Soc. Entomol. Fr, Vol.104 (4) ; pp.425-427
    4. Cho P S (1969) Illustrated encyclopedia of fauna & plora of Korea (X): Insecta (II), Samhwa Press,
    5. Folmer O , Black M , Hoeh W , Lutz R , Vrijenhoek R (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates , Mol. Mar. Bio. Biotech, Vol.3 ; pp.294-299
    6. Goldstein P Z , Desalle R (2003) Calibrating phylogenetic species formation in a threatened insect using DNA from historical specimens , Mol. Ecol, Vol.12 ; pp.1993-1998
    7. Gusakov A A (2002) Scarabaeid beetles of the subfamily Trichiinae (Coleoptera, Scarabaeidae) in the Fauna of Russia. Species of the genus Gymnodus Kirby, 1827 , Eidos, Kaluga,
    8. Hajibabaei M , Smith MA , Janzen DH , Rodriguez JJ , Whitfield JB , Hebert PDN (2006) A minimalist barcode can identify a specimen whose DNA is degraded , Mol. Ecol. Not, Vol.6 ; pp.959-964
    9. Han T , Kang T , Jeong J , Lee Y , Chung H , Park S , Lee S , Kim K , Park H (2012) Pseudocryptic speciation of Chrysochroa fulgidissima (Coleoptera: Buprestidae) with two new species fromKorea, China and Vietnam , Zoo. J. Linn. Soc, Vol.164 ; pp.71-98
    10. Han T , Lee W , Lee Y , Kim N , Lee S , Park H (2014) Barcoding old Korean lepidopteran specimens using newly designed specific primer pairs , J. Asia-Pacific Entomol, Vol.17 ; pp.679-684
    11. Hebert PDN , Cywinska A , Ball SL , deWaard JR (2003) Biological identifications through DNA barcodes , Proc. R. Soc. B Biol. Sci, Vol.270 ; pp.313-321
    12. Kim CW (1978) Distribution atlas of insects of Korea, Korea University Press,
    13. Kim JI (1998) Decreasing state of Coleoptera population in Korea , J. Kor. Biota, Vol.3 ; pp.351-361
    14. Kim J I (2001) Scarabaeoidea II , Ins. Kor., Suppl, Vol.17 ; pp.1-197
    15. Kim J I (2011) Insect fauna of Korea. Pleurosticti, Scarabaeoidea, Coleoptera, Insecta, Arthropoda. Insect fauna of Korea , Nat. Inst. Bio. Res, Vol.12 (1) ; pp.1-263
    16. Kimura M (1980) A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences , J. Mol. Evol, Vol.16 ; pp.111-120
    17. Landvik M , Wahlberg N , Roslin T (2013) The identity of the Finnish Osmoderma (Coleoptera: Scarabaeidae, Cetoniinae) population established by COI sequencing , Entomol. Fennica, Vol.8 ; pp.147-155
    18. Lee S P (2013) Red data book of endangered insects in Korea II. Red Data Book 8, Nat. Inst. Bio. Res,
    19. Lee GE , Han T , Jeong J , Kim SH , Park IG , Park H (2015) Molecular phylogeny of the genus Dicronocephalus (Coleoptera, Scarabaeidae, Cetoniinae) based on mtCOI and 16S rRNA genes , Zookeys, Vol.501 ; pp.63-87
    20. Lewis G (1887) Notes on a new species of Osmoderma and a Trichius from Japan , Wie. Entomol. Zeit, Vol.6 (2) ; pp.49
    21. Motschulsky V de (1845) Remarques sur la collection de Coléoptères russes de Victor de Motschulsky , Bull. Soc. Imp . Nat. Moscou, Vol.18 ; pp.1-127
    22. Nardi G , Micó E (2010) Osmoderma cristinae , The IUCN Red List of Threatened Species. Version 2014.3, http://www.uicnredlist.org (accessed on 15 March 2016)
    23. Nieto A , Alexander KNA (2010) European red list of saproxylic beetles, Publications Office of the European Union,
    24. Nei M Miyamoto M M , Cracraft J (1991) Relative efficiencies of different tree-making methods for molecular data , Phylogenetic Analysis of DNA Sequences, Oxford University Press, ; pp.90-128
    25. Nei M (1996) Phylogenetic analysis in molecular evolutionary genetics , Ann. Rev. Gen, Vol.30 ; pp.371-403
    26. Ranius T (2002a) Osmoderma eremita as an indicator of species richness of beetles in tree hollows , Biodiv. Conserv, Vol.11 ; pp.931-941
    27. Ranius T (2002b) Influence of stand size and quality of tree hollows on saproxylic beetles in Sweden , Bio. Conserv, Vol.103 ; pp.85-91
    28. Saitou N , Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees , Mol. Biol. Evol, Vol.4 ; pp.406-425
    29. Shin YH , Yoon IB (1994) Check List of Insects from Korea, Kunkuk University Press,
    30. Stebnicka Z (1980) Scarabaeoides of the Democratic People's Republic of Korea , Acta Zool. Cracov, Vol.24 ; pp.191-297
    31. Tamura K , Peterson D , Peterson N , Stecher G , Nei M , Kumar S (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods , Mol. Biol. Evol, Vol.28 ; pp.2731-2739
    32. Tauzin P (1994) Le genre Osmoderma Le Peletier et Audinet-Serville 1828 (Coleoptera, Cetoniidae, Trichiinae, Osmodermatini) , Systématique, biologie et distribution (Première partie). L’Entomologiste, Vol.50 ; pp.195-214
    33. Tauzin P (2006) Nouveau statut pour la variété castaneum de Osmoderma barnabita Motschoulsky, 1845 (Coleoptera, Scarabaeoidea, Cetoniidae, Trichiinae, Osmodermatini) , Cetoniimania, Vol.1-2 ; pp.39-49
    34. Tauzin P (2013) Présence de Osmoderma opicum Lewis, 1887 en Corée du Nord (Coleoptera, Scarabaeoidae, Cetoniinae, Trichiini) , Lambillionea CXIII, Vol.2 ; pp.174-178
    35. Wandeler P , Hoeck PEA , Keller LF (2007) Back to the future: museum specimens in population genetics , Trend. Ecol. Evol, Vol.22 ; pp.634-642
    36. Untergrasser A , Cutcutache I , Koressaar TYJ , Faircloth B C , Remm M , Rozen SG (2012) Primer3 - new capabilities and interfaces , Nucl. Acids Res, Vol.40 ; pp.e115