Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Manipulation of the seagrass-associated microbiome reduces disease severity

Host-associated microbes influence host health and function and can be a first line of defence against infections. While research increasingly shows that terrestrial plant microbiomes contribute to bacterial, fungal, and oomycete disease resistance, no comparable experimental work has investigated marine plant microbiomes or more diverse disease agents. We test the hypothesis that the eelgrass (Zostera marina) leaf microbiome increases resistance to seagrass wasting disease. From field eelgrass with paired diseased and asymptomatic tissue, 16S rRNA gene amplicon sequencing revealed that bacterial composition and richness varied markedly between diseased and asymptomatic tissue in one of the two years. This suggests that the influence of disease on eelgrass microbial communities may vary with environmental conditions. We next experimentally reduced the eelgrass microbiome with antibiotics and bleach, then inoculated plants with Labyrinthula zosterae, the causative agent of wasting disease. We detected significantly higher disease severity in eelgrass with a native microbiome than an experimentally reduced microbiome. Our results over multiple experiments do not support a protective role of the eelgrass microbiome against L. zosterae. Further studies of these marine host-microbe-pathogen relationships may continue to show new relationships between plant microbiomes and diseases.

Combined transcriptome and proteome analysis reveal the key physiological processes in seed germination stimulated by decreased salinity in the seagrass Zostera marina L

BACKGROUND

Zostera marina L., or eelgrass, is the most widespread seagrass species throughout the temperate northern hemisphere. Unlike the dry seeds of terrestrial plants, eelgrass seeds must survive in water, and salinity is the key factor influencing eelgrass seed germination. In the present study, transcriptome and proteome analysis were combined to investigate the mechanisms via which eelgrass seed germination was stimulated by low salinity, in addition to the dynamics of key metabolic pathways under germination.

RESULTS

According to the results, low salinity stimulated the activation of Ca2+ signaling and phosphatidylinositol signaling, and further initiated various germination-related physiological processes through the MAPK transduction cascade. Starch, lipids, and storage proteins were mobilized actively to provide the energy and material basis for germination; abscisic acid synthesis and signal transduction were inhibited whereas gibberellin synthesis and signal transduction were activated, weakening seed dormancy and preparing for germination; cell wall weakening and remodeling processes were activated to provide protection for cotyledon protrusion; in addition, multiple antioxidant systems were activated to alleviate oxidative stress generated during the germination process; ERF transcription factor has the highest number in both stages suggested an active role in eelgrass seed germination.

CONCLUSION

In summary, for the first time, the present study investigated the mechanisms by which eelgrass seed germination was stimulated by low salinity and analyzed the transcriptomic and proteomic features during eelgrass seed germination comprehensively. The results of the present study enhanced our understanding of seagrass seed germination, especially the molecular ecology of seagrass seeds.

Environmental DNA metabarcoding reveals distinct fish assemblages supported by seagrass (Zostera marina and Zostera pacifica) beds in different geographic settings in Southern California

Seagrass beds are disappearing at a record pace despite their known value to our oceans and coastal communities. Simultaneously, our coastlines are under the constant pressure of climate change which is impacting their chemical, physical and biological characteristics. It is thus pertinent to evaluate and record habitat use so we can understand how these different environments contribute to local biodiversity. This study evaluates the assemblages of fish found at five Zostera beds in Southern California using environmental DNA (eDNA) metabarcoding. eDNA is a powerful biodiversity monitoring tool that offers key advantages to conventional monitoring. Results from our eDNA study found 78 species of fish that inhabit these five beds around Southern California representing embayment, open coastal mainland and open coastal island settings. While each bed had the same average number of species found throughout the year, the composition of these fish assemblages was strongly site dependent. There were 35 fish that were found at both open coast and embayment seagrass beds, while embayment seagrass sites had 20 unique fish and open coast sites had 23 unique fish. These results demonstrate that seagrass fish assemblages are heterogenous based on their geographic positioning and that marine managers must take this into account for holistic conservation and restoration efforts.

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

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.

Chloroplast genomic comparison provides insights into the evolution of seagrasses

BACKGROUND

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have evolved to live entirely submerged in marine waters. Thus, these species are ideal for studying plant adaptation to marine environments. Herein, we sequenced the chloroplast (cp) genomes of two seagrass species (Zostera muelleri and Halophila ovalis) and performed a comparative analysis of them with 10 previously published seagrasses, resulting in various novel findings.

RESULTS

The cp genomes of the seagrasses ranged in size from 143,877 bp (Zostera marina) to 178,261 bp (Thalassia hemprichii), and also varied in size among different families in the following order: Hydrocharitaceae > Cymodoceaceae > Ruppiaceae > Zosteraceae. The length differences between families were mainly related to the expansion and contraction of the IR region. In addition, we screened out 2,751 simple sequence repeats and 1,757 long repeat sequence types in the cp genome sequences of the 12 seagrass species, ultimately finding seven hot spots in coding regions. Interestingly, we found nine genes with positive selection sites, including two ATP subunit genes (atpA and atpF), three ribosome subunit genes (rps4, rps7, and rpl20), one photosystem subunit gene (psbH), and the ycf2, accD, and rbcL genes. These gene regions may have played critical roles in the adaptation of seagrasses to diverse environments. In addition, phylogenetic analysis strongly supported the division of the 12 seagrass species into four previously recognized major clades. Finally, the divergence time of the seagrasses inferred from the cp genome sequences was generally consistent with previous studies.

CONCLUSIONS

In this study, we compared chloroplast genomes from 12 seagrass species, covering the main phylogenetic clades. Our findings will provide valuable genetic data for research into the taxonomy, phylogeny, and species evolution of seagrasses.

Overexpression of Zostera japonica heat shock protein gene ZjHsp70 enhances the thermotolerance of transgenic Arabidopsis

BACKGROUND

Heat shock protein 70s (Hsp70s) are major members of the heat shock protein family and play a variety of roles to protect plants against stress. Plant Hsp70s are a conserved and widely expressed family of heat shock proteins. They have two main functional regions: N-terminal nucleic acid binding region and C-terminal substrate binding region.

METHODS AND RESULTS

In this study, we cloned the Hsp70 gene of Zostera japonica (ZjHsp70) based on the sequence obtained by transcriptome sequencing. The transcriptional levels of ZjHsp70 increased significantly at 1 h after heat treatment. ZjHsp70 was located in the cytoplasm and nucleus. The overexpression of ZjHsp70 in Arabidopsis resulted in increased heat tolerance, lower contents of malondialdehyde and higher antioxidant enzyme activity than in the wild type. ZjHsp70 may achieve this goal by maintaining highly active antioxidant enzymes.

CONCLUSIONS

We show that ZjHsp70 can improve plant heat tolerance by maintaining high antioxidant enzyme activity under high temperature stress. This study provided a basis to study the role of ZjHsp70 in thermotolerance in more detail.

Profiling the cell walls of seagrasses from A (Amphibolis) to Z (Zostera)

BACKGROUND

The polyphyletic group of seagrasses shows an evolutionary history from early monocotyledonous land plants to the marine environment. Seagrasses form important coastal ecosystems worldwide and large amounts of seagrass detritus washed on beaches might also be valuable bioeconomical resources. Despite this importance and potential, little is known about adaptation of these angiosperms to the marine environment and their cell walls.

