A curated DNA barcode reference library for parasitoids of northern European cyclically outbreaking geometrid moths.

Large areas of forests are annually damaged or destroyed by outbreaking insect pests. Understanding the factors that trigger and terminate such population eruptions has become crucially important, as plants, plant-feeding insects, and their natural enemies may respond differentially to the ongoing changes in the global climate. In northernmost Europe, climate-driven range expansions of the geometrid moths Epirrita autumnata and Operophtera brumata have resulted in overlapping and increasingly severe outbreaks. Delayed density-dependent responses of parasitoids are a plausible explanation for the 10-year population cycles of these moth species, but the impact of parasitoids on geometrid outbreak dynamics is unclear due to a lack of knowledge on the host ranges and prevalences of parasitoids attacking the moths in nature. To overcome these problems, we reviewed the literature on parasitism in the focal geometrid species in their outbreak range and then constructed a DNA barcode reference library for all relevant parasitoid species based on reared specimens and sequences obtained from public databases. The combined recorded parasitoid community of E. autumnata and O. brumata consists of 32 hymenopteran species, all of which can be reliably identified based on their barcode sequences. The curated barcode library presented here opens up new opportunities for estimating the abundance and community composition of parasitoids across populations and ecosystems based on mass barcoding and metabarcoding approaches. Such information can be used for elucidating the role of parasitoids in moth population control, possibly also for devising methods for reducing the extent, intensity, and duration of outbreaks.

Tick-borne pathogens in Finland: comparison of Ixodes ricinus and I. persulcatus in sympatric and parapatric areas.

BACKGROUND: Almost 3500 tick samples, originally collected via a nationwide citizen science campaign in 2015, were screened to reveal the prevalence and distribution of a wide spectrum of established and putative tick-borne pathogens vectored by Ixodes ricinus and I. persulcatus in Finland. The unique geographical distribution of these two tick species in Finland allowed us to compare pathogen occurrence between an I. ricinus-dominated area (southern Finland), an I. persulcatus-dominated area (northern Finland), and a sympatric area (central Finland). RESULTS: Of the analysed ticks, almost 30% carried at least one pathogen and 2% carried more than one pathogen. A higher overall prevalence of tick-borne pathogens was observed in I. ricinus than in I. persulcatus: 30.0% (604/2014) versus 24.0% (348/1451), respectively. In addition, I. ricinus were more frequently co-infected than I. persulcatus: 2.4% (49/2014) versus 0.8% (12/1451), respectively. Causative agents of Lyme borreliosis, i.e. bacterial genospecies in Borrelia burgdorferi (sensu lato) group, were the most prevalent pathogens (overall 17%). "Candidatus Rickettsia tarasevichiae" was found for the first time in I. ricinus ticks and in Finnish ticks in general. Moreover, Babesia divergens, B. venatorum and "Candidatus Neoehrlichia mikurensis" were reported for the first time from the Finnish mainland. CONCLUSIONS: The present study provides valuable information on the prevalence and geographical distribution of various tick-borne pathogens in I. ricinus and I. persulcatus ticks in Finland. Moreover, this comprehensive subset of ticks revealed the presence of rare and potentially dangerous pathogens. The highest prevalence of infected ticks was in the I. ricinus-dominated area in southern Finland, while the prevalence was essentially equal in sympatric and I. persulcatus-dominated areas. However, the highest infection rates for both species were in areas of their dominance, either in south or north Finland.

Crowdsourcing-based nationwide tick collection reveals the distribution of Ixodes ricinus and I. persulcatus and associated pathogens in Finland.

A national crowdsourcing-based tick collection campaign was organized in 2015 with the objective of producing novel data on tick distribution and tick-borne pathogens in Finland. Nearly 20 000 Ixodes ticks were collected. The collected material revealed the nationwide distribution of I. persulcatus for the first time and a shift northwards in the distribution of I. ricinus in Finland. A subset of 2038 tick samples containing both species was screened for Borrelia burgdorferi sensu lato (the prevalence was 14.2% for I. ricinus and 19.8% for I. persulcatus), B. miyamotoi (0.2% and 0.4%, respectively) and tick-borne encephalitis virus (TBEV; 0.2% and 3.0%, respectively). We also report new risk areas for TBEV in Finland and, for the first time, the presence of B. miyamotoi in ticks from mainland Finland. Most importantly, our study demonstrates the overwhelming power of citizen science in accomplishing a collection effort that would have been impossible with the scientific community alone.

