A long noncoding RNA at the cortex locus controls adaptive coloration in butterflies.

Evolutionary variation in the wing pigmentation of butterflies and moths offers striking examples of adaptation by crypsis and mimicry. The cortex locus has been independently mapped as the locus controlling color polymorphisms in 15 lepidopteran species, suggesting that it acts as a genomic hotspot for the diversification of wing patterns, but functional validation through protein-coding knockouts has proven difficult to obtain. Our study unveils the role of a long noncoding RNA (lncRNA) which we name ivory, transcribed from the cortex locus, in modulating color patterning in butterflies. Strikingly, ivory expression prefigures most melanic patterns during pupal development, suggesting an early developmental role in specifying scale identity. To test this, we generated CRISPR mosaic knock-outs in five nymphalid butterfly species and show that ivory mutagenesis yields transformations of dark pigmented scales into white or light-colored scales. Genotyping of Vanessa cardui germline mutants associates these phenotypes to small on-target deletions at the conserved first exon of ivory. In contrast, cortex germline mutant butterflies with confirmed null alleles lack any wing phenotype and exclude a color patterning role for this adjacent gene. Overall, these results show that a lncRNA gene acts as a master switch of color pattern specification and played key roles in the adaptive diversification of wing patterns in butterflies.

The ivory lncRNA regulates seasonal color patterns in buckeye butterflies.

Long noncoding RNAs (lncRNAs) are transcribed elements increasingly recognized for their roles in regulating gene expression. Thus far, however, we have little understanding of how lncRNAs contribute to evolution and adaptation. Here, we show that a conserved lncRNA, ivory, is an important color patterning gene in the buckeye butterfly Junonia coenia. ivory overlaps with cortex, a locus linked to multiple cases of crypsis and mimicry in Lepidoptera. Along with a companion paper by Livraghi et al., we argue that ivory, not cortex, is the color pattern gene of interest at this locus. In J. coenia, a cluster of cis-regulatory elements (CREs) in the first intron of ivory are genetically associated with natural variation in seasonal color pattern plasticity, and targeted deletions of these CREs phenocopy seasonal phenotypes. Deletions of different ivory CREs produce other distinct phenotypes as well, including loss of melanic eyespot rings, and positive and negative changes in overall wing pigmentation. We show that the color pattern transcription factors Spineless, Bric-a-brac, and Ftz-f1 bind to the ivory promoter during wing pattern development, suggesting that they directly regulate ivory. This case study demonstrates how cis-regulation of a single noncoding RNA can exert diverse and nuanced effects on the evolution and development of color patterns, including modulating seasonally plastic color patterns.

Sex-linked gene traffic underlies the acquisition of sexually dimorphic UV color vision in Heliconius butterflies.

The acquisition of novel sexually dimorphic traits poses an evolutionary puzzle: How do new traits arise and become sex-limited? Recently acquired color vision, sexually dimorphic in animals like primates and butterflies, presents a compelling model for understanding how traits become sex-biased. For example, some Heliconius butterflies uniquely possess UV (ultraviolet) color vision, which correlates with the expression of two differentially tuned UV-sensitive rhodopsins, UVRh1 and UVRh2. To discover how such traits become sexually dimorphic, we studied Heliconius charithonia, which exhibits female-specific UVRh1 expression. We demonstrate that females, but not males, discriminate different UV wavelengths. Through whole-genome shotgun sequencing and assembly of the H. charithonia genome, we discovered that UVRh1 is present on the W chromosome, making it obligately female-specific. By knocking out UVRh1, we show that UVRh1 protein expression is absent in mutant female eye tissue, as in wild-type male eyes. A PCR survey of UVRh1 sex-linkage across the genus shows that species with female-specific UVRh1 expression lack UVRh1 gDNA in males. Thus, acquisition of sex linkage is sufficient to achieve female-specific expression of UVRh1, though this does not preclude other mechanisms, like cis-regulatory evolution from also contributing. Moreover, both this event, and mutations leading to differential UV opsin sensitivity, occurred early in the history of Heliconius. These results suggest a path for acquiring sexual dimorphism distinct from existing mechanistic models. We propose a model where gene traffic to heterosomes (the W or the Y) genetically partitions a trait by sex before a phenotype shifts (spectral tuning of UV sensitivity).

