Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

Climate change is impacting locally adapted species such as the keystone tree species cork oak (Quercus suber L.). Quantifying the importance of environmental variables in explaining the species distribution can help build resilient populations in restoration projects and design forest management strategies. Using landscape genomics, we investigated the population structure and ecological adaptation of this tree species across the Mediterranean Basin. We applied genotyping by sequencing and derived 2,583 single nucleotide polymorphism markers genotyped from 81 individuals across 17 sites in the studied region. We implemented an approach based on the nearest neighbour haplotype ‘coancestry’ and uncovered a weak population structure along an east–west climatic gradient across the Mediterranean region. We identified genomic regions potentially involved in local adaptation and predicted differences in the genetic composition across the landscape under current and future climates. Variants associated with temperature and precipitation variables were detected, and we applied a nonlinear multivariate association method, gradient forest, to project these gene–environment relationships across space. The model allowed the identification of geographic areas within the western Mediterranean region most sensitive to climate change: south‐western Iberia and northern Morocco. Our findings provide a preliminary assessment towards a potential management strategy for the conservation of cork oak in the Mediterranean Basin.

We used a population genomic approach to unravel the population structure, genetic differentiation, and genetic diversity of three widespread wild bee species across the Iberian Peninsula, Andrena agilissima, Andrena flavipes and Lasioglossum malachurum. Our results demonstrated that genetic lineages in the Ebro River valley or near the Pyrenees mountains are different from the rest of Iberia. This relatively congruent pattern across species once more supports the hypothesis of “refugia within refugia” in the Iberian Peninsula. The results for A. flavipes and A. agilissima showed an unexpected pattern of genetic differentiation, with the generalist polylectic A. flavipes having lower levels of genetic diversity (Ho = 0.0807, He = 0.2883) and higher differentiation (FST = 0.5611), while the specialist oligolectic A. agilissima had higher genetic diversity (Ho = 0.2104, He = 0.3282) and lower differentiation values (FST = 0.0957). For L. malachurum, the smallest and the only social species showed the lowest inbreeding coefficient (FIS = 0.1009) and the lowest differentiation level (FST = 0.0663). Overall, our results, suggest that this pattern of population structure and genetic diversity could be explained by the combined role of past climate changes and the life-history traits of the species (i.e., size, sociality and host-plant specialization), supporting the role of the Iberian refugia as a biodiversity hotspot.