NADPH oxidase-dependent H2O2 production mediates salicylic acid-induced salt tolerance in mangrove plant Kandelia obovata by regulating Na+/K+ and redox homeostasis.

Salicylic acid (SA) is known to enhance salt tolerance in plants. However, the mechanism of SA-mediated response to high salinity in halophyte remains unclear. Using electrophysiological and molecular biological methods, we investigated the role of SA in response to high salinity in mangrove species, Kandelia obovata, a typical halophyte. Exposure of K. obovata roots to high salinity resulted in a rapid increase in endogenous SA produced by phenylalanine ammonia lyase pathway. The application of exogenous SA improved the salt tolerance of K. obovata, which depended on the NADPH oxidase-mediated H2O2. Exogenous SA and H2O2 increased Na+ efflux and reduced K+ loss by regulating the transcription levels of Na+ and K+ transport-related genes, thus reducing the Na+/K+ ratio in the salt-treated K. obovata roots. In addition, exogenous SA-enhanced antioxidant enzyme activity and its transcripts, and the expressions of four genes related to AsA-GSH cycle as well, then alleviated oxidative damages in the salt-treated K. obovata roots. However, the above effects of SA could be reversed by diphenyleneiodonium chloride (the NADPH oxidase inhibitor) and paclobutrazol (a SA biosynthesis inhibitor). Collectively, our results demonstrated that SA-induced salt tolerance of K. obovata depends on NADPH oxidase-generated H2O2 that affects Na+/K+ and redox homeostasis in response to high salinity.

Comparative transcriptomics of the chilling stress response in two Asian mangrove species, Bruguiera gymnorhiza and Rhizophora apiculata.

Low temperatures largely determine the geographic limits of plant species by reducing survival and growth. Inter-specific differences in the geographic distribution of mangrove species have been associated with cold tolerance, with exclusively tropical species being highly cold-sensitive and subtropical species being relatively cold-tolerant. To identify species-specific adaptations to low temperatures, we compared the chilling stress response of two widespread Indo-West Pacific mangrove species from Rhizophoraceae with differing latitudinal range limits-Bruguiera gymnorhiza (L.) Lam. ex Savigny (subtropical range limit) and Rhizophora apiculata Blume (tropical range limit). For both species, we measured the maximum photochemical efficiency of photosystem II (Fv/Fm) as a proxy for the physiological condition of the plants and examined gene expression profiles during chilling at 15 and 5 °C. At 15 °C, B. gymnorhiza maintained a significantly higher Fv/Fm than R. apiculata. However, at 5 °C, both species displayed equivalent Fv/Fm values. Thus, species-specific differences in chilling tolerance were only found at 15 °C, and both species were sensitive to chilling at 5 °C. At 15 °C, B. gymnorhiza downregulated genes related to the light reactions of photosynthesis and upregulated a gene involved in cyclic electron flow regulation, whereas R. apiculata downregulated more RuBisCo-related genes. At 5 °C, both species repressed genes related to CO2 assimilation. The downregulation of genes related to light absorption and upregulation of genes related to cyclic electron flow regulation are photoprotective mechanisms that likely contributed to the greater photosystem II photochemical efficiency of B. gymnorhiza at 15 °C. The results of this study provide evidence that the distributional range limits and potentially the expansion rates of plant species are associated with differences in the regulation of photosynthesis and photoprotective mechanisms under low temperatures.

Brassinosteroid enhances salt tolerance via S-nitrosoglutathione reductase and nitric oxide signaling pathway in mangrove Kandelia obovata.

