Comparative transcriptomics for mangrove species: an expanding resource

Positive selection and heat-response transcriptomes reveal adaptive features of the Brassicaceae desert model, Anastatica hierochuntica

Plant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts. We generated an A. hierochuntica reference transcriptome and identified extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat, and detecting positively selected genes in A. hierochuntica. The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat 'stress-ready' state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype. We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.

De novo genome assembly and in natura epigenomics reveal salinity-induced DNA methylation in the mangrove tree Bruguiera gymnorhiza

Mangroves are adapted to harsh environments, such as high ultraviolet (UV) light, low nutrition, and fluctuating salinity in coastal zones. However, little is known about the transcriptomic and epigenomic basis of the resilience of mangroves due to limited available genome resources. We performed a de novo genome assembly and in natura epigenome analyses of the mangrove Bruguiera gymnorhiza, one of the dominant mangrove species. We also performed the first genome-guided transcriptome assembly for mangrove species. The 309 Mb of the genome is predicted to encode 34 403 genes and has a repeat content of 48%. Depending on its growing environment, the natural B. gymnorhiza population showed drastic morphological changes associated with expression changes in thousands of genes. Moreover, high-salinity environments induced genome-wide DNA hypermethylation of transposable elements (TEs) in the B. gymnorhiza. DNA hypermethylation was concurrent with the transcriptional regulation of chromatin modifier genes, suggesting robust epigenome regulation of TEs in the B. gymnorhiza genome under high-salinity environments. The genome and epigenome data in this study provide novel insights into the epigenome regulation of mangroves and a better understanding of the adaptation of plants to fluctuating, harsh natural environments.

Characterization of Laguncularia racemosa transcriptome and molecular response to oil pollution

Mangroves are ecosystems of economic and ecological importance. Laguncularia racemosa (Combretaceae), popularly known as white mangrove, is a species that greatly contributes to the community structure of neotropical and West African mangrove forests. Despite the significance of these ecosystems, they have been destroyed by oil spills that can cause yellowing of leaves, increased sensitivity to other stresses and death of trees. However, the molecular response of plants to oil stress is poorly known. In this work, Illumina reads were de novo assembled into 46,944 transcripts of L. racemosa roots and leaves, including putative isoform variants. In addition to improving the genomic information available for mangroves, the L. racemosa assembled transcriptome allowed us to identify reference genes to normalize quantitative real-time PCR (qPCR) expression data from oil-stressed mangrove plants, which were used in RNASeq validation. The analysis of expression changes induced by the oil exposure revealed 310 and 286 responsive transcripts of leaves and roots, respectively, mainly up-regulated. Enriched GO categories related to chloroplasts and photosynthesis were found among both leaf and root oil-responsive transcripts, while "response to heat" and "response to hypoxia" were exclusively enriched in leaves and roots, respectively. The comparison of L. racemosa 12-h-oil-stressed leaf expression profile to previous Arabidopsis heat-stress studies and co-expression evidence also pointed to similarities between the heat and oil responses, in which the HSP-coding genes seem to play a key role. A subset of the L. racemosa oil-responsive root genes exhibited similar up-regulation profiles to their Arabidopsis homologs involved in hypoxia responses, including the HRA1 and LBD41 TF-coding genes. Genes linked to the ethylene pathway such as those coding for ERF TFs were also modulated during the L. racemosa root response to oil stress. Taken together, these results show that oil contamination affects photosynthesis, protein metabolism, hypoxia response and the ethylene pathway in L. racemosa 12-h-oil-exposed leaves and roots.

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.

Comparative Analysis of Transcriptomes in Rhizophoraceae Provides Insights into the Origin and Adaptive Evolution of Mangrove Plants in Intertidal Environments

Mangroves are woody plants that grow at the interface between land and sea in tropical and subtropical latitudes, where they exist in conditions of high salinity, extreme tides, strong winds, high temperatures, and muddy, anaerobic soils. Rhizophoraceae is a key mangrove family, with highly developed morphological and physiological adaptations to extreme conditions. It is an ideal system for the study of the origin and adaptive evolution of mangrove plants. In this study, we characterized and comprehensively compared the transcriptomes of four mangrove species, from all four mangrove genera, as well as their closest terrestrial relative in Rhizophoraceae, using RNA-Seq. We obtained 41,936-48,845 unigenes with N50 values of 982-1,185 bp and 61.42-69.48% annotated for the five species in Rhizophoraceae. Orthology annotations of Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Clusters of Orthologous Groups revealed overall similarities in the transcriptome profiles among the five species, whereas enrichment analysis identified remarkable genomic characteristics that are conserved across the four mangrove species but differ from their terrestrial relative. Based on 1,816 identified orthologs, phylogeny analysis and divergence time estimation revealed a single origin for mangrove species in Rhizophoraceae, which diverged from the terrestrial lineage ~56.4 million years ago (Mya), suggesting that the transgression during the Paleocene-Eocene Thermal Maximum may have been responsible for the entry of the mangrove lineage of Rhizophoraceae into intertidal environments. Evidence showed that the ancestor of Rhizophoraceae may have experienced a whole genome duplication event ~74.6 Mya, which may have increased the adaptability and survival chances of Rhizophoraceae during and following the Cretaceous-Tertiary extinction. The analysis of positive selection identified 10 positively selected genes from the ancestor branch of Rhizophoraceae mangroves, which were mainly associated with stress response, embryo development, and regulation of gene expression. Positive selection of these genes may be crucial for increasing the capability of stress tolerance (i.e., defense against salt and oxidative stress) and development of adaptive traits (i.e., vivipary) of Rhizophoraceae mangroves, and thus plays an important role in their adaptation to the stressful intertidal environments.

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.

