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

The chromosome-level genome of the submerged plant Cryptocoryne crispatula provides insights into the terrestrial-freshwater transition in Araceae

Plant terrestrialization (i.e. the transition to a terrestrial environment) is a significant evolutionary event that has been intensively studied. While certain plant lineages, particularly in angiosperms, have re-adapted to freshwater habitats after colonizing terrene, however, the molecular mechanism of the terrestrial-freshwater (T-F) transition remains limited. Here, the basal monocot Araceae was selected as the study object to explore the T-F transition adaptation mechanism by comparative genomic analysis. Our findings revealed that the substitution rates significantly increased in the lineage of freshwater Araceae, which may promote their adaptation to the freshwater habitat. Additionally, 20 gene sets across all four freshwater species displayed signs of positive selection contributing to tissue development and defense responses in freshwater plants. Comparative synteny analysis showed that genes specific to submerged plants were enriched in cellular respiration and photosynthesis. In contrast, floating plants were involved in regulating gene expression, suggesting that gene and genome duplications may provide the original material for plants to adapt to the freshwater environment. Our study provides valuable insights into the genomic aspects of the transition from terrestrial to aquatic environments in Araceae, laying the groundwork for future research in the angiosperm.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

A trait-based framework for seagrass ecology: Trends and prospects

In the last three decades, quantitative approaches that rely on organism traits instead of taxonomy have advanced different fields of ecological research through establishing the mechanistic links between environmental drivers, functional traits, and ecosystem functions. A research subfield where trait-based approaches have been frequently used but poorly synthesized is the ecology of seagrasses; marine angiosperms that colonized the ocean 100M YA and today make up productive yet threatened coastal ecosystems globally. Here, we compiled a comprehensive trait-based response-effect framework (TBF) which builds on previous concepts and ideas, including the use of traits for the study of community assembly processes, from dispersal and response to abiotic and biotic factors, to ecosystem function and service provision. We then apply this framework to the global seagrass literature, using a systematic review to identify the strengths, gaps, and opportunities of the field. Seagrass trait research has mostly focused on the effect of environmental drivers on traits, i.e., "environmental filtering" (72%), whereas links between traits and functions are less common (26.9%). Despite the richness of trait-based data available, concepts related to TBFs are rare in the seagrass literature (15% of studies), including the relative importance of neutral and niche assembly processes, or the influence of trait dominance or complementarity in ecosystem function provision. These knowledge gaps indicate ample potential for further research, highlighting the need to understand the links between the unique traits of seagrasses and the ecosystem services they provide.

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.

Comparative analysis of mitochondrial genomes reveals marine adaptation in seagrasses

BACKGROUND

Seagrasses are higher marine flowering plants that evolved from terrestrial plants, but returned to the sea during the early evolution of monocotyledons through several separate lineages. Thus, they become a good model for studying the adaptation of plants to the marine environment. Sequencing of the mitochondrial (mt) genome of seagrasses is essential for understanding their evolutionary characteristics.

RESULTS

In this study, we sequenced the mt genome of two endangered seagrasses (Zostera japonica and Phyllospadix iwatensis). These data and data on previously sequenced mt genomes from monocotyledons provide new evolutionary evidence of genome size reduction, gene loss, and adaptive evolution in seagrasses. The mt genomes of Z. japonica and P. iwatensis are circular. The sizes of the three seagrasses (including Zostera marine) that have been sequenced to date are smaller than that of other monocotyledons. Additionally, we found a large number of repeat sequences in seagrasses. The most abundant long repeat sequences were 31-40 bp repeats. Our study also found that seagrass species lost extensive ribosomal protein genes during evolution. The rps7 gene and the rpl16 gene of P. iwatensis are exceptions to this trend. The phylogenetic analysis based on the mt genome strongly supports the previous results. Furthermore, we identified five positive selection genes (atp8, nad3, nad6, ccmFn, and matR) in seagrasses that may be associated with their adaptation to the marine environment.

CONCLUSIONS

In this study, we sequenced and annotated the mt genomes of Z. japonica and P. iwatensis and compared them with the genome of other monocotyledons. The results of this study will enhance our understanding of seagrass adaptation to the marine environment and can inform further investigations of the seagrass mt genome.

De novo assembly and annotation of the transcriptome of the endangered seagrass Zostera capensis: Insights from differential gene expression under thermal stress

