Pushed waves, trailing edges, and extreme events: Eco-evolutionary dynamics of a geographic range shift in the owl limpet, Lottia gigantea

As climatic variation re-shapes global biodiversity, understanding eco-evolutionary feedbacks during species range shifts is of increasing importance. Theory on range expansions distinguishes between two different forms: "pulled" and "pushed" waves. Pulled waves occur when the source of the expansion comes from low-density peripheral populations, while pushed waves occur when recruitment to the expanding edge is supplied by high-density populations closer to the species' core. How extreme events shape pushed/pulled wave expansion events, as well as trailing-edge declines/contractions, remains largely unexplored. We examined eco-evolutionary responses of a marine invertebrate (the owl limpet, Lottia gigantea) that increased in abundance during the 2014-2016 marine heatwaves near the poleward edge of its geographic range in the northeastern Pacific. We used whole-genome sequencing from 19 populations across >11 degrees of latitude to characterize genomic variation, gene flow, and demographic histories across the species' range. We estimated present-day dispersal potential and past climatic stability to identify how contemporary and historical seascape features shape genomic characteristics. Consistent with expectations of a pushed wave, we found little genomic differentiation between core and leading-edge populations, and higher genomic diversity at range edges. A large and well-mixed population in the northern edge of the species' range is likely a result of ocean current anomalies increasing larval settlement and high-dispersal potential across biogeographic boundaries. Trailing-edge populations have higher differentiation from core populations, possibly driven by local selection and limited gene flow, as well as high genomic diversity likely as a result of climatic stability during the Last Glacial Maximum. Our findings suggest that extreme events can drive poleward range expansions that carry the adaptive potential of core populations, while also cautioning that trailing-edge extirpations may threaten unique evolutionary variation. This work highlights the importance of understanding how both trailing and leading edges respond to global change and extreme events.

Genomic responses to parallel temperature gradients in the eelgrass Zostera marina in adjacent bays

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

Species' attributes predict the relative magnitude of ecological and genetic recovery following mass mortality

Theoretically, species' characteristics should allow estimation of dispersal potential and, in turn, explain levels of population genetic differentiation. However, a mismatch between traits and genetic patterns is often reported for marine species, and interpreted as evidence that life-history traits do not influence dispersal. Here, we couple ecological and genomic methods to test the hypothesis that species with attributes favouring greater dispersal potential-e.g., longer pelagic duration, higher fecundity and larger population size-have greater realized dispersal overall. We used a natural experiment created by a large-scale and multispecies mortality event which created a "clean slate" on which to study recruitment dynamics, thus simplifying a usually complex problem. We surveyed four species of differing dispersal potential to quantify the abundance and distribution of recruits and to genetically assign these recruits to probable parental sources. Species with higher dispersal potential recolonized a broader extent of the impacted range, did so more quickly and recovered more genetic diversity than species with lower dispersal potential. Moreover, populations of taxa with higher dispersal potential exhibited more immigration (71%-92% of recruits) than taxa with lower dispersal potential (17%-44% of recruits). By linking ecological with genomic perspectives, we demonstrate that a suite of interacting life-history and demographic attributes do influence species' realized dispersal and genetic neighbourhoods. To better understand species' resilience and recovery in this time of global change, integrative eco-evolutionary approaches are needed to more rigorously evaluate the effect of dispersal-linked attributes on realized dispersal and population genetic differentiation.

From the individual to the colony: Marine invertebrates as models to understand levels of biological organization

