Testing the role of ecology and life history in structuring genetic variation across a landscape: a trait‐based phylogeographic approach

Hypotheses to explain phylogeographic structure traditionally invoke geographic features, but often fail to provide a general explanation for spatial patterns of genetic variation. Organisms' intrinsic characteristics might play more important roles than landscape features in determining phylogeographic structure. We developed a novel comparative approach to explore the role of ecological and life‐history variables in determining spatial genetic variation and tested it on frog communities in Panama. We quantified spatial genetic variation within 31 anuran species based on mitochondrial DNA sequences, for which hierarchical approximate Bayesian computation analyses rejected simultaneous divergence over a common landscape. Regressing ecological variables, on genetic divergence allowed us to test the importance of individual variables revealing that body size, current landscape resistance, geographic range, biogeographic origin and reproductive mode were significant predictors of spatial genetic variation. Our results support the idea that phylogeographic structure represents the outcome of an interaction between organisms and their environment, and suggest a conceptual integration we refer to as trait‐based phylogeography.     

Environmental Variables Explain Genetic Structure in a Beetle-Associated Nematode

The distribution of a species is a complex expression of its ecological and evolutionary history and integrating population genetic, environmental, and ecological data can provide new insights into the effects of the environment on the population structure of species. Previous work demonstrated strong patterns of genetic differentiation in natural populations of the hermaphroditic nematode Pristionchus pacificus in its La Réunion Island habitat, but gave no clear understanding of the role of the environment in structuring this variation. Here, we present what is to our knowledge the first study to statistically evaluate the role of the environment in shaping the structure and distribution of nematode populations. We test the hypothesis that genetic structure in P. pacificus is influenced by environmental variables, by combining population genetic analyses of microsatellite data from 18 populations and 370 strains, with multivariate statistics on environmental data, and species distribution modelling. We assess and quantify the relative importance of environmental factors (geographic distance, altitude, temperature, precipitation, and beetle host) on genetic variation among populations. Despite the fact that geographic populations of P. pacificus comprise vast genetic diversity sourced from multiple ancestral lineages, we find strong evidence for local associations between environment and genetic variation. Further, we show that significantly more genetic variation in P. pacificus populations is explained by environmental variation than by geographic distances. This supports a strong role for environmental heterogeneity vs. genetic drift in the divergence of populations, which we suggest may be influenced by adaptive forces.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Environmental Heterogeneity Explains the Genetic Structure of Continental and Mediterranean Populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.     

Geographic variation in genetic and demographic performance: new insights from an old biogeographical paradigm

The ‘centre–periphery hypothesis’ (CPH) is a long‐standing postulate in ecology that states that genetic variation and demographic performance of a species decrease from the centre to the edge of its geographic range. This hypothesis is based on an assumed concordance between geographical peripherality and ecological marginality such that environmental conditions become harsher towards the limits of a species range. In this way, the CPH sets the stage for understanding the causes of distribution limits. To date, no study has examined conjointly the consistency of these postulates. In an extensive literature review we discuss the birth and development of the CPH and provide an assessment of the CPH by reviewing 248 empirical studies in the context of three main themes. First, a decrease in species occurrence towards their range limits was observed in 81% of studies, while only 51% demonstrated reduced abundance of individuals. A decline in genetic variation, increased differentiation among populations and higher rates of inbreeding were demonstrated by roughly one in two studies (47, 45 and 48%, respectively). However, demographic rates, size and population performance less often followed CPH expectations (20–30% of studies). We highlight the impact of important methodological, taxonomic, and biogeographical biases on such validation rates. Second, we found that geographic and ecological marginality gradients are not systematically concordant, which casts doubt on the reliability of a main assumption of the CPH. Finally, we attempt to disentangle the relative contribution of geographical, ecological and historical processes on the spatial distribution of genetic and demographic parameters. While ecological marginality gradients explain variation in species' demographic performance better than geographic gradients, contemporary and historical factors may contribute interactively to spatial patterns of genetic variation. We thereby propose a framework that integrates species' ecological niche characteristics together with current and past range structure to investigate spatial patterns of genetic and demographic variation across species ranges.     

