Population Connectivity Research Articles

A Larval “Recruitment Kernel” to Predict Hatching Locations and Quantify Recruitment Patterns

AbstractLarval recruitment, a critical component of population connectivity, has been under investigated compared to larval dispersal. We developed a backward‐in‐time Lagrangian particle tracking model to predict larval hatching locations and proposed a larval recruitment kernel, to quantify recruitment patterns. Combining field data and a hydrodynamic model, our backtracking model predicted Lake Whitefish (Coregonus clupeaformis) hatching locations in Lake Erie. We found a strong linear correlation (r = 0.95–0.98) between travel distance (i.e., distance along a trajectory) and pelagic larval duration (PLD), and a moderate correlation (r = 0.66–0.68) between linear distance (i.e., displacement) and PLD. This questions the wide use of PLD as a proxy for dispersal distance. We defined the recruitment kernel using the probability density function of the linear recruitment distance. Characteristics of the recruitment kernel, such as theoretical self‐recruitment, median‐recruitment distance, long‐distance recruitment, and openness convey significant information about population connectivity that are distinct from those derived using the well‐known dispersal kernel (e.g., theoretical local retention).

Genetic predictors of population resilience: A case study of native Brook Trout in headwater streams

AbstractObjectivePopulations of eastern Brook Trout Salvelinus fontinalis face threats from several sources, such as habitat fragmentation, climate change, and competition with introduced salmonids. As a native species, understanding how these populations will respond to disturbances is paramount to their management and effective conservation. A population's ability to respond to disturbance, its resilience, is influenced by several factors. One such group of factors is population genetics.MethodsWe calculated population resilience metrics based on transient dynamics using population projection matrix models. Long‐term demographic data from 23 headwater stream Brook Trout populations were used to parameterize models. Genetic data were collected, and genetic indices were calculated. Partial redundancy analysis was then used to evaluate relationships between resilience metrics and genetic indices.ResultInbreeding coefficient, rarefied allelic richness, pairwise genetic differentiation (FST), and effective population size were all found to be important variables in predicting resilience.ConclusionOur results suggest that genetic isolation may increase the demographic resilience in Brook Trout through faster generation times and higher juvenile survival, but this likely comes at the cost of increased extinction risk and truncated size structures. Genetic indices can provide insight into gene flow between populations, thus the relationship between population connectivity and resilience. Given the importance of connectivity to population resilience, restoring and maintaining movement corridors could affect resilience in headwater Brook Trout populations.

Landscape features influencing gene flow and connectivity of an endangered passerine.

Dispersal of individuals and gene flow are crucial aspects to maintain genetic diversity and viability of populations, especially in the case of threatened species. Landscape composition and structure may facilitate or limit individual movement within and among populations. We used a landscape genetics approach to assess the connectivity patterns of the threatened Dupont's lark (Chersophilus duponti subsp. duponti), considering their genetic patterns and the landscape features associated with its gene flow in Spain. We analysed the genetic relatedness based on 11 species-specific polymorphic microsatellites on 416 Dupont's lark individuals sampled across peninsular Spain between 2017 and 2019, covering most of the European distribution of the species. To assess the relationship between the landscape composition and the species gene flow, we estimated genetic distance at the individual level (Dps). Next, we built a set of environmental surfaces from two time periods (years 1990 and 2018), based on factors such as land use and topography, influencing individuals' movement. We then obtained resistance surfaces from an optimization process on landscape variables. Landscape genetics analyses were done for single and composite surface models for each year separately. Our findings from both time points show that scatter or mosaic-structured vegetation composed by low agricultural and tree cover and high presence of sclerophyllous shrubs favoured Dupont's lark dispersal, while dense and continuous tree cover, as well as areas of intensive agriculture, were limiting factors. Our results suggest the importance of steppe habitat patches for the species' establishment and dispersal. In addition, our results provide key information to develop conservation measures, including conserving and restoring steppe habitats as scattered and/or mosaic-structured vegetation that could warrant the connectivity and persistence of Dupont's lark populations.

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia.

Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan. We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM). Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM. This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

