Evolutionary and ecological causes of the latitudinal diversity gradient in hylid frogs: treefrog trees unearth the roots of high tropical diversity

Why are there more species in the tropics than in temperate regions? In recent years, this long-standing question has been addressed primarily by seeking environmental correlates of diversity. But to understand the ultimate causes of diversity patterns, we must also examine the evolutionary and biogeographic processes that directly change species numbers (i.e., speciation, extinction, and dispersal). With this perspective, we dissect the latitudinal diversity gradient in hylid frogs. We reconstruct a phylogeny for 124 hylid species, estimate divergence times and diversification rates for major clades, reconstruct biogeographic changes, and use ecological niche modeling to identify climatic variables that potentially limit dispersal. We find that hylids originated in tropical South America and spread to temperate regions only recently (leaving limited time for speciation). There is a strong relationship between the species richness of each region and when that region was colonized but not between the latitudinal positions of clades and their rates of diversification. Temperature seasonality seemingly limits dispersal of many tropical clades into temperate regions and shows significant phylogenetic conservatism. Overall, our study illustrates how two general principles (niche conservatism and the time-for-speciation effect) may help explain the latitudinal diversity gradient as well as many other diversity patterns across taxa and regions.     

Multilocus tests of Pleistocene refugia and ancient divergence in a pair of Atlantic Forest antbirds (Myrmeciza)

The Atlantic Forest (AF) harbours one of the most diverse vertebrate faunas of the world, including 199 endemic species of birds. Understanding the evolutionary processes behind such diversity has become the focus of many recent, primarily single locus, phylogeographic studies. These studies suggest that isolation in forest refugia may have been a major mechanism promoting diversification, although there is also support for a role of riverine and geotectonic barriers, two sets of hypotheses that can best be tested with multilocus data. Here we combined multilocus data (one mtDNA marker and eight anonymous nuclear loci) from two species of parapatric antbirds, Myrmeciza loricata and M. squamosa, and Approximate Bayesian Computation to determine whether isolation in refugia explains current patterns of genetic variation and their status as independent evolutionary units. Patterns of population structure, differences in intraspecific levels of divergence and coalescent estimates of historical demography fit the predictions of a recently proposed model of refuge isolation in which climatic stability in the northern AF sustains higher diversity and demographic stability than in the southern AF. However, a pre‐Pleistocene divergence associated with their abutting range limits in a region of past tectonic activity also suggests a role for rivers or geotectonic barriers. Little or no gene flow between these species suggests the development of reproductive barriers or competitive exclusion. Our results suggests that limited marker sampling in recent AF studies may compromise estimates of divergence times and historical demography, and we discuss the effects of such sampling on this and other studies.     

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon-specific attributes.     

Intra- and interspecific morphological variation in the field mouse species Apodemus argenteus and A. speciosus in the Japanese archipelago: the role of insular isolation and biogeographic gradients

Two species of field mice, Apodemus argenteus and A. speciosus, occur in sympatry across the Japanese archipelago. The inter- and intraspecific patterns of morphological differentiation have been evaluated, using a Fourier analysis of the mandible outline. The relative importance of the effect of insular isolation and latitudinal climatic gradient on the size and shape of the two species was assessed by a comparison of the populations from the large island of Honshu and the surrounding small-island populations. The size variation in A. argenteus is correlated with the climatic gradient whilst the shape variation corresponds mainly to a random differentiation of the small-island populations from a Honshu-like basic morphological pattern. A. speciosus displays increased size on small islands, and its shape variation is related to both the climatic gradient and insularity. Finally, the two species are differentiated by both the size and shape of the mandible across the Japanese archipelago, suggesting that interspecific competition between both species is reduced via niche partitioning. Our results emphasize the importance of insular isolation on shape differentiation, but a part of the morphological differentiation is also related to the latitudinal climatic gradient. Isolation on small islands could have favoured such a response to environmental factors by lowering the gene flow that prevents almost any significant differentiation within Honshu populations.     

Evolutionary diversification of the lizard genus Bassiana (Scincidae) across Southern Australia

