Molecular Ecology

Di Lecce, Irene; Sudyka, Joanna; Perrier, Charles; Szulkin, Marta

doi: 10.1111/mec.17481pmid: 39044486

Urbanisation has been increasing worldwide in recent decades, driving environmental change and exerting novel selective pressures on wildlife. Phenotypic differences between urban and rural individuals have been widely documented in several taxa. However, the extent to which urbanisation impacts mating strategies is less known. Here, we investigated extra‐pair paternity variation in great tits (Parus major) and blue tits (Cyanistes caeruleus) breeding in nestboxes set in a gradient of urbanisation in Warsaw, Poland, over three breeding seasons. Urbanisation was quantified as the amount of light pollution, noise pollution, impervious surface area (ISA) and tree cover within a 100‐m radius around each nestbox. We obtained genotypes for 1213 great tits at 7344 SNP markers and for 1299 blue tits at 9366 SNP markers with a genotyping‐by‐sequencing method, and inferred extra‐pair paternity by computing a genomewide relatedness matrix. We report higher extra‐pair paternity in blue tits breeding in more urbanised areas, for example, with higher light pollution and ISA, and lower tree cover. However, no such trend was found in great tits. Late‐stage survival of individual nestlings in both species was not associated with paternity or urbanisation proxies, thus we were not able to detect fitness benefits or drawbacks of being an extra‐pair offspring in relation to urbanisation. Our results contribute to the growing body of knowledge reporting on the effects of urbanisation on avian ecology and behaviour, and confirm species‐specific and population‐specific patterns of extra‐pair paternity variation.

Giglio, Rachael M.; Bowden, Courtney F.; Brook, Ryan K.; Piaggio, Antoinette J.; Smyser, Timothy J.

doi: 10.1111/mec.17489pmid: 39148259

Globalization has led to the frequent movement of species out of their native habitat. Some of these species become highly invasive and capable of profoundly altering invaded ecosystems. Feral swine (Sus scrofa × domesticus) are recognized as being among the most destructive invasive species, with populations established on all continents except Antarctica. Within the United States (US), feral swine are responsible for extensive crop damage, the destruction of native ecosystems, and the spread of disease. Purposeful human‐mediated movement of feral swine has contributed to their rapid range expansion over the past 30 years. Patterns of deliberate introduction of feral swine have not been well described as populations may be established or augmented through small, undocumented releases. By leveraging an extensive genomic database of 18,789 samples genotyped at 35,141 single nucleotide polymorphisms (SNPs), we used deep neural networks to identify translocated feral swine across the contiguous US. We classified 20% (3364/16,774) of sampled animals as having been translocated and described general patterns of translocation using measures of centrality in a network analysis. These findings unveil extensive movement of feral swine well beyond their dispersal capabilities, including individuals with predicted origins >1000 km away from their sampling locations. Our study provides insight into the patterns of human‐mediated movement of feral swine across the US and from Canada to the northern areas of the US. Further, our study validates the use of neural networks for studying the spread of invasive species.

Croft, Lilly; Matheson, Paige; Butterworth, Nathan J.; McGaughran, Angela

doi: 10.1111/mec.17492pmid: 39136044

Invasive species often undergo demographic bottlenecks that cause a decrease in genetic diversity and associated reductions in population fitness. Despite this, they manage to thrive in novel environments. Investigating the effects of inbreeding and genetic bottlenecks on population fitness for invasive species is, therefore, key to understanding how they may survive in new environments. We used the blowfly Calliphora vicina (Sciences, Mathématiques et Physique, 1830, 2, 1), which is native to Europe and was introduced to Australia and New Zealand, to examine the effects of genetic diversity on population fitness. We first collected 59 samples from 15 populations across New Zealand and one in Australia, and used 20,501 biallelic SNPs to investigate population genomic diversity, structure and admixture. We then explored the impacts of repeated experimental bottlenecks on population fitness by creating inbred and outbred lines of C. vicina and measuring a variety of fitness traits. In wild‐caught samples, we found low overall genetic diversity, signals of genetic admixture and limited (<3%) genetic differentiation between North and South Island populations, with genetic links between the South Island and Australia. Following experimental bottlenecks, we found significant reductions in fitness for inbred lines. However, fitness effects were not felt equally across all phenotypic traits. Moreover, they were not enough to cause population collapse in any experimental line, suggesting that C. vicina (when under relaxed selection, as in laboratory settings) may be able to compensate for population bottlenecks even when highly inbred. Our results demonstrate the value of a tractable experimental system for investigating processes that may facilitate or hamper biological invasion.

Pallarés, Susana; Ortego, Joaquín; Carbonell, José Antonio; Franco‐Fuentes, Eduardo; Bilton, David T.; Millán, Andrés; Abellán, Pedro