RESULTS

We investigated polysaccharide composition of nine seagrass species from the Mediterranean, Red Sea and eastern Indian Ocean. Sequential extraction revealed a similar seagrass cell wall polysaccharide composition to terrestrial angiosperms: arabinogalactans, pectins and different hemicelluloses, especially xylans and/or xyloglucans. However, the pectic fractions were characterized by the monosaccharide apiose, suggesting unusual apiogalacturonans are a common feature of seagrass cell walls. Detailed analyses of four representative species identified differences between organs and species in their constituent monosaccharide composition and lignin content and structure. Rhizomes were richer in glucosyl units compared to leaves and roots. Enhalus had high apiosyl and arabinosyl abundance, while two Australian species of Amphibolis and Posidonia, were characterized by high amounts of xylosyl residues. Interestingly, the latter two species contained appreciable amounts of lignin, especially in roots and rhizomes whereas Zostera and Enhalus were lignin-free. Lignin structure in Amphibolis was characterized by a higher syringyl content compared to that of Posidonia.

CONCLUSIONS

Our investigations give a first comprehensive overview on cell wall composition across seagrass families, which will help understanding adaptation to a marine environment in the evolutionary context and evaluating the potential of seagrass in biorefinery incentives.

Identification of the genome sequence of Zostera associated varicosavirus 1, a novel negative-sense RNA virus, in the common eelgrass (Zostera marina) transcriptome

Varicosaviruses (the genus Varicosavirus) are bipartite, negative-sense, single-stranded RNA viruses that infect plants. We analyzed a transcriptome dataset isolated from the common eelgrass (Zostera marina) and identified a novel varicosavirus named Zostera associated varicosavirus 1 (ZaVV1). The ZaVV1 genome consists of two genomic segments: RNA1 (6,632-nt) has an open reading frame (ORF) encoding a large multi-functional polymerase protein (L), while RNA2 (4,304-nt) has four ORFs: one for a nucleocapsid protein and three for proteins with unknown functions (P2, P3, and P4). Sequence comparison and phylogenetic analysis using L proteins showed that ZaVV1 is a novel member of the genus Varicosavirus of the family Rhabdoviridae. The conserved regulatory elements involved in transcription termination/polyadenylation and transcription initiation were identified in the ZaVV1 gene-junction regions with the consensus sequence 3'-UAUUAUUCUUUUUGCUCU-5' (in the negative-sense genome). The ZaVV1 genome sequence may be useful for studying the phylogenetic relationships of varicosaviruses and genome evolution of rhabdoviruses. Keywords: Zostera associated varicosavirus 1, Varicosavirus, Rhabdoviridae, common eelgrass, Zostera marina.

New insights into physiological effects of anoxia under darkness on the iconic seagrass Zostera marina based on a combined analysis of transcriptomics and metabolomics

Coastal hypoxia/anoxia is a major emerging threat to global coastal ecosystems. Macroalgae blooms of tens of kilometers are often observed in open waters. These blooms not only cause a lack of oxygen, but also benthic light limitation. We explored the physiological responses of Zostera marina L. to anoxia under darkness. After exposing Z. marina to anoxia under darkness for 72 h, we measured the elongation of leaves and the decrease in maximal quantum yield of photosystem II (Fv/Fm), and investigated the transcriptomic and metabolomic responses to anoxic stress based on RNA-sequencing and liquid chromatography-mass spectrometry (LC-MS) technology. The results showed that anoxic stress significantly reduced the leaf Fv/Fm, and had a significant negative effect on the photosynthesis and growth of Z. marina. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of up-regulated differentially expressed genes (DEGs) showed that glycolysis was the most significant enrichment pathway (p < 0.001), and most of the important products in glycolysis were significantly up-regulated. This indicated that the glycolysis process of anaerobic respiration is promoted under anoxia. The metabolite results also showed that glyceraldehyde 3-phosphate in the glycolysis pathway was significantly up-regulated. Moreover, three genes encoding sucrose synthase (gene-ZOSMA_310G00150, gene-ZOSMA_81G00980, and gene-ZOSMA_8G00730) and one gene encoding alpha-amylase (gene-ZOSMA_95G00270) were significantly up-regulated, providing the sugar basis for the subsequent increase in glycolysis. Furthermore, gene-encoding oxoglutarate dehydrogenase, the rate-limiting step of the tricarboxylic acid (TCA) cycle, was significantly down-regulated, indicating that this cycle was inhibited under anoxia. Metabolomic results showed that L-tryptophan, L-phenylalanine, and DL-leucine were significantly up-regulated. Only significantly decreased glutamate and non-significantly decreased glutamine, substances consumed in alanine and γ-aminobutyric acid (GABA) shunt mechanisms, were detected in the leaves, while GABA and alanine were not detected. The results of this study show that anoxic stress induces a programmed transcriptomic and metabolomic response in seagrass, most likely reflecting a complex strategy of acclimation and adaptation in seagrass to resist anoxic stress.

Heat-stress induced expression of stress-inducible nucleotide exchange factor Fes1 in seagrass Zostera japonica

Seagrass meadows are among the four most productive marine natural ecosystems in the world. Zostera japonica (Z. japonica) is the most widely distributed species of seagrass in China. Nucleotide exchange factors (NEFs) promote the release of ADP during heat stress, accelerating the rate-limiting step of Heat shock protein 70 (Hsp70). Although NEFs play an important role in abiotic stress tolerance of plants, NEFs in seagrass have not been studied. In this study, we cloned Fes1 from Z. japonica (ZjFes1) by rapid amplification of the cDNA ends using RACE, and full length ZjFes1 was 1171 bp. It contained an 81 bp 5'-terminal untranslated region (UTR), 109 bp 3'-UTR and 981 bp open reading frame (ORF). The ORF (ZjFes1) was predicted to encode a polypeptide of 326 amino acids with theoretical molecular weight (MW) of 36.10 kDa and pI of 5.22. ZjFes1 shared 89% amino acid identity with Fes1 from Zostera marina (Z. marina). The transcriptional levels of ZjFes1 increased significantly 1 h after heat treatment. ZjFes1 was localized to the cytoplasm. Taken together, we found that ZjFes1 was a stress-inducible gene that may be involved in heat stress response. This study lays the foundation for further studies on the role of ZjFes1 in heat resistance.

Integrating genetics, biophysical, and demographic insights identifies critical sites for seagrass conservation