Otolith biogeochemistry reveals possible impacts of extreme climate events on population connectivity of a highly migratory fish, Japanese Spanish mackerel Scomberomorus niphonius

Otolith biogeochemistry reveals possible impacts of extreme climate events on population connectivity of a highly migratory fish, Japanese Spanish mackerel Scomberomorus niphonius

Genetic history of esophageal cancer group in southwestern China revealed by Y‐chromosome STRs and genomic evolutionary connection analysis

Genetic and environmental factors play crucial roles in the development of esophageal cancer (EC) and contribute uniquely or cooperatively to human cancer susceptibility. Sichuan is located in the interior of southwestern China, and the northern part of Sichuan is one of the regions with a high occurrence of EC. However, the factors influencing the high incidence rate of EC in the Sichuan Han Chinese population and its corresponding genetic background and origins are still poorly understood. Here, we utilized genome-wide single nucleotide polymorphisms (SNPs) and Y-chromosome short tandem repeats (Y-STRs) to characterize the genetic structure, connection, and origin of cancer groups and general populations. We generated Y-STR-based haplotype data from 214 Sichuan individuals, including the Han Chinese EC population and a control group of Han Chinese individuals. Our results, obtained from Y-STR-based population statistical methods (analysis of molecular variance (AMOVA), principal component analysis (PCA), and phylogenetic analysis), demonstrated that there was a genetic substructure difference between the EC population in the high-incidence area of northern Sichuan Province and the control population. Additionally, there was a strong genetic relationship between the EC population in the northern Sichuan high-incidence area and those at high risk in both the Fujian and Chaoshan areas. In addition, we obtained high-density SNP data from saliva samples of 60 healthy Han Chinese individuals from three high-prevalence areas of EC in China: Sichuan Nanchong, Fujian Quanzhou, and Henan Xinxiang. As inferred from the allele frequency of SNPs and sharing patterns of haplotype segments, the evolutionary history and admixture events suggested that the Han population from Nanchong in northern Sichuan Province shared a close genetic relationship with the Han populations from Xinxiang in Henan Province and Quanzhou in Fujian Province, both of which are regions with a high prevalence of EC. Our study illuminated the genetic profile and connection of the Northern Sichuan Han population and enriched the genomic resources and features of the Han Chinese populations in China, especially for the Y-STR genetic data of the Han Chinese EC population. Populations living in different regions with high incidences of EC may share similar genetic backgrounds, which offers new insights for further exploring the genetic mechanisms underlying EC.

Combining acoustic telemetry with archival tagging to investigate the spatial dynamic of the understudied pollack, Pollachius pollachius.

Combining fish tracking methods is a promising way of leveraging the strengths of each approach while mitigating their individual weaknesses. Acoustic telemetry provides presence information as the fish move within receiver range, eliminating the need for tag recovery. Archival tags, on the other hand, record environmental variables on tag retrieval, enabling continuous path reconstruction of a fish beyond coastal regions. This study capitalizes on the combination of both methods for geolocating pollack, Pollachius pollachius, an understudied species of the northeast Atlantic, where declining stocks are raising concern. Essential knowledge of population structure and connectivity between essential habitats is critically lacking and could help inform stock assessment and management. The aims of the study were (1) to evaluate the feasibility of double-tagging pollack, known for being challenging to tag, and (2) to track seasonal movements across the Channel to gain first insights into pollack spatial ecology. In 2022, an extensive network of acoustic receivers was been deployed in the Channel along the French, English, and Belgian coasts as part of the Fish Intel project. We tagged 83 pollack with acoustic transmitters, among which 48 were double-tagged with data storage tags. Post-tagging survival assessment, conducted on a subset of 35 individuals, revealed a successful procedure with a 97% short-term survival rate. By October 2023, the acoustic telemetry network detected 30 out of 83 pollack at least once, with no large-scale movements observed across the Channel. Presence in the network fluctuates seasonally, peaking in summer, particularly among immature fish. Integrating acoustic detections with temperature and depth time series in a geolocation model enabled trajectory reconstruction of 10 recaptured pollack, seven of which were detected by the network. This combined tracking approach revealed coastal movements along the coast of Brittany in France, highlighting the ecological significance of the Iroise Sea for pollack throughout the year, particularly in summer. The geolocation model also suggested movements towards the entrance of the western Channel. This study highlights the complementarity of acoustic telemetry and archival tagging in reconstructing fish movements in their natural environment. As data accumulate, these innovative tracking methods promise to continually unveil new insights into the spatial ecology of the understudied pollack, which is essential for the species' management.

