- Research article
- Open Access
Analysis of a normalised expressed sequence tag (EST) library from a key pollinator, the bumblebee Bombus terrestris
© Sadd et al; licensee BioMed Central Ltd. 2010
- Received: 30 November 2009
- Accepted: 15 February 2010
- Published: 15 February 2010
The bumblebee, Bombus terrestris (Order Hymenoptera), is of widespread importance. This species is extensively used for commercial pollination in Europe, and along with other Bombus spp. is a key member of natural pollinator assemblages. Furthermore, the species is studied in a wide variety of biological fields. The objective of this project was to create a B. terrestris EST resource that will prove to be valuable in obtaining a deeper understanding of this significant social insect.
A normalised cDNA library was constructed from the thorax and abdomen of B. terrestris workers in order to enhance the discovery of rare genes. A total of 29'428 ESTs were sequenced. Subsequent clustering resulted in 13'333 unique sequences. Of these, 58.8 percent had significant similarities to known proteins, with 54.5 percent having a "best-hit" to existing Hymenoptera sequences. Comparisons with the honeybee and other insects allowed the identification of potential candidates for gene loss, pseudogene evolution, and possible incomplete annotation in the honeybee genome. Further, given the focus of much basic research and the perceived threat of disease to natural and commercial populations, the immune system of bumblebees is a particularly relevant component. Although the library is derived from unchallenged bees, we still uncover transcription of a number of immune genes spanning the principally described insect immune pathways. Additionally, the EST library provides a resource for the discovery of genetic markers that can be used in population level studies. Indeed, initial screens identified 589 simple sequence repeats and 854 potential single nucleotide polymorphisms.
The resource that these B. terrestris ESTs represent is valuable for ongoing work. The ESTs provide direct evidence of transcriptionally active regions, but they will also facilitate further functional genomics, gene discovery and future genome annotation. These are important aspects in obtaining a greater understanding of this key pollinator species.
- Unique Sequence
- Repeat Length
- Custom Perl Script
- Population Level Study
- Slim Category
Social insects, especially those belonging to the order Hymenoptera, have been an extremely successful group. They take on a wide diversity of roles around the world, and while some are considered as pests, others are seen as vital for ecosystem functioning and even commercial prosperity. Bumblebees (the genus Bombus) fall into the latter beneficial groups, being important members of natural pollinator assemblages and providing commercial pollination services .
Bumblebees have been an object of diverse scientific study due to their importance, behaviour, social life, and a number of other fascinating traits. To name but a few, bumblebees have been utilised in research on social evolution [2, 3] and organisation [4, 5], development [6, 7], plant-pollinator interactions [8–10], learning [11, 12], invasion biology , host-parasite ecology [14–17], ecological immunology [18–20], and community ecology [21, 22]. While not officially domesticated, bumblebees still present an important agricultural resource with a significant economic load [23, 24]. Bumblebees in both Europe (mainly B. terrestris) and North America (mainly B. impatiens) are bred commercially by the hundred thousands for this reason, with the colonies being used in the pollination of a variety of standard glasshouse crops . Many have considered their economic importance to be on the rise given the recent declines in honeybee numbers . Bumblebees not only provide valuable pollination services in the agricultural sector, but the pollination they carry out is an integral part of many natural ecosystems [26–28]. The importance of bumblebees comes sharply into focus when we consider the threat that pollinators currently face worldwide. This pollinator crisis has been epitomised by the recent honeybee colony collapse disorder , but the bumblebee has also suffered population declines in many areas over a number of decades [29–31].
Given their importance, a good knowledge of bumblebees from the genomic level up is of great interest. An improved set of genomic resources for bumblebees would facilitate further studies in the bumblebee, and comparisons with other bees such as the honeybee, Apis mellifera, where the genome is already sequenced . There is currently a limited set of resources available for genetic investigation, and at the time of writing, for B. terrestris, there were 61 cDNA sequences and 288 nucleotide sequences (172 of which represent microsatellite markers) in NCBI Genbank (http://www.ncbi.nlm.nih.gov, November 2009). Further tools available include a linkage map for B. terrestris and a BAC library .
In order to greatly expand the genomic resources available for B. terrestris we constructed a normalised expressed sequence tag (EST) library from thorax and abdomen tissue of workers. ESTs are short, unedited, randomly selected, single pass reads from cDNA libraries . They provide evidence of transcriptionally active regions in an organism and are an excellent resource for upstream work including gene discovery, functional genomics and marker discovery [35, 36]. Standard cDNA libraries can have problems with redundancy, and rare transcripts are often under-represented. Therefore, we employed normalisation in order to increase the number of unique transcripts discovered .
Here we describe a collection of 29,428 ESTs, which represent a valuable resource for further ecological and evolutionary studies in Bombus spp. and for comparative studies with other Hymenoptera and insects. Sequences have been deposited in the GenBank, EMBL, and DDBJ nucleotide sequence databases under the Accession numbers FN611035 through to FN640462. For further material, including consensus sequences of assembled contigs, please contact the authors.
Sequence assembly and features
A total of 29,428 quality controlled EST sequences were assembled using TGI clustering tools (TGICL). The assembly generated 4,682 contigs and 8,651 singletons. The average number of sequences per contig was 4.43 with the maximum being 35. Contigs averaged 960 bases in length with the longest contig being 3,542 bases. In total 13,333 unique sequences were produced. A table showing all EST and contig statistics can be found as an additional file (Additional file 1).
Taxonomic distribution of the "best-hits" for the 13,333 post-assembly Bombus terrestris EST sequences (BLASTX, cutoff = 1e-05).
Number of assembled ESTs
Existing Bombus terrestris sequences
Other Bombus spp.
Gene ontology (GO) terms
Comparison of the percentage distribution of Bombus terrestris Gene Ontology terms (overview GO slim terms) based on the reported EST library and the distribution of Gene Ontology terms for all proteins predicted in each of the Apis mellifera and Tribolium castaneum genomes.
GO slim term*
translation regulator activity
signal transducer activity
transcription regulator activity
enzyme regulator activity
structural molecule activity
response to stimulus
cellular amino acid and derivative metabolic process
nucleobase, nucleoside, nucleotide and nucleic acid metabolic process
extracellular structure organization
multicellular organismal development
regulation of biological process
Comparisons to the honeybee and other insects
The GC content of the B. terrestris unique sequences was calculated to be 36 percent (omitting ambiguous bases). Gene regions in the A. mellifera genome have a GC content of 29 percent, while in D. melanogaster and A. gambiae it is 44 and 47 percent, respectively . The GC percentage based on these B. terrestris ESTs suggests that lower GC content of gene regions is not restricted to Apis, and may be a wider phenomenon. However, calculations of GC-content in Nasonia vitripennis mRNA (version dated 9th July 2007) at 46 percent hint that it is not Hymenopteran wide and further species would need to be studied to determine the exact range of this low GC-content.
The 7,216 "hits" of unique sequences to A. mellifera proteins were to 4,480 unique Apis proteins. This level of redundancy suggests that 62 percent of the assembled sequences are actually derived from unique genes, yielding an estimation of 8,278 unique genes represented by this EST library. In addition to a comparison with the predicted proteins in the A. mellifera genome, the unique sequences derived from the B. terrestris EST library were also compared to the entire A. mellifera genome sequence (BLASTN, cutoff = 1e-08) . A total of 10,024 unique sequences had "hits" to the genome, 6,749 of which already had "hits" to predicted A. mellifera proteins, but 3,275 which did not. These unique sequences with "hits" to the genome were distributed on the Apis linkage groups as follows, with the figure in brackets being the percent distribution of actual Apis protein coding sequences: 9.98 (9.57) percent on LG1, 5.89 (5.12) on LG2, 4.94 (4.61) on LG3, 4.94 (4.52) on LG4, 6.74 (5.77) on LG5, 3.73 (3.55) on LG6, 3.5 (3.56) on LG7, 5.25 (4.79) on LG8, 3.79 (3.76) on LG9, 5.04 (4.44) on LG10, 6.24 (5.97) on LG11, 3.20 (2.80) on LG12, 3.25 (3.03) on LG13, 4.03 (4.13) on LG14, 4.35 (4.61) on LG 15, 2.77 (2.30) on LG16, 0.02 (0.14) on the Mitochondria and 22.38 (27.34) on unassembled regions.
Unique sequences from the 20 longest ORFs that have "hits" to other insect proteins and also "hits" away from protein coding regions in the Apis genome, but are not part of the Apis protein set.
Project Identification (bom001no-)
Potential reason for absence from A. mellifera protein set
similar to F-box protein 28
similar to p53-like protein
similar to ENSANGP00000017887
similar to GA18228-PA
P. humanus corporis
nuclear pore complex protein nup98
Immune genes and pathways
While the bumblebee workers used to produce this EST library were not immune challenged, the vast amount of work on immunity and host-parasite interactions in these insects [19, 43–45] makes this category of genes nonetheless worthy of further interest. Furthermore, parasites have been implicated in declines of natural and commercial pollinators, such as bumblebees, worldwide [25, 46, 47], and thus knowledge of genes and pathways involved in immunity and parasite defence is highly important.
Based on the annotations of the unique sequences, and in comparison with proposed A. mellifera immune genes , we find that 134 unique sequences have "best hits" to proteins from genes characterised as immune related (Additional file 6). These 107 unique sequences in Bombus terrestris match to 67 Apis genes (a similar level of redundancy to that shown with all A. mellifera protein "hits"). In addition, based on the classification by Sackton et al., 101 unique sequences have "best-hits" in D. melanogaster to proteins of 58 immune genes. However, a core set of only 32 unique sequences had both immune classified "hits" in D. melanogaster and A. mellifera, with "hits" to 24 and 25 unique genes, respectively. This could be due differential annotation of immune genes in the two organisms, or different divergence from B. terrestris among different genes.
Potential EST-derived genetic markers
Traditional methods of isolating and characterising genetic markers can be expensive and time-consuming. However, ESTs can be a valuable resource for the initial identification of genetic markers that can subsequently be used in studies of molecular ecology . Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) allow population level studies both within the target species and related species.
Using MISA , a number of SSRs were identified in the unique sequences. We restricted our search for SSRs to di-, tri-, tetra-, penta-, and hexa- nucleotide motifs that were repeated at least six times for di-nucleotide and five times for all other motifs. We found a total of 589 markers found in 507 unique sequences. 397 di- (mean repeat length = 7.65, range = 6-17), 154 tri- (mean repeat length = 5.79, range = 5-17), 26 tetra- (mean repeat length = 6, range = 5-11), 6 penta- (mean repeat length = 5.83, range = 5-7), and 6 hexa- (mean repeat length = 6.33, range = 5-9) nucleotide repeats were found. 61 of the SSRs were found in compound formation with another (maximum of 100 nucleotides separation). While polymorphism remains to be verified for these markers with appropriate laboratory screens, previous levels of polymorphism in EST derived SSRs suggest that a large number of these markers will be variable .
These identified potential SSRs together with existing markers from B. terrestris[52–54] make up a significant resource for molecular ecology studies. While EST derived SSRs tend to show a lower level of polymorphism than those identified from genomic libraries, they can more often be used across related species [36, 51, 55]. Given the pervasive importance of species of the genus Bombus in temperate pollinator communities, this is a highly advantageous feature. In addition to SSRs, contigs built from a number of EST reads can be useful sources of potential SNPs. While the normalization procedure used in this work reduces such redundancy and thus the likelihood of finding SNPs, 1,973 contigs contain at least four sequences, a minimum for reliable identification of SNPs . These contigs cover 2,124,958 bases. Using SNPidentifier  with minimum minor allele frequency of 0.1 and a minimum occurrence of 2, we were able to find 854 potential SNPs occuring in 573 contigs. These potential SNPs provide an additional set of genetic markers for subsequent population studies.
The resource provided by the presented EST library will serve researchers in a diverse range of biological fields. The detected markers will facilitate further population level studies, and the sequences themselves can be used in investigations of particular genes, and also in functional genomic approaches. Further, as demonstrated here, the sequences allow genomic comparisons with the honeybee A. mellifera and other insects, informing on such things as potential gene loss or pseudogene evolution. This EST library will also be a practical resource in the assembly and annotation of the impending B. terrestris genome http://www.hgsc.bcm.tmc.edu. These elements, coupled with the major commercial and ecosystem functions of B. terrestris and related hymenoptera, mean that this EST library is a useful resource for ongoing research.
Colonies were set-up with field caught queens collected from Northern Switzerland in the spring of 2007 (Aesch, Basel), and were kept under red-light illumination at 27 ± 1°C with sugar water (ApiInvert®) and pollen provided ad libitum. Worker bees used for the EST library originated from four colonies (two workers each) that had been confirmed as Bombus terrestris by visual checks and based on a species-specific length polymorphism marker in the region between mitochondrial CO1 and CO2 genes. The colonies were also confirmed to be free of common parasites by microscopic investigation of faeces. Worker bees were age controlled, and one bee from each colony was 7 and one bee 14 days post adult eclosion at the time of sacrifice. Bees were sacrificed by snap freezing in liquid nitrogen, and were subsequently stored at -80°C. Tissue from both the thorax and abdomen was used as material for the EST library.
RNA isolation and cDNA library construction
Total RNA was extracted with Solution D . cDNA was synthesized using the SMART approach (Mint-Universal cDNA synthesis kit, Evrogen, Russia), subsequently normalized using duplex-specific nuclease (Trimmer kit, Evrogen, Russia) according to manufacturer's instructions, and directionally cloned into pAL32 plasmid vector (Evrogen, Russia). Plasmids were transferred via electroporation into E. coli DH10B (Invitrogen, U.S.A.).
EST sequencing, quality control and assembly
Plasmids from the normalized libraries were 5' end sequenced using the pALforward primer (5'-CTCGGGAAGCGCGCCATT-3') and Big Dye Terminator chemistry (ABI). Collected reads from 3730XL capillary sequencers (ABI) were base-called using Phred http://www.phrap.org and subsequently quality and vector clipped using Lucy http://compbio.dfci.harvard.edu/tgi/software with standard parameters. For clustering and assembly, the TGI tools developed at TIGR http://compbio.dfci.harvard.edu/tgi/software were used.
Annotation and comparisons
Once clusters were created, the resulting unique sequences were used for various comparisons. BLASTX searches with a cut-off of 1e-05 were run against each of the nr database in Genbank (version dated 19th October 2009), and databases comprising of all predicted proteins in the genomes of A. mellifera, N. vitripennis, D. melanogaster and T. castaneum.
Redundancy in the unique sequences was determined by comparison with A. mellifera. This was calculated as the number of unique A. mellifera proteins that were "best-hits" in the BLASTX search divided by the number of unique B. terrestris sequences with hits. GC-content was calculated using a custom Perl script that amalgamated all unique sequences and calculated the proportion of G and C bases in all unambiguous bases.
Annotation of GO terms  to the current B. terrestris EST derived unique sequences were carried out using the Perl based annot8r . A database containing SwissProt and Trembl (Uniprot Knowledgebase Release 15.9) with manually curated GO terms was used. Terms that had been derived electronic annotation were omitted to improve reliability of the annotations. BLASTP against this database was carried out with a cut-off of 1e-05. In addition, predicted proteins from the entire genomes of A. mellifera and T. castaneum were compared with the same database for comparisons of the distribution of GO terms. For statistical comparisons, GO terms were matched to higher-level GO Slim categories. Pairwise comparisons between B. terrestris and both A. mellifera and T. castaneum were then carried out for the distribution of sequences in each GO Slim category using Fisher's exact tests in R2.92 for Mac .
The EST unique sequences were compared with the A. mellifera genome  using BLASTN, cut-off 1e-08. A custom Perl script was then used to extract positional information of all the "best-hits". Linkage group and position on the linkage group was then compared with the position of predicted protein coding regions within the genome to obtain a relative position of the "hits" (inside a protein-coding region; spanning a protein-coding region boundary; outside a protein coding region). For those outside a predicted protein-coding region, the distance to the closest was calculated. For those B. terrestris sequences with hits greater than 4'375 bases (the largest A. mellifera UTR, obtained from the UCSC genome browser ) away from a protein coding region, a custom Perl script was used to interface with ORF Finder http://www.ncbi.nlm.nih.gov/gorf to obtain ORF information. Of these sequences, the 20 with the largest ORFs were investigated manually. BLASTP searches were carried out against the nr database, and if the "best-hit" was to another insect protein (cut-off = 1e-20), the "best-hit" was taken and blasted (TBLASTN) against the A. mellifera genome. It was confirmed that the "hit" location was the same as for the EST derived unique sequence and the Blast alignments were manually examined for potential features inhibiting protein production.
SSR and SNP marker identification
With MISA the configuration file was adjusted to identify motifs of two, three, four, five or six nucleotides, that repeat at least six times for di-nucleotide repeats and five times for all others. MISA was then run on all the post-assembly unique sequences. With SNP identifier we adjusted the minimum minor allele frequency to 0.1 and occurrence to 2. This meant that post-assembly sequences needed to be made up of at least four ESTs to be considered. SNPidentifier was run on the alignments of the EST sequences.
We acknowledge Monique Roenick and Janina Thiel (Berlin) for technical assistance, and Alfred Beck (Berlin) for computational support. This study was supported by the Max Planck Society, the Swiss National Science Foundation, CCES (GEDIHAP) and the Genetic Diversity Centre of ETH Zurich.
- Goulson D: Bumblebees - their behaviour and ecology. 2003, New York, USA: Oxford University PressGoogle Scholar
- Pamilo P: Evolution of Colony Characteristics in Social Insects .1. Sex Allocation. Am Nat. 1991, 137 (1): 83-107. 10.1086/285147.View ArticleGoogle Scholar
- Page RE, Metcalf RA: Multiple Mating, Sperm Utilization, and Social Evolution. Am Nat. 1982, 119 (2): 263-281. 10.1086/283907.View ArticleGoogle Scholar
- Jandt JM, Dornhaus A: Spatial organization and division of labour in the bumblebee Bombus impatiens. Anim Behav. 2009, 77 (3): 641-651. 10.1016/j.anbehav.2008.11.019.View ArticleGoogle Scholar
- Schmid-Hempel P, Schmid-Hempel R: Transmission of a pathogen in Bombus terrestris, with a note on division of labor in social insects. Behavioral Ecology and Sociobiology. 1993, 33 (5): 319-327. 10.1007/BF00172930.View ArticleGoogle Scholar
- Yerushalmi S, Bodenhaimer S, Bloch G: Developmentally determined attenuation in circadian rhythms links chronobiology to social organization in bees. J Exp Biol. 2006, 209 (6): 1044-1051. 10.1242/jeb.02125.PubMedView ArticleGoogle Scholar
- Hartfelder K, Cnaani J, Hefetz A: Caste-specific differences in ecdysteroid titers in early larval stages of the bumblebee Bombus terrestris. J Insect Physiol. 2000, 46 (11): 1433-1439. 10.1016/S0022-1910(00)00067-6.PubMedView ArticleGoogle Scholar
- Cooley AM, Carvallo G, Willis JH: Is floral diversification associated with pollinator divergence? Flower shape, flower colour and pollinator preference in chilean Mimulus. Annals of Botany. 2008, 101 (5): 641-650. 10.1093/aob/mcn014.PubMed CentralPubMedView ArticleGoogle Scholar
- Suzuki K, Dohzono I, Hiei K: Evolution of pollinator generalization in bumblebee-pollinated plants. Plant Species Biology. 2007, 22 (3): 141-159. 10.1111/j.1442-1984.2007.00187.x.View ArticleGoogle Scholar
- Gegear RJ, Burns JG: The birds, the bees, and the virtual flowers: Can pollinator Behavior drive ecological speciation in flowering plants?. Am Nat. 2007, 170 (4): 551-566. 10.1086/521230.PubMedView ArticleGoogle Scholar
- Ings TC, Raine NE, Chittka L: A population comparison of the strength and persistence of innate colour preference and learning speed in the bumblebee Bombus terrestris. Behavioral Ecology and Sociobiology. 2009, 63 (8): 1207-1218. 10.1007/s00265-009-0731-8.View ArticleGoogle Scholar
- Leadbeater E, Chittka L: The dynamics of social learning in an insect model, the bumblebee (Bombus terrestris). Behavioral Ecology and Sociobiology. 2007, 61 (11): 1789-1796. 10.1007/s00265-007-0412-4.View ArticleGoogle Scholar
- Schmid-Hempel P, Schmid-Hempel R, Brunner PC, Seeman OD, Allen GR: Invasion success of the bumblebee, Bombus terrestris, despite a drastic genetic bottleneck. Heredity. 2007, 99 (4): 414-422. 10.1038/sj.hdy.6801017.PubMedView ArticleGoogle Scholar
- Otti O, Schmid-Hempel P: Nosema bombi: A pollinator parasite with detrimental fitness effects. J Invertebr Pathol. 2007, 96 (2): 118-124. 10.1016/j.jip.2007.03.016.PubMedView ArticleGoogle Scholar
- Baer B, Schmid-Hempel P: Bumblebee workers from different sire groups vary in susceptibility to parasite infection. Ecol Lett. 2003, 6 (2): 106-110. 10.1046/j.1461-0248.2003.00411.x.View ArticleGoogle Scholar
- Durrer S, Schmid-Hempel P: Shared Use of Flowers Leads to Horizontal Pathogen Transmission. Proc R Soc Lond Ser B-Biol Sci. 1994, 258 (1353): 299-302. 10.1098/rspb.1994.0176.View ArticleGoogle Scholar
- Schmid-Hempel P, Reber Funk C: The distribution of genotypes of the trypanosome parasite, Crithidia bombi, in populations of its host, Bombus terrestris. Parasitology. 2004, 129: 147-158. 10.1017/S0031182004005542.PubMedView ArticleGoogle Scholar
- Doums C, Moret Y, Benelli E, Schmid-Hempel P: Senescence of immune defence in Bombus workers. Ecol Entomol. 2002, 27 (2): 138-144. 10.1046/j.1365-2311.2002.00388.x.View ArticleGoogle Scholar
- Moret Y, Schmid-Hempel P: Survival for immunity: The price of immune system activation for bumblebee workers. Science. 2000, 290 (5494): 1166-1168. 10.1126/science.290.5494.1166.PubMedView ArticleGoogle Scholar
- Sadd BM, Schmid-Hempel P: Facultative but persistent transgenerational immunity via the mother's eggs in bumblebees. Current Biology. 2007, 17 (24): R1046-R1047. 10.1016/j.cub.2007.11.007.PubMedView ArticleGoogle Scholar
- Hatfield RG, LeBuhn G: Patch and landscape factors shape community assemblage of bumble bees, Bombus spp. (Hymenoptera: Apidae), in montane meadows. Biological Conservation. 2007, 139 (1-2): 150-158. 10.1016/j.biocon.2007.06.019.View ArticleGoogle Scholar
- Williams PH, Araujo MB, Rasmont P: Can vulnerability among British bumblebee (Bombus) species be explained by niche position and breadth?. Biological Conservation. 2007, 138 (3-4): 493-505. 10.1016/j.biocon.2007.06.001.View ArticleGoogle Scholar
- Velthuis HHW, van Doorn A: A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie. 2006, 37 (4): 421-451. 10.1051/apido:2006019.View ArticleGoogle Scholar
- Palma G, Quezada-Euan JJG, Reyes-Oregel V, Melendez V, Moo-Valle H: Production of greenhouse tomatoes (Lycopersicon esculentum) using Nannotrigona perilampoides, Bombus impatiens and mechanical vibration (Hym.: Apoidea). Journal of Applied Entomology. 2008, 132 (1): 79-85. 10.1111/j.1439-0418.2007.01246.x.View ArticleGoogle Scholar
- Cox-Foster D, Vanengelsdorp D: Saving the HONEYBEE. Scientific American. 2009, 300 (4): 40-10.1038/scientificamerican0409-40.PubMedView ArticleGoogle Scholar
- Bingham RA, Ort AR: Efficient pollination of alpine plants. Nature. 1998, 391 (6664): 238-239. 10.1038/34564.View ArticleGoogle Scholar
- Ono A, Dohzono I, Sugawara T: Bumblebee pollination and reproductive biology of Rhododendron semibarbatum (Ericaceae). Journal of Plant Research. 2008, 121 (3): 319-327. 10.1007/s10265-008-0155-y.PubMedView ArticleGoogle Scholar
- Forup ML, Henson KSE, Craze PG, Memmott J: The restoration of ecological interactions: plant-pollinator networks on ancient and restored heathlands. Journal of Applied Ecology. 2008, 45 (3): 742-752. 10.1111/j.1365-2664.2007.01390.x.View ArticleGoogle Scholar
- Biesmeijer JC, Roberts SPM, Reemer M, Ohlemuller R, Edwards M, Peeters T, Schaffers AP, Potts SG, Kleukers R, Thomas CD: Parallel declines in pollinators and insect-pollinated plants in Britain and the Netherlands. Science. 2006, 313 (5785): 351-354. 10.1126/science.1127863.PubMedView ArticleGoogle Scholar
- Grixti JC, Wong LT, Cameron SA, Favret C: Decline of bumble bees (Bombus) in the North American Midwest. Biological Conservation. 2009, 142 (1): 75-84. 10.1016/j.biocon.2008.09.027.View ArticleGoogle Scholar
- Fitzpatrick U, Murray TE, Paxton RJ, Breen J, Cotton D, Santorum V, Brown MJF: Rarity and decline in bumblebees - A test of causes and correlates in the Irish fauna. Biological Conservation. 2007, 136 (2): 185-194. 10.1016/j.biocon.2006.11.012.View ArticleGoogle Scholar
- Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443 (7114): 931-949. 10.1038/nature05260.View ArticleGoogle Scholar
- Wilfert L, Gadau J, Schmid-Hempel P: A core linkage map of the bumblebee Bombus terrestris. Genome. 2006, 49 (10): 1215-1226. 10.1139/G06-075.PubMedView ArticleGoogle Scholar
- Wilfert L, Torres MM, Reber-Funk C, Schmid-Hempel R, Tomkins J, Gadau J, Schmid-Hempel P: Construction and characterization of a BAC-library for a key pollinator, the bumblebee Bombus terrestris L. Insectes Sociaux. 2009, 56 (1): 44-48. 10.1007/s00040-008-1034-1.View ArticleGoogle Scholar
- Nagaraj SH, Gasser RB, Ranganathan S: A hitchhiker's guide to expressed sequence tag (EST) analysis. Briefings in Bioinformatics. 2007, 8 (1): 6-21. 10.1093/bib/bbl015.PubMedView ArticleGoogle Scholar
- Bouck A, Vision T: The molecular ecologist's guide to expressed sequence tags. Molecular Ecology. 2007, 16 (5): 907-924. 10.1111/j.1365-294X.2006.03195.x.PubMedView ArticleGoogle Scholar
- Bonaldo MDF, Lennon G, Soares MB: Normalization and subtraction: Two approaches to facilitate gene discovery. Genome Research. 1996, 6 (9): 791-806. 10.1101/gr.6.9.791.PubMedView ArticleGoogle Scholar
- Schmid R, Blaxter M: annot8r: GO, EC and KEGG annotation of EST datasets. BMC Bioinformatics. 2008, 9 (1): 180-10.1186/1471-2105-9-180.PubMed CentralPubMedView ArticleGoogle Scholar
- Tribolium Genome Sequencing Consortium: The genome of the model beetle and pest Tribolium castaneum. Nature. 2008, 452 (7190): 949-955. 10.1038/nature06784.View ArticleGoogle Scholar
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Research. 2006, 16 (11): 1334-1338. 10.1101/gr.5204306.PubMed CentralPubMedView ArticleGoogle Scholar
- Sharkey MJ: Phylogeny and classification of Hymenoptera. Zootaxa. 2007, 521-548. 1668Google Scholar
- Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, Kent WJ: The UCSC Table Browser data retrieval tool. Nucleic Acids Res. 2004, 32: D493-D496. 10.1093/nar/gkh103.PubMed CentralPubMedView ArticleGoogle Scholar
- Baer B, Schmid-Hempel P: Experimental variation in polyandry affects parasite loads and fitness in a bumble-bee. Nature. 1999, 397 (6715): 151-154. 10.1038/16451.View ArticleGoogle Scholar
- Schmid-Hempel P: On the evolutionary ecology of host-parasite interactions: addressing the question with regard to bumblebees and their parasites. Naturwissenschaften. 2001, 88 (4): 147-158. 10.1007/s001140100222.PubMedView ArticleGoogle Scholar
- Sadd BM, Schmid-Hempel P: Insect immunity shows specificity in protection upon secondary pathogen exposure. Current Biology. 2006, 16 (12): 1206-1210. 10.1016/j.cub.2006.04.047.PubMedView ArticleGoogle Scholar
- Colla SR, Otterstatter MC, Gegear RJ, Thomson JD: Plight of the bumble bee: Pathogen spillover from commercial to wild populations. Biological Conservation. 2006, 129 (4): 461-467. 10.1016/j.biocon.2005.11.013.View ArticleGoogle Scholar
- Thorp RW: Species Profile: Bombus franklini. Red list of pollinator insects of North America. Edited by: Shepherd M, Vaughan, DM, Black SH. 2005, Portland, Oregon, USA: The Xerces Society for Invertebrate ConservationGoogle Scholar
- Evans JD, Aronstein K, Chen YP, Hetru C, Imler JL, Jiang H, Kanost M, Thompson GJ, Zou Z, Hultmark D: Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol Biol. 2006, 15 (5): 645-656. 10.1111/j.1365-2583.2006.00682.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39 (12): 1461-1468. 10.1038/ng.2007.60.PubMedView ArticleGoogle Scholar
- MISA - MIcroSAtellite identification tool. [http://pgrc.ipk-gatersleben.de/misa/]
- Ellis JR, Burke JM: EST-SSRs as a resource for population genetic analyses. Heredity. 2007, 99 (2): 125-132. 10.1038/sj.hdy.6801001.PubMedView ArticleGoogle Scholar
- Funk CR, Schmid-Hempel R, Schmid-Hempel P: Microsatellite loci for Bombus spp. Molecular Ecology Notes. 2006, 6 (1): 83-86. 10.1111/j.1471-8286.2005.01147.x.View ArticleGoogle Scholar
- Stolle E, Rohde M, Vautrin D, Solignac M, Schmid-Hempel P, Schmid-Hempel R, Moritz RFA: Novel microsatellite DNA loci for Bombus terrestris (Linnaeus, 1758). Molecular Ecology Resources. 2009, 9 (5): 1345-1352. 10.1111/j.1755-0998.2009.02610.x.PubMedView ArticleGoogle Scholar
- Estoup A, Solignac M, Cornuet JM, Goudet J, Scholl A: Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe. Molecular Ecology. 1996, 5 (1): 19-31. 10.1111/j.1365-294X.1996.tb00288.x.PubMedView ArticleGoogle Scholar
- Pashley CH, Ellis JR, McCauley DE, Burke JM: EST databases as a source for molecular markers: Lessons from Helianthus. Journal of Heredity. 2006, 97 (4): 381-388. 10.1093/jhered/esl013.PubMedView ArticleGoogle Scholar
- Long AD, Beldade P, Macdonald SJ: Estimation of population heterozygosity and library construction-induced mutation rate from expressed sequence tag collections. Genetics. 2007, 176 (1): 711-714. 10.1534/genetics.106.063610.PubMed CentralPubMedView ArticleGoogle Scholar
- Gorbach DM, Hu ZL, Du ZQ, Rothschild MF: SNP discovery in Litopenaeus vannamei with a new computational pipeline. Animal Genetics. 2009, 40 (1): 106-109. 10.1111/j.1365-2052.2008.01792.x.PubMedView ArticleGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA Isolation by acid guanidinium thiocyanate phenol chloroform extraction. Anal Biochem. 1987, 162 (1): 156-159. 10.1016/0003-2697(87)90021-2.PubMedView ArticleGoogle Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JGR, Korf I, Lapp H: The Bioperl Toolkit: Perl Modules for the Life Sciences. Genome Research. 2002, 12 (10): 1611-1618. 10.1101/gr.361602.PubMed CentralPubMedView ArticleGoogle Scholar
- R: A language and environment for statistical computing. [http://www.R-project.org]
- Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.