Developmental Transcriptomics Reveals a Gene Network Driving Mimetic Color Variation in a Bumble Bee

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

The diversity, evolution, and development of setal morphologies in bumble bees (Hymenoptera: Apidae: Bombus spp.)

Bumble bees are characterized by their thick setal pile that imparts aposematic color patterns often used for species-level identification. Like all bees, the single-celled setae of bumble bees are branched, an innovation thought important for pollen collection. To date no studies have quantified the types of setal morphologies and their distribution on these bees, information that can facilitate understanding of their adaptive ecological function. This study defines several major setal morphotypes in the common eastern bumble bee Bombus impatiens Cresson, revealing these setal types differ by location across the body. The positions of these types of setae are similar across individuals, castes, and sexes within species. We analyzed the distribution of the two most common setal types (plumose and spinulate) across the body dorsum of half of the described bumble bee species. This revealed consistently high density of plumose (long-branched) setae across bumble bees on the head and mesosoma, but considerable variation in the amount of metasomal plumosity. Variation on the metasoma shows strong phylogenetic signal at subgeneric and smaller group levels, making it a useful trait for species delimitation research, and plumosity has increased from early Bombus ancestors. The distribution of these setal types suggests these setae may serve several functions, including pollen-collecting and thermoregulatory roles, and probable mechanosensory functions. This study further examines how and when setae of the pile develop, evidence for mechanosensory function, and the timing of pigmentation as a foundation for future genetic and developmental research in these bees.

Mito-nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern "Rocky Mountain'' color form with ferruginous mid-abdominal segments (B. m. melanopygus) and a southern "Pacific'' form with black mid-abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color-controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito-nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.