Combining population genomics and biophysical modelling to assess connectivity patterns in an Antarctic fish.

Connectivity is a fundamental process of population dynamics in marine ecosystems. In the last decade, with the emergence of new methods, combining different approaches to understand the patterns of connectivity among populations and their regulation has become increasingly feasible. The Western Antarctic Peninsula (WAP) is characterized by complex oceanographic dynamics, where local conditions could act as barriers to population connectivity. Here, the notothenioid fish Harpagifer antarcticus, a demersal species with a complex life cycle (adults with poor swim capabilities and pelagic larvae), was used to assess connectivity along the WAP by combining biophysical modelling and population genomics methods. Both approaches showed congruent patterns. Areas of larvae retention and low potential connectivity, observed in the biophysical model output, coincide with four genetic groups within the WAP: (1) South Shetland Islands, (2) Bransfield Strait, (3) the central and (4) the southern area of WAP (Marguerite Bay). These genetic groups exhibited limited gene flow between them, consistent with local oceanographic conditions, which would represent barriers to larval dispersal. The joint effect of geographic distance and larval dispersal by ocean currents had a greater influence on the observed population structure than each variable evaluated separately. The combined effect of geographic distance and a complex oceanographic dynamic would be generating limited levels of population connectivity in the fish H. antarcticus along the WAP. Based on this, population connectivity estimations and priority areas for conservation were discussed, considering the marine protected area proposed for this threatened region of the Southern Ocean.

Appalachians to the Andes: Potential population connectivity and loop migration of Canada Warblers (Cardellina canadensis) revealed by light-level geolocators

Appalachians to the Andes: Potential population connectivity and loop migration of Canada Warblers (Cardellina canadensis) revealed by light-level geolocators

Broad-scale genetic assessment of Southern Ground-Hornbills (Bucorvus leadbeateri) to inform population management

The Southern Ground-hornbill (SGH) (Bucorvus leadbeateri) is considered an umbrella species for biodiversity conservation in savannah biomes since they require large territories and significant protection measures that help to conserve a wide range of biodiversity with similar savanna and grassland requirements. Declines of the species are attributed to low reproductive rates coupled with multiple anthropogenic threats, including secondary poisoning, and persecution. Little is known about connectivity and population structure of SGH populations in Africa, south of the equator. Knowledge of population differentiation is needed to ensure that targeted conservation management plans can be implemented to slow population declines and ensure survival of the species. To inform a long-term conservation strategy, we investigated the broad-scale population structure of Southern Ground-hornbill across their sub-equatorial range. Our study based on 16 microsatellite loci identified moderate variation (average of 5.889 alleles per locus and a mean observed heterozygosity of 0.546) similar to other long-lived avian species. In contrast, mitochondrial DNA sequences analysis identified low diversity (Hd = 0.3313, π = 0.0015). A Bayesian assignment approach, principal component analysis, analysis of molecular variance and phylogenetic analysis identified weak to moderate population structuring across long distances and mitochondrial data showed a shallow phylogeny. Restriction to long-distance dispersal was detected that could not be attributed to isolation by distance, suggesting that other factors, such as their dispersal biology, are shaping the observed genetic differentiation. Although our study does not support the designation of populations as independent conservation units, we advocate that population management should continue to follow the Precautionary Principle (mixing founders from the same range state, rather than allowing mixing of founders from the extremes of the range) until there is scientific certainty. Following further research, if no independent conservation units are detected, then the global captive population can contribute to reintroductions across the range. In the wild, populations at the edge of the species range may need additional management strategies and gene flow should be promoted between neighbouring populations.

