27Parallel adaptation results from independent evolution of similar traits between closely related 28 lineages and allows testing to which extent evolution is repeatable. Parallel adaptation often 29 involves similar gene expression changes but the identity of genes shaped by parallel selection 30 and the causes of expression parallelism remains largely unknown. By comparing genomes 31 and transcriptomes of four independent foothill-alpine population pairs across four 32 treatments, we addressed genetic basis, plasticity and functional consequences of gene 33 expression parallelism in alpine adaptation. Seeds of four population pairs of Arabidopsis 34 arenosa from distinct mountain regions were raised under four treatments that differed in 35 temperature and irradiance, factors varying strongly with elevation. Parallelism in gene 36 expression was quantified by RNA-seq in leaves of young plants. By manipulating 37temperature and irradiance, we also tested for parallelism in plasticity (gene-by-environment 38 interaction [GEI]). We found significant parallelism in differential gene expression across 39 four independently recruited alpine ecotypes with an overrepresentation of genes involved in 40 biotic stress response. In addition, we demonstrated significant parallelism in GEI indicating 41 shared response to environmental variation in our foothill and alpine populations. Fraction of 42 genes showing expression parallelism also encompassed genomic variants showing outlying 43 differentiation, with greater enrichment of such variants in cis-regulatory elements. In 44 summary, our results suggest frequent evolutionary repeatability in shaping expression 45 difference associated with colonization of an alpine environment and support the hypothesis 46 of an important role of genetic variation in cis-regulatory elements in gene expression 47 parallelism.48 49