Forces driving transposable element load variation during Arabidopsis range expansion.

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Arabidopsis EXECUTER1 interacts with WRKY transcription factors to mediate plastid-to-nucleus singlet oxygen signaling.

Chloroplasts produce singlet oxygen (1O2), which causes changes in nuclear gene expression through plastid-to-nucleus retrograde signaling to increase plant fitness. However, the identity of this 1O2-triggered pathway remains unclear. Here, we identify mutations in GENOMES UNCOUPLED4 (GUN4) and GUN5 as suppressors of phytochrome-interacting factor1 (pif1) pif3 in regulating the photo-oxidative response in Arabidopsis thaliana. GUN4 and GUN5 specifically interact with EXECUTER1 (EX1) and EX2 in plastids, and this interaction is alleviated by treatment with Rose Bengal (RB) or white light. Impaired expression of GUN4, GUN5, EX1, or EX2 leads to insensitivity to excess light and overexpression of EX1 triggers photo-oxidative responses. Strikingly, upon light irradiation or RB treatment, EX1 transiently accumulates in the nucleus and the nuclear fraction of EX1 shows a similar molecular weight as the plastid-located protein. Point mutagenesis analysis indicated that nuclear localization of EX1 is required for its function. EX1 acts as a transcriptional co-activator and interacts with the transcription factors WRKY18 and WRKY40 to promote the expression of 1O2-responsive genes. This study suggests that EX1 may act in plastid-to-nucleus signaling and establishes a 1O2-triggered retrograde signaling pathway that allows plants adapt to changing light environments during chloroplast development.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa.

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

COG3 confers the chilling tolerance to mediate OsFtsH2-D1 module in rice.

The chilling stress induced by the global climate change harms rice production, especially at seedling and booting stage, which feed half the population of the world. Although there are key quantitative trait locus genes identified in the individual stage, few genes have been reported and functioned at both stages. Utilizing chromosome segment substitution lines (CSSLs) and a combination of map-based cloning and phenotypes of the mutants and overexpression lines, we identified the major gene Chilling-tolerance in Geng/japonica rice 3 (COG3) of q chilling-tolerance at the booting and seedling stage 11 (qCTBS11) conferred chilling tolerance at both seedling and booting stages. COG3 was significantly upregulated in Nipponbare under chilling treatment compared with its expression in 93-11. The loss-of-function mutants cog3 showed a reduced chilling tolerance. On the contrary, overexpression enhanced chilling tolerance. Genome evolution and genetic analysis suggested that COG3 may have undergone strong selection in temperate japonica during domestication. COG3, a putative calmodulin-binding protein, physically interacted with OsFtsH2 at chloroplast. In cog3-1, OsFtsH2-mediated D1 degradation was impaired under chilling treatment compared with wild-type. Our results suggest that COG3 is necessary for maintaining OsFtsH2 protease activity to regulate chilling tolerance at the booting and seedling stage.

Polygenic adaptation of rosette growth in Arabidopsis thaliana.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its evolution in natural populations. Here, we investigate evolutionary and environmental factors shaping variation in the growth rate of Arabidopsis thaliana. We used plant diameter as a proxy to monitor plant growth over time in environments that mimicked latitudinal differences in the intensity of natural light radiation, across a set of 278 genotypes sampled within four broad regions, including an outgroup set of genotypes from China. A field experiment conducted under natural conditions confirmed the ecological relevance of the observed variation. All genotypes markedly expanded their rosette diameter when the light supply was decreased, demonstrating that environmental plasticity is a predominant source of variation to adapt plant size to prevailing light conditions. Yet, we detected significant levels of genetic variation both in growth rate and growth plasticity. Genome-wide association studies revealed that only 2 single nucleotide polymorphisms associate with genetic variation for growth above Bonferroni confidence levels. However, marginally associated variants were significantly enriched among genes with an annotated role in growth and stress reactions. Polygenic scores computed from marginally associated variants confirmed the polygenic basis of growth variation. For both light regimes, phenotypic divergence between the most distantly related population (China) and the various regions in Europe is smaller than the variation observed within Europe, indicating that the evolution of growth rate is likely to be constrained by stabilizing selection. We observed that Spanish genotypes, however, reach a significantly larger size than Northern European genotypes. Tests of adaptive divergence and analysis of the individual burden of deleterious mutations reveal that adaptive processes have played a more important role in shaping regional differences in rosette growth than maladaptive evolution.

