Metabolic changes in tomato plants caused by psychrotolerant Antarctic endophytic bacteria might be implicated in cold stress mitigation.

Climate change is responsible for mild winters and warm springs that can induce premature plant development, increasing the risk of exposure to cold stress with a severe reduction in plant growth. Tomato plants are sensitive to cold stress and beneficial microorganisms can increase their tolerance. However, scarce information is available on mechanisms stimulated by bacterial endophytes in tomato plants against cold stress. This study aimed to clarify metabolic changes stimulated by psychrotolerant endophytic bacteria in tomato plants exposed to cold stress and annotate compounds possibly associated with cold stress mitigation. Tomato seeds were inoculated with two bacterial endophytes isolated from Antarctic Colobanthus quitensis plants (Ewingella sp. S1.OA.A_B6 and Pseudomonas sp. S2.OTC.A_B10) or with Paraburkholderia phytofirmans PsJN, while mock-inoculated seeds were used as control. The metabolic composition of tomato plants was analyzed immediately after cold stress exposure (4°C for seven days) or after two and four days of recovery at 25°C. Under cold stress, the content of malondialdehyde, phenylalanine, ferulic acid, and p-coumaric acid was lower in bacterium-inoculated compared to mock-inoculated plants, indicating a reduction of lipid peroxidation and the stimulation of phenolic compound metabolism. The content of two phenolic compounds, five putative phenylalanine-derived dipeptides, and three further phenylalanine-derived compounds was higher in bacterium-inoculated compared to mock-inoculated samples under cold stress. Thus, psychrotolerant endophytic bacteria can reprogram polyphenol metabolism and stimulate the accumulation of secondary metabolites, like 4-hydroxybenzoic and salicylic acid, which are presumably involved in cold stress mitigation, and phenylalanine-derived dipeptides possibly involved in plant stress responses.

Tailored Mass Spectral Data Exploration Using the SpecXplore Interactive Dashboard.

Untargeted metabolomics promises comprehensive characterization of small molecules in biological samples. However, the field is hampered by low annotation rates and abstract spectral data. Despite recent advances in computational metabolomics, manual annotations and manual confirmation of in-silico annotations remain important in the field. Here, exploratory data analysis methods for mass spectral data provide overviews, prioritization, and structural hypothesis starting points to researchers facing large quantities of spectral data. In this research, we propose a fluid means of dealing with mass spectral data using specXplore, an interactive Python dashboard providing interactive and complementary visualizations facilitating mass spectral similarity matrix exploration. Specifically, specXplore provides a two-dimensional t-distributed stochastic neighbor embedding embedding as a jumping board for local connectivity exploration using complementary interactive visualizations in the form of partial network drawings, similarity heatmaps, and fragmentation overview maps. SpecXplore makes use of state-of-the-art ms2deepscore pairwise spectral similarities as a quantitative backbone while allowing fast changes of threshold and connectivity limitation settings, providing flexibility in adjusting settings to suit the localized node environment being explored. We believe that specXplore can become an integral part of mass spectral data exploration efforts and assist users in the generation of structural hypotheses for compounds of interest.

Untargeted Plant Metabolomics: Evaluation of Lyophilization as a Sample Preparation Technique.

Lyophilization is a common method used for stabilizing biological samples prior to storage or to concentrate extracts. However, it is possible that this process may alter the metabolic composition or lead to the loss of metabolites. In this study, the performance of lyophilization is investigated in the example of wheat roots. To this end, native and 13C-labelled, fresh or already lyophilized root samples, and (diluted) extracts with dilution factors up to 32 and authentic reference standards were investigated. All samples were analyzed using RP-LC-HRMS. Results show that using lyophilization for the stabilization of plant material altered the metabolic sample composition. Overall, 7% of all wheat metabolites detected in non-lyophilized samples were not detected in dried samples anymore, and up to 43% of the remaining metabolites exhibited significantly increased or decreased abundances. With respect to extract concentration, less than 5% of the expected metabolites were completely lost by lyophilization and the recovery rates of the remaining metabolites were slightly reduced with increasing concentration factors to an average of 85% at an enrichment factor of 32. Compound annotation did not indicate specific classes of wheat metabolites to be affected.

