EFFECTS OF POSTGLACIAL RANGE EXPANSION ON ALLOZYME AND QUANTITATIVE GENETIC VARIATION OF THE PITCHER-PLANT MOSQUITO, WYEOMYIA SMITHII

Evolutionary potential varies across populations and traits in the neotropical oak Quercus oleoides

Heritable variation in polygenic (quantitative) traits is critical for adaptive evolution and is especially important in this era of rapid climate change. In this study, we examined the levels of quantitative genetic variation of populations of the tropical tree Quercus oleoides Cham. and Schlect. for a suite of traits related to resource use and drought resistance. We tested whether quantitative genetic variation differed across traits, populations and watering treatments. We also tested potential evolutionary factors that might have shaped such a pattern: selection by climate and genetic drift. We measured 15 functional traits on 1322 1-year-old seedlings of 84 maternal half-sib families originating from five populations growing under two watering treatments in a greenhouse. We estimated the additive genetic variance, coefficient of additive genetic variation and narrow-sense heritability for each combination of traits, populations and treatments. In addition, we genotyped a total of 119 individuals (with at least 20 individuals per population) using nuclear microsatellites to estimate genetic diversity and population genetic structure. Our results showed that gas exchange traits and growth exhibited strikingly high quantitative genetic variation compared with traits related to leaf morphology, anatomy and photochemistry. Quantitative genetic variation differed between populations even at geographical scales as small as a few kilometers. Climate was associated with quantitative genetic variation, but only weakly. Genetic structure and diversity in neutral markers did not relate to coefficient of additive genetic variation. Our study demonstrates that quantitative genetic variation is not homogeneous across traits and populations of Q. oleoides. More importantly, our findings suggest that predictions about potential responses of species to climate change need to consider population-specific evolutionary characteristics.

When do trade-offs occur? The roles of energy constraints and trait flexibility in bushcricket populations

In many animal species, the expression of sexually selected traits is negatively correlated with traits associated with survival such as immune function, a relationship termed a 'trade-off'. But an alternative in which sexually selected traits are positively correlated with survival traits is also widespread. We propose that the nature of intertrait relationships is largely determined by overall energy expenditure, energy availability and trait flexibility, with trade-offs expected when individuals are subject to energy constraints. We tested this hypothesis in Ephippiger diurnus, a European bushcricket in which males are distinguished by two prominent sexually selected traits, acoustic calls and a large spermatophore transferred to the female at mating, and where immune function may be critical in survival. Ephippiger diurnus are distributed as small, isolated populations that are differentiated genetically and behaviourally. We analysed songs, spermatophores and the immune function in male individuals from eight populations spanning a range of song types. As predicted, we only found trade-offs in those populations that expended more energy on song and were less flexible in their ability to adjust that expenditure. Ultimately, energy constraints and resulting trade-offs may limit the evolution of song exaggeration in E. diurnus populations broadcasting long calls comprised of multiple 'syllables'.

Evolutionary transition from blood feeding to obligate nonbiting in a mosquito

The evolutionary transformation from a blood-feeding to an obligate nonbiting lifestyle is occurring uniquely within the genetic background of a single species of mosquito, Wyeomyia smithii, as a product of selection in nature. Associated genetic changes in metabolic pathways indicate a high anticipatory metabolic investment prior to consuming blood, presumably balanced by the reproductive benefits from an imminent blood meal. This evolutionary transformation provides a starting point for determining pivotal upstream genetic changes between biters and nonbiters and for identifying universal nonbiting genes or pathways in mosquitoes. If there is no bite, there is no transmission of pathogens; hence W. smithii offers a different approach to investigate control of blood-feeding vectors of human diseases.