Weak transgenerational effects of ancestral nitrogen and phosphorus availabilities on offspring phenotypes in Arabidopsis thaliana.

Nutrient availability significantly regulates plant growth and metabolic functions, but whether and how the long-term exposure of ancestral plants to contrasting nutrient environments influences offspring phenotypic performance (i.e., transgenerational plasticity) remain poorly addressed. Here we conducted experimental manipulations using Arabidopsis thaliana with the ancestral plants grown in different nitrogen (N) and phosphorus (P) availabilities over eleven consecutive generations, and then examined the offspring phenotypic performance under the interactive effects of current and ancestral nutrient environments. We found that current rather than ancestral nutrient environments dominantly explained the variations in offspring plant traits (i.e., flowering time, aboveground biomass and biomass allocation fractions), suggesting the relatively weak transgenerational effects of ancestral N and P availabilities on offspring phenotypes. In contrast, increasing N and P availabilities in the offspring generation remarkably shortened the flowering time, increased the aboveground biomass, and altered biomass allocation fractions differentially among organs. Despite the overall weak transgenerational phenotypic plasticity, under the low nutrient environment, the offspring of ancestral plants from the low nutrient environment had a significantly higher fruit mass fraction than those from the suitable nutrient environment. Taken together, our findings suggest that A. thaliana exhibits a much stronger within- than trans-generational trait plasticity under contrasting nutrient availabilities, and may provide important insights into the understanding of plant adaptation and evolutionary processes under changing nutrient environments.

Transcriptional repression specifies the central cell for double fertilization.

Double fertilization is a key innovation for the evolutionary success of angiosperms by which the two fertilized female gametes, the egg cell and central cell, generate the embryo and endosperm, respectively. The female gametophyte (embryo sac) enclosed in the sporophyte is derived from a one-celled haploid cell lineage. It undergoes successive events of mitotic divisions, cellularization, and cell specification to give rise to the mature embryo sac, which contains the two female gametes accompanied by two types of accessory cells, namely synergids and antipodals. How the cell fate of the central cell is specified has long been equivocal and is further complicated by the structural diversity of female gametophyte across plant taxa. Here, MADS-box protein AGL80 was verified as a transcriptional repressor that directly suppresses the expression of accessory cell-specific genes to specify the central cell. Further genetic rescue and phylogenetic assay of the AGL80 orthologs revealed a possible conserved mechanism in the Brassicaceae family. Results from this study provide insight into the molecular determination of the second female gamete cell in Brassicaceae.

Genomes shed light on the evolution of Begonia, a mega-diverse genus.

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

Scaling-up and proteomic analysis reveals photosynthetic and metabolic insights toward prolonged H2 photoproduction in Chlamydomonas hpm91 mutant lacking proton gradient regulation 5 (PGR5).

