Structure, function, and evolution of plant O-methyltransferases

The biosynthetic pathway of the hallucinogen mescaline and its heterologous reconstruction

Mescaline, among the earliest identified natural hallucinogens, holds great potential in psychotherapy treatment. Nonetheless, despite the existence of a postulated biosynthetic pathway for more than half a century, the specific enzymes involved in this process are yet to be identified. In this study, we investigated the cactus Lophophora williamsii (Peyote), the largest known natural producer of the phenethylamine mescaline. We employed a multi-faceted approach, combining de novo whole-genome and transcriptome sequencing with comprehensive chemical profiling, enzymatic assays, molecular modeling, and pathway engineering for pathway elucidation. We identified four groups of enzymes responsible for the six catalytic steps in the mescaline biosynthetic pathway, and an N-methyltransferase enzyme that N-methylates all phenethylamine intermediates, likely modulating mescaline levels in Peyote. Finally, we reconstructed the mescaline biosynthetic pathway in both Nicotiana benthamiana plants and yeast cells, providing novel insights into several challenges hindering complete heterologous mescaline production. Taken together, our study opens up avenues for exploration of sustainable production approaches and responsible utilization of mescaline, safeguarding this valuable natural resource for future generations.

Identification of O-Methyltransferases Potentially Contributing to the Structural Diversity of 2-(2-Phenylethyl)chromones in Agarwood

2-(2-Phenylethyl)chromones (PECs) are the primary constituents responsible for the promising pharmacological activities and unique fragrance of agarwood. However, the O-methyltransferases (OMTs) involved in the formation of diverse methylated PECs have not been reported. In this study, we identified one Mg2+-dependent caffeoyl-CoA-OMT subfamily enzyme (AsOMT1) and three caffeic acid-OMT subfamily enzymes (AsOMT2-4) from NaCl-treated Aquilaria sinensis calli. AsOMT1 not only converts caffeoyl-CoA to feruloyl-CoA but also performs nonregioselective methylation at either the 6-OH or 7-OH position of 6,7-dihydroxy-PEC. On the other hand, AsOMT2-4 preferentially utilizes PECs as substrates to produce structurally diverse methylated PECs. Additionally, AsOMT2-4 also accepts nonPEC-type substrates such as caffeic acid and apigenin to generate methylated products. Protein structure prediction and site-directed mutagenesis revealed that residues of L313 and I318 in AsOMT3, as well as S292 and F313 in AsOMT4 determine the distinct regioselectivity of these two OMTs toward apigenin. These findings provide important biochemical evidence of the remarkable structural diversity of PECs in agarwood.

Specialized metabolite modifications in Brassicaceae seeds and plants: diversity, functions and related enzymes

Covering: up to 2023Specialized metabolite (SM) modifications and/or decorations, corresponding to the addition or removal of functional groups (e.g. hydroxyl, methyl, glycosyl or acyl group) to SM structures, contribute to the huge diversity of structures, activities and functions of seed and plant SMs. This review summarizes available knowledge (up to 2023) on SM modifications in Brassicaceae and their contribution to SM plasticity. We give a comprehensive overview on enzymes involved in the addition or removal of these functional groups. Brassicaceae, including model (Arabidopsis thaliana) and crop (Brassica napus, Camelina sativa) plant species, present a large diversity of plant and seed SMs, which makes them valuable models to study SM modifications. In this review, particular attention is given to the environmental plasticity of SM and relative modification and/or decoration enzymes. Furthermore, a spotlight is given to SMs and related modification enzymes in seeds of Brassicaceae species. Seeds constitute a large reservoir of beneficial SMs and are one of the most important dietary sources, providing more than half of the world's intake of dietary proteins, oil and starch. The seed tissue- and stage-specific expressions of A. thaliana genes involved in SM modification are presented and discussed in the context of available literature. Given the major role in plant phytochemistry, biology and ecology, SM modifications constitute a subject of study contributing to the research and development in agroecology, pharmaceutical, cosmetics and food industrial sectors.

Discovery of a novel flavonol O-methyltransferase possessing sequential 4'- and 7-O-methyltransferase activity from Camptotheca acuminata Decne

The biosynthetic route for flavonol in Camptotheca acuminata has been recently elucidated from a chemical point of view. However, the genes involved in flavonol methylation remain unclear. It is a critical step for fully uncovering the flavonol metabolism in this ancient plant. In this study, the multi-omics resource of this plant was utilized to perform flavonol O-methyltransferase-oriented mining and screening. Two genes, CaFOMT1 and CaFOMT2 are identified, and their recombinant CaFOMT proteins are purified to homogeneity. CaFOMT1 exhibits strict substrate and catalytic position specificity for quercetin, and selectively methylates only the 4'-OH group. CaFOMT2 possesses sequential O-methyltransferase activity for the 4'-OH and 7-OH of quercetin. These CaFOMT genes are enriched in the leaf and root tissues. The catalytic dyad and critical substrate-binding sites of the CaFOMTs are determined by molecular docking and further verified through site-mutation experiments. PHE181 and MET185 are designated as the critical sites for flavonol substrate selectivity. Genomic environment analysis indicates that CaFOMTs evolved independently and that their ancestral genes are different from that of the known Ca10OMT. This study provides molecular insights into the substrate-binding pockets of two new CaFOMTs responsible for flavonol metabolism in C. acuminata.

Integrative Dissection of Lignin Composition in Tartary Buckwheat Seed Hulls for Enhanced Dehulling Efficiency

The rigid hull encasing Tartary buckwheat seeds necessitates a laborious dehulling process before flour milling, resulting in considerable nutrient loss. Investigation of lignin composition is pivotal in understanding the structural properties of tartary buckwheat seeds hulls, as lignin is key determinant of rigidity in plant cell walls, thus directly impacting the dehulling process. Here, the lignin composition of seed hulls from 274 Tartary buckwheat accessions is analyzed, unveiling a unique lignin chemotype primarily consisting of G lignin, a common feature in gymnosperms. Furthermore, the hardness of the seed hull showed a strong negative correlation with the S lignin content. Genome-wide detection of selective sweeps uncovered that genes governing the biosynthesis of S lignin, specifically two caffeic acid O-methyltransferases (COMTs) and one ferulate 5-hydroxylases, are selected during domestication. This likely contributed to the increased S lignin content and decreased hardness of seed hulls from more domesticated varieties. Genome-wide association studies identified robust associations between FtCOMT1 and the accumulation of S lignin in seed hull. Transgenic Arabidopsis comt1 plants expressing FtCOMT1 successfully reinstated S lignin content, confirming its conserved function across plant species. These findings provide valuable metabolic and genetic insights for the potential redesign of Tartary buckwheat seed hulls.

A Novel O- and S-Methyltransferase from Pleurotus sapidus Is Involved in Flavor Formation

Increasing consumer aversion to non-natural flavoring substances is prompting a heightened interest in enzymatic processes for flavor production. This includes methylation reactions, which are often performed by using hazardous chemicals. By correlation of aroma profile data and transcriptomic analysis, a novel O-methyltransferase (OMT) catalyzing a respective reaction within the formation of p-anisaldehyde was identified in the mushroom Pleurotus sapidus. Heterologous expression in E. coli followed by purification allowed for further characterization of the enzyme. Besides p-hydroxybenzaldehyde, the proposed precursor of p-anisaldehyde, the enzyme catalyzed the methylation of further hydroxylated aromatic compounds at the meta- and para-position. The Km values determined for p-hydroxybenzaldehyde and S-adenosyl-l-methionine were 80 and 107 μM, respectively. Surprisingly, the studied enzyme enabled the transmethylation of thiol-nucleophiles, as indicated by the formation of 2-methyl-3-(methylthio)furan from 2-methyl-3-furanthiol. Moreover, the enzyme was crystallized at a resolution of 2.0 Å, representing the first published crystal structure of a basidiomycetous OMT.

Genetic diversity and candidate genes for transient waterlogging tolerance in mungbean at the germination and seedling stages

