RNA viruses of Crithidia bombi, a parasite of bumblebees.

Leishbuviridae (Bunyavirales) are a diverse monophyletic group of negative-sense single-stranded RNA virus infecting parasitic flagellates of the family Trypanosomatidae. The presence of RNA viruses in trypanosomatids can influence the virulence of the latter. Here, we performed a screening for viruses in Crithidia bombi - a common parasite of important pollinators Bombus spp. (bumblebees) that negatively affects its host in stressful conditions. The majority (8/10) of C. bombi isolates collected in Europe and North America were positive for a virus that we named Crithidia bombi leishbuvirus 1 with high conservation of amino acid sequences between isolates. The results of our comparative phylogenetic analyses of the trypanosomatids and their viruses suggest that the high mobility of bumblebees and frequent coinfections by different strains of C. bombi determine an extensive viral exchange between the latter.

Impact of climate change on parasite infection of an important pollinator depends on host genotypes.

Climate change is predicted to affect host-parasite interactions, and for some hosts, parasite infection is expected to increase with rising temperatures. Global population declines of important pollinators already have been attributed to climate change and parasitism. However, the role of climate in driving parasite infection and the genetic basis for pollinator hosts to respond often remain obscure. Based on decade-long field data, we investigated the association between climate and Nosema bombi (Microsporidia) infection of buffed-tailed bumblebees (Bombus terrestris), and whether host genotypes play a role. For this, we genotyped 876 wild bumblebee queens and screened for N. bombi infection of those queens between 2000 and 2010. We recorded seven climate parameters during those 11 years and tested for correlations between climate and infection prevalence. Here we show that climatic factors drive N. bombi infection and that the impact of climate depends on mitochondrial DNA cytochrome oxidase I (COI) haplotypes of the host. Infection prevalence was correlated with climatic variables during the time when queens emerge from hibernation. Remarkably, COI haplotypes best predict this association between climatic factors and infection. In particular, two host haplotypes ("A" and "B") displayed phenotypic plasticity in response to climatic variation: Temperature was positively correlated with infection of host haplotype B, but not haplotype A. The likelihood of infection of haplotype A was associated with moisture, conferring greater resistance to parasite infection during wetter years. In contrast, infection of haplotype B was unrelated to moisture. To the best of our knowledge, this is the first study that identifies specific host genotypes that confer differential parasite resistance under variable climatic conditions. Our results underscore the importance of mitochondrial haplotypes to ward off parasites in a changing climate. More broadly, this also suggests that COI may play a pertinent role in climate change adaptations of insect pollinators.

Genetic variation and microbiota in bumble bees cross-infected by different strains of C. bombi.

The bumblebee Bombus terrestris is commonly infected by a trypanosomatid gut parasite Crithidia bombi. This system shows a striking degree of genetic specificity where host genotypes are susceptible to different genotypes of parasite. To a degree, variation in host gene expression underlies these differences, however, the effects of standing genetic variation has not yet been explored. Here we report on an extensive experiment where workers of twenty colonies of B. terrestris were each infected by one of twenty strains of C. bombi. To elucidate the host's genetic bases of susceptibility to infection (measured as infection intensity), we used a low-coverage (~2 x) genome-wide association study (GWAS), based on angsd, and a standard high-coverage (~15x) GWAS (with a reduced set from a 8 x 8 interaction matrix, selected from the full set of twenty). The results from the low-coverage approach remained ambiguous. The high-coverage approach suggested potentially relevant genetic variation in cell surface and adhesion processes. In particular, mucin, a surface mucoglycoprotein, potentially affecting parasite binding to the host gut epithelia, emerged as a candidate. Sequencing the gut microbial community of the same bees showed that the abundance of bacterial taxa, such as Gilliamella, Snodgrassella, or Lactobacillus, differed between 'susceptible' and 'resistant' microbiota, in line with earlier studies. Our study suggests that the constitutive microbiota and binding processes at the cell surface are candidates to affect infection intensity after the first response (captured by gene expression) has run its course. We also note that a low-coverage approach may not be powerful enough to analyse such complex traits. Furthermore, testing large interactions matrices (as with the full 20 x 20 combinations) for the effect of interaction terms on infection intensity seems to blur the specific host x parasite interaction effects, likely because the outcome of an infection is a highly non-linear process dominated by variation in individually different pathways of host defence (immune) responses.

