Apparent ‘sympatric’ speciation in ecologically similar herbivorous beetles facilitated by multiple colonizations of an island

Predicting the dispersal and invasion dynamics of ambrosia beetles through demographic reconstruction and process-explicit modeling

Evaluating potential routes of invasion of pathogens and vectors of sanitary importance is essential for planning and decision-making at multiple scales. An effective tool are process-explicit models that allow coupling environmental, demographic and dispersal information to evaluate population growth and range dynamics as a function of the abiotic conditions in a region. In this work we simulate multiple dispersal/invasion routes in Mexico that could be taken by ambrosia beetles and a specific symbiont, Harringtonia lauricola, responsible for a severe epiphytic of Lauraceae in North America. We used Xyleborus bispinatus Eichhoff 1868 as a study subject and estimated its demography in the laboratory in a temperature gradient (17, 20, 26, 29, 35 °C), which we then used to parameterize a process-based model to estimate its metapopulation dynamics. The maximum intrinsic growth rate of X. bispinatus is 0.13 with a thermal optimum of 26.2 °C. The models suggest important regions for the establishment and dispersal the states of Veracruz, Chiapas and Oaxaca (high host and secondary vectors diversity), the Isthmus of Tehuantepec (connectivity region), and Michoacán and Jalisco (important avocado plantations). The use of hybrid process-based models is a promising tool to refine the predictions applied to the study of biological invasions and species distributions.

Two sympatric lineages of Australian Cnestus solidus share Ambrosiella symbionts but not Wolbachia

Sympatric lineages of inbreeding species provide an excellent opportunity to investigate species divergence patterns and processes. Many ambrosia beetle lineages (Curculionidae: Scolytinae) reproduce by predominant inbreeding through sib mating in nests excavated in woody plant parts wherein they cultivate symbiotic ambrosia fungi as their sole source of nutrition. The Xyleborini ambrosia beetle species Cnestus solidus and Cnestus pseudosolidus are sympatrically distributed across eastern Australia and have overlapping morphological variation. Using multilocus sequencing analysis of individuals collected from 19 sites spanning their sympatric distribution, we assessed their phylogenetic relationships, taxonomic status and microbial symbionts. We found no genetic differentiation between individuals morphologically identified as C. solidus and C. pseudosolidus confirming previous suggestions that C. pseudosolidus is synonymous to C. solidus. However, within C. solidus we unexpectedly discovered the sympatric coexistence of two morphologically indistinguishable but genetically distinct lineages with small nuclear yet large mitochondrial divergence. At all sites except one, individuals of both lineages carried the same primary fungal symbiont, a new Ambrosiella species, indicating that fungal symbiont differentiation may not be involved in lineage divergence. One strain of the maternally inherited bacterial endosymbiont Wolbachia was found at high prevalence in individuals of the more common lineage but not in the other, suggesting that it may influence host fitness. Our data suggest that the two Australian Cnestus lineages diverged allopatrically, and one lineage then acquired Wolbachia. Predominant inbreeding and Wolbachia infection may have reinforced reproductive barriers between these two lineages after their secondary contact contributing to their current sympatric distribution.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Host and geography together drive early adaptive radiation of Hawaiian planthoppers

The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.

Revision of the Hawaiian psyllid genus Swezeyana, with descriptions of seven new species (Hemiptera, Psylloidea, Triozidae)

The endemic Hawaiian genus Swezeyana Caldwell, 1940 is highly distinctive due to the extremely long genal processes. In addition, some of the immatures are ornamented with extraordinary tubercles and tentacles. Two Swezeyana species are redescribed, and seven new species are described, bringing the total number of species in the genus to nine. All species are hosted by a single, endemic host plant, Planchonellasandwicensis (Sapotaceae), which is distributed across all major islands in the archipelago. The majority of Swezeyana species are single island endemics. A sister taxon pair is found sympatrically on the same individual plants on Kauai, and putative sister or at least closely related species are also found sympatrically on Oahu and Hawaii, suggesting these taxa may have diversified in sympatry. However, there is no observed ecological niche partitioning, despite some striking morphological diversity, as all Swezeyana species have free-living immatures that are found on the leaf surface, and therefore no apparent biological shifts are coincident with occupying the same host plant. Two species groups are represented by strikingly different female terminalia structure and endoskeletal development, although ovipositor structure is very similar between the two groups. Mitochondrial DNA barcodes (COI and cytB) are provided for eight of the nine species. A phylogenetic analysis of the mitochondrial barcode regions indicates species relationships within Swezeyana and provides a comparison of genetic divergence with other Hawaiian endemic genera.