The Afrotropical breeding grounds of the Palearctic-African migratory painted lady butterflies (Vanessa cardui).

Migratory insects are key players in ecosystem functioning and services, but their spatiotemporal distributions are typically poorly known. Ecological niche modeling (ENM) may be used to predict species seasonal distributions, but the resulting hypotheses should eventually be validated by field data. The painted lady butterfly (Vanessa cardui) performs multigenerational migrations between Europe and Africa and has become a model species for insect movement ecology. While the annual migration cycle of this species is well understood for Europe and northernmost Africa, it is still unknown where most individuals spend the winter. Through ENM, we previously predicted suitable breeding grounds in the subhumid regions near the tropics between November and February. In this work, we assess the suitability of these predictions through i) extensive field surveys and ii) two-year monitoring in six countries: a large-scale monitoring scheme to study butterfly migration in Africa. We document new breeding locations, year-round phenological information, and hostplant use. Field observations were nearly always predicted with high probability by the previous ENM, and monitoring demonstrated the influence of the precipitation seasonality regime on migratory phenology. Using the updated dataset, we built a refined ENM for the Palearctic-African range of V. cardui. We confirm the relevance of the Afrotropical region and document the missing natural history pieces of the longest migratory cycle described in butterflies.

Ultrabithorax modifies a regulatory network of genes essential for butterfly eyespot development in a wing sector-specific manner.

Nymphalid butterfly species often have a different number of eyespots in forewings and hindwings, but how the hindwing identity gene Ultrabithorax (Ubx) drives this asymmetry is not fully understood. We examined a three-gene regulatory network for eyespot development in the hindwings of Bicyclus anynana butterflies and compared it with the same network previously described for forewings. We also examined how Ubx interacts with each of these three eyespot-essential genes. We found similar genetic interactions between the three genes in fore- and hindwings, but we discovered three regulatory differences: Antennapedia (Antp) merely enhances spalt (sal) expression in the eyespot foci in hindwings, but is not essential for sal activation, as in forewings; Ubx upregulates Antp in all hindwing eyespot foci but represses Antp outside these wing regions; and Ubx regulates sal in a wing sector-specific manner, i.e. it activates sal expression only in the sectors that have hindwing-specific eyespots. We propose a model for how the regulatory connections between these four genes evolved to produce wing- and sector-specific variation in eyespot number.

Unusual magnetotransport in twisted bilayer graphene.

SignificanceWhen two sheets of graphene are twisted to the magic angle of 1.1∘, the resulting flat moiré bands can host exotic correlated electronic states such as superconductivity and ferromagnetism. Here, we show transport properties of a twisted bilayer graphene device at 1.38∘, far enough above the magic angle that we do not expect exotic correlated states. Instead, we see several unusual behaviors in the device's resistivity upon tuning both charge carrier density and perpendicular magnetic field. We can reproduce these behaviors with a surprisingly simple model based on Hofstadter's butterfly. These results shed light on the underlying properties of twisted bilayer graphene.

Testing the selective sequestration hypothesis: Monarch butterflies preferentially sequester plant defences that are less toxic to themselves while maintaining potency to others.

Herbivores that sequester toxins are thought to have cracked the code of plant defences. Nonetheless, coevolutionary theory predicts that plants should evolve toxic variants that also negatively impact specialists. We propose and test the selective sequestration hypothesis, that specialists preferentially sequester compounds that are less toxic to themselves while maintaining toxicity to enemies. Using chemically distinct plants, we show that monarch butterflies sequester only a subset of cardenolides from milkweed leaves that are less potent against their target enzyme (Na+ /K+ -ATPase) compared to several dominant cardenolides from leaves. However, sequestered compounds remain highly potent against sensitive Na+ /K+ -ATPases found in most predators. We confirmed this differential toxicity with mixtures of purified cardenolides from leaves and butterflies. The genetic basis of monarch adaptation to sequestered cardenolides was also confirmed with transgenic Drosophila that were CRISPR-edited with the monarch's Na+ /K+ -ATPase. Thus, the monarch's selective sequestration appears to reduce self-harm while maintaining protection from enemies.

Got milkweed? Genetic assimilation as potential source for the evolution of nonmigratory monarch butterfly wing shape.

