Post-pollination barriers do not explain the persistence of two distinct Antirrhinum subspecies with parapatric distribution

Empirical studies of post-pollination barriers to gene flow between recently diverged plant species are important to understand ecological processes underlying speciation. Using greenhouse and common garden experiments, we investigated the strength of post-pollination barriers that restrict or prevent gene flow between two subspecies of Antirrhinum: Antirrhinum majus pseudomajus and A. m. striatum. The two are distributed parapatrically but share the same major pollinators (bumblebees), and form narrow hybrid zone in many areas of southern France and northern Spain where they come into close contact. We assessed the strength and symmetry of mating barriers and their homogeneity among populations by comparing fruit set, seed set and offspring performances between intra- and inter-subspecific crosses performed in parental and hybrid populations. Although all populations showed high levels of self-incompatibility, we found very little evidence for barriers to gene flow once pollen had been transferred to the stigma. We suggest that reproductive isolation in this system mostly involves barriers related to processes occurring before pollination, with little or no role of post-pollination barriers.     

Locally asymmetric introgressions between subspecies suggest circular range expansion at the Antirrhinum majus global scale

Assessing processes of geographic expansion in contact zones is a crucial step towards an accurate prediction of the evolution of species genetic diversity. The geographic distribution of cytonuclear discordance often reflects genetic introgression patterns across a species geographic range. Antirrhinum majus pseudomajus and A. m. striatum are two interfertile subspecies that occupy nonoverlapping areas but enter in contact in many locations at the margin of their geographic distribution. We found that genetic introgression between both subspecies was asymmetric at the local scale and geographically oriented in opposite directions at both ends of their contact zone perimeter in the Pyrenees. Our results suggest that the geographic expansion of A. majus subspecies was circular around the perimeter of their contact zone and pinpoint the need to integrate different spatial scales to unravel complex patterns of species geographic expansion.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.     

The Influence of Prior Learning Experience on Pollinator Choice: An Experiment Using Bumblebees on Two Wild Floral Types of Antirrhinum majus

Understanding how pollinator behavior may influence pollen transmission across floral types is a major challenge, as pollinator decision depends on a complex range of environmental cues and prior experience. Here we report an experiment using the plant Antirrhinum majus and the bumblebee Bombus terrestris to investigate how prior learning experience may affect pollinator preferences between floral types when these are presented together. We trained naive bumblebees to forage freely on flowering individuals of either A. majus pseudomajus (magenta flowers) or A. majus striatum (yellow flowers) in a flight cage. We then used a Y-maze device to expose trained bumblebees to a dual choice between the floral types. We tested the influence of training on their choice, depending on the type of plant signals available (visual signals, olfactory signals, or both). Bumblebees had no innate preference for either subspecies. Bumblebees trained on the yellow-flowered subspecies later preferred the yellow type, even when only visual or only olfactory signals were available, and their preference was not reinforced when both signal types were available. In contrast, bumblebees trained on the magenta-flowered subspecies showed no further preference between floral types and took slightly more time to make their choice. Since pollinator constancy has been observed in wild populations of A. majus with mixed floral types, we suggest that such constancy likely relies on short-term memory rather than acquired preference through long-term memory induced by prior learning.     

The role of climate for the range limits of parapatric European land salamanders

Abrupt range limits of parapatric species may serve as a model system to understand the factors that determine species’ range borders. Theory suggests that parapatric range limits can be caused by abiotic conditions along environmental gradients, biotic interactions or a combination of both. Geographic ranges of the parapatric salamanders, Salamandra salamandra and S. atra, meet in small contact zones in the European Alps and to date, the cause of parapatry and the restricted range of S. atra remain elusive. We combine multivariate approaches and climatic data analysis to explore niche differentiation among the two salamanders with respect to the available climatic environment at their contact zones. Our purpose is to evaluate whether climatic conditions explain the species’ sharp range limits or if biotic interactions may play a role for range delimitation. Analyses were carried out in three contact zones in Switzerland to assess possible geographic variation. Our results indicate that both species occur at localities with different climatic conditions as well as the presence of a strong climatic gradient across the species’ range limits. Although the species’ climatic niches differ moderately (with a wider niche breadth for S. atra), interspecific niche overlap is found. Comparisons among the contact zones confirm geographic variation in the species’ climatic niches as well as in the conditions within the geographically available space. Our results suggest that the change in climatic conditions along the recognized gradient represents a determining factor for species’ range limits within contact zones. However, our analyses of geographic variation in climatic conditions reveal that both salamander species can occur in a much wider range of conditions than observed within contact zones. This finding and the interspecific climatic niche overlap within each contact zone provides indirect evidence that biotic interactions (likely competition) between the two species may also determine their range limits.     

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long‐toed salamander (Ambystoma macrodactylum Baird)