Clean and sustainable H2 production is crucial to a carbon-neutral world. H2 generation by Chlamydomonas reinhardtii is an attractive approach for solar-H2 from H2O. However, it is currently not large-scalable because of lacking desirable strains with both optimal H2 productivity and sufficient knowledge of underlying molecular mechanism. We hereby carried out extensive and in-depth investigations of H2 photoproduction of hpm91 mutant lacking PGR5 (Proton Gradient Regulation 5) toward its up-scaling and fundamental mechanism issues. We show that hpm91 is at least 100-fold scalable (up to 10 L) with continuous H2 collection of 7287 ml H2/10L-HPBR in averagely 26 days under sulfur deprivation. Also, we show that hpm91 is robust and active during sustained H2 photoproduction, most likely due to decreased intracellular ROS relative to wild type. Moreover, we obtained quantitative proteomic profiles of wild type and hpm91 at four representing time points of H2 evolution, leading to 2229 and 1350 differentially expressed proteins, respectively. Compared to wild type, major proteome alterations of hpm91 include not only core subunits of photosystems and those related to anti-oxidative responses but also essential proteins in photosynthetic antenna, C/N metabolic balance, and sulfur assimilation toward both cysteine biosynthesis and sulfation of metabolites during sulfur-deprived H2 production. These results reveal not only new insights of cellular and molecular basis of enhanced H2 production in hpm91 but also provide additional candidate gene targets and modules for further genetic modifications and/or in artificial photosynthesis mimics toward basic and applied research aiming at advancing solar-H2 technology.

Adaptation and Phenotypic Diversification in Arabidopsis through Loss-of-Function Mutations in Protein-Coding Genes.

According to the less-is-more hypothesis, gene loss is an engine for evolutionary change. Loss-of-function (LoF) mutations resulting in the natural knockout of protein-coding genes not only provide information about gene function but also play important roles in adaptation and phenotypic diversification. Although the less-is-more hypothesis was proposed two decades ago, it remains to be explored on a large scale. In this study, we identified 60,819 LoF variants in 1071 Arabidopsis (Arabidopsis thaliana) genomes and found that 34% of Arabidopsis protein-coding genes annotated in the Columbia-0 genome do not have any LoF variants. We found that nucleotide diversity, transposable element density, and gene family size are strongly correlated with the presence of LoF variants. Intriguingly, 0.9% of LoF variants with minor allele frequency larger than 0.5% are associated with climate change. In addition, in the Yangtze River basin population, 1% of genes with LoF mutations were under positive selection, providing important insights into the contribution of LoF mutations to adaptation. In particular, our results demonstrate that LoF mutations shape diverse phenotypic traits. Overall, our results highlight the importance of the LoF variants for the adaptation and phenotypic diversification of plants.

Transposable elements drive rapid phenotypic variation in Capsella rubella.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3' UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Natural variation in the HAN1 gene confers chilling tolerance in rice and allowed adaptation to a temperate climate.

Rice (Oryza sativa L.) is a chilling-sensitive staple crop that originated in subtropical regions of Asia. Introduction of the chilling tolerance trait enables the expansion of rice cultivation to temperate regions. Here we report the cloning and characterization of HAN1, a quantitative trait locus (QTL) that confers chilling tolerance on temperate japonica rice. HAN1 encodes an oxidase that catalyzes the conversion of biologically active jasmonoyl-L-isoleucine (JA-Ile) to the inactive form 12-hydroxy-JA-Ile (12OH-JA-Ile) and fine-tunes the JA-mediated chilling response. Natural variants in HAN1 diverged between indica and japonica rice during domestication. A specific allele from temperate japonica rice, which gained a putative MYB cis-element in the promoter of HAN1 during the divergence of the two japonica ecotypes, enhances the chilling tolerance of temperate japonica rice and allows it to adapt to a temperate climate. The results of this study extend our understanding of the northward expansion of rice cultivation and provide a target gene for the improvement of chilling tolerance in rice.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara).

BACKGROUND: Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown. RESULTS: Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding-noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding-noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding-noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare. CONCLUSION: We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Parallel Evolution of Common Allelic Variants Confers Flowering Diversity in Capsella rubella.

