Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F-ST) between the two introduced populations exceeds that of the naturally sympatric populations (F-ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ (pi over bar approximate to 0.002 and pi over bar approximate to 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.

Insects rely on their innate immune system to successfully mediate complex interactions with their microbiota, as well as the microbes present in the environment. Previous work has shown that components of the canonical immune gene repertoire evolve rapidly and have evolutionary characteristics originating from interactions with fast-evolving microorganisms. Although these interactions are likely to vary among populations, there is a poor understanding of the microevolutionary dynamics of immune genes, especially in non-Dipteran insects. Here, we use the full set of canonical insect immune genes to investigate microevolutionary dynamics acting on these genes between and among populations by comparing three allopatric populations of the green-veined white butterfly, Pieris napi (Linne; Lepidoptera, Pieridae). Immune genes showed increased genetic diversity compared to genes from the rest of the genome and various functional categories exhibited different types of signatures of selection, at different evolutionary scales, presenting a complex pattern of selection dynamics. Signatures of balancing selection were identified in 10 genes, and 17 genes appear to be under positive selection. Genes involved with the cellular arm of the immune response as well as the Toll pathway appear to be enriched among our outlier loci, regardless of functional category. This suggests that the targets of selection might focus upon an entire pathway, rather than functional subsets across pathways. Our microevolutionary results are similar to previously observed macroevolutionary patterns from diverse taxa, suggesting that either the immune system is robust to dramatic differences in life history and microbial communities, or that diverse microbes exert similar selection pressures.

Interactions between herbivorous insects and their host plants are a central component of terrestrial food webs and a critical topic in agriculture, where a substantial fraction of potential crop yield is lost annually to pests. Important insights into plant-insect interactions have come from research on specific plant defences and insect detoxification mechanisms. Yet, much remains unknown about the molecular mechanisms that mediate plant-insect interactions. Here we use multiple genome-wide approaches to map the molecular basis of herbivory from both plant and insect perspectives, focusing on butterflies and their larval host plants. Parallel genome-wide association studies in the cabbage white butterfly, Pieris rapae, and its host plant, Arabidopsis thaliana, pinpointed a small number of butterfly and plant genes that influenced herbivory. These genes, along with much of the genome, were regulated in a dynamic way over the time course of the feeding interaction. Comparative analyses, including diverse butterfly/plant systems, showed a variety of genome-wide responses to herbivory, as well as a core set of highly conserved genes in butterflies as well as their host plants. These results greatly expand our understanding of the genomic causes and evolutionary consequences of ecological interactions across two of nature's most diverse taxa, butterflies and flowering plants.

While large-scale genomic approaches are increasingly revealing the genetic basis of polymorphic phenotypes such as colour morphs, such approaches are almost exclusively conducted in species with high-quality genomes and annotations. Here, we use Pool-Seq data for both genome assembly and SNP frequency estimation, followed by scanning for F-ST outliers to identify divergent genomic regions. Using paired-end, short-read sequencing data from two groups of individuals expressing divergent phenotypes, we generate a de novo rough-draft genome, identify SNPs and calculate genomewide F-ST differences between phenotypic groups. As genomes generated by Pool-Seq data are highly fragmented, we also present an approach for super-scaffolding contigs using existing protein-coding data sets. Using this approach, we reanalysed genomic data from two recent studies of birds and butterflies investigating colour pattern variation and replicated their core findings, demonstrating the accuracy and power of a Pool-Seq-only approach. Additionally, we discovered new regions of high divergence and new annotations that together suggest novel parallels between birds and butterflies in the origins of their colour pattern variation.

The global loss of biodiversity continues at an alarming rate. Genomic approaches have been suggested as a promising tool for conservation practice as scaling up to genome-wide data can improve traditional conservation genetic inferences and provide qualitatively novel insights. However, the generation of genomic data and subsequent analyses and interpretations remain challenging and largely confined to academic research in ecology and evolution. This generates a gap between basic research and applicable solutions for conservation managers faced with multifaceted problems. Before the real-world conservation potential of genomic research can be realized, we suggest that current infrastructures need to be modified, methods must mature, analytical pipelines need to be developed, and successful case studies must be disseminated to practitioners.

Numerous studies have found amino acid variation at the phosphoglucose isomerase (PGI) gene associated with organismal performance and fitness. Here we focus upon recent advances in the study of this gene, highlighting novel species being studied, new tools being used, and emerging insights into the evolutionary dynamics acting on this gene. Our synthesis highlights questions that are coming into focus, as well as the need for attention in specific areas, such as manipulative experiments to establish mechanistic insights and a causative role of allelic variation.