Genetic isolation through time: allochronic differentiation of a phenologically atypical population of the pine processionary moth

Seasonal differences in the timing of flight between the invasive winter moth and native Bruce spanworm promotes reproductive isolation

The European winter moth, Operophtera brumata L. (Lepidoptera: Geometridae), was accidentally introduced to North America on at least 4 separate occasions, where it has been hybridizing with the native Bruce spanworm, O. bruceata Hulst, at rates up to 10% per year. Both species are known to respond to the same sex pheromones and to produce viable offspring, but whether they differ in the seasonal timing of their mating flights is unknown. Therefore, we collected adult male moths weekly along 2 transects in the northeastern United States and genotyped individuals using polymorphic microsatellite markers as males of these 2 species cannot be differentiated morphologically. Along each transect, we then estimated the cumulative proportions (i.e., the number of individuals out of the total collected) of each species on each calendar day. Our results indicate that there are significant differences between the species regarding their seasonal timing of flight, and these allochronic differences likely are acting to promote reproductive isolation between these 2 species. Lastly, our results suggest that the later flight observed by winter moth compared to Bruce spanworm may be limiting its inland spread in the northeastern United States because of increased exposure to extreme winter events.

Consequences of coupled barriers to gene flow for the build-up of genomic differentiation

Theory predicts that when different barriers to gene flow become coincident, their joint effects enhance reproductive isolation and genomic divergence beyond their individual effects, but empirical tests of this "coupling" hypothesis are rare. Here, we analyze patterns of gene exchange among populations of European corn borer moths that vary in the number of acting barriers, allowing for comparisons of genomic variation when barrier traits or loci are in coincident or independent states. We find that divergence is mainly restricted to barrier loci when populations differ by a single barrier, whereas the coincidence of temporal and behavioral barriers is associated with divergence of two chromosomes harboring barrier loci. Furthermore, differentiation at temporal barrier loci increases in the presence of behavioral divergence and differentiation at behavioral barrier loci increases in the presence of temporal divergence. Our results demonstrate how the joint action of coincident barrier effects leads to levels of genomic differentiation that far exceed those of single barriers acting alone, consistent with theory arguing that coupling allows indirect selection to combine with direct selection and thereby lead to a stronger overall barrier to gene flow. Thus, the state of barriers-independent or coupled-strongly influences the accumulation of genomic differentiation.

Patterns of genomic and allochronic strain divergence in the fall armyworm, Spodoptera frugiperda (J.E. Smith)

Speciation is the process through which reproductive isolation develops between distinct populations. Because this process takes time, speciation studies often necessarily examine populations within a species that are at various stages of divergence. The fall armyworm, Spodoptera frugiperda (J.E. Smith), is comprised of two strains (R = Rice & C = Corn) that serve as a novel system to explore population divergence in sympatry. Here, we use ddRADSeq data to show that fall armyworm strains in the field are largely genetically distinct, but some interstrain hybridization occurs. Although we detected F1 hybrids of both R‐ and C‐strain maternal origin, only hybrids with R‐strain mtDNA were found to contribute to subsequent generations, possibly indicating a unidirectional barrier to gene flow. Although these strains have been previously defined as “host plant‐associated,” we recovered an equal proportion of R‐ and C‐strain moths in fields dominated by C‐strain host plants. As an alternative to host‐associated divergence, we tested the hypothesis that differences in nightly activity patterns could account for reproductive isolation by genotyping temporally collected moths. Our data indicates that strains exhibit a significant shift in the timing of their nightly activities in the field. This divergence in phenology creates a prezygotic reproductive barrier that likely maintains the genetic isolation between strains. Thus, we conclude that it may be ecologically inaccurate to refer to the C‐ and R‐ strain as “host‐associated” and they should more appropriately be considered “allochronic strains.”

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.

Spatial orientation of social caterpillars is influenced by polarized light

Processionary caterpillars of Thaumetopoea pityocampa (in Europe) and Ochrogaster lunifer (in Australia) (Lepidoptera: Notodontidae) form single files of larvae crawling head-to-tail when moving to feeding and pupation sites. We investigated if the processions are guided by polarization vision. The heading orientation of processions could be manipulated with linear polarizing filters held above the leading caterpillar. Exposure to changes in the angle of polarization around the caterpillars resulted in corresponding changes in heading angles. Anatomical analysis indicated specializations for polarization vision of stemma I in both species. Stemma I has a rhabdom with orthogonal and aligned microvilli, and an opaque and rugged surface, which are optimizations for skylight polarization vision, similar to the dorsal rim of adult insects. Stemmata II-VI have a smooth and shiny surface and lobed rhabdoms with non-orthogonal and non-aligned microvilli; they are thus optimized for general vision with minimal polarization sensitivity. Behavioural and anatomical evidence reveal that polarized light cues are important for larval orientation and can be robustly detected with a simple visual system.

