Re-Engineering Fungal Nonribosomal Peptide Synthetases by Module Dissection and Duplicated Thiolation Domains.

Unnatural product (uNP) nonribosomal peptides promise to be a valuable source of pharmacophores for drug discovery. However, the extremely large size and complexity of the nonribosomal peptide synthetase (NRPS) enzymes pose formidable challenges to the production of such uNPs by combinatorial biosynthesis and synthetic biology. Here we report a new NRPS dissection strategy that facilitates the engineering and heterologous production of these NRPSs. This strategy divides NRPSs into "splitting units", each forming an enzyme subunit that contains catalytically independent modules. Functional collaboration between the subunits is then facilitated by artificially duplicating, at the N-terminus of the downstream subunit, the linker - thiolation domain - linker fragment that is resident at the C-terminus of the upstream subunit. Using the suggested split site that follows a conserved motif in the linker connecting the adenylation and the thiolation domains allows cognate or chimeric splitting unit pairs to achieve productivities that match, and in many cases surpass those of hybrid chimeric enzymes, and even those of intact NRPSs, upon production in a heterologous chassis. Our strategy provides facile options for the rational engineering of fungal NRPSs and for the combinatorial reprogramming of nonribosomal peptide production.

A DNA demethylase reduces seed size by decreasing the DNA methylation of AT-rich transposable elements in soybean.

Understanding how to increase soybean yield is crucial for global food security. The genetic and epigenetic factors influencing seed size, a major crop yield determinant, are not fully understood. We explore the role of DNA demethylase GmDMEa in soybean seed size. Our research indicates that GmDMEa negatively correlates with soybean seed size. Using CRISPR-Cas9, we edited GmDMEa in the Dongnong soybean cultivar, known for small seeds. Modified plants had larger seeds and greater yields without altering plant architecture or seed nutrition. GmDMEa preferentially demethylates AT-rich transposable elements, thus activating genes and transcription factors associated with the abscisic acid pathway, which typically decreases seed size. Chromosomal substitution lines confirm that these modifications are inheritable, suggesting a stable epigenetic method to boost seed size in future breeding. Our findings provide insights into epigenetic seed size control and suggest a strategy for improving crop yields through the epigenetic regulation of crucial genes. This work implies that targeted epigenetic modification has practical agricultural applications, potentially enhancing food production without compromising crop quality.

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana.

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

Onygenaleosides A-F, 6/5 Bicyclic Ring Skeleton Triterpene Glycosides with Insecticidal Activity from Onygenales sp. YX1425.

Six new squalene derived polyether glycosides, onygenaleosides A-F (1-6), that possess a 6/5 bicyclic fused ring skeleton were isolated from the cultures of Onygenales sp. YX1425, along with two known analogues (7 and 8). The planar structures of the new compounds were elucidated based on analysis of NMR and MS spectroscopy data, and the absolute configuration of 1 was determined by the advanced Mosher method and quantum chemical calculations. Compound 2 was active against Spodoptera frugiperda with an LC50 value of 193.4 ± 1.1 µg/mL.

Heritable epigenetic modification of BpPIN1 is associated with leaf shapes in Betula pendula.

The new variety Betula pendula 'Dalecarlica', selected from Betula pendula, shows high ornamental value owing to its lobed leaf shape. In this study, to identify the genetic components of leaf shape formation, we performed bulked segregant analysis and molecular marker-based fine mapping to identify the causal gene responsible for lobed leaves in B. pendula 'Dalecarlica'. The most significant variations associated with leaf shape were identified within the gene BpPIN1 encoding a member of the PIN-FORMED family, responsible for the auxin efflux carrier. We further confirmed the hypomethylation at the promoter region promoting the expression level of BpPIN1, which causes stronger and longer veins and lobed leaf shape in B. pendula 'Dalecarlica'. These results indicated that DNA methylation at the BpPIN1 promoter region is associated with leaf shapes in B. pendula. Our findings revealed an epigenetic mechanism of BpPIN1 in the regulation of leaf shape in Betula  Linn. (birch), which could help in the molecular breeding of ornamental traits.

Climate-responsive DNA methylation is involved in the biosynthesis of lignin in birch.