Mungbean [Vigna radiata var. radiata (L.) Wilczek] production in Asia is detrimentally affected by transient soil waterlogging caused by unseasonal and increasingly frequent extreme precipitation events. While mungbean exhibits sensitivity to waterlogging, there has been insufficient exploration of germplasm for waterlogging tolerance, as well as limited investigation into the genetic basis for tolerance to identify valuable loci. This research investigated the diversity of transient waterlogging tolerance in a mini-core germplasm collection of mungbean and identified candidate genes for adaptive traits of interest using genome-wide association studies (GWAS) at two critical stages of growth: germination and seedling stage (i.e., once the first trifoliate leaf had fully-expanded). In a temperature-controlled glasshouse, 292 genotypes were screened for tolerance after (i) 4 days of waterlogging followed by 7 days of recovery at the germination stage and (ii) 8 days of waterlogging followed by 7 days of recovery at the seedling stage. Tolerance was measured against drained controls. GWAS was conducted using 3,522 high-quality DArTseq-derived SNPs, revealing five significant associations with five phenotypic traits indicating improved tolerance. Waterlogging tolerance was positively correlated with the formation of adventitious roots and higher dry masses. FGGY carbohydrate kinase domain-containing protein was identified as a candidate gene for adventitious rooting and mRNA-uncharacterized LOC111241851, Caffeoyl-CoA O-methyltransferase At4g26220 and MORC family CW-type zinc finger protein 3 and zinc finger protein 2B genes for shoot, root, and total dry matter production. Moderate to high broad-sense heritability was exhibited for all phenotypic traits, including seed emergence (81%), adventitious rooting (56%), shoot dry mass (81%), root dry mass (79%) and SPAD chlorophyll content (70%). The heritability estimates, marker-trait associations, and identification of sources of waterlogging tolerant germplasm from this study demonstrate high potential for marker-assisted selection of tolerance traits to accelerate breeding of climate-resilient mungbean varieties.

Whole genome identification, molecular docking and expression analysis of enzymes involved in the selenomethionine cycle in Cardamine hupingshanensis

BACKGROUND

The selenomethionine cycle (SeMTC) is a crucial pathway for the metabolism of selenium. The basic bioinformatics and functions of four enzymes involved in the cycle including S-adenosyl-methionine synthase (MAT), SAM-dependent methyltransferase (MTase), S-adenosyl-homocysteine hydrolase (SAHH) and methionine synthase (MTR), have been extensively reported in many eukaryotes. The identification and functional analyses of SeMTC genes/proteins in Cardamine hupingshanensis and their response to selenium stress have not yet been reported.

RESULTS

In this study, 45 genes involved in SeMTC were identified in the C. hupingshanensis genome. Phylogenetic analysis showed that seven genes from ChMAT were clustered into four branches, twenty-seven genes from ChCOMT were clustered into two branches, four genes from ChSAHH were clustered into two branches, and seven genes from ChMTR were clustered into three branches. These genes were resided on 16 chromosomes. Gene structure and homologous protein modeling analysis illustrated that proteins in the same family are relatively conserved and have similar functions. Molecular docking showed that the affinity of SeMTC enzymes for selenium metabolites was higher than that for sulfur metabolites. The key active site residues identified for ChMAT were Ala269 and Lys273, while Leu221/231 and Gly207/249 were determined as the crucial residues for ChCOMT. For ChSAHH, the essential active site residues were found to be Asn87, Asp139 and Thr206/207/208/325. Ile204, Ser111/329/377, Asp70/206/254, and His329/332/380 were identified as the critical active site residues for ChMTR. In addition, the results of the expression levels of four enzymes under selenium stress revealed that ChMAT3-1 genes were upregulated approximately 18-fold, ChCOMT9-1 was upregulated approximately 38.7-fold, ChSAHH1-2 was upregulated approximately 11.6-fold, and ChMTR3-2 genes were upregulated approximately 28-fold. These verified that SeMTC enzymes were involved in response to selenium stress to varying degrees.

CONCLUSIONS

The results of this research are instrumental for further functional investigation of SeMTC in C. hupingshanensis. This also lays a solid foundation for deeper investigations into the physiological and biochemical mechanisms underlying selenium metabolism in plants.

O- and N-Methyltransferases in benzylisoquinoline alkaloid producing plants

BACKGROUND

Secondary metabolites such as benzylisoquinoline alkaloids (BIA) have attracted considerable attention because of their pharmacological properties and potential therapeutic applications. Methyltransferases (MTs) can add methyl groups to alkaloid molecules, altering their physicochemical properties and bioactivity, stability, solubility, and recognition by other cellular components. Five types of O-methyltransferases and two types of N-methyltransferases are involved in BIA biosynthesis.

OBJECTIVE

Since MTs may be the source for the discovery and development of novel biomedical, agricultural, and industrial compounds, we performed extensive molecular and phylogenetic analyses of O- and N-methyltransferases in BIA-producing plants.

METHODS

MTs involved in BIA biosynthesis were isolated from transcriptomes of Berberis koreana and Caulophyllum robustum. We also mined the methyltransferases of Coptis japonica, Papaver somniferum, and Nelumbo nucifera from the National Center for Biotechnology Information protein database. Then, we analyzed the functional motifs and phylogenetic analysis.

RESULT

We mined 42 O-methyltransferases and 8 N-methyltransferases from the five BIA-producing plants. Functional motifs for S-adenosyl-L-methionine-dependent methyltransferases were retained in most methyltransferases, except for the three O-methyltransferases from N. nucifera. Phylogenetic analysis revealed that the methyltransferases were grouped into four clades, I, II, III and IV. The clustering patterns in the phylogenetic analysis suggested a monophyletic origin of methyltransferases and gene duplication within species. The coexistence of different O-methyltransferases in the deep branch subclade might support some cases of substrate promiscuity.

CONCLUSIONS

Methyltransferases may be a source for the discovery and development of novel biomedical, agricultural, and industrial compounds. Our results contribute to further understanding of their structure and reaction mechanisms, which will require future functional studies.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation

BACKGROUND

Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family.

FINDINGS

The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species.

CONCLUSIONS

The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Genome-wide identification of the OMT gene family in Cucumis melo L. and expression analysis under abiotic and biotic stress

Background

O-methyltransferase (OMT)-mediated O-methylation is a frequent modification that occurs during natural product biosynthesis, and it increases the diversity and stability of secondary metabolites. However, detailed genome-wide identification and expression analyses of OMT gene family members have not been performed in melons. In this study, we aimed to perform the genome-wide identification of OMT gene family members in melon to identify and clarify their actions during stress.

Methods

Genome-wide identification of OMT gene family members was performed using data from the melon genome database. The Cucumis melo OMT genes (CmOMTs) were then compared with the genes from two representative monocotyledons and three representative dicotyledons. The basic information, cis-regulatory elements in the promoter, predicted 3-D-structures, and GO enrichment results of the 21 CmOMTs were analyzed.

Results

In our study, 21 CmOMTs (named CmOMT1-21) were obtained by analyzing the melon genome. These genes were located on six chromosomes and divided into three groups composed of nine, six, and six CmOMTs based on phylogenetic analysis. Gene structure and motif descriptions were similar within the same classes. Each CmOMT gene contains at least one cis-acting element associated with hormone transport regulation. Analysis of cis-acting elements illustrated the potential role of CmOMTs in developmental regulation and adaptations to various abiotic and biotic stresses. The RNA-seq and quantitative real-time PCR (qRT-PCR) results indicated that NaCl stress significantly induced CmOMT6/9/14/18 and chilling and high temperature and humidity (HTH) stresses significantly upregulated CmOMT14/18. Furthermore, the expression pattern of CmOMT18 may be associated with Fusarium oxysporum f. sp. melonis race 1.2 (FOM1.2) and powdery mildew resistance. Our study tentatively explored the biological functions of CmOMT genes in various stress regulation pathways and provided a conceptual basis for further detailed studies of the molecular mechanisms.

Diverse O-methyltransferases catalyze the biosynthesis of floral benzenoids that repel aphids from the flowers of waterlily Nymphaea prolifera

Nymphaea is a key genus of the ANA grade (Amborellales, Nymphaeales, and Austrobaileyales) of basal flowering plants, which serve as a key model to study the early evolution of floral traits. In this study, we comprehensively investigated the emission, biosynthesis, and biological function of the floral scent in a night-blossoming waterlily Nymphaea prolifera. The headspace volatile collection combined with GC-MS analysis showed that the floral scent of N. prolifera is predominately comprised by methylated benzenoids including anisole, veratrole, guaiacol, and methoxyanisole. Moreover, the emission of these floral benzenoids in N. prolifera exhibited temporal and spatial pattern with circadian rhythm and tissue specificity. By creating and mining transcriptomes of N. prolifera flowers, 12 oxygen methyltransferases (NpOMTs) were functionally identified. By in vitro enzymatic assay, NpOMT3, 6, and 7 could produce anisole and NpOMT5, 7, 9, produce guaiacol, whereas NpOMT3, 6, 9, 11 catalyzed the formation of veratrole. Methoxyanisole was identified as the universal product of all NpOMTs. Expression patterns of NpOMTs provided implication for their roles in the production of the respective benzenoids. Phylogenetic analysis of OMTs suggested a Nymphaea-specific expansion of the OMT family, indicating the evolution of lineage-specific functions. In bioassays, anisole, veratrole, and guaiacol in the floral benzenoids were revealed to play the critical role in repelling waterlily aphids. Overall, this study indicates that the basal flowering plant N. prolifera has evolved a diversity and complexity of OMT genes for the biosynthesis of methylated benzenoids that can repel insects from feeding the flowers. These findings provide new insights into the evolutional mechanism and ecological significance of the floral scent from early-diverged flowering plants.

