Population structure and genetic diversity among eelgrass (Zostera marina) beds and depths in San Francisco Bay

Facultative Annual Life Cycles in Seagrasses

Plant species usually have either annual or perennial life cycles, but facultative annual species have annual or perennial populations depending on their environment. In terrestrial angiosperms, facultative annual species are rare, with wild rice being one of the few examples. Our review shows that in marine angiosperms (seagrasses) facultative annual species are more common: six (of 63) seagrass species are facultative annual. It concerns Zostera marina, Z. japonica, Halophila decipiens, H. beccarii, Ruppia maritima, and R. spiralis. The annual populations generally produce five times more seeds than their conspecific perennial populations. Facultative annual seagrass species occur worldwide. Populations of seagrasses are commonly perennial, but the facultative annual species had annual populations when exposed to desiccation, anoxia-related factors, shading, or heat stress. A system-wide 'experiment' (closure of two out of three connected estuaries for large-scale coastal protection works) showed that the initial annual Z. marina population could shift to a perennial life cycle within 5 years, depending on environmental circumstances. We discuss potential mechanisms and implications for plant culture. Further exploration of flexible life histories in plant species, and seagrasses in particular, may aid in answering questions about trade-offs between vegetative and sexual reproduction, and preprogrammed senescence.

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.

Reciprocal Field Transplant Experiment and Comparative Transcriptome Analysis Provide Insights Into Differences in Seed Germination Time of Two Populations From Different Geographic Regions of Zostera marina L

Seagrasses are the only submerged marine higher plants, which can colonize the sea through sexual (via seeds) reproduction. The transition between seed dormancy and germination is an important ecological trait and a key stage in the life cycle of higher plants. According to our observations, the seeds of Zostera marina L. (eelgrass) in Swan Lake (SL) and Qingdao Bay (QB) in northern China have the same maturation time (summer) but different germination time. To investigate this phenomenon, we further carried out reciprocal transplantation experiment and transcriptome analysis. Results revealed that differences in the seed germination time between the two sites do exist and are determined by internal molecular mechanisms as opposed to environmental factors. Furthermore, we conducted comparative transcriptome analysis of seeds at the mature and early germination stages in both locations. The results that the number of genes related to energy, hormone and cell changes was higher in SL than in QB, could account for that the dormancy depth of seeds in SL was deeper than that in QB; consequently, the seeds in SL needed to mobilize more related genes to break dormancy and start germination. The results could have important practical implications for seagrass meadow restoration via seeds and provide in-depth and comprehensive data for understanding the molecular mechanisms related to seagrass seed germination.

Global Diversity and Biogeography of the Zostera marina Mycobiome

Seagrasses are marine flowering plants that provide critical ecosystem services in coastal environments worldwide. Marine fungi are often overlooked in microbiome and seagrass studies, despite terrestrial fungi having critical functional roles as decomposers, pathogens, or endophytes in global ecosystems. Here, we characterize the distribution of fungi associated with the seagrass Zostera marina, using leaves, roots, and rhizosphere sediment from 16 locations across its full biogeographic range. Using high-throughput sequencing of the ribosomal internal transcribed spacer (ITS) region and 18S rRNA gene, we first measured fungal community composition and diversity. We then tested hypotheses of neutral community assembly theory and the degree to which deviations suggested that amplicon sequence variants (ASVs) were plant selected or dispersal limited. Finally, we identified a core mycobiome and investigated the global distribution of differentially abundant ASVs. We found that the fungal community is significantly different between sites and that the leaf mycobiome follows a weak but significant pattern of distance decay in the Pacific Ocean. Generally, there was evidence for both deterministic and stochastic factors contributing to community assembly of the mycobiome, with most taxa assembling through stochastic processes. The Z. marina core leaf and root mycobiomes were dominated by unclassified Sordariomycetes spp., unclassified Chytridiomycota lineages (including Lobulomycetaceae spp.), unclassified Capnodiales spp., and Saccharomyces sp. It is clear from the many unclassified fungal ASVs and fungal functional guilds that knowledge of marine fungi is still rudimentary. Further studies characterizing seagrass-associated fungi are needed to understand the roles of these microorganisms generally and when associated with seagrasses.

