Evolutionary Relationship Between Platycerus Stag Beetles and Their Mycangium-Associated Yeast Symbionts

Investigation of Fungal Community Structure in the Gut of the Stag Beetle Dorcus hopei (Coleoptera; Lucanidae): Comparisons Among Developmental Stages

Stag beetles, recognized as common saproxylic insects, are valued for their vibrant coloration and distinctive morphology. These beetles play a crucial ecological role in decomposition and nutrient cycling, serving as a vital functional component in ecosystem functioning. Although previous studies have confirmed that stag beetles are predominantly fungivores, the fluctuations in their intestinal fungal communities at different developmental stages remain poorly understood. In the current study, high-throughput sequencing was employed to investigate the dynamic changes within intestinal fungal communities at various developmental stages in the stag beetle Dorcus hopei. Results showed that microbial diversity was higher during the larval stage than during the pupal and adult stages. Furthermore, significant differences were identified in the composition of the intestinal fungal communities across the larval, pupal, and adult stages, suggesting that developmental transitions may be crucial factors contributing to variations in fungal community composition and diversity. Dominant genera included Candida, Scheffersomyces, Phaeoacremonium, and Trichosporon. Functional predictions indicated a greater diversity and relative abundance of endosymbiotic fungi in the larval gut, suggesting a potential dependency of larvae on beneficial gut fungi for nutrient acquisition. Additionally, the application of abundance-based β-null deviation and niche width analyses revealed that the adult gut exerted a stronger selection pressure on its fungal community, favoring certain taxa. This selection process culminates in a more robust co-occurrence network of fungal communities within the adult gut, thereby enhancing their adaptability to environmental fluctuations. This study advances our understanding of the intestinal fungal community structure in stag beetles, providing a crucial theoretical foundation for the development of saproxylic beetle resources, biomass energy utilization, plastic degradation strategies, and beetle conservation efforts.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae)

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Presence of non-symbiotic yeasts in a symbiont-transferring organ of a stag beetle that lacks yeast symbionts found in other stag beetles

Dispersal from wood to wood is essential for wood-inhabiting fungi and wood-inhabiting insects play an important role in the dispersal success of such fungi. However, it is poorly understood whether wood-inhabiting insects which change the habitats from wood to non-wood environments can contribute to the fungal dispersal. Larvae of most stag beetles (Coleoptera: Lucanidae) are wood feeders, while adults are sap feeders. Female adults of lulcanids possess specialized organs (mycetangia) for transportation of fungal symbionts and harbor specific yeasts (e.g., Scheffersomyces spp.) within. Here, we report that the lucanid Aegus subnitidus harbors non-specific yeasts facultatively in mycetangia. We conducted yeast isolation from mycetangia and hindguts of female adults, in a larval gallery in wood-associated materials, and in female-visiting fermented tree sap using culture-dependent methods. Less than half of the females carried a total of 20 yeast species with small amounts using mycetangia and a female harbored up to five species, suggesting the absence of an intimate association with specific yeasts that are found in other lucanids. Yeast species compositions markedly differed between the larval gallery and sap. Most yeasts from the larval galley exhibited xylose-assimilation abilities, while few yeasts from sap did. Mycetangial yeasts comprised a combination from both sources. In hindguts, most yeasts were found in sap (> 70%) with no yeasts in the larval gallery. Sap-associated yeasts in each female mycetangium were also obtained from the female-visiting sap patch, while mycetangial, larval gallery-associated yeasts were absent from the patch, suggesting the survival success of larval gallery-associated yeasts in mycetangia through sap patches. Therefore, wood-inhabiting insects that possess mycetangia can potentially act as vectors of non-symbiotic wood-inhabiting yeasts dispersing from wood to wood via other environments.

Morphological adaptation for ectosymbiont maintenance and transmission during metamorphosis in Lagria beetles

