Ester A. Serrão

Standardizing methods to address clonality in population studies

Although clonal species are dominant in many habitats, from unicellular organisms to plants and animals, ecological and particularly evolutionary studies on clonal species have been strongly limited by the difficulty in assessing the number, size and longevity of genetic individuals within a population. The development of molecular markers has allowed progress in this area, and although allozymes remain of limited use due to their typically low level of polymorphism, more polymorphic markers have been discovered during the last decades, supplying powerful tools to overcome the problem of clonality assessment. However, population genetics studies on clonal organisms lack a standardized framework to assess clonality, and to adapt conventional data analyses to account for the potential bias due to the possible replication of the same individuals in the sampling. Moreover, existing studies used a variety of indices to describe clonal diversity and structure such that comparison among studies is difficult at best. We emphasize the need for standardizing studies on clonal organisms, and particularly on clonal plants, in order to clarify the way clonality is taken into account in sampling designs and data analysis, and to allow further comparison of results reported in distinct studies. In order to provide a first step towards a standardized framework to address clonality in population studies, we review, on the basis of a thorough revision of the literature on population structure of clonal plants and of a complementary revision on other clonal organisms, the indices and statistics used so far to estimate genotypic or clonal diversity and to describe clonal structure in plants. We examine their advantages and weaknesses as well as various conceptual issues associated with statistical analyses of population genetics data on clonal organisms. We do so by testing them on results from simulations, as well as on two empirical data sets of microsatellites of the seagrasses Posidonia oceanica and Cymodocea nodosa. Finally, we also propose a selection of new indices and methods to estimate clonal diversity and describe clonal structure in a way that should facilitate comparison between future studies on clonal plants, most of which may be of interest for clonal organisms in general.

North Atlantic phylogeography and large-scale population differentiation of the seagrass Zostera marina L.

As the most widespread seagrass in temperate waters of the Northern Hemisphere, Zostera marina provides a unique opportunity to investigate the extent to which the historical legacy of the last glacial maximum (LGM18 000–10 000 years bp) is detectable in modern population genetic structure. We used sequences from the nuclear rDNA–internal transcribed spacer (ITS) and chloroplast matK‐intron, and nine microsatellite loci to survey 49 populations (> 2000 individuals) from throughout the species’ range. Minimal sequence variation between Pacific and Atlantic populations combined with biogeographical groupings derived from the microsatellite data, suggest that the trans‐Arctic connection is currently open. The east Pacific and west Atlantic are more connected than either is to the east Atlantic. Allelic richness was almost two‐fold higher in the Pacific. Populations from putative Atlantic refugia now represent the southern edges of the distribution and are not genetically diverse. Unexpectedly, the highest allelic diversity was observed in the North Sea–Wadden Sea–southwest Baltic region. Except for the Mediterranean and Black Seas, significant isolation‐by‐distance was found from ~150 to 5000 km. A transition from weak to strong isolation‐by‐distance occurred at ~150 km among northern European populations suggesting this scale as the natural limit for dispersal within the metapopulation. Links between historical and contemporary processes are discussed in terms of the projected effects of climate change on coastal marine plants. The identification of a high genetic diversity hotspot in Northern Europe provides a basis for restoration decisions.

Genetic structure at range edge: low diversity and high inbreeding in Southeast Asian mangrove ( Avicennia marina ) populations

