Transcriptional homeostasis of a mangrove species, Ceriops tagal, in saline environments, as revealed by microarray analysis

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.

Mangrove Forests: Natural Laboratories for Studying Epigenetic and Climate Changes

In the adaptation market, plants cash in the changes of their DNA (either genetic or epigenetic) to purchase fitness. Fitness is negatively affected by stressful conditions caused by climate change and well-designed studies are required to investigate the fine-tuning cooperation of epigenetic and genetic changes in response to those stresses. Mangrove trees are promising model systems for studying climate change because the effects of environmental changes are already evident in their natural habitats where they are exposed to different salinity levels ranging from saltwater to freshwater. In addition, as mangrove species are known to have very low genetic diversity caused by their stressful living conditions, epigenetic variation is likely to be a vital source for them to respond to environmental changes. This mini review aims to provide an overview of available studies on epigenetic regulation and adaptation of mangroves.

Epigenetic and transcriptional responses underlying mangrove adaptation to UV-B

Tropical plants have adapted to strong solar ultraviolet (UV) radiation. Here we compare molecular responses of two tropical mangroves Avecennia marina and Rhizophora apiculata to high-dose UV-B. Whole-genome bisulfate sequencing indicates that high UV-B induced comparable hyper- or hypo-methylation in three sequence contexts (CG, CHG, and CHH, where H refers to A, T, or C) in A. marina but mainly CHG hypomethylation in R. apiculata. RNA and small RNA sequencing reveals UV-B induced relaxation of transposable element (TE) silencing together with up-regulation of TE-adjacent genes in R. apiculata but not in A. marina. Despite conserved upregulation of flavonoid biosynthesis and downregulation of photosynthesis genes caused by high UV-B, A. marina specifically upregulated ABC transporter and ubiquinone biosynthesis genes that are known to be protective against UV-B-induced damage. Our results point to divergent responses underlying plant UV-B adaptation at both the epigenetic and transcriptional level.

Gene Body Methylation Confers Transcription Robustness in Mangroves During Long-Term Stress Adaptation

Whether induced epigenetic changes contribute to long-term adaptation remains controversial. Recent studies indicate that environmentally cued changes in gene body methylation (gbM) can facilitate acclimatization. However, such changes are often associated with genetic variation and their contribution to long-term stress adaptation remains unclear. Using whole-genome bisulfite sequencing, we examined evolutionary gains and losses of gbM in mangroves that adapted to extreme intertidal environments. We treated mangrove seedlings with salt stress, and investigated expression changes in relation with stress-induced or evolutionarily-acquired gbM changes. Evolution and function of gbM was compared with that of genetic variation. Mangroves gained much more gbM than their terrestrial relatives, mainly through convergent evolution. Genes that convergently gained gbM during evolution are more likely to become methylated in response to salt stress in species where they are normally not marked. Stress-induced and evolutionarily convergent gains of gbM both correlate with reduction in expression variation, conferring genome-wide expression robustness under salt stress. Moreover, convergent gbM evolution is uncoupled with convergent sequence evolution. Our findings suggest that transgenerational inheritance of acquired gbM helps environmental canalization of gene expression, facilitating long-term stress adaptation of mangroves in the face of a severe reduction in genetic diversity.

Molecular adaptation to salinity fluctuation in tropical intertidal environments of a mangrove tree Sonneratia alba

BACKGROUND

Mangroves have adapted to intertidal zones - the interface between terrestrial and marine ecosystems. Various studies have shown adaptive evolution in mangroves at physiological, ecological, and genomic levels. However, these studies paid little attention to gene regulation of salt adaptation by transcriptome profiles.

RESULTS

We sequenced the transcriptomes of Sonneratia alba under low (fresh water), medium (half the seawater salinity), and high salt (seawater salinity) conditions and investigated the underlying transcriptional regulation of salt adaptation. In leaf tissue, 64% potential salinity-related genes were not differentially expressed when salinity increased from freshwater to medium levels, but became up- or down-regulated when salt concentrations further increased to levels found in sea water, indicating that these genes are well adapted to the medium saline condition. We inferred that both maintenance and regulation of cellular environmental homeostasis are important adaptive processes in S. alba. i) The sulfur metabolism as well as flavone and flavonol biosynthesis KEGG pathways were significantly enriched among up-regulated genes in leaves. They are both involved in scavenging ROS or synthesis and accumulation of osmosis-related metabolites in plants. ii) There was a significantly increased percentage of transcription factor-encoding genes among up-regulated transcripts. High expressions of salt tolerance-related TF families were found under high salt conditions. iii) Some genes up-regulated in response to salt treatment showed signs of adaptive evolution at the amino acid level and might contribute to adaptation to fluctuating intertidal environments.