Genome analysis of Cephalotrichum gorgonifer and identification of the biosynthetic pathway for rasfonin, an inhibitor of KRAS dependent cancer.

BACKGROUND: Fungi are important sources for bioactive compounds that find their applications in many important sectors like in the pharma-, food- or agricultural industries. In an environmental monitoring project for fungi involved in soil nitrogen cycling we also isolated Cephalotrichum gorgonifer (strain NG_p51). In the course of strain characterisation work we found that this strain is able to naturally produce high amounts of rasfonin, a polyketide inducing autophagy, apoptosis, necroptosis in human cell lines and showing anti-tumor activity in KRAS-dependent cancer cells. RESULTS: In order to elucidate the biosynthetic pathway of rasfonin, the strain was genome sequenced, annotated, submitted to transcriptome analysis and genetic transformation was established. Biosynthetic gene cluster (BGC) prediction revealed the existence of 22 BGCs of which the majority was not expressed under our experimental conditions. In silico prediction revealed two BGCs with a suite of enzymes possibly involved in rasfonin biosynthesis. Experimental verification by gene-knock out of the key enzyme genes showed that one of the predicted BGCs is indeed responsible for rasfonin biosynthesis. CONCLUSIONS: This study identified a biosynthetic gene cluster containing a key-gene responsible for rasfonin production. Additionally, molecular tools were established for the non-model fungus Cephalotrichum gorgonifer which allows strain engineering and heterologous expression of the BGC for high rasfonin producing strains and the biosynthesis of rasfonin derivates for diverse applications.

Identification of a UDP-glucosyltransferase conferring deoxynivalenol resistance in Aegilops tauschii and wheat.

Aegilops tauschii is the diploid progenitor of the wheat D subgenome and a valuable resource for wheat breeding, yet, genetic analysis of resistance against Fusarium head blight (FHB) and the major Fusarium mycotoxin deoxynivalenol (DON) is lacking. We treated a panel of 147 Ae. tauschii accessions with either Fusarium graminearum spores or DON solution and recorded the associated disease spread or toxin-induced bleaching. A k-mer-based association mapping pipeline dissected the genetic basis of resistance and identified candidate genes. After DON infiltration nine accessions revealed severe bleaching symptoms concomitant with lower conversion rates of DON into the non-toxic DON-3-O-glucoside. We identified the gene AET5Gv20385300 on chromosome 5D encoding a uridine diphosphate (UDP)-glucosyltransferase (UGT) as the causal variant and the mutant allele resulting in a truncated protein was only found in the nine susceptible accessions. This UGT is also polymorphic in hexaploid wheat and when expressed in Saccharomyces cerevisiae only the full-length gene conferred resistance against DON. Analysing the D subgenome helped to elucidate the genetic control of FHB resistance and identified a UGT involved in DON detoxification in Ae. tauschii and hexaploid wheat. This resistance mechanism is highly conserved since the UGT is orthologous to the barley UGT HvUGT13248 indicating descent from a common ancestor of wheat and barley.

Gramiketides, Novel Polyketide Derivatives of Fusarium graminearum, Are Produced during the Infection of Wheat.