Background

Relatively recent (Plio-Pleistocene) climatic variations had strong impacts on the fauna and flora of temperate-zone North America and Europe; genetic analyses suggest that many lineages were restricted to unglaciated refuges during this time, and have expanded their ranges since then. Temperate-zone Australia experienced less severe glaciation, suggesting that patterns of genetic structure among species may reflect older (aridity-driven) divergence events rather than Plio-Pleistocene (thermally-mediated) divergences. The lizard genus Bassiana (Squamata, Scincidae) contains three species that occur across a wide area of southern Australia (including Tasmania), rendering them ideally-suited to studies on the impact of past climatic fluctuations.

Methodology/Principal Findings

We performed molecular phylogenetic and dating analyses using two partial mitochondrial genes (ND2 and ND4) of 97 samples of Bassiana spp. Our results reveal a pattern of diversification beginning in the Middle Miocene, with intraspecific diversification arising from 5.7 to 1.7 million years ago in the Upper Miocene-Lower Pleistocene.

Conclusions/Significance

In contrast to the temperate-zone Northern Hemisphere biota, patterns of evolutionary diversification within southern Australian taxa appear to reflect geologically ancient events, mostly relating to east-west discontinuities imposed by aridity rather than (as is the case in Europe and North America) relatively recent recolonisation of northern regions from unglaciated refugia to the south.     

武夷山典型植被类型土壤动物群落的结构特征

选择武夷山不同海拔高度的4个典型植被群落(EVB,COF,DWF和ALM),对其土壤动物群落的结构和变化进行研究,结果表明,植被沿海拔梯度变化对土壤动物类群数、密度和多样性产生了显著的影响。不同海拔高度样地中,密度-类群DG多样性指数和土壤动物数量沿海拔高度上升呈显著递减趋势;数量占优势的土壤动物类群膜翅目(F=6.35,P0.001)、蜱螨目(F=9.25,P0.01)、弹尾目(F=5.93,P0.05)和鞘翅目(F=4.42,P0.05)的密度在不同海拔样地之间存在显著的差异。除高山草甸外,土壤动物群落在土层中的垂直分布呈明显"表聚现象"。但高山草甸的弹尾目、膜翅目、双翅目(幼虫)和等足目等类群较多地聚集在10-25cm土层。不同海拔高度土壤动物的类群总数和总个体数量表现为2种季节变化模式:(1)常绿阔叶林和针叶林秋季大量发生;(2)亚高山矮林和高山草甸夏季大量发生且数量接近。     

Barriers,rather than refugia,underlie the origin of diversity in toads endemic to the Brazilian Atlantic Forest

In this study, we investigated the relative contribution of geographic barriers and Pleistocene refuges in the diversification of the Rhinella crucifer species complex, a group of endemic toads with a widespread distribution in the Brazilian Atlantic Forest (AF). We used intensive sampling and multilocus DNA sequence data to compare nucleotide diversity between refuge and nonrefuge areas, investigate regional demographic patterns, estimate demographic parameters related to genetic breaks and test refuge versus barrier scenarios of diversification using approximate Bayesian computation. We did not find higher levels of genetic diversity in putative refuge areas, either at regional or biome scale. Rather, the demographic history of the species complex supports regional differences with moderate population growth in the north and central regions and stability in southern AF. Genetic breaks were dated to the Plio–Pleistocene; however, our analyses rejected the role of refuges in creating a northern and central divergence, supporting a recent colonization scenario at a smaller scale within the central AF. Overall, our data rule out massive climatically driven fragmentation and large‐scale recolonization events for populations across the biome. We confirmed the importance of geographic barriers in creating main divergences and underscored the importance of searching for cryptic discontinuities in the landscape. Comparison of our results with those of other AF taxa indicates organismal specific responses to moderate shifts in habitat and that multiple refuges may constitute a more realistic model for diversification of Atlantic Forest biota.     

Molecular substitution rate increases with latitude in butterflies: evidence for a trans‐glacial latitudinal layering of populations?

