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.

The spatial signature of microevolutionary processes structuring genetic variation may play an important role in the detection of loci under selection. However, the spatial location of samples has not yet been used to quantify this. Here, we present a new two-step method of spatial outlier detection at the individual and deme levels using the power spectrum of Moran eigenvector maps (MEM). The MEM power spectrum quantifies how the variation in a variable, such as the frequency of an allele at a SNP locus, is distributed across a range of spatial scales defined by MEM spatial eigenvectors. The first step (Moran spectral outlier detection: MSOD) uses genetic and spatial information to identify outlier loci by their unusual power spectrum. The second step uses Moran spectral randomization (MSR) to test the association between outlier loci and environmental predictors, accounting for spatial autocorrelation. Using simulated data from two published papers, we tested this two-step method in different scenarios of landscape configuration, selection strength, dispersal capacity and sampling design. Under scenarios that included spatial structure, MSOD alone was sufficient to detect outlier loci at the individual and deme levels without the need for incorporating environmental predictors. Follow-up with MSR generally reduced (already low) false-positive rates, though in some cases led to a reduction in power. The results were surprisingly robust to differences in sample size and sampling design. Our method represents a new tool for detecting potential loci under selection with individual-based and population-based sampling by leveraging spatial information that has hitherto been neglected.

A wealth of plant species used by humans for different purposes, but mainly as food, originated and domesticated in the Mesoamerican region. Papaya (Carica papaya) is the third most cultivated tropical crop worldwide, and it has been hypothesized that Mesoamerica is the most likely center of its origin and domestication. In support of it, many wild populations of papaya occur throughout Mesoamerica and hence represent the gene pool of genetic variability for further evolution and future crop management. Despite its importance, a dearth of information exists regarding the status of wild populations of papaya, as compared to the extent of knowledge, and interest, on domesticated varieties. We review the evidence on the extant wild populations of papaya, as well as its origin and distribution. Also, we synthetize what is known on the domestication history of the species, including the domestication syndrome that distinguishes wild and domesticated papayas. Moreover, we make an account of the use of genetic markers to assess genetic diversity of wild and domesticated papaya, and discuss the importance of papaya as the first species with a transgenic cultivar to be released for human consumption, and one that has its complete genome sequenced. Evidence from different disciplines strongly suggest that papaya originated and was domesticated in Mesoamerica, and that wild populations in the region possess, still, high genetic diversity compared to the domesticated papaya. Finally, we outline papaya as an excellent model species for genomic studies that will help gain insight into the domestication process and improvement of papaya and other tropical crops.

• Premise of the study: Wild populations of domesticated species constitute a genetic reservoir and are fundamental to the evolutionary potential of species. Wild papaya (Carica papaya) is a rare, short-lived, gap-colonizing, dioecious tree that persists in the forest by continuous dispersal. Theoretically, these life-history characteristics render wild papaya highly susceptible to habitat fragmentation, with anticipated negative effects on its gene pool. Further, species dioecy may cause founder effects to generate local biases in sex ratio, decreasing effective population size.• Methods: We contrasted the genetic diversity and structure of C. papaya between wild populations from rainforest fragments and continuous forest at Los Tuxtlas, Mexico. We evaluated recent migration rates among populations as well as landscape resistance to gene flow. Finally, we calculated the sex ratio of the populations in both habitats.• Key results: Populations of wild papaya in rainforest fragments showed lower genetic diversity and higher population differentiation than populations in continuous rainforest. Estimates of recent migration rates showed a higher percentage of migrants moving from the continuous forest to the forest fragments than in the opposite direction. Agricultural land and cattle pasture were found to be the most resistant matrices to gene flow. Finally, biased sex ratios were seen to affect the effective population size in both habitats.• Conclusions: The mating system, rarity, and short life cycle of C. papaya are exacerbating the effects of rainforest fragmentation on its genetic diversity, threatening the persistence of its natural populations in the proposed place of origin as well as its genetic reservoir.

Background and aims Few studies have evaluated the genetic structure and evolutionary history of wild varieties of important crop species. The wild papaya (Carica papaya) is a key element of early successional tropical and sub-tropical forests in Mexico, and constitutes the genetic reservoir for evolutionary potential of the species. In this study we aimed to determine how diverse and structured is the genetic variability of wild populations of C. papaya in Northern Mesoamerica. Moreover, we assessed if genetic structure and evolutionary history coincide with hypothetized (1) pre-Pleistocene events (Isthmus of Tehuantepec sinking), (2) Pleistocene refugia or (3) recent patterns.Methods We used six nuclear and two chloroplast (cp) DNA markers to assess the genetic diversity and phylogeographical structure of 19 wild populations of C. papaya in its natural distribution in Northern Mesoamerica.Key Results We found high genetic diversity (H o ¼ 0Á681 for nuclear markers, and h ¼ 0Á701 for cpDNA markers) and gene flow between populations of C. papaya (migration r up to 420 km). A lack of phylogeographical structure was found with the cpDNA markers (NST < GST), whereas a recent population structure was inferred with the nuclear markers. Evidence indicates that pre-Pleistocene events or refugia did not play an important role in the genetic structuring of wild papaya.Conclusions Because of its life history characteristics and lack of an ancient phylogeographical structure found with the cpDNA markers, we suggest that C. papaya was dispersed throughout the lowland rain forests of Mexico (along the coastal plains and foothills of Sierras). This scenario supports the hypothesis that tropical forests in Northern Mesoamerica did not experience important climate fluctuations during the Pleistocene, and that the life history of C. papaya could have promoted long-distance dispersal and rapid colonization of lowland rainforests. Moreover, the results obtained with the nuclear markers suggest recent human disturbances. The fragmentation of tropical habitats in Northern Mesoamerica appears to be the main driver of genetic structuring, and the major threat to the dispersion and survival of the species in the wild.

