GENETIC ARCHITECTURE, EVOLUTION AND MAINTENANCE OF A COLOR CLINE AT A CONTINENTAL SCALE

On a geological time scale the conditions on earth are very variable and biological patterns (for example the distributions of species) are very dynamic. Understanding large scale patterns of variation observed today thus requires a deep understanding of the historical factors that drove their evolution. In this thesis, we reevaluated the evolution and maintenance of a continental color cline observed in the European barn owl (Tyto alba) using population genetic tools. The colour cline spans from south-est Europe where most individual have pure white underparts to north and east Europe where most individuals have rufous-brown underparts. Our results globally showed that the old scenario, stipulating that the color cline evolved by secondary contact of two color morphs (white and rufous) that evolved in allopatry during the last ice age has to be revised. We collected samples of about 700 barn owls from the Western Palearctic to establish the first population genetic data set for this species. Individuals were genotyped at 22 microsatellites markers, at one mitochondrial gene, and at a candidate color gene. The color of each individuals was assessed and their sex determined by molecular methods. We first showed that the genetic variation in Western Europe is very limited compared to the heritable color variation. We found no evidences of different glacial lineages, and showed that selection must be involved in the maintenance of the color cline (chapter 1). Using computer simulations, we demonstrated that the post-glacial colonization of Europe occurred from the Iberian Peninsula and that the color cline could not have evolved by neutral demographic processes during this colonization (chapter 2). Finally we reevaluated the whole history of the establishment of the Western Palearctic variation of the barn owl (chapter 3): This study showed that all Western European barn owls descend from white barn owls phenotypes from the Middle East that colonized the Iberian Peninsula via North-Africa. Following the end of the last ice age (20'000 years ago), these white barn owls colonized Western Europe and under selection a novel rufous phenotype evolved (during or after the colonization). An important part of the color variation could be explained by a single mutation in the melanocortin-1-receptor (MC1R) gene that appeared during or after the colonization. The colonization of Europe reached until Greece, where the rufous birds encountered white ones (which reached Greece from the Middle East over the Bosporus) in a secondary contact zone. Our analyses show that white and rufous barn owls in Greece interbreed only to a limited extent. This suggests that barn owls are at the verge of becoming two species in Greece and demonstrates that European barn owls represent an incipient ring species around the Mediterranean. The revisited history of the establishment of the European barn owl color cline makes this model system remarkable for several aspects. It is a very clear example of strong local adaptation that can be achieved despite high gene flow (strong color and MC1R differentiation despite almost no neutral genetic differentiation). It also offers a wonderful model system to study the interactions between colonization processes and selection processes which have, for now, been remarkably understudied despite their potentially ubiquitous importance. Finally it represents a very interesting case in the speciation continuum and appeals for further studying the amount of gene flow that occurs between the color morphs in Greece. -- Sur l'échelle des temps géologiques, les conditions sur terre sont très variables et les patrons biologiques (telle que la distribution des espèces) sont très dynamiques. Si l'on veut comprendre des patrons que l'on peut observer à large échelle aujourd'hui, il est nécessaire de d'abord comprendre les facteurs historiques qui ont gouverné leur établissement. Dans cette thèse, nous allons réévaluer, grâce à des outils modernes de génétique des populations, l'évolution et la maintenance d'un cline de couleur continental observé chez l'effraie des clochers européenne (Tyto alba). Globalement, nos résultats montrent que le scenario accepté jusqu'à maintenant, qui stipule que le cline de couleur a évolué à partir du contact secondaire de deux morphes de couleur (blanches et rousses) ayant évolué en allopatrie durant les dernières glaciations, est à revoir. Afin de constituer le premier jeu de données de génétique des populations pour cette espèce, nous avons récolté des échantillons d'environ 700 effraies de l'ouest Paléarctique. Nous avons génotypé tous les individus à 22 loci microsatellites, sur un gène mitochondrial et sur un autre gène participant au déterminisme de la couleur. Nous avons aussi mesuré la couleur de tous les individus et déterminé leur sexe génétiquement. Nous avons tout d'abord pu montrer que la variation génétique neutre est négligeable en comparaison avec la variation héritable de couleur, qu'il n'existe qu'une seule lignée européenne et que de la sélection doit être impliquée dans le maintien du cline de couleur (chapitre 1). Grâce à des simulations informatiques, nous avons démontré que l'ensemble de l'Europe de l'ouest a été recolonisé depuis la Péninsule Ibérique après les dernières glaciations et que le cline de couleur ne peut pas avoir évolué par des processus neutre durant cette colonisation (chapitre 2). Finalement, nous avons réévalué l'ensemble de l'histoire postglaciaire de l'espèce dans l'ouest Paléarctique (chapitre 3): l'ensemble des effraies du Paléarctique descendent d'effraie claire du Moyen-Orient qui ont colonisé la péninsule ibérique en passant par l'Afrique du nord. Après la fin de la dernière glaciation (il y a 20'000 ans), ces effraies claires ont colonisé l'Europe de l'ouest et ont évolués par sélection le phénotype roux (durant ou après la colonisation). Une part importante de la variation de couleur peut être expliquée par une mutation sur le gène MC1R qui est apparue durant ou juste après la colonisation. Cette vague de colonisation s'est poursuivie jusqu'en Grèce où ces effraies rousses ont rencontré dans une zone de contact secondaire des effraies claires (qui sont remontées en Grèce depuis le Moyen-Orient via le Bosphore). Nos analyses montrent que le flux de gènes entre effraies blanches et rousses est limité en Grèce, ce qui suggère qu'elles sont en passe de former deux espèces et ce qui montre que les effraies constituent un exemple naissant de spéciation en anneaux autour de la Méditerranée. L'histoire revisitée des effraies des clochers de l'ouest Paléarctique en fait un système modèle remarquable pour plusieurs aspects. C'est un exemple très claire de forte adaptation locale maintenue malgré un fort flux de gènes (différenciation forte de couleur et sur le gène MC1R malgré presque aucune structure neutre). Il offre également un très bon système pour étudier l'interaction entre colonisation et sélection, un thème ayant été remarquablement peu étudié malgré son importance. Et il offre finalement un cas très intéressant dans le « continuum de spéciation » et il serait très intéressant d'étudier plus en détail l'importance du flux de gènes entre les morphes de couleur en Grèce.

Multi-scale support vector algorithms for hot spot detection and modelling

The algorithmic approach to data modelling has developed rapidly these last years, in particular methods based on data mining and machine learning have been used in a growing number of applications. These methods follow a data-driven methodology, aiming at providing the best possible generalization and predictive abilities instead of concentrating on the properties of the data model. One of the most successful groups of such methods is known as Support Vector algorithms. Following the fruitful developments in applying Support Vector algorithms to spatial data, this paper introduces a new extension of the traditional support vector regression (SVR) algorithm. This extension allows for the simultaneous modelling of environmental data at several spatial scales. The joint influence of environmental processes presenting different patterns at different scales is here learned automatically from data, providing the optimum mixture of short and large-scale models. The method is adaptive to the spatial scale of the data. With this advantage, it can provide efficient means to model local anomalies that may typically arise in situations at an early phase of an environmental emergency. However, the proposed approach still requires some prior knowledge on the possible existence of such short-scale patterns. This is a possible limitation of the method for its implementation in early warning systems. The purpose of this paper is to present the multi-scale SVR model and to illustrate its use with an application to the mapping of Cs137 activity given the measurements taken in the region of Briansk following the Chernobyl accident.

Providing more informative projections of climate change impact on plant distribution in a mountain environment

Summary Due to their conic shape and the reduction of area with increasing elevation, mountain ecosystems were early identified as potentially very sensitive to global warming. Moreover, mountain systems may experience unprecedented rates of warming during the next century, two or three times higher than that records of the 20th century. In this context, species distribution models (SDM) have become important tools for rapid assessment of the impact of accelerated land use and climate change on the distribution plant species. In my study, I developed and tested new predictor variables for species distribution models (SDM), specific to current and future geographic projections of plant species in a mountain system, using the Western Swiss Alps as model region. Since meso- and micro-topography are relevant to explain geographic patterns of plant species in mountain environments, I assessed the effect of scale on predictor variables and geographic projections of SDM. I also developed a methodological framework of space-for-time evaluation to test the robustness of SDM when projected in a future changing climate. Finally, I used a cellular automaton to run dynamic simulations of plant migration under climate change in a mountain landscape, including realistic distance of seed dispersal. Results of future projections for the 21st century were also discussed in perspective of vegetation changes monitored during the 20th century. Overall, I showed in this study that, based on the most severe A1 climate change scenario and realistic dispersal simulations of plant dispersal, species extinctions in the Western Swiss Alps could affect nearly one third (28.5%) of the 284 species modeled by 2100. With the less severe 61 scenario, only 4.6% of species are predicted to become extinct. However, even with B1, 54% (153 species) may still loose more than 80% of their initial surface. Results of monitoring of past vegetation changes suggested that plant species can react quickly to the warmer conditions as far as competition is low However, in subalpine grasslands, competition of already present species is probably important and limit establishment of newly arrived species. Results from future simulations also showed that heavy extinctions of alpine plants may start already in 2040, but the latest in 2080. My study also highlighted the importance of fine scale and regional. assessments of climate change impact on mountain vegetation, using more direct predictor variables. Indeed, predictions at the continental scale may fail to predict local refugees or local extinctions, as well as loss of connectivity between local populations. On the other hand, migrations of low-elevation species to higher altitude may be difficult to predict at the local scale. Résumé La forme conique des montagnes ainsi que la diminution de surface dans les hautes altitudes sont reconnues pour exposer plus sensiblement les écosystèmes de montagne au réchauffement global. En outre, les systèmes de montagne seront sans doute soumis durant le 21ème siècle à un réchauffement deux à trois fois plus rapide que celui mesuré durant le 20ème siècle. Dans ce contexte, les modèles prédictifs de distribution géographique de la végétation se sont imposés comme des outils puissants pour de rapides évaluations de l'impact des changements climatiques et de la transformation du paysage par l'homme sur la végétation. Dans mon étude, j'ai développé de nouvelles variables prédictives pour les modèles de distribution, spécifiques à la projection géographique présente et future des plantes dans un système de montagne, en utilisant les Préalpes vaudoises comme zone d'échantillonnage. La méso- et la microtopographie étant particulièrement adaptées pour expliquer les patrons de distribution géographique des plantes dans un environnement montagneux, j'ai testé les effets d'échelle sur les variables prédictives et sur les projections des modèles de distribution. J'ai aussi développé un cadre méthodologique pour tester la robustesse potentielle des modèles lors de projections pour le futur. Finalement, j'ai utilisé un automate cellulaire pour simuler de manière dynamique la migration future des plantes dans le paysage et dans quatre scénarios de changement climatique pour le 21ème siècle. J'ai intégré dans ces simulations des mécanismes et des distances plus réalistes de dispersion de graines. J'ai pu montrer, avec les simulations les plus réalistes, que près du tiers des 284 espèces considérées (28.5%) pourraient être menacées d'extinction en 2100 dans le cas du plus sévère scénario de changement climatique A1. Pour le moins sévère des scénarios B1, seulement 4.6% des espèces sont menacées d'extinctions, mais 54% (153 espèces) risquent de perdre plus 80% de leur habitat initial. Les résultats de monitoring des changements de végétation dans le passé montrent que les plantes peuvent réagir rapidement au réchauffement climatique si la compétition est faible. Dans les prairies subalpines, les espèces déjà présentes limitent certainement l'arrivée de nouvelles espèces par effet de compétition. Les résultats de simulation pour le futur prédisent le début d'extinctions massives dans les Préalpes à partir de 2040, au plus tard en 2080. Mon travail démontre aussi l'importance d'études régionales à échelle fine pour évaluer l'impact des changements climatiques sur la végétation, en intégrant des variables plus directes. En effet, les études à échelle continentale ne tiennent pas compte des micro-refuges, des extinctions locales ni des pertes de connectivité entre populations locales. Malgré cela, la migration des plantes de basses altitudes reste difficile à prédire à l'échelle locale sans modélisation plus globale.

Effect of biogeographic history on population vulnerability in European amphibians.

The genetic diversity of populations, which contributes greatly to their adaptive potential, is negatively affected by anthropogenic habitat fragmentation and destruction. However, continental-scale losses of genetic diversity also resulted from the population expansions that followed the end of the last glaciation, an element that is rarely considered in a conservation context. We addressed this issue in a meta-analysis in which we compared the spatial patterns of vulnerability of 18 widespread European amphibians in light of phylogeographic histories (glacial refugia and postglacial routes) and anthropogenic disturbances. Conservation statuses significantly worsened with distances from refugia, particularly in the context of industrial agriculture; human population density also had a negative effect. These findings suggest that features associated with the loss of genetic diversity in post-glacial amphibian populations (such as enhanced fixation load or depressed adaptive potential) may increase their susceptibility to current threats (e.g., habitat fragmentation and pesticide use). We propose that the phylogeographic status of populations (i.e., refugial vs. post-glacial) should be considered in conservation assessments for regional and national red lists.

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus, a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus, yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.

21st century climate change threatens mountain flora unequally across Europe

Continental-scale assessments of 21st century global impacts of climate change on biodiversity have forecasted range contractions for many species. These coarse resolution studies are, however, of limited relevance for projecting risks to biodiversity in mountain systems, where pronounced microclimatic variation could allow species to persist locally, and are ill-suited for assessment of species-specific threat in particular regions. Here, we assess the impacts of climate change on 2632 plant species across all major European mountain ranges, using high-resolution (ca. 100 m) species samples and data expressing four future climate scenarios. Projected habitat loss is greater for species distributed at higher elevations; depending on the climate scenario, we find 36-55% of alpine species, 31-51% of subalpine species and 19-46% of montane species lose more than 80% of their suitable habitat by 2070-2100. While our high-resolution analyses consistently indicate marked levels of threat to cold-adapted mountain florae across Europe, they also reveal unequal distribution of this threat across the various mountain ranges. Impacts on florae from regions projected to undergo increased warming accompanied by decreased precipitation, such as the Pyrenees and the Eastern Austrian Alps, will likely be greater than on florae in regions where the increase in temperature is less pronounced and rainfall increases concomitantly, such as in the Norwegian Scandes and the Scottish Highlands. This suggests that change in precipitation, not only warming, plays an important role in determining the potential impacts of climate change on vegetation.

A concentric analysis of the impact of urbanization on the threatened European tree frog in an agricultural landscape

Pond-breeding amphibians are affected by site-specific factors and regional and landscape-scale patterns of land use. Recent anthropogenic landscape modifications (drainage, agriculture intensification, larger road networks, and increased traffic) affect species by reducing the suitable habitat area and fragmenting remaining populations. Using a robust concentric approach based on permutation tests, we evaluated the impact of recent landscape changes on the presence of the endangered European tree frog (Hyla arborea.) in wetlands. We analyzed the frequency of 1 traffic and 14 land-use indices at 20 circular ranges (from 100-m up to 2-km radii) around 76 ponds identified in western Switzerland. Urban areas and road surfaces had a strong adverse effect on tree frog presence even at relatively great distances (from 100 m up to 1 km). When traffic measurements were considered instead of road surfaces, the effect increased, suggesting a negative impact due to a vehicle-induced effect. Altogether, our results indicate that urbanization and traffic must be taken into account when pond creation is an option in conservation management plans, as is the case for the European tree frog in western Switzerland. We conclude that our easy-to-use and robust concentric method of analysis can successfully assist managers in identifying potential sites for pond creation, where probability of the presence of tree frogs is maximized.

'Next generation' biogeography: towards understanding the drivers of species diversification and persistence

The drivers of species diversification and persistence are of great interest to current biogeography, especially in those global biodiversity hotspots' harbouring most of Earth's animal and plant life. Classical multispecies biogeographical work has yielded fascinating insights into broad-scale patterns of diversification, and DNA-based intraspecific phylogeographical studies have started to complement this picture at much finer temporal and spatial scales. The advent of novel next-generation sequencing (NGS) technologies provides the opportunity to greatly scale up the numbers of individuals, populations and species sampled, potentially merging intraspecific and interspecific approaches to biogeographical inference. Here, we outline these prospects and issues by using the example of an undisputed hotspot, the Cape of southern Africa. We outline the current state of knowledge on the biogeography of species diversification within the Cape, review the literature for phylogeographical evidence of its likely drivers and mechanisms, and suggest possible ways forward based on NGS approaches. We demonstrate the potential of these methods and current bioinformatic issues with the help of restriction-site-associated DNA (RAD) sequencing data for three highly divergent species of the Restionaceae, an important plant radiation in the Cape. A thorough understanding of the mechanisms that facilitate species diversification and persistence in spatially structured, species-rich environments will require the adoption of novel genomic and bioinformatic tools in biogeographical studies.

Late Pan-African magmatism in the Himalaya: new geochronological and geochemical data from the Ordovician Tso Morari metagranites (Ladakh, NW India)

Two granitic plutons, the Tso Morari gneiss and the Rupshu metagranite, crop out in the Tso Morari area. The Polokongka La granite, classically interpreted as a young intrusion in the Tso Morari gneiss, has been recognized as the undeformed facies of the latter. Conventional isotope dilution U-Pb zircon dating on single-grain and small multi-grain fractions yielded magmatic ages of 479 +/- 2 Ma for the Tso Morari gneiss and the Polokongka La granite, and 482.5 +/- 1 Ma for the Rupshu granite. There is a great difference in zircon morphology between the Tso Morari gneiss (peraluminous type) and the Rupshu granite (alkaline type). This difference is confirmed by whole-rock chemistry. The Tso Morari gneiss is a typical deformed S-type granite, resulting from crustal anatexis. On the other hand, the Rupshu granite is an essentially metaluminous alkali-calcic intrusion derived from a different source material. Data compilation from other Himalayan Cambro-Ordovician granites reveals huge and widespread magmatic activity all along and beyond the northern Indian plate between 570 and 450 Ma, with a peak at 500-480 Ma. A major, continental-scale tectonic event is required to generate such a large magmatic belt; it has been tentatively compared to the Variscan post-orogenic extensional regime of Western Europe, as a late evolution stage of a Pan-African orogenic event.

Shield morphology of curly overhair in 22 genera of Soricidae (Insectivora, Mammalia)

The structure of the shield of curly overhairs was studied in 22 genera of shrews, using scanning electron microscopy. The cross section of the shield is quadriconcave with two sides showing particular scale patterns or a relief which can be grouped into 4 morphological types: 1) a smooth type with, at most, shallow U-shaped notches; 2) a type with uniserial, V-shaped tiled notches; 3) a type with a groove and irregular notches; 4) a type with a deep ridged groove. Myosorex, which occupies a basal phyletic position, shows type 3. This type is therefore interpreted as an ancestral character state. It is found also in Feroculus and in some Sylvisorex. In the Crocidurinae, this type evolved into type 2 (Scutisorex and some Indomalayan Crocidura), and finally into type 1 (Suncus and most Crocidura). In the Soricinae, it evolved into type 4, which is common to all genera except Megasorex and Notiosorex. These two genera possess type 1 which is interpreted as being a synapomorphous character of these genera. The regression of the complex structure under dry climatic conditions supports the hypothesis that the function of the grooved form of hair is water repulsion.