The plant pathogen Fusarium graminearum is a proficient producer of mycotoxins and other in part still unknown secondary metabolites, some of which might act as virulence factors on wheat. The PKS15 gene is expressed only in planta, so far hampering the identification of an associated metabolite. Here we combined the activation of silent gene clusters by chromatin manipulation (kmt6) with blocking the metabolic flow into the competing biosynthesis of the two major mycotoxins deoxynivalenol and zearalenone. Using an untargeted metabolomics approach, two closely related metabolites were found in triple mutants (kmt6 tri5 pks4,13) deficient in production of the major mycotoxins deoxynivalenol and zearalenone, but not in strains with an additional deletion in PKS15 (kmt6 tri5 pks4,13 pks15). Characterization of the metabolites, by LC-HRMS/MS in combination with a stable isotope-assisted tracer approach, revealed that they are likely hybrid polyketides comprising a polyketide part consisting of malonate-derived acetate units and a structurally deviating part. We propose the names gramiketide A and B for the two metabolites. In a biological experiment, both gramiketides were formed during infection of wheat ears with wild-type but not with pks15 mutants. The formation of the two gramiketides during infection correlated with that of the well-known virulence factor deoxynivalenol, suggesting that they might play a role in virulence.

Heterogeneous multimeric metabolite ion species observed in LC-MS based metabolomics data sets.

Covalent or non-covalent heterogeneous multimerization of molecules associated with extracts from biological samples analyzed via LC-MS are quite difficult to recognize/annotate and therefore the prevalence of multimerization remains largely unknown. In this study, we utilized 13C labeled and unlabeled Pichia pastoris extracts to recognize heterogeneous multimers. More specifically, between 0.8% and 1.5% of the biologically-derived features detected in our experiments were confirmed to be heteromers, about half of which we could successfully annotate with monomeric partners. Interestingly, we found specific chemical classes such as nucleotides to disproportionately contribute to heteroadducts. Furthermore, we compiled these compounds into the first MS/MS library that included data from heteromultimers to provide a starting point for other labs to improve the annotation of such ions in other metabolomics data sets. Then, the detected heteromers were also searched in publicly accessible LC-MS datasets available in Metabolights, Metabolomics WB and GNPS/MassIVE to demonstrate that these newly annotated ions are also relevant to other public datasets. Furthermore, in additional datasets (Triticum aestivum, Fusarium graminearum, and Trichoderma reesei) our developed workflow also detected 0.5%-4.9% of metabolite features to originate from heterodimers, demonstrating heteroadducts to be present in metabolomics studies at a low percentage.

Polaramycin B, and not physical interaction, is the signal that rewires fungal metabolism in the Streptomyces-Aspergillus interaction.

Co-culturing the bacterium Streptomyces rapamycinicus and the ascomycete Aspergillus nidulans has previously been shown to trigger the production of orsellinic acid (ORS) and its derivates in the fungal cells. Based on these studies it was assumed that direct physical contact is a prerequisite for the metabolic reaction that involves a fungal amino acid starvation response and activating chromatin modifications at the biosynthetic gene cluster (BGC). Here we show that not physical contact, but a guanidine containing macrolide, named polaramycin B, triggers the response. The substance is produced constitutively by the bacterium and above a certain concentration, provokes the production of ORS. In addition, several other secondary metabolites were induced by polaramycin B. Our genome-wide transcriptome analysis showed that polaramycin B treatment causes downregulation of fungal genes necessary for membrane stability, general metabolism and growth. A compensatory genetic response can be observed in the fungus that included upregulation of BGCs and genes necessary for ribosome biogenesis, translation and membrane stability. Our work discovered a novel chemical communication, in which the antifungal bacterial metabolite polaramycin B leads to the production of antibacterial defence chemicals and to the upregulation of genes necessary to compensate for the cellular damage caused by polaramycin B.

Towards a broader view of the metabolome: untargeted profiling of soluble and bound polyphenols in plants.