The eelgrass Zostera marina is an important foundation species of coastal areas in the Northern Hemisphere, but is continuing to decline, despite management actions. The development of new management tools is therefore urgent in order to prioritize limited resources for protecting meadows most vulnerable to local extinctions and identifying most valuable present and historic meadows to protect and restore, respectively. We assessed 377 eelgrass meadows along the complex coastlines of two fjord regions on the Swedish west coast-one is currently healthy and the other is substantially degraded. Shoot dispersal for all meadows was assessed with Lagrangian biophysical modeling (scale: 100-1,000 m) and used for barrier analysis and clustering; a subset (n = 22) was also assessed with population genetic methods (20 microsatellites) including diversity, structure, and network connectivity. Both approaches were in very good agreement, resulting in seven subpopulation groupings or management units (MUs). The MUs correspond to a spatial scale appropriate for coastal management of "waterbodies" used in the European Water Framework Directive. Adding demographic modeling based on the genetic and biophysical data as a third approach, we are able to assess past, present, and future metapopulation dynamics to identify especially vulnerable and valuable meadows. In a further application, we show how the biophysical approach, using eigenvalue perturbation theory (EPT) and distribution records from the 1980s, can be used to identify lost meadows where restoration would best benefit the present metapopulation. The combination of methods, presented here as a toolbox, allows the assessment of different temporal and spatial scales at the same time, as well as ranking of specific meadows according to key genetic, demographic and ecological metrics. It could be applied to any species or region, and we exemplify its versatility as a management guide for eelgrass along the Swedish west coast.

The highly efficient NDH-dependent photosystem I cyclic electron flow pathway in the marine angiosperm Zostera marina

Zostera marina, a fully submerged marine angiosperm with a unique evolutionary history associated with its terrestrial origin, has distinct photochemical characteristics caused by its oxygen-evolving complex (OEC) being prone to deactivation in visible light. Based on the present phylogenetic analysis, the chloroplast NADPH dehydrogenase-like (NDH) complex was found to be completed in of Z. marina, unlike other marine plants, suggesting its crucial role. Thus, the responses of electron transport to irradiation were investigated through multiple chlorophyll fluorescence techniques and Western blot analysis. Moreover, the respective contribution of the two photosystem I cyclic electron flow (PSI-CEF) pathways to the generation of trans-thylakoid proton gradient (∆pH) was also examined using inhibitors. The contributions of the two PSI-CEF pathways to ∆pH were similar; furthermore, there was a trade-off between the two pathways under excess irradiation: the PGR5/L1-dependent PSI-CEF decreased gradually following its activation during the initial illumination, while NDH-dependent PSI-CEF was activated gradually with exposure duration. OEC inactivation was continuously under excess irradiation, which exhibits a positive linear correlation with the activation of NDH-dependent PSI-CEF. We suggest that PGR5/L1-dependent PSI-CEF was preferentially activated to handle the excess electron caused by the operation of OEC during the initial illumination. Subsequently, the increasing OEC inactivation with exposure duration resulted in a deficit of electrons. Limited electrons from PSI might preferentially synthesize NADPH, which could support the function of NDH-dependent PSI-CEF to generate ∆pH and ATP via reducing ferredoxin, thereby maintaining OEC stability.

Overexpression of Seagrass Nucleotide Exchange Factor Gene ZjFes1 Enhances Heat Tolerance in Transgenic Arabidopsis

Nucleotide exchange factors (NEFs) play an important role in plant abiotic stress response, but their characteristics and functions in seagrass have not been studied. Zostera japonica (Z. japonica) is one of the most widely distributed seagrass species in China and are distributed in subtropical and temperate coastal areas. Z. japonica is intertidal seagrass, which often undergoes heat stress during summer when the tide ebbs. Overexpression of ZjFes1 in Arabidopsis results in an increase in heat tolerance. We found that ZjFes1 associates with ZjHsp70 in vivo by yeast two-hybrid and bimolecular fluorescence complementarity (BiFC). Overexpression of ZjFes1 leads to selective reduction of Hsp70 transcription and an increase in Hsp101. In conclusion, our results suggest that ZjFes1 may be an active regulator of heat tolerance.

Transcriptome profiling analysis of the seagrass, Zostera muelleri under copper stress

Copper (Cu) in an essential trace metal but it can also contaminate coastal waters at high concentrations mainly from agricultural run-off and mining activities which are detrimental to marine organisms including seagrasses. The molecular mechanisms driving Cu toxicity in seagrasses are not clearly understood yet. Here, we investigated the molecular responses of the Australian seagrass, Z. muelleri at the whole transcriptomic level after 7 days of exposure to 250 μg Cu L^-1 and 500 μg Cu L^-1. The leaf-specific whole transcriptome results showed a concentration-dependent disturbance in chloroplast function, regulatory stress responses and defense mechanisms. This study provided new insights into the responses of seagrasses to trace metal stress and reports possible candidate genes which can be considered as biomarkers to improve conservation and management of seagrass meadows.

Effect of reduced irradiance on 13C uptake, gene expression and protein activity of the seagrass Zostera muelleri

Photosynthesis in the seagrass Zostera muelleri remains poorly understood. We investigated the effect of reduced irradiance on the incorporation of ¹³C, gene expression of photosynthetic, photorespiratory and intermediates recycling genes as well as the enzymatic content and activity of Rubisco and PEPC within Z. muelleri. Following 48 h of reduced irradiance, we found that i) there was a ∼7 fold reduction in ¹³C incorporation in above ground tissue, ii) a significant down regulation of photosynthetic, photorespiratory and intermediates recycling genes and iii) no significant difference in enzyme activity and content. We propose that Z. muelleri is able to alter its physiology in order to reduce the amount of C lost through photorespiration to compensate for the reduced carbon assimilation as a result of reduced irradiance. In addition, the first estimated rate constant (K_cat) and maximum rates of carboxylation (V_cmax) of Rubisco is reported for the first time for Z. muelleri.

Molecular Characterization of ZosmaNRT2, the Putative Sodium Dependent High-Affinity Nitrate Transporter of Zostera marina L

One of the most important adaptations of seagrasses during sea colonization was the capacity to grow at the low micromolar nitrate concentrations present in the sea. In contrast to terrestrial plants that use H⁺ symporters for high-affinity NO₃⁻ uptake, seagrasses such as Zostera marina L. use a Na⁺-dependent high-affinity nitrate transporter. Interestingly, in the Z. marina genome, only one gene (Zosma70g00300.1; NRT2.1) is annotated to this function. Analysis of this sequence predicts the presence of 12 transmembrane domains, including the MFS domains of the NNP transporter family and the “nitrate signature” that appears in all members of the NNP family. Phylogenetic analysis shows that this sequence is more related to NRT2.5 than to NRT2.1, sharing a common ancestor with both monocot and dicot plants. Heterologous expression of ZosmaNRT2-GFP together with the high-affinity nitrate transporter accessory protein ZosmaNAR2 (Zosma63g00220.1) in Nicotiana benthamiana leaves displayed four-fold higher fluorescence intensity than single expression of ZosmaNRT2-GFP suggesting the stabilization of NRT2 by NAR2. ZosmaNRT2-GFP signal was present on the Hechtian-strands in the plasmolyzed cells, pointing that ZosmaNRT2 is localized on the plasma membrane and that would be stabilized by ZosmaNAR2. Taken together, these results suggest that Zosma70g00300.1 would encode a high-affinity nitrate transporter located at the plasma membrane, equivalent to NRT2.5 transporters. These molecular data, together with our previous electrophysiological results support that ZosmaNRT2 would have evolved to use Na⁺ as a driving ion, which might be an essential adaptation of seagrasses to colonize marine environments.

Excess copper promotes photoinhibition and modulates the expression of antioxidant-related genes in Zostera muelleri

