Development and characterization of microsatellite loci for the haploid-diploid red seaweed Gracilaria vermiculophylla

Phenology and thallus size in a non-native population of Gracilaria vermiculophylla

Phenology, or seasonal variation in life cycle events, is poorly described for many macroalgal species. We describe the phenology of a non-native population of Gracilaria vermiculophylla whose thalli are free-living or anchored by decorating polychaetes to tube caps. At a site in South Carolina, USA, we sampled 100 thalli approximately every month from January 2014 to January 2015. We assessed the reproductive state and measured thallus size based on wet weight, thallus length, and thallus surface area from herbarium mounts. Because life cycle stage cannot be assigned using morphology, we implemented a PCR assay to determine the life cycle stage-tetrasporophyte, female gametophyte, or male gametophyte-of each thallus. Tetrasporophytes dominated throughout the year, making up 81%-100% of thalli sampled per month. Reproductive tetrasporophytes varied between 0% and 65% of monthly samples and were most common in warm summer months (July through September) when thalli also tended to be larger. The vast majority of the reproductive thalli were worm-anchored and not fixed to hard substratum via a holdfast. Thus, free-living thalli can be reproductive and potentially seed new non-native populations. Given G. vermiculophylla reproduction seems tied closely to temperature, our work suggests phenology may change with climate-related changes in seawater temperatures. We also highlight the importance of understanding the natural history of macroalgae to better understand the consequence of range expansions on population dynamics.

Founder effects shape linkage disequilibrium and genomic diversity of a partially clonal invader

The genomic variation of an invasive species may be affected by complex demographic histories and evolutionary changes during the invasion. Here, we describe the relative influence of bottlenecks, clonality, and population expansion in determining genomic variability of the widespread red macroalga Agarophyton vermiculophyllum. Its introduction from mainland Japan to the estuaries of North America and Europe coincided with shifts from predominantly sexual to partially clonal reproduction and rapid adaptive evolution. A survey of 62,285 SNPs for 351 individuals from 35 populations, aligned to 24 chromosome-length scaffolds indicate that linkage disequilibrium (LD), observed heterozygosity (H_o ), Tajima's D, and nucleotide diversity (Pi) were greater among non-native than native populations. Evolutionary simulations indicate LD and Tajima's D were consistent with a severe population bottleneck. Also, the increased rate of clonal reproduction in the non-native range could not have produced the observed patterns by itself but may have magnified the bottleneck effect on LD. Elevated marker diversity in the genetic source populations could have contributed to the increased H_o and Pi observed in the non-native range. We refined the previous invasion source region to a ~50 km section of northeastern Honshu Island. Outlier detection methods failed to reveal any consistently differentiated loci shared among invaded regions, probably because of the complex A. vermiculophyllum demographic history. Our results reinforce the importance of demographic history, specifically founder effects, in driving genomic variation of invasive populations, even when localized adaptive evolution and reproductive system shifts are observed.

Using RAD-seq to develop sex-linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double-digest restriction site-associated DNA sequencing (ddRAD-seq) to identify sex-linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD-seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non-native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female-linked and 19 putatively male-linked sequences. Four female- and eight male-linked markers amplified in all three life cycle stages. Using one female- and one male-linked marker that were sex-specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non-native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD-seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex-linked markers in other haplodiplontic macroalgae for which genomes are lacking.

The Contribution of Clonality to Population Genetic Structure in the Sea Anemone, Diadumene lineata

Ecological and evolutionary processes differ depending on how genetic diversity is organized in space. For clonal organisms, the organization of both genetic and genotypic diversity can influence the fitness effects of competition, the mating system, and reproductive mode, which are key drivers of life cycle evolution. Understanding how individual reproductive behavior contributes to population genetic structure is essential for disentangling these forces, particularly in species with complex and plastic life cycles. The widespread sea anemone, Diadumene lineata, exhibits temperature-dependent fission, which contributes to predictable variation in clonal rate along the Atlantic coast of the United States, part of its non-native range. Because warmer conditions lead to higher rates of clonality, we expected to find lower genotypic and genetic diversity in lower versus higher latitude populations. We developed primers for 11 microsatellite loci and genotyped 207 anemones collected from 8 sites ranging from Florida to Massachusetts. We found clonal influence at all sites, and as predicted, the largest clones were found at lower latitude sites. We also found genetic signatures of sex in the parts of the range where gametogenesis is most common. Evidence of sex outside the native range is novel for this species and provides insights into the dynamics of this successful invader. Our findings also illustrate challenges that partially clonal taxa pose for eco-evolutionary studies, such as difficulty sampling statistically robust numbers of genets and interpretating common population genetic metrics. For example, we found high among-locus variation in FIS, which makes the meaning of mean multilocus FIS unclear.