A reference quality genome assembly for the jewel scarab Chrysina gloriosa.

The jewel scarab Chrysina gloriosa is one of the most charismatic beetles in the United States and is found from the mountains of West Texas to the Southeastern Arizona sky islands. This species is highly sought by professional and amateur collectors worldwide due to its gleaming metallic coloration. However, the impact of the large-scale collection of this beetle on its populations is unknown, and there is a limited amount of genetic information available to make informed decisions about its conservation. As a first step, we present the genome of C. gloriosa, which we reconstructed using a single female specimen sampled from our ongoing effort to document population connectivity and the demographic history of this beetle. Using a combination of long-read sequencing and Omni-C data, we reconstructed the C. gloriosa genome at a near-chromosome level. Our genome assembly consisted of 454 scaffolds spanning 642 MB, with the 10 largest scaffolds capturing 98% of the genome. The scaffold N50 was 72 MB, and the BUSCO score was 95.5%. This genome assembly will be an essential tool to accelerate understanding C. gloriosa biology and help make informed decisions for the conservation of Chrysina and other species with similar distributions in this region. This genome assembly will further serve as a community resource for comparative genomic analysis.

Phenotypic and genetic divergence in a cold‐adapted grasshopper may lead to lineage‐specific responses to rapid climate change

AbstractAimSpecies responses to global warming will depend on intraspecific diversity, yet studies of factors governing biogeographic patterns of variability are scarce. Here, we investigate the evolutionary processes underlying genetic and phenotypic diversity in the flightless and cold‐adapted grasshopper Sigaus piliferus, and project its suitable space in time.LocationTe Ika‐a‐Māui Aotearoa—North Island of New Zealand.MethodsWe used mitochondrial sequences to investigate population connectivity and demographic trends using phylogeographic tools and neutrality statistics. Metric data were used to document phenotypic variation using naïve clustering. We used niche metrics to assess intraspecific niche variation, and niche modelling to investigate suitability under past and future scenarios. Multiple matrix regressions with randomization explored the processes contributing to phenotypic differentiation among grasshopper populations.ResultsNiche models and demographic analyses suggest suitable space for this grasshopper was more restricted during glacial than interglacial stages. Genealogical relationships among ND2 haplotypes revealed a deep north–south split partly concordant with phenotypic and niche variation, suggesting two ecotypes that have mixed during recolonisation of the central volcanic region. Multiple matrix regressions with randomization indicate a link between climate and phenotypic differentiation inferred from leg and pronotum dimensions but not pronotum shape. Niche projections predict severe habitat reduction due to climate warming.Main conclusionsThe current distribution and intraspecific diversity of S. piliferus reflect complex biogeographical scenarios consistent with Quaternary climates and volcanism. Phenotypic divergence appears adaptive. Current levels of genetic and phenotypic variation suggest adaptive potential, yet the pace of anthropogenic warming over the next 50 years could result in small populations that may collapse before adapting. Differences in niche features between diverging intraspecific lineages suggest distinct responses to climate change, and this has implications for prioritising conservation actions and management strategies.

Spatio-temporal connectivity and dispersal seasonal patterns in the Adriatic Sea using a retention clock approach

Hydrodynamic features play a key role in determining the dispersal and connectivity of fish populations, especially in highly energetic areas determined by currents, river flow, and meteorologically induced fluctuations. Understanding how species interact with these physical processes is essential for managing vulnerable populations and identifying areas that require effective conservation efforts. This study examines the hydrodynamics that regulate connectivity in the Adriatic Sea, a shallow and semi-enclosed basin that is widely recognized as one of the most important areas in the Mediterranean Sea for protection. A high-resolution hydrodynamic model coupled with a lagrangian tracking module serves as the numerical tool. Lagrangian particles, representing eggs and larvae with typical biological characteristics of generic marine organisms inhabiting the region, are released throughout the basin at different times during a test year to identify the most likely pathways of individual dispersal. The temporal component of connectivity is highlighted using a previously developed retention clock matrix over different larval durations. Seasonality is a critical factor in dispersal, with greater variability and reduced efficiency in winter compared to summer. The potential implications of the results for improved assessment and management of high value marine species in the basin are discussed.