The diversity and success of holometabolous insects is partly driven by metamorphosis, which allows for the exploitation of different niches and decouples growth and tissue differentiation from reproduction. Despite its benefits, metamorphosis comes with the cost of temporal vulnerability during pupation and challenges associated with tissue reorganizations. These rearrangements can also affect the presence, abundance, and localization of beneficial microbes in the host. However, how symbionts are maintained or translocated during metamorphosis and which adaptations are necessary from each partner during this process remains unknown for the vast majority of symbiotic systems. Here, we show that Lagria beetles circumvent the constraints of metamorphosis by maintaining defensive symbionts on the surface in specialized cuticular structures. The symbionts are present in both sexes throughout larval development and during the pupal phase, in line with a protective role during the beetle’s immature stages. By comparing symbiont titer and morphology of the cuticular structures between sexes using qPCR, fluorescence in situ hybridization, and micro-computed tomography, we found that the organs likely play an important role as a symbiont reservoir for transmission to female adults, since symbiont titers and structures are reduced in male pupae. Using symbiont-sized fluorescent beads, we demonstrate transfer from the region of the dorsal symbiont-housing organs to the opening of the reproductive tract of adult females, suggesting that symbiont relocation on the outer surface is possible, even without specialized symbiont adaptations or motility. Our results illustrate a strategy for holometabolous insects to cope with the challenge of symbiont maintenance during metamorphosis via an external route, circumventing problems associated with internal tissue reorganization. Thereby, Lagria beetles keep a tight relationship with their beneficial partners during growth and metamorphosis.

Local ecological divergence of two closely related stag beetles based on genetic, morphological, and environmental analyses

The process of phenotypic adaptation to the environments is widely recognized. However, comprehensive studies integrating phylogenetic, phenotypic, and ecological approaches to assess this process are scarce. Our study aims to assess whether local adaptation may explain intraspecific differentiation by quantifying multidimensional differences among populations in closely related lucanid species, Platycerus delicatulus and Platycerus kawadai, which are endemic saproxylic beetles in Japan. First, we determined intraspecific analysis units based on nuclear and mitochondrial gene analyses of Platycerus delicatulus and Platycerus kawadai under sympatric and allopatric conditions. Then, we compared differences in morphology and environmental niche between populations (analysis units) within species. We examined the relationship between morphology and environmental niche via geographic distance. P. kawadai was subdivided into the “No introgression” and “Introgression” populations based on mitochondrial COI gene – nuclear ITS region discordance. P. delicatulus was subdivided into “Allopatric” and “Sympatric” populations. Body length differed significantly among the populations of each species. For P. delicatulus, character displacement was suggested. For P. kawadai, the morphological difference was likely caused by geographic distance or genetic divergence rather than environmental differences. The finding showed that the observed mitochondrial–nuclear discordance is likely due to historical mitochondrial introgression following a range of expansion. Our results show that morphological variation among populations of P. delicatulus and P. kawadai reflects an ecological adaptation process based on interspecific interactions, geographic distance, or genetic divergence. Our results will deepen understanding of ecological specialization processes across the distribution and adaptation of species in natural systems.

Diverse Heat Tolerance of the Yeast Symbionts of Platycerus Stag Beetles in Japan