Phenylalanine (Phe) is a central precursor for numerous secondary plant metabolites with a multitude of biological functions. Recent studies on the fungal disease Fusarium head blight in wheat showed numerous Phe-derived defence metabolites to be induced in the presence of the pathogen. These studies also suggest a partial incorporation of Phe-derived secondary metabolites into the cell wall. To broaden the view of the metabolome to bound Phe derivatives, an existing approach using 13C-labelled Phe as tracer was extended. The developed workflow consists of three successive extractions with an acidified acetonitrile-methanol-water mixture to remove the soluble plant metabolites, followed by cell wall hydrolysis with 4M aqueous NaOH, acidification with aqueous HCl, and liquid-liquid extraction of the hydrolysate with ethyl acetate. The untargeted screening of Phe-derived metabolites revealed 156 soluble compounds and 90 compounds in the hydrolysed samples including known cell wall constituents like ferulic acid, coumaric acid, and tricin. Forty-nine metabolites were found exclusively in the hydrolysate. The average cumulative extraction yield of the soluble metabolites was 99.6%, with a range of 91.8 to 100%. Repeatability coefficients of variation of the protocol ranged from 10.5 to 25.9%, with a median of 16.3%. To demonstrate the suitability of the proposed method for a typical metabolomics application, mock-treated and Fusarium graminearum-treated wheat samples were compared. The study revealed differences between the hydrolysates of the two sample types, confirming the differential incorporation of Phe-derived metabolites into the cell wall under infection conditions.

PeakBot: machine-learning-based chromatographic peak picking.

MOTIVATION: Chromatographic peak picking is among the first steps in data processing workflows of raw LC-HRMS datasets in untargeted metabolomics applications. Its performance is crucial for the holistic detection of all metabolic features as well as their relative quantification for statistical analysis and metabolite identification. Random noise, non-baseline separated compounds and unspecific background signals complicate this task. RESULTS: A machine-learning-based approach entitled PeakBot was developed for detecting chromatographic peaks in LC-HRMS profile-mode data. It first detects all local signal maxima in a chromatogram, which are then extracted as super-sampled standardized areas (retention-time versus m/z). These are subsequently inspected by a custom-trained convolutional neural network that forms the basis of PeakBot's architecture. The model reports if the respective local maximum is the apex of a chromatographic peak or not as well as its peak center and bounding box. In training and independent validation datasets used for development, PeakBot achieved a high performance with respect to discriminating between chromatographic peaks and background signals (accuracy of 0.99). For training the machine-learning model a minimum of 100 reference features are needed to learn their characteristics to achieve high-quality peak-picking results for detecting such chromatographic peaks in an untargeted fashion. PeakBot is implemented in python (3.8) and uses the TensorFlow (2.5.0) package for machine-learning related tasks. It has been tested on Linux and Windows OSs. AVAILABILITY AND IMPLEMENTATION: The package is available free of charge for non-commercial use (CC BY-NC-SA). It is available at https://github.com/christophuv/PeakBot. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Fungal Melanin Biosynthesis Pathway as Source for Fungal Toxins.

Contamination of food and feed with toxin-producing fungi is a major threat in agriculture and for human health. The filamentous fungus Alternaria alternata is one of the most widespread postharvest contaminants and a weak plant pathogen. It produces a large variety of secondary metabolites with alternariol and its derivatives as characteristic mycotoxin. Other important phyto- and mycotoxins are perylene quinones (PQs), some of which have anticancer properties. Here, we discovered that the PQ altertoxin (ATX) biosynthesis shares most enzymes with the 1,8-dihydroxynaphthalene (1,8-DHN) melanin pathway. However, melanin was formed in aerial hyphae and spores, and ATXs were synthesized in substrate hyphae. This spatial separation is achieved through the promiscuity of a polyketide synthase, presumably producing a pentaketide (T4HN), a hexaketide (AT4HN), and a heptaketide (YWA1) as products. T4HN directly enters the altertoxin and DHN melanin pathway, whereas AT4HN and YWA1 can be converted only in aerial hyphae, which probably leads to a higher T4HN concentration, favoring 1,8-DHN melanin formation. Whereas the production of ATXs was strictly dependent on the CmrA transcription factor, melanin could still be produced in the absence of CmrA to some extent. This suggests that different cues regulate melanin and toxin formation. Since DHN melanin is produced by many fungi, PQs or related compounds may be produced in many more fungi than so far assumed. IMPORTANCE Mycotoxins are a major threat for human health. Food safety control relies on the identification of the toxins or the detection of the expression of the respective genes. The latter method, however, relies on the knowledge of the biosynthetic pathway and the key genes. Alternaria alternata is a major food contaminant and produces many different mycotoxins with altertoxins and other perylene quinones as prominent examples. Here, we discovered that the biosynthetic pathway for altertoxins shares most of the enzymes with the dihydroxynaphthalene (DHN) melanin pathway. Because the DHN melanin pathway is widespread among fungi, the production of mycotoxins of the perylene quinone class could be more widespread than so far anticipated.