A well‐documented consequence of repeated global ice ages is the negative relationship between latitude and intraspecific genetic diversity. However, little is known about additional effects of such major climatic events on population genetic structure. Here we studied the phylogeographic structure of five lycaenid butterfly species with varied ecological adaptations, sampled across a latitudinal gradient in the Holarctic region. We found a positive correlation between latitude and substitution rate of mitochondrial DNA sequences in all species investigated. We propose that this result is the signal of increased genetic drift and founder effects during post‐glacial recolonization of northern populations. Given that phylogenetic branch length is the result of a cumulative process over evolutionary time, we hypothesize that a latitudinal layering of populations has generally been maintained during repeated cycles of glaciation, possibly due to a neutral spatial effect and/or local adaptive advantage. This trans‐glacial latitudinal layering could be viewed as a particular case of the more general phenomenon of intraspecific structuring that is created and maintained in a fluctuating environmental gradient.     

Latitudinal‐diversity gradients can be shaped by biotic processes: new insights from an eco‐evolutionary model

The processes involved in shaping latitudinal‐diversity gradients (LDGs) have been a longstanding source of debate and research. Climatic, historical and evolutionary factors have all been shown to contribute to the formation of LDGs. However, meta‐analyses have shown that different clades have LDG slopes that may vary in more than one order of magnitude. Such large variation cannot be explained solely by climatic or historical factors (e.g. difference in surface area between temperate and tropical zones) given that all clades within a geographic region are subject to the same conditions. Therefore, biotic processes intrinsic to each taxonomic group could be relevant in explaining rate differences in diversity decline across latitudinal gradients among groups. In this study, we developed a model simulating multiple competing species subjected (or not) to a demographic Allee effect. We simulated the range expansion of these species across an environmental gradient to show how these two overlooked factors (competition and Allee effects) are capable of modulating LDGs. Allee effects resulted in a steeper LDG given a higher probability of local extinction and lower colonization capacity compared to species without Allee effects. Likewise, stronger competition also led to a steeper decline in species diversity compared to scenarios with weaker species antagonistic interactions. This pattern occurred mostly due to the strength of priority effects, wherein scenarios with strong competition, species that dispersed earlier in the landscape were able to secure many patches whereas late‐arriving species were progressively precluded from expanding their ranges. Overall, our results suggest that the effect of biotic processes in shaping macroecological patterns could be more important than it is currently appreciated.     

Molecular phylogenetics and biogeography of the Neocellia Series of Anopheles mosquitoes in the Oriental Region

Molecular studies of population divergence and speciation across the Oriental Region are sparse, despite the region’s high biodiversity and extensive Pliocene and Pleistocene environmental change. A molecular phylogenetic study of the Neocellia Series of Anopheles mosquitoes was undertaken to identify patterns of diversification across the Oriental Region and to infer the role of Pleistocene and Pliocene climatic change. A robust phylogeny was constructed using CO2 and ND5 mitochondrial genes and ITS2 and D3 nuclear ribosomal markers. Bayesian analysis of mitochondrial genes was used to date divergence events. The repeated contraction and expansion of forest habitat resulting from Pleistocene climatic fluctuations appears to have had a substantial impact on intraspecific diversification, but has not driven speciation within this group. Primarily early to mid Pliocene speciation was detected within the Annularis Group, whereas speciation within the Maculatus and Jamesii Groups occurred during the mid and late Pliocene. Both allopatric divergence driven by late Pliocene environmental changes and ecological adaptation, involving altitudinal replacement and seasonality, are likely to have influenced speciation in the Maculatus Group.     

Metabolic Rate and Climatic Fluctuations Shape Continental Wide Pattern of Genetic Divergence and Biodiversity in Fishes