The genus Platycerus (Coleoptera: Lucanidae) is a small stag beetle group, which is adapted to cool-temperate deciduous broad-leaved forests in East Asia. Ten Platycerus species in Japan form a monophyletic clade endemic to Japan and inhabit species-specific climatic zones. They are reported to have co-evolutionary associations with their yeast symbionts of the genus Sheffersomyces based on host cytochrome oxidase subunit I (COI) and yeast intergenic spacer (IGS) phylogenies. Here we examined the heat tolerances of the yeast colonies isolated from the mycangia of 37 females belonging ten Japanese Platycerus species. The upper limits of growth and survival temperatures of each colony were decided by cultivating it at ten temperature levels between 17.5 and 40°C. Although both temperatures varied during 25.0–31.25°C, the maximum survival temperatures (MSTs) were a little higher than the maximum growth temperatures (MGTs) in 16 colonies. Pearson’s correlations between these temperatures and environmental factors (elevation and 19 bioclimatic variables from Worldclim database) of host beetle collection sites were calculated. These temperatures were significantly correlated with elevation negatively, the maximum temperature of the warmest month (Bio5) positively, and some precipitative variables, especially in the warm season (Bio12, 13, 16, 18) negatively. Sympatric Platycerus kawadai and Platycerus albisomni share the same lineage of yeast symbionts that exhibit the same heat tolerance, but the elevational lower range limit of P. kawadai is higher than that of P. albisomni. Based on the field survey in their sympatric site, the maximum temperature of host wood of P. kawadai larvae is higher about 2–3°C than that of P. albisomni larvae in the summer, which may restrict the elevational range of P. kawadai to higher area. In conclusion, it is suggested that the heat tolerance of yeast symbionts restricts the habitat range of their host Platycerus species or/and that the environmental condition that host Platycerus species prefers affect the heat tolerance of its yeast symbionts.

Lateral Transmission of Yeast Symbionts Among Lucanid Beetle Taxa

To deepen understanding the evolutionary process of lucanid–yeast association, the lateral transmission process of yeast symbionts among stag beetle genera Platycerus and Prismognathus around the border between Japan and South Korea was estimated based on molecular analyses and species distribution modelings. Phylogenetic analyses were based on yeast ITS and IGS sequences and beetle COI sequences using Prismognathus dauricus from the Tsushima Islands and Pr. angularis from Kyushu, Japan, as well as other sequence data from our previous studies. The range overlap based on the species distribution model (SDM) and differentiation in ecological space were analyzed. Based on the IGS sequences, Clade II yeast symbionts were shared by Platycerus hongwonpyoi and Pr. dauricus in South Korea and the Tsushima Islands, and Platycerus viridicuprus in Japan. Clade III yeasts were shared by Pr. dauricus from the Tsushima Islands and Pr. angularis in Japan. During the Last Interglacial period when the land bridge between Japan and the Korean Peninsula existed, range overlap was predicted to occur between Pl. viridicuprus and Pr. dauricus in Kyushu and between Pr. dauricus and Pr. angularis in Kyushu and the Tsushima Islands. The ecological space of Pl. hongwonpyoi was differentiated from that of Pl. viridicuprus and Pr. angularis. We demonstrated the paleogeographical lateral transmission process of Scheffersomyces yeast symbionts among lucanid genera and species: putative transmission of yeasts from Pr. dauricus to Pl. viridicuprus in Kyushu and from Pr. angularis to Pr. dauricus in Kyushu or the Tsushima Islands. We also found that the yeast symbionts are likely being replaced in Pr. dauricus on the Tsushima Islands. We present novel estimates of the lateral transmission process of microbial symbionts based on phylogenetic, SDM and environmental analyses among lucanid beetles.

Integrating intraspecific variation in species distribution models by quantifying niche differentiation

Intraspecific variation provides insight into heterogeneous demography and adaptive history across distribution ranges of organisms. Most species distribution models assume that species respond to the environment as a single undifferentiated entity across their entire distribution. However, populations occupying different niches might differ in their ability to cope with climate change. Here, we assessed potential climatic niche differentiation at the species and subspecies levels and identified the palaeoclimatic range of three cold-adapted, low-dispersal beetle species: Carabus irregularis, Platycerus albisomni and Platycerus takakuwai. Our results showed that: (1) MaxEnt models incorporating information derived from intraspecific variation outperformed the species-level models; (2) tests of niche similarity revealed niche conservatism in most subspecies, except for two subspecies of C. irregularis, C. i. irregularis and C. i. bucephalus; and (3) historical predictions suggested substantial shifts within species ranges, with multiple glacial refugia during the Last Glacial Maximum. In conclusion, we recommend closer examination of intraspecific variation when predicting species distributions, in order to obtain more accurate generalizations regarding range shifts under climate change.