The spread of blood-borne pathogens by mosquitoes relies on their taking a blood meal; if there is no bite, there is no disease transmission. Although many species of mosquitoes never take a blood meal, identifying genes that distinguish blood feeding from obligate nonbiting is hampered by the fact that these different lifestyles occur in separate, genetically incompatible species. There is, however, one unique extant species with populations that share a common genetic background but blood feed in one region and are obligate nonbiters in the rest of their range: Wyeomyia smithii. Contemporary blood-feeding and obligate nonbiting populations represent end points of divergence between fully interfertile southern and northern populations. This divergence has undoubtedly resulted in genetic changes that are unrelated to blood feeding, and the challenge is to winnow out the unrelated genetic factors to identify those related specifically to the evolutionary transition from blood feeding to obligate nonbiting. Herein, we determine differential gene expression resulting from directional selection on blood feeding within a polymorphic population to isolate genetic differences between blood feeding and obligate nonbiting. We show that the evolution of nonbiting has resulted in a greatly reduced metabolic investment compared with biting populations, a greater reliance on opportunistic metabolic pathways, and greater reliance on visual rather than olfactory sensory input. W. smithii provides a unique starting point to determine if there are universal nonbiting genes in mosquitoes that could be manipulated as a means to control vector-borne disease.

Natural Variation and Genetics of Photoperiodism in Wyeomyia smithii

Seasonal change in the temperate and polar regions of Earth determines how the world looks around us and, in fact, how we live our day-to-day lives. For biological organisms, seasonal change typically involves complex physiological and metabolic reorganization, the majority of which is regulated by photoperiodism. Photoperiodism is the ability of animals and plants to use day length or night length, resulting in life-historical transformations, including seasonal development, migration, reproduction, and dormancy. Seasonal timing determines not only survival and reproductive success but also the structure and organization of complex communities and, ultimately, the biomes of Earth. Herein, a small mosquito, Wyeomyia smithii, that lives only in the water-filled leaves of a carnivorous plant over a wide geographic range, is used to explore the genetic and evolutionary basis of photoperiodism. Photoperiodism in W. smithii is considered in the context of its historical biogeography in nature to examine the startling finding that recent rapid climate change can drive genetic change in plants and animals at break-neck speed, and to challenge the ponderous 80+ year search for connections between daily and seasonal time-keeping mechanisms. Finally, a model is proposed that reconciles the seemingly disparate 24-h daily clock driven by the invariant rotation of Earth about its axis with the evolutionarily flexible seasonal timer orchestrated by variable seasonality driven by the rotation of Earth about the Sun.

Female preference functions drive interpopulation divergence in male signalling: call diversity in the bushcricket Ephippiger diurnus

Female preferences play a major role in the elaboration and diversification of male traits: as a selective pressure on males, variation in female preferences can generate population divergence and ultimately, speciation. We studied how interpopulation differences in the shape of female mate preference functions may have shaped male advertisement signals in the bushcricket Ephippiger diurnus. This species is distributed as geographically isolated populations with striking interpopulation variation in male acoustic signals, most notably in the number of syllables per call. Here, we asked whether differences in the shape of preference functions exist among populations and whether those differences may have driven male signal evolution resulting in the observed differences in syllable numbers. Our results reveal fundamental differences in female preferences among populations, with differences in the overall preference function shape corresponding to differences in male signals. These differences in female preferences best explain the major differences in male signals among populations. The interpopulation variation in signals and preferences potentially reflects the evolutionary history of the species and may contribute to further divergence among populations and subsequent speciation.

Evolutionary Divergence of Circadian and Photoperiodic Phenotypes in the Pitcher-Plant Mosquito, Wyeomyia smithii

For decades, chronobiologists have investigated the relationship between the circadian clock that mediates daily activities and the photoperiodic timer that mediates seasonal activities. The main experiment used to infer a circadian basis for photoperiodic time measurement is the Nanda-Hamner protocol (NH). Herein, the authors compare additive and nonadditive (dominance and epistasis) genetic effects that lead to the divergence of populations of the pitcher-plant mosquito, Wyeomyia smithii, for critical photoperiod (CPP) and amplitude of the rhythmic response to NH for 3 temporal-geographic scales: 1) Over geological time between populations in northern and southern clades, 2) over millennial time between populations within the northern clade, and 3) over generational time between lines selected for long and short CPP from within a single population. The authors show that the pattern of additive, dominance, and epistatic effects depends on the time scale over which populations or lines have diverged. Patterns for genetic differences between populations for CPP and response to NH reveal similarities over geological and millennial time scales but differences over shorter periods of evolution. These results, and the observation that neither the period nor amplitude of the NH rhythm are significantly correlated with CPP among populations, lead the authors to conclude that the rhythmic response to NH has evolved independently of photoperiodic response in populations of W. smithii. The implication is that in this species, genetic modification of the circadian clock has not been the basis for the adaptive modification of photoperiodic time measurement over the climatic gradient of North America.