Lignin is one of the most important secondary metabolites and essential to the formation of cell walls. Changes in lignin biosynthesis have been reported to be associated with environmental variations and can influence plant fitness and their adaptation to abiotic stresses. However, the molecular mechanisms underlying this association remain unclear. In this study, we evaluated the relations between the lignin biosynthesis and environmental factors and explored the role of epigenetic modification (DNA methylation) in contributing to these relations if any in natural birch. Significantly negative correlations were observed between the lignin content and temperature ranges. Analyzing the transcriptomes of birches in two habitats with different temperature ranges showed that the expressions of genes and transcription factors (TFs) involving lignin biosynthesis were significantly reduced at higher temperature ranges. Whole-genome bisulfite sequencing revealed that promoter DNA methylation of two NAC-domain TFs, BpNST1/2 and BpSND1, may be involved in the inhibition of these gene expressions, and thereby reduced the content of lignin. Based on these results we proposed a DNA methylation-mediated lignin biosynthesis model which responds to environmental factors. Overall, this study suggests the possibility of environmental signals to induce epigenetic variations that result in changes in lignin content, which can aid to develop resilient plants to combat ongoing climate changes or to manipulate secondary metabolite biosynthesis for agricultural, medicinal, or industrial values.

CRISPR/Cas9-Mediated Multiple Knockouts in Abscisic Acid Receptor Genes Reduced the Sensitivity to ABA during Soybean Seed Germination.

Abscisic acid (ABA) is an important plant hormone that regulates numerous functions in plant growth, development, and stress responses. Several proteins regulate the ABA signal transduction mechanism in response to environmental stress. Among them, the PYR1/PYL/RCAR family act as ABA receptors. This study used the CRISPR/Cas9 gene-editing system with a single gRNA to knock out three soybean PYL genes: GmPYL17, GmPYL18, and GmPYL19. The gRNA may efficiently cause varying degrees of deletion of GmPYL17, GmPYL18, and GmPYL19 gene target sequences, according to the genotyping results of T0 plants. A subset of induced alleles was successfully transferred to progeny. In the T2 generation, we obtained double and triple mutant genotypes. At the seed germination stage, CRISPR/Cas9-created GmPYL gene knockout mutants, particularly gmpyl17/19 double mutants, are less susceptible to ABA than the wild type. RNA-Seq was used to investigate the differentially expressed genes related to the ABA response from germinated seedlings under diverse treatments using three biological replicates. The gmpyl17/19-1 double mutant was less susceptible to ABA during seed germination, and mutant plant height and branch number were higher than the wild type. Under ABA stress, the GO enrichment analysis showed that certain positive germination regulators were activated, which reduced ABA sensitivity and enhanced seed germination. This research gives a theoretical basis for a better understanding of the ABA signaling pathway and the participation of the key component at their molecular level, which helps enhance soybean abiotic stress tolerance. Furthermore, this research will aid breeders in regulating and improving soybean production and quality under various stress conditions.

Transgenerationally Transmitted DNA Demethylation of a Spontaneous Epialleles Using CRISPR/dCas9-TET1cd Targeted Epigenetic Editing in Arabidopsis.

CRISPR/dCas9 is an important DNA modification tool in which a disarmed Cas9 protein with no nuclease activity is fused with a specific DNA modifying enzyme. A previous study reported that overexpression of the TET1 catalytic domain (TET1cd) reduces genome-wide methylation in Arabidopsis. A spontaneous naturally occurring methylation region (NMR19-4) was identified in the promoter region of the PPH (Pheophytin Pheophorbide Hydrolase) gene, which encodes an enzyme that can degrade chlorophyll and accelerate leaf senescence. The methylation status of NMR19-4 is associated with PPH expression and leaf senescence in Arabidopsis natural accessions. In this study, we show that the CRISPR/dCas9-TET1cd system can be used to target the methylation of hypermethylated NMR19-4 region to reduce the level of methylation, thereby increasing the expression of PPH and accelerating leaf senescence. Furthermore, hybridization between transgenic demethylated plants and hypermethylated ecotypes showed that the demethylation status of edited NMR19-4, along with the enhanced PPH expression and accelerated leaf senescence, showed Mendelian inheritance in F1 and F2 progeny, indicating that spontaneous epialleles are stably transmitted trans-generationally after demethylation editing. Our results provide a rational approach for future editing of spontaneously mutated epialleles and provide insights into the epigenetic mechanisms that control plant leaf senescence.

Genome-Wide Identification, Expression and Interaction Analysis of ARF and AUX/IAA Gene Family in Soybean.