Distinct heat response molecular mechanisms emerge in cassava vasculature compared to leaf mesophyll tissue under high temperature stress

With growing concerns over global warming, cultivating heat-tolerant crops has become paramount to prepare for the anticipated warmer climate. Cassava (Manihot esculenta Crantz), a vital tropical crop, demonstrates exceptional growth and productivity under high-temperature (HT) conditions. Yet, studies elucidating HT resistance mechanisms in cassava, particularly within vascular tissues, are rare. We dissected the leaf mid-vein from leaf, and did the comparative transcriptome profiling between mid-vein and leaf to figure out the cassava vasculature HT resistance molecular mechanism. Anatomical microscopy revealed that cassava leaf veins predominantly consisted of vasculature. A thermal imaging analysis indicated that cassava experienced elevated temperatures, coinciding with a reduction in photosynthesis. Transcriptome sequencing produced clean reads in total of 89.17G. Using Venn enrichment, there were 65 differentially expressed genes (DEGs) and 93 DEGs had been found highly specifically expressed in leaf and mid-vein. Further investigation disclosed that leaves enhanced pyruvate synthesis as a strategy to withstand high temperatures, while mid-veins fortified themselves by bolstering lignin synthesis by comprehensive GO and KEGG analysis of DEGs. The identified genes in these metabolic pathways were corroborated through quantity PCR (QPCR), with results aligning with the transcriptomic data. To verify the expression localization of DEGs, we used in situ hybridization experiments to identify the expression of MeCCoAMT(caffeoyl-coenzyme A-3-O-methyltransferase) in the lignin synthesis pathway in cassava leaf veins xylem. These findings unravel the disparate thermotolerance mechanisms exhibited by cassava leaves and mid-veins, offering insights that could potentially inform strategies for enhancing thermotolerance in other crops.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

Biosynthesis of kratom opioids

Mitragynine, an analgesic alkaloid from the plant Mitragyna speciosa (kratom), offers a safer alternative to clinical opioids such as morphine, owing to its more favorable side effect profile. Although kratom has been traditionally used for stimulation and pain management in Southeast Asia, the mitragynine biosynthesis pathway has remained elusive. We embarked on a search for mitragynine biosynthetic genes from the transcriptomes of kratom and other members of the Rubiaceae family. We studied their functions in vitro and in vivo. Our investigations led to the identification of several reductases and an enol methyltransferase that forms a new clade within the SABATH methyltransferase family. Furthermore, we discovered a methyltransferase from Hamelia patens (firebush), which catalyzes the final step. With the tryptamine 4-hydroxylase from the psychedelic mushroom Psilocybe cubensis, we accomplished the four-step biosynthesis for mitragynine and its stereoisomer, speciogynine in both yeast and Escherichia coli when supplied with tryptamine and secologanin. Although we have yet to pinpoint the authentic hydroxylase and methyltransferase in kratom, our discovery completes the mitragynine biosynthesis. Through these breakthroughs, we achieved the microbial biosynthesis of kratom opioids for the first time. The remarkable enzyme promiscuity suggests the possibility of generating derivatives and analogs of kratom opioids in heterologous systems.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using F_ST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.

Deciphering the genomic regions governing major agronomic traits and blast resistance using genome wide association mapping in finger millet

Finger millet (Eleusine coracana L.) is climate resilient minor millet of Asia and Africa with wide adaptation and unparallel nutritional profile. To date, genomic resources available in finger millet are scanty and genetic control of agronomic traits remains elusive. Here, a collection of eco-geographically diverse 186 genotypes was quantified for variation in 13 agronomic traits and reaction to blast to identify marker-trait associations (MTAs) using genotyping-by-sequencing (GBS) and genome-wide association study (GWAS). GBS generated 2977 high quality single nucleotide polymorphism (SNPs) markers and identified three subpopulations with varying admixture levels. General linear and mixed model approaches of GWAS to correct for population structure and genetic relatedness identified 132 common MTAs for agronomic traits across the years. The phenotypic variance explained by the makers varied from 4.8% (TP692389-flag leaf width) to 20% (TP714446-green fodder weight). Of these, 26 MTAs showed homology with candidate genes having role in plant growth, development and photosynthesis in the genomes of foxtail millet, rice, maize, wheat and barley. We also found 4 common MTAs for neck blast resistance, which explained 5.9-15.1% phenotypic variance. Three MTAs for neck blast resistance showed orthologues in related genera having putative functions in pathogen defense in plants. The results of this work lay a foundation for understanding the genetic architecture of agronomic traits and blast resistance in finger millet and provide a framework for genomics assisted breeding.

Two O-Methyltransferases from Phylogenetically Unrelated Cow Parsley (Anthriscus sylvestris) and Hinoki-Asunaro (Thujopsis dolabrata var. hondae) as a Signature of Lineage-Specific Evolution in Lignan Biosynthesis

O-Methyltransferases (OMTs) play important roles in antitumor lignan biosynthesis. To date, six OMTs catalyzing the methylation of dibenzylbutyrolactone lignans as biosynthetic precursors of antitumor lignans have been identified. However, there is still no systematic understanding of the diversity and regularity of the biosynthetic mechanisms among various plant lineages. Herein, we report the characterization of two OMTs from Anthriscus sylvestris and Thujopsis dolabrata var. hondae [designated as AsSecoNorYatein (SNY) OMT and TdSNYOMT] together with the six known OMTs to evaluate their diversity and regularity. Although A. sylvestris 5-O-methylthujaplicatin (SecoNorYatein) and 4-O-demethylyatein (NorYatein) OMT (AsSNYOMT) and TdSNYOMT accept 5-O-methylthujaplicatin and 4-O-demethylyatein as substrates, phylogenetic analysis indicated that these two OMTs shared low amino acid sequence identity, 33.8%, indicating a signature of parallel evolution. The OMTs and the six previously identified OMTs were found to be diverse in terms of their substrate specificity, regioselectivity and amino acid sequence identity, indicating independent evolution in each plant species. Meanwhile, two-entropy analysis detected four amino acid residues as being specifically acquired by dibenzylbutyrolactone lignan OMTs. Site-directed mutation of AsSNYOMT indicated that two of them contributed specifically to 5-O-methylthujaplicatin methylation. The results provide a new example of parallel evolution and the diversity and regularity of OMTs in plant secondary (specialized) metabolism.

Discovery of specific catalytic activity toward IAA/FA by LaSABATHs based on genome-wide phylogenetic and enzymatic analysis of SABATH gene family from Larix kaempferi

The SABATH methyltransferases catalyze methylation of small-molecule metabolites, which participate in plant growth, development and defense response. Given lack of genome-wide studies on gymnosperms SABATH family, the formation and functional differentiation mechanism of the Larix kaempferi SABATH gene family was systematically and exhaustively explored by analyzing gene sequence characteristics, phylogenetic relationship, expression pattern, and enzyme activities. Phylogenetic analysis showed that 247 SABATH genes from 14 land plants were divided into 4 clades, and lineage-specific gene duplication events were important factors that contributed to the evolution of the SABATH gene family in gymnosperms and angiosperms. Substrate specificity analysis of 18 Larix SABATH proteins showed that LaSABATHs could catalyze O-methylation of indole-3-acetic acid (IAA) and farnesic acid (FA), N-methylation of theobromine, and S-methylation of thiobenzoic acid. Furthermore, only LaSABATH2 and LaSABATH29 could catalyze O-methylation of FA, and only LaSABATH30 could catalyze O-methylation of IAA. Homology modeling and molecular docking studies showed the hydrogen bond formed between the His188 of LaSABATH30 and IAA and the noticeable hydrophobic IAA-binding pocket may be helpful for IAA methylation. In this study, identification of proteins with significant specific catalytic activity toward FA and IAA provided high-quality candidate genes for forest genetics and breeding.

Transcriptome-wide modulation by Sargassum vulgare and Acanthophora spicifera extracts results in a prime-triggered plant signalling cascade in tomato and sweet pepper