Copper (Cu) is an essential micronutrient for plants and as such is vital to many metabolic processes. Nevertheless, when present at elevated concentrations, Cu can exert toxic effects on plants by disrupting protein functions and promoting oxidative stress. Due to their proximity to the urbanised estuaries, seagrasses are vulnerable to chemical contamination via industrial runoff, waste discharges and leachates. Zostera muelleri is a common seagrass species that forms habitats in the intertidal areas along the temperate coast of Australia. Previous studies have shown the detrimental effects of Cu exposure on photosynthetic efficiency of Z. muelleri. The present study focuses on the impacts of sublethal Cu exposure on the physiological and molecular responses. By means of a single addition, plants were exposed to 250 and 500 μg Cu L^-1 (corresponding to 3.9 and 7.8 μM, respectively) as well as uncontaminated artificial seawater (control) for 7 days. Chlorophyll fluorescence parameters, measured as the effective quantum yield (ϕPSII), the maximum quantum yield (Fv/Fm) and non-photochemical quenching (NPQ) were assessed daily, while Cu accumulation in leaf tissue, total reactive oxygen species (ROS) and the expression of genes involved in antioxidant activities and trace metal binding were determined after 1, 3 and 7 days of exposure. Z. muelleri accumulated Cu in the leaf tissue in a concentration-dependent manner and the bioaccumulation was saturated by day 3. Cu exposure resulted in an acute suppression of ϕPSII and Fv/Fm. These two parameters also showed a concentration- and time-dependent decline. NPQ increased sharply during the first few days before subsequently decreasing towards the end of the experiment. Cu accumulation induced oxidative stress in Z. muelleri as an elevated level of ROS was detected on day 7. Lower Cu concentration promoted an up-regulation of genes encoding Cu/Zn-superoxide dismutase (Cu/Zn-sod), ascorbate peroxidase (apx), catalase (cat) and glutathione peroxidase (gpx), whereas no significant change was detected with higher Cu concentration. Exposure to Cu at any concentration failed to induce regulation in the expression level of genes encoding metallothionein type 2 (mt2), metallothionein type 3 (mt3) and cytochrome c oxidase copper chaperone (cox17). It is concluded that chlorophyll fluorescence parameters provide timely probe of the status of photosynthetic machinery under Cu stress. In addition, when exposed to a moderate level of Cu, Z. muelleri mitigates any induced oxidative stress by up-regulating transcripts coding for antioxidant enzymes.

Individual-based genetic analyses support asexual hydrochory dispersal in Zostera noltei

Dispersal beyond the local patch in clonal plants was typically thought to result from sexual reproduction via seed dispersal. However, evidence for the separation, transport by water, and re-establishment of asexual propagules (asexual hydrochory) is mounting suggesting other important means of dispersal in aquatic plants. Using an unprecedented sampling size and microsatellite genetic identification, we describe the distribution of seagrass clones along tens of km within a coastal lagoon in Southern Portugal. Our spatially explicit individual-based sampling design covered 84 km2 and collected 3 185 Zostera noltei ramets from 803 sites. We estimated clone age, assuming rhizome elongation as the only mechanism of clone spread, and contrasted it with paleo-oceanographic sea level change. We also studied the association between a source of disturbance and the location of large clones. A total of 16 clones were sampled more than 10 times and the most abundant one was sampled 59 times. The largest distance between two samples from the same clone was 26.4 km and a total of 58 and 10 clones were sampled across more than 2 and 10 km, respectively. The number of extremely large clone sizes, and their old ages when assuming the rhizome elongation as the single causal mechanism, suggests other processes are behind the span of these clones. We discuss how the dispersal of vegetative fragments in a stepping-stone manner might have produced this pattern. We found higher probabilities to sample large clones away from the lagoon inlet, considered a source of disturbance. This study corroborates previous experiments on the success of transport and re-establishment of asexual fragments and supports the hypothesis that asexual hydrochory is responsible for the extent of these clones.

Genomic comparison of two independent seagrass lineages reveals habitat-driven convergent evolution

Seagrasses are marine angiosperms that live fully submerged in the sea. They evolved from land plant ancestors, with multiple species representing at least three independent return-to-the-sea events. This raises the question of whether these marine angiosperms followed the same adaptation pathway to allow them to live and reproduce under the hostile marine conditions. To compare the basis of marine adaptation between seagrass lineages, we generated genomic data for Halophila ovalis and compared this with recently published genomes for two members of Zosteraceae, as well as genomes of five non-marine plant species (Arabidopsis, Oryza sativa, Phoenix dactylifera, Musa acuminata, and Spirodela polyrhiza). Halophila and Zosteraceae represent two independent seagrass lineages separated by around 30 million years. Genes that were lost or conserved in both lineages were identified. All three species lost genes associated with ethylene and terpenoid biosynthesis, and retained genes related to salinity adaptation, such as those for osmoregulation. In contrast, the loss of the NADH dehydrogenase-like complex is unique to H. ovalis. Through comparison of two independent return-to-the-sea events, this study further describes marine adaptation characteristics common to seagrass families, identifies species-specific gene loss, and provides molecular evidence for convergent evolution in seagrass lineages.

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.

Molecular physiology reveals ammonium uptake and related gene expression in the seagrass Zostera muelleri

Seagrasses are important marine foundation species, which are presently threatened by coastal development and global change worldwide. The molecular mechanisms that drive seagrass responses to anthropogenic stresses, including elevated levels of nutrients such as ammonium, remains poorly understood. Despite the evidence that seagrasses can assimilate ammonium by using glutamine synthetase (GS)/glutamate synthase (glutamine-oxoglutarate amidotransferase or GOGAT) cycle, the regulation of this fundamental metabolic pathway has never been studied at the gene expression level in seagrasses so far. Here, we combine (i) reverse transcription quantitative real-time PCR (RT-qPCR) to measure expression of key genes involved in the GS/GOGAT cycle, and (ii) stable isotope labelling and mass spectrometry to investigate ¹⁵N-ammonium assimilation in the widespread Australian species Zostera muelleri subsp. capricorni (Z. muelleri). We demonstrate that exposure to a pulse of ammonium in seawater can induce changes in GS gene expression of Z. muelleri, and further correlate these changes in gene expression with ¹⁵N-ammonium uptake rate in above- and below-ground tissue.

The Genome of a Southern Hemisphere Seagrass Species (Zostera muelleri)

Seagrasses are marine angiosperms that evolved from land plants but returned to the sea around 140 million years ago during the early evolution of monocotyledonous plants. They successfully adapted to abiotic stresses associated with growth in the marine environment, and today, seagrasses are distributed in coastal waters worldwide. Seagrass meadows are an important oceanic carbon sink and provide food and breeding grounds for diverse marine species. Here, we report the assembly and characterization of the Zostera muelleri genome, a southern hemisphere temperate species. Multiple genes were lost or modified in Z. muelleri compared with terrestrial or floating aquatic plants that are associated with their adaptation to life in the ocean. These include genes for hormone biosynthesis and signaling and cell wall catabolism. There is evidence of whole-genome duplication in Z. muelleri; however, an ancient pan-commelinid duplication event is absent, highlighting the early divergence of this species from the main monocot lineages.

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 relative importance of trait vs. genetic differentiation for the outcome of interactions among plant genotypes

