Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Loading...
Thumbnail Image

Authors

Yu, Lei
Khachaturyan, Marina
Matschiner, Michael
Healey, Adam
Bauer, Diane
Cameron, Brenda
Cusson, Mathieu
Duffy, J. Emmett
Fodrie, F. Joel
Gill, Diana

Journal Title

Journal ISSN

Volume Title

Publisher

Nature Research

Abstract

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Description

DATA AVAILABILITY : Genome data have been deposited in Genbank (short read archive, Supplementary Data Table 3). Coding sequences of Z. japonica and Z. marina for the ASTRAL analysis can be found on figshare (https://doi.org/10.6084/m9.figshare.21626327.v1). VCF files of the 11,705 core SNPs can be accessed at https://doi.org/10.6084/m9.figshare.21629471.v1. Source data for Fig. 1b,c are given, as well as statistics of sequencing coverage, mapping rate and further specifications of each sequenced library (Supplementary Tables 1–3). Source data are provided with this paper.
CODE AVAILABILITY : Custom-made scripts are deposited on GitHub for SNP filtering (github.com/leiyu37/populationGenomics_ZM.git), for clone mate detection (github.com/leiyu37/Detecting-clonemates.git), for heterozygote and nucleotide diversity quantification (github.com/leiyu37/populationGenomics_ZM.git) and to prepare SplitsTree input files (https://github.com/leiyu37/populationGenomics_ZM/blob/main/10_SplitsTree/vcf2alignment.py) and SNAPP input files (github.com/mmatschiner/snapp_prep). Scripts for calculating D-statistics are available at github.com/mmatschiner/tutorials/blob/master/analysis_of_introgression_with_snp_data/src/plot_d.rb. Scripts to prepare the gene presence/absence analysis are deposited on https://github.com/leiyu37/populationGenomics_ZM/tree/main/gene_presense_absence_analysis. Further software code for the MSMC analysis are found at http://lh3lh3.users.sourceforge.net/snpable.shtml (generation of mappability mask file for each of six chromosomes using SNPable) and at https://github.com/stschiff/msmc-tools (generation of ramet-specific mask file based on a bam file using bamCaller.py).
SUPPLEMENTARY MATERIAL : Supplementary Notes 1–8, Tables 1–6 and Figs. 1–12.
SUPPLEMENTARY DATA : Table 1: Sequence coverage. Supplementary Data Table 2: Mapping rate. Supplementary Data Table 3: Accession number of each library.
SOURCE DATA : Fig. 1b,c.

Keywords

Marine biology, Plant evolution, Population genetics, SDG-14: Life below water

Sustainable Development Goals

Life below water

Citation

Yu, L., Khachaturyan, M., Matschiner, M. et al. Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina). Nature Plants 9, 1207–1220 (2023). https://doi.org/10.1038/s41477-023-01464-3.