CONCLUSIONS

This study first elucidates the mechanism of high-salt adaptation in mangroves at the whole-transcriptome level by salt gradient experimental treatments. It reveals that several candidate genes (including salt-related genes, TF-encoding genes, and PSGs) and major pathways are involved in adaptation to high-salt environments. Our study also provides a valuable resource for future investigation of adaptive evolution in extreme environments.

Identification and characterization of evolutionarily conserved alternative splicing events in a mangrove genus Sonneratia

Alternative splicing (AS), which produces multiple mRNA transcripts from a single gene, plays crucial roles in plant growth, development and environmental stress responses. Functional significances of conserved AS events among congeneric species have not been well characterized. In this study, we performed transcriptome sequencing to characterize AS events in four common species of Sonneratia, a mangrove genus excellently adaptive to intertidal zones. 7,248 to 12,623 AS events were identified in approximately 25% to 35% expressed genes in the roots of the four species. The frequency of AS events in Sonneratia was associated with genomic features, including gene expression level and intron/exon number and length. Among the four species, 1,355 evolutionarily conserved AS (ECAS) events were identified from 1,170 genes. Compared with non-ECAS events, ECAS events are of shorter length and less possibility to introduce premature stop codons (PTCs) and frameshifts. Functional annotations of the genes containing ECAS events showed that four of the 26 enriched Gene Ontology (GO) terms are involved in proton transport, signal transduction and carbon metabolism, and 60 genes from another three GO terms are implicated in responses to osmotic, oxidative and heat stresses, which may contribute to the adaptation of Sonneratia species to harsh intertidal environments.

Convergent adaptive evolution in marginal environments: unloading transposable elements as a common strategy among mangrove genomes

Several clades of mangrove trees independently invade the interface between land and sea at the margin of woody plant distribution. As phenotypic convergence among mangroves is common, the possibility of convergent adaptation in their genomes is quite intriguing. To study this molecular convergence, we sequenced multiple mangrove genomes. In this study, we focused on the evolution of transposable elements (TEs) in relation to the genome size evolution. TEs, generally considered genomic parasites, are the most common components of woody plant genomes. Analyzing the long terminal repeat-retrotransposon (LTR-RT) type of TE, we estimated their death rates by counting solo-LTRs and truncated elements. We found that all lineages of mangroves massively and convergently reduce TE loads in comparison to their nonmangrove relatives; as a consequence, genome size reduction happens independently in all six mangrove lineages; TE load reduction in mangroves can be attributed to the paucity of young elements; the rarity of young LTR-RTs is a consequence of fewer births rather than access death. In conclusion, mangrove genomes employ a convergent strategy of TE load reduction by suppressing element origination in their independent adaptation to a new environment.

Transcriptomics analysis of salt stress tolerance in the roots of the mangrove Avicennia officinalis

Salinity affects growth and development of plants, but mangroves exhibit exceptional salt tolerance. With direct exposure to salinity, mangrove roots possess specific adaptations to tolerate salt stress. Therefore, studying the early effects of salt on mangrove roots can help us better understand the tolerance mechanisms. Using two-month-old greenhouse-grown seedlings of the mangrove tree Avicennia officinalis subjected to NaCl treatment, we profiled gene expression changes in the roots by RNA-sequencing. Of the 6547 genes that were differentially regulated in response to salt treatment, 1404 and 5213 genes were significantly up- and down-regulated, respectively. By comparative genomics, 93 key salt tolerance-related genes were identified of which 47 were up-regulated. Upon placing all the differentially expressed genes (DEG) in known signaling pathways, it was evident that most of the DEGs involved in ethylene and auxin signaling were up-regulated while those involved in ABA signaling were down-regulated. These results imply that ABA-independent signaling pathways also play a major role in salt tolerance of A. officinalis. Further, ethylene response factors (ERFs) were abundantly expressed upon salt treatment and the Arabidopsis mutant aterf115, a homolog of AoERF114 is characterized. Overall, our results would help in understanding the possible molecular mechanism underlying salt tolerance in plants.

The origin, diversification and adaptation of a major mangrove clade (Rhizophoreae) revealed by whole-genome sequencing