Functional trait differences and genetic distance are increasingly used as metrics to predict the. outcome of species interactions and the maintenance of diversity. We apply these ideas to intraspecific diversity for the seagrass Zostera marina (eelgrass), by explicitly testing the influence of trait distance and genetic relatedness on the outcome of pairwise interactions among eelgrass genotypes. Increasing trait distance (but not relatedness) between eelgrass genotypes decreased the likelihood that both would persist over a year-long field experiment, contrary to our expectations based on niche partitioning. In plots in which one genotype excluded another, the biomass and growth of the remaining genotype increased with the trait distance and genetic relatedness of the initial pair, presumably due to a legacy of past interactions. Together these results suggest that sustained competition among functionally similar genotypes did not produce a clear winner, but rapid exclusion occurred among genotypes with distinct trait combinations. Borrowing from coexistence theory, we argue that fitness differences between genotypes with distinct traits overwhelmed any stabilizing effects of niche differentiation. Previously observed effects of eelgrass genetic diversity on performance may rely on nonadditive interactions among multiple genotypes or sufficient environmental heterogeneity to increase stabilizing forces and/or interactions.

Genome-wide survey of the seagrass Zostera muelleri suggests modification of the ethylene signalling network

Seagrasses are flowering plants which grow fully submerged in the marine environment. They have evolved a range of adaptations to environmental challenges including light attenuation through water, the physical stress of wave action and tidal currents, high concentrations of salt, oxygen deficiency in marine sediment, and water-borne pollination. Although, seagrasses are a key stone species of the costal ecosystems, many questions regarding seagrass biology and evolution remain unanswered. Genome sequence data for the widespread Australian seagrass species Zostera muelleri were generated and the unassembled data were compared with the annotated genes of five sequenced plant species (Arabidopsis thaliana, Oryza sativa, Phoenix dactylifera, Musa acuminata, and Spirodela polyrhiza). Genes which are conserved between Z. muelleri and the five plant species were identified, together with genes that have been lost in Z. muelleri. The effect of gene loss on biological processes was assessed on the gene ontology classification level. Gene loss in Z. muelleri appears to influence some core biological processes such as ethylene biosynthesis. This study provides a foundation for further studies of seagrass evolution as well as the hormonal regulation of plant growth and development.

Metabarcoding is powerful yet still blind: a comparative analysis of morphological and molecular surveys of seagrass communities

In the context of the sixth wave of extinction, reliable surveys of biodiversity are increasingly needed to infer the cause and consequences of species and community declines, identify early warning indicators of tipping points, and provide reliable impact assessments before engaging in activities with potential environmental hazards. DNA metabarcoding has emerged as having potential to provide speedy assessment of community structure from environmental samples. Here we tested the reliability of metabarcoding by comparing morphological and molecular inventories of invertebrate communities associated with seagrasses through estimates of alpha and beta diversity, as well as the identification of the most abundant taxa. Sediment samples were collected from six Zostera marina seagrass meadows across Brittany, France. Metabarcoding surveys were performed using both mitochondrial (Cytochrome Oxidase I) and nuclear (small subunit 18S ribosomal RNA) markers, and compared to morphological inventories compiled by a long-term benthic monitoring network. A sampling strategy was defined to enhance performance and accuracy of results by preventing the dominance of larger animals, boosting statistical support through replicates, and using two genes to compensate for taxonomic biases. Molecular barcodes proved powerful by revealing a remarkable level of diversity that vastly exceeded the morphological survey, while both surveys identified congruent differentiation of the meadows. However, despite the addition of individual barcodes of common species into taxonomic reference databases, the retrieval of only 36% of these species suggest that the remaining were either not present in the molecular samples or not detected by the molecular screening. This finding exemplifies the necessity of comprehensive and well-curated taxonomic reference libraries and multi-gene surveys. Overall, results offer methodological guidelines and support for metabarcoding as a powerful and repeatable method of characterizing communities, while also presenting suggestions for improvement, including implementation of pilot studies prior to performing full "blind" metabarcoding assessments to optimize sampling and amplification protocols.

De novo assembly and characterization of the transcriptome of seagrass Zostera marina using Illumina paired-end sequencing

Background

The seagrass Zostera marina is a monocotyledonous angiosperm belonging to a polyphyletic group of plants that can live submerged in marine habitats. Zostera marina L. is one of the most common seagrasses and is considered a cornerstone of marine plant molecular ecology research and comparative studies. However, the mechanisms underlying its adaptation to the marine environment still remain poorly understood due to limited transcriptomic and genomic data.

Principal Findings

Here we explored the transcriptome of Z. marina leaves under different environmental conditions using Illumina paired-end sequencing. Approximately 55 million sequencing reads were obtained, representing 58,457 transcripts that correspond to 24,216 unigenes. A total of 14,389 (59.41%) unigenes were annotated by blast searches against the NCBI non-redundant protein database. 45.18% and 46.91% of the unigenes had significant similarity with proteins in the Swiss-Prot database and Pfam database, respectively. Among these, 13,897 unigenes were assigned to 57 Gene Ontology (GO) terms and 4,745 unigenes were identified and mapped to 233 pathways via functional annotation against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). We compared the orthologous gene family of the Z. marina transcriptome to Oryza sativa and Pyropia yezoensis and 11,667 orthologous gene families are specific to Z. marina. Furthermore, we identified the photoreceptors sensing red/far-red light and blue light. Also, we identified a large number of genes that are involved in ion transporters and channels including Na⁺ efflux, K⁺ uptake, Cl⁻ channels, and H⁺ pumping.

Conclusions

Our study contains an extensive sequencing and gene-annotation analysis of Z. marina. This information represents a genetic resource for the discovery of genes related to light sensing and salt tolerance in this species. Our transcriptome can be further utilized in future studies on molecular adaptation to abiotic stress in Z. marina.

Current European Labyrinthula zosterae are not virulent and modulate seagrass (Zostera marina) defense gene expression

Pro- and eukaryotic microbes associated with multi-cellular organisms are receiving increasing attention as a driving factor in ecosystems. Endophytes in plants can change host performance by altering nutrient uptake, secondary metabolite production or defense mechanisms. Recent studies detected widespread prevalence of Labyrinthula zosterae in European Zostera marina meadows, a protist that allegedly caused a massive amphi-Atlantic seagrass die-off event in the 1930's, while showing only limited virulence today. As a limiting factor for pathogenicity, we investigated genotype × genotype interactions of host and pathogen from different regions (10-100 km-scale) through reciprocal infection. Although the endophyte rapidly infected Z. marina, we found little evidence that Z. marina was negatively impacted by L. zosterae. Instead Z. marina showed enhanced leaf growth and kept endophyte abundance low. Moreover, we found almost no interaction of protist × eelgrass-origin on different parameters of L. zosterae virulence/Z. marina performance, and also no increase in mortality after experimental infection. In a target gene approach, we identified a significant down-regulation in the expression of 6/11 genes from the defense cascade of Z. marina after real-time quantitative PCR, revealing strong immune modulation of the host's defense by a potential parasite for the first time in a marine plant. Nevertheless, one gene involved in phenol synthesis was strongly up-regulated, indicating that Z. marina plants were probably able to control the level of infection. There was no change in expression in a general stress indicator gene (HSP70). Mean L. zosterae abundances decreased below 10% after 16 days of experimental runtime. We conclude that under non-stress conditions L. zosterae infection in the study region is not associated with substantial virulence.

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.