BACKGROUND: The plant hormones auxin affects most aspects of plant growth and development. The auxin transport and signaling are regulated by different factors that modulate plant morphogenesis and respond to external environments. The modulation of gene expression by Auxin Response Factors (ARFs) and inhibiting Auxin/Indole-3-Acetic Acid (Aux/IAA) proteins are involved in auxin signaling pathways. These components are encoded by gene families with numerous members in most flowering plants. METHODS: However, there is no genome-wide analysis of the expression profile and the structural and functional properties of the ARF and Aux/IAA gene families in soybean. Using various online tools to acquire of genomic and expression data, and analyzing them to differentiate the selected gene family's expression, interaction, and responses in plant growth and development. RESULTS: Here, we discovered 63 GmIAAs and 51 GmARFs in a genome-wide search for soybean and analyzed the genomic, sequential and structural properties of GmARFs and GmIAAs. All of the GmARFs found have the signature B3 DNA-binding (B3) and ARF (Aux rep) domains, with only 23 possessing the C-terminal PB1 (Phox and Bem1) domain (Aux/IAA), according to domain analysis. The number of exons in GmARFs and GmIAAs genes varies from two to sixteen, indicating that the gene structure of GmARFs and GmIAAs is highly variable. Based on phylogenetic analysis, the 51 GmARFs and 63 GmIAAs were classified into I-V and I-VII groups. The expression pattern of GmARFs and GmIAAs revealed that the GmARF expression is more specific to a particular part of the plant; for example, ARF 2, 7, and 11 are highly expressed in the root. In contrast, GmIAAs expression has occurred in various parts of the plants. The interaction of ARF with functional genes showed extensive interactions with genes involved in auxin transport which helps to control plant growth and development. Furthermore, we also elaborate on the DNA-protein interaction of ARFs by identifying the residues involved in DNA recognition. CONCLUSIONS: This study will improve our understanding of the auxin signaling system and its regulatory role in plant growth and development.

Chemometrics and genome mining reveal an unprecedented family of sugar acid-containing fungal nonribosomal cyclodepsipeptides.

Xylomyrocins, a unique group of nonribosomal peptide secondary metabolites, were discovered in Paramyrothecium and Colletotrichum spp. fungi by employing a combination of high-resolution tandem mass spectrometry (HRMS/MS)-based chemometrics, comparative genome mining, gene disruption, stable isotope feeding, and chemical complementation techniques. These polyol cyclodepsipeptides all feature an unprecedented d-xylonic acid moiety as part of their macrocyclic scaffold. This biosynthon is derived from d-xylose supplied by xylooligosaccharide catabolic enzymes encoded in the xylomyrocin biosynthetic gene cluster, revealing a novel link between carbohydrate catabolism and nonribosomal peptide biosynthesis. Xylomyrocins from different fungal isolates differ in the number and nature of their amino acid building blocks that are nevertheless incorporated by orthologous nonribosomal peptide synthetase (NRPS) enzymes. Another source of structural diversity is the variable choice of the nucleophile for intramolecular macrocyclic ester formation during xylomyrocin chain termination. This nucleophile is selected from the multiple available alcohol functionalities of the polyol moiety, revealing a surprising polyspecificity for the NRPS terminal condensation domain. Some xylomyrocin congeners also feature N-methylated amino acid residues in positions where the corresponding NRPS modules lack N-methyltransferase (M) domains, providing a rare example of promiscuous methylation in the context of an NRPS with an otherwise canonical, collinear biosynthetic program.

A Putative Plasma Membrane Na+/H+ Antiporter GmSOS1 Is Critical for Salt Stress Tolerance in Glycine max.

Soybean (Glycine max) is a staple crop and a major source of vegetable protein and vegetable oil. The growth of soybean is dramatically inhibited by salt stress, especially by the excessive toxic Na+. Salt Overly Sensitive 1 (SOS1) is the only extensively characterized Na+ efflux transporter in multiple plant species so far. However, the role of GmSOS1 in soybean salt stress responses remains unclear. Herein, we created three gmsos1 mutants using the CRISPR-Cas9 system in soybean. We found a significant accumulation of Na+ in the roots of the gmsos1 mutants, resulting in the imbalance of Na+ and K+, which links to impaired Na+ efflux and increased K+ efflux in the roots of the gmsos1 mutants under salt stress. Compared to the wild type, our RNA-seq analysis revealed that the roots of the gmsos1-1 showed preferential up and downregulation of ion transporters under salt stress, supporting impaired stress detection or an inability to develop a comprehensive response to salinity in the gmsos1 mutants. Our findings indicate that the plasma membrane Na+/H+ exchanger GmSOS1 plays a critical role in soybean salt tolerance by maintaining Na+ homeostasis and provides evidence for molecular breeding to improve salt tolerance in soybean and other crops.

