LIFE-HISTORY VARIATION AND COMPARATIVE PHYLOGEOGRAPHY OF DARTERS (PISCES: PERCIDAE) FROM THE NORTH AMERICAN CENTRAL HIGHLANDS

Phylogeography seeks to evaluate the relationship between genetic variation and geographic distribution of a species to examine the influence of historical events on divergence. Congruent phylogeographic patterns in codistributed species indicate historical association of the taxa being compared, and the uniform action of biogeographic events in shaping genetic variation. We sought to evaluate the congruence of patterns of genetic variation of five closely related fish species across a well-defined biogeographic boundary. We gathered allozymic and mitochondrial DNA (mtDNA) sequence data for five species of darters (Percidae: Etheostomatini) from populations distributed among biogeographic regions in the Ozark and Ouachita highlands of the south-central United States. Comparisons among species revealed noncongruence in the magnitude of genetic divergence in both allozymes and mtDNA sequences. We hypothesized that noncongruence resulted, in part, from differences in life histories of the species in our comparison. To address this hypothesis, we evaluated the association between gene flow (measured by N_em) and variation in body size and fecundity variables because they have been shown to influence gene flow in fishes. Correlation analysis revealed an association between gene flow and fecundity (r = 0.88), but not with body size (r = 0.36) or reproductive investment per individual (r = –0.23). The result was similar when independent contrasts of the original variables were used in correlation analyses. Phylogeographic analysis of mtDNA sequence data indicated the importance of history, evident in gene trees of Percina nasuta and Percina phoxocephala. Divergence rates between these two taxa may differ because of historically persistent differences in population sizes, reflected in present-day abundance and fecundity differences. Conversely, Percina caprodes showed little evidence of divergence in mtDNA sequences and yielded the highest mean N_em values from allozyme data. Comparisons among closely related, codistributed taxa may help discriminate among the possibilities for noncongruence in biogeographic studies by allowing evaluation of the influence of life history on patterns of gene flow.     

Genetic and phenotypic variation in a colourful treefrog across five geographic barriers

Aim Understanding the patterns and processes underlying phenotype in a polytypic species provides key insights into microevolutionary mechanisms of diversification. The red‐eyed treefrog, Agalychnis callidryas, exhibits strong regional differentiation in colour pattern, corresponding to five admixed mitochondrial DNA clades. We evaluated spatial diversity patterns across multiple, putative barriers to examine the fine‐scale processes that mediate phenotypic divergence between some regions while maintaining homogeneity between others. Location We examined patterns of phenotypic diversification among 17 sites that span five putative biogeographic barriers in lower Central America (Costa Rica and Panama). Methods We tested the extent to which genetic distance (F_ST) derived from six multilocus nuclear genotypes covaried with measures of phenotypic distance (leg coloration) within and between biogeographic regions. We used linear regression analyses to determine the role of geographic and genetic factors in structuring spatial patterns of phenotypic diversity. Results The factors that best explained patterns of phenotypic diversity varied among biogeographic regions. We identified one geographic barrier that impeded gene exchange and resulted in concordant phenotypic divergence across the Continental Divide, isolating Caribbean and Pacific populations. Across Caribbean Costa Rican populations, one barrier structured phenotypic but not genetic diversity patterns, indicating a role for selection. In other regions, the putative barriers had no determining effect on either genetic or leg colour structure. Main conclusions The processes mediating the distribution and diversification of colour pattern in this polytypic, wide‐ranging treefrog varied among biogeographic regions. Spatially varying selection combined with the isolating effects of geographic factors probably resulted in the patchy distribution of colour diversity across Costa Rican and Panamanian populations.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Discordance in body size, colour pattern, and advertisement call across genetically distinct populations in a Neotropical anuran (Dendropsophus ebraccatus)

Patterns of intraspecific geographic variation in morphology and behaviour, when examined in a phylogenetic context, can provide insight into the microevolutionary processes driving population divergence and ultimately speciation. In the present study, we quantified behavioural and phenotypic variation among populations from genetically divergent regions in the Central American treefrog, Dendropsophus ebraccatus . Our fine-scale population comparisons demonstrated regional divergence in body size, colour pattern frequencies, and male advertisement call. None of the characters covaried with phylogenetic history or geographic proximity among sampled populations, indicating the importance of highly localized selection pressures and genetic drift in shaping character divergence among isolated regions. The study underscores how multiple phenotypic characters can evolve independently across relatively small spatial scales.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 298–313.     

Rapid postglacial diversification and long‐term stasis within the songbird genus Junco: phylogeographic and phylogenomic evidence