The relative importance of ecological vs. nonecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages – such as the distinct genetic groups that make up the ranges of many northern species – represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long‐toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long‐toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non‐ecological factors to play a more important role in the maintenance of northern taxa.     

The use of geostatistics and GIS for evolutionary history studies: the case of the nose‐horned viper (Vipera ammodytes) in the Balkan Peninsula

Geostatistics and geographical information system (GIS) procedures are novel techniques helpful for the identification of environmental correlates sustaining contact zones among subspecies or closely related species. In this paper, we tried to infer evolutionary scenarios for Vipera ammodytes across the European part of its distribution area using geostatistics and ecological niche‐based models, hence trying to solve several biogeographical questions that remained unclear after the application of classical morphological tools and genetic analyses. Eleven morphological traits from 871 vipers were analysed with geostatistics and ecological niche‐based modelling. Interpolation by kriging was used to generate surfaces of morphological variation, which were combined with spatial principal components analysis (SPCA). SPCA maps were used to test putative morphological differentiated groups with discriminant function analysis (DFA). Maximum entropy modelling and seven environmental variables were used to identify factors limiting the distribution of groups and areas for the potential occurrence of such groups. Three patterns of morphological variation were observed: a north‐west/south‐west cline, transition zones with steep clines of variation in a west–east arc, and particular character traits that disturbed the general cline. SPCA identified between three and nine putative population groups, of which three were supported by DFA. Areas of potential occurrence of these groups were coherent with the range of the three subspecies of V. ammodytes currently recognized. The distribution of all subspecies was mostly related to precipitation in the driest month. Areas of probable sympatry between subspecies are generally small and restricted. The main patterns of geographic variation of morphological characters for V. ammodytes were similar to the patterns obtained for Vipera latastei and Vipera monticola; the same environmental factors limit the distribution of differentiated groups of vipers in the Balkans and the Iberian Peninsula. The influence of humidity on the variation of morphological traits in spatially separated viper taxa from the two European peninsulas coincides with their phylogenetic relatedness. Geostatistics and GIS procedures were successful in the identification of environmental correlates sustaining contact zones among V. ammodytes subspecies in the Balkans. The same techniques should be applied for studying other parapatric forms and refugia regions. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 651–666.     

Ecological Niche Transferability Using Invasive Species as a Case Study

Species distribution modeling is widely applied to predict invasive species distributions and species range shifts under climate change. Accurate predictions depend upon meeting the assumption that ecological niches are conserved, i.e., spatially or temporally transferable. Here we present a multi-taxon comparative analysis of niche conservatism using biological invasion events well documented in natural history museum collections. Our goal is to assess spatial transferability of the climatic niche of a range of noxious terrestrial invasive species using two complementary approaches. First we compare species’ native versus invasive ranges in environmental space using two distinct methods, Principal Components Analysis and Mahalanobis distance. Second we compare species’ native versus invaded ranges in geographic space as estimated using the species distribution modeling technique Maxent and the comparative index Hellinger’s I. We find that species exhibit a range of responses, from almost complete transferability, in which the invaded niches completely overlap with the native niches, to a complete dissociation between native and invaded ranges. Intermediate responses included expansion of dimension attributable to either temperature or precipitation derived variables, as well as niche expansion in multiple dimensions. We conclude that the ecological niche in the native range is generally a poor predictor of invaded range and, by analogy, the ecological niche may be a poor predictor of range shifts under climate change. We suggest that assessing dimensions of niche transferability prior to standard species distribution modeling may improve the understanding of species’ dynamics in the invaded range.     

The Role of Geographical and Ecological Factors on Population Divergence of the Neotropical otter <Emphasis Type="Italic">Lontra longicaudis</Emphasis> (Carnivora,Mustelidae)

The geographic distribution of the populations of a species are influenced by the spatial structure of the ecosystems, the environmental factors and the presence of geographic barriers. The Neotropical otter, Lontra longicaudis, is widely distributed throughout the Americas, where a wide range of environmental conditions and geographical features could promote genetic and morphological variation on the three currently recognized subspecies. In this study, we combined phylogeographic, morphometric and environmental niche modelling analyses to examine whether: (1) genetic variation is associated with the presence of barriers to gene flow and/or hydrography; (2) genetic and morphologic variation are associated with environmental variation; and (3) the observed variation in L. longicaudis populations corresponds to the previously defined subspecies. We found strong phylogeographic structure between the northern (L. l. annectens) and the two-southern subspecies (L. l. longicaudis and L. l. enudris), and although shallower, we also detected genetic differentiation between the two South American subspecies. Such genetic differentiation corresponds to the hydrography and to the geographical barriers characteristic of the distributional area of the species. We found a correlation between the shape of the skull and mandible with the environmental variation through the distribution of the species, and we rejected the hypothesis of niche equivalency and similarity between the three identified genetic lineages, suggesting adaptations to different environmental conditions. Our results support that the variation in environmental conditions, in concert with geographical barriers to gene flow and hydrography, have led to population divergence of L. longicaudis along the Neotropics. These results have important taxonomic implications for the species and its conservation.     

