Montane Grasshoppers Research Articles

High-Elevation Populations of Montane Grasshoppers Exhibit Greater Developmental Plasticity in Response to Seasonal Cues.

Populations of insects can differ in how sensitive their development, growth, and performance are to environmental conditions such as temperature and daylength. The environmental sensitivity of development can alter phenology (seasonal timing) and ecology. Warming accelerates development of most populations. However, high-elevation and season-limited populations can exhibit developmental plasticity to either advance or prolong development depending on conditions. We examine how diurnal temperature variation and daylength interact to shape growth, development, and performance of several populations of the montane grasshopper, Melanoplus boulderensis, along an elevation gradient. We then compare these experimental results to observed patterns of development in the field. Although populations exhibited similar thermal sensitivities of development under long-day conditions, development of high-elevation populations was more sensitive to temperature under short-day conditions. This developmental plasticity resulted in rapid development of high elevation populations in short-day conditions with high temperature variability, consistent with their observed capacity for rapid development in the field when conditions are permissive early in the season. Notably, accelerated development generally did not decrease body size or alter body shape. Developmental conditions did not strongly influence thermal tolerance but altered the temperature dependence of performance in difficult-to-predict ways. In sum, the high-elevation and season-limited populations exhibited developmental plasticity that enables advancing or prolonging development consistent with field phenology. Our results suggest these patterns are driven by the thermal sensitivity of development increasing when days are short early in the season compared to when days are long later in the season. Developmental plasticity will shape phenological responses to climate change with potential implications for community and ecosystem structure.

Glacial connectivity and current population fragmentation in sky islands explain the contemporary distribution of genomic variation in two narrow‐endemic montane grasshoppers from a biodiversity hotspot

AbstractAimCold‐adapted biotas from mid‐latitudes often show small population sizes, harbour low levels of local genetic diversity and are highly vulnerable to extinction due to ongoing climate warming and the progressive shrinking of montane and alpine ecosystems. In this study, we use a suite of analytical approaches to infer the demographic processes that have shaped contemporary patterns of genomic variation in Omocestus bolivari and Omocestus femoralis, two narrow‐endemic and red‐listed Iberian grasshoppers forming highly fragmented populations in the sky island archipelago of the Baetic System.LocationSouth‐eastern Iberia.MethodsWe quantified genomic variation in the two focal taxa and coupled ecological niche models and a spatiotemporally explicit simulation approach based on coalescent theory to determine the relative statistical support of a suite of competing demographic scenarios representing contemporary population isolation (i.e. a predominant role of genetic drift) versus historical connectivity and post‐glacial colonization of sky islands (i.e. pulses of gene flow and genetic drift linked to Pleistocene glacial cycles).ResultsInference of spatial patterns of genetic structure, environmental niche modelling and statistical evaluation of alternative species‐specific demographic models within an approximate Bayesian computation framework collectively supported genetic admixture during glacial periods and post‐glacial colonization of sky islands, rather than long‐term population isolation, as the scenario best explaining the current distribution of genomic variation in the two focal taxa. Moreover, our analyses revealed that isolation in sky islands has also led to extraordinary genetic fragmentation and contributed to reduce local levels of genetic diversity.Main conclusionsThis study exemplifies the potential of integrating genomic and environmental niche modelling data across biological and spatial replicates to determine whether organisms with similar habitat requirements have experienced concerted/idiosyncratic responses to Quaternary climatic oscillations, which can ultimately help to reach more general conclusions about the vulnerability of mountain biodiversity hotspots to ongoing climate warming.

Demographic consequences of dispersal-related trait shift in two recently diverged taxa of montane grasshoppers.

Although the pervasiveness of intraspecific wing-size polymorphism and transitions to flightlessness have long captivated biologists, the demographic outcomes of shifts in dispersal ability are not yet well understood and have been seldom studied at early stages of diversification. Here, we use genomic data to infer the consequences of dispersal-related trait variation in the taxonomically controversial short-winged (Chorthippus corsicus corsicus) and long-winged (Chorthippus corsicus pascuorum) Corsican grasshoppers. Our analyses revealed lack of contemporary hybridization between sympatric long- and short-winged forms and phylogenomic reconstructions supported their taxonomic distinctiveness, rejecting the hypothesis of intraspecific wing polymorphism. Statistical evaluation of alternative models of speciation strongly supported a scenario of Pleistocene divergence (<1.5 Ma) with ancestral gene flow. According to neutral expectations from differences in dispersal capacity, historical effective migration rates from the long- to the short-winged taxon were threefold higher than in the opposite direction. Although populations of the two taxa present a marked genetic structure and have experienced parallel demographic histories, our coalescent-based analyses suggest that reduced dispersal has fueled diversification in the short-winged C. c. corsicus. Collectively, our study illustrates how dispersal reduction can speed up geographical diversification and increase the opportunity for allopatric speciation in topographically complex landscapes.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers.

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200m), and one narrower, K. tristis (1600-2000m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using FST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n=130; K. tristis: 7363 SNPs, n=135). Stronger genetic structure was found in the more narrowly distributed K.tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n=3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.

