What one genus of showy moths can say about migration, adaptation, and wing pattern.

The Ornate Moth, Utetheisa ornatrix, has served as a model species in chemical ecology studies for decades. Like in the widely publicized stories of the Monarch and other milkweed butterflies, the Ornate Moth and its relatives are tropical insects colonizing whole continents assisted by their chemical defenses. With the recent advances in genomic techniques and evo-devo research, it is becoming a model for studies in other areas, from wing pattern development to phylogeography, from toxicology to epigenetics. We used a genomic approach to learn about Utetheisa's evolution, detoxification, dispersal abilities, and wing pattern diversity. We present an evolutionary genomic analysis of the worldwide genus Utetheisa, then focusing on U. ornatrix. Our reference genome of U. ornatrix reveals gene duplications in the regions possibly associated with detoxification abilities, which allows them to feed on toxic food plants. Finally, comparative genomic analysis of over 100 U. ornatrix specimens from the museum with apparent differences in wing patterns suggest the potential roles of cortex and lim3 genes in wing pattern formation of Lepidoptera and the utility of museum-preserved collection specimens for wing pattern research.

Massive seasonal high-altitude migrations of nocturnal insects above the agricultural plains of East China.

Long-distance migrations of insects contribute to ecosystem functioning but also have important economic impacts when the migrants are pests or provide ecosystem services. We combined radar monitoring, aerial sampling, and searchlight trapping, to quantify the annual pattern of nocturnal insect migration above the densely populated agricultural lands of East China. A total of ~9.3 trillion nocturnal insect migrants (15,000 t of biomass), predominantly Lepidoptera, Hemiptera, and Diptera, including many crop pests and disease vectors, fly at heights up to 1 km above this 600 km-wide region every year. Larger migrants (>10 mg) exhibited seasonal reversal of movement directions, comprising northward expansion during spring and summer, followed by southward movements during fall. This north-south transfer was not balanced, however, with southward movement in fall 0.66× that of northward movement in spring and summer. Spring and summer migrations were strongest when the wind had a northward component, while in fall, stronger movements occurred on winds that allowed movement with a southward component; heading directions of larger insects were generally close to the track direction. These findings indicate adaptations leading to movement in seasonally favorable directions. We compare our results from China with similar studies in Europe and North America and conclude that ecological patterns and behavioral adaptations are similar across the Northern Hemisphere. The predominance of pests among these nocturnal migrants has severe implications for food security and grower prosperity throughout this heavily populated region, and knowledge of their migrations is potentially valuable for forecasting pest impacts and planning timely management actions.

A polygenic explanation for Haldane's rule in butterflies.