Seaweed extracts (SWEs) are becoming integrated into crop production systems due to their multiple beneficial effects including growth promotion and induction of defence mechanisms. However, the comprehensive molecular mechanisms of these effects are yet to be elucidated. The current study investigated the transcriptomic changes induced by SWEs derived from Sargassum vulgare and Acanthophora spicifera on tomato and sweet pepper plants. Tomato and sweet pepper plants were subjected to foliar treatment with alkaline extracts prepared from the above seaweeds. Transcriptome changes in the plants were assessed 72 h after treatments using RNA sequencing. The treated plants were also analysed for defence enzyme activities, nutrient composition and phytohormonal profiles. The results showed the significant enrichment of genes associated with several growth and defence processes including photosynthesis, carbon and nitrogen metabolism, plant hormone signal transduction, plant-pathogen interaction, secondary metabolite metabolism, MAPK signalling and amino acid biosynthesis. Activities of defence enzymes were also significantly increased in SWE-treated plants. Plant nutrient profiling showed significant increases in calcium, potassium, nitrogen, sulphur, boron, copper, iron, manganese, zinc and phosphorous levels in SWE-treated plants. Furthermore, the levels of auxins, cytokinins and gibberellins were also significantly increased in the treated plants. The severity of bacterial leaf spot and early blight incidence in plants treated with SWE was significantly reduced, in addition to other effects like an increase in chlorophyll content, plant growth, and fruit yield. The results demonstrated the complex effect of S. vulgare and A. spicifera extracts on the plants' transcriptome and provided evidence of a strong role of these extracts in increasing plant growth responses while priming the plants against pathogenic attack simultaneously. The current study contributes to the understanding of the molecular mechanisms of SWEs in plants and helps their usage as a viable organic input for sustainable crop production.

Structural diversification of bioactive bibenzyls through modular co-culture leading to the discovery of a novel neuroprotective agent

Bibenzyls, a kind of important plant polyphenols, have attracted growing attention for their broad and remarkable pharmacological activities. However, due to the low abundance in nature, uncontrollable and environmentally unfriendly chemical synthesis processes, these compounds are not readily accessible. Herein, one high-yield bibenzyl backbone-producing Escherichia coli strain was constructed by using a highly active and substrate-promiscuous bibenzyl synthase identified from Dendrobium officinale in combination with starter and extender biosynthetic enzymes. Three types of efficiently post-modifying modular strains were engineered by employing methyltransferases, prenyltransferase, and glycosyltransferase with high activity and substrate tolerance together with their corresponding donor biosynthetic modules. Structurally different bibenzyl derivatives were tandemly and/or divergently synthesized by co-culture engineering in various combination modes. Especially, a prenylated bibenzyl derivative (12) was found to be an antioxidant that exhibited potent neuroprotective activity in the cellular and rat models of ischemia stroke. RNA-seq, quantitative RT-PCR, and Western-blot analysis demonstrated that 12 could up-regulate the expression level of an apoptosis-inducing factor, mitochondria associated 3 (Aifm3), suggesting that Aifm3 might be a new target in ischemic stroke therapy. This study provides a flexible plug-and-play strategy for the easy-to-implement synthesis of structurally diverse bibenzyls through a modular co-culture engineering pipeline for drug discovery.

Two O-Methyltransferases Mediate Multiple Methylation Steps in the Biosynthesis of Coumarins in Cnidium monnieri

Coumarins with methoxy groups such as osthole (1), xanthotoxin (2), bergapten (3), and isopimpinellin (4) are typical bioactive ingredients of many medicinal plants. The methylation steps remain widely unknown. Herein, we report the discovery of two methyltransferases in the biosynthesis of O-methyl coumarins in Cnidium monnieri by transcriptome mining, heterologous expression, and in vitro enzymatic assays. The results reveal that (i) CmOMT1 catalyzes the methylation of osthenol (8) as the final step in the biosynthesis of 1, (ii) CmOMT2 shows the highest efficiency and preference for methylating xanthotoxol (11) to form 2, and (iii) CmOMT1 and CmOMT2 also efficiently transform bergaptol (10) and 8-hydroxybergapten (13) into 3 or 4, suggesting the CmOMTs mediate multistep methylations in the biosynthesis of linear furanocoumarins in C. monnieri.

Infection by Pseudocercospora musae leads to an early reprogramming of the Musa paradisiaca defense transcriptome

Deep sequencing technologies such as RNA sequencing can help unravel mechanisms governing defense or resistance responses in plant-pathogen interactions. Several studies have been carried out to investigate the transcriptomic changes in Musa germplasm against Yellow Sigatoka disease, but the defense response of Musa paradisiaca has not been investigated so far. We carried out transcriptome sequencing of M. paradisiaca var. Kachkal infected with the pathogen Pseudocercospora musae and found that a vast set of genes were upregulated while many genes were downregulated in the resistant cultivar as a result of infection. After transcriptome assembly and differential gene expression analysis, 429 upregulated and 156 downregulated genes were filtered out (considering fold change ± 2, p < 0.01). Functional annotation of the differentially expressed genes (DEGs) enriched the upregulated genes into 49 gene ontology (GO) classes of biological processes (BP), 20 classes of molecular function (MF) and 9 classes of cellular component (CC). Similarly, the downregulated genes were classified into 35 GO classes of BP, 28 classes of MF and 6 classes of CC. The KEGG enrichment analysis revealed that the upregulated genes were most highly represented in 'metabolic' and 'biosynthesis of secondary metabolites' pathways. Additionally, 'plant hormone signal transduction', 'plant-pathogen interaction' and 'phenylpropanoid biosynthesis' pathways were also significantly enriched indicating their involvement in resistance responses against the pathogen. The RNA-seq analysis also depicts that a range of important defense-related genes are modulated as a result of infection, all of which are responsible for either mediating or activating resistance responses in the host. We studied and validated the expression profiles of ten important defense-related genes potentially involved in conferring resistance to the pathogen through qRT-PCR. Almost all the selected defense-related genes were found to be highly and significantly upregulated within 24 h post inoculation (hpi) and for some genes, the expression remained consistently high till the later time point of 72 hpi. These results, thus, indicate that the infection by P. musae leads to a rapid reprogramming of the defense transcriptome of the resistant banana cultivar. The defense-related genes identified to be modulated in response to infection are important not only for pathogen recognition and perception but also for activation and persistence of defense in the host.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03245-9.

Transcriptome Analysis of Soursop (Annona muricata L.) Fruit under Postharvest Storage Identifies Genes Families Involved in Ripening

Soursop (Annona muricata L.) is climacteric fruit with a short ripening period and postharvest shelf life, leading to a rapid softening. In this study, transcriptome analysis of soursop fruits was performed to identify key gene families involved in ripening under postharvest storage conditions (Day 0, Day 3 stored at 28 ± 2 °C, Day 6 at 28 ± 2 °C, Day 3 at 15 ± 2 °C, Day 6 at 15 ± 2 °C, Day 9 at 15 ± 2 °C). The transcriptome analysis showed 224,074 transcripts assembled clustering into 95, 832 unigenes, of which 21, 494 had ORF. RNA-seq analysis showed the highest number of differentially expressed genes on Day 9 at 15 ± 2 °C with 9291 genes (4772 up-regulated and 4519 down-regulated), recording the highest logarithmic fold change in pectin-related genes. Enrichment analysis presented significantly represented GO terms and KEGG pathways associated with molecular function, metabolic process, catalytic activity, biological process terms, as well as biosynthesis of secondary metabolites, plant hormone signal, starch, and sucrose metabolism, plant–pathogen interaction, plant–hormone signal transduction, and MAPK-signaling pathways, among others. Network analysis revealed that pectinesterase genes directly regulate the loss of firmness in fruits stored at 15 ± 2 °C.

A new approach to zip-lignin: 3,4-dihydroxybenzoate is compatible with lignification

Renewed interests in the development of bioenergy, biochemicals, and biomaterials have elicited new strategies for engineering the lignin of biomass feedstock plants. This study shows, for the first time, that 3,4-dihydroxybenzoate (DHB) is compatible with the radical coupling reactions that assemble polymeric lignin in plants. We introduced a bacterial 3-dehydroshikimate dehydratase into hybrid poplar (Populus alba × grandidentata) to divert carbon flux away from the shikimate pathway, which lies upstream of lignin biosynthesis. Transgenic poplar wood had up to 33% less lignin with p-hydroxyphenyl units comprising as much as 10% of the lignin. Mild alkaline hydrolysis of transgenic wood released fewer ester-linked p-hydroxybenzoate groups than control trees, and revealed the novel incorporation of cell-wall-bound DHB, as well as glycosides of 3,4-dihydroxybenzoic acid (DHBA). Two-dimensional nuclear magnetic resonance (2D-NMR) analysis uncovered DHBA-derived benzodioxane structures suggesting that DHB moieties were integrated into the lignin polymer backbone. In addition, up to 40% more glucose was released from transgenic wood following ionic liquid pretreatment and enzymatic hydrolysis. This work highlights the potential of diverting carbon flux from the shikimate pathway for lignin engineering and describes a new type of 'zip-lignin' derived from the incorporation of DHB into poplar lignin.