Grasshopper phenological responses to climate gradients, variability, and change

AbstractSpecies have responded to recent climate change via pronounced, but often variable, shifts in seasonal timing (phenology). Can this variability be accounted for by comparing responses along an elevational gradient and among species differing in seasonal timing and traits? We resurveyed montane grasshopper communities along a 1300‐m elevational gradient to examine phenological responses to ~50 yr of climate warming. Phenological responses were elevation dependent, influenced by relative warming (which was more pronounced at higher elevations) and by substantial interannual variability in available growing degree‐days (GDDs). The total number of GDDs available for development within a season declines with elevation, and correspondingly high‐elevation populations required fewer GDDs to reach adulthood. Although phenological advancements were particularly pronounced at higher elevations, consistent with a greater magnitude of climate warming, species as a whole initiated and completed development earlier in warmer years. During warm years, phenological advances were most pronounced early in development, and the timing of adulthood tended to occur earlier (at a similar or greater accumulation of GDDs). Phenological advancement was also more pronounced in early season species, which exhibit developmental plasticity to capitalize on warming. Interactions with elevation were also strongest for the early season species. Our study shows the importance of incorporating developmental information to understand the environmental and organismal mechanisms that drive variable phenological responses and suggests that the phenological impacts of warming will be concentrated early in development, among early season species, and at high elevations. Considering these drivers of phenological variability can aid understanding and predicting climate‐induced changes to community and ecosystem structure.

Genomic data reveal deep genetic structure but no support for current taxonomic designation in a grasshopper species complex.

Taxonomy has traditionally relied on morphological and ecological traits to interpret and classify biological diversity. Over the last decade, technological advances and conceptual developments in the field of molecular ecology and systematics have eased the generation of genomic data and changed the paradigm of biodiversity analysis. Here we illustrate how traditional taxonomy has led to species designations that are supported neither by high throughput sequencing data nor by the quantitative integration of genomic information with other sources of evidence. Specifically, we focus on Omocestus antigai and Omocestus navasi, two montane grasshoppers from the Pyrenean region that were originally described based on quantitative phenotypic differences and distinct habitat associations (alpine vs. Mediterranean-montane habitats). To validate current taxonomic designations, test species boundaries, and understand the factors that have contributed to genetic divergence, we obtained phenotypic (geometric morphometrics) and genome-wide SNP data (ddRADSeq) from populations covering the entire known distribution of the two taxa. Coalescent-based phylogenetic reconstructions, integrative Bayesian model-based species delimitation, and landscape genetic analyses revealed that populations assigned to the two taxa show a spatial distribution of genetic variation that do not match with current taxonomic designations and is incompatible with ecological/environmental speciation. Our results support little phenotypic variation among populations and a marked genetic structure that is mostly explained by geographic distances and limited population connectivity across the abrupt landscapes characterizing the study region. Overall, this study highlights the importance of integrative approaches to identify taxonomic units and elucidate the evolutionary history of species.

Distributional shifts – not geographic isolation – as a probable driver of montane species divergence

As biodiversity hotspots, montane regions have been a focus of research to understand the divergence process. Like their oceanic counterparts, the diversity of the ‘sky islands’ might be ascribed to geographic isolation of mountaintops. However, because the sky islands, and especially those in northern latitudes, are subject to extreme climatic events such as the glacial cycles that drove both altitudinal and geographical shifts in species’ distributions, the dynamic colonization process is also a possible factor driving divergence. Here we test these two hypotheses (i.e. isolation versus colonization) in a flightless montane grasshopper, Melanoplus oregonensis, which is a member of a diverse group that radiated across the Rocky Mountains of North America. Using approximate Bayesian computation (ABC) and spatially explicit simulations that account for spatial heterogeneity and temporal shifts in species distributions, we show that a colonization model of the sky islands from refugial populations provides a significantly better fit to the empirical genetic data than a model of the geographic isolation among sky islands. Moreover, support for the colonization model holds irrespective of whether the movement of individuals was modeled as a diffusion process or was informed by differences in habitat suitabilities across the landscape. With validation analyses to confirm the models provide a good fit to the data, as well as general power and quality analyses, the research not only adds to a growing body of work on the complex dynamics underlying montane biodiversity, but it also provides much needed evaluation of competing hypotheses based on explicit models of the divergence process, as opposed to inferences about diversification drivers from species diversity patterns.