Two robust rules have been discovered about animal hybrids: Heterogametic hybrids are more unfit (Haldane's rule), and sex chromosomes are disproportionately involved in hybrid incompatibility (the large-X/Z effect). The exact mechanisms causing these rules in female heterogametic taxa such as butterflies are unknown but are suggested by theory to involve dominance on the sex chromosome. We investigate hybrid incompatibilities adhering to both rules in Papilio and Heliconius butterflies and show that dominance theory cannot explain our data. Instead, many defects coincide with unbalanced multilocus introgression between the Z chromosome and all autosomes. Our polygenic explanation predicts both rules because the imbalance is likely greater in heterogametic females, and the proportion of introgressed ancestry is more variable on the Z chromosome. We also show that mapping traits polygenic on a single chromosome in backcrosses can generate spurious large-effect QTLs. This mirage is caused by statistical linkage among polygenes that inflates estimated effect sizes. By controlling for statistical linkage, most incompatibility QTLs in our hybrid crosses are consistent with a polygenic basis. Since the two genera are very distantly related, polygenic hybrid incompatibilities are likely common in butterflies.

Meiotic drive against chromosome fusions in butterfly hybrids.

Species frequently differ in the number and structure of chromosomes they harbor, but individuals that are heterozygous for chromosomal rearrangements may suffer from reduced fitness. Chromosomal rearrangements like fissions and fusions can hence serve as a mechanism for speciation between incipient lineages, but their evolution poses a paradox. How can rearrangements get fixed between populations if heterozygotes have reduced fitness? One solution is that this process predominantly occurs in small and isolated populations, where genetic drift can override natural selection. However, fixation is also more likely if a novel rearrangement is favored by a transmission bias, such as meiotic drive. Here, we investigate chromosomal transmission distortion in hybrids between two wood white (Leptidea sinapis) butterfly populations with extensive karyotype differences. Using data from two different crossing experiments, we uncover that there is a transmission bias favoring the ancestral chromosomal state for derived fusions, a result that shows that chromosome fusions actually can fix in populations despite being counteracted by meiotic drive. This means that meiotic drive not only can promote runaway chromosome number evolution and speciation, but also that it can be a conservative force acting against karyotypic change and the evolution of reproductive isolation. Based on our results, we suggest a mechanistic model for why chromosome fusion mutations may be opposed by meiotic drive and discuss factors contributing to karyotype evolution in Lepidoptera.

Dosage compensation in Bombyx mori is achieved by partial repression of both Z chromosomes in males.

SignificanceGenes on sex chromosomes (i.e. human chX) are regulated differently in males and females to balance gene expression levels between sexes (XY vs. XX). This sex-specific regulation is called dosage compensation (DC). DC is achieved by altering the shape and compaction of sex chromosomes specifically in one sex. In this study, we use Oligopaints to examine DC in silkworms. This study visualizes this phenomenon in a species with ZW sex chromosomes, which evolved independently of XY. Our data support a long-standing model for how DC mechanisms evolved across species, and we show potential similarity between DC in silkworms and nematodes, suggesting that this type of DC may have emerged multiple independent times throughout evolution.

Genome sequencing and assembly of Indian golden silkmoth, Antheraea assamensis Helfer (Saturniidae, Lepidoptera).

Muga silkworm (Antheraea assamensis), one of the economically important wild silkmoths, is unique among saturniid silkmoths. It is confined to the North-eastern part of India. Muga silk has the highest value among the other silks. Unlike other silkmoths, A. assamensis has a low chromosome number (n = 15), and ZZ/ZO sex chromosome system. Here, we report the first high-quality draft genome of A. assamensis, assembled by employing the Illumina and PacBio sequencing platforms. The assembled genome of A. assamensis is 501.18 Mb long, with 2697 scaffolds and an N50 of 683.23 Kb. The genome encompasses 18,385 protein-coding genes, 86.29% of which were functionally annotated. Phylogenetic analysis of A. assamensis revealed its divergence from other Antheraea species approximately 28.7 million years ago. Moreover, an investigation into detoxification-related gene families, CYP450, GST, and ABC-transporter, revealed a significant expansion in A. assamensis as compared to the Bombyx mori. This expansion is comparable to Spodoptera litura, suggesting adaptive responses linked to the polyphagous behavior observed in these insects. This study provides valuable insights into the molecular basis of evolutionary divergence and adaptations in muga silkmoth. The genome assembly reported in this study will significantly help in the functional genomics studies on A. assamensis and other Antheraea species along with comparative genomics analyses of Bombycoidea insects.

Toward an Integrated Understanding of the Lepidoptera Microbiome.

Research over the past 30 years has led to a widespread acceptance that insects establish widespread and diverse associations with microorganisms. More recently, microbiome research has been accelerating in lepidopteran systems, leading to a greater understanding of both endosymbiont and gut microorganisms and how they contribute to integral aspects of the host. Lepidoptera are associated with a robust assemblage of microorganisms, some of which may be stable and routinely detected in larval and adult hosts, while others are ephemeral and transient. Certain microorganisms that populate Lepidoptera can contribute significantly to the hosts' performance and fitness, while others are inconsequential. We emphasize the context-dependent nature of the interactions between players. While our review discusses the contemporary literature, there are major avenues yet to be explored to determine both the fundamental aspects of host-microbe interactions and potential applications for the lepidopteran microbiome; we describe these avenues after our synthesis.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae.

BACKGROUND: The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops. RESULTS: The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae. CONCLUSIONS: Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

Opsin Gene Duplication in Lepidoptera: Retrotransposition, Sex Linkage, and Gene Expression.

Color vision in insects is determined by signaling cascades, central to which are opsin proteins, resulting in sensitivity to light at different wavelengths. In certain insect groups, lineage-specific evolution of opsin genes, in terms of copy number, shifts in expression patterns, and functional amino acid substitutions, has resulted in changes in color vision with subsequent behavioral and niche adaptations. Lepidoptera are a fascinating model to address whether evolutionary change in opsin content and sequence evolution are associated with changes in vision phenotype. Until recently, the lack of high-quality genome data representing broad sampling across the lepidopteran phylogeny has greatly limited our ability to accurately address this question. Here, we annotate opsin genes in 219 lepidopteran genomes representing 33 families, reconstruct their evolutionary history, and analyze shifts in selective pressures and expression between genes and species. We discover 44 duplication events in opsin genes across ∼300 million years of lepidopteran evolution. While many duplication events are species or family specific, we find retention of an ancient long-wavelength-sensitive (LW) opsin duplication derived by retrotransposition within the speciose superfamily Noctuoidea (in the families Nolidae, Erebidae, and Noctuidae). This conserved LW retrogene shows life stage-specific expression suggesting visual sensitivities or other sensory functions specific to the early larval stage. This study provides a comprehensive order-wide view of opsin evolution across Lepidoptera, showcasing high rates of opsin duplications and changes in expression patterns.

A time course analysis through diapause reveals dynamic temporal patterns of microRNAs associated with endocrine regulation in the butterfly Pieris napi.

Organisms inhabiting highly seasonal environments must cope with a wide range of environmentally induced challenges. Many seasonal challenges require extensive physiological modification to survive. In winter, to survive extreme cold and limited resources, insects commonly enter diapause, which is an endogenously derived dormant state associated with minimized cellular processes and low energetic expenditure. Due to the high degree of complexity involved in diapause, substantial cellular regulation is required, of which our understanding primarily derives from the transcriptome via messenger RNA expression dynamics. Here we aim to advance our understanding of diapause by investigating microRNA (miRNA) expression in diapausing and direct developing pupae of the butterfly Pieris napi. We identified coordinated patterns of miRNA expression throughout diapause in both head and abdomen tissues of pupae, and via miRNA target identification, found several expression patterns to be enriched for relevant diapause-related physiological processes. We also identified two candidate miRNAs, miR-14-5p and miR-2a-3p, that are likely involved in diapause progression through their activity in the ecdysone pathway, a critical regulator of diapause termination. miR-14-5p targets phantom, a gene in the ecdysone synthesis pathway, and is upregulated early in diapause. miR-2a-3p has been found to be expressed in response to ecdysone, and is upregulated during diapause termination. Together, the expression patterns of these two miRNAs match our current understanding of the timing of hormonal regulation of diapause in P. napi and provide interesting candidates to further explore the mechanistic role of microRNAs in diapause regulation.

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem.

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera).

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus.

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Transcriptomic and proteomic elucidation of Z chromosome dosage compensation in Helicoverpa armigera.

