Population genetic structure of annual and perennial populations of Zostera marina L. along the Pacific coast of Baja California and the Gulf of California

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.

Environment predicts seagrass genotype, phenotype, and associated biodiversity in a temperate ecosystem

Coastal vegetative ecosystems are among the most threatened in the world, facing multiple anthropogenic stressors. A good example of this is seagrass, which supports carbon capture, coastal stabilization, and biodiversity, but is declining globally at an alarming rate. To understand the causes and consequences of changes to these ecosystems, we need to determine the linkages between different biotic and abiotic components. We used data on the seagrass, Zostera marina, collected by citizen scientists across 300 km of the south coast of the United Kingdom as a case study. We assembled data on seagrass genotype, phenotype, infauna, and associated bathymetry, light, sea surface temperature, and wave and current energy to test hypotheses on the distribution and diversity of this temperate sub-tidal ecosystem. We found spatial structure in population genetics, evident through local assortment of genotypes and isolation by distance across a broader geographic scale. By integrating our molecular data with information on seagrass phenotype and infauna, we demonstrate that these ecosystem components are primarily linked indirectly through the effects of shared environmental factors. It is unusual to examine genotypic, phenotypic, and environmental data in a single study, but this approach can inform both conservation and restoration of seagrass, as well as giving new insights into a widespread and important ecosystem.

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.

Using next generation sequencing of alpine plants to improve fecal metabarcoding diet analysis for Dall's sheep

Objectives

Dall’s sheep (Ovis dalli dalli) are important herbivores in the mountainous ecosystems of northwestern North America, and recent declines in some populations have sparked concern. Our aim was to improve capabilities for fecal metabarcoding diet analysis of Dall’s sheep and other herbivores by contributing new sequence data for arctic and alpine plants. This expanded reference library will provide critical reference sequence data that will facilitate metabarcoding diet analysis of Dall’s sheep and thus improve understanding of plant-animal interactions in a region undergoing rapid climate change.

Data description

We provide sequences for the chloroplast rbcL gene of 16 arctic-alpine vascular plant species that are known to comprise the diet of Dall’s sheep. These sequences contribute to a growing reference library that can be used in diet studies of arctic herbivores.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05590-z.

Population genomics reveals a mismatch between management and biological units in green abalone (Haliotis fulgens)

Effective fishery management strategies should be based on stock delimitation and knowledge of the spatial scale at which species are distributed. However, a mismatch often occurs between biological and management units of fishery resources. The green abalone (Haliotis fulgens) supports an important artisanal fishery in the west coast of the Baja California Peninsula (BCP), Mexico, which has shown a declining tendency despite the several management measures. Thus, the aim of this study was to characterize the spatial patterns of neutral genomic variation of green abalone along the BCP to test whether the genomic structure patterns support the current green abalone management areas. To test this hypothesis, a set of 2,170 putative neutral single nucleotide polymorphisms discovered by a double digest restriction-site associated DNA approach was used on 10 locations along the BCP. The results revealed a population structure with three putative groups: Guadalupe Island and northern and southern BCP locations. The contemporary gene flow might be explained by local oceanographic features, where it is bidirectional within the southern region but with a predominant southward flow from the northern region. These findings indicated that the administrative areas did not match the biological units of H. fulgens fishery; hence, the stock assessment and management areas should be revised.

Shared genomic outliers across two divergent population clusters of a highly threatened seagrass

The seagrass, Zostera capensis, occurs across a broad stretch of coastline and wide environmental gradients in estuaries and sheltered bays in southern and eastern Africa. Throughout its distribution, habitats are highly threatened and poorly protected, increasing the urgency of assessing the genomic variability of this keystone species. A pooled genomic approach was employed to obtain SNP data and examine neutral genomic variation and to identify potential outlier loci to assess differentiation across 12 populations across the ∼9,600 km distribution of Z. capensis. Results indicate high clonality and low genomic diversity within meadows, which combined with poor protection throughout its range, increases the vulnerability of this seagrass to further declines or local extinction. Shared variation at outlier loci potentially indicates local adaptation to temperature and precipitation gradients, with Isolation-by-Environment significantly contributing towards shaping spatial variation in Z. capensis. Our results indicate the presence of two population clusters, broadly corresponding to populations on the west and east coasts, with the two lineages shaped only by frequency differences of outlier loci. Notably, ensemble modelling of suitable seagrass habitat provides evidence that the clusters are linked to historical climate refugia around the Last Glacial Maxi-mum. Our work suggests a complex evolutionary history of Z. capensis in southern and eastern Africa that will require more effective protection in order to safeguard this important ecosystem engineer into the future.