Seagrasses are important marine ecosystem engineers but anthropogenic impacts and climate change have led to numerous population declines globally. In South Africa, Zostera capensis is endangered due to fragmented populations and heavy anthropogenic pressures on estuarine ecosystems that house the core of the populations. Addressing questions of how pressures such as climate change affect foundational species, including Z. capensis are crucial to supporting their conservation and underpin restoration efforts. Here we use ecological transcriptomics to study key functional responses of Z. capensis through quantification of gene expression after thermal stress and present the first reference transcriptome of Z. capensis. Four de novo reference assemblies (Trinity, IDBA-tran, RNAspades, SOAPdenovo) filtered through the EvidentialGene pipeline resulted in 153,755 transcripts with a BUSCO score of 66.1% for completeness. Differential expression analysis between heat stressed (32 °C for three days) and pre-warming plants identified genes involved in photosynthesis, oxidative stress, translation, metabolic and biosynthetic processes in the Z. capensis thermal stress response. This reference transcriptome is a significant contribution to the limited available genomic resources for Z. capensis and represents a vital tool for addressing questions around the species restoration and potential functional responses to warming marine environments.

Tissue-specific transcriptome profiles identify functional differences key to understanding whole plant response to life in variable salinity

Plants endure environmental stressors via adaptation and phenotypic plasticity. Studying these mechanisms in seagrasses is extremely relevant as they are important primary producers and functionally significant carbon sinks. These mechanisms are not well understood at the tissue level in seagrasses. Using RNA-seq, we generated transcriptome sequences from tissue of leaf, basal leaf meristem and root organs of Posidonia australis, establishing baseline in situ transcriptomic profiles for tissues across a salinity gradient. Samples were collected from four P. australis meadows growing in Shark Bay, Western Australia. Analysis of gene expression showed significant differences between tissue types, with more variation among leaves than meristem or roots. Gene ontology enrichment analysis showed the differences were largely due to the role of photosynthesis, plant growth and nutrient absorption in leaf, meristem and root organs, respectively. Differential gene expression of leaf and meristem showed upregulation of salinity regulation processes in higher salinity meadows. Our study highlights the importance of considering leaf meristem tissue when evaluating whole-plant responses to environmental change.

This article has an associated First Person interview with the first author of the paper.

Summary: Differences in seagrass leaf, meristem and root transcriptomes across variable salinities are due to tissue-specific processes. Leaf meristem contained the broadest process range, indicating preferential use for inferring plant-wide activity.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptations to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptations, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the three main clades that was best explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light- or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and WGDs using phylogenomics.

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.

Ecological and evolutionary diversification of sulphated polysaccharides in diverse photosynthetic lineages: A review

Sulphated polysaccharides (SPs) are carbohydrate macromolecules with sulphate esters that are found among marine algae, seagrasses, mangroves and some terrestrial plants. The sulphate concentration in the ocean (28 mM) since ancient time could have driven the production of SPs in marine algae. SPs have a gelatinous property that can protect marine algae against desiccation and salinity stress. Agar and carrageenan are red algal SPs that are widely used as gelling agents in the food and pharmaceutical industries. The information on the SPs from freshwater and land plants are limited. In this review, we reviewed the taxonomic distribution and composition of SPs in different photosynthetic lineages, and explored the association of SP production in these diversified photosynthetic organisms with evolution history and environmental stresses. We also reviewed the genes/proteins involved in SP biosynthesis. Insights into SP biosynthetic machinery may shed light on the evolution that accompanied adaptation to life on earth.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass)

Background: Seagrasses (Alismatales) are the only fully marine angiosperms.  Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass  Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings.  

Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. 

Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. 

Conclusions: As an important marine angiosperm, the improved  Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Seagrasses are marine flowering plants found in tropical and sub-tropical areas that live in coastal regions between the sea and land. All seagrass species evolved from terrestrial monocotyledons, providing the opportunity to study plant adaptation to sea environments. Here, we sequenced the chloroplast genomes (cpGenomes) of three Zostera species, then analyzed and compared their cpGenome structures and sequence variations. We also performed a phylogenetic analysis using published seagrass chloroplasts and calculated the selection pressure of 17 species within seagrasses and nine terrestrial monocotyledons, as well as estimated the number of shared genes of eight seagrasses. The cpGenomes of Zosteraceae species ranged in size from 143,877 bp (Zostera marina) to 152,726 bp (Phyllospadix iwatensis), which were conserved and displayed similar structures and gene orders. Additionally, we found 17 variable hotspot regions as candidate DNA barcodes for Zosteraceae species, which will be helpful for studying the phylogenetic relationships and interspecies differences between seagrass species. Interestingly, nine genes had positive selection sites, including two ATP subunit genes (atpA and atpF), two ribosome subunit genes (rps4 and rpl20), two DNA-dependent RNA polymerase genes (rpoC1 and rpoC2), as well as accD, clpP, and ycf2. These gene regions may have played key roles in the seagrass adaptation to diverse environments. The Branch model analysis showed that seagrasses had a higher rate of evolution than terrestrial monocotyledons, suggesting that seagrasses experienced greater environmental pressure. Moreover, a branch-site model identified positively selected sites (PSSs) in ccsA, suggesting their involvement in the adaptation to sea environments. These findings are valuable for further investigations on Zosteraceae cpGenomes and will serve as an excellent resource for future studies on seagrass adaptation to sea environments.