The developmental and evolutionary principles of coloniality in marine animals remain largely unexplored. Although many common traits have evolved independently in different groups of colonial animals, questions about their significance for colonial life histories remain unanswered. In 2018 (Nov. 25 - Dec. 8), the inaugural course on the Evolution of Coloniality and Modularity took place at the Center for Marine Biology of the University of São Paulo (CEBIMAR-USP), Brazil. During the intensive two-week graduate-level course, we addressed some of the historical ideas about animal coloniality by focal studies in bryozoans, tunicates, cnidarians, and sponges. We discussed many historical hypotheses and ways to test these using both extant and paleontological data, and we carried direct observations of animal colonies in the different phyla to address questions about coloniality. We covered topics related to multi-level selection theory and studied colonial traits, including modular miniaturization, polymorphism, brooding, and allorecognition. Course participants carried out short research projects using local species of animals to address questions on allorecognition and regeneration in ascidians and sponges, fusion and chimerism in anthoathecate hydrozoans, and evolution of polymorphism in bryozoans. Although many questions remain unanswered, this course served as a foundation to continue to develop a developmental and evolutionary synthesis of clonal and modular development in colonial marine organisms.

Inbreeding shapes the evolution of marine invertebrates

Inbreeding is a potent evolutionary force shaping the distribution of genetic variation within and among populations of plants and animals. Yet, our understanding of the forces shaping the expression and evolution of nonrandom mating in general, and inbreeding in particular, remains remarkably incomplete. Most research on plant mating systems focuses on self-fertilization and its consequences for automatic selection, inbreeding depression, purging, and reproductive assurance, whereas studies of animal mating systems have often assumed that inbreeding is rare, and that natural selection favors traits that promote outbreeding. Given that many sessile and sedentary marine invertebrates and marine macroalgae share key life history features with seed plants (e.g., low mobility, modular construction, and the release of gametes into the environment), their mating systems may be similar. Here, we show that published estimates of inbreeding coefficients (F_IS ) for sessile and sedentary marine organisms are similar and at least as high as noted in terrestrial seed plants. We also found that variation in F_IS within invertebrates is related to the potential to self-fertilize, disperse, and choose mates. The similarity of F_IS for these organismal groups suggests that inbreeding could play a larger role in the evolution of sessile and sedentary marine organisms than is currently recognized. Specifically, associations between traits of marine invertebrates and F_IS suggest that inbreeding could drive evolutionary transitions between hermaphroditism and separate sexes, direct development and multiphasic life cycles, and external and internal fertilization.

Genetic distance predicts trait differentiation at the subpopulation but not the individual level in eelgrass, Zostera marina

Ecological studies often assume that genetically similar individuals will be more similar in phenotypic traits, such that genetic diversity can serve as a proxy for trait diversity. Here, we explicitly test the relationship between genetic relatedness and trait distance using 40 eelgrass (Zostera marina) genotypes from five sites within Bodega Harbor, CA. We measured traits related to nutrient uptake, morphology, biomass and growth, photosynthesis, and chemical deterrents for all genotypes. We used these trait measurements to calculate a multivariate pairwise trait distance for all possible genotype combinations. We then estimated pairwise relatedness from 11 microsatellite markers. We found significant trait variation among genotypes for nearly every measured trait; however, there was no evidence of a significant correlation between pairwise genetic relatedness and multivariate trait distance among individuals. However, at the subpopulation level (sites within a harbor), genetic (FST) and trait differentiation were positively correlated. Our work suggests that pairwise relatedness estimated from neutral marker loci is a poor proxy for trait differentiation between individual genotypes. It remains to be seen whether genomewide measures of genetic differentiation or easily measured "master" traits (like body size) might provide good predictions of overall trait differentiation.

The developmental transcriptomes of two sea biscuit species with differing larval types

BACKGROUND

Larval developmental patterns are extremely varied both between and within phyla, however the genetic mechanisms leading to this diversification are poorly understood. We assembled and compared the developmental transcriptomes for two sea biscuit species (Echinodermata: Echinoidea) with differing patterns of larval development, to provide a resource for investigating the evolution of alternate life cycles. One species (Clypeaster subdepressus) develops via an obligately feeding larva which metamorphoses 3-4 weeks after fertilization; the other (Clypeaster rosaceus) develops via a rare, intermediate larval type-facultative feeding- and can develop through metamorphosis entirely based on egg provisioning in under one week.