Tick-borne bacterial pathogens in southwestern Finland.

BACKGROUND: Ixodes ricinus and Ixodes persulcatus are the main vectors of Lyme borreliosis spirochetes and several other zoonotic bacteria in northern Europe and Russia. However, few studies screening bacterial pathogens in Finnish ticks have been conducted. Therefore, reports on the occurrence and prevalence of several bacterial pathogens detected from ticks elsewhere in Europe and Russia are altogether missing from Finland. The main aim of the current study was to produce novel data on the occurrence and prevalence of several tick-borne bacterial pathogens in ticks collected from southwestern Finland. METHODS: Ticks were collected in 2013-2014 by blanket dragging from 25 localities around southwestern Finland, and additionally from a dog in Lempäälä. Collected ticks were molecularly identified and screened for Borrelia burgdorferi s.l., Borrelia miyamotoi, Rickettsia, Bartonella and Candidatus Neoehrlichia mikurensis using quantitative PCR. Furthermore, detected Rickettsia spp. were sequenced using conventional PCR to determine species. RESULTS: A total of 3169 ticks in 1174 DNA samples were screened for the listed pathogens. The most common bacteria detected was B. burgdorferi (s.l.) (18.5 % nymphal and 23.5 % adult ticks), followed by Rickettsia spp. (1.1 %; 5.1 %) and B. miyamotoi (0.51 %; 1.02 %). B. miyamotoi and Rickettsia spp. were also detected in larval samples (minimum infection rates 0.31 % and 0.21 %, respectively). Detected Rickettsia spp. were identified by sequencing as R. helvetica and R. monacensis. All screened samples were negative for Bartonella spp. and Ca. N. mikurensis. CONCLUSIONS: In the current study we report for the first time the presence of Rickettsia in Finnish ticks. Furthermore, Rickettsia spp. and B. miyamotoi were found from larval tick samples, emphasizing the importance they may have as vectors of these pathogens. Comparisons of tick density estimates and B. burgdorferi (s.l.) prevalence made between the current study and a previous study conducted in 2000 in ten out of the 25 study localities suggest that an increase in tick abundance and B. burgdorferi (s.l.) prevalence has occurred in at least some of the study localities.

Assessing the abundance, seasonal questing activity, and Borrelia and tick-borne encephalitis virus (TBEV) prevalence of Ixodes ricinus ticks in a Lyme borreliosis endemic area in Southwest Finland.

Studies have revealed that Ixodes ricinus (Acari: Ixodidae) have become more abundant and their geographical distribution extended northwards in some Nordic countries during the past few decades. However, ecological data of tick populations in Finland are sparse. In the current study, I. ricinus abundance, seasonal questing activity, and their Borrelia spp. and tick-borne encephalitis virus (TBEV) prevalence were evaluated in a Lyme borreliosis endemic area in Southwest Finland, Seili Island, where a previous study mapping tick densities was conducted 12 years earlier. A total of 1940 ticks were collected from five different biotopes by cloth dragging during May-September 2012. The overall tick density observed was 5.2 ticks/100m(2) for nymphs and adults. Seasonal questing activity of ticks differed between biotopes and life stages: bimodal occurrences were observed especially for nymphal and adult ticks in forested biotopes, while larvae in pastures exhibited mostly unimodal occurrence. Prevalence of Borrelia and TBEV in ticks was evaluated using conventional and real-time PCR. All samples were negative for TBEV. Borrelia prevalence was 25.0% for adults (n=44) and the minimum infection rate (MIR) 5.6% for pooled nymph samples (191 samples, 1-14 individuals per sample; 30/191 positive). No Borrelia were detected in pooled larval samples (63 samples, 1-139 individuals per sample). Five species of Borrelia were identified from the samples: B. afzelii, B. burgdorferi s.s., B. garinii, B. valaisiana and B. miyamotoi. In Finland, B. valaisiana and B. miyamotoi have previously been reported from the Åland Islands but not from the mainland or inner archipelago. The results of the present study suggest an increase in I. ricinus abundance on the island.

Delayed density-dependent parasitism of eggs and pupae as a contributor to the cyclic population dynamics of the autumnal moth.