Entangled effects of allelic and clonal (genotypic) richness in the resistance and resilience of experimental populations of the seagrass Zostera noltii to diatom invasion

BACKGROUND

The relationship between species diversity and components of ecosystem stability has been extensively studied, whilst the influence of the genetic component of biodiversity remains poorly understood. Here we manipulated both genotypic and allelic richness of the seagrass Zostera noltii, in order to explore their respective influences on the resistance of the experimental population to stress. Thus far intra-specific diversity was seldom taken into account in management plans, and restoration actions showed very low success. Information is therefore needed to understand the factors affecting resistance and resilience of populations.

RESULTS

Our results show a positive influence of both allelic and genotypic richness on the resistance of meadows to environmental perturbations. They also show that at the low genotypic (i.e. clonal) richness levels used in prior experimental approaches, the effects of genotypic and allelic richness could not be disentangled and allelic richness was a likely hidden treatment explaining at least part of the effects hitherto attributed to genotypic richness.

CONCLUSIONS

Altogether, these results emphasize the need to acknowledge and take into account the interdependency of both genotypic and allelic richness in experimental designs attempting to estimate their importance alone or in combination. A positive influence of allelic richness on resistance to perturbations, and of allelic richness combined with genotypic richness on the recovery (resilience) of the experimental populations is supported by differential mortality. These results, on the key species structuring of one of the most threatened coastal ecosystem worldwide, seagrass meadows, support the need to better take into account the distinct compartments of clonal and genetic diversity in management strategies, and in possible restoration plans in the future.

Genetic diversity enhances restoration success by augmenting ecosystem services

Disturbance and habitat destruction due to human activities is a pervasive problem in near-shore marine ecosystems, and restoration is often used to mitigate losses. A common metric used to evaluate the success of restoration is the return of ecosystem services. Previous research has shown that biodiversity, including genetic diversity, is positively associated with the provision of ecosystem services. We conducted a restoration experiment using sources, techniques, and sites similar to actual large-scale seagrass restoration projects and demonstrated that a small increase in genetic diversity enhanced ecosystem services (invertebrate habitat, increased primary productivity, and nutrient retention). In our experiment, plots with elevated genetic diversity had plants that survived longer, increased in density more quickly, and provided more ecosystem services (invertebrate habitat, increased primary productivity, and nutrient retention). We used the number of alleles per locus as a measure of genetic diversity, which, unlike clonal diversity used in earlier research, can be applied to any organism. Additionally, unlike previous studies where positive impacts of diversity occurred only after a large disturbance, this study assessed the importance of diversity in response to potential environmental stresses (high temperature, low light) along a water-depth gradient. We found a positive impact of diversity along the entire depth gradient. Taken together, these results suggest that ecosystem restoration will significantly benefit from obtaining sources (transplants or seeds) with high genetic diversity and from restoration techniques that can maintain that genetic diversity.

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.

Isolation and characterization of a SEPALLATA-like gene, ZjMADS1, from marine angiosperm Zostera japonica

In flowering plants, floral homeotic MADS-box genes, which constitute a large multigene family, play important roles in the specification of floral organs as defined by the ABCDE model. In this study, a MADS-box gene, ZjMADS1, was isolated and characterized from the marine angiosperm Zostera japonica. The predicted length of the ZjMADS1 protein was 246 amino acids (AA), and the AA sequence was most similar to those of the SEPALLATA (SEP) subfamily, corresponding to E-function genes. Southern blot analysis suggested the presence of two SEP3-like genes in the Z. japonica genome. ZjMADS1 mRNA levels were extremely high in the spadices, regardless of the developmental stage, compared to other organs from the reproductive and vegetative shoots. These results suggest that the ZjMADS1 gene may be involved in spadix development in Z. japonica and act as an E-function gene in floral organ development in marine angiosperms.

A set of microsatellite primers for Zostera japonica (Zosteraceae)

PREMISE OF THE STUDY

Polymorphic microsatellite primers were developed in the seagrass Zostera japonica to investigate genetic variation and to identify clonal structure. •

METHODS AND RESULTS

Thirteen polymorphic loci and 23 monomorphic loci were developed in Z. japonica. Two to 13 alleles per locus were observed at the polymorphic loci across 57 individuals of two Z. japonica populations. The observed and expected heterozygosities within populations ranged from 0.0000 to 1.0000 and from 0.0000 to 0.8542, respectively. •

CONCLUSIONS

Our study showed high-level polymorphism at the polymorphic loci in Z. japonica. These primers would be a powerful tool to study genetic variation, clonal structure, and mating systems.

Identification of seagrasses in the gut of a marine herbivorous fish using DNA barcoding and visual inspection techniques

Traditional visual diet analysis techniques were compared with DNA barcoding in juvenile herbivorous rabbitfish Siganus fuscescens collected in Moreton Bay, Australia, where at least six species of seagrass occur. The intergenic spacer trnH-psbA, suggested as the optimal gene for barcoding angiosperms, was used for the first time to identify the seagrass in fish guts. Four seagrass species and one alga were identified visually from gut contents; however, there was considerable uncertainty in visual identification with 38 of 40 fish having unidentifiable plant fragments in their gut. PCR and single-strand conformational polymorphism (SSCP) were able to discriminate three seagrass families from visually cryptic gut contents. While effective in identifying cryptic gut content to family level, this novel method is likely to be most efficient when paired with visual identification techniques.

An evaluation of small-scale genetic diversity and the mating system in Zostera noltii on an intertidal sandflat in the Wadden Sea

BACKGROUND AND AIMS

The dwarf eelgrass, Zostera noltii, is a predominant inhabitant of soft-bottom intertidal regions along the coasts of northern Europe. It is a monoecious, protogynous angiosperm in which the potential for self-fertilization and inbreeding are high, especially if clone sizes exceed pollen dispersal distances. The aim of the present study was to determine the relationship between mating system and clonal structure, examine the relative roles of geitonogamous selfing and biparental inbreeding, measure pollen availability (multiple paternities) and estimate pollen dispersal.

METHODS

A 100-m(2) plot was established in a large, intertidal Z. noltii meadow on the island of Sylt in the German Wadden Sea. A total of 256 adult shoots was sampled: one from the centre of 100 fixed 1-m(2) quadrats (large scale resolution) and an additional 156 from within eight randomly selected 1-m(2) sub-quadrats (small-scale resolution). DNA was extracted from seeds and leaf tissue of all samples and genotyped with nine microsatellite loci.

KEY RESULTS

Mating system analysis revealed high multilocus and single locus outcrossing rates. Average pollen dispersal distance was nearly the same as the average genet (clone) size. Multiple paternity was common and 20-30 % of mature seeds originated from matings within the plot. Among inbred seeds, most resulted from geitonogamy rather than biparental inbreeding.

CONCLUSIONS

Moderate disturbances intrinsic to the intertidal habitat appear to facilitate seed recruitment by gap formation. Pollen dispersal distances are sufficient to maintain outcrossing and high clonal diversity.

