On the distribution and diversity of tissue-specific somatic mutations in honey bee (Apis mellifera) drones

Somatic mutations originate from both exogenous (e.g. UV radiation, chemical agents) and endogenous (e.g., DNA replication, defective DNA repair) sources and can have significant impacts on an animal’s reproductive success. This may be especially true for haploid organisms that are susceptible to any deleterious alleles inherited from their parent and any that arise over their lifetime. Unfortunately, little is known about the rate of somatic mutation accumulation across individuals and tissues of haplodiploid animal populations, the functional processes through which they arise, and their distribution across tissues and the genome. Here, we generated short-read whole-genome sequencing data for four tissues of haploid honey bee males. We paired this with estimates of telomere length and tissue-specific DNA content to address three major questions: is there variance in somatic mutational load across haploid individuals and specific tissues therein, does increased DNA content in a tissue contribute to somatic mutational load, and does telomere length correlate with mutational load? Our results suggest that variance in somatic mutational load is better captured across individuals than across tissues, that tissue-specific DNA content is not associated with somatic mutation load, and that variance in telomere length does not correlate with somatic mutation loads across tissues. To our knowledge, this is the first observational study on somatic mutational load in Apoidea and likely Hymenoptera. It serves as a useful advent for additional studies understanding the processes through which haploids tolerate or repair somatic mutations.