Gene body DNA methylation in seagrasses: inter- and intraspecific differences and interaction with transcriptome plasticity under heat stress

The role of DNA methylation and its interaction with gene expression and transcriptome plasticity is poorly understood, and current insight comes mainly from studies in very few model plant species. Here, we study gene body DNA methylation (gbM) and gene expression patterns in ecotypes from contrasting thermal environments of two marine plants with contrasting life history strategies in order to explore the potential role epigenetic mechanisms could play in gene plasticity and responsiveness to heat stress. In silico transcriptome analysis of CpGO/E ratios suggested that the bulk of Posidonia oceanica and Cymodocea nodosa genes possess high levels of intragenic methylation. We also observed a correlation between gbM and gene expression flexibility: genes with low DNA methylation tend to show flexible gene expression and plasticity under changing conditions. Furthermore, the empirical determination of global DNA methylation (5-mC) showed patterns of intra and inter-specific divergence that suggests a link between methylation level and the plants' latitude of origin and life history. Although we cannot discern whether gbM regulates gene expression or vice versa, or if other molecular mechanisms play a role in facilitating transcriptome responsiveness, our findings point to the existence of a relationship between gene responsiveness and gbM patterns in marine plants.

Seagrasses in an era of ocean warming: a review

Seagrasses are valuable sources of food and habitat for marine life and are one of Earth's most efficient carbon sinks. However, they are facing a global decline due to ocean warming and eutrophication. In the last decade, with the advent of new technology and molecular advances, there has been a dramatic increase in the number of studies focusing on the effects of ocean warming on seagrasses. Here, we provide a comprehensive review of the future of seagrasses in an era of ocean warming. We have gathered information from published studies to identify potential commonalities in the effects of warming and the responses of seagrasses across four distinct levels: molecular, biochemical/physiological, morphological/population, and ecosystem/planetary. To date, we know that although warming strongly affects seagrasses at all four levels, seagrass responses diverge amongst species, populations, and over depths. Furthermore, warming alters seagrass distribution causing massive die-offs in some seagrass populations, whilst also causing tropicalization and migration of temperate species. In this review, we evaluate the combined effects of ocean warming with other environmental stressors and emphasize the need for multiple-stressor studies to provide a deeper understanding of seagrass resilience. We conclude by discussing the most significant knowledge gaps and future directions for seagrass research.

Phenotypic plasticity under rapid global changes: The intrinsic force for future seagrasses survival

Coastal oceans are particularly affected by rapid and extreme environmental changes with dramatic consequences for the entire ecosystem. Seagrasses are key ecosystem engineering or foundation species supporting diverse and productive ecosystems along the coastline that are particularly susceptible to fast environmental changes. In this context, the analysis of phenotypic plasticity could reveal important insights into seagrasses persistence, as it represents an individual property that allows species' phenotypes to accommodate and react to fast environmental changes and stress. Many studies have provided different definitions of plasticity and related processes (acclimation and adaptation) resulting in a variety of associated terminology. Here, we review different ways to define phenotypic plasticity with particular reference to seagrass responses to single and multiple stressors. We relate plasticity to the shape of reaction norms, resulting from genotype by environment interactions, and examine its role in the presence of environmental shifts. The potential role of genetic and epigenetic changes in underlying seagrasses plasticity in face of environmental changes is also discussed. Different approaches aimed to assess local acclimation and adaptation in seagrasses are explored, explaining strengths and weaknesses based on the main results obtained from the most recent literature. We conclude that the implemented experimental approaches, whether performed with controlled or field experiments, provide new insights to explore the basis of plasticity in seagrasses. However, an improvement of molecular analysis and the application of multi-factorial experiments are required to better explore genetic and epigenetic adjustments to rapid environmental shifts. These considerations revealed the potential for selecting the best phenotypes to promote assisted evolution with fundamental implications on restoration and preservation efforts.

The Genetic Component of Seagrass Restoration: What We Know and the Way Forwards

