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

Growth, Morphology and Respiratory Cost Responses to Salinity in the Mangrove Plant Rhizophora Stylosa Depend on Growth Temperature

Mangrove plants, which have evolved to inhabit tidal flats, may adjust their physiological and morphological traits to optimize their growth in saline habitats. Furthermore, the confined distribution of mangroves within warm regions suggests that warm temperature is advantageous to their growth in saline environments. We analyzed growth, morphology and respiratory responses to moderate salinity and temperature in a mangrove species, Rhizophora stylosa. The growth of R. stylosa was accelerated in moderate salinity compared with its growth in fresh water. Under warm conditions, the increased growth is accompanied by increased specific leaf area (SLA) and specific root length. Low temperature resulted in a low relative growth rate due to a low leaf area ratio and small SLA, regardless of salinity. Salinity lowered the ratio of the amounts of alternative oxidase to cytochrome c oxidase in the mitochondrial respiratory chain in leaves. Salinity enhanced the leaf respiration rate for maintenance, but under warm conditions this enhancement was compensated by a low leaf respiration rate for growth. In contrast, salinity enhanced overall leaf respiration rates at low temperature. Our results indicate that under moderate saline conditions R. stylosa leaves require warm temperatures to grow with a high rate of resource acquisition without enhancing respiratory cost.

Genome assembly and multiomic analyses reveal insights into flower and bark colors of Lagerstroemia excelsa

Lagerstroemia excelsa is a unique plant species from China, holds a significant aesthetic and economic value, and plays a crucial role in landscape architecture and horticulture. Thus far, there is little genetic and genomic information available about this species, which limits its use in development of new cultivars. In this study, a high-quality genome map of L. excelsa was obtained via whole-genome sequencing. Results showed that its genome size is about 330.4 Mb and a scaffold mapping rate is approximately 97.20%, resulting in 24 pseudochromosomes. L. excelsa might have undergone a recent whole-genome triplication event and diverged from the pomegranate about 32.3 million years ago (MYA). Subsequently, the divergence time between L. indica and L. excelsa was around 5.9 MYA. The transcriptomic and metabolomic analyses of L. excelsa and L. indica indicated that the chalcone synthase pathway may play a key role in regulating flower color differentiation between the two species. Additionally, a transcription factor LeMYB103 may be involved in regulating anthocyanin synthesis by interacting with LeMYB66, resulting in the accumulation of anthocyanins in the stem bark. This study is the first step toward genomic analysis of L. excelsa, which may provide a foundation for further molecular investigation of this species and offer valuable insights into the molecular mechanisms underlying the flower and stem bark colors in L. excelsa, two important ornamental traits in Lagerstroemia breeding.

Convergent Isobilateral Leaves Increase the Risk for Mangroves Facing Human-Induced Rapid Environmental Changes

Understanding plant adaptations in extreme environments is crucial, as these adaptations often confer advantages for survival. However, a significant gap exists regarding the genetic mechanisms underlying these adaptations and their responses to human-induced rapid environmental change (HIREC). This study addresses the question of whether genetic convergence occurs among plants with similar adaptive features, specifically focusing on isobilateral leaves in mangrove species. Here, we analyse the genetic convergence of isobilateral leaves in mangroves that have independently adapted to coastal intertidal zones. Our findings reveal that genetic convergence is evident in gene families involved in leaf abaxial and adaxial development, with strong selection pressures identified in photosynthesis and leaf polarity pathways. Despite these adaptations, mangrove species with isobilateral leaves occupy narrower ecological niches and face diminishing suitable habitat areas projected under various HIREC scenarios. These results indicate that while convergent traits enhance local adaptation, they may also increase vulnerability to ongoing environmental changes. This research provides valuable insight into the interplay between genetic adaptation and environmental resilience, underscoring the necessity for targeted biodiversity conservation strategies that safeguard specific adaptive traits amid rapid environmental shifts.

A nearly telomere-to-telomere diploid genome assembly of Firmiana kwangsiensis, a threatened species in China

Firmiana kwangsiensis is a tree species of high ornamental value. The species is critically endangered in the wild, and is listed as a first-class national protected wild plant in China, and a Plant Species with Extremely Small Populations in need of urgent protection. We have assembled a chromosome-scale, haplotype-resolved genome for F. kwangsiensis using a combination of PacBio HiFi sequencing, ONT sequencing, and Hi-C sequencing. The final assembled genome is 2.3 G in size and comprises 2n = 40 chromosomes. All chromosomal ends contain telomeric characteristic motifs (TTTAGGG), and there are only 2 gaps within the rDNA regions, both close to a T2T genome assembly. Two complete sets of haplotypes are present, Haplotype A (1169.19 Mb) and Haplotype B (1157.87 Mb), with contig N50 lengths of 58.37 Mb and 57.27 Mb, respectively. The genome contains a total of 67,527 coding genes, with 62,351 genes functionally annotated here. This is the first report of the genome of F. kwangsiensis, and lays the foundation for future conservation genomics research into this species.

MangroveDB: A Comprehensive Online Database for Mangroves Based on Multi-Omics Data

