Wolbachia strain diversity in a complex group of sympatric cryptic parasitoid wasp species.

BACKGROUND: Maternally-inherited symbionts can induce pre-mating and/or post-mating reproductive isolation between sympatric host lineages, and speciation, by modifying host reproductive phenotypes. The large parasitoid wasp genus Cotesia (Braconidae) includes a diversity of cryptic species, each specialized in parasitizing one to few related Lepidoptera host species. Here, we characterized the infection status of an assemblage of 21 Cotesia species from 15 countries by several microbial symbionts, as a first step toward investigating whether symbionts may provide a barrier to gene flow between these parasitoid host lineages. RESULTS: The symbiotic microbes Arsenophonus, Cardinium, Microsporidium and Spiroplasma were not detected in the Cotesia wasps. However, the endosymbiotic bacterium Wolbachia was present in at least eight Cotesia species, and hence we concentrated on it upon screening additional DNA extracts and SRAs from NCBI. Some of the closely related Cotesia species carry similar Wolbachia strains, but most Wolbachia strains showed patterns of horizontal transfer between phylogenetically distant host lineages. CONCLUSIONS: The lack of co-phylogenetic signal between Wolbachia and Cotesia suggests that the symbiont and hosts have not coevolved to an extent that would drive species divergence between the Cotesia host lineages. However, as the most common facultative symbiont of Cotesia species, Wolbachia may still function as a key-player in the biology of the parasitoid wasps. Its precise role in the evolution of this complex clade of cryptic species remains to be experimentally investigated.

Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies.

Sex chromosomes play an outsized role in adaptation and speciation, and thus deserve particular attention in evolutionary genomics. In particular, fusions between sex chromosomes and autosomes can produce neo-sex chromosomes, which offer important insights into the evolutionary dynamics of sex chromosomes. Here, we investigate the evolutionary origin of the previously reported Danaus neo-sex chromosome within the tribe Danaini. We assembled and annotated genomes of Tirumala septentrionis (subtribe Danaina), Ideopsis similis (Amaurina), Idea leuconoe (Euploeina) and Lycorea halia (Itunina) and identified their Z-linked scaffolds. We found that the Danaus neo-sex chromosome resulting from the fusion between a Z chromosome and an autosome corresponding to the Melitaea cinxia chromosome (McChr) 21 arose in a common ancestor of Danaina, Amaurina and Euploina. We also identified two additional fusions as the W chromosome further fused with the synteny block McChr31 in I. similis and independent fusion occurred between ancestral Z chromosome and McChr12 in L. halia. We further tested a possible role of sexually antagonistic selection in sex chromosome turnover by analysing the genomic distribution of sex-biased genes in I. leuconoe and L. halia. The autosomes corresponding to McChr21 and McChr31 involved in the fusions are significantly enriched in female- and male-biased genes, respectively, which could have hypothetically facilitated fixation of the neo-sex chromosomes. This suggests a role of sexual antagonism in sex chromosome turnover in Lepidoptera. The neo-Z chromosomes of both I. leuconoe and L. halia appear fully compensated in somatic tissues, but the extent of dosage compensation for the ancestral Z varies across tissues and species.

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.

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.

Active around the year: Butterflies and moths adapt their life cycles to a warming world.

Living in a warming world requires adaptations to altered annual temperature regimes. In Europe, spring is starting earlier, and the vegetation period is ending later in the year. These climatic changes are leading not only to shifts in distribution ranges of flora and fauna, but also to phenological shifts. Using long-term observation data of butterflies and moths collected during the past decades across northern Austria, we test for phenological shifts over time and changes in the number of generations. On average, Lepidoptera adults emerged earlier in the year and tended to extend their flight periods in autumn. Many species increased the annual number of generations. These changes were more pronounced at lower altitudes than at higher altitudes, leading to an altered phenological zonation. Our findings indicate that climate change does not only affect community composition but also the life history of insects. Increased activity and reproductive periods might alter Lepidoptera-host plant associations and food webs.

