Phylogenetic composition of host plant communities drives plant-herbivore food web structure.

Insects tend to feed on related hosts. The phylogenetic composition of host plant communities thus plays a prominent role in determining insect specialization, food web structure, and diversity. Previous studies showed a high preference of insect herbivores for congeneric and confamilial hosts suggesting that some levels of host plant relationships may play more prominent role that others. We aim to quantify the effects of host phylogeny on the structure of quantitative plant-herbivore food webs. Further, we identify specific patterns in three insect guilds with different life histories and discuss the role of host plant phylogeny in maintaining their diversity. We studied herbivore assemblages in three temperate forests in Japan and the Czech Republic. Sampling from a canopy crane, a cherry picker and felled trees allowed a complete census of plant-herbivore interactions within three 0·1 ha plots for leaf chewing larvae, miners, and gallers. We analyzed the effects of host phylogeny by comparing the observed food webs with randomized models of host selection. Larval leaf chewers exhibited high generality at all three sites, whereas gallers and miners were almost exclusively monophagous. Leaf chewer generality dropped rapidly when older host lineages (5-80 myr) were collated into a single lineage but only decreased slightly when the most closely related congeneric hosts were collated. This shows that leaf chewer generality has been maintained by feeding on confamilial hosts while only a few herbivores were shared between more distant plant lineages and, surprisingly, between some congeneric hosts. In contrast, miner and galler generality was maintained mainly by the terminal nodes of the host phylogeny and dropped immediately after collating congeneric hosts into single lineages. We show that not all levels of host plant phylogeny are equal in their effect on structuring plant-herbivore food webs. In the case of generalist guilds, it is the phylogeny of deeper plant lineages that drives the food web structure whereas the terminal relationships play minor roles. In contrast, the specialization and abundance of monophagous guilds are affected mainly by the terminal parts of the plant phylogeny and do not generally reflect deeper host phylogeny.

Mycobiomes of two distinct clades of ambrosia gall midges (Diptera: Cecidomyiidae) are species-specific in larvae but similar in nutritive mycelia.

IMPORTANCE: Ambrosia gall midges are endophagous insect herbivores whose larvae live enclosed within a single gall for their entire development period. They may exhibit phytomycetophagy, a remarkable feeding mode that involves the consumption of plant biomass and mycelia of their cultivated gall symbionts. Thus, AGMs are ideal model organisms for studying the role of microorganisms in the evolution of host specificity in insects. However, compared to other fungus-farming insects, insect-fungus mutualism in AGMs has been neglected. Our study is the first to use DNA metabarcoding to characterize the complete mycobiome of the entire system of the gall-forming insects as we profiled gall surfaces, nutritive mycelia, and larvae. Interestingly, larval mycobiomes were significantly different from their nutritive mycelia, although Botryosphaeria dothidea dominated the nutritive mycelia, regardless of the evolutionary separation of the tribes studied. Therefore, we confirmed a long-time hypothesized paradigm for the important evolutionary association of this fungus with AGMs.

Concentration-dependent effect of plant secondary metabolites on bacterial and fungal microbiomes in caterpillar guts.

IMPORTANCE: The caterpillar gut is an excellent model system for studying host-microbiome interactions, as it represents an extreme environment for microbial life that usually has low diversity and considerable variability in community composition. Our study design combines feeding caterpillars on a natural and artificial diet with controlled levels of plant secondary metabolites and uses metabarcoding and quantitative PCR to simultaneously profile bacterial and fungal assemblages, which has never been performed. Moreover, we focus on multiple caterpillar species and consider diet breadth. Contrary to many previous studies, our study suggested the functional importance of certain microbial taxa, especially bacteria, and confirmed the previously proposed lower importance of fungi for caterpillar holobiont. Our study revealed the lack of differences between monophagous and polyphagous species in the responses of microbial assemblages to plant secondary metabolites, suggesting the limited role of the microbiome in the plasticity of the herbivore diet.

Seasonal Shifts in Bacterial and Fungal Microbiomes of Leaves and Associated Leaf-Mining Larvae Reveal Persistence of Core Taxa Regardless of Diet.