Mangroves are dominant flora of intertidal zones along tropical and subtropical coastline around the world that offer important ecological and economic value. Recently, the genomes of mangroves have been decoded, and massive omics data were generated and deposited in the public databases. Reanalysis of multi-omics data can provide new biological insights excluded in the original studies. However, the requirements for computational resource and lack of bioinformatics skill for experimental researchers limit the effective use of the original data. To fill this gap, we uniformly processed 942 transcriptome data, 386 whole-genome sequencing data, and provided 13 reference genomes and 40 reference transcriptomes for 53 mangroves. Finally, we built an interactive web-based database platform MangroveDB (https://github.com/Jasonxu0109/MangroveDB), which was designed to provide comprehensive gene expression datasets to facilitate their exploration and equipped with several online analysis tools, including principal components analysis, differential gene expression analysis, tissue-specific gene expression analysis, GO and KEGG enrichment analysis. MangroveDB not only provides query functions about genes annotation, but also supports some useful visualization functions for analysis results, such as volcano plot, heatmap, dotplot, PCA plot, bubble plot, population structure, and so on. In conclusion, MangroveDB is a valuable resource for the mangroves research community to efficiently use the massive public omics datasets.

Phylogenomics reveals a new generic delimitation of Asian Millettia, with reinstatements of Pongamia and Otosema (Fabaceae: Millettieae)

The generic relationships within the core Millettieae have not been well resolved to date. Millettia is considered one of the most problematic genera to classify within the core Millettieae, with approximately 150 species found in Asia and Africa. Recent molecular evidence has indicated that Millettia is highly polyphyletic. In this study, we utilized plastomes and single-copy nuclear genes from genome skimming sequencing to reconstruct the generic relationships and estimate the divergence time for the core Millettieae, focusing on the Asian species of Millettia. Our results revealed a broader circumscription of the core Millettieae with the inclusion of the Phaseoleae subtribe Ophrestiinae, consisting of 7 clades and 14 subclades within two main groups (canavanine and non-canavanine). The Asian species of Millettia were grouped into three distinct, well-supported subclades, leading us to propose a much narrower generic concept of Millettia, along with the reinstatements of the genera Pongamia and Otosema. Millettia s.str. was restricted to include only seven species. Pongamia was recognized as a medium-sized genus that includes Fordia and Ibatiria, with approximately 56 species. Otosema, usually treated as a subgenus or section of Millettia, was resurrected as a distinct genus with three species. Our findings also indicated that the sampled African Millettia species belong to three subclades distinct from the Asian Millettia subclades, warranting further investigation in future studies. The dating result showed that the early rapid diversification of the core Millettieae during the Mid-Miocene Climatic Optimum involved several major problematic genera, including Derris s.lato, Lonchocarpus s.lato, and Millettia s.lato, which may explain why these plant groups are difficult to classify. These results offer new insights into the phylogeny and taxonomy of the core Millettieae.

Comparative Plastome Analysis Between Endangered Mangrove Species Acanthus ebracteatus and Acanthus Relatives Provides Insights into Its Origin and Adaptive Evolution

Acanthus ebracteatus is a typical true mangrove species with great ecological and medicinal values. However, it has become endangered in China. Moreover, because of the similar morphology and distribution, it is commonly confused with the congeneric mangrove species, A. ilicifolius, which poses challenges to the protection and proper medicinal utilization of A. ebracteatus. Plastomes provide a solution for molecular identification and adaptive evolution investigation of plants. In this study, we dissected the complete plastome for A. ebracteatus and performed comparative analysis to A. ilicifolius and three non-mangrove relatives (A. montanus, A. leucostachyus and A. mollis). Both plastome sequences and structure are highly conserved between the two mangrove species, while less similar between mangrove and non-mangrove species. Phylogenetic analysis suggested that the mangrove species were divergent from the non-mangrove groups at approximately 15.15 million years ago (Mya), where early to middle Miocene global warming and high sea level might act as one of the main forces driving the mangrove lineage entering into intertidal environments. Furthermore, 12 single nucleotide polymorphisms (SNPs) and 10 insertions/deletions (indels) were detected between the plastomes of A. ebracteatus and A. ilicifolius. PCR validation further demonstrated the effectiveness of the plastid marker in distinguishing the two sibling mangrove species. Taken together, our study broadens the understanding of the origin and evolution of Acanthus mangrove plants, and provided valuable information on the correct identification and protection of endangered mangrove species A. ebracteatus.

CYP76BK1 orthologs catalyze furan and lactone ring formation in clerodane diterpenoids across the mint family

The Lamiaceae (mint family) is the largest known source of furanoclerodanes, a subset of clerodane diterpenoids with broad bioactivities including insect antifeedant properties. The Ajugoideae subfamily, in particular, accumulates significant numbers of structurally related furanoclerodanes. The biosynthetic capacity for formation of these diterpenoids is retained across most Lamiaceae subfamilies, including the early-diverging Callicarpoideae which forms a sister clade to the rest of Lamiaceae. VacCYP76BK1, a cytochrome P450 monooxygenase from Vitex agnus-castus, was previously found to catalyze the formation of the proposed precursor to furan and lactone-containing labdane diterpenoids. Through transcriptome-guided pathway exploration, we identified orthologs of VacCYP76BK1 in Ajuga reptans and Callicarpa americana. Functional characterization demonstrated that both could catalyze the oxidative cyclization of clerodane backbones to yield a furan ring. Subsequent investigation revealed a total of 10 CYP76BK1 orthologs across six Lamiaceae subfamilies. Through analysis of available chromosome-scale genomes, we identified four CYP76BK1 members as syntelogs within a conserved syntenic block across divergent subfamilies. This suggests an evolutionary lineage that predates the speciation of the Lamiaceae. Functional characterization of the CYP76BK1 orthologs affirmed conservation of function, as all catalyzed furan ring formation. Additionally, some orthologs yielded two novel lactone ring moieties. The presence of the CYP76BK1 orthologs across Lamiaceae subfamilies closely overlaps with the distribution of reported furanoclerodanes. Together, the activities and distribution of the CYP76BK1 orthologs identified here support their central role in furanoclerodane biosynthesis within the Lamiaceae family. Our findings lay the groundwork for biotechnological applications to harness the economic potential of this promising class of compounds.

Genomic convergence in terrestrial root plants through tandem duplication in response to soil microbial pressures

Despite increasing reports of convergent adaptation, evidence for genomic convergence across diverse species worldwide is lacking. Here, our study of 205 Archaeplastida genomes reveals evidence of genomic convergence through tandem duplication (TD) across different lineages of root plants despite their genomic diversity. TD-derived genes, notably prevalent in trees with developed root systems embedded in soil, are enriched in enzymatic catalysis and biotic stress responses, suggesting adaptations to environmental pressures. Correlation analyses suggest that many factors, particularly those related to soil microbial pressures, are significantly associated with TD dynamics. Conversely, flora transitioned to aquatic, parasitic, halophytic, or carnivorous lifestyles-reducing their interaction with soil microbes-exhibit a consistent decline in TD frequency. This trend is further corroborated in mangroves that independently adapted to hypersaline intertidal soils, characterized by diminished microbial activity. Our findings propose TD-driven genomic convergence as a widespread adaptation to soil microbial pressures among terrestrial root plants.

The diversity, drivers, consequences and management of plant invasions in the mangrove ecosystems

Mangrove ecosystems, which occupy intertidal environments across tropical and subtropical regions, provide crucial ecosystem services, such as protecting the coastal areas by reducing the impact of cyclones, storms, and tidal waves. Anthropogenic activities such as human settlements, deforestation, pollution, and climate change have increased the risk of biological invasions in mangrove habitats. Plant species can be introduced to mangrove habitats via anthropogenic means, such as trade and transportation, urbanisation, and agriculture, as well as through natural processes like wind, floods, cyclones, and animal-assisted seed dispersal. Additionally, some native plants can become invasive due to the changes in the mangrove ecosystem. Invasive species can significantly affect coastal ecosystems by out-competing native flora for resources, thereby altering fundamental properties, functions, and ecosystem services of the mangrove forests. The successful establishment of invasive species depends on a complex interplay of factors involving the biological attributes of the invading species and the ecological dynamics of the invaded habitat. This review focuses on exploring the mechanisms of invasion, strategies used by invasive plants, the effects of invasive plants on mangrove habitats and their possible management strategies. Based on the literature, managing invasive species is possible by biological, chemical, or physical methods. Some non-native mangrove species introduced through restoration activities can often become more intrusive than native species. Therefore, restoration activities should prioritise avoiding the use of non-native plant species.

Plant reproduction: Mangroves use 'child-care' in extreme environments

Using population genomics and molecular biological methods, a recent study reveals the molecular mechanisms of vivipary in mangroves, a term designating all woody plants of the tropical intertidal zones.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Methanolic Extract of Rhizophora mangle (Rhizophoraceae) Leaves: Phytochemical Characterization and Anthelmintic Evaluation against Schistosoma mansoni

Rhizophora mangle is commonly used in traditional medicine to treat infections, reduce inflammation, and promote healing. This study aimed to analyze the phytochemical profile of the methanolic extract of R. mangle leaves (MELRm) and evaluate its in vitro schistosomicidal activity against Schistosoma mansoni as well as its cytotoxicity. Plant material was collected in Itamaracá City, Pernambuco, Brazil. The extract was analyzed using UV/Vis spectrophotometry and high-performance liquid chromatography (HPLC). The motility, mortality, and cell viability of adult worms were assessed in a schistosomicidal assay, while cytotoxicity was evaluated through a colorimetric assay with MTT on RAW 264.7 cells. The primary compounds identified in MELRm were phenolic compounds. In the schistosomicidal assay, all concentrations of MELRs induced changes in the motility of adult worms. At a concentration of 400 μg/mL, MELRs resulted in 56.25% mortality after 72 h of incubation. After 120 h, mortality rates of 75%, 62.5%, and 50% were observed at MELRm concentrations of 400, 200, and 100 μg/mL, respectively. No eggs were detected at any MELRm concentration. MELRs did not show cytotoxicity towards RAW 264.7 cells at the concentrations tested. These results indicate that MELRs demonstrate schistosomicidal activity in vitro, suggesting they are promising candidates for in vivo studies.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes

Spartina alterniflora is an exo-recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high-quality, chromosome-level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high-throughput chromosome conformation capture (Hi-C) technology and Illumina-based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole-genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species-specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA-seq analysis identified several ion transporter genes including the high-affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up-regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ-ECS, GSH-S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte.

CYP76BK1 orthologs catalyze furan and lactone ring formation in clerodane diterpenoids across the mint family

The Lamiaceae (mint family) is the largest known source of furanoclerodanes, a subset of clerodane diterpenoids with broad bioactivities including insect antifeedant properties. The Ajugoideae subfamily, in particular, accumulates significant numbers of structurally related furanoclerodanes. The biosynthetic capacity for formation of these diterpenoids is retained across most Lamiaceae subfamilies, including the early-diverging Callicarpoideae which forms a sister clade to the rest of Lamiaceae. VacCYP76BK1, a cytochrome P450 monooxygenase from Vitex agnus-castus, was previously found to catalyze the formation of the proposed precursor to furan and lactone-containing labdane diterpenoids. Through transcriptome-guided pathway exploration, we identified orthologs of VacCYP76BK1 in Ajuga reptans and Callicarpa americana. Functional characterization demonstrated that both could catalyze the oxidative cyclization of clerodane backbones to yield a furan ring. Subsequent investigation revealed a total of ten CYP76BK1 orthologs across six Lamiaceae subfamilies. Through analysis of available chromosome-scale genomes, we identified four CYP76BK1 members as syntelogs within a conserved syntenic block across divergent subfamilies. This suggests an evolutionary lineage that predates the speciation of the Lamiaceae. Functional characterization of the CYP76BK1 orthologs affirmed conservation of function, as all catalyzed furan ring formation. Additionally, some orthologs yielded two novel lactone ring moieties. The presence of the CYP76BK1 orthologs across Lamiaceae subfamilies closely overlaps with the distribution of reported furanoclerodanes. Together, the activities and distribution of the CYP76BK1 orthologs identified here support their central role in furanoclerodane biosynthesis within the Lamiaceae family. Our findings lay the groundwork for biotechnological applications to harness the economic potential of this promising class of compounds.

Global analysis of key post-transcriptional regulation in early leaf development of Limonium bicolor identifies a long non-coding RNA that promotes salt gland development and salt resistance

Limonium bicolor, known horticulturally as sea lavender, is a typical recretohalophyte with salt glands in its leaf epidermis that secrete excess Na+ out of the plant. Although many genes have been proposed to contribute to salt gland initiation and development, a detailed analysis of alternative splicing, alternative polyadenylation patterns, and long non-coding RNAs (lncRNAs) has been lacking. Here, we applied single-molecule long-read mRNA isoform sequencing (Iso-seq) to explore the complexity of the L. bicolor transcriptome in leaves during salt gland initiation (stage A) and salt gland differentiation (stage B) based on the reference genome. We identified alternative splicing events and the use of alternative poly(A) sites unique to stage A or stage B, leading to the hypothesis that they might contribute to the differentiation of salt glands. Based on the Iso-seq data and RNA in situ hybridization of candidate genes, we selected the lncRNA Btranscript_153392 for gene editing and virus-induced gene silencing to dissect its function. In the absence of this transcript, we observed fewer salt glands on the leaf epidermis, leading to diminished salt secretion and salt tolerance. Our data provide transcriptome resources for unraveling the mechanisms behind salt gland development and furthering crop transformation efforts towards enhanced survivability in saline soils.

Datasets of fungal diversity and pseudo-chromosomal genomes of mangrove rhizosphere soil in China

With climate change and anthropic influence on the coastal ecosystems, mangrove ecosystems are disappearing at an alarming rate. Accordingly, it becomes important to track, study, record and store the mangrove microbial community considering their ecological importance and potential for biotechnological applications. Here, we provide information on mangrove fungal community composition and diversity in mangrove ecosystems with different plant species and from various locations differing in relation to anthropic influences. We describe twelve newly assembled genomes, including four chromosomal-level genomes of fungal isolates from the mangrove ecosystems coupled with functional annotations. We envisage that these data will be of value for future studies including comparative genome analysis and large-scale temporal and/or spatial research to elucidate the potential mechanisms by which mangrove fungal communities assemble and evolve. We further anticipate that the genomes represent valuable resources for bioprospecting related to industrial or clinical uses.

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.

Asteraceae genome database: a comprehensive platform for Asteraceae genomics

Asteraceae, the largest family of angiosperms, has attracted widespread attention for its exceptional medicinal, horticultural, and ornamental value. However, researches on Asteraceae plants face challenges due to their intricate genetic background. With the continuous advancement of sequencing technology, a vast number of genomes and genetic resources from Asteraceae species have been accumulated. This has spurred a demand for comprehensive genomic analysis within this diverse plant group. To meet this need, we developed the Asteraceae Genomics Database (AGD; http://cbcb.cdutcm.edu.cn/AGD/). The AGD serves as a centralized and systematic resource, empowering researchers in various fields such as gene annotation, gene family analysis, evolutionary biology, and genetic breeding. AGD not only encompasses high-quality genomic sequences, and organelle genome data, but also provides a wide range of analytical tools, including BLAST, JBrowse, SSR Finder, HmmSearch, Heatmap, Primer3, PlantiSMASH, and CRISPRCasFinder. These tools enable users to conveniently query, analyze, and compare genomic information across various Asteraceae species. The establishment of AGD holds great significance in advancing Asteraceae genomics, promoting genetic breeding, and safeguarding biodiversity by providing researchers with a comprehensive and user-friendly genomics resource platform.

Utilities of environmental radioactivity tracers in assessing sequestration potential of carbon in the coastal wetland ecosystems

Demand for accurate estimation of coastal blue carbon sequestration rates in a regular interval has recently surged due to the increasing awareness of nature-based climate solutions to alleviate adverse impacts stemming from the recent global warming. The robust estimation method is, however, far from well-established. The international community requires, moreover, to quantify its effect of "management." This article tries to provide the environmental isotope community with basic biophysical features of coastal blue carbon ecosystems to identify a suitable set of environmental isotopes for promoting coastal ocean-based climate solutions. This article reviews (i) the primary biophysical characteristics of coastal blue carbon ecosystems and hydrology, (ii) their consequential impact on the accumulation and preservation of organic carbon (OC) in the sediment column, (iii) suitable environmental isotopes to quantifying the sedimentary organic carbon accumulation, outwelling of the carbon-containing byproducts of decomposition of biogenic organic matter and acid neutralizing alkalinity produced in situ sediment to the offshore. Above-ground biomass is not cumulative over the years except for mangrove forests within coastal blue carbon systems. Non-gaseous carbon sequestration and loss occur mainly as a form of sediment organic carbon (SOC) and dissolved carbon in an intertidal and subtidal bottom sediment body in a slow, patchy, and dispersive way, on which this article focuses. Investigating environmental radionuclides is probably the most cost-effective effort to contribute to defining the offshore spatial extent of coastal blue carbon systems except for seagrass beds (e.g., Ra isotopes), to quantify millimeter per year scale carbon accretion and loss within the systems (e.g., 7Be, 210Pb) and a liter per meter of coastline per a day scale water movement from the systems (Ra isotopes). A millimeter-scale spatial and an annual (or less) time-scale resolution offered by the use of environmental isotopes would equip us with a novel tool to enhance the carbon storage capacity of the coastal blue carbon system.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Origin and evolution of a new tetraploid mangrove species in an intertidal zone

Polyploidy is a major factor in the evolution of plants, yet we know little about the origin and evolution of polyploidy in intertidal species. This study aimed to identify the evolutionary transitions in three true-mangrove species of the genus Acanthus distributed in the Indo-West Pacific region. For this purpose, we took an integrative approach that combined data on morphology, cytology, climatic niche, phylogeny, and biogeography of 493 samples from 42 geographic sites. Our results show that the Acanthus ilicifolius lineage distributed east of the Thai-Malay Peninsula possesses a tetraploid karyotype, which is morphologically distinct from that of the lineage on the west side. The haplotype networks and phylogenetic trees for the chloroplast genome and eight nuclear genes reveal that the tetraploid species has two sub-genomes, one each from A. ilicifolius and A . ebracteatus, the paternal and maternal parents, respectively. Population structure analysis also supports the hybrid speciation history of the new tetraploid species. The two sub-genomes of the tetraploid species diverged from their diploid progenitors during the Pleistocene. Environmental niche models revealed that the tetraploid species not only occupied the near-entire niche space of the diploids, but also expanded into novel environments. Our findings suggest that A. ilicifolius species distributed on the east side of the Thai-Malay Peninsula should be regarded as a new species, A. tetraploideus, which originated from hybridization between A. ilicifolius and A. ebracteatus, followed by chromosome doubling. This is the first report of a true-mangrove allopolyploid species that can reproduce sexually and clonally reproduction, which explains the long-term adaptive potential of the species.

Comparative analysis of lipid and flavonoid biosynthesis between Pongamia and soybean seeds: genomic, transcriptional, and metabolic perspectives

BACKGROUND

Soybean (Glycine max) is a vital oil-producing crop. Augmenting oleic acid (OA) levels in soybean oil enhances its oxidative stability and health benefits, representing a key objective in soybean breeding. Pongamia (Pongamia pinnata), known for its abundant oil, OA, and flavonoid in the seeds, holds promise as a biofuel and medicinal plant. A comparative analysis of the lipid and flavonoid biosynthesis pathways in Pongamia and soybean seeds would facilitate the assessment of the potential value of Pongamia seeds and advance the genetic improvements of seed traits in both species.

RESULTS

The study employed multi-omics analysis to systematically compare differences in metabolite accumulation and associated biosynthetic genes between Pongamia seeds and soybean seeds at the transcriptional, metabolic, and genomic levels. The results revealed that OA is the predominant free fatty acid in Pongamia seeds, being 8.3 times more abundant than in soybean seeds. Lipidomics unveiled a notably higher accumulation of triacylglycerols (TAGs) in Pongamia seeds compared to soybean seeds, with 23 TAG species containing OA. Subsequently, we identified orthologous groups (OGs) involved in lipid biosynthesis across 25 gene families in the genomes of Pongamia and soybean, and compared the expression levels of these OGs in the seeds of the two species. Among the OGs with expression levels in Pongamia seeds more than twice as high as in soybean seeds, we identified one fatty acyl-ACP thioesterase A (FATA) and two stearoyl-ACP desaturases (SADs), responsible for OA biosynthesis, along with two phospholipid:diacylglycerol acyltransferases (PDATs) and three acyl-CoA:diacylglycerol acyltransferases (DGATs), responsible for TAG biosynthesis. Furthermore, we observed a significantly higher content of the flavonoid formononetin in Pongamia seeds compared to soybean seeds, by over 2000-fold. This difference may be attributed to the tandem duplication expansions of 2,7,4'-trihydroxyisoflavanone 4'-O-methyltransferases (HI4'OMTs) in the Pongamia genome, which are responsible for the final step of formononetin biosynthesis, combined with their high expression levels in Pongamia seeds.

CONCLUSIONS

This study extends beyond observations made in single-species research by offering novel insights into the molecular basis of differences in lipid and flavonoid biosynthetic pathways between Pongamia and soybean, from a cross-species comparative perspective.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Mitochondrial metagenomics reveal the independent colonization of the world's coasts by intertidal oribatid mites (Acari, Oribatida, Ameronothroidea)

Oribatid mites are an ancient group that already roamed terrestrial ecosystems in the early and middle Devonian. The superfamily of Ameronothroidea, a supposedly monophyletic lineage, represents the only group of oribatid mites that has successfully invaded the marine coastal environment. By using mitogenome data and nucleic ribosomal RNA genes (18S, 5.8S, 28S), we show that Ameronothroidea are a paraphyletic assemblage and that the land-to-sea transition happened three times independently. Common ancestors of the tropical Fortuyniidae and Selenoribatidae were the first to colonize the coasts and molecular calibration of our phylogeny dates this event to a period in the Triassic and Jurassic era (225-146 mya), whereas present-day distribution indicates that this event might have happened early in this period during the Triassic, when the supercontinent Pangaea still existed. The cold temperate northern hemispheric Ameronothridae colonized the marine littoral later in the late Jurassic-Early Cretaceous and had an ancient distribution on Laurasian coasts. The third and final land-to-sea transition happened in the same geological period, but approx. 30 my later when ancestors of Podacaridae invaded coastal marine environments of the Gondwanan landmasses.

Chromosome-level genome assembly of the threatened resource plant Cinnamomum chago

Cinnamomum chago is a tree species endemic to Yunnan province, China, with potential economic value, phylogenetic importance, and conservation priority. We assembled the genome of C. chago using multiple sequencing technologies, resulting in a high-quality, chromosomal-level genome with annotation information. The assembled genome size is approximately 1.06 Gb, with a contig N50 length of 92.10 Mb. About 99.92% of the assembled sequences could be anchored to 12 pseudo-chromosomes, with only one gap, and 63.73% of the assembled genome consists of repeat sequences. In total, 30,497 genes were recognized according to annotation, including 28,681 protein-coding genes. This high-quality chromosome-level assembly and annotation of C. chago will assist us in the conservation and utilization of this valuable resource, while also providing crucial data for studying the evolutionary relationships within the Cinnamomum genus, offering opportunities for further research and exploration of its diverse applications.

De novo assembly and analysis of Sonneratia ovata genome and population analysis

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

Stomatal dynamics are regulated by leaf hydraulic traits and guard cell anatomy in nine true mangrove species

Stomatal regulation is critical for mangroves to survive in the hyper-saline intertidal zone where water stress is severe and water availability is highly fluctuant. However, very little is known about the stomatal sensitivity to vapour pressure deficit (VPD) in mangroves, and its co-ordination with stomatal morphology and leaf hydraulic traits. We measured the stomatal response to a step increase in VPD in situ, stomatal anatomy, leaf hydraulic vulnerability and pressure-volume traits in nine true mangrove species of five families and collected the data of genome size. We aimed to answer two questions: (1) Does stomatal morphology influence stomatal dynamics in response to a high VPD in mangroves? with a consideration of possible influence of genome size on stomatal morphology; and (2) do leaf hydraulic traits influence stomatal sensitivity to VPD in mangroves? We found that the stomata of mangrove plants were highly sensitive to a step rise in VPD and the stomatal responses were directly affected by stomatal anatomy and hydraulic traits. Smaller, denser stomata was correlated with faster stomatal closure at high VPD across the species of Rhizophoraceae, and stomata size negatively and vein density positively correlated with genome size. Less negative leaf osmotic pressure at the full turgor (πo) was related to higher operating steady-state stomatal conductance (gs); and a higher leaf capacitance (Cleaf) and more embolism resistant leaf xylem were associated with slower stomatal responses to an increase in VPD. In addition, stomatal responsiveness to VPD was indirectly affected by leaf morphological traits, which were affected by site salinity and consequently leaf water status. Our results demonstrate that mangroves display a unique relationship between genome size, stomatal size and vein packing, and that stomatal responsiveness to VPD is regulated by leaf hydraulic traits and stomatal morphology. Our work provides a quantitative framework to better understand of stomatal regulation in mangroves in an environment with high salinity and dynamic water availability.

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.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Genome-wide identification of Avicennia marina aquaporins reveals their role in adaptation to intertidal habitats and their relevance to salt secretion and vivipary

Aquaporins (AQPs) regulate the transport of water and other substrates, aiding plants in adapting to stressful environments. However, the knowledge of AQPs in salt-secreting and viviparous Avicennia marina is limited. In this study, 46 AmAQPs were identified in A. marina genome, and their subcellular localisation and function in transporting H2 O2 and boron were assessed through bioinformatics analysis and yeast transformation. Through analysing their expression patterns via RNAseq and real-time quantitative polymerase chain reaction, we found that most AmAQPs were downregulated in response to salt and tidal flooding. AmPIP (1;1, 1;7, 2;8, 2;9) and AmTIP (1;5, 1;6) as salt-tolerant candidate genes may contribute to salt secretion together with Na+ /H+ antiporters. AmPIP2;1 and AmTIP1;5 were upregulated during tidal flooding and may be regulated by anaerobic-responsive element and ethylene-responsive element cis-elements, aiding in adaptation to tidal inundation. Additionally, we found that the loss of the seed desiccation and dormancy-related TIP3 gene, and the loss of the seed dormancy regulator DOG1 gene, or DOG1 protein lack heme-binding capacity, may be genetic factors contributing to vivipary. Our findings shed light on the role of AQPs in A. marina adaptation to intertidal environments and their relevance to salt secretion and vivipary.

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.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Diazotrophic nitrogen fixation through aerial roots occurs in Avicennia marina: implications for adaptation of mangrove plant growth to low-nitrogen tidal flats

Nitrogen limitation of primary production is common in coastal ecosystems. Mangrove trees maintain high levels of nitrogen fixation around their roots. The interior aerial space of mangrove roots, in which atmospheric gas is supplied through lenticels, could be efficient sites for nitrogen fixation. We measured tidal variations of partial pressure of N2 in root aerenchyma and conducted field experiments using 15 N2 as a tracer to track N2 movement through aerial roots of Avicennia marina. We used the acetylene reduction assay to identify the root parts harboring diazotrophs. The nitrogenase activity and estimated nitrogen fixation through aerenchyma were higher in pneumatophores and absorbing roots than in cable roots. Positive correlations between root nitrogen contents and turnover rates of root nitrogen derived from N2 through aerenchyma suggested that the internal supply of N2 to diazotrophs could be the main source for nitrogen assimilation by A. marina roots. Our results confirmed that N2 is supplied to diazotrophs through aerial roots and that nitrogen fixation occurs in A. marina roots. The aerial root structures, which occur across families of mangrove plants, could be an adaptation to survival in not only low-oxygen environments but also tidal flats with little plant-available nitrogen.

The transmembrane protein LbRSG from the recretohalophyte Limonium bicolor enhances salt gland development and salt tolerance

Salt glands are the unique epidermal structures present in recretohalophytes, plants that actively excrete excess Na+ by salt secretory structures to avoid salt damage. Here, we describe a transmembrane protein that localizes to the plasma membrane of the recretohalophyte Limonium bicolor. As virus-induced gene silencing of the corresponding gene LbRSG in L. bicolor decreased the number of salt glands, we named the gene Reduced Salt Gland. We detected LbRSG transcripts in salt glands by in situ hybridization and transient transformation. Overexpression and silencing of LbRSG in L. bicolor pointed to a positive role in salt gland development and salt secretion by interacting with Lb3G16832. Heterologous LbRSG expression in Arabidopsis enhanced salt tolerance during germination and the seedling stage by alleviating NaCl-induced ion stress and osmotic stress after replacing or deleting the (highly) negatively charged region of extramembranous loop. After screened by immunoprecipitation-mass spectrometry and verified using yeast two-hybrid, PGK1 and BGLU18 were proposed to interact with LbRSG to strengthen salt tolerance. Therefore, we identified (highly) negatively charged regions in the extramembrane loop that may play an essential role in salt tolerance, offering hints about LbRSG function and its potential to confer salt resistance.

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.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata)

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

Invasive Trends of Spartina alterniflora in the Southeastern Coast of China and Potential Distributional Impacts on Mangrove Forests

Mangrove forests are one of the most productive and seriously threatened ecosystems in the world. The widespread invasion of Spartina alterniflora has seriously imperiled the security of mangroves as well as coastal mudflat ecosystems. Based on a model evaluation index, we selected RF, GBM, and GLM as a predictive model for building a high-precision ensemble model. We used the species occurrence records combined with bioclimate, sea-land topography, and marine environmental factors to predict the potentially suitable habitats of mangrove forests and the potentially suitable invasive habitats of S. alterniflora in the southeastern coast of China. We then applied the invasion risk index (IRI) to assess the risk that S. alterniflora would invade mangrove forests. The results show that the suitable habitats for mangrove forests are mainly distributed along the coastal provinces of Guangdong, Hainan, and the eastern coast of Guangxi. The suitable invasive habitats for S. alterniflora are mainly distributed along the coast of Zhejiang, Fujian, and relatively less in the southern provinces. The high-risk areas for S. alterniflora invasion of mangrove forests are concentrated in Zhejiang and Fujian. Bioclimate variables are the most important variables affecting the survival and distribution of mangrove forests and S. alterniflora. Among them, temperature is the most important environmental variable determining the large-scale distribution of mangrove forests. Meanwhile, S. alterniflora is more sensitive to precipitation than temperature. Our results can provide scientific insights and references for mangrove forest conservation and control of S. alterniflora.

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.

Genetic structure of an endangered species Ormosia henryi in southern China, and implications for conservation

BACKGROUND

The evergreen broadleaved forest (EBLF) is an iconic vegetation type of East Asia, and it contributes fundamentally to biodiversity-based ecosystem functioning and services. However, the native habitat of EBLFs keeps on decreasing due to anthropogenic activities. Ormosia henryi is a valuable rare woody species in EBLFs that is particularly sensitive to habitat loss. In this study, ten natural populations of O. henryi in southern China were sampled, and then genotyping by sequencing (GBS) was applied to elucidate the standing genetic variation and population structure of this endangered species.

RESULTS

In ten O. henryi populations, 64,158 high-quality SNPs were generated by GBS. Based on these markers, a relatively low level of genetic diversity was found with the expected heterozygosity (He) ranging from 0.2371 to 0.2901. Pairwise F_ST between populations varied from 0.0213 to 0.1652, indicating a moderate level of genetic differentiation. However, contemporary gene flow between populations were rare. Assignment test and principal component analysis (PCA) both supported that O. henryi populations in southern China could be divided into four genetic groups, and prominent genetic admixture was found in those populations located in southern Jiangxi Province. Mantel tests and multiple matrix regression with randomization (MMRR) analyses suggested that isolation by distance (IBD) could be the possible reason for describing the current population genetic structure. In addition, the effective population size (Ne) of O. henryi was extremely small, and showed a continuous declining trend since the Last Glacial Period.

CONCLUSIONS

Our results indicate that the endangered status of O. henryi is seriously underestimated. Artificial conservation measures should be applied as soon as possible to prevent O. henryi from the fate of extinction. Further studies are needed to elucidate the mechanism that leading to the continuous loss of genetic diversity in O. henryi and help to develop a better conservation strategy.

Chromosomal-Scale Genome Assemblies of Two Coastal Plant Species, Scaevola taccada and S. hainanensis-Insight into Adaptation Outside of the Common Range

While most of the species in Goodeniaceae family, excluding the Scaevola genus, are endemic to Australasia, S. taccada and S. hainanensis have expanded their distribution range to the tropical coastlines of the Atlantic and Indian Oceans. S. taccada appears to be highly adapted to coastal sandy lands and cliffs, and it has become invasive in places. S. hainanensis is found mainly in salt marshes near mangrove forests, and is at risk of extinction. These two species provide a good system to investigate adaptive evolution outside the common distribution range of this taxonomic group. Here, we report their chromosomal-scale genome assemblies with the objective of probing their genomic mechanisms related to divergent adaptation after leaving Australasia. The scaffolds were assembled into eight chromosome-scale pseudomolecules, which covered 90.12% and 89.46% of the whole genome assembly for S. taccada and S. hainanensis, respectively. Interestingly, unlike many mangroves, neither species has undergone whole-genome duplication. We show that private genes, specifically copy-number expanded genes are essential for stress response, photosynthesis, and carbon fixation. The gene families that are expanded in S. hainanensis and contracted in S. taccada might have facilitated adaptation to high salinity in S. hainanensis. Moreover, the genes under positive selection in S. hainanensis have contributed to its response to stress and its tolerance of flooding and anoxic environments. In contrast, compared with S. hainanensis, the more drastic copy number expansion of FAR1 genes in S. taccada might have facilitated its adaptation to the stronger light radiation present in sandy coastal lands. In conclusion, our study of the chromosomal-scale genomes of S. taccada and S. hainanensis provides novel insights into their genomic evolution after leaving Australasia.

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.

Diverse mangroves deviate from other angiosperms in their genome size, leaf cell size and cell packing density relationships

BACKGROUND AND AIMS

While genome size limits the minimum sizes and maximum numbers of cells that can be packed into a given leaf volume, mature cell sizes can be substantially larger than their meristematic precursors and vary in response to abiotic conditions. Mangroves are iconic examples of how abiotic conditions can influence the evolution of plant phenotypes.

METHODS

Here, we examined the coordination between genome size, leaf cell sizes, cell packing densities and leaf size in 13 mangrove species across four sites in China. Four of these species occurred at more than one site, allowing us to test the effect of climate on leaf anatomy.

RESULTS

We found that genome sizes of mangroves were very small compared to other angiosperms, but, like other angiosperms, mangrove cells were always larger than the minimum size defined by genome size. Increasing mean annual temperature of a growth site led to higher packing densities of veins (Dv) and stomata (Ds) and smaller epidermal cells but had no effect on stomatal size. In contrast to other angiosperms, mangroves exhibited (1) a negative relationship between guard cell size and genome size; (2) epidermal cells that were smaller than stomata; and (3) coordination between Dv and Ds that was not mediated by epidermal cell size. Furthermore, mangrove epidermal cell sizes and packing densities covaried with leaf size.

CONCLUSIONS

While mangroves exhibited coordination between veins and stomata and attained a maximum theoretical stomatal conductance similar to that of other angiosperms, the tissue-level tradeoffs underlying these similar relationships across species and environments were markedly different, perhaps indicative of the unique structural and physiological adaptations of mangroves to their stressful environments.

Meliaceae genomes provide insights into wood development and limonoids biosynthesis

Meliaceae is a useful plant family owing to its high-quality timber and its many limonoids that have pharmacological and biological activities. Although some genomes of Meliaceae species have been reported, many questions regarding their unique family features, namely wood quality and natural products, have not been answered. In this study, we provide the whole-genome sequence of Melia azedarach comprising 237.16 Mb with a contig N50 of 8.07 Mb, and an improved genome sequence of Azadirachta indica comprising 223.66 Mb with a contig N50 of 8.91 Mb. Moreover, genome skimming data, transcriptomes and other published genomes were comprehensively analysed to determine the genes and proteins that produce superior wood and valuable limonoids. Phylogenetic analysis of chloroplast genomes, single-copy gene families and single-nucleotide polymorphisms revealed that Meliaceae should be classified into two subfamilies: Cedreloideae and Melioideae. Although the Meliaceae species did not undergo additional whole-genome duplication events, the secondary wall biosynthetic genes of the woody Cedreloideae species, Toona sinensis, expanded significantly compared to those of A. indica and M. azedarach, especially in downstream transcription factors and cellulose/hemicellulose biosynthesis-related genes. Moreover, expanded special oxidosqualene cyclase catalogues can help diversify Sapindales skeletons, and the clustered genes that regulate terpene chain elongation, cyclization and modification would support their roles in limonoid biosynthesis. The expanded clans of terpene synthase, O-methyltransferase and cytochrome P450, which are mainly derived from tandem duplication, are responsible for the different limonoid classes among the species. These results are beneficial for further investigations of wood development and limonoid biosynthesis.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.