doi: 10.1111/mec.17487pmid: 39108249

An intricate interplay between evolutionary and demographic processes has frequently resulted in complex patterns of genetic and phenotypic diversity in alpine lineages, posing serious challenges to species delimitation and biodiversity conservation planning. Here we integrate genomic data, geometric morphometric analyses and thermal tolerance experiments to explore the role of Pleistocene climatic changes and adaptation to alpine environments on patterns of genomic and phenotypic variation in diving beetles from the taxonomically complex Agabus bipustulatus species group. Genetic structure and phylogenomic analyses revealed the presence of three geographically cohesive lineages, two representing trans‐Palearctic and Iberian populations of the elevation‐generalist A. bipustulatus and another corresponding to the strictly‐alpine A. nevadensis, a narrow‐range endemic taxon from the Sierra Nevada mountain range in southeastern Iberia. The best‐supported model of lineage divergence, along with the existence of pervasive genetic introgression and admixture in secondary contact zones, is consistent with a scenario of population isolation and connectivity linked to Quaternary climatic oscillations. Our results suggest that A. nevadensis is an alpine ecotype of A. bipustulatus, whose genotypic, morphological and physiological differentiation likely resulted from an interplay between population isolation and local altitudinal adaptation. Remarkably, within the Iberian Peninsula, such ecotypic differentiation is unique to Sierra Nevada populations and has not been replicated in other alpine populations of A. bipustulatus. Collectively, our study supports fast ecotypic differentiation and incipient speciation processes within the study complex and points to Pleistocene glaciations and local adaptation along elevational gradients as key drivers of biodiversity generation in alpine environments.

Kuo, Wen‐Hsi; Zhong, Limei; Wright, Sara J.; Goad, David M.; Olsen, Kenneth M.

doi: 10.1111/mec.17484pmid: 39072878

Species that repeatedly evolve phenotypic clines across environmental gradients have been highlighted as ideal systems for characterizing the genomic basis of local environmental adaptation. However, few studies have assessed the importance of observed phenotypic clines for local adaptation: conspicuous traits that vary clinally may not necessarily be the most critical in determining local fitness. The present study was designed to fill this gap, using a plant species characterized by repeatedly evolved adaptive phenotypic clines. White clover is naturally polymorphic for its chemical defence cyanogenesis (HCN release with tissue damage); climate‐associated cyanogenesis clines have evolved throughout its native and introduced range worldwide. We performed landscape genomic analyses on 415 wild genotypes from 43 locations spanning much of the North American species range to assess the relative importance of cyanogenesis loci vs. other genomic factors in local climatic adaptation. We find clear evidence of local adaptation, with temperature‐related climatic variables best describing genome‐wide differentiation between sampling locations. The same climatic variables are also strongly correlated with cyanogenesis frequencies and gene copy number variations (CNVs) at cyanogenesis loci. However, landscape genomic analyses indicate no significant contribution of cyanogenesis loci to local adaptation. Instead, several genomic regions containing promising candidate genes for plant response to seasonal cues are identified — some of which are shared with previously identified QTLs for locally adaptive fitness traits in North American white clover. Our findings suggest that local adaptation in white clover is likely determined primarily by genes controlling the timing of growth and flowering in response to local seasonal cues. More generally, this work suggests that caution is warranted when considering the importance of conspicuous phenotypic clines as primary determinants of local adaptation.

doi: 10.1111/mec.17016pmid: N/A

Zhou, Qiujie; Karunarathne, Piyal; Andersson‐Li, Lili; Chen, Chen; Opgenoorth, Lars; Heer, Katrin; Piotti, Andrea; Vendramin, Giovanni Giuseppe; Nakvasina, Elena; Lascoux, Martin; Milesi, Pascal

doi: 10.1111/mec.17495pmid: 39148357

Most tree species underwent cycles of contraction and expansion during the Quaternary. These cycles led to an ancient and complex genetic structure that has since been affected by extensive gene flow and by strong local adaptation. The extent to which hybridization played a role in this multi‐layered genetic structure is important to be investigated. To study the effect of hybridization on the joint population genetic structure of two dominant species of the Eurasian boreal forest, Picea abies and P. obovata, we used targeted resequencing and obtained around 480 K nuclear SNPs and 87 chloroplast SNPs in 542 individuals sampled across most of their distribution ranges. Despite extensive gene flow and a clear pattern of Isolation‐by‐Distance, distinct genetic clusters emerged, indicating the presence of barriers and corridors to migration. Two cryptic refugia located in the large hybrid zone between the two species played a critical role in shaping their current distributions. The two species repeatedly hybridized during the Pleistocene and the direction of introgression depended on latitude. Our study suggests that hybridization helped both species to overcome main shifts in their distribution ranges during glacial cycles and highlights the importance of considering whole species complex instead of separate entities to retrieve complex demographic histories.