Structural and Functional Characterization of Allatostatin Receptor Type-C of Thaumetopoea pityocampa, a Potential Target for Next-Generation Pest Control Agents

Insect neuropeptide receptors, including allatostatin receptor type C (AstR-C), a G protein-coupled receptor, are among the potential targets for designing next-generation pesticides that despite their importance in offering a new mode-of-action have been overlooked. Focusing on AstR-C of Thaumetopoea pityocampa, a common pest in Mediterranean countries, by employing resonance energy transfer-based methods, we showed Gαi/o coupling and β-arrestin recruitment of the receptor at sub-nanomolar and nanomolar ranges of the endogenous ligand, AST-C, respectively. Molecular docking and molecular dynamics simulation studies revealed the importance of extracellular loop 2 in AstRC/AST-C interaction, and a combination of in silico and in vitro approaches showed the substantial role of Q271^6.55 in G protein-dependent activation of AstR-C possibly via contributing to the flexibility of the receptor's structure. The functional and structural insights obtained on T. pit AstR-C positively assist future efforts in developing environmentally friendly pest control agents that are needed urgently.

Resource competition shapes biological rhythms and promotes temporal niche differentiation in a community simulation

Competition for resources often contributes strongly to defining an organism's ecological niche. Endogenous biological rhythms are important adaptations to the temporal dimension of niches, but how other organisms influence such temporal niches has not been much studied, and the role of competition in particular has been even less examined. We investigated how interspecific competition and intraspecific competition for resources shape an organism's activity rhythms.To do this, we simulated communities of one or two species in an agent-based model. Individuals in the simulation move according to a circadian activity rhythm and compete for limited resources. Probability of reproduction is proportional to an individual's success in obtaining resources. Offspring may have variance in rhythm parameters, which allow for the population to evolve over time.We demonstrate that when organisms are arrhythmic, one species will always be competitively excluded from the environment, but the existence of activity rhythms allows niche differentiation and indefinite coexistence of the two species. Two species which are initially active at the same phase will differentiate their phase angle of entrainment over time to avoid each other. When only one species is present in an environment, competition within the species strongly selects for niche expansion through arrhythmicity, but the addition of an interspecific competitor facilitates evolution of increased rhythmic amplitude when combined with additional adaptations for temporal specialization. Finally, if individuals preferentially mate with others who are active at similar times of day, then disruptive selection by intraspecific competition can split one population into two reproductively isolated groups separated in activity time.These simulations suggest that biological rhythms are an effective method to temporally differentiate ecological niches and that competition is an important ecological pressure promoting the evolution of rhythms and sleep. This is the first study to use ecological modeling to examine biological rhythms.

Does hunger lead to hybridization in a genus of sexually cannibalistic insects (Orthoptera: Prophalangopsidae)?

Allochronic isolation can be a strong mechanism for reproductive isolation and speciation. However, imperfect allochrony and the expression of phenological plasticity can erode temporal barriers to gene flow and result in hybridization between divergent lineages. Here, we combine behavioural ecology and genomics to investigate this scenario in two closely related species of grigs in the genus Cyphoderris. These species exhibit a unique mating system whereby females feed on the fleshy hind wings of the male during copulation, and copulation with conspecific males is more likely in food-restricted females than in well-fed females. In western Canada, Cyphoderris buckelli and Cyphoderris monstrosa are sympatric but largely allochronically separated, with C. buckelli breeding earlier. However, their breeding seasons can overlap, leading to potential for older C. buckelli females to mate with young C. monstrosa males to obtain resources via sexual cannibalism. We used behavioural assays to test whether female feeding status affects the propensity for interspecific mating between C. buckelli females and C. monstrosa males. We then tested for hybridization and gene exchange in wild populations of both species, using morphology, mitochondrial DNA and genome-wide nuclear markers. We found that interspecific courtship and mating can occur, but the relationship between food restriction and increased propensity for hybridization was not significant. Although we observed intraspecific population genetic structure in both species, we found no signatures of hybridization in the morphological or genetic datasets, which suggests that postmating reproductive barriers might be preventing successful hybridization in the wild.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Phenology-dependent cold exposure and thermal performance of Ostrinia nubilalis ecotypes

Background