Transcriptomic data have been used to study sex chromosome dosage compensation (SCDC) in approximately 10 Lepidoptera ZW species, yielding a consensus compensation pattern of Z ≈ ZZ < AA . $$ \approx \mathrm{ZZ}<\mathrm{AA}. $$ It remains unclear whether this compensation pattern holds when examining more Lepidoptera ZW species and/or using proteomic data to analyse SCDC. Here we combined transcriptomic and proteomic data as well as transcriptional level of six individual Z genes to reveal the SCDC pattern in Helicoverpa armigera, a polyphagous lepidopteran pest of economic importance. Transcriptomic analysis showed that the Z chromosome expression of H. armigera was balanced between male and female but substantially reduced relative to autosome expression, exhibiting an SCDC pattern of Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ . When using H. amigera midgut proteomic data, the SCDC pattern of this species changed from Z ≈ ZZ < AA $$ \approx \mathrm{ZZ}<\mathrm{AA} $$ at transcriptomic level to Z = ZZ = AA at the proteomic level. RT-qPCR analysis of transcript abundance of six Z genes found that compensation for each Z gene could vary from no compensation to overcompensation, depending on the individual genes and tissues tested. These results demonstrate for the first time the existence of a translational compensation mechanism, which is operating in addition to a translational mechanism, such as has been reported in other lepidopteran species. And the transcriptional compensation mechanism functions to accomplish Z chromosome dosage balance between the sexes (M = F on the Z chromosome), whereas the translation compensation mechanism operates to achieve dosage compensation between Z chromosome and autosome (Z = AA).