Flowering time is an adaptive life history trait. Capsella rubella, a close relative of Arabidopsis thaliana and a young species, displays extensive variation for flowering time but low standing genetic variation due to an extreme bottleneck event, providing an excellent opportunity to understand how phenotypic diversity can occur with a limited initial gene pool. Here, we demonstrate that common allelic variation and parallel evolution at the FLC locus confer variation in flowering time in C. rubella. We show that two overlapping deletions in the 5' untranslated region (UTR) of C. rubella FLC, which are associated with local changes in chromatin conformation and histone modifications, reduce its expression levels and promote flowering. We further show that these two pervasive variants originated independently in natural C. rubella populations after speciation and spread to an intermediate frequency, suggesting a role of this parallel cis-regulatory change in adaptive evolution. Our results provide an example of how parallel mutations in the same 5' UTR region can shape phenotypic evolution in plants.

Parallel Speciation of Wild Rice Associated with Habitat Shifts.

The occurrence of parallel speciation strongly implies the action of natural selection. However, it is unclear how general a phenomena parallel speciation is since it was only shown in a small number of animal species. In particular, the adaptive process and mechanisms underlying the process of parallel speciation remain elusive. Here, we used an integrative approach incorporating population genomics, common garden, and crossing experiments to investigate parallel speciation of the wild rice species Oryza nivara from O. rufipogon. We demonstrated that O. nivara originated multiple times from different O. rufipogon populations and revealed that different O. nivara populations have evolved similar phenotypes under divergent selection, a reflection of recurrent local adaptation of ancient O. rufipogon populations to dry habitats. Almost completed premating isolation was detected between O. nivara and O. rufipogon in the absence of any postmating barriers between and within these species. These results suggest that flowering time is a "magic" trait that contributes to both local adaptation and reproductive isolation in the origin of wild rice species. Our study thus demonstrates a convincing case of parallel ecological speciation as a consequence of adaptation to new environments.

Synthesis, physicochemical properties, and health aspects of structured lipids: A review.

Structured lipids (SLs) refer to a new type of functional lipids obtained by chemically, enzymatically, or genetically modifying the composition and/or distribution of fatty acids in the glycerol backbone. Due to the unique physicochemical characteristics and health benefits of SLs (for example, calorie reduction, immune function improvement, and reduction in serum triacylglycerols), there is increasing interest in the research and application of novel SLs in the food industry. The chemical structures and molecular architectures of SLs define mainly their physicochemical properties and nutritional values, which are also affected by the processing conditions. In this regard, this holistic review provides coverage of the latest developments and applications of SLs in terms of synthesis strategies, physicochemical properties, health aspects, and potential food applications. Enzymatic synthesis of SLs particularly with immobilized lipases is presented with a short introduction to the genetic engineering approach. Some physical features such as solid fat content, crystallization and melting behavior, rheology and interfacial properties, as well as oxidative stability are discussed as influenced by chemical structures and processing conditions. Health-related considerations of SLs including their metabolic characteristics, biopolymer-based lipid digestion modulation, and oleogelation of liquid oils are also explored. Finally, potential food applications of SLs are shortly introduced. Major challenges and future trends in the industrial production of SLs, physicochemical properties, and digestion behavior of SLs in complex food systems, as well as further exploration of SL-based oleogels and their food application are also discussed.

Rice FLUORESCENT1 Is Involved in the Regulation of Chlorophyll.

Chlorophyll biosynthesis plays essential roles in photosynthesis and plant growth in response to environmental conditions. The accumulation of excess chlorophyll biosynthesis intermediates under light results in the production of reactive oxygen species and oxidative stress. In this study, we identified a rice (Oryza sativa) mutant, oxidation under photoperiod (oxp), that displayed photobleached lesions on its leaves, reduced growth and decreased chlorophyll content during light/dark cycles or following a dark-to-light transition. The oxp mutant accumulated more chlorophyll precursors (5-aminolevulinic acid and protochlorophyllide) than the wild type in the dark, and more singlet oxygen following light exposure. Several singlet-oxygen-responsive genes were greatly upregulated in oxp, whereas the expression patterns of OsPORA and OsPORB, two genes encoding the chlorophyll biosynthesis enzyme NADPH:protochlorop hyllide oxidoreductase, were altered in de-etiolated oxp seedlings. Molecular and complementation studies revealed that oxp is a loss-of-function mutant in LOC_Os01g32730, a homolog of FLUORESCENT (FLU) in Arabidopsis thaliana. Rice PHYTOCHROME-INTERACTING FACTOR-LIKE14 (OsPIL14) transcription factor directly bound to the OsFLU1 promoter and activated its expression. Dark-grown transgenic rice seedlings overexpressing OsPIL14 accumulated more chlorophyll and turned green faster than the wild type upon light illumination. Thus, OsFLU1 is an important regulator of chlorophyll biosynthesis in rice.