Continent wide response of mountain vegetation to climate change

Climate impact studies have indicated ecological fingerprints of recent global warming across a wide range of habitats. Whereas these studies have shown responses from various local case studies, a coherent large-scale account on temperature-driven changes of biotic communities has been lacking. Here we use 867 vegetation samples above the treeline from 60 summit sites in all major European mountain systems to show that ongoing climate change gradually transforms mountain plant communities. We provide evidence that the more cold-adapted species decline and the more warm-adapted species increase, a process described here as thermophilisation. At the scale of individual mountains this general trend may not be apparent, but at the¦larger, continental scale we observed a significantly higher abundance of thermophilic species in 2008, compared with 2001. Thermophilisation of mountain plant communities mirrors the degree of recent warming and is more pronounced in areas where the temperature increase has been higher. In view of the projected climate change the observed transformation suggests a progressive decline of cold mountain habitats and their biota.

Genetic consequences of population expansions and contractions in the common hippopotamus (Hippopotamus amphibius) since the Late Pleistocene.

Over the past two decades, an increasing amount of phylogeographic work has substantially improved our understanding of African biogeography, in particular the role played by Pleistocene pluvial-drought cycles on terrestrial vertebrates. However, still little is known on the evolutionary history of semi-aquatic animals, which faced tremendous challenges imposed by unpredictable availability of water resources. In this study, we investigate the Late Pleistocene history of the common hippopotamus (Hippopotamus amphibius), using mitochondrial and nuclear DNA sequence variation and range-wide sampling. We documented a global demographic and spatial expansion approximately 0.1-0.3 Myr ago, most likely associated with an episode of massive drainage overflow. These events presumably enabled a historical continent-wide gene flow among hippopotamus populations, and hence, no clear continental-scale genetic structuring remains. Nevertheless, present-day hippopotamus populations are genetically disconnected, probably as a result of the mid-Holocene aridification and contemporary anthropogenic pressures. This unique pattern contrasts with the biogeographic paradigms established for savannah-adapted ungulate mammals and should be further investigated in other water-associated taxa. Our study has important consequences for the conservation of the hippo, an emblematic but threatened species that requires specific protection to curtail its long-term decline.

