In this research, the Bayesian quantile regression model is applied to investigate the teleconnections between large oceanic–atmospheric indices and drought standardized precipitation index (SPI) in Iran. The 12-month SPI time series from 138 synoptic stations for 1952–2014 were selected as the drought index. Three oceanic–atmospheric indices, the North Atlantic Oscillation (NAO), the Southern Oscillation Index (SOI) and the Multivariate El Niño/Southern Oscillation Index (MEI), were selected as covariates. The results show that NAO has the weakest impact on drought in different quantiles and different regions in Iran. La Niña conditions amplified droughts through all SPI quantiles in western, Caspian Sea coastal regions and southern regions. The positive phase of MEI significantly modulates low SPI quantiles (i.e. drought conditions) throughout the Zagros region, Caspian Sea coastal regions and southern regions. The study shows that the effect of large oceanic–atmospheric indices have heterogeneous impacts on extreme dry and wet conditions.