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

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

Niche divergence promotes rapid diversification of East African sky island white‐eyes (Aves: Zosteropidae)

The Eastern Afromontane biodiversity hotspot composed of highly fragmented forested highlands (sky islands) harbours exceptional diversity and endemicity, particularly within birds. To explain their elevated diversity within this region, models founded on niche conservatism have been offered, although detailed phylogeographic studies are limited to a few avian lineages. Here, we focus on the recent songbird genus Zosterops, represented by montane and lowland members, to test the roles of niche conservatism versus niche divergence in the diversification and colonization of East Africa's sky islands. The species‐rich white‐eyes are a typically homogeneous family with an exceptional colonizing ability, but in contrast to their diversity on oceanic islands, continental diversity is considered depauperate and has been largely neglected. Molecular phylogenetic analysis of ~140 taxa reveals extensive polyphyly among different montane populations of Z. poliogastrus. These larger endemic birds are shown to be more closely related to taxa with divergent habitat types, altitudinal distributions and dispersal abilities than they are to populations of restricted endemics that occur in neighbouring montane forest fragments. This repeated transition between lowland and highland habitats over time demonstrate that diversification of the focal group is explained by niche divergence. Our results also highlight an underestimation of diversity compared to morphological studies that has implications for their taxonomy and conservation. Molecular dating suggests that the spatially extensive African radiation arose exceptionally rapidly (1–2.5 Ma) during the fluctuating Plio‐Pleistocene climate, which may have provided the primary driver for lineage diversification.     

Patterns of speciation are similar across mountainous and lowland regions for a Neotropical plant radiation (Costaceae: Costus)

High species richness and endemism in tropical mountains are recognized as major contributors to the latitudinal diversity gradient. The processes underlying mountain speciation, however, are largely untested. The prevalence of steep ecogeographic gradients and the geographic isolation of populations by topographic features are predicted to promote speciation in mountains. We evaluate these processes in a species-rich Neotropical genus of understory herbs that range from the lowlands to montane forests and have higher species richness in topographically complex regions. We ask whether climatic niche divergence, geographic isolation, and pollination shifts differ between mountain-influenced and lowland Amazonian sister pairs inferred from a 756-gene phylogeny. Neotropical Costus ancestors diverged in Central America during a period of mountain formation in the last 3 million years with later colonization of Amazonia. Although climatic divergence, geographic isolation, and pollination shifts are prevalent in general, these factors do not differ between mountain-influenced and Amazonian sister pairs. Despite higher climatic niche and species diversity in the mountains, speciation modes in Costus appear similar across regions. Thus, greater species richness in tropical mountains may reflect differences in colonization history, diversification rates, or the prevalence of rapidly evolving plant life forms, rather than differences in speciation mode.     

Genomic insight into “sky island” species diversification in a mountainous biodiversity hotspot

“Sky island” species diversification contributes greatly to mountainous biodiversity. However, the underlying genomic divergence and the inferred drivers remain largely unknown. In this study, we examined the diversification history of five diploid species with three exclusively endemic to the sky islands (mountains) of the Himalaya–Hengduan Mountains biodiversity hotspot. All of them together comprise a clade of the genus Eutrema (Brassicaceae). We resequenced genomes of multiple individuals of the found populations for each species. We recovered the inconsistent phylogenetic relationships between plastome and nuclear‐genome trees for one species. Based on nuclear population genomic data, we detected high genetic divergence between five species with limited gene flow. Four species seemed to diverge mainly through geographical isolation, whereas one arose through hybrid origin. The origins of the sampled five species were dated to within the late Miocene when mountains were uplifted and climates oscillated. All species decreased their population sizes since the inferred origin of each species initially, but only two of them expanded after the Last Glacial Maximum. Together, these findings suggest that geographic isolation plays an important role in driving the sky island species diversification of the sampled species in addition to the occasional gene flow that might have led to the hybrid origin of some sky island species, similar to the species diversification of sea islands.     

Diversification of montane species via elevation shifts: the case of the Kaçkar cricket Phonochorion (Orthoptera)

