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

NO enhances the adaptability to high-salt environments by regulating osmotic balance, antioxidant defense, and ion homeostasis in eelgrass based on transcriptome and metabolome analysis

Introduction

Eelgrass is a typical marine angiosperm that exhibits strong adaptability to high-salt environments. Previous studies have shown that various growth and physiological indicators were significantly affected after the nitrate reductase (NR) pathway for nitric oxide (NO) synthesis in eelgrass was blocked.

Methods

To analyze the molecular mechanism of NO on the adaptability to high-salt environment in eelgrass, we treated eelgrass with artificial seawater (control group) and artificial seawater with 1 mM/L Na2WO4 (experimental group). Based on transcriptomics and metabolomics, we explored the molecular mechanism of NO affecting the salt tolerance of eelgrass.

Results

We obtained 326, 368, and 859 differentially expressed genes (DEGs) by transcriptome sequencing in eelgrass roots, stems, and leaves, respectively. Meanwhile, we obtained 63, 52, and 36 differentially accumulated metabolites (DAMs) by metabolomics in roots, stems, and leaves, respectively. Finally, through the combined analysis of transcriptome and metabolome, we found that the NO regulatory mechanism of roots and leaves of eelgrass is similar to that of terrestrial plants, while the regulatory mechanism of stems has similar and unique features.

Discussion

NO in eelgrass roots regulates osmotic balance and antioxidant defense by affecting genes in transmembrane transport and jasmonic acid-related pathways to improve the adaptability of eelgrass to high-salt environments. NO in eelgrass leaves regulates the downstream antioxidant defense system by affecting the signal transduction of plant hormones. NO in the stems of eelgrass regulates ion homeostasis by affecting genes related to ion homeostasis to enhance the adaptability of eelgrass to high-salt environments. Differently, after the NO synthesis was inhibited, the glyoxylate and dicarboxylate metabolism, as well as the tricarboxylic acid (TCA) cycle, was regulated by glucose metabolism as a complementary effect to cope with the high-salt environment in the stems of eelgrass. These are studies on the regulatory mechanism of NO in eelgrass, providing a theoretical basis for the study of the salt tolerance mechanism of marine plants and the improvement of terrestrial crop traits. The key genes discovered in this study can be applied to increase salt tolerance in terrestrial crops through cloning and molecular breeding methods in the future.

Evolution of ethylene as an abiotic stress hormone in streptophytes

All land plants modulate their growth and physiology through intricate signaling cascades. The majority of these are at least modulated-and often triggered-by phytohormones. Over the past decade, it has become apparent that some phytohormones have an evolutionary origin that runs deeper than plant terrestrialization-many emerged in the streptophyte algal progenitors of land plants. Ethylene is such a case. Here we synthesize the current knowledge on the evolution of the phytohormone ethylene and speculate about its deeply conserved role in adjusting stress responses of streptophytes for more than half a billion years of evolution.

Advances in understanding multilevel responses of seagrasses to hypersalinity

Human- and nature-induced hypersaline conditions in coastal systems can lead to profound alterations of the structure and vitality of seagrass meadows and their socio-ecological benefits. In the last two decades, recent research efforts (>50 publications) have contributed significantly to unravel the physiological basis underlying the seagrass-hypersalinity interactions, although most (∼70%) are limited to few species (e.g. Posidonia oceanica, Zostera marina, Thalassia testudinum, Cymodocea nodosa). Variables related to photosynthesis and carbon metabolism are among the most prevalent in the literature, although other key metabolic processes such as plant water relations and responses at molecular (i.e. gene expression) and ultrastructure level are attracting attention. This review emphasises all these latest insights, offering an integrative perspective on the interplay among biological responses across different functional levels (from molecular to clonal structure), and their interaction with biotic/abiotic factors including those related to climate change. Other issues such as the role of salinity in driving the evolutionary trajectory of seagrasses, their acclimation mechanisms to withstand salinity increases or even the adaptive properties of populations that have historically lived under hypersaline conditions are also included. The pivotal role of the costs and limits of phenotypic plasticity in the successful acclimation of marine plants to hypersalinity is also discussed. Finally, some lines of research are proposed to fill the remaining knowledge gaps.

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.

