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.

Vicariance is the simplest explanation for divergence between sister lineages separated by a potential barrier, and the northern Andes would seem to provide an ideal example of a vicariant driver of divergence. We evaluated the potential role of the uplift of the Eastern Cordillera (EC) of the Colombian Andes and the Mérida Andes (MA) of Venezuela as drivers of vicariance between lowland populations co-distributed on both flanks. We synthesized published geological data and provided a new reconstruction showing that the EC-MA grew from north to south, reaching significant heights and separating drainages and changing sediment composition by 38–33 million years ago (Ma). A few lowland passes across the EC-MA may have reached their current heights (~1,900 m a.s.l.) at 3–5 Ma. We created a comparative phylogeographic data set for 37 lineages of lowland tetrapods. Based on molecular phylogenetic analyses, most divergences between sister populations or species across the EC-MA occurred during Pliocene and the Quaternary and a few during the latest Miocene, and coalescent simulations rejected synchronous divergence for most groups. Divergence times were on average slightly but significantly more recent in homeotherms relative to poikilotherms. Because divergence ages are mostly too recent relative to the geological history and too asynchronous relative to each other, divergence across the northern Andes may be better explained by organism-environment interactions concomitant with climate oscillations during the Pleistocene, and/or dispersal across portals through the Andes.

AimTo evaluate the potential role of the orogeny of the Eastern Cordillera (EC) of the Colombian Andes and the Mérida Andes (MA) of Venezuela as drivers of vicariance between populations of 37 tetrapod lineages co-distributed on both flanks, through geological reconstruction and comparative phylogeographic analyses.LocationNorthwestern South AmericaMethodsWe first reviewed and synthesized published geological data on the timing of uplift for the EC-MA. We then combined newly generated mitochondrial DNA sequence data with published datasets to create a comparative phylogeographic dataset for 37 independent tetrapod lineages. We reconstructed time-calibrated molecular phylogenies for each lineage under Bayesian inference to estimate divergence times between lineages located East and West of the Andes. We performed a comparative phylogeographic analysis of all lineages within each class of tetrapod using hierarchical approximate Bayesian computation (hABC) to test for synchronous vicariance across the EC-MA. To evaluate the potential role of life history in explaining variation in divergence times among lineages, we evaluated 13 general linear models (GLM) containing up to six variables each (maximum elevation, range size, body length, thermoregulation, type of dispersal, and taxonomic class).ResultsOur synthesis of geological evidence suggested that the EC-MA reached significant heights by 38–33 million years ago (Ma) along most of its length, and we reject the oft-cited date of 2–5 Ma. Based on mtDNA divergence from 37 lineages, however, the median estimated divergence time across the EC-MA was 3.26 Ma (SE = 2.84) in amphibians, 2.58 Ma (SE = 1.81) in birds, 2.99 Ma (SE = 4.68) in reptiles and 1.43 Ma (SE = 1.23) in mammals. Using Bayes Factors, the hypothesis for a single temporal divergence interval containing synchronous divergence events was supported for mammals and but not supported for amphibians, non-avian reptiles, or birds. Among the six life-history variables tested, only thermoregulation successfully explained variation in divergence times (minimum AICc, R2 0.10), with homeotherms showing more recent divergence relative to poikilotherms.Main conclusionsOur results reject the hypothesis of the rise Andean Cordillera as driver of vicariance of lowland population because divergence dates are too recent and too asynchronous. We discuss alternative explanations, including dispersal through mountain passes, and suggest that changes in the climatic conditions during the Pliocene and Pleistocene interacted with tetrapod physiology, promoting older divergences in amphibians and reptiles relative to mammals and birds on an already established orogen.