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

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.

Genetic diversity and population structure of natural provenances of Sonneratia caseolaris in Vietnam

BACKGROUND

Sommeratia caseolaris is considered the most important mangrove species for reforestation and conservation programs. Therefore, the knowledge of genetic diversity and the population structure of the species has important implications both for the conservation of existing genetic resources and development programs. In the present study, the genetic diversity and structure population of eight populations of S. caseolaris from the Northern to the Southern Coast of Vietnam were determined using nine ISSR molecular markers.

RESULTS

Eight populations of the mangrove species Sonneratia caseolaris were sampled across the natural range in Vietnam to evaluate the genetic diversity of the species. Nine ISSR markers were used to analyse 30 individuals from each population. There were moderate to high levels of genetic diversity (I = 0.447; h = 0.300). PCoA analysis gave very similar results to UPGMA dendrogram construction with the eight populations clustered into three genetic groups which mostly aligned with geographical distances among them. AMOVA analysis results indicated that most (81 %) of the genetic variation was within populations.

CONCLUSION

The current study also indicates the high level of genetic variation existing among and within the natural population of S. caseolaris in Vietnam. These results open new perspectives towards the conservation of the species' genetic resources and their future use in conservation and reforestation programs.

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.

Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree

BACKGROUND AND AIMS

In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado.

METHODS

The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated.

KEY RESULTS

High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction.

CONCLUSIONS

The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Manifestation of Triploid Heterosis in the Root System after Crossing Diploid and Autotetraploid Energy Willow Plants

Successful use of woody species in reducing climatic and environmental risks of energy shortage and spreading pollution requires deeper understanding of the physiological functions controlling biomass productivity and phytoremediation efficiency. Targets in the breeding of energy willow include the size and the functionality of the root system. For the combination of polyploidy and heterosis, we have generated triploid hybrids (THs) of energy willow by crossing autotetraploid willow plants with leading cultivars (Tordis and Inger). These novel Salix genotypes (TH3/12, TH17/17, TH21/2) have provided a unique experimental material for characterization of Mid-Parent Heterosis (MPH) in various root traits. Using a root phenotyping platform, we detected heterosis (TH3/12: MPH 43.99%; TH21/2: MPH 26.93%) in the size of the root system in soil. Triploid heterosis was also recorded in the fresh root weights, but it was less pronounced (MPH%: 9.63-19.31). In agreement with root growth characteristics in soil, the TH3/12 hybrids showed considerable heterosis (MPH: 70.08%) under in vitro conditions. Confocal microscopy-based imaging and quantitative analysis of root parenchyma cells at the division-elongation transition zone showed increased average cell diameter as a sign of cellular heterosis in plants from TH17/17 and TH21/2 triploid lines. Analysis of the hormonal background revealed that the auxin level was seven times higher than the total cytokinin contents in root tips of parental Tordis plants. In triploid hybrids, the auxin-cytokinin ratios were considerably reduced in TH3/12 and TH17/17 roots. In particular, the contents of cytokinin precursor, such as isopentenyl adenosine monophosphate, were elevated in all three triploid hybrids. Heterosis was also recorded in the amounts of active gibberellin precursor, GA19, in roots of TH3/12 plants. The presented experimental findings highlight the physiological basics of triploid heterosis in energy willow roots.

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.

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.

Comparative transcriptome analysis unveiling reactive oxygen species scavenging system of Sonneratia caseolaris under salinity stress

