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

Mitogenome selection shaped the terrestrial adaptation of Grapsidae (Decapoda: Brachyura)

The colonization of aquatic to terrestrial habitats by brachyuran crabs requires genetic innovations as well as morphological adaptations to adapt to terrestrial environments. The genetic basis of such adaptive evolution, however, is largely unknown. This study focuses on terrestrialization in Geograpsus (Grapsidae) the only highly terrestrial genus in this family, which represents a notable example of terrestrial adaptive radiation. Here, we sequenced the mitogenomes of two Geograpsus species and used the mitogenomes of 215 representative crabs to construct phylogenetic and time frameworks that we used to infer terrestrial origins and evolution. Using mitochondrial genomic data, we demonstrated that marine crab ancestors began to settle on land during the early Eocene. Ocean acidification, the Paleocene-Eocene Thermal Maximum (PETM), and mangrove expansion at that time may have driven the diversification and ecological expansion of these terrestrial crabs. Evolutionary analyses reveal strong positive selection signals on monophyletic lineages of Grapsidae, especially the terrestrial species of Geograpsus. Positively selected sites in functionally important regions of ND5 and ND4 may imply enhanced energy metabolism in Grapsidae compared to other crabs, and may have played an important role in their terrestrial adaptation. Overall, our work provides valuable resources and opportunities to reveal the adaptation of crabs to complex terrestrial environments.

Molecular mechanism of vivipary as revealed by the genomes of viviparous mangroves and non-viviparous relatives

Vivipary is a prominent feature of mangroves, allowing seeds to complete germination while attached to the mother plant, and equips propagules to endure and flourish in challenging coastal intertidal wetlands. However, vivipary-associated genetic mechanisms remain largely elusive. Genomes of two viviparous mangrove species and a non-viviparous inland relative were sequenced and assembled at the chromosome level. Comparative genomic analyses between viviparous and non-viviparous genomes revealed that DELAY OF GERMINATION 1 (DOG1) family genes (DFGs), the proteins from which are crucial for seed dormancy, germination, and reserve accumulation, are either lost or dysfunctional in the entire lineage of true viviparous mangroves but are present and functional in their inland, non-viviparous relatives. Transcriptome dynamics at key stages of vivipary further highlighted the roles of phytohormonal homeostasis, proteins stored in mature seeds, and proanthocyanidins in vivipary under conditions lacking DFGs. Population genomic analyses elucidate dynamics of syntenic regions surrounding the missing DFGs. Our findings demonstrated the genetic foundation of constitutive vivipary in Rhizophoraceae mangroves.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

From swamp to field: how genes from mangroves and its associates can enhance crop salinity tolerance

Salinity stress is a critical challenge in crop production and requires innovative strategies to enhance the salt tolerance of plants. Insights from mangrove species, which are renowned for their adaptability to high-salinity environments, provides valuable genetic targets and resources for improving crops. A significant hurdle in salinity stress is the excessive uptake of sodium ions (Na+) by plant roots, causing disruptions in cellular balance, nutrient deficiencies, and hampered growth. Specific ion transporters and channels play crucial roles in maintaining a low Na+/K+ ratio in root cells which is pivotal for salt tolerance. The family of high-affinity potassium transporters, recently characterized in Avicennia officinalis, contributes to K+ homeostasis in transgenic Arabidopsis plants even under high-salt conditions. The salt overly sensitive pathway and genes related to vacuolar-type H+-ATPases hold promise for expelling cytosolic Na+ and sequestering Na+ in transgenic plants, respectively. Aquaporins contribute to mangroves' adaptation to saline environments by regulating water uptake, transpiration, and osmotic balance. Antioxidant enzymes mitigate oxidative damage, whereas genes regulating osmolytes, such as glycine betaine and proline, provide osmoprotection. Mangroves exhibit increased expression of stress-responsive transcription factors such as MYB, NAC, and CBFs under high salinity. Moreover, genes involved in various metabolic pathways, including jasmonate synthesis, triterpenoid production, and protein stability under salt stress, have been identified. This review highlights the potential of mangrove genes to enhance salt tolerance of crops. Further research is imperative to fully comprehend and apply these genes to crop breeding to improve salinity resilience.

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.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Rapid diversification of grey mangroves (Avicennia marina) driven by geographic isolation and extreme environmental conditions in the Arabian Peninsula

Biological systems occurring in ecologically heterogeneous and spatially discontinuous habitats provide an ideal opportunity to investigate the relative roles of neutral and selective factors in driving lineage diversification. The grey mangroves (Avicennia marina) of Arabia occur at the northern edge of the species' range and are subject to variable, often extreme, environmental conditions, as well as historic large fluctuations in habitat availability and connectivity resulting from Quaternary glacial cycles. Here, we analyse fully sequenced genomes sampled from 19 locations across the Red Sea, the Arabian Sea and the Persian/Arabian Gulf (PAG) to reconstruct the evolutionary history of the species in the region and to identify adaptive mechanisms of lineage diversification. Population structure and phylogenetic analyses revealed marked genetic structure correlating with geographic distance and highly supported clades among and within the seas surrounding the Arabian Peninsula. Demographic modelling showed times of divergence consistent with recent periods of geographic isolation and low marine connectivity during glaciations, suggesting the presence of (cryptic) glacial refugia in the Red Sea and the PAG. Significant migration was detected within the Red Sea and the PAG, and across the Strait of Hormuz to the Arabian Sea, suggesting gene flow upon secondary contact among populations. Genetic-environment association analyses revealed high levels of adaptive divergence and detected signs of multi-loci local adaptation driven by temperature extremes and hypersalinity. These results support a process of rapid diversification resulting from the combined effects of historical factors and ecological selection and reveal mangrove peripheral environments as relevant drivers of lineage diversity.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Evolution of sea-surfing plant propagule as revealed by the genomes of Heritiera mangroves