RESULTS

Overall, the two transcriptomes are highly similar, containing largely orthologous contigs with similar functional annotation. However, we found distinct differences in gene expression patterns between the two species. Larvae from C. rosaceus, the facultative planktotroph, turned genes on at earlier stages and had less differentiation in gene expression between larval stages, whereas, C. subdepressus showed a higher degree of stage-specific gene expression.

CONCLUSION

This study is the first genetic analysis of a species with facultatively feeding larvae. Our results are consistent with known developmental differences between the larval types and raise the question of whether earlier onset of developmental genes is a key step in the evolution of a reduced larval period. By publishing a transcriptome for this rare, intermediate, larval type, this study adds developmental breadth to the current genetic resources, which will provide a valuable tool for future research on echinoderm development as well as studies on the evolution of development in general.

Multiple dimensions of intraspecific diversity affect biomass of eelgrass and its associated community

Genetic diversity within key species can play an important role in the functioning of entire communities. However, the extent to which different dimensions of diversity (e.g., the number of genotypes vs. the extent of genetic differentiation among those genotypes) best predicts functioning is unknown and may yield clues into the different mechanisms underlying diversity effects. We explicitly test the relative influence of genotypic richness and genetic relatedness on eelgrass productivity, biomass, and the diversity of associated invertebrate grazers in a factorial field experiment using the seagrass species, Zostera marina (eelgrass). Genotypic richness had the strongest effect on eelgrass biomass accumulation, such that plots with more genotypes at the end of the experiment attained a higher biomass. Genotypic diversity (richness + evenness) was a stronger predictor of biomass than richness alone, and both genotype richness and diversity were positively correlated with trait diversity. The relatedness of genotypes in a plot reduced eelgrass biomass independently of richness. Plots containing eelgrass with greater trait diversity also had a higher abundance of invertebrate grazers, while the diversity and relatedness of eelgrass genotypes had little effect on invertebrate abundance or richness. Our work extends previous findings by explicitly relating genotypic diversity to trait diversity, thus mechanistically connecting genotypic diversity to plot-level yields. We also show that other dimensions of diversity, namely relatedness, influence eelgrass performance independent of trait differentiation. Ultimately, richness and relatedness captured fundamentally different components of intraspecific variation and should be treated as complementary rather than competing dimensions of biodiversity affecting ecosystem functioning.

Rarity and persistence

Rarity is a population characteristic that is usually associated with a high risk of extinction. We argue here, however, that chronically rare species (those with low population densities over many generations across their entire ranges) may have individual-level traits that make populations more resistant to extinction. The major obstacle to persistence at low density is successful fertilisation (union between egg and sperm), and chronically rare species are more likely to survive when (1) fertilisation occurs inside or close to an adult, (2) mate choice involves long-distance signals, (3) adults or their surrogate gamete dispersers are highly mobile, or (4) the two sexes are combined in a single individual. In contrast, external fertilisation and wind- or water-driven passive dispersal of gametes, or sluggish or sedentary adult life habits in the absence of gamete vectors, appear to be incompatible with sustained rarity. We suggest that the documented increase in frequency of these traits among marine genera over geological time could explain observed secular decreases in rates of background extinction. Unanswered questions remain about how common chronic rarity actually is, which traits are consistently associated with chronic rarity, and how chronically rare species are distributed among taxa, and among the world's ecosystems and regions.

LIFE-HISTORY VARIATION WITHIN A POPULATION OF THE COLONIAL ASCIDIAN BOTRYLLUS SCHLOSSERI. I. THE GENETIC AND ENVIRONMENTAL CONTROL OF SEASONAL VARIATION