Ion homeostasis for salinity tolerance in plants: a molecular approach

Soil salinity is one of the major environmental stresses faced by the plants. Sodium chloride is the most important salt responsible for inducing salt stress by disrupting the osmotic potential. Due to various innate mechanisms, plants adapt to the sodic niche around them. Genes and transcription factors regulating ion transport and exclusion such as salt overly sensitive (SOS), Na⁺ /H⁺ exchangers (NHXs), high sodium affinity transporter (HKT) and plasma membrane protein (PMP) are activated during salinity stress and help in alleviating cells of ion toxicity. For salt tolerance in plants signal transduction and gene expression is regulated via transcription factors such as NAM (no apical meristem), ATAF (Arabidopsis transcription activation factor), CUC (cup-shaped cotyledon), Apetala 2/ethylene responsive factor (AP2/ERF), W-box binding factor (WRKY) and basic leucine zipper domain (bZIP). Cross-talk between all these transcription factors and genes aid in developing the tolerance mechanisms adopted by plants against salt stress. These genes and transcription factors regulate the movement of ions out of the cells by opening various membrane ion channels. Mutants or knockouts of all these genes are known to be less salt-tolerant compared to wild-types. Using novel molecular techniques such as analysis of genome, transcriptome, ionome and metabolome of a plant, can help in expanding the understanding of salt tolerance mechanism in plants. In this review, we discuss the genes responsible for imparting salt tolerance under salinity stress through transport dynamics of ion balance and need to integrate high-throughput molecular biology techniques to delineate the issue.

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.

Transcriptome sequencing revealed molecular mechanisms underlying tolerance of Suaeda salsa to saline stress

The halophyte Suaeda salsa displayed strong resistance to salinity. Up to date, molecular mechanisms underlying tolerance of S. salsa to salinity have not been well understood. In the present study, S. salsa seedlings were treated with 30‰ salinity and then leaves and roots were subjected to Illumina sequencing. Compared with the control, 68,599 and 77,250 unigenes were significantly differentially expressed in leaves and roots in saline treatment, respectively. KEGG enrichment analyses indicated that photosynthesis process, carbohydrate, lipid and amino acid metabolisms were all downregulated in saline treatment, which should inhibit growth of S. salsa. Expression levels of Na+/H+ exchanger, V-H+ ATPase, choline monooxygenase, potassium and chloride channels were upregulated in saline treatment, which could relieve reduce over-accumulation of Na+ and Cl-. Fe-SOD, glutathione, L-ascorbate and flavonoids function as antioxidants in plants. Genes in relation to them were all upregulated, suggesting that S. salsa initiated various antioxidant mechanisms to tolerate high salinity. Besides, plant hormones, especially auxin, ethylene and jasmonic acid signaling transduction pathways were all upregulated in response to saline treatment, which were important to gene regulations of ion transportation and antioxidation. These changes might comprehensively contribute to tolerance of S. salsa to salinity. Overall, the present study provided new insights to understand the mechanisms underlying tolerance to salinity in halophytes.

Transcriptomic analysis reveals the molecular adaptation to NaCl stress in Zostera marina L

The seagrass Zostera marina L. shows optimal growth in marine water and reduced growth under low salinity conditions. However, little is known about the molecular mechanisms underlying its adaptation to high salinity in Z. marina. In this study, transcriptomic analyses were performed using RNA-seq of the following two groups with different NaCl content: the CK group (seagrasses grown in the absence of NaCl) and the NaCl group (seagrasses grown in the presence of 400 mM NaCl for 6 h). Approximately 316 million high-quality reads were generated, and 87.9% of the data were mapped to the reference genome. Moreover, differentially expressed genes between the CK and NaCl groups were identified. According to a functional analysis, the up-regulated genes after the NaCl treatment were significantly enriched in nitrogen metabolism, calcium signalling and DNA replication while the down-regulated genes were significantly enriched in photosynthesis. A comparative transcriptomic analysis detected many differentially expressed genes and pathways required for adaptation to NaCl stress, providing a foundation for future studies investigating the molecular mechanisms of salt adaptation in Z. marina. We discuss how molecular changes in these processes may have contributed to the NaCl adaptation.

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.

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.

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.