Coastal forests, such as mangroves, protect much of the tropical and subtropical coasts. Long-distance dispersal via sea-surfing propagules is essential for coastal plants, but the genomic and molecular basis of sea-surfing plant propagule evolution remains unclear. Heritiera fomes and Heritiera littoralis are two coastal plants with typical buoyant fruits. We de novo sequenced and assembled their high-quality genomes. Our phylogenomic analysis indicates H. littoralis and H. fomes originated (at ~6.08 Mya) just before the start of Quaternary sea-level fluctuations. Whole-genome duplication occurred earlier, permitting gene copy gains in the two species. Many of the expanded gene families are involved in lignin and flavonoid biosynthesis, likely contributing to buoyant fruit emergence. It is repeatedly revealed that one duplicated copy to be under positive selection while the other is not. By examining H. littoralis fruits at three different developmental stages, we found that gene expression levels remain stable from young to intermediate. However, ~1000 genes are up-regulated and ~ 3000 genes are down-regulated as moving to mature. Particularly in fruit epicarps, the upregulation of WRKY12 and E2Fc likely constrains the production of p-Coumaroyl-CoA, the key internal substrate for lignin biosynthesis. Hence, to increase fruit impermeability, methylated lignin biosynthesis is shut down by down-regulating the genes CCoAOMT, F5H, COMT, and CSE, while unmethylated lignins are preferentially produced by upregulating CAD and CCR. Similarly, cutin polymers and cuticular waxes accumulate with high levels before maturation in epicarps. Overall, our genome assemblies and analyses uncovered the genomic evolution and temporal transcriptional regulation of sea-surfing propagule.

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.

Marine-montane transitions coupled with gill and genetic convergence in extant crustacean

Marine-terrestrial transition represents an important aspect of organismal evolution that requires numerous morphological and genetic innovations and has been hypothesized to be caused by geological changes. We used talitrid crustaceans with marine-coastal-montane extant species at a global scale to investigate the marine origination and terrestrial adaptation. Using genomic data, we demonstrated that marine ancestors repeatedly colonized montane terrestrial habitats during the Oligocene to Miocene. Biological transitions were well correlated with plate collisions or volcanic island formation, and top-down cladogenesis was observed on the basis of a positive relationship between ancestral habitat elevation and divergence time for montane lineages. We detected convergent variations of convoluted gills and convergent evolution of SMC3 associated with montane transitions. Moreover, using CRISPR-Cas9 mutagenesis, we proposed that SMC3 potentially regulates the development of exites, such as talitrid gills. Our results provide a living model for understanding biological innovations and related genetic regulatory mechanisms associated with marine-terrestrial transitions.

The role of oceanic currents in the dispersal and connectivity of the mangrove Rhizophora mangle on the Southwest Atlantic region

Dispersal is a crucial mechanism to living beings, allowing them to reach new resources such that populations and species can occupy new environments. However, directly observing the dispersal mechanisms of widespread species can be costly or even impractical, which is the case for mangrove trees. The influence of ocean currents on mangrove dispersal is increasingly evident; however, few studies mechanistically relate the patterns of population distribution with the dispersal by oceanic currents under an integrated framework. Here, we evaluate the role of oceanic currents on connectivity of Rhizophora mangle along the Southwest Atlantic. We inferred population genetic structure and migration rates, simulated the displacement of propagules and tested our hypotheses with Mantel tests and redundancy analysis. We observed populations structured in two major groups, north and south, which is corroborated by other studies with Rhizophora and other coastal plants. Inferred recent migration rates do not indicate ongoing gene flow between sites. Conversely, long-term migration rates were low across groups and contrasting dispersal patterns within each one, which is consistent with long-distance dispersal events. Our hypothesis tests suggest that both isolation by distance and isolation by oceanography (derived from the oceanic currents) can explain the neutral genetic variation of R. mangle in the region. Our findings expand current knowledge of mangrove connectivity and highlight how the association of molecular methods with oceanographic simulations improve the interpretation of the dispersal process. This integrative approach is a cost- and time-efficient strategy to include dispersal and connectivity data into marine protected areas planning and management.

Chloroplast genome analysis of three Acanthus species reveal the adaptation of mangrove to intertidal habitats

