Nelson C. Coelho

Isolation by oceanographic distance explains genetic structure for Macrocystis pyrifera in the Santa Barbara Channel.

Ocean currents are expected to be the predominant environmental factor influencing the dispersal of planktonic larvae or spores; yet, their characterization as predictors of marine connectivity has been hindered by a lack of understanding of how best to use oceanographic data. We used a high‐resolution oceanographic model output and Lagrangian particle simulations to derive oceanographic distances (hereafter called transport times) between sites studied for Macrocystis pyrifera genetic differentiation. We build upon the classical isolation‐by‐distance regression model by asking how much additional variability in genetic differentiation is explained when adding transport time as predictor. We explored the extent to which gene flow is dependent upon seasonal changes in ocean circulation. Because oceanographic transport between two sites is inherently asymmetric, we also compare the explanatory power of models using the minimum or the mean transport times. Finally, we compare the direction of connectivity as estimated by the oceanographic model and genetic assignment tests. We show that the minimum transport time had higher explanatory power than the mean transport time, revealing the importance of considering asymmetry in ocean currents when modelling gene flow. Genetic assignment tests were much less effective in determining asymmetry in gene flow. Summer‐derived transport times, in particular for the month of June, which had the strongest current speed, greatest asymmetry and highest spore production, resulted in the best‐fit model explaining twice the variability in genetic differentiation relative to models that use geographic distance or habitat continuity. The best overall model also included habitat continuity and explained 65% of the variation in genetic differentiation among sites.

Habitat continuity and geographic distance predict population genetic differentiation in giant kelp

Isolation by distance (IBD) models are widely used to predict levels of genetic connectivity as a function of Euclidean distance, and although recent studies have used GIS-landscape ecological approaches to improve the predictability of spatial genetic structure, few if any have addressed the effect of habitat continuity on gene flow. Landscape effects on genetic connectivity are even less understood in marine populations, where habitat mapping is particularly challenging. In this study, we model spatial genetic structure of a habitat-structuring species, the giant kelp Macrocystis pyrifera, using highly variable microsatellite markers. GIS mapping was used to characterize habitat continuity and distance between sampling sites along the mainland coast of the Santa Barbara Channel, and their roles as predictors of genetic differentiation were evaluated. Mean dispersal distance (sigma) and effective population size (Ne) were estimated by comparing our IBD slope with those from simulations incorporating habitat continuity and spore dispersal characteristics of the study area. We found an allelic richness of 7-50 alleles/locus, which to our knowledge is the highest reported for macroalgae. The best regression model relating genetic distance to habitat variables included both geographic distance and habitat continuity, which were respectively, positively and negatively related to genetic distance. Our results provide strong support for a dependence of gene flow on both distance and habitat continuity and elucidate the combination of Ne and a that explained genetic differentiation.

High and distinct range-edge genetic diversity despite local bottlenecks.

The genetic consequences of living on the edge of distributional ranges have been the subject of a largely unresolved debate. Populations occurring along persistent low latitude ranges (rear-edge) are expected to retain high and unique genetic diversity. In contrast, currently less favourable environmental conditions limiting population size at such range-edges may have caused genetic erosion that prevails over past historical effects, with potential consequences on reducing future adaptive capacity. The present study provides an empirical test of whether population declines towards a peripheral range might be reflected on decreasing diversity and increasing population isolation and differentiation. We compare population genetic differentiation and diversity with trends in abundance along a latitudinal gradient towards the peripheral distribution range of Saccorhiza polyschides , a large brown seaweed that is the main structural species of kelp forests in SW Europe. Signatures of recent bottleneck events were also evaluated to determine whether the recently recorded distributional shifts had a negative influence on effective population size. Our findings show decreasing population density and increasing spatial fragmentation and local extinctions towards the southern edge. Genetic data revealed two well supported groups with a central contact zone. As predicted, higher differentiation and signs of bottlenecks were found at the southern edge region. However, a decrease in genetic diversity associated with this pattern was not verified. Surprisingly, genetic diversity increased towards the edge despite bottlenecks and much lower densities, suggesting that extinctions and recolonizations have not strongly reduced diversity or that diversity might have been even higher there in the past, a process of shifting genetic baselines.

Genes Left Behind: Climate Change Threatens Cryptic Genetic Diversity in the Canopy-Forming Seaweed Bifurcaria bifurcata