Differences between ecological niches in northern and southern populations of Angolan black and white colobus monkeys (Colobus angolensis palliatus and Colobus angolensis sharpei) throughout Kenya and Tanzania

Ecological niche models can be useful for clarifying relationships between environmental factors and a species’ geographic distribution. In this study, we use presence‐only data and environmental layers to create an ecological niche model to better understand the distribution of the East African Angolan black and white colobus monkey, Colobus angolensis palliatus, and to assess whether the model supports considering the population as two separate subspecies, Colobus angolensis sharpei and C. a. palliatus. We found the range of the predicted distribution for suitable habitat of C. a. palliatus as currently classified to be only 12.4% of that shown in the International Union for Conservation of Nature Red List range map and to be fragmented. As C. angolensis is considered a “Least Concern” species, this difference suggests that generalized maps may lead to understating the species’ extinction risk. When presence points were divided into two previously proposed subspecies —C. a. palliatus (Kenya and Northern Tanzania) and C. a. sharpei (Southern Tanzania)—we found significant environmental differences between the distributions. The most important ecological variable for C. a. palliatus was predominantly precipitation of the driest month (69.1%) whereas for C. a. sharpei annual precipitation (44.8%) and land cover (normalized difference vegetation index, 16.4%) were the most important. When comparing suitable ranges for the separate distributions, we found only a 1.2% geographical overlap. These differences are consistent with previous subspecies delineations of C. a. palliatus and C. a. sharpei based upon morphology, pelage, and genetics. Our study suggests that extirpation of C. a. palliatus in suitable habitat areas and occurrence of this subspecies in anthropogenic environments, warrant further consideration for conservation actions.     

Quantifying the climatic niche of symbiont partners in a lichen symbiosis indicates mutualist‐mediated niche expansions

The large distributional areas and ecological niches of many lichenized fungi may in part be due to the plasticity in interactions between the fungus (mycobiont) and its algal or cyanobacterial partners (photobionts). On the one hand, broad‐scale phylogenetic analyses show that partner compatibility in lichens is rather constrained and shaped by reciprocal selection pressures and codiversification independent of ecological drivers. On the other hand, sub‐species‐level associations among lichen symbionts appear to be environmentally structured rather than phylogenetically constrained. In particular, switching between photobiont ecotypes with distinct environmental preferences has been hypothesized as an adaptive strategy for lichen‐forming fungi to broaden their ecological niche. The extent and direction of photobiont‐mediated range expansions in lichens, however, have not been examined comprehensively at a broad geographic scale. Here we investigate the population genetic structure of Lasallia pustulata symbionts at sub‐species‐level resolution across the mycobiont's Europe‐wide range, using fungal MCM7 and algal ITS rDNA sequence markers. We show that variance in occurrence probabilities in the geographic distribution of genetic diversity in mycobiont‐photobiont interactions is closely related to changes in climatic niches. Quantification of niche extent and overlap based on species distribution modeling and construction of Hutchinsonian climatic hypervolumes revealed that combinations of fungal–algal interactions change at the sub‐species level along latitudinal temperature gradients and in Mediterranean climate zones. Our study provides evidence for symbiont‐mediated niche expansion in lichens. We discuss our results in the light of symbiont polymorphism and partner switching as potential mechanisms of environmental adaptation and niche evolution in mutualisms.     

Niche Divergence in a Brown Lemur (<Emphasis Type="Italic">Eulemur</Emphasis> spp.) Hybrid Zone: Using Ecological Niche Models to Test Models of Stability

Endogenous selection is often implicated in the maintenance of stability of natural hybrid zones. Environmental conditions often vary across these zones, suggesting that local adaptation to ecological conditions could also play a role in this process. We used niche modeling to investigate these alternatives in a hybrid zone between two species of brown lemur (Eulemur rufifrons and E. cinereiceps) in southeastern Madagascar. We produced ecological niche models (ENMs) for parental and hybrid populations and compared values of niche overlap to null expectations using identity and background tests. All three taxonomic groups had nonequivalent ENMs with limited spatial overlap, supporting a role for niche divergence and local adaptation in the maintenance of this zone. However, values of niche overlap between ENMs were not greater than null expectations controlling for background environmental differences. These results could suggest that taxa in this hybrid zone inhabit portions of their environments that are more similar to their backgrounds, i.e., niche conservatism. Nevertheless, we did find evidence of niche divergence when using background tests that examined environmental variables separately. Although we could not rule out models indicating selection against hybrids, most lines of evidence were consistent with predictions for the bounded superiority model of hybrid zone stability. This study thus provides support that exogenous, environmental selection may be responsible for maintaining the hybrid zone, and may be implicated in the evolutionary divergence of these taxa.     

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.     