Acanthus is a distinctive genus that covers three species with different ecological niches including Acanthus mollis (arid terrestrial), Acanthus leucostachyus (damp forest) and Acanthus ilicifolius (coastal intertidal). It is an intriguing question how these species evolved from terrestrial to coastal intertidal. In the present study, we assembled chloroplast genomes of A. ilicifolius, A. leucostachyus and A. mollis, which exhibited typical quadripartite structures. The sizes were 150,758, 154,686 and 150,339 bp that comprised a large single copy (LSC, 82,963, 86,461 and 82,612 bp), a small single copy (SSC, 17,191, 17,511 and 17,019 bp), and a pair of inverted repeats (IRs, 25,302, 25,357 and 25,354 bp), respectively. Gene annotation revealed that A. ilicifolius, A. leucostachyus and A. mollis contained 113, 112 and 108 unique genes, each of which contained 79, 79 and 74 protein-coding genes, 30, 29 and 30 tRNAs, and 4 rRNA genes, respectively. Differential gene analysis revealed plenty of ndhs gene deletions in the terrestrial plant A. mollis. Nucleotide diversity analysis showed that the psbK, ycf1, ndhG, and rpl22 have the highest nucleotide variability. Compared to A. leucostachyus and A. mollis, seven genes in A. ilicifolius underwent positive selection. Among them, the atpF gene showed a strong positive selection throughout terrestrial to marine evolution and was important for adaptation to coastal intertidal habitats. Phylogenetic analysis indicated that A. ilicifolius has a closer genetic relationship with A. leucostachyus than A. mollis which further confirmed the evolutionary direction of Acanthus going from terrestrial to coastal intertidal zones.

Expansion and adaptive evolution of the WRKY transcription factor family in Avicennia mangrove trees

Mangroves are adapted to intertidal zones, which present extreme environmental conditions. WRKYs are among the most prominent transcription factors (TFs) in higher plants and act through various interconnected networks to regulate responses to multiple abiotic stressors. Here, based on omic data, we investigated the landscape and evolutionary patterns of WRKYs in the main mangrove genus Avicennia. We found that both the number and the proportion of TFs and WRKYs in Avicennia species exceeded their inland relatives, indicating a significant expansion of WRKYs in Avicennia. We identified 109 WRKY genes in the representative species Avicennia marina. Comparative genomic analysis showed that two recent whole-genome duplication (WGD) events played a critical role in the expansion of WRKYs, and 88% of Avicennia marina WRKYs (AmWRKYs) have been retained following these WGDs. Applying comparative transcriptomics on roots under experimental salt gradients, we inferred that there is high divergence in the expression of WGD-retained AmWRKYs. Moreover, we found that the expression of 16 AmWRKYs was stable between freshwater and moderately saline water but increased when the trees were exposed to high salinity. In particular, 14 duplicates were retained following the two recent WGD events, indicating potential neo- and sub-functionalization. We also found that WRKYs could interact with other upregulated genes involved in signalling pathways and natural antioxidant biosynthesis to enhance salt tolerance, contributing to the adaptation to intertidal zones. Our omic data of the WRKY family in A. marina broadens the understanding of how a TF family relates to the adaptive evolution of mangroves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00177-y.

Genome wide identification and characterization of MATE family genes in mangrove plants

Multidrug and Toxic Compound Extrusion (MATE) proteins are essential transporters that extrude metabolites and participate in plant development and cellular detoxification. MATE transporters, which play crucial roles in the survival of mangrove plants under highly challenged environments, by specialized salt extrusion mechanisms, are mined from their genomes and reported here for the first time. Through homology search and domain prediction in the genome assemblies of Avicennia marina, Bruguiera sexangula, Ceriops zippeliana, Kandelia obovata, Rhizophora apiculata and Ceriops tagal, 74, 68, 66, 66, 63 and 64 MATE proteins, respectively were identified. The phylogenetic analysis divided the identified proteins into five major clusters and following the clustering pattern of the functionally characterized proteins, functions of the transporters in each cluster were predicted. Amino acid sequences, exon-intron structure, motif details and subcellular localization pattern for all the 401 proteins are described. The custom designed repeat masking libraries generated for each of these genomes, which will be of extensive use for the researchers worldwide, are also provided in this paper. This is the first study on the MATE genes in mangroves and the results provide comprehensive information on the molecular mechanisms enabling the survival of mangroves under hostile conditions.

Assessment of the Genetic Diversity and Population Structure of Rhizophora mucronata along Coastal Areas in Thailand

Unique and biodiverse, mangrove ecosystems provide humans with benefits and contribute to coastal protection. Rhizophora mucronata, a member of the Rhizophoraceae family, is prevalent in the mangrove forests of Thailand. R. mucronata's population structure and genetic diversity have received scant attention. Here, we sequenced the entire genome of R. mucronata using 10× Genomics technology and obtained an assembly size of 219 Mb with the N50 length of 542,540 bases. Using 2857 single nucleotide polymorphism (SNP) markers, this study investigated the genetic diversity and population structure of 80 R. mucronata accessions obtained from the mangrove forests in Thailand. The genetic diversity of R. mucronata was moderate (I = 0.573, Ho = 0.619, He = 0.391). Two subpopulations were observed and confirmed from both population structure and principal component analysis (PCA). Analysis of molecular variance (AMOVA) showed that there was more variation within populations than between them. Mean pairwise genetic differentiation (FST = 0.09) showed that there was not much genetic difference between populations. Intriguingly, the predominant clustering pattern in the R. mucronata population did not correspond to the Gulf of Thailand and the Andaman Sea, which are separated by the Malay Peninsula. Several factors could have influenced the R. mucronata genetic pattern, such as hybridization and anthropogenic factors. This research will provide important information for the future conservation and management of R. mucronata in Thailand.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Evolution of woody plants to the land-sea interface - The atypical genomic features of mangroves with atypical phenotypic adaptation