The high cost of movement in an arid working landscape for an endangered amphibian.

Connectivity is essential for the maintenance of genetic diversity and stability of wildlife populations. Drought and changing precipitation regimes have caused natural aquatic amphibian breeding habitats to disappear or become isolated and have led to the replacement of natural surface water with artificial livestock water tanks. Terrestrial movement is the only means of responding to aquatic threats in arid landscapes and to allow population connectivity. Aridity may present an impenetrable barrier in hydrologically fragmented environments. We used a facultatively paedomorphic and federally endangered salamander to assess the challenges of movement across arid working lands. Sonoran tiger salamanders (Ambystoma mavortium stebbinsi) are endemic to the San Rafael Valley of southeastern Arizona, United States of America, where they depend on livestock water tanks as breeding habitat. The ecology of this species' metamorphs outside of stock tanks is virtually unknown. To assess survival on the landscape during terrestrial movement we used radio-transmitters to track 78 adult metamorphosed salamanders over 2 years. Sonoran tiger salamanders moved up to 1 km from the tank edge, and average distances moved of over 400 m were higher than most Ambystoma species. However, during the study period, none reached neighboring stock tanks. We found high mortality due to predation and desiccation. Individuals that dispersed to terrestrial habitat in summer survived longer than individuals that dispersed in spring. High mortality suggests terrestrial movement is exceptionally risky and may contribute to isolated subpopulations and elevated levels of inbreeding. Conservation actions that improve and maintain artificial aquatic habitats as well as increase connectivity may improve long-term management for pond-breeding amphibians in arid regions.

Prioritizing conservation areas to mitigate connectivity loss and local extinction risk of a small carnivore (Leopardus guttulus) in South America

Effective conservation management depends on the maintenance of key areas that allow population connectivity across the landscape. However, the lack of knowledge of how habitat conversion affects species movement hinders the identification of these areas. Here, we analyzed the impact of habitat fragmentation on landscape connectivity for Leopardus guttulus, a small Neotropical felid threatened by the high habitat fragmentation across the Atlantic Forest, and mapped and ranked the most important core areas and corridors for conservation actions. We also estimated genetic diversity indices and predicted the viability of the current core areas in the future. Our analyses suggest that L. guttulus populations are fragmented, and connectivity links between populations are few and weak. We predict that due to their size, estimated density and low connectivity, some current core areas may not maintain viable populations in the long-term. Also, ongoing land-use changes may further isolate remaining populations, leading to progressive reductions in the populations they support. In this study, we spatially prioritize the most critical areas for L. guttulus conservation and highlighted the urge that exists in the adoption of management measures for its conservation.

Topographical and cell type-specific connectivity of rostral and caudal forelimb corticospinal neuron populations

Corticospinal neurons (CSNs) synapse directly on spinal neurons, a diverse assortment of cells with unique structural and functional properties necessary for body movements. CSNs modulating forelimb behavior fractionate into caudal forelimb area (CFA) and rostral forelimb area (RFA) motor cortical populations. Despite their prominence, the full diversity of spinal neurons targeted by CFA and RFA CSNs is uncharted. Here, we use anatomical and RNA sequencing methods to show that CSNs synapse onto a remarkably selective group of spinal cell types, favoring inhibitory populations that regulate motoneuron activity and gate sensory feedback. CFA and RFA CSNs target similar spinal neuron types, with notable exceptions that suggest that these populations differ in how they influence behavior. Finally, axon collaterals of CFA and RFA CSNs target similar brain regions yet receive highly divergent inputs. These results detail the rules of CSN connectivity throughout the brain and spinal cord for two regions critical for forelimb behavior.

The widespread keeping of wild pets in the Neotropics: An overlooked risk for human, livestock and wildlife health