Comparative DNA Methylome of Phytoplasma Associated Retrograde Metamorphosis in Sesame (Sesamum indicum L.)

Phytoplasma-associated diseases such as phyllody and little leaf are critical threats to sesame cultivation worldwide. The mechanism of the dramatic conversion of flowers to leafy structures leading to yield losses and the drastic reduction in leaf size due to Phytoplasma infection remains yet to be identified. Cytosine methylation profiles of healthy and infected sesame plants studied using Whole Genome Bisulfite Sequencing (WGBS) and Quantitative analysis of DNA methylation with the real-time PCR (qAMP) technique revealed altered DNA methylation patterns upon infection. Phyllody was associated with global cytosine hypomethylation, though predominantly in the CHH (where H = A, T or C) context. Interestingly, comparable cytosine methylation levels were observed between healthy and little leaf-affected plant samples in CG, CHG and CHH contexts. Among the different genomic fractions, the highest number of differentially methylated Cytosines was found in the intergenic regions, followed by promoter, exonic and intronic regions in decreasing order. Further, most of the differentially methylated genes were hypomethylated and were mainly associated with development and defense-related processes. Loci for STOREKEEPER protein-like, a DNA-binding protein and PP2-B15, an F-Box protein, responsible for plugging sieve plates to maintain turgor pressure within the sieve tubes were found to be hypomethylated by WGBS, which was confirmed by methylation-dependent restriction digestion and qPCR. Likewise, serine/threonine-protein phosphatase-7 homolog, a positive regulator of cryptochrome signaling involved in hypocotyl and cotyledon growth and probable O-methyltransferase 3 locus were determined to be hypermethylated. Phytoplasma infection-associated global differential methylation as well as the defense and development-related loci reported here for the first time significantly elucidate the mechanism of phytoplasma-associated disease development.

Characterization and Structural Analysis of Emodin-O-Methyltransferase from Aspergillus terreus

All O-methylated derivatives of emodin, including physcion, questin, and 1-O-methylemodin, show potential antifungal activities. Notably, emodin and questin are two pivotal intermediates of geodin biosynthesis in Aspergillus terreus. Although most of the geodin biosynthetic steps have been investigated, the key O-methyltransferase (OMT) responsible for the O-methylation of emodin to generate questin has remained unidentified. Herein, through phylogenetic tree analysis and in vitro biochemical assays, the long-sought class II emodin-O-methyltransferase GedA has been functionally characterized. Additionally, the catalytic mechanism and key residues at the catalytic site of GedA were elucidated by enzyme-substrate-methyl donor analogue ternary complex crystal structure determination and site-directed mutagenesis. As we demonstrate, GedA adopts a typical general acid/base (E446/H373)-mediated transmethylation mechanism. In particular, residue D374 is also crucial for efficient catalysis through blocking the formation of intramolecular hydrogen bonds in emodin. This study will facilitate future engineering of GedA for the production of physcion or other site-specific O-methylated anthraquinone derivatives with potential applications as biopesticides.

Molecular Characterisation of Flavanone O-methylation in Eucalyptus

Flavonoids are ubiquitous polyphenolic compounds in plants, long recognised for their health-promoting properties in humans. Methylated flavonoids have received increasing attention due to the potential of methylation to enhance medicinal efficacy. Recently, Eucalyptus species with high levels of the O-methylated flavanone pinostrobin have been identified. Pinostrobin has potential commercial value due to its numerous pharmacological and functional food benefits. Little is known about the identity or mode of action of the enzymes involved in methylating flavanones. This study aimed to identify and characterise the methyltransferase(s) involved in the regiospecific methylation of pinostrobin in Eucalyptus and thereby add to our limited understanding of flavanone biosynthesis in plants. RNA-seq analysis of leaf tips enabled the isolation of a gene encoding a flavanone 7-O-methyltransferase (EnOMT1) in Eucalyptus. Biochemical characterisation of its in vitro activity revealed a range of substrates upon which EnOMT1 acts in a regiospecific manner. Comparison to a homologous sequence from a Eucalyptus species lacking O-methylated flavonoids identified critical catalytic amino acid residues within EnOMT1 responsible for its activity. This detailed molecular characterisation identified a methyltransferase responsible for chemical ornamentation of the core flavanone structure of pinocembrin and helps shed light on the mechanism of flavanone biosynthesis in Eucalyptus.

Diversification of Chemical Structures of Methoxylated Flavonoids and Genes Encoding Flavonoid-O-Methyltransferases

The O-methylation of specialized metabolites in plants is a unique decoration that provides structural and functional diversity of the metabolites with changes in chemical properties and intracellular localizations. The O-methylation of flavonoids, which is a class of plant specialized metabolites, promotes their antimicrobial activities and liposolubility. Flavonoid O-methyltransferases (FOMTs), which are responsible for the O-methylation process of the flavonoid aglycone, generally accept a broad range of substrates across flavones, flavonols and lignin precursors, with different substrate preferences. Therefore, the characterization of FOMTs with the physiology roles of methoxylated flavonoids is useful for crop improvement and metabolic engineering. In this review, we summarized the chemodiversity and physiology roles of methoxylated flavonoids, which were already reported, and we performed a cross-species comparison to illustrate an overview of diversification and conserved catalytic sites of the flavonoid O-methyltransferases.

Genomic insights into the evolution of Echinochloa species as weed and orphan crop

As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.

Evolution of DIMBOA-Glc O-Methyltransferases from Flavonoid O-Methyltransferases in the Grasses

O-Methylated benzoxazinoids (BXs) and flavonoids are widespread defenses against herbivores and pathogens in the grasses (Poaceae). Recently, two flavonoid O-methyltransferases (FOMTs), ZmFOMT2 and ZmFOMT3, have been reported to produce phytoalexins in maize (Zea mays). ZmFOMT2 and ZmFOMT3 are closely related to the BX O-methyltransferases (OMTs) ZmBX10-12 and ZmBX14, suggesting a common evolutionary origin in the Poaceae. Here, we studied the evolution and enzymatic requirements of flavonoid and BX O-methylation activities in more detail. Using BLAST searches and phylogenetic analyses, we identified enzymes homologous to ZmFOMT2 and ZmFOMT3, ZmBX10-12, and ZmBX14 in several grasses, with the most closely related candidates found almost exclusively in species of the Panicoideae subfamily. Biochemical characterization of candidate enzymes from sorghum (Sorghum bicolor), sugar cane (Saccharum spp.), and teosinte (Zea nicaraguensis) revealed either flavonoid 5-O-methylation activity or DIMBOA-Glc 4-O-methylation activity. However, DIMBOA-Glc 4-OMTs from maize and teosinte also accepted flavonols as substrates and converted them to 3-O-methylated derivatives, suggesting an evolutionary relationship between these two activities. Homology modeling, sequence comparisons, and site-directed mutagenesis led to the identification of active site residues crucial for FOMT and BX OMT activity. However, the full conversion of ZmFOMT2 activity into BX OMT activity by switching these residues was not successful. Only trace O-methylation of BXs was observed, indicating that amino acids outside the active site cavity are also involved in determining the different substrate specificities. Altogether, the results of our study suggest that BX OMTs have evolved from the ubiquitous FOMTs in the PACMAD clade of the grasses through a complex series of amino acid changes.

Systematic analysis reveals O-methyltransferase gene family members involved in flavonoid biosynthesis in grape

O-methyltransferases (OMTs) are an important group of enzymes involved in the methylation of various secondary metabolites, including flavonoids. However, the features and functions of OMTs have not been comprehensively studied in grape (Vitis vinifera), a rich source of methylated flavonoids. Here, 47 OMT members were identified in grape genome. They were unevenly distributed on grape chromosomes and some genes were tandem duplicated, indicating the role of duplication processes in the expansion of this gene family. Based on the phylogenetic relationship, these OMTs were clustered into CCoAOMT and COMT subclades, which were further supported by the results of conserved motif and gene structure analysis. Correlation analysis revealed that three members (VvCCoAOMT1, VvCCoAOMT4, and VvCOMT1) were potentially involved in the synthesis of most methylated flavonoids in the berry skins. Expression profiling based on RNA-seq data and qRT-PCR experiments indicated that VvCCoAOMT1 and VvCCoAOMT4 had specific and high expression in berry skins, and responded to abscisic acid and high temperature treatments; and that VvCOMT1 expression was significantly induced during berry development and UVC treatment. Cis-regulatory element analysis suggested important roles of OMTs in growth, development, and defense against stresses. We further demonstrated the transcriptional regulation of VvCCoAOMT4 by VvMYBA1, a master regulator of grape berry anthocyanin, and verified the protein localization of VvCCoAOMT4 in membrane and nucleus. These findings facilitate a better understanding of the characteristics of OMT gene family, especially of the potential members involved in the formation of O-methylated flavonoids in grape.

