Mitochondrial DNA variation and GIS analysis confirm a secondary origin of geographical variation in the bushcricket Ephippiger ephippiger (Orthoptera: Tettigonioidea), and resurrect two subspecies

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'.

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.

Association between the geographic distribution during the last glacial maximum of Asian wild rice, Oryza rufipogon (Poaceae), and its current genetic variation

PREMISE OF THE STUDY

The combination of traditional population genetic studies and species distribution modeling (SDM) provides many new insights in detecting phylogeographic signals. In Asian wild rice (Oryza rufipogon), the progenitor of cultivated Asian rice, geographical subdivision has been documented in many genetic studies although the root cause of this subdivision remains unknown. Surprisingly, environmental factors associated with the spatial and temporal distribution of O. rufipogon have rarely been examined. The aim of this study is to understand the historical distribution of O. rufipogon and its relationship to the current geographical pattern of genetic variation.

METHODS

We used SDM to examine the present, past, and future distribution of O. rufipogon. The estimated distribution during the Last Glacial Maximum was then compared with genetic data from our previous work.

KEY RESULTS

The predicted paleodistribution of O. rufipogon at the Last Glacial Maximum was separated into disconnected east and west ranges. This past distribution is consistent with the current geographic pattern of genetic variation, with two genetic groups that intergrade. Annual precipitation is the single factor that contributes most to SDM estimates. SDM predictions for 2080 indicate a general trend of increasing probability of presence and range expansion.

CONCLUSIONS

(1) The historically disjunct distribution potentially contributes to the current genetic subdivision of O. rufipogon. (2) Water availability is an important factor that limits the distribution of O. rufipogon. (3) Global warming is a lesser threat than other human-mediated factors to the conservation of this endangered species.

Integrating GIS-based environmental data into evolutionary biology

Many evolutionary processes are influenced by environmental variation over space and time, including genetic divergence among populations, speciation and evolutionary change in morphology, physiology and behaviour. Yet, evolutionary biologists have generally not taken advantage of the extensive environmental data available from geographic information systems (GIS). For example, studies of phylogeography, speciation and character evolution often ignore or use only crude proxies for environmental variation (e.g. latitude and distance between populations). Here, we describe how the integration of GIS-based environmental data, along with new spatial tools, can transform evolutionary studies and reveal new insights into the ecological causes of evolutionary patterns.

The past and future influence of geographic information systems on hybrid zone, phylogeographic and speciation research

Over the past two decades geographers have developed an increasingly sophisticated technology termed a geographic information system (GIS). A GIS has the ability to store, map and analyse spatial data. The powerful analytical capabilities of a GIS could serve to enhance our understanding of the spatial component of the evolutionary process. In particular, phylogeographers, hybrid zone and speciation researchers could benefit enormously from incorporating this sophisticated technology from the discipline of geography, as they have done so readily from other disciplines (e.g. genetics). Indeed, an increasing number of researchers in these fields are beginning to include GIS analyses into their research programmes. Some of this integration has taken the form of analysing the spatial relationship between populations and hybrid zones. Several other researchers have also begun to incorporate GIS into their work through the use of GIS-based niche models. These models estimate a multidimensional niche for a species using known geo-referenced populations and digital climate maps. Here, I review the recent integration of GIS and GIS-based predictive niche models into the above evolutionary sub-disciplines. I also describe evolutionary analyses that could be further enhanced through the implementation of GIS.

CONTRASTING PATTERNS OF HYBRIDIZATION IN LARGE HOUSE SPIDERS (TEGENARIA ATRICA GROUP, AGELENIDAE)

The integrity of species is not fixed and may vary geographically. Here we investigate the geographic distributions and interactions of species in the Tegenaria atrica group (Araneae: Agelenidae). Detailed mapping of T. saeva and T. gigantea in England and Wales shows them to be broadly allopatric in southern England with a tightly defined, and possibly long-standing, narrow zone of parapatry in central southern England. In the north of England (Yorkshire), by contrast, the species are broadly sympatric as a result of recent range expansions. GIS techniques are used to map the species distributions and to quantify, we believe for the first time, the intimacy of interspecific interactions. The extent and nature of hybridization in these two areas is examined through regression and multivariate analyses of morphology. We show that the relative incidence of hybridization is much greater in Yorkshire than within the parapatric zone in the south. Clear patterns of asymmetric introgression are observed in both northern and southern England, with a greater impact of T. gigantea on T. saeva than vice versa. We find no sign of morphological reproductive character displacement at the zone of parapatry that might indicate reinforcement, although we cannot exclude more subtle effects, for example via cuticular pheromones. The integrity of these two species seems to be breaking down in northern England, a process that might gain momentum as the gene pools become more similar.

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.

Morphological and genetic divergence of intralacustrine stickleback morphs in Iceland: a case for selective differentiation?

The evolutionary processes involved in population divergence and local adaptation are poorly understood. Theory predicts that divergence of adjacent populations is possible but depends on several factors including gene flow, divergent selection, population size and the number of genes involved in divergence and their distribution on the genome. We analyse variation in neutral markers, markers linked to putative quantitative trait loci and morphological traits in a recent (<10000 years) zone of primary divergence between stickleback morphs in Lake Thingvallavatn, Iceland. Environmental factors, especially predation, are clearly implicated in reducing gene flow between morphs. There is continuous morphological and genetic variation between habitats with a zone centre similar to secondary contact zones. Individual microsatellite loci are implicated as being linked to adaptive variation by direct tests as well as by differences in cline shape. Patterns of linkage disequilibria indicate that the morphs have diverged at several loci. This divergence shows parallels and differences with the well-studied limnetic-benthic stickleback morphs, both in phenotypic divergence and at the genomic level.

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.