The global redistribution of biodiversity will intensify in the coming decades of climate change, making projections of species range shifts and of associated genetic losses important components of conservation planning. Highly-structured marine species, notably brown seaweeds, often harbor unique genetic variation at warmer low-latitude rear edges and thus are of particular concern. Here, a combination of Ecological Niche Models (ENMs) and molecular data is used to forecast the potential near-future impacts of climate change for a warm-temperate, canopy forming seaweed, Bifurcaria bifurcata. ENMs for B. bifurcata were developed using marine and terrestrial climatic variables, and its range projected for 2040-50 and 2090-2100 under two greenhouse emission scenarios. Geographical patterns of genetic diversity were assessed by screening 18 populations spawning the entire distribution for two organelle genes and 6 microsatellite markers. The southern limit of B. bifurcata was predicted to shift northwards to central Morocco by the mid-century. By 2090-2100, depending on the emission scenario, it could either retreat further north to western Iberia or be relocated back to Western Sahara. At the opposing margin, B. bifurcata was predicted to expand its range to Scotland or even Norway. Microsatellite diversity and endemism were highest in Morocco, where a unique and very restricted lineage was also identified. Our results imply that B. bifurcata will maintain a relatively broad latitudinal distribution. Although its persistence is not threatened, the predicted extirpation of a unique southern lineage or even the entire Moroccan diversity hotspot will erase a rich evolutionary legacy and shrink global diversity to current (low) European levels. NW Africa and similarly understudied southern regions should receive added attention if expected range changes and diversity loss of warm-temperate species is not to occur unnoticed.

Reproductive strategies and isolation-by-demography in a marine clonal plant along an eutrophication gradient.

Genetic diversity in clonal organisms includes two distinct components, (i) the diversity of genotypes or clones (i.e. genotypic richness) in a population and (ii) that of the alleles (i.e. allelic and gene diversity within populations, and differentiation between populations). We investigated how population differentiation and genotypic components are associated across a gradient of eutrophication in a clonal marine plant. To that end, we combined direct measurements of sexual allocation (i.e. flower and seed counts) and genotypic analyses, which are used as an estimator of effective sexual reproduction across multiple generations. Genetic differentiation across sites was also modelled according to a hypothesis here defined as isolation‐by‐demography, in which we use population‐specific factors, genotypic richness and eutrophication that are hypothesized to affect the source‐sink dynamics and thus influence the genetic differentiation between a pair of populations. Eutrophic populations exhibited lower genotypic richness, in agreement with lower direct measurements of sexual allocation and contemporaneous gene flow. Genetic differentiation, while not explained by distance, was best predicted by genotypic richness and habitat quality. A multiple regression model using these two predictors was considered the best model (R2 = 0.43). In this study, the relationship between environment and effective sexual–asexual balance is not simply (linearly) predicted by direct measurements of sexual allocation. Our results indicate that population‐specific factors and the isolation‐by‐demography model should be used more often to understand genetic differentiation.

Looking into the black box: simulating the role of self-fertilization and mortality in the genetic structure of Macrocystis pyrifera

Patterns of spatial genetic structure (SGS), typically estimated by genotyping adults, integrate migration over multiple generations and measure the effective gene flow of populations. SGS results can be compared with direct ecological studies of dispersal or mating system to gain additional insights. When mismatches occur, simulations can be used to illuminate the causes of these mismatches. Here, we report a SGS and simulation‐based study of self‐fertilization in Macrocystis pyrifera, the giant kelp. We found that SGS is weaker than expected in M. pyrifera and used computer simulations to identify selfing and early mortality rates for which the individual heterozygosity distribution fits that of the observed data. Only one (of three) population showed both elevated kinship in the smallest distance class and a significant negative slope between kinship and geographical distance. All simulations had poor fit to the observed data unless mortality due to inbreeding depression was imposed. This mortality could only be imposed for selfing, as these were the only simulations to show an excess of homozygous individuals relative to the observed data. Thus, the expected data consistently achieved nonsignificant differences from the observed data only under models of selfing with mortality, with best fits between 32% and 42% selfing. Inbreeding depression ranged from 0.70 to 0.73. The results suggest that density‐dependent mortality of early life stages is a significant force in structuring Macrocystis populations, with few highly homozygous individuals surviving. The success of these results should help to validate simulation approaches even in data‐poor systems, as a means to estimate otherwise difficult‐to‐measure life cycle parameters.

Characterization of 12 polymorphic microsatellite markers in the sugar kelp Saccharina latissima

Saccharina latissima is an ecologically and economically important kelp species native to the coastal regions of the Northern Hemisphere. This species has considerable phylogeographic structure and morphological plasticity, but lack of resolution of available genetic markers prevents a finer characterization of its genetic diversity and structure. Here, we describe 12 microsatellite loci identified in silico in a genomic library, and assess their polymorphism in three distant populations. Allelic richness at the species level was relatively high (5–23 alleles per locus), as was gene diversity within populations (0.42 < HE < 0.62). In addition, individuals readily form distinct genotypic clusters matching their populations of origin. The variation detected confirms the great potential of these markers to investigate the biogeography and population dynamics of S. latissima, and to better characterize its genetic resources for the establishing farming industry.