Evolutionary divergence of core and post-translational circadian clock genes in the pitcher-plant mosquito, Wyeomyia smithii

BACKGROUND

Internal circadian (circa, about; dies, day) clocks enable organisms to maintain adaptive timing of their daily behavioral activities and physiological functions. Eukaryotic clocks consist of core transcription-translation feedback loops that generate a cycle and post-translational modifiers that maintain that cycle at about 24 h. We use the pitcher-plant mosquito, Wyeomyia smithii (subfamily Culicini, tribe Sabethini), to test whether evolutionary divergence of the circadian clock genes in this species, relative to other insects, has involved primarily genes in the core feedback loops or the post-translational modifiers. Heretofore, there is no reference transcriptome or genome sequence for any mosquito in the tribe Sabethini, which includes over 375 mainly circumtropical species.

METHODS

We sequenced, assembled and annotated the transcriptome of W. smithii containing nearly 95 % of conserved single-copy orthologs in animal genomes. We used the translated contigs and singletons to determine the average rates of circadian clock-gene divergence in W. smithii relative to three other mosquito genera, to Drosophila, to the butterfly, Danaus, and to the wasp, Nasonia.

RESULTS

Over 1.08 million cDNA sequence reads were obtained consisting of 432.5 million nucleotides. Their assembly produced 25,904 contigs and 54,418 singletons of which 62 % and 28 % are annotated as protein-coding genes, respectively, sharing homology with other animal proteomes.

DISCUSSION

The W. smithii transcriptome includes all nine circadian transcription-translation feedback-loop genes and all eight post-translational modifier genes we sought to identify (Fig. 1). After aligning translated W. smithii contigs and singletons from this transcriptome with other insects, we determined that there was no significant difference in the average divergence of W. smithii from the six other taxa between the core feedback-loop genes and post-translational modifiers.

CONCLUSIONS

The characterized transcriptome is sufficiently complete and of sufficient quality to have uncovered all of the insect circadian clock genes we sought to identify (Fig. 1). Relative divergence does not differ between core feedback-loop genes and post-translational modifiers of those genes in a Sabethine species (W. smithii) that has experienced a continual northward dispersal into temperate regions of progressively longer summer day lengths as compared with six other insect taxa. An associated microarray platform derived from this work will enable the investigation of functional genomics of circadian rhythmicity, photoperiodic time measurement, and diapause along a photic and seasonal geographic gradient.

Spatial genetic structure patterns of phenotype-limited and boundary-limited expanding populations: a simulation study

Range expansions may create a unique spatial genetic pattern characterized by alternate genetically homogeneous domains and allele frequency clines. Previous attempts to model range expansions have mainly focused on the loss of genetic diversity during expansions. Using individual-based models, we examined spatial genetic patterns under two expansion scenarios, boundary-limited range expansions (BLRE) and phenotype-limited range expansions (PhLRE). Our simulation revealed that the genetic diversity within populations lost quickly during the range expansion, while the genetic difference accumulated between populations. Consequently, accompanying the expansions, the overall diversity featured a slow decrease. Specifically, during BLREs, high speed of boundary motion facilitated the maintenance of total genetic diversity and sharpened genetic clines. Very slight constraints on boundary motion of BLREs drastically narrowed the homogeneous domains and increased the allele frequency fluctuations from those levels exhibited by PhLREs. Even stronger constraints, however, surprisingly brought the width of homogeneous domains and the allele frequency fluctuations back to the normal levels of PhLREs. Furthermore, high migration rates maintained a higher total genetic diversity than low ones did during PhLREs. Whereas, the total genetic diversities during BLREs showed a contrary pattern: higher when migration was low than those when migration was high. Besides, the increase of migration rates helped maintain a greater number of homogeneous domains during PhLREs, but their effects on the number of homogeneous domains during BLREs were not monotonous. Previous studies have showed that the homogenous domains can merge to form a few broad domains as the expansion went on, leading to fewer homogeneous domains. Our simulations, meanwhile, revealed that the range expansions could also rebuild homogeneous domains from the clines during the range expansion. It is possible that that the number of homogeneous domains was determined by the interaction of merging and newly emerging homogeneous domains.