Microorganisms are key mediators of interactions between insect herbivores and their host plants. Despite a substantial interest in studying various aspects of these interactions, temporal variations in microbiomes of woody plants and their consumers remain understudied. In this study, we investigated shifts in the microbiomes of leaf-mining larvae (Insecta: Lepidoptera) and their host trees over one growing season in a deciduous temperate forest. We used 16S and ITS2 rRNA gene metabarcoding to profile the bacterial and fungal microbiomes of leaves and larvae. We found pronounced shifts in the leaf and larval microbiota composition and richness as the season progressed, and bacteria and fungi showed consistent patterns. The quantitative similarity between leaf and larval microbiota was very low for bacteria (~9%) and decreased throughout the season, whereas fungal similarity increased and was relatively high (~27%). In both leaves and larvae, seasonality, along with host taxonomy, was the most important factor shaping microbial communities. We identified frequently occurring microbial taxa with significant seasonal trends, including those more prevalent in larvae (Streptococcus, Candida sake, Debaryomyces prosopidis, and Neoascochyta europaea), more prevalent in leaves (Erwinia, Seimatosporium quercinum, Curvibasidium cygneicollum, Curtobacterium, Ceramothyrium carniolicum, and Mycosphaerelloides madeirae), and frequent in both leaves and larvae (bacterial strain P3OB-42, Methylobacterium/Methylorubrum, Bacillus, Acinetobacter, Cutibacterium, and Botrytis cinerea). Our results highlight the importance of considering seasonality when studying the interactions between plants, herbivorous insects, and their respective microbiomes, and illustrate a range of microbial taxa persistent in larvae, regardless of their occurrence in the diet. IMPORTANCE Leaf miners are endophagous insect herbivores that feed on plant tissues and develop and live enclosed between the epidermis layers of a single leaf for their entire life cycle. Such close association is a precondition for the evolution of more intimate host-microbe relationships than those found in free-feeding herbivores. Simultaneous comparison of bacterial and fungal microbiomes of leaves and their tightly linked consumers over time represents an interesting study system that could fundamentally contribute to the ongoing debate on the microbial residence of insect gut. Furthermore, leaf miners are ideal model organisms for interpreting the ecological and evolutionary roles of microbiota in host plant specialization. In this study, the larvae harbored specific microbial communities consisting of core microbiome members. Observed patterns suggest that microbes, especially bacteria, may play more important roles in the caterpillar holobiont than generally presumed.

Individual variability in habitat selection by aquatic insects is driven by taxonomy rather than specialisation.

Habitat selection, the choice of a habitat based on its perceived quality, is a key mechanism structuring freshwater communities. To date, individual variability in habitat selection has been neglected, and specialisation has never been considered in this type of studies. We examined the individual differences in the habitat selection of backswimmers (Notonectidae) and diving beetles (Dytiscidae). From each family, we selected one habitat generalist able to coexist with fish (Notonecta glauca, Dytiscus marginalis), and one species specialised to fishless habitats (Notonecta obliqua, Acilius sulcatus). We performed a mesocosm experiment quantifying the consistency in individuals' decisions in response to fish and vegetation structure, in relation to sex and specialisation. Neither the overall pattern of preferences nor consistency in individuals' decisions differed between specialists and generalists or between the sexes, but both were consistent within families. At the population level, backswimmers preferred fishless pools with submersed and floating macrophytes, while diving beetles showed no clear preferences. Individual decisions of backswimmers were consistent and likely driven by conspecific/heterospecific attraction. In diving beetles, individual decisions were primarily density-dependent. Our results reinforce the significance of habitat selectivity for aquatic community assembly, while suggesting a range of mechanisms driving variability in individual behaviour.

Fungi are more transient than bacteria in caterpillar gut microbiomes.