Role of oceanography in shaping the genetic structure in the North Pacific hake Merluccius productus

Determining the relative influence of biotic and abiotic factors on genetic connectivity among populations remains a major challenge in evolutionary biology and in the management and conservation of species. North Pacific hake (Merluccius productus) inhabits upwelling regions in the California Current ecosystem from the Gulf of California to the Gulf of Alaska. In this study, we examined mitochondrial DNA (mtDNA) and microsatellite variation to estimate levels of genetic differentiation of M. productus in relation to the role of oceanographic features as potential barriers to gene flow. Samples were obtained from nine sites spanning a large part of the geographic range of the species, from Puget Sound, Washington to Costa Rica. The microsatellite results revealed three genetically discrete populations: one spanning the eastern Pacific coast, and two apparently resident populations circumscribed to the Puget Sound and the northern Gulf of California (F_ST = 0.032, p = 0.036). Cytochrome b sequence data indicated that isolation between the Puget Sound and northern Gulf of California populations from the coastal Pacific were recent phenomena (18.5 kyr for Puget Sound and 40 kyr for the northern Gulf of California). Oceanographic data obtained from the Gulf of California support the hypothesis that permanent fronts within the region, and strong gradients at the entrance to the Gulf of California act as barriers to gene flow. A seascape genetics approach found significant genetic–environment associations, where the daytime sea surface temperature and chlorophyll concentrations were the best predictive variables for the observed genetic differentiation. Considering the potential causes of genetic isolation among the three populations, e.g. spawning areas in different latitudes associated with upwelling processes, oceanographic barriers, asymmetric migration and specialized diet, oceanographic barriers appear to be a likely mechanism restricting gene flow.

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.

The Kuroshio current influences genetic diversity and population genetic structure of a tropical seagrass, Enhalus acoroides

Information on genetic diversity and differentiation of seagrass populations is essential for the conservation of coastal ecosystems. However, little is known about the seagrasses in the Indo-West Pacific Ocean, where the world's highest diversity of seagrasses occurs. The influence of sea currents on these populations is also unknown. We estimated the genetic diversity and population genetic structure and identified reproductive features in Enhalus acoroides populations from the Yaeyama Islands, Hainan Island and the Philippines. The Philippines are situated at the centre of the E. acoroides range, Yaeyama and Hainan are peripheral populations, and the Yaeyama population is at the northern limit of the species range. The powerful Kuroshio Current flows from the Philippines to Yaeyama. Genetic analyses using nine microsatellite markers indicated that reproduction of E. acoroides is mostly sexual. Clonal diversity does not decrease in northern populations, although genetic diversity does. However, the genetic diversity of the Yaeyama populations is greater than that of the Hainan populations. Significant genetic differentiation among most populations was evident; however, the Yaeyama and north-east Philippines populations were genetically similar, despite being separated by ~1100 km. An assignment test suggested that recruitment occurs from the north-east Philippines to Yaeyama. The strong current in this region is probably responsible for the extant genetic diversity and recruitment patterns.

Edge Effects and the Population Structure of Humboldt Bay, California, Eelgrass (Zostera marinaL.)

The physical structure of a habitat (e.g., fragmentation) may affect the distribution of genetic diversity within a population, and genetic diversity may alter ecological function. This study investigates the population genetics of Humboldt Bay, California, eelgrass in light of the habitat’s fragmented physical structure. Historical and ongoing dredging operations in the bay maintain channels between fragmented tidal mudflats that are inhabited by eelgrass. A sample of 469 individuals, from 11 discrete fragments, was genotyped using microsatellite analysis. Ramets were collected at multiple points located at the edge and in the interior of each fragment. The sampling scheme was designed to detect evidence of population structure within and between fragments, as well as between edge and interior bed positions. Genetic diversity is not shown to differ among fragments and is not consistent with widespread clonality. Heterozygosity levels indicate that Hardy-Weinberg equilibrium predominates across loci in most instances. Comparisons of edge and interior bed positions also reveal similarities between positions for the metrics described above. Results suggest efficient gene flow between fragments and between bed positions. Such findings, which indicate the unstructured nature of the population, can guide resource managers in making evolutionarily informed conservation decisions by discouraging assumptions of low genetic diversity and insignificant sexual reproduction.