Molecular cloning and characterization of two distinct caffeoyl CoA O-methyltransferases (CCoAOMTs) from the liverwort Marchantia paleacea

Caffeoyl CoA O-methyltransferases (CCoAOMTs) catalyze the transfer of a methyl group from S-adenosylmethionine to a hydroxyl moiety of caffeoyl-CoA as part of the lignin biosynthetic pathway. CCoAOMT-like proteins also catalyze to a variety of flavonoids, coumarins, and phenylpropanoids. Several CCoAOMTs that prefer flavonoids as substrates have been characterized from liverworts. Here, we cloned two CCoAOMT genes, MpalOMT2 and MpalOMT3, from the liverwort Marchantia paleacea. MpalOMT3 has a second ATG codon downstream and the truncated version that lacks 11 amino acids was named MpalOMT3-Tr. Phylogenetic analysis placed MpalOMT3 at the root of the clade with true CCoAOMTs from vascular plants and placed MpalOMT2 between the CCoAOMT and CCoAOMT-like proteins. Recombinant OMTs methylated caffeoyl CoA, phenylpropanoids, and flavonoids containing two or three vicinal hydroxyl groups. MpalOMT3 showed higher catalytic activity for phenylpropanoids than MpalOMT2, but MpalOMT2 showed more promiscuous towards eriodictyol and myricetin. The lignin content in Arabidopsis thaliana stems increased with constitutive heterologous expression of MpalOMT3-Tr, but not MpalOMT2. Subcellular localization experiments indicated that the N-terminus of MpalOMT3 probably served as a chloroplast transit peptide and inhibited its enzymatic activity. Combining the phylogenetic analysis and functional characterization, we conclude that the liverwort M. paleacea harbors true CCoAOMT and CCoAOMT-like genes.

Genome-Wide Identification and Characterization of GhCOMT Gene Family during Fiber Development and Verticillium Wilt Resistance in Cotton

Caffeic acid O-methyltransferases (COMTs) play an essential role in lignin synthesis procession, especially in the plant’s phenylalanine metabolic pathway. The content of COMT genes in cotton and the relationship between their expression patterns have not been studied clearly in cotton. In this study, we have identified 190 COMT genes in cotton, which were classified into three groups (I, II and III), and mapped on the cotton chromosomes. In addition, we found that 135 of the 190 COMT genes result from dispersed duplication (DSD) and whole-genome duplication (WGD), indicating that DSD and WGD were the main forces driving COMT gene expansion. The Ka/Ks analysis showed that GhCOMT43 and GhCOMT41 evolved from GaCOMT27 and GrCOMT14 through positive selection. The results of qRT-PCR showed that GhCOMT13, GhCOMT28, GhCOMT39 and GhCOMT55 were related to lignin content during the cotton fiber development. GhCOMT28, GhCOMT39, GhCOMT55, GhCOMT56 and GhCOMT57 responded to Verticillium Wilt (VW) and maybe related to VW resistance through lignin synthesis. Conclusively, this study found that GhCOMTs were highly expressed in the secondary wall thickening stage and VW. These results provide a clue for studying the functions of GhCOMTs in the development of cotton fiber and VW resistance and could lay a foundation for breeding cotton cultivates with higher quantity and high resistance to VW.

Genome sequencing of turmeric provides evolutionary insights into its medicinal properties

Curcuma longa, or turmeric, is traditionally known for its immense medicinal properties and has diverse therapeutic applications. However, the absence of a reference genome sequence is a limiting factor in understanding the genomic basis of the origin of its medicinal properties. In this study, we present the draft genome sequence of C. longa, belonging to Zingiberaceae plant family, constructed using 10x Genomics linked reads and Oxford Nanopore long reads. For comprehensive gene set prediction and for insights into its gene expression, transcriptome sequencing of leaf tissue was also performed. The draft genome assembly had a size of 1.02 Gbp with ~70% repetitive sequences, and contained 50,401 coding gene sequences. The phylogenetic position of C. longa was resolved through a comprehensive genome-wide analysis including 16 other plant species. Using 5,388 orthogroups, the comparative evolutionary analysis performed across 17 species including C. longa revealed evolution in genes associated with secondary metabolism, plant phytohormones signaling, and various biotic and abiotic stress tolerance responses. These mechanisms are crucial for perennial and rhizomatous plants such as C. longa for defense and environmental stress tolerance via production of secondary metabolites, which are associated with the wide range of medicinal properties in C. longa.

The role of watermelon caffeic acid O-methyltransferase (ClCOMT1) in melatonin biosynthesis and abiotic stress tolerance

Melatonin is a pleiotropic signaling molecule that regulates plant growth and responses to various abiotic stresses. The last step of melatonin synthesis in plants can be catalyzed by caffeic acid O-methyltransferase (COMT), a multifunctional enzyme reported to have N-acetylserotonin O-methyltransferase (ASMT) activity; however, the ASMT activity of COMT has not yet been characterized in nonmodel plants such as watermelon (Citrullus lanatus). Here, a total of 16 putative O-methyltransferase (ClOMT) genes were identified in watermelon. Among them, ClOMT03 (Cla97C07G144540) was considered a potential COMT gene (renamed ClCOMT1) based on its high identities (60.00-74.93%) to known COMT genes involved in melatonin biosynthesis, expression in almost all tissues, and upregulation under abiotic stresses. The ClCOMT1 protein was localized in the cytoplasm. Overexpression of ClCOMT1 significantly increased melatonin contents, while ClCOMT1 knockout using the CRISPR/Cas-9 system decreased melatonin contents in watermelon calli. These results suggest that ClCOMT1 plays an essential role in melatonin biosynthesis in watermelon. In addition, ClCOMT1 expression in watermelon was upregulated by cold, drought, and salt stress, accompanied by increases in melatonin contents. Overexpression of ClCOMT1 enhanced transgenic Arabidopsis tolerance against such abiotic stresses, indicating that ClCOMT1 is a positive regulator of plant tolerance to abiotic stresses.

Melatonin synthesis genes N-acetylserotonin methyltransferases evolved into caffeic acid O-methyltransferases and both assisted in plant terrestrialization

Terrestrialization is one of the most momentous events in the history of plant life, which leads to the subsequent evolution of plant diversity. The transition species, in this process, had to acquire a range of adaptive mechanisms to cope with the harsh features of terrestrial environments compared to that of aquatic habitat. As an ancient antioxidant, a leading regulator of ROS signaling or homeostasis, and a presumed plant master regulator, melatonin likely assisted plants transition to land and their adaption to terrestrial ecosystems. N-acetylserotonin methyltransferases (ASMT) and caffeic acid O-methyltransferases (COMT), both in the O-methyltransferase (OMT) family, catalyze the core O-methylation reaction in melatonin biosynthesis. How these two enzymes with close relevance evolved in plant evolutionary history and whether they participated in plant terrestrialization remains unknown. Using combined phylogenetic evidence and protein structure analysis, it is revealed that COMT likely evolved from ASMT by gene duplication and subsequent divergence. Newly emergent COMT gained a significantly higher ASMT activity to produce greater amounts of melatonin for immobile plants to acclimate to the stressful land environments after evolving from the more environmentally-stable aquatic conditions. The COMT genes possess more conserved substrate-binding sites at the amino acid level and more open protein conformation compared to ASMT, and getting a new function to catalyze the lignin biosynthesis. This development directly contributed to the dominance of vascular plants among the Earth's flora and prompted plant colonization of land. Thus, ASMT, together with its descendant COMT, might play key roles in plant transition to land. The current study provides new insights into plant terrestrialization with gene duplication contributing to this process along with well-known horizontal gene transfer.

Secondary Metabolites of Pseudomonas aeruginosa LV Strain Decrease Asian Soybean Rust Severity in Experimentally Infected Plants

Asian Soybean Rust (ASR), a disease caused by Phakopsora pachyrhizi, causing yield losses up to 90%. The control is based on the fungicides which may generate resistant fungi. The activation of the plant defense system, should help on ASR control. In this study, secondary metabolites of Pseudomonas aeruginosa LV strain were applied on spore germination and the expression of defense genes in infected soybean plants. The F4A fraction and the pure metabolites were used. In vitro, 10 µg mL⁻¹ of F4A reduced spore germination by 54%, while 100 µg mL⁻¹ completely inhibited. Overexpression of phenylalanine ammonia lyase (PAL), O-methyltransferase (OMT) and pathogenesis related protein-2 (PR-2; glucanases) defense-related genes were detected 24 and 72 h after soybean sprouts were sprayed with an organocopper antimicrobial compound (OAC). Under greenhouse conditions, the best control was observed in plants treated with 60 µg mL⁻¹ of PCA, which reduced ASR severity and lesion frequency by 75% and 43%, respectively. Plants sprayed with 2 and 20 µg mL⁻¹ of F4A also decreased severity (41%) and lesion frequency (32%). The significant reduction in spore germination ASR in plant suggested that the strain of these metabolites are effective against P. pachyrhizi, and they can be used for ASR control.

