Multifactorial studies assessing the cumulative effects of natural and anthropogenic stressors on individual stress response are crucial to understand how organisms and populations cope with environmental change. We tested direct and indirect causal pathways through which environmental stressors affect the stress response of wild gilthead seabream in Mediterranean costal lagoons using an integrative PLS-PM approach. We integrated information on 10 environmental variables and 36 physiological variables into seven latent variables reflecting lagoons features and fish health. These variables concerned fish lipid reserves, somatic structure, inorganic contaminant loads, and individual trophic and stress response levels. This modelling approach allowed explaining 30 % of the variance within these 46 variables considered. More importantly, 54 % of fish stress response was explained by the dependent lagoon features, fish age, fish diet, fish reserve, fish structure and fish contaminant load latent variables included in our model. This integrative study sheds light on how individuals deal with contrasting environments and multiple ecological pressures.