Genome-Wide Identification and Characterization of Caffeic Acid O-Methyltransferase Gene Family in Soybean.

Soybean is one of the most important legumes, providing high-quality protein for humans. The caffeic acid O-methyltransferase (COMT) gene has previously been demonstrated to be a critical gene that regulates lignin production in plant cell walls and plays an important function in plant growth and development. However, the COMT gene family has not been studied in soybeans. In this study, 55 COMT family genes in soybean were identified by phylogenetic analysis and divided into two groups, I and II. The analysis of conserved domains showed that all GmCOMTs genes contained Methyltransferase-2 domains. Further prediction of cis-acting elements showed that GmCOMTs genes were associated with growth, light, stress, and hormonal responses. Eventually, based on the genomic data of soybean under different stresses, the results showed that the expression of GmCOMTs genes was different under different stresses, such as salt and drought stress. This study has identified and characterized the COMT gene family in soybean, which provides an important theoretical basis for further research on the biological functions of COMT genes and promotes revealing the role of GmCOMTs genes under stress resistance.

Loss of chromatin remodeler DDM1 causes segregation distortion in Arabidopsis thaliana.

MAIN CONCLUSION: In ddm1 mutants, the DNA methylation is primarily affected in the heterochromatic region of the chromosomes, which is associated with the segregation distortion of SNPs in the F2 progenies. Segregation distortion (SD) is common in most genetic mapping experiments and a valuable resource to determine how gene loci induce deviation. Meiotic DNA crossing over and SD are under the control of several types of epigenetic modifications. DNA methylation is an important regulatory epigenetic modification that is inherited across generations. In the present study, we investigated the relationship between SD and DNA methylation. The ecotypes Col-0/C24 and chromatin remodeler mutants ddm1-10/Col and ddm1-15/C24 were reciprocally crossed to obtain F2 generations. A total of 300 plants for each reciprocally crossed plant in the F2 generations were subjected to next-generation sequencing to detect the single-nucleotide polymorphisms (SNPs) as DNA markers. All SNPs were analyzed using the Chi-square test method to determine their segregation ratio in F2 generations. Through the segregation ratio, whole-genome SNPs were classified into 16 classes. In class 10, the SNPs in the reciprocal crosses of wild type showed the expected Mendelian ratio of 1:2:1, while those in the reciprocal crosses of ddm1 mutants showed distortion. In contrast, all SNPs in class 16 displayed a normal 1:2:1 ratio, and class 1 showed SD, regardless of wild type or mutants, as assessed using CAPS (cleaved amplified polymorphic sequences) marker analysis to confirm the next-generation sequencing. In ddm1 mutants, the DNA methylation is highly reduced throughout the whole genome and more significantly in the heterochromatic regions of chromosomes. Our results showed that the ddm1 mutants exhibit low levels of DNA methylation, which facilitates the SD of SNPs primarily located in the heterochromatic region of chromosomes by reducing the heterozygous ratio. The present study will provide a strong base for future research focusing on the impact of DNA methylation on trait segregation and plant evolution.

Role of epigenetic modifications in the development of crops essential traits.

Epigenetic modification refers to the chemical modifications of chromosomal DNA and histones, mainly including DNA methylation, histone modifications and non-coding RNAs. Without altering the DNA sequence, these heritable modifications can affect gene expression profiles by changing the chromatin state and play an important role in regulating the growth and development of plants. When the specific epigenetic modifications are changed, crops can obtain excellent phenotypes and stronger environmental adaptability. Therefore, artificially changing the epigenetic modifications are expected to obtain high-quality germplasm resources more suitable for agricultural production. In this review, we summarize the main types of plant epigenetic modifications, highlight the research progresses of functional plant epigenetic modifications on the important traits and responses to environmental stress, and identify the main problems that need be solved in the application of epigenetics in crop improvement, thereby providing new insights for the functional epigenetic modifications on crop breeding and improvement.

Genome-wide identification, classification, and expression analysis of the JmjC domain-containing histone demethylase gene family in birch.

