Do group dynamics affect colour morph clines during a range shift?

Funded by Strategic Research Area Biodiversity and Ecosystem Services in a Changing Climate (BECC) Lund and Gothenburg Universities Wenner-Gren Foundation EU FP7 Swedish Research Council Royal Swedish Academy of Sciences (KVA) Stiftelsen Anna-Greta and Holger Crafoords Fund Crafoord FoundationPeer reviewedPostprin

Latitudinal clines in sexual selection, sexual size dimorphism, and sex-specific genetic dispersal during a poleward range expansion

Acknowledgements This work was supported by Macquarie University (to AC, RYD), an EU FP7, Marie Curie International Incoming Fellowship (to RYD, BH; project code ‘MOVE2ADAPT’), a WennerGren Foundation Postdoctoral Stipend (to RYD, BH), the Oscar and Lili Lamm Foundation (to RYD, BH), Biodiversity and Ecosystem Services in a Changing Climate (BECC; a joint LundGothenburg University initiative) (to LL, BH), the Swedish Research Council (to EIS, BH (2014-5222, 2016-689)), the Crafoord Foundation, “Stina Werners Stiftelse” and “Erik Philip Sörensens Stiftelse” (to EIS). We thank Hanna Bensch, John Waller, Paul Caplat and Martin Andersson for field and lab assistance, the Grimsö Research Station and Mikael Åkesson for field support and Sonu Yadav for analysis advice. We thank Julian Catchen, Martin Stervander and Dag Ahren for bioinformatics advice and Maren Wellenreuther for helpful discussion.Peer reviewedPostprin

Genomic divergence and a lack of recent introgression between commercial and wild bumblebees (Bombus terrestris)

The global movement of bees for agricultural pollination services can affect local pollinator populations via hybridization. When commercial bumblebees are of the same species but of different geographic origin, intraspecific hybridization may result in beneficial integration of new genetic variation, or alternatively may disrupt locally adapted gene complexes. However, neither the existence nor the extent of genomic introgression and evolutionary divergence between wild and commercial bumblebees is fully understood. We obtained whole-genome sequencing data from wild and commercial Bombus terrestris collected from sites in Southern Sweden with and without long-term use of commercially imported B. terrestris. We search for evidence of introgression, dispersal and genome-wide differentiation in a comparative genomic analysis of wild and commercial bumblebees. Commercial B. terrestris were found in natural environments near sites where commercial bumblebees were used, as well as drifting wild B. terrestris in commercial bumblebee colonies. However, we found no evidence for widespread, recent genomic introgression of commercial B. terrestris into local wild conspecific populations. We found that wild B. terrestris had significantly higher nucleotide diversity (Nei's pi, pi), while the number of segregating sites (Watterson's theta, theta w) was higher in commercial B. terrestris. A highly divergent region on chromosome 11 was identified in commercial B. terrestris and found to be enriched with structural variants. The genes present in this region are involved in flight muscle contraction and structure and pathogen immune response, providing evidence for differing evolutionary processes operating in wild and commercial B. terrestris. We did not find evidence for recent introgression, suggesting that co-occurring commercial B. terrestris have not disrupted evolutionary processes in wild B. terrestris populations