Understanding adaptation involves establishing connections between selective agents and beneficial population responses. However, relatively little attention has been paid to seasonal adaptation, in part, because it requires complex and integrative knowledge about seasonally fluctuating environmental factors, the effects of variable phenology on exposure to those factors, and evidence for temporal specialization. In the European corn borer moth, Ostrinia nubilalis, sympatric pheromone strains exploit the same host plant (Zea mays) but may genetically differ in phenology and be reproductively “isolated by time.” Z strain populations in eastern North America have been shown to have a prolonged larval diapause and produce one annual mating flight (July), whereas E strain populations complete an earlier (June) and a later (August) mating flight by shortening diapause duration. Here, we find evidence consistent with seasonal “adaptation by time” between these ecotypes.

Results

We use 12 years of field observation of adult seasonal abundance to estimate phenology of ecotype life cycles and to quantify life-stage specific climatic conditions. We find that the observed reduction of diapause duration in the E strain leads their non-diapausing, active life stages to experience a ~ 4 °C colder environment compared to the equivalent life stages in the Z strain. For a representative pair of populations under controlled laboratory conditions, we compare life-stage specific cold tolerance and find non-diapausing, active life stages in the E strain have as much as a 60% greater capacity to survive rapid cold shock. Enhanced cold hardiness appears unrelated to life-stage specific changes in the temperature at which tissues freeze.

Conclusions

Our results suggest that isolation by time and adaptation by time may both contribute to population divergence, and they argue for expanded study in this species of allochronic populations in nature experiencing the full spectrum of seasonal environments. Cyclical selective pressures are inherent properties of seasonal habitats. Diverse fluctuating selective agents across each year (temperature, predation, competition, precipitation, etc.) may therefore be underappreciated drivers of biological diversity.

Shifted phenology in the pine processionary moth affects the outcome of tree-insect interaction

In the Mediterranean and temperate regions, an increase in the frequency and intensity of drought events has been recorded, probably due to climate change. In consequence, trees will more frequently experience hydric stress, a condition that can be expected to affect insect-tree interactions, while adaptation mechanisms may be further in course. The effect of tree water stress on the performance of two allochronic populations of Thaumetopoea pityocampa was here studied. Namely, we compared a unique population of this insect, in which the larvae develop in the summer (SP), with the typical population having winter larval development (WP), to test the adaptation hypothesis to host plant status. Larvae of each population were fed on needles of young potted Pinus pinaster plants under two water supply regimes: (i) well-watered (control) and (ii) subjected to 3 months of drought stress. Compared to control, stressed plants had higher amounts of soluble sugars, phenols, and higher C/N ratio, whereas water content and chlorophylls concentrations were lower. In general, T. pityocampa larvae had lower performances on water-stressed plants, as shown by lower survival rates, lower needle consumption, and longer development times. Yet, the detrimental effects of tree stress were only significant for the WP larvae, while SP larvae were able to overcome such conditions. Results demonstrate that tree water stress can negatively affect T. pityocampa populations. Furthermore, the evidence is also provided that responses to the physiological condition of the host trees may occur at the population level, as a result of adaptation mechanisms driven by climate change.

Temporal isolation between sympatric host plants cascades across multiple trophic levels of host-associated insects

Phenological differences between host plants can promote temporal isolation among host-associated populations of insects with life cycles tightly coupled to plant phenology. Divergence in the timing of spring budbreak between two sympatric sister oak species has been shown to promote temporal isolation between host plants and their host-associated populations of a cynipid gall wasp. Here, we examined the generality of this mechanism by testing the hypothesis of cascading temporal isolation for five additional gall-formers and three natural enemy species associated with these same oak species. The timing of adult emergence from galls differed significantly between host-associated populations for all nine species and parallels the direction of the phenological differences between host plants. Differences in emergence timing can reduce gene flow between host-associated populations by diminishing mating opportunities and/or reducing the fitness of immigrants due to differences in the availability of ephemeral resources. Our study suggests that cascading temporal isolation could be a powerful 'biodiversity generator' across multiple trophic levels in tightly coupled plant-insect systems.

Environmental Adaptations, Ecological Filtering, and Dispersal Central to Insect Invasions

Insect invasions, the establishment and spread of nonnative insects in new regions, can have extensive economic and environmental consequences. Increased global connectivity accelerates rates of introductions, while climate change may decrease the barriers to invader species' spread. We follow an individual-level insect- and arachnid-centered perspective to assess how the process of invasion is influenced by phenotypic heterogeneity associated with dispersal and stress resistance, and their coupling, across the multiple steps of the invasion process. We also provide an overview and synthesis on the importance of environmental filters during the entire invasion process for the facilitation or inhibition of invasive insect population spread. Finally, we highlight important research gaps and the relevance and applicability of ongoing natural range expansions in the context of climate change to gain essential mechanistic insights into insect invasions.