Roadside habitat: Boon or bane for pollinating insects?

Pollinators, which provide vital services to wild ecosystems and agricultural crops, are facing global declines and habitat loss. As undeveloped land becomes increasingly scarce, much focus has been directed recently to roadsides as potential target zones for providing floral resources to pollinators. Roadsides, however, are risky places for pollinators, with threats from vehicle collisions, toxic pollutants, mowing, herbicides, and more. Although these threats have been investigated, most studies have yet to quantify the costs and benefits of roadsides to pollinators and, therefore, do not address whether the costs outweigh the benefits for pollinator populations using roadside habitats. In this article, we address how, when, and under what conditions roadside habitats may benefit or harm pollinators, reviewing existing knowledge and recommending practical questions that managers and policymakers should consider when planning pollinator-focused roadside management.

Regional differences in thermoregulation between two European butterfly communities.

Understanding how different organisms cope with changing temperatures is vital for predicting future species' distributions and highlighting those at risk from climate change. As ectotherms, butterflies are sensitive to temperature changes, but the factors affecting butterfly thermoregulation are not fully understood. We investigated which factors influence thermoregulatory ability in a subset of the Mediterranean butterfly community. We measured adult thoracic temperature and environmental temperature (787 butterflies; 23 species) and compared buffering ability (defined as the ability to maintain a consistent body temperature across a range of air temperatures) and buffering mechanisms to previously published results from Great Britain. Finally, we tested whether thermoregulatory ability could explain species' demographic trends in Catalonia. The sampled sites in each region differ climatically, with higher temperatures and solar radiation but lower wind speeds in the Catalan sites. Both butterfly communities show nonlinear responses to temperature, suggesting a change in behaviour from heat-seeking to heat avoidance at approximately 22°C. However, the communities differ in the use of buffering mechanisms, with British populations depending more on microclimates for thermoregulation compared to Catalan populations. Contrary to the results from British populations, we did not find a relationship between region-wide demographic trends and butterfly thermoregulation, which may be due to the interplay between thermoregulation and the habitat changes occurring in each region. Thus, although Catalan butterfly populations seem to be able to thermoregulate successfully at present, evidence of heat avoidance suggests this situation may change in the future.

Influence of Host Plants and Tending Ants on the Cuticular Hydrocarbon Profile of a Generalist Myrmecophilous Caterpillar.

In myrmecophilous organisms, which live in symbiosis with ants, cuticular hydrocarbons (CHCs) play a pivotal role in interspecific communication and defense against chemical-oriented predators. Although these interactions form complex information webs, little is known about the influence of biotic environmental factors on the CHC profiles of myrmecophiles. Here, we analyzed the effect of different host plants and tending ants on the larval CHC profile of Synargis calyce (Lepidoptera: Riodinidae), a polyphagous species with facultative myrmecophily. Groups of caterpillars were fed individually with three host plant species (without tending ants), and with two tending ant species. Through gas chromatography analysis, we compared the cuticular profiles of treatments and found a high similarity between plants and caterpillars (65-82%), but a low similarity between caterpillars and their tending ants (30 - 25%). Cluster analysis showed that caterpillars, ants, and plants form distinct groups, indicating that S. calyce caterpillars have their own chemical profile. These results are similar to those observed for Lycaenidae caterpillars indicating that there is functional convergence in the chemical strategies used by myrmecophilous caterpillar species with similar ecology. Also, the results suggest that the cuticular compounds of S. calyce are primarily influenced by their host plants rather than their tending ants. Thus, we propose that these caterpillars present a trade-off between camouflage and directly informing their presence to ants, maintaining their unique chemical profile, though slightly affected by biotic environmental factors.

The decline of butterflies in Europe: Problems, significance, and possible solutions.