Abstract Zoonoses constitute a major risk to human health. Comprehensive assessments on the potential emergence of novel disease outbreaks are essential to ensure the effectiveness of sanitary controls and to establish mitigating actions. Through a continental‐scale survey of rural human settlements conducted over 13 years in 15 Neotropical countries, we document the vast extent of poaching to meet the local demand for pets, resulting in thousands of families living with ca. 275 species of wild animals without any sanitary controls. Parrots account for ca. 80% of wild pets, dying mostly from diseases at an average age of 1 year. This culturally rooted tradition, which dates back to pre‐Columbian times, may lead to health risks by bringing wild animals prone to carrying parasites and pathogens into close contact with humans and their exotic pets and livestock. Although animal pathogens and parasites have been transmitted to humans for centuries, the current trend of human population growth and connectivity can increase the risk of zoonotic outbreaks spreading at an unprecedented pace. Similarly, disease transmission from humans and poultry to wild animals is also expected to be facilitated via wild pets, leading to conservation problems. Several studies have highlighted the risk posed by wildlife city markets for cross‐species disease transmission, ignoring the risk associated with widespread pet ownership of wild animals poached locally in rural areas. Given its geographic and social dimensions, a holistic approach is required to reduce this illegal activity as well as to strengthen health surveillance of seized individuals and people in close contact with poached pets, which would benefit both people and wildlife. Read the free Plain Language Summary for this article on the Journal blog.

Pollinator-mediated connectivity in fragmented urban green spaces—tracking pollen grain movements in the city center

Cities are expanding rapidly and emerging as complex ecosystems. This expansion of urbanized areas has resulted in the alteration and fragmentation of several natural habitats. Recent studies have shown that cities support a significant level of biodiversity. This unexpected role of cities holds promise for advancing knowledge of urban ecosystems and promoting their conservation. In particular, the connectivity of plant populations through pollen transfer is crucial for the long-term persistence of insect-pollinated plant species. In this study, we tracked pollen movement areas in four isolated patches of urban greenery in urbanized area using quantum dots. We studied Fritillaria imperialis (spring) and Hemerocallis sp. (mid-summer). Our research revealed frequent pollen transfer between small, isolated flowering patches, even when these locations were not connected by green corridors. Common elements found in urban ecosystems, such as streets, trails, and pavements, do not stop the dispersion of pollen grains. Moreover, the migration pathways of pollen grains vary for each species according to different factors. For F. imperialis, we identified the proportion of green areas in proximity to the study location as a key factor in shaping pollen transfer. For Hemerocallis, we discovered that the proportion of green areas, the distance between study sites, and the frequency of pollinator visits are influential factors. Our study demonstrates that small isolated plant populations exchange pollen due to pollinator movement, indicating that these small populations may serve as stepping stones for pollinators among larger populations.

Genomic analysis of population history for Hawaiian monk seals

The Hawaiian monk seal Neomonachus schauinslandi, one of the world’s most endangered pinnipeds, has faced decades of declines and has been the focus of intensive conservation efforts. A myriad of conservation threats has led to range-wide population declines, but population trends among islands can vary widely in response to heterogeneous threats. Populations in the Northwestern Hawaiian Islands have been declining, whereas Main Hawaiian Islands numbers are expanding. Molecular data can provide information to disentangle population structure and dynamics; however, previous studies have yielded insufficient resolution in such a genetically depauperate species. Advances in genomic technology and affordability offer a novel opportunity to revisit questions about Hawaiian monk seal trends with high-resolution markers that provide better discrimination ability in low-diversity species. Here, we investigated region- and island-level population structuring and connectivity. We used BestRAD sequencing on 169 seals from 14 islands that span the archipelago to estimate genetic diversity, genetic differentiation, population structure, and migration rates. We did not find robust evidence for island-level population structure. For the first time, our data set provided resolution to differentiate regional populations with low but significant genetic differentiation. Further, DAPC illustrated population structure with evidence for connectivity, which mirrored our migration rate estimates. Future conservation decisions will need to consider the balance of maintaining connectivity between regions while not homogenizing and losing valuable, yet rare, regional unique variation.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.