DOES NICHE CONSERVATISM PROMOTE SPECIATION? A CASE STUDY IN NORTH AMERICAN SALAMANDERS

Recent speciation research has generally focused on how lineages that originate in allopatry evolve intrinsic reproductive isolation, or how ecological divergence promotes nonallopatric speciation. However, the ecological basis of allopatric isolation, which underlies the most common geographic mode of speciation, remains poorly understood and largely unstudied. Here, we explore the ecological and evolutionary factors that promote speciation in Desmognathus and Plethodon salamanders from temperate eastern North America. Based on published molecular phylogenetic estimates and the degree of geographic range overlap among extant species, we find strong evidence for a role for geographic isolation in speciation. We then examine the relationship between climatic variation and speciation in 16 sister-taxon pairs using geographic information system maps of climatic variables, new methods for modeling species' potential geographic distributions, and data on geographic patterns of genetic variation. In contrast to recent studies in tropical montane regions, we found no evidence for parapatric speciation along climatic gradients. Instead, many montane sister taxa in the Appalachian Highlands inhabit similar climatic niches and seemingly are allopatric because they are unable to tolerate the climatic conditions in the intervening lowlands. This temporal and spatial-ecological pattern suggests that niche conservatism, rather than niche divergence, plays the primary role in promoting allopatric speciation and montane endemism in this species-rich group of vertebrates. Our results demonstrate that even the relatively subtle climatic differences between montane and lowland habitats in eastern North America may play a key role in the origin of new species.     

Assortative mating in poison-dart frogs based on an ecologically important trait

The origin of new species can be influenced by both deterministic and stochastic factors. Mate choice and natural selection may be important deterministic causes of speciation (as opposed to the essentially stochastic factors of geographic isolation and genetic drift). Theoretical models predict that speciation is more likely when mate choice depends on an ecologically important trait that is subject to divergent natural selection, although many authors have considered such mating/ecology pleiotropy, or "magic-traits" to be unlikely. However, phenotypic signals are important in both mate choice and ecological processes such as avoiding predation. In chemically defended species, it may be that the phenotypic characteristics influencing mate choice are the same signals being used to transmit a warning to potential predators, although few studies have demonstrated this in wild populations. We tested for assortative mating between two color morphs of the Strawberry Poison-Dart Frog, Dendrobates pumilio, a group with striking geographic variation in aposematic color patterns. We found that females significantly prefer individuals of their own morph under two different light treatments, indicating strong assortative mating based on multiple coloration cues that are also important ecological signals. This study provides a rare example of one phenotypic trait affecting both ecological viability and nonrandom mating, indicating that mating/ecology pleiotropy is plausible in wild populations, particularly for organisms that are aposematically colored and visually orienting.     

Structural equation modelling analysis of evolutionary and ecological patterns in Australian Banksia

Evolutionary history of species, their geographic ranges, ecological ranges, genetic diversity, and resistance to pathogen infection, have been viewed as being mutually linked through a complex network of interactions. Previous studies have described simple correlations between pairs of these factors, while rarely separated the direct effects among multiple interacting factors. This study was to separate the effect of multiple interacting factors, to reveal the strength of the interactions among these factors, and to explore the mechanisms underlying the ecological and evolutionary processes shaping the geographic range, genetic diversity and fitness of species. I assembled comparative data on evolutionary history, geographic range, ecological range, genetic diversity, and resistance to pathogen infection for thirteen Banksia species from Australia. I used structural equation modelling to test multivariate hypotheses involving evolutionary history, geographic range, genetic diversity and fitness. Key results are: (1) Species with longer evolutionary times tend to occupy larger geographic ranges; (2) higher genetic diversity is directly associated with longer flowering duration in Banksia; and (3) species with higher genotypic diversity have higher level of resistance to infection caused by the pathogen Phytophthora cinnamomi, whereas heterozygosity has the opposite relationship with capacity of resistance to the infections caused by this pathogen. These results revealed a mutually linked and complex network of interactions among gene, species, environment and pathogen in evolutionary and ecological scales. These findings also have great practical significance and help to provide preemptive management options in pathogen control.     

Phylogenetic analyses reveal extensive cryptic speciation and host specialization in an economically important mite taxon