The role of environment and core-margin effects on range-wide phenotypic variation in a montane grasshopper.

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range-wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal-related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120000years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among-population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal-related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range-wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.

Hierarchical genetic structure shaped by topography in a narrow-endemic montane grasshopper.

BackgroundUnderstanding the underlying processes shaping spatial patterns of genetic structure in free-ranging organisms is a central topic in evolutionary biology. Here, we aim to disentangle the relative importance of neutral (i.e. genetic drift) and local adaptation (i.e. ecological divergence) processes in the evolution of spatial genetic structure of the Morales grasshopper (Chorthippus saulcyi moralesi), a narrow-endemic taxon restricted to the Central Pyrenees. More specifically, we analysed range-wide patterns of genetic structure and tested whether they were shaped by geography (isolation-by-distance, IBD), topographic complexity and present and past habitat suitability models (isolation-by-resistance, IBR), and environmental dissimilarity (isolation-by-environment, IBE).ResultsDifferent clustering analyses revealed a deep genetic structure that was best explained by IBR based on topographic complexity. Our analyses did not reveal a significant role of IBE, a fact that may be due to low environmental variation among populations and/or consequence of other ecological factors not considered in this study are involved in local adaptation processes. IBR scenarios informed by current and past climate distribution models did not show either a significant impact on genetic differentiation after controlling for the effects of topographic complexity, which may indicate that they are not capturing well microhabitat structure in the present or the genetic signal left by dispersal routes defined by habitat corridors in the past.ConclusionsOverall, these results indicate that spatial patterns of genetic variation in our study system are primarily explained by neutral divergence and migration-drift equilibrium due to limited dispersal across abrupt reliefs, whereas environmental variation or spatial heterogeneity in habitat suitability associated with the complex topography of the region had no significant effect on genetic discontinuities after controlling for geography. Our study highlights the importance of considering a comprehensive suite of potential isolating mechanisms and analytical approaches in order to get robust inferences on the processes promoting genetic divergence of natural populations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-016-0663-7) contains supplementary material, which is available to authorized users.

Exploring the population genetic consequences of the colonization process with spatio‐temporally explicit models: insights from coupled ecological, demographic and genetic models in montane grasshoppers

Understanding the genetic consequences of shifting species distributions is critical for evaluating the impact of climate-induced distributional changes. However, the demographic expansion associated with the colonization process typically takes place across a heterogeneous environment, with population sizes and migration rates varying across the landscape. Here we describe an approach for coupling ecological-niche models (ENMs) with demographic and genetic models to explore the genetic consequences of distributional shifts across a heterogeneous landscape. Analyses of a flightless grasshopper from the sky islands of the Rocky Mountains of North America are used to show how biologically informed predictions can be generated about the genetic consequences of a colonization process across a spatially and temporally heterogeneous landscape (i.e. the suitability of habitats for the montane species differs across the landscape and is itself not static, with the displacement of contemporary populations into glacial refugia). By using (i) ENMs for current climatic conditions and the last glacial maximum to (ii) parameterize a demographic model of the colonization process, which then (iii) informs coalescent simulations, a set of models can be generated that capture different processes associated with distributional shifts. We discuss how the proposed approach for model generation can be integrated into a statistical framework for estimating key demographic parameters and testing hypotheses about the conditions for which distributional shifts may (or may not) enhance species divergence, including the importance of habitat stability, past gene-flow among currently isolated populations, and maintenance of refugial populations in multiple geographic regions.

Importance of Ecological Scale in Montane Grasshopper (Orthoptera: Acrididae) Species Structure in Similar Habitat between Differing Soil Textures and Dominant Vegetative Canopy Coverage

Discrepancies in scale dynamics often make cohesive structural conclusions difficult, especially when dealing with ecological variance. We studied presence and abundance of grasshopper species in similar, yet distinct, montane habitat of southeast Wyoming and northern Colorado, USA. By limiting ecological variance, grasshopper species structure at two behavioral scales (grouped species dynamics and individual species interactions) was maximized with regard to soil texture and vegetative canopy coverage. Combining univariate and multivariate statistical methods, we note montane grasshopper species interactions to be scale dependent.

Coupling Genetic and Ecological-Niche Models to Examine How Past Population Distributions Contribute to Divergence

