Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

Production and characterization of Trichoderma asperellum chitinases and their use in synergy with Bacillus thuringiensis for lepidopteran control

BACKGROUND

Despite their known negative effects on ecosystems and human health, synthetic pesticides are still largely used to control crop insect pests. Currently, the biopesticide market for insect biocontrol mainly relies on the entomopathogenic bacterium Bacillus thuringiensis (Bt). New biocontrol tools for crop protection might derive from fungi, in particular from Trichoderma spp., which are known producers of chitinases and other bioactive compounds able to negatively affect insect survival.

RESULTS

In this study, we first developed an environmentally sustainable production process for obtaining chitinases from Trichoderma asperellum ICC012. Then, we investigated the biological effects of this chitinase preparation - alone or in combination with a Bt-based product - when orally administered to two lepidopteran species. Our results demonstrate that T. asperellum efficiently produces a multi-enzymatic cocktail able to alter the chitin microfibril network of the insect peritrophic matrix, resulting in delayed development and larval death. The co-administration of T. asperellum chitinases and sublethal concentrations of Bt toxins increased larval mortality. This synergistic effect was likely due to the higher amount of Bt toxins that passed the damaged peritrophic matrix and reached the target receptors on the midgut cells of chitinase-treated insects.

CONCLUSION

Our findings may contribute to the development of an integrated pest management technology based on fungal chitinases that increase the efficacy of Bt-based products, mitigating the risk of Bt-resistance development. © 2024 Society of Chemical Industry.

De Novo Genome Assembly and Annotation of Leptosia nina Provide New Insights into the Evolutionary Dynamics of Genes Involved in Host-Plant Adaptation of Pierinae Butterflies

In interactions between plants and herbivorous insects, the traits enabling phytophagous insects to overcome chemical defenses of their host plants have evolved multiple times. A prominent example of such adaptive key innovations in herbivorous insects is nitrile specifier proteins (NSPs) that enabled Pierinae butterflies to colonize Brassicales host plants that have a glucosinolate-myrosinase defense system. Although the evolutionary aspects of NSP-encoding genes have been studied in some Pierinae taxa (especially among Pieris butterflies), the ancestral evolutionary state of NSPs is unclear due to the limited genomic information available for species within Pierinae. Here, we generate a high-quality genome assembly and annotation of Leptosia nina, a member of a small tribe, Leptosiaini. L. nina uses as its main host Capparaceae plants, one of the ancestral hosts within Pierinae. By using ∼90-fold coverage of Oxford Nanopore long reads and Illumina short reads for subsequent polishing and error correction, we constructed a final genome assembly that consisted of 286 contigs with a total of 225.8 Mb and an N50 of 10.7 Mb. Genome annotation with transcriptome hints predicted 16,574 genes and covered 98.3% of BUSCO genes. A typical NSP gene is composed of three tandem domains found in Pierinae butterflies; unexpectedly, we found a new NSP-like gene in Pierinae composed of only two tandem domains. This newly found NSP-like gene in L. nina provides important insights into the evolutionary dynamics of domain and gene duplication events relating to host-plant adaptation in Pierinae butterflies.

Host genetics and larval host plant modulate microbiome structure and evolution underlying the intimate insect-microbe-plant interactions in Parnassius species on the Qinghai-Tibet Plateau

Insects harbor a remarkable diversity of gut microbiomes critical for host survival, health, and fitness, but the mechanism of this structured symbiotic community remains poorly known, especially for the insect group consisting of many closely related species that inhabit the Qinghai-Tibet Plateau. Here, we firstly analyzed population-level 16S rRNA microbial dataset, comprising 11 Parnassius species covering 5 subgenera, from 14 populations mostly sampled in mountainous regions across northwestern-to-southeastern China, and meanwhile clarified the relative importance of multiple factors on gut microbial community structure and evolution. Our findings indicated that both host genetics and larval host plant modulated gut microbial diversity and community structure. Moreover, the effect analysis of host genetics and larval diet on gut microbiomes showed that host genetics played a critical role in governing the gut microbial beta diversity and the symbiotic community structure, while larval host plant remarkably influenced the functional evolution of gut microbiomes. These findings of the intimate insect-microbe-plant interactions jointly provide some new insights into the correlation among the host genetic background, larval host plant, the structure and evolution of gut microbiome, as well as the mechanisms of high-altitude adaptation in closely related species of this alpine butterfly group.

Integrating learning into animal range dynamics under rapid human-induced environmental change