Understanding the genetic composition and mating systems of edge populations provides important insights into the environmental and demographic factors shaping species’ distribution ranges. We analysed samples of the mangrove Avicennia marina from Vietnam, northern Philippines and Australia, with microsatellite markers. We compared genetic diversity and structure in edge (Southeast Asia, and Southern Australia) and core (North and Eastern Australia) populations, and also compared our results with previously published data from core and southern edge populations. Comparisons highlighted significantly reduced gene diversity and higher genetic structure in both margins compared to core populations, which can be attributed to very low effective population size, pollinator scarcity and high environmental pressure at distribution margins. The estimated level of inbreeding was significantly higher in northeastern populations compared to core and southern populations. This suggests that despite the high genetic load usually associated with inbreeding, inbreeding or even selfing may be advantageous in margin habitats due to the possible advantages of reproductive assurance, or local adaptation. The very high level of genetic structure and inbreeding show that populations of A. marina are functioning as independent evolutionary units more than as components of a metapopulation system connected by gene flow. The combinations of those characteristics make these peripheral populations likely to develop local adaptations and therefore to be of particular interest for conservation strategies as well as for adaptation to possible future environmental changes.

Network analysis identifies weak and strong links in a metapopulation system

The identification of key populations shaping the structure and connectivity of metapopulation systems is a major challenge in population ecology. The use of molecular markers in the theoretical framework of population genetics has allowed great advances in this field, but the prime question of quantifying the role of each population in the system remains unresolved. Furthermore, the use and interpretation of classical methods are still bounded by the need for a priori information and underlying assumptions that are seldom respected in natural systems. Network theory was applied to map the genetic structure in a metapopulation system by using microsatellite data from populations of a threatened seagrass, Posidonia oceanica, across its whole geographical range. The network approach, free from a priori assumptions and from the usual underlying hypotheses required for the interpretation of classical analyses, allows both the straightforward characterization of hierarchical population structure and the detection of populations acting as hubs critical for relaying gene flow or sustaining the metapopulation system. This development opens perspectives in ecology and evolution in general, particularly in areas such as conservation biology and epidemiology, where targeting specific populations is crucial.

Ecological genetics in the North Atlantic: environmental gradients and adaptation at specific loci.

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

Within‐population spatial genetic structure, neighbourhood size and clonal subrange in the seagrass Cymodocea nodosa

The extent of clonality within populations strongly influences their spatial genetic structure (SGS), yet this is hardly ever thoroughly analysed. We employed spatial autocorrelation analysis to study effects of sexual and clonal reproduction on dispersal of the dioecious seagrass Cymodocea nodosa. Analyses were performed both at genet level (i.e. excluding clonal repeats) and at ramet level. Clonal structure was characterized by the clonal subrange, a spatial measure of the linear limits where clonality still affects SGS. We show that the clonal subrange is equivalent to the distance where the probability of clonal identity approaches zero. This combined approach was applied to two meadows with different levels of disturbance, Cadiz (stable) and Alfacs (disturbed). Genotypic richness, the proportion of the sample representing distinct genotypes, was moderate (0.38 Cadiz, 0.46 Alfacs) mostly due to dominance of a few clones. Expected heterozygosities were comparable to those found in other clonal plants. SGS analyses at the genet level revealed extremely restricted gene dispersal in Cadiz (Sp = 0.052, a statistic reflecting the decrease of pairwise kinship with distance), the strongest SGS found for seagrass species, comparable only to values for selfing herbaceous land plants. At Cadiz the clonal subrange extended across shorter distances (20–25 m) than in Alfacs (30–35 m). Comparisons of sexual and vegetative components of gene dispersal suggest that, as a dispersal vector within meadows, clonal spread is at least as important as sexual reproduction. The restricted dispersal and SGS pattern in both meadows indicates that the species follows a repeated seedling recruitment strategy.

Implications of Extreme Life Span in Clonal Organisms: Millenary Clones in Meadows of the Threatened Seagrass Posidonia oceanica