Seagrasses are marine flowering plants providing key ecological services and functions in coasts and estuaries across the globe. Increased environmental changes fueled by human activities are affecting their existence, compromising natural habitats and ecosystems’ biodiversity and functioning. In this context, restoration of disturbed seagrass environments has become a worldwide priority to reverse ecosystem degradation and to recover ecosystem functionality and associated services. Despite the proven importance of genetic research to perform successful restoration projects, this aspect has often been overlooked in seagrass restoration. Here, we aimed to provide a comprehensive perspective of genetic aspects related to seagrass restoration. To this end, we first reviewed the importance of studying the genetic diversity and population structure of target seagrass populations; then, we discussed the pros and cons of different approaches used to restore and/or reinforce degraded populations. In general, the collection of genetic information and the development of connectivity maps are critical steps for any seagrass restoration activity. Traditionally, the selection of donor population preferred the use of local gene pools, thought to be the best adapted to current conditions. However, in the face of rapid ocean changes, alternative approaches such as the use of climate-adjusted or admixture genotypes might provide more sustainable options to secure the survival of restored meadows. Also, we discussed different transplantation strategies applied in seagrasses and emphasized the importance of long-term seagrass monitoring in restoration. The newly developed information on epigenetics as well as the application of assisted evolution strategies were also explored. Finally, a view of legal and ethical issues related to national and international restoration management is included, highlighting improvements and potential new directions to integrate with the genetic assessment. We concluded that a good restoration effort should incorporate: (1) a good understanding of the genetic structure of both donors and populations being restored; (2) the analysis of local environmental conditions and disturbances that affect the site to be restored; (3) the analysis of local adaptation constraints influencing the performances of donor populations and native plants; (4) the integration of distribution/connectivity maps with genetic information and environmental factors relative to the target seagrass populations; (5) the planning of long-term monitoring programs to assess the performance of the restored populations. The inclusion of epigenetic knowledge and the development of assisted evolution programs are strongly hoped for the future.

Evolutionary Origins of Drought Tolerance in Spermatophytes

It is commonly known that drought stress is a major constraint limiting crop production. Drought stress and associated drought tolerance mechanisms are therefore under intense investigation with the view to future production of drought tolerant crops. With an ever-growing population and variable climate, novel approaches need to be considered to sustainably feed future generations. In this context, definitions of drought tolerance are highly variable, which poses a major challenge for the systematic assessment of this trait across the plant kingdom. Furthermore, drought tolerance is a polygenic trait and understanding the evolution of this complex trait may inform us about patterns of gene gain and loss in relation to diverse drought adaptations. We look at the transition of plants from water to land, and the role of drought tolerance in enabling this transition, before discussing the first drought tolerant plant and common drought responses amongst vascular plants. We reviewed the distribution of a combined "drought tolerance" trait in very broad terms to encompass different experimental systems and definitions used in the current literature and assigned a binary trait "tolerance vs. sensitivity" in 178 extant plant species. By simplifying drought responses of plants into this "binary" trait we were able to explore the evolution of drought tolerance across the wider plant kingdom, compared to previous studies. We show how this binary "drought tolerance/sensitivity" trait has evolved and discuss how incorporating this information into an evolutionary genomics framework could provide insights into the molecular mechanisms underlying extreme drought adaptations.

The Cell Wall of Seagrasses: Fascinating, Peculiar and a Blank Canvas for Future Research

Seegrasses are a polyphyletic group of angiosperm plants, which evolved from early monocotyledonous land plants and returned to the marine environment around 140 million years ago. Today, seagrasses comprise the five families Zosteraceae, Hydrocharitaceae, Posidoniaceae, Cymodoceaceae, and Ruppiaceae and form important coastal ecosystems worldwide. Despite of this ecological importance, the existing literature on adaption of these angiosperms to the marine environment and especially their cell wall composition is limited up to now. A unique feature described for some seagrasses is the occurrence of polyanionic, low-methylated pectins mainly composed of galacturonic acid and apiose (apiogalacturonans). Furthermore, sulfated galactans have been detected in some species. Recently, arabinogalactan-proteins (AGPs), highly glycosylated proteins of the cell wall of land plants, have been isolated for the first time from a seagrass of the baltic sea. Obviously, seagrass cell walls are characterized by new combinations of structural polysaccharide and glycoprotein elements known from macroalgae and angiosperm land plants. In this review, current knowledge on cell walls of seagrasses is summarized and suggestions for future investigations are given.

m6A RNA Methylation in Marine Plants: First Insights and Relevance for Biological Rhythms

Circadian regulations are essential for enabling organisms to synchronize physiology with environmental light-dark cycles. Post-transcriptional RNA modifications still represent an understudied level of gene expression regulation in plants, although they could play crucial roles in environmental adaptation. N⁶-methyl-adenosine (m⁶A) is the most prevalent mRNA modification, established by "writer" and "eraser" proteins. It influences the clockwork in several taxa, but only few studies have been conducted in plants and none in marine plants. Here, we provided a first inventory of m⁶A-related genes in seagrasses and investigated daily changes in the global RNA methylation and transcript levels of writers and erasers in Cymodocea nodosa and Zostera marina. Both species showed methylation peaks during the dark period under the same photoperiod, despite exhibiting asynchronous changes in the m⁶A profile and related gene expression during a 24-h cycle. At contrasting latitudes, Z. marina populations displayed overlapping daily patterns of the m⁶A level and related gene expression. The observed rhythms are characteristic for each species and similar in populations of the same species with different photoperiods, suggesting the existence of an endogenous circadian control. Globally, our results indicate that m⁶A RNA methylation could widely contribute to circadian regulation in seagrasses, potentially affecting the photo-biological behaviour of these plants.