Many mangrove forests have undergone major changes as a result of human activity and global climate change. Sonneratia caseolaris is a common tree located in inner mangroves, and its range extends inland along tidal creeks, as far as the influence of salinity extends. This study investigated the physiological and molecular response mechanisms of S. caseolaris by analyzing its antioxidant defense capacity, including its differentially expressed genes (DEGs) under similar salt stress conditions. Salt treatment significantly affected the osmoprotectants and lipid peroxidation in S. caseolaris seedlings, which increased proline (Pro) content by 31.01-54.90% during all sample periods and decreased malonaldehyde (MDA) content by 12.81 and 18.17% at 25 and 40 days under 3.0% NaCl treatment. Antioxidant enzyme activities increased significantly following 3.0% NaCl treatment. Transcriptome analysis following De novo assembly showed 26,498 matched unigenes. The results showed that 1,263 DEGs responded to transcription factors (TFs) and plant phytohormones and mediated oxidoreductase activity to scavenge reactive oxygen species (ROS) in the control vs. 3.0% NaCl comparison. In addition, the transcription levels of genes associated with auxin and ethylene signal transduction also changed. Under salt stress, ROS scavenging genes (POD, CAT, and APX) and part of AP2, MYB, NAC, C2C2, bHLH, and WRKY TFs were upregulated. This study identified important pathways and candidate genes involved in S. caseolaris salinity tolerance and provided suggestions for further research into the mechanisms of salt tolerance in S. caseolaris.

Integrated Full-Length Transcriptome and MicroRNA Sequencing Approaches Provide Insights Into Salt Tolerance in Mangrove (Sonneratia apetala Buch.-Ham.)

MicroRNAs (miRNAs) are small RNA molecules that serve as key players in plant stress responses. Although stress-regulated miRNAs have been explored in various plants, they are not well studied in mangroves. Herein, we combined PacBio isoform sequencing (Iso-Seq) with BGISEQ short-read RNA-seq to probe the role of miRNAs in the salt stress response of the mangrove plant, Sonneratia apetala Buch.-Ham. A total of 1,702,463 circular consensus sequencing reads were generated that produced 295,501 nonredundant full-length transcripts from the leaves of a 1-year-old S. apetala. After sequencing nine small RNA libraries constructed from control and 1- and 28-day 300 mM NaCl treatments, we identified 143 miRNAs (114 known and 29 novel) from a total of >261 million short reads. With the criteria of |log2FC| ≥ 1 and q-value < 0.05, 42 and 70 miRNAs were differentially accumulated after 1- and 28-day salt treatments, respectively. These differential accumulated miRNAs potentially targeted salt-responsive genes encoding transcription factors, ion homeostasis, osmotic protection, and detoxificant-related proteins, reminiscent of their responsibility for salinity adaptation in S. apetala. Particularly, 62 miRNAs were Sonneratia specific under salt stress, of which 34 were co-expressed with their 131 predicted targets, thus producing 140 miRNA-target interactions. Of these, 82 miRNA-target pairs exhibited negative correlations. Eighteen miRNA targets were categorized for the 'environmental information processing' during KEGG analysis and were related to plant hormone signal transduction (ko04075), MAPK signaling pathway-plant (ko04016), and ABC transporters (ko02010). These results underscored miRNAs as possible contributors to mangrove success in severe environments and offer insights into an miRNA-mediated regulatory mechanism of salt response in S. apetala.

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.

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.

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.

Depth-dependent variability of biological nitrogen fixation and diazotrophic communities in mangrove sediments

BACKGROUND

Nitrogen-fixing prokaryotes (diazotrophs) contribute substantially to nitrogen input in mangrove sediments, and their structure and nitrogen fixation rate (NFR) are significantly controlled by environmental conditions. Despite the well-known studies on diazotrophs in surficial sediments, the diversity, structure, and ecological functions of diazotrophic communities along environmental gradients of mangrove sediment across different depths are largely unknown. Here, we investigated how biological nitrogen fixation varied with the depth of mangrove sediments from the perspectives of both NFR and diazotrophic communities.

RESULTS