BiP-overexpressing soybean plants display accelerated hypersensitivity response (HR) affecting the SA-dependent sphingolipid and flavonoid pathways

Biotic and abiotic environmental stresses have limited the increase in soybean productivity. Overexpression of the molecular chaperone BiP in transgenic plants has been associated with the response to osmotic stress and drought tolerance by maintaining cellular homeostasis and delaying hypersensitive cell death. Here, we evaluated the metabolic changes in response to the hypersensitivity response (HR) caused by the non-compatible bacteria Pseudomonas syringae pv. tomato in BiP-overexpressing plants. The HR-modified metabolic profiles in BiP-overexpressing plants were significantly distinct from the wild-type untransformed. The transgenic plants displayed a lower abundance of HR-responsive metabolites as amino acids, sugars, carboxylic acids and signal molecules, including p-aminobenzoic acid (PABA) and dihydrosphingosine (DHS), when compared to infected wild-type plants. In contrast, salicylic acid (SA) biosynthetic and signaling pathways were more stimulated in transgenic plants, and both pathogenesis-related genes (PRs) and transcriptional factors controlling the SA pathway were more induced in the BiP-overexpressing lines. Furthermore, the long-chain bases (LCBs) and ceramide biosynthetic pathways showed alterations in gene expression and metabolite abundance. Thus, as a protective pathway against pathogens, HR regulation by sphingolipids and SA may account at least in part by the enhanced resistance of transgenic plants. GmNAC32 transcriptional factor was more induced in the transgenic plants and it has also been reported to regulate flavonoid synthesis in response to SA. In fact, the BiP-overexpressing plants showed an increase in flavonoids, mainly prenylated isoflavones, as precursors for phytoalexins. Our results indicate that the BiP-mediated acceleration in the hypersensitive response may be a target for metabolic engineering of plant resistance against pathogens.

Isoflavone malonyl-CoA acyltransferase GmMaT2 is involved in nodulation of soybean by modifying synthesis and secretion of isoflavones

Malonyl-CoA:flavonoid acyltransferases (MaTs) modify isoflavones, but only a few have been characterized for activity and assigned to specific physiological processes. Legume roots exude isoflavone malonates into the rhizosphere, where they are hydrolyzed into isoflavone aglycones. Soybean GmMaT2 was highly expressed in seeds, root hairs, and nodules. GmMaT2 and GmMaT4 recombinant enzymes used isoflavone 7-O-glucosides as acceptors and malonyl-CoA as an acyl donor to generate isoflavone glucoside malonates. GmMaT2 had higher activity towards isoflavone glucosides than GmMaT4. Overexpression in hairy roots of GmMaT2 and GmMaT4 produced more malonyldaidzin, malonylgenistin, and malonylglycitin, and resulted in more nodules than control. However, only GmMaT2 knockdown (KD) hairy roots showed reduced levels of malonyldaidzin, malonylgenistin, and malonylglycitin, and, likewise, reduced nodule numbers. These were consistent with the up-regulation of only GmMaT2 by rhizobial infection, and higher expression levels of early nodulation genes in GmMaT2- and GmMaT4-overexpressing roots, but lower only in GmMaT2-KD roots compared with control roots. Higher malonyl isoflavonoid levels in transgenic hairy roots were associated with higher levels of isoflavones in root exudates and more nodules, and vice versa. We suggest that GmMaT2 participates in soybean nodulation by catalyzing isoflavone malonylation and affecting malonyl isoflavone secretion for activation of Nod factor and nodulation.

Identification and expression analysis of caffeoyl-coenzyme A O-methyltransferase family genes related to lignin biosynthesis in tea plant (Camellia sinensis)

Tea plant, an economically important crop, is used in producing tea, which is a non-alcoholic beverage. Lignin, the second most abundant component of the cell wall, reduces the tenderness of tea leaves and affects tea quality. Caffeoyl-coenzyme A O-methyltransferase (CCoAOMT) involved in lignin biosynthesis affects the efficiency of lignin synthesis and lignin composition. A total of 10 CsCCoAOMTs were identified based on tea plant genome. Systematic analysis of CCoAOMTs was conducted for its physicochemical properties, phylogenetic relationships, conserved motifs, gene structure, and promoter cis-element prediction. Phylogenetic analysis suggested that all the CsCCoAOMT proteins can be categorized into three clades. The promoters of six CsCCoAOMT genes possessed lignin-specific cis-elements, indicating they are possibly essential for lignin biosynthesis. According to the distinct tempo-spatial expression profiles, five genes were substantially expressed in eight tested tissues. Most CsCCoAOMT genes were expressed in stems and leaves in three tea plant cultivars 'Longjing 43,' 'Anjibaicha,' and 'Fudingdabai' by RT-qPCR detection and analysis. The expression levels of two genes (CsCCoAOMT5 and CsCCoAOMT6) were higher than those of the other genes. The expression levels of most CsCCoAOMT genes in 'Longjing 43' were significantly higher than that those in 'Anjibaicha' and 'Fudingdabai.' Correlation analysis revealed that only the expression levels of CsCCoAOMT6 were positively correlated with lignin content in the leaves and stems. These results lay a foundation for the future exploration of the roles of CsCCoAOMTs in lignin biosynthesis in tea plant.

Biochemical Characterization of a Flavonoid O-methyltransferase from Perilla Leaves and Its Application in 7-Methoxyflavonoid Production

Methylation is a common structural modification that can alter and improve the biological activities of natural compounds. O-Methyltransferases (OMTs) catalyze the methylation of a wide array of secondary metabolites, including flavonoids, and are potentially useful tools for the biotechnological production of valuable natural products. An OMT gene (PfOMT3) was isolated from perilla leaves as a putative flavonoid OMT (FOMT). Phylogenetic analysis and sequence comparisons showed that PfOMT3 is a class II OMT. Recombinant PfOMT3 catalyzed the methylation of flavonoid substrates, whereas no methylated product was detected in PfOMT3 reactions with phenylpropanoid substrates. Structural analyses of the methylation products revealed that PfOMT3 regiospecifically transfers a methyl group to the 7-OH of flavonoids. These results indicate that PfOMT3 is an FOMT that catalyzes the 7-O-methylation of flavonoids. PfOMT3 methylated diverse flavonoids regardless of their backbone structure. Chrysin, naringenin and apigenin were found to be the preferred substrates of PfOMT3. Recombinant PfOMT3 showed moderate OMT activity toward eriodictyol, luteolin and kaempferol. To assess the biotechnological potential of PfOMT3, the biotransformation of flavonoids was performed using PfOMT3-transformed Escherichia coli. Naringenin and kaempferol were successfully bioconverted to the 7-methylated products sakuranetin and rhamnocitrin, respectively, by E. coli harboring PfOMT3.

A thorough analysis and categorization of bacterial interrupted adenylation domains, including previously unidentified families

Interrupted adenylation (A) domains are key to the immense structural diversity seen in the nonribosomal peptide (NRP) class of natural products (NPs). Interrupted A domains are A domains that contain within them the catalytic portion of another domain, most commonly a methylation (M) domain. It has been well documented that methylation events occur with extreme specificity on either the backbone (N-) or side chain (O- or S-) of the amino acid (or amino acid-like) building blocks of NRPs. Here, through taxonomic and phylogenetic analyses as well as multiple sequence alignments, we evaluated the similarities and differences between interrupted A domains. We probed their taxonomic distribution amongst bacterial organisms, their evolutionary relatedness, and described conserved motifs of each type of M domain found to be embedded in interrupted A domains. Additionally, we categorized interrupted A domains and the M domains within them into a total of seven distinct families and six different types, respectively. The families of interrupted A domains include two new families, 6 and 7, that possess new architectures. Rather than being interrupted between the previously described a2–a3 or a8–a9 of the ten conserved A domain sequence motifs (a1–a10), family 6 contains an M domain between a6–a7, a previously unknown interruption site. Family 7 demonstrates that di-interrupted A domains exist in Nature, containing an M domain between a2–a3 as well as one between a6–a7, displaying a novel arrangement. These in-depth investigations of amino acid sequences deposited in the NCBI database highlighted the prevalence of interrupted A domains in bacterial organisms, with each family of interrupted A domains having a different taxonomic distribution. They also emphasized the importance of utilizing a broad range of bacteria for NP discovery. Categorization of the families of interrupted A domains and types of M domains allowed for a better understanding of the trends of naturally occurring interrupted A domains, which illuminated patterns and insights on how to harness them for future engineering studies.