Brassinosteroid (BR) has been shown to modulate plant tolerance to various stresses. S-nitrosoglutathione reductase (GSNOR) is involved in the plant response to environment stress by fine-turning the level of nitric oxide (NO). However, whether GSNOR is involved in BR-regulated Na+ /K+ homeostasis to improve the salt tolerance in halophyte is unknown. Here, we firstly reported that high salinity increases the expression of BR-biosynthesis genes and the endogenous levels of BR in mangrove Kandelia obovata. Then, salt-induced BR triggers the activities and gene expressions of GSNOR and antioxidant enzymes, thereafter decrease the levels of malondialdehyde, hydrogen peroxide. Subsequently, BR-mediated GSNOR negatively regulates NO contributions to the reduction of reactive oxygen species generation and induction of the gene expression related to Na+ and K+ transport, leading to the decrease of Na+ /K+ ratio in the roots of K. obovata. Finally, the applications of exogenous BR, NO scavenger, BR biosynthetic inhibitor and GSNOR inhibitor further confirm the function of BR. Taken together, our result provides insight into the mechanism of BR in the response of mangrove K. obovata to high salinity via GSNOR and NO signaling pathway by reducing oxidative damage and modulating Na+ /K+ homeostasis.

Comparative Transcriptomics and Metabolomics Analyses of Avicennia marina and Kandelia obovata under Chilling Stress during Seedling Stage.

One of the most productive ecosystems in the world, mangroves are susceptible to cold stress. However, there is currently insufficient knowledge of the adaptation mechanisms of mangrove plants in response to chilling stress. This study conducted a comparative analysis of transcriptomics and metabolomics to investigate the adaptive responses of Kandelia obovata (chilling-tolerant) and Avicennia marina (chilling-sensitive) to 5 °C. The transcriptomics results revealed that differentially expressed genes (DEGs) were mostly enriched in signal transduction, photosynthesis-related pathways, and phenylpropanoid biosynthesis. The expression pattern of genes involved in photosynthesis-related pathways in A. marina presented a downregulation of most DEGs, which correlated with the decrease in total chlorophyll content. In the susceptible A. marina, all DEGs encoding mitogen-activated protein kinase were upregulated. Phenylpropanoid-related genes were observed to be highly induced in K. obovata. Additionally, several metabolites, such as 4-aminobutyric acid, exhibited higher levels in K. obovata than in A. marina, suggesting that chilling-tolerant varieties regulated more metabolites in response to chilling. The investigation defined the inherent distinctions between K. obovata and A. marina in terms of signal transduction gene expression, as well as phenylpropanoid and flavonoid biosynthesis, during exposure to low temperatures.

Genome-Wide Identification, Characterization, and Expression Analysis of the Copper-Containing Amine Oxidase Gene Family in Mangrove Kandelia obovata.

Copper-containing amine oxidases (CuAOs) are known to have significant involvement in the process of polyamine catabolism, as well as serving crucial functions in plant development and response to abiotic stress. A genome-wide investigation of the CuAO protein family was previously carried out in sweet orange (Citrus sinensis) and sweet cherry (Prunus avium L.). Six CuAO (KoCuAO1-KoCuAO6) genes were discovered for the first time in the Kandelia obovata (Ko) genome through a genome-wide analysis conducted to better understand the key roles of the CuAO gene family in Kandelia obovata. This study encompassed an investigation into various aspects of gene analysis, including gene characterization and identification, subcellular localization, chromosomal distributions, phylogenetic tree analysis, gene structure analysis, motif analysis, duplication analysis, cis-regulatory element identification, domain and 3D structural variation analysis, as well as expression profiling in leaves under five different treatments of copper (CuCl2). Phylogenetic analysis suggests that these KoCuAOs, like sweet cherry, may be subdivided into three subgroups. Examining the chromosomal location revealed an unequal distribution of the KoCuAO genes across four out of the 18 chromosomes in Kandelia obovata. Six KoCuAO genes have coding regions with 106 and 159 amino acids and exons with 4 and 12 amino acids. Additionally, we discovered that the 2.5 kb upstream promoter region of the KoCuAOs predicted many cis elements linked to phytohormones and stress responses. According to the expression investigations, CuCl2 treatments caused up- and downregulation of all six genes. In conclusion, our work provides a comprehensive overview of the expression pattern and functional variety of the Kandelia obovata CuAO gene family, which will facilitate future functional characterization of each KoCuAO gene.

Genome-Wide Analysis of bHLH Family Genes and Identification of Members Associated with Cold/Drought-Induced Photoinhibition in Kandelia obovata.