Through acetylene reduction assay, nifH gene amplicon and metagenomic sequencing, we found that the NFR increased but the diversity of diazotrophic communities decreased with the depth of mangrove sediments. The structure of diazotrophic communities at different depths was largely driven by salinity and exhibited a clear divergence at the partitioning depth of 50 cm. Among diazotrophic genera correlated with NFR, Agrobacterium and Azotobacter were specifically enriched at 50-100 cm sediments, while Anaeromyxobacter, Rubrivivax, Methylocystis, Dickeya, and Methylomonas were more abundant at 0-50 cm. Consistent with the higher NFR, metagenomic analysis demonstrated the elevated abundance of nitrogen fixation genes (nifH/D/K) in deep sediments, where nitrification genes (amoA/B/C) and denitrification genes (nirK and norB) became less abundant. Three metagenome-assembled genomes (MAGs) of diazotrophs from deep mangrove sediments indicated their facultatively anaerobic and mixotrophic lifestyles as they contained genes for low-oxygen-dependent metabolism, hydrogenotrophic respiration, carbon fixation, and pyruvate fermentation.

CONCLUSIONS

This study demonstrates the depth-dependent variability of biological nitrogen fixation in terms of NFR and diazotrophic communities, which to a certain extent relieves the degree of nitrogen limitation in deep mangrove sediments. Video Abstract.

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.

Weak Effect of Gypsy Retrotransposon Bursts on Sonneratia alba Salt Stress Gene Expression

Transposable elements (TEs) are an important source of genetic diversity and can be co-opted for the regulation of host genes. However, to what extent the pervasive TE colonization of plant genomes has contributed to stress adaptation remains controversial. Plants inhabiting harsh environments in nature provide a unique opportunity to answer this question. We compared TE compositions and their evolutionary dynamics in the genomes of two mangrove species: the pioneer Sonneratia alba and its less salt-tolerant relative S. caseolaris. Age distribution, strength of purifying selection and the removal rate of LTR (long terminal repeat) retrotransposons were estimated. Phylogenetic analysis of LTR retrotransposons and their distribution in the genome of S. alba were surveyed. Small RNA sequencing and whole-genome bisulfite sequencing was conducted using leaves of S. alba. Expression pattern of LTR retrotransposons and their nearby genes were examined using RNA-seq data of S. alba under different salt treatments. S. alba possesses more TEs than S. caseolaris. Particularly, many more young Gypsy LTR retrotransposons have accumulated in S. alba than in S. caseolaris despite an increase in purifying selection against TE insertions. The top two most abundant Gypsy families in S. alba preferentially insert in gene-poor regions. They are under relaxed epigenetic repression, probably due to the presence of CHROMO domains in their 3'-ends. Although a considerable number of TEs in S. alba showed differential expression under salt stress, only four copies were significantly correlated with their nearby genes in expression levels. One such TE-gene pair involves Abscisic acid 8'-hydroxylase 3 functioning in abscisic acid catabolism. This study sheds light on the evolutionary dynamics and potential function of TEs in an extremophile. Our results suggest that the conclusion on co-option of TEs should be cautious even though activation of TEs by stress might be prevalent.

Genomic Convergence in the Adaptation to Extreme Environments

Convergent evolution is especially common in plants that have independently adapted to the same extreme environments (i.e., extremophile plants). The recent burst of omics data has alleviated many limitations that have hampered molecular convergence studies of non-model extremophile plants. In this review, we summarize cases of genomic convergence in these taxa to examine the extent and type of genomic convergence during the process of adaptation to extreme environments. Despite being well studied by candidate gene approaches, convergent evolution at individual sites is rare and often has a high false-positive rate when assessed in whole genomes. By contrast, genomic convergence at higher genetic levels has been detected during adaptation to the same extreme environments. Examples include the convergence of biological pathways and changes in gene expression, gene copy number, amino acid usage, and GC content. Higher convergence levels play important roles in the adaptive evolution of extremophiles and may be more frequent and involve more genes. In several cases, multiple types of convergence events have been found to co-occur. However, empirical and theoretical studies of this higher level convergent evolution are still limited. In conclusion, both the development of powerful approaches and the detection of convergence at various genetic levels are needed to further reveal the genetic mechanisms of plant adaptation to extreme environments.