The maximum size and age that clonal organisms can reach remains poorly known, although we do know that the largest natural clones can extend over hundreds or thousands of metres and potentially live for centuries. We made a review of findings to date, which reveal that the maximum clone age and size estimates reported in the literature are typically limited by the scale of sampling, and may grossly underestimate the maximum age and size of clonal organisms. A case study presented here shows the occurrence of clones of slow-growing marine angiosperm Posidonia oceanica at spatial scales ranging from metres to hundreds of kilometres, using microsatellites on 1544 sampling units from a total of 40 locations across the Mediterranean Sea. This analysis revealed the presence, with a prevalence of 3.5 to 8.9%, of very large clones spreading over one to several (up to 15) kilometres at the different locations. Using estimates from field studies and models of the clonal growth of P. oceanica, we estimated these large clones to be hundreds to thousands of years old, suggesting the evolution of general purpose genotypes with large phenotypic plasticity in this species. These results, obtained combining genetics, demography and model-based calculations, question present knowledge and understanding of the spreading capacity and life span of plant clones. These findings call for further research on these life history traits associated with clonality, considering their possible ecological and evolutionary implications.

Population genetics of dwarf eelgrass Zostera noltii throughout its biogeographic range

The marine angiosperm Zostera noltii (dwarf eelgrass), an important facilitator species and food source for invertebrates and waterfowl, predominantly inhabits intertidal habitats along eastern Atlantic shores from Mauritania to southern Norway/Kattegat Sea and throughout the Mediterranean, Black and Azov seas. We used 9 microsatellite loci to characterize population struc- ture at a variety of spatial scales among 33 populations from 11 localities throughout the entire biogeographic range. Isolation by distance analysis suggested a panmictic genetic neighborhood of 100 to 150 km. At the global scale, a neighbor-joining tree based on Reynolds distances revealed strongly-supported groups corresponding to northern Europe, Mauritania and the Black/Azov Sea; separate Mediterranean and Atlantic-Iberian groups were poorly supported. Clones (genets with multiple ramets) were present in most populations but were generally small (ca. <3 m 2 ). Exceptions were found in Mauritania (ca. 29 m in length), the Azov Sea (ca. 40 m in length) and the Black Sea (ca. 50 m in length). Although genetic diversity and allelic richness generally decreased from Mauri- tania to Denmark, the putative post-glacial recolonization route, both were unexpectedly high among populations from the German Wadden Sea.

Intriguing asexual life in marginal populations of the brown seaweed Fucus vesiculosus

Reproduction of attached large brown algae is known to occur only by sexual zygotes. Using microsatellites we show evolution of asexual reproduction in the bladder wrack promoting population persistence in the brackish water Baltic Sea (< 6 psu). Here a dwarf morph of Fucus vesiculosus is dominated by a single clone but clonal reproduction is also present in the common form of the species. We describe a possible mechanism for vegetative reproduction of attached algae, and conclude that clonality plays an important role in persistence and dispersal of these marginal populations, in which sexual reproduction is impaired by low salinity.

Genetic structure in the Mediterranean seagrass Posidonia oceanica: Disentangling past vicariance events from contemporary patterns of gene flow

The Mediterranean Sea is a two‐basin system, with the boundary zone restricted to the Strait of Sicily and the narrow Strait of Messina. Two main population groups are recognized in the Mediterranean endemic seagrass Posidonia oceanica, corresponding to the Western and the Eastern basins. To address the nature of the East–West cleavage in P. oceanica, the main aims of this study were: (i) to define the genetic structure within the potential contact zone (i.e. the Strait of Sicily) and clarify the extent of gene flow between the two population groups, and (ii) to investigate the role of present water circulation patterns vs. past evolutionary events on the observed genetic pattern. To achieve these goals, we utilized SSR markers and we simulated, with respect to current regime, the possible present‐day dispersal pattern of Posidonia floating fruits using 28‐day numerical Lagrangian trajectories. The results obtained confirm the presence of the two main population groups, without any indices of reproductive isolation, with the break zone located at the level of the Southern tip of Calabria. The populations in the Strait of Sicily showed higher affinity with Western than with Eastern populations. This pattern of genetic structure probably reflects historical avenues of recolonization from relict glacial areas and past vicariance events, but seems to persist as a result of the low connectivity among populations via marine currents, as suggested by our dispersal simulation analysis.