Monarch butterflies (Danaus plexippus) are well studied for their annual long-distance migration from as far north as Canada to their overwintering grounds in Central Mexico. At the end of the cold season, monarchs start to repopulate North America through short-distance migration over the course of multiple generations. Interestingly, some populations in various tropical and subtropical islands do not migrate and exhibit heritable differences in wing shape and size, most likely an adaptation to island life. Less is known about forewing differences between long- and short-distance migrants in relation to island populations. Given their different migratory behaviors, we hypothesized that these differences would be reflected in wing morphology. To test this, we analyzed forewing shape and size of three different groups: nonmigratory, lesser migratory (migrate short-distances), and migratory (migrate long-distances) individuals. Significant differences in shape appear in all groups using geometric morphometrics. As variation found between migratory and lesser migrants has been shown to be caused by phenotypic plasticity, and lesser migrants develop intermediate forewing shapes between migratory and nonmigratory individuals, we suggest that genetic assimilation might be an important mechanism to explain the heritable variation found between migratory and nonmigratory populations. Additionally, our research confirms previous studies which show that forewing size is significantly smaller in nonmigratory populations when compared to both migratory phenotypes. Finally, we found sexual dimorphism in forewing shape in all three groups, but for size in nonmigratory populations only. This might have been caused by reduced constraints on forewing size in nonmigratory populations.

A butterfly egg-killing hypersensitive response in Brassica nigra is controlled by a single locus, PEK, containing a cluster of TIR-NBS-LRR receptor genes.

Knowledge of plant recognition of insects is largely limited to a few resistance (R) genes against sap-sucking insects. Hypersensitive response (HR) characterizes monogenic plant traits relying on R genes in several pathosystems. HR-like cell death can be triggered by eggs of cabbage white butterflies (Pieris spp.), pests of cabbage crops (Brassica spp.), reducing egg survival and representing an effective plant resistance trait before feeding damage occurs. Here, we performed genetic mapping of HR-like cell death induced by Pieris brassicae eggs in the black mustard Brassica nigra (B. nigra). We show that HR-like cell death segregates as a Mendelian trait and identified a single dominant locus on chromosome B3, named PEK (Pieris  egg- killing). Eleven genes are located in an approximately 50 kb region, including a cluster of genes encoding intracellular TIR-NBS-LRR (TNL) receptor proteins. The PEK locus is highly polymorphic between the parental accessions of our mapping populations and among B. nigra reference genomes. Our study is the first one to identify a single locus potentially involved in HR-like cell death induced by insect eggs in B. nigra. Further fine-mapping, comparative genomics and validation of the PEK locus will shed light on the role of these TNL receptors in egg-killing HR.

Heat waves induce milkweed resistance to a specialist herbivore via increased toxicity and reduced nutrient content.

Over the last decade, a large effort has been made to understand how extreme climate events disrupt species interactions. Yet, it is unclear how these events affect plants and herbivores directly, via metabolic changes, and indirectly, via their subsequent altered interaction. We exposed common milkweed (Asclepias syriaca) and monarch caterpillars (Danaus plexippus) to control (26:14°C, day:night) or heat wave (HW) conditions (36:24°C, day:night) for 4 days and then moved each organism to a new control or HW partner to disentangle the direct and indirect effects of heat exposure on each organism. We found that the HW directly benefited plants in terms of growth and defence expression (increased latex exudation and total cardenolides) and insect her'bivores through faster larval development. Conversely, indirect HW effects caused both plant latex and total cardenolides to decrease after subsequent herbivory. Nonetheless, increasing trends of more toxic cardenolides and lower leaf nutritional quality after herbivory by HW caterpillars likely led to reduced plant damage compared to controls. Our findings reveal that indirect impacts of HWs may play a greater role in shaping plant-herbivore interactions via changes in key physiological traits, providing valuable understanding of how ecological interactions may proceed in a changing world.

The Impact of Extensive Surgical Resection of Butterfly Glioblastomas on Outcomes in the Presence of TERT Mutation and EGFR Amplification: A Retrospective Cohort Study.