Human-induced rapid environmental change (HIREC) is creating environments deviating considerably from natural habitats in which species evolved. Concurrently, climate warming is pushing species' climatic envelopes to geographic regions that offer novel ecological conditions. The persistence of species is likely affected by the interplay between the degree of ecological novelty and phenotypic plasticity, which in turn may shape an organism's range-shifting ability. Current modelling approaches that forecast animal ranges are characterized by a static representation of the relationship between habitat use and fitness, which may bias predictions under conditions imposed by HIREC. We argue that accounting for dynamic species-resource relationships can increase the ecological realism of range shift predictions. Our rationale builds on the concepts of ecological fitting, the process whereby individuals form successful novel biotic associations based on the suite of traits they carry at the time of encountering the novel condition, and behavioural plasticity, in particular learning. These concepts have revolutionized our view on fitness in novel ecological settings, and the way these processes may influence species ranges under HIREC. We have integrated them into a model of range expansion as a conceptual proof of principle highlighting the potentially substantial role of learning ability in range shifts under HIREC.

The angiosperm radiation played a dual role in the diversification of insects and insect pollinators

Interactions with angiosperms have been hypothesised to play a crucial role in driving diversification among insects, with a particular emphasis on pollinator insects. However, support for coevolutionary diversification in insect-plant interactions is weak. Macroevolutionary studies of insect and plant diversities support the hypothesis that angiosperms diversified after a peak in insect diversity in the Early Cretaceous. Here, we used the family-level fossil record of insects as a whole, and insect pollinator families in particular, to estimate diversification rates and the role of angiosperms on insect macroevolutionary history using a Bayesian process-based approach. We found that angiosperms played a dual role that changed through time, mitigating insect extinction in the Cretaceous and promoting insect origination in the Cenozoic, which is also recovered for insect pollinator families only. Although insects pollinated gymnosperms before the angiosperm radiation, a radiation of new pollinator lineages began as angiosperm lineages increased, particularly significant after 50 Ma. We also found that global temperature, increases in insect diversity, and spore plants were strongly correlated with origination and extinction rates, suggesting that multiple drivers influenced insect diversification and arguing for the investigation of different explanatory variables in further studies.

Evolution of realized niche breadth diversity driven by community dynamics

Why many herbivorous insects are host plant specialists, with non-negligible exceptions, is a conundrum of evolutionary biology, especially because the host plants are not necessarily optimal larval diets. Here, I present a novel model of host plant preference evolution of two insect species. Because habitat preference evolution is contingent upon demographic dynamics, I integrate the evolutionary framework with the modern coexistence theory. The results show that the two insect species can evolve into a habitat specialist and generalist, when they experience both negative and positive frequency-dependent community dynamics. This happens because the joint action of positive and negative frequency dependence creates multiple (up to nine) eco-evolutionary equilibria. Furthermore, initial condition dependence due to positive frequency dependence allows specialization to poor habitats. Thus, evolved habitat preferences do not necessarily correlate with the performances. The model provides explanations for counterintuitive empirical patterns and mechanistic interpretations for phenomenological models of niche breadth evolution.

Taxonomic Structure and Wing Pattern Evolution in the Parnassius mnemosyne Species Complex (Lepidoptera, Papilionidae)

In our study, using the analysis of DNA barcodes and morphology (wing color, male genitalia, and female sphragis shape), we show that the group of species close to P. mnemosyne comprises the western and eastern phylogenetic lineages. The eastern lineage includes P. stubbendorfii, P. glacialis, and P. hoenei. The western lineage includes three morphologically similar species: P. mnemosyne (Western Eurasia), P. turatii (southwestern Europe), and P. nubilosusstat. nov. (Turkmenistan and NE Iran), as well as the morphologically differentiated P. ariadne (Altai). The latter species differs from the rest of the group in the presence of red spots on the wings. Parnassius mnemosyne s.s. is represented by four differentiated mitochondrial clusters that show clear association with specific geographic regions. We propose to interpret them as subspecies: P. mnemosyne mnemosyne (Central and Eastern Europe, N Caucasus, N Turkey), P. mnemosyne adolphi (the Middle East), P. mnemosyne falsa (Tian Shan), and P. mnemosyne gigantea (Gissar-Alai in Central Asia). We demonstrate that in P. ariadne, the red spots on the wing evolved as a reversion to the ancestral wing pattern. This reversion is observed in Altai, where the distribution areas of the western lineage, represented by P. ariadne, and the eastern lineage, represented by P. stubbendorfii, overlap. These two species hybridize in Altai, and we hypothesize that the color change in P. ariadne is the result of reinforcement of prezygotic isolation in the contact zone. The lectotype of Parnassius mnemosyne var. nubilosus Christoph, 1873, is designated.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Parasite-host network analysis provides insights into the evolution of two mistletoe lineages (Loranthaceae and Santalaceae)