A total evidence phylogeny for the processionary moths of the genus Thaumetopoea (Lepidoptera: Notodontidae: Thaumetopoeinae)

Processionary moths belong to a group of about 100 species well known for their social behaviour and their urticating systems. The genus Thaumetopoea s.l. includes about 15 species and has been divided into three genera (Helianthocampa, Thaumetopoea s.s., and Traumatocampa) in the last revision, based on key morphological features of the adults and on the host plants of the larvae. We performed a total evidence approach to resolve the phylogeny of the genus Thaumetopoea s.l., analysing all valid taxa included in this group, plus a broad array of close relatives. Thaumetopoea was monophyletic and supported by several apomorphies. Further subclades corroborated by synapomorphies were identified. Our phylogeny suggests that Thaumetopoea must be regarded as a single genus. The mapping of key life history traits on the total evidence tree allowed us to sketch a plausible identikit of the Thaumetopoea ancestor and to track the evolution of the genus. The ancestor originated in the eastern Mediterranean area, and used broadleaved host plants. Subsequently, a switch to conifers occurred, just once, in a large subclade. The ancestor pupated in the soil, like several current species, but in a few taxa this trait was lost, together with the related morphological adaptations.

Effect of heat waves on embryo mortality in the pine processionary moth

Extreme climate events such as heat waves are predicted to become more frequent with climate change, representing a challenge for many organisms. The pine processionary moth Thaumetopoea pityocampa is a Mediterranean pine defoliator, which typically lays eggs during the summer. We evaluated the effects of heat waves on egg mortality of three populations with different phenologies: a Portuguese population with a classical life cycle (eggs laid in summer), an allochronic Portuguese population reproducing in spring, and a Tunisian population from the extreme southern limit of T. pityocampa distribution range, in which eggs are laid in fall. We tested the influence of three consecutive hot days on egg survival and development time, using either constant (CT) or daily cycling temperatures (DT) with equivalent mean temperatures. Maximum temperatures (T max) used in the experiment ranged from 36 to 48°C for DT and from 30 to 42°C for CT. Heat waves had a severe negative effect on egg survival when T max reached 42°C for all populations. No embryo survived above this threshold. At high mean temperatures (40°C), significant differences were observed between populations and between DT and CT regimes. Heat waves further increased embryo development time. The knowledge we gained about the upper lethal temperature to embryos of this species will permit better prediction of the potential expansion of this insect under different climate warming scenarios.

Sympatric, temporally isolated populations of the pine white butterfly Neophasia menapia, are morphologically and genetically differentiated

Temporal isolation remains an understudied, and potentially under-appreciated, mechanism of reproductive isolation. Phenological differences have been discovered in populations of the pine white butterfly (Neophasia menapia), a typically univoltine species found throughout western North America. At two locations in the Coast Range of California there are two periods of adult emergence per year, one in early summer (July) and one in late summer/autumn (September/October). Differences in flight time are accompanied by differences in wing shape and pigmentation. Here we use a combination of population genomics and morphological analyses to assess the extent to which temporal isolation is able to limit gene flow between sympatric early and late flights. Not only did we detect both genetic and morphological differences between early and late flights at the two sites, we also found that the patterns of differentiation between the two flights were different at each location, suggesting an independent origin for the two sympatric flights. Additionally, we found no evidence that these sympatric flights originated via colonization from any of the other sampled localities. We discuss several potential hypotheses about the origin of these temporally isolated sympatric flights.

Experimental evidence for heritable reproductive time in 2 allochronic populations of pine processionary moth

Phenology allows organisms to overcome seasonally variable conditions through life-cycle adjustment. Changes in phenology can drastically modify the evolutionary trajectory of a population, while a shift in the reproductive time may cause allochronic differentiation. The hypothesis of heritable reproductive time was experimentally tested, by studying a unique population of the pine processionary moth Thaumetopoea pityocampa (Den. & Schiff.) which has a shifted phenology, and however co-occurs with the typical population following the classical life cycle. When populations of both types were reared under controlled conditions, the reproductive time was maintained asynchronous, as observed in the field. The shifted population was manipulated in the laboratory to reproduce later than usual, yet the offspring emerged in the next year at the expected dates thus "coming back" to the usual cycle. Hybrids from crosses performed between the 2 populations showed an intermediate phenology. From the emergence times of parents and offspring, a high heritability of the reproductive time (h = 0.76) was observed. The offspring obtained from each type of cross was genetically characterized using microsatellite markers. Bayesian clustering analysis confirmed that hybrids can be successfully identified and separated from the parental genetic classes by genotyping. Findings support the hypothesis that, for this particular population, incipient allochronic speciation is due to a heritable shift in the reproductive time that further causes assortative mating and might eventually cause ecological adaptation/maladaptation in response to environmental changes.