We review changes in the status of butterflies in Europe, focusing on long-running population data available for the United Kingdom, the Netherlands, and Belgium, based on standardized monitoring transects. In the United Kingdom, 8% of resident species have become extinct, and since 1976 overall numbers declined by around 50%. In the Netherlands, 20% of species have become extinct, and since 1990 overall numbers in the country declined by 50%. Distribution trends showed that butterfly distributions began decreasing long ago, and between 1890 and 1940, distributions declined by 80%. In Flanders (Belgium), 20 butterflies have become extinct (29%), and between 1992 and 2007 overall numbers declined by around 30%. A European Grassland Butterfly Indicator from 16 European countries shows there has been a 39% decline of grassland butterflies since 1990. The 2010 Red List of European butterflies listed 38 of the 482 European species (8%) as threatened and 44 species (10%) as near threatened (note that 47 species were not assessed). A country level analysis indicates that the average Red List rating is highest in central and mid-Western Europe and lowest in the far north of Europe and around the Mediterranean. The causes of the decline of butterflies are thought to be similar in most countries, mainly habitat loss and degradation and chemical pollution. Climate change is allowing many species to spread northward while bringing new threats to susceptible species. We describe examples of possible conservation solutions and a summary of policy changes needed to conserve butterflies and other insects.

A switch to feeding on cycads generates parallel accelerated evolution of toxin tolerance in two clades of Eumaeus caterpillars (Lepidoptera: Lycaenidae).

We assembled a complete reference genome of Eumaeus atala, an aposematic cycad-eating hairstreak butterfly that suffered near extinction in the United States in the last century. Based on an analysis of genomic sequences of Eumaeus and 19 representative genera, the closest relatives of Eumaeus are Theorema and Mithras We report natural history information for Eumaeus, Theorema, and Mithras Using genomic sequences for each species of Eumaeus, Theorema, and Mithras (and three outgroups), we trace the evolution of cycad feeding, coloration, gregarious behavior, and other traits. The switch to feeding on cycads and to conspicuous coloration was accompanied by little genomic change. Soon after its origin, Eumaeus split into two fast evolving lineages, instead of forming a clump of close relatives in the phylogenetic tree. Significant overlap of the fast evolving proteins in both clades indicates parallel evolution. The functions of the fast evolving proteins suggest that the caterpillars developed tolerance to cycad toxins with a range of mechanisms including autophagy of damaged cells, removal of cell debris by macrophages, and more active cell proliferation.

Mutualistic interactions shape global spatial congruence and climatic niche evolution in Neotropical mimetic butterflies.

Understanding the mechanisms underlying species distributions and coexistence is both a priority and a challenge for biodiversity hotspots such as the Neotropics. Here, we highlight that Müllerian mimicry, where defended prey species display similar warning signals, is key to the maintenance of biodiversity in the c. 400 species of the Neotropical butterfly tribe Ithomiini (Nymphalidae: Danainae). We show that mimicry drives large-scale spatial association among phenotypically similar species, providing new empirical evidence for the validity of Müller's model at a macroecological scale. Additionally, we show that mimetic interactions drive the evolutionary convergence of species climatic niche, thereby strengthening the co-occurrence of co-mimetic species. This study provides new insights into the importance of mutualistic interactions in shaping both niche evolution and species assemblages at large spatial scales. Critically, in the context of climate change, our results highlight the vulnerability to extinction cascades of such adaptively assembled communities tied by positive interactions.

Widespread Gene Expression Divergence in Butterfly Sensory Tissues Plays a Fundamental Role During Reproductive Isolation and Speciation.

Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation.

Additive and interactive effects of anthropogenic stressors on an insect herbivore.

The pressures of global change acting on wild plants and animals include exposure to environmental toxins, the introduction of non-native species, and climate change. Relatively few studies have been reported in which these three main classes of stressors have been examined simultaneously, allowing for the possibility of synergistic effects in an experimental context. In this study, we exposed caterpillars of the Melissa blue butterfly (Lycaeides melissa) to three concentrations of chlorantraniliprole, under three experimental climates, on a diet of a native or a non-native host plant throughout larval development in a fully factorial experiment. We find that high pesticide exposure and a non-native diet exhibit strong negative effects on caterpillars, resulting in 62% and 42% reduction in survival, respectively, while interactive effects tend to be weaker, ranging from 15% to 22% reduction in survival. Interactive effects have been shown to be strong in other contexts, but do not appear to be universal; however, our study shows that the cumulative effects of stressors acting in isolation (additively) are sufficiently strong to severely reduce survival and by extension population persistence in the wild.