The wheat curl mite (WCM) is a major pest in cereal crops around the world and the vector of at least four known pathogens capable of reducing yields in crops such as wheat, corn, barley, oats, millet and rye. Current taxonomy recognizes WCM as a single species, Aceria tosichella; however, recent genetic, physiological and ecological studies have shown that WCM is likely to be a species complex. In this study we assessed genetic variation and phylogenetic relationships among WCM from four continents and a wide range of host plants using DNA sequence data from one mitochondrial gene, one nuclear gene and a single nuclear intergenic spacer region. Phylogenetic analyses revealed 11 unique mite lineages associated with specific plant hosts including wheat and barley. Host associations were consistent across continents, often with a single haplotype dominating a host plant regardless of geographic origin. The genetic and ecological differences identified in this study support the notion that WCM is a species complex in need of major taxonomic revision. These findings have implications for control of WCM globally, particularly within the context of identifying plants that form ‘green bridge’ refuges, assessing disease transmission risk, and identifying resistance in cereal genotypes to WCM and associated pathogens.     

RELATIONSHIP OF GENETIC DIVERSITY AND NICHE CENTRALITY: A SURVEY AND ANALYSIS

The distribution of genetic diversity within and among populations in relation to species’ geographic ranges is important to understanding processes of evolution, speciation, and biogeography. One hypothesis predicts that natural populations at geographic range margins will have lower genetic diversity relative to those located centrally in species’ distributions owing to a link between geographic and environmental marginality; alternatively, genetic variation may be unrelated with geographic marginality via decoupling of geographic and environmental marginality. We investigate the predictivity of geographic patterns of genetic variation based on geographic and environmental marginality using published genetic diversity data for 40 species (insects, plants, birds, mammals, worms). Only about half of species showed positive relationships between geographic and environmental marginality. Three analyses (sign test, multiple linear regression, and meta‐analysis of correlation effect sizes) showed a negative relationship between genetic diversity and distance to environmental niche centroid, but no consistent relationship of genetic diversity with distance to geographic range center.     

Genetic and phylogeographic structures of the symbiotic fungus Tuber magnatum

The quality and market price of truffles vary with the species and, traditionally, the place of origin. The premium species Tuber magnatum produces white truffles and has a patchy distribution restricted to Italy and some Balkan areas. We used polymorphic microsatellites to evaluate 316 specimens grouped into 26 populations sampled across the species' geographic range to determine if natural populations of T. magnatum are genetically differentiated. We found that the southernmost and the northwesternmost populations were significantly differentiated from the rest of the populations. The simple sequence repeat data also could be used to make inferences about the postglacial T. magnatum expansion pattern. This study is the first to identify a genetic and phylogeographic structure in T. magnatum. The presence of a genetic structure can be of practical interest in tracing truffle populations according to their geographic origin for marketing strategies. Evidence for extensive outcrossing in field populations of T. magnatum also is provided for the first time.     

Speciation along a latitudinal gradient: The origin of the Neotropical cycad sister pair Dioon sonorense–D. vovidesii (Zamiaceae)

Latitude is correlated with environmental components that determine the distribution of biodiversity. In combination with geographic factors, latitude‐associated environmental variables are expected to influence speciation, but empirical evidence on how those factors interplay is scarce. We evaluated the genetic and environmental variation among populations in the pair of sister species Dioon sonorense–D. vovidesii, two cycads distributed along a latitudinal environmental gradient in northwestern Mexico, to reveal their demographic histories and the environmental factors involved in their divergence. Using genome‐wide loci data, we determined the species delimitation, estimated the gene flow, and compared multiple demographic scenarios of divergence. Also, we estimated the variation of climatic variables among populations and used ecological niche models to test niche overlap between species. The effect of geographic and environmental variables on the genetic variation among populations was evaluated using linear models. Our results suggest the existence of a widespread ancestral population that split into the two species ~829 ky ago. The geographic delimitation along the environmental gradient occurs in the absence of major geographic barriers, near the 28th parallel north, where a zonation of environmental seasonality exists. The northern species, D. vovidesii, occurs in more seasonal environments but retains the same niche of the southern species, D. sonorense. The genetic variation throughout populations cannot be solely explained by stochastic processes; the latitudinal‐associated seasonality has been an additive factor that strengthened the species divergence. This study represents an example of how speciation can be achieved by the effect of the latitude‐associated factors on the genetic divergence among populations.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