Genomewide transcriptional reprogramming in the seagrass Cymodocea nodosa under experimental ocean acidification

Here, we report the first use of massive-scale RNA-sequencing to explore seagrass response to CO₂ -driven ocean acidification (OA). Large-scale gene expression changes in the seagrass Cymodocea nodosa occurred at CO₂ levels projected by the end of the century. C. nodosa transcriptome was obtained using Illumina RNA-Seq technology and de novo assembly, and differential gene expression was explored in plants exposed to short-term high CO₂ /low pH conditions. At high pCO₂ , there was a significant increased expression of transcripts associated with photosynthesis, including light reaction functions and CO₂ fixation, and also to respiratory pathways, specifically for enzymes involved in glycolysis, in the tricarboxylic acid cycle and in the energy metabolism of the mitochondrial electron transport. The upregulation of respiratory metabolism is probably supported by the increased availability of photosynthates and increased energy demand for biosynthesis and stress-related processes under elevated CO₂ and low pH. The upregulation of several chaperones resembling heat stress-induced changes in gene expression highlighted the positive role these proteins play in tolerance to intracellular acid stress in seagrasses. OA further modifies C. nodosa secondary metabolism inducing the transcription of enzymes related to biosynthesis of carbon-based secondary compounds, in particular the synthesis of polyphenols and isoprenoid compounds that have a variety of biological functions including plant defence. By demonstrating which physiological processes are most sensitive to OA, this research provides a major advance in the understanding of seagrass metabolism in the context of altered seawater chemistry from global climate change.

Carbon-concentrating mechanisms in seagrasses

Seagrasses are unique angiosperms that carry out growth and reproduction submerged in seawater. They occur in at least three families of the Alismatales. All have chloroplasts mainly in the cells of the epidermis. Living in seawater, the supply of inorganic carbon (Ci) to the chloroplasts is diffusion limited, especially under unstirred conditions. Therefore, the supply of CO2 and bicarbonate across the diffusive boundary layer on the outer side of the epidermis is often a limiting factor. Here we discuss the evidence for mechanisms that enhance the uptake of Ci into the epidermal cells. Since bicarbonate is plentiful in seawater, a bicarbonate pump might be expected; however, the evidence for such a pump is not strongly supported. There is evidence for a carbonic anhydrase outside the outer plasmalemma. This, together with evidence for an outward proton pump, suggests the possibility that local acidification leads to enhanced concentrations of CO2 adjacent to the outer tangential epidermal walls, which enhances the uptake of CO2, and this could be followed by a carbon-concentrating mechanism (CCM) in the cytoplasm and/or chloroplasts. The lines of evidence for such an epidermal CCM are discussed, including evidence for special 'transfer cells' in some but not all seagrass leaves in the tangential inner walls of the epidermal cells. It is concluded that seagrasses have a CCM but that the case for concentration of CO2 at the site of Rubisco carboxylation is not proven.

Transcriptome Analysis of the Response to NaCl in Suaeda maritima Provides an Insight into Salt Tolerance Mechanisms in Halophytes

Although salt tolerance is a feature representative of halophytes, most studies on this topic in plants have been conducted on glycophytes. Transcriptome profiles are also available for only a limited number of halophytes. Hence, the present study was conducted to understand the molecular basis of salt tolerance through the transcriptome profiling of the halophyte Suaeda maritima, which is an emerging plant model for research on salt tolerance. Illumina sequencing revealed 72,588 clustered transcripts, including 27,434 that were annotated using BLASTX. Salt application resulted in the 2-fold or greater upregulation of 647 genes and downregulation of 735 genes. Of these, 391 proteins were homologous to proteins in the COGs (cluster of orthologous groups) database, and the majorities were grouped into the poorly characterized category. Approximately 50% of the genes assigned to MapMan pathways showed homology to S. maritima. The majority of such genes represented transcription factors. Several genes also contributed to cell wall and carbohydrate metabolism, ion relation, redox responses and G protein, phosphoinositide and hormone signaling. Real-time PCR was used to validate the results of the deep sequencing for the most of the genes. This study demonstrates the expression of protein kinase C, the target of diacylglycerol in phosphoinositide signaling, for the first time in plants. This study further reveals that the biochemical and molecular responses occurring at several levels are associated with salt tolerance in S. maritima. At the structural level, adaptations to high salinity levels include the remodeling of cell walls and the modification of membrane lipids. At the cellular level, the accumulation of glycinebetaine and the sequestration and exclusion of Na+ appear to be important. Moreover, this study also shows that the processes related to salt tolerance might be highly complex, as reflected by the salt-induced enhancement of transcription factor expression, including hormone-responsive factors, and that this process might be initially triggered by G protein and phosphoinositide signaling.

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.