Despite an increasing number of studies on caterpillar (Insecta: Lepidoptera) gut microbiota, bacteria have been emphasized more than fungi. Therefore, we lack data on whether fungal microbiota is resident or transient and shaped by factors similar to those of bacteria. We sampled nine polyphagous caterpillar species from several tree species at multiple sites to determine the factors shaping leaf and gut bacterial and fungal microbiota as well as the extent to which caterpillars acquire microbiota from their diet. We performed 16S and ITS2 DNA metabarcoding of the leaves and guts to determine the composition and richness of the respective microbiota. While spatial variables shaped the bacterial and fungal microbiota of the leaves, they only affected fungi in the guts, whereas the bacteria were shaped primarily by caterpillar species, with some species harboring more specific bacterial consortia. Leaf and gut microbiota significantly differed; in bacteria, this difference was more pronounced. The quantitative similarity between leaves and guts significantly differed among caterpillar species in bacteria but not fungi, suggesting that some species have more transient bacterial microbiota. Our results suggest the complexity of the factors shaping the gut microbiota, while highlighting interspecific differences in microbiota residency within the same insect functional group.

Specialization directs habitat selection responses to a top predator in semiaquatic but not aquatic taxa.

Habitat selectivity has become an increasingly acknowledged mechanism shaping the structure of freshwater communities; however, most studies have focused on the effect of predators and competitors, neglecting habitat complexity and specialization. In this study, we examined the habitat selection of semiaquatic (amphibians: Bufonidae; odonates: Libellulidae) and aquatic organisms (true bugs: Notonectidae; diving beetles: Dytiscidae). From each family, we selected one habitat generalist species able to coexist with fish (Bufo bufo, Sympetrum sanguineum, Notonecta glauca, Dytiscus marginalis) and one species specialized in fishless habitats (Bufotes viridis, Sympetrum danae, Notonecta obliqua, Acilius sulcatus). In a mesocosm experiment, we quantified habitat selection decisions in response to the non-consumptive presence of fish (Carassius auratus) and vegetation structure mimicking different successional stages of aquatic habitats (no macrophytes; submerged and floating macrophytes; submerged, floating, and littoral-emergent macrophytes). No congruence between habitat specialists and generalists was observed, but a similar response to fish and vegetation structure defined both semiaquatic and aquatic organisms. While semiaquatic generalists did not distinguish between fish and fishless pools, specialists avoided fish-occupied pools and had a preferred vegetation structure. In aquatic taxa, predator presence affected habitat selection only in combination with vegetation structure, and all species preferred fishless pools with floating and submerged macrophytes. Fish presence triggered avoidance only in the generalist bug N. glauca. Our results highlight the significance of habitat selectivity for structuring freshwater ecosystems and illustrate how habitat selection responses to a top predator are dictated by specialization and life history.

Wing wettability gradient in a damselfly Lestes sponsa (Odonata: Lestidae) reflects the submergence behaviour during underwater oviposition.

The phenomenon of hydrophobicity of insect cuticles has received great attention from technical fields due to its wide applicability to industry or medicine. However, in an ecological/evolutionary context such studies remain scarce. We measured spatial differences in wing wettability in Lestes sponsa (Odonata: Lestidae), a damselfly species that can submerge during oviposition, and discussed the possible functional significance. Using dynamic contact angle (CA) measurements together with scanning electron microscopy (SEM), we investigated differences in wettability among distal, middle and proximal wing regions, and in surface nanostructures potentially responsible for observed differences. As we moved from distal towards more proximal parts, mean values of advancing and receding CAs gradually increased from 104° to 149°, and from 67° to 123°, respectively, indicating that wing tips were significantly less hydrophobic than more proximal parts. Moreover, values of CA hysteresis for the respective wing parts decreased from 38° to 26°, suggesting greater instability of the structure of the wing tips. Accordingly, compared with more proximal parts, SEM revealed higher damage of the wax nanostructures at the distal region. The observed wettability gradient is well explained by the submergence behaviour of L. sponsa during underwater oviposition. Our study thus proposed the existence of species-dependent hydrophobicity gradient on odonate wings caused by different ovipositional strategies.

Effect of Bacterial and Fungal Microbiota Removal on the Survival and Development of Bryophagous Beetles.

Insect microbiota may play a wide range of roles in host physiology. Among others, microbiota can be involved in diet processing or protection against pathogens, both of which are potentially important in bryophagous (moss-feeding) insects, which survive on extreme diets and live in the stable environment of moss clumps suitable for the growth of fungi and bacteria. We treated Cytilus sericeus (Forster, 1771) (Coleoptera: Byrrhidae) as a model organism with bactericides and fungicides to test the effect of bacterial and fungal removal on egg hatching and larval development. Furthermore, we supplied larvae with adult feces to determine whether feces is a source of beneficial microbiota or pathogens. Bactericides had a positive effect, but fungicides had a negative effect on beetle fitness, both of which manifested during egg hatching. The feces did not play a positive role. Our conclusions indicate the presence of beneficial fungal microbiota associated with eggs but not transmitted through feces. Based on preliminary cultivation and fungicide tests, Fusarium or Penicillium may be important for suppressing pathogens, but their exact role needs to be further studied.