BACKGROUND: Histone methylation occurs primarily on lysine residues and requires a set of enzymes capable of reading, writing, and erasing to control its establishment and deletion, which is essential for maintaining chromatin structure and gene expression. Histone methylation and demethylation are contributed to plant growth and development, and are involved in adapting to environmental stresses. The JmjC domain-containing proteins are extensively studied for their function in histone lysine demethylation in plants, and play a critical role in sustaining histone methylation homeostasis. RESULTS: In this study, a total of 21 JmjC domain-containing histone demethylase proteins (JHDMs) in birch were identified and classified into five subfamilies based on structural characteristics and phylogenetic relationships among Arabidopsis, rice, maize, and birch. Although the BpJMJ genes displayed significant schematic variation, their distribution on the chromosomes is relatively uniform. Additionally, the BpJMJ genes in birch have never experienced a tandem-duplication event proved by WGD analysis and were remaining underwent purifying selection (Ka/Ks < < 1). A typical JmjC domain was found in all BpJMJ genes, some of which have other essential domains for their functions. In the promoter regions of BpJMJ genes, cis-acting elements associated with hormone and abiotic stress responses were overrepresented. Under abiotic stresses, the transcriptome profile reveals two contrasting expression patterns within 21 BpJMJ genes. Furthermore, it was established that most BpJMJ genes had higher expression in young tissues under normal conditions, with BpJMJ06/16 having the highest expression in germinating seeds and participating in the regulation of BpGA3ox1/2 gene expression. Eventually, BpJMJ genes were found to directly interact with genes involved in the "intracellular membrane" in respond to cold stress. CONCLUSIONS: The present study will provide a foundation for future experiments on histone demethylases in birch and a theoretical basis for epigenetic research on growth and development in response to abiotic stresses.

Harzianic Acid from Trichoderma afroharzianum Is a Natural Product Inhibitor of Acetohydroxyacid Synthase.

Acetohydroxyacid synthase (AHAS) is the first enzyme in the branched-chain amino acid biosynthetic pathway and is a validated target for herbicide and fungicide development. Here we report harzianic acid (HA, 1) produced by the biocontrol fungus Trichoderma afroharzianum t-22 (Tht22) as a natural product inhibitor of AHAS. The biosynthetic pathway of HA was elucidated with heterologous reconstitution. Guided by a putative self-resistance enzyme in the genome, HA was biochemically demonstrated to be a selective inhibitor of fungal AHAS, including those from phytopathogenic fungi. In addition, HA can inhibit a common resistant variant of AHAS in which the active site proline is mutated. Structural analysis of AHAS complexed with HA revealed the molecular basis of competitive inhibition, which differs from all known commercial AHAS inhibitors. The alternative binding mode also rationalizes the selectivity of HA, as well as effectiveness toward resistant mutants. A proposed role of HA biosynthesis by Tht22 in the rhizosphere is discussed based on the data.

Genome-wide identification and expression analysis of the AT-hook Motif Nuclear Localized gene family in soybean.

BACKGROUND: Soybean is an important legume crop and has significant agricultural and economic value. Previous research has shown that the AT-Hook Motif Nuclear Localized (AHL) gene family is highly conserved in land plants, playing crucial roles in plant growth and development. To date, however, the AHL gene family has not been studied in soybean. RESULTS: To investigate the roles played by the AHL gene family in soybean, genome-wide identification, expression patterns and gene structures were performed to analyze. We identified a total of 63 AT-hook motif genes, which were characterized by the presence of the AT-hook motif and PPC domain in soybean. The AT-hook motif genes were distributed on 18 chromosomes and formed two distinct clades (A and B), as shown by phylogenetic analysis. All the AHL proteins were further classified into three types (I, II and III) based on the AT-hook motif. Type-I was belonged to Clade-A, while Type-II and Type-III were belonged to Clade-B. Our results also showed that the main type of duplication in the soybean AHL gene family was segmented duplication event. To discern whether the AHL gene family was involved in stress response in soybean, we performed cis-acting elements analysis and found that AHL genes were associated with light responsiveness, anaerobic induction, MYB and gibberellin-responsiveness elements. This suggest that AHL genes may participate in plant development and mediate stress response. Moreover, a co-expression network analysis showed that the AHL genes were also involved in energy transduction, and the associated with the gibberellin pathway and nuclear entry signal pathways in soybean. Transcription analysis revealed that AHL genes in Jack and Williams82 have a common expression pattern and are mostly expressed in roots, showing greater sensitivity under drought and submergence stress. Hence, the AHL gene family mainly reacts on mediating stress responses in the roots and provide comprehensive information for further understanding of the AT-hook motif gene family-mediated stress response in soybean. CONCLUSION: Sixty-three AT-hook motif genes were identified in the soybean genome. These genes formed into two distinct phylogenetic clades and belonged to three different types. Cis-acting elements and co-expression network analyses suggested that AHL genes participated in significant biological processes. This work provides important theoretical basis for the understanding of AHLs biological functions in soybean.