Phylogenetic analysis of the genus Laparocerus, with comments on colonisation and diversification in Macaronesia (Coleoptera, Curculionidae, Entiminae)

The flightless Entiminae weevil genus Laparocerus is the species-richest genus, with 237 species and subspecies, inhabiting Macaronesia (Madeira archipelago, Selvagens, Canary Islands) and the continental 'Macaronesian enclave' in Morocco (one single polytypic species). This is the second contribution to gain insight of the genus and assist in its systematic revision with a mitochondrial phylogenetic analysis. It centres on the Canarian clade, adding the 12S rRNA gene to the combined set of COII and 16S rRNA used in our first contribution on the Madeiran clade (here re-analysed). The nuclear 28S rRNA was also used to produce an additional 4-gene tree to check coherency with the 3-gene tree. A total of 225 taxa (95%) has been sequenced, mostly one individual per taxa. Plausible explanations for incoherent data (mitochondrial introgressions, admixture, incomplete lineage sorting, etc.) are discussed for each of the monophyletic subclades that are coincident with established subgenera, or are restructured or newly described. The overall mean genetic divergence (p-distance) among species is 8.2%; the mean divergence within groups (subgenera) ranks from 2.9 to 7.0% (average 4.6%), and between groups, from 5.4% to 12.0% (average 9.2%). A trustful radiation event within a young island (1.72 Ma) was used to calibrate and produce a chronogram using the software RelTime. These results confirm the monophyly of both the Madeiran (36 species and subspecies) and the Canarian (196 species and subspecies) clades, which originated ca. 11.2 Ma ago, and started to radiate in their respective archipelagos ca. 8.5 and 7.7 Ma ago. The Madeiran clade seems to have begun in Porto Santo, and from there it jumped to the Desertas and to Madeira, with additional radiations. The Canarian clade shows a sequential star-shape radiation process generating subclades with a clear shift from East to West in coherence with the decreasing age of the islands. Laparocerus garretai from the Selvagens belongs to a Canarian subclade, and Laparocerus susicus from Morocco does not represent the ancestral continental lineage, but a back-colonisation from the Canaries to Africa. Dispersal processes, colonisation patterns, and ecological remarks are amply discussed. Diversification has been adaptive as well as non-adaptive, and the role of 'geological turbulence' is highlighted as one of the principal drivers of intra-island allopatric speciation. Based on the phylogenetic results, morphological features and distribution, five new monophyletic subgenera are described: Aridotroxsubg. n., Belicariussubg. n., Bencomiussubg. n., Canariotroxsubg. n., and Purpuraniussubg. n., totalling twenty subgenera in Laparocerus.

Evolutionary consequences of climate-induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short- and long-term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short-term processes acting on reproductive barriers between species and (ii) long-term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate-induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.

Fission and fusion in island taxa--serendipity, or something to be expected?

A well-used metaphor for oceanic islands is that they act as 'natural laboratories' for the study of evolution. But how can islands or archipelagos be considered analogues of laboratories for understanding the evolutionary process itself? It is not necessarily the case that just because two or more related species occur on an island or archipelago, somehow, this can help us understand more about their evolutionary history. But in some cases, it can. In this issue of Molecular Ecology, Garrick et al. () use population-level sampling within closely related taxa of Galapagos giant tortoises to reveal a complex demographic history of the species Chelonoidis becki - a species endemic to Isabela Island, and geographically restricted to Wolf Volcano. Using microsatellite genotyping and mitochondrial DNA sequencing, they provide a strong case for C. becki being derived from C. darwini from the neighbouring island of Santiago. But the interest here is that colonization did not happen only once. Garrick et al. () reveal C. becki to be the product of a double colonization event, and their data reveal these two founding lineages to be now fusing back into one. Their results are compelling and add to a limited literature describing the evolutionary consequences of double colonization events. Here, we look at the broader implications of the findings of Garrick et al. () and suggest genomic admixture among multiple founding populations may be a characteristic feature within insular taxa.

Mitochondrial capture misleads about ecological speciation in the Daphnia pulex complex