Arabinogalactan-proteins of Zostera marina L. contain unique glycan structures and provide insight into adaption processes to saline environments

Seagrasses evolved from monocotyledonous land plants that returned to the marine habitat. This transition was accomplished by substantial changes in cell wall composition, revealing habitat-driven adaption to the new environment. Whether arabinogalactan-proteins (AGPs), important signalling molecules of land plants, are present in seagrass cell walls is of evolutionary and plant development interest. AGPs of Zostera marina L. were isolated and structurally characterised by analytical and bioinformatics methods as well as by ELISA with different anti-AGP antibodies. Calcium-binding capacity of AGPs was studied by isothermal titration calorimetry (ITC) and microscopy. Bioinformatic searches of the Z. marina proteome identified 9 classical AGPs and a large number of chimeric AGPs. The glycan structures exhibit unique features, including a high degree of branching and an unusually high content of terminating 4-O-methyl-glucuronic acid (4-OMe GlcA) residues. Although the common backbone structure of land plant AGPs is conserved in Z. marina, the terminating residues are distinct with high amounts of uronic acids. These differences likely result from the glycan-active enzymes (glycosyltransferases and methyltransferases) and are essential for calcium-binding properties. The role of this polyanionic surface is discussed with regard to adaption to the marine environment.

Transcriptional Landscapes of Divergent Sporophyte Development in Two Mosses, Physcomitrium (Physcomitrella) patens and Funaria hygrometrica

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

SMRT sequencing of full-length transcriptome of seagrasses Zostera japonica

Seagrass meadows are among the four most productive marine ecosystems in the world. Zostera japonica (Z. japonica) is the most widely distributed species of seagrass in China. However, there is no reference genome or transcriptome available for Z. japonica, impeding progress in functional genomic and molecular ecology studies in this species. Temperature is the main factor that controls the distribution and growth of seagrass around the world, yet how seagrass responds to heat stress remains poorly understood due to the lack of genomic and transcriptomic data. In this study, we applied a combination of second- and third-generation sequencing technologies to sequence full-length transcriptomes of Z. japonica. In total, we obtained 58,134 uniform transcripts, which included 46,070 high-quality full-length transcript sequences. We identified 15,411 simple sequence repeats, 258 long non-coding RNAs and 28,038 open reading frames. Exposure to heat elicited a complex transcriptional response in genes involved in posttranslational modification, protein turnover and chaperones. Overall, our study provides the first large-scale full-length trascriptome in Zostera japonica, allowing for structural, functional and comparative genomics studies in this important seagrass species. Although previous studies have focused specifically on heat shock proteins, we found that examination of other heat stress related genes is important for studying response to heat stress in seagrass. This study provides a genetic resource for the discovery of genes related to heat stress tolerance in this species. Our transcriptome can be further utilized in future studies to understand the molecular adaptation to heat stress in Zostera japonica.

De novo Assembly of the Pokeweed Genome Provides Insight Into Pokeweed Antiviral Protein (PAP) Gene Expression

Ribosome-inactivating proteins (RIPs) are RNA glycosidases thought to function in defense against pathogens. These enzymes remove purine bases from RNAs, including rRNA; the latter activity decreases protein synthesis in vitro, which is hypothesized to limit pathogen proliferation by causing host cell death. Pokeweed antiviral protein (PAP) is a RIP synthesized by the American pokeweed plant (Phytolacca americana). PAP inhibits virus infection when expressed in crop plants, yet little is known about the function of PAP in pokeweed due to a lack of genomic tools for this non-model species. In this work, we de novo assembled the pokeweed genome and annotated protein-coding genes. Sequencing comprised paired-end reads from a short-insert library of 83X coverage, and our draft assembly (N50 = 42.5 Kb) accounted for 74% of the measured pokeweed genome size of 1.3 Gb. We obtained 29,773 genes, 73% of which contained known protein domains, and identified several PAP isoforms. Within the gene models of each PAP isoform, a long 5' UTR intron was discovered, which was validated by RT-PCR and sequencing. Presence of the intron stimulated reporter gene expression in tobacco. To gain further understanding of PAP regulation, we complemented this genomic resource with expression profiles of pokeweed plants subjected to stress treatments [jasmonic acid (JA), salicylic acid, polyethylene glycol, and wounding]. Cluster analysis of the top differentially expressed genes indicated that some PAP isoforms shared expression patterns with genes involved in terpenoid biosynthesis, JA-mediated signaling, and metabolism of amino acids and carbohydrates. The newly sequenced promoters of all PAP isoforms contained cis-regulatory elements associated with diverse biotic and abiotic stresses. These elements mediated response to JA in tobacco, based on reporter constructs containing promoter truncations of PAP-I, the most abundant isoform. Taken together, this first genomic resource for the Phytolaccaceae plant family provides new insight into the regulation and function of PAP in pokeweed.

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.