Many empirical analyses of life-history tactics are based on the assumption that demographic variation ought to be greatest among populations or species living in different environments. However, in a single population of the sessile colonial sea squirt Botryllus schlosseri, there are two discrete life-history morphs. Semelparous colonies are characterized by a) death immediately following the production of a single clutch, b) early age at first reproduction, c) rapid growth to first reproduction, and d) high reproductive effort. In contrast, iteroparous colonies a) produce at least three clutches before dying, b) postpone sexual reproduction until they are nearly twice the age of semelparous colonies, c) grow at about half the rate of semelparous colonies, and d) invest roughly 75% less in reproductive effort than semelparous colonies. Semelparous colonies numerically dominate the population through midsummer; later in the summer, iteroparous colonies are most numerous. Field and laboratory common-garden experiments, along with breeding studies, indicate that the demographic differences between the morphs are genetically determined. Consequently, the seasonal switch from dominance by semelparous colonies to dominance by iteroparous colonies may be an evolved response to a seasonally changing environment. On theoretical grounds, temporal variation in selection is thought to play a relatively unimportant role in maintaining genetic polymorphism; nonetheless, the seasonally recurrent life-history polymorphism shown in this study indicates that temporal variation in selection can lead to the maintenance of genetic polymorphism for traits strongly affecting fitness.

EVOLUTIONARY GENETICS OF ALLORECOGNITION IN THE COLONIAL HYDROID HYDRACTINIA SYMBIOLONGICARPUS

Many sedentary, clonal marine invertebrates compete intensively with conspecifics for habitable space. Allorecognition systems mediate the nature and outcome of these intraspecific competitive interactions, such that the initiation of agonistic behavior and the potential for intergenotypic fusion depend strongly on the relatedness of the contestants. The dependence of these behaviors on relatedness, along with the extraordinary precision with which self can be discriminated from nonself, suggest that allorecognition systems are highly polymorphic genetically. However, allotypic specificity of this sort could be produced by any number of genetic scenarios, ranging from relatively few loci with abundant allelic variation to numerous loci with relatively few alleles per locus. At this point, virtually nothing is known of the formal genetics of allorecognition in marine invertebrates; consequently, the evolutionary dynamics of such systems remain poorly understood. In this paper, we characterize the formal genetics of allorecognition in the marine hydrozoan Hydractinia symbiolongicarpus. Hydractinia symbiolongicarpus colonizes gastropod shells occupied by hermit crabs. When two or more individuals grow into contact, one of three outcomes ensues: fusion (compatibility), transitory fusion (a temporary state of compatibility), and rejection (incompatibility, often accompanied by the production of agonistic structures termed hyperplastic stolons). Observed patterns of compatibility between unrelated, half-sib pairs, and full-sib pairs show that unrelated and half-sib pairs under laboratory culture have a very low probability of being fusible, whereas full sibs have a roughly 30% rate of fusion in experimental pairings. The genetic simulations indicate that roughly five loci, with 5-7 alleles per locus, confer specificity in this species. In ecological terms, the reproductive ecology of H. symbiolongicarpus should promote the cosettlement of kin, some of which should be full sibs, and some half sibs. Thus, there is potential for kin selection to play a major role in the evolution of the H. symbiolongicarpus allorecognition system. In genetic terms, this system conforms to theoretical predictions for a recognition system selected to distinguish among classes of kin, in addition to self from nonself.

Life-history predicts past and present population connectivity in two sympatric sea stars

Life‐history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic‐developing sea star, Parvulastra exigua, has lower levels of within‐ and among‐population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic‐developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic‐developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post‐Pleistocene introgression.

Temporal stability in patterns of genetic diversity and structure of a marine foundation species (Zostera marina)