Mistletoes are ecologically important parasitic plants, with > 1600 species from five lineages worldwide. Mistletoe lineages exhibit distinct patterns of species diversification and host specificity, however, the mechanisms underlying these differences are poorly understood. In this study, we analysed a comprehensive parasite-host network, including 280 host species from 60 families and 22 mistletoe species from two lineages (Santalaceae and Loranthaceae) in Xishuangbanna, located in a biodiversity hotspot of tropical Asia. We identified the factors that predict the infection strength of mistletoes. We also detected host specificity and the phylogenetic signal of mistletoes and their hosts. We found that this interaction network could be largely explained by a model based on the relative abundance of species. Host infection was positively correlated with diameter at breast height and tree coverage, but negatively correlated with wood density. Overall, closely related mistletoe species tend to interact more often with similar hosts. However, the two lineages showed a significantly different network pattern. Rates of host generality were higher in Loranthaceae than in Santalaceae, although neither lineage showed phylogenetic signal for host generality. This study demonstrates that the neutral interaction hypothesis provides suitable predictions of the mistletoe-host interaction network, and mistletoe species show significant phylogenetic signals for their hosts. Our findings also indicate that high species diversification in Loranthaceae may be explained by high rates of host generality and the evolutionary history shared by Loranthaceae species with diverse host plants in the tropics.

Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

A global phylogeny of butterflies reveals their evolutionary history, ancestral hosts and biogeographic origins

Butterflies are a diverse and charismatic insect group that are thought to have evolved with plants and dispersed throughout the world in response to key geological events. However, these hypotheses have not been extensively tested because a comprehensive phylogenetic framework and datasets for butterfly larval hosts and global distributions are lacking. We sequenced 391 genes from nearly 2,300 butterfly species, sampled from 90 countries and 28 specimen collections, to reconstruct a new phylogenomic tree of butterflies representing 92% of all genera. Our phylogeny has strong support for nearly all nodes and demonstrates that at least 36 butterfly tribes require reclassification. Divergence time analyses imply an origin ~100 million years ago for butterflies and indicate that all but one family were present before the K/Pg extinction event. We aggregated larval host datasets and global distribution records and found that butterflies are likely to have first fed on Fabaceae and originated in what is now the Americas. Soon after the Cretaceous Thermal Maximum, butterflies crossed Beringia and diversified in the Palaeotropics. Our results also reveal that most butterfly species are specialists that feed on only one larval host plant family. However, generalist butterflies that consume two or more plant families usually feed on closely related plants.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Polyhydroxy Acids as Fabaceous Plant Components Induce Oviposition of the Common Grass Yellow Butterfly, Eurema Mandarina

The common grass yellow butterfly, Eurema mandarina is a Fabaceae-feeding species, the females of which readily oviposit on Albizia julibrissin and Lespedeza cuneata in mainland Japan. We previously demonstrated that the methanolic leaf extracts of these plants, and their highly polar aqueous fractions strongly elicit female oviposition. Furthermore, the three subfractions obtained by ion-exchange chromatographic separation of the aqueous fraction have been found to be less effective alone, but synergistically stimulate female oviposition when combined. This indicates that female butterflies respond to multiple compounds with different acidity. We have previously identified d-pinitol from the neutral/amphoteric subfractions and glycine betaine from the basic subfractions as oviposition stimulants of E. mandarina. The present study aimed to identify active compounds in the remaining acidic subfractions of A. julibrissin and L. cuneata leaf extracts. GC-MS analyses of trimethylsilyl-derivatized samples revealed the presence of six compounds in the acidic subfractions. In bioassays using these authentic chemicals, erythronic acid (EA) and threonic acid (TA) were moderately active in eliciting oviposition responses in E. mandarina, with their d-isomers showing slightly higher activity than their l-isomers. Female responsiveness differed between d-EA and l-TA, the major isomers of these compounds in plants, with the response to d-EA reaching a plateau at concentrations above 0.005% and that to l-TA peaking at a concentration of 0.01%. The natural concentrations of d-EA and l-TA in fresh A. julibrissin and L. cuneata leaves were sufficient to stimulate oviposition. Furthermore, mixing 0.001% d-EA or 0.001% l-TA, to which females are mostly unresponsive, with 0.1% d-pinitol resulted in a synergistic enhancement of the oviposition response. These findings demonstrate that E. mandarina females utilize both polyhydroxy acids, EA and TA, as chemical cues for oviposition.

Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry

The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of  480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.

Systematics of Opsiphanes Doubleday, [1849] (Lepidoptera: Nymphalidae, Satyrinae, Brassolini): an integrative approach