New microsatellite markers for examining genetic variation in peripheral and core populations of the coastal giant salamander (Dicamptodon tenebrosus

Abstract The Coastal Giant Salamander (Dicamptodon tenebrosus) is classified as threatened at the northern periphery of its range in British Columbia (BC), Canada, primarily due to forestry practices and habitat fragmentation. Characterising dispersal behaviour and population connectivity is therefore a priority for this region, while genetic differentiation in core versus peripheral locations remains unstudied in this wide-ranging species. We present seven new polymorphic microsatellite markers for use in population genetic analyses of D. tenebrosus. We examine locus characteristics and genetic variation in 12 streams at the species' northern range limit in BC, and within two regions representing sub-peripheral (North Cascades) and core localities (South Cascades) in Washington State, United States. In BC, the number of alleles per locus ranged from 2-5 and observed heterozygosity ranged from 0.044-0.825. Genetic differentiation was highest between BC and the South Cascades, and intermediate between BC and the North Cascades. Across loci, mean allelic richness was similar across regions, while private allelic richness was highest in the core locality (corrected for sample size). These new microsatellite loci will be a valuable addition to existing markers for detailed landscape and population genetic analyses of D. tenebrosus across its range

Current and Historical Drivers of Landscape Genetic Structure Differ in Core and Peripheral Salamander Populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management

Genomic comparisons reveal biogeographic and anthropogenic impacts in the koala (Phascolarctos cinereus): a dietary-specialist species distributed across heterogeneous environments

The Australian koala is an iconic marsupial with highly specific dietary requirements distributed across heterogeneous environments, over a large geographic range. The distribution and genetic structure of koala populations has been heavily influenced by human actions, specifically habitat modification, hunting and translocation of koalas. There is currently limited information on population diversity and gene flow at a species-wide scale, or with consideration to the potential impacts of local adaptation. Using species-wide sampling across heterogeneous environments, and high-density genome-wide markers (SNPs and PAVs), we show that most koala populations display levels of diversity comparable to other outbred species, except for those populations impacted by population reductions. Genetic clustering analysis and phylogenetic reconstruction reveals a lack of support for current taxonomic classification of three koala subspecies, with only a single evolutionary significant unit supported. Furthermore, similar to 70% of genetic variance is accounted for at the individual level. The Sydney Basin region is highlighted as a unique reservoir of genetic diversity, having higher diversity levels (i.e., Blue Mountains region; AvHe(corr)-0.20, PL% = 68.6). Broad-scale population differentiation is primarily driven by an isolation by distance genetic structure model (49% of genetic variance), with clinal local adaptation corresponding to habitat bioregions. Signatures of selection were detected between bioregions, with no single region returning evidence of strong selection. The results of this study show that although the koala is widely considered to be a dietary-specialist species, this apparent specialisation has not limited the koala's ability to maintain gene flow and adapt across divergent environments as long as the required food source is available

Genetic and morphological divergence in island and mainland birds: Informing conservation priorities

Evolutionary processes can complicate conservation efforts for species with uncertain taxonomic classifications and discrete geographic populations. Discordant morphological and genetic patterns across the geographic range of species further calls for the identification of evolutionary significant units for conservation. Using island and mainland populations of a small Australian passerine (the superb fairy-wren, Malurus cyaneus), we examine the relationship between morphological and genetic divergence among two subspecies, M. c. ashbyi (Kangaroo Island, South Australia) and M. c. leggei (South Australia, mainland), using eight microsatellite markers. Island birds showed clear evidence for morphological divergence, with a larger body size and thinner bill compared to mainland birds. Two genetic clusters were found using Bayesian methods, comprising mainland and island regions. Estimates of recent migration rates between all sites were very low (<2%). Morphological and genetic differentiation between island and mainland sites correlated significantly, but not when controlling for isolation by distance. Genetic and morphological substructure was evident with three distinct genetic clusters in each region. Males, the highly sedentary sex, appeared to drive correlations between morphological and genetic differentiation. Our study provides evidence that the subspecies classification of M. cyaneus in island and mainland regions encapsulates two independently diverging populations that can be recognised in conservation planning.Rachael Y. Dudaniec, Beth E. Schlotfeldt, Terry Bertozzi, Stephen C. Donnellan and Sonia Kleindorferhttp://www.journals.elsevier.com/biological-conservation

Selection despite low genetic diversity and high gene flow in a rapid island invasion of the bumblebee, Bombus terrestris

Invasive species are predicted to adjust their morphological, physiological and life-history traits to adapt to their non-native environments. Although a loss of genetic variation during invasion may restrict local adaptation, introduced species often thrive in novel environments. Despite being founded by just a few individuals, Bombus terrestris (Hymenoptera: Apidae) has in less than 30 years successfully spread across the island of Tasmania (Australia), becoming abundant and competitive with native pollinators. We use RADseq to investigate what neutral and adaptive genetic processes associated with environmental and morphological variation allow B. terrestris to thrive as an invasive species in Tasmania. Given the widespread abundance of B. terrestris, we expected little genetic structure across Tasmania and weak signatures of environmental and morphological selection. We found high gene flow with low genetic diversity, although with significant isolation-by-distance and spatial variation in effective migration rates. Restricted migration was evident across the mid-central region of Tasmania, corresponding to higher elevations, pastural land, low wind speeds and low precipitation seasonality. Tajima's D indicated a recent population expansion extending from the south to the north of the island. Selection signatures were found for loci in relation to precipitation, wind speed and wing loading. Candidate loci were annotated to genes with functions related to cuticle water retention and insect flight muscle stability. Understanding how a genetically impoverished invasive bumblebee has rapidly adapted to a novel island environment provides further understanding about the evolutionary processes that determine successful insect invasions, and the potential for invasive hymenopteran pollinators to spread globally

Dealing with uncertainty in landscape genetic resistance models: a case of three co-occurring marsupials

Landscape genetics lacks explicit methods for dealing with the uncertainty in landscape resistance estimation, which is particularly problematic when sample sizes of individuals are small. Unless uncertainty can be quantified, valuable but small data sets may be rendered unusable for conservation purposes. We offer a method to quantify uncertainty in landscape resistance estimates using multimodel inference as an improvement over single model-based inference. We illustrate the approach empirically using co-occurring, woodland-preferring Australian marsupials within a common study area: two arboreal gliders (Petaurus breviceps, and Petaurus norfolcensis) and one ground-dwelling antechinus (Antechinus flavipes). First, we use maximum-likelihood and a bootstrap procedure to identify the best-supported isolation-by-resistance model out of 56 models defined by linear and non-linear resistance functions. We then quantify uncertainty in resistance estimates by examining parameter selection probabilities from the bootstrapped data. The selection probabilities provide estimates of uncertainty in the parameters that drive the relationships between landscape features and resistance. We then validate our method for quantifying uncertainty using simulated genetic and landscape data showing that for most parameter combinations it provides sensible estimates of uncertainty. We conclude that small data sets can be informative in landscape genetic analyses provided uncertainty can be explicitly quantified. Being explicit about uncertainty in landscape genetic models will make results more interpretable and useful for conservation decision-making, where dealing with uncertainty is critical

P.brev_GenepopInput

Microsatellite datafile in genepop format - P.brevicep