Taxonomically exhaustive and continent wide patterns of genetic divergence within and between species have rarely been described and the underlying evolutionary causes shaping biodiversity distribution remain contentious. Here, we show that geographic patterns of intraspecific and interspecific genetic divergence among nearly all of the North American freshwater fish species (>750 species) support a dual role involving both the late Pliocene-Pleistocene climatic fluctuations and metabolic rate in determining latitudinal gradients of genetic divergence and very likely influencing speciation rates. Results indicate that the recurrent glacial cycles caused global reduction in intraspecific diversity, interspecific genetic divergence, and species richness at higher latitudes. At the opposite, longer geographic isolation, higher metabolic rate increasing substitution rate and possibly the rapid accumulation of genetic incompatibilities, led to an increasing biodiversity towards lower latitudes. This indicates that both intrinsic and extrinsic factors similarly affect micro and macro evolutionary processes shaping global patterns of biodiversity distribution. These results also indicate that factors favouring allopatric speciation are the main drivers underlying the diversification of North American freshwater fishes.     

Natural selection drives altitudinal divergence at the albumin locus in deer mice, Peromyscus maniculatus

In populations that are distributed across steep environmental gradients, the potential for local adaptation is largely determined by the spatial scale of fitness variation relative to dispersal distance. Since altitudinal gradients are generally characterized by dramatic ecological transitions over relatively short linear distances, adaptive divergence across such gradients will typically require especially strong selection to counterbalance the homogenizing effect of gene flow. Here we report the results of a study that was designed to test for evidence of adaptive divergence across an altitudinal gradient in a natural population of deer mice, Peromyscus maniculatus. We conducted a multilocus survey of allozyme variation across a steep altitudinal gradient in the southern Rocky Mountains that spanned several distinct biomes, from prairie grassland to alpine tundra. As a control for the effects of altitude, we also surveyed the same loci in mice sampled along a latitudinal transect through the prairie grassland that ran perpendicular to the east-west altitudinal transect. We used a coalescent-based simulation model to identify loci that deviated from neutral expectations, and we then assessed whether locus-specific patterns of variation were nonrandom with respect to altitude. Results indicated that the albumin locus (Alb) reflects a history of diversifying selection across the altitudinal gradient. This conclusion is supported by two main lines of evidence: (1) Alb was characterized by levels of divergence across the altitudinal transect that exceeded neutral expectations in two consecutive years of sampling (in contrast to the spatial pattern of variation across the latitudinal transect), and (2) levels of divergence at the Alb locus exhibited a positive association with altitudinal distance in both years (in contrast to the pattern observed at unlinked loci). We conclude that clinal variation at the Alb locus reflects a balance between gene flow and diversifying selection that results from elevational changes in fitness rankings among alternative genotypes.     

Adaptive latitudinal trends in the mandible shape of Apodemus wood mice

Aim Size and shape of the mandible are investigated across the latitudinal range of the European wood mouse (Apodemus sylvaticus), in order to address the relative importance of genetic structure, insularity, and geographical gradient in patterning morphological variation. Results are compared with those on two Asiatic species of wood mice, A. argenteus and A. speciosus. Location The European wood mouse is sampled by a set of trapping localities including both, islands and mainland populations, as well as the four genetic groups identified in previous studies. The localities cover a latitudinal gradient from 55° N to 36° N. Methods Different Fourier methods are applied to the outlines of mandibles and their results compared in the case of A. sylvaticus. All provide similar results and allow a quantification of the size and shape variations across the geographical range of the European wood mouse. Using the method allowing for the best reduction of the informative data set, a comparison of the European wood mouse with the two Asiatic species was performed. Results Within the European wood mouse A. sylvaticus, a strong latitudinal gradient in mandible shape overrides the influence of insularity and genetic structure. Yet, random morphological divergence in insular conditions can be identified as a secondary process of shape differentiation. Size displays no obvious pattern of variation, neither with insularity or latitude. A comparison with two other species of wood mice suggests that a similar latitudinal gradient in mandible shape exists in different species, mandibles being flatter in the north and wider in the south. Main conclusion The latitudinal gradient in mandible shape observed in the three species of wood mice is interpreted as an intraspecific adaptive response to gradual changes in feeding behaviour.     

Morphological disparity opposes latitudinal diversity gradient in lacertid lizards