Opsiphanes Doubleday, [1849] (Lepidoptera: Nymphalidae: Satyrinae: Brassolini) is an exclusively Neotropical genus, occurring from Argentina to Mexico. Until the present study, Opsiphanes was considered to contain 14 species, 60 subspecies, and 38 synonyms. The considerable phenotypic variation of species and subspecies of the genus has affected the taxonomy of the group by causing the proliferation of several names that have been proposed to represent their diversity, taxa that have often not been adequately described and/or delimited. The present study analyzed information on the immature stages and morphology, with molecular data and distribution data, in order to provide revised taxonomic hypotheses for Opsiphanes species and subspecies. These analyses of approximately 5,500 specimens and all species known for the genus made it possible to define two groups: "cassiae" and "quiteria". The "quiteria" group was subdivided into seven subgroups: "boisduvallii", "camena", "zelotes", "sallei", "quiteria", "fabricii", and "invirae". The statuses of three species and two subespecies are reinstated: Opsiphanes badius Stichel, 1902 stat. rest., Opsiphanes quirinus Godman & Salvin, 1881 stat. rest., Opsiphanes merianae Stichel, 1902 stat. rest., Opsiphanes bogotanus castaneus Stichel, 1904 stat. rest. and Opsiphanes badius cauca Röber, 1906 stat. rest. Six subspecies are here treated as species: Opsiphanes mexicana Bristow, 1991 stat. nov., Opsiphanes zelus Stichel, 1908 stat. nov., Opsiphanes farrago Stichel, 1904 stat. nov., Opsiphanes barkeri Bristow, 1991 stat. nov., Opsiphanes caliensis Bristow, 1991 stat. nov., and Opsiphanes cuspidatus Stichel, 1904 stat. nov. One subjective synonym is treated as a valid subspecies: Opsiphanes invirae pernambucoensis Bristow, 1991 stat. rev. One species is treated as a subspecies: Opsiphanes cassiae tamarindi C. Felder & R. Felder, 1861 stat. nov. Eight new statuses are proposed: Opsiphanes cassiae incolumis Stichel, 1904 stat. nov., Opsiphanes cassiae tamarindi C. Felder & R. Felder, 1861 stat. nov., Opsiphanes badius angostura Bristow, 1979 stat. nov., Opsiphanes fabricii camposi Bristow, 1991 stat. nov., Opsiphanes fabricii numatius Fruhstorfer, 1912 stat. nov., Opsiphanes merianae notanda Stichel, 1904 stat. nov., Opsiphanes periphetes Fruhstorfer, 1912 stat. nov., and Opsiphanes cuspidatus relucens Fruhstorfer, 1907 stat. nov. Seven subjective synonyms are reinstated: Opsiphanes crameri C. Felder & R. Felder, 1862 syn. rest. of Opsiphanes cassiae cassiae (Linnaeus, 1758); Opsiphanes tamarindi latifascia Rothschild, 1916 syn. rest. of Opsiphanes cassiae incolumis Stichel, 1904 stat. nov.; Opsiphanes erebus Röber, 1927 syn. rest. of Opsiphanes quiteria quirinalis Staudinger, 1887; Opsiphanes cassina aequatorialis Stichel, 1902 syn. rest., Opsiphanes invirae pseudophilon Fruhstorfer, 1907 syn. rest., Opsiphanes invirae remoliatus Fruhstorfer, 1907 syn. rest., and Opsiphanes invirae agasthenes Fruhstorfer, 1907 syn. rest. of Opsiphanes invirae invirae (Hübner, [1808]). Twenty-five new synonyms are proposed: Pavonia Godart [1824] syn. nov. of Bia Hübner, [1819]; Opsiphanes bogotanus phrataphernes Fruhstorfer, 1912 syn. nov., and Opsiphanes bogotanus blandini Bristow, 1991 syn. nov. of Opsiphanes bogotanus bogotanus Distant, 1875; Opsiphanes cassiae alajuela Bristow, 1991 syn. nov. of Opsiphanes bogotanus castaneus Stichel, 1904 stat. rest.; Opsiphanes cassiae rubigatus Stichel, 1904 syn. nov., Opsiphanes cassiae strophios Fruhstorfer, 1907 syn. nov., and Opsiphanes tamarindi xiphos Fruhstorfer, 1907 syn. nov. of Opsiphanes cassiae cassiae (Linnaeus, 1758); Opsiphanes tamarindi corrosus Stichel, 1904 syn. nov., and Opsiphanes tamarindi kleisthenes Fruhstorfer, 1912 syn. nov. of Opsiphanes cassiae tamarindi C. Felder & R. Felder, 1861 stat. nov.; Opsiphanes mutatus parodizi Bristow, 1991 syn. nov. of Opsiphanes farrago Stichel, 1904 stat. nov.; Opsiphanes sallei kennerleyi Bristow, 1991 syn. nov. of Opsiphanes sallei colombiana Bristow, 1991; Opsiphanes quiteria talamancensis Bristow, 1991 syn. nov. of Opsiphanes quirinus Godman & Salvin, 1881 stat. rest.; Opsiphanes quiteria quaestor Stichel, 1902 syn. nov., Opsiphanes quiteria bolivianus Stichel, 1902 syn. nov., and Opsiphanes quiteria cardenasi Bristow, 1991 syn. nov. of Opsiphanes quiteria quiteria (Stoll, 1780); Opsiphanes quiteria phylas Fruhstorfer, 1912 syn. nov. of Opsiphanes quiteria quirinalis Staudinger, 1887; Opsiphanes cassina chiriquensis Stichel, 1902 syn. nov. of Opsiphanes fabricii fabricii (Boisduval, 1870); Opsiphanes cassina milesi Bristow, 1991 syn. nov. of Opsiphanes merianae notanda Stichel, 1904 stat. nov.; Opsiphanes cassina aucotti Bristow, 1991 syn. nov. of Opsiphanes periphetes Fruhstorfer, 1912 stat. nov.; Opsiphanes cassina C. Felder & R. Felder, 1862 syn. nov., Opsiphanes invirae intermedius Stichel, 1902 syn. nov., Opsiphanes invirae amplificatus Stichel, 1904 syn. nov., Opsiphanes sticheli Röber, 1906 syn. nov., Opsiphanes invirae roraimaensis Bristow, 1991 syn. nov., and Opsiphanes invirae sieberti Bristow, 1991 syn. nov. of Opsiphanes invirae invirae (Hübner, [1808]). To ensure unambiguous identification of names, nine neotypes were designated for: Opsiphanes bogotanus Distant, 1875, Opsiphanes aurivillii Röber, 1906, Papilio glycerie Fabricius, 1787, Opsiphanes zelotes zelus Stichel, 1908, Opsiphanes badius var. cauca Röber, 1906, Opsiphanes erebus Röber, 1927, Potamis invirae Hübner, [1808], Opsiphanes sticheli Röber, 1906, and Opsiphanes invirae ledon Fruhstorfer, 1912; and nine lectotypes for: Opsiphanes bogotanus phrataphernes Fruhstorfer, 1912, Opsiphanes tamarindi cherocles Fruhstorfer, 1912, Caligo tamarindi Boisduval, 1870, Opsiphanes sallei nicandrus Fruhstorfer, 1912, Opsiphanes quiteria augeias Fruhstorfer, 1912, Opsiphanes quirinus Godman & Salvin, 1881, Opsiphanes quiteria var. meridionalis Staudinger, 1887, Opsiphanes quiteria oresbios Fruhstorfer, 1912, and Opsiphanes quiteria phylas Fruhstorfer, 1912, . The present taxonomic scheme proposed for Opsiphanes includes 23 species, 23 subspecies, and 69 synonyms.