Localized up‐down altitudinal shifts and subsequent isolation–admixture of montane species in response to glacial cycles has been proposed as a mechanism for the high diversity along Anatolian mountains. However, specific predictions of the proposed mechanism (the elevation shift model) have yet to be tested. Here, we provide a first assessment of this model for promoting inter‐ and intraspecific genetic diversity in the bush‐cricket genus Phonochorion endemic to the West Lesser Caucasus hotspot. Mitochondrial genes were analysed by Bayesian Markov Chain Monte Carlo inferences and coalescent simulations. Timing of diversification was estimated using a multispecies coalescent model. Divergence with gene flow was tested using an isolation with migration model. Population genetic parameters and genetic structuring were determined using Bayesian coalescent methods and spatial analysis. Demographic history was assessed using mismatch distributions and extended Bayesian skyline plots. Speciation events corresponded both to the Miocene and Pleistocene while intraspecific divergence was Pleistocene based. There was evidence for moderate levels of gene flow between species during diversification; however, incomplete lineage sorting could explain the data as well as gene flow. Overall diversification patterns within the genus Phonochorion agree with the predictions of the elevations shift model. Genetic patterns of diversification were driven mainly by Pleistocene glacial cycles and reflected the nature and distribution of sky islands. There was also some albeit weak evidence of demographic expansions coinciding with glacial cooling. However, evidence for divergence with gene flow was inconclusive.     

Accounting for data heterogeneity in patterns of biodiversity: an application of linear mixed effect models to the oceanic island biogeography of spore‐producing plants

The general dynamic model of oceanic island biogeography describes the evolution of species diversity properties, including species richness (SR), through time. We investigate the hypothesis that SR in organisms with high dispersal capacities is better predicted by island area and elevation (as a surrogate of habitat diversity) than by time elapsed since island emergence and geographic isolation. Linear mixed effect models (LMMs) subjected to information theoretic model selection were employed to describe moss and liverwort SR patterns from 67 oceanic islands across 12 archipelagos. Random effects, which are used to modulate model parameters to take differences among archipelagos into account, included only a random intercept in the best‐fit model for liverworts and in one of the two best‐fit models for mosses. In this case, the other coefficients are constant across archipelagos, and we interpret the intercept as a measure of the intrinsic carrying capacity of islands within each archipelago, independently of their size, age, elevation and geographic isolation. The contribution of area and elevation to the models was substantially higher than that of time, with the least contribution made by measures of geographic isolation. This reinforces the idea that oceanic barriers are not a major impediment for migration in bryophytes and, together with the almost complete absence of in situ insular diversification, explains the comparatively limited importance of time in the models. We hence suggest that time per se has little independent role in explaining bryophyte SR and principally features as a variable accounting for the changing area and topographic complexity during the life‐cycle of oceanic islands. Simple area models reflecting habitat availability and diversity might hence prevail over more complex temporal models reflecting in‐situ speciation and dispersal (time, geographic connectivity) in explaining patterns of biodiversity for exceptionally mobile organisms.     

Docosahexaenoic and eicosapentaenoic acids increase prion formation in neuronal cells

Background

Genetic tests of paleoecological hypotheses have been rare, partly because recent genetic divergence is difficult to detect and time. According to fossil plant data, continuous woodland in the southwestern USA and northern Mexico became fragmented during the last 10,000 years, as warming caused cool-adapted species to retreat to high elevations. Most genetic studies of resulting 'sky islands' have either failed to detect recent divergence or have found discordant evidence for ancient divergence. We test this paleoecological hypothesis for the region with intraspecific mitochondrial DNA and microsatellite data from sky-island populations of a sedentary bird, the Mexican jay (Aphelocoma ultramarina). We predicted that populations on different sky islands would share common, ancestral alleles that existed during the last glaciation, but that populations on each sky island, owing to their isolation, would contain unique variants of postglacial origin. We also predicted that divergence times estimated from corrected genetic distance and a coalescence model would post-date the last glacial maximum.

Results

Our results provide multiple independent lines of support for postglacial divergence, with the predicted pattern of shared and unique mitochondrial DNA haplotypes appearing in two independent sky-island archipelagos, and most estimates of divergence time based on corrected genetic distance post-dating the last glacial maximum. Likewise, an isolation model based on multilocus gene coalescence indicated postglacial divergence of five pairs of sky islands. In contrast to their similar recent histories, the two archipelagos had dissimilar historical patterns in that sky islands in Arizona showed evidence for older divergence, suggesting different responses to the last glaciation.

Conclusion