While global variation in taxonomic diversity is strongly linked to latitude, the extent to which morphological disparity follows geographical gradients is less well known. We estimated patterns of lineage diversification, morphological disparity and rates of phenotypic evolution in the Old World lizard family Lacertidae, which displays a nearly inverse latitudinal diversity gradient with decreasing species richness towards the tropics. We found that lacertids exhibit relatively constant rates of lineage accumulation over time, although the majority of morphological variation appears to have originated during recent divergence events, resulting in increased partitioning of disparity within subclades. Among subclades, tropical arboreal taxa exhibited the fastest rates of shape change while temperate European taxa were the slowest, resulting in an inverse relationship between latitudinal diversity and rates of phenotypic evolution. This pattern demonstrates a compelling counterexample to the ecological opportunity theory of diversification, suggesting an uncoupling of the processes generating species diversity and morphological differentiation across spatial scales.     

Deep Phylogeographic Structure and Environmental Differentiation in the Carnivorous Plant Sarracenia alata

We collected ～29 kb of sequence data using Roche 454 pyrosequencing in order to estimate the timing and pattern of diversification in the carnivorous pitcher plant Sarracenia alata. Utilizing modified protocols for reduced representation library construction, we generated sequence data from 86 individuals across 10 populations from throughout the range of the species. We identified 76 high-quality and high-coverage loci (containing over 500 SNPs) using the bioinformatics pipeline PRGmatic. Results from a Bayesian clustering analysis indicate that populations are highly structured, and are similar in pattern to the topology of a population tree estimated using *BEAST. The pattern of diversification within Sarracenia alata implies that riverine barriers are the primary factor promoting population diversification, with divergence across the Mississippi River occurring more than 60,000 generations before present. Further, significant patterns of niche divergence and the identification of several outlier loci suggest that selection may contribute to population divergence. Our results demonstrate the feasibility of using next-generation sequencing to investigate intraspecific genetic variation in nonmodel species.     

EVIDENCE FOR CLIMATE‐DRIVEN DIVERSIFICATION? A CAUTION FOR INTERPRETING ABC INFERENCES OF SIMULTANEOUS HISTORICAL EVENTS

Approximate Bayesian computation (ABC) is rapidly gaining popularity in population genetics. One example, msBayes , infers the distribution of divergence times among pairs of taxa, allowing phylogeographers to test hypotheses about historical causes of diversification in co‐distributed groups of organisms. Using msBayes , we infer the distribution of divergence times among 22 pairs of populations of vertebrates distributed across the Philippine Archipelago. Our objective was to test whether sea‐level oscillations during the Pleistocene caused diversification across the islands. To guide interpretation of our results, we perform a suite of simulation‐based power analyses. Our empirical results strongly support a recent simultaneous divergence event for all 22 taxon pairs, consistent with the prediction of the Pleistocene‐driven diversification hypothesis. However, our empirical estimates are sensitive to changes in prior distributions, and our simulations reveal low power of the method to detect random variation in divergence times and bias toward supporting clustered divergences. Our results demonstrate that analyses exploring power and prior sensitivity should accompany ABC model selection inferences. The problems we identify are potentially mitigable with uniform priors over divergence models (rather than classes of models) and more flexible prior distributions on demographic and divergence‐time parameters.     

Environmental harshness is positively correlated with intraspecific divergence in mammals and birds

Life on Earth is conspicuously more diverse in the tropics. Although this intriguing geographical pattern has been linked to many biotic and abiotic factors, their relative importance and potential interactions are still poorly understood. The way in which latitudinal changes in ecological conditions influence evolutionary processes is particularly controversial, as there is evidence for both a positive and a negative latitudinal gradient in speciation rates. Here, we identify and address some methodological issues (how patterns are analysed and how latitude is quantified) that could lead to such conflicting results. To address these issues, we assemble a comprehensive data set of the environmental correlates of latitude (including climate, net primary productivity and habitat heterogeneity) and combine it with biological, historical and molecular data to explore global patterns in recent divergence events (subspeciation). Surprisingly, we find that the harsher conditions that typify temperate habitats (lower primary productivity, decreased rainfall and more variable and unpredictable temperatures) are positively correlated with greater subspecies richness in terrestrial mammals and birds. Thus, our findings indicate that intraspecific divergence is greater in regions with lower biodiversity, a pattern that is robust to both sampling variation and latitudinal biases in taxonomic knowledge. We discuss possible causal mechanisms for the link between environmental harshness and subspecies richness (faster rates of evolution, greater likelihood of range discontinuities and more opportunities for divergence) and conclude that this pattern supports recent indications that latitudinal gradients of diversity are maintained by simultaneously higher potentials for both speciation and extinction in temperate than tropical regions.     