A phylogenetic study to assess the link between biome specialization and diversification in swallowtail butterflies

The resource-use hypothesis, proposed by E.S. Vrba, states that habitat fragmentation caused by climatic oscillations would affect particularly biome specialists (species inhabiting only one biome), which might show higher speciation and extinction rates than biome generalists. If true, lineages would accumulate biome-specialist species. This effect would be particularly exacerbated for biomes located at the periphery of the global climatic conditions, namely, biomes that have high/low precipitation and high/low temperature such as rainforest (warm-humid), desert (warm-dry), steppe (cold-dry) and tundra (cold-humid). Here, we test these hypotheses in swallowtail butterflies, a clade with more than 570 species, covering all the continents but Antarctica, and all climatic conditions. Swallowtail butterflies are among the most studied insects, and they are a model group for evolutionary biology and ecology studies. Continental macroecological rules are normally tested using vertebrates, this means that there are fewer examples exploring terrestrial invertebrate patterns at global scale. Here, we compiled a large Geographic Information System database on swallowtail butterflies' distribution maps and used the most complete time-calibrated phylogeny to quantify diversification rates (DRs). In this paper, we aim to answer the following questions: (1) Are there more biome-specialist swallowtail butterflies than biome generalists? (2) Is DR related to biome specialization? (3) If so, do swallowtail butterflies inhabiting extreme biomes show higher DRs? (4) What is the effect of species distribution area? Our results showed that swallowtail family presents a great number of biome specialists which showed substantially higher DRs compared to generalists. We also found that biome specialists are unevenly distributed across biomes. Overall, our results are consistent with the resource-use hypothesis, species climatic niche and biome fragmentation as key factors promoting isolation.

What Goes in Must Come Out? The Metabolic Profile of Plants and Caterpillars, Frass, And Adults of Asota (Erebidae: Aganainae) Feeding on Ficus (Moraceae) in New Guinea

