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.

The protein-coding regions of eukaryotic genes are fragmented into exons that, like the genes within which they are situated, can be duplicated, deleted, or reorganized. Cataloging and organizing within-gene exon similarities is necessary for a systematic study of exon evolution and its consequences. To facilitate the study of exon duplications, we present Exonize, a computational tool that identifies and classifies coding exon duplications in annotated genomes. Exonize implements a graph-based framework to handle clusters of related exons resulting from repeated rounds of exon duplication. The interdependence between duplicated exons or groups of exons across transcripts is classified. By identifying duplication events between exonic and intronic regions, Exonize can detect unannotated or degenerate exons. To aid in data parsing and downstream analysis, the Python module exonize_analysis is provided. The application of Exonize to 20 eukaryote genomes identifies full-exon duplications in at least 4% of vertebrate genes, with more than 900 human genes having a full-exon duplication event.

To date, the use of CRISPR/Cas9 technology in ecological-model species for validating genotype to phenotype connections has focused primarily on visual phenotypes using G0 mutations, coupled with analyses of resulting mosaic phenotypes. However, studies of physiological phenotypes necessitate germline mutations in order to assess non-visible phenotypic effects, and thus, dedicated efforts to develop efficient germline mutations in ecological model species are needed. Here, we applied the CRISPR/Cas9 technology to an ecological model species, the speckled wood butterfly (Pararge aegeria). We targeted yellow-y, which is required for the production of black melanin, as yellow-y loss of function (LOF) mutations are not lethal and easy to phenotype, affording efficient assessment of G0 and germline mutations. To explore what factors may affect the efficiency of transformation, we employed four alternative treatments, including variation in sgRNAs and their concentrations. Color changes in the head capsule of first larval instars, as well as adult wing color, were used as indicators of successful knockouts. Individuals with wings that were at least 50% transformed were mated, with their G1 offspring assessed for the presence of germline mutations. Our CRISPR/Cas9 technique was highly efficient at generating LOF mutations in yellow-y. Across all treatments, nearly 80% of adults exhibited mosaic LOF phenotypes, of which nearly 30% appeared to have 100% LOF phenotypes. Crosses between adults exhibiting at least 50% LOF phenotypes resulted in fully transformed offspring, revealing a high incidence of germline LOF mutations in yellow-y. We provide a detailed protocol on how to obtain high germline LOF mutation efficiency in order to advance the study of genotype–phenotype connections for non-visible physiological traits across natural populations of this and other model ecological species.

Parasitoid wasps are major causes of mortality of many species, making host immune defences a common target of adaptive evolution, though such targets outside model species are poorly understood. In this study, we used two tests of positive selection to compare across three closely related Galerucella leaf beetles that show substantial differences in their phenotypic response to the shared parasitoid wasp Asecodes parviclava, their main natural enemy. Using a codon-based test, which detects excess amino acid fixations per locus along each species’ lineage, we found more evidence of positive selection on parasitoid-relevant immune genes in the species with the strongest immunocompetence (G. pusilla) compared with the species having weaker immunocompetence (G. tenella and G. calmariensis). Moreover, genes coding for the early phases in the immune response cascade were predominantly among the positively selected immune genes, providing targets for future functional genomic study to pin-point connections between genotypic and phenotypic differences in defences towards a parasitoid wasp. In contrast, genome-wide analyses of the haplotype frequency spectrum, which quantify selection over recent evolutionary time scales, revealed similar signatures of positive selection on immune genes across species. These results advance the field of host-parasitoid dynamics by providing novel insights into the tempo and mode of insect host evolutionary dynamics, and offering a framework for making genotype to phenotype connections for immunocompetence phenotypes.

Evolutionary changes in geographic distribution and larval host plants may promote the rapid diversification of montane insects, but this scenario has been rarely investigated. We studied the rapid radiation of the butterfly genus Colias, which has diversified in mountain ecosystems in Eurasia, Africa, and the Americas. Based on a data set of 150 nuclear protein-coding genetic loci and mitochondrial genomes, we constructed a time-calibrated phylogenetic tree of Colias species with broad taxon sampling. We then inferred their ancestral geographic ranges, historical diversification rates, and the evolution of host use. We found that the most recent common ancestor of Colias was likely geographically widespread and originated ~3.5 Ma. The group subsequently diversified in different regions across the world, often in tandem with geographic expansion events. No aspect of elevation was found to have a direct effect on diversification. The genus underwent a burst of diversification soon after the divergence of the Neotropical lineage, followed by an exponential decline in diversification rate toward the present. The ancestral host repertoire included the legume genera Astragalus and Trifolium but later expanded to include a wide range of Fabaceae genera and plants in more distantly related families, punctuated with periods of host range expansion and contraction. We suggest that the widespread distribution of the ancestor of all extant Colias lineages set the stage for diversification by isolation of populations that locally adapted to the various different environments they encountered, including different host plants. In this scenario, elevation is not the main driver but might have accelerated diversification by isolating populations.