Replicate phylogenies and post-glacial range expansion of the pitcher-plant mosquito, Wyeomyia smithii, in North America

Herein we tested the repeatability of phylogenetic inference based on high throughput sequencing by increased taxon sampling using our previously published techniques in the pitcher-plant mosquito, Wyeomyia smithii in North America. We sampled 25 natural populations drawn from different localities nearby 21 previous collection localities and used these new data to construct a second, independent phylogeny, expressly to test the reproducibility of phylogenetic patterns. Comparison of trees between the two data sets based on both maximum parsimony and maximum likelihood with Bayesian posterior probabilities showed close correspondence in the grouping of the most southern populations into clear clades. However, discrepancies emerged, particularly in the middle of W. smithii's current range near the previous maximum extent of the Laurentide Ice Sheet, especially concerning the most recent common ancestor to mountain and northern populations. Combining all 46 populations from both studies into a single maximum parsimony tree and taking into account the post-glacial historical biogeography of associated flora provided an improved picture of W. smithii's range expansion in North America. In a more general sense, we propose that extensive taxon sampling, especially in areas of known geological disruption is key to a comprehensive approach to phylogenetics that leads to biologically meaningful phylogenetic inference.

Footprints in time: comparative quantitative trait loci mapping of the pitcher-plant mosquito, Wyeomyia smithii

Identifying regions of the genome contributing to phenotypic evolution often involves genetic mapping of quantitative traits. The focus then turns to identifying regions of 'major' effect, overlooking the observation that traits of ecological or evolutionary relevance usually involve many genes whose individual effects are small but whose cumulative effect is large. Herein, we use the power of fully interfertile natural populations of a single species of mosquito to develop three quantitative trait loci (QTL) maps: one between two post-glacially diverged populations and two between a more ancient and a post-glacial population. All demonstrate that photoperiodic response is genetically a highly complex trait. Furthermore, we show that marker regressions identify apparently 'non-significant' regions of the genome not identified by composite interval mapping, that the perception of the genetic basis of adaptive evolution is crucially dependent upon genetic background and that the genetic basis for adaptive evolution of photoperiodic response is highly variable within contemporary populations as well as between anciently diverged populations.

Island biology and morphological divergence of the Skyros wall lizard Podarcis gaigeae: a combined role for local selection and genetic drift on color morph frequency divergence?

BACKGROUND

Patterns of spatial variation in discrete phenotypic traits can be used to draw inferences about the adaptive significance of traits and evolutionary processes, especially when compared to patterns of neutral genetic variation. Population divergence in adaptive traits such as color morphs can be influenced by both local ecology and stochastic factors such as genetic drift or founder events. Here, we use quantitative color measurements of males and females of Skyros wall lizard, Podarcis gaigeae, to demonstrate that this species is polymorphic with respect to throat color, and the morphs form discrete phenotypic clusters with limited overlap between categories. We use divergence in throat color morph frequencies and compare that to neutral genetic variation to infer the evolutionary processes acting on islet- and mainland populations.

RESULTS