BACKGROUND AND OBJECTIVES: This study aimed to assess if extensive surgical resection enhances outcomes in wild-type Isocitrate Dehydrogenase (IDH) butterfly glioblastoma (B-GBM) patients, despite the presence of Telomerase Reverse Transcriptase (TERT) mutation and Epidermal Growth Factor Receptor (EGFR) amplification. METHODS: The study, retrospectively conducted from 2014 to 2022, involved 723 GBM patients, 41 of whom met the criteria for IDH wild-type B-GBM. Exclusion criteria comprised prior diagnoses or treatments for low-grade glial tumors. Surgeons, employing two approaches-partial and extensive surgery-categorized patients based on age, sex, tumor location, corpus callosum involvement, and genetic characteristics. The interval between initial surgery and tumor recurrence/tumor-free period (TR/TFP) and overall survival (OS) were recorded and compared between the partial and extensive resection groups, analyzing the impact of resection width on TR/TFP and OS. Preoperative assessments utilized thin-section cranial computed tomography (CT) and contrast-enhanced magnetic resonance imaging (MRI). Intraoperatively, tumor excision was guided by sodium fluorescein, and margins were delineated via neuronavigation. Genetic alterations (TERT mutations and EGFR amplifications) were correlated with surgical type, TR/TFP, and OS. Karnofsky Performance Scale (KPS) evaluations were performed pre- and post-operatively and at key intervals, comparing outcomes between surgical groups. Standard radiotherapy and chemotherapy regimens were administered to all patients. RESULTS: Extensive resection yielded significantly longer TR/TFP compared to partial resection, despite TERT gene mutation and EGFR amplification being linked to shorter TR/TFP and OS. Its impact on OS, however, was not significant. KPS scores indicated a superior quality of life after extensive resection, with sustained improvement upon recurrence. CONCLUSIONS: Extensive resection of B-GBM, even in the presence of adverse genetic alterations, may prolong TR/TFP, offering patients a period of improved comfort with minimal distress.

High toxin concentration in pollen may deter collection by bees in butterfly-pollinated Rhododendron molle.

BACKGROUNDS AND AIMS: The hypothesis that plants evolve features that protect accessible pollen from consumption by flower visitors remains poorly understood. METHODS: To explore potential chemical defence against pollen consumption, we examined the pollinator assemblage, foraging behaviour, visitation frequency and pollen transfer efficiency in Rhododendron molle, a highly toxic shrub containing rhodojaponin III. Nutrient (protein and lipid) and toxic components in pollen and other tissues were measured. KEY RESULTS: Overall in the five populations studied, floral visits by butterflies and bumblebees were relatively more frequent than visits by honeybees. All foraged for nectar but not pollen. Butterflies did not differ from bumblebees in the amount of pollen removed per visit, but deposited more pollen per visit. Pollination experiments indicated that R. molle was self-compatible, but both fruit and seed production were pollen-limited. Our analysis indicated that the pollen was not protein-poor and had a higher concentration of the toxic compound rhodojaponin III than petals and leaves, this compound was undetectable in nectar. CONCLUSION: Pollen toxicity in Rhododendron flowers may discourage pollen robbers (bees) from taking the freely accessible pollen grains, while the toxin-free nectar rewards effective pollinators, promoting pollen transfer. This preliminary study supports the hypothesis that chemical defence in pollen would be likely to evolve in species without physical protection from pollinivores.

Potent Nitrogen-containing Milkweed Toxins are Differentially Regulated by Soil Nitrogen and Herbivore-induced Defense.

Theories have been widely proposed and tested for impacts of soil nitrogen (N) on phytochemical defenses. Among the hundreds of distinct cardenolide toxins produced by milkweeds (Asclepias spp.), few contain N, yet these appear to be the most toxic against specialist herbivores. Because N- and non-N-cardenolides coexist in milkweed leaves and likely have distinct biosynthesis, they present an opportunity to address hypotheses about drivers of toxin expression. We tested effects of soil N and herbivore-damage on cardenolide profiles of two milkweed species differing in life-history strategies (Asclepias syriaca and A. curassavica), and the toxicity of their leaves. In particular leaf extracts were tested against the target enzymes (Na+/K+-ATPase extracted from neural tissue) from both monarch butterflies (Danaus plexippus) as well as less cardenolide-resistant queen butterflies, D. gilippus. Increasing soil N enhanced biomass of Asclepias syriaca but had weak effects on cardenolides, including causing a significant reduction in the N-cardenolide labriformin; feeding by monarch caterpillars strongly induced N-cardenolides (labriformin), its precursors, and total cardenolides. Conversely, soil N had little impact on A. curassavica biomass, but was the primary driver of increasing N-cardenolides (voruscharin, uscharin and their precursors); caterpillar induction was weak. Butterfly enzyme assays revealed damage-induced cardenolides substantially increased toxicity of both milkweeds to both butterflies, swamping out effects of soil N on cardenolide concentration and composition. Although these two milkweed species differentially responded to soil N with allocation to growth and specific cardenolides, leaf toxicity to butterfly Na+/K+-ATPases was primarily driven by herbivore-induced defense. Thus, both biotic and abiotic factors shape the composition of phytochemical defense expression, and their relative importance may be dictated by plant life-history differences.