Characteristic reflection patterns in the southeast Canadian cordillera, northern Appalachians and Swiss Alps

There are some striking similarities and some differences between the seismic reflection sections recorded across the fold and thrust belts of the southeast Canadian Cordillera, Quebec-Maine Appalachians and Swiss Alps. In the fold and thrust belts of all three mountain ranges, seismic reflection surveys have yielded high-quality images of. (1) nappes (thin thrust sheets) stacked on top of ancient continental margins; (2) ramp anticlines in the hanging walls of faults that have ramp-flat or listric geometries; (3) back thrusts and back folds that developed during the terminal phases of orogeny; and (4) tectonic wedges and regional decollements. A principal result of the Cordilleran and Appalachian deep crustal studies has been the recognition of master decollements along which continental margin strata have been transported long distances, whereas a principal result of the Swiss Alpine deep crustal program has been the identification of the Adriatic indenter, a crustal-scale wedge that caused delamination of the European lithosphere. Significant crustal roots are observed beneath the fold and thrust belts of the Alps, southeast Canadian Cordillera and parts of the southern Appalachians, but such structures beneath the northern Appalachians have probably been removed by post-orogenic collapse and/or crustal attenuation associated with the Mesozoic opening of the Atlantic Ocean.

A framework for assessing the scale of influence of environmental factors on ecological patterns