Insect herbivores have evolved a broad spectrum of adaptations in response to the diversity of chemical defences employed by plants. Here we focus on two species of New Guinean Asota and determine how these specialist moths deal with the leaf alkaloids of their fig (Ficus) hosts. As each focal Asota species is restricted to one of three chemically distinct species of Ficus, we also test whether these specialized interactions lead to similar alkaloid profiles in both Asota species. We reared Asota caterpillars on their respective Ficus hosts in natural conditions and analyzed the alkaloid profiles of leaf, frass, caterpillar, and adult moth samples using UHPLC-MS/MS analyses. We identified 43 alkaloids in our samples. Leaf alkaloids showed various fates. Some were excreted in frass or found in caterpillars and adult moths. We also found two apparently novel indole alkaloids-likely synthesized de novo by the moths or their microbiota-in both caterpillar and adult tissue but not in leaves or frass. Overall, alkaloids unique or largely restricted to insect tissue were shared across moth species despite feeding on different hosts. This indicates that a limited number of plant compounds have a direct ecological function that is conserved among the studied species. Our results provide evidence for the importance of phytochemistry and metabolic strategies in the formation of plant-insect interactions and food webs in general. Furthermore, we provide a new potential example of insects acquiring chemicals for their benefit in an ecologically relevant insect genus.

Host association, environment, and geography underlie genomic differentiation in a major forest pest

Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high‐throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance‐based redundancy analyses (dbRDA), we show that differentiation of thousands of genome‐wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host‐associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.

Highly Stretchable Hydrogels as Wearable and Implantable Sensors for Recording Physiological and Brain Neural Signals

Recording electrophysiological information such as brain neural signals is of great importance in health monitoring and disease diagnosis. However, foreign body response and performance loss over time are major challenges stemming from the chemomechanical mismatch between sensors and tissues. Herein, microgels are utilized as large crosslinking centers in hydrogel networks to modulate the tradeoff between modulus and fatigue resistance/stretchability for producing hydrogels that closely match chemomechanical properties of neural tissues. The hydrogels exhibit notably different characteristics compared to nanoparticles reinforced hydrogels. The hydrogels exhibit relatively low modulus, good stretchability, and outstanding fatigue resistance. It is demonstrated that the hydrogels are well suited for fashioning into wearable and implantable sensors that can obtain physiological pressure signals, record the local field potentials in rat brains, and transmit signals through the injured peripheral nerves of rats. The hydrogels exhibit good chemomechanical match to tissues, negligible foreign body response, and minimal signal attenuation over an extended time, and as such is successfully demonstrated for use as long-term implantable sensory devices. This work facilitates a deeper understanding of biohybrid interfaces, while also advancing the technical design concepts for implantable neural probes that efficiently obtain physiological information.

Phylogeny and Biogeographic History of Parnassius Butterflies (Papilionidae: Parnassiinae) Reveal Their Origin and Deep Diversification in West China

Simple Summary

Butterflies of the genus Parnassius are distributed in the mountains across the Northern Hemisphere. Studies have shown that this genus originated in the regions of West China–Central Asia. To further explore the spatiotemporal pattern and driving mechanism of Parnassius diversification, we reconstructed the phylogeny and biogeographic history of Parnassius based on 45 species. Ancestral area reconstruction obtained by using the statistical dispersal–extinction cladogenesis model indicated that Parnassius originated in West China (Qinghai–Tibet Plateau and Xinjiang) during the Middle Miocene. Paleoenvironment changes and host plants were probably influenced by the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America. Furthermore, ancient gene introgression might have contributed to the spread of Parnassius butterflies from the high mountains of the Qinghai–Tibet Plateau to the low-altitude areas of Central East China. This study will provide an understanding of the phylogeny and biogeographic history of the genus Parnassius.

Abstract

We studied 239 imagoes of 12 Parnassius species collected from the mountains of the Qinghai–Tibet Plateau (QTP) and its neighbouring areas in China. We selected three mitochondrial gene (COI, ND1, and ND5) sequences, along with the homologous gene sequences of other Parnassius species from GenBank, to reconstruct the phylogenetic tree and biogeographic history of this genus. Our results show that Parnassius comprises eight monophyletic subgenera, with subgenus Parnassius at the basal position; the genus crown group originated during the Middle Miocene (ca. 16.99 Ma), and species diversification continued during sustained cooling phases after the Middle Miocene Climate Optimum (MMCO) when the QTP and its neighbouring regions experienced rapid uplift and extensive orogeny. A phylogenetic network analysis based on transcriptomes from GenBank suggests that ancient gene introgression might have contributed to the spread of the Parnassius genus to different altitudes. Ancestral area reconstruction indicates that Parnassius most likely originated in West China (QTP and Xinjiang) and then spread to America in two dispersal events as subgenera Driopa and Parnassius, along with their host plants Papaveraceae and Crassulaceae, respectively. Our study suggests that extensive mountain-building processes led to habitat fragmentation in the QTP, leading to the early diversification of Parnassius, and climate cooling after MMCO was the driving mechanism for the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Expanding the Menu: Are Polyphagy and Gene Family Expansions Linked across Lepidoptera?