Identification of salt gland-associated genes and characterization of a dehydrin from the salt secretor mangrove Avicennia officinalis

BACKGROUND

Salt stress is a major challenge for growth and development of plants. The mangrove tree Avicennia officinalis has evolved salt tolerance mechanisms such as salt secretion through specialized glands on its leaves. Although a number of structural studies on salt glands have been done, the molecular mechanism of salt secretion is not clearly understood. Also, studies to identify salt gland-specific genes in mangroves have been scarce.

RESULTS

By subtractive hybridization (SH) of cDNA from salt gland-rich cell layers (tester) with mesophyll tissues as the driver, several Expressed Sequence Tags (ESTs) were identified. The major classes of ESTs identified include those known to be involved in regulating metabolic processes (37%), stress response (17%), transcription (17%), signal transduction (17%) and transport functions (12%). A visual interactive map generated based on predicted functional gene interactions of the identified ESTs suggested altered activities of hydrolase, transmembrane transport and kinases. Quantitative Real-Time PCR (qRT-PCR) was carried out to validate the expression specificity of the ESTs identified by SH. A Dehydrin gene was chosen for further experimental analysis, because it is significantly highly expressed in salt gland cells, and dehydrins are known to be involved in stress remediation in other plants. Full-length Avicennia officinalis Dehydrin1 (AoDHN1) cDNA was obtained by Rapid Amplification of cDNA Ends. Phylogenetic analysis and further characterization of this gene suggested that AoDHN1 belongs to group II Late Embryogenesis Abundant proteins. qRT-PCR analysis of Avicennia showed up-regulation of AoDHN1 in response to salt and drought treatments. Furthermore, some functional insights were obtained by growing E. coli cells expressing AoDHN1. Growth of E. coli cells expressing AoDHN1 was significantly higher than that of the control cells without AoDHN1 under salinity and drought stresses, suggesting that the mangrove dehydrin protein helps to mitigate the abiotic stresses.

CONCLUSIONS

Thirty-four ESTs were identified to be enriched in salt gland-rich tissues of A. officinalis leaves. qRT-PCR analysis showed that 10 of these were specifically enriched in the salt gland-rich tissues. Our data suggest that one of the selected genes, namely, AoDHN1 plays an important role to mitigate salt and drought stress responses.

Identification and analysis of red sea mangrove (Avicennia marina) microRNAs by high-throughput sequencing and their association with stress responses

Although RNA silencing has been studied primarily in model plants, advances in high-throughput sequencing technologies have enabled profiling of the small RNA components of many more plant species, providing insights into the ubiquity and conservatism of some miRNA-based regulatory mechanisms. Small RNAs of 20 to 24 nucleotides (nt) are important regulators of gene transcript levels by either transcriptional or by posttranscriptional gene silencing, contributing to genome maintenance and controlling a variety of developmental and physiological processes. Here, we used deep sequencing and molecular methods to create an inventory of the small RNAs in the mangrove species, Avicennia marina. We identified 26 novel mangrove miRNAs and 193 conserved miRNAs belonging to 36 families. We determined that 2 of the novel miRNAs were produced from known miRNA precursors and 4 were likely to be species-specific by the criterion that we found no homologs in other plant species. We used qRT-PCR to analyze the expression of miRNAs and their target genes in different tissue sets and some demonstrated tissue-specific expression. Furthermore, we predicted potential targets of these putative miRNAs based on a sequence homology and experimentally validated through endonucleolytic cleavage assays. Our results suggested that expression profiles of miRNAs and their predicted targets could be useful in exploring the significance of the conservation patterns of plants, particularly in response to abiotic stress. Because of their well-developed abilities in this regard, mangroves and other extremophiles are excellent models for such exploration.

Life at the extreme: lessons from the genome

Extremophile plants thrive in places where most plant species cannot survive. Recent developments in high-throughput technologies and comparative genomics are shedding light on the evolutionary mechanisms leading to their adaptation.

Transcriptome characterization and sequencing-based identification of salt-responsive genes in Millettia pinnata, a semi-mangrove plant

Semi-mangroves form a group of transitional species between glycophytes and halophytes, and hold unique potential for learning molecular mechanisms underlying plant salt tolerance. Millettia pinnata is a semi-mangrove plant that can survive a wide range of saline conditions in the absence of specialized morphological and physiological traits. By employing the Illumina sequencing platform, we generated ~192 million short reads from four cDNA libraries of M. pinnata and processed them into 108,598 unisequences with a high depth of coverage. The mean length and total length of these unisequences were 606 bp and 65.8 Mb, respectively. A total of 54,596 (50.3%) unisequences were assigned Nr annotations. Functional classification revealed the involvement of unisequences in various biological processes related to metabolism and environmental adaptation. We identified 23,815 candidate salt-responsive genes with significantly differential expression under seawater and freshwater treatments. Based on the reverse transcription-polymerase chain reaction (RT-PCR) and real-time PCR analyses, we verified the changes in expression levels for a number of candidate genes. The functional enrichment analyses for the candidate genes showed tissue-specific patterns of transcriptome remodelling upon salt stress in the roots and the leaves. The transcriptome of M. pinnata will provide valuable gene resources for future application in crop improvement. In addition, this study sets a good example for large-scale identification of salt-responsive genes in non-model organisms using the sequencing-based approach.