IMPORTANCE Fungi have important functional roles when associated with land plants, yet very little is known about the roles of fungi associated with marine plants, like seagrasses. In this study, we report the results of a global effort to characterize the fungi associated with the seagrass Zostera marina across its full biogeographic range. Although we defined a putative global core fungal community, it is apparent from the many fungal sequences and predicted functional guilds that had no matches to existing databases that general knowledge of seagrass-associated fungi and marine fungi is lacking. This work serves as an important foundational step toward future work investigating the functional ramifications of fungi in the marine ecosystem.

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.

New Insights into Different Reproductive Effort and Sexual Recruitment Contribution between Two Geographic Zostera marina L. Populations in Temperate China

Seagrasses are important components of global coastal ecosystems, and the eelgrass Zostera marina L. is widely distributed along the Atlantic and Pacific coasts in the temperate northern hemisphere, but limited datum related to the contribution of sexual reproduction to population recruitment have been reported. This study aimed to understand eelgrass sexual reproduction and population recruitment in Swan Lake (SLL), and Huiquan Bay (HQB) was included for comparison. Random sampling, permanent quadrats or cores and laboratory seed germination-based experimental methods were employed. The flowering, seed production, seed banks, seed germination, seedling survival, and seedling growth of eelgrass were investigated from July 2014 to December 2015 to evaluate the contribution of sexual reproduction to population recruitment. Results indicated a dominant role of asexual reproduction in HQB, while sexual reproduction played a relatively important role in SLL. The highest flowering shoot density in SLL was 517.27 ± 504.29 shoots m^-2 (June) and represented 53.34% of the total shoots at the center site. The potential seed output per reproductive shoot and per unit area in SLL were 103.67 ± 37.95 seeds shoot^-1 and 53,623.66 ± 19,628.11 seeds m^-2, respectively. The maximum seed bank density in SLL was 552.21 ± 204.94 seeds m^-2 (October). Seed germination mainly occurred from the middle of March to the end of May, and the highest seedling density was 296.88 ± 274.27 seedlings m^-2 in April. The recruitment from seedlings accounted for 41.36% of the Z. marina population recruitment at the center site, while the sexual recruitment contribution at the patch site (50.52%) was greater than that at the center site. Seeds in SLL were acclimated to spring germination, while in HQB, they were acclimated to autumn germination (early October-late November). Seed bank density in HQB was very low, with a value of 254.35 ± 613.34 seeds m^-2 (early October). However, seeds in HQB were significantly larger and heavier than those in SLL (size: P = 0.004; weight: P < 0.001). The recruitment from seedlings accounted for as low as 2.53% of the Z. marina population recruitment in HQB. Our laboratory seed germination experiment, which was conducted in autumn, showed that the seed germination percent in HQB was significantly greater than in SLL at optimal germination temperatures (10 and 15°C; P < 0.001). A laboratory seed germination test at suitable temperature may be a potential novel approach to identify the ecological differences among different geographic populations. It is suggested that the Z. marina population recruitment may have different strategies and adapt to specific local conditions, such as in SLL and HQB, and the temperature regime may control morphological and phonological variations.

Genotypic Diversity and Short-term Response to Shading Stress in a Threatened Seagrass: Does Low Diversity Mean Low Resilience?

Seagrasses that are predominantly clonal often have low levels of genetic variation within populations and predicting their response to changing conditions requires an understanding of whether genetic variation confers increased resistance to environmental stressors. A higher level of genetic diversity is assumed to benefit threatened species due to the increased likelihood of those populations having genotypes that can persist under environmental change. To test this idea, we conducted an in situ shading experiment with six geographically distinct meadows of the threatened seagrass Posidonia australis that vary in genetic diversity. Different genotypes within meadows varied widely in their physiological and growth responses to reduced light during a simulated short-term turbidity event. The majority of meadows were resistant to the sudden reduction in light availability, but a small subset of meadows with low genotypic diversity were particularly vulnerable to the early effects of shading, showing substantially reduced growth rates after only 3 weeks. Using the photosynthetic performance (maximum quantum yield) of known genotypes, we simulated meadows of varying genetic diversity to show that higher diversity can increase meadow resilience to stress by ensuring a high probability of including a high-performing genotype. These results support the hypothesis that complementarity among genotypes enhances the adaptive capacity of a population, and have significant implications for the conservation of declining P. australis meadows close to the species range edge on the east coast of Australia, where the genotypic diversity is low.

