Research note | Open | Published:
Identification of 24 new microsatellite loci in the sweat bee Lasioglossum malachurum (Hymenoptera: Halictidae)
BMC Research Notesvolume 10, Article number: 753 (2017)
The objective here is to identify highly polymorphic microsatellite loci for the Palaearctic sweat bee Lasioglossum malachurum. Sweat bees (Hymenoptera: Halictidae) are widespread pollinators that exhibit an unusually large range of social behaviours from non-social, where each female nests alone, to eusocial, where a single queen reproduces while the other members of the colony help to rear her offspring. They thus represent excellent models for understanding social evolution.
24 new microsatellite loci were successfully optimized. When amplified across 23–40 unrelated females, the number of alleles per locus ranged from 3 to 17 and the observed heterozygosities 0.45 to 0.95. Only one locus showed evidence of significant deviation from Hardy–Weinberg equilibrium. No evidence of linkage disequilibrium was found. These 24 loci will enable researchers to gain greater understanding of colony relationships within this species, an important model for the study of eusociality. Furthermore, 22 of the same loci were also successfully amplified in L. calceatum, suggesting that these loci may be useful for investigating the ecology and evolution of sweat bees in general.
Sweat bees (Hymenoptera: Halictidae) are widespread pollinators which exhibit an unusually large range of social behaviours from non-social, where each female nests alone, to eusocial, where a single queen reproduces while the other members of the colony help to rear her offspring . Sweat bees are also unusual in that social and non-social species are often closely related, with multiple evolutionary transitions having occurred between sociality and non-sociality . Sweat bees thus represent excellent models for understanding social evolution [1, 2]. Here we present a new set of microsatellite loci developed from Lasioglossum malachurum (Kirby, 1802), a haplodiploid eusocial species that has been particularly well studied, mainly because it is widely distributed in the Western Palaearctic and because it often occurs in large, dense nesting aggregations that facilitate behavioural research [3,4,5,6]. Microsatellite markers are widely used in social evolution research, for example to investigate population structure, estimate genetic relatedness and assign offspring to parents [7,8,9]. Microsatellite loci have been developed for this species previously [3, 10] but most of them have comparatively low heterozygosities and are difficult to combine into multiplex reactions because of highly specific annealing temperatures and polymerase chain reaction (PCR) mixes. Here, we report 24 new microsatellite markers developed for L. malachurum, 14 of which have been efficiently amplified in two multiplex sets. These markers should substantially aid future studies on sweat bee behaviour and ecology.
Lasioglossum malachurum females were sampled from a field site at Denton in East Sussex, UK in 2015. Genomic DNA was extracted from head, abdomen and/or legs using an ammonium acetate extraction method [11, 12]. DNA concentration was quantified using a Fluostar Optima fluorimeter and its quality assessed using gel electrophoresis. DNA from one foundress (female M4) from Denton was digested using MboI and the fragments enriched for dinucleotide and tetranucleotide repeat motifs (following ). An Illumina paired end library was then compiled using this repeat-enriched genomic DNA. The NEBNext Ultra library preparation kit (New England Biolabs Inc. Cat. No. E7370S) protocol was followed and DNA sequencing was conducted using a MiSeq Benchtop Sequencer (Illumina). Primer sets were designed from 53 microsatellite sequences using PRIMER3 v0.4.0 . Sequences were confirmed to be unique using BLAST software .
Each 2 µl PCR contained approximately 10 ng of air-dried genomic DNA, 0.2 µM of each primer and 1 µl QIAGEN Multiplex PCR mix (QIAGEN Inc. Cat. No. 20614) following . As we required loci that could be reliably multiplexed together for efficient use we designed primers with very similar melting temperatures (± 2 °C) enabling these to be amplified at the same annealing temperature (57 °C). The following PCR profile was used: 95 °C for 15 min, followed by 44 cycles of 94 °C for 30 s, 57 °C for 90 s, 72 °C for 90 s and finally 60 °C for 30 min. PCR amplification was performed using a DNA Engine Tetrad ®Thermal Cycler (MJ Research, Bio-Rad, Hemel Hempstead, Herts, UK). PCR products were genotyped on an ABI 3730 48-well capillary DNA Analyser using the LIZ size standard (Applied Biosystems Inc. Cat. No. 4322682). Alleles were scored using GENEMAPPERv3.7 software (Applied Biosystems Inc.). Of the 53 markers, 24 could be scored reliably across the test sample (23–40 females all from the same field site at Denton) (Table 1). The remaining 29 were found to be either monomorphic or unreliable following our PCR methodology (Table 2). It is possible that with more specific optimization, some of these could be used in future studies. We successfully incorporated 14 of the optimized markers into two multiplex panels (using the above PCR reagents and concentrations) with no dropout or artifacts produced (Table 1).
The numbers of alleles and heterozygosities were calculated for each of the 24 loci using CERVUS v3.0.6 and with the sample sizes shown in Table 1 . Tests for deviation from Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were conducted using GENEPOP web version 4.2 . To correct p-values in multiple tests, the Q Value was applied to LD p-values. The q value is a measure of the significance in terms of false discovery rate, rather than conventional Bonferroni correction which attempts to measure significance in terms of false positives only . Observed levels of heterozygosity ranged from 0.45 to 0.95 with 3–17 alleles per locus (Table 1). Only Lma31 deviated from HWE (p = 0.049). No groups of loci displayed LD, providing no evidence of physical linkage based on the individuals genotyped.
These loci are likely to be useful for investigating the ecology and behaviour of L. malachurum and also potentially that of other sweat bees. Indeed, we have successfully amplified 22 of the 24 loci in L. calceatum (Scopoli) individuals sampled in the UK; only Lma20 and Lma21 failed to amplify and 17 of the 22 loci that did amplify were polymorphic (Table 1; Davison & Field, in prep.).
Due to the relatively short read length of the MiSeq Benchtop Sequencing system we were unable to design primer sets to amplify greater than 300 bases. This may however be somewhat fortuitous; the incorporation of larger markers into multiplex panels often proves problematic, since they are generally harder to amplify than markers with smaller products and are more susceptible to dropout .
polymerase chain reaction
Schwarz MP, Richards MH, Danforth BN. Changing paradigms in insect social evolution: insights from halictine and allodapine bees. Annu Rev Entomol. 2007;52:127–50.
Danforth BN. Evolution of sociality in a primitively eusocial lineage of bees. Proc Natl Acad Sci USA. 2002;99:286–90.
Paxton RJ, Arevalo E, Field J. Microsatellite loci for the eusocial Lasioglossum malachurum and other sweat bees (Hymenoptera, Halictidae). Mol Ecol Notes. 2003;3:82–4.
Strohm E, Bordon-Hauser A. Advantages and disadvantages of large colony size in a halictid bee: the queen’s perspective. Behav Ecol. 2003;14:546–53.
Wyman LM, Richards MH. Colony social organization of Lasioglossum malachurum Kirby (Hymenoptera, Halictidae) in southern Greece. Insect Soc. 2003;50:201–11.
Zobel MU, Paxton RJ. Is big the best? Queen size, usurpation and nest closure in a primitively eusocial sweat bee (Lasioglossum malachurum). Behav Ecol Sociobiol. 2007;61:435–47.
Paxton RJ, Ayasse M, Field J, Soro A. Complex sociogenetic organization and reproductive skew in a primitively eusocial sweat bee, Lasioglossum malachurum, as revealed by microsatellites. Mol Ecol. 2002;11:2405–16.
Field J, Paxton RJ, Soro A, Bridge C. Cryptic plasticity underlies a major evolutionary transition. Curr Biol. 2010;20:2028–31.
Soro A, Field J, Bridge C, Cardinal SC, Paxton RJ. Genetic differentiation across the social transition in a socially polymorphic sweat bee, Halictus rubicundus. Mol Ecol. 2010;19:3351–63.
Soro A, Ayasse M, Zobel MU, Paxton RJ. Complex sociogenetic organization and the origin of unrelated workers in a eusocial sweat bee, Lasioglossum malachurum. Insect Soc. 2009;56:55–63.
Nicholls JA, Double MC, Rowell DM, Magrath RD. The evolution of cooperative and pair breeding in thornbills Acanthiza (Pardalotidae). J Avian Biol. 2000;31:165–76.
Richardson DS, Jury FL, Blaakmeer K, Komdeur J, Burke T. Parentage assignment and extra-group paternity in a cooperative breeder: the Seychelles warbler (Acrocephalus sechellensis). Mol Ecol. 2001;10:2263–73.
Armour JAL, Neumann R, Gobert S, Jeffreys AJ. Isolation of human simple repeat loci by hybridization selection. Hum Mol Gen. 1994;3:599–605.
Rozen S, Skaletsky HJ. PRIMER version 3 on the WWW for general users and for biologist programmers. In: Krawetz S, Misener S, editors. Bioinformatics methods and protocols: methods in molecular biology. New Jersey: Humana Press; 2000. p. 365–86.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and psi–BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Kenta T, Gratten J, Haigh NS, Hinten GN, Slate J, Butlin RK, et al. Multiplex SNP-SCALE:a cost-effective medium-throughput single nucleotide polymorphism genotyping method. Mol Ecol Res. 2008;8:1230–8.
Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16:1099–106.
Rousset F. Genepop’007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Res. 2008;8:103–10.
Storey J. A direct approach to false discovery rates. J R Stat Soc Series B Stat Methodol. 2002;64:479–98.
Toouli CD, Turner DR, Grist SA, Morley AA. The effect of cycle number and target size on polymerase chain reaction amplification of polymorphic repetitive sequences. Anal Biochem. 2000;280:324–6.
PP: Tested and optimised the primers, scored loci, performed the analysis, co-wrote the paper. CC: Collected samples, scored loci and assisted with the analysis. GJH: prepared the MiSeq library and designed primer sets. DAD: Assisted with marker development, primer design, analysis and manuscript preparation. JF: Wrote the grant application, collected samples and co-wrote the paper. All authors read and approved the final manuscript.
We thank Paul Davison for assisting with sample collection and the NERC Biomolecular Analysis Facility—Sheffield, UK for their facilities.
The authors declare that they have no competing interests.
Availability of data and materials
The sequences acquired herein have been submitted to Genbank, Accession Numbers MG273262–MG273314.
Consent for publication
Ethics approval and consent to participate
This work was funded by the UK Natural Environment Research Council (NERC Grant NE/M003191/1 to J.F.). This project has also received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant Agreement No. 695744 to JF).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.