Testing the genetic predictions of a biogeographical model in a dominant endemic Eastern Pacific coral (Porites panamensis) using a genetic seascape approach

The coral fauna of the Eastern Tropical Pacific (ETP) is depauperate and peripheral; hence, it has drawn attention to the factors allowing its survival. Here, we use a genetic seascape approach and ecological niche modeling to unravel the environmental factors correlating with the genetic variation of Porites panamensis, a hermatypic coral endemic to the ETP. Specifically, we test if levels of diversity and connectivity are higher among abundant than among depauperate populations, as expected by a geographically relaxed version of the Abundant Center Hypothesis (rel-ACH). Unlike the original ACH, referring to a geographical center of distribution of maximal abundance, the rel-ACH refers only to a center of maximum abundance, irrespective of its geographic position. The patterns of relative abundance of P. panamensis in the Mexican Pacific revealed that northern populations from Baja California represent its center of abundance; and southern depauperate populations along the continental margin are peripheral relative to it. Genetic patterns of diversity and structure of nuclear DNA sequences (ribosomal DNA and a single copy open reading frame) and five alloenzymatic loci partially agreed with rel-ACH predictions. We found higher diversity levels in peninsular populations and significant differentiation between peninsular and continental colonies. In addition, continental populations showed higher levels of differentiation and lower connectivity than peninsular populations in the absence of isolation by distance in each region. Some discrepancies with model expectations may relate to the influence of significant habitat discontinuities in the face of limited dispersal potential. Environmental data analyses and niche modeling allowed us to identify temperature, water clarity, and substrate availability as the main factors correlating with patterns of abundance, genetic diversity, and structure, which may hold the key to the survival of P. panamensis in the face of widespread environmental degradation.

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

The seagrass Zostera marina is widely distributed in coastal regions throughout much of the northern hemisphere, forms the foundation of an important ecological habitat, and is suffering population declines. Studies in the Atlantic and Pacific oceans indicate that the degree of population genetic differentiation is location dependent. San Francisco Bay, California, USA, is a high-current, high-wind environment where rafting of seed-bearing shoots has the potential to enhance genetic connectivity among Z. marina populations. We tested Z. marina from six locations, including one annual population, within the bay to assess population differentiation and to compare levels of within-population genetic diversity. Using 7 microsatellite loci, we found significant differentiation among all populations. The annual population had significantly higher clonal diversity than the others but showed no detectible differences in heterozygosity or allelic richness. There appears to be sufficient input of genetic variation through sexual reproduction or immigration into the perennial populations to prevent significant declines in the number and frequency of alleles. In additional depth comparisons, we found differentiation among deep and shallow portions in 1 of 3 beds evaluated. Genetic drift, sweepstakes recruitment, dispersal limitation, and possibly natural selection may have combined to produce genetic differentiation over a spatial scale of 3-30 km in Z. marina. This implies that the scale of genetic differentiation may be smaller than expected for seagrasses in other locations too. We suggest that populations in close proximity may not be interchangeable for use as restoration material.

Discordant distribution of populations and genetic variation in a sea star with high dispersal potential

Patiria miniata, a broadcast-spawning sea star species with high dispersal potential, has a geographic range in the intertidal zone of the northeast Pacific Ocean from Alaska to California that is characterized by a large range gap in Washington and Oregon. We analyzed spatial genetic variation across the P. miniata range using multilocus sequence data (mtDNA, nuclear introns) and multilocus genotype data (microsatellites). We found a strong phylogeographic break at Queen Charlotte Sound in British Columbia that was not in the location predicted by the geographical distribution of the populations. However, this population genetic discontinuity does correspond to previously described phylogeographic breaks in other species. Northern populations from Alaska and Haida Gwaii were strongly differentiated from all southern populations from Vancouver Island and California. Populations from Vancouver Island and California were undifferentiated with evidence of high gene flow or very recent separation across the range disjunction between them. The surprising and discordant spatial distribution of populations and alleles suggests that historical vicariance (possibly caused by glaciations) and contemporary dispersal barriers (possibly caused by oceanographic conditions) both shape population genetic structure in this species.