Caterpillar gut and host plant phylloplane mycobiomes differ: a new perspective on fungal involvement in insect guts.

Compared with the highly diverse microbiota of leaves, herbivorous insects exhibit impoverished gut microbial communities. Research to date has focused on the bacterial component of these gut microbiomes, neglecting the fungal component. As caterpillar gut bacterial microbiomes are derived strongly from their diet, we hypothesized that their mycobiomes would reflect the host leaf mycobiomes. Using the ITS2 rDNA and V5-V6 16S rRNA gene regions for DNA metabarcoding of caterpillar gut and host leaf sample pairs we compared their mycobiome genus diversity and compositions and identified genera associated with caterpillar guts. Leaves and caterpillar guts harbored different mycobiomes with quite low qualitative similarity (Jaccard index = 38.03%). The fungal genera most significantly associated with the caterpillar gut included Penicillium, Mucor and unidentified Saccharomycetales, whereas leaf-associated genera included Holtermanniella, Gibberella (teleomorph of Fusarium) and Seimatosporium. Although caterpillar gut and leaf mycobiomes had similar genus richness overall, this indicator was not correlated for individual duplets. Moreover, as more samples entered the analysis, mycobiome richness increased more rapidly in caterpillar guts than in leaves. The results suggest that the mycobiota of the caterpillar gut differs from that of their feeding substrate; further, the mycobiomes appear to be richer than the well-studied bacterial microbiotas.

Plant phylogeny drives arboreal caterpillar assemblages across the Holarctic.

Assemblages of insect herbivores are structured by plant traits such as nutrient content, secondary metabolites, physical traits, and phenology. Many of these traits are phylogenetically conserved, implying a decrease in trait similarity with increasing phylogenetic distance of the host plant taxa. Thus, a metric of phylogenetic distances and relationships can be considered a proxy for phylogenetically conserved plant traits and used to predict variation in herbivorous insect assemblages among co-occurring plant species.Using a Holarctic dataset of exposed-feeding and shelter-building caterpillars, we aimed at showing how phylogenetic relationships among host plants explain compositional changes and characteristics of herbivore assemblages.Our plant-caterpillar network data derived from plot-based samplings at three different continents included >28,000 individual caterpillar-plant interactions. We tested whether increasing phylogenetic distance of the host plants leads to a decrease in caterpillar assemblage overlap. We further investigated to what degree phylogenetic isolation of a host tree species within the local community explains abundance, density, richness, and mean specialization of its associated caterpillar assemblage.The overlap of caterpillar assemblages decreased with increasing phylogenetic distance among the host tree species. Phylogenetic isolation of a host plant within the local plant community was correlated with lower richness and mean specialization of the associated caterpillar assemblages. Phylogenetic isolation had no effect on caterpillar abundance or density. The effects of plant phylogeny were consistent across exposed-feeding and shelter-building caterpillars.Our study reveals that distance metrics obtained from host plant phylogeny are useful predictors to explain compositional turnover among hosts and host-specific variations in richness and mean specialization of associated insect herbivore assemblages in temperate broadleaf forests. As phylogenetic information of plant communities is becoming increasingly available, further large-scale studies are needed to investigate to what degree plant phylogeny structures herbivore assemblages in other biomes and ecosystems.

Quantitative assessment of plant-arthropod interactions in forest canopies: A plot-based approach.