Understanding the impact of climate-induced distributional shifts on species divergence, like those accompanying the Pleistocene glacial cycles [1, 2], requires tools that explicitly incorporate the geographic configuration of past distributions into analyses of genetic differentiation. Depending on the historical distribution of species, genetic differences may accumulate among ancestral source populations, but there is long-standing debate whether displacements into glacial refugia promoted divergence. Here we integrate coalescent-based genetic models [3, 4] with ecological-niche modeling [5, 6] to generate expectations for patterns of genetic variation based on an inferred past distribution of a species. Reconstruction of the distribution of a montane grasshopper species during the last glacial maximum suggests that Melanoplus marshalli populations from the sky islands of Colorado and Utah were likely colonized from multiple ancestral source populations. The genetic analyses provide compelling evidence that the historical distribution of M. marshalli-namely, spatial separation of multiple refugia-was conducive to genetic differentiation. The coupling of genetic and ecological-niche modeling provides a new and flexible tool for integrating paleoenvironmental details into species-specific predictions of population structure that can increase our understanding of why the glacial cycles promoted speciation in some taxa and yet inhibited diversification in others [7, 8].

ESTIMATING A GEOGRAPHICALLY EXPLICIT MODEL OF POPULATION DIVERGENCE

Patterns of genetic variation can provide valuable insights for deciphering the relative roles of different evolutionary processes in species differentiation. However, population-genetic models for studying divergence in geographically structured species are generally lacking. Since these are the biogeographic settings where genetic drift is expected to predominate, not only are population-genetic tests of hypotheses in geographically structured species constrained, but generalizations about the evolutionary processes that promote species divergence may also be potentially biased. Here we estimate a population-divergence model in montane grasshoppers from the sky islands of the Rocky Mountains. Because this region was directly impacted by Pleistocene glaciation, both the displacement into glacial refugia and recolonization of montane habitats may contribute to differentiation. Building on the tradition of using information from the genealogical relationships of alleles to infer the geography of divergence, here the additional consideration of the process of gene-lineage sorting is used to obtain a quantitative estimate of population relationships and historical associations (i.e., a population tree) from the gene trees of five anonymous nuclear loci and one mitochondrial locus in the broadly distributed species Melanoplus oregonensis. Three different approaches are used to estimate a model of population divergence; this comparison allows us to evaluate specific methodological assumptions that influence the estimated history of divergence. A model of population divergence was identified that significantly fits the data better compared to the other approaches, based on per-site likelihood scores of the multiple loci, and that provides clues about how divergence proceeded in M. oregonensis during the dynamic Pleistocene. Unlike the approaches that either considered only the most recent coalescence (i.e., information from a single individual per population) or did not consider the pattern of coalescence in the gene genealogies, the population-divergence model that best fits the data was estimated by considering the pattern of gene lineage coalescence across multiple individuals, as well as loci. These results indicate that sampling of multiple individuals per population is critical to obtaining an accurate estimate of the history of divergence so that the signal of common ancestry can be separated from the confounding influence of gene flow-even though estimates suggest that gene flow is not a predominant factor structuring patterns of genetic variation across these sky island populations. They also suggest that the gene genealogies contain information about population relationships, despite the lack of complete sorting of gene lineages. What emerges from the analyses is a model of population divergence that incorporates both contemporary distributions and historical associations, and shows a latitudinal and regional structuring of populations reminiscent of population displacements into multiple glacial refugia. Because the population-divergence model itself is built upon the specific events shaping the history of M. oregonensis, it provides a framework for estimating additional population-genetic parameters relevant to understanding the processes governing differentiation in geographically structured species and avoids the problems of relying on overly simplified and inaccurate divergence models. The utility of these approaches, as well as the caveats and future improvements, for estimating population relationships and historical associations relevant to genetic analyses of geographically structured species are discussed.

Variable nuclear markers for Melanoplus oregonensis identified from the screening of a genomic library

AbstractWe report the isolation of seven single‐copy nuclear loci from the montane grasshopper Melanoplus oregonensis. With an average length of approximately 930 bp and between five and 67 variable sites, these loci are useful for population genetic analyses within M. oregonensis. Amplification of nuclear loci in additional species of western Melanoplus grasshoppers suggests that they will also be useful for a variety of population genetic, phylogeographic and phylogenetic studies.

Importance of genetic drift during Pleistocene divergence as revealed by analyses of genomic variation

