Identifying sources of genetic variation and reconstructing invasion routes for non-native introduced species is central to understanding the circumstances under which they may evolve increased invasiveness. In this study, we used genome-wide single nucleotide polymorphisms to study the colonization history of Centaurea solstitialis in its native range in Eurasia and invasions into the Americas. We leveraged this information to pinpoint key evolutionary shifts in plant size, a focal trait associated with invasiveness in this species. Our analyses revealed clear population genomic structure of potential source populations in Eurasia, including deep differentiation of a lineage found in the southern Apennine and Balkan Peninsulas and divergence among populations in Asia, eastern Europe, and western Europe. We found strongest support for an evolutionary scenario in which western European populations were derived from an ancient admixture event between populations from eastern Europe and Asia, and subsequently served as the main genetic ‘bridgehead’ for introductions to the Americas. Introductions to California appear to be from a single source region, and multiple, independent introductions of divergent genotypes likely occurred into the Pacific Northwest. Plant size has evolved significantly at three points during range expansion, including a large size increase in the lineage responsible for the aggressive invasion of California’s interior. These results reveal a long history of colonization, admixture, and trait evolution in C. solstitialis, and suggest routes for improving evidence-based management decisions for one of the most ecologically and economically damaging invasive species in the western United States.

The G-matrix occupies an important position in evolutionary biology both as a summary of the inheritance of quantitative traits and as an ingredient in predicting how those traits will respond to selection and drift. Consequently, the stability of G has an important bearing on the accuracy of predicted evolutionary trajectories. Furthermore, G should evolve in response to stable features of the adaptive landscape and their trajectories through time. Although the stability and evolution of G might be predicted from knowledge of selection in natural populations, most empirical comparisons of G-matrices have been made in the absence of such a priori predictions. We present a theoretical argument that within-sex G-matrices should be more stable than between-sex B-matrices because they are more powerfully exposed to multivariate stabilizing selection. We tested this conjecture by comparing estimates of B-and within-sex G-matrices among three populations of the garter snake Thamnophis elegans. Matrix comparisons using Flury's hierarchical approach revealed that within-sex G-matrices had four principal components in common (full CPC), whereas B-matrices had only a single principal component in common and eigenvalues that were more variable among populations. These results suggest that within-sex G is more stable than B, as predicted by our theoretical argument. K E Y W O R D S : Flury hierarchy, genetic covariance matrix, sexual dimorphism, Thamnophis elegans.The additive genetic variance-covariance matrix, or G-matrix, plays a central role in evolutionary theory in predicting deterministic responses to selection as well as the stochastic consequences of finite population size (Lande 1979). In particular, the deterministic evolution of the phenotypic mean is affected both by selection and the G-matrix. Unfortunately, however, we lack an analytical framework for predicting how G itself will evolve in response to selection in populations of finite size. Likewise, no equations have been derived that specify the stability of G in evolving populations. In the absence of an analytical framework, our knowledge of how G responds to selection and finite population size is based on simulation and empirical studies .Correlational selection should have profound consequences for G-matrix evolution and stability. Indeed, simulations have shown that G stability depends on trait-specific issues of mutation and selection, as well as on population size. As expected from theoretical considerations (Lande 1980a), the stability of the principal components of G (i.e., the angle of its eigenvectors) is enhanced by strong correlational selection, as well as by large population size (Jones et al.

Quaternary climatic oscillations caused changes in sea level that altered the size, number and degree of isolation of islands, particularly in land-bridge archipelagoes. Elucidating the demographic effects of these oscillations increases our understanding of the role of climate change in shaping evolutionary processes in archipelagoes. The Puerto Rican Bank (PRB) (Puerto Rico and the Eastern Islands, which comprise Vieques, Culebra, the Virgin Islands and associated islets) in the eastern Caribbean Sea periodically coalesced during glaciations and fragmented during interglacial periods of the quaternary. To explore population-level consequences of sea level changes, we studied the phylogeography of the frog Eleutherodactylus antillensis across the archipelago. We tested hypotheses encompassing vicariance and dispersal narratives by sequencing mtDNA (c. 552 bp) of 285 individuals from 58 localities, and four nuDNA introns (totalling c. 1633 bp) from 173 of these individuals. We found low support for a hypothesis of divergence of the Eastern Islands populations prior to the start of the penultimate interglacial c. 250 kya, and higher support for a hypothesis of colonization of the Eastern Islands from sources in eastern Puerto Rico during the penultimate and last glacial period, when a land bridge united the PRB. The Río Grande de Loíza Basin in eastern Puerto Rico delineates a phylogeographic break. Haplotypes shared between the PRB and St. Croix (an island c. 105 km south-east of this archipelago) likely represent human-mediated introductions. Our findings illustrate how varying degrees of connectivity and isolation influence the evolution of tropical island organisms.

The margins of an expanding range are predicted to be challenging environments for adaptation. Marginal populations should often experience low effective population sizes (Ne) where genetic drift is high due to demographic expansion and/or census population size is low due to unfavourable environmental conditions. Nevertheless, invasive species demonstrate increasing evidence of rapid evolution and potential adaptation to novel environments encountered during colonization, calling into question whether significant reductions in Ne are realized during range expansions in nature. Here we report one of the first empirical tests of the joint effects of expansion dynamics and environment on effective population size variation during invasive range expansion. We estimate contemporary values of Ne using rates of linkage disequilibrium among genome‐wide markers within introduced populations of the highly invasive plant Centaurea solstitialis (yellow starthistle) in North America (California, USA), and within native Eurasian populations. As predicted, we find that Ne within the invaded range is positively correlated with both expansion history (time since founding) and habitat quality (abiotic climate). History and climate had independent additive effects with similar effect sizes, indicating an important role for both factors in this invasion. These results support theoretical expectations for the population genetics of range expansion, though whether these processes can ultimately arrest the spread of an invasive species remains an unanswered question.