Rapid Metabolome and Bioactivity Profiling of Fungi Associated with the Leaf and Rhizosphere of the Baltic Seagrass Zostera marina

Zostera marina (eelgrass) is a marine foundation species with key ecological roles in coastal habitats. Its bacterial microbiota has been well studied, but very little is known about its mycobiome. In this study, we have isolated and identified 13 fungal strains, dominated by Penicillium species (10 strains), from the leaf and the root rhizosphere of Baltic Z. marina. The organic extracts of the fungi that were cultured by an OSMAC (One-Strain–Many-Compounds) regime using five liquid culture media under both static and shaking conditions were investigated for their chemical and bioactivity profiles. All extracts showed strong anti-quorum sensing activity, and the majority of the Penicillium extracts displayed antimicrobial or anti-biofilm activity against Gram-negative environmental marine and human pathogens. HPLC-DAD-MS-based rapid metabolome analyses of the extracts indicated the high influence of culture conditions on the secondary metabolite (SM) profiles. Among 69 compounds detected in all Penicillium sp. extracts, 46 were successfully dereplicated. Analysis of SM relatedness in culture conditions by Hierarchical Cluster Analysis (HCA) revealed generally low similarity and showed a strong effect of medium selection on chemical profiles of Penicillium sp. This is the first study assessing both the metabolite and bioactivity profile of the fungi associated with Baltic eelgrass Z. marina.

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

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

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Use of a draft genome of coffee (Coffea arabica) to identify SNPs associated with caffeine content

Arabica coffee (Coffea arabica) has a small gene pool limiting genetic improvement. Selection for caffeine content within this gene pool would be assisted by identification of the genes controlling this important trait. Sequencing of DNA bulks from 18 genotypes with extreme high- or low-caffeine content from a population of 232 genotypes was used to identify linked polymorphisms. To obtain a reference genome, a whole genome assembly of arabica coffee (variety K7) was achieved by sequencing using short read (Illumina) and long-read (PacBio) technology. Assembly was performed using a range of assembly tools resulting in 76 409 scaffolds with a scaffold N50 of 54 544 bp and a total scaffold length of 1448 Mb. Validation of the genome assembly using different tools showed high completeness of the genome. More than 99% of transcriptome sequences mapped to the C. arabica draft genome, and 89% of BUSCOs were present. The assembled genome annotated using AUGUSTUS yielded 99 829 gene models. Using the draft arabica genome as reference in mapping and variant calling allowed the detection of 1444 nonsynonymous single nucleotide polymorphisms (SNPs) associated with caffeine content. Based on Kyoto Encyclopaedia of Genes and Genomes pathway-based analysis, 65 caffeine-associated SNPs were discovered, among which 11 SNPs were associated with genes encoding enzymes involved in the conversion of substrates, which participate in the caffeine biosynthesis pathways. This analysis demonstrated the complex genetic control of this key trait in coffee.

The Tropical Invasive Seagrass, Halophila stipulacea, Has a Superior Ability to Tolerate Dynamic Changes in Salinity Levels Compared to Its Freshwater Relative, Vallisneria americana

The tropical seagrass species, Halophila stipulacea, originated from the Indian Ocean and the Red Sea, subsequently invading the Mediterranean and has recently established itself in the Caribbean Sea. Due to its invasive nature, there is growing interest in understanding this species' capacity to adapt to new conditions. One approach to understanding the natural tolerance of a plant is to compare the tolerant species with a closely related non-tolerant species. We compared the physiological responses of H. stipulacea exposed to different salinities, with that of its nearest freshwater relative, Vallisneria americana. To achieve this goal, H. stipulacea and V. americana plants were grown in dedicated microcosms, and exposed to the following salt regimes: (i) H. stipulacea: control (40 PSU, practical salinity units), hyposalinity (25 PSU) and hypersalinity (60 PSU) for 3 weeks followed by a 4-week recovery phase (back to 40 PSU); (ii) V. americana: control (1 PSU), and hypersalinity (12 PSU) for 3 weeks, followed by a 4-week recovery phase (back to 1 PSU). In H. stipulacea, leaf number and chlorophyll content showed no significant differences between control plants and plants under hypo and hypersalinities, but a significant decrease in leaf area under hypersalinity was observed. In addition, compared with control plants, H. stipulacea plants exposed to hypo and hypersalinity were found to have reduced below-ground biomass and C/N ratios, suggesting changes in the allocation of resources in response to both stresses. There was no significant effect of hypo/hypersalinity on dark-adapted quantum yield of photosystem II (Fv/Fm) suggesting that H. stipulacea photochemistry is resilient to hypo/hypersalinity stress. In contrast to the seagrass, V. americana exposed to hypersalinity displayed significant decreases in above-ground biomass, shoot number, leaf number, blade length and Fv/Fm, followed by significant recoveries of all these parameters upon return of the plants to non-saline control conditions. These data suggest that H. stipulacea shows remarkable tolerance to both hypo and hypersalinity. Resilience to a relatively wide range of salinities may be one of the traits explaining the invasive nature of this species in the Mediterranean and Caribbean Seas.