This study is one of the first to provide explicit support from genetic data for a postglacial divergence scenario predicted by one of the best paleoecological records in the world. Our results demonstrate that sky islands act as generators of genetic diversity at both recent and historical timescales and underscore the importance of thorough sampling and the use of loci with fast mutation rates to studies that test hypotheses concerning recent genetic divergence.     

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.     

Alpine biogeography of Parnassian butterflies during Quaternary climate cycles in North America

Growth of alpine glaciers during the Pleistocene had profound effects on montane landscapes in North America and the organisms now inhabiting alpine ecosystems. Biogeography of this region has often been viewed as a system of sky islands despite the fact that species richness patterns deviate from a strict island biogeographic model. One explanation is that alpine species are not in equilibrium because of late Quaternary geographic range shifts. Genetic data can provide evidence of nonequilibrium dynamics and the distributional shifts that occur during glaciation events in alpine landscapes. Using mitochondrial and nuclear sequence data, we examine the evolutionary history of butterflies in the Parnassius phoebus complex. We test explicit, alternative models of the biogeographic history of Parnassius smintheus and Parnassius behrii , including an equilibrium island model, ancestral radiation and fragmentation, an expanding alpine archipelago and an alpine archipelago refuge model. Our results support the alpine archipelago refuge model, in which alpine butterflies undergo population contraction during glacial climates followed by population expansion during interglacial phases. While butterflies can disperse between distant mountain ranges during glacial periods, gene flow is rare. We find evidence of recent connectivity between California and Colorado, population expansion events following deglaciation ∼20 000 years B.P., and small population sizes during the last glacial period. An analysis of lineage splitting suggests that morphological differences in P. smintheus and P. behrii are the result of late Pleistocene divergence (∼48 000 years B.P.) with limited gene flow. Our results demonstrate that spatially complex and nonequilibrium population dynamics influence alpine diversity patterns.     

CLIMATE‐DRIVEN DIVERSIFICATION AND PLEISTOCENE REFUGIA IN PHILIPPINE BIRDS: EVIDENCE FROM PHYLOGEOGRAPHIC STRUCTURE AND PALEOENVIRONMENTAL NICHE MODELING

Avian diversification in oceanic archipelagos is largely attributed to isolation across marine barriers. During glacial maxima, lowered sea levels resulted in repeated land connections between islands joined by shallow seas. Consequently, such islands are not expected to show endemism. However, if climate fluctuations simultaneously caused shifts in suitable environmental conditions, limiting populations to refugia, then occurrence on and dispersal across periodic land bridges are not tenable. To assess the degree to which paleoclimate barriers, rather than marine barriers, drove avian diversification in the Philippine Archipelago, we produced ecological niche models for current‐day, glacial maxima, and interglacial climate scenarios to infer potential Pleistocene distributions and paleoclimate barriers. We then tested marine and paleoclimate barriers for correspondence to geographic patterns of population divergence, inferred from DNA sequences from eight codistributed bird species. In all species, deep‐water channels corresponded to zones of genetic differentiation, but six species exhibited deeper divergence associated with a periodic land bridge in the southern Philippines. Ecological niche models for these species identified a common paleoclimate barrier that coincided with deep genetic structure among populations. Although dry land connections joined southern Philippine islands during low sea level stands, unfavorable environmental conditions limited populations within landmasses, resulting in long‐term isolation and genetic differentiation. These results highlight the complex nature of diversification in archipelagos: marine barriers, changes in connectivity due to sea level change, and climate‐induced refugia acted in concert to produce great species diversity and endemism in the Philippines.     

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.     

Pliocene forest dynamics as a primary driver of African bird speciation