Genomic Variation among Strains of Crithidia bombi and C. expoeki.

In this study, we sequenced and analyzed the genomes of 40 strains, in addition to the already-reported two type strains, of two Crithidia species infecting bumblebees in Alaska and Central Europe and demonstrated that different strains of Crithidia bombi and C. expoeki vary considerably in terms of single nucleotide polymorphisms and gene copy number. Based on the genomic structure, phylogenetic analyses, and the pattern of copy number variation, we confirmed the status of C. expoeki as a separate species. The Alaskan populations appear to be clearly separated from those of Central Europe. This pattern fits a scenario of rapid host-parasite coevolution, where the selective advantage of a given parasite strain is only temporary. This study provides helpful insights into possible scenarios of selection and diversification of trypanosomatid parasites.IMPORTANCE A group of trypanosomatid flagellates includes several well-studied medically and economically important parasites of vertebrates and plants. Nevertheless, the vast majority of trypanosomatids infect only insects (mostly flies and true bugs) and, because of that, has attracted little research attention in the past. Of several hundred trypanosomatid species, only four can infect bees (honeybees and bumblebees). Because of such scarcity, these parasites are severely understudied. We analyzed whole-genome information for a total of 42 representatives of bee-infecting trypanosomatids collected in Central Europe and Alaska from a population genetics point of view. Our data shed light on the evolution, selection, and diversification in this important group of trypanosomatid parasites.

The genomes of Crithidia bombi and C. expoeki, common parasites of bumblebees.

Trypanosomatids (Trypanosomatidae, Kinetoplastida) are flagellated protozoa containing many parasites of medical or agricultural importance. Among those, Crithidia bombi and C. expoeki, are common parasites in bumble bees around the world, and phylogenetically close to Leishmania and Leptomonas. They have a simple and direct life cycle with one host, and partially castrate the founding queens greatly reducing their fitness. Here, we report the nuclear genome sequences of one clone of each species, extracted from a field-collected infection. Using a combination of Roche 454 FLX Titanium, Pacific Biosciences PacBio RS, and Illumina GA2 instruments for C. bombi, and PacBio for C. expoeki, we could produce high-quality and well resolved sequences. We find that these genomes are around 32 and 34 MB, with 7,808 and 7,851 annotated genes for C. bombi and C. expoeki, respectively-which is somewhat less than reported from other trypanosomatids, with few introns, and organized in polycistronic units. A large fraction of genes received plausible functional support in comparison primarily with Leishmania and Trypanosoma. Comparing the annotated genes of the two species with those of six other trypanosomatids (C. fasciculata, L. pyrrhocoris, L. seymouri, B. ayalai, L. major, and T. brucei) shows similar gene repertoires and many orthologs. Similar to other trypanosomatids, we also find signs of concerted evolution in genes putatively involved in the interaction with the host, a high degree of synteny between C. bombi and C. expoeki, and considerable overlap with several other species in the set. A total of 86 orthologous gene groups show signatures of positive selection in the branch leading to the two Crithidia under study, mostly of unknown function. As an example, we examined the initiating glycosylation pathway of surface components in C. bombi, finding it deviates from most other eukaryotes and also from other kinetoplastids, which may indicate rapid evolution in the extracellular matrix that is involved in interactions with the host. Bumble bees are important pollinators and Crithidia-infections are suspected to cause substantial selection pressure on their host populations. These newly sequenced genomes provide tools that should help better understand host-parasite interactions in these pollinator pathogens.

Large-scale cultivation of the bumblebee gut microbiota reveals an underestimated bacterial species diversity capable of pathogen inhibition.