Highly Disturbed Populations of Seagrass Show Increased Resilience but Lower Genotypic Diversity

The response of seagrass systems to a severe disturbance provides an opportunity to quantify the degree of resilience in different meadows, and subsequently to test whether there is a genetic basis to resilience. We used existing data on levels of long-standing disturbance from poor water quality, and the responses of seagrass (Zostera muelleri) after an extreme flood event in Moreton Bay, Queensland, Australia. Sites were grouped into high and low disturbance categories, in which seagrass showed high and low resilience, respectively, as determined by measuring rates of key feedback processes (nutrient removal, suppression of sediment resuspension, and algal grazing), and physiological and morphological traits. Theoretically, meadows with higher genotypic diversity would be expected to have greater resilience. However, because the more resilient meadows occur in areas historically exposed to high disturbance, the alternative is also possible, that selection will have resulted in a narrower, less diverse subset of genotypes than in less disturbed meadows. Levels of genotypic and genetic diversity (allelic richness) based on 11 microsatellite loci, were positively related (R² = 0.58). Genotypic diversity was significantly lower at highly disturbed sites (R = 0.49) than at less disturbed sites (R = 0.61). Genotypic diversity also showed a negative trend with two morphological characteristics known to confer resilience on seagrass in Moreton Bay, leaf chlorophyll concentrations and seagrass biomass. Genetic diversity did not differ between disturbed and undisturbed sites. We postulate that the explanation for these results is historical selection for genotypes that confer protection against disturbance, reducing diversity in meadows that contemporarily show greater resilience.

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.

The Huperzia selago Shoot Tip Transcriptome Sheds New Light on the Evolution of Leaves

Lycopodiophyta—consisting of three orders, Lycopodiales, Isoetales and Selaginellales, with different types of shoot apical meristems (SAMs)—form the earliest branch among the extant vascular plants. They represent a sister group to all other vascular plants, from which they differ in that their leaves are microphylls—that is, leaves with a single, unbranched vein, emerging from the protostele without a leaf gap—not megaphylls. All leaves represent determinate organs originating on the flanks of indeterminate SAMs. Thus, leaf formation requires the suppression of indeterminacy, that is, of KNOX transcription factors. In seed plants, this is mediated by different groups of transcription factors including ARP and YABBY.

We generated a shoot tip transcriptome of Huperzia selago (Lycopodiales) to examine the genes involved in leaf formation. Our H. selago transcriptome does not contain any ARP homolog, although transcriptomes of Selaginella spp. do. Surprisingly, we discovered a YABBY homolog, although these transcription factors were assumed to have evolved only in seed plants.

The existence of a YABBY homolog in H. selago suggests that YABBY evolved already in the common ancestor of the vascular plants, and subsequently was lost in some lineages like Selaginellales, whereas ARP may have been lost in Lycopodiales. The presence of YABBY in the common ancestor of vascular plants would also support the hypothesis that this common ancestor had a simplex SAM. Furthermore, a comparison of the expression patterns of ARP in shoot tips of Selaginella kraussiana (Harrison CJ, etal. 2005. Independent recruitment of a conserved developmental mechanism during leaf evolution. Nature 434(7032):509–514.) and YABBY in shoot tips of H. selago implies that the development of microphylls, unlike megaphylls, does not seem to depend on the combined activities of ARP and YABBY. Altogether, our data show that Lycopodiophyta are a diverse group; so, in order to understand the role of Lycopodiophyta in evolution, representatives of Lycopodiales, Selaginellales, as well as of Isoetales, have to be examined.

Identification of Two Novel Amalgaviruses in the Common Eelgrass (Zostera marina) and in Silico Analysis of the Amalgavirus +1 Programmed Ribosomal Frameshifting Sites