Effects of nitrogen and phosphorus supply on stoichiometry of six elements in leaves of Arabidopsis thaliana.

BACKGROUND AND AIMS: Plant elemental composition is of fundamental importance for plant growth and metabolic functions. However, knowledge of how multi-elemental stoichiometry responds to varying nitrogen (N) and phosphorus (P) availabilities remains limited. METHODS: We conducted experimental manipulations with nine repeat experiments to investigate the effects of N and P supply on the concentrations and variability of six elements, carbon (C), N, P, potassium (K), calcium (Ca) and magnesium (Mg), in leaves of Arabidopsis thaliana. KEY RESULTS: N supply increased the concentrations of N, K and Mg, decreased the concentration of P, but exerted little influence on the concentrations of C and Ca in green leaves. P supply increased the concentrations of P and Ca, decreased the concentration of C, initially increased and then decreased the concentration of K, but showed little influence on the concentrations of N and Mg in green leaves. Multivariate patterns among the concentrations of these six elements in green leaves was influenced by the type of nutrient supply (i.e. N or P). Elemental variability decreased with increasing elemental concentrations in green leaves at the intraspecific level, supporting the Stability of Limiting Elements Hypothesis that was originally proposed from a meta-analysis of pooled data across species or communities. Compared with green leaves, the senesced leaves showed greater variability in C, N, P, K and Mg concentrations but lower variability in Ca concentration. CONCLUSIONS: N and P supplies exerted differential influences on the concentrations of C, N, P, K, Ca and Mg in green leaves. The specific C content should be considered when assessing C cycling under global nutrient changes. Stage-dependent patterns of leaf stoichiometric homeostasis differed among elements with various chemical characteristics. These findings can help to improve our understanding of plant eco-physiological responses and acclimation under global nutrient changes from the stoichiometric perspective of multiple elements.

A Novel N-Methyltransferase in Arabidopsis Appears to Feed a Conserved Pathway for Nicotinate Detoxification among Land Plants and Is Associated with Lignin Biosynthesis.

The Preiss-Handler pathway, which salvages nicotinate (NA) for NAD synthesis, is an indispensable biochemical pathway in land plants. Various NA conjugations (mainly methylation and glycosylation) have been detected and have long been proposed for NA detoxification in plants. Previously, we demonstrated that NA O-glucosylation functions as a mobilizable storage form for NAD biosynthesis in the Brassicaceae. However, little is known about the functions of other NA conjugations in plants. In this study, we first found that N-methylnicotinate is a ubiquitous NA conjugation in land plants. Furthermore, we functionally identified a novel methyltransferase (At3g53140; NANMT), which is mainly responsible for N-methylnicotinate formation, from Arabidopsis (Arabidopsis thaliana). We also established that trigonelline is a detoxification form of endogenous NA in plants. Combined phylogenetic analysis and enzymatic assays revealed that NA N-methylation activity was likely derived from the duplication and subfunctionalization of an ancestral caffeic acid O-methyltransferase (COMT) gene in the course of land plant evolution. COMT enzymes, which function in S-lignin biosynthesis, also have weak NANMT activity. Our data suggest that NA detoxification conferred by NANMT and COMT might have facilitated the retention of the Preiss-Handler pathway in land plants.