Aim Montane tropics are areas of high endemism, and mechanisms driving this endemism have been receiving increasing attention at a global scale. A general trend is that climatic factors do not explain the species richness of species with small to medium‐sized geographic ranges, suggesting that geological and evolutionary processes must be considered. On the African continent, several hypotheses including both refugial and geographic uplift models have been advanced to explain avian speciation and diversity in the lowland forest and montane regions of central and eastern Africa; montane regions in particular are recognized as hotspots of vertebrate endemism. Here, we examine the possible role of these models in driving speciation in a clade of African forest robins. Location Africa. Methods We constructed the first robustly supported molecular phylogenetic hypothesis of forest robins. On this phylogeny, we reconstructed habitat‐based distributions and geographic distributions relative to the Albertine Rift. We also estimated the timing of lineage divergences via a molecular clock. Results Robust estimates of phylogenetic relationships and clock‐based divergences reject Miocene tectonic uplift and Pleistocene forest refugia as primary drivers of speciation in forest robins. Instead, our data suggest that most forest robin speciation took place in the Late Pliocene, from 3.2 to 2.2 Ma. Distributional patterns are complex, with the Albertine Rift region serving as a general east–west break across the group. Montane distributions are inferred to have evolved four times. Main conclusions Phylogenetic divergence dates coincide with a single period of lowland forest retraction in the late Pliocene, suggesting that most montane speciation resulted from the rapid isolation of populations in montane areas, rather than montane areas themselves being drivers of speciation. This conclusion provides additional evidence that Pliocene climate change was a major driver of speciation in broadly distributed African animal lineages. We further show that lowland forest robins are no older than their montane relatives, suggesting that lowland areas are not museums which house ‘ancient’ taxa; rather, for forest robins, montane areas should be viewed as living museums of a late Pliocene diversification event. A forest refugial pattern is operating in Africa, but it is not constrained to the Pleistocene.     

Islands within islands: two montane palaeo‐endemic birds impacted by recent anthropogenic fragmentation

Anthropogenic habitat fragmentation of species that live in naturally patchy metapopulations such as mountaintops or sky islands experiences two levels of patchiness. Effects of such multilevel patchiness on species have rarely been examined. Metapopulation theory suggests that patchy habitats could have varied impacts on persistence, dependent on differential migration. It is not known whether montane endemic species, evolutionarily adapted to natural patchiness, are able to disperse between anthropogenic fragments at similar spatial scales as natural patches. We investigated historic and contemporary gene flow between natural and anthropogenic patches across the distribution range of a Western Ghats sky‐island‐endemic bird species complex. Data from 14 microsatellites for 218 individuals detected major genetic structuring by deep valleys, including one hitherto undescribed barrier. As expected, we found strong effects of historic genetic differentiation across natural patches, but not across anthropogenic fragments. Contrastingly, contemporary differentiation (D_PS) was higher relative to historic differentiation (F_ST) in anthropogenic fragments, despite the species’ ability to historically traverse shallow valleys. Simulations of recent isolation resulted in high D_PS/F_ST values, confirming recent isolation in Western Ghats anthropogenic fragments and also suggesting that this ratio can be used to identifying recent fragmentation in the context of historic connectedness. We suggest that in this landscape, in addition to natural patchiness affecting population connectivity, anthropogenic fragmentation additionally impacts connectivity, making anthropogenic fragments akin to islands within natural islands of montane habitat, a pattern that may be recovered in other sky‐island systems.     

Coupling genetic and ecological-niche models to examine how past population distributions contribute to divergence

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

Lineage diversification and historical demography of a sky island salamander,Plethodon ouachitae,from the Interior Highlands

Sky islands provide ideal opportunities for understanding how climatic changes associated with Pleistocene glacial cycles influenced species distributions, genetic diversification, and demography. The salamander Plethodon ouachitae is largely restricted to high‐elevation, mesic forest on six major mountains in the Ouachita Mountains. Because these mountains are separated by more xeric, low‐elevation valleys, the salamanders appear to be isolated on sky islands where gene flow among populations on different mountains may be restricted. We used DNA sequence data along with ecological niche modelling and coalescent simulations to test several hypotheses related to diversifications in sky island habitats. Our results revealed that P. ouachitae is composed of seven well‐supported lineages structured across six major mountains. The species originated during the Late Pliocene, and lineage diversification occurred during the Middle Pleistocene in a stepping stone fashion with a cyclical pattern of dispersal to a new mountain followed by isolation and divergence. Diversification occurred primarily on an east–west axis, which is likely related to the east–west orientation of the Ouachita Mountains and the more favourable cooler and wetter environmental conditions on north slopes compared to south‐facing slopes and valleys. All non‐genealogical coalescent methods failed to detect significant population expansion in any lineages. Bayesian skyline plots showed relatively stable population sizes over time, but indicated a slight to moderate amount of population growth in all lineages starting approximately 10 000–12 000 years ago. Our results provide new insight into sky island diversifications from a previously unstudied region, and further demonstrate that climatic changes during the Pleistocene had profound effects on lineage diversification and demography, especially in species from environmentally sensitive habitats in montane regions.     

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.     

Sky islands as foci for divergence of fig trees and their pollinators in southwest China