Comparative analysis of expressed sequence tag (EST) libraries in the seagrass Zostera marina subjected to temperature stress

Global warming is associated with increasing stress and mortality on temperate seagrass beds, in particular during periods of high sea surface temperatures during summer months, adding to existing anthropogenic impacts, such as eutrophication and habitat destruction. We compare several expressed sequence tag (EST) in the ecologically important seagrass Zostera marina (eelgrass) to elucidate the molecular genetic basis of adaptation to environmental extremes. We compared the tentative unigene (TUG) frequencies of libraries derived from leaf and meristematic tissue from a control situation with two experimentally imposed temperature stress conditions and found that TUG composition is markedly different among these conditions (all P < 0.0001). Under heat stress, we find that 63 TUGs are differentially expressed (d.e.) at 25 degrees C compared with lower, no-stress condition temperatures (4 degrees C and 17 degrees C). Approximately one-third of d.e. eelgrass genes were characteristic for the stress response of the terrestrial plant model Arabidopsis thaliana. The changes in gene expression suggest complex photosynthetic adjustments among light-harvesting complexes, reaction center subunits of photosystem I and II, and components of the dark reaction. Heat shock encoding proteins and reactive oxygen scavengers also were identified, but their overall frequency was too low to perform statistical tests. In all conditions, the most abundant transcript (3-15%) was a putative metallothionein gene with unknown function. We also find evidence that heat stress may translate to enhanced infection by protists. A total of 210 TUGs contain one or more microsatellites as potential candidates for gene-linked genetic markers. Data are publicly available in a user-friendly database at http://www.uni-muenster.de/Evolution/ebb/Services/zostera .

Eelgrass meadows in the California Channel Islands and adjacent coast reveal a mosaic of two species, evidence for introgression and variable clonality

BACKGROUND AND AIMS

Seagrasses are important facilitator species in shallow, soft-bottom marine environments worldwide and, in many places, are threatened by coastal development and eutrophication. One narrow-leaved species (Zostera marina) and one wide-leaved species, variously designated as Z. marina, Z. pacifica or Z. asiatica, are found off the California Channel Islands and adjacent California-Mexico coast. The aim of the present study was to confirm species identification genetically and to link patterns of genetic diversity, connectivity and hybridization among and within the populations with historical sea levels (Ice Age) or the contemporary environment.

METHODS

Samples (n = 11-100) were collected from 28 sites off five California Channel Islands and six sites off the adjacent coast of southern California and Baja California, Mexico. DNA polymorphisms of the rDNA-ITS (internal transcribed spacer) cistron (nuclear), the matK intron (chloroplast) and nine microsatellite loci (nuclear) were examined in a population genetic and phylogeographic context.

KEY RESULTS

All wide-leaved individuals were Z. pacifica, whereas narrow-leaved forms were Z. marina. Microsatellite genotypes were consistent with hybridization between the two species in three populations. The present distribution of Z. pacifica follows a glacial age land mass rather than present oceanographic regimes, but no link was observed between the present distribution of Z. marina and past or present environments. Island populations of Z. marina often were clonal and characterized by low genotypic diversity compared with populations along the Baja California coast. The high level of clonal connectivity around Santa Catalina Island indicated the importance of dispersal and subsequent re-establishment of vegetative fragments.

CONCLUSIONS

The pristine environmental conditions of offshore islands do not guarantee maximum genetic diversity. Future restoration and transplantation efforts of seagrasses must recognize cryptic species and consider the degree of both genetic and genotypic variation in candidate donor populations.

Direct and indirect effects of the herbicides Glyphosate, Bentazone and MCPA on eelgrass (Zostera marina)

Eelgrass beds are important habitats for many organisms, but there has been a decline in the area covered by eelgrass during the last decades due to increased eutrophication resulting in increased shading from phytoplankton. The use of herbicides in terrestrial agriculture has also increased over the last century, and while the effects of herbicides on non-target organisms have been well studied in freshwater they are overlooked in coastal waters. It is not known if herbicides have any effect on the distribution of eelgrass (Zostera marina), or how natural phytoplankton communities respond to the same herbicides. Direct and indirect effects of the herbicides Glyphosate, Bentazone and MCPA both as single toxicants and as mixtures, on the eelgrass plants were investigated in this study. The direct effects on eelgrass were examined by measuring the four different endpoints; the relatively growth rate as length and weight, the chlorophyll a and b ratio, as well as the RNA-DNA ratio, at the end of a 3 days exposure period. The indirect effect was investigated by measuring the effect on a natural phytoplankton production from Roskilde Fjord six times during 1 year. The results showed that the chlorophyll a-b and RNA-DNA ratios were the most sensitive endpoints in single herbicide experiments. The effects of herbicide mixtures on eelgrass were much larger compared to the single substances. Nearly a halving was found for both the relatively growth rate in length and weight, and the RNA-DNA as well as the chlorophyll a-b ratios were also significant reduced. This indicates a possible synergistic effect, and calculations based on the concentration addition model indicate that the low concentrations mixture has a synergistic effect, whereas the high concentration mixture has an antagonistic effect on eelgrass (Z. marina). The low concentrations mixture is the one with the highest relevance for coastal areas. The effect on phytoplankton showed some variation over the year but mainly in a stimulation in primary production. This means that herbicides affect eelgrass both directly and indirectly.

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.

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

The Baja California peninsula represents a biogeographical boundary contributing to regional differentiation among populations of marine animals. We investigated the genetic characteristics of perennial and annual populations of the marine angiosperm, Zostera marina, along the Pacific coast of Baja California and in the Gulf of California, respectively. Populations of Z. marina from five coastal lagoons along the Pacific coast and four sites in the Gulf of California were studied using nine microsatellite loci. Analyses of variance revealed significant interregional differentiation, but no subregional differentiation. Significant spatial differentiation, assessed using theta(ST) values, was observed among all populations within the two regions. Z. marina populations along the Pacific coast are separated by more than 220 km and had the greatest theta(ST) (0.13-0.28) values, suggesting restricted gene flow. In contrast, lower but still significant genetic differentiation was observed among populations within the Gulf of California (theta(ST) = 0.04-0.18), even though populations are separated by more than 250 km. This suggests higher levels of gene flow among Gulf of California populations relative to Pacific coast populations. Direction of gene flow was predominantly southward among Pacific coast populations, whereas no dominant polarity in the Gulf of California populations was observed. The test for isolation by distance (IBD) showed a significant correlation between genetic and geographical distances in Gulf of California populations, but not in Pacific coast populations, perhaps because of shifts in currents during El Niño Southern Oscillation (ENSO) events along the Pacific coast.

Genetic neighbourhood of clone structures in eelgrass meadows quantified by spatial autocorrelation of microsatellite markers