Natural systems composed of closely related taxa that vary in the degree of phenotypic divergence and geographic isolation provide an opportunity to investigate the rate of phenotypic diversification and the relative roles of selection and drift in driving lineage formation. The genus Junco (Aves: Emberizidae) of North America includes parapatric northern forms that are markedly divergent in plumage pattern and colour, in contrast to geographically isolated southern populations in remote areas that show moderate phenotypic divergence. Here, we quantify patterns of phenotypic divergence in morphology and plumage colour and use mitochondrial DNA genes, a nuclear intron, and genomewide SNPs to reconstruct the demographic and evolutionary history of the genus to infer relative rates of evolutionary divergence among lineages. We found that geographically isolated populations have evolved independently for hundreds of thousands of years despite little differentiation in phenotype, in sharp contrast to phenotypically diverse northern forms, which have diversified within the last few thousand years as a result of the rapid postglacial recolonization of North America. SNP data resolved young northern lineages into reciprocally monophyletic lineages, indicating low rates of gene flow even among closely related parapatric forms, and suggesting a role for strong genetic drift or multifarious selection acting on multiple loci in driving lineage divergence. Juncos represent a compelling example of speciation in action, where the combined effects of historical and selective factors have produced one of the fastest cases of speciation known in vertebrates.     

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.     

The influence of life‐history strategy on genetic differentiation and lineage divergence in darters (Percidae: Etheostomatinae)

Recent studies determined that darters with specialized breeding strategies can exhibit deep lineage divergence over fine geographic scales without apparent physical barriers to gene flow. However, the extent to which intrinsic characteristics interact with extrinsic factors to influence population divergence and lineage diversification in darters is not well understood. This study employed comparative phylogeographic and population genetic methods to investigate the influence of life history on gene flow, dispersal ability, and lineage divergence in two sympatric sister darters with differing breeding strategies. Our results revealed highly disparate phylogeographic histories, patterns of genetic structure, and dispersal abilities between the two species suggesting that life history may contribute to lineage diversification in darters, especially by limiting dispersal among large river courses. Both species also showed striking differences in demographic history, indicating that extrinsic factors differentially affected each species during the Pleistocene. Collectively, our results indicate that intrinsic and extrinsic factors have influenced levels of gene flow among populations within both species examined. However, we suggest that life‐history strategy may play a more important role in lineage diversification in darters than previously appreciated, a finding that has potentially important implications for understanding diversification of the rich North American freshwater fish fauna.     

Genetic divergence among sympatric colour morphs of the Dalmatian wall lizard (<Emphasis Type="Italic">Podarcis melisellensis</Emphasis>)

If alternative phenotypes in polymorphic populations do not mate randomly, they can be used as model systems to study adaptive diversification and possibly the early stages of sympatric speciation. In this case, non random mating is expected to support genetic divergence among the different phenotypes. In the present study, we use population genetic analyses to test putatively neutral genetic divergence (of microsatellite loci) among three colour morphs of the lizard Podarcis melisellensis, which is associated with differences in male morphology, performance and behaviour. We found weak evidence of genetic divergence, indicating that gene flow is somewhat restricted among morphs and suggesting possible adaptive diversification.     

Testing Pleistocene refugia theory: phylogeographical analysis of Desmognathus wrighti,a high-elevation salamander in the southern Appalachians

During the colder climates of the Pleistocene, the ranges of high-elevation species in unglaciated areas may have expanded, leading to increased gene flow among previously isolated populations. The phylogeography of the pygmy salamander, Desmognathus wrighti, an endemic species restricted to the highest mountain peaks of the southern Appalachians, was examined to test the hypothesis that the range of D. wrighti expanded along with other codistributed taxa during the Pleistocene. Analyses of genetic variation at 14 allozymic loci and of the 12S rRNA gene in the mtDNA genome was conducted on individuals sampled from 14 population isolates throughout the range of D. wrighti. In contrast to the genetic patterns of many other high-elevation animals and plants, genetic distances derived from both molecular markers showed significant isolation by distance and genetic structuring of populations, suggesting long-term isolation of populations. Phylogeographical analyses revealed four genetically distinct population clusters that probably remained fragmented during the Pleistocene, although there was also evidence supporting recent gene flow among some population groups. Support for isolation by distance is rare among high-elevation species in unglaciated areas of North and Middle America, although not uncommon among Plethodontid Salamanders, and this pattern suggests that populations of D. wrighti did not expand entirely into suitable habitat during the Pleistocene. We propose that intrinsic barriers to dispersal, such as species interactions with other southern Appalachian plethodontid salamanders, persisted during the Pleistocene to maintain the fragmented distribution of D. wrighti and allow for significant genetic divergence of populations by restricting gene flow.     

Phylogeography of the yucca moth Tegeticula maculata: the role of historical biogeography in reconciling high genetic structure with limited speciation

