A polygenic explanation for Haldane's rule in butterflies.

Two robust rules have been discovered about animal hybrids: Heterogametic hybrids are more unfit (Haldane's rule), and sex chromosomes are disproportionately involved in hybrid incompatibility (the large-X/Z effect). The exact mechanisms causing these rules in female heterogametic taxa such as butterflies are unknown but are suggested by theory to involve dominance on the sex chromosome. We investigate hybrid incompatibilities adhering to both rules in Papilio and Heliconius butterflies and show that dominance theory cannot explain our data. Instead, many defects coincide with unbalanced multilocus introgression between the Z chromosome and all autosomes. Our polygenic explanation predicts both rules because the imbalance is likely greater in heterogametic females, and the proportion of introgressed ancestry is more variable on the Z chromosome. We also show that mapping traits polygenic on a single chromosome in backcrosses can generate spurious large-effect QTLs. This mirage is caused by statistical linkage among polygenes that inflates estimated effect sizes. By controlling for statistical linkage, most incompatibility QTLs in our hybrid crosses are consistent with a polygenic basis. Since the two genera are very distantly related, polygenic hybrid incompatibilities are likely common in butterflies.

Cell death and wing reduction during the metamorphosis of sex-specific flightless morphs in winter geometrid moths.

Winter geometrid moths show striking sexual dimorphism by having female-specific flightless morphs. The evolutionary grades of wing reduction in winter geometrid moths vary and range from having short wings, vestigial wings, to being wingless. Although the ontogenetic processes underlying the development of the wingless or short-wing morphs in Lepidoptera has been well studied, the mechanisms underlying the development of vestigial wing morphs in winter geometrid moths during metamorphosis are poorly understood. In the winter geometrid moth Sebastosema bubonaria Warren, 1896, the males have functional wings, but the females have vestigial wings. Here, we studied the ontogenetic processes affecting wing reduction in the winter geometrid moth S. bubonaria using light microscopy and transmission electron microscopy, and compared the ontogenetic process of wing reduction in this species with that in another species of the wingless-female winter moth that we investigated previously. Our results showed that, in the vestigial-wing morphs, the loss of pupal wing epithelium was terminated in the middle of the wing degeneration process, whereas in the wingless morph, the pupal wing epithelium disappeared almost completely and the final appearance of the wings differed slightly among flightless morphs. We propose that the extent and location of cell death in the pupal wing play an important role in the various patterns of reduced wings that are observed in flightless moths.

Simple and complex, sexually dimorphic retinal mosaic of fritillary butterflies.

Butterflies have variable sets of spectral photoreceptors that underlie colour vision. The photoreceptor organization may be optimized for the detection of body coloration. Fritillaries (Argynnini) are nymphalid butterflies exhibiting varying degrees of sexual dimorphism in wing coloration. In two sister species, the females have orange (Argynnis paphia) and dark wings (Argynnis sagana), respectively, while the males of both species have orange wings with large patches of pheromone-producing androconia. In spite of the differences in female coloration, the eyes of both species exhibit an identical sexual dimorphism. The female eyeshine is uniform yellow, while the males have a complex retinal mosaic with yellow and red-reflecting ommatidia. We found the basic set of ultraviolet-, blue- and green-peaking photoreceptors in both sexes. Males additionally have three more photoreceptor classes, peaking in green, yellow and red, respectively. The latter is the basal R9, indirectly measured through hyperpolarizations in the green-peaking R1-2. In many nymphalid tribes, including the closely related Heliconiini, the retinal mosaic is complex in both sexes. We hypothesize that the simple mosaic of female Argynnini is a secondary reduction, possibly driven by the use of olfaction for intraspecific recognition, whereas vision remains the primary sense for the task in the males. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

Admixture of evolutionary rates across a butterfly hybrid zone.

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

First Annotated Genome of a Mandibulate Moth, Neomicropteryx cornuta, Generated Using PacBio HiFi Sequencing.

We provide a new, annotated genome assembly of Neomicropteryx cornuta, a species of the so-called mandibulate archaic moths (Lepidoptera: Micropterigidae). These moths belong to a lineage that is thought to have split from all other Lepidoptera more than 300 Ma and are consequently vital to understanding the early evolution of superorder Amphiesmenoptera, which contains the order Lepidoptera (butterflies and moths) and its sister order Trichoptera (caddisflies). Using PacBio HiFi sequencing reads, we assembled a highly contiguous genome with a contig N50 of nearly 17 Mb. The assembled genome length of 541,115,538 bp is about half the length of the largest published Amphiesmenoptera genome (Limnephilus lunatus, Trichoptera) and double the length of the smallest (Papilio polytes, Lepidoptera). We find high recovery of universal single copy orthologs with 98.1% of BUSCO genes present and provide a genome annotation of 15,643 genes aided by resolved isoforms from PacBio IsoSeq data. This high-quality genome assembly provides an important resource for studying ecological and evolutionary transitions in the early evolution of Amphiesmenoptera.

Zoogeographical barriers causing discontinuous osteometrical variations in the northern treeshrew skulls.

Morphological variation of the skull was examined in the northern treeshrew (Tupaia belangeri) from various localities across Southeast Asia. Through a multivariate analysis, the treeshrews from South Vietnam exhibited distinct morphological characteristics compared to other populations from Thailand and Laos, and Malaysia. The plots of the specimens of North Vietnam are not randomly mixed with Thailand plots segregation in the scatteregrams of canonical discriminant analysis. Since the skulls of the population from North Vietnam were morphologically similar to those form central Laos and northern and northeastern Thailand, the zoogeographical barrier effect of Mekong River was not clearly confirmed. The population of the Kanchanaburi in western Thailand is clearly smaller in size compared to the other populations. The southern border of the distribution of this species is determined by the Isthmus of Kra or Kangar-Pattani Line. In the northern treeshrew, which is distributed from southern China to Bangladesh and southern Thailand, we have detected osteometrical geographical variation driven by geography. These results indicate that the skull morphology in the Tupaia glis-belangeri complex distinctively differs in South Vietnam, western Thailand, and southern Thailand. The zoogeographical barrier and factor separating these districts are expected to clarify in the future.

Molecular phylogeny and historical biogeography of Parnara butterflies (Lepidoptera: Hesperiidae).

The butterfly genus Parnara (Hesperiinae: Baorini), of which some are major pests of economic crops (e.g., rice, wild rice stems and sugarcane), currently consists of 10 species and several subspecies and has a highly disjunct distribution in Australia, Africa, and Asia. We determined the systematic relationships and biogeographical history of the genus by reconstructing the phylogeny based on eight genes and 101 specimens representing all 10 recognized species. Four species delimitation methods (ABGD, bPTP, GMYC and BPP) were also employed to assess the taxonomic status of each species. Based on these results and analyses, we recognize 11 extant species in the genus. The status of the taxon P. naso poutieri (Boisduval, 1833) from Madagascar is revised as a distinct species, Parnara poutieri (Boisduval, 1833) stat. rev. The subspecies P. guttata mangala (Moore, 1866) syn. nov. is synonymized with P. guttata guttata (Bremer & Grey, 1853), while P. bada (Moore, 1878) is provisionally treated as a complex of two species, namely P. bada and P. apostata (Snellen, 1886). The monophyly of Parnara is strongly supported, with the following relationships: P. amalia + ((P. monasi + (P. poutieri + P. naso)) + ((P. kawazoei + P. bada complex) + (P. ganga + (P. ogasawarensis + (P. guttata + P. batta))))). Divergence time and ancestral range estimates indicate that the common ancestor of Parnara originated in an implausible area of Australia, Africa, and Oriental region in the mid-Oligocene and then differentiated in the late Miocene-late Pliocene. Dispersal and range expansion have played an important role in diversification of the genus in Asia and Afica. Relatively stable geotectonic plates at the time when most extant lineages appeared during the late Miocene-early Pliocene might have been the factor responsible for the relatively constant low dynamic rate of diversification within the group.

Anchored phylogenomics illuminates the skipper butterfly tree of life.

BACKGROUND: Butterflies (Papilionoidea) are perhaps the most charismatic insect lineage, yet phylogenetic relationships among them remain incompletely studied and controversial. This is especially true for skippers (Hesperiidae), one of the most species-rich and poorly studied butterfly families. METHODS: To infer a robust phylogenomic hypothesis for Hesperiidae, we sequenced nearly 400 loci using Anchored Hybrid Enrichment and sampled all tribes and more than 120 genera of skippers. Molecular datasets were analyzed using maximum-likelihood, parsimony and coalescent multi-species phylogenetic methods. RESULTS: All analyses converged on a novel, robust phylogenetic hypothesis for skippers. Different optimality criteria and methodologies recovered almost identical phylogenetic trees with strong nodal support at nearly all nodes and all taxonomic levels. Our results support Coeliadinae as the sister group to the remaining skippers, the monotypic Euschemoninae as the sister group to all other subfamilies but Coeliadinae, and the monophyly of Eudaminae plus Pyrginae. Within Pyrginae, Celaenorrhinini and Tagiadini are sister groups, the Neotropical firetips, Pyrrhopygini, are sister to all other tribes but Celaenorrhinini and Tagiadini. Achlyodini is recovered as the sister group to Carcharodini, and Erynnini as sister group to Pyrgini. Within the grass skippers (Hesperiinae), there is strong support for the monophyly of Aeromachini plus remaining Hesperiinae. The giant skippers (Agathymus and Megathymus) once classified as a subfamily, are recovered as monophyletic with strong support, but are deeply nested within Hesperiinae. CONCLUSIONS: Anchored Hybrid Enrichment sequencing resulted in a large amount of data that built the foundation for a new, robust evolutionary tree of skippers. The newly inferred phylogenetic tree resolves long-standing systematic issues and changes our understanding of the skipper tree of life. These resultsenhance understanding of the evolution of one of the most species-rich butterfly families.

Species richness of Eurasian Zephyrus hairstreaks (Lepidoptera: Lycaenidae: Theclini) with implications on historical biogeography: An NDM/VNDM approach.

AIM: A database based on distributional records of Eurasian Zephyrus hairstreaks (Lepidoptera: Lycaenidae: Theclini) was compiled to analyse their areas of endemism (AoEs), species richness and distribution patterns, to explore their locations of past glacial refugia and dispersal routes. METHODS: Over 2000 Zephyrus hairstreaks occurrences are analysed using the NDM/VNDM algorithm, for the recognition of AoEs. Species richness was calculated by using the option 'Number of different classes' to count the different classes of a variable presented in each 3.0°×3.0° grid cell, and GIS software was used to visualize distribution patterns of endemic species. RESULTS: Centres of species richness of Zephyrus hairstreaks are situated in the eastern Qinghai-Tibet Plateau (EQTP), Hengduan Mountain Region (HDMR) and the Qinling Mountain Region (QLMR). Latitudinal gradients in species richness show normal distribution with the peak between 25° N and 35° N in the temperate zone, gradually decreasing towards the poles. Moreover, most parts of central and southern China, especially the area of QLMR-EQTP-HDMR, were identified as AoEs that may have played a significant role as refugia during Quaternary global cooling. There are four major distributional patterns of Zephyrus hairstreaks in Eurasia: Sino-Japanese, Sino-Himalayan, high-mountain and a combined distribution covering all three patterns. CONCLUSIONS: Zephyrus hairstreaks probably originated at least 23-24 Myr ago in E. Asia between 25° N to 35° N in the temperate zone. Cenozoic orogenies caused rapid speciation of this tribe and extrusion of the Indochina block resulted in vicariance between the Sino-Japanese and the Sino-Himalayan patterns. The four distribution patterns provided two possible dispersal directions: Sino-Japanese dispersal and Sino-Himalayan dispersal.

Curious Oviposition Behavior in Phyllium westwoodii (Phasmatodea: Phylliidae): Preliminary Observations.

We report that in a leaf insect, Phyllium westwoodii Wood-Mason (Phasmatodea: Phylliidae), two differing apertures can be used for oviposition, the color of eggs being affected by which aperture is used. Eggs which are forcibly propelled from the internal space within the valvulae of the abdomen are brown, whereas white eggs emerge slowly from the opening between the eighth sternite and the valvulae, and are deposited close to the ventral surface of the female. This unusual oviposition system does not appear to have been previously reported in phasmatids or in other insects.

Theclini butterflies from Weixi, China, with description of two new species (Lepidoptera: Lycaenidae).

Thirteen species belonging to nine genera of the tribe Theclini are recorded from Weixi county, Yunnan province, China. Among which Noseozephyrus Koiwaya and Shirozuozephyrus kirbariensis machimurai Koiwaya are reported as a new record genus and a new record species for China, respectively, and two new species, Noseozephyrus lisus sp. n. and Shirozuozephyrus fibonacci sp. n., are described. All the species and their male genitalia are illustrated. A key for the genus of Shirozuozephyrus is provided.

Phylogeography of Fischer's blue, Tongeia fischeri, in Japan: Evidence for introgressive hybridization.

The widespread lycaenid butterfly Tongeia fischeri is distributed from eastern Europe to northeastern Asia and represented by three geographically isolated populations in Japan. In order to clarify the phylogeographic history of the species, we used sequences of three mitochondrial (COI, Cyt b and ND5) and two nuclear (Rpl5 and Ldh) genes of 207 individuals collected from 55 sites throughout Japan and five sites on the Asian continent. Phylogenetic trees and the median-joining network revealed six evolutionary mitochondrial haplotype clades, which corresponded to the geographic distribution of the species. Common ancestors of Japanese T. fischeri might have come to Japan during the mid-Pleistocene by multiple dispersals of continental populations, probably via a land bridge or narrow channel between western Japan and the Korean Peninsula. The geographical patterns of variation of mitochondrial and nuclear markers are discordant in northeastern Kyushu, possibly as a result of introgressive hybridization during the ancient contact between the Kyushu and Shikoku populations in the last glacial maximum. The phylogeographic pattern of T. fischeri in Japan are probably related to the geological history, Pleistocene climatic oscillations and distribution of the host plant.

Molecular phylogeny, laboratory rearing, and karyotype of the bombycid moth, Trilocha varians.

This study describes the molecular phylogeny, laboratory rearing, and karyotype of a bombycid moth, Trilocha varians (F. Walker) (Lepidoptera: Bombycidae), which feeds on leaves of Ficus spp. (Rosales: Moraceae). The larvae of this species were collected in Taipei city, Taiwan, and the Ryukyu Archipelago (Ishigaki and Okinawa Islands, Japan). Molecular phylogenetic analyses revealed that T. varians belongs to the subfamily Bombycinae, thus showing a close relationship to the domesticated silkworm Bombyx mori (L.), a lepidopteran model insect. A laboratory method was developed for rearing T. varians and the time required for development from the embryo to adult was determined. From oviposition to adult emergence, the developmental zero was 10.47 °C and total effective temperature was 531.2 day-degrees, i.e., approximately 30 days for one generation when reared at 28 °C. The haploid of T. varians consisted of n = 26 chromosomes. In highly polyploid somatic nuclei, females showed a large heterochromatin body, indicating that the sex chromosome system in T. varians is WZ/ZZ (female/male). The results of the present study should facilitate the utilization of T. varians as a reference species for B. mori, thereby leading to a greater understanding of the ecology and evolution of bombycid moths.

Female sex pheromone and male behavioral responses of the bombycid moth Trilocha varians: comparison with those of the domesticated silkmoth Bombyx mori.

Analysis of female sex pheromone components and subsequent field trap experiments demonstrated that the bombycid moth Trilocha varians uses a mixture of (E,Z)-10,12-hexadecadienal (bombykal) and (E,Z)-10,12-hexadecadienyl acetate (bombykyl acetate) as a sex pheromone. Both of these components are derivatives of (E,Z)-10,12-hexadecadienol (bombykol), the sex pheromone of the domesticated silkmoth Bombyx mori. This finding prompted us to compare the antennal and behavioral responses of T. varians and B. mori to bombykol, bombykal, and bombykyl acetate in detail. The antennae of T. varians males responded to bombykal and bombykyl acetate but not to bombykol, and males were attracted only when lures contained both bombykal and bombykyl acetate. In contrast, the antennae of B. mori males responded to all the three components. Behavioral analysis showed that B. mori males responded to neither bombykal nor bombykyl acetate. Meanwhile, the wing fluttering response of B. mori males to bombykol was strongly inhibited by bombykal and bombykyl acetate, thereby indicating that bombykal and bombykyl acetate act as behavioral antagonists for B. mori males. T. varians would serve as a reference species for B. mori in future investigations into the molecular mechanisms underlying the evolution of sex pheromone communication systems in bombycid moths.