Population genetic structure of eelgrass (Zostera marina) on the Korean coast: Current status and conservation implications for future management

Seagrasses provide numerous ecosystem services for coastal and estuarine environments, such as nursery functions, erosion protection, pollution filtration, and carbon sequestration. Zostera marina (common name “eelgrass”) is one of the seagrass bed-forming species distributed widely in the northern hemisphere, including the Korean Peninsula. Recently, however, there has been a drastic decline in the population size of Z. marina worldwide, including Korea. We examined the current population genetic status of this species on the southern coast of Korea by estimating the levels of genetic diversity and genetic structure of 10 geographic populations using eight nuclear microsatellite markers. The level of genetic diversity was found to be significantly lower for populations on Jeju Island [mean allelic richness (AR) = 1.92, clonal diversity (R) = 0.51], which is located approximately 155 km off the southernmost region of the Korean Peninsula, than for those in the South Sea (mean AR = 2.69, R = 0.82), which is on the southern coast of the mainland. South Korean eelgrass populations were substantially genetically divergent from one another (F_ST = 0.061–0.573), suggesting that limited contemporary gene flow has been taking place among populations. We also found weak but detectable temporal variation in genetic structure within a site over 10 years. In additional depth comparisons, statistically significant genetic differentiation was observed between shallow (or middle) and deep zones in two of three sites tested. Depleted genetic diversity, small effective population sizes (N_e) and limited connectivity for populations on Jeju Island indicate that these populations may be vulnerable to local extinction under changing environmental conditions, especially given that Jeju Island is one of the fastest warming regions around the world. Overall, our work will inform conservation and restoration efforts, including transplantation for eelgrass populations at the southern tip of the Korean Peninsula, for this ecologically important species.

Limited effects of source population identity and number on seagrass transplant performance

Global declines in coastal foundation species highlight the importance of effective restoration. In this study, we examined the effects of source population identity and diversity (one vs. three sources per plot) on seagrass (Zostera marina) transplant success. The field experiment was replicated at two locations in Massachusetts with adjacent natural Zostera marina beds to test for local adaptation and source diversity effects on shoot density. We also collected morphological and genetic data to characterize variation within and among source populations, and evaluate whether they were related to performance. Transplants grew and expanded until six months post-transplantation, but then steadily declined at both sites. Prior to declines, we observed variation in performance among source populations at one site that was related to morphological traits: the populations with the longest leaves had the highest shoot densities, whereas the population with the shortest leaves performed the worst at six months post-transplantation. In addition, multiple source plots at this same transplant site consistently had similar or higher shoot densities than single source plots, and shoots from weak-performing populations showed improved performance in multiple source plots. We found no evidence for home site advantage or benefits of population-level genetic variation in early transplant performance at either site. Our results show limited effects of source population on early transplant performance and suggest that factors (e.g., morphology) other than home site advantage and population genetic variation serve a role. Based on our overall findings that transplant success varied among source populations and that population diversity at the plot level had positive but limited effects on individual and plot performance, we support planting shoots from multiple source sites in combination to enhance transplant success, particularly in the absence of detailed information on individual source characteristics.

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.

The Structure of Genetic Diversity in Eelgrass (Zostera marina L.) along the North Pacific and Bering Sea Coasts of Alaska