An everlasting pioneer: the story of Antirrhinum research

Despite the tremendous success of Arabidopsis thaliana, no single model can represent the vast range of form that is seen in the approximately 250,000 existing species of flowering plants (angiosperms). Here, we consider the history and future of an alternative angiosperm model--the snapdragon Antirrhinum majus. We ask what made Antirrhinum attractive to the earliest students of variation and inheritance, and how its use led to landmark advances in plant genetics and to our present understanding of plant development. Finally, we show how the wide diversity of Antirrhinum species, combined with classical and molecular genetics--the two traditional strengths of Antirrhinum--provide an opportunity for developmental, evolutionary and ecological approaches. These factors make A. majus an ideal comparative angiosperm.     

Deep Genetic Divergence between Disjunct Refugia in the Arctic-Alpine King’s Crown,Rhodiola integrifolia (Crassulaceae)

Despite the strength of climatic variability at high latitudes and upper elevations, we still do not fully understand how plants in North America that are distributed between Arctic and alpine areas responded to the environmental changes of the Quaternary. To address this question, we set out to resolve the evolutionary history of the King’s Crown, Rhodiola integrifolia using multi-locus population genetic and phylogenetic analyses in combination with ecological niche modeling. Our population genetic analyses of multiple anonymous nuclear loci revealed two major clades within R. integrifolia that diverged from each other ~ 700 kya: one occurring in Beringia to the north (including members of subspecies leedyi and part of subspecies integrifolia), and the other restricted to the Southern Rocky Mountain refugium in the south (including individuals of subspecies neomexicana and part of subspecies integrifolia). Ecological niche models corroborate our hypothesized locations of refugial areas inferred from our phylogeographic analyses and revealed some environmental differences between the regions inhabited by its two subclades. Our study underscores the role of geographic isolation in promoting genetic divergence and the evolution of endemic subspecies in R. integrifolia. Furthermore, our phylogenetic analyses of the plastid spacer region trnL-F demonstrate that among the native North American species, R. integrifolia and R. rhodantha are more closely related to one another than either is to R. rosea. An understanding of these historic processes lies at the heart of making informed management decisions regarding this and other Arctic-alpine species of concern in this increasingly threatened biome.     

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.     

Application of epidermal characters to the taxonomy of European species of Antirrhinum (Schrophulariaceae)

Epidermal characters provide additional features in the assessment of taxonomic relationships among the European species of Antirrhinum. Much of this study confirms the classification derived by Webb but important differences have been observed, viz. A. mollissimum may be considered as identical to A. molle whilst A. boissieri, A. ambiguum and A. rupestre may be assigned as distinct species. Evidence from epidermal characters supports Rothmaler's view that A. majus spp. litigiosum is a subspecies of A. majus and should not be amalgamated with A. barrelieri as proposed by Webb. A. × huteri is proposed to be a hybrid between A. hispanicum and A. boissieri. Thus it is proposed that the section Antirrhinum comprises 20 species with A. majus divided into five subspecies.     

Lack of congruence of genetic and niche divergence in Podarcis hispanicus complex

Niche divergence among closely related lineages can be informative on the ecological and evolutionary processes involved in differentiation, particularly in the case of cryptic species complexes. Here we compared phylogenetic relationships and niche similarity between pairs of lineages included in the Podarcis hispanicus complex to examine patterns of niche divergence and its role in the organization of current diversity patterns, as allopatric, parapatric, and sympatric lineages occur in the Western Mediterranean Basin. First, we used ecological niche models to characterize the realized climatic niche of each Podarcis hispanicus complex lineage based on topographic and climatic variables, to identify important variables, and to test for niche conservatism or divergence between pairs of lineages. Variables related to precipitation generally exhibited the highest contribution to niche models, highlighting the importance of rainfall levels in shaping distributions of Podarcis wall lizards. We found that most forms have significant differences in realized climatic niches that do not follow the pattern of mitochondrial divergence. These results lend support to the hypothesis that genetic divergence across Podarcis hispanicus complex most likely occurred in allopatric conditions, mostly with significant niche divergence. Competition after secondary contact is also suggested by the common occurrence of niche overlap between lineages that exhibit strictly parapatric distribution. The almost continuous distribution of Podarcis lizards in the study area appears to be a result of a combination of complementary suitable niches and competition, which seem two important mechanisms limiting geographic distributions and restricting the existence of extensive contact zones.