Limited dispersal distances in plant populations frequently cause local genetic structure, which can be quantified by spatial autocorrelation. In clonal plants, three levels of spatial organization can contribute to positive autocorrelation; namely, the neighbourhood of (a) ramets, (b) clone fragments and (c) entire clones. Here we use data from an exhaustive sampling scheme on a clonal plant to measure the contribution of the neighbourhoods of each distinct clonal structure to total spatial autocorrelation. Four plots (256 grid points each) within dense meadows of the marine clonal plant Zostera marina (eelgrass) were sampled for clone structure with nine microsatellite markers ( approximately 80 alleles). We found significant coancestry (f(ij)), at all three levels of spatial organization, even when not allowing for joins between samples of identical genets. In addition, absolute values of f(ij) and the maximum distance with significant positive f(ij) decreased with the progressive exclusion of joins between alike genotypes. The neighbourhood of this clonal plant thus consists of three levels of organization, which are reflected in different kinship structures. Each of these kinship structures may affect the level of biparental inbreeding and the physical distance between flowering shoots and their outcrossing neighbourhood. These results also emphasize the notion that spatial autocorrelation crucially depends on the scale and intensity of sampling.

Genetic diversity enhances the resistance of a seagrass ecosystem to disturbance

Motivated by recent global reductions in biodiversity, empirical and theoretical research suggests that more species-rich systems exhibit enhanced productivity, nutrient cycling, or resistance to disturbance or invasion relative to systems with fewer species. In contrast, few data are available to assess the potential ecosystem-level importance of genetic diversity within species known to play a major functional role. Using a manipulative field experiment, we show that increasing genotypic diversity in a habitat-forming species (the seagrass Zostera marina) enhances community resistance to disturbance by grazing geese. The time required for recovery to near predisturbance densities also decreases with increasing eelgrass genotypic diversity. However, there is no effect of diversity on resilience, measured as the rate of shoot recovery after the disturbance, suggesting that more rapid recovery in diverse plots is due solely to differences in disturbance resistance. Genotypic diversity did not affect ecosystem processes in the absence of disturbance. Thus, our results suggest that genetic diversity, like species diversity, may be most important for enhancing the consistency and reliability of ecosystems by providing biological insurance against environmental change.

Flexible mating: cross-pollination affects sex-expression in a marine clonal plant

Many functionally hermaphroditic plants have evolved mechanisms to reduce interference between the sex functions and to optimize reproductive output. In addition to physical mechanisms such as the spatial (herkogamy) and temporal (dichogamy) separation of male and female functions, plasticity in sex expression by means of mate-recognition (flexible mating) could be important in plants with variable access to cross-pollen. This applies particularly to clonal plants because of their modular growth form. We experimentally tested for the effects of pollen source and vegetative neighbourhood on instantaneous sex ratio and seed production in the self-compatible clonal marine angiosperm Zostera marina L. To this end, we exposed the (monoecious) flowering shoots to self and cross-pollen and to neighbourhoods of their own and a mix of foreign vegetative shoots. Flowering shoots that had been exposed to cross-pollen showed (1) a significantly lower female/male ratio at peak flowering, evidence for mate-recognition, and (2) a significantly higher seed set by the end of the season. Both effects were independent of the genetic composition of their vegetative neighbourhood. The results suggest that Z. marina maintains a cryptic self-incompatibility system not previously described for angiosperms with sub-aqueous pollination. In Z. marina, and possibly other self-compatible clonal plant species, mate-recognition could be a means of increasing the out-crossing probability for flowering shoots with central positions within their clone.

Phylogenetic relationships in the genera Zostera and Heterozostera (Zosteraceae) based on matK sequence data

Phylogenetic analysis of the plastid (chloroplast) DNA matK gene of Zosteraceae species was undertaken. A molecular phylogenetic tree based on matK sequence data showed the monophyly of Heterozostera tasmanica and subgenus Zosterella and did not support the separation of Heterozostera from the genus Zostera. The tree based on matK supported the monophyly of the subgenus Zostera, and showed that Zosteraceae consist of three main groups: Phyllospadix, which is clearly defined by being dioecious; the subgenus Zosterella and Heterozostera; and the subgenus Zostera. Character-state reconstruction of chromosome number and geographic distribution for our molecular phylogenetic tree showed that 2 n=12 is a plesiomorphic character for Zostera and Heterozostera, that the chromosome number was doubled or tripled in two lineages, and that the initial speciation of Zostera and Heterozostera occurred in the Northern Hemisphere. The matK tree showed the close affinity of Z. noltii and Z. japonica, which have disjunct distributions. Zostera marina, which is the only widely distributed species in the subgenus Zostera, also occurring in the northern Atlantic, was shown to be embedded within other subgenus Zostera species.

Inbreeding depression influences genet size distribution in a marine angiosperm

Although inbreeding depression is a major genetic phenomena influencing individual fitness, it is difficult to measure in wild populations. An alternative approach is to correlate heterozygosity, measured using highly polymorphic markers, with a fitness-correlated trait. In clonal plants, genet size is predicted to be fitness correlated. Here we test the prediction that the genet size distribution of the marine clonal plant Zostera marina (eelgrass) is influenced by inbreeding depression. We used nine polymorphic microsatellite markers to access the fine scale clonal structure and to measure individual heterozygosity within 4 plots (each corresponds to 256 m2, sampled at 1-m intervals) in two populations along the German Baltic Coast. The same plots were also sampled for flowering and vegetative shoots to obtain estimates for sexual reproductive output at the level of the genetic individual. We found substantial differences in the genet size distribution between the two populations that may be explained by different disturbance frequency. In both populations, clone size was significantly positively correlated with the total number of flowering shoots, indicating that larger clones have a higher reproductive output. Individual heterozygosity was significantly positively associated with clone size. The effect was much stronger in Falkenstein (low disturbance) than in Maasholm (high disturbance). The results indicate that in a low disturbance population the relatively outbred clones occupy a higher proportion of the available space, possibly because they outcompete relatively inbred neighbours.

Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: Implications for the origin and diversification of seagrasses in Japanese waters

Seagrasses are composed of four families belonging to angiosperms and they are thought to become adaptive to aquatic life independently. Zosteraceae is one such family and because of the relatively high species diversity around Japan and Korea coast areas, the family might have arisen therefrom. To elucidate the origin and evolution of Zosteraceae which consists of three genera, Phyllospadix, Zostera, and Heterozostera, 2.8 kb nucleotide sequences of rbcL and matK genes in the chloroplast genome were examined for various species, including cosmopolitan Z. marina and endemic Z. caulescens. The phylogenetic analysis reveals the following three features. First, based on the synonymous nucleotide substitution rate of the rice chloroplast genome, we estimated the divergence times between Zosteraceae and its closest relative, Potamogetonaceae, and between different genera, Zostera and Phyllospadix, as approximately 100 million years (myr) and 36 myr, respectively, suggesting that Zosteraceae emerged somewhere in the period from 36 myr ago to 100 myr ago. Second, two subgenera of Zostera, Zostera and Zosterella, exhibit their reciprocal monophyly and appear to have differentiated from each other approximately 33 myr ago. However, the third genus Heterozostera branched off only 5 myr ago from the stem lineage leading to Zosterella and this seems too recent in comparison with the ancient divergence of the two subgenera. Third, we estimated the most recent common ancestor of subgenus Zostera as 6 myr. In Z. marina four haplotypes were found in the sample and have diversified in the past 1.5 myr. One haplotype is shared by both sides of the Japan Archipelago and its closely related haplotypes occur also in eastern Pacific Ocean. Based on these phylogeographic analyses, we propose a provisional age related classification of Zosteraceae to argue the origin and evolution.