The North American ecological species Daphniapulicaria and Daphniapulex are thought to have diverged from a common ancestor by adaptation to sympatric but ecologically distinct lake and pond habitats respectively. Based on mtDNA relationships, European D. pulicaria is considered a different species only distantly related to its North American counterpart, but both species share a lactate dehydrogenase (Ldh) allele F supposedly involved in lake adaptation in North America, and the same allele is also carried by the related Holarctic Daphniatenebrosa. The correct inference of the species' ancestral relationships is therefore critical for understanding the origin of their adaptive divergence. Our species tree inferred from unlinked nuclear loci for D. pulicaria and D. pulex resolved the European and North American D. pulicaria as sister clades, and we argue that the discordant mtDNA gene tree is best explained by capture of D. pulex mtDNA by D. pulicaria in North America. The Ldh gene tree shows that F-class alleles in D. pulicaria and D. tenebrosa are due to common descent (as opposed to introgression), with D. tenebrosa alleles paraphyletic with respect to D. pulicaria alleles. That D. tenebrosa still segregates the ancestral and derived amino acids at the two sites distinguishing the pond and lake alleles suggests that D. pulicaria inherited the derived states from the D. tenebrosa ancestry. Our results suggest that some adaptations restricting the gene flow between D. pulicaria and D. pulex might have evolved in response to selection in ancestral environments rather than in the species' current sympatric habitats. The Arctic (D. tenebrosa) populations are likely to provide important clues about these issues.

Ecological speciation in dynamic landscapes

Although verbal theories of speciation consider landscape changes, ecological speciation is usually modelled in a fixed geographical arrangement. Yet landscape changes occur, at different spatio-temporal scales, due to geological, climatic or ecological processes, and these changes result in repeated divisions and reconnections of populations. We examine the effect of such landscape dynamics on speciation. We use a stochastic, sexual population model with polygenic inheritance, embedded in a landscape dynamics model (allopatry-sympatry oscillations). We show that, under stabilizing selection, allopatry easily generates diversity, but species coexistence is evolutionarily unsustainable. Allopatry produces refuges whose persistence depends on the characteristic time scales of the landscape dynamics. Under disruptive selection, assuming that sympatric speciation is impossible due to Mendelian inheritance, allopatry is necessary for ecological differentiation. The completion of reproductive isolation, by reinforcement, then requires several sympatric phases. These results demonstrate that the succession of past, current and future geographical arrangements considerably influence the speciation process.

Allopatric origin of cryptic butterfly species that were discovered feeding on distinct host plants in sympatry

Surveys of tropical insects are increasingly uncovering cryptic species - morphologically similar yet reproductively isolated taxa once thought to comprise a single interbreeding entity. The vast majority of such species are described from a single location. This leaves us with little information on geographic range and intraspecific variation and limits our ability to infer the forces responsible for generating such diversity. For example, in herbivorous and parasitic insects, multiple specialists are often discovered within what were thought to be single more generalized species. Host shifts are likely to have contributed to speciation in these cases. But when and where did those shifts occur, and were they facilitated by geographic isolation? We attempted to answer these questions for two cryptic species within the butterfly Cymothoe egesta that were recently discovered on different host plants in central Cameroon. We first used mtDNA markers to separate individuals collected on the two hosts within Cameroon and then extended our analysis to incorporate individuals collected across the entire pan-Afrotropical range of the original taxon. To our surprise, we found that the species are almost entirely allopatric, dividing the original range and overlapping only in the narrow zone of West-Central Africa where they were first discovered in sympatry. This finding, combined with analyses of genetic variation within each butterfly species, strongly suggests that speciation occurred in allopatry, probably during the Pleistocene. We discuss the implications of our results for understanding speciation among other cryptic species recently discovered in the tropics and argue that more work is needed on geographic patterns and host usage in such taxa.

Effects of gallery density and species ratio on the fitness and fecundity of two sympatric bark beetles (Coleoptera: Curculionidae)

Interspecific interactions among tree-killing bark beetle species may have ecologically important consequences on beetle population dynamics. Using two tree-killing beetle species (Dendroctonus brevicomis and D. frontalis), we performed observational and experimental studies to verify cross-attraction and co-colonization under field conditions in northern Arizona and test the effects of gallery density and species ratio on response variables of average gallery length, offspring size (progeny fitness), and offspring production per centimeter gallery (fecundity). Our results show that both D. frontalis and D. brevicomis aggregate to pheromones synthesized de novo by D. brevicomis under field conditions and that galleries of both D. brevicomis and D. frontalis occurred together in the same region of a single host tree with significant frequency. In experimental manipulations of species ratios, the presence of conspecific beetles in the gallery environment strongly mediated fecundity, but D. frontalis was the only species that suffered negative impacts from the presence of heterospecific beetles in the gallery environment. Interactions did not result in any apparent fitness effects for progeny of either species, which suggests that multispecies aggregations and co-colonization may be a dominant ecological strategy in the region and result in niche sharing.

Does a predatory insect contribute to the divergence between cave- and surface-adapted fish populations?