De novo transcriptomic analysis of the female and male adults of the blood fluke Schistosoma turkestanicum

Background

Schistosoma turkestanicum is a parasite of considerable veterinary importance as an agent of animal schistosomiasis in many countries, including China. The S. turkestanicum cercariae can also infect humans, causing cercarial dermatitis in many countries and regions of the world. In spite of its significance as a pathogen of animals and humans, there is little transcriptomic and genomic data in the public databases.

Methods

Herein, we performed the transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females of S. turkestanicum and de novo transcriptome assembly.

Results

Approximately 81.1 (female) and 80.5 (male) million high-quality clean reads were obtained and then 29,526 (female) and 41,346 (male) unigenes were assembled. A total of 34,624 unigenes were produced from S. turkestanicum females and males, with an average length of 878 nucleotides (nt) and N50 of 1480 nt. Of these unigenes, 25,158 (72.7 %) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 21,995 (63.5 %), 22,189 (64.1 %) and 13,754 (39.7 %) of the unigenes had significant similarity in the NCBI non-redundant protein (NR), non-redundant nucleotide (NT) and Swiss-Prot databases, respectively. In addition, 3150 unigenes were identified to be expressed specifically in females and 1014 unigenes were identified to be expressed specifically in males. Interestingly, several pathways associated with gonadal development and sex maintenance were found, including the Wnt signaling pathway (103; 2 %) and progesterone-mediated oocyte maturation (77; 1.5 %).

Conclusions

The present study characterized and compared the transcriptomes of adult female and male blood fluke, S. turkestanicum. These results will not only serve as valuable resources for future functional genomics studies to understand the molecular aspects of S. turkestanicum, but also will provide essential information for ongoing whole genome sequencing efforts on this pathogenic blood fluke.

De Novo Sequencing and Analysis of Lemongrass Transcriptome Provide First Insights into the Essential Oil Biosynthesis of Aromatic Grasses

Aromatic grasses of the genus Cymbopogon (Poaceae family) represent unique group of plants that produce diverse composition of monoterpene rich essential oils, which have great value in flavor, fragrance, cosmetic, and aromatherapy industries. Despite the commercial importance of these natural aromatic oils, their biosynthesis at the molecular level remains unexplored. As the first step toward understanding the essential oil biosynthesis, we performed de novo transcriptome assembly and analysis of C. flexuosus (lemongrass) by employing Illumina sequencing. Mining of transcriptome data and subsequent phylogenetic analysis led to identification of terpene synthases, pyrophosphatases, alcohol dehydrogenases, aldo-keto reductases, carotenoid cleavage dioxygenases, alcohol acetyltransferases, and aldehyde dehydrogenases, which are potentially involved in essential oil biosynthesis. Comparative essential oil profiling and mRNA expression analysis in three Cymbopogon species (C. flexuosus, aldehyde type; C. martinii, alcohol type; and C. winterianus, intermediate type) with varying essential oil composition indicated the involvement of identified candidate genes in the formation of alcohols, aldehydes, and acetates. Molecular modeling and docking further supported the role of identified protein sequences in aroma formation in Cymbopogon. Also, simple sequence repeats were found in the transcriptome with many linked to terpene pathway genes including the genes potentially involved in aroma biosynthesis. This work provides the first insights into the essential oil biosynthesis of aromatic grasses, and the identified candidate genes and markers can be a great resource for biotechnological and molecular breeding approaches to modulate the essential oil composition.