How plants adapt and diverge in extreme environments is a key question of plant evolution and ecology. Mangrove invasion of intertidal environments is facilitated by adaptive phenotypes such as aerial roots, salt-secreting leaf, and viviparity, and genomic mechanisms including whole genome duplication and transposable element number reduction. However, a number of mangroves lack these typical phenotypes. The question we ask is whether these phenotypically atypical mangroves also have distinct genomic features? The sibling mangrove species Lumnitzera littorea and Lumnitzera racemosa provide a model to study this question. We sequenced and assembled their genomes to chromosome level, together with a closely related species Combretum micranthum. While most mangroves have small genomes, the genomes of both Lumnitzera species are large (1443 and 1317 Mb) and carry a high proportion of repeat sequences (~75%). Moreover, Lumnitzera species have not undergone post-gamma whole-genome duplications. Their genome size increased mainly due to the expansion of repeat sequences in their ancestors. However, Lumnitzera genomes have reduced transposable elements by constraining the proliferation of new LTR-RTs. Meanwhile, the two species have more gene families contracted than expanded, and some gene families with reversed size change may underlie their differentiation in root morphology and local distribution. We identified 86 chromosomal inversions, five of which are measured between 6.5 and 12.8 megabases. A number of genes located in these inversions function in pigment biosynthesis, a process likely involved in flower colour differentiation between the Lumnitzera species. We conclude that the mangroves with atypical phenotypes also have atypical genomic evolution.

In silico analysis of NAC gene family in the mangrove plant Avicennia marina provides clues for adaptation to intertidal habitats

NAC (NAM, ATAF1/2, CUC2) transcription factors (TFs) constitute a plant-specific gene family. It is reported that NAC TFs play important roles in plant growth and developmental processes and in response to biotic/abiotic stresses. Nevertheless, little information is known about the functional and evolutionary characteristics of NAC TFs in mangrove plants, a group of species adapting coastal intertidal habitats. Thus, we conducted a comprehensive investigation for NAC TFs in Avicennia marina, one pioneer species of mangrove plants. We totally identified 142 NAC TFs from the genome of A. marina. Combined with NAC proteins having been functionally characterized in other organisms, we built a phylogenetic tree to infer the function of NAC TFs in A. marina. Gene structure and motif sequence analyses suggest the sequence conservation and transcription regulatory regions-mediated functional diversity. Whole-genome duplication serves as the driver force to the evolution of NAC gene family. Moreover, two pairs of NAC genes were identified as positively selected genes of which AmNAC010/040 may be imposed on less constraint toward neofunctionalization. Quite a few stress/hormone-related responsive elements were found in promoter regions indicating potential response to various external factors. Transcriptome data revealed some NAC TFs were involved in pneumatophore and leaf salt gland development and response to salt, flooding and Cd stresses. Gene co-expression analysis found a few NAC TFs participates in the special biological processes concerned with adaptation to intertidal environment. In summary, this study provides detailed functional and evolutionary information about NAC gene family in mangrove plant A. marina and new perspective for adaptation to intertidal habitats.

Uninterrupted embryonic growth leading to viviparous propagule formation in woody mangrove

Vivipary is a rare sexual reproduction phenomenon where embryos germinate directly on the maternal plants. However, it is a common genetic event of woody mangroves in the Rhizophoraceae family. The ecological benefits of vivipary in mangroves include the nurturing of seedlings in harsh coastal and saline environments, but the genetic and molecular mechanisms of vivipary remain unclear. Here we investigate the viviparous embryo development and germination processes in mangrove Kandelia obovata by a transcriptomic approach. Many key biological pathways and functional genes were enriched in different tissues and stages, contributing to vivipary. Reduced production of abscisic acid set a non-dormant condition for the embryo to germinate directly. Genes involved in the metabolism of and response to other phytohormones (gibberellic acid, brassinosteroids, cytokinin, and auxin) are expressed precociously in the axis of non-vivipary stages, thus promoting the embryo to grow through the seed coat. Network analysis of these genes identified the central regulatory roles of LEC1 and FUS3, which maintain embryo identity in Arabidopsis. Moreover, photosynthesis related pathways were significantly up-regulated in viviparous embryos, and substance transporter genes were highly expressed in the seed coat, suggesting a partial self-provision and maternal nursing. We conclude that the viviparous phenomenon is a combinatorial result of precocious loss of dormancy and enhanced germination potential during viviparous seed development. These results shed light on the relationship between seed development and germination, where the continual growth of the embryo replaces a biphasic phenomenon until a mature propagule is established.

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.

Extensive gene flow in secondary sympatry after allopatric speciation