Rising minimum temperatures contribute to 50 years of occupancy decline among cold-adapted Arctic and boreal butterflies in North America.

Global climate change has been identified as a potential driver of observed insect declines, yet in many regions, there are critical data gaps that make it difficult to assess how communities are responding to climate change. Poleward regions are of particular interest because warming is most rapid while biodiversity data are most sparse. Building on recent advances in occupancy modeling of presence-only data, we reconstructed 50 years (1970-2019) of butterfly occupancy trends in response to rising minimum temperatures in one of the most under-sampled regions of North America. Among 90 modeled species, we found that cold-adapted species are far more often in decline compared with their warm-adapted, more southernly distributed counterparts. Furthermore, in a post hoc analysis using species' traits, we find that species' range-wide average annual temperature is the only consistent predictor of occupancy changes. Species with warmer ranges were most likely to be increasing in occupancy. This trend results in the majority of butterflies increasing in occupancy probability over the last 50 years. Our results provide the first look at macroscale butterfly biodiversity shifts in high-latitude North America. These results highlight the potential of leveraging the wealth of presence-only data, the most abundant source of biodiversity data, for inferring changes in species distributions.

Substantial urbanization-driven declines of larval and adult moths in a subtropical environment.

Recent work has shown the decline of insect abundance, diversity and biomass, with potential implications for ecosystem services. These declines are especially pronounced in regions with high human activity, and urbanization is emerging as a significant contributing factor. However, the scale of these declines and the traits that determine variation in species-specific responses remain less well understood, especially in subtropical and tropical regions, where insect diversity is high and urban footprints are rapidly expanding. Here, we surveyed moths across an entire year in protected forested sites across an urbanization gradient to test how caterpillar and adult life stages of subtropical moths (Lepidoptera) are impacted by urbanization. Specifically, we assess how urban development affects the total biomass of caterpillars, abundance of adult moths and quantify how richness and phylogenetic diversity of macro-moths are impacted by urban development. Additionally, we explore how life-history traits condition species' responses to urban development. At the community level, we find that urban development decreases caterpillar biomass and adult moth abundance. We also find sharp declines of adult macro-moths in response to urban development across the phylogeny, leading to a decrease in species richness and phylogenetic diversity in more urban sites. Finally, our study found that smaller macro-moths are less impacted by urban development than larger macro-moths in subtropical environments, perhaps highlighting the tradeoffs of metabolic costs of urban heat favoring smaller moths over the relative benefits of dispersal for larger moths. In summary, our research underscores the far-reaching consequences of urbanization on moths and provides compelling evidence that urban forests alone may not be sufficient to safeguard biodiversity in cities.

Diet influences resource allocation in chemical defence but not melanin synthesis in an aposematic moth.

