Summary1. Functional variability (FV) of populations can be decomposed into three main features: the individual variability of multiple traits, the strength of correlations between those traits and the main direction of these correlations, the latter two being known as 'phenotypic integration'. Evolutionary biology has long recognized that FV in natural populations is key to determining potential evolutionary responses, but this topic has been little studied in functional ecology. 2. Here, we focus on the arctico-alpine perennial plant species Polygonum viviparum L.. We used a comprehensive sampling of seven functional traits in 29 wild populations covering the whole environmental niche of the species. The niche of the species was captured by a temperature gradient, which separated alpine stressful habitats from species-rich, competitive subalpine ones. We sought to assess the relative roles of abiotic stress and biotic interactions in shaping different aspects of functional variation within and among populations, that is, the multi-trait variability, the strength of correlations between traits and the main directions of functional trade-offs. 3. Populations with the highest extent of functional variability were found in the warm end of the gradient, whereas populations exhibiting the strongest degree of phenotypic integration were located in sites with intermediate temperatures. This could reveal both the importance of environmental filtering and population demography in structuring FV. Interestingly, we found that the main axes of multivariate functional variation were radically different within and across population. 4. Although the proximate causes of FV structure remain uncertain, our study presents a robust methodology for the quantitative study of functional variability in connection with species' niches. It also opens up new perspectives for the conceptual merging of intraspecific functional patterns with community ecology.