Lineage origin and expansion of a Neotropical migrant songbird after recent glaciation events

Birds of the Northern Hemisphere often harbour the genetic signature of postglaciation expansion but analyses identifying the location of refugia and the directionality of expansions are rare. Here we explore the evolutionary history of yellow warbler lineages, focusing on how these lineages recolonized their current range. We genotyped samples from 696 yellow warblers via direct sequencing of a 333-bp control region I mitochondrial DNA fragment or lineage-specific genotyping. Phylogenetic analysis revealed two monophyletic clades: a highly migratory group including previously identified eastern and western lineages and a less migratory group including a lineage consisting of tropical residents and a new 'southern' lineage localized in southwest United States. We then modelled the expansion of the eastern and western lineages, identified the location of potential refugia and assessed the importance of migration as a historical factor promoting gene flow. The expansion of the eastern lineage proceeded from a main refugia in the eastern United States, with possible contribution of an additional local refugia. In the western lineage, the expansion proceeded from a single refugia possibly located in western United States. Because two lineages overlapped to varying degrees in central North America, we suggest that the Canadian Prairies offered a bridge of riparian habitats where the lineages met after glacier retreat, while the US Central Great Plains acted as a barrier that limited secondary contact. Finally, gene flow was more important along the north-south axis of migration than away from it, suggesting spring migration played a role in the dispersal of lineages.

Gis-based niche models reveal unifying climatic mechanisms that maintain the location of avian hybrid zones in a North American suture zone

The existence of suture zones in North America has recently been verified, yet the environmental factors responsible for the maintenance of their structure and position have remained undetermined. The Great Plains suture zone in the US is perhaps the most significant and broadly studied in North America. Numerous avian hybrid zones that cluster in this region have been extensively studied over the last half-century. A primary result of this work is that exogenous factors present in this region have promoted and maintained avian species divergence. Yet, to date the variables most important in the generation and the clustered positioning of these hybrid zones in relation to one another have not been determined. The present study aims to advance our understanding of this North American suture zone by using niche modelling to quantify which environmental variables are most important in its formation and maintenance. Through the generation of fundamental niches for four hybrids and their parental species, similarities linking the hybrid zones were uncovered. The results show that temperature holds this suture zone in position and also suggest that temperature played a primary role in promoting niche differentiation among these sister species pairs during Pleistocene glacial cycles.

DO SUTURE ZONES EXIST?

Remington (1968) argued that 13 suture zones exist in North America. Remington defined a suture zone as, "a band of geographic overlap between major biotic assemblages, including some pairs of species or semispecies which hybridize in the zone" (p. 322). Although initially controversial, the idea that suture zones exist has picked up momentum over the past decade, due largely to the phylogeographic work of Hewitt, Avise, and their colleagues. Nevertheless, the reality of suture zones has not yet been subjected to rigorous analysis using statistical and geographic information system (GIS) approaches. To test for the existence of Remington's suture zones, we first identified 117 terrestrial hybrid zones in Canada and the United States through a literature search for the key words "cline," "contact zone," "hybrid zone," and "hybridization" in articles published between 1970 and 2002. The 117 hybrid zones were mapped using a GIS approach and compared with a digitized version of Remington's original suture zone map. Overall, there does appear to be an association between hybrid zones and suture zones, but this association is largely attributable to clustering of hybrid zones in only two of the 13 suture zones recognized by Remington. The results suggest that evolutionary biologists should retain some skepticism toward Remington's suture zones.

BARRIERS TO SYMPATRY BETWEEN AVIAN SIBLING SPECIES (PARIDAE: BAEOLOPHUS) IN LOCAL SECONDARY CONTACT

Range limits and secondary contact zones often occur at ecotones between major associations of habitat and climate. Therefore, understanding processes that limit sympatry between species in such areas provides an important framework for testing biogeographic and evolutionary hypotheses. Theoretical and empirical work has shown that the evolution of species borders is influenced by a complexity of factors, including gene flow from central to peripheral populations and the ability of species to adapt locally to environmental conditions. However, few studies have used bioclimatic models, combined with molecular and morphological data, to predict geographic range limits in the context of gene flow across a secondary contact zone. In this study, I applied these methods to test specific hypotheses about barriers to sympatry between closely related species where they approach and contact each other. Specifically, I examined the importance of historical isolation, local adaptation, and symmetry of gene flow in limiting sympatry and range expansion of ecologically distinct species across environmental gradients. Molecular (mitochondrial DNA, allozymes), morphological, and bioclimatic data were obtained for two avian sibling species (Baeolophus inornatus and B. ridgwayi) that exist in recent, narrow secondary contact in northern California. These species are broadly allopatric and occupy rangewide associations of oak and pinyon-juniper woodlands, respectively, although B. inornatus also inhabits mixed or juniper woodlands locally. Patterns of molecular variation generally were congruent with morphological and bioclimatic data, and support prior evidence for a history of isolation, adaptation, and divergence in distinctive, species-specific vegetation-climate associations. However, molecular and morphological clines fall east of the limit of oaks, and individuals of B. inornatus in this juniper-associated contact zone experience bioclimates that are more similar to B. ridgwayi than to B. inornatus in oak habitat. Thus, B. inornatus is able to adapt and expand locally into the range of its close relative, but not vice versa. These data support the hypothesis that gene flow is asymmetrical where peripheral populations meet at range boundaries. Physiological differences between species may play an important role in influencing these patterns. Empirical studies that highlight the importance of local adaptation and patterns of gene flow in which closely related species contact across ecotones are central to understanding limits on geographic ranges, sympatry, and introgression-a cornerstone of biogeographic and speciation theory.