EgMIXTA1, a MYB-Type Transcription Factor, Promotes Cuticular Wax Formation in Eustoma grandiflorum Leaves.

In the aerial plant organs, cuticular wax forms a hydrophobic layer that can protect cells from dehydration, repel pathogen attacks, and prevent organ fusion during development. The MIXTA gene encodes an MYB-like transcription factor, which is associated with epicuticular wax biosynthesis to increase the wax load on the surface of leaves. In this study, the AmMIXTA-homologous gene EgMIXTA1 was functionally characterized in the Eustoma grandiflorum. EgMIXTA1 was ubiquitously, but highly, expressed in leaves and buds. We identified the Eustoma MIXTA homolog and developed the plants for overexpression. EgMIXTA1-overexpressing plants had more wax crystal deposition on the leaf surface compared to wild-type and considerably more overall cuticular wax. In the leaves of the overexpression line, the cuticular transpiration occurred more slowly than in those of non-transgenic plants. Analysis of gene expression indicated that several genes, such as EgCER3, EgCER6, EgCER10, EgKCS1, EgKCR1, and EgCYP77A6, which are known to be involved in wax biosynthesis, were induced by EgMIXTA1-overexpression lines. Expression of another gene, WAX INDUCER1/SHINE1, encoding a transcription factor that stimulates the production of cutin, was also significantly higher in the overexpressors than in wild-type. However, the expression of a lipid-related gene, EgABCG12, did not change relative to the wild-type. These results suggest that EgMIXTA1 is involved in the biosynthesis of cuticular waxes.

Combinatorial Biosynthesis of Sulfated Benzenediol Lactones with a Phenolic Sulfotransferase from Fusarium graminearum PH-1.

Total biosynthesis or whole-cell biocatalytic production of sulfated small molecules relies on the discovery and implementation of appropriate sulfotransferase enzymes. Although fungi are prominent biocatalysts and have been used to sulfate drug-like phenolics, no gene encoding a sulfotransferase enzyme has been functionally characterized from these organisms. Here, we identify a phenolic sulfotransferase, FgSULT1, by genome mining from the plant-pathogenic fungus Fusarium graminearum PH-1. We expressed FgSULT1 in a Saccharomyces cerevisiae chassis to modify a broad range of benzenediol lactones and their nonmacrocyclic congeners, together with an anthraquinone, with the resulting unnatural natural product (uNP) sulfates displaying increased solubility. FgSULT1 shares low similarity with known animal and plant sulfotransferases. Instead, it forms a sulfotransferase family with putative bacterial and fungal enzymes for phase II detoxification of xenobiotics and allelochemicals. Among fungi, putative FgSULT1 homologues are encoded in the genomes of Fusarium spp. and a few other genera in nonsyntenic regions, some of which may be related to catabolic sulfur recycling. Computational structure modeling combined with site-directed mutagenesis revealed that FgSULT1 retains the key catalytic residues and the typical fold of characterized animal and plant sulfotransferases. Our work opens the way for the discovery of hitherto unknown fungal sulfotransferases and provides a synthetic biological and enzymatic platform that can be adapted to produce bioactive sulfates, together with sulfate ester standards and probes for masked mycotoxins, precarcinogenic toxins, and xenobiotics.IMPORTANCE Sulfation is an expedient strategy to increase the solubility, bioavailability, and bioactivity of nutraceuticals and clinically important drugs. However, chemical or biological synthesis of sulfoconjugates is challenging. Genome mining, heterologous expression, homology structural modeling, and site-directed mutagenesis identified FgSULT1 of Fusarium graminearum PH-1 as a cytosolic sulfotransferase with the typical fold and active site architecture of characterized animal and plant sulfotransferases, despite low sequence similarity. FgSULT1 homologues are sparse in fungi but form a distinct clade with bacterial sulfotransferases. This study extends the functionally characterized sulfotransferase superfamily to the kingdom Fungi and demonstrates total biosynthetic and biocatalytic synthetic biological platforms to produce unnatural natural product (uNP) sulfoconjugates. Such uNP sulfates may be utilized for drug discovery in human and veterinary medicine and crop protection. Our synthetic biological methods may also be adapted to generate masked mycotoxin standards for food safety and environmental monitoring applications and to expose precarcinogenic xenobiotics.