Mangroves invade some very marginal habitats for woody plants-at the interface between land and sea. Since mangroves anchor tropical coastal communities globally, their origin, diversification and adaptation are of scientific significance, particularly at a time of global climate change. In this study, a combination of single-molecule long reads and the more conventional short reads are generated from Rhizophora apiculata for the de novo assembly of its genome to a near chromosome level. The longest scaffold, N50 and N90 for the R. apiculata genome, are 13.3 Mb, 5.4 Mb and 1.0 Mb, respectively. Short reads for the genomes and transcriptomes of eight related species are also generated. We find that the ancestor of Rhizophoreae experienced a whole-genome duplication ~70 Myrs ago, which is followed rather quickly by colonization and species diversification. Mangroves exhibit pan-exome modifications of amino acid (AA) usage as well as unusual AA substitutions among closely related species. The usage and substitution of AAs, unique among plants surveyed, is correlated with the rapid evolution of proteins in mangroves. A small subset of these substitutions is associated with mangroves' highly specialized traits (vivipary and red bark) thought to be adaptive in the intertidal habitats. Despite the many adaptive features, mangroves are among the least genetically diverse plants, likely the result of continual habitat turnovers caused by repeated rises and falls of sea level in the geologically recent past. Mangrove genomes thus inform about their past evolutionary success as well as portend a possibly difficult future.

A multilevel investigation to discover why Kandelia candel thrives in high salinity

The halophilic mangrove species Kandelia candel is an excellent model for understanding why halophytes thrive in high salinity. Preliminary transcriptomic analyses revealed that genes involved in diverse functions, such as in phenylpropanoid and amino acid metabolisms, and those in DNA replication and damage repair were highly responsive to salt stress. Proteomic analyses revealed that the proteins involved in light reaction of photosynthesis, amino acid and carbohydrate metabolisms, secondary metabolite biosynthesis and posttranslational modification showed increased levels in response to salt stress. The metabolisms of phenylpropanoids and amino acids under salt stress were systematically examined based on the preliminary omics analyses. The activities of phenylpropanoid biosynthetic enzymes and the contents of phenols, flavonoids, anthocyanins and lignins were significantly increased under salt stress. In the free amino acid pool, glutamate was the most abundant. Together with γ-aminobutyric acid, glutamate levels further increased, while proline levels remained unchanged in response to salt stress. These findings point to the potential importance of phenylpropanoids and free amino acids in salt tolerance of K. candel that have been observed, but not systemically investigated at the levels of gene expression, enzyme activity and metabolite accumulation in glycophytes and non-tree halophytes.

Small RNA transcriptomes of mangroves evolve adaptively in extreme environments

MicroRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs) are key players in plant stress responses. Here, we present the sRNA transcriptomes of mangroves Bruguiera gymnorrhiza and Kandelia candel. Comparative computational analyses and target predictions revealed that mangroves exhibit distinct sRNA regulatory networks that differ from those of glycophytes. A total of 32 known and three novel miRNA families were identified. Conserved and mangrove-specific miRNA targets were predicted; the latter were widely involved in stress responses. The known miRNAs showed differential expression between the mangroves and glycophytes, reminiscent of the adaptive stress-responsive changes in Arabidopsis. B. gymnorrhiza possessed highly abundant but less conserved TAS3 trans-acting siRNAs (tasiRNAs) in addition to tasiR-ARFs, with expanded potential targets. Our results indicate that the evolutionary alteration of sRNA expression levels and the rewiring of sRNA-regulatory networks are important mechanisms underlying stress adaptation. We also identified sRNAs that are involved in salt and/or drought tolerance and nutrient homeostasis as possible contributors to mangrove success in stressful environments.

Transcriptome analysis of the Holly mangrove Acanthus ilicifolius and its terrestrial relative, Acanthus leucostachyus, provides insights into adaptation to intertidal zones

BACKGROUND

Acanthus is a unique genus consisting of both true mangrove and terrestrial species; thus, it represents an ideal system for studying the origin and adaptive evolution of mangrove plants to intertidal environments. However, little is known regarding the two respects of mangrove species in Acanthus. In this study, we sequenced the transcriptomes of the pooled roots and leaves tissues for a mangrove species, Acanthus ilicifolius, and its terrestrial congener, A. leucostachyus, to illustrate the origin of the mangrove species in this genus and their adaptive evolution to harsh habitats.

RESULTS