No abstract

IntroductionTheory predicts that habitat fragmentation, by reducing population size and increasing isolation among remnant populations, can alter their genetic diversity and structure. A cascade of effects is expected: genetic drift and inbreeding after a population bottleneck, changes in biotic interactions that may affect, as in the case of plants, pollen dynamics, mating system, reproductive success. The detection of the effects of contemporary habitat fragmentation on the genetic structure of populations are conditioned by the magnitude of change, given the few number of generations since the onset of fragmentation, especially for long-lived organisms. However, the present-day genetic structure of populations may bear the signature of past demography events. Here, we examine the effects of rainforest fragmentation on the genetic diversity, population structure, mating system (outcrossing rate), indirect gene flow and contemporary pollen dynamics in the understory herb Aphelandra aurantiaca. Also, we assessed its present-day genetic structure under different past demographic scenarios.MethodsTwelve populations of A. aurantiaca were sampled in large (4), medium (3), and small (5) forest fragments in the lowland tropical rainforest at Los Tuxtlas region. Variation at 11 microsatellite loci was assessed in 28–30 reproductive plants per population. In two medium- and two large-size fragments we estimated the density of reproductive plants, and the mating system by analyzing the progeny of different mother plants per population.ResultsDespite prevailing habitat fragmentation, populations of A. aurantiaca possess high genetic variation (He = 0.61), weak genetic structure (Rst = 0.037), and slight inbreeding in small fragments. Effective population sizes (Ne) were large, but slightly lower in small fragments. Migrants derive mostly from large and medium size fragments. Gene dispersal is highly restricted but long distance gene dispersal events were detected. Aphelandra aurantiaca shows a mixed mating system (tm = 0.81) and the outcrossing rate have not been affected by habitat fragmentation. A strong pollen pool structure was detected due to few effective pollen donors (Nep) and low distance pollen movement, pointing that most plants received pollen from close neighbors. Past demographic fluctuations may have affected the present population genetic structure as Bayesian coalescent analysis revealed the signature of past population expansion, possibly during warmer conditions after the last glacial maximum.DiscussionHabitat fragmentation has not increased genetic differentiation or reduced genetic diversity of A. aurantiaca despite dozens of generations since the onset of fragmentation in the region of Los Tuxtlas. Instead, past population expansion is compatible with the lack of observed genetic structure. The predicted negative effects of rainforest fragmentation on genetic diversity and population structure of A. aurantiaca seem to have been buffered owing to its large effective populations and long-distance disper...

We evaluated the role of gene flow and wild-crop introgression on the structure and genetic diversity of Lima bean (Phaseolus lunatus) in the Yucatan Peninsula, an important Mesoamerican diversity area for this crop, using a genotyping-by-sequencing approach (15,168 SNP markers) and two scales. At the local scale, STRUCTURE and NGSEP analyses showed predominantly crop-to-wild introgression, but also evidence of a bidirectional gene flow in the two wild-weedy-crop complexes studied (Itzinté and Dzitnup). The ABBA-BABA tests showed a higher introgression in Itzinté (the older complex) than in Dzitnup (the younger one); at the allelic level, the wild-crop introgression in Itzinté was similar in both directions, in Dzitnup it was higher from crop-to-wild; and at the chromosomal level, introgression in Itzinté was from wild-to-crop, whereas in Dzitnup it occured in the opposite direction. Also, we foundHEvalues slightly higher in the domesticated accessions than in the wild ones, in both complexes (Itzinté: wild = 0.31, domesticated = 0.34; Dzinup: wild = 0.27, domesticated = 0.36), but %Pandπestimators were higher in the wild accessions than in the domesticated ones. At a regional scale, STRUCTURE and MIGRATE showed a low gene flow, predominantly from crop-to-wild; and STRUCTURE, Neighbor-Joining and PCoA analyses indicated the existence of two wild groups and one domesticated group, with a marked genetic structure based in the existence of domesticated MI and wild MII gene pools. Also, at the regional scale, we found a higher genetic diversity in the wild accessions than in the domesticated ones, in all estimators used (e.g.,HE= 0.27 andHE= 0.17, respectively). Our results indicate that gene flow and introgression are playing an important role at the local scale, but its consequences on the structure and genetic diversity of the Lima bean are not clearly reflected at the regional scale, where diversity patterns between wild and domesticated populations could be reflecting historical events.