http://lod.bco-dmo.org/id/dataset/924852 eng; USA utf8 dataset Highest level of data collection, from a common set of sensors or instrumentation, usually within the same research project Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact 2024-04-10 ISO 19115-2 Geographic Information - Metadata - Part 2: Extensions for Imagery and Gridded Data ISO 19115-2:2009(E) Sample collection and sequence accession information for Zostera marina whole genome resequencing from specimens collected at 16 geographic locations worldwide in 2017 2024-04-10 publication 2024-04-10 revision Marine Biological Laboratory/Woods Hole Oceanographic Institution Library (MBLWHOI DLA) 2024-07-02 publication https://doi.org/10.26008/1912/bco-dmo.924852.1 John J. Stachowicz University of California-Davis principalInvestigator Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. publisher Cite this dataset as: Stachowicz, J. J. (2024) Sample collection and sequence accession information for Zostera marina whole genome resequencing from specimens collected at 16 geographic locations worldwide in 2017. Biological and Chemical Oceanography Data Management Office (BCO-DMO). (Version 1) Version Date 2024-04-10 [if applicable, indicate subset used]. doi:10.26008/1912/bco-dmo.924852.1 [access date] Zostera marina whole genome resequencing Dataset Description: Methods and Sampling: &lt;p&gt;We conducted a range-wide sampling collection of 190 &lt;em&gt;Zostera marina&lt;/em&gt; specimens from 16 geographic locations (Yu et al. 2023 Nature Plants). The chosen populations feature a mix of sexual and vegetative reproduction with the exception of mostly vegetative reproduction at the sites PO and NN, apparent through extended clones. Chosen locations were a subset of the Zostera Experimental Network sites that were previously analysed using 24 microsatellite loci. Although a sampling distance of &amp;gt;2 meters was maintained to reduce the likelihood of collecting the same genet/clone twice, this was not always successful and thus provided an estimate of local clonal diversity. Plant tissue was selected from the basal meristematic part of the shoot after peeling away the leaf sheath to minimize epiphytes (bacteria and diatoms), frozen in liquid nitrogen and stored at -80 degrees Celsius (°C) until DNA extraction. Quality control was performed following Joint Genome Institute guidelines (&lt;a href=&quot;https://jgi.doe.gov/wp-content/uploads/2013/11/Genomic-DNA-Sample-QC.pdf&quot; target=&quot;_blank&quot;&gt;https://jgi.doe.gov/wp-content/uploads/2013/11/Genomic-DNA-Sample-QC.pdf&lt;/a&gt;). Plate-based DNA library preparation for Illumina sequencing was performed on the PerkinElmer Sciclone NGS robotic liquid handling system using Kapa Biosystems library preparation kit. About 200 nanograms (ng) of sample DNA was sheared to a length of around 600 base pairs (bp) using a Covaris LE220 focused ultrasonicator. Selected fragments were end-repaired, A-tailed, and ligated with sequencing adaptors containing a unique molecular index barcode. Libraries were quantified using KAPA Biosystems' next-generation sequencing library qPCR-kit on a Roche LightCycler 480 real-time PCR instrument. Quantified libraries were then pooled together and prepared for sequencing on the Illumina HiSeq2500 sequencer using TruSeq SBS sequencing kits (v4) following a 2 × 150 bp indexed run recipe to a targeted depth of approximately 40x coverage. The quality of the raw reads was assessed by FastQC (&lt;a href=&quot;https://www.bioinformatics.babraham.ac.uk/projects/fastqc/&quot; target=&quot;_blank&quot;&gt;https://www.bioinformatics.babraham.ac.uk/projects/fastqc/&lt;/a&gt;) and visualized by MultiQC58. BBDuk (&lt;a href=&quot;https://jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/&quot; target=&quot;_blank&quot;&gt;https://jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/&lt;/a&gt;) was used to remove adapters and for quality filtering, discarding sequence reads (1) with more than one 'N' (maxns = 1), (2) shorter than 50 bp after trimming (minlength = 50), and (3) with average quality &amp;lt;10 after trimming (maq = 10). FastQC and MultiQC were used for a second round of quality check for the clean reads. Sequencing coverage and mapping rate were calculated for each sample (see Yu et al. 2023 for details).&lt;/p&gt; Funding provided by NSF Division of Ocean Sciences (NSF OCE) Award Number: OCE-1829976 Award URL: https://www.nsf.gov/awardsearch/show-award?AWD_ID=1829976 completed John J. Stachowicz University of California-Davis 530-752-1113 Department of Evolution and Ecology 1 Shields Ave Davis CA 95616 USA jjstachowicz@ucdavis.edu pointOfContact asNeeded Dataset Version: 1 Unknown SampleID BioProject Library_Code Sample_Description Location_Name Lat Long SRA_accession Library_type Sampled_tissue Instrument reads Number_of_bases_Raw_fastq Coverage_Raw_fastq Number_of_bases_Clean_fastq Coverage_Clean_fastq Mapped_pcnt Properly_paired_pcnt Illumina HiSeq2500 sequencer theme None, User defined sample identification NCBI BioProject accession seq_descrip sample description site latitude longitude NCBI SRA study accession (SRP) instrument number featureType BCO-DMO Standard Parameters Automated DNA Sequencer instrument BCO-DMO Standard Instruments otherRestrictions otherRestrictions Access Constraints: none. Use Constraints: Please follow guidelines at: http://www.bco-dmo.org/terms-use Distribution liability: Under no circumstances shall BCO-DMO be liable for any direct, incidental, special, consequential, indirect, or punitive damages that result from the use of, or the inability to use, the materials in this data submission. If you are dissatisfied with any materials in this data submission your sole and exclusive remedy is to discontinue use. Using genomics to link traits to ecosystem function in the eelgrass Zostera marina https://osprey.bco-dmo.org/project/777102 Using genomics to link traits to ecosystem function in the eelgrass Zostera marina <p><em>NSF Award Abstract:</em><br /> Seagrass ecosystems provide important services to coastal regions, including primary production, carbon storage, nutrient cycling, habitat for fisheries species, and erosion control. At the same time, eelgrass is threatened by direct destruction, pollution, and other human impacts on the environment. We know that genetic diversity in eelgrass enhances seagrass bed growth and persistence, but application of this knowledge to restoration and conservation is limited. This work will guide restoration programs by considering what specific aspects of diversity are important to conservation and restoration of seagrass ecosystems, helping to guide the selection of source material to improve restoration success (which is often low). The project integrates the effects of multiple components of diversity and clarifies the extent to which genetic and ecological uniqueness can predict ecosystem functions.</p> <p>Intellectual Merit: Genetic diversity as measured by the number of genetically distinct individuals (genets) in an assemblage influences critical ecosystem functions in a wide range of ecosystems. Functional diversity, the presence of key traits, or population flexibility to respond to environmental change are all potential mechanisms underlying these patterns, but distinguishing among them requires a clear link between genetic diversity and the phenotypes present in an assemblage. The investigators, and others, have previously demonstrated that genet diversity in eelgrass (Zostera marina) increases stand productivity, animal community diversity, and resilience to environmental change. These genet diversity effects are associated with increases in genetically determined trait diversity. Predicting trait diversity without having to measure traits of every genet remains a major barrier to wider application of functional diversity approaches in restoration and management. In this project, the investigators assess the association between Single Nucleotide Polymorphisms (SNPs) across the genome and performance-related traits that we will measure at the individual, population, and seascape-scale. They also assess environmental correlates of trait differentiation from field sampling. Finally, the research team will compare the predictive power of genomic SNP diversity versus other metrics of intraspecific diversity for the functioning (productivity, invertebrate abundance) of field planted eelgrass assemblages. If genomic variation can reliably be used to predict functional traits, then the value of genomic sequencing efforts for informing management will be greatly enhanced. Broader Impacts: Seagrass restoration and mitigation is currently of major interest in California and elsewhere and the project results will inform current initiatives regarding eelgrass management in California through the state's Ocean Protection Council. In addition to recruiting individual students from diverse backgrounds to work on the project, the project broadens participation of students in STEM fields through its partnership with three existing outreach/training programs at UC Davis.</p> <p>This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.</p> ZosteraEcoGenomics largerWorkCitation project eng; USA oceans -117.225474 7.493273 32.713756 67.267997 2017-08-15 2017-10-30 In Zostera marina beds worldwide, including western and eastern margins of both the Atlantic and Pacific Oceans. Project centered in Bodega Bay, CA 38.31 N; 123.059 W 0 BCO-DMO catalogue of parameters from Sample collection and sequence accession information for Zostera marina whole genome resequencing from specimens collected at 16 geographic locations worldwide in 2017 Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact http://lod.bco-dmo.org/id/dataset-parameter/931661.rdf Name: SampleID Units: unitless Description: <p>Sample identifier</p> http://lod.bco-dmo.org/id/dataset-parameter/931662.rdf Name: BioProject Units: unitless Description: <p>NCBI BioProject number</p> http://lod.bco-dmo.org/id/dataset-parameter/931663.rdf Name: Library_Code Units: unitless Description: <p>Library Code</p> http://lod.bco-dmo.org/id/dataset-parameter/931664.rdf Name: Sample_Description Units: unitless Description: <p>Sample Description</p> http://lod.bco-dmo.org/id/dataset-parameter/931665.rdf Name: Location_Name Units: unitless Description: <p>Location name</p> http://lod.bco-dmo.org/id/dataset-parameter/931666.rdf Name: Lat Units: decimal degrees Description: <p>Site latitude</p> http://lod.bco-dmo.org/id/dataset-parameter/931667.rdf Name: Long Units: decimal degrees Description: <p>Site longitude; negative values = West</p> http://lod.bco-dmo.org/id/dataset-parameter/931668.rdf Name: SRA_accession Units: unitless Description: <p>NCBI SRA Study Accession identifier</p> http://lod.bco-dmo.org/id/dataset-parameter/931669.rdf Name: Library_type Units: unitless Description: <p>Library type</p> http://lod.bco-dmo.org/id/dataset-parameter/931670.rdf Name: Sampled_tissue Units: unitless Description: <p>Description of sampled tissue</p> http://lod.bco-dmo.org/id/dataset-parameter/931671.rdf Name: Instrument Units: unitless Description: <p>Type of sequencing instrument</p> http://lod.bco-dmo.org/id/dataset-parameter/931672.rdf Name: reads Units: unitless Description: <p>read length</p> http://lod.bco-dmo.org/id/dataset-parameter/931673.rdf Name: Number_of_bases_Raw_fastq Units: count Description: <p>Total number of base pairs sequenced in raw files</p> http://lod.bco-dmo.org/id/dataset-parameter/931674.rdf Name: Coverage_Raw_fastq Units: genomes Description: <p>Genome Coverage (times covered) - Raw</p> http://lod.bco-dmo.org/id/dataset-parameter/931675.rdf Name: Number_of_bases_Clean_fastq Units: count Description: <p>Genome Coverage (times covered) - Cleaned</p> http://lod.bco-dmo.org/id/dataset-parameter/931676.rdf Name: Coverage_Clean_fastq Units: genomes Description: <p>Fraction of the genome covered by cleaned sequencing</p> http://lod.bco-dmo.org/id/dataset-parameter/931677.rdf Name: Mapped_pcnt Units: percent Description: <p>Percent of reads mapped to the genome</p> http://lod.bco-dmo.org/id/dataset-parameter/931678.rdf Name: Properly_paired_pcnt Units: percent Description: <p>Percent properly paired</p> GB/NERC/BODC > British Oceanographic Data Centre, Natural Environment Research Council, United Kingdom Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact 44854 https://darchive.mblwhoilibrary.org/server/api/core/bitstreams/b8d8fe78-ba17-4a9a-8984-f2fe87114ddf/content download https://doi.org/10.26008/1912/bco-dmo.924852.1 download onLine dataset &lt;p&gt;We conducted a range-wide sampling collection of 190 &lt;em&gt;Zostera marina&lt;/em&gt; specimens from 16 geographic locations (Yu et al. 2023 Nature Plants). The chosen populations feature a mix of sexual and vegetative reproduction with the exception of mostly vegetative reproduction at the sites PO and NN, apparent through extended clones. Chosen locations were a subset of the Zostera Experimental Network sites that were previously analysed using 24 microsatellite loci. Although a sampling distance of &amp;gt;2 meters was maintained to reduce the likelihood of collecting the same genet/clone twice, this was not always successful and thus provided an estimate of local clonal diversity. Plant tissue was selected from the basal meristematic part of the shoot after peeling away the leaf sheath to minimize epiphytes (bacteria and diatoms), frozen in liquid nitrogen and stored at -80 degrees Celsius (°C) until DNA extraction. Quality control was performed following Joint Genome Institute guidelines (&lt;a href=&quot;https://jgi.doe.gov/wp-content/uploads/2013/11/Genomic-DNA-Sample-QC.pdf&quot; target=&quot;_blank&quot;&gt;https://jgi.doe.gov/wp-content/uploads/2013/11/Genomic-DNA-Sample-QC.pdf&lt;/a&gt;). Plate-based DNA library preparation for Illumina sequencing was performed on the PerkinElmer Sciclone NGS robotic liquid handling system using Kapa Biosystems library preparation kit. About 200 nanograms (ng) of sample DNA was sheared to a length of around 600 base pairs (bp) using a Covaris LE220 focused ultrasonicator. Selected fragments were end-repaired, A-tailed, and ligated with sequencing adaptors containing a unique molecular index barcode. Libraries were quantified using KAPA Biosystems' next-generation sequencing library qPCR-kit on a Roche LightCycler 480 real-time PCR instrument. Quantified libraries were then pooled together and prepared for sequencing on the Illumina HiSeq2500 sequencer using TruSeq SBS sequencing kits (v4) following a 2 × 150 bp indexed run recipe to a targeted depth of approximately 40x coverage. The quality of the raw reads was assessed by FastQC (&lt;a href=&quot;https://www.bioinformatics.babraham.ac.uk/projects/fastqc/&quot; target=&quot;_blank&quot;&gt;https://www.bioinformatics.babraham.ac.uk/projects/fastqc/&lt;/a&gt;) and visualized by MultiQC58. BBDuk (&lt;a href=&quot;https://jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/&quot; target=&quot;_blank&quot;&gt;https://jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/&lt;/a&gt;) was used to remove adapters and for quality filtering, discarding sequence reads (1) with more than one 'N' (maxns = 1), (2) shorter than 50 bp after trimming (minlength = 50), and (3) with average quality &amp;lt;10 after trimming (maq = 10). FastQC and MultiQC were used for a second round of quality check for the clean reads. Sequencing coverage and mapping rate were calculated for each sample (see Yu et al. 2023 for details).&lt;/p&gt; Specified by the Principal Investigator(s) &lt;p&gt;To analyse genetic loci present throughout the global distribution range of eelgrass, we focused on identifying core genes that are present in genomes of all individuals. To do so, each of the 190 ramets were de novo assembled using HipMer (k = 51). To categorize, extract and compare core and variable (shell and cloud) genes, primary transcript sequences (21,483 gene models) from the &lt;em&gt;Z. marina&lt;/em&gt; reference (V3.1; ref. 19) were aligned using BLAT using default parameters to each de novo assembly. Genes were considered present if the transcript aligned with either (1) &amp;gt;60% identity and &amp;gt;60% coverage from a single alignment or (2) &amp;gt;85% identity and &amp;gt;85% coverage split across three or fewer scaffolds. Individual presence–absence-variation calls were combined into a matrix to classify genes into core, cloud and shell categories based on their observation across the population.&lt;/p&gt; &lt;p&gt;The total number of genes considered was 20,100. Because identical genotypes and fragmented, low-quality assemblies can bias and skew presence–absence-variation analyses, only 141 single representatives of clones and ramets with greater than 17,500 genes were kept to ensure that only unique, high-quality assemblies were retained. Genes were classified using discriminant analysis of principal components into cloud, shell and core gene clusters based on their frequency. Core genes were the largest category, with 18,717 genes that were on average observed in 97% of ramets.&lt;/p&gt; &lt;p&gt;SNP mapping, calling and filtering The quality-filtered reads were mapped against the chromosome-level &lt;em&gt;Z. marina&lt;/em&gt; reference genome V3.1 using BWA MEM62. The alignments were converted to BAM format and sorted using Samtools. The Mark- Duplicates module in GATK4 was used to identify and tag duplicate reads in the BAM files. The mapping rate for each genotype was calculated using Samtools (Supplementary Data Table 2 of Yu et al. 2023). HaplotypeCaller (GATK4) was used to generate a Genomic Variant Call Format (GVCF) file for each sample, and all the GVCF files were combined by CombineGVCFs (GATK4). GenotypeGVCFs (GATK4) was used to call genetic variants. BCFtools64 was used to remove SNPs within 20 base pairs of an indel or other variant type, as these variant types may cause erroneous SNPs calls. VariantsToTable (GATK4) was used to extract INFO annotations. SNPs meeting one or more than one of the following criteria were marked by VariantFiltration (GATK4): MQ &amp;lt; 40.0; FS &amp;gt; 60.0; QD &amp;lt; 10.0; MQRandSum &amp;gt; 2.5 or MQRandSum &amp;lt; −2.5; Read- PosRandSum &amp;lt; −2.5; ReadPosRandSum &amp;gt; 2.5; SOR &amp;gt; 3.0; DP &amp;gt; 10,804.0 (2 × average DP). Those SNPs were excluded by SelectVariants (GATK4). A total of 3,975,407 SNPs were retained. VCFtools65 was used to convert individual genotypes to missing data when GQ &amp;lt; 30 or DP &amp;lt; 10. Individual homozygous reference calls with one or more reads supporting the variant allele, and individual homozygous variant calls with ≥1 read supporting the reference, were set as missing data. Only bi-allelic SNPs were kept (3,892,668 SNPs). To avoid the reference-genome-related biases, due to the large Pacific–Atlantic genomic divergence, we focused on the 18,717 core genes that were on average observed in 97% of ramets. Bedtools66 was used to find overlap between the SNPs and the core genes, and only those SNPs were kept (ZM_HQ_SNPs, 763,580 SNPs). Genotypes that were outside our custom quality criteria were represented as missing data.&lt;/p&gt; &lt;p&gt;Possible parent-descendant pairs under selfing as well as clonemates were detected based on the shared heterozygosity. To ensure that all genotypes assessed originated by random mating, ten ramets showing evidence for selfing were excluded. Seventeen multiple sampled clonemates were also excluded. After the exclusion of 37 samples owing to missing data, selfing, or clonality, 153 samples were left for further analyses.&lt;/p&gt; &lt;p&gt;The chloroplast genome was de novo assembled by NOVOPlasty. The chloroplast genome of &lt;em&gt;Z. marina&lt;/em&gt; was represented by a circular molecule of 143,968 bp with a classic quadripartite structure: two identical inverted repeats (IRa and IRb) of 24,127 bp each, a large single-copy region of 83,312 bp, and a small single-copy region of 12,402 bp. All regions were equally taken into SNP calling analysis except for 9,818 bp encoding 23S and 16S ribosomal RNAs due to bacterial contamination in some samples. The raw Illumina reads of each individual were aligned by BWA MEM to the assembled chloroplast genome. The alignments were converted to BAM format and then sorted using Samtools. Genomic sites were called as variable positions when the frequency of variant reads was &amp;gt;50% and the total coverage of the position was &amp;gt;30% of the median coverage (174 variable positions). Then 11 positions likely related to microsatellites and 12 positions reflecting minute inversions caused by hairpin structures were removed from the final set of variable positions for the haplotype reconstruction (151 SNPs). For the phylogenetic tree reconstruction, we further selected 108 SNPs that represent parsimony-informative sites (that is, no singletons).&lt;/p&gt; Specified by the Principal Investigator(s) - Created a lookup table to match SRA accession number to BioProject by running script: https://gist.github.com/adyork/07ae4060e3017ff8cddd141dbe2ffe3a#file-sra_to_bioproject_all-ipynb. - Saved this lookup table as &quot;srp_and_bioproject.csv&quot;. - Imported original file &quot;Stachowicz_AccessionNumbers.csv&quot; into the BCO-DMO system. - Imported file &quot;srp_and_bioproject.csv&quot;. - Joined the two files together, adding the BioProject numbers to the primary data table, matching on the SRA accession number. - Renamed fields to comply with BCO-DMO naming conventions. - Saved the final file as &quot;924852_v1_z_marina_whole_genome_resequencing.csv&quot;. Specified by BCO-DMO Data Managers asNeeded 7.x-1.1 Biological and Chemical Oceanography Data Management Office (BCO-DMO) Unavailable 508-289-2009 WHOI MS#36 Woods Hole MA 02543 USA info@bco-dmo.org http://www.bco-dmo.org Monday - Friday 8:00am - 5:00pm For questions regarding this resource, please contact BCO-DMO via the email address provided. pointOfContact Illumina HiSeq2500 sequencer Illumina HiSeq2500 sequencer PI Supplied Instrument Name: Illumina HiSeq2500 sequencer Instrument Name: Automated DNA Sequencer Instrument Short Name:Automated Sequencer Instrument Description: A DNA sequencer is an instrument that determines the order of deoxynucleotides in deoxyribonucleic acid sequences.