Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Little is known about how alternative splicing and other post-transcriptional molecular mechanisms impact plasticity. Steward et al. use transcriptomic and genomic data from the butterfly Bicyclus anynana, finding that splicing plasticity is likely to experience genetic constraints.

Vertically transmitted rhabdoviruses are found across three insect families and have dynamic interactions with their hosts

A small number of free-living viruses have been found to be obligately vertically transmitted, but it remains uncertain how widespread vertically transmitted viruses are and how quickly they can spread through host populations. Recent metagenomic studies have found several insects to be infected with sigma viruses (Rhabdoviridae). Here, we report that sigma viruses that infect Mediterranean fruit flies (Ceratitis capitata), Drosophila immigrans, and speckled wood butterflies (Pararge aegeria) are all vertically transmitted. We find patterns of vertical transmission that are consistent with those seen in Drosophila sigma viruses, with high rates of maternal transmission, and lower rates of paternal transmission. This mode of transmission allows them to spread rapidly in populations, and using viral sequence data we found the viruses in D. immigrans and C. capitata had both recently swept through host populations. The viruses were common in nature, with mean prevalences of 12% in C. capitata, 38% in D. immigrans and 74% in P. aegeria We conclude that vertically transmitted rhabdoviruses may be widespread in a broad range of insect taxa, and that these viruses can have dynamic interactions with their hosts.

Differentiation in putative male sex pheromone components across and within populations of the African butterfly Bicyclus anynana as a potential driver of reproductive isolation

Sexual traits are often the most divergent characters among closely related species, suggesting an important role of sexual traits in speciation. However, to prove this, we need to show that sexual trait differences accumulate before or during the speciation process, rather than being a consequence of it. Here, we contrast patterns of divergence among putative male sex pheromone (pMSP) composition and the genetic structure inferred from variation in the mitochondrial cytochrome oxidase 1 and nuclear CAD loci in the African butterfly Bicyclus anynana (Butler, 1879) to determine whether the evolution of "pheromonal dialects" occurs before or after the differentiation process. We observed differences in abundance of some shared pMSP components as well as differences in the composition of the pMSP among B. anynana populations. In addition, B. anynana individuals from Kenya displayed differences in the pMSP composition within a single population that appeared not associated with genetic differences. These differences in pMSP composition both between and within B. anynana populations were as large as those found between different Bicyclus species. Our results suggest that "pheromonal dialects" evolved within and among populations of B. anynana and may therefore act as precursors of an ongoing speciation process.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Butterflies and moths (Lepidoptera) vary in chromosome number. Here, the authors sequence the genome of the Glanville fritillary butterfly, Melitaea cinxia, show it has the ancestral lepidopteran karyotype and provide insight into how chromosomal fusions have shaped karyotype evolution in butterflies and moths.

Mitochondrial DNA signature for range-wide populations of Bicyclus anynana suggests a rapid expansion from recent refugia

This study investigates the genetic diversity, population structure and demographic history of the afrotropical butterfly Bicyclus anynana using mitochondrial DNA (mtDNA). Samples from six wild populations covering most of the species range from Uganda to South Africa were compared for the cytochrome c oxidase subunit gene (COI). Molecular diversity indices show overall high mtDNA diversity for the populations, but low nucleotide divergence between haplotypes. Our results indicate relatively little geographic population structure among the southern populations, especially given the extensive distributional range and an expectation of limited gene flow between populations. We implemented neutrality tests to assess signatures of recent historical demographic events. Tajima's D test and Fu's F(S) test both suggested recent population growth for the populations. The results were only significant for the southernmost populations when applying Tajima's D, but Fu's F(S) indicated significant deviations from neutrality for all populations except the one closest to the equator. Based on our own findings and those from pollen and vegetation studies, we hypothesize that the species range of B. anynana was reduced to equatorial refugia during the last glacial period, and that the species expanded southwards during the past 10.000 years. These results provide crucial background information for studies of phenotypic and molecular adaptation in wild populations of B. anynana.

Translating environmental gradients into discontinuous reaction norms via hormone signalling in a polyphenic butterfly

Polyphenisms-the expression of discrete phenotypic morphs in response to environmental variation-are examples of phenotypic plasticity that may potentially be adaptive in the face of predictable environmental heterogeneity. In the butterfly Bicyclus anynana, we examine the hormonal regulation of phenotypic plasticity that involves divergent developmental trajectories into distinct adult morphs for a suite of traits as an adaptation to contrasting seasonal environments. This polyphenism is induced by temperature during development and mediated by ecdysteroid hormones. We reared larvae at separate temperatures spanning the natural range of seasonal environments and measured reaction norms for ecdysteroids, juvenile hormones (JHs) and adult fitness traits. Timing of peak ecdysteroid, but not JH titres, showed a binary response to the linear temperature gradient. Several adult traits (e.g. relative abdomen mass) responded in a similar, dimorphic manner, while others (e.g. wing pattern) showed a linear response. This study demonstrates that hormone dynamics can translate a linear environmental gradient into a discrete signal and, thus, that polyphenic differences between adult morphs can already be programmed at the stage of hormone signalling during development. The range of phenotypic responses observed within the suite of traits indicates both shared regulation and independent, trait-specific sensitivity to the hormone signal.