Geographically close islet- and mainland populations of the Skyros wall lizard exhibit strong divergence in throat color morph frequencies. Population variation in throat color morph frequencies between islets was higher than that between mainland populations, and the effective population sizes on the islets were small (Ne:s < 100). Population divergence (FST) for throat color morph frequencies fell within the neutral FST-distribution estimated from microsatellite markers, and genetic drift could thus not be rejected as an explanation for the pattern. Moreover, for both comparisons among mainland-mainland population pairs and between mainland-islet population pairs, morph frequency divergence was significantly correlated with neutral divergence, further pointing to some role for genetic drift in divergence also at the phenotypic level of throat color morphs.

CONCLUSIONS

Genetic drift could not be rejected as an explanation for the pattern of population divergence in morph frequencies. In spite of an expected stabilising selection, throat color frequencies diverged in the islet populations. These results suggest that there is an interaction between selection and genetic drift causing divergence even at a phenotypic level in these small, subdivided populations.

Resolving postglacial phylogeography using high-throughput sequencing

The distinction between model and nonmodel organisms is becoming increasingly blurred. High-throughput, second-generation sequencing approaches are being applied to organisms based on their interesting ecological, physiological, developmental, or evolutionary properties and not on the depth of genetic information available for them. Here, we illustrate this point using a low-cost, efficient technique to determine the fine-scale phylogenetic relationships among recently diverged populations in a species. This application of restriction site-associated DNA tags (RAD tags) reveals previously unresolved genetic structure and direction of evolution in the pitcher plant mosquito, Wyeomyia smithii, from a southern Appalachian Mountain refugium following recession of the Laurentide Ice Sheet at 22,000-19,000 B.P. The RAD tag method can be used to identify detailed patterns of phylogeography in any organism regardless of existing genomic data, and, more broadly, to identify incipient speciation and genome-wide variation in natural populations in general.

Evolution of photoperiodic time measurement is independent of the circadian clock in the pitcher-plant mosquito, Wyeomyia smithii

For over 70 years, researchers have debated whether the ability to use day length as a cue for the timing of seasonal events (photoperiodism) is related to the endogenous circadian clock that regulates the timing of daily events. Models of photoperiodism include two components: (1) a photoperiodic timer that measures the length of the day, and (2) a photoperiodic counter that elicits the downstream photoperiodic response after a threshold number of days has been counted. Herein, we show that there is no geographical pattern of genetic association between the expression of the circadian clock and the photoperiodic timer or counter. We conclude that the photoperiodic timer and counter have evolved independently of the circadian clock in the pitcher-plant mosquito Wyeomyia smithii and hence, the evolutionary modification of photoperiodism throughout the range of W. smithii has not been causally mediated by a corresponding evolution of the circadian clock.

Quantitative trait loci associated with photoperiodic response and stage of diapause in the pitcher-plant mosquito, Wyeomyia smithii

A wide variety of temperate animals rely on length of day (photoperiodism) to anticipate and prepare for changing seasons by regulating the timing of development, reproduction, dormancy, and migration. Although the molecular basis of circadian rhythms regulating daily activities is well defined, the molecular basis for the photoperiodic regulation of seasonal activities is largely unknown. We use geographic variation in the photoperiodic control of diapause in the pitcher-plant mosquito Wyeomyia smithii to create the first QTL map of photoperiodism in any animal. For critical photoperiod (CPP), we detect QTL that are unique, a QTL that is sex linked, QTL that overlap with QTL for stage of diapause (SOD), and a QTL that interacts epistatically with the circadian rhythm gene, timeless. Results presented here confirm earlier studies concluding that CPP is under directional selection over the climatic gradient of North America and that the evolution of CPP is genetically correlated with SOD. Despite epistasis between timeless and a QTL for CPP, timeless is not located within any detectable QTL, indicating that it plays an ancillary role in the evolution of photoperiodism in W. smithii. Finally, we highlight one region of the genome that includes loci contributing to CPP, SOD, and hormonal regulation of development.

An unusual phylogeography in the bushcricket Ephippiger ephippiger from Southern France