A total of 1940 isolates from gut samples of 60 bumblebees representing Bombus pascuorum, Bombus terrestris, Bombus lucorum and Bombus lapidarius was collected and identified through state-of the-art taxonomic methods. The bacterial species diversity in these Bombus species exceeded that suggested by phylotype analysis through 16S rRNA amplicon sequencing, and revealed that B. pascuorum and B. terrestris had a unique microbiota composition, each. Representatives of most phylotypes reported earlier and detected in the present study were effectively isolated, and included several novel bacterial taxa and species reported for the first time in the bumblebee gut. Isolates were screened in pectin degradation assays and growth inhibition assays against the honeybee pathogens Paenibacillus larvae, Melissococcus plutonius and Ascosphaera apis and the bumblebee parasite Crithidia bombi. While inhibitory activity against each of these pathogens was observed, only one single culture was able to degrade pectin and polygalacturonic acid in vitro. The availability of accurately identified microbial isolates will facilitate future evaluation of the functional potential of the bumblebee gut microbiota.

Parasite infection of specific host genotypes relates to changes in prevalence in two natural populations of bumblebees.

The antagonistic relationship between parasites and their hosts is strongly influenced by genotype-by-genotype interactions. Defense against parasitism is commonly studied in the context of immune system-based mechanisms and, thus, the focus in the search for candidate genes in host-parasite interactions is often on immune genes. In this study, we investigated the association between prevalence of parasite infection and host mitochondrial DNA (mtDNA) haplotypes in two natural populations of bumblebees (Bombus terrestris). The two most common haplotypes of the host populations, termed A and B, differ by a single nonsynonymous nucleotide substitution within the coding region of cytochrome oxidase I, an important player in metabolic pathways. We screened infection by Nosema bombi, a common endoparasite of bumblebees, and the corresponding host mtDNA-haplotype frequencies in over 1400 bumblebees between 2000 and 2010. The island population of Gotland showed lower mtDNA diversity compared to the mainland population in Switzerland. Over time, we observed large fluctuations in infection prevalence, as well as variation in host haplotype frequencies in both populations. Our long-term observation revealed that N. bombi infection of specific host genotypes is transient: We found that with increasing infection prevalence, proportionally more individuals with haplotype B, but fewer individuals with haplotype A were infected. This suggests that the presence of N. bombi in specific host genotypes relates to infection prevalence. This may be a result of parasite competition, or differential resilience of host types to ward off infections. The findings highlight the important role of host mtDNA haplotypes in the interaction with parasites.

The Bee Microbiome: Impact on Bee Health and Model for Evolution and Ecology of Host-Microbe Interactions.

As pollinators, bees are cornerstones for terrestrial ecosystem stability and key components in agricultural productivity. All animals, including bees, are associated with a diverse community of microbes, commonly referred to as the microbiome. The bee microbiome is likely to be a crucial factor affecting host health. However, with the exception of a few pathogens, the impacts of most members of the bee microbiome on host health are poorly understood. Further, the evolutionary and ecological forces that shape and change the microbiome are unclear. Here, we discuss recent progress in our understanding of the bee microbiome, and we present challenges associated with its investigation. We conclude that global coordination of research efforts is needed to fully understand the complex and highly dynamic nature of the interplay between the bee microbiome, its host, and the environment. High-throughput sequencing technologies are ideal for exploring complex biological systems, including host-microbe interactions. To maximize their value and to improve assessment of the factors affecting bee health, sequence data should be archived, curated, and analyzed in ways that promote the synthesis of different studies. To this end, the BeeBiome consortium aims to develop an online database which would provide reference sequences, archive metadata, and host analytical resources. The goal would be to support applied and fundamental research on bees and their associated microbes and to provide a collaborative framework for sharing primary data from different research programs, thus furthering our understanding of the bee microbiome and its impact on pollinator health.

Large scale patterns of abundance and distribution of parasites in Mexican bumblebees.