Many populations of forest Lepidoptera exhibit 10-year cycles in densities, with impressive outbreaks across large regions. Delayed density-dependent interactions with natural enemies are recognized as key factors driving these cyclic population dynamics, but emphasis has typically been on the larval stages. Eggs, pupae and adults also suffer mortality from predators, parasitoids and pathogens, but little is known about possible density relationships between mortality factors and these non-feeding life stages. In a long-term field study, we experimentally deployed autumnal moth (Epirrita autumnata) eggs and pupae to their natural enemies yearly throughout the 10-year population cycle in northern Norway. The abundance of another geometrid, the winter moth (Operophtera brumata), increased in the study area, permitting comparisons between the two moth species in predation and parasitism. Survival of autumnal moth eggs and pupae was related to the moth abundance in an inverse and delayed manner. Egg and pupal parasitoids dominated as density-dependent mortality factors and predicted the subsequent growth rate of the host population size. In contrast, effects of egg and pupal predators were weakly density dependent, and generally predation remained low. Parasitism rates did not differ between the autumnal and winter moth pupae, whereas predators preferred winter moth pupae over those of the autumnal moth. We conclude that parasitism of the autumnal moth by egg and pupal parasitoids can be related to the changes of the moth density in a delayed density-dependent manner. Furthermore, egg and pupal parasitoids cannot be overlooked as causal factors for the population cycles of forest Lepidoptera in general.

Multi-objective optimization shapes ecological variation.

Ecological systems contain a huge amount of quantitative variation between and within species and locations, which makes it difficult to obtain unambiguous verification of theoretical predictions. Ordinary experiments consider just a few explanatory factors and are prone to providing oversimplified answers because they ignore the complexity of the factors that underlie variation. We used multi-objective optimization (MO) for a mechanistic analysis of the potential ecological and evolutionary causes and consequences of variation in the life-history traits of a species of moth. Optimal life-history solutions were sought for environmental conditions where different life stages of the moth were subject to predation and other known fitness-reducing factors in a manner that was dependent on the duration of these life stages and on variable mortality rates. We found that multi-objective optimal solutions to these conditions that the moths regularly experience explained most of the life-history variation within this species. Our results demonstrate that variation can have a causal interpretation even for organisms under steady conditions. The results suggest that weather and species interactions can act as underlying causes of variation, and MO acts as a corresponding adaptive mechanism that maintains variation in the traits of organisms.

Experimental test of parasitism hypothesis for population cycles of a forest lepidopteran.

Population cycles of herbivores are thought to be driven by trophic interaction mechanisms, either between food plant and herbivore or between the herbivorous prey and its natural enemies. Observational data have indicated that hymenopteran parasitoids cause delayed density-dependent mortality in cyclic autumnal moth (Epirrita autumnata) populations. We experimentally tested the parasitism hypothesis of moth population cycles by establishing a four-year parasitoid-exclusion experiment, with parasitoid-proof exclosures, parasitoid-permeable exclosures, and control plots. The exclusion of parasitoids led to high autumnal moth abundances, while the declining abundance in both the parasitoid-permeable exclosures and the control plots paralleled the naturally declining density in the study area and could be explained by high rates of parasitism. Our results provide firm experimental support for the hypothesis that hymenopteran parasitoids have a causal relationship with the delayed density-dependent component required in the generation of autumnal moth population cycles.

The evolution of female flightlessness among Ennominae of the Holarctic forest zone (Lepidoptera, Geometridae).