Evolutionary expansions and contractions of gene families are often correlated with key innovations and/or ecological characteristics. In butterflies and moths (Lepidoptera), expansions of gene families involved in detoxification of plant specialized metabolites are hypothesized to facilitate a polyphagous feeding style. However, analyses supporting this hypothesis are mostly based on a limited number of lepidopteran species. We applied a phylogenomics approach, using 37 lepidopteran genomes, to analyze if gene family evolution (gene gain and loss) is associated with the evolution of polyphagy. Specifically, we compared gene counts and evolutionary gene gain and loss rates of gene families involved in adaptations with plant feeding. We correlated gene evolution to host plant family range (phylogenetic diversity) and specialized metabolite content of plant families (functional metabolite diversity). We found a higher rate for gene loss than gene gain in Lepidoptera, a potential consequence of genomic rearrangements and deletions after (potentially small-scale) duplication events. Gene family expansions and contractions varied across lepidopteran families, and were associated to host plant use and specialization levels. Within the family Noctuidae, a higher expansion rate for gene families involved in detoxification can be related to the large number of polyphagous species. However, gene family expansions are observed in both polyphagous and monophagous lepidopteran species and thus seem to be species-specific in the taxa sampled. Nevertheless, a significant positive correlation of gene counts of the carboxyl- and choline esterase and glutathione-S-transferase detoxification gene families with the level of polyphagy was identified across Lepidoptera.

Lifestyle Evolution Analysis by Binary-State Speciation and Extinction (BiSSE) Model

Phylogenetic comparative methods (PCMs) combine statistics and evolutionary models to infer the dynamics of trait evolution and diversification that underlie the observed phylogeny. While PCMs have been used to study macro-evolutionary processes and evolutionary transitions of macroorganisms, their application to microbes is still limited. With the abundance of publicly available genomic and trait character data for diverse microbes nowadays, applications of PCMs on these data can provide insights into the fundamental principles that govern microbial evolution. Here, we introduce the Binary-State Speciation and Extinction (BiSSE) model, which is a relatively simple yet powerful approach for analyzing trait evolution. We begin by explaining the theoretical background and intuition behind the BiSSE model. Then, R commands for running the BiSSE model are presented. Finally, we introduce a case study that successfully applied the BiSSE model to investigate generalist and specialist microbial lifestyle evolution.

Food Plant Shifts Drive the Diversification of Sack-Bearer Moths

AbstractLepidoptera are a highly diverse group of herbivorous insects; however, some superfamilies have relatively few species. Two alternative hypotheses for drivers of Lepidoptera diversity are shifts in food plant use or shifts from concealed to external feeding as larvae. Many studies address the former hypothesis but with bias toward externally feeding taxa. One of the most striking examples of species disparity between sister lineages in Lepidoptera is between the concealed-feeding sack-bearer moths (Mimallonoidea), which contain about 300 species, and externally feeding Macroheterocera, which have over 74,000 species. We provide the first dated tree of Mimallonidae to understand the diversification dynamics of these moths in order to fill a knowledge gap pertaining to drivers of diversity within an important concealed-feeding clade. We find that Mimallonidae is an ancient Lepidoptera lineage that originated in the Cretaceous ∼105 million years ago and has had a close association with the plant order Myrtales for the past 40 million years. Diversification dynamics are tightly linked with food plant usage in this group. Reliance on Myrtales may have influenced diversification of Mimallonidae because clades that shifted away from the ancestral condition of feeding on Myrtales have the highest speciation rates in the family.

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

Neofunctionalization of an ancient domain allows parasites to avoid intraspecific competition by manipulating host behaviour

Intraspecific competition is a major force in mediating population dynamics, fuelling adaptation, and potentially leading to evolutionary diversification. Among the evolutionary arms races between parasites, one of the most fundamental and intriguing behavioural adaptations and counter-adaptations are superparasitism and superparasitism avoidance. However, the underlying mechanisms and ecological contexts of these phenomena remain underexplored. Here, we apply the Drosophila parasite Leptopilina boulardi as a study system and find that this solitary endoparasitic wasp provokes a host escape response for superparasitism avoidance. We combine multi-omics and in vivo functional studies to characterize a small set of RhoGAP domain-containing genes that mediate the parasite's manipulation of host escape behaviour by inducing reactive oxygen species in the host central nervous system. We further uncover an evolutionary scenario in which neofunctionalization and specialization gave rise to the novel role of RhoGAP domain in avoiding superparasitism, with an ancestral origin prior to the divergence between Leptopilina specialist and generalist species. Our study suggests that superparasitism avoidance is adaptive for a parasite and adds to our understanding of how the molecular manipulation of host behaviour has evolved in this system.