Genetic differentiation of the pine processionary moth at the southern edge of its range: contrasting patterns between mitochondrial and nuclear markers

The pine processionary moth (Thaumetopoea pityocampa) is an important pest of coniferous forests at the southern edge of its range in Maghreb. Based on mitochondrial markers, a strong genetic differentiation was previously found in this species between western (pityocampa clade) and eastern Maghreb populations (ENA clade), with the contact zone between the clades located in Algeria. We focused on the moth range in Algeria, using both mitochondrial (a 648 bp fragment of the tRNA-cox2) and nuclear (11 microsatellite loci) markers. A further analysis using a shorter mtDNA fragment and the same microsatellite loci was carried out on a transect in the contact zone between the mitochondrial clades. Mitochondrial diversity showed a strong geographical structure and a well-defined contact zone between the two clades. In particular, in the pityocampa clade, two inner subclades were found whereas ENA did not show any further structure. Microsatellite analysis outlined a different pattern of differentiation, with two main groups not overlapping with the mitochondrial clades. The inconsistency between mitochondrial and nuclear markers is probably explained by sex-biased dispersal and recent afforestation efforts that have bridged isolated populations.

Temporally isolated lineages of pink salmon reveal unique signatures of selection on distinct pools of standing genetic variation

A species' genetic diversity bears the marks of evolutionary processes that have occurred throughout its history. However, robust detection of selection in wild populations is difficult and often impeded by lack of replicate tests. Here, we investigate selection in pink salmon (Oncorhynchus gorbuscha) using genome scans coupled with inference from a haploid-assisted linkage map. Pink salmon have a strict 2-year semelparous life history which has resulted in temporally isolated (allochronic) lineages that remain sympatric through sharing of spawning habitats in alternate years. The lineages differ in a range of adaptive traits, suggesting different genetic backgrounds. We used genotyping by sequencing of haploids to generate a high-density linkage map with 7035 loci and screened an existing panel of 8036 loci for signatures of selection. The linkage map enabled identification of novel genomic regions displaying signatures of parallel selection shared between lineages. Furthermore, 24 loci demonstrated divergent selection and differences in genetic diversity between lineages, suggesting that adaptation in the 2 lineages has arisen from different pools of standing genetic variation. Findings have implications for understanding asynchronous population abundances as well as predicting future ecosystem impacts from lineage-specific responses to climate change.

Contemporary gene flow between wild An. gambiae s.s. and An. arabiensis

Background

In areas where the morphologically indistinguishable malaria mosquitoes Anopheles gambiae Giles and An. arabiensis Patton are sympatric, hybrids are detected occasionally via species-diagnostic molecular assays. An. gambiae and An. arabiensis exhibit both pre- and post-reproductive mating barriers, with swarms largely species-specific and male F1 (first-generation) hybrids sterile. Consequently advanced-stage hybrids (back-crosses to parental species), which would represent a route for potentially-adaptive introgression, are expected to be very rare in natural populations. Yet the use of one or two physically linked single-locus diagnostic assays renders them indistinguishable from F1 hybrids and levels of interspecific gene flow are unknown.

Methods

We used data from over 350 polymorphic autosomal SNPs to investigate post F1 gene flow via patterns of genomic admixture between An. gambiae and An. arabiensis from eastern Uganda. Simulations were used to investigate the statistical power to detect hybrids with different levels of crossing and to identify the hybrid category significantly admixed genotypes could represent.

Results

A range of admixture proportions were detected for 11 field-collected hybrids identified via single-locus species-diagnostic PCRs. Comparison of admixture data with simulations indicated that at least seven of these hybrids were advanced generation crosses, with backcrosses to each species identified. In addition, of 36 individuals typing as An. gambiae or An. arabiensis that exhibited outlying admixture proportions, ten were identified as significantly mixed backcrosses, and at least four of these were second or third generation crosses.

Conclusions

Our results show that hybrids detected using standard diagnostics will often be hybrid generations beyond F1, and that in our study area around 5% (95% confidence intervals 3%-9%) of apparently ‘pure’ species samples may also be backcrosses. This is likely an underestimate because of rapidly-declining detection power beyond the first two backcross generations. Post-F1 gene flow occurs at a far from inconsequential rate between An. gambiae and An. arabiensis, and, especially for traits under strong selection, could readily lead to adaptive introgression of genetic variants relevant for vector control.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-345) contains supplementary material, which is available to authorized users.