The TOR kinase pathway is relevant for nitrogen signaling and antagonism of the mycoparasite Trichoderma atroviride.

Trichoderma atroviride (Ascomycota, Sordariomycetes) is a well-known mycoparasite applied for protecting plants against fungal pathogens. Its mycoparasitic activity involves processes shared with plant and human pathogenic fungi such as the production of cell wall degrading enzymes and secondary metabolites and is tightly regulated by environmental cues. In eukaryotes, the conserved Target of Rapamycin (TOR) kinase serves as a central regulator of cellular growth in response to nutrient availability. Here we describe how alteration of the activity of TOR1, the single and essential TOR kinase of T. atroviride, by treatment with chemical TOR inhibitors or by genetic manipulation of selected TOR pathway components affected various cellular functions. Loss of TSC1 and TSC2, that are negative regulators of TOR complex 1 (TORC1) in mammalian cells, resulted in altered nitrogen source-dependent growth of T. atroviride, reduced mycoparasitic overgrowth and, in the case of Δtsc1, a diminished production of numerous secondary metabolites. Deletion of the gene encoding the GTPase RHE2, whose mammalian orthologue activates mTORC1, led to rapamycin hypersensitivity and altered secondary metabolism, but had an only minor effect on vegetative growth and mycoparasitic overgrowth. The latter also applied to mutants missing the npr1-1 gene that encodes a fungus-specific kinase known as TOR target in yeast. Genome-wide transcriptome analysis confirmed TOR1 as a regulatory hub that governs T. atroviride metabolism and processes associated to ribosome biogenesis, gene expression and translation. In addition, mycoparasitism-relevant genes encoding terpenoid and polyketide synthases, peptidases, glycoside hydrolases, small secreted cysteine-rich proteins, and G protein coupled receptors emerged as TOR1 targets. Our results provide the first in-depth insights into TOR signaling in a fungal mycoparasite and emphasize its importance in the regulation of processes that critically contribute to the antagonistic activity of T. atroviride.

Luteapyrone, a Novel Æ´-Pyrone Isolated from the Filamentous Fungus Metapochonia lutea.

In the process of screening for new bioactive microbial metabolites we found a novel Æ´-pyrone derivative for which we propose the trivial name luteapyrone, in a recently described microscopic filamentous fungus, Metapochonia lutea BiMM-F96/DF4. The compound was isolated from the culture extract of the fungus grown on modified yeast extract sucrose medium by means of flash chromatography followed by preparative HPLC. The chemical structure was elucidated by NMR and LC-MS. The new compound was found to be non-cytotoxic against three mammalian cell lines (HEK 263, KB-3.1 and Caco-2). Similarly, no antimicrobial activity was observed in tested microorganisms (gram positive and negative bacteria, yeast and fungi).

Polyphasic Approach Utilized for the Identification of Two New Toxigenic Members of Penicillium Section Exilicaulis, P. krskae and P. silybi spp. nov.