Genetic diversity and population structure reflect complex interactions among a diverse set of processes that may vary temporally, limiting their potential to predict ecological and evolutionary outcomes. Yet, the stability of these patterns is rarely tested. We resampled eelgrass (Zostera marina) meadows from published studies to determine variability in genetic diversity and structure within and between meadows over 5-12 years. The meadows sampled (San Francisco, Tomales and Bodega Bays in California and the Virginia coastal bays) represent a range of life histories (annual vs perennial), age (well-established vs restored) and environments (rural vs urbanized). In all of these systems, neither diversity nor differentiation (F_ST) changed over time. Differences among tidal heights within Bodega Bay were also remarkably consistent, with the high intertidal being more diverse than the subtidal, and tidal height differentiation being modest but significant at both time points. Historical studies used only a few microsatellite loci; therefore, our temporal comparisons were based on 4-5 loci. However, analysis of the current data using a set of 12 loci show that 4-5 loci are sufficient to describe diversity and differentiation patterns in this system. This temporal consistency was not because of the resampling of large clones, underscoring the feasibility and relevance of understanding drivers of the differences. Because seagrasses are declining at rapid rates, restoration and conservation are increasingly a coastal management priority. Our results argue that surveys of eelgrass genetic structure and diversity at decadal scales can provide accurate depictions of populations, increasing the utility of published genetic data for restoration and designing networks of reserves.

Patterns of Mass Mortality among Rocky Shore Invertebrates across 100 km of Northeastern Pacific Coastline

Mass mortalities in natural populations, particularly those that leave few survivors over large spatial areas, may cause long-term ecological perturbations. Yet mass mortalities may remain undocumented or poorly described due to challenges in responding rapidly to unforeseen events, scarcity of baseline data, and difficulties in quantifying rare or patchily distributed species, especially in remote or marine systems. Better chronicling the geographic pattern and intensity of mass mortalities is especially critical in the face of global changes predicted to alter regional disturbance regimes. Here, we couple replicated post-mortality surveys with preceding long-term surveys and historical data to describe a rapid and severe mass mortality of rocky shore invertebrates along the north-central California coast of the northeastern Pacific Ocean. In late August 2011, formerly abundant intertidal populations of the purple sea urchin (Strongylocentrotus purpuratus, a well-known ecosystem engineer), and the predatory six-armed sea star (Leptasterias sp.) were functionally extirpated from ~100 km of coastline. Other invertebrates, including the gumboot chiton (Cryptochiton stelleri) the ochre sea star (Pisaster ochraceus), and subtidal populations of purple sea urchins also exhibited elevated mortality. The pattern and extent of mortality suggest the potential for long-term population, community, and ecosystem consequences, recovery from which may depend on the different dispersal abilities of the affected species.

Gene co-expression modules underlying polymorphic and monomorphic zooids in the colonial hydrozoan, Hydractinia symbiolongicarpus

Advances in sequencing technology have forced a quantitative revolution in Evolutionary Biology. One important feature of this renaissance is that comprehensive genomic resources can be obtained quickly for almost any taxon, thus speeding the development of new model organisms. Here, we analyze 20 RNA-seq libraries from morphologically, sexually, and genetically distinct polyp types from the gonochoristic colonial hydrozoan, Hydractinia symbiolongicarpus (Cnidaria). Analyses of these data using weighted gene co-expression networks highlight deeply conserved genetic elements of animal spermatogenesis and demonstrate the utility of these methods in identifying modules of genes that correlate with different zooid types across various statistical contrasts. RNA-seq data and analytical scripts described here are deposited in publicly available databases.

Kinship and the evolution of social behaviours in the sea

Until recently, little attention has been paid to the existence of kin structure in the sea, despite the fact that many marine organisms are sessile or sedentary. This lack of attention to kin structure, and its impacts on social evolution, historically stems from the pervasive assumption that the dispersal of gametes and larvae is almost always sufficient to prevent any persistent associations of closely related offspring or adults. However, growing evidence, both theoretical and empirical, casts doubt on the generality of this assumption, not only in species with limited dispersal, but also in species with long dispersive phases. Moreover, many marine organisms either internally brood their progeny or package them in nurseries, both of which provide ample opportunities for kinship to influence the nature and outcomes of social interactions among family members. As the evidence for kin structure within marine populations mounts, it follows that kin selection may play a far greater role in the evolution of both behaviours and life histories of marine organisms than is presently appreciated.