How do microbiota associated with an invasive seaweed vary across scales?

Communities are shaped by scale dependent processes. To study the diversity and variation of microbial communities across scales, the invasive and widespread seaweed Agarophyton vermiculophyllum presents a unique opportunity. We characterized pro- and eukaryotic communities associated with this holobiont across its known distribution range, which stretches over the northern hemisphere. Our data reveal that community composition and diversity in the holobiont vary at local but also larger geographic scales. While processes acting at the local scale (i.e., within population) are the main structuring drivers of associated microbial communities, changes in community composition also depend on processes acting at larger geographic scales. Interestingly, the largest analysed scale (i.e., native and non-native ranges) explained variation in the prevalence of predicted functional groups, which could suggest a functional shift in microbiota occurred over the course of the invasion process. While high variability in microbiota at the local scale supports A. vermiculophyllum to be a generalist host, we also identified a number of core taxa. These geographically independent holobiont members imply that cointroduction of specific microbiota may have additionally promoted the invasion process.

What's ploidy got to do with it? Understanding the evolutionary ecology of macroalgal invasions necessitates incorporating life cycle complexity

Biological invasions represent grave threats to terrestrial, aquatic, and marine ecosystems, but our understanding of the role of evolution during invasions remains rudimentary. In marine environments, macroalgae account for a large percentage of invaders, but their complicated life cycles render it difficult to move methodologies and predictions wholesale from species with a single, free-living ploidy stage, such as plants or animals. In haplodiplontic macroalgae, meiosis and fertilization are spatiotemporally separated by long-lived, multicellular haploid and diploid stages, and gametes are produced by mitosis, not meiosis. As a consequence, there are unique eco-evolutionary constraints that are not typically considered in invasions. First, selfing can occur in both monoicious (i.e., hermaphroditic) and dioicious (i.e., separate sexes) haplodiplontic macroalgae. In the former, fertilization between gametes produced by the same haploid thallus results in instantaneous, genome-wide homozygosity. In the latter, cross-fertilization between separate male and female haploids that share the same diploid parent is analogous to selfing in plants or animals. Separate sexes, therefore, cannot be used as a proxy for outcrossing. Second, selfing likely facilitates invasions (i.e., Baker's law) and the long-lived haploid stage may enable purging of deleterious mutations, further contributing to invasion success. Third, asexual reproduction will result in the dominance of one ploidy and/or sex and the loss of the other(s). Whether or not sexual reproduction can be recovered depends on which stage is maintained. Finally, fourth, haplodiplontic life cycles are predicted to be maintained through niche differentiation in the haploid and diploid stages. Empirical tests are rare, but fundamental to our understanding of macroalgal invasion dynamics. By highlighting these four phenomena, we can build a framework with which to empirically and theoretically address important gaps in the literature on marine evolutionary ecology, of which biological invasions can serve as unnatural laboratories.

Nonnative Gracilaria vermiculophylla tetrasporophytes are more difficult to debranch and are less nutritious than gametophytes

Theory predicts that the maintenance of haplodiplontic life cycles requires ecological differences between the haploid gametophytes and diploid sporophytes, yet evidence of such differences remain scarce. The haplodiplontic red seaweed Gracilaria vermiculophylla has invaded the temperate estuaries of the Northern Hemisphere, where it commonly modifies detrital and trophic pathways. In native populations, abundant hard substratum enables spore settlement, and gametophyte:tetrasporophyte ratios are ~40:60. In contrast, many non-native populations persist in soft-sediment habitats without abundant hard substratum, and can be 90%-100% tetrasporophytic. To test for ecologically relevant phenotypic differences, we measured thallus morphology, protein content, organic content, "debranching resistance" (i.e., tensile force required to remove a branch from its main axis node), and material properties between male gametophytes, female gametophytes, and tetrasporophytes from a single, nonnative site in Charleston Harbor, South Carolina, USA in 2015 and 2016. Thallus length and surface area to volume ratio differed between years, but were not significantly different between ploidies. Tetrasporophytes had lower protein content than gametophytes, suggesting the latter may be more attractive to consumers. More force was required to pull a branch from the main axis of tetrasporophytes relative to gametophytes. A difference in debranching resistance may help to maintain tetrasporophyte thallus durability relative to gametophytes, providing a potential advantage in free-floating populations. These data may shed light on the invasion ecology of an important ecosystem engineer, and may advance our understanding of life cycle evolution and the maintenance of life cycle diversity.

Combining niche shift and population genetic analyses predicts rapid phenotypic evolution during invasion

The rapid evolution of non-native species can facilitate invasion success, but recent reviews indicate that such microevolution rarely yields expansion of the climatic niche in the introduced habitats. However, because some invasions originate from a geographically restricted portion of the native species range and its climatic niche, it is possible that the frequency, direction, and magnitude of phenotypic evolution during invasion have been underestimated. We explored the utility of niche shift analyses in the red seaweed Gracilaria vermiculophylla, which expanded its range from the northeastern coastline of Japan to North America, Europe, and northwestern Africa within the last 100 years. A genetically informed climatic niche shift analysis indicates that native source populations occur in colder and highly seasonal habitats, while most non-native populations typically occur in warmer, less seasonal habitats. This climatic niche expansion predicts that non-native populations evolved greater tolerance for elevated heat conditions relative to native source populations. We assayed 935 field-collected and 325 common-garden thalli from 40 locations, and as predicted, non-native populations had greater tolerance for ecologically relevant extreme heat (40°C) than did Japanese source populations. Non-native populations also had greater tolerance for cold and low-salinity stresses relative to source populations. The importance of local adaptation to warm temperatures during invasion was reinforced by evolution of parallel clines: Populations from warmer, lower-latitude estuaries had greater heat tolerance than did populations from colder, higher-latitude estuaries in both Japan and eastern North America. We conclude that rapid evolution plays an important role in facilitating the invasion success of this and perhaps other non-native marine species. Genetically informed ecological niche analyses readily generate clear predictions of phenotypic shifts during invasions and may help to resolve debate over the frequency of niche conservatism versus rapid adaptation during invasion.

Genetic identification of source and likely vector of a widespread marine invader

The identification of native sources and vectors of introduced species informs their ecological and evolutionary history and may guide policies that seek to prevent future introductions. Population genetics provides a powerful set of tools to identify origins and vectors. However, these tools can mislead when the native range is poorly sampled or few molecular markers are used. Here, we traced the introduction of the Asian seaweed Gracilaria vermiculophylla (Rhodophyta) into estuaries in coastal western North America, the eastern United States, Europe, and northwestern Africa by genotyping more than 2,500 thalli from 37 native and 53 non‐native sites at mitochondrial cox1 and 10 nuclear microsatellite loci. Overall, greater than 90% of introduced thalli had a genetic signature similar to thalli sampled from the coastline of northeastern Japan, strongly indicating this region served as the principal source of the invasion. Notably, northeastern Japan exported the vast majority of the oyster Crassostrea gigas during the 20th century. The preponderance of evidence suggests G. vermiculophylla may have been inadvertently introduced with C. gigas shipments and that northeastern Japan is a common source region for estuarine invaders. Each invaded coastline reflected a complex mix of direct introductions from Japan and secondary introductions from other invaded coastlines. The spread of G. vermiculophylla along each coastline was likely facilitated by aquaculture, fishing, and boating activities. Our ability to document a source region was enabled by a robust sampling of locations and loci that previous studies lacked and strong phylogeographic structure along native coastlines.

The importance of effective sampling for exploring the population dynamics of haploid-diploid seaweeds

The mating system partitions genetic diversity within and among populations and the links between life history traits and mating systems have been extensively studied in diploid organisms. As such most evolutionary theory is focused on species for which sexual reproduction occurs between diploid male and diploid female individuals. However, there are many multicellular organisms with biphasic life cycles in which the haploid stage is prolonged and undergoes substantial somatic development. In particular, biphasic life cycles are found across green, brown and red macroalgae. Yet, few studies have addressed the population structure and genetic diversity in both the haploid and diploid stages in these life cycles. We have developed some broad guidelines with which to develop population genetic studies of haploid-diploid macroalgae and to quantify the relationship between power and sampling strategy. We address three common goals for studying macroalgal population dynamics, including haploid-diploid ratios, genetic structure and paternity analyses.