In order to facilitate the study of the evolution of female flightlessness among the geometrid subfamily Ennominae (Lepidoptera, Geometridae), we carried out a phylogenetic analysis based on a morphological data matrix, and DNA sequences. We used seven nuclear gene fragments, elongation factor 1alpha (EF-1alpha), wingless (wgl), isocitrate dehydrogenase (IDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein S5 (RpS5) and segments D1 and D2 of the 28S rRNA gene, and one mitochondrial gene fragment, cytochrome oxidase subunit I (COI). Sampling included 55 species of Ennominae covering all tribes with flightless females of the Holarctic boreal zone, and some other geometrids used as outgroups. Our results clearly confirmed that Ennominae (including Alsophila of the traditional subfamily Alsophilinae) is a monophyletic group, as well as supported the previously established morphology-based division of Ennominae into "ennomine" and "boarmiine" groups of genera. A number of taxonomic ambiguities were resolved but the monophyly of the traditionally recognised tribe Bistonini, comprising a number of flightless species, remained ambiguous. Bistonini is thus suggested to be subsumed to the tribe Boarmiini in the broad sense. Indeed, an analysis of timing of divergence suggested that Boarmiini s. lat. rapidly diversified in the late Oligocene/early Miocene. Within the Ennominae, seven independent origins of female flightlessness were revealed facilitating phylogenetic comparative analyses to be performed in search of causes and consequences of this phenomenon. The present phylogenetic hypothesis supports the conclusions of the "adaptive story", a hypothesis of the sequence of evolutionary events leading to flightlessness, we have presented earlier (Snäll et al., 2007). In particular, in the "boarmiine" group, the tribe Boarmiini s. lat. clearly represents a group of geometrids in which female flightlessness has evolved more frequently than in any other tribes, suggesting that this clade has likely been predisposed to evolutionary events leading to the manifestation of female flightlessness. The ancestor of the wing-reduced Ennominae has likely been a winged but slow flying forest moth feeding polyphagously on deciduous trees.

Fecundity of the autumnal moth depends on pooled geometrid abundance without a time lag: implications for cyclic population dynamics.

1. The abundance and fecundity-related body size variation of the cyclic autumnal moth Epirrita autumnata were monitored from the early increase phase and throughout the outbreak to the end of the density decline in northernmost Norway during 1999-2006. Another geometrid, the winter moth Operophtera brumata, did not increase in density until the autumnal moth had its post-peak in 2004, and was at its own peak concurrent with the steeply declining autumnal moth abundance in 2005-06. 2. The body size variables measured (forewing lengths of males and females and hind femur lengths of males) of the autumnal moth showed a similar density-dependent response, i.e. increasing density was associated with decreasing body size and fecundity. Nevertheless, regression analyses clearly ranked the pooled geometrid abundance without a time lag as the best predictor for the body size variation, ahead of the abundance of the autumnal moth or past abundance of all geometrids. 3. Nondelayed effects of lowered food quality and absolute shortage of foliage under congested conditions are the most plausible reasons for reduced body size. 4. Two most commonly proposed causal factors of the autumnal moth population cycle, i.e. delayed inducible resistance of the host plant (mountain birch Betula pubescens czerepanovii) and delayed density-dependent parasitism by specialized hymenopteran parasitoids, cannot easily explain the diverging population trends between the autumnal and winter moths. 5. We suggest that either the inducible resistance of the host tree or the host utilization of the most important parasitoids and/or pathogens have to be strictly species-specific between these closely related moth species to produce the population dynamics observed. That fecundity of the autumnal moth was best related to the pooled geometrid abundance weakens support for the former hypothesis, while our lack of host-specific information limits conclusions about the role of natural enemies.

Tree species diversity influences herbivore abundance and damage: meta-analysis of long-term forest experiments.

Plant monocultures are commonly believed to be more susceptible to herbivore attacks than stands composed of several plant species. However, few studies have experimentally tested the effects of tree species diversity on herbivory. In this paper, we present a meta-analysis of uniformly collected data on insect herbivore abundance and damage on three tree species (silver birch, black alder and sessile oak) from seven long-term forest diversity experiments in boreal and temperate forest zones. Our aim was to compare the effects of forest diversity on herbivores belonging to different feeding guilds and inhabiting different tree species. At the same time we also examined the variation in herbivore responses due to tree age and sampling period within the season, the effects of experimental design (plot size and planting density) and the stability of herbivore responses over time. Herbivore responses varied significantly both among insect feeding guilds and among host tree species. Among insect feeding guilds, only leaf miner densities were consistently lower and less variable in mixed stands as compared to tree monocultures regardless of the host tree species. The responses of other herbivores to forest diversity depended largely on host tree species. Insect herbivory on birch was significantly lower in mixtures than in birch monocultures, whereas insect herbivory on oak and alder was higher in mixtures than in oak and alder monocultures. The effects of tree species diversity were also more pronounced in older trees, in the earlier part of the season, at larger plots and at lower planting density. Overall our results demonstrate that forest diversity does not generally and uniformly reduce insect herbivory and suggest instead that insect herbivore responses to forest diversity are highly variable and strongly dependent on the host tree species and other stand characteristics as well as on the type of the herbivore.

Reduction in size and fecundity of the autumnal moth, Epirrita autumnata, in the increase phase of a population cycle.

Increasing fecundity with increasing density has been observed for many cyclic herbivore populations, including some forest Lepidoptera. We monitored population density, body size and reproductive capacity of the cyclic lepidopteran, the autumnal moth (Epirrita autumnata, Geometridae), from the early increase phase to the devastating outbreak density in northernmost Norway. Larval density of the species increased exponentially from 1998 to 2002 and remained at the outbreak level also in 2003. Within the same period, the body size and fecundity of individuals reduced as analysed from several parallel datasets on larvae, pupae and adults. In another study area in northernmost Finland, the density increase of the autumnal moth was moderate only, and true outbreak density was not attained during the study. Despite that, a reduction was again detected in the size and fecundity of individuals. Possible factors responsible for the reduced size and fecundity of individuals in the Norwegian population were quantitative shortage of foliage, rapid and delayed inducible resistances of the host, mountain birch (Betula pubescens ssp. czerepanovii), as well as crowding-induced responses of larvae. These factors likely acted in concert, although non-delayed responses to the density were emphasized. Our findings did not support the hypotheses of climatic release, inducible susceptibility of the host tree and mast depression (i.e. lowered chemical defence of the host tree after its mast seeding) as promoters of the fecundity-based density increase of the autumnal moth, since the reduced fecundity in relation to increased density was strongly against the predictions of these hypotheses. Therefore, we suggest that the density increase of autumnal moth populations is promoted by high survival rather than exceptionally high fecundity.

Performance of the cyclic autumnal moth, Epirrita autumnata, in relation to birch mast seeding.

The mast depression hypothesis has been put forward to explain the 9- to 10-year population cycle of the autumnal moth (Epirrita autumnata; Lepidoptera: Geometridae) in northern Fennoscandia. We analysed long-term data from Finnish Lapland in order to evaluate the critical assumption of the mast depression hypothesis: that better individual performance of herbivores, followed by high annual growth rate of populations, occurs in the year following mast seeding of the host, the mountain birch ( Betula pubescens ssp. czerepanovii). Since mast seeding has been suggested to occur at the expense of chemical defence against herbivores, we bioassayed the quality of birch leaves from the same trees by means of yearly growth trials with autumnal moth larvae. We also measured the size of wild adults as a determinant of potential fecundity of the species in different years. The relative growth rate of larvae was poorer in post-mast years compared to other years, rather than better as assumed by the hypothesis. Conversely, a slight indication of the increase in potential fecundity was observed due to the somewhat larger adult size in post-mast years. Population growth rate estimates, however, showed that the increase in fecundity would have to be much higher to facilitate population increase towards a cyclic peak with outbreak density. Accordingly, our two data sets do not support the assumption of a higher annual growth rate in autumnal moth populations subsequent to mast seeding of the host, thereby contradicting the predictions of the mast depression hypothesis. Temperatures, when indexed by the North Atlantic Oscillation and accumulated thermal sums, were observed to correlate with the abundance or rate of population change of the autumnal moth. The factors underlying the regular population cycles of the autumnal moth, however, remain unidentified. Overall, we suggest that the causal agents in cyclic insect population dynamics should be clarified by field experimentation, since trophic interactions are complex and are further modified by abiotic factors such as climate.

Spatial responses of two herbivore groups to a geometrid larva on mountain birch.

Direct or plant-mediated interactions between herbivores may modify their spatial distribution among and within plants. In this study, we examined the effect of a leaf-chewing geometrid, the autumnal moth (Epirrita autumnata), on two different herbivore groups, leaf rolling Deporaus betulae weevils and Eriocrania spp. leafminers, both feeding on mountain birch (Betula pubescens ssp. czerepanovii). The exact locations of herbivores within tree canopies were mapped during three successive summers. In the first 2 years, some trees were artificially colonized by eggs of the autumnal moth to induce both rapid and delayed resistance in the foliage. The natural infection levels of the pathogenic rust fungus (Melampsoridium betulinum), potentially involved in species interactions, were also recorded. At the level of the whole tree, the density of D. betulae leaf rolls was lower in trees infested by the autumnal moth in the same year. However, the feeding locations within trees were partly segregated: D. betulae favoured shadier branches, while E. autumnata preferred the sunny parts of the canopy. The autumnal moth did not affect current- or following-year density of leafminers at the tree or branch level. Trees infected by rust had fewer leafminers in the same summer than noninfected trees. There were no interaction effects between defoliation by the autumnal moth and rust infection, and no delayed effects on the abundance of other herbivores the following year. Taken together, these findings suggest that the autumnal moth has a negative, partially plant-mediated impact on D. betulae, and can reduce the extent of current-year defoliation caused by D. betulae. This may be beneficial for the mountain birch, since the greater part of D. betulae damage occurs around or after the end of the larval period of the autumnal moth, which may be a critical time for tree recovery after moth outbreaks.

Delayed induced changes in the biochemical composition of host plant leaves during an insect outbreak.

In birch, Betula pubescens, herbivore-induced delayed induced resistance (DIR) of defoliated trees may cause a strong reduction in the potential fecundity of a geometrid folivore Epirrita autumnata. In this study, we examined the biochemical basis of DIR in birch leaves during a natural outbreak of E. autumnata. A set of experimental trees was defoliated at four sites by wild larvae in the peak year of the outbreak, whereas control trees were protected from defoliation by spraying with an insecticide. The biochemical composition of leaves was analysed in the following year and, although the DIR response was weak during this outbreak, causing less than a 20% reduction in the potential fecundity of E. autumnata, some consistent relationships between defoliation, biochemistry and pupal mass of E. autumnata suggested a general biochemical basis for the defoliation-induced responses in birch leaves. Total concentrations of nitrogen, sugars and acetone-insoluble residue (e.g. cell wall polysaccharides, cell-wall-bound phenolics, protein, starch, lignin and hemicellulose) were consistently lower, and total concentrations of phenolics, especially of gallotannins and soluble proanthocyanidins, were higher in the leaves of trees defoliated in the previous year than in those protected from defoliation. The capacity of tannins to precipitate proteins correlated with contents of gallotannins, and was highest in defoliated trees. The pupal mass of E. autumnata showed a strong, positive correlation with concentrations of nitrogen and sugars, and a negative correlation with the acetone-insoluble residue and gallotannins in foliage. Correlations with other measured biochemical traits were weak. The correlation coefficients between biochemical traits and pupal mass consistently had similar signs for both defoliated and insecticide-sprayed trees, suggesting that variation in leaf quality due to defoliation in the previous year was based on similar biochemical traits as variation for other reasons. We suggest that DIR is associated with reduced growth activity of leaves, and may be seen as a delay in the biochemical maturation of leaves in defoliated trees. This explains the high concentration of gallotannins in defoliated trees, a characteristic feature of young leaves. However, the lower content of nitrogen and the higher content of soluble proanthocyanidins in defoliated trees are traits usually characterising mature, not young, leaves, indicating defoliation-induced changes in chemistry in addition to modified leaf age. Our results emphasise the importance of understanding the natural changes in chemistry during leaf maturation when interpreting defoliation-induced changes in leaf biochemistry.

Consequences of herbivory in the mountain birch (Betula pubescens ssp tortuosa): importance of the functional organization of the tree.

Three types of experiments indicate that the functional organization of the mountain birch may influence the ways in which the tree responds to simulated or natural herbivory. The first experiment showed that herbivory to both short and long shoot leaves affects plant development but, because growth largely proceeds by resources of the previous year, is manifested only in the year following the damage. The second experiment showed that even partial damage to a single long shoot leaf caused the axillary bud of that leaf to produce a shorter shoot the next year. Therefore, the value of a leaf depends also on the organ which it is subtending. In the third experiment we manipulated the apical dominance of shoots in ramets and caused improvement to leaf quality in extant shoots. Ramets within a tree responded individually, probably mediated by disturbance of the hormonal control because removal of apical buds elicited the response although removal of the same number of basal buds did not. Induced amelioration is a different response to induced resistance. The two responses are triggered by different cues and may occur in the same plant. By altering hormonal balance of shoots it is potentially possible for herbivores to induce amelioration of food quality. The ways in which herbivory is simulated may explain variability of results obtained when herbivory-induced responses in plants have been studied.