Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean

Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean View ORCID Profile Ekaterina Bocharova , View ORCID Profile Aleksandr Volkov , View ORCID Profile Solenn Stoeckel doi: https://doi.org/10.1101/2025.04.24.650399 Ekaterina Bocharova 1 Faculty of Life Sciences, University of Vienna , Vienna, Austria 2 Vavilov Institute of General Genetics of Russian Academy of Sciences , Moscow, Russia Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Ekaterina Bocharova For correspondence: bocharova.ekaterina{at}gmail.com solenn.stoeckel{at}inrae.fr Aleksandr Volkov 3 Branch for the freshwater fisheries of VNIRO (VNIIPRKH) , Moscow region, Russia Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Aleksandr Volkov Solenn Stoeckel 4 DECOD, L’Institut Agro , IFREMER, INRAE, 35000, Rennes, France 5 IGEPP, INRAE, Institut Agro, Université de Rennes , Le Rheu, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Solenn Stoeckel Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract Reproductive mode is a key factor shaping genetic diversity, evolutionary potential, and the processes of dispersal and colonization. Clonality is particularly common in harsh environments and at the margins of species ranges, where it supports persistence, enables rapid growth, and promotes the maintenance of locally adapted genotypes. In the rapidly changing Arctic, increasing ecological connectivity is eroding historical barriers for sessile species. Evaluating genetic diversity in this context, before global change further alters Arctic ocean, is essential for understanding evolutionary dynamics during range expansion and for informing conservation strategies. Aulactinia stella is a circumpolar sea anemone with physiological characteristics in laboratory conditions suggesting a potential for clonal reproduction. In this study, we investigated its reproductive modes in natural populations across the Arctic ocean, from the northern Pacific to the Atlantic, and examined how genetic diversity is structured between adults and juveniles at five sampled sites. Across all study sites, we observed only females or individuals lacking gonads, with the exception of Kamchatka, where males were also present. Genetic indices and changes in genotype frequencies between adults and juveniles confirmed that this species reproduces partially by parthenogenesis. Populations on the Atlantic side were highly clonal with clonal rates ( c ) estimated at 80-99%, whereas populations on the Pacific side reproduced more sexually (c around 50%). Allelic diversity was twice as high in Kamchatka and Kuril populations, suggesting North Pacific coasts being the main last glacial refugia of A. stella . We found a stepping-stone pattern of genetic structure from Kamchatka to Atlantic populations, consistent with contemporary ocean currents and melted summer sea ice. Only a subset of the juvenile genetic diversity, mostly of local origin, was found in the established adults, while juveniles exhibited lower levels of genetic differentiation across the Arctic Ocean. Our findings underscore the need for further ecological and behavioral investigations to elucidate the mechanisms allowing the current possibilities of dispersal of this species across the Arctic Ocean. Download figure Open in new tab Graphical Abstract (photo: Ekaterina Bocharova) Introduction The way individuals reproduce is the most influential trait in governing genetic diversity and potentially evolutionary response to environmental changes ( Duminil et al., 2007 ; Pierre et al., 2022 ). Reproductive mode governs the persistence of DNA variations through generations. It has therefore major impacts on the level of genome-wide genetic diversity and its structuration within and among populations ( Ellegren & Galtier 2016 ). Reproductive modes also profoundly affect other biological traits that play important roles for ecological processes, including population growth, dispersal, colonization of new areas and the genetic basis for adaptation ( Yu et al., 2016 ; He et al., 2024 ). Clonality is a dominant reproductive mode found in populations colonizing new area, especially when enduring harsh conditions at the edge of their living range ( Cadotte et al., 2006 ; Burke & Bonduriansky, 2018 ; Roiloa, 2019 ; Franklin et al., 2021 ) and in recently invaded, resource-poor environments to which individuals are poorly adapted (geographic parthenogenesis; Peck et al., 1998 ). All groups and kingdoms of Eukaryota domain, including most of the species that are essential for ecosystem functions, have species that reproduce using partial clonality ( Avise, 2015 ; Schön et al., 2009 ; Tibayrenc et al., 2015 ; Orive & Krueger-Hadfield, 2021 ). In this reproductive mode, parents can generate offspring that are genetically identical to themselves at all loci on their genomes, with the exception of rare recombination events and occasional somatic mutations ( De Meeûs et al., 2007 ; Stoeckel et al., 2021 a). In such species, the proportion of descendants produced by clonality, i.e. the rate of clonality, is a key biological factor to understand their genetic and trait evolution ( Yu et al., 2016 ; Becheler et al., 2017 ). Clonality often increases in populations expanding at range margins, as it allows pioneers to establish and persist in the absence of mates ( Baker & Stebbins, 1965 ; Liu et al., 2006 ; Wang et al., 2019 ). Moreover, clonality enables locally successful genotypes to persist across generations without recombining with less adapted ones, thereby accelerating population growth ( Fouqueau & Roze, 2021 ; Pierre et al., 2022 ). It also decreases the loss of allele diversity and of heterozygosity by genetic drift, especially in small founding populations (Navascues et al., 2010; Stoeckel & Masson, 2014 ; Reichel et al., 2016 ; Stoeckel et al., 2021 b). However, the absence of recombination in clonal reproduction constrains selection by preventing the combination of adaptive alleles and hindering the removal of deleterious variants ( de Visser & Elena, 2007 ; López-Villavicencio et al., 2013; Pierre et al. 2022 ). The Arctic Ocean has been previously considered as an insurmountable barrier for sessile and semi-sessile marine species dispersal due to limited access, harsh environmental conditions and inadequate food resources ( Chan et al., 2019 ). Climate change is opening new routes across the Arctic Ocean, linking northern Pacific and Atlantic biomes ( Vermeij & Roopnarine, 2008 ). The recent acceleration of summer sea-ice loss, now surpassing levels recorded since the Holocene Climatic Optimum (∼8,000 years before present, estimated at +0.3 °C relative warming above pre-industrial levels; compared to +1.34 °C in the Arctic in 2024; Laakkonen et al., 2021 ; E.U. Copernicus Global Climate Highlights, 2024), is likely to facilitate renewed secondary contact between species and populations from the northern Atlantic and Pacific biomes ( Brandt et al., 2023 ). Understanding genetic diversity and population structure in such species provide useful insights into the ecological and evolutionary processes that occur during range expansion under real and complex conditions ( Sax et al., 2007 ; Moran & Alexender, 2014; Bridle & Hoffmann, 2022 ). It also provides a frame of reference to managers and policymakers to identify priorities and mechanisms for mitigating adverse effects of global changes on intraspecific diversity and for controlling future invasive populations ( Chan et al., 2019 ). Sea anemones are semi-sessile organisms reproducing using partial clonality ( Bocharova & Kozevich, 2011 ; Bocharova, 2016 ). Aulactinia stella (Cnidaria, Anthozoa) is an intertidal and sublittoral species belonging to Actiniidae family. It was first observed and described on the eastern coast of the USA in the 19th century ( Verrill, 1864 ), probably in the southernmost part of its distribution area. Its geographic range is rather circumpolar, along the coasts of the Arctic Ocean and of the northern parts of Atlantic and Pacific Oceans (Grebelnyi, 1980; Sanamyan & Sanamyan, 1998 , 2006 ). First studies of A. stella populations in the Arctic Ocean, on the Murman Coast, were carried on about fifty years ago ( Loseva, 1974 ). Adults of A. stella mainly live fixed to rocky ground, even if it can occasionally be found buried in sandy or muddy soils, hidden in crevices or fixed to debris ( Bocharova, 2012 ). It feed predating smaller mollusks, crustaceans, fish larvae and other animals ( Bocharova, 2012 ). The population density of A. stella ranges from 1 to 18 individuals per square meter in the littoral zone of the White Sea, and from 1 to 14 individuals per square meter along the Barents Sea coasts ( Bocharova, 2012 ). In laboratory conditions, adults can survive for several decades, but their natural lifespan in the wild still remains unknown. The life cycle of A. stella includes a swimming larval stage, which can be brooded by the parental anemone until it develops into a polyp. The polyp stage is benthic, but it retains the ability to be transported passively by water currents until it settles as adult on a solid substrate. In response to unfavorable abiotic or biotic conditions, many sea anemone species are capable of detaching from the substrate and temporarily drifting with tidal and coastal currents, a behavior that is particularly common when winter conditions become extreme (Riemann-Zürneck, 2008). Some circum-Antarctic and Arctic sea anemones present a distinct post-metamorphic juvenile stage, following the larval phase, during which they drift within pelagic currents (Riemann-Zürneck, 2008). To date, no comparable pelagic stage has been documented or investigated for A. stella despite the presence of a free-living polyp stage within the benthic fauna with similar potential ( Bocharova, 2012 ). Reproduction occurs annually from May to October, with phenological variation depending on local environmental conditions, but a peak period between June and August ( Bocharova, 2012 ). Along the Kamchatka coasts during this time, males release gametes passively into the water column, that are later ingested by females. Females brood the developing larvae internally ( Bocharova, 2012 ). As for many circumpolar marine species, direct observations of reproductive and dispersal processes still remain unrealistic. Consequently, genetic analyses remain the most effective approach for inferring key biological and ecological traits of this species in natural conditions. Developmental and histological evidence suggests that this species may reproduce through parthenogenesis ( Bocharova, 2012 ). Preliminary analyses based on a limited number of mitochondrial and ribosomal sequences hypothesized that these populations are likely highly clonal and exhibit strong genetic structuring across the Arctic Ocean despite the region’s increasing ice melt over recent decades ( Bocharova & Mugue, 2012 ; Bocharova, 2015 ). While these sequences were suitable for phylogenetic analyses, they lacked the resolution necessary required to confirm clonality and accurately characterize the genetic structure of A. stella populations. The recent development of microsatellite primer pairs now enables investigating the biology of this circumpolar sea anemone more precisely, using approaches based on genetic diversity and its spatial patterns ( Bocharova et al., 2018 ). In this study, we aimed to confirm the occurrence of clonal reproduction in natural populations of A. stella , quantify its extent and variability across the species’ range, and use high-resolution, polymorphic markers to accurately assess population structure throughout the Arctic Ocean. To investigate clonal reproduction in this species, we analyzed genotypic and genetic indices across five geographically distinct sites spanning the Arctic Ocean, including locations in the northern Atlantic and northern Pacific. We leveraged previous sampling campaigns (2012–2017), which included fewer sites and greater temporal gaps in the Pacific. While this may influence some interpretations, our results still provide a robust overview of the species’ reproductive variation and current biogeography, especially given its laboratory lifespan of at least more than one decade. A total of 314 individuals, representing two age cohorts, juveniles and adults, were sampled to estimate the contribution of clonal reproduction and its variation across the species’ range. By examining temporal changes in genotypic frequencies within each site, we aimed to infer the prevalence of clonal reproduction. We hypothesized that the limited energetic and physiological cost of clonal reproduction may facilitate the establishment and persistence of populations under harsh polar coastal conditions ( Peck et al., 1998 ; Klimeš, 2008 ; Wang et al., 2018 ). Conversely, we also considered that the genetic recombination associated with sexual reproduction may promote the emergence of genotypes better adapted to local environments ( Pierre et al., 2022 ). Our findings on reproductive modes are interpreted in the context of the specific ecological challenges faced by species expanding into previously unfavorable habitats. In addition, we used the genotyped populations to assess genetic differentiation and identity in state between individuals and sampled sites. Sampling across the full distribution range allowed us to pinpoint regions with higher genetic diversity, which may correspond to historical refugia. By comparing genetic differentiation between juveniles and adults, we aimed to identify present-day or historical corresponding dispersal routes. Although A. stella was first described along the North Atlantic coasts, suggesting a possible Atlantic origin, previous limited sequence data instead indicated higher genetic diversity in North Pacific populations. If confirmed using more polymorphic and resolutive markers that are microsatellites, this may suggest a more recent expansion into the North Atlantic via the Arctic Ocean. Considering the free-living polyp stage which can disperse passively via marine currents, as well as the potential for semi-sessile adults to be transported on drifting debris or through human-mediated movement of materials, we expect a limited but ongoing level of gene flow facilitated by favorable summer sea conditions. Under this scenario, we would expect to observe shared probabilities of identity in state among populations and a stepping-stone pattern of genetic differentiation across the Arctic, potentially with preferential dispersal routes connecting sampled populations. Alternatively, considering populations evolving isolated due to current Arctic conditions for ∼8,000 years (Holocene Thermal Maximum) or even ∼125,000 years (last interglacial Eemian period), we would expect very distinct and strongly structured genetic diversity with even fixed and private alleles between populations ( Slatkin & Takahata, 1985 ), out of artefactual homoplasy. Ultimately, our findings will provide a baseline framework for understanding the genetic structure of this species across its Arctic distribution, offering critical insights at a time of profound environmental change in this region ( Laakkonen et al., 2021 ; Brandt et al., 2023 ). Materials and methods Sample collection In September 2012, we sampled at the low tide in the littoral zone 23 Aulactinia stella in the White Sea Biological Station of Moscow State University (WSBS), Kandalaksha Bay of the White Sea and 30 individuals at Liinkhamari, in the Varanger Fjord of the Barents Sea (see Bocharova 2015 ). We also sampled using SCUBA diving at 5-18 m depth in the sublittoral zone 29 individuals at the Cape Broyle in the Admirals Cove area, Newfoundland, Canada, in June 2017. Along Kamchatka coast of the Pacific Ocean, we collected 5 individuals in the Avachinsky Bay in June-July 2014 and 9 individuals in the Viluchinsky Bay in June-July 2016 at 5-21 m depth. We used these 14 Kamchatka adult individuals for mitochondrial DNA sequencing, but only 6 of them brooded juveniles. Therefore, for microsatellite analysis, we additionally analysed 3 females with their brooded offspring collected in the Avachinsky Bay in July 2012 (2 samples) and June 2013 (1 sample). At low tide, we sampled 20 individuals in the intertidal zone of the Matua Island of the Kuril Islands in August 2016 (see map of the sampled sites on Figure S1 ). We checked that the phenotypic and biometric appearance of each individual belonged to A. stella species using descriptions in Loseva (1974) , Fautin Dunn et al. (1980) , Sanamyan and Sanamyan (2008), Sanamyan et al. (2019 , 2020 ) and WoRMS ( Rodríguez et al., 2025 ). In addition, we used nematocyst identification from tentacles, pharynx, body wall and mesentery ( Manuel, 1981 ), and mitochondrial DNA sequence analysis to confirm the species. To analyze genetic diversity in two successive generations, we purposely genotyped large females and the juveniles they brooded in their gastrovascular cavity. The quantity of juveniles in such females varied from 1 to 39 per female and we sequenced them all. DNA extraction, mtDNA sequencing and genotyping by microsatellite loci We extracted the total DNA of each sampled adult and juvenile polyp. For each adult, we used a part of the adult pedal disk measuring approximately 3x3x3 mm; And for juveniles, we extracted DNA from the whole juvenile polyps which measured from 1 to 5 mm. For each sample, we desiccated and grinded the sampled tissues. We then used the Wizard SV Genomic DNA Purification System (Promega, USA) following the manufacturer’s protocol. Extracted DNA was preserved in 96% ethanol before genotyping. For mitochondrial gene (12S rRNA, 16S rRNA, cytochrome oxidase III) sequencing, we only used adult subsamples: 23 individuals from the White Sea, 30 from the Barents Sea, 29 from Newfoundland, 14 from Kamchatka, and 20 from the Kuril Islands. We applied the primer sets and methodology thoroughly described for 16S rRNA and COIII ( Geller & Walton, 2001 ) and 12S rRNA mitochondrial fragments ( Bocharova 2015 ). Sequencing was performed by 3500 Genetic Analyzer (Applied Biosystems, USA) using BigDye Terminator v3.1 Cycle Sequencing Kit. The sequences were analyzed by Geneious Prime software. By concatenating these three mtDNA sequences (12S, 16S, COIII), we grouped the samples into haplotypes that were aligned and perfectly matched the five previously described haplotypes ( Bocharova, 2015 ). Including these mitochondrial haplotypes (Pacific 1-5), we used geneHapR ( Zhang et al., 2023 ) package to produce the haplotype network ( Figure 2 ). We then used the extracted nuclear DNA to amplify 41 microsatellites using primer pairs previously developed in Bocharova et. al. (2018) for this species and another sea anemone species of Actiniidae family, Cribrinopsis albopunctata ( Sanamyan & Sanamyan 2006 , Sanamyan et al. 2019 ) with the corresponding PCR conditions. PCR products were visualized in 3500 Genetic Analyzer (Applied Biosystems, USA) using POP7 gel polymer. We analyzed the obtained chromatograms using GenMapper Software (ThermoFisher Scientific, USA). Out of these 41 primer pairs, only ten pairs (Act007, Act011, Act028, Act067, Act173, Act177, Act235, Act238, Act249, Act386,) successfully amplified and were polymorphic in our samples. We then multiplexed the ten microsatellites into four PCR pools using three specific dyes (TAM, ROX, FAM) for visualization: Pool 1) Act007 (FAM), Act011 (TAM); Pool 2) Act238 (FAM), Act249 (TAM), Act173 (ROX); Pool 3) Act067 (FAM), Act177 (ROX), Act386 (TAM); Pool 4) Act28 (TAM), Act235 (ROX). To estimate if this marker set was powerful enough to distinguish between individuals, we computed the overall observed unbiased probability of identity in state under panmixia (hereafter overall pidu) and the overall observed unbiased probability of identity in state between sibs (overall pids) as recommended by Waits et al. (2001) . To achieve allele dosage within each individual per locus, we used the dose-effect method with some modifications (Rodzen & Mаy, 2002; Drauch Schreier et аl., 2011). In cases of more than 1 allele, dose-effect is counted following the instruction ( Figure S2 ): two alleles of almost the same height were considered as AABB ( Figure S2A ); One high allele and two low alleles as AABC ( Figure S2B ); One very high allele and one low allele as AAAB ( Figure S2C ); Tree alleles of the same height were considered as ABC + null allele ( Figure S2D ); One high allele and one allele half lower were considered as AAB + null allele ( Figure S2E ). Analyses of genetic diversity within spatio-temporal samples Parents that produce descendants by clonality can result to observe multiple times the same Multi-Locus Genotype (hereafter MLG) among the genotyped samples ( Halkett et al., 2005 ; Arnaud-Haond et al., 2007 ). However, these repeated MLGs may also come from genetic marker lacking sufficient resolution to disentangle between descendants produced by sexual reproduction, either under panmixia (the lower probability to observe two identical sexual descendants) or between mating sibs (the higher probability to observe two identical sexual descendants). To estimate the power of our marker set to distinguish different individuals and assess the upper and lower bounds for the probability of observing identical MLG drawing a pair of individuals sampled in a same population, we therefore computed the probabilities of identity under panmixia ( pidu ) and considering only sib mating ( pids ; Jacquard, 1968 , 1970 ; Waits et al., 2001 ) using formula adapted for tetraploid populations implemented in GenAPoPop software ( Stoeckel et al., 2024 ). To assess the level of genetic diversity within populations and compare between populations, for each age cohort in one site, we reported the sample size ( N ) and computed the average number of observed alleles and its variance Var(Na) over the ten microsatellites, the average effective number of alleles and its variance Var( A E ), the average gene diversity and its variance Var( H E ), the average observed heterozygosity and its variance Var( H 0 ) and the average value of the inbreeding coefficient and its variance Var( F IS ) over all loci. We also computed genotypic indices including the number of unique genotypes (G), the unbiased genotypic richness R ( Dorken & Eckert, 2001 ), the genotypic evenness as Pareto β to assess the potential uneven contribution of some clonal MLGs to the genotypic composition of populations ( Arnaud-Haond et al., 2007 ; Bailleul et al., 2016 ). We searched for clonal signatures in genetic diversity by analyzing and its variance Var( F IS ) and the overall estimate of the linkage disequilibrium ( , Agapow & Burt, 2001 ) as commonly interpreted using one-shoot samples (Halket et al., 2005; Arnaud-Haond et al., 2020 ; Stoeckel et al., 2021 a). Finally, to assess rates of clonality within spatial sites, we used the Bayesian estimate of rates of clonality using genotype frequency changes between adult and juvenile at the same sample site (ClonEstiMate, Becheler et al., 2017 adapted to autpolyploids in GenAPoPop Stoeckel et al., 2024 ). To assess the genetic structure between populations and generations, we computed temporal and spatial pST ( Ronfort et al., 1998 ) using GenAPoPop, and plotted the unrooted minimum spanning tree between individuals using their pairwise number of shared alleles (identity-in-state) using the classical equal-angle algorithm ( Felsenstein, 2004 ) to picture the coancestry between individuals across populations at our marker set. To assess the dynamics (speed and direction) of genetic diversity along sampling sites, we compared juvenile and adult subpopulations sampled at the same sites and between sites computing the corresponding pST . We also used a Discriminant Analysis of Principal Components (DAPC, Jombart et al., 2010 ) to investigate the genetic structure of our five sites and tens cohorts. We chose to retain 80 principal component eigenvalues. We used a knee locator approach to find the minimal number of clusters best structuring the differentiation between individuals and groups. We verified this value by increasing and decreasing the number of clusters and checking that we obtained the same groups and distribution of clusters among sites and cohorts. Statistical analyses To test if distribution of population genetic indices were significantly different between Arctic-Atlantic and Pacific populations and between juveniles and adults, we used Mann-Whitney test and reported the median values of the two groups ( median group1 versus median group2 ) in the order they are mentioned, the sum of the rank of the first group ( U ), the number of measures in each group ( n1=4 for the cohorts sampled on the Pacific side of the Arctic Ocean and n2=6 for the cohorts sampled on the Atlantic side) and the probability value ( p ) of the hypothesis that the two groups have identical distribution of their genetic indices. To compare pST between juveniles and adults among pairs of the same sites, we used Wilcoxon signed rank test. We reported the median pST value of the two groups ( median group1 versus median group2 ) in the order they are mentioned, the sum of the ranks of the differences above zero ( W ) and the probability value ( p ) of the hypothesis that the two age cohorts have identical distribution of pST . All these tests were performed using python 3.11 and Scipy.stats 1.13 ( Virtanen et al., 2020 ). Results We genotyped 314 individuals of Aulactinia stella sampled on five sites, including a total of 75 adults, ranging from 9 to 21 with a median of 15 adults per site, and a total of 239 juveniles, ranging from 6 to 93 with a median of 49 juveniles per site ( Table 1 ). We found both females and males in the Kamchatka population, whereas we found only females in the Kuril Islands, Newfoundland, the Barents and the White Seas. For subsequent genetic analyses, we only used the adult females from each site. View this table: View inline View popup Download powerpoint Table 1: Summary of population genetic indices in ten spatio-temporal populations of Aulactinia stella . For each age cohort at one site, are reported the sample size N , number of genotypes G , genotypic richness R , genotypic evenness Paretoβ , overall estimate of linkage disequilibrium , average number of observed alleles Na and its variance V(Na) over the ten microsatellites, average effective number of alleles and its variance V( A E ), average observed heterozygosity and its variance V( H 0 ), average gene diversity and its variance V( H E ), average value of the inbreeding coefficient and its variance V( F IS ) over all loci, unbiased probability of identity under panmixia pidu and between sibs pids . We found 15 brooding females in White Sea, 21 in Barents Sea, 14 in Newfoundland, 9 in Kamchatka and 19 in Kuril Islands sampled sites. The oral disc diameter of mature females varied from 10 to 45 mm while the overall size of the juveniles released from the gastrovascular cavities of the adults varied from 1 to 5 mm. Gene, allele and genotypic diversities in nuclear microsatellites Gene, allele and genotypic diversities in nuclear microsatellites were consistent over generations ( Table 1 ). We found a maximum of four different alleles for each locus in one individual over all our samples and concluded that A. stella must be at least an autotetraploid species. Using our ten microsatellites, the mean numbers of alleles per locus per site and developmental stage ranged from 1.2 to 5.7, with a median of 2.75. The mean effective numbers of alleles per locus per site and developmental stage ranged from 1.19 to 2.53 with a median of 1.69. Over all samples, probability of identities under panmixia was 8.3×10 -9 and probability of identity between sibs was 4.5×10 -3 . These low probabilities that two individuals generated by sexual reproduction would result into similar genotypes advocated that the ten microsatellites we used were powerful enough to genetically distinguish between sexual descendants, even resulting from a sib mating ( Waits et al. 2001 ). We thus considered afterward that individuals with identical MLGs were individuals of the same clone (i.e., ramets of a same genet). Over all sites and cohorts, we found 154 unique MLGs and 26 of these MLGs were found repeated in more than one sampled individual, i.e. genets with more than one ramet ( Table S1 ). In Kuril Islands, we genotyped no repeated MLG. In Kamchatka, we found three individuals of the same clone. In Barents Sea, we found three repeated MLGs, including respectively two, five and thirteen ramets. In White Sea, we found eight repeated MLGs involving two to forty ramets, for a total of 88 ramets. In Newfoundland, we found fourteen repeated MLGs involving two to eleven ramets, for a total of 49 ramets. All ramets of repeated MLGs were exclusively sampled on the same site. In all populations with repeated MLGs, repeated MLGs were found both in adults and in juveniles. Some of these genets were even shared by adults and juveniles, except in the Newfoundland samples. We found higher genetic and genotypic diversity, and lower probabilities of identity on the Pacific side sites (i.e., Kuril and Kamchatka) than found on the Atlantic side of the Arctic Ocean (i.e., Newfoundland, Barents and White Seas), even if differences among sites within oceans still remained on some indices ( Table 1 ). The four cohorts sampled in Kuril Islands and in Kamchatka showed nearly twice more average gene diversity, ranging from 0.372 to 0.460, than found in the six cohorts sampled in Newfoundland, Barents and White Seas, that ranged from 0.097 to 0.266 (0.427 versus 0.174, U=24, n1=4, n2=6, p=0.009), while they both showed similar standard deviation of gene diversity over all loci (U=15, n1=4, n2=6, p=0.610; Table 1 ). Cohorts in the Pacific sites also showed significantly higher mean numbers of alleles per locus (4.9 versus 2.0, U=23.5, n1=4, n2=6, p=0.019) and higher effective number of alleles (2.43 versus 1.43, U=24, n1=4, n2=6, p=0.009). Kamchatka and Kuril samples showed higher genotypic richness measured as R than populations on the Atlantic side (1.00 versus 0.37, U=24, n1=4, n2=6, p=0.009) and higher genotypic evenness measured as Pareto β values (3.64 versus 0.45, U=23.5, n1=4, n2=6, p=0.019). Over all loci, Kamchatka and Kuril samples showed lower probabilities of identity under panmixia (<10 -6 versus 0.003, U=0, n1=4, n2=6, p=0.009) and between sibs (0.008 versus 0.118, U=0, n1=4, n2=6, p=0.009) than found in samples from the Atlantic side. Heterozygosity and deviation from the genotypic proportions expected under Hardy-Weinberg hypotheses We found no difference in observed heterozygosity between cohorts on the Pacific and Atlantic sides (0.241 versus 0.181, U=15, n1=4, n2=6, p=0.61), though variance of observed heterozygosity between loci was nearly two times higher in cohorts on the Atlantic side (0.389 versus 0.701, U=1, n1=4, n2=6, p=0.019). Mean Fis over all loci were significantly lower in Newfoundland, the Barents and the White Seas than in Kamchatka and Kurils (0.406 versus - 0.066, U=24, n1=4, n2=6, p=0.009) and even negative in the juvenile cohort of Barents Sea, and in adults and juveniles of White Sea. Variances of Fis values between loci were significantly higher, three to ten times, in cohorts on the Atlantic side (1.08 versus 6.41, U=24, n1=4, n2=6, p=0.009), than found in populations on the Pacific side. Though linkage disequilibrium seemed higher in cohorts on the Atlantic side, they were not significantly different from values measured in cohorts on the Pacific side (0.051versus 0.356, U=5, n1=4, n2=6, p=0.171). Estimated rates of clonality using temporal genotype frequencies Credible intervals including 99% of posterior probabilities of rates of clonality estimated from temporal samples were quite sharp, spanning around their peaks on a range of 0.06 to maximum 0.25 ( Figure 1 ). Clonality was intermediate in the Kuril Islands and Kamchatka, but very high in the Barents Sea, the White Sea, and Newfoundland. Download figure Open in new tab Figure 1: Posterior density probabilities of estimated rates of clonality using changes of genotype frequencies over one generation for each sampled site using ClonEstiMate method for autopolyploids in GenAPoPop software (Becheler et al., 2007; Stoeckel et al., 2024 ) on the ten nuclear microsatellites. Kuril Islands in blue, Kamchatka in green, Newfoundland in red, White Sea in pink and Barents Sea in gold. Genetic structure between sites using mitochondrial gene sequences Using three gene sequences (12S, 16S, COIII), we found that all the 82 adults sampled on the Atlantic side (Newfoundland, the White and the Barents Seas) shared the same mitochondrial haplotype, Pacific-3 (classification according to Bocharova, 2015 ; GenBank Accession Numbers: 12S: KT310190 , 16S: KT310200 , COIII: KT310212 ). On the Pacific side, all the 20 individuals sampled in the Kuril Islands belonged to another mitochondrial haplotype, Pacific-1 (GenBank Accession Numbers: 12S: KT310188 , 16S: KT310198 , COIII: KT310210 ) while Kamchatka’s samples distributed among five mitochondrial haplotypes ( Figure 2 ). Download figure Open in new tab Figure 2: Haplotype network using three mitochondrial gene sequences (12S, 16S, COIII). The proportion of each sampled site per haplotype is represented in color. The circle sizes of each haplotype are proportional to their number in our dataset. We used haplotype names initially identified in Bocharova (2015) , thus Pacific 1 to 5. We found no new haplotypes. Genetic differentiation between sites and cohorts using nuclear microsatellites Among adults, genetic differentiation between sites ranged from 0.155 to 0.612 while among juveniles, genetic differentiation between sites ranged from 0.105 to 0.530 ( Table S2 ). Within each sampled site, pST between adults and juveniles ranged from 0.009 to 0.088, therefore between six to eleven times lower than between sites. Between sites, adults were more genetically structured than juveniles (0.366 versus 0.245, W=43, p=0.065; Figures 3 , S1 , S3 ). Considering the adults and juveniles together by site, genetic differentiation between sites ranged from 0.096 to 0.521 ( Table S2a ). Globally, in juveniles, pST between sites showed very close values and followed ranks of geographical and along-the-coast distances ( Table S2b , Figure S1 ). In contrast, adults of White Sea appeared highly genetically isolated from the other sites, excepted from juveniles in Newfoundland. pST genetic differentiation hierarchically structured the samples into one Pacific side group and one Atlantic side group then subdivided by sites and then by age cohorts ( Table S2b , Figure 3 ). Download figure Open in new tab Figure 3: Oriented graph of pST values among sites and cohorts (vertices), based on ten nuclear microsatellite loci. Orange arrows (edges) and the pST values in boxes illustrate pairwise genetic differentiation between sites, allowing inference of the dynamics and direction of genetic spread. Comparisons of pST values among adults (blue boxes) and juveniles (green boxes) between site pairs reveal temporal genetic differentiation, even if apparent increases or decreases in genetic differentiation may arise by chance due to sampling variability within the same underlying distribution, or as a result of homoplasy or gene flow between populations. For instance, adults from the Barents Sea are twice as identical in state with Newfoundland juveniles as adults from Newfoundland are with Barents Sea juveniles, suggesting preferential gene flow from the Barents Sea toward Newfoundland. Moreover, genetic differentiation between Newfoundland and Barents Sea adults is approximately twice that observed between juveniles, indicating either enhanced gene flow among juveniles or processes promoting endemism in established adults. This hierarchical genetic structure was congruently supported by the DAPC results ( Figure 4 ). Interestingly, the Bayesian information criteria of the DAPC decreased continuously with the increasing number of clusters, as expected in a spatial stepping-stone model ( Figure S4 ). The minimal number of clusters maximizing the distance between groups was 4. Whatever the number of clusters, DAPC grouped together individuals of Barents Sea and Newfoundland into one common and compact group, while the other sites each stood for their own genetic groups with couples to handful individuals assigned to other groups. Looking at the genetic structure between cohorts and sites, juveniles from the Kuril Islands shared more genetic diversity with adults from Kamchatka than the inverse ( Figure 2 ). Adults of Barents Sea shared more similar genetic diversity with juveniles of Newfoundland and of White Sea than the inverse. DAPC proportion of each individual assigned to one DAPC cluster also supported this dynamic tendency ( Figure 4 ): one descendant sampled in Kuril was assigned to Kamchatka genetic cluster; Three adults sampled in Kamchatka were assigned to Barents-Newfoundland genetic cluster and one adult sampled in Barents Sea was assigned to Kamchatka genetic cluster; Four adults sampled in Barents Sea were assigned to white sea genetic cluster, and one adult and ten descendants sampled in White Sea were assigned to Barents-Newfoundland genetic cluster. Download figure Open in new tab Figure 4: Discriminant Analysis of Principal Components of the genetic structure of Aulactinia stella across the Arctic Ocean using the ten nuclear microsatellites. We found four clusters being the minimal number of clusters maximizing the distance between groups (see Figure S4 ). Individuals are represented as columns corresponding to the membership probability to the four supported clusters (in color). Sampled sites are represented successively, separated by a large white bar. Within each site, adults are on the left side of the dashed vertical black line and juveniles on the right. Similarly, the minimum spanning tree of pairwise genetic distances between individuals showed two main groups, one including Kurils and Kamchatka (Pacific side group), and one including Newfoundland, Barents and White Seas (Atlantic side group, Figure 5 ). Atlantic side group showed more compact genetic distances between individuals than found in the Pacific side group. Atlantic side group also showed more star and chaplet topological subgraphs, commonly observed in clonal lineages ( Stoeckel et al., 2024 ). In accordance with pST , juveniles were more centrally connected in this network and less genetically distant to the other populations than the adults within each site. Some pairwise genetic distances between individuals within site, like in Kamchatka, were larger than between two individuals from different sites ( Figure 5 , e.g. Kamchatka’s individuals linked to White and Barents seas individuals). Download figure Open in new tab Figure 5: Minimum spanning tree of the genetic differentiation between sampled sites and cohorts of A.stella using pairwise p ST using the ten nuclear microsatellites. Kamchatka adults and juveniles respectively in green and pale green; Kuril adults and juveniles respectively in dark blue and blue; Barents sea adults and juveniles respectively in gold and pale gold colors; White sea adults and juveniles respectively in pink and pale pink; Newfoundland adults and juveniles respectively in dark red and red. Discussion We investigated reproductive modes, genetic diversity and differentiation in five populations of the autotetraploid sea anemone Aulactinia stella , sampled from North Pacific, Arctic, and North Atlantic coasts using three mitochondrial markers and ten nuclear microsatellites. A. stella is broadly distributed throughout the northern hemisphere where it is an important semi-sessile predator in shallow-water communities and prey for higher trophic levels ( Ottaway, 1977 ). We genotyped two successive age cohorts from each of the five sites. We benefited from past sampling campaigns (from 2012 to 2017) to genotype individuals, with fewer sites and greater temporal dispersion in the Pacific Ocean: This could conceivably influence our results. The ten nuclear microsatellites we used offered sufficient power to distinguish clones from inbred offspring. Despite our standardized and reproducible protocol for allele dosage inference ( Figure S2 ), tetraploid microsatellite datasets remain sensitive to technical biases, particularly those arising from null alleles and amplification biases, which may affect quantitative estimates prior to downstream analyses. Genetic differentiation in A. stella revealed a hierarchical structure: first separating into a Pacific and an Arctic-Atlantic group, each of which was further subdivided by sampling site, and then by age cohort. This type of structure closely resembles patterns seen in the intertidal macroalga Fucus distichus and in other Arctic and peri-Arctic sessile coastal marine species ( Coyer et al., 2011 ; Bringloe et al., 2020 ). Environmental factors, rather than species-specific biological or behavioral traits, may thus likely be the primary drivers of genetic differentiation in this region. We therefore examined reproductive modes and assessed genetic diversity and differentiation across these hierarchical groups. Reproductive modes in Aulactinia stella populations We found strong genetic evidence that A. stella reproduces clonally in natural populations, consistent with the absence of males in some sites and with previous observations of embryos developing parthenogenetically within parental mesenteries ( Bocharova, 2012 , 2015 ; Kaliszewicz et al., 2012 ). In the North Pacific, repeated multilocus genotypes (MLGs) were rare; however, genetic indices, including interlocus FIS variances close to 1 and low levels of linkage disequilibrium, indicated intermediate levels of clonality (Navascues et al., 2010; Stoeckel et al., 2021 a,b, 2024 ). It was further supported by changes in genotype frequencies between adults and juveniles, which yielded clonality estimates of 0.37 to 0.65 ( Becheler et al., 2017 ; Stoeckel et al., 2024 ). The high clonal evenness in these populations likely reduced the probability of resampling identical MLGs ( Arnaud-Haond et al., 2007 , Stoeckel et al., 2021 a). In contrast, Arctic and North Atlantic populations, characterized by lower allelic diversity, exhibited extensive clonality. Many individuals shared identical MLGs, and population-level genetic indices, including negative mean FIS, large variance of FIS across loci, and significant linkage disequilibrium, indicated high rates of clonal reproduction ( Halkett et al., 2005 ; De Meeûs et al., 2006 ; Navascues et al., 2010; Stoeckel et al., 2021 a). Consistently, clonality rates inferred from changes in genotype frequencies between adults and juveniles ranged from 0.74 to 1.0 in Arctic and North Atlantic populations, nearly twice the rates estimated for the North Pacific. The advantages of clonal reproduction at the limits of species’ ranges High levels of uniparental reproduction, including clonality, are commonly observed in species expanding into new or marginal environments (Barrett, 2015; Pannel et al., 2015; Burke & Bonduriansky, 2019 , Pereyra et al., 2023 ), particularly when enduring harsh conditions ( Klimeš 2008 ) and when heading towards the poles ( Peck et al., 1998 ). The post-glacial Baltic sea provides several examples, including the emblematic case of Fucus vesiculosus ( Johannesson & André 2006 ; Pereyra et al., 2009 ; 2023 ). In plant species, the persistence of connection between ramets (mainly by roots) offers ecological advantages that may alone explain the prevalence of clonality at range margins, including regeneration after injury, distribution of physiological functions among ramets, translocation of resources between interconnected ramets (clonal integration) and space monopolization using different strategies (including phalanx versus guerrilla) to address intra- and inter-specific competitions depending on the environmental context ( Cadotte et al., 2006 ; Gonzalvez et al., 2008; Roiloa, 2019 ; Franklin et al., 2021 ). In contrast, dispersed clones in marine species such as A. stella do not benefit from such clonal integration. In this context, clonality should only favor kin cooperation (e.g., through brooding or cooperative behavior Bocharova, 2012 ; Rudin & Briffa, 2012 ) but more importantly it should favor the persistence of local populations under harsh conditions by avoiding the energetic and biological costs of sexual reproduction ( Klimeš, 2008 ; Pierre et al., 2022 ). Baker (1955) postulated that uniparental reproduction (including clonality) provides a critical advantage for colonization: individuals arriving in new areas can reproduce clonally to establish a new area. First, clonality allows a population to maintain itself locally until compatible mates become available ( Baker & Stebbins, 1965 , Daly et al., 2023 ). This hypothesis is particularly relevant for gonochoristic (or dioecious) sessile organisms that require close proximity for successful sexual reproduction (Barrett, 2015), such as in sea anemones. Second, in new harsh environments, clonal individuals that manage to survive minimize the risks of genetic reassortment associated with sexual reproduction for their offspring ( Peniston et al., 2021 ; Pierre et al., 2022 ). Third, clonality maintains higher effective population sizes than in identical more sexual populations ( Reichel et al. 2016 ; Stoeckel et al., 2021 b). Clonality limits the negative effects of genetic drift during founder events, including the fixation of deleterious alleles in small populations (Navascues et al., 2010; He et al., 2024 ). Our results in A. stella , which indicated a twofold increase in clonality at the species’ range edge, suggest that clonal reproduction may play an important role in supporting population persistence and maintaining higher effective population sizes and the more adapted genotypes under unfavorable conditions during range expansion (Baker’s conjecture, Pannel et al., 2015; Rafajlović et al., 2017 ; geographic parthenogenesis, Peck et al. 1998 ), without necessarily relying on mechanisms such as clonal integration (as proposed in plants) or habitat modification (as suggested for ecological engineering species) to explain clonal success. Aulactinia stella spreading across the Arctic Ocean from north Pacific? Genetic diversity was nearly twice as high in North Pacific populations compared to Atlantic ones. Kamchatka and Kuril Islands harbored the greatest allelic and genotypic diversity, whereas Newfoundland, Barents Sea, and White Sea populations showed more compact identities in state. All Arctic and North Atlantic individuals carried the Pacific-3 haplotype, rare in the Pacific, while Kuril populations were fixed for Pacific-1, the dominant haplotype in Kamchatka. This distribution, considering our samples, suggests a North Pacific origin, with Kamchatka as a likely glacial refugium ( He et al., 2024 ). In contrast, Arctic and North Atlantic populations may be more recently established (since the last glaciation) or alternatively, may have originated from smaller, more localized glacial refugia, evolving relatively isolated since the last warmer interglacial period (Eemian, minimum ∼116,000 years ago). But the DAPC results ( Figure S4 ) and pairwise genetic distance ( Figure 5 ) supported a stepping-stone genetic structure between sites and genetic groups rather than isolated as expected in an island model with no migration ( Jombart et al., 2010 ). Interestingly, the path of the lowest pST genetic structure between sites align with major summer currents (e.g., Anadyr, Beaufort gyre, transpolar drift), areas of low ice cover, and coastal distances ( Rudels & Carmack, 2022 ; Timmermans & Marshall, 2020 ; Constantin & Johnson, 2023 ) and increasing Arctic maritime traffic that may further accelerate dispersal, either through ballast water transport or passive enhancement by ship wakes ( Melia et al., 2016 ; Boylan, 2021 ; PAME, 2024 ). Alterations in Arctic Ocean conditions, together with increasing human activities (including anthropogenic contributions to floating debris, ship-generated wake currents, and maritime transport of materials) may facilitate connectivity among A. stella populations, potentially via stepping-stone dispersal along coasts during the free-living polyp stage. Since at least 1993, temperatures in the upper 10 m of the water column have commonly reached 8–10 °C in coastal and sheltered bays, and 2–3 °C in the open sea, with higher temperatures in coastal and tidal currents along the Eurasian Arctic coasts, making this migration physiologically feasible ( Locarnini et al. 2024 ; E.U. Copernicus Marine Service, 2024 ). These maritime and tidal currents are estimated to transport approximately 1 Sverdrup, a volume comparable to the combined mean discharge of all the world’s rivers, but substantially lower than the Gulf Stream, which rather transports 60–90 Sverdrups Constantin & Johnson, 2023 ). Taken together, these factors may explain the reduced genetic differentiation observed among juveniles within the established stepping-stone connectivity pattern of A. stella populations. Similar genetic patterns along with evidence for human-mediated dispersal have been reported for the coastal sea anemone Diadumene lineata ( Glon et al., 2020 ), which, however, is more fecund and lacks the brooding behavior observed in A. stella . Nevertheless, the biology of A. stella , like that of many Arctic marine organisms, remains poorly understood once individuals move beyond nearshore habitats and become inaccessible to direct observation. In particular, dedicated studies are needed to determine whether this Arctic sea anemone can persist while attached to debris or as free-living polyp or adult drifters within coastal and tidal currents, potentially enabling effective migration across the Arctic Ocean. Juveniles less genetically structured than adults: methodological artifacts or gene flow with local endemism? Comparisons of genetic diversity between adults and juveniles revealed slightly lower levels of genetic differentiation among juveniles ( Figure S3 ). Several non-exclusive hypotheses could explain this pattern, requiring further research to disentangle their relative contributions. Methodological artefacts must first be ruled out through more geographically balanced and systematically designed sampling. Using polymorphic markers with lower mutation rates could also help distinguish genuine shared ancestry from homoplasy that has higher chance to occur in microsatellite loci across distant populations. Clarifying the life cycle of A. stella , particularly the dynamics of clonal interference and mortality during recruitment, will be essential to assess whether reduced genetic structure in juveniles and the absence of shared multilocus genotypes (MLGs) across sites reflect local genetic endemism, gregarious behavior promoting shared ancestry (as found in other sea anemone species: Lane et al., 2020 ; Ryan et al., 2021 ), or persistent founder effects with limited gene flow. A. stella occupies a range of habitats, such as the Barents and White Seas, which differ significantly in salinity despite their proximity (∼100 km). Environmental factors, including temperature, salinity, depth, tidal exposure, and substrate disturbance, are known to influence reproductive modes and genetic differentiation ( Eriksson & Johansson, 2005 ; Krueger-Hadfield et al., 2013 ; Berković et al., 2014 ; Dudgeon et al., 2017 ; McMahon et al., 2017 ; Ryan & Miller, 2019 ; Williams et al., 2024 ). Similar habitat-driven patterns are observed in other sea anemones, such as Actinia tenebrosa and Diadumene lineata , where localized clonal distributions and endemic lineages correlate with temperature, salinity, UV exposure, and community composition ( Sherman & Ayre, 2008 ; Ryan et al., 2021 ). It may explain the higher genetic differentiation found in adults between sites. Restricted dispersal or local adaptation in Aulactinia stella ? No MLGs were shared between sampling sites, even over relatively shorter distances (e.g., Barents vs. White Sea, ∼360 km; Kuril Islands vs. Kamchatka, ∼680 km). The presence of local, site-specific clones and high levels of genetic differentiation between adult populations is common in cold-water sea anemones, including A. tenebrosa , Anthothoe albocincta , Diadumene lineata , and several North Pacific species ( Veale & Lavery, 2012 ; Sherman et al., 2007 ; Ryan et al., 2021 ; Edmands & Potts, 1997 ), which contrasts with tropical species that often share genets across hundreds of kilometers ( Russo et al., 1994 ; Emms et al., 2020 ). Such patterns would suggest that Arctic conditions may constrain and challenge the dispersal of sea anemones ( Chan et al., 2019 ), but other passively dispersing taxa in cold-water habitats such as macroalgae, seagrasses, and bivalves, show weaker population structure ( Coolen et al., 2020 ; Pereyra et al., 2023 ; Krueger-Hadfield et al., 2021 ). In fact, polyp stage commonly allows locally-active and passive transport by tidal flows and circumpolar maritime currents, facilitating long-distance dispersal ( Riemann-Zürneck, 1998 ). For instance, Urticina eques exhibits very low differentiation over hundreds of kilometers along the British coast ( Cava-Solé et al., 1994 ). In the case of A. stella , Arctic populations may have only narrow temporal windows of suitable thermal conditions for successful sexual reproduction, particularly in the intertidal and shallow sublittoral zones. These habitats are also characterized by extreme variability in salinity and desiccation at low tide. In such habitats, A. stella is often found hiding under Fucus algae, in rock crevices, or in tide pools, behaviors that indicate exposure to environmental stress. Local environmental filtering through selection and post-settlement mortality (e.g., Fucus vesiculosus : Johannesson & André, 2006 ; Pereyra et al., 2009 ; Fucus distichus : Laughinghouse et al., 2015 ; Actinia tenebrosa : Sherman et al., 2008 ; seagrasses: Kendrick et al., 2012 ; Acropora millepora : van Oppen et al., 2011 ), as well as individual habitat preferences and gregarious settlement (e.g., Diadumene lineata : Ryan et al., 2021 ), may explain the roughly twofold higher genetic differentiation observed among adults compared to juveniles in A. stella . These patterns highlight the need for further ecological and behavioral studies to clarify the potential mechanisms and temporal dynamics of dispersal in this species in the Arctic Ocean. Data Availability All data (including adults, juveniles, MLG genotypes, sampling GPS locations, and clone sizes) as well as the scripts used to reproduce the results and figures, are openly available on Zenodo https://doi.org/10.5281/zenodo.17193174 Supplementary Material Supplementary material provides two tables and three figures. Table S1 reports distributions of repeated MLGs within sites and MLGs shared between adults and juveniles; Table S2 reports pairwise pST between sites and between age cohorts; Figure S1 reports a map of the sample sites with the minimum spanning trees of genetic differentiation of adults, juveniles and overall among sites; Figure S2 reports how allele dosage of the microsatellites was achieved; Figure S3 reports the change of the Bayesian information criteria of the DAPC with increasing number of genetic clusters arguing for a stepping-stone genetic structure. Funding This work was funded by the Russian Foundation for Basic Research (Grant No. 16-04-01685 to E.S. Bocharova), by the French National Research Agency (project Clonix2D ANR-18-CE32-0001 to S. Stoeckel) and the French Embassy in the Russian Federation (Metchnikov Stipend 2019 to E.S. Bocharova). Author contributions EB and SS laid the foundation of this work while EB visited SS lab in 2018 and 2019. EB and SS conceived the study and were responsible for funding applications. EB was responsible for data acquisition. She supervised and achieved sample collection, determined the species and sexed the sampled individuals. She also conceived, developed the sequences and microsatellites, supervised and achieved sequencing and genotyping of all the samples. From sequencing machines, she genotyped the samples. AV contributed to methodology, supervision and administration of data acquisition. EB performed data analyses and produced figure on mitochondrial sequences. SS performed data analyses and produced tables and figures on microsatellites. SS and EB tracked and managed the bibliography, interpreted results, wrote the manuscript and managed its edition. All authors read and approved this manuscript. Conflict of interest disclosure The authors declare that they comply with the PCI rule of having no financial conflicts of interest in relation to the content of the article. Supplementary information Download figure Open in new tab Download figure Open in new tab Figure S1: Minimum spanning tree of A). the pST of adults, B). the pST of juveniles and C). the pST of all age cohorts over all loci between sampled sites. Download figure Open in new tab Figure S2: Counting of allele dosage of the microsatellite locus 249TAM (as an example) in Aulactinia stella (within each individual per locus): A). 129, 129, 133, 133 (AABB), B). 129, 129, 137, 141 (AABC), C). 137, 137, 137, 141 (AAAB), D). 129, 133, 137 (ABC+null allele), E). 133, 133, 137 (AAB+null allele). Download figure Open in new tab Figure S3: Distributions of pST between pairs of sites among juveniles and adults (light blue violin plots). We reported the probability that difference of pST between juveniles and adults would be due to a random (sampling, locus) effects as a Wilcoxon test (significance bar), the pST values in each age cohort (Black dots), the extremum and mean values (dark blue bars). Download figure Open in new tab Figure S4: Change of the Bayesian information criteria (BIC) with increasing number of genetic clusters used in the Discriminant Analysis of Principal Components to investigate the genetic differentiation between individuals and among sampled sites and cohorts. The continuously decreasing pattern observed is typical of a stepping-stone structure of the studied demes. View this table: View inline View popup Download powerpoint Table S1: Distribution of repeated MLGs within sites and MLGs shared between adults and juveniles. Genet ID, the number of ramets within the genet and if the considered genet has been found in adults, juveniles or both are reported. No repeated MLG was found in Kuril islands sampled genotypes. View this table: View inline View popup Download powerpoint Table S2: Pairwise pST between sites (A) and between age cohorts within sites and among sites (B). Acknowledgements The authors thank N.P. Sanamyan, K.E. Sanamyan, B.V. Osadchenko and I.A. Kosevich for providing material and assistance in species identifications, N.S. Mugue, T.V. Neretina, A.E. Barmintseva and S.D. Grebelnyi for professional consultations, A. Mercier, J.-F. Hamel, E. Montgomery and the Ocean Sciences Field Services of Memorial University for their assistance with the animal collections. We thank Patricia Nadan and Pascale Le Névé (INRAE IGEPP) for their administrative support along this project. A preprint version of this article has been peer-reviewed and recommended by PCI Evol Biol ( https://doi.org/10.24072/pci.evolbiol.100866 ). We warmly thank Benoit Nabholz, Kerstin Johannesson, Malgorzata Lagisz and one anonymous reviewer for their insightful comments that helped improving our manuscript and its data availability and reproducibility. Funder Information Declared French National Research Agency , project Clonix2D ANR-18-CE32-0001 French Embassy in the Russian Federation , Metchnikov Stipend 2019 Russian Foundation for Basic Research , Grant No. 16-04-01685 Footnotes This is the final accepted, peer-reviewed and recommended by PCI Evol Biol. “A preprint version of this article has been peer-reviewed and recommended by PCI Evol Biol ( https://doi.org/10.24072/pci.evolbiol.100866 ).” https://doi.org/10.5281/zenodo.15102245 References ↵ Agapow P-M , Burt A ( 2001 ) Indices of multilocus linkage disequilibrium . Molecular Ecology Notes , 1 , 101 – 102 . doi: 10.1046/j.1471-8278.2000.00014.x OpenUrl CrossRef Web of Science ↵ Arnaud-Haond S , Duarte CM , Alberto F , Serrão EA ( 2007 ) Standardizing methods to address clonality in population studies . Molecular Ecology , 16 , 5115 – 5139 . doi: 10.1111/j.1365-294X.2007.03535.x OpenUrl CrossRef PubMed Web of Science ↵ Arnaud-Haond S , Stoeckel S , Bailleul D ( 2020 ) New insights into the population genetics of partially clonal organisms: When seagrass data meet theoretical expectations . Molecular Ecology , 29 , 3248 – 3260 . doi: 10.1111/mec.15532 OpenUrl CrossRef ↵ Avise JC ( 2015 ) Evolutionary perspectives on clonal reproduction in vertebrate animals . Proceedings of the National Academy of Sciences , 112 , 8867 – 8873 . doi: 10.1073/pnas.1501820112 OpenUrl Abstract / FREE Full Text ↵ Bailleul D , Stoeckel S , Arnaud-Haond S ( 2016 ) RClone: a package to identify MultiLocus Clonal Lineages and handle clonal data sets in R . Methods in Ecology and Evolution , 7 , 966 – 970 . doi: 10.1111/2041-210X.12550 OpenUrl CrossRef ↵ Baker HG ( 1955 ) Self-compatibility and establishment after ‘“long-distance” dispersal . Evolution , 9 , 347 – 349 . doi: 10.1111/j.1558-5646.1955.tb01544.x OpenUrl CrossRef ↵ Baker HG , Stebbins GL ( 1965 ) The genetics of colonizing species . New York. Barrett SCH (2015) Foundations of invasion genetics: the Baker and Stebbins legacy. Molecular Ecology , 24 , 1927 – 1941 . doi: 10.1111/mec.13014 OpenUrl CrossRef ↵ Becheler R , Masson J-P , Arnaud-Haond S , Halkett F , Mariette S , Guillemin M-L , Valero M , Destombe C , Stoeckel S ( 2017 ) ClonEstiMate, a Bayesian method for quantifying rates of clonality of populations genotyped at two-time steps . Molecular Ecology Resources , 17 , e251 – e267 . doi: 10.1111/1755-0998.12698 OpenUrl CrossRef PubMed ↵ Berković B , Cabaço S , Barrio JM , Santos R , Serrão EA , Alberto F ( 2014 ) Extending the life history of a clonal aquatic plant: Dispersal potential of sexual and asexual propagules of Zostera noltii . Aquatic Botany , 113 , 123 – 129 . doi: 10.1016/j.aquabot.2013.10.007 OpenUrl CrossRef ↵ Bocharova ES ( 2012 ) Osobennosti razmnozheniya Aulactinia stella ( Verrill, 1864 ) (Hexacorallia, Actiniidae) v Belom i Barentsevom moryakh [Reproduction characters of Aulactinia stella ( Verrill, 1864 ) (Hexacorallia, Actiniidae) in the White and the Barents Seas]. PhD thesis, Moscow State University , Moscow . ↵ Bocharova ES ( 2015 ) Reproductive biology and genetic diversity of the Sea Anemone Aulactinia stella ( Verrill, 1864 ). Hydrobiologia , 759 , 27 – 38 . doi: 10.1007/s10750-015-2250-0 OpenUrl CrossRef ↵ Goffredo S , Dubinsky Z Bocharova E ( 2016 ) Reproduction of Sea Anemones and Other Hexacorals . In: The Cnidaria, Past, Present and Future (eds Goffredo S , Dubinsky Z ), pp. 239 – 248 . Springer International Publishing , Cham . doi: 10.1007/978-3-319-31305-4_15 OpenUrl CrossRef ↵ Bocharova ES , Kozevich IA ( 2011 ) Modes of reproduction in sea anemones (Cnidaria , Anthozoa). Biology Bulletin , 38 , 849 – 860 . doi: 10.1134/S1062359011090020 OpenUrl CrossRef ↵ Bocharova ES , Mugue NS ( 2012 ) Sea anemones Aulactinia stella ( Verrill, 1864 ) (Hexacorallia, Actiniidae) can brood offspring from other individuals of the same species. Doklady Biological Sciences , 444 , 173 – 175 . doi: 10.1134/S0012496612030040 OpenUrl CrossRef PubMed ↵ Bocharova ES , Sergeev AA , Volkov AA ( 2018 ) Identification of microsatellite loci in sea anemones Aulactinia stella and Cribrinopsis albopunctata (family Actiniidae) . F1000Research , 7 , 232 . doi: 10.12688/f1000research.13724.1 OpenUrl CrossRef ↵ Boylan BM ( 2021 ) Increased maritime traffic in the Arctic: Implications for governance of Arctic sea routes . Marine Policy , 131 , 104566 . doi: 10.1016/j.marpol.2021.104566 OpenUrl CrossRef ↵ Brandt S , Wassmann P , Piepenburg D ( 2023 ) Revisiting the footprints of climate change in Arctic marine food webs: An assessment of knowledge gained since 2010 . Frontiers in Marine Science , 10 . doi: 10.3389/fmars.2023.1096222 OpenUrl CrossRef ↵ Bridle J , Hoffmann A ( 2022 ) Understanding the biology of species’ ranges: when and how does evolution change the rules of ecological engagement? Philosophical Transactions of the Royal Society B: Biological Sciences , 377 , 20210027 . doi: 10.1098/rstb.2021.0027 OpenUrl CrossRef PubMed ↵ Bringloe TT , Verbruggen H , Saunders GW ( 2020 ) Unique biodiversity in Arctic marine forests is shaped by diverse recolonization pathways and far northern glacial refugia . Proceedings of the National Academy of Sciences , 117 , 22590 – 22596 . doi: 10.1073/pnas.2002753117 OpenUrl Abstract / FREE Full Text ↵ Burke NW , Bonduriansky R ( 2018 ) The geography of sex: sexual conflict, environmental gradients and local loss of sex in facultatively parthenogenetic animals . Philosophical Transactions of the Royal Society B: Biological Sciences , 373 , 20170422 . doi: 10.1098/rstb.2017.0422 OpenUrl CrossRef PubMed ↵ Burke NW , Bonduriansky R ( 2019 ) The paradox of obligate sex: The roles of sexual conflict and mate scarcity in transitions to facultative and obligate asexuality . Journal of Evolutionary Biology , 32 , 1230 – 1241 . doi: 10.1111/jeb.13523 OpenUrl CrossRef PubMed ↵ Cadotte MW , Murray BR , Lovett-Doust J ( 2006 ) Ecological patterns and biological invasions: using regional species inventories in Macroecology . Biological Invasions , 8 , 809 – 821 . doi: 10.1007/s10530-005-3839-4 OpenUrl CrossRef Web of Science ↵ Cava-Solé AM , Thorpe JP , Todd CD ( 1994 ) High genetic similarity between geographically distant populations in a sea anemone with low dispersal capabilities . Journal of the Marine Biological Association of the United Kingdom , 74 , 895 – 902 . doi: 10.1017/S0025315400090123 OpenUrl CrossRef ↵ Chan FT , Stanislawczyk K , Sneekes AC , Dvoretsky A , Gollasch S , Minchin D , David M , Jelmert A , Albretsen J , Bailey SA ( 2019 ) Climate change opens new frontiers for marine species in the Arctic: Current trends and future invasion risks . Global Change Biology , 25 , 25 – 38 . doi: 10.1111/gcb.14469 OpenUrl CrossRef PubMed ↵ Constantin A , Johnson RS ( 2023 ) On the dynamics of the near-surface currents in the Arctic Ocean . Nonlinear Analysis: Real World Applications , 73 , 103894 . doi: 10.1016/j.nonrwa.2023.103894 OpenUrl CrossRef ↵ Coolen JWP , Boon AR , Crooijmans R , van Pelt H , Kleissen F , Gerla D , Beermann J , Birchenough SNR , Becking LE , Luttikhuizen PC ( 2020 ) Marine stepping-stones: Connectivity of Mytilus edulis populations between offshore energy installations . Molecular Ecology , 29 , 686 – 703 . doi: 10.1111/mec.15364 OpenUrl CrossRef ↵ Coyer JA , Hoarau G , Van Schaik J , Luijckx P , Olsen JL ( 2011 ) Trans-Pacific and trans-Arctic pathways of the intertidal macroalga Fucus distichus L. reveal multiple glacial refugia and colonizations from the North Pacific to the North Atlantic . Journal of Biogeography , 38 , 756 – 771 . doi: 10.1111/j.1365-2699.2010.02437.x OpenUrl CrossRef Web of Science ↵ Daly EZ , Chabrerie O , Massol F , Facon B , Hess MCM , Tasiemski A , Grandjean F , Chauvat M , Viard F , Forey E , Folcher L , Buisson E , Boivin T , Baltora-Rosset S , Ulmer R , Gibert P , Thiébaut G , Pantel JH , Heger T , Richardson DM , Renault D ( 2023 ) A synthesis of biological invasion hypotheses associated with the introduction–naturalisation–invasion continuum . Oikos , 2023 , e09645 . doi: 10.1111/oik.09645 OpenUrl CrossRef ↵ De Meeûs T , Lehmann L , Balloux F ( 2006 ) Molecular epidemiology of clonal diploids: A quick overview and a short DIY (do it yourself) notice. Infection , Genetics and Evolution , 6 , 163 – 170 . doi: 10.1016/j.meegid.2005.02.004 OpenUrl CrossRef PubMed Web of Science ↵ De Meeûs T , Prugnolle F , Agnew P ( 2007 ) Asexual reproduction: Genetics and evolutionary aspects . Cellular and Molecular Life Sciences , 64 , 1355 – 1372 . doi: 10.1007/s00018-007-6515-2 OpenUrl CrossRef PubMed ↵ De Visser JAGM , Elena SF ( 2007 ) The evolution of sex: empirical insights into the roles of epistasis and drift . Nature Reviews Genetics , 8 , 139 – 149 . doi: 10.1038/nrg1985 OpenUrl CrossRef PubMed Web of Science ↵ Dorken ME , Eckert CG ( 2001 ) Severely reduced sexual reproduction in northern populations of a clonal plant, Decodonverticillatus (Lythraceae) . Journal of Ecology , 89 , 339 – 350 . doi: 10.1046/j.1365-2745.2001.00558.x OpenUrl CrossRef Drauch Schreier A , Gille D , Mahardja B , May B ( 2011 ) Neutral markers confirm the octoploid origin and reveal spontaneous autopolyploidy in white sturgeon, Acipenser transmontanus : Neutral markers and ploidy in white sturgeon . Journal of Applied Ichthyology , 27 , 24 – 33 . doi: 10.1111/j.1439-0426.2011.01873.x OpenUrl CrossRef ↵ Dudgeon S , Kübler JE , West JA , Kamiya M , Krueger-Hadfield SA ( 2017 ) Asexuality and the cryptic species problem . Perspectives in Phycology , 47 – 59 . doi: 10.1127/pip/2017/0070 OpenUrl CrossRef ↵ Duminil J , Fineschi S , Hampe A , Jordano P , Salvini D , Vendramin GG , Petit RJ ( 2007 ) Can population genetic structure be predicted from life-history traits? The American Naturalist , 169 , 662 – 672 . doi: 10.1086/513490 OpenUrl CrossRef PubMed Web of Science ↵ Edmands S , Potts DC ( 1997 ) Population genetic structure in brooding sea anemones ( Epiactis spp. ) with contrasting reproductive modes . Marine Biology , 127 , 485 – 498 . doi: 10.1007/s002270050037 OpenUrl CrossRef ↵ Ellegren H , Galtier N ( 2016 ) Determinants of genetic diversity . Nature Reviews Genetics , 17 , 422 – 433 . doi: 10.1038/nrg.2016.58 OpenUrl CrossRef PubMed ↵ Emms MA , Saenz-Agudelo P , Giles EC , Gatins R , Nanninga GB , Scott A , Hobbs J-PA , Frisch AJ , Mills SC , Beldade R , Berumen ML ( 2020 ) Comparative phylogeography of three host sea anemones in the Indo-Pacific . Journal of Biogeography , 47 , 487 – 500 . doi: 10.1111/jbi.13775 OpenUrl CrossRef ↵ Eriksson BK , Johansson G ( 2005 ) Effects of sedimentation on macroalgae: species-specific responses are related to reproductive traits . Oecologia , 143 , 438 – 448 . doi: 10.1007/s00442-004-1810-1 OpenUrl CrossRef PubMed Web of Science E.U. Copernicus Climate Change Service ( 2025 ) Global Climate Highlights 2024. 17 January 2025. https://climate.copernicus.eu/global-climate-highlights-2024 ↵ E.U. Copernicus Marine Service ( 2024 ) Global Ocean Ensemble Physics Reanalysis . 26 June 2024. doi: 10.48670/moi-00024 OpenUrl CrossRef ↵ Fautin Dunn DF , Devaney DM , Roth B ( 1980 ) Stylobates: A shell-forming sea anemone (Coelenterata, Anthozoa , Actiniidae). Pacific Science , 34 , 379 – 388 . OpenUrl ↵ Felsenstein J ( 2004 ) Inferring phylogenies . Sinauer Associates, Sunderland , Mass . ↵ Fouqueau L , Roze D ( 2021 ) The evolution of sex along an environmental gradient . Evolution , 75 , 1334 – 1347 . doi: 10.1111/evo.14237 OpenUrl CrossRef PubMed ↵ Franklin S , Alpert P , Salguero-Gómez R , Janovský Z , Herben T , Klimešová J , Douhovnikoff V ( 2021 ) Next-gen plant clonal ecology . Perspectives in Plant Ecology, Evolution and Systematics , 49 , 125601 . doi: 10.1016/j.ppees.2021.125601 OpenUrl CrossRef ↵ Geller JB , Walton ED ( 2001 ) Breaking up and getting together: evolution of symbiosis and cloning by fission in sea anemones (genus Anthopleura ) . Evolution , 55 , 1781 – 1794 . doi: 10.1111/j.0014-3820.2001.tb00827.x OpenUrl CrossRef PubMed Web of Science ↵ Glon H , Daly M , Carlton JT , Flenniken MM , Currimjee Z ( 2020 ) Mediators of invasions in the sea: life history strategies and dispersal vectors facilitating global sea anemone introductions . Biological Invasions , 22 , 3195 – 3222 . doi: 10.1007/s10530-020-02321-6 OpenUrl CrossRef PubMed Grebelnyi SD Distribution of the sea anemones in the Arctic . In: The theoretical and practical significance of Coelenterata , pp. 20 – 33 . Publishing House of the Zoological Institute of the USSR Academy of Sciences , Leningrad . ↵ Halkett F , Simon J-C , Balloux F ( 2005 ) Tackling the population genetics of clonal and partially clonal organisms . Trends in Ecology & Evolution , 20 , 194 – 201 . doi: 10.1016/j.tree.2005.01.001 OpenUrl CrossRef PubMed Web of Science ↵ He Z-Z , Stotz GC , Liu X , Liu J-J , Wang Y-G , Yang J , Li L-F , Zhang W-J , Nan P , Song Z-P ( 2024 ) A global synthesis of the patterns of genetic diversity in endangered and invasive plants . Biological Conservation , 291 , 110473 . doi: 10.1016/j.biocon.2024.110473 OpenUrl CrossRef ↵ Jacquard A ( 1968 ) Panmixie et consanguinité. Quelques précisions de langage . Population (French Edition) , 23 , 1065 . doi: 10.2307/1528645 OpenUrl CrossRef ↵ Jacquard A ( 1970 ) Structure génétiques des populations . Masson & cie , Paris . ↵ Johannesson K , André C ( 2006 ) Invited review: Life on the margin: genetic isolation and diversity loss in a peripheral marine ecosystem, the Baltic Sea . Molecular Ecology , 15 , 2013 – 2029 . doi: 10.1111/j.1365-294X.2006.02919.x OpenUrl CrossRef PubMed Web of Science ↵ Jombart T , Devillard S , Balloux F ( 2010 ) Discriminant analysis of principal components: a new method for the analysis of genetically structured populations . BMC Genetics , 11 , 94 . doi: 10.1186/1471-2156-11-94 OpenUrl CrossRef PubMed ↵ Kaliszewicz A , Panteleeva N , Olejniczak I , Boniecki P , Sawicki M ( 2012 ) Internal brooding affects the spatial structure of intertidal sea anemones in the Arctic-boreal region . Polar Biology , 35 , 1911 – 1919 . doi: 10.1007/s00300-012-1232-1 OpenUrl CrossRef ↵ Kendrick GA , Waycott M , Carruthers TJB , Cambridge ML , Hovey R , Krauss SL , Lavery PS , Les DH , Lowe RJ , Vidal OM i, Ooi JLS , Orth RJ , Rivers DO , Ruiz-Montoya L , Sinclair EA , Statton J , van Dijk JK , Verduin JJ ( 2012 ) The central role of dispersal in the maintenance and persistence of seagrass populations . BioScience , 62 , 56 – 65 . doi: 10.1525/bio.2012.62.1.10 OpenUrl CrossRef Web of Science ↵ Klimeš L ( 2008 ) Clonal splitters and integrators in harsh environments of the Trans-Himalaya . Evolutionary Ecology , 22 , 351 – 367 . doi: 10.1007/s10682-007-9195-3 OpenUrl CrossRef ↵ Krueger-Hadfield SA , Guillemin M-L , Destombe C , Valero M , Stoeckel S ( 2021 ) Exploring the genetic consequences of clonality in haplodiplontic taxa . Journal of Heredity , 112 , 92 – 107 . doi: 10.1093/jhered/esaa063 OpenUrl CrossRef PubMed ↵ Krueger-Hadfield SA , Kübler JE , Dudgeon SR ( 2013 ) Reproductive effort of Mastocarpus papillatus (Rhodophyta) along the California coast . Journal of Phycology , 49 , 271 – 281 . doi: 10.1111/jpy.12034 OpenUrl CrossRef PubMed ↵ Laakkonen HM , Hardman M , Strelkov P , Väinölä R ( 2021 ) Cycles of trans-Arctic dispersal and vicariance, and diversification of the amphi-boreal marine fauna . Journal of Evolutionary Biology , 34 , 73 – 96 . doi: 10.1111/jeb.13674 OpenUrl CrossRef ↵ Lane SM , Wilson AJ , Briffa M ( 2020 ) Analysis of direct and indirect genetic effects in fighting sea anemones . Behavioral Ecology , 31 , 540 – 547 . doi: 10.1093/beheco/arz217 OpenUrl CrossRef PubMed ↵ Laughinghouse HDI , Müller KM , Adey WH , Lara Y , Young R , Johnson G ( 2015 ) Evolution of the northern rockweed, Fucus distichus , in a regime of glacial cycling: implications for benthic algal phylogenetics . PLOS ONE , 10 , e0143795 . doi: 10.1371/journal.pone.0143795 OpenUrl CrossRef PubMed ↵ Liu J , Dong M , Miao SL , Li ZY , Song MH , Wang RQ ( 2006 ) Invasive alien plants in China: role of clonality and geographical origin . Biological Invasions , 8 , 1461 – 1470 . doi: 10.1007/s10530-005-5838-x OpenUrl CrossRef Web of Science ↵ Locarnini RA , Mishonov AV , Baranova OK , Reagan JR , Boyer TP , Seidov D , Wang Z , Garcia HE , Bouchard C , Cross SL , Paver CR , Dukhovskoy D ( 2024 ) World Ocean Atlas 2023, Volume 1: Temperature . doi: 10.25923/54bh-1613 OpenUrl CrossRef ↵ Loseva LM ( 1974 ) Polovoye sostoyanie Bunodactis stella [Sex composition of Bunodactis stella ] . In: Biologiya razvitiya morskikh organizmov (Developmental Biology of Marine Organisms) , pp. 133 – 142 . Akademiya Nauk SSSR , Apatity . ↵ Manuel RL ( 1981 ) British Anthozoa: Keys and notes for the identification of the species . Acad. Pr , London u.a . ↵ McMahon KM , Evans RD , van Dijk K , Hernawan U , Kendrick GA , Lavery PS , Lowe R , Puotinen M , Waycott M ( 2017 ) disturbance is an important driver of clonal richness in tropical seagrasses . Frontiers in Plant Science , 8 . doi: 10.3389/fpls.2017.02026 OpenUrl CrossRef ↵ Melia N , Haines K , Hawkins E ( 2016 ) Sea ice decline and 21st century trans-Arctic shipping routes . Geophysical Research Letters , 43 , 9720 – 9728 . doi: 10.1002/2016GL069315 OpenUrl CrossRef Moran EV , Alexander JM ( 2014 ) Evolutionary responses to global change: lessons from invasive species . Ecology Letters , 17 , 637 – 649 . doi: 10.1111/ele.12262 OpenUrl CrossRef PubMed Navascués M , Stoeckel S , Mariette S ( 2010 ) Genetic diversity and fitness in small populations of partially asexual, self-incompatible plants . Heredity , 104 , 482 – 492 . doi: 10.1038/hdy.2009.159 OpenUrl CrossRef PubMed ↵ Orive ME , Krueger-Hadfield SA ( 2021 ) Sex and asex: a clonal lexicon . Journal of Heredity , 112 , 1 – 8 . doi: 10.1093/jhered/esaa058 OpenUrl CrossRef PubMed ↵ Ottaway JR ( 1977 ) Predators of sea anemones . Tuatara , 22 , 213 – 220 . OpenUrl Pannell JR , Auld JR , Brandvain Y , Burd M , Busch JW , Cheptou P-O , Conner JK , Goldberg EE , Grant A-G , Grossenbacher DL , Hovick SM , Igic B , Kalisz S , Petanidou T , Randle AM , de Casas RR , Pauw A , Vamosi JC , Winn AA ( 2015 ) The scope of Baker’s law . New Phytologist , 208 , 656 – 667 . doi: 10.1111/nph.13539 OpenUrl CrossRef PubMed ↵ Peck JR , Yearsley JM , Waxman D ( 1998 ) Explaining the geographic distributions of sexual and asexual populations . Nature , 391 , 889 – 892 . doi: 10.1038/36099 OpenUrl CrossRef ↵ Peniston JH , Barfield M , Holt RD , Orive ME ( 2021 ) Environmental fluctuations dampen the effects of clonal reproduction on evolutionary rescue . Journal of Evolutionary Biology , 34 , 710 – 722 . doi: 10.1111/jeb.13778 OpenUrl CrossRef PubMed ↵ Pereyra RT , Bergström L , Kautsky L , Johannesson K ( 2009 ) Rapid speciation in a newly opened postglacial marine environment, the Baltic Sea . BMC Evolutionary Biology , 9 , 70 . doi: 10.1186/1471-2148-9-70 OpenUrl CrossRef PubMed ↵ Pereyra RT , Rafajlović M , De Wit P , Pinder M , Kinnby A , Töpel M , Johannesson K ( 2023 ) Clones on the run: The genomics of a recently expanded partially clonal species . Molecular Ecology , 32 , 4209 – 4223 . doi: 10.1111/mec.16996 OpenUrl CrossRef ↵ Pierre J-S , Stoeckel S , Wajnberg E ( 2022 ) The advantage of sex: Reinserting fluctuating selection in the pluralist approach . PLOS ONE , 17 , e0272134 . doi: 10.1371/journal.pone.0272134 OpenUrl CrossRef ↵ PAME ( 2024 ) Arctic Ship Traffic Data (ASTD) . Protection of the Arctic Marine Environment, Working Group of the Arctic Council . https://pame.is/ourwork/arctic-shipping/current-shipping-projects/astd/ ↵ Rafajlović M , Kleinhans D , Gulliksson C , Fries J , Johansson D , Ardehed A , Sundqvist L , Pereyra RT , Mehlig B , Jonsson PR , Johannesson K ( 2017 ) Neutral processes forming large clones during colonization of new areas . Journal of Evolutionary Biology , 30 , 1544 – 1560 . doi: 10.1111/jeb.13124 OpenUrl CrossRef PubMed ↵ Reichel K , Masson J-P , Malrieu F , Arnaud-Haond S , Stoeckel S ( 2016 ) Rare sex or out of reach equilibrium? The dynamics of F IS in partially clonal organisms . BMC Genetics , 17 , 76 . doi: 10.1186/s12863-016-0388-z OpenUrl CrossRef PubMed ↵ Riemann-Zürneck K ( 1998 ) How sessile are sea anemones? a review of free-living forms in the Actiniaria cnidaria: Anthozoa . Marine Ecology , 19 , 247 – 261 . doi: 10.1111/j.1439-0485.1998.tb00466.x OpenUrl CrossRef ↵ Rodríguez E , Fautin D , Daly M ( 2025 ) World list of Actiniaria . doi: 10.14284/568 OpenUrl CrossRef Rodzen JA , May B ( 2002 ) Inheritance of microsatellite loci in the white sturgeon ( Acipenser transmontanus ) . Genome , 45 , 1064 – 1076 . doi: 10.1139/g02-083 OpenUrl CrossRef PubMed ↵ Roiloa SR ( 2019 ) Clonal traits and plant invasiveness: The case of Carpobrotus N . E.Br. (Aizoaceae). Perspectives in Plant Ecology, Evolution and Systematics , 40 , 125479 . doi: 10.1016/j.ppees.2019.125479 OpenUrl CrossRef ↵ Ronfort J , Jenczewski E , Bataillon T , Rousset F ( 1998 ) Analysis of Population Structure in Autotetraploid Species . Genetics , 150 , 921 – 930 . doi: 10.1093/genetics/150.2.921 OpenUrl Abstract / FREE Full Text ↵ Rudels B , Carmack E ( 2022 ) Arctic ocean water mass structure and circulation . Oceanography , 35 , 52 – 65 . doi: 10.5670/oceanog.2022.116 OpenUrl CrossRef ↵ Rudin FS , Briffa M ( 2012 ) Is boldness a resource-holding potential trait? Fighting prowess and changes in startle response in the sea anemone , Actinia equina. Proceedings of the Royal Society B: Biological Sciences , 279 , 1904 – 1910 . doi: 10.1098/rspb.2011.2418 OpenUrl CrossRef PubMed ↵ Russo CAM , Solé-Cava AM , Thorpe JP ( 1994 ) Population structure and genetic variation in two tropical sea anemones (Cnidaria , Actinidae) with different reproductive strategies. Marine Biology , 119 , 267 – 276 . doi: 10.1007/BF00349566 OpenUrl CrossRef ↵ Ryan WH , Aida J , Krueger-Hadfield SA ( 2021 ) The contribution of clonality to population genetic structure in the sea anemone, Diadumene lineata . Journal of Heredity , 112 , 122 – 139 . doi: 10.1093/jhered/esaa050 OpenUrl CrossRef PubMed ↵ Ryan WH , Miller TE ( 2019 ) Reproductive strategy changes across latitude in a clonal sea anemone . Marine Ecology Progress Series , 611 , 129 – 141 . doi: 10.3354/meps12862 OpenUrl CrossRef ↵ Sanamyan NP , Sanamyan KE ( 1998 ) Some Actiniaria from the Commander Islands (Cnidaria: Anthozoa) . Zoosystematica Rossica , 1 , 1 – 8 . OpenUrl ↵ Sanamyan NP , Sanamyan KE ( 2006 ) The genera Urticina and Cribrinopsis (Anthozoa: Actiniaria) from the north-western Pacific . Journal of Natural History , 40 , 359 – 393 . doi: 10.1080/00222930600703532 OpenUrl CrossRef ↵ Sanamyan NP , Sanamyan KE , McDaniel N , Martynov AV , Korshunova TA , Bocharova ES ( 2019 ) A revision of sea anemones of the genus Cribrinopsis Carlgren, 1921(Actiniaria: Actiniidae) from British Columbia with the description of a new species . Marine Biodiversity , 49 , 1951 – 1969 . doi: 10.1007/s12526-019-00956-w OpenUrl CrossRef ↵ Sanamyan NP , Sanamyan KE , Mercier A , Hamel J-F , Bocharova ES ( 2020 ) Morphological and molecular assessment of large sea anemones (Actiniaria: Actiniidae) in Newfoundland (eastern Canada) . Polar Biology , 43 , 495 – 509 . doi: 10.1007/s00300-020-02652-z OpenUrl CrossRef ↵ Sax DF , Stachowicz JJ , Brown JH , Bruno JF , Dawson MN , Gaines SD , Grosberg RK , Hastings A , Holt RD , Mayfield MM , O’Connor MI , Rice WR ( 2007 ) Ecological and evolutionary insights from species invasions . Trends in Ecology & Evolution , 22 , 465 – 471 . doi: 10.1016/j.tree.2007.06.009 OpenUrl CrossRef PubMed Web of Science ↵ Schön I , Martens K , Dijk P (Eds.) ( 2009 ) Lost Sex: The Evolutionary Biology of Parthenogenesis . Springer Netherlands , Dordrecht . doi: 10.1007/978-90-481-2770-2 OpenUrl CrossRef ↵ Sherman CDH , Ayre DJ ( 2008 ) Fine-scale adaptation in a clonal sea anemone . Evolution , 62 , 1373 – 1380 . doi: 10.1111/j.1558-5646.2008.00375.x OpenUrl CrossRef PubMed Web of Science ↵ Sherman CDH , Hunt A , Ayre DJ ( 2008 ) Is life history a barrier to dispersal? Contrasting patterns of genetic differentiation along an oceanographically complex coast . Biological Journal of the Linnean Society , 95 , 106 – 116 . doi: 10.1111/j.1095-8312.2008.01044.x OpenUrl CrossRef Web of Science ↵ Sherman CDH , Peucker AJ , Ayre DJ ( 2007 ) Do reproductive tactics vary with habitat heterogeneity in the intertidal sea anemone Actinia tenebrosa ? Journal of Experimental Marine Biology and Ecology , 340 , 259 – 267 . doi: 10.1016/j.jembe.2006.09.016 OpenUrl CrossRef ↵ Slatkin M , Takahata N ( 1985 ) The average frequency of private alleles in a partially isolated population . Theoretical Population Biology , 28 , 314 – 331 . doi: 10.1016/0040-5809(85)90032-2 OpenUrl CrossRef ↵ Stoeckel S , Arnaud-Haond S , Krueger-Hadfield SA ( 2021 ) The combined effect of haplodiplonty and partial clonality on genotypic and genetic diversity in a finite mutating population . Journal of Heredity , 112 , 78 – 91 . doi: 10.1093/jhered/esaa062 OpenUrl CrossRef PubMed ↵ Stoeckel S , Becheler R , Bocharova E , Barloy D ( 2024 ) GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms . Molecular Ecology Resources , 24 , e13886 . doi: 10.1111/1755-0998.13886 OpenUrl CrossRef ↵ Stoeckel S , Masson J-P ( 2014 ) The exact distributions of F IS under partial asexuality in small finite populations with mutation . PLOS ONE , 9 , e85228 . doi: 10.1371/journal.pone.0085228 OpenUrl CrossRef Stoeckel S , Porro B , Arnaud-Haond S ( 2021 ) The discernible and hidden effects of clonality on the genotypic and genetic states of populations: Improving our estimation of clonal rates . Molecular Ecology Resources , 21 , 1068 – 1084 . doi: 10.1111/1755-0998.13316 OpenUrl CrossRef PubMed ↵ Tibayrenc M , Avise JC , Ayala FJ ( 2015 ) In the light of evolution IX: Clonal reproduction: Alternatives to sex . Proceedings of the National Academy of Sciences , 112 , 8824 – 8826 . doi: 10.1073/pnas.1508087112 OpenUrl FREE Full Text ↵ Timmermans M-L , Marshall J ( 2020 ) understanding arctic ocean circulation: a review of ocean dynamics in a changing climate . Journal of Geophysical Research: Oceans , 125 , e2018J C 014378 . doi: 10.1029/2018JC014378 OpenUrl CrossRef ↵ Van Oppen MJH , Peplow LM , Kininmonth S , Berkelmans R ( 2011 ) Historical and contemporary factors shape the population genetic structure of the broadcast spawning coral, Acropora millepora , on the Great Barrier Reef . Molecular Ecology , 20 , 4899 – 4914 . doi: 10.1111/j.1365-294X.2011.05328.x OpenUrl CrossRef PubMed ↵ Veale A , Lavery S ( 2012 ) The population genetic structure of the waratah anemone ( Actinia tenebrosa ) around New Zealand . New Zealand Journal of Marine and Freshwater Research , 46 , 523 – 536 . doi: 10.1080/00288330.2012.730053 OpenUrl CrossRef ↵ Vermeij GJ , Roopnarine PD ( 2008 ) The coming Arctic invasion . Science , 321 , 780 – 781 . doi: 10.1126/science.1160852 OpenUrl Abstract / FREE Full Text ↵ Verrill AE ( 1864 ) List of the polyps and corals sent by the Museum of Comparative Zoology to other institutions in exchange, with annotations . Bulletin of the Museum of Comparative Zoology , 1 , 29 – 60 . OpenUrl ↵ Virtanen P , Gommers R , Oliphant TE , Haberland M , Reddy T , Cournapeau D , Burovski E , Peterson P , Weckesser W , Bright J , Van Der Walt SJ , Brett M , Wilson J , Millman KJ , Mayorov N , Nelson ARJ , Jones E , Kern R , Larson E , Carey CJ , Polat İ , Feng Y , Moore EW , VanderPlas J , Laxalde D , Perktold J , Cimrman R , Henriksen I , Quintero EA , Harris CR , Archibald AM , Ribeiro AH , Pedregosa F , Van Mulbregt P , SciPy 1.0 Contributors , Vijaykumar A , Bardelli AP , Rothberg A , Hilboll A , Kloeckner A , Scopatz A , Lee A , Rokem A , Woods CN , Fulton C , Masson C , Häggström C , Fitzgerald C , Nicholson DA , Hagen DR , Pasechnik DV , Olivetti E , Martin E , Wieser E , Silva F , Lenders F , Wilhelm F , Young G , Price GA , Ingold G-L , Allen GE , Lee GR , Audren H , Probst I , Dietrich JP , Silterra J , Webber JT , Slavič J , Nothman J , Buchner J , Kulick J , Schönberger JL , De Miranda Cardoso JV , Reimer J , Harrington J , Rodríguez JLC , Nunez-Iglesias J , Kuczynski J , Tritz K , Thoma M , Newville M , Kümmerer M , Bolingbroke M , Tartre M , Pak M , Smith NJ , Nowaczyk N , Shebanov N , Pavlyk O , Brodtkorb PA , Lee P , McGibbon RT , Feldbauer R , Lewis S , Tygier S , Sievert S , Vigna S , Peterson S , More S , Pudlik T , Oshima T , Pingel TJ , Robitaille TP , Spura T , Jones TR , Cera T , Leslie T , Zito T , Krauss T , Upadhyay U , Halchenko YO , Vázquez-Baeza Y ( 2020 ) SciPy 1.0: fundamental algorithms for scientific computing in Python . Nature Methods , 17 , 261 – 272 . doi: 10.1038/s41592-019-0686-2 OpenUrl CrossRef PubMed ↵ Waits LP , Luikart G , Taberlet P ( 2001 ) Estimating the probability of identity among genotypes in natural populations: cautions and guidelines . Molecular Ecology , 10 , 249 – 256 . doi: 10.1046/j.1365-294X.2001.01185.x OpenUrl CrossRef PubMed Web of Science ↵ Wang Y-J , Chen D , Yan R , Yu F-H , van Kleunen M ( 2019 ) Invasive alien clonal plants are competitively superior over co-occurring native clonal plants . Perspectives in Plant Ecology, Evolution and Systematics , 40 , 125484 . doi: 10.1016/j.ppees.2019.125484 OpenUrl CrossRef ↵ Wang Z , Xie L , Prather CM , Guo H , Han G , Ma C ( 2018 ) What drives the shift between sexual and clonal reproduction of Caragana stenophylla along a climatic aridity gradient? BMC Plant Biology , 18 , 91 . doi: 10.1186/s12870-018-1313-6 OpenUrl CrossRef PubMed ↵ Williams TM , Krueger-Hadfield SA , Hill-Spanik KM , Kosaki RK , Stoeckel S , Spalding HL ( 2024 ) The reproductive system of the cryptogenic alga Chondria tumulosa (Florideophyceae) at Manawai, Papahānaumokuākea Marine National Monument . Phycologia , 63 , 36 – 44 . doi: 10.1080/00318884.2023.2284074 OpenUrl CrossRef ↵ Yu F-H , Roiloa SR , Alpert P ( 2016 ) Editorial: Global Change, Clonal Growth, and Biological Invasions by Plants . Frontiers in Plant Science , 7 . doi: 10.3389/fpls.2016.01467 OpenUrl CrossRef ↵ Zhang R , Jia G , Diao X ( 2023 ) geneHapR: an R package for gene haplotypic statistics and visualization . BMC Bioinformatics , 24 , 199 . doi: 10.1186/s12859-023-05318-9 OpenUrl CrossRef PubMed View the discussion thread. Back to top Previous Next Posted November 13, 2025. Download PDF Data/Code Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean Ekaterina Bocharova , Aleksandr Volkov , Solenn Stoeckel bioRxiv 2025.04.24.650399; doi: https://doi.org/10.1101/2025.04.24.650399 Share This Article: Copy Citation Tools Increased clonality, decreased allele diversity and high genetic structure in tetraploid sea anemone Aulactinia stella populations from North Pacific to Atlantic across the Arctic Ocean Ekaterina Bocharova , Aleksandr Volkov , Solenn Stoeckel bioRxiv 2025.04.24.650399; doi: https://doi.org/10.1101/2025.04.24.650399 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Evolutionary Biology Subject Areas All Articles Animal Behavior and Cognition (7634) Biochemistry (17690) Bioengineering (13892) Bioinformatics (41935) Biophysics (21451) Cancer Biology (18587) Cell Biology (25499) Clinical Trials (138) Developmental Biology (13375) Ecology (19899) Epidemiology (2067) Evolutionary Biology (24316) Genetics (15609) Genomics (22505) Immunology (17736) Microbiology (40393) Molecular Biology (17180) Neuroscience (88598) Paleontology (666) Pathology (2832) Pharmacology and Toxicology (4824) Physiology (7641) Plant Biology (15151) Scientific Communication and Education (2045) Synthetic Biology (4294) Systems Biology (9824) Zoology (2271)