Pleistocene glaciations have played a major role in species divergence. The bushcricket Ephippiger ephippiger shows unusual patterns of intraspecific variation in multiple traits across Southern Europe. This is centred in Southern France, and evidence implies that it results from secondary contact after differentiation in Pleistocene refugia. However, the possible time scales involved, locations of the refugia and patterns of expansion remain obscure. This study sequenced the COII (507 BP) and cyt b (428 BP) mitochondrial genes to examine the intraspecific phylogeography of Western European samples of E. ephippiger. A minimum evolution tree revealed little resolution between described subspecies of E. ephippiger. Strikingly, populations from the Pyrenees and Mediterranean coastal region contained a complex genetic structure corresponding to major river valleys, independent of the traditional taxonomy. Samples of the subspecies E. e. vitium formed a distinct clade, perhaps supporting their taxonomic status. However, other forms (cruciger and cunii) were not genetically distinct, which is surprising given differences in their morphology and behaviour. The extent of the genetic divergence between Pyreneen valleys is unexpectedly deep, with average Tamura-Nei distances of around 14% (net distances of 11%) separating the main clades of coding COII sequences. Cyt b showed a similar pattern, but was confounded by some non-coding probable pseudogenes. If a conventional insect molecular clock is applied, these cryptic clades must pre-date the Pleistocene, and hypotheses for their history are discussed. However, mtDNA divergence in Ephippiger is not evolving in a clock-like manner, because a likelihood ratio test rejects clock assumptions for the COII sequences.

Phylogeography of postglacial range expansion in Nigronia serricornis Say (Megaloptera: Corydalidae)

We have examined the effects of post-Wisconsinan glacial range expansion on the phylogeography of the saw-combed fishfly, Nigronia serricornis Say (Megaloptera: Corydalidae) because aquatic insects are under-represented in postglacial studies (and in phylogeography in general), and because the effects of ecological degradation on the population genetics of environmental indicator species like N. serricornis cannot be measured unless the underlying phylogeography is understood. Sequence data from a 630-base fragment of the mitochondrial cytochrome oxidase I (COI) gene were subjected to amova and nested clade analysis for 30 populations (n = 344) of N. serricornis. Both the amova and nested clade analysis revealed substantial population structure; 44.4% of the variance occured among populations. Three northward migrations are apparent: one from Tennessee into Illinois, Indiana, Wisconsin, Michigan and Ontario, a second that radiated eastward from Pennsylvania, and a third that moved along the coast from North Carolina into Connecticut, Massachusetts, Maine and then into New York. The latter two of these migrations were the result of contiguous range expansions, while the former expansion, out of Tennessee, appears to have been rapid with little gene flow from the source population. Additional clades included a group of haplotypes in central Kentucky that appear to have expanded along preglacial drainages, and clades in North Carolina and Georgia that have remained centrally located. Haplotype diversity decreased from south to north, a pattern that has been widely reported for animal and plant populations that expanded with the retreat of the Wisconsinan glaciation.

The conversion of variance and the evolutionary potential of restricted recombination

Genetic recombination is usually considered to facilitate adaptive evolution. However, recombination prevents the reliable cotransmission of interacting gene combinations and can disrupt complexes of coadapted genes. If interactions between genes have important fitness effects, restricted recombination may lead to evolutionary responses that are different from those predicted from a purely additive model and could even aid adaptation. Theory and data have demonstrated that phenomena that limit the effectiveness of recombination via increasing homozygosity, such as inbreeding and population subdivision and bottlenecks, can temporarily increase the additive genetic variance available to these populations. This effect has been attributed to the conversion of nonadditive to additive genetic variance. Analogously, phenomena such as chromosomal inversions and apomictic parthenogenesis that physically restrict recombination in part or all of the genome may also result in a release of additive variance. Here, we review and synthesize literature concerning the evolutionary potential of populations with effectively or physically restricted recombination. Our goal is to emphasize the common theme of increased short-term access to additive genetic variance in all of these situations and to motivate research directed towards a more complete characterization of the relevance of the conversion of variance to the evolutionary process.

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.

Epistasis underlying a fitness trait within a natural population of the pitcher-plant mosquito, Wyeomyia smithii