Plant basic helix-loop-helix (bHLH) transcription factors play pivotal roles in responding to stress, including cold and drought. However, it remains unclear how bHLH family genes respond to these stresses in Kandelia obovata. In this study, we identified 75 bHLH members in K. obovata, classified into 11 subfamilies and unevenly distributed across its 18 chromosomes. Collineation analysis revealed that segmental duplication primarily drove the expansion of KobHLH genes. The KobHLH promoters were enriched with elements associated with light response. Through RNA-seq, we identified several cold/drought-associated KobHLH genes. This correlated with decreased net photosynthetic rates (Pn) in the leaves of cold/drought-treated plants. Weighted gene co-expression network analysis (WGCNA) confirmed that 11 KobHLH genes were closely linked to photoinhibition in photosystem II (PS II). Among them, four Phytochrome Interacting Factors (PIFs) involved in chlorophyll metabolism were significantly down-regulated. Subcellular localization showed that KobHLH52 and KobHLH30 were located in the nucleus. Overall, we have comprehensively analyzed the KobHLH family and identified several members associated with photoinhibition under cold or drought stress, which may be helpfulfor further cold/drought-tolerance enhancement and molecular breeding through genetic engineering in K. obovata.

Comparative chloroplast genome analyses provide insights into evolutionary history of Rhizophoraceae mangroves.

Background: The Rhizophoraceae family comprises crucial mangrove plants that inhabit intertidal environments. In China, eight Rhizophoraceae mangrove species exist. Although complete chloroplast (Cp) genomes of four Rhizophoraceae mangrove plants have been reported, the Cp genomes of the remaining four species remain unclear, impeding a comprehensive understanding of the evolutionary history of this family. Methods: Illumina high-throughput sequencing was employed to obtain the DNA sequences of Rhizophoraceae species. Cp genomes were assembled by NOVOPlasty and annotated using CpGAVAS software. Phylogenetic and divergence time analyses were conducted using MEGA and BEAST 2 software. Results: Four novel Cp genomes of Rhizophoraceae mangrove species (Bruguiera sexangula, Bruguiera gymnorrhiza, Bruguiera × rhynchopetala and Rhizophora apiculata) were successfully assembled. The four Cp genomes ranged in length from 163,310 to 164,560 bp, with gene numbers varying from 124 to 128. The average nucleotide diversity (Pi) value of the eight Rhizophoraceae Cp genomes was 0.00596. Phylogenetic trees constructed based on the complete Cp genomes supported the monophyletic origin of Rhizophoraceae. Divergence time estimation based on the Cp genomes of representative species from Malpighiales showed that the origin of Rhizophoraceae occurred at approximately 58.54-50.02 million years ago (Mya). The divergence time within the genus Rhizophora (∼4.51 Mya) was much earlier than the divergence time within the genus Bruguiera (∼1.41 Mya), suggesting recent speciation processes in these genera. Our data provides new insights into phylogenetic relationship and evolutionary history of Rhizophoraceae mangrove plants.

Genome-Wide Identification and Expression Analysis of the Copper Transporter (COPT/Ctr) Gene Family in Kandelia obovata, a Typical Mangrove Plant.

The copper transporter (COPT/Ctr) gene family plays a critical part in maintaining the balance of the metal, and many diverse species depend on COPT to move copper (Cu) across the cell membrane. In Arabidopsis thaliana, Oryza sativa, Medicago sativa, Zea mays, Populus trichocarpa, Vitis vinifera, and Solanum lycopersicum, a genome-wide study of the COPT protein family was performed. To understand the major roles of the COPT gene family in Kandelia obovata (Ko), a genome-wide study identified four COPT genes in the Kandelia obovata genome for the first time. The domain and 3D structural variation, phylogenetic tree, chromosomal distributions, gene structure, motif analysis, subcellular localization, cis-regulatory elements, synteny and duplication analysis, and expression profiles in leaves and Cu were all investigated in this research. Structural and sequence investigations show that most KoCOPTs have three transmembrane domains (TMDs). According to phylogenetic research, these KoCOPTs might be divided into two subgroups, just like Populus trichocarpa. KoCOPT gene segmental duplications and positive selection pressure were discovered by universal analysis. According to gene structure and motif analysis, most KoCOPT genes showed consistent exon-intron and motif organization within the same group. In addition, we found five hormones and four stress- and seven light-responsive cis-elements in the KoCOPTs promoters. The expression studies revealed that all four genes changed their expression levels in response to copper (CuCl2) treatments. In summary, our study offers a thorough overview of the Kandelia obovata COPT gene family's expression pattern and functional diversity, making it easier to characterize each KoCOPT gene's function in the future.