Comparative analysis of two phenologically divergent populations of the pine processionary moth (Thaumetopoea pityocampa) by de novo transcriptome sequencing

The pine processionary moth Thaumetopoea pityocampa is a Mediterranean lepidopteran defoliator that experiences a rapid range expansion towards higher latitudes and altitudes due to the current climate warming. Its phenology - the time of sexual reproduction - is certainly a key trait for the local adaptation of the processionary moth to climatic conditions. Moreover, an exceptional case of allochronic differentiation was discovered ca. 15 years ago in this species. A population with a shifted phenology (the summer population, SP) co-exists near Leiria, Portugal, with a population following the classical cycle (the winter population, WP). The existence of this population is an outstanding opportunity to decipher the genetic bases of phenology. No genomic resources were so far available for T. pityocampa. We developed a high-throughput sequencing approach to build a first reference transcriptome, and to proceed with comparative analyses of the sympatric SP and WP. We pooled RNA extracted from whole individuals of various developmental stages, and performed a transcriptome characterisation for both populations combining Roche 454-FLX and traditional Sanger data. The obtained sequences were clustered into ca. 12,000 transcripts corresponding to 9265 unigenes. The mean transcript coverage was 21.9 reads per bp. Almost 70% of the de novo assembled transcripts displayed significant similarity to previously published proteins and around 50% of the transcripts contained a full-length coding region. Comparative analyses of the population transcriptomes allowed to investigate genes specifically expressed in one of the studied populations only, and to identify the most divergent homologous SP/WP transcripts. The most divergent pairs of transcripts did not correspond to obvious phenology-related candidate genes, and 43% could not be functionally annotated. This study provides the first comprehensive genome-wide resource for the target species T. pityocampa. Many of the assembled genes are orthologs of published Lepidoptera genes, which allows carrying out gene-specific re-sequencing. Data mining has allowed the identification of SNP loci that will be useful for population genomic approaches and genome-wide scans of population differentiation to identify signatures of selection.

Phenotypic divergence in reproductive traits of a moth population experiencing a phenological shift

Allochrony that is reproductive isolation by time may further lead to divergence of reproductive adaptive traits in response to different environmental pressures over time. A unique "summer" population of the pine processionary moth Thaumetopoea pityocampa, reproductively isolated from the typical winter populations by allochronic differentiation, is here analyzed. This allochronically shifted population reproduces in the spring and develops in the summer, whereas "winter" populations reproduce in the late summer and have winter larval development. Both summer and winter populations coexist in the same pine stands, yet they face different climatic pressures as their active stages are present in different seasons. The occurrence of significant differences between the reproductive traits of the summer population and the typical winter populations (either sympatric or allopatric) is thus hypothesized. Female fecundity, egg size, egg covering, and egg parasitism were analyzed showing that the egg load was lower and that egg size was higher in the summer population than in all the studied winter populations. The scales that cover the egg batches of T. pityocampa differed significantly between populations in shape and color, resulting in a looser and darker covering in the summer population. The single specialist egg parasitoid species of this moth was almost missing in the summer population, and the overall parasitism rates were lower than in the winter population. Results suggest the occurrence of phenotypic differentiation between the summer population and the typical T. pityocampa winter populations for the life-history traits studied. This work provides an insight into how ecological divergence may follow the process of allochronic reproductive isolation.

Are heat waves susceptible to mitigate the expansion of a species progressing with global warming?

A number of organisms, especially insects, are extending their range in response of the increasing trend of warmer temperatures. However, the effects of more frequent climatic anomalies on these species are not clearly known. The pine processionary moth, Thaumetopoea pityocampa, is a forest pest that is currently extending its geographical distribution in Europe in response to climate warming. However, its population density largely decreased in its northern expansion range (near Paris, France) the year following the 2003 heat wave. In this study, we tested whether the 2003 heat wave could have killed a large part of egg masses. First, the local heat wave intensity was determined. Then, an outdoor experiment was conducted to measure the deviation between the temperatures recorded by weather stations and those observed within sun-exposed egg masses. A second experiment was conducted under laboratory conditions to simulate heat wave conditions (with night/day temperatures of 20/32°C and 20/40°C compared to the control treatment 13/20°C) and measure the potential effects of this heat wave on egg masses. No effects were noticed on egg development. Then, larvae hatched from these egg masses were reared under mild conditions until the third instar and no delayed effects on the development of larvae were found. Instead of eggs, the 2003 heat wave had probably affected directly or indirectly the young larvae that were already hatched when it occurred. Our results suggest that the effects of extreme climatic anomalies occurring over narrow time windows are difficult to determine because they strongly depend on the life stage of the species exposed to these anomalies. However, these effects could potentially reduce or enhance the average warming effects. As extreme weather conditions are predicted to become more frequent in the future, it is necessary to disentangle the effects of the warming trend from the effects of climatic anomalies when predicting the response of a species to climate change.