In cases of recurrent colonizations of similar habitats from the same base population, it is commonly expected that repeated phenotypic adaptation is caused by parallel changes in genetic variation. However, it is becoming increasingly clear that similar phenotypic variation may also evolve by alternative genetic pathways. Here, we explore the repeated evolution of photoperiodic plasticity for diapause induction across Swedish populations of the speckled wood butterfly, Pararge aegeria. This species has colonized Scandinavia at least twice, and population genomic results show that one of the candidate regions associated with spatial variation in photoperiodism is situated on the Z-chromosome. Here, we assay hybrid crosses between several populations that differ in photoperiodic plasticity for sex-linked inheritance of the photoperiodic reaction norm. We find that while a cross between more distantly related populations from the two different colonization events shows strong sex-dependent inheritance of photoperiodic plasticity, a cross between two more closely related populations within the oldest colonization range shows no such effect. We conclude that the genotype–phenotype map for photoperiodic plasticity varies across these populations and that similar local phenotypic adaptation has evolved during recurrent colonization events by partly non-parallel genetic changes.

Quantifying the tempo and mode via modern phylogenetic comparative methods can provide key insights into how selection and constraints shape trait evolution on a macroevolutionary time scale. Here, we elucidate the evolution of hibernation (winter) diapause, a complex and defining life-history trait that allows temporal escape from harsh winters in temperate regions for many insects, including our model system, butterflies. Butterflies can diapause in all major life stages, and the availability of global-scale phylogenies makes them an ideal model system for studying diapause evolution. First, using a thorough literature survey, we scored the developmental stage of hibernation diapause (egg, larva, pupa, adult) vs. absence of diapause. We find that larval diapause is most common, while pupal, egg, and adult diapause are relatively rare. Next, we determined that the loss of diapause occurred at a much higher rate and that gains primarily occurred from the non-diapause state. While ancestral state estimation at deeper nodes remained uncertain, we found consistent patterns for some families and strong evidence for extensive convergence in diapause evolution. Contrary to expectations, we find no support for increased gain of diapause during the Eocene–Oligocene glaciation (~35 million years ago). Overall, the evolution of diapause in butterflies has a complex history, has evolved convergently, and has likely predated the major glaciation event consistent with the deep history of diapause evolution in insects. This study advances our understanding of the evolution of a complex and important life-history trait and establishes a macroevolutionary foundation for future studies on the ultimate and proximate basis of diapause evolution.

The neritid snail Theodoxus fluviatilis is found across habitats differing in salinity, from shallow waters along the coast of the Baltic Sea to lakes throughout Europe. Living close to the water surface makes this species vulnerable to changes in salinity in their natural habitat, and the lack of a free-swimming larval stage limits this species’ dispersal. Together, these factors have resulted in a patchy distribution of quite isolated populations differing in their salinity tolerances. In preparation for investigating the mechanisms underlying the physiological differences in osmoregulation between populations that cannot be explained solely by phenotypic plasticity, we present here an annotated draft genome assembly for T. fluviatilis, generated using PacBio long reads, Illumina short reads, and transcriptomic data. While the total assembly size (1045 kb) is similar to those of related species, it remains highly fragmented (N scaffolds = 35,695; N50 = 74 kb) though moderately high in complete gene content (BUSCO single copy complete: 74.3%, duplicate: 2.6%, fragmented: 10.6%, missing: 12.5% using metazoa n = 954). Nevertheless, we were able to generate gene annotations of 21,220 protein-coding genes (BUSCO single copy complete: 65.1%, duplicate: 16.7%, fragmented: 9.1%, missing: 9.1% using metazoa n = 954). Not only will this genome facilitate comparative evolutionary studies across Gastropoda, as this is the first genome assembly for the basal snail family Neritidae, it will also greatly facilitate the study of salinity tolerance in this species. Additionally, we discuss the challenges of working with a species where high molecular weight DNA isolation is very difficult. 

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Comparative analyses of gene birth–death dynamics have the potential to reveal gene families that played an important role in the evolution of morphological, behavioral, or physiological variation. Here, we used whole genomes of 30 species of butterflies and moths to identify gene birth–death dynamics among the Lepidoptera that are associated with specialist or generalist feeding strategies. Our work advances this field using a uniform set of annotated proteins for all genomes, investigating associations while correcting for phylogeny, and assessing all gene families rather than a priori subsets. We discovered that the sizes of several important gene families (e.g. those associated with pesticide resistance, xenobiotic detoxification, and/or protein digestion) are significantly correlated with diet breadth. We also found 22 gene families showing significant shifts in gene birth–death dynamics at the butterfly (Papilionoidea) crown node, the most notable of which was a family of pheromone receptors that underwent a contraction potentially linked with a shift to visual-based mate recognition. Our findings highlight the importance of uniform annotations, phylogenetic corrections, and unbiased gene family analyses in generating a list of candidate genes that warrant further exploration.