In the conventional view, species are separate gene pools delineated by reproductive isolation (RI). In an alternative view, species may also be delineated by a small set of 'speciation genes' without full RI, a view that has gained broad acceptance. A recent survey, however, suggested that the extensive literature on 'speciation with gene flow' is mostly (if not all) about exchanges in the early stages of speciation. There is no definitive evidence that the observed gene flow actually happened after speciation is completed. Here, we wish to know whether 'good species' (defined by the 'secondary sympatry' test) do continue to exchange genes and, importantly, under what conditions such exchanges can be observed. De novo whole-genome assembly and re-sequencing of individuals across the range of two closely related mangrove species (Rhizophora mucronata and R. stylosa) reveal the genomes to be well delineated in allopatry. They became sympatric in northeastern Australia but remain distinct species. Nevertheless, their genomes harbor ∼4000-10 000 introgression blocks averaging only about 3-4 Kb. These fine-grained introgressions indicate continual gene flow long after speciation as non-introgressable 'genomic islets,' ∼1.4 Kb in size, often harbor diverging genes of flower or gamete development. The fine-grained introgression in secondary sympatry may help settle the debate about sympatric vs. micro-allopatric speciation. In conclusion, true 'good species' may often continue to exchange genes but the opportunity for detection is highly constrained.

Establishing the relationship between non-human primates and mangrove forests at the global, national, and local scales

Global and spatially explicit information about the interaction between habitat and wildlife species is critical to enhancing conservation efforts. Despite the recognized importance of mangrove forests to non-human primates, the relationship between the two lacks understanding. To counter this, we created the MangPrim-21 database to map and measure the locations of interactions between all non-human primates and all mangrove forests globally. We report our findings across the global, national, and local scales for all inventoried non-human primates and all inventoried mangrove forests. Globally, we find that half of all non-primates potentially use mangrove forests, and more than half of the global mangrove forest falls within the delineated range of at least one non-human primate species. Nationally, we find that Indonesia, Madagascar, Brazil, Cameroon, and Malaysia likely have the most non-human primate and mangrove forest interactions. At the subnational level, we find that several discrete locations in Kalimantan are critical to both mangrove forests and non-human primates. The MangPrim-21 database provides a globally consistent and locally applicable database of non-human primate and mangrove forest interactions. The results presented have broader implications for non-human primate and mangrove conservation and global actions to protect both. Additionally, our results raise questions about the idea that non-human primates primarily use mangrove forests as a refuge from human encroachment and habitat degradation.

Assessment of the Genetic Diversity and Population Structure of Rhizophora apiculata Blume (Rhizophoraceae) in Thailand

Rhizophora apiculata is one of the most widespread and economically important mangrove trees in the Indo-West Pacific region. Knowledge of the genetic variation of R. apiculata in Thailand is limited. Here, we generated a whole-genome sequence of R. apiculata using the 10× Genomics technology. R. apiculata genome assembly was 230.47 Mb. Based on its genome, 2640 loci of high-quality biallelic SNPs were identified from 82 R. apiculata accessions collected from 17 natural mangrove forests in Thailand to assess the genetic diversity and population structure among them. A moderate level of genetic diversity of R. apiculata was observed. The average observed heterozygosity (Ho = 0.48) was higher than the average expected heterozygosity (He = 0.36). Two subpopulations were observed and confirmed from three approaches: population structure, PCA, and phylogenetic analyses. They corresponded to the Gulf of Thailand and the Andaman Sea separated by the Malay Peninsula. AMOVA analyses indicated that genetic variation was attributable to 76.22% within populations and 23.78% among populations. A high level of genetic differentiation between the two subpopulations (FST = 0.24, p < 0.001) was observed. This study evaluated the genetic diversity and population structure of R. apiculata, providing useful information for sustainable mangrove management in Thailand.

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.

Adaptation to a new environment with pre-adaptive genomic features - Evidence from woody plants colonizing the land-sea interface

Adaptation to new environments is a key evolutionary process which presumably involves complex genomic changes. Mangroves, a collection of approximately 80 woody plants that have independently invaded intertidal zones >20 times, are ideal for studying this process. We assembled near-chromosome-scale genomes of three Xylocarpus species as well as an outgroup species using single-molecule real-time sequencing. Phylogenomic analysis reveals two separate lineages, one with the mangrove Xylocarpus granatum and the other comprising a mangrove Xylocarpus moluccensis and a terrestrial Xylocarpus rumphii. In conjunction with previous studies, we identified several genomic features associated with mangroves: (i) signals of positive selection in genes related to salt tolerance and root development; (ii) genome-wide elevated ratios of non-synonymous to synonymous substitution relative to terrestrial relatives; and (iii) active elimination of long terminal repeats. These features are found in the terrestrial X. rumphii in addition to the two mangroves. These genomic features, not being strictly mangrove-specific, are hence considered pre-adaptive. We infer that the coastal but non-intertidal habitat of X. rumphii may have predisposed the common ancestor to invasion of true mangrove habitats. Other features including the preferential retention of duplicated genes and intolerance to pseudogenization are not found in X. rumphii and are likely true adaptive features in mangroves. In conclusion, by studying adaptive shift and partial shifts among closely related species, we set up a framework to study genomic features that are acquired at different stages of the pre-adaptation and adaptation to new environments.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Evolution of coastal forests based on a full set of mangrove genomes