SPECIATION AND POPULATION GENETIC STRUCTURE IN TROPICAL PACIFIC SEA URCHINS

Unlike populations of many terrestrial species, marine populations often are not separated by obvious, permanent barriers to gene flow. When species have high dispersal potential and few barriers to gene flow, allopatric divergence is slow. Nevertheless, many marine species are of recent origin, even in taxa with high dispersal potential. To understand the relationship between genetic structure and recent species formation in high dispersal taxa, we examined population genetic structure among four species of sea urchins in the tropical Indo-West Pacific that have speciated within the past one to three million years. Despite high potential for gene flow, mtDNA sequence variation among 200 individuals of four species in the urchin genus Echinometra shows a signal of strong geographic effects. These effects include (1) substantial population heterogeneity; (2) lower genetic variation in peripheral populations; and (3) isolation by distance. These geographic patterns are especially strong across scales of 5000-10,000 km, and are weaker over scales of 2500-5000 km. As a result, strong geographic patterns would not have been readily visible except over the wide expanse of the tropical Pacific. Surface currents in the Pacific do not explain patterns of gene flow any better than do patterns of simple spatial proximity. Finally, populations of each species tend to group into large mtDNA regions with similar mtDNA haplotypes, but these regional boundaries are not concordant in different species. These results show that all four species have accumulated mtDNA differences over similar spatial and temporal scales but that the precise geographic pattern of genetic differentiation varies for each species. These geographic patterns appear much less deterministic than in other well-known coastal marine systems and may be driven by chance and historical accident.     

Distribution, ecology, life history, genetic variation, and risk of extinction of nonhuman primates from Costa Rica

We examined the association between geographic distribution, ecological traits, life history, genetic diversity, and risk of extinction in nonhuman primate species from Costa Rica. All of the current nonhuman primate species from Costa Rica are included in the study; spider monkeys (Ateles geoffroyi), howling monkeys (Alouatta palliata), capuchins (Cebus capucinus), and squirrel monkeys (Saimiri oerstedii). Geographic distribution was characterized accessing existing databases. Data on ecology and life history traits were obtained through a literature review. Genetic diversity was characterized using isozyme electrophoresis. Risk of extinction was assessed from the literature. We found that species differed in all these traits. Using these data, we conducted a Pearson correlation between risk of extinction and ecological and life history traits, and genetic variation, for widely distributed species. We found a negative association between risk of extinction and population birth and growth rates; indicating that slower reproducing species had a greater risk of extinction. We found a positive association between genetic variation and risk of extinction; i.e., species showing higher genetic variation had a greater risk of extinction. The relevance of these traits for conservation efforts is discussed.     

Genetic response to climatic change: insights from ancient DNA and phylochronology

Understanding how climatic change impacts biological diversity is critical to conservation. Yet despite demonstrated effects of climatic perturbation on geographic ranges and population persistence, surprisingly little is known of the genetic response of species. Even less is known over ecologically long time scales pertinent to understanding the interplay between microevolution and environmental change. Here, we present a study of population variation by directly tracking genetic change and population size in two geographically widespread mammal species (Microtus montanus and Thomomys talpoides) during late-Holocene climatic change. We use ancient DNA to compare two independent estimates of population size (ecological and genetic) and corroborate our results with gene diversity and serial coalescent simulations. Our data and analyses indicate that, with population size decreasing at times of climatic change, some species will exhibit declining gene diversity as expected from simple population genetic models, whereas others will not. While our results could be consistent with selection, independent lines of evidence implicate differences in gene flow, which depends on the life history strategy of species.     

Does geography or ecology best explain 'cultural' variation among chimpanzee communities?