Augustijnen, Hannah; Lucek, Kay

doi: 10.1111/mec.17488pmid: 39119885

Replicated secondary contact zones can provide insights into the barriers to gene flow that are important during speciation and can reveal to which degree secondary contact may result in similar evolutionary outcomes. Here, we studied two secondary contact zones between highly differentiated Alpine butterflies of the genus Erebia using whole‐genome resequencing data. We assessed the genomic relationships between populations and species and found hybridization to be rare, with no to little current or historical introgression in either contact zone. There are large similarities between contact zones, consistent with an allopatric origin of interspecific differentiation, with no indications for ongoing reinforcing selection. Consistent with expected reduced effective population size, we further find that scaffolds related to the Z‐chromosome show increased differentiation compared to the already high levels across the entire genome, which could also hint towards a contribution of the Z chromosome to species divergence in this system. Finally, we detected the presence of the endosymbiont Wolbachia, which can cause reproductive isolation between its hosts, in all E. cassioides, while it appears to be fully or largely absent in contact zone populations of E. tyndarus. We discuss how this rare pattern may have arisen and how it may have affected the dynamics of speciation upon secondary contact.

Gómez Quijano, María José; Gross, Briana L.; Etterson, Julie R.

doi: 10.1111/mec.17483pmid: 39056407

Adaptive differentiation of traits and underlying loci can occur at a small geographical scale if natural selection is stronger than countervailing gene flow and drift. We investigated this hypothesis using coupled quantitative genetic and genomic approaches for a wind‐pollinated tree species, Quercus rubra, along the steep, narrow gradient of the Lake Superior coast that encompasses four USDA Hardiness Zones within 100 km. For the quantitative genetic component of this study, we examined phenotypic differentiation among eight populations in a common garden, measuring seed mass, germination, height, stem diameter, leaf number, specific leaf area and survival. For the genomic component, we quantified genetic differentiation for 26 populations from the same region using RAD‐seq. Because hybridisation with Quercus ellipsoidalis occurs in other parts of the species' range, we included two populations of this congener for comparison. In the common garden study, we found a strong signal of population differentiation that was significantly associated with at least one climate factor for nine of 10 measured traits. In contrast, we found no evidence of genomic differentiation among populations based on FST or any other measures. However, both distance‐based and genotype–environment association analyses identified loci showing the signature of selection, with one locus in common across five analyses. This locus was associated with the minimum temperature of the coldest month, a factor that defines the climate zones and was also significant in the common garden analyses. In addition, we documented introgression from Q. ellipsoidalis into Q. rubra, with rates of introgression correlated with the climate gradient. In sum, this study reveals signatures of selection at the quantitative trait and genomic level consistent with climate adaptation, a pattern that is more often documented at a much broader geographical scale, especially in long‐lived wind‐pollinated species.

Vidal‐Villarejo, Mireia; Dößelmann, Bianca; Kogler, Benedikt; Hammerschmidt, Michael; Oppliger, Barbara; Oppliger, Hans; Schmid, Karl

doi: 10.1111/mec.17482pmid: 39082382

The spread and adaptation of fungal plant pathogens in agroecosystems are facilitated by environmental homogeneity. Metagenomic sequencing of infected tissues allowed us to monitor eco‐evolutionary dynamics and interactions between host, pathogen and plant microbiome. Exserohilum turcicum, the causal agent of northern corn leaf blight (NCLB) in maize, is distributed in multiple clonal lineages throughout Europe. To characterize regional pathogen diversity, we conducted metagenomic DNA sequencing on 241 infected leaf samples from the highly susceptible Swiss maize landrace Rheintaler Ribelmais, collected over 3 years (2016–2018) from an average of 14 agricultural farms within the Swiss Rhine Valley. All major European clonal lineages of E. turcicum were identified. Lineages differ by their mating types which indicates potential for sexual recombination and rapid evolution of new pathogen strains, although we found no evidence of recent recombination. The associated eukaryotic and prokaryotic leaf microbiome exhibited variation in taxonomic diversity between years and locations and is likely influenced by local weather conditions. A network analysis revealed distinct clusters of eukaryotic and prokaryotic taxa that correlates with the frequency of E. turcicum sequencing reads, suggesting causal interactions. Notably, the yeast genus Metschnikowia exhibited a strongly negative association with E. turcicum, supporting its known potential as biological control agent against fungal pathogens. Our findings show that metagenomic sequencing is a useful tool for analysing the role of environmental factors and potential pathogen–microbiome interactions in shaping pathogen dynamics and evolution, suggesting their potential for effective pathogen management strategies.