Immigrant inviability, where individuals from foreign, ecologically divergent habitats are less likely to survive, can restrict gene flow among diverging populations and result in speciation. I investigated whether a predatory aquatic insect (Belostoma sp.) selects against migrants between cave and surface populations of a fish (Poecilia mexicana). Cavefish were more susceptible to attacks in the light, whereas surface fish were more susceptible in darkness. Environmentally dependent susceptibility to attacks may thus contribute to genetic and phenotypic differentiation between the populations. This study highlights how predation-in this case in conjunction with differences in other environmental factors-can be an important driver in speciation.

Genetic diversity in two sibling species of the Anopheles punctulatus group of mosquitoes on Guadalcanal in the Solomon Islands

BACKGROUND

The mosquito Anopheles irenicus, a member of the Anopheles punctulatus group, is geographically restricted to Guadalcanal in the Solomon Islands. It shows remarkable morphological similarities to one of its sibling species, An. farauti sensu stricto (An. farauti s.s.), but is dissimilar in host and habitat preferences. To infer the genetic variations between these two species, we have analyzed mitochondrial cytochrome oxidase subunit II (COII) and nuclear ribosomal internal transcribed spacer 2 (ITS2) sequences from Guadalcanal and from one of its nearest neighbours, Malaita, in the Solomon Islands.

RESULTS

An. farauti s.s. was collected mostly from brackish water and by the human bait method on both islands, whereas An. irenicus was only collected from fresh water bodies on Guadalcanal Island. An. irenicus is distributed evenly with An. farauti s.s. (Phi SC = 0.033, 0.38%) and its range overlaps in three of the seven sampling sites. However, there is a significant population genetic structure between the species (Phi CT = 0.863, P < 0.01; Phi ST = 0.865, P < 0.01 and FST = 0.878, P < 0.01). Phylogenetic analyses suggest that An. irenicus is a monophyletic species, not a hybrid, and is closely related to the An. farauti s.s. on Guadalcanal. The time estimator suggests that An. irenicus diverged from the ancestral An. farauti s.s. on Guadalcanal within 29,000 years before present (BP). An. farauti s.s. expanded much earlier on Malaita (texp = 24,600 BP) than the populations on Guadalcanal (texp = 16,800 BP for An. farauti s.s. and 14,000 BP for An. irenicus).

CONCLUSION

These findings suggest that An. irenicus and An. farauti s.s. are monophyletic sister species living in sympatry, and their populations on Guadalcanal have recently expanded. Consequently, the findings further suggest that An. irenicus diverged from the ancestral An. farauti s.s. on Guadalcanal.

Nuclear and mitochondrial sequences confirm complex colonization patterns and clear species boundaries for flightless weevils in the Galápagos archipelago

Nuclear sequence data were collected from endemic Galápagos species and an introduced close relative, and contrasted with mitochondrial DNA sequences, continuing investigation into the colonization history and modes of diversification in the weevil genus Galapaganus. The current combined phylogeny together with previously published penalized likelihood age estimates builds a complex picture of the archipelago's colonization history. The present reconstruction relies on submerged platforms to explain the early divergence of the young southern Isabela endemics or the Española or San Cristobal populations. Diversity is later built through inter-island divergence starting on older islands and continuing on two simultaneous tracks towards younger islands. The amount of diversity generated through intra-island processes is skewed towards older islands, suggesting that island age significantly influences diversity. Phylogenetic concordance between nuclear and mitochondrial datasets and well-supported monophyletic species in mitochondrial derived topologies appear to reject the possibility of inter-species hybridization. These clear species boundaries might be related to the tight host associations of adult weevils in discrete ecological zones. If shared hosts facilitate hybridization, then host- or habitat-promoted divergences could prevent it, even in the case of species that share islands, since the altitudinal partitioning of habitats minimizes range overlap.

Phylogenetic analysis of community assembly and structure over space and time

Evolutionary ecologists are increasingly combining phylogenetic data with distributional and ecological data to assess how and why communities of species differ from random expectations for evolutionary and ecological relatedness. Of particular interest have been the roles of environmental filtering and competitive interactions, or alternatively neutral effects, in dictating community composition. Our goal is to place current research within a dynamic framework, specifically using recent phylogenetic studies from insular environments to provide an explicit spatial and temporal context. We compare communities over a range of evolutionary, ecological and geographic scales that differ in the extent to which speciation and adaptation contribute to community assembly and structure. This perspective allows insights into the processes that can generate community structure, as well as the evolutionary dynamics of community assembly.

Evaluating nuclear protein-coding genes for phylogenetic utility in beetles