A modified method for efficient RNA isolation from mangrove root tissues rich in secondary metabolites.

Secondary metabolites in mangroves often interfere with RNA extraction yielding poor concentration and quality, which is unsuitable for downstream applications. As existing protocols yielded low-quality RNA from root tissues of Kandelia candel (L.) Druce and Rhizophora mucronata Lam., an optimized method was developed for improving the quality and yield of RNA. Compared with three other methods, this optimized protocol gave better RNA yield and purity for both species. The absorbance ratios were ≥1.9 for A260/280 and A260/230, while RNA integrity number values ranged from 7.5 to 9.6. Results show that our modified method is efficient in obtaining high-quality RNA from mangrove roots and is suitable for downstream experiments such as cDNA synthesis, real-time quantitative PCR and next-generation sequencing.

Description and genomic characterization of Nocardioides bruguierae sp. nov., isolated from Bruguiera gymnorhiza.

Strains BSK12Z-3T and BSK12Z-4, two Gram-stain-positive, aerobic, non-spore-forming strains, were isolated from Shankou Mangrove Nature Reserve, Guangxi Zhuang Autonomous Region, China. The diagnostic diamino acid in the cell-wall peptidoglycan of strain BSK12Z-3T was LL-diaminopimelic acid and MK-8(H4) was the predominant menaquinone. The polar lipids comprised diphosphatidylglycerol (DPG), phosphatidylglycerol (PG) and phospholipid (PL). The major fatty acids was iso-C16:0. Phylogenetic analysis based on 16S rRNA gene sequences suggested that the two strains fell within the genus Nocardioides, appearing most closely related to Nocardioides ginkgobilobae KCTC 39594T (97.5-97.6 % sequence similarity) and Nocardioides marinus DSM 18248T (97.4-97.6 %). Genome-based phylogenetic analysis confirmed that strains BSK12Z-3T and BSK12Z-4 formed a distinct phylogenetic cluster within the genus Nocardioides. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of strains BSK12Z-3T, BSK12Z-4 with their most related species N. marinus DSM18248T were within the ranges of 77.2-77.3 % and 21.3-21.4 %, respectively, clearly indicated that strains BSK12Z-3T, BSK12Z-4 represented novel species. Strains BSK12Z-3T and BSK12Z-4 exhibited 99.9 % 16S rRNA gene sequence similarity. The ANI and dDDH values between the two strains were 97.8 % and 81.1 %, respectively, suggesting that they belong to the same species. However, DNA fingerprinting discriminated that they were not from one clonal origin. Based on phylogenomic and phylogenetic analyses coupled with phenotypic and chemotaxonomic characterizatons, strains BSK12Z-3T and BSK12Z-4 could be classified as a novel species of the genus Nocardioides, for which the name Nocardioides bruguierae sp. nov., is proposed. The type strain is BSK12Z-3T (=CGMCC 4.7709T = JCM 34554T).

Bacterial communities vary across populations and tissue type in red mangroves (Rhizophora mangle, Rhizophoraceae) along an expanding front.

Plant-associated microbial communities may be important sources of functional diversity and genetic variation that influence host evolution. Bacteria provide benefits for their hosts, yet in most plant systems we know little about their taxonomic composition or variation across tissues and host range. Red Mangrove (Rhizophora mangle L.) is a vital coastal plant species that is currently expanding poleward and with it, perhaps, its microbiome. We explored variability in bacterial communities across tissues, individuals, and populations. We collected samples from six sample types from 5 to 10 individuals at each of three populations and used 16S rRNA gene (iTag) sequencing to describe their bacterial communities. Core community members and dominant bacterial classes were determined for each sample type. Pairwise PERMANOVA of Bray-Curtis dissimilarity and Indicator Species Analysis revealed significant differences in bacterial communities between sample types and populations. We described the previously unexplored microbiome of the reproductive tissues of R. mangle. Populations and most sample types were associated with distinct communities. Bacterial communities associated with R. mangle are influenced by host geography and sample type. Our study provides a foundation for future work exploring the functional roles of these microbes and their relevance to biogeochemical cycling.