Determining what factors affect the structuring of genetic variation is key to deciphering the relative roles of different evolutionary processes in species differentiation. Such information is especially critical to understanding how the frequent shifts and fragmentation of species distributions during the Pleistocene translates into species differences, and why the effect of such rapid climate change on patterns of species diversity varies among taxa. Studies of mitochondrial DNA (mtDNA) have detected significant population structure in many species, including those directly impacted by the glacial cycles. Yet, understanding the ultimate consequence of such structure, as it relates to how species divergence occurs, requires demonstration that such patterns are also shared with genomic patterns of differentiation. Here we present analyses of amplified fragment length polymorphisms (AFLPs) in the montane grasshopper Melanoplus oregonensis to assess the evolutionary significance of past demographic events and associated drift-induced divergence as inferred from mtDNA. As an inhabitant of the sky islands of the northern Rocky Mountains, this species was subject to repeated and frequent shifts in species distribution in response to the many glacial cycles. Nevertheless, significant genetic structuring of M. oregonensis is evident at two different geographic and temporal scales: recent divergence associated with the recolonization of the montane meadows in individual sky islands, as well as older divergence associated with displacements into regional glacial refugia. The genomic analyses indicate that drift-induced divergence, despite the lack of long-standing geographic barriers, has significantly contributed to species divergence during the Pleistocene. Moreover, the finding that divergence associated with past demographic events involves the repartitioning of ancestral variation without significant reductions of genomic diversity has intriguing implications - namely, the further amplification of drift-induced divergence by selection.

Association of genetic lineages with ecological features in a polyphagous montane grasshopper species

Abstract Understanding the relationships between intraspecific diversity and habitat has been greatly enhanced with the advent of molecular markers. Such understanding may help to elucidate complex microevolutionary processes. Melanoplus alpinus, a montane/alpine grasshopper species found in the central and northern Rocky Mountains, has a disjunct distribution with highly divergent mtDNA lineages, structured among meadows and drainages within mountain ranges. Previous analyses showed that genetic differentiation is not distance driven, so habitat factors may account for genetic structuring and diversity within the species. In this study, 3 mtDNA lineages were analyzed for their association with ecological variables, and PCR-RFLP haplotype genetic diversities were analyzed for relationships with elevation, elevation-latitude index, and meadow size. Chi-square analysis revealed mtDNA lineage predominance in particular soil textures, indicating the potential of soil to function as a limiting environmental fa...

Tergoceracris, a new genus and six new species of montane grasshoppers (Orthoptera: Acrididae: Ommatolampinae) from Dominican Republic and Puerto Rico

Tergoceracris, a new genus and six new species of montane grasshoppers (Orthoptera: Acrididae: Ommatolampinae) from Dominican Republic and Puerto Rico

Tergoceracris, a new genus and six new species of montane grasshoppers (Orthoptera: Acrididae: Ommatolampinae) from Dominican Republic and Puerto Rico

The new genus Tergoceracris consists of six new species; all are brachypterous, semi-arboreal and montane grasshoppers. They are distinguished from all other known neotropical acridoids by having exaggerated and vertically projecting furculae in the males. Two new species are from the Dominican Republic (Tergoceracris ocampensis n. sp. and Tergoceracris ebanoverde n. sp.) and four from Puerto Rico (Tergoceracris guajataca n. sp., Tergoceracris cerropunta n. sp., Tergoceracris cayey n. sp., and Tergoceracris luquillensis n. sp.).

Did the pleistocene glaciations promote divergence? Tests of explicit refugial models in montane grasshopprers.

There is a long-standing debate over whether or not the Pleistocene glaciations promoted speciation. While some models predict that extensive mixing of populations during interglacial expansion would have inhibited divergence, others postulate that divergence among allopatric glacial refuges or founder events during recolonization of previously glaciated areas would have promoted differentiation. Using a combination of traditional and coalescent based population genetic approaches, this study finds that the glaciations did not inhibit divergence among populations of the grasshopper Melanoplus oregonensis. Instead, drift associated with recolonization of previously glaciated areas, as well as divergence among multiple allopatric glacial refugia, have both contributed to differentiation in this montane grasshopper from the 'sky islands' of the northern Rocky Mountains. Significant population structure was detected by phylogenetic and FST analyses, including significant FST values among individual pairs of sky-island populations. In addition to clustering of haplotypes within populations, there is some evidence of regional phylogeographic structure, although none of the 'regional groups' form a monophyletic clade and there is a lack of concordance between the genealogical and geographical positions of some haplotypes. However, coalescent simulations confirm there is significant regional phylogeographic structure that most likely reflects divergence among multiple ancestral refugial populations, and indicate that it is very unlikely that the observed gene tree could have been produced by the fragmentation of a single widespread ancestral population. Thus, rather than inhibiting differentiation, the glaciations appear to have promoted population divergence in M. oregonensis, suggesting that they may have contributed to the radiation of Melanoplus species during the Pleistocene.