Tegeticula maculata is one of the most ancient and morphologically variable lineages within the yucca moths, yet has apparently undergone little diversification in comparison with much younger yucca moth lineages that have rapidly diversified. A phylogeographic approach was used to determine the number of independent lineages within T. maculata and to examine whether these patterns corresponded with morphological differences between its subspecies maculata and extranea. Phylogenetic analysis of mitochondrial DNA sequence variation indicated that the two subspecies are in separate clades, but there was also an equally deep split within subspecies maculata. There was no evidence for gene flow among regions and there was considerable substructure within clades. The phylogeographic structure of moth populations among and within subspecies can be explained in part by historical biogeographic boundaries and increasingly patchy postglacial distribution of the exclusive host plant, Hesperoyucca whipplei. Local specialization and co-adaptation would be possible in the absence of apparent gene flow, yet gross morphological divergence is limited to the very old split between the subspecies. Sorting of ancient mitochondrial lineages followed by local genetic differentiation may explain the pattern of high genetic structure with limited speciation.     

Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.     

Historical Divergence and Gene Flow in the Genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence—and ultimately speciation—remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticate Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last ~150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.     

Biogeographical patterns of genetic differentiation in dung beetles of the genus Trypocopris (Coleoptera, Geotrupidae) inferred from mtDNA and AFLP analyses

Aim To examine the phylogeography and population structure of three dung beetle species of the genus Trypocopris (Coleoptera, Geotrupidae). We wanted to test whether genetic differences and genealogies among populations were in accordance with morphologically described subspecies and we aimed to establish times of divergence among subspecies to depict the appropriate temporal framework of their phylogeographical differentiation. We also wished to investigate the historical demographic events and the relative influences of gene flow and drift on the distribution of genetic variability of the different populations. Location Europe (mostly Italy). Methods We collected adult males from dung pats from 15 Italian localities over the period 2000–2002. For sequence analysis, some dried specimens from Albania, Croatia, Slovakia and Spain were also used. We applied cytochrome oxidase I mitochondrial DNA sequencing and the amplified fragment length polymorphism (AFLP) technique to determine whether phylogeographical patterns within the three species support the proposed hypotheses of subspecies designations, and to detect further structure among populations that might mediate diversification. Results and main conclusions The results show a high concordance between the distribution of mtDNA variation and the main morphological groups recognized as subspecies, which thus may represent independent evolutionary units. The degree of mitochondrial divergence suggests that speciation events occurred during the Pliocene, while diversification of the main subspecific lineages took place in the Pleistocene, from c. 0.3 to 1.5 Ma. Mitochondrial and nuclear data also reveal that there is phylogeographical structuring among populations within each of the main groups and that both contemporary and historical processes determined this pattern of genetic structure. Geographical populations form monophyletic clades in both phylogenetic and network reconstructions. Despite the high levels of intrapopulational diversity, F_ST values indicate moderate but significant genetic differentiation among populations, and a Bayesian clustering analysis of the AFLP data clearly separates the geographical populations. Nucleotide and gene diversity estimates reveal interspecific differences in the degree of diversification among populations that may be related to the different ecological requirements of the three species.     

Comparative phylogeography of three codistributed stomatopods: origins and timing of regional lineage diversification in the Coral Triangle

The Indonesian-Australian Archipelago is the center of the world's marine biodiversity. Although many biogeographers have suggested that this region is a "center of origin," criticism of this theory has focused on the absence of processes promoting lineage diversification in the center. In this study we compare patterns of phylogeographic structure and gene flow in three codistributed, ecologically similar Indo-West Pacific stomatopod (mantis shrimp) species. All three taxa show evidence for limited gene flow across the Maluku Sea with deep genetic breaks between populations from Papua and Northern Indonesia, suggesting that limited water transport across the Maluku Sea may limit larval dispersal and gene flow across this region. All three taxa also show moderate to strong genetic structure between populations from Northern and Southern Indonesia, indicating limited gene flow across the Flores and Java Seas. Despite the similarities in phylogeographic structure, results indicate varied ages of the genetic discontinuities, ranging from the middle Pleistocene to the Pliocene. Concordance of genetic structure across multiple taxa combined with temporal discordance suggests that regional genetic structures have arisen from the action of common physical processes operating over extended time periods. The presence in all three species of both intraspecific genetic structure as well as deeply divergent lineages that likely represent cryptic species suggests that these processes may promote lineage diversification within the Indonesian-Australian Archipelago, providing a potential mechanism for the center of origin. Efforts to conserve biodiversity in the Coral Triangle should work to preserve both existing biodiversity as well as the processes creating the biodiversity.     

Speciation in Chitaura grasshoppers (Acrididae: Oxyinae) on the island of Sulawesi: colour patterns, morphology and contact zones