Research on canopy arthropods has progressed from species inventories to the study of their interactions and networks, enhancing our understanding of how hyper-diverse communities are maintained. Previous studies often focused on sampling individual tree species, individual trees or their parts. We argue that such selective sampling is not ideal when analyzing interaction network structure, and may lead to erroneous conclusions. We developed practical and reproducible sampling guidelines for the plot-based analysis of arthropod interaction networks in forest canopies. Our sampling protocol focused on insect herbivores (leaf-chewing insect larvae, miners and gallers) and non-flying invertebrate predators (spiders and ants). We quantitatively sampled the focal arthropods from felled trees, or from trees accessed by canopy cranes or cherry pickers in 53 0.1 ha forest plots in five biogeographic regions, comprising 6,280 trees in total. All three methods required a similar sampling effort and provided good foliage accessibility. Furthermore, we compared interaction networks derived from plot-based data to interaction networks derived from simulated non-plot-based data focusing either on common tree species or a representative selection of tree families. All types of non-plot-based data showed highly biased network structure towards higher connectance, higher web asymmetry, and higher nestedness temperature when compared with plot-based data. Furthermore, some types of non-plot-based data showed biased diversity of the associated herbivore species and specificity of their interactions. Plot-based sampling thus appears to be the most rigorous approach for reconstructing realistic, quantitative plant-arthropod interaction networks that are comparable across sites and regions. Studies of plant interactions have greatly benefited from a plot-based approach and we argue that studies of arthropod interactions would benefit in the same way. We conclude that plot-based studies on canopy arthropods would yield important insights into the processes of interaction network assembly and dynamics, which could be maximised via a coordinated network of plot-based study sites.

Phenotypic plasticity in specialists: How long-spined larval Sympetrum depressiusculum (Odonata: Libellulidae) responds to combined predator cues.

Phenotypic plasticity is a common defensive strategy in species experiencing variable predation risk, such as habitat generalists. Larvae of generalist dragonflies can elongate their abdominal spines in environments with fish, but long spines render larvae susceptible to invertebrate predators. Long-spined specialists adapted to fish-heavy habitats are not expected to have phenotypic plasticity in this defence trait, but no empirical studies have been undertaken. Moreover, in comparison to prey responding to multiple predators that induce similar phenotypes, relatively little is known regarding how species react to combinations of predators that favour opposing traits. We examined plasticity of larval dragonfly Sympetrum depressiusculum, a long-spined habitat specialist. In a rearing experiment, larvae were exposed to four environments: (i) no predator control, (ii) fish cues (Carassius auratus), (iii) invertebrate cues (Anax imperator), as well as (iv) a combination of (ii) and (iii). Compared with the control, fish but not invertebrate cues resulted in longer spines for two (one lateral, one dorsal) of the six spines measured. Interestingly, the combined-cue treatment led to the elongation of all four dorsal spines compared with the fish treatment alone, whereas lateral spines showed no response. Our experiment provided evidence of morphological plasticity in a long-spined specialist dragonfly. We showed that nearly all spines can elongate, but also react differently under specific predator settings. Therefore, while spine plasticity evolved in direct response to a single predator type (fish), plasticity was maintained against invertebrate predators as long as fish were also present. Selective spine induction under the combined condition suggests that S. depressiusculum can successfully survive in environments with both predators. Therefore, phenotypic plasticity may be an effective strategy for habitat generalists and specialists. Although more studies are necessary to fully understand how selection shapes the evolution of phenotypic plasticity, we demonstrated that in dragonflies, presence or absence of a specific predator is not the only factor that determines plastic defence responses.

Vertical canopy gradient shaping the stratification of leaf-chewer-parasitoid interactions in a temperate forest.

Knowledge about herbivores and their parasitoids in forest canopies remains limited, despite their diversity and ecological importance. Thus, it is important to understand the factors that shape the herbivore-parasitoid community structure, particularly the effect of vertical gradient. We investigated a quantitative community dataset of exposed and semiconcealed leaf-chewing larvae and their parasitoids along a vertical canopy gradient in a temperate forest. We sampled target insects using an elevated work platform in a 0.2 ha broadleaf deciduous forest plot in the Czech Republic. We analyzed the effect of vertical position among three canopy levels (first [lowest], second [middle], and third [highest]) and tree species on community descriptors (density, diversity, and parasitism rate) and food web structure. We also analyzed vertical patterns in density and parasitism rate between exposed and semiconcealed hosts, and the vertical preference of the most abundant parasitoid taxa in relation to their host specificity. Tree species was an important determinant of all community descriptors and food web structure. Insect density and diversity varied with the vertical gradient, but was only significant for hosts. Both host guilds were most abundant in the second level, but only the density of exposed hosts declined in the third level. Parasitism rate decreased from the first to third level. The overall parasitism rate did not differ between guilds, but semiconcealed hosts suffered lower parasitism in the third level. Less host-specific taxa (Ichneumonidae, Braconidae) operated more frequently lower in the canopy, whereas more host-specific Tachinidae followed their host distribution. The most host-specific Chalcidoidea preferred the third level. Vertical stratification of insect density, diversity, and parasitism rate was most pronounced in the tallest tree species. Therefore, our study contradicts the general paradigm of weak arthropod stratification in temperate forest canopies. However, in the network structure, vertical variation might be superseded by variation among tree species.