Drift by drift: effective population size is limited by advection

BACKGROUND

Genetic estimates of effective population size often generate surprising results, including dramatically low ratios of effective population size to census size. This is particularly true for many marine species, and this effect has been associated with hypotheses of "sweepstakes" reproduction and selective hitchhiking.

RESULTS

Here we show that in advective environments such as oceans and rivers, the mean asymmetric transport of passively dispersed reproductive propagules will act to limit the effective population size in species with a drifting developmental stage. As advection increases, effective population size becomes decoupled from census size as the persistence of novel genetic lineages is restricted to those that arise in a small upstream portion of the species domain.

CONCLUSION

This result leads to predictions about the maintenance of diversity in advective systems, and complements the "sweepstakes" hypothesis and other hypotheses proposed to explain cases of low allelic diversity in species with high fecundity. We describe the spatial extent of the species domain in which novel allelic diversity will be retained, thus determining how large an appropriately placed marine reserve must be to allow the persistence of endemic allelic diversity.

SOUTHERN HOSPITALITY: A LATITUDINAL GRADIENT IN GENE FLOW IN THE MARINE ENVIRONMENT

In recent years population genetics and phylogeographic studies have become increasingly valuable tools for inferring both historical and present-day genetic patterns within marine species. Here, we take a comparative approach to population-level study, analyzing original mitochondrial DNA data from 969 individuals representing 28 chiton (Mollusca: Polyplacophora) species to uncover large-scale genetic patterns along the Pacific coast of North America. The data reveal a distinct latitudinal connectivity gradient among chitons: species that exist at lower latitudes tend to have more isolated populations. This trend appears to be a product of between-species differences; within species, no significant gradient in connectivity is observed. Lower average annual sea surface temperatures are hypothesized to contribute to longer larval duration (and by extension, greater connectivity) among lecithotrophic species, providing a mechanism for the observed positive correlation between gene flow and latitude. Because increased isolation among populations may lead to speciation, a latitudinal trend in gene flow may contribute to the increased species diversity observed at lower latitudes.

Does disturbance enhance genotypic diversity in clonal organisms? A field test in the marine angiosperm Zostera marina

Physical disturbance has often been invoked to control genotypic diversity in sessile clonal organisms, yet experimental evidence is lacking. I studied the effects of physical disturbance on genet dynamics and genotypic diversity in a clonal marine angiosperm, Zostera marina (eelgrass). In replicated plots of 1 m2, the vegetation canopy was removed in gaps of zero (control), 25%, 50% and 75% of the area (n = 6 replicates). Before removal and during two consecutive years, the genotypic composition was determined using genetic markers (DNA microsatellites) in a 5 x 5 pixel grid per plot. An aggregate index of genet dynamics summarizing recruitment, increase, loss and decrease of clones was maximal at intermediate disturbance levels (quadratic polynomial P = 0.02). Physical disturbance also increased the occurrence of new genotypes, possibly reflecting recruitment (linear model, P < 0.05). Contrary to expectations, there was no competitive advantage of more heterozygous genotypes over less heterozygous ones. In the absence of disturbance, in particular, clones with lower individual heterozygosity were more likely to increase in area over a 1-year time period than more heterozygous ones, while there was no such correlation in plots with disturbance (logistic model, P(disturbance x heterozygosity) = 0.036). Undisturbed plots revealed background recruitment independent of canopy gaps, suggesting that Z. marina exhibits a strategy of continual recruitment. Effects of experimental disturbance (linear or quadratic) on clonal diversity were not detectable. Instead, initial (pretreatment) clonal diversity accounted for between 68% and 91% of the variance in diversity, indicating remarkable resilience of genotypic diversity in the face of physical disturbance.