Characterization of ten highly polymorphic microsatellite loci for the intertidal mussel Perna perna, and cross species amplification within the genus

BackgroundThe brown mussel Perna perna (Linnaeus, 1758) is a dominant constituent of intertidal communities and a strong invader with multiple non-native populations distributed around the world. In a previous study, two polymorphic microsatellite loci were developed and used to determine population-level genetic diversity in invasive and native P. perna populations. However, higher number of microsatellite markers are required for reliable population genetic studies.In this context, in order to understand P. perna origins and history of invasion and to compare population genetic structure in native versus invaded areas, we developed 10 polymorphic microsatellite markers.FindingsDescribed microsatellite markers were developed from an enriched genomic library. Analyses and characterization of loci using 20 individuals from a population in Western Sahara revealed on average 11 alleles per locus (range: 5–27) and mean gene diversity of 0.75 (range: 0.31 - 0.95). One primer pair revealed possible linkage disequilibrium while heterozygote deficiency was significant at four loci. Six of these markers cross-amplified in P. canaliculus (origin: New Zealand).ConclusionsDeveloped markers will be useful in addressing a variety of questions concerning P. perna, including dispersal scales, genetic variation and population structure, in both native and invaded areas.

Microsatellite markers for the Arctic copepod Calanus glacialis and cross-amplification with C. finmarchicus

Calanus glacialis is a major component of Arctic zooplankton and a keystone species in Arctic marine ecosystems. Due to the observed climate warming, its numbers are being reduced to the advantage of a sibling Atlantic species Calanus finmarchicus. We developed and characterized the first set of microsatellite markers in this species to investigate its population genetic structure and dispersal capabilities. Nine polymorphic loci displayed an average of 7.3 alleles (range between 2 and 13) and the levels of expected heterozygosity ranged from 0.039 to 0.806. These provide a valuable tool to understand present connectivity patterns across Arctic regions, look for signatures of past climate effects and predict the response to future climate-driven environmental changes. Additionally, due to the cross-amplification with C. finmarchicus, the markers can be used to discriminate between these sibling species.

Cryptic diversity, geographical endemism and allopolyploidy in NE Pacific seaweeds

BackgroundMolecular markers are revealing a much more diverse and evolutionarily complex picture of marine biodiversity than previously anticipated. Cryptic and/or endemic marine species are continually being found throughout the world oceans, predominantly in inconspicuous tropical groups but also in larger, canopy-forming taxa from well studied temperate regions. Interspecific hybridization has also been found to be prevalent in many marine groups, for instance within dense congeneric assemblages, with introgressive gene-flow being the most common outcome. Here, using a congeneric phylogeographic approach, we investigated two monotypic and geographically complementary sister genera of north-east Pacific intertidal seaweeds (Hesperophycus and Pelvetiopsis), for which preliminary molecular tests revealed unexpected conflicts consistent with unrecognized cryptic diversity and hybridization.ResultsThe three recovered mtDNA clades did not match a priori species delimitations. H. californicus was congruent, whereas widespread P. limitata encompassed two additional narrow-endemic species from California - P. arborescens (here genetically confirmed) and P. hybrida sp. nov. The congruence between the genotypic clusters and the mtDNA clades was absolute. Fixed heterozygosity was apparent in a high proportion of loci in P. limitata and P. hybrida, with genetic analyses showing that the latter was composed of both H. californicus and P. arborescens genomes. All four inferred species could be distinguished based on their general morphology.ConclusionsThis study confirmed additional diversity and reticulation within NE Pacific Hesperophycus/Pelvetiopsis, including the validity of the much endangered, modern climatic relict P. arborescens, and the identification of a new, stable allopolyploid species (P. hybrida) with clearly discernable ancestry (♀ H. californicus x ♂ P. arborescens), morphology, and geographical distribution. Allopolyploid speciation is otherwise completely unknown in brown seaweeds, and its unique occurrence within this genus (P. limitata possibly representing a second example) remains enigmatic. The taxonomic separation of Hesperophycus and Pelvetiopsis is not supported and the genera should be synonymized; we retain only the latter. The transitional coastline between Point Conception and Monterey Bay represented a diversity hotspot for the genus and the likely sites of extraordinary evolutionary events of allopolyploid speciation at sympatric range contact zones. This study pinpoints how much diversity (and evolutionary processes) potentially remains undiscovered even on a conspicuous seaweed genus from the well-studied Californian intertidal shores let alone in other, less studied marine groups and regions/depths.