For animals that synthesise their chemical compounds de novo, resources, particularly proteins, can influence investment in chemical defences and nitrogen-based wing colouration such as melanin. Competing for the same resources often leads to trade-offs in resource allocation. We manipulated protein availability in the larval diet of the wood tiger moth, Arctia plantaginis, to test how early life resource availability influences relevant life history traits, melanin production and chemical defences. We expected higher dietary protein to result in more effective chemical defences in adult moths and a higher amount of melanin in the wings. According to the resource allocation hypothesis, we also expected individuals with less melanin to have more resources to allocate to chemical defences. We found that protein-deprived moths had a slower larval development, and their chemical defences were less unpalatable for bird predators, but the expression of melanin in their wings did not differ from that of moths raised on a high-protein diet. The amount of melanin in the wings, however, unexpectedly correlated positively with chemical defences. Our findings demonstrate that the resources available in early life have an important role in the efficacy of chemical defences, but melanin-based warning colours are less sensitive to resource variability than other fitness-related traits.

Local climate change velocities and evolutionary history explain multidirectional range shifts in a North American butterfly assemblage.

Species are often expected to shift their distributions either poleward or upslope to evade warming climates and colonise new suitable climatic niches. However, from 18-years of fixed transect monitoring data on 88 species of butterfly in the midwestern United States, we show that butterflies are shifting their centroids in all directions, except towards regions that are warming the fastest (southeast). Butterflies shifted their centroids at a mean rate of 4.87 km year-1. The rate of centroid shift was significantly associated with local climate change velocity (temperature by precipitation interaction), but not with mean climate change velocity throughout the species' ranges. Species tended to shift their centroids at a faster rate towards regions that are warming at slower velocities but increasing in precipitation velocity. Surprisingly, species' thermal niche breadth (range of climates butterflies experience throughout their distribution) and wingspan (often used as metric for dispersal capability) were not correlated with the rate at which species shifted their ranges. We observed high phylogenetic signal in the direction species shifted their centroids. However, we found no phylogenetic signal in the rate species shifted their centroids, suggesting less conserved processes determine the rate of range shift than the direction species shift their ranges. This research shows important signatures of multidirectional range shifts (latitudinal and longitudinal) and uniquely shows that local climate change velocities are more important in driving range shifts than the mean climate change velocity throughout a species' entire range.

Selfish herd effects in aggregated caterpillars and their interaction with warning signals.

Larval Lepidoptera gain survival advantages by aggregating, especially when combined with aposematic warning signals, yet reductions in predation risk may not be experienced equally across all group members. Hamilton's selfish herd theory predicts that larvae that surround themselves with their group mates should be at lower risk of predation, and those on the periphery of aggregations experience the greatest risk, yet this has rarely been tested. Here, we expose aggregations of artificial 'caterpillar' targets to predation from free-flying, wild birds to test for marginal predation when all prey are equally accessible and for an interaction between warning coloration and marginal predation. We find that targets nearer the centre of the aggregation survived better than peripheral targets and nearby targets isolated from the group. However, there was no difference in survival between peripheral and isolated targets. We also find that grouped targets survived better than isolated targets when both are aposematic, but not when they are non-signalling. To our knowledge, our data provide the first evidence to suggest that avian predators preferentially target peripheral larvae from aggregations and that prey warning signals enhance predator avoidance of groups.

Haemolymph transfusions transfer heritable learned novel odour preferences to naive larvae of Bicyclus anynana butterflies.

The mechanisms whereby environmental experiences of parents are transmitted to their offspring to impact their behaviour and fitness are poorly understood. Previously, we showed that naive Bicyclus anynana butterfly larvae, whose parents fed on a normal plant feed but coated with a novel odour, inherited an acquired preference towards that odour, which had initially elicited avoidance in the naive parents. Here, we performed simple haemolymph transfusions from odour-fed and control-fed larvae to naive larval recipients. We found that larvae injected with haemolymph from odour-fed donors stopped avoiding the novel odour, and their naive offspring preferred the odour more, compared to the offspring of larvae injected with control haemolymph. These results indicate that factors in the haemolymph, potentially the odour molecule itself, play an important role in odour learning and preference transmission across generations. Furthermore, this mechanism of odour preference inheritance, mediated by the haemolymph, bypasses the peripheral odour-sensing mechanisms taking place in the antennae, mouthparts or legs, and may mediate food plant switching and diversification in Lepidoptera or more broadly across insects.