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.

Significance Factors influencing Wolbachia transfer into new species remain poorly understood. This is important as Wolbachia can influence speciation and is being developed as a novel arthropod-borne disease control approach. We show the native microbiota of Anopheles impede vertical transmission of Wolbachia . Antibiotic microbiome perturbation enables Wolbachia transmission in two Anopheles species. Mosquitoes with altered microbiomes do not exhibit blood meal-induced mortality associated with Wolbachia infection, suggesting that mosquitoes are killed by interactions between Wolbachia and other bacteria present in the mosquito. We identified Asaia as the bacterium responsible for inhibiting Wolbachia transmission, and partially responsible for blood meal-induced mortality. These results suggest that microbial interactions profoundly affect the host, and that microbiome incompatibility may influence distribution of Wolbachia in arthropods.

Male-killing bacteria are widespread in arthropods, and can profoundly alter the reproductive biology of their host species. Here we detail the first case of complete suppression of a male killer. The nymphalid butterfly Hypolimnas bolina is infected with a strain of the bacterium Wolbachia, wBol1, which kills male host embryos in Polynesian populations, but does not do so in many areas of Southeast Asia, where both males and female adults are naturally infected, and wBol1-infected females produce a 1:1 sex ratio. We demonstrate that absence of male killing by wBol1 is associated with dominant zygotic suppression of the action of the male killer. Simulations demonstrate host suppressors of male-killer action can spread very rapidly, and historical data indicating the presence of male killing in Southeast Asia in the very recent past suggests suppressor spread has been a very recent occurrence. Thus, male killer/host interactions are much more dynamic than previously recognised, with rapid and dramatic loss of the phenotype. Our results also indicate that suppression can render male killers completely quiescent, leading to the conclusion that some species that do not currently express a male killer may have done so in the past, and thus that more species have had their biology affected by these parasites than previously believed.

Microbial interactions are an underappreciated force in shaping insect microbiome communities. Although pairwise patterns of symbiont interactions have been identified, we have a poor understanding regarding the scale and the nature of co-occurrence and co-exclusion interactions within the microbiome. To characterize these patterns in mosquitoes, we sequenced the bacterial microbiome of Aedes aegypti, Ae. albopictus, and Culex quinquefasciatus caught in the field or reared in the laboratory and used these data to generate interaction networks. For collections, we used traps that attracted host-seeking or ovipositing female mosquitoes to determine how physiological state affects the microbiome under field conditions. Interestingly, we saw few differences in species richness or microbiome community structure in mosquitoes caught in either trap. Co-occurrence and co-exclusion analysis identified 116 pairwise interactions substantially increasing the list of bacterial interactions observed in mosquitoes. Networks generated from the microbiome of Ae. aegypti often included highly interconnected hub bacteria. There were several instances where co-occurring bacteria co-excluded a third taxa, suggesting the existence of tripartite relationships. Several associations were observed in multiple species or in field and laboratory-reared mosquitoes indicating these associations are robust and not influenced by environmental or host factors. To demonstrate that microbial interactions can influence colonization of the host, we administered symbionts to Ae. aegypti larvae that either possessed or lacked their resident microbiota. We found that the presence of resident microbiota can inhibit colonization of particular bacterial taxa. Our results highlight that microbial interactions in mosquitoes are complex and influence microbiome composition.

The ratio of males to females in a species is often considered to be relatively constant, at least over ecological time. Hamilton noted that the spread of "selfish" sex ratio-distorting elements could be rapid and produce a switch to highly biased population sex ratios. Selection against a highly skewed sex ratio should promote the spread of mutations that suppress the sex ratio distortion. We show that in the butterfly Hypolimnas bolina the suppression of sex biases occurs extremely fast, with a switch from a 100:1 population sex ratio to 1:1 occurring in fewer than 10 generations.

Sex-ratio distorters are found in numerous species and can reach high frequencies within populations. Here, we address the compelling, but poorly tested, hypothesis that the sex ratio bias caused by such elements profoundly alters their host's mating system. We compare aspects of female and male reproductive biology between island populations of the butterfly Hypolimnas bolina that show varying degrees of female bias, because of a male-killing Wolbachia infection. Contrary to expectation, female bias leads to an increase in female mating frequency, up to a point where male mating capacity becomes limiting. We show that increased female mating frequency can be explained as a facultative response to the depleted male mating resources in female biased populations. In other words, this system is one where male-killing bacteria trigger a vicious circle of increasing male fatigue and female promiscuity.

Background: The interaction between the Blue Moon butterfly, Hypolimnas bolina, and Wolbachia has attracted interest because of the high prevalence of male-killing achieved within the species, the ecological consequences of this high prevalence, the intensity of selection on the host to suppress the infection, and the presence of multiple Wolbachia infections inducing different phenotypes. We examined diversity in the co-inherited marker, mtDNA, and the partitioning of this between individuals of different infection status, as a means to investigate the population biology and evolutionary history of the Wolbachia infections.

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.