In-depth study of intriguing bacterial interrupted adenylation domains from seven distinct families and six different types.

Characterization of O-methyltransferases involved in the biosynthesis of tetrandrine in Stephania tetrandra

Tetrandrine is the most effective small molecule that has been found to inhibit the Ebola virus. It is a typical bisbenzylisoquinoline alkaloid and is the main active ingredient in Stephania tetrandra. Metabolic engineering and synthetic biology are potential methods for efficient and rapid acquisition of tetrandrine. S-adenosyl-L-methionine: (S)-norcoclaurine-6-O-methyltransferase (6OMT) is a rate-limiting step involved in the biosynthesis of tetrandrine. In this study, we identify S-adenosyl-L-methionine: (S)-norcoclaurine-6-O-methyltransferase from S. tetrandra, which catalyzes the conversion of (S)-norcoclaurine to (S)-coclaurine. Four 6OMT-like genes were cloned from S. tetrandra. An in vitro enzyme assay showed that St6OMT1 could catalyze the conversion of (S)-norcoclaurine to produce (S)-coclaurine. St6OMT2 can catalyze the production of very few (S)-coclaurine molecules, accompanied by more by-products with m/z 300, compared to St6OMT1. The newly discovered 6OMTs will provide an optional genetic component for benzylisoquinoline alkaloid (BIA) synthetic biology research. This work will lay the foundation for the analysis of the biosynthetic pathway of tetrandrine in S. tetrandra.

Biotransformation Changes Bioaccumulation and Toxicity of Diclofenac in Aquatic Organisms

Biotransformation plays a crucial role in regulating the bioaccumulation potential and toxicity of organic compounds in organisms but is, in general, poorly understood for emerging contaminants. Here, we have used diclofenac as a model compound to study the impact of biotransformation on the bioaccumulation potential and toxicity in two keystone aquatic invertebrates: Gammarus pulex and Hyalella azteca. In both species, diclofenac was transformed into several oxidation products and conjugates, including two novel products, that is, diclofenac taurine conjugate (DCF-M403) and unexpected diclofenac methyl ester (DCF-M310.03). The ratios of biotransformation products to parent compound were 12-17 for DCF-M403 and 0.01-0.7 for DCF-M310.03 after 24 h exposure. Bioconcentration factors (BCFs) of diclofenac were 0.5 and 3.2 L kg_ww^-1 in H. azteca and G. pulex, respectively, whereas BCFs of DCF-M310.03 was 164.5 and 104.7 L kg_ww^-1, respectively, representing a 25- to 110-fold increase. Acute toxicity of DCF-M310.03 was also higher than the parent compound in both species, which correlated well with the increased bioconcentration potential. The LC₅₀ of diclofenac in H. azteca was 216 mg L^-1, while that of metabolite DCF-M310.03 was reduced to only 0.53 mg L^-1, representing a 430-fold increase in acute toxicity compared to diclofenac. DCF-M403 is less toxic than its parent compound toward H. azteca, which may be linked to its slightly lower hydrophobicity. Furthermore, the transformation of diclofenac to its methyl ester derivative was explored in crude invertebrate extracts spiked with an S-adenosylmethionine cofactor, revealing possible catalysis by an S-adenosylmethionine-dependent carboxylic acid methyltransferase. Methylation of diclofenac was further detected in fish hepatocytes and human urine, indicating a broader relevance. Therefore, potentially methylated metabolites of polar contaminants should be considered for a comprehensive risk assessment in the future.

Expression Atlas of Selaginella moellendorffii Provides Insights into the Evolution of Vasculature, Secondary Metabolism, and Roots

Selaginella moellendorffii is a representative of the lycophyte lineage that is studied to understand the evolution of land plant traits such as the vasculature, leaves, stems, roots, and secondary metabolism. However, only a few studies have investigated the expression and transcriptional coordination of Selaginella genes, precluding us from understanding the evolution of the transcriptional programs behind these traits. We present a gene expression atlas comprising all major organs, tissue types, and the diurnal gene expression profiles for S. moellendorffii We show that the transcriptional gene module responsible for the biosynthesis of lignocellulose evolved in the ancestor of vascular plants and pinpoint the duplication and subfunctionalization events that generated multiple gene modules involved in the biosynthesis of various cell wall types. We demonstrate how secondary metabolism is transcriptionally coordinated and integrated with other cellular pathways. Finally, we identify root-specific genes and show that the evolution of roots did not coincide with an increased appearance of gene families, suggesting that the development of new organs does not coincide with increased fixation of new gene functions. Our updated database at conekt.plant.tools represents a valuable resource for studying the evolution of genes, gene families, transcriptomes, and functional gene modules in the Archaeplastida kingdom.

Transcriptomic and metabolomic analysis provides insights into anthocyanin and procyanidin accumulation in pear

BACKGROUND

Pear is one of the most important fruit crops worldwide. Anthocyanins and procyanidins (PAs) are important secondary metabolites that affect the appearance and nutritive quality of pear. However, few studies have focused on the molecular mechanism underlying anthocyanin and PA accumulation in pear.

RESULTS

We conducted metabolome and transcriptome analyses to identify candidate genes involved in anthocyanin and PA accumulation in young fruits of the pear cultivar 'Clapp Favorite' (CF) and its red mutation cultivar 'Red Clapp Favorite' (RCF). Gene-metabolite correlation analyses revealed a 'core set' of 20 genes that were strongly correlated with 10 anthocyanin and seven PA metabolites. Of these, PcGSTF12 was confirmed to be involved in anthocyanin and PA accumulation by complementation of the tt19-7 Arabidopsis mutant. Interestingly, PcGSTF12 was found to be responsible for the accumulation of procyanidin A3, but not petunidin 3, 5-diglucoside, opposite to the function of AtGSTs in Arabidopsis. Transformation with PcGSTF12 greatly promoted or repressed genes involved in anthocyanin and PA biosynthesis, regulation, and transport. Electrophoretic mobility shift and luciferase reporter assays confirmed positive regulation of PcGSTF12 by PcMYB114.

CONCLUSION

These findings identify a core set of genes for anthocyanin and PA accumulation in pear. Of these, PcGSTF12, was confirmed to be involved in anthocyanin and PA accumulation. Our results also identified an important anthocyanin and PA regulation node comprising two core genes, PcGSTF12 and PcMYB114. These results provide novel insights into anthocyanin and PA accumulation in pear and represent a valuable data set to guide future functional studies and pear breeding.

Isolation and characterization of two O-methyltransferases involved in benzylisoquinoline alkaloid biosynthesis in sacred lotus (Nelumbo nucifera)

Benzylisoquinoline alkaloids (BIAs) are a major class of plant metabolites with many pharmacological benefits. Sacred lotus (Nelumbo nucifera) is an ancient aquatic plant of medicinal value because of antiviral and immunomodulatory activities linked to its constituent BIAs. Although more than 30 BIAs belonging to the 1-benzylisoquinoline, aporphine, and bisbenzylisoquinoline structural subclasses and displaying a predominant R-enantiomeric conformation have been isolated from N. nucifera, its BIA biosynthetic genes and enzymes remain unknown. Herein, we report the isolation and biochemical characterization of two O-methyltransferases (OMTs) involved in BIA biosynthesis in sacred lotus. Five homologous genes, designated NnOMT1-5 and encoding polypeptides sharing >40% amino acid sequence identity, were expressed in Escherichia coli Functional characterization of the purified recombinant proteins revealed that NnOMT1 is a regiospecific 1-benzylisoquinoline 6-O-methyltransferase (6OMT) accepting both R- and S-substrates, whereas NnOMT5 is mainly a 7-O-methyltransferase (7OMT), with relatively minor 6OMT activity and a strong stereospecific preference for S-enantiomers. Available aporphines were not accepted as substrates by either enzyme, suggesting that O-methylation precedes BIA formation from 1-benzylisoquinoline intermediates. Km values for NnOMT1 and NnOMT5 were 20 and 13 μm for (R,S)-norcoclaurine and (S)-N-methylcoclaurine, respectively, similar to those for OMTs from other BIA-producing plants. Organ-based correlations of alkaloid content, OMT activity in crude extracts, and OMT gene expression supported physiological roles for NnOMT1 and NnOMT5 in BIA metabolism, occurring primarily in young leaves and embryos of sacred lotus. In summary, our work identifies two OMTs involved in BIA metabolism in the medicinal plant N. nucifera.