Performance of DNA metabarcoding, standard barcoding, and morphological approach in the identification of host-parasitoid interactions.

Understanding interactions between herbivores and parasitoids is essential for successful biodiversity protection and monitoring and for biological pest control. Morphological identifications employ insect rearing and are complicated by insects' high diversity and crypsis. DNA barcoding has been successfully used in studies of host-parasitoid interactions as it can substantially increase the recovered real host-parasitoid diversity distorted by overlooked species complexes, or by species with slight morphological differences. However, this approach does not allow the simultaneous detection and identification of host(s) and parasitoid(s). Recently, high-throughput sequencing has shown high potential for surveying ecological communities and trophic interactions. Using mock samples comprising insect larvae and their parasitoids, we tested the potential of DNA metabarcoding for identifying individuals involved in host-parasitoid interactions to different taxonomic levels, and compared it to standard DNA barcoding and morphological approaches. For DNA metabarcoding, we targeted the standard barcoding marker cytochrome oxidase subunit I using highly degenerate primers, 2*300 bp sequencing on a MiSeq platform, and RTAX classification using paired-end reads. Additionally, using a large host-parasitoid dataset from a Central European floodplain forest, we assess the completeness and usability of a local reference library by confronting the number of Barcoding Index Numbers obtained by standard barcoding with the number of morphotypes. Overall, metabarcoding recovery was high, identifying 92.8% of the taxa present in mock samples, and identification success within individual taxonomic levels did not significantly differ among metabarcoding, standard barcoding, and morphology. Based on the current local reference library, 39.4% parasitoid and 90.7% host taxa were identified to the species level. DNA barcoding estimated higher parasitoid diversity than morphotyping, especially in groups with high level of crypsis. This study suggests the potential of metabarcoding for effectively recovering host-parasitoid diversity, together with more accurate identifications obtained from building reliable and comprehensive reference libraries, especially for parasitoids.

Fungal Associates of the Xylosandrus compactus (Coleoptera: Curculionidae, Scolytinae) Are Spatially Segregated on the Insect Body.

Studies of symbioses have traditionally focused on explaining one-to-one interactions between organisms. In reality, symbioses are often much more dynamic. They can involve many interacting members, and change depending on context. In studies of the ambrosia symbiosis-the mutualism between wood borer beetles and fungi-two variables have introduced uncertainty when explaining interactions: imprecise symbiont identification, and disregard for anatomical complexity of the insects. The black twig borer, Xylosandrus compactus Eichhoff, is a globally invasive ambrosia beetle that infests >200 plant species. Despite many studies on this beetle, reports of its primary symbionts are conflicting. We sampled adult X. compactus and infested plant material in central Florida to characterize the fungal symbiont community using dilution series, beetle partitioning, and DNA-based identification. X. compactus was consistently associated with two fungal taxa, Fusarium spp. and Ambrosiella xylebori Multivariate analyses revealed that A. xylebori was strongly associated with the beetle mycangium while Fusarium spp. were associated with the abdomen and external surfaces. The Fusarium spp. carried by X. compactus are not members of the Ambrosia Fusarium Clade, and are probably not mutualists. Fungal community composition of the mycangium was less variable than external body surfaces, thus providing a more consistent fungal inoculum. This is the first report of spatial partitioning as a mechanism for maintenance of a multimember ambrosia fungus community. Our results provide an explanation for discrepancies among previous reports, and suggest that conflicting results are not due to differences in symbiont communities, but due to inconsistent and incomplete sampling.