Bumblebees are highly valued for their pollination services in natural ecosystems as well as for agricultural crops. These precious pollinators are known to be declining worldwide, and one major factor contributing to this decline are infections by parasites. Knowledge about parasites in wild bumblebee populations is thus of paramount importance for conservation purposes. We here report the geographical distribution of Crithidia and Nosema, two common parasites of bumblebees, in a yet poorly investigated country: Mexico. Based on sequence divergence of the Cytochrome b and Glycosomal glyceraldehyde phosphate deshydrogenase (gGPDAH) genes, we discovered the presence of a new Crithidia species, which is mainly distributed in the southern half of the country. It is placed by Bayesian inference as a sister species to C. bombi. We suggest the name Crithidia mexicana for this newly discovered organism. A population of C. expoeki was encountered concentrated on the flanks of the dormant volcanic mountain, Iztaccihuatl, and microsatellite data showed evidence of a bottleneck in this population. This study is the first to provide a large-scale insight into the health status of endemic bumblebees in Mexico, based on a large sample size (n=3,285 bees examined) over a variety of host species and habitats.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim.

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.

A depauperate immune repertoire precedes evolution of sociality in bees.

BACKGROUND: Sociality has many rewards, but can also be dangerous, as high population density and low genetic diversity, common in social insects, is ideal for parasite transmission. Despite this risk, honeybees and other sequenced social insects have far fewer canonical immune genes relative to solitary insects. Social protection from infection, including behavioral responses, may explain this depauperate immune repertoire. Here, based on full genome sequences, we describe the immune repertoire of two ecologically and commercially important bumblebee species that diverged approximately 18 million years ago, the North American Bombus impatiens and European Bombus terrestris. RESULTS: We find that the immune systems of these bumblebees, two species of honeybee, and a solitary leafcutting bee, are strikingly similar. Transcriptional assays confirm the expression of many of these genes in an immunological context and more strongly in young queens than males, affirming Bateman's principle of greater investment in female immunity. We find evidence of positive selection in genes encoding antiviral responses, components of the Toll and JAK/STAT pathways, and serine protease inhibitors in both social and solitary bees. Finally, we detect many genes across pathways that differ in selection between bumblebees and honeybees, or between the social and solitary clades. CONCLUSIONS: The similarity in immune complement across a gradient of sociality suggests that a reduced immune repertoire predates the evolution of sociality in bees. The differences in selection on immune genes likely reflect divergent pressures exerted by parasites across social contexts.

The invasion of southern South America by imported bumblebees and associated parasites.

The Palaearctic Bombus ruderatus (in 1982/1983) and Bombus terrestris (1998) have both been introduced into South America (Chile) for pollination purposes. We here report on the results of sampling campaigns in 2004, and 2010-2012 showing that both species have established and massively expanded their range. Bombus terrestris, in particular, has spread by some 200 km year(-1) and had reached the Atlantic coast in Argentina by the end of 2011. Both species, and especially B. terrestris, are infected by protozoan parasites that seem to spread along with the imported hosts and spillover to native species. Genetic analyses by polymorphic microsatellite loci suggest that the host population of B. terrestris is genetically diverse, as expected from a large invading founder population, and structured through isolation by distance. Genetically, the populations of the trypanosomatid parasite, Crithidia bombi, sampled in 2004 are less diverse, and distinct from the ones sampled later. Current C. bombi populations are highly heterozygous and also structured through isolation by distance correlating with the genetic distances of B. terrestris, suggesting the latter's expansion to be a main structuring factor for the parasite. Remarkably, wherever B. terrestris spreads, the native Bombus dahlbomii disappears although the reasons remain unclear. Our ecological and genetic data suggest a major invasion event that is currently unfolding in southern South America with disastrous consequences for the native bumblebee species.

Probing mixed-genotype infections I: extraction and cloning of infections from hosts of the trypanosomatid Crithidia bombi.

We here present an efficient, precise and reliable method to isolate and cultivate healthy and viable single Crithidia bombi cells from bumblebee faeces using flow cytometry. We report a precision of >99% in obtaining single trypanosomatid cells for further culture and analysis ("cloning"). In the study, we have investigated the use of different liquid media to cultivate C. bombi and present an optimal medium for obtaining viable clones from all tested, infected host donors. We show that this method can be applied to genotype a collection of clones from natural infections. Furthermore, we show how to cryo-preserve C. bombi cells to be revived to become infective clones after at least 4 years of storage. Considering the high prevalence of infections in natural populations, our method provides a powerful tool in studying the level and diversity of these infections, and thus enriches the current methodology for the studies of complex host-parasite interactions.