Juvenile habitat partitioning and relative productivity in allochronically isolated sockeye salmon (Oncorhynchus nerka)

Allochronic divergence, like spatial isolation, may contribute to population diversity and adaptation, however the challenges for tracking habitat utilization in shared environments are far greater. Adult Klukshu River (Yukon, Canada) sockeye salmon, Oncorhynchus nerka, return as genetically distinct “early” and “late” runs. Early and late adult spawning populations (1999 and 2000) and their subsequent fry (sampled at 7 sites in 2000 and at 8 sites in 2001 throughout Klukshu Lake and River) were genotyped at eight microsatellite loci. Bayesian assignment was used to determine the spatial distribution of early versus late fry; although intermixed, the distribution of fry significantly differed in Klukshu Lake and in the Klukshu River in 2001, based on crosstab analyses. Late-run fry predominated in Klukshu Lake at all sites, while early-run fry were most common in the north and south of Klukshu Lake and in Klukshu River. Early-run spawners had significantly higher relative productivity (early life survival) than late-run fish (2.9 times more fry produced per early-run adult in 2000, and 9.2 times more in 2001). This study demonstrates spatial habitat partitioning and differences in the contribution of allochronically isolated populations to fry abundance, and highlights annual variability that likely contributes to recruitment variation.

Parallel allochronic divergence in a winter moth due to disruption of reproductive period by winter harshness

The disruption of reproductive timing by climatic harshness may result in the temporal isolation of conspecific populations and, ultimately, in speciation. However, whether temporal isolation alone can act as the force initiating speciation and how often the same type of climatic disruption results in the divergence of allochronic populations in a lineage are largely unknown. The reproductive period of the winter geometrid moth Inurois punctigera is separated into early and late winter in habitats with severe winters, but not in habitats with mild winters, suggesting that the reproductive season is disrupted by the harshness of the mid-winter period. Here, we show that sympatric pairs of early- and late-winter populations that differ in origin exist in different regions, suggesting a parallel divergence of reproductive timing. In each region, significant genetic differentiation exists between these early- and late-winter populations, suggesting that the temporal reproductive isolation has persisted. Moreover, we demonstrate that the temporal isolation, in comparison with geographic isolation, contributes greatly to the genetic differentiation among geographic and temporal populations by an analysis of molecular variance and by a comparison of genetic differentiations (F(ST) ) between geographic populations with and without difference in reproductive season. Our results suggest that adaptive divergence of allochronically reproducing populations has occurred independently in different regions, implying the generality of the role of temporal isolation in initiating speciation in a winter moth lineage.

Developmental plasticity as a cohesive evolutionary process between sympatric alternate-year insect cohorts

Many species, particularly insects, pass through a series of distinct phases during their life history, with the developmental timing directed towards appropriate resources. Any factor that creates variation in developmental timing may partition a population into discrete populations-or 'cohorts'. Where there is continued failure to recruit outside the natal cohort then alternate cohorts will have their own internal dynamics, eventually leading to independent demographic and evolutionary trajectories. By contrast, continued variation in development rates within a cohort-cohort splitting-may homogenise otherwise independent demographic units. Using a panel of 14 microsatellite loci, we quantify the genetic signature of apparent demographic isolation between coexisting, but alternate, semivoltine cohorts of the damselfly Coenagrion mercuriale at locations that span its distribution in the UK. We find consistently low levels of genetic divergence between sympatric cohorts of C. mercuriale, indicative of developmental plasticity during the larval stage (unregulated development) whereby some individuals complete their development outside the predominant 2-year (semivoltine) period. Thus, individuals that alter their developmental rate successfully recruit to a different cohort. Despite maintaining contrasting population sizes, gene flow between alternate cohorts broadly is sufficient to place them on a similar evolutionary trajectory and also buffers against loss of genetic diversity. Such flexible larval development permits a response to local conditions and may facilitate response to environmental change.

Genetically regulated temporal variation of novel courtship elements in the Hawaiian cricket genus Laupala