The genome sequences of two novel monopartite RNA viruses were identified in a common eelgrass (Zostera marina) transcriptome dataset. Sequence comparison and phylogenetic analyses revealed that these two novel viruses belong to the genus Amalgavirus in the family Amalgaviridae. They were named Zostera marina amalgavirus 1 (ZmAV1) and Zostera marina amalgavirus 2 (ZmAV2). Genomes of both ZmAV1 and ZmAV2 contain two overlapping open reading frames (ORFs). ORF1 encodes a putative replication factory matrix-like protein, while ORF2 encodes a RNA-dependent RNA polymerase (RdRp) domain. The fusion protein (ORF1+2) of ORF1 and ORF2, which mediates RNA replication, was produced using the +1 programmed ribosomal frameshifting (PRF) mechanism. The +1 PRF motif sequence, UUU_CGN, which is highly conserved among known amalgaviruses, was also found in ZmAV1 and ZmAV2. Multiple sequence alignment of the ORF1+2 fusion proteins from 24 amalgaviruses revealed that +1 PRF occurred only at three different positions within the 13-amino acid-long segment, which was surrounded by highly conserved regions on both sides. This suggested that the +1 PRF may be constrained by the structure of fusion proteins. Genome sequences of ZmAV1 and ZmAV2, which are the first viruses to be identified in common eelgrass, will serve as useful resources for studying evolution and diversity of amalgaviruses.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

A new mechanistic understanding of light-limitation in the seagrass Zostera muelleri

In this study we investigated the effect of light-limitation (∼20 μmol photons m^-2 s^-1) on the southern hemisphere seagrass, Zostera muelleri. RNA sequencing, chlorophyll fluorometry and HPLC techniques were used to investigate how the leaf-specific transcriptome drives changes in photosynthesis and photo-pigments in Z. muelleri over 6 days. 1593 (7.51%) genes were differentially expressed on day 2 and 1481 (6.98%) genes were differentially expressed on day 6 of the experiment. Differential gene expression correlated with significant decreases in rETR_Max, I_k, an increase in Yi (initial photosynthetic quantum yield of photosystem II), and significant changes in pigment composition. Regulation of carbohydrate metabolism was observed along with evidence that abscisic acid may serve a role in the low-light response of this seagrass. This study provides a novel understanding of how Z. muelleri responds to light-limitation in the marine water column and provides potential molecular markers for future conservation monitoring efforts.

SeagrassDB: An open-source transcriptomics landscape for phylogenetically profiled seagrasses and aquatic plants

Seagrasses and aquatic plants are important clades of higher plants, significant for carbon sequestration and marine ecological restoration. They are valuable in the sense that they allow us to understand how plants have developed traits to adapt to high salinity and photosynthetically challenged environments. Here, we present a large-scale phylogenetically profiled transcriptomics repository covering seagrasses and aquatic plants. SeagrassDB encompasses a total of 1,052,262 unigenes with a minimum and maximum contig length of 8,831 bp and 16,705 bp respectively. SeagrassDB provides access to 34,455 transcription factors, 470,568 PFAM domains, 382,528 prosite models and 482,121 InterPro domains across 9 species. SeagrassDB allows for the comparative gene mining using BLAST-based approaches and subsequent unigenes sequence retrieval with associated features such as expression (FPKM values), gene ontologies, functional assignments, family level classification, Interpro domains, KEGG orthology (KO), transcription factors and prosite information. SeagrassDB is available to the scientific community for exploring the functional genic landscape of seagrass and aquatic plants at: http://115.146.91.129/index.php.

Tissue-specific transcriptomic profiling provides new insights into the reproductive ecology and biology of the iconic seagrass species Posidonia oceanica

Seagrasses form extensive meadows in shallow coastal waters and are among the world's most productive ecosystems. Seagrasses can produce both clonally and sexually, and flowering has long been considered infrequent, but important for maintaining genetically diverse stands. Here we investigate the molecular mechanisms involved in flowering of the seagrass Posidonia oceanica, an iconic species endemic to the Mediterranean. We generated a de novo transcriptome of this non-model species for leaf, male and female flower tissue of three individuals, and present molecular evidence for genes that may be involved in the flowering process and on the reproductive biology of the species. We present evidence that suggests that P. oceanica exhibits a strategy of protogyny, where the female part of the hermaphroditic flower develops before the male part, in order to avoid self-fertilization. We found photosynthetic genes to be up-regulated in the female flower tissues, indicating that this may be capable of photosynthesis. Finally, we detected a number of interesting genes, previously known to be involved in flowering pathways responding to light and temperature cues and in pathways involved in anthocyanin and exine synthesis. This first comparative transcriptomic approach of leaf, male and female tissue provides a basis for functional genomics research on flower development in P. oceanica and other seagrass species.