BackgroundThe border ocellus, or eyespot, is a conspicuous color pattern element in butterfly wings. For two decades, it has been hypothesized that transcription factors such as Distal-less (Dll) are responsible for eyespot pattern development in butterfly wings, based on their expression in the prospective eyespots. In particular, it has been suggested that Dll is a determinant for eyespot size. However, functional evidence for this hypothesis has remained incomplete, due to technical difficulties.ResultsHere, we show that ectopically expressed Dll induces ectopic elemental color patterns in the adult wings of the blue pansy butterfly, Junonia orithya (Lepidoptera, Nymphalidae). Using baculovirus-mediated gene transfer, we misexpressed Dll protein fused with green fluorescent protein (GFP) in pupal wings, resulting in ectopic color patterns, but not the formation of intact eyespots. Induced changes included clusters of black and orange scales (a basic feature of eyespot patterns), black and gray scales, and inhibition of cover scale development. In contrast, ectopic expression of GFP alone did not induce any color pattern changes using the same baculovirus-mediated gene transfer system.ConclusionsThese results suggest that Dll plays an instructive role in the development of color pattern elements in butterfly wings, although Dll alone may not be sufficient to induce a complete eyespot. This study thus experimentally supports the hypothesis of Dll function in eyespot development.

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. Population genetics theory has been developed and tested to explain the rapid evolution of sex chromosomes and sex determination. However, we do not know why the sex chromosomes are divergent in some taxa and conserved in others. Addressing this question requires comparing closely related taxa with conserved and divergent sex chromosomes to identify biological features that could explain these differences. Cytological karyotypes suggest that muscid flies (e.g., house fly) and blow flies are such a taxonomic pair. The sex chromosomes appear to differ across muscid species, whereas they are conserved across blow flies. Despite the cytological evidence, we do not know the extent to which muscid sex chromosomes are independently derived along different evolutionary lineages. To address that question, we used genomic and transcriptomic sequence data to identify young sex chromosomes in two closely related muscid species, horn fly (Haematobia irritans) and stable fly (Stomoxys calcitrans). We provide evidence that the nascent sex chromosomes of horn fly and stable fly were derived independently from each other and from the young sex chromosomes of the closely related house fly (Musca domestica). We present three different scenarios that could have given rise to the sex chromosomes of horn fly and stable fly, and we describe how the scenarios could be distinguished. Distinguishing between these scenarios in future work could identify features of muscid genomes that promote sex chromosome divergence.

It is often desirable but difficult to retrieve information on the mature phenotype of an immature tissue sample that has been subjected to gene expression analysis. This problem cannot be ignored when individual variation within a species is large. To circumvent this problem in the butterfly wing system, we developed a new surgical method for removing a single forewing from a pupa using Junonia orithya; the operated pupa was left to develop to an adult without eclosion. The removed right forewing was subjected to gene expression analysis, whereas the non-removed left forewing was examined for color patterns. As a test case, we focused on Distal-less (Dll), which likely plays an active role in inducing elemental patterns, including eyespots. The Dll expression level in forewings was paired with eyespot size data from the same individual. One third of the operated pupae survived and developed wing color patterns. Dll expression levels were significantly higher in males than in females, although male eyespots were smaller in size than female eyespots. Eyespot size data showed weak but significant correlations with the Dll expression level in females. These results demonstrate that a single-wing removal method was successfully applied to the butterfly wing system and suggest the weak and non-exclusive contribution of Dll to eyespot size determination in this butterfly. Our novel methodology for establishing correspondence between gene expression and phenotype can be applied to other candidate genes for color pattern development in butterflies. Conceptually similar methods may also be applicable in other developmental systems.

Selection pressures can vary within localized areas and across massive geographical scales. Temperature is one of the best studied ecologically variable abiotic factors that can affect selection pressures across multiple spatial scales. Organisms rely on physiological (thermal tolerance) and behavioral (thermal preference) mechanisms to thermoregulate in response to environmental temperature. In addition, spatial heterogeneity in temperatures can select for local adaptation in thermal tolerance, thermal preference, or both. However, the concordance between thermal tolerance and preference across genotypes and sexes within species and across populations is greatly understudied. The house fly, Musca domestica, is a well‐suited system to examine how genotype and environment interact to affect thermal tolerance and preference. Across multiple continents, house fly males from higher latitudes tend to carry the male‐determining gene on the Y chromosome, whereas those from lower latitudes usually have the male determiner on the third chromosome. We tested whether these two male‐determining chromosomes differentially affect thermal tolerance and preference as predicted by their geographical distributions. We identify effects of genotype and developmental temperature on male thermal tolerance and preference that are concordant with the natural distributions of the chromosomes, suggesting that temperature variation across the species range contributes to the maintenance of the polymorphism. In contrast, female thermal preference is bimodal and largely independent of congener male genotypes. These sexually dimorphic thermal preferences suggest that temperature‐dependent mating dynamics within populations could further affect the distribution of the two chromosomes. Together, the differences in thermal tolerance and preference across sexes and male genotypes suggest that different selection pressures may affect the frequencies of the male‐determining chromosomes across different spatial scales.

Sex determination, the developmental process by which sexually dimorphic phenotypes are established, evolves fast. Evolutionary turnover in a sex determination pathway may occur via selection on alleles that are genetically linked to a new master sex determining locus on a newly formed proto‐sex chromosome. Species with polygenic sex determination, in which master regulatory genes are found on multiple different proto‐sex chromosomes, are informative models to study the evolution of sex determination and sex chromosomes. House flies are such a model system, with male determining loci possible on all six chromosomes and a female‐determiner on one of the chromosomes as well. The two most common male‐determining proto‐Y chromosomes form latitudinal clines on multiple continents, suggesting that temperature variation is an important selection pressure responsible for maintaining polygenic sex determination in this species. Temperature‐dependent fitness effects could be manifested through temperature‐dependent gene expression differences across proto‐Y chromosome genotypes. These gene expression differences may be the result of cis regulatory variants that affect the expression of genes on the proto‐sex chromosomes, or trans effects of the proto‐Y chromosomes on genes elswhere in the genome. We used RNA‐seq to identify genes whose expression depends on proto‐Y chromosome genotype and temperature in adult male house flies. We found no evidence for ecologically meaningful temperature‐dependent expression differences of sex determining genes between male genotypes, but we were probably not sampling an appropriate developmental time‐point to identify such effects. In contrast, we identified many other genes whose expression depends on the interaction between proto‐Y chromosome genotype and temperature, including genes that encode proteins involved in reproduction, metabolism, lifespan, stress response, and immunity. Notably, genes with genotype‐by‐temperature interactions on expression were not enriched on the proto‐sex chromosomes. Moreover, there was no evidence that temperature‐dependent expression is driven by chromosome‐wide cis‐regulatory divergence between the proto‐Y and proto‐X alleles. Therefore, if temperature‐dependent gene expression is responsible for differences in phenotypes and fitness of proto‐Y genotypes across house fly populations, these effects are driven by a small number of temperature‐dependent alleles on the proto‐Y chromosomes that may have trans effects on the expression of genes on other chromosomes.