Two new species, Penicillium krskae (isolated from the air as a lab contaminant in Tulln (Austria, EU)) and Penicillium silybi (isolated as an endophyte from asymptomatic milk thistle (Silybum marianum) stems from Josephine County (Oregon, USA)) are described. The new taxa are well supported by phenotypic (especially conidial ornamentation under SEM, production of red exudate and red pigments), physiological (growth at 37 °C, response to cycloheximide and CREA), chemotaxonomic (production of specific extrolites), and multilocus phylogenetic analysis using RNA-polymerase II second largest subunit (RPB2), partial tubulin (benA), and calmodulin (CaM). Both new taxa are resolved within the section Exilicaulis in series Restricta and show phylogenetic affiliation to P. restrictum sensu stricto. They produce a large spectrum of toxic anthraquinoid pigments, namely, monomeric anthraquinones related to emodic and chloremodic acids and other interesting bioactive extrolites (i.e., endocrocin, paxilline, pestalotin, and 7-hydroxypestalotin). Of note, two bianthraquinones (i.e., skyrin and oxyskyrin) were detected in a culture extract of P. silybi. Two new chloroemodic acid derivatives (2-chloro-isorhodoptilometrin and 2-chloro-desmethyldermoquinone) isolated from the exudate of P. krskae ex-type culture were analyzed by nuclear magnetic resonance (NMR) and liquid chromatography-mass spectrometry (LC-MS).

Enhanced Metabolome Coverage and Evaluation of Matrix Effects by the Use of Experimental-Condition-Matched 13C-Labeled Biological Samples in Isotope-Assisted LC-HRMS Metabolomics.

Stable isotope-assisted approaches can improve untargeted liquid chromatography-high resolution mass spectrometry (LC-HRMS) metabolomics studies. Here, we demonstrate at the example of chemically stressed wheat that metabolome-wide internal standardization by globally 13C-labeled metabolite extract (GLMe-IS) of experimental-condition-matched biological samples can help to improve the detection of treatment-relevant metabolites and can aid in the post-acquisition assessment of putative matrix effects in samples obtained upon different treatments. For this, native extracts of toxin- and mock-treated (control) wheat ears were standardized by the addition of uniformly 13C-labeled wheat ear extracts that were cultivated under similar experimental conditions (toxin-treatment and control) and measured with LC-HRMS. The results show that 996 wheat-derived metabolites were detected with the non-condition-matched 13C-labeled metabolite extract, while another 68 were only covered by the experimental-condition-matched GLMe-IS. Additional testing is performed with the assumption that GLMe-IS enables compensation for matrix effects. Although on average no severe matrix differences between both experimental conditions were found, individual metabolites may be affected as is demonstrated by wrong decisions with respect to the classification of significantly altered metabolites. When GLMe-IS was applied to compensate for matrix effects, 272 metabolites showed significantly altered levels between treated and control samples, 42 of which would not have been classified as such without GLMe-IS.

The Lipoxygenase Lox1 Is Involved in Light- and Injury-Response, Conidiation, and Volatile Organic Compound Biosynthesis in the Mycoparasitic Fungus Trichoderma atroviride.

The necrotrophic mycoparasite Trichoderma atroviride is a biological pest control agent frequently applied in agriculture for the protection of plants against fungal phytopathogens. One of the main secondary metabolites produced by this fungus is 6-pentyl-α-pyrone (6-PP). 6-PP is an organic compound with antifungal and plant growth-promoting activities, whose biosynthesis was previously proposed to involve a lipoxygenase (Lox). In this study, we investigated the role of the single lipoxygenase-encoding gene lox1 encoded in the T. atroviride genome by targeted gene deletion. We found that light inhibits 6-PP biosynthesis but lox1 is dispensable for 6-PP production as well as for the ability of T. atroviride to parasitize and antagonize host fungi. However, we found Lox1 to be involved in T. atroviride conidiation in darkness, in injury-response, in the production of several metabolites, including oxylipins and volatile organic compounds, as well as in the induction of systemic resistance against the plant-pathogenic fungus Botrytis cinerea in Arabidopsis thaliana plants. Our findings give novel insights into the roles of a fungal Ile-group lipoxygenase and expand the understanding of a light-dependent role of these enzymes.