Diversification linked to larval host plant in the butterfly Eumedonia eumedon.

It is widely accepted that the relationship between phytophagous insects and their host plants influences insect diversification. However, studies addressed at documenting host-associated genetic differentiation (HAD) and the mechanisms that drive reproductive isolation in host-associated lineages (or host races) are still scarce relative to insect diversity. To uncover further evidence on the HAD processes in Lepidoptera, we investigated the genetic structure of the geranium argus butterfly (Eumedonia eumedon) and tested for isolation by ecology (IBE) vs. isolation by distance (IBD). Genomic data revealed an array of host races (three of them in the same mountain range, the Cantabrian Mountains, northern Iberia) at apparently distinct levels of reproductive isolation. We found a pattern of IBE mediated by HAD at both local and European scales, in which genetic differentiation between populations and individuals correlated significantly with the taxonomic relatedness of the host plants. IBD was significant only when considered at the wider European scale. We hypothesize that, locally, HAD between Geranium-feeding populations was caused (at least partially) by allochrony, that is via adaptation of adult flight time to the flowering period of each host plant species. Nevertheless, the potential reproductive isolation between populations using Erodium and populations using Geranium cannot be explained by allochrony or IBD, and other mechanisms are expected to be at play.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae).

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Mechanisms underpinning community stability along a latitudinal gradient: Insights from a niche-based approach.

At large scales, the mechanisms underpinning stability in natural communities may vary in importance due to changes in species composition, mean abundance, and species richness. Here we link species characteristics (niche positions) and community characteristics (richness and abundance) to evaluate the importance of stability mechanisms in 156 butterfly communities monitored across three European countries and spanning five bioclimatic regions. We construct niche-based hierarchical structural Bayesian models to explain first differences in abundance, population stability, and species richness between the countries, and then explore how these factors impact community stability both directly and indirectly (via synchrony and population stability). Species richness was partially explained by the position of a site relative to the niches of the species pool, and species near the centre of their niche had higher average population stability. The differences in mean abundance, population stability, and species richness then influenced how much variation in community stability they explained across the countries. We found, using variance partitioning, that community stability in Finnish communities was most influenced by community abundance, whereas this aspect was unimportant in Spain with species synchrony explaining most variation; the UK was somewhat intermediate with both factors explaining variation. Across all countries, the diversity-stability relationship was indirect with species richness reducing synchrony which increased community stability, with no direct effects of species richness. Our results suggest that in natural communities, biogeographical variation observed in key drivers of stability, such as population abundance and species richness, leads to community stability being limited by different factors and that this can partially be explained due to the niche characteristics of the European butterfly assemblage.

Challenges and opportunities for using natural history collections to estimate insect population trends.

Natural history collections (NHC) provide a wealth of information that can be used to understand the impacts of global change on biodiversity. As such, there is growing interest in using NHC data to estimate changes in species' distributions and abundance trends over historic time horizons when contemporary survey data are limited or unavailable. However, museum specimens were not collected with the purpose of estimating population trends and thus can exhibit spatiotemporal and collector-specific biases that can impose severe limitations to using NHC data for evaluating population trajectories. Here we review the challenges associated with using museum records to track long-term insect population trends, including spatiotemporal biases in sampling effort and sparse temporal coverage within and across years. We highlight recent methodological advancements that aim to overcome these challenges and discuss emerging research opportunities. Specifically, we examine the potential of integrating museum records and other contemporary data sources (e.g. collected via structured, designed surveys and opportunistic citizen science programs) in a unified analytical framework that accounts for the sampling biases associated with each data source. The emerging field of integrated modelling provides a promising framework for leveraging the wealth of collections data to accurately estimate long-term trends of insect populations and identify cases where that is not possible using existing data sources.