The Kandelia obovata transcription factor KoWRKY40 enhances cold tolerance in transgenic Arabidopsis.

BACKGROUND: WRKY transcription factors play key roles in plant development processes and stress response. Kandelia obovata is the most cold-resistant species of mangrove plants, which are the important contributors to coastal marine environment. However, there is little known about the WRKY genes in K. obovata. RESULTS: In this study, a WRKY transcription factor gene, named KoWRKY40, was identified from mangrove plant K. obovata. The full-length cDNA of KoWRKY40 gene was 1420 nucleotide bases, which encoded 318 amino acids. The KoWRKY40 protein contained a typical WRKY domain and a C2H2 zinc-finger motif, which were common signatures to group II of WRKY family. The three-dimensional (3D) model of KoWRKY40 was formed by one α-helix and five ß-strands. Evolutionary analysis revealed that KoWRKY40 has the closest homology with a WRKY protein from another mangrove plant Bruguiera gymnorhiza. The KoWRKY40 protein was verified to be exclusively located in nucleus of tobacco epidermis cells. Gene expression analysis demonstrated that KoWRKY40 was induced highly in the roots and leaves, but lowly in stems in K. obovata under cold stress. Overexpression of KoWRKY40 in Arabidopsis significantly enhanced the fresh weight, root length, and lateral root number of the transgenic lines under cold stress. KoWRKY40 transgenic Arabidopsis exhibited higher proline content, SOD, POD, and CAT activities, and lower MDA content, and H2O2 content than wild-type Arabidopsis under cold stress condition. Cold stress affected the expression of genes related to proline biosynthesis, antioxidant system, and the ICE-CBF-COR signaling pathway, including AtP5CS1, AtPRODH1, AtMnSOD, AtPOD, AtCAT1, AtCBF1, AtCBF2, AtICE1, AtCOR47 in KoWRKY40 transgenic Arabidopsis plants. CONCLUSION: These results demonstrated that KoWRKY40 conferred cold tolerance in transgenic Arabidopsis by regulating plant growth, osmotic balance, the antioxidant system, and ICE-CBF-COR signaling pathway. The study indicates that KoWRKY40 is an important regulator involved in the cold stress response in plants.

Chromosome-level genome assembly of Indian mangrove (Ceriops tagal) revealed a genome-wide duplication event predating the divergence of Rhizophoraceae mangrove species.

Mangrove ecosystems are unique, highly diverse, provide benefits to humans, and aid in coastal protection. The Indian mangrove, or spurred mangrove, [Ceriops tagal (Perr.) C. B. Rob.] is a member of the Rhizophoraceae family and is commonly found along the intertidal zones in tropical regions in Southeast Asia, southern Asia, and Africa. Here, we present the first high-quality reference genome assembly of the Ceriops species. A preliminary draft assembly, generated from the 10× Genomics linked-read library, was scaffolded using the proximity ligation chromatin contact mapping technique (Hi-C) to obtain a chromosome-scale assembly of 231,919,005 bases with an N50 length of 11,408,429 bases. The benchmarking universal single-copy orthologs (BUSCO) analysis revealed that C. tagal gene predictions recovered 95.8% of the highly conserved orthologs. Phylogenetic analyses suggested that C. tagal diverged from the last common ancestor of flat-leaf spurred mangrove [C. decandra (Griff.) Ding Hou] and C. zippeliana Blume ∼10.4 million yr ago (MYA), and the last common ancestor of genera Ceriops, Kandelia, and Rhizophora diverged from that of genus Bruguiera ∼49.4 MYA. In addition, our analysis of the transversion rate at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a recent whole-genome duplication in C. tagal. The STRUCTURE and principal component analyses illustrated that C. tagal individuals investigated in this study were the admixture of two subpopulations, the genetic background of which was influenced primarily by location. The availability of genomic and transcriptomic resources and biodiversity data reported in this work will be useful for future studies that may shed light on adaptive evolutions of mangrove species.