Preparation of uniformly labelled 13C- and 15N-plants using customised growth chambers.

BACKGROUND: Stable isotopically labelled organisms have found wide application in life science research including plant physiology, plant stress and defense as well as metabolism related sciences. Therefore, the reproducible production of plant material enriched with stable isotopes such as 13C and 15N is of considerable interest. A high degree of enrichment (> 96 atom %) with a uniformly distributed isotope (global labelling) is accomplished by a continuous substrate supply during plant growth/cultivation. In the case of plants, 13C-labelling can be achieved by growth in 13CO2(g) atmosphere while global 15N-labelling needs 15N- containing salts in the watering/nutrient solution. Here, we present a method for the preparation of 13C and 15N-labelled plants by the use of closed growth chambers and hydroponic nutrient supply. The method is exemplified with durum wheat. RESULTS: In total, 330 g of globally 13C- and 295 g of 15N-labelled Triticum durum wheat was produced during 87 cultivation days. For this, a total of 3.88 mol of 13CO2(g) and 58 mmol of 15N were consumed. The degree of enrichment was determined by LC-HRMS and ranged between 96 and 98 atom % for 13C and 95-99 atom % for 15N, respectively. Additionally, the isotopically labelled plant extracts were successfully used for metabolome-wide internal standardisation of native T.durum plants. Application of an isotope-assisted LC-HRMS workflow enabled the detection of 652 truly wheat-derived metabolites out of which 143 contain N. CONCLUSION: A reproducible cultivation which makes use of climate chambers and hydroponics was successfully adapted to produce highly enriched, uniformly 13C- and 15N-labelled wheat. The obtained plant material is suitable to be used in all kinds of isotope-assisted research. The described technical equipment and protocol can easily be applied to other plants to produce 13C-enriched biological samples when the necessary specific adaptations e.g. temperature and light regime, as well as nutrient supply are considered. Additionally, the 15N-labelling method can also be carried out under regular glasshouse conditions without the need for customised atmosphere.

Stable Isotope-Assisted Plant Metabolomics: Combination of Global and Tracer-Based Labeling for Enhanced Untargeted Profiling and Compound Annotation.

Untargeted approaches and thus biological interpretation of metabolomics results are still hampered by the reliable assignment of the global metabolome as well as classification and (putative) identification of metabolites. In this work we present an liquid chromatography-mass spectrometry (LC-MS)-based stable isotope assisted approach that combines global metabolome and tracer based isotope labeling for improved characterization of (unknown) metabolites and their classification into tracer derived submetabolomes. To this end, wheat plants were cultivated in a customized growth chamber, which was kept at 400 ± 50 ppm 13CO2 to produce highly enriched uniformly 13C-labeled sample material. Additionally, native plants were grown in the greenhouse and treated with either 13C9-labeled phenylalanine (Phe) or 13C11-labeled tryptophan (Trp) to study their metabolism and biochemical pathways. After sample preparation, liquid chromatography-high resolution mass spectrometry (LC-HRMS) analysis and automated data evaluation, the results of the global metabolome- and tracer-labeling approaches were combined. A total of 1,729 plant metabolites were detected out of which 122 respective 58 metabolites account for the Phe- and Trp-derived submetabolomes. Besides m/z and retention time, also the total number of carbon atoms as well as those of the incorporated tracer moieties were obtained for the detected metabolite ions. With this information at hand characterization of unknown compounds was improved as the additional knowledge from the tracer approaches considerably reduced the number of plausible sum formulas and structures of the detected metabolites. Finally, the number of putative structure formulas was further reduced by isotope-assisted annotation tandem mass spectrometry (MS/MS) derived product ion spectra of the detected metabolites. A major innovation of this paper is the classification of the metabolites into submetabolomes which turned out to be valuable information for effective filtering of database hits based on characteristic structural subparts. This allows the generation of a final list of true plant metabolites, which can be characterized at different levels of specificity.