The distribution of living organisms, habitats and ecosystems is primarily driven by abiotic environmental factors that are spatially structured. Assessing the spatial structure of environmental factors, e.g., through spatial autocorrelation analyses (SAC), can thus help us understand their scale of influence on the distribution of organisms, habitats, and ecosystems. Yet SAC analyses of environmental factors are still rarely performed in biogeographic studies. Here, we describe a novel framework that combines SAC and statistical clustering to identify scales of spatial patterning of environmental factors, which can then be interpreted as the scales at which those factors influence the geographic distribution of biological and ecological features. We illustrate this new framework with datasets at different spatial or thematic resolutions. This framework is conceptually and statistically robust, providing a valuable approach to tackle a wide range of issues in ecological and environmental research and particularly when building predictors for ecological models. The new framework can significantly promote fundamental research on all spatially-structured ecological patterns. It can also foster research and application in such fields as global change ecology, conservation planning, and landscape management.

Horizontal, but not vertical, biotic interactions affect fine-scale plant distribution patterns in a low-energy system

Studies of species range determinants have traditionally focused on abiotic variables (typically climatic conditions), and therefore the recent explicit consideration of biotic interactions represents an important advance in the field. While these studies clearly support the role of biotic interactions in shaping species distributions, most examine only the influence of a single species and/or a single interaction, failing to account for species being subject to multiple concurrent interactions. By fitting species distribution models (SDMs), we examine the influence of multiple vertical (i.e., grazing, trampling, and manuring by mammalian herbivores) and horizontal (i.e., competition and facilitation; estimated from the cover of dominant plant species) interspecific interactions on the occurrence and cover of 41 alpine tundra plant species. Adding plant-plant interactions to baseline SDMs (using five field-quantified abiotic variables) significantly improved models' predictive power for independent data, while herbivore-related variables had only a weak influence. Overall, abiotic variables had the strongest individual contributions to the distribution of alpine tundra plants, with the importance of horizontal interaction variables exceeding that of vertical interaction variables. These results were consistent across three modeling techniques, for both species occurrence and cover, demonstrating the pattern to be robust. Thus, the explicit consideration of multiple biotic interactions reveals that plant-plant interactions exert control over the fine-scale distribution of vascular species that is comparable to abiotic drivers and considerably stronger than herbivores in this low-energy system.