Eelgrass (Zostera marina) populations occupying coastal waters of Alaska are separated by a peninsula and island archipelago into two Large Marine Ecosystems (LMEs). From populations in both LMEs, we characterize genetic diversity, population structure, and polarity in gene flow using nuclear microsatellite fragment and chloroplast and nuclear sequence data. An inverse relationship between genetic diversity and latitude was observed (heterozygosity: R2 = 0.738, P < 0.001; allelic richness: R2 = 0.327, P = 0.047), as was significant genetic partitioning across most sampling sites (θ = 0.302, P < 0.0001). Variance in allele frequency was significantly partitioned by region only in cases when a population geographically in the Gulf of Alaska LME (Kinzarof Lagoon) was instead included with populations in the Eastern Bering Sea LME (θp = 0.128-0.172; P < 0.003), suggesting gene flow between the two LMEs in this region. Gene flow among locales was rarely symmetrical, with notable exceptions generally following net coastal ocean current direction. Genetic data failed to support recent proposals that multiple Zostera species (i.e. Z. japonica and Z. angustifolia) are codistributed with Z. marina in Alaska. Comparative analyses also failed to support the hypothesis that eelgrass populations in the North Atlantic derived from eelgrass retained in northeastern Pacific Last Glacial Maximum refugia. These data suggest northeastern Pacific populations are derived from populations expanding northward from temperate populations following climate amelioration at the terminus of the last Pleistocene glaciation.

Conservation of eelgrass (Zostera marina) genetic diversity in a mesocosm-based restoration experiment

Eelgrass (Zostera marina) forms the foundation of an important shallow coastal community in protected estuaries and bays. Widespread population declines have stimulated restoration efforts, but these have often overlooked the importance of maintaining the evolutionary potential of restored populations by minimizing the reduction in genetic diversity that typically accompanies restoration. In an experiment simulating a small-scale restoration, we tested the effectiveness of a buoy-deployed seeding technique to maintain genetic diversity comparable to the seed source populations. Seeds from three extant source populations in San Francisco Bay were introduced into eighteen flow-through baywater mesocosms. Following seedling establishment, we used seven polymorphic microsatellite loci to compare genetic diversity indices from 128 shoots to those found in the source populations. Importantly, allelic richness and expected heterozygosity were not significantly reduced in the mesocosms, which also preserved the strong population differentiation present among source populations. However, the inbreeding coefficient F_IS was elevated in two of the three sets of mesocosms when they were grouped according to their source population. This is probably a Wahlund effect from confining all half-siblings within each spathe to a single mesocosm, elevating F_IS when the mesocosms were considered together. The conservation of most alleles and preservation of expected heterozygosity suggests that this seeding technique is an improvement over whole-shoot transplantation in the conservation of genetic diversity in eelgrass restoration efforts.

Scaling of processes shaping the clonal dynamics and genetic mosaic of seagrasses through temporal genetic monitoring

Theoretically, the dynamics of clonal and genetic diversities of clonal plant populations are strongly influenced by the competition among clones and rate of seedling recruitment, but little empirical assessment has been made of such dynamics through temporal genetic surveys. We aimed to quantify 3 years of evolution in the clonal and genetic composition of Zostera marina meadows, comparing parameters describing clonal architecture and genetic diversity at nine microsatellite markers. Variations in clonal structure revealed a decrease in the evenness of ramet distribution among genets. This illustrates the increasing dominance of some clonal lineages (multilocus lineages, MLLs) in populations. Despite the persistence of these MLLs over time, genetic differentiation was much stronger in time than in space, at the local scale. Contrastingly with the short-term evolution of clonal architecture, the patterns of genetic structure and genetic diversity sensu stricto (that is, heterozygosity and allelic richness) were stable in time. These results suggest the coexistence of (i) a fine grained (at the scale of a 20 × 30 m quadrat) stable core of persistent genets originating from an initial seedling recruitment and developing spatial dominance through clonal elongation; and (ii) a local (at the scale of the meadow) pool of transient genets subjected to annual turnover. This simultaneous occurrence of initial and repeated recruitment strategies highlights the different spatial scales at which distinct evolutionary drivers and mating systems (clonal competition, clonal growth, propagule dispersal and so on) operate to shape the dynamics of populations and the evolution of polymorphism in space and time.