The dynamics of populations and their divergence over time have shaped current levels of biodiversity and in the case of the “sky islands” of mountainous southwest (SW) China have resulted in an area of exceptional botanical diversity. Ficus tikoua is a prostrate fig tree subendemic to the area that displays unique intraspecific diversity, producing figs typical of different pollination modes in different parts of its range. By combining climate models, genetic variation in populations of the tree's obligate fig wasp pollinators and distributions of the different plant phenotypes, we examined how this unusual situation may have developed. We identified three genetically distinct groups of a single Ceratosolen pollinator species that have largely parapatric distributions. The complex topography of the region contributed to genetic divergence among the pollinators by facilitating geographical isolation and providing refugia. Migration along elevations in response to climate oscillations further enhanced genetic differentiation of the three pollinator groups. Their distributions loosely correspond to the distributions of the functionally significant morphological differences in the male figs of their host plants, but postglacial expansion of one group has not been matched by spread of its associated plant phenotype, possibly due to a major river barrier. The results highlight how interplay between the complex topography of the “sky island” complex and climate change has shaped intraspecies differentiation and relationships between the plant and its pollinator. Similar processes may explain the exceptional botanical diversity of SW China.     

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi‐locus and RNA‐seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ?a?i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ?a?i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).     

Cryptic diversity within the endemic prehensile‐tailed gecko Urocotyledon inexpectata across the Seychelles Islands: patterns of phylogeographical structure and isolation at the multilocus level

The deciphering of the process of genetic differentiation of species with insular distributions is relevant for biogeographical and conservation reasons. Despite their importance as old gondwanic islands and part of the western Indian Ocean biodiversity hotspot, little is known about the genetic structure of taxa from the Seychelles Islands. We have examined the patterns of structure and isolation within Urocotyledon inexpectata (Reptilia: Geckkonidae), an endemic species from this archipelago. Genetic diversity was screened from populations across the archipelago for both mitochondrial and nuclear genes. Gene genealogies and model‐based inference were used to explore patterns and timings of isolation between the main lineages. High levels of genetic diversity were found for the mitochondrial and some of the nuclear markers. This species harbours at least two highly differentiated lineages, exclusively distributed across the northern and southern groups of the islands. The main split between these was dated back to the Miocene–late Pliocene, but isolation events throughout the Pliocene and Pleistocene were also inferred. Migration between groups of islands was apparently nonexistent, except for one case. The low dispersal capabilities of this species, together with the intrinsic fragmented nature of its geographical distribution, seem to have resulted in highly structured populations, despite the cyclic periods of contact between the different island groups. These populations may currently represent more than one species, making U. inexpectata another example of a morphologically cryptic lineage with deep genetic divergence within gekkonids. The observed patterns suggest a hypothetical biogeographic scenario (of a main north–south phylogeographic break) for the Seychelles that can be further tested with the exploration of the phylogeographic structure of other Seychellois taxa. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 177–191.     

Spatial processes that maintain biodiversity in plant communities

The composition of communities of sessile organisms, and the change in species diversity with time, is a spatially explicit phenomenon. Three spatial factors clearly affect diversity: (1) the structure and heterogeneity of the landscape that limits species immigration and ultimate community size; (2) neighborhood interactions that determine colonization and extinction rates and influence residence times of local populations; and (3) disturbances that open spatially contiguous areas for recolonization by less abundant species. The importance of these three factors was first reviewed and then examined with a spatially explicit, multi-species model of plant dispersal, competition and establishment, with an assumption of neutrality (all species had equivalent life histories) that reduced the initial dimensionality of the problem. The simulations assumed that the probability of immigration was a linear function of mainland abundance and distance to islands, similar to the equilibrium theory of island biogeography and the unified neutral theory of biodiversity. The rate of increase in species richness was not constant across island sizes, declining as island area became very large. This pattern was explained by the spatial dynamics of colonization and establishment, a non-random process that cannot be explained by passive sampling alone. Simulations showed that population establishment depended critically on rare long-distance dispersal events while population persistence was achieved by the formation of aggregated species distributions that developed through restricted dispersal and local competitive interactions. Nevertheless, species richness always declined to a single species in the absence of disturbances, while up to 40 species could persist to 10,000 years when spatially dependent mortality was added. Further explorations with spatially explicit models will be required to fully appreciate the consequence of land use change and altered disturbance regimes on patterns of species distribution and the maintenance of diversity.