The evolution of coevolution in the study of species interactions

The study of reciprocal adaptation in interacting species has been an active and inspiring area of evolutionary research for nearly 60 years. Perhaps owing to its great natural history and potential consequences spanning population divergence to species diversification, coevolution continues to capture the imagination of biologists. Here we trace developments following Ehrlich and Raven's classic paper, with a particular focus on the modern influence of two studies by Dr. May Berenbaum in the 1980s. This series of classic work presented a compelling example exhibiting the macroevolutionary patterns predicted by Ehrlich and Raven and also formalized a microevolutionary approach to measuring selection, functional traits, and understanding reciprocal adaptation between plants and their herbivores. Following this breakthrough was a wave of research focusing on diversifying macroevolutionary patterns, mechanistic chemical ecology, and natural selection on populations within and across community types. Accordingly, we breakdown coevolutionary theory into specific hypotheses at different scales: reciprocal adaptation between populations within a community, differential coevolution among communities, lineage divergence, and phylogenetic patterns. We highlight progress as well as persistent gaps, especially the link between reciprocal adaptation and diversification.

Discordant Phylogenomic Placement of Hydnoraceae and Lactoridaceae Within Piperales Using Data From All Three Genomes

Phylogenetic relationships within the magnoliid order Piperales have been studied extensively, yet the relationships of the monotypic family Lactoridaceae and the holoparasitic Hydnoraceae to the remainder of the order remain a matter of debate. Since the first confident molecular phylogenetic placement of Hydnoraceae among Piperales, different studies have recovered various contradictory topologies. Most phylogenetic hypotheses were inferred using only a few loci and have had incomplete taxon sampling at the genus level. Based on these results and an online survey of taxonomic opinion, the Angiosperm Phylogeny Group lumped both Hydnoraceae and Lactoridaceae in Aristolochiaceae; however, the latter family continues to have unclear relationships to the aforementioned taxa. Here we present extensive phylogenomic tree reconstructions based on up to 137 loci from all three subcellular genomes for all genera of Piperales. We infer relationships based on a variety of phylogenetic methods, explore instances of phylogenomic discordance between the subcellular genomes, and test alternative topologies. Consistent with these phylogenomic results and a consideration of the principles of phylogenetic classification, we propose to exclude Hydnoraceae and Lactoridaceae from the broad circumscription of Aristolochiaceae, and instead favor recognition of four monophyletic and morphologically well circumscribed families in the perianth-bearing Piperales: Aristolochiaceae, Asaraceae, Hydnoraceae, and Lactoridaceae, with a total of six families in the order.

Two complete mitochondrial genomes of Papilio butterflies obtained from historical specimens (Lepidoptera: Papilionidae)

Museum specimens are collected for education, exhibition, and various multiple scientific purposes. However, millions of specimens remain in their collection boxes for years without being analyzed. Historical specimens have been known to contain low-quality DNA; hence, it is difficult to utilize their sequence information in phylogenetic studies. However, recent advances in high-throughput sequencing (HTS) make these collections amenable to phylogenomic studies. In this study, two historical specimens (Papilio xuthus Linnaeus, 1767, and Papilio thoas Linnaeus, 1771) were sampled and DNA extracted for HTS via the Miseq platform. Two complete mitogenomes were assembled, even though the DNA quality of those specimens was highly fragmented, below 250 bp in length. The 37 genes of 60 mitogenomes were aligned and used for inferring the phylogenetic relationships of Papilioninae. These two newly sequenced mitogenomes are correctly grouped in the genus Papilio, and this result indicates that historical specimens show great potential for phylogenetic studies with HTS technology.

Sharing and double-dating in the lichen world

Historic and modern efforts to understand lichen diversity and evolution have overwhelmingly concentrated on that of the fungal partner, which represents one of the most taxonomically diverse nutritional modes among the Fungi. But what about the algal and cyanobacterial symbionts? An explosion of studies on these cryptic symbionts over the past 20+ years has facilitated a richer understanding of their diversity, patterns of association, and the symbiosis itself. In a From the Cover article in this issue of Molecular Ecology, Dal Forno et al. (2021) provide new insight into one of the most fascinating lichen symbioses. By sequencing cyanobacterial symbionts from over 650 specimens, they reveal the presence of overlooked cyanobacterial diversity, evidence for symbiont sharing among distantly related fungi, and utilize a comparative dating framework to demonstrate temporal discordance among interacting fungal and cyanobacterial lineages.