Probing mixed-genotype infections II: high multiplicity in natural infections of the trypanosomatid, Crithidia bombi, in its host, Bombus spp.

Mixed-genotype infections have major consequences for many essential elements of host-parasite interactions. With genetic exchange between co-infecting parasite genotypes increased diversity among parasite offspring and the emergence of novel genotypes from infected hosts is possible. We here investigated mixed- genotype infections using the host, Bombus spp. and its trypanosome parasite Crithidia bombi as our study case. The natural infections of C. bombi were genotyped with a novel method for a representative sample of workers and spring queens in Switzerland. We found that around 60% of all infected hosts showed mixed-genotype infections with an average of 2.47±0.22 (S.E.) and 3.65±1.02 genotypes per worker or queen, respectively. Queens, however, harboured up to 29 different genotypes. Based on the genotypes of co-infecting strains, these could be putatively assigned to either 'primary' and 'derived' genotypes - the latter resulting from genetic exchange among the primary genotypes. High genetic relatedness among co-infecting derived but not primary genotypes supported this scenario. Co-infection in queens seems to be a major driver for the diversity of genotypes circulating in host populations.

Diversity of Nosema associated with bumblebees (Bombus spp.) from China.

Bumblebees (Bombus spp.) are important pollinators of many economically important crops and microsporidia are among the most important infections of these hosts. Using molecular markers, we screened a large sample (n=1,009 bees) of workers of 27 different Bombus spp. from China (Sichuan, Qinghai, Inner Mongolia, and Gansu provinces). The results showed that 62 individuals representing 12 Bombus spp. were infected by microsporidia with an overall prevalence of 6.1%. Based on the haplotypes (ssrRNA sequences), we confirmed the presence of Nosema bombi, Nosema ceranae and (likely) Nosema thomsoni. In addition, four new putatively novel taxa were identified by phylogenetic reconstruction: Nosema A, Nosema B-complex, Nosema C-complex and Nosema D-complex. In many cases, hosts were infected by more than one Nosema taxon. Possible caveats of sequence analyses are discussed.

Genetic exchange and emergence of novel strains in directly transmitted trypanosomatids.

The breeding structure of protozoan infections, i.e. whether and how frequently parasites exchange genes ("sexual reproduction"), is a crucially important parameter for many important questions; it also matters for how new virulent strains might emerge. Whether protozoan parasites are clonal or sexual is therefore a hotly debated issue. For trypanosomatids, few experimental tests of breeding structure exist to date and are limited to the vector-borne human diseases Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major. We infected the natural host (Bombus terrestris) of the monoxenous parasite Crithidia bombi (Trypanosomatida) either with a single strain of the parasite or in mixed infections and tested for genetic exchange among co-infecting strains using microsatellite markers. We show that strains regularly exchange genetic material, with occasional self-crossing during mixed infections. Most offspring clones fit the expected allelic pattern from a standard Mendelian segregation. In some cases, alleles are lost or gained, leading to an entirely new genotype different from either parent. Genetic exchange in C. bombi therefore does occur and the process also leads to allelic loss or gain that could result from slippage during recombination. The majority of novel offspring types correspond to a recombination of parental alleles. The case of C. bombi demonstrates that directly transmitted, monoxenic trypanosomatids can also exchange genes. Sex therefore seems to be found in very different lineages of the trypanosomatids. Furthermore, the data allowed estimating a frequency at which C. bombi shows genetic exchange in populations.

A second generation genetic map of the bumblebee Bombus terrestris (Linnaeus, 1758) reveals slow genome and chromosome evolution in the Apidae.

