Population genetics at three spatial scales of a rare sponge living in fragmented habitats

  • Andrea Blanquer1Email author and

    Affiliated with

    • Maria J Uriz1

      Affiliated with

      BMC Evolutionary Biology201010:13

      DOI: 10.1186/1471-2148-10-13

      Received: 15 July 2009

      Accepted: 14 January 2010

      Published: 14 January 2010

      Abstract

      Background

      Rare species have seldom been studied in marine habitats, mainly because it is difficult to formally assess the status of rare species, especially in patchy benthic organisms, for which samplings are often assumed to be incomplete and, thus, inappropriate for establishing the real abundance of the species. However, many marine benthic invertebrates can be considered rare, due to the fragmentation and rarity of suitable habitats. Consequently, studies on the genetic connectivity of rare species in fragmented habitats are basic for assessing their risk of extinction, especially in the context of increased habitat fragmentation by human activities. Sponges are suitable models for studying the intra- and inter-population genetic variation of rare invertebrates, as they produce lecitotrophic larvae and are often found in fragmented habitats.

      Results

      We investigated the genetic structure of a Mediterranean sponge, Scopalina lophyropoda (Schmidt), using the allelic size variation of seven specific microsatellite loci. The species can be classified as "rare" because of its strict habitat requirements, the low number of individuals per population, and the relatively small size of its distribution range. It also presents a strong patchy distribution, philopatric larval dispersal, and both sexual and asexual reproduction. Classical genetic-variance-based methods (AMOVA) and differentiation statistics revealed that the genetic diversity of S. lophyropoda was structured at the three spatial scales studied: within populations, between populations of a geographic region, and between isolated geographic regions, although some stochastic gene flow might occur among populations within a region. The genetic structure followed an isolation-by-distance pattern according to the Mantel test. However, despite philopatric larval dispersal and fission events in the species, no single population showed inbreeding, and the contribution of clonality to the population makeup was minor (only ca. 4%).

      Conclusions

      The structure of the S. lophyropoda populations at all spatial scales examined confirms the philopatric larval dispersal that has been reported. Asexual reproduction does not seem to play a relevant role in the populations. The heterozygote excess and the lack of inbreeding could be interpreted as a hitherto unknown outcrossing strategy of the species. The envisaged causes for this strategy are sperm dispersal, a strong selection against the mating of genetically related individuals to avoid inbreeding depression or high longevity of genets combined with stochastic recruitment events by larvae from other populations. It should be investigated whether this strategy could also explain the genetic diversity of many other patchy marine invertebrates whose populations remain healthy over time, despite their apparent rarity.

      Background

      Many marine invertebrates are sessile as adults, have limited dispersal abilities, and their populations are in general well-structured (e.g., [18]), despite the fact that marine environments have long been considered to be more open than terrestrial environments [9]. Genetic exchange (i.e., gene flow) among populations of sessile invertebrates is restricted to the dispersal of gametes, and asexual and/or sexual propagula (often lecitotrophic, short-living, small larvae).

      The particular rocky habitats of most marine sessile invertebrates (i.e., sponges, cnidarians, brooding bryozoans and ascidians) are often fragmented in space. Habitat fragmentation has been reported to account for population isolation in invertebrates with poor dispersal capabilities through the hampering of genetic exchange among populations [10]. Moreover, at long temporal scales, fragmented habitats can promote species rarity [11]. Consequently, fragmented habitats may either put the survival of the species at risk or trigger genetic, biological and/or ecological strategies to overcome the problems that arise from isolation. At evolutionary time-scales, fragmentation of habitats can also favor speciation (e.g., [12, 13]).

      Rare species are considered to be strongly vulnerable in fragmented habitats [14], especially when gene flow among patches is restricted or absent. Most definitions of the rarity of species include ideas about habitat requirements, the number of individuals, and the size of distribution ranges, which can vary over temporal and spatial scales [15]. The status of "rare" species is especially difficult to assess appropriately for marine invertebrates, as there is a lack of quantitative data in these respects [11, 14]. Accordingly, many cryptic, patchily-distributed marine invertebrates are not considered rare because of assumed inappropriate sampling and the incorrect assumption that marine distant populations are always connected through the exchange of pelagically dispersing propagula [16, 17]. Since habitat fragmentation can be aggravated by human perturbations, studies on the genetic connectivity of rare species in fragmented habitats are particularly necessary for the assessment of their extinction risk, and this should be taken into account when developing management rules or for conservation issues.

      Sponges are key, widespread elements in benthic ecosystems (e.g., [18, 19]). Thus, they represent suitable models for studying inter-population connectivity, the colonization of new habitats, and intra- and inter-population genetic variation in sessile marine invertebrates. However, despite their ecological importance, very little is known about the genetic issues of sponge populations such as gene flow, genetic variation, and effective dispersal. Therefore, genetic studies on sponges are still largely lacking.

      Scopalina lophyropoda (Schmidt, 1862) has only been recorded from a few Western Mediterranean localities [2024]. Recent reports from Southern and Eastern Mediterranean localities [25, 26] need genetic confirmation, because cryptic species have been recently reported [13, 27]. The species only dwells on the lower zones of north-facing rocky walls and blocks (i.e., close to the bottom), in areas with a certain amount of suspended organic material (e.g., close to a river mouth) [13, 28]. The species habitat is strongly fragmented along the species distribution area. After performing exhaustive samplings at distances shorter than 100 km, with sub-samplings at each sampling site, we can state that S. lophyropoda can be classed as a "rare species" according to the definition of Chapman [11], owing to its strict habitat requirements, the low number of individuals per population, and the relatively small size of its distribution range: the well-known Western Mediterranean sublittoral [13].

      Several known reproduction traits of S. lophyropoda allowed us to predict some aspects of its population genetics. First, the larvae do not move far away from the release point, but remain close to their parents in a vertical posture, suggesting philopatry [22, 28, 29]. Second, asexual reproduction by fission [30] indicates that clonality may play some role in determining the genetic makeup of populations. Third, the release of asexual propagula containing several mature larvae [31] may represent a stochastic dispersal event. A previous study on the species at a small spatial scale (from centimeters to tens of meters) using molecular markers confirmed the restricted dispersal of philopatric larvae, but revealed that the extent of clonality was low [32].

      The main objective of this study was to gain knowledge on the genetic connectivity among populations of a marine rare sponge species inhabiting fragmented habitats (S. lophyropoda) in order to better understand the maintenance and growth of its populations, which have been monitored for the last 20 years. The study was performed at three spatial-scales: within populations, between populations of a geographic region, and between isolated regions. We also assessed the relative extent of sexual and asexual reproduction in the populations.

      The knowledge acquired on this sponge species may be of interest in studies dealing with those rare modular invertebrates with limited dispersal, and whose habitats are disrupted, which are frequent in benthic ecosystems [11].

      Results

      Microsatellite polymorphism across loci

      All of the amplifications were visualized on agarose gels and produced positive bands. However, despite good amplifications, double checking, and repetitions, the sizes of the alleles could not be estimated properly from the electropherograms for five individuals (ca. 2%) for locus Scol_d, three individuals (ca. 1%) for Scol_n and one (ca. 0.5%) for both, Scol_p and Scol_c. These unsized alleles were not included in the analyses, but cannot be considered null alleles. All microsatellite loci were polymorphic (Table 1). A total of 46 alleles were detected over the seven loci, with a range from 4 (Scol_c, Scol_i and Scol_p) to 12 (Scol_l) alleles per locus. There were private alleles in all populations except FE (probably due to the low number of individuals in this population; Table 2). A third of all of the linkage disequilibrium tests revealed linkages between pairs of loci, but these linkages were not consistent across populations. In each population, all loci behaved as neutral according to the Ewens-Watterson test. All loci, but Scol_l and Scol_n, showed negative values for the inbreeding coefficient (FIS; Table 1).
      Table 1

      Summary of statistics for the Scopalina lophyropoda microsatellites

      Locus

      A

      HO

      HE

      FIS

      Scol_c

      4

      0.4932

      0.4417

      -0.173*

      Scol_d

      4

      0.7834

      0.808

      -0.076

      Scol_i

      8

      0.1486

      0.1411

      -0.071

      Scol_l

      12

      0.6036

      0.6578

      0.013

      Scol_n

      8

      0.4155

      0.5698

      0.112

      Scol_p

      4

      0.1991

      0.1862

      -0.098*

      Scol_r

      6

      0.3378

      0.3187

      -0.169*

      A = the total number of alleles per locus; the observed (HO) and expected (HE) heterozygosities and the inbreeding coefficients (FIS; *P ≤ 0.05) for each locus.

      Genetic diversity at different spatial scales

      All the "a priori" defined populations of Scopalina lophyropoda, were genetically diverse (Table 2), with mean gene diversities (HE) ranging from 0.322 (GA) to 0.469 (SA). Single locus exact tests for the Hardy-Weinberg equilibrium indicated significant deviations (P < 0.05) in 9 out of 42 comparisons. In the MU population there was heterozygote deficit, whereas in all other populations the departures were due to heterozygote excesses (Table 2). Multi-locus fixation indexes, FIS, ranged from 0.017 in MU (the only positive value) to -0.099 in SA, and both were in the Blanes region. Multi-locus FIS values for all populations, but that in MU, were negative, although only values for CG, SA and GA were significantly different from 0 (P < 0.05), revealing outcrossing.
      Table 2

      Summary of statistics for the Scopalina lophyropoda microsatellites for each region and population site

      Region

      Location

      N

      Na

      Pa

      HO

      HE

      FIS

      CG

       

      32

      5.1

      5

      0.509

      0.467

      -0.088*

       

      CG

      32

      5.1

      5

      0.509

      0.467

      -0.088*

      BL

       

      97

      5

      2

      0.469

      0.465

      -0.023

       

      MU

      40

      4.4

      1

      0.443

      0.450

      0.017

       

      SA

      43

      4.3

      1

      0.515

      0.469

      -0.099*

       

      FE

      14

      3

      0

      0.432

      0.406

      -0.069

        

      93

           

      PC

        

      4.1

      3

      0.343

      0.344

      -0.016

       

      BA

      40

      3

      1

      0.375

      0.368

      -0.018

       

      GA

      53

      3.9

      2

      0.346

      0.322

      -0.075*

      N = the sample size; Na = the mean number of alleles per locus; Pa = the number of private alleles; the mean observed (HO) and expected heterozygosities (HE), and the inbreeding coefficients (FIS; *P ≤ 0.05).

      The global estimates of FST and RST, over all populations, were significantly (P < 0.05) different from 0 (FST = 0.122 and RST = 0.076). Multi-locus pairwise FST indicated significant differences among all the population pairs (Table 3A). Regression of FST/(1-FST) values against the logarithm of the geographic distances revealed significant isolation by distance according to the Mantel test (R2 = 0.4579; P < 0.05, Figure 1).
      Table 3

      Pairwise multi-locus FST values between Scopalina lophyropoda population pairs (A) and between geographical region pairs (B) (* P < 0.05).

       

      CG

      MU

      SA

      FE

      BA

      GA

       

      CG

      BL

      PC

      CG

      -

           

      CG

      -

        

      Mu

      0.195*

      -

          

      BL

      0.150*

      -

       

      SA

      0.127*

      0.021*

      -

         

      PC

      0.179*

      0.058*

      -

      FE

      0.144*

      0.094*

      0.066*

      -

            

      BA

      0.141*

      0.094*

      0.050*

      0.051*

      -

       

      B

         

      GA

      0.216*

      0.105*

      0.100*

      0.053*

      0.060*

      -

          

      A

                
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-10-13/MediaObjects/12862_2009_Article_1229_Fig1_HTML.jpg
      Figure 1

      Genetic isolation by distance of all theScopalina lophyropodapopulations as inferred using multi-locus estimates of F ST /(1-F ST ) and the logarithm of the geographical distance (the Mantel test). The correlation coefficient R2 = 0.45, p < 0.05.

      As for the geographical regions (Table 2), Cabo de Gata Natural Park had the highest value of mean number of alleles (5.1) whereas Port Cros National Park had the lowest value (4.1). Mean expected heterozygosities were similar for Cabo de Gata and Blanes regions (0.467 vs 0.465), but were slightly lower for Port Cros (0.344). Multi-locus pairwise FST indicated significant differences among the geographical regions (Table 3B). The global FST and RST values among the regions (FST = 0.104; RST = 0.042) revealed the presence of significant genetic structure (P < 0.05).

      Overall, the hierarchical AMOVA revealed significant genetic structuring of the analyzed samples at the three geographical levels (7.57% among regions, 4.66% among populations within regions; and 87.77% within populations; Table 4). The analyses revealed that most of the genetic variation was found within the populations.
      Table 4

      The results of the analysis of molecular variance (AMOVA) on Scopalina lophyropoda microsatellites.

      Source of variation

      d. f.

      Sum of squares

      Variance components

      Percentage of variation

      Fixation indices

      Pvalue

      Among regions

      2

      46.955

      0.1207

      7.57

      0.076 FCT

      P < 0.001

      Among populations within regions

      3

      19.827

      0.0744

      4.66

      0.050 FSC

      P < 0.001

      Within populations

      438

      613.166

      139.992

      87.77

      0.123 FST

      P < 0.05

      Total

      433

      679.948

      159.499

         
      Bayesian clustering analyses detected the highest likelihood for the model with three genetically homogeneous groups of individuals (K), as shown by the highest peak in the ΔK for K = 3 (Figure 2), according to Evanno et al. (2005).
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-10-13/MediaObjects/12862_2009_Article_1229_Fig2_HTML.jpg
      Figure 2

      Values of ΔKfor each number of genetically homogeneous populations (K) considering all theScopalina lophyropodaindividuals.

      The Wilcoxon test detected recent bottlenecks in the CG, MU and GA populations (P two-tails values of 0.016, 0.0078 and 0.039, respectively), under the Stepwise Mutation Model (SMM). The tests performed, assuming the Two Phase model (TPM) and the Infinite Allele model (IAM), rendered non-significant results. With respect to the geographical regions, the Wilcoxon test again detected recent bottlenecks in each of the three regions (CG, BL, PC) using the SMM model (P two-tails values of 0.016, 0.023 and 0,039, respectively) but not using the TPM and IAM models.

      Contribution of clonal reproduction to the population structure

      The number of multi locus genotypes (MLGs) occurring more than once was variable in the populations studied (Table 5). Individuals sharing MLGs were found in all populations except FE and SA. In Port Cros National Park, one MLG was shared by individuals from different populations: two from the BA population and one from the GA population (Table 5).
      Table 5

      Identical multi-locus genotypes (MLGs) that were detected in the Scopalina lophyropoda populations.

          

      Genotype

       

      Region

      Pop.

      MLG

      N

      Scol_c

      Scol_d

      Scol_i

      Scol_l

      Scol_n

      Scol_p

      Scol_r

      MLGsim

      CG

      CG

                
        

      1

      2

      146/148

      313/319

      153/153

      165/169

      199/201

      149/149

      118/124

      *

        

      2

      2

      144/148

      317/319

      153/153

      165/169

      199/199

      149/149

      124/124

      *

      Bl

                 
       

      MU

      1

      2

      146/148

      317/317

      153/153

      169/169

      193/193

      149/149

      124/126

      ns

       

      SA

      0

               
       

      FE

      0

               

      PC

                 
       

      BA

      1

      2

      148/148

      313/313

      153/153

      171/171

      193/193

      149/149

      124/124

      ns

        

      2

      2

      148/148

      313/319

      153/153

      169/169

      193/193

      149/149

      124/124

      ns

        

      3

      2

      146/148

      313/319

      153/153

      169/171

      193/215

      149/149

      124/124

      ns

        

      4

      2

      148/148

      315/319

      153/153

      169/171

      193/193

      149/149

      124/124

      ns

        

      5

      2

      146/148

      313/319

      153/153

      169/171

      193/199

      149/149

      124/124

      ns

        

      6

      3

      148/148

      319/319

      153/153

      171/171

      193/193

      149/149

      124/124

      *

        

      7

      4

      146/148

      313/319

      153/153

      169/169

      193/199

      149/149

      124/124

      *

       

      GA

      1

      2

      148/148

      311/313

      153/153

      169/169

      193/193

      149/149

      124/124

      ns

        

      2

      2

      148/148

      311/313

      153/153

      169/171

      199/215

      149/149

      124/124

      ns

        

      3

      2

      148/148

      315/315

      153/153

      171/171

      193/193

      149/149

      120/124

      *

        

      4

      2

      148/148

      313/313

      152/153

      169/169

      193/217

      149/149

      124/124

      *

        

      5

      3

      148/148

      311/313

      153/153

      171/171

      193/193

      149/149

      124/124

      ns

       

      BA+GA

      1

      2+1

      146/148

      313/319

      153/153

      169/171

      193/193

      149/149

      124/124

      ns

      The genotypes are characterized using the allele sizes; N, the number of individuals sharing MLG; MLGsim, (* indicates significantly small probabilities of being the result of random mating; and ns, non-significance at P < 0.05). Numbers typed in bold indicate the allele with the highest frequency in the population.

      The probabilities of finding identical MLGs by random mating were significantly small (P < 0.05) for the ramets from CG, for the groups of three and four ramets sharing MLGs in BA, and for two groups of two ramets presenting identical MLGs in GA (Table 5). A sexual origin for the two individuals in BA and the one in GA that shared MLG could not be rejected (P < 0.05). The probability of finding the latter MLG in three different genets was estimated by considering the allele frequencies in the whole Port Cros National Park region (the BA and GA populations).

      Discussion

      Microsatellites have proven their use as powerful markers for population and clonality studies in the sponge Scopalina lophyropoda. The differentiation statistics used revealed that the genetic diversity was spatially structured in this rare species at the three spatial scales analyzed (as revealed by the FST values and the AMOVA). Most diversity was found within the populations, but differentiation among populations and among geographic regions was also significant. The Bayesian clustering method (STRUCTURE), however, turned out to be more conservative and only detected three genetically independent groups, which correspond to the three geographical regions. However, when the polymorphism shown by the selected markers is limited, as it occurs in sponge microsatellites, or there is a strong isolation by distance pattern, the STRUCTURE algorithm is limited to resolve population clusters [33], and has been reported to provide a lower resolution than differentiation statistics [3436]. Thus, considering all other results (Mantel test, AMOVA, FST values, and the presence of private alleles), we can conclude that the sampled populations are genetically well structured, though some gene flow may take place within regions.

      The positive regression between genetic differentiation and geographic distance between populations of this species revealed distance isolation and short-range effective dispersal. This pattern agrees with the philopatric behavior reported for S. lophyropoda larvae [28, 29], which, in combination with the strict habitat requirement of the species [13], could contribute to its extremely fragmented distribution. Environmental stress is likely to be stronger in small fragmented populations than in large populations [37], and it has been reported to lead to rapid genetic divergence [38]. The later studies have suggested that the intensity and direction of the selection could differ among small populations [37] and lead to speciation. The presence of several cryptic Scopalina species in the distribution area of S. lophyropoda [27] could be interpreted in that way.

      Gene flow can be maintained by dispersal of sexual and/or asexual propagules. In S. lophyropoda, the contribution of clonality to the populations was minor. Only 4% (9 out of 222) of the ramets analyzed can be assumed to be the result of asexual reproduction. A plausible explanation is that fissions are balanced with fusions between clone-mates [30], since clones remain attached to the rocky substrate and can contact each other when growing.

      Since asexual dispersal appears to be unimportant in the species and the few clones identified use to be found at smaller distances than 0.5 m [32], gene flow between the populations depends on the dispersal of sexual products: either larvae, or sperm. Unfortunately, little is known about the dispersal of sponge sperm in general, and nothing is known about S. lophyropoda sperm in particular but the reported short-range dispersal [22, 28, 29, 32] of S. lophyropoda larvae, fully agrees with the strong genetic structure in the species. The strong genetic structuring (AMOVA, FST values and Mantel test results) indicates that there is a current absence of connectivity among the regions but STRUCTURE results point to some, although restricted, gene flow among the populations within a region, which may occur through sporadic stochastic dispersal events. A combined sexual-asexual dispersal mechanism consisting in asexual fragments of ripe individuals containing larvae has been speculated [31], which, if it occurs in the field, could enhance gene flow over longer distances than those suggested by the observed larval behavior in the field.

      Bottleneck tests under IAM and TPM did not reveal any recent reduction in size of the populations, which is in agreement with the growth detected during the last 20 years in the several populations studied (authors' obs.). Only under the less probable mutation model for the microsatellites (SMM) [3942], the bottleneck analysis indicated recent reductions in size in some populations. Although we do not have information on the mutation model of S. lophyropoda microsatellites, some studies in other taxa (e.g. humans, fish, insects) reported that microsatellites fit the mixed TPM or the IAM, but not the SMM [3942]. Thus, field monitoring, the best models for microsatellite mutation, and the presence of private alleles in all populations (rare alleles are the first that disappear after a bottleneck), indicate that if any size reduction has occurred, it was not recently.

      Most marine invertebrates with non-feeding and short dispersed larvae present genetically structured populations and inbreeding, whatever the Phylum they belong to [2, 4345]. As for sponges, only one study of population genetics using microsatellites is available up to now (Crambe crambe; [2]). Some population traits of S. lophyropoda are shared with those reported for C. crambe, whereas other aspects are contrasting. The populations are genetically structured in both species, despite their respective patchy vs. continuous distributions and the differences in larval behavior [28]. Restricted dispersal seems to be the main factor responsible for the isolation by distance that is found in the two sponges. However, the contribution of clonality to the population makeup was rather insignificant in S. lophyropoda (4% on average) while it reached 23% in C. crambe [2, 28]. There are two possible explanations for the differences found between both species. First, in the study on C. crambe, several individuals did not amplify for each locus, which may contribute to a failure in differentiation of some individuals that were not real clones. On the other hand, fusion rates between ramets (also much more frequent in S. lophyropoda than in C. crambe) are more balanced in the former with fission rates and thus could account for the low extent of clonality in our "snapshot" sampling.

      Moreover, all C. crambe populations showed inbreeding [2], as expected in long-lived, low-dispersing organisms due to temporary Wahlund effects [46, 47]. Furthermore, individuals in small populations usually present inbreeding, which can increase their extinction probability by decreasing their potential for adaptation [39], especially when fragmented [14]. However, contrary to most predictions, we show here that S. lophyropoda has no inbreeding depression and exhibits heterozygote excess for most of the loci studied. In fact, our results contrast with theory and experimental evidences, which suggest that the effective size of small populations, and its response to selection, decreases with time [48]. The unexpected genetic results in the S. lophyropoda populations, however, support studies that predict that inbreeding should not be that common in marine invertebrates, despite their restricted dispersal, strong population structure, and consanguinity (e.g.,[4951]). Our study provides indirect evidence on the existence of genetic strategies for maintaining the high genetic diversity even when the populations are isolated by distance and there is restricted gene flow among the populations. Selection against homozygote and/or for outcrossing, extreme genet longevity combined with stochastic recruitment events from other populations [51] are among the strategies that may allow S. lophyropoda and other rare, marine, benthic invertebrates to persist in a fragmented scenario.

      Conclusions

      Scopalina lophyropoda, as other sessile marine invertebrates, a priori considered rare, lives without apparent constrains in fragmented populations with natural low abundances and restricted dispersal abilities. Mechanisms that are as yet poorly understood may be acting in those populations to maintain the necessary genetic diversity for the species to adapt to changing environments. This study has shown the results of these potential strategies, which would act in S. lophyropoda to avoid inbreeding, despite the species biological features. Efforts should be concentrated on the study of these mechanisms, which are responsible for the persistence of most marine sessile invertebrates with small, scattered populations and which have restricted dispersal, for their implications in the maintenance of marine biodiversity. Results on these mechanisms should be considered in the implementation of conservation directives, as rare species inhabiting fragmented habitats may require particular protection rules.

      Methods

      Sampling

      Scopalina lophyropoda populations were sampled by SCUBA diving from several localities in the Northwestern Mediterranean. The species identity was confirmed through amplification of a fragment of mtDNA cytochrome c oxidase gen, given the presence of several cryptic species [13, 27]. Due to its rarity, after an exhaustive search, only 8 populations were found and sampled along the western Mediterranean. Samplings were performed every 100 km from the Adriatic to the South of the Iberian Peninsula and at the Atlantic Canary Islands. Populations from Banyuls-sur-Mer and Marseilles that yielded only 4 and 8 individuals, respectively, were genetically characterized, but excluded from the data analyses to avoid artifacts from the use of extremely unbalanced data. Consequently, the study included fragments from 222 individuals collected at 6 locations belonging to 3 distinct geographic regions (Figure 3A): Cabo de Gata Natural Park, Spain (CG), 36°51'N, 2°10'W; Port Cros National Park, France (PC, Figure 3B), 43°01'N, 6°23'E (Ile de Bagaud -BA- and Ilot de la Gabinière -GA-); and Blanes, Catalonia (BL, Figure 3C), 41°40'N, 2°47'E, (Muntanyeta -MU-, Santa Anna Niells -SA-, and Fenals -FE-); Table 2). Cabo de Gata Natural Park and Port Cros National Park are marine protected areas (MPA) and are therefore under special management conditions. Sponge fragments were preserved in absolute ethanol, which was changed three times, before storing at -20°C until DNA extraction.
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-10-13/MediaObjects/12862_2009_Article_1229_Fig3_HTML.jpg
      Figure 3

      Location of theScopalina lophyropodapopulations. (A). Sampling sites at the Port Cros Park region (Ile de Bagaud, BA and Ilot de la Gabinière, GA; B), and at the Blanes region (Muntanyeta, MU, Santa Anna Niells, SA and Fenals, FE; C).

      Microsatellite genotyping

      The samples were cleaned from foreign tissues. The DNA was extracted using the DNeasy Tissue Kit (Qiagen) and screened for the seven microsatellite loci that had been previously isolated from S. lophyropoda [52], under the amplifying conditions described in that paper. Alleles were sized on an ABI 3700 automated sequencer (Scientific and Technical Services of the University of Barcelona) using GENESCAN™ and GENOTYPER™ software for the genotyping.

      Identical multi-locus genotypes

      Individuals sharing multi-locus genotypes (MLGs) could have resulted either from sampling individuals resulting from asexual reproduction (ramets), or from sampling two or more individuals resulting from sexual reproduction (genets). The probability of the latter depends on the allele frequencies in the population and the observed heterozygosities. To test this probability, we used the program MLGsim (10,000 simulations, P < 0.05; [53]). This program, using a simulation approach, calculates significance values for the likelihood that a MLG observed more than once in a population is the result of sexual reproduction.

      Population data analyses

      The TRANSFORMER-3 program [54] was used to generate the input files for the software packages used. Within each population, the analyses were performed at the genet level (i.e., excluding those individuals that presented significantly low probability (P < 0.05) of having a sexual origin). Allele frequencies, the observed and mean number of alleles per locus [55], and the observed and expected heterozygosities were estimated using POPGENE [56]. Wright's [57] inbreeding coefficients, FIS, as in [58] and departures from the Hardy-Weinberg equilibrium (HWE), using Bonferroni corrections, were estimated using FSTAT version 2.9.3 [59]. Ewens-Watterson tests of neutrality (10,000 simulated samples) and linkage disequilibria among all pairs of loci were also tested with POPGENE [56] and GENEPOP [60], respectively.

      Differentiation statistics based on microsatellite allele data are widely used for studying moderately-structured populations [61]. Wright's [57] FST values (as in [58]) were estimated using ARLEQUIN [62] and Slatkin's [63] global RST value (as in [64] Michalakis & Excoffier 1996) was estimated using GENEPOP [60]. Isolation by distance was tested using the correlation between pairwise FST/(1-FST) values and the logarithm of the geographical distances between the populations (as calculated from the minimum distance by sea), using a Mantel test. The analysis of molecular variance (AMOVA) is a testing procedure based on permutation analysis for studying molecular variation within species. It was performed using ARLEQUIN software [62]. The number of genetically homogeneous groups (K) was inferred using a Bayesian algorithm (STRUCTURE) [65]. K was inferred using the ad hoc statistic ΔK as in Evanno et al. [66]. The probability of recent effective population size reductions from allele data frequencies was tested using BOTTLENECK software [67]. This program is based on the assumption that populations, which have experienced a recent reduction of their effective population size, exhibit both a reduction of the allelic diversity and heterozygosity, although the allele numbers are reduced faster than the heterozygosity. The software computes for each population sample and for each locus, the distribution of the heterozygosity expected from the number of alleles, given the sample size under the assumption of mutation-drift equilibrium.

      Declarations

      Acknowledgements

      Authors thank the authorities in Cabo de Gata Natural Park and Port Cros National Park. We are also grateful to T. Perez and C. Lejeusne from Centre d'Océanologie de Marseille for the logistics and field assistance in Port Cros National Park; and to M. Pascual from University of Barcelona for laboratory facilities. This research was supported by an I3P predoctoral fellowship from CSIC to AB and funded by project CMT2007-66635-C02-01 of the Spanish Government (CICYT) to MJU.

      Authors’ Affiliations

      (1)
      Departament d'Ecologia Marina, Centre d'Estudis Avançats de Blanes, CSIC

      References

      1. Duran S, Giribet G, Turon X: Phylogeographical history of the sponge Crambe crambe (Porifera, Poecilosclerida): range expansion and recent invasion of the Macaronesian islands from the Mediterranean Sea. Mol Ecol 2004,13(1):109–122.PubMedView Article
      2. Duran S, Pascual M, Estoup A, Turon X: Strong population structure in the marine sponge Crambe crambe (Poecilosclerida) as revealed by microsatellite markers. Mol Ecol 2004,13(3):511–522.PubMedView Article
      3. Gutierrez-Rodriguez C, Lasker HR: Microsatellite variation reveals high levels of genetic variability and population structure in the gorgonian coral Pseudopterogorgia elisabethae across the Bahamas. Mol Ecol 2004,13(8):2211–2221.PubMedView Article
      4. Le Goff-Vitry MC, Pybus OG, Rogers AD: Genetic structure of the deep-sea coral Lophelia pertusa in the northeast Atlantic revealed by microsatellites and internal transcribed spacer sequences. Mol Ecol 2004,13(3):537–549.PubMedView Article
      5. Baums IB, Miller MW, Hellberg ME: Regionally isolated populations of an imperiled Caribbean coral, Acropora palmata. Mol Ecol 2005,14(5):1377–1390.PubMedView Article
      6. Lopez-Legentil S, Turon X, Planes S: Genetic structure of the star sea squirt, Botryllus schlosseri, introduced in southern European harbours. Mol Ecol 2006,15(13):3957–3967.PubMedView Article
      7. Vollmer SV, Palumbi SR: Restricted gene flow in the Caribbean staghorn coral Acropora cervicomis: Implications for the recovery of endangered reefs. J Hered 2007,98(1):40–50.PubMedView Article
      8. Pemberton AJ, Hansson LJ, Craig SF, Hughes RN, Bishop JDD: Microscale genetic differentiation in a sessile invertebrate with cloned larvae: investigating the role of polyembryony. Mar Biol 2007,153(1):71–82.View Article
      9. Caley MJ, Buckley KA, Jones GP: Separating ecological effects of habitat fragmentation, degradation, and loss on coral commensals. Ecology 2001,82(12):3435–3448.View Article
      10. Fahrig L: Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol Syst 2003, 34:487–515.View Article
      11. Chapman MG: Are there adequate data to assess how well theories of rarity apply to marine invertebrates? Biodivers Conserv 1999,8(10):1295–1318.View Article
      12. Lejeusne C, Chevaldonne P: Brooding crustaceans in a highly fragmented habitat: the genetic structure of Mediterranean marine cave-dwelling mysid populations. Mol Ecol 2006,15(13):4123–4140.PubMedView Article
      13. Blanquer A, Uriz MJ: 'A posteriori' searching for phenotypic characters to describe new cryptic species of sponges revealed by molecular markers (Dictyonellidae: Scopalina). Invertebr Syst 2008,22(5):489–502.View Article
      14. Sanderson W: Rarity of marine benthic species in Great Britain: development and application of assessment criteria. Aquat Conserv 1996,6(4):245–256.View Article
      15. Gaston K: Rarity. London: Chapman & Hall; 1994.
      16. Underwood AJ, Fairweather PG: Supply-side ecology and benthic marine assemblages. Trends Ecol Evol 1989,4(1):16–20.PubMedView Article
      17. Ward RD, Woodwark M, Skibinski DOF: A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes. J Fish Biol 1994,44(2):213–232.View Article
      18. Diaz M, Soest Rv: Sponges in time and space. Rotterdam: Balkema; 1994.
      19. Wulff J: Assessing and monitoring coral reef sponges: Why and how? Bullf Mar Sci 2001,69(2):831–846.
      20. Topsent E: Éponges observees dans le parages de Monaco. 1er Partie. Bull Inst Océanogr Monaco 1934, 650:1–42.
      21. Boury-Esnault N: Spongiaires de la zone rocheuse de Banyuls-sur-mer. II. Systématique. Vie et Milieu ser B 1971, 22:287–350.
      22. Uriz M: Morfología y comportamiento de la larva parenquímula de Scopalina lophyropoda Schmidt 1862 (demospongia, Halichondrida) y la formación del rhagon. Investigación Pesquera 1982, 42:313–322.
      23. Martinez A: Estudio sistemático de los poriferos del sureste ibérico. Murcia: University of Murcia; 1991.
      24. Cocito S, Bedulli D, Sgorbini S: Distribution patterns of the sublittoral epibenthic assemblages on a rocky shoal in the Ligurian Sea (NW Mediterranean). Sci Mar 2002,66(2):175–181.View Article
      25. Mustapha K, Zarrouk S, Souissi A, Abed A: Diversite des demosponges tunisiennes. Bull Inst Natl Sci Tech Salammbó 2003, 30:55–77.
      26. Voultsiadou E: Demosponge distribution in the eastern Mediterranean: a NW-SE gradient. Helgoland Mar Res 2005,59(3):237–251.View Article
      27. Blanquer A, Uriz MJ: Cryptic speciation in marine sponges evidenced by mitochondrial and nuclear genes: A phylogenetic approach. Mol Phylogenet Evol 2007,45(1):392–397.PubMedView Article
      28. Uriz MJ, Maldonado M, Turon X, Marti R: How do reproductive output, larval behaviour, and recruitment contribute to adult spatial patterns in Mediterranean encrusting sponges? Mar Ecol Prog Ser 1998, 167:137–148.View Article
      29. Mariani S, Uriz MJ, Turon X, Alcoverro T: Dispersal strategies in sponge larvae: integrating the life history of larvae and the hydrologic component. Oecologia 2006,149(1):174–184.PubMedView Article
      30. Blanquer A, Uriz MJ, Agell G: Hidden diversity in sympatric sponges: adjusting life-history dynamics to share substrate. Mar Ecol Prog Ser 2008, 371:109–115.View Article
      31. Maldonado M, Uriz MJ: Sexual propagation by sponge fragments. Nature 1999,398(6727):476–476.View Article
      32. Blanquer A, Uriz MJ, Caujapé-Castells J: Small-scale spatial genetic structure in Scopalina lophyropoda , an encrusting sponge with philopatric larval dispersal and frequent fission and fusion events. Mar Ecol Prog Ser 2009, 380:95–102.View Article
      33. Pritchard J, Wen W: Documentation for STRUCTURE software. Version 2. [http://​pritch.​bsd.​uchicago.​edu] 2003.
      34. Sequeira F, Alexandrino J, Weiss S, Ferrand N: Documenting the advantages and limitations of different classes of molecular markers in a well-established phylogeographic context: lessns from the Iberian endemic Golden-striped salamander, Chioglossa lusitanica (Caudata: Salamandridae). Biological Journal of the Linnean Society 2008, 95:371–387.View Article
      35. Worley K, Strobeck C, Arthur S, Carey J, Schwantje H, Veitch A, Coltman DW: Population genetic structure of North American thinhorn sheep (Ovis dalli). Molecular Ecology 2004,13(9):2545–2556.PubMedView Article
      36. Zulliger D, Tanner S, Ruch M, Ribi G: Genetic structure of the high dispersal Atlanto-Mediterranean sea star astropecten aranciacus revealed by mitochondrial DNA sequences and microsatellite loci. Mar Biol 2009, 156:597–610.View Article
      37. Willi Y, Van Buskirk J, Hoffmann AA: Limits to the adaptive potential of small populations. Annu Rev Ecol Evol Syst 2006, 37:433–458.View Article
      38. Antonovics J: The nature of limts to natural selection. Ann Mo Bot Gard 1976,63(2):224–247.View Article
      39. Di Rienzo A, Peterson A, Garza J, Valdes A, Slatkin M, Freimer N: Mutational processes of simple-sequence repeat loci in human populations. Proceedings of the National Academy of Sciences of the USA 1994, 91:3166–3170.PubMedView Article
      40. Estoup A, Cornuet JM, Goldstein D, Schlötterer C: Microsatellite evolution: inferences from population data. In Microsatellites: Evolution and Applications. Amsterdam: NATO Press; 1999:55–86.
      41. Ellegren H: Microsatellite mutations in the germline: implications for evolutionary inference. Trends in Genetics 2000, 16:551–558.PubMedView Article
      42. Schlötterer C: Evolutionary dynamic of microsatellite DNA. Chromosoma 2000, 109:365–371.PubMedView Article
      43. Goffredo S, Mezzomonaco L, Zaccanti F: Genetic differentiation among populations of the Mediterranean hermaphroditic coral Balanophyllia europaea (Scleratinia: Dendrophyllidae). Mar Biol 2004, 145:1075–1083.View Article
      44. Constantini F, Fauvelot C, Abbiati M: Fine-scale enetic structuring in Corallium rubrum : evidence of inbreeding and limited effective larval dispersal. Mar Ecol Prog Ser 2007, 340:109–119.View Article
      45. Pérez-Portela R, Turon X: Cryptic divergence and strong population structure in the colnial invertebrate Pycnoclavella communis (Ascidiacea) inferred from molecular data. Zoology 2008, 111:163–178.PubMedView Article
      46. Knowlton N, Jackson J, Thornhill N: Inbreeding and outbreeding in marine invertebrates. In The natural history of inbreeding and outbreeding: theoretical and empirical perspectives. Chicago: University of Chicago Press; 1993:200–249.
      47. Ayre DJ, Hughes TP: Genotypic diversity and gene flow in brooding and spawning corals along the Great Barrier Reef, Australia. Evolution 2000,54(5):1590–1605.PubMed
      48. Swindell WR, Bouzat JL: Modeling the adaptive potential of isolated populations: Experimental simulations using Drosophila. Evolution 2005,59(10):2159–2169.PubMed
      49. Ayre DJ, Miller KJ: Where do clonal coral larvae go? Adult genotypic diversity conflicts with reproductive effort in the brooding coral Pocillopora damicornis. Mar Ecol Prog Ser 2004, 277:95–105.View Article
      50. Addison JA, Hart MW: Spawning, copulation and inbreeding coefficients in marine invertebrates. Biol Letters 2005,1(4):450–453.View Article
      51. McFadden C: Contributions of sexual and asexual reproduction to population structure in the clonal soft coral Alcyonium rudyi . Evolution 1997,51(1):112–126.View Article
      52. Blanquer A, Uriz MJ, Pascual M: Polymorphic microsatellite loci isolated from the marine sponge Scopalina lophyropoda (Demospongiae: Halichondrida). Mol Ecol Notes 2005,5(3):466–468.View Article
      53. Stenberg P, Lundmark M, Saura A: MLGsim: a program for detecting clones using a simulation approach. Mol Ecol Notes 2003,3(2):329–331.View Article
      54. Caujapé-Castells J, Baccarani-Rosas M: Transformer-3: a program for the analysis of molecular population genetic data. [http://​exegen.​org/​es/​software.​php] EXEGEN software & Jardín Botánico Canario "Viera y Clavijo"; 2005.
      55. Kimura M, Crow J: The number of alleles that can be maintained in a finite population. Genetics 1964, 49:725–738.PubMed
      56. Yeh F, Yang R, Boyle T, Ye Z, Mao J: POPGENE, the user-friendly shareware for population genetic analysis. [http://​www.​ualberta.​ca/​~fyeh/​] Molecular Biology and Biotechnology Centre University of Alberta. Canada; 1997.
      57. Wright S: The genetical structure of populations. Ann Eugenics 1951, 15:323–354.
      58. Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution 1984,38(6):1358–1370.View Article
      59. Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered 1995,86(6):485–486.
      60. Raymond M, Rousset F: GENEPOP (Version1.2) - Population genetics software for exact tests and ecumenicism. J Hered 1995,86(3):248–249.
      61. Balloux F, Lugon-Moulin N: The estimation of population differentiation with microsatellite markers. Mol Ecol 2002,11(2):155–165.PubMedView Article
      62. Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 2005, 1:47–50.PubMed
      63. Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics 1995,139(1):457–462.PubMed
      64. Michalakis Y, Excoffier L: A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 1996,142(3):1061–1064.PubMed
      65. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics 2000,155(2):945–959.PubMed
      66. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 2005,14(8):2611–2620.PubMedView Article
      67. Cornuet JM, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 1996,144(4):2001–2014.PubMed

      Copyright

      © Blanquer and Uriz. 2010

      This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

      Advertisement