Environmental drivers of annual population fluctuations in a trans-Saharan insect migrant.

Many latitudinal insect migrants including agricultural pests, disease vectors, and beneficial species show huge fluctuations in the year-to-year abundance of spring immigrants reaching temperate zones. It is widely believed that this variation is driven by climatic conditions in the winter-breeding regions, but evidence is lacking. We identified the environmental drivers of the annual population dynamics of a cosmopolitan migrant butterfly (the painted lady Vanessa cardui) using a combination of long-term monitoring and climate and atmospheric data within the western part of its Afro-Palearctic migratory range. Our population models show that a combination of high winter NDVI (normalized difference vegetation index) in the Savanna/Sahel of sub-Saharan Africa, high spring NDVI in the Maghreb of North Africa, and frequent favorably directed tailwinds during migration periods are the three most important drivers of the size of the immigration to western Europe, while our atmospheric trajectory simulations demonstrate regular opportunities for wind-borne trans-Saharan movements. The effects of sub-Saharan vegetative productivity and wind conditions confirm that painted lady populations on either side of the Sahara are linked by regular mass migrations, making this the longest annual insect migration circuit so far known. Our results provide a quantification of the environmental drivers of large annual population fluctuations of an insect migrant and hold much promise for predicting invasions of migrant insect pests, disease vectors, and beneficial species.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Cardenolides, toxicity, and the costs of sequestration in the coevolutionary interaction between monarchs and milkweeds.

For highly specialized insect herbivores, plant chemical defenses are often co-opted as cues for oviposition and sequestration. In such interactions, can plants evolve novel defenses, pushing herbivores to trade off benefits of specialization with costs of coping with toxins? We tested how variation in milkweed toxins (cardenolides) impacted monarch butterfly (Danaus plexippus) growth, sequestration, and oviposition when consuming tropical milkweed (Asclepias curassavica), one of two critical host plants worldwide. The most abundant leaf toxin, highly apolar and thiazolidine ring-containing voruscharin, accounted for 40% of leaf cardenolides, negatively predicted caterpillar growth, and was not sequestered. Using whole plants and purified voruscharin, we show that monarch caterpillars convert voruscharin to calotropin and calactin in vivo, imposing a burden on growth. As shown by in vitro experiments, this conversion is facilitated by temperature and alkaline pH. We next employed toxin-target site experiments with isolated cardenolides and the monarch's neural Na+/K+-ATPase, revealing that voruscharin is highly inhibitory compared with several standards and sequestered cardenolides. The monarch's typical >50-fold enhanced resistance to cardenolides compared with sensitive animals was absent for voruscharin, suggesting highly specific plant defense. Finally, oviposition was greatest on intermediate cardenolide plants, supporting the notion of a trade-off between benefits and costs of sequestration for this highly specialized herbivore. There is apparently ample opportunity for continued coevolution between monarchs and milkweeds, although the diffuse nature of the interaction, due to migration and interaction with multiple milkweeds, may limit the ability of monarchs to counteradapt.

In vivo visualization of butterfly scale cell morphogenesis in Vanessa cardui.