Stable Isotope-Assisted Plant Metabolomics: Investigation of Phenylalanine-Related Metabolic Response in Wheat Upon Treatment With the Fusarium Virulence Factor Deoxynivalenol.

The major Fusarium mycotoxin deoxynivalenol (DON) is a virulence factor in wheat and has also been shown to induce defense responses in host plant tissue. In this study, global and tracer labeling with 13C were combined to annotate the overall metabolome of wheat spikes and to evaluate the response of phenylalanine-related pathways upon treatment with DON. At anthesis, spikes of resistant and susceptible cultivars as well as two related near isogenic wheat lines (NILs) differing in the presence/absence of the major resistance QTL Fhb1 were treated with 1 mg DON or water (control), and samples were collected at 0, 12, 24, 48, and 96 h after treatment (hat). A total of 172 Phe-derived wheat constituents were detected with our untargeted approach employing 13C-labeled phenylalanine and subsequently annotated as flavonoids, lignans, coumarins, benzoic acid derivatives, hydroxycinnamic acid amides (HCAAs), as well as peptides. Ninety-six hours after the DON treatment, up to 30% of the metabolites biosynthesized from Phe showed significantly increased levels compared to the control samples. Major metabolic changes included the formation of precursors of compounds implicated in cell wall reinforcement and presumed antifungal compounds. In addition, also dipeptides, which presumably are products of proteolytic degradation of truncated proteins generated in the presence of the toxin, were significantly more abundant upon DON treatment. An in-depth comparison of the two NILs with correlation clustering of time course profiles revealed some 70 DON-responsive Phe derivatives. While several flavonoids had constitutively different abundance levels between the two NILs differing in resistance, other Phe-derived metabolites such as HCAAs and hydroxycinnamoyl quinates were affected differently in the two NILs after treatment with DON. Our results suggest a strong activation of the general phenylpropanoid pathway and that coumaroyl-CoA is mainly diverted towards HCAAs in the presence of Fhb1, whereas the metabolic route to monolignol(-conjugates), lignans, and lignin seems to be favored in the absence of the Fhb1 resistance quantitative trait loci.

Methanol Generates Numerous Artifacts during Sample Extraction and Storage of Extracts in Metabolomics Research.

Many metabolomics studies use mixtures of (acidified) methanol and water for sample extraction. In the present study, we investigated if the extraction with methanol can result in artifacts. To this end, wheat leaves were extracted with mixtures of native and deuterium-labeled methanol and water, with or without 0.1% formic acid. Subsequently, the extracts were analyzed immediately or after storage at 10 °C, -20 °C or -80 °C with an HPLC-HESI-QExactive HF-Orbitrap instrument. Our results showed that 88 (8%) of the >1100 detected compounds were derived from the reaction with methanol and either formed during sample extraction or short-term storage. Artifacts were found for various substance classes such as flavonoids, carotenoids, tetrapyrrols, fatty acids and other carboxylic acids that are typically investigated in metabolomics studies. 58 of 88 artifacts were common between the two tested extraction variants. Remarkably, 34 of 73 (acidified extraction solvent) and 33 of 73 (non-acidified extraction solvent) artifacts were formed de novo as none of these meth(ox)ylated metabolites were found after extraction of native leaf samples with CD3OH/H2O. Moreover, sample extracts stored at 10 °C for several days, as can typically be the case during longer measurement sequences, led to an increase in both the number and abundance of methylated artifacts. In contrast, frozen sample extracts were relatively stable during a storage period of one week. Our study shows that caution has to be exercised if methanol is used as the extraction solvent as the detected metabolites might be artifacts rather than natural constituents of the biological system. In addition, we recommend storing sample extracts in deep freezers immediately after extraction until measurement.