Influence of drainage divides versus arid corridors on genetic structure and demography of a widespread freshwater turtle,Emydura macquarii krefftii,from Australia

The influence of Pleistocene climatic cycles on Southern Hemisphere biotas is not yet well understood. Australia's eastern coastal margin provides an ideal setting for examining the relative influence of landscape development, sea level fluctuation, and cyclic climatic aridity on the evolution of freshwater biodiversity. We examined the impact of climatic oscillations and physical biogeographic barriers on the evolutionary history of the wide‐ranging Krefft's river turtle (Emydura macquarii krefftii), using range‐wide sampling (649 individuals representing 18 locations across 11 drainages) and analysis of mitochondrial sequences (~1.3‐kb control region and ND4) and nuclear microsatellites (12 polymorphic loci). A range of phylogeographic (haplotype networks, molecular dating), demographic (neutrality tests, mismatch distributions), and population genetic analyses (pairwise F_ST, analysis of molecular variance, Bayesian clustering analysis) were implemented to differentiate between competing demographic (local persistence vs. range expansion) and biogeographic (arid corridor vs. drainage divide) scenarios. Genetic data reveal population genetic structure in Krefft's river turtles primarily reflects isolation across drainage divides. Striking north‐south regional divergence (2.2% ND4 p‐distance; c. 4.73 Ma, 95% higher posterior density (HPD) 2.08–8.16 Ma) was consistent with long‐term isolation across a major drainage divide, not an adjacent arid corridor. Ancient divergence among regional lineages implies persistence of northern Krefft's populations despite the recurrent phases of severe local aridity, but with very low contemporary genetic diversity. Stable demography and high levels of genetic diversity are inferred for southern populations, where aridity was less extreme. Range‐wide genetic structure in Krefft's river turtles reflects contemporary and historical drainage architecture, although regional differences in the extent of Plio–Pleistocene climatic aridity may be reflected in current levels of genetic diversity.     

Idiosyncratic responses to climate‐driven forest fragmentation and marine incursions in reed frogs from Central Africa and the Gulf of Guinea Islands

Organismal traits interact with environmental variation to mediate how species respond to shared landscapes. Thus, differences in traits related to dispersal ability or physiological tolerance may result in phylogeographic discordance among co‐distributed taxa, even when they are responding to common barriers. We quantified climatic suitability and stability, and phylogeographic divergence within three reed frog species complexes across the Guineo‐Congolian forests and Gulf of Guinea archipelago of Central Africa to investigate how they responded to a shared climatic and geological history. Our species‐specific estimates of climatic suitability through time are consistent with temporal and spatial heterogeneity in diversification among the species complexes, indicating that differences in ecological breadth may partly explain these idiosyncratic patterns. Likewise, we demonstrated that fluctuating sea levels periodically exposed a land bridge connecting Bioko Island with the mainland Guineo‐Congolian forest and that habitats across the exposed land bridge likely enabled dispersal in some species, but not in others. We did not find evidence that rivers are biogeographic barriers across any of the species complexes. Despite marked differences in the geographic extent of stable climates and temporal estimates of divergence among the species complexes, we recovered a shared pattern of intermittent climatic suitability with recent population connectivity and demographic expansion across the Congo Basin. This pattern supports the hypothesis that genetic exchange across the Congo Basin during humid periods, followed by vicariance during arid periods, has shaped regional diversity. Finally, we identified many distinct lineages among our focal taxa, some of which may reflect incipient or unrecognized species.