Assembly of a hybrid mangrove, Bruguiera hainesii, and its two ancestral contributors, Bruguiera cylindrica and Bruguiera gymnorhiza.

Mangroves are plants that live in tropical and subtropical coastal regions of the world, they are adapted to high salt environments and cyclic tidal flooding. Mangroves play important ecological roles, including acting as breeding grounds for many fish species and to prevent coastal erosion. The genomes of three mangrove species, Bruguiera gymnorhiza, Bruguiera cylindrica, and a hybrid of the two, Bruguiera hainesii were sequenced, assembled and annotated. The two progenitor species, B. gymnorhiza and B. cylindrica, were found to be highly similar to each other and sufficiently similar to B. parviflora to allow it to be used for reference based scaffolding to generate chromosome level scaffolds. The two subgenomes of B. hainesii were independently assembled and scaffolded. Analysis of B. hainesii confirms that it is a hybrid and the hybridisation event was estimated at 2.4 to 3.5 million years ago using a Bayesian Relaxed Molecular Clock approach.

A de novo reference assembly of the yellow mangrove Ceriops zippeliana genome.

Mangroves are of great ecological and economical importance, providing shelters for a wide range of species and nursery habitats for commercially important marine species. Ceriops zippeliana (yellow mangrove) belongs to Rhizophoraceae family and is commonly distributed in the tropical and subtropical coastal communities. In this study, we present a high-quality assembly of the C. zippeliana genome. We constructed an initial draft assembly of 240,139,412 bases with an N50 contig length of 564,761 bases using the 10x Genomics linked-read technology. This assembly was further scaffolded with RagTag using a chromosome-scale assembly of a closely related Ceriops species as a reference. The final assembly contained 243,228,612 bases with an N50 scaffold length of 10,559,178 Mb. The size of the final assembly was close to those estimated using DNA flow cytometry (248 Mb) and the k-mer distribution analysis (246 Mb). We predicted a total of 23,474 gene models and 21,724 protein-coding genes in the C. zippeliana genome, of which 16,002 were assigned gene ontology terms. We recovered 97.1% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs analysis. The phylogenetic analysis based on single-copy orthologous genes illustrated that C. zippeliana and Ceriops tagal diverged approximately 10.2 million years ago (MYA), and their last common ancestor and Kandelia obovata diverged approximately 29.9 MYA. The high-quality assembly of C. zippeliana presented in this work provides a useful genomic resource for studying mangroves' unique adaptations to stressful intertidal habitats and for developing sustainable mangrove forest restoration and conservation programs.

De Novo Reference Assembly of the Upriver Orange Mangrove (Bruguiera sexangula) Genome.

Upriver orange mangrove (Bruguiera sexangula) is a member of the most mangrove-rich taxon (Rhizophoraceae family) and is commonly distributed in the intertidal zones in tropical and subtropical latitudes. In this study, we employed the 10× Genomics linked-read technology to obtain a preliminary de novo assembly of the B. sexangula genome, which was further scaffolded to a pseudomolecule level using the Bruguiera parviflora genome as a reference. The final assembly of the B. sexangula genome contained 260 Mb with an N50 scaffold length of 11,020,310 bases. The assembly comprised 18 pseudomolecules (corresponding to the haploid chromosome number in B. sexangula), covering 204,645,832 bases or 78.6% of the 260-Mb assembly. We predicted a total of 23,978 protein-coding sequences, 17,598 of which were associated with gene ontology terms. Our gene prediction recovered 96.6% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. The chromosome-level assembly presented in this work provides a valuable genetic resource to help strengthen our understanding of mangroves' physiological and morphological adaptations to the intertidal zones.

A chromosome-scale reference genome assembly of yellow mangrove (Bruguiera parviflora) reveals a whole genome duplication event associated with the Rhizophoraceae lineage.