Genomic studies are now poised to explore whole communities of species. The ~70 species of woody plants that anchor the coastal ecosystems of the tropics, collectively referred to as mangroves, are particularly suited to this exploration. In this study, we de novo sequenced the whole genomes of 32 mangroves, which we combined with other sequences of 30 additional species, comprising almost all mangroves globally. These community-wide genomic data will be valuable for ecology, evolution and biodiversity research. While the data revealed 27 independent origins of mangroves, the total phylogeny shows only modest increases in species number, even in coastal areas of active speciation, suggesting that mangrove extinction is common. A possible explanation for common extinction is the frequent sea-level rises and falls (SLRs and SLFs) documented in the geological record. Indeed, near-extinctions of species with extremely small population size (N) often happened during periods of rapid SLR, as revealed by the genome-wide heterozygosity of almost all mangroves. Reduction in N has possibly been further compounded by population fragmentation and the subsequent accumulation of deleterious mutations, thus pushing mangroves even closer to extinction. Crucially, the impact of the next SLR will be exacerbated by human encroachment into these mangrove habitats, potentially altering the ecosystems of tropical coasts irreversibly.

Identification, Characterization and Function of Orphan Genes Among the Current Cucurbitaceae Genomes

Orphan genes (OGs) that are missing identifiable homologs in other lineages may potentially make contributions to a variety of biological functions. The Cucurbitaceae family consists of a wide range of fruit crops of worldwide or local economic significance. To date, very few functional mechanisms of OGs in Cucurbitaceae are known. In this study, we systematically identified the OGs of eight Cucurbitaceae species using a comparative genomics approach. The content of OGs varied widely among the eight Cucurbitaceae species, ranging from 1.63% in chayote to 16.55% in wax gourd. Genetic structure analysis showed that OGs have significantly shorter protein lengths and fewer exons in Cucurbitaceae. The subcellular localizations of OGs were basically the same, with only subtle differences. Except for aggregation in some chromosomal regions, the distribution density of OGs was higher near the telomeres and relatively evenly distributed on the chromosomes. Gene expression analysis revealed that OGs had less abundantly and highly tissue-specific expression. Interestingly, the largest proportion of these OGs was significantly more tissue-specific expressed in the flower than in other tissues, and more detectable expression was found in the male flower. Functional prediction of OGs showed that (1) 18 OGs associated with male sterility in watermelon; (2) 182 OGs associated with flower development in cucumber; (3) 51 OGs associated with environmental adaptation in watermelon; (4) 520 OGs may help with the large fruit size in wax gourd. Our results provide the molecular basis and research direction for some important mechanisms in Cucurbitaceae species and domesticated crops.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments

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

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 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.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

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.

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.

Mangrove diversity is more than fringe deep

Mangroves form coastal tropical forests in the intertidal zone and are an important component of shoreline protection. In comparison to other tropical forests, mangrove stands are thought to have relatively low genetic diversity with population genetic structure gradually increasing with distance along a coastline. We conducted genetic analyses of mangrove forests across a range of spatial scales; within a 400 m² parcel comprising 181 Rhizophora mangle (red mangrove) trees, and across four sites ranging from 6-115 km apart in Honduras. In total, we successfully genotyped 269 R. mangle trees, using a panel of 677 SNPs developed with 2b-RAD methodology. Within the 400 m² parcel, we found two distinct clusters with high levels of genetic differentiation (F_ST = 0.355), corresponding to trees primarily located on the seaward fringe and trees growing deeper into the forest. In contrast, there was limited genetic differentiation (F_ST = 0.027-0.105) across the sites at a larger scale, which had been predominantly sampled along the seaward fringe. Within the 400 m² parcel, the cluster closest to the seaward fringe exhibited low genetic differentiation (F_ST = 0.014-0.043) with the other Honduran sites, but the cluster further into the forest was highly differentiated from them (F_ST = 0.326-0.414). These findings contradict the perception that genetic structure within mangroves forests occurs mainly along a coastline and highlights that there is greater genetic structure at fine spatial scales.

Genetic and molecular mechanisms underlying mangrove adaptations to intertidal environments

Mangroves are halophytic plants belonging to diverse angiosperm families that are adapted to highly stressful intertidal zones between land and sea. They are special, unique, and one of the most productive ecosystems that play enormous ecological roles and provide a large number of benefits to the coastal communities. To thrive under highly stressful conditions, mangroves have innovated several key morphological, anatomical, and physio-biochemical adaptations. The evolution of the unique adaptive modifications might have resulted from a host of genetic and molecular changes and to date we know little about the nature of these genetic and molecular changes. Although slow, new information has accumulated over the last few decades on the genetic and molecular regulation of the mangrove adaptations, a comprehensive review on it is not yet available. This review provides up-to-date consolidated information on the genetic, epigenetic, and molecular regulation of mangrove adaptive traits.

Whole-Genome Sequencing of Acer catalpifolium Reveals Evolutionary History of Endangered Species