Although nuclear protein-coding genes have proven broadly useful for phylogenetic inference, relatively few such genes are regularly employed in studies of Coleoptera, the most diverse insect order. We increase the number of loci available for beetle systematics by developing protocols for three genes previously unused in beetles (alpha-spectrin, RNA polymerase II and topoisomerase I) and by refining protocols for five genes already in use (arginine kinase, CAD, enolase, PEPCK and wingless). We evaluate the phylogenetic performance of each gene in a Bayesian framework against a presumably known test phylogeny. The test phylogeny covers 31 beetle specimens and two outgroup taxa of varying age, including three of the four extant beetle suborders and a denser sampling in Adephaga and in the carabid genus Bembidion. All eight genes perform well for Cenozoic divergences and accurately separate closely related species within Bembidion, but individual genes differ markedly in accuracy over the older Mesozoic and Permian divergences. The concatenated data reconstruct the test phylogeny with high support in both Bayesian and parsimony analyses, indicating that combining data from multiple nuclear loci will be a fruitful approach for assembling the beetle tree of life.

The hidden side of invasions: massive introgression by local genes

Despite hundreds of reports involving both plants and animals, the mechanisms underlying introgression remain obscure, even if some form of selection is frequently invoked. Introgression has repeatedly been reported in species that have recently colonized a new habitat, suggesting that demographic processes should be given more attention for understanding the mechanisms of introgression. Here we show by spatially explicit simulations that massive introgression of neutral genes takes place during the invasion of an occupied territory if interbreeding is not severely prevented between the invading and the local species. We also demonstrate that introgression occurs almost exclusively from the local to the invading species, especially for populations located far away from the source of the invasion, and this irrespective of the relative densities of the two species. This pattern is strongest at markers experiencing reduced gene flow, in keeping with the observation that organelle genes are often preferentially introgressed across species boundaries. A survey of the literature shows that a majority of published empirical studies of introgression during range expansions, in animals and in plants, follow the predictions of our model. Our results imply that speciation genes can be identified by comparing genomes of interfertile native and invading species pairs.

Colonization history, ecological shifts and diversification in the evolution of endemic Galápagos weevils

Mitochondrial DNA sequence data were obtained for eight species of flightless Galapaganus endemic weevils and one winged close relative in order to study their colonization history and modes of diversification in the Galápagos Archipelago. Contrary to most other insular radiations, the phylogeny estimates we recovered for Galapaganus do not follow the progression rule of island biogeography. The penalized likelihood age estimates of colonization of the archipelago exceed the age of the emerged islands and underscore the potential role of now sunken seamounts for the early evolution of Galapaganus. The phylogeny proposes one intra-island origin for Galapaganus endemics, but monophyly tests suggest a larger contribution of in-situ speciation on older islands. Generalist habitat preferences were reconstructed as ancestral while shifts to highland habitats were reconstructed as having evolved independently on different islands. Magnitudes and patterns of diversification rate were found to differ between older and younger islands. Our analyses reveal that the colonization sequence of islands and timing of colonization of Galapaganus could be linked with the geological and volcanic history of the islands in a rather complex scenario. Even though most islands appear to have been colonized soon after their emergence, there are notable deviations from the pattern of sequential colonization expected under the progression rule when considering only the extant emerged islands. Patterns of diversification rate variation on older and younger islands correspond to the volcanic activity or remnants of such activity, while the pattern of independent evolution of restricted habitat preferences in different islands suggests that habitat shifts could also have contributed to species diversity in Galapaganus.

DNA based cladograms augment the discovery of a new Ips species from China (Coleoptera: Curculionidae: Scolytinae)

The implementation of DNA in taxonomic study is in its infancy because the association of the amount and type of nucleotide change with species boundaries has not been fully examined for most taxa. Mitochondrial cytochrome c oxidase I (COI) nucleotide data is currently the most popular molecular marker for delimiting species boundaries and a standard pair-wise nucleotide divergence between groups of individuals has been suggested for the recognition of new species. It is unlikely that such a standard would be applicable across animal species, but the association of the amount and type of nucleotide change with species boundaries could help with the establishment of a taxon-specific DNA taxonomy. This study utilizes DNA data from nuclear and mitochondrial genes to improve the taxonomy of an important forest beetle pest, Ips. Amount and type of nucleotide difference are associated with monophyletic species based on a cladistic analysis of these data. As a result, a new species from China is described for a clade of beetles whose nucleotide differences exceeded the amount of evolutionary change observed within currently recognized species. The COI data are analyzed independently with an expanded taxon data set, including pair-wise nucleotide differences between recognized sister species. The wide range of average intraspecific pair-wise nucleotide difference (0-10.0%) suggests limitations to the application of a standard percent nucleotide difference as a means to identify species boundaries. At most, average COI nucleotide intraspecific difference provides an informal guide to identify potential clades that may warrant further systematic investigation.

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.