Mangrove forest ecosystems support a diverse flora and fauna of marine and terrestrial species and have important direct and indirect economic, ecological and social values to mankind. Yellow mangrove (Bruguiera parviflora) belongs to the Rhizophoraceae family and is widely distributed in the intertidal zones along sheltered coastal areas in tropical latitudes. Here, we present a high-quality, chromosome-level assembly of the B. parviflora genome. We employed the 10x Genomics linked-read technology to obtain a preliminary assembly, which was subsequently scaffolded using the long-range chromatin contact mapping technique (HiC) to obtain a final assembly containing 213,026,782 bases in 10,045 scaffolds with an N50 length of 10,906,948 bases. Our gene prediction recovered 96.5% of the highly conserved orthologues in the Embryophyta lineage based on the Benchmarking Universal Single-Copy Orthologues (BUSCO) analysis. We analysed the transversion rate at fourfold-degenerate sites from orthologous gene pairs and discovered evidence supporting a recent whole-genome duplication event in B. parviflora and other Rhizophoreae members. Comparative studies based on single-copy orthologous genes indicated that B. parviflora and Bruguiera gymnorrhiza diverged approximately 24.1 million years ago. The population structure analysis revealed that 63 B. parviflora accessions from different geographical regions in Thailand were an admixture of two subpopulations. The examination of alternative splicing events in B. parviflora showed that the most prevalent splicing mechanism was intron retention. This high-quality genome assembly together with the genetic diversity information obtained from the germplasm provide useful genomic resources for future studies on comparative phylogenetics and evolution of adaptive traits in mangrove species.

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.

The First De Novo Transcriptome Assembly and Transcriptomic Dynamics of the Mangrove Tree Rhizophora stylosa Griff. (Rhizophoraceae).

Mangroves are salt-tolerant plant species that grow in coastal saline water and are adapted to harsh environmental conditions. In this study, we de novo assembled and functionally annotated the transcriptome of Rhizophora stylosa, the widely distributed mangrove from the largest mangrove family (Rhizophoraceae). The final transcriptome consists of 200,491 unigenes with an average length, and N50 of 912.7 and 1334 base pair, respectively. We then compared the genome-wide expression profiles between the two morphologically distinct natural populations of this species growing under different levels of salinity depending on their distance from the ocean. Among the 200,491 unigenes, 40,253 were identified as differentially expressed between the two populations, while 15,741 and 24,512 were up- and down-regulated, respectively. Functional annotation assigned thousands of upregulated genes in saline environment to the categories related to abiotic stresses such as response to salt-, osmotic-, and oxidative-stress. Validation of those genes may contribute to a better understanding of adaptation in mangroves species. This study reported, for the first time, the transcriptome of R. stylosa, and the dynamic of it in response to salt stress and provided a valuable resource for elucidation of the molecular mechanism underlying the salt stress response in mangroves and other plants that live under stress.

Relict inland mangrove ecosystem reveals Last Interglacial sea levels.

Climatic oscillations during the Pleistocene played a major role in shaping the spatial distribution and demographic dynamics of Earth's biota, including our own species. The Last Interglacial (LIG) or Eemian Period (ca. 130 to 115 thousand years B.P.) was particularly influential because this period of peak warmth led to the retreat of all ice sheets with concomitant changes in global sea level. The impact of these strong environmental changes on the spatial distribution of marine and terrestrial ecosystems was severe as revealed by fossil data and paleogeographic modeling. Here, we report the occurrence of an extant, inland mangrove ecosystem and demonstrate that it is a relict of the LIG. This ecosystem is currently confined to the banks of the freshwater San Pedro Mártir River in the interior of the Mexico-Guatemala El Petén rainforests, 170 km away from the nearest ocean coast but showing the plant composition and physiognomy typical of a coastal lagoon ecosystem. Integrating genomic, geologic, and floristic data with sea level modeling, we present evidence that this inland ecosystem reached its current location during the LIG and has persisted there in isolation ever since the oceans receded during the Wisconsin glaciation. Our study provides a snapshot of the Pleistocene peak warmth and reveals biotic evidence that sea levels substantially influenced landscapes and species ranges in the tropics during this period.