BACKGROUND: The bumblebee Bombus terrestris is an ecologically and economically important pollinator and has become an important biological model system. To study fundamental evolutionary questions at the genomic level, a high resolution genetic linkage map is an essential tool for analyses ranging from quantitative trait loci (QTL) mapping to genome assembly and comparative genomics. We here present a saturated linkage map and match it with the Apis mellifera genome using homologous markers. This genome-wide comparison allows insights into structural conservations and rearrangements and thus the evolution on a chromosomal level. RESULTS: The high density linkage map covers ~ 93% of the B. terrestris genome on 18 linkage groups (LGs) and has a length of 2'047 cM with an average marker distance of 4.02 cM. Based on a genome size of ~ 430 Mb, the recombination rate estimate is 4.76 cM/Mb. Sequence homologies of 242 homologous markers allowed to match 15 B. terrestris with A. mellifera LGs, five of them as composites. Comparing marker orders between both genomes we detect over 14% of the genome to be organized in synteny and 21% in rearranged blocks on the same homologous LG. CONCLUSIONS: This study demonstrates that, despite the very high recombination rates of both A. mellifera and B. terrestris and a long divergence time of about 100 million years, the genomes' genetic architecture is highly conserved. This reflects a slow genome evolution in these bees. We show that data on genome organization and conserved molecular markers can be used as a powerful tool for comparative genomics and evolutionary studies, opening up new avenues of research in the Apidae.

Molecular divergence defines two distinct lineages of Crithidia bombi (Trypanosomatidae), parasites of bumblebees.

This study provides, for the first time, sequence data for the protozoan flagellates Crithidia bombi and Crithidia mellificae (Kinetoplastea: Trypanosomatidae). We amplified the partial sequences of the small subunit ribosomal RNA (SSU rRNA), glycosomal glyceraldehyde phosphate dehydrogenase (gGAPDH), cytochrome b (Cyt b), and the complete internal transcribed spacer region 1 (ITS1) of the ribosomal RNA gene region for 66 clones of C. bombi from Switzerland and Alaska. Furthermore, we sequenced the same stretch of SSU rRNA and gGAPDH for one isolate of C. mellificae from Switzerland. All four molecular markers classified the C. bombi samples into two distinct lineages A and B. Both lineages were found in the two sampling locations. Variation within lineages was small or non-existing. Sequence differences between lineages were 1.64% for SSU rRNA, 4.36% for gGAPDH, and 12.02% for Cyt b. The ITS1-sequences of lineages A and B have diverged so much that no alignment was possible. With regard to ITS1, we additionally found fragment length polymorphism (variation in microsatellite repeat numbers) as well as nucleotide diversity within each lineage. Furthermore, the sequences of SSU rRNA and gGAPDH of C. mellificae were different from both lineages of C. bombi. The separation of lineages A and B, based on sequence differences and phylogenetic reconstruction, is so pronounced as to characterize two species of "C. bombi." We propose to retain C. bombi for the more common lineage A and suggest the name Crithidia expoeki n. sp. for lineage B.

Novel microsatellite DNA loci for Bombus terrestris (Linnaeus, 1758).

We present details and characteristics of 123 novel polymorphic microsatellite DNA loci for Bombus terrestris. Thirty-four of these loci have been tested in nine other Bombus species and 25 of them showed polymorphisms in at least one species. These microsatellite DNA loci together with the already established 60 loci will be useful for characterizing wild and managed populations of B. terrestris and other Bombus species as well as for detailed genetic studies in including mapping studies and genome annotations.

Heritability of sperm length in the bumblebee Bombus terrestris.

Sperm length is highly variable, both between and within species, but the evolutionary significance of this variation is poorly understood. Sexual selection on sperm length requires a significant additive genetic variance, but few studies have actually measured this. Here we present the first estimates of narrow sense heritability of sperm length in a social insect, the bumblebee Bombus terrestris. In spite of a balanced and straightforward rearing design of colonies, and the possibility to replicate measurements of sperm within single males nested within colonies, the analysis proved to be complex. Several appropriate statistical models were derived, each depending on different assumptions. The heritability estimates obtained ranged from h (2) = 0.197 +/- 0.091 to h (2) = 0.429 +/- 0.154. All our estimates were substantially lower than previous estimates of sperm length heritability in non-social insects and vertebrates.