During metamorphosis, the wings of a butterfly sprout hundreds of thousands of scales with intricate microstructures and nano-structures that determine the wings' optical appearance, wetting characteristics, thermodynamic properties, and aerodynamic behavior. Although the functional characteristics of scales are well known and prove desirable in various applications, the dynamic processes and temporal coordination required to sculpt the scales' many structural features remain poorly understood. Current knowledge of scale growth is primarily gained from ex vivo studies of fixed scale cells at discrete time points; to fully understand scale formation, it is critical to characterize the time-dependent morphological changes throughout their development. Here, we report the continuous, in vivo, label-free imaging of growing scale cells of Vanessa cardui using speckle-correlation reflection phase microscopy. By capturing time-resolved volumetric tissue data together with nanoscale surface height information, we establish a morphological timeline of wing scale formation and gain quantitative insights into the underlying processes involved in scale cell patterning and growth. We identify early differences in the patterning of cover and ground scales on the young wing and quantify geometrical parameters of growing scale features, which suggest that surface growth is critical to structure formation. Our quantitative, time-resolved in vivo imaging of butterfly scale development provides the foundation for decoding the processes and biomechanical principles involved in the formation of functional structures in biological materials.

Outcomes of endoscopic modified cartilage-perichondrium tympanoplasty on different sized perforations.

OBJECTIVES: To evaluate the anatomic and functional outcomes of type1 tympanoplasty with endoscopic modified butterfly cartilage-perichondrium technique. METHODS: In our modification, perichondrium was elevated circumferentially till the attached part of the composite graft was approximately same size and shape of the perforation, cartilage was trimmed based on the perforation but 0.5 mm larger. Cartilage portion of the graft was placed medial to the edge of the perforation, then perichondrium was rolled out and draped on the circumferential raw surface of remaining tympanic membrane around. RESULTS: At 4 months postop, the anatomic integrity rate of the tympanic membrane perforation for small & medium sized perforation and large sized perforation group were 100 % and 94 % (p > 0.05). For the small & medium perforation group, the mean pre and 4 months postop ACs were 30 ± 8 dB and 18 ± 6.4 dB (p < 0.01). The mean pre and 4 months postop ABGs were 19 ± 11 dB and 9 ± 3 dB (p < 0.01). For the large perforation group, the mean pre and 4 months postop ACs were 43 ± 12.5 dB and 21.5 ± 7 dB (p < 0.01). The mean pre and 4 months postop ABGs were 34 ± 8.5 dB and 12.5 ± 6 dB (p < 0.01). The differences of mean 4 months postop ACs and mean 4 months postop ABGs between the two groups were not significant (p > 0.05). CONCLUSIONS: Compared to the conventional inlay butterfly cartilage tympanoplasty technique, large or marginal perforations can be sealed more securely by this modification.

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe.

BACKGROUND: Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored. RESULTS: Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies. CONCLUSIONS: Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

aristaless1 has a dual role in appendage formation and wing color specification during butterfly development.

BACKGROUND: Highly diverse butterfly wing patterns have emerged as a powerful system for understanding the genetic basis of phenotypic variation. While the genetic basis of this pattern variation is being clarified, the precise developmental pathways linking genotype to phenotype are not well understood. The gene aristaless, which plays a role in appendage patterning and extension, has been duplicated in Lepidoptera. One copy, aristaless1, has been shown to control a white/yellow color switch in the butterfly Heliconius cydno, suggesting a novel function associated with color patterning and pigmentation. Here we investigate the developmental basis of al1 in embryos, larvae, and pupae using new antibodies, CRISPR/Cas9, RNAi, qPCR assays of downstream targets, and pharmacological manipulation of an upstream activator. RESULTS: We find that Al1 is expressed at the distal tips of developing embryonic appendages consistent with its ancestral role. In developing wings, we observe Al1 accumulation within developing scale cells of white H. cydno during early pupation while yellow scale cells exhibit little Al1 at this time point. Reduced Al1 expression is also associated with yellow scale development in al1 knockouts and knockdowns. We propose that Al1 expression in future white scales might be related to an observed downregulation of the enzyme Cinnabar and other genes that synthesize and transport the yellow pigment, 3-hydroxykynurenine (3-OHK). Finally, we provide evidence that Al1 activation is under the control of Wnt signaling. CONCLUSIONS: We propose a model in which high levels of Al1 during early pupation, which are mediated by Wnt, are important for melanic pigmentation and specifying white portions of the wing while reduced levels of Al1 during early pupation promote upregulation of proteins needed to move and synthesize 3-OHK, promoting yellow pigmentation. In addition, we discuss how the ancestral role of aristaless in appendage extension may be relevant in understanding the cellular mechanism behind color patterning in the context of the heterochrony hypothesis.