Acer catalpifolium is an endangered species restricted to remote localities of West China. Understanding the genomic content and evolution of A. catalpifolium is essential to conservation efforts of this rare and ecologically valuable plant. Here, we report a high-quality genome of A. catalpifolium consisting of ∼654 Mbp and ∼35,132 protein-coding genes. We detected 969 positively selected genes in two Acer genomes compared with four other eudicots, 65 of which were transcription factors. We hypothesize that these positively selected mutations in transcription factors might affect their function and thus contribute to A. catalpifolium’s decline-type population. We also identified 179 significantly expanded gene families compared with 12 other eudicots, some of which are involved in stress responses, such as the FRS–FRF family. We inferred that A. catalpifolium has experienced gene family expansions to cope with environmental stress in its evolutionary history. Finally, 109 candidate genes encoding key enzymes in the lignin biosynthesis pathway were identified in A. catalpifolium; of particular note were the large range and high copy number of cinnamyl alcohol dehydrogenase genes. The chromosome-level genome of A. catalpifolium presented here may serve as a fundamental genomic resource for better understanding endangered Acer species, informing future conservation efforts.

Global diversification and evolutionary history of onchidiid slugs (Gastropoda, Pulmonata)

Many marine species are specialized to specific parts of a habitat. In a mangrove forest, for instance, species may be restricted to the mud surface, the roots and trunks of mangrove trees, or rotting logs, which can be regarded as distinct microhabitats. Shifts to new microhabitats may be an important driver of sympatric speciation. However, the evolutionary history of these shifts is still poorly understood in most groups of marine organisms, because it requires a well-supported phylogeny with relatively complete taxon sampling. Onchidiid slugs are an ideal case study for the evolutionary history of habitat and microhabitat shifts because onchidiid species are specialized to different tidal zones and microhabitats in mangrove forests and rocky shores, and the taxonomy of the family in the Indo-West Pacific has been recently revised in a series of monographs. Here, DNA sequences for onchidiid species from the North and East Pacific, the Caribbean, and the Atlantic are used to reconstruct phylogenetic relationships among Onchidella species, and are combined with new data for Indo-West Pacific species to reconstruct a global phylogeny of the family. The phylogenetic relationships of onchidiid slugs are reconstructed based on three mitochondrial markers (COI, 12S, 16S) and three nuclear markers (28S, ITS2, H3) and nearly complete taxon sampling (all 13 genera and 62 of the 67 species). The highly-supported phylogeny presented here suggests that ancestral onchidiids most likely lived in the rocky intertidal, and that a lineage restricted to the tropical Indo-West Pacific colonized new habitats, including mudflats, mangrove forests, and high-elevation rainforests. Many onchidiid species in the Indo-West Pacific diverged during the Miocene, around the same time that a high diversity of mangrove plants appears in the fossil record, while divergence among Onchidella species occurred earlier, likely beginning in the Eocene. It is demonstrated that ecological specialization to microhabitats underlies the divergence between onchidiid genera, as well as the diversification through sympatric speciation in the genera Wallaconchis and Platevindex. The geographic distributions of onchidiid species also indicate that allopatric speciation played a key role in the diversification of several genera, especially Onchidella and Peronia. The evolutionary history of several morphological traits (penial gland, rectal gland, dorsal eyes, intestinal loops) is examined in relation to habitat and microhabitat evolutionary transitions and that the rectal gland of onchidiids is an adaptation to high intertidal and terrestrial habitats.

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.

Genomic insights into molecular adaptation to intertidal environments in the mangrove Aegiceras corniculatum

Mangroves have colonised extreme intertidal environments characterised by high salinity, hypoxia and other abiotic stresses. Aegiceras corniculatum, a pioneer mangrove species that has evolved two specialised adaptive traits (salt secretion and crypto-vivipary) is an attractive ecological model to investigate molecular mechanisms underlying adaptation to intertidal environments. We assembled de novo a high-quality reference genome of A. corniculatum and performed comparative genomic and transcriptomic analyses to investigate molecular mechanisms underlying adaptation to intertidal environments. We provide evidence that A. corniculatum experienced a whole-genome duplication (WGD) event c. 35 Ma. We infer that maintenance of cellular environmental homeostasis is an important adaptive process in A. corniculatum. The 14-3-3 and H⁺ -ATPase protein-coding genes, essential for the salt homeostasis, were preferentially retained after the recent WGD event. Using comparative transcriptomics, we show that genes upregulated under high-salt conditions are involved in salt transport and ROS scavenging. We also found that all homologues of DELAY OF GERMINATION1 (DOG1) had lost their heme-binding ability in A. corniculatum, and that this may contribute to crypto-vivipary. Our study provides insight into the genomic correlates of phenotypic adaptation to intertidal environments. This could contribute not only within the genomics community, but also to the field of plant evolution.

Is speciation an unrelenting march to reproductive isolation?

Speciation is often portrayed as an "incomplete" or "incipient" process if two groups of organisms, technically distinguishable either by morphology or genetics, can exchange genes. The ultimate outcome of diversification, given this perspective, is complete reproductive isolation. But an increasing amount of evidence suggests that speciation is rarely complete and inter-fertility between different taxonomically accepted species is consistently maintained. In this issue of Molecular Ecology, Linan et al. (2021) provide results that bridge evolutionary processes from populations to phylogenies that indicate suites of closely related tree species in the Mascarene Islands actively exchange genes, evolving as a nested set of syngameons with a hierarchical pattern of interfertility. The deep insight into diversification provided by this study is particularly powerful because of the genomic scale of the data and the complete taxonomic sampling of an island clade evolving in situ. The prevalence of syngameon dynamics in a broad range of organisms indicates that we should adopt a fluid and comprehensive approach to defining evolutionary units for conservation and research. We should move beyond focusing on single endangered species in evolutionary and ecological isolation from other species but consider the entire network of potentially interfertile species and the potential for future adaptation and innovation, particularly in a human dominated world.