The Hawaiian cricket genus Laupala (Gryllidae: Trigonidiinae) has undergone rapid and extensive speciation, with divergence in male song and female acoustic preference playing a role in maintaining species boundaries. Recent study of interspecific differences in the diel rhythmicity of singing and mating, suggests that temporal variation in behavior may reduce gene flow between species. In addition, Laupala perform an elaborate and protracted courtship, providing potential for further temporal variation. However, whether these behavioral differences have a genetic basis or result from environmental variation is unknown. We observed courtship and mating in a common garden study of the sympatric species, Laupala cerasina and Laupala paranigra. We document interspecific differences in the onset and duration of courtship, spermatophore production rate, and diel mating rhythmicity. Our study demonstrates a genetic contribution to interspecific behavioral differences, and suggests an evolutionary pathway to the origins of novel timing phenotypes.

Short- and long-day responses in the pre-adult developmental duration of two species of Camponotus ants

We assessed the effect of different day/night lengths on the pre-adult developmental time of two species of Camponotus ants that normally develop in dark underground nests. We assayed larval (egg-to-pupal formation), pupal (pupal formation-to-adult emergence), and pre-adult (egg-to-adult emergence) durations in these ants under three different light/dark (LD) cycles of 12:12 h, 10:14 h, and 14:10 h. We observed that the pre-adult development time of ants under these day lengths was significantly different. Although both species developed fastest under 12:12 h LD, when asymmetric LD cycles were compared, night-active species (Camponotus compressus) developed faster under short days (10:14 h) and day-active species (C. paria) developed faster under long days (14:10 h). This day/night-length-mediated difference in pre-adult developmental duration was mostly due to modulation of larval duration; however, in day-active species it was also via altered pupal duration. These results thus indicate that the two species of Camponotus ants respond differently to short and long days, suggesting that seasonal timers regulate pre-adult development time in tropical ant species living in dark underground nests.

Incipient allochronic speciation by climatic disruption of the reproductive period

Disruptive selection of life-cycle timing may cause temporal isolation directly and, ultimately, allochronic speciation. Despite the fact that segregation of the reproductive period among related species has been broadly observed across taxa, it remains controversial whether temporal isolation can function as the primary process of speciation. In the Japanese winter geometrid moth Inurois punctigera, allochronic divergence has resulted from climatic disruption of the reproductive period. In habitats with severe midwinter, two sympatric groups of moth reproduce allochronically in early and late winter. These groups are genetically diverging sister lineages and now co-occur allochronically throughout Japan. By contrast, in habitats with milder midwinter these lineages form a continuous adult period and gene flow has been facilitated between the lineages. These results, together with the fact that there is no difference in larval host use, indicate that temporal isolation has been the sole mechanism for allochronic isolation in colder habitats and that allochrony is not a by-product of other adaptations. Thus, the allochronic divergence of sympatric I. punctigera populations represents an incipient speciation process driven by midwinter disruption of the reproductive period.

Allochronic seasonal peak activities of Dermacentor and Haemaphysalis spp. under continental climate in Hungary

Hard ticks (Acari: Ixodidae) were collected from the vegetation at monthly intervals in 2008 on six places of three different biotopes in Hungary. Except for Haemaphysalis concinna and H. punctata, predominance of females was observed in the questing population. From among the species with biphasic activity period, metastriate ticks (Dermacentor spp. and H. inermis) had higher prevalence of host-seeking males in the autumn than in the spring, as opposed to prostriate Ixodes ricinus. In comparison with the permanently mild weather experienced at the beginning of 2007, sharply rising temperatures during consecutive winter days in January and February of 2008 appeared to be more efficient in triggering a 1-2 month earlier spring peak activity of hard ticks, except for I. ricinus. Regarding seasonality, D. marginatus was most numerous in February and March, whereas D. reticulatus in September and October; H. punctata showed peak activity in March-April, H. concinna in May, and H. inermis in November-December. Within these genera such a temporal difference (allochrony) between seasonal peak activities of sympatric species under continental climate is described for the first time.

Incipient allochronic speciation due to non-selective assortative mating by flowering time, mutation and genetic drift

We model the evolution of flowering time using a multilocus quantitative genetic model with non-selective assortative mating and mutation to investigate incipient allochronic speciation in a finite population. For quantitative characters with evolutionary parameters satisfying empirical observations and two approximate inequalities that we derived, disjunct clusters in the population flowering phenology originated within a few thousand generations in the absence of disruptive natural or sexual selection. Our simulations and the conditions we derived showed that cluster formation was promoted by limited population size, high mutational variance of flowering time, short individual flowering phenology and a long flowering season. By contrast, cluster formation was hindered by inbreeding depression, stabilizing selection and pollinator limitation. Our results suggest that incipient allochronic speciation in populations of limited size (satisfying two inequalities) could be a common phenomenon.