Wings are not perfect: increased wing asymmetry in a tropical butterfly as a response to forest fragmentation.

Habitat fragmentation and ecosystem changes have the potential to affect animal populations in different ways. To effectively monitor these changes, biomonitoring tools have been developed and applied to detect changes in population structure and/or individual traits that reflect such changes. Fluctuating asymmetry (FA) represents random deviations from perfect symmetry in bilateral traits from perfect symmetry in response to genetic and/or environmental stresses. In this study, we evaluated the use of FA as a tool to monitor stress caused by forest fragmentation and edge formation, using the tropical butterfly M. helenor (Nymphalidae) as a model species. We collected adult butterflies from three fragments of Atlantic Forest in Brazil encompassing both edge and interior habitats. Four wing traits (wing length, wing width, ocelli area, and ocelli diameter) were evaluated. Butterflies captured at edge sites exhibited higher FA values for wing length and wing width compared to those captured at interior sites, whereas traits related to ocelli did not show differences between the two habitat types. Our results suggest that the differences in abiotic and biotic conditions between forest interior and edges can act as a source of stress, impacting the symmetry of flight-related traits. On the other hand, as ocelli are crucial for butterfly camouflage and counter-predator strategies, our results indicate that this trait may be more conserved. By employing FA, we identified trait-specific responses to habitat fragmentation, thus suggesting its potential as a biomarker for environmental stress that can be used in butterflies to monitor habitat quality and change.

The Basis of Navigation Across Species.

Animals navigate a wide range of distances, from a few millimeters to globe-spanning journeys of thousands of kilometers. Despite this array of navigational challenges, similar principles underlie these behaviors across species. Here, we focus on the navigational strategies and supporting mechanisms in four well-known systems: the large-scale migratory behaviors of sea turtles and lepidopterans as well as navigation on a smaller scale by rats and solitarily foraging ants. In lepidopterans, rats, and ants we also discuss the current understanding of the neural architecture which supports navigation. The orientation and navigational behaviors of these animals are defined in terms of behavioral error-reduction strategies reliant on multiple goal-directed servomechanisms. We conclude by proposing to incorporate an additional component into this system: the observation that servomechanisms operate on oscillatory systems of cycling behavior. These oscillators and servomechanisms comprise the basis for directed orientation and navigational behaviors.

Infrared optical and thermal properties of microstructures in butterfly wings.

While surface microstructures of butterfly wings have been extensively studied for their structural coloration or optical properties within the visible spectrum, their properties in infrared wavelengths with potential ties to thermoregulation are relatively unknown. The midinfrared wavelengths of 7.5 to 14 µm are particularly important for radiative heat transfer in the ambient environment, because of the overlap with the atmospheric transmission window. For instance, a high midinfrared emissivity can facilitate surface cooling, whereas a low midinfrared emissivity can minimize heat loss to surroundings. Here we find that the midinfrared emissivity of butterfly wings from warmer climates such as Archaeoprepona demophoon (Oaxaca, Mexico) and Heliconius sara (Pichincha, Ecuador) is up to 2 times higher than that of butterfly wings from cooler climates such as Celastrina echo (Colorado) and Limenitis arthemis (Florida), using Fourier-transform infrared (FTIR) spectroscopy and infrared thermography. Our optical computations using a unit cell approach reproduce the spectroscopy data and explain how periodic microstructures play a critical role in the midinfrared. The emissivity spectrum governs the temperature of butterfly wings, and we demonstrate that C. echo wings heat up to 8 °C more than A. demophoon wings under the same sunlight in the clear sky of Irvine, CA. Furthermore, our thermal computations show that butterfly wings in their respective habitats can maintain a moderate temperature range through a balance of solar absorption and infrared emission. These findings suggest that the surface microstructures of butterfly wings potentially contribute to thermoregulation and provide an insight into butterflies' survival.