We investigate colour pattern and morphological variation in Chitaura grasshoppers on the Indonesian island of Sulawesi, and examine the relationship between divergence in these two sets of characters and population history, as reflected by variation in mitochondrial DNA. Analysis of colour pattern variation identifies a total of 21 distinct Operational Taxonomic Units in Sulawesi, the majority of which have parapatric distributions. Patterns of phenotypic variation at two contact zones in north Sulawesi suggest genetic independence between three of these colour forms, indicating that speciation has occurred. Despite this, colour pattern divergence is only coincident with morphological differentiation at one of these contact zones. In addition, neither type of phenotypic divergence is associated with geographical structuring in mitochondrial DNA, suggesting that historical isolation has had little influence on evolutionary diversification. Instead, divergence in colour pattern and morphology appears to have occurred rapidly or under conditions of continued gene flow, possibly in response to spatially variable natural selection. This result has implications for the identification of conservation units based solely on molecular markers.     

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae)

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.     

Using regional comparative phylogeographic data from snake lineages to infer historical processes in Middle America

Understanding how historical processes have either similarly, or differentially, shaped the evolution of lineages or biotic assemblages is important for a broad spectrum of fields. Gaining such understanding can be particularly challenging, however, especially for regions that have a complex geologic and biological history. In this study we apply a broad comparative approach to distill such regional biogeographic perspectives, by characterizing sets of divergence times for major biogeographic boundaries estimated from multiple codistributed lineages of snakes. We use a large combined (mitochondrial gene sequence) phylogeographic/phylogenetic dataset containing several clades of snakes that range across Middle America – the tropical region between Mexico and northwestern South America. This region is known for its complex tectonic history, and poorly understood historical biogeography. Based on our results, we highlight how phylogeographic transition zones between Middle and South America and the Nicaragua Depression appear to have undergone multiple episodes of diversification in different lineages. This is in contrast to other examples we find where apparently a single vicariant period is shared across multiple lineages. We specifically evaluate the distributions of divergence time estimates across multiple lineages and estimate the number of temporal periods of lineage diversification per biogeographic break. Overall, our results highlight a great deal of shared temporal divergence, and provide important hypotheses for yet unstudied lineages. These multi‐lineage comparisons across multiple spatial and temporal scales provide excellent predictive power for identifying the roles of geology, climate, ecology and natural history in shaping regional biodiversity.     

Divergence genetics analysis reveals historical population genetic processes leading to contrasting phylogeographic patterns in co-distributed species

Coalescent samplers are computational time machines for inferring the historical demographic genetic processes that have given rise to observable patterns of spatial genetic variation among contemporary populations. We have used traditional characterizations of population structure and coalescent‐based inferences about demographic processes to reconstruct the population histories of two co‐distributed marine species, the frilled dog whelk, Nucella lamellosa, and the bat star, Patiria miniata. Analyses of population structure were consistent with previous work in both species except that additional samples of N. lamellosa showed a larger regional genetic break on Vancouver Island (VI) rather than between the southern Alexander Archipelago as in P. miniata. Our understanding of the causes, rather than just the patterns, of spatial genetic variation was dramatically improved by coalescent analyses that emphasized variation in population divergence times. Overall, gene flow was greater in bat stars (planktonic development) than snails (benthic development) but spatially homogeneous within species. In both species, these large phylogeographic breaks corresponded to relatively ancient divergence times between populations rather than regionally restricted gene flow. Although only N. lamellosa shows a large break on VI, population separation times on VI are congruent between species, suggesting a similar response to late Pleistocene ice sheet expansion. The absence of a phylogeographic break in P. miniata on VI can be attributed to greater gene flow and larger effective population size in this species. Such insights put the relative significance of gene flow into a more comprehensive historical biogeographic context and have important implications for conservation and landscape genetic studies that emphasize the role of contemporary gene flow and connectivity in shaping patterns of population differentiation.     

The influence of gene flow and drift on genetic and phenotypic divergence in two species of Zosterops in Vanuatu

Colonization of an archipelago sets the stage for adaptive radiation. However, some archipelagos are home to spectacular radiations, while others have much lower levels of diversification. The amount of gene flow among allopatric populations is one factor proposed to contribute to this variation. In island colonizing birds, selection for reduced dispersal ability is predicted to produce changing patterns of regional population genetic structure as gene flow-dominated systems give way to drift-mediated divergence. If this transition is important in facilitating phenotypic divergence, levels of genetic and phenotypic divergence should be associated. We consider population genetic structure and phenotypic divergence among two co-distributed, congeneric (Genus: Zosterops) bird species inhabiting the Vanuatu archipelago. The more recent colonist, Z. lateralis, exhibits genetic patterns consistent with a strong influence of distance-mediated gene flow. However, complex patterns of asymmetrical gene flow indicate variation in dispersal ability or inclination among populations. The endemic species, Z. flavifrons, shows only a partial transition towards a drift-mediated system, despite a long evolutionary history on the archipelago. We find no strong evidence that gene flow constrains phenotypic divergence in either species, suggesting that levels of inter-island gene flow do not explain the absence of a radiation across this archipelago.     

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.