Much attention has been paid to geographic variation in chimpanzee behavior, but few studies have applied quantitative techniques to explain this variation. Here, we apply methods typically utilized in macroecology to explain variation in the putative cultural traits of chimpanzees. We analyzed published data containing 39 behavioral traits from nine chimpanzee communities. We used a canonical correspondence analysis to examine the relative importance of environmental characteristics and geography, which may be a proxy for inter-community gene flow and/or social transmission, for explaining geographic variation in chimpanzee behavior. We found that geography, and longitude in particular, was the best predictor of behavioral variation. Chimpanzee communities in close longitudinal proximity to each other exhibit similar behavioral repertoires, independent of local ecological factors. No ecological variables were significantly related to behavioral variation. These results support the idea that inter-community dispersal patterns have played a major role in structuring behavioral variation. We cannot be certain whether behavioral variation has a genetic basis, is the result of innovation and diffusion, or a combination of the two.     

History matters more when explaining genetic diversity within the context of the core–periphery hypothesis

The core–periphery hypothesis (CPH) predicts that populations located at the periphery of a species' range should have lower levels of genetic variation than those at the centre of the range. However, most of the research on the CPH focuses on geographic distance and not on ecological distance, or uses categorical definitions of core and periphery to explain the distribution of genetic diversity. We use current climate data and historical climate data from the last glacial maxima to develop quantitative estimates of contemporary and historical ecological suitability using ecological niche models. We analysed genetic diversity using 12 polymorphic microsatellites to estimate changes in heterozygosity, allelic richness and population differentiation in 31 populations of the wood frog (Lithobates sylvaticus) spanning the species’ entire eastern clade (33^o to 45^o latitude) from Alabama, USA, to Nova Scotia, Canada. Our data support predictions based on the CPH. Populations showed significant differences in genetic diversity across the range, with lower levels of genetic variation at the geographic range edge and in areas with lower levels of historical and contemporary ecological suitability. However, history and geography (not current ecological suitability) best explain the patterns. This study highlights the importance of examining more than just geography when assessing the CPH, and the importance of historical ecological suitability in the maintenance of genetic diversity and population differentiation.     

Ecological specialization correlates with genotypic differentiation in sympatric host-populations of the pea aphid

The pea aphid, Acyrthosiphon pisum, encompasses distinct host races specialized on various Fabaceae species, but the extent of genetic divergence associated with ecological specialization varies greatly depending on plant and geographic origins of aphid populations. Here, we studied the genetic structure of French sympatric pea aphid populations collected on perennial (pea and faba bean) and annual (alfalfa and red clover) hosts using 14 microsatellite loci. Classical and Bayesian population genetics analyses consistently identified genetic clusters mostly related to plant origin: the pea/faba bean cluster was highly divergent from the red clover and the alfalfa ones, indicating they represent different stages along the continuum of genetic differentiation. Some genotypes were assigned to a cluster differing from the one expected from their plant origin while others exhibited intermediate genetic characteristics. These results suggest incomplete barriers to gene flow. However, this limited gene flow seems insufficient to prevent ecological specialization and genetic differentiation in sympatry.     

Accounting for regional niche variation in habitat suitability models

Ecological niche modeling has become an increasingly important tool to address issues in many fields of basic and applied ecology. The ecological niche space occupied across the geographic range, particularly for wide-ranging species, may vary for a variety of evolutionary and non-evolutionary reasons. However, ecological niche models are often applied over large geographic areas without regard for the potential effects of regional variation in adaptation, environmental conditions and their interactions, and species responses, thus significantly reducing their accuracy and utility. We develop regionally partitioned ecological niche models, using GARP, for the wide-ranging North American tree Gleditsia triacanthos (Fabaceae) . Models were constructed based on known tree occurrences at peripheral and range-centre locations, as well as across the geographic range as a whole.
Our results suggest that the niche space occupied by G. triacanthos varies regionally and that between some regions in particular there may be a complete absence of niche overlap. In particular, while there is some overlap between the niche space occupied by trees in the western and central regions of the range, there appears to be virtually no overlap in the niche space occupied by trees in the south of the range and that occupied by western and central trees. This lack of overlap appears to be driven primarily by regional differences in abiotic conditions, rather than regional adaptation per se. The results of our study have several important implications for the future development of habitat suitability models over large geographic areas. Spatial partitioning of data is clearly necessary to improve predictions of models where regional niche variation occurs. For wide-ranging species in particular, regional differences in ecological characteristics may cause apparent niche variation.