We obtained 73,039 and 69,580 contigs with N50 values of 741 and 1557 bp for A. ilicifolius and A. leucostachyus, respectively. Phylogenetic analyses based on four nuclear segments and three chloroplast fragments revealed that mangroves and terrestrial species in Acanthus fell into different clades, indicating a single origin of the mangrove species in Acanthus. Based on 6634 orthologs, A. ilicifolius and A. leucostachyus were found to be highly divergent, with a peak of synonymous substitution rate (Ks) distribution of 0.145 and an estimated divergence time of approximately 16.8 million years ago (MYA). The transgression in the Early to Middle Miocene may be the major reason for the entry of the mangrove lineage of Acanthus into intertidal environments. Gene ontology (GO) classifications of the full transcriptomes did not show any apparent differences between A. ilicifolius and A. leucostachyus, suggesting the absence of gene components specific to the mangrove transcriptomes. A total of 99 genes in A. ilicifolius were identified with signals of positive selection. Twenty-three of the 99 positively selected genes (PSGs) were found to be involved in salt, heat and ultraviolet stress tolerance, seed germination and embryo development under periodic inundation. These stress-tolerance related PSGs may be crucial for the adaptation of the mangrove species in this genus to stressful marine environments and may contribute to speciation in Acanthus.

CONCLUSIONS

We characterized the transcriptomes of one mangrove species of Acanthus, A. ilicifolius, and its terrestrial relative, A. leucostachyus, and provided insights into the origin of the mangrove Acanthus species and their adaptive evolution to abiotic stresses in intertidal environments.

Transcriptome assembly, profiling and differential gene expression analysis of the halophyte Suaeda fruticosa provides insights into salt tolerance

BACKGROUND

Improvement of crop production is needed to feed the growing world population as the amount and quality of agricultural land decreases and soil salinity increases. This has stimulated research on salt tolerance in plants. Most crops tolerate a limited amount of salt to survive and produce biomass, while halophytes (salt-tolerant plants) have the ability to grow with saline water utilizing specific biochemical mechanisms. However, little is known about the genes involved in salt tolerance. We have characterized the transcriptome of Suaeda fruticosa, a halophyte that has the ability to sequester salts in its leaves. Suaeda fruticosa is an annual shrub in the family Chenopodiaceae found in coastal and inland regions of Pakistan and Mediterranean shores. This plant is an obligate halophyte that grows optimally from 200-400 mM NaCl and can grow at up to 1000 mM NaCl. High throughput sequencing technology was performed to provide understanding of genes involved in the salt tolerance mechanism. De novo assembly of the transcriptome and analysis has allowed identification of differentially expressed and unique genes present in this non-conventional crop.

RESULTS

Twelve sequencing libraries prepared from control (0 mM NaCl treated) and optimum (300 mM NaCl treated) plants were sequenced using Illumina Hiseq 2000 to investigate differential gene expression between shoots and roots of Suaeda fruticosa. The transcriptome was assembled de novo using Velvet and Oases k-45 and clustered using CDHIT-EST. There are 54,526 unigenes; among these 475 genes are downregulated and 44 are upregulated when samples from plants grown under optimal salt are compared with those grown without salt. BLAST analysis identified the differentially expressed genes, which were categorized in gene ontology terms and their pathways.

CONCLUSIONS

This work has identified potential genes involved in salt tolerance in Suaeda fruticosa, and has provided an outline of tools to use for de novo transcriptome analysis. The assemblies that were used provide coverage of a considerable proportion of the transcriptome, which allows analysis of differential gene expression and identification of genes that may be involved in salt tolerance. The transcriptome may serve as a reference sequence for study of other succulent halophytes.

Proteomic analysis of salt-responsive proteins in the leaves of mangrove Kandelia candel during short-term stress

Salt stress is a major abiotic stress that limits crop productivity in many regions of the world. A comparative proteomic approach to identify salt stress-responsive proteins and to understand the molecular mechanisms was carried out in the woody halophyte Kandelia candel. Four-leaf-old K. candel seedlings were exposed to 150 (control), 300, 450, and 600 mM NaCl for 3 days. Proteins extracted from the leaves of K. candel seedlings were separated by two-dimensional gel electrophoresis (2-DE). More than 900 protein spots were detected on each gel, and 53 differentially expressed protein spots were located with at least two-fold differences in abundance on 2-DE maps, of which 48 were identified by matrix-assisted laser desorption ionization time-of-flight/time-of-flight mass spectrometry (MALDI-TOF-TOF/MS). The results showed that K. candel could withstand up to 450 mM NaCl stress by up-regulating proteins that are mainly involved in photosynthesis, respiration and energy metabolism, Na(+) compartmentalization, protein folding and assembly, and signal transduction. Physiological data, including superoxide dismutase (SOD) and dehydroascorbate reductase (DHAR) activities, hydrogen peroxide (H2O2) and superoxide anion radicals (O2(-)) contents, as well as Na(+) content and K(+)/Na(+) ratios all correlated well with our proteomic results. This study provides new global insights into woody halophyte salt stress responses. Identification of differentially expressed proteins promotes better understanding of the molecular basis for salt stress reduction in K. candel.