19The conserved fatty acyl reductase (FAR) family is involved in biosynthesis of fatty alcohols that 20 serve a range of biological roles. In moths, butterflies (Lepidoptera), and bees (Hymenoptera), 21FARs biosynthesize fatty alcohol pheromones participating in mate-finding strategies. Using a 22 combination of next-generation sequencing, analysis of transposable elements (TE) in the 23 genomic environment of FAR genes, and functional characterization of FARs from Bombus 24 lucorum, B. lapidarius, and B. terrestris, we uncovered a massive expansion of the FAR gene 25 family in Hymenoptera, presumably facilitated by TEs. Expansion occurred in the common 26 ancestor of bumblebees (Bombini) and stingless bees (Meliponini) after their divergence from the 27 honeybee lineage. We found that FARs from the expanded FAR-A orthology group contributed to 28 the species-specific male marking pheromone composition. Our results indicate that TE-mediated 29 expansion and functional diversification of the FAR gene family played a key role in the 30 evolution of pheromone communication in the crown group of Hymenoptera. 31 32 33 34 35 Abbreviations: MMP: male marking pheromone, FA: fatty acid, FAME: fatty acid methyl ester, 36 FAR: fatty acyl reductase, LG: labial gland, FB: fat body, TE: transposable element. 37 38 (Lepidoptera) are the most well-studied model of insect pheromone biosynthesis and have been 64 the subject of substantial research effort related to FARs. Variation in FAR enzymatic specificities 65 is a source of sex pheromone signal diversity among moths in the genus Ostrinia 19 and is also 66 responsible for the distinct pheromone composition in two reproductively isolated races of the 67European corn borer Ostrinia nubilalis 20 . Divergence in pheromone biosynthesis can potentially 68 install or strengthen reproductive barriers, ultimately leading to speciation 21 . However, the 69 biological significance of a large number of insect FAR paralogs remains unclear, as all FARs 70 implicated in moth and butterfly sex pheromone biosynthesis are restricted to a single clade, 71indicating that one FAR group was exclusively recruited for pheromone biosynthesis 20,22-24 . 72 While more than 20 FARs have been experimentally characterized from 23 moth and butterfly 73 (Lepidoptera) species 25 , FARs from other insect orders have received far less attention. Single 74 FAR genes have been isolated and experimentally characterized from Drosophila (Diptera) 14 , the 75 European honeybee (Hymenoptera) 26 and the scale insect Ericeus pela (Hemiptera) 27 . Our 76 limited knowledge about FAR function prevents us from drawing inferences about the biological 77 significance of the FAR gene family expansion in insects. 78 Bumblebees (Hymenoptera: Apidae) are a convenient experimental model to study 79 insect FAR evolution because the majority of bumblebee species produces fatty alcohols as 80 species-specific components of male marking pheromones (MMPs) 28 , which are presumed to be 81 biosynthesized by some of the numerous bumbl...