Identification, characterization and expression analysis of lineage-specific genes within mangrove species Aegiceras corniculatum

Lineage-specific genes (LSGs) are the genes that have no recognizable homology to any sequences in other species, which are important drivers for the generation of new functions, phenotypic changes, and facilitating species adaptation to environment. Aegiceras corniculatum is one of major mangrove plant species adapted to waterlogging and saline conditions, and the exploration of aegiceras-specific genes (ASGs) is important to reveal its adaptation to the harsh environment. Here, we performed a systematic analysis on ASGs, focusing on their sequence characterization, origination and expression patterns. Our results reveal that there are 4823 ASGs in the genome, approximately 11.84% of all protein-coding genes. High proportion (45.78%) of ASGs originate from gene duplication, and the time of gene duplication of ASGs is consistent with the timing of two genome-wide replication (WGD) events that occurred in A. corniculatum, and also coincides with a short period of global warming during the Paleocene-Eocene Maximum (PETM, 55.5 million years ago). Gene structure analysis showed that ASGs have shorter protein lengths, fewer exons, and higher isoelectric point. Expression patterns analysis showed that ASGs had low levels of expression and more tissue-specific expression. Weighted gene co-expression network analysis (WGCNA) revealed that 86 ASGs co-expressed gene modules were primarily involved in pathways related to adversity stress, including plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, peroxisome and pentose phosphate pathway. This study provides a comprehensive analysis of the characteristics and potential functions of ASGs and identifies key candidate genes, which will contribute to the subsequent further investigation of the adaptation of A. corniculatum to intertidal coastal wetland habitats.

A reference-grade genome identifies salt-tolerance genes from the salt-secreting mangrove species Avicennia marina

Water scarcity and salinity are major challenges facing agriculture today, which can be addressed by engineering plants to grow in the boundless seawater. Understanding the mangrove plants at the molecular level will be necessary for developing such highly salt-tolerant agricultural crops. With this objective, we sequenced the genome of a salt-secreting and extraordinarily salt-tolerant mangrove species, Avicennia marina, that grows optimally in 75% seawater and tolerates >250% seawater. Our reference-grade ~457 Mb genome contains 31 scaffolds corresponding to its chromosomes. We identified 31,477 protein-coding genes and a salinome consisting of 3246 salinity-responsive genes and homologs of 614 experimentally validated salinity tolerance genes. The salinome provides a strong foundation to understand the molecular mechanisms of salinity tolerance in plants and breeding crops suitable for seawater farming.

Global phylogeography of a pantropical mangrove genus Rhizophora

Rhizophora is a key genus for revealing the formation process of the pantropical distribution of mangroves. In this study, in order to fully understand the historical scenario of Rhizophora that achieved pantropical distribution, we conducted phylogeographic analyses based on nucleotide sequences of chloroplast and nuclear DNA as well as microsatellites for samples collected worldwide. Phylogenetic trees suggested the monophyly of each AEP and IWP lineages respectively except for R. samoensis and R. × selala. The divergence time between the two lineages was 10.6 million years ago on a dated phylogeny, and biogeographic stochastic mapping analyses supported these lineages separated following a vicariant event. These data suggested that the closure of the Tethys Seaway and the reduction in mangrove distribution followed by Mid-Miocene cooling were key factors that caused the linage diversification. Phylogeographic analyses also suggested the formation of the distinctive genetic structure at the AEP region across the American continents around Pliocene. Furthermore, long-distance trans-pacific dispersal occurred from the Pacific coast of American continents to the South Pacific and formed F1 hybrid, resulting in gene exchange between the IWP and AEP lineages after 11 million years of isolation. Considering the phylogeny and phylogeography with divergence time, a comprehensive picture of the historical scenario behind the pantropical distribution of Rhizophora is updated.

SMARTdenovo: a de novo assembler using long noisy reads

Long-read single-molecule sequencing has revolutionized de novo genome assembly and enabled the automated reconstruction of reference-quality genomes. It has also been widely used to study structural variants, phase haplotypes and more. Here, we introduce the assembler SMARTdenovo, a single-molecule sequencing (SMS) assembler that follows the overlap-layout-consensus (OLC) paradigm. SMARTdenovo (RRID: SCR_017622) was designed to be a rapid assembler, which, unlike contemporaneous SMS assemblers, does not require highly accurate raw reads for error correction. It has performed well in the evaluation of congeneric assemblers and has been successfully users for various assembly projects. It is compatible with Canu for assembling high-quality genomes, and several of the assembly strategies in this program have been incorporated into subsequent popular assemblers. The assembler has been in use since 2015; here we provide information on the development of SMARTdenovo and how to implement its algorithms into current projects.