We selected on divergent photoperiodic response in three separate lines from a natural population of the pitcher-plant mosquito, Wyeomyia smithii. Line crosses reveal that there exists within a population, diverse epistatic variation for a fitness trait that could contribute to adaptive potential following founder events or rapid climate change.

The population genetics of a biological control introduction: mitochondrial DNA and microsatellie variation in native and introduced populations of Aphidus ervi, a parisitoid wasp

Introductions of biological control agents may cause bottlenecks in population size despite efforts to avoid them. We examined the population genetics of Aphidius ervi (Hymenoptera: Braconidae), a parasitoid that was introduced to North America from Western Europe in 1959 to control pea aphids. To explore the phylogeographical relationships of A. ervi we sequenced 1249 bp of mitochondrial DNA (mtDNA) from 27 individuals from the native range and 51 individuals from the introduced range. Most individuals from Western Europe, the Middle East and North America shared one of two common haplotypes, consistent with the known history of the introduction. However, some A. ervi from the Pacific Northwest have a haplotype that is most similar to haplotypes found in Japan, raising the possibility of a second accidental introduction. To examine population structure and assess whether a bottleneck occurred upon introduction to North America, we assayed variation at 5 microsatellite loci in 62 individuals from 2 native populations and 230 individuals from 6 introduced populations. Introduced samples had fewer rare alleles than native samples (F1,34 = 13.5, P = 0.0008), but heterozygosity did not differ significantly. These results suggest that a mild bottleneck occurred in spite of the introduction of over 1000 individuals. Using a hierarchical Bayesian approach, the founding population size was estimated to be 245 individuals. amova showed significant genetic differentiation between the European and North American samples, and a Bayesian assignment approach clustered individuals into four groups, with most European individuals in one group and most North American individuals in the other three. These results highlight that genetic changes are associated with founder events in rapidly growing natural populations, even when the founding population size is relatively large.

Out of Anatolia: longitudinal gradients in genetic diversity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae

Many studies have addressed the latitudinal gradients in intraspecific genetic diversity of European taxa generated during postglacial range expansion from southern refugia. Although Asia Minor is known to be a centre of diversity for many taxa, relatively few studies have considered its potential role as a Pleistocene refugium or a potential source for more ancient westward range expansion into Europe. Here we address these issues for an oak gallwasp, Andricus quercustozae (Hymenoptera: Cynipidae), whose distribution extends from Morocco along the northern coast of the Mediterranean through Turkey to Iran. We use sequence data for a fragment of the mitochondrial gene cytochrome b and allele frequency data for 12 polymorphic allozyme loci to answer the following questions: (1). which regions represent current centres of genetic diversity for A. quercustozae? Do eastern populations represent one refuge or several discrete glacial refugia? (2). Can we infer the timescale and sequence of the colonization processes linking current centres of diversity? Our results suggest that A. quercustozae was present in five distinct refugia (Iberia, Italy, the Balkans, southwestern Turkey and northeastern Turkey) with recent genetic exchange between Italy and Hungary. Genetic diversity is greatest in the Turkish refugia, suggesting that European populations are either (a). derived from Asia Minor, or (b). subject to more frequent population bottlenecks. Although Iberian populations show the lowest diversity for putatively selectively neutral markers, they have colonized a new oak host and represent a genetically and biologically discrete entity within the species.

The genetic legacy of the Quaternary ice ages

Global climate has fluctuated greatly during the past three million years, leading to the recent major ice ages. An inescapable consequence for most living organisms is great changes in their distribution, which are expressed differently in boreal, temperate and tropical zones. Such range changes can be expected to have genetic consequences, and the advent of DNA technology provides most suitable markers to examine these. Several good data sets are now available, which provide tests of expectations, insights into species colonization and unexpected genetic subdivision and mixture of species. The genetic structure of human populations may be viewed in the same context. The present genetic structure of populations, species and communities has been mainly formed by Quaternary ice ages, and genetic, fossil and physical data combined can greatly help our understanding of how organisms were so affected.