Antiparasitic Ovalicin Derivatives from Pseudallescheria boydii, a Mutualistic Fungus of French Guiana Termites

Social insects are in mutualism with microorganisms, contributing to their resistance against infectious diseases. The fungus Pseudallescheria boydii SNB-CN85 isolated from termites produces ovalicin derivatives resulting from the esterification of the less hindered site of the ovalicin epoxide by long-chain fatty acids. Their structures were elucidated using spectroscopic analysis and semisynthesis from ovalicin. For ovalicin, these compounds displayed antiprotozoal activities against Plasmodium falciparum and Trypanosoma brucei, with IC50 values of 19.8 and 1.1 µM, respectively, for the most active compound, i.e., ovalicin linoleate. In parallel, metabolomic profiling of a collection of P. boydii strains associated with termites made it possible to highlight this class of compounds together with tyroscherin derivatives in all strains. Finally, the complete genome of P. boydii strains was obtained by sequencing, and the cluster of potential ovalicin and ovalicin biosynthesis genes was annotated. Through these metabolomic and genomic analyses, a new ovalicin derivative named boyden C, in which the 6-membered ring of ovalicin was opened by oxidative cleavage, was isolated and structurally characterized.


Introduction
Insects might be the most represented and diverse living creatures on earth [1]. In nature, insects live in symbiosis with bacteria and fungi to protect themselves from pathogens [2]. These microorganisms play a key role in the ecological success of their host [3], providing them antipathogenic metabolites that could eventually find applications in therapeutics [4]. Insect/microorganism associations are well described for social insects, such as ants, bees, and some wasp species [5,6]. However, little work has been conducted thus far on termite/microorganism interactions (outside the context of trophobioses) [7][8][9][10].
Termites (Blattodea: Isoptera) are eusocial insects comprising more than 3000 species, mostly distributed in tropical and subtropical regions. All termites live in nests containing 1 to 1.5 million individuals divided into castes (reproductive individuals, workers, and soldiers) that cooperatively perform colony tasks. These insects play an important role in ecosystems through the decomposition of plant lignocellulose, having a positive impact in nutrient cycling and soil-water dynamics through their nest construction and tunnel networks in the soil [11]. Termites are vital members of old-growth tropical forests, being perhaps the main decomposers of dead plant material at all stages of humification (decay). Therefore, termite colonies viability is considered a marker of human influence on the environment [12,13]. Nests protect termites from aggressions of biotic and abiotic origin In parallel, new methodologies for dereplication of complex mixtures have been developed to accelerate and expand the discovery of new bioactive natural molecules. In order to avoid focusing research efforts on already known molecules, dereplication strategies have been developed [29,30]. Dereplication allows for the annotation of known molecules in a mixture but also allows for targeting and focusing on potentially novel molecules to be isolated and characterized. Mass spectrometry-based molecular networking approaches can be used to overcome these difficulties, allowing the targeting of specific metabolites through the tandem mass spectrometry data set organization they confer [29][30][31][32][33][34][35]. In this study, a molecular network approach allowed us to navigate the chemical space of P. boydii metabolomes and was used to annotate the chemical diversity of these fungi. A new compound was identified by this method and was isolated and fully characterized. Then, it was chosen to compare the metabolomic and genomic data of all P. boydii in our collection. For this purpose, the whole genome sequence of the four P. boydii strains was established, and the cluster of potential genes for biosynthesis of ovalicin and its fatty acid derivatives was annotated.

Results
The study of the P. boydii SNB-CN85 ethyl acetate extract led to the isolation of ovalicin (1) and the five ovalicin analogs 2-6 ( Figure 1).
Preliminary inspection of the 1 H and 13 C NMR spectra confirmed that Compounds 2-5 were ovalicin analogs (Figures S1-S29 in the Supplementary Materials, Table 1, Tables S1-S5 in the Supplementary Materials). In fact, Compounds 2-5 exhibited the same 1 H and 13 C NMR shifts as ovalicin, and the only difference was observed at position 14 with protons and carbons at δH 4.11 and 4.17/δC 69.4-69.5 ppm instead of δH 2.63 and 2.97/δC 51.0 ppm for ovalicin. For Compound 6, the same chemical shift but also the loss of the Omethyl group at position 7 was observed (Table S5 in    Preliminary inspection of the 1 H and 13 C NMR spectra confirmed that Compounds 2-5 were ovalicin analogs (Figures S1-S29 in the Supplementary Materials, Table 1, Tables S1-S5 in the Supplementary Materials). In fact, Compounds 2-5 exhibited the same 1 H and 13 C NMR shifts as ovalicin, and the only difference was observed at position 14 with protons and carbons at δ H 4.11 and 4.17/δ C 69.4-69.5 ppm instead of δ H 2.63 and 2.97/δ C 51.0 ppm for ovalicin. For Compound 6, the same chemical shift but also the loss of the O-methyl group at position 7 was observed (  13 C, HSQC spectra and 1 H-1 H COSY and 1 H-13 C HMBC correlations indicated the presence of an octadecadienoate chain, as shown by the presence of four olefinic methines (4H at δ H 5.33-5.36), 12 methylenes (δ H 1.28-2.78), and one triplet methyl (δ H 0.91). The chain was also determined following the sequence of 1 [36,37]. Due to too small quantities, this procedure was not possible for Compounds 5 and 6. To confirm the structures and configurations of Compounds 2, 3, and 4, ovalicin esters were synthesized from Compound 1 by esterification with linoleic acid, oleic acid, and stearic acid, respectively ( Figure S34 in the Supplementary Materials). The comparison of the NMR spectra obtained between the isolated and synthetic compounds allowed us to confirm the identification and configurations of fungal metabolites 2-4 ( Figures S35-S43 in the Supplementary Materials). The absolute configurations Compounds 2-4 were also confirmed first by comparison with circular dichroism data of 1 ( Figure S33 in the Supplementary Materials) for the ovalicin moiety. Furthermore, all these compounds showed negative optical rotation. Secondly, for the fatty acid moiety, our 13 C NMR analyses were in accordance with literature data concerning the specific chemical shifts of methylene at C-11 position between two double bonds ( Table S6 in the Supplementary Materials) [38,39]. lowing the sequence of 1 H-1 H COSY cross-peaks H-8′/H-9′/H-10′/H-11′/H-12′/H-13′. The fatty acid chain was linked to C-14 thanks to the H-14/C-1′ correlation ( Figure 2). Compounds 3-4 and 6 presented similar spectral data. The only difference was observed in the fatty acid chains. The exact mass of each molecule, the presence of olefinic protons between δH 5.33 and 5.36 ppm, and the 1 H-1 H COSY correlations confirmed the presence of two double bonds for 2 and 6, one double bond for 3 and 5, and none for 4 ( Figure 2).   [36,37]. Due to too small quantities, this procedure was not possible for Compounds 5 and 6. To confirm the structures and configurations of Compounds 2, 3, and 4, ovalicin esters were synthesized from Compound 1 by esterification with linoleic acid, oleic acid, and stearic acid, respectively ( Figure S34 in the Supplementary Materials). The comparison of the NMR spectra obtained between the isolated and synthetic compounds allowed us to confirm the identification and configurations of fungal metabolites 2-4 ( Figures S35-S43 in the Supplementary Materials). The absolute configurations Compounds 2-4 were also confirmed first by comparison with circular dichroism data of 1 ( Figure S33 in the Supplementary Materials) for the ovalicin moiety. Furthermore, all these compounds showed negative optical rotation. Secondly, for the fatty acid moiety, our 13 C NMR analyses were in accordance with literature data concerning the specific chemical shifts of methylene at C-11′ position between two double bonds ( Table S6 in the Supplementary Materials) [38,39].
Compounds 1-4 antiplasmodial activity against Plasmodium falciparum and antitrypanosomal activity against Trypanosoma brucei brucei are summarized in Table 2. Ovalicin 1 and fatty acid ester of ovalicin 2 and 4 showed high inhibition percentages at IC50 values Compounds 1-4 antiplasmodial activity against Plasmodium falciparum and antitrypanosomal activity against Trypanosoma brucei brucei are summarized in Table 2. Ovalicin 1 and fatty acid ester of ovalicin 2 and 4 showed high inhibition percentages at IC 50 values below 20.6 µM against Plasmodium falciparum. Against Trypanosoma brucei brucei, all compounds exhibited IC 50 values below 5 µM, but 1 to 2 were the only compounds to show activity (IC 50 < 1.1 µM). Compounds 5 and 6 were isolated in too small amounts and could not be tested. To the best of our knowledge, only ovalicin is known to show antiangiogenic properties [40] and potent antiplasmodial activity [41]. Consequently, the isolated compounds were also tested for their cytotoxicity on normal primary human umbilical vein endothelial cells (HUVECs). All compounds demonstrated low cytotoxicity on the tested cell lines (IC 50 > 20 µM), contrary to reported literature for ovalicin with an IC 50 of 13 µM [42,43]. Finally, isolated compounds were also tested for their potential antimicrobial activities against C. albicans. The minimal inhibitory concentrations (MICs) observed for molecules 1-3 were greater than 128 µg/mL, showing no significant inhibition. As it has been previously shown that strains of P. boydii were able to biosynthesize derivatives of tyrocherin and ovalicin [20][21][22], an evaluation was conducted of whether our different strains were also capable of biosynthesizing the fatty acid derivatives of ovalicin isolated from SNB-CN85. For this purpose, EtOAc extracts of four P. boydii strains, SNB-CN71, -CN73, -CN81, and -CN85, were analyzed by reversed-phase liquid chromatography-high-resolution mass spectrometry (RPPLC-HRMS/MS) using the datadependent acquisition mode. The resulting spectral data were further investigated by the molecular network approach (MN) [29][30][31]. Data were processed by MZmine 2.51 [32,33], and an MN was generated by Metgem [34]. The resulting network grouped 3021 nodes into 83 clusters ( Figure 3 and Figure S44 in the Supplementary Materials). Strain mapping was applied by attributing a given color code to each strain in the collection, allowing us to visually navigate through the MN and check for the presence and distribution of specific compound classes within the collection. Relative quantification of the ions was represented using pie charts, with their proportions based on the respective areas of the peaks detected on extracted ion chromatography (XIC).
Molecules 2022, 27, x FOR PEER REVIEW  8 of 21 tions, although they were grown under the same conditions. Compound 6 was not detected in the cluster, likely because this compound is a demethylated derivative of ovalicin but mostly because it seems to be present in very small quantities in the crude extract. Following these results, the biosynthetic pathway of ovalicin was investigated. For this purpose, whole genome sequencing of the four P. boydii strains SNB-CN71, -CN73, -CN81, and -CN85 was performed by combining short and long sequences to achieve de novo hybrid assemblies. The quality of the four 42-43 Mb genomes was assessed by the BUSCO score (95.9-96.1% completeness) ( Table S6 in the Supplementary Materials) [44]. It was considered that ovalicin is a known fumagillol analog, an intermediate in the fumagillin biosynthetic pathway [45,46]. While the fumagillin biosynthetic pathway has been described [47], this is not the case for ovalicin and its fatty ester derivatives. Furthermore, β-trans-bergamotene is a biosynthetic intermediate that ovalicin shares with fumagillin and fumagillol [45,46]. When our four genomes were analyzed using the antiSMASH pipeline [48], a cluster of biosynthetic genes common to all four strains, comprising enzymes responsible for producing ovalicin and its fatty acid esters, was identified. This cluster, composed of 18 genes conserved in the four P. boydii strains, was intertwined with the pseurotin A biosynthetic gene cluster (6 genes out of 18). Pseurotin A is another specialized metabolite produced by P. boydii strains whose protonated species were identified in the molecular networks ( Figure 3, Figure S44 in the Supplementary Materials). These two gene clusters appear to be frequently intertwined, as the pseurotin/fumagillin gene , were present in a large relative proportion in SNB-CN71, -CN73, and -CN81. N-Methyltyrocherin and pseurotin A were not detected in SNB-CN85 and SNB-CN71, but they were previously isolated from these strains. Therefore, pseudallicins A-D isolated previously from SNB-CN85 extract strains were not found during this analysis. This may be because the concentration of these compounds in the crude extracts appeared to be very low, whereas they had been concentrated in the previously studied fractions. Sodium-cationized ovalicin and Compounds 1-5 were detected in a specific cluster present in all extracts. By studying the fragmentation spectra of the compounds grouped in the ovalicin fatty acid cluster, it was possible to determine the presence of ovalicin analogs bearing C18 fatty acid chains with three double bonds (Compound 9), one to two epoxides (Compounds 10, 11 and 12), or with a C16 fatty acid chain (Compounds 5 and 8). Compound 11 exhibited both an epoxide and a double bond ( Figure S45 in the Supplementary Materials). It is interesting to note that this collection of fatty acid esters of ovalicin is quite large, and some of these specialized metabolites are present in the different species of the collection. Finally, not all strains produced the different derivatives in the same proportions, although they were grown under the same conditions. Compound 6 was not detected in the cluster, likely because this compound is a demethylated derivative of ovalicin but mostly because it seems to be present in very small quantities in the crude extract.
Following these results, the biosynthetic pathway of ovalicin was investigated. For this purpose, whole genome sequencing of the four P. boydii strains SNB-CN71, -CN73, -CN81, and -CN85 was performed by combining short and long sequences to achieve de novo hybrid assemblies. The quality of the four 42-43 Mb genomes was assessed by the BUSCO score (95.9-96.1% completeness) ( Table S6 in the Supplementary Materials) [44]. It was considered that ovalicin is a known fumagillol analog, an intermediate in the fumagillin biosynthetic pathway [45,46]. While the fumagillin biosynthetic pathway has been described [47], this is not the case for ovalicin and its fatty ester derivatives. Furthermore, β-trans-bergamotene is a biosynthetic intermediate that ovalicin shares with fumagillin and fumagillol [45,46]. When our four genomes were analyzed using the antiSMASH pipeline [48], a cluster of biosynthetic genes common to all four strains, comprising enzymes responsible for producing ovalicin and its fatty acid esters, was identified. This cluster, composed of 18 genes conserved in the four P. boydii strains, was intertwined with the pseurotin A biosynthetic gene cluster (6 genes out of 18). Pseurotin A is another specialized metabolite produced by P. boydii strains whose protonated species were identified in the molecular networks ( Figure 3, Figure S44 in the Supplementary Materials). These two gene clusters appear to be frequently intertwined, as the pseurotin/fumagillin gene cluster has also been reported in Aspergillus fumigatus [49]. The eight genes coding for enzymes involved in the ovalicin biosynthetic pathway were named Pbo_Ova_A-H ( Figure 4, Table S7 in the Supplementary Materials), and the six genes coding for pseurotin A were named Pbo_Pse_A-F (Table S8 in the Supplementary Materials). Their functions were predicted using the Pfam and UniProt databases [50,51]. Like fumagillin, the ovalicin biosynthetic pathway is initiated by the action of the terpene synthase Pbo_Ova_A, which converts farnesyl pyrophosphate to β-trans-bergamoten (Btb). A cytochrome P450 monooxygenase Pbo_Ova_B whose af510 is the closest homolog (85% similarity) was predicted to catalyze the formation of intermediate (Btb1) Figure S46 in the Supplementary Materials). Finally, fatty acids were likely added to dihydroovalicin via an acyltransferase, the most likely candidate of which was Pbo_Ova_G, whose gene was present at the edge of the Pbo_Ova/Pse biosynthetic gene cluster.  This biosynthetic scheme was fully consistent with experimental data proposed previously, concerning configurations of compounds 2-6 and sheds light on subunits' ovalicin absolute configurations.
To more deeply explore the chemical diversity of the P. boydii metabolome, we looked for nodes with intense signals but could not annotate them by searching standards and analogs in MetGem databases. It was chosen to isolate the specialized metabolites showing an ion signal at m/z 323.1433 (Figure 3 and Figure S44 in the Supplementary Materials), which could be a potential original derivative of ovalicin in view of the observed mass. The molecular formula of Compound 7 was determined to be C 15 H 24 O 6 by ESI-HRMS (m/z 323.1473 [M + Na] + , calcd. for C 15 H 24 O 6 Na + , 323.1471, err 1.3 ppm), indicating the presence of two additional oxygen and the lack of one carbon compared with the molecular formula of ovalicin (1). The 1 H NMR spectra of Compounds 1 and 7 displayed some similarities, suggesting that 7 was an analog of 1. The remarkable difference between Compounds 1 and 7 in 1 H NMR was the absence, in 7, of the singlet methoxy signal integrating for 3H at δ H 3.49 (Table 3).
However, Compound 7 exhibited characteristic dimethylallyl signals at δ H 5.24 (m, 1H), δ H 1.73 (m, 3H), and 1.67 (m, 3H). The examination of the COSY correlations allowed us to link this dimethylallyl moiety to H-8 at 3.26, (t, J = 6.5 Hz, 1H) [8]. The complete proton assignment of Compound 7 was deduced from the careful examination of COSY, HSQC, and HMBC correlations ( Figure 5). HMBC analysis indicated that methylene protons at position 2 at δ H 2.22 ppm correlated with a carbonyl carbon of a carboxylic acid at δ C 178.0 ppm. The analysis of the HSQC correlations showed that Compound 7 had additional oxygenated methylene proton signals at δ H 3.64 and 3.60 ppm. Consequently, we deduced that Compound 7 was an analog of demethylated ovalicin in which the C-1-C-7 bond was cleaved by oxidation or hydrolysis. The relative and absolute stereochemistry at C-7 was defined as drawn in Figure 5 based on biosynthetic considerations (Figure 4). Compound 7 was named boyden C. Unfortunately, as this compound was isolated in small quantities, it could not be biologically tested.

General Experimental Procedures
Optical rotations ([α]D) were measured on an Anton Paar MCP 300 polarimeter in a 100 mm long 350 µL cell (Anton Paar, Graz, Austria). IR spectra were recorded using a Perkin Elmer Spectrometer BX FT-IR (Thermo Fisher scientific, Les Ulis, France). Nuclear magnetic resonance (NMR) spectra were recorded on a Bruker 600 MHz spectrometer equipped with a 3 mm inverse detection cryoprobe or a Bruker 500 and 800 MHz spectrometer equipped with a 5 mm inverse detection probe (Bruker, Rheinstetten, Germany). Chemical shifts (δ) are reported as ppm based on the TMS signal. HR-ESI-MS measurements were performed using a Waters Acquity UPLC system with column bypass coupled to a Waters LCT Premier time-of-flight mass spectrometer equipped with an electrospray interface (ESI) (Waters, Manchester, England). Flash chromatography was performed on a Grace Reveleris system with dual UV and ELSD detection equipped with a 40 g C18 column. Analytical and preparative HPLCs were conducted with a Gilson system equipped with a 322 pump device, a GX-271 fraction collector, a 171 diode array detector, and a prepELSII electrospray nebulizer detector (Gilson, Middelton, USA). Columns used for these experiments included a Phenomenex Luna C18 5 µm 4.6 × 250 mm analytical column and a Phenomenex Luna C18 5 µm 21.2 × 250 mm preparative column. The flow rate was set to 1 or 21 mL/min using a linear gradient of H2O mixed with an increasing proportion of CH3CN. Both solvents were HPLC grade and modified with 0.1% formic acid. Potato dextrose agar (PDA) was purchased from Fluka Analytical (Janssen Pharma-

General Experimental Procedures
Optical rotations ([α]D) were measured on an Anton Paar MCP 300 polarimeter in a 100 mm long 350 µL cell (Anton Paar, Graz, Austria). IR spectra were recorded using a Perkin Elmer Spectrometer BX FT-IR (Thermo Fisher scientific, Les Ulis, France). Nuclear magnetic resonance (NMR) spectra were recorded on a Bruker 600 MHz spectrometer equipped with a 3 mm inverse detection cryoprobe or a Bruker 500 and 800 MHz spectrometer equipped with a 5 mm inverse detection probe (Bruker, Rheinstetten, Germany). Chemical shifts (δ) are reported as ppm based on the TMS signal. HR-ESI-MS measurements were performed using a Waters Acquity UPLC system with column bypass coupled to a Waters LCT Premier time-of-flight mass spectrometer equipped with an electrospray interface (ESI) (Waters, Manchester, England). Flash chromatography was performed on a Grace Reveleris system with dual UV and ELSD detection equipped with a 40 g C18 column. Analytical and preparative HPLCs were conducted with a Gilson system equipped with a 322 pump device, a GX-271 fraction collector, a 171 diode array detector, and a prepELSII electrospray nebulizer detector (Gilson, Middelton, WI, USA). Columns used for these experiments included a Phenomenex Luna C18 5 µm 4.6 × 250 mm analytical column and a Phenomenex Luna C18 5 µm 21.2 × 250 mm preparative column. The flow rate was set to 1 or 21 mL/min using a linear gradient of H 2 O mixed with an increasing proportion of CH 3 CN. Both solvents were HPLC grade and modified with 0.1% formic acid. Potato dextrose agar (PDA) was purchased from Fluka Analytical (Janssen Pharmaceuticalaan, Geel, Belgium). Molecular analyses were performed externally by BACTUP, Saint-Priest, France. MS/MS experiments were performed by direct injection analysis (DIA, 10 µL/min) on an Agilent Q-ToF 6540 (APCI interface) mass spectrometer (Agilent, Santa Clara, CA, USA). Compounds were dissolved in CH 2 Cl 2 at 2 mg/mL and then diluted in CH 3 CN at 200 µg/mL. The corona current was fixed to 2 µA, while the desolvation gas flow was set to 8 L/min (heated at 300 • C). The capillary, fragmentor, and skimmer voltages were set to 3000 V, 100 V, and 45 V, respectively. A collision energy of 15 eV was used.

Collection and Identification of Pseudallescheria Boydii
Termite workers were collected from an aerial termite nest located in Rémire-Montjoly, French Guiana, in July 2011. The termites were identified as Termes cf. hispaniolae and Nasutitermes corniger by Prof. Reginaldo Constantino (University of Brasilia, Brasilia, Brazil). The workers were surface-sterilized by successive soakings in 70% EtOH (2 min), 5% NaOCl (2 min), and sterile water. The termite was subsequently placed in a Petri dish containing solid PDA medium. After 1 week at 25 • C, the first fungal hyphae to emerge from the insect were sampled and transferred into other Petri dishes. One of the microbial colonies consisted of a pure fungus, which was saved in triplicate at −80 • C in H 2 O−glycerol (50/50). A sample was submitted for amplification of the nuclear ribosomal internal transcribed spacer region ITS. Sequencing allowed for strain identification by NCBI sequence comparison (Blastn ® , Rockville, MD, USA). The sequences have been registered in the NCBI GenBank database (http://www.ncbi.nlm.nih.gov, accessed on 20 December 2021) under the registry numbers KJ023743 (SNB-CN85), KC684883 (SNB-CN73), KJ023735 (SNB-CN71), and KJ023720 (SNB-CN81).

Culture, Extraction, and Isolation
The P. boydii strains were cultivated on PDA at 26 • C for 15 days, initially on a small scale. The fungus and culture medium were then transferred into a large container and macerated with EtOAc for 24 h. The organic solvent was then collected by filtration, washed with H 2 O in a separatory funnel and evaporated, yielding the four crude extracts.
The P. boydii strain SNB-CN85 was cultivated on a large scale in 150 × 14 cm Petri dishes, yielding 3.57 g extract. A portion of the extract (1.13 g) was subjected to reversedphase flash chromatography using a gradient of H 2 O mixed with an increasing proportion of CH 3 CN to afford 6 fractions (A to E). The column was then eluted with CH 3 CN/CH 2 Cl 2 1:1 (Fraction F). Fraction D (300 mg, eluted with H 2 O/CH 3 CN, 2:8) was purified by reversedphase flash chromatography using the same solvent system as mentioned above with a modified gradient. The major compound was collected in pure form and was identified as ovalicin (1) Figure S34).

Esterification
The free fatty acids were dissolved in MeOH (500 µL). Each solution was cooled down to 0 • C. SOCl 2 (20 µL, 275.4 µmol) was added, and the solution was stirred for 1 h at 50 • C. The solution was concentrated under reduced pressure, and H 2 O (500 µL) was added. The aqueous solution was extracted with a mixture of Et 2 O/heptane, 5:95 (2 mL). The organic layers were filtered on a Na 2 SO 4 /SiO 2 column and concentrated to afford fatty acid methyl esters.  HCl (0.1 mL), the aqueous solution was extracted with Et 2 O (2 × 2.0 mL). The combined organic layers were washed with water (2 mL), filtered, and evaporated to afford a crude mixture (32.2 mg), which was subjected to preparative HPLC (Luna C18, isocratic 5:95 elution over 20 min) to afford synthetic ovalicin linoleate (2)

Identification of the Double-Bond Position in Fatty Acid Methyl Esters by Liquid Chromatography/Atmospheric Pressure Chemical Ionization Mass Spectrometry
The double-bond positions in the fatty acid methyl esters synthetized from ovalicin linoleate (2) and ovalicin oleate (3)  Compounds were eluted at a flow rate of 0.4 mL·min −1 with a gradient from 5% B to 100% B in 20 min and then 100% B for 3 min. The injection volume was fixed at 5 µL for all analyses. Mass spectra were recorded in electrospray ionization positive ion mode with the following parameters: gas temperature 325 • C, drying gas flow rate 10 L·min −1 , nebulizer pressure 30 psi, sheath gas temperature 350 • C, sheath gas flow rate 10 L·min −1 , capillary voltage 3500 V, nozzle voltage 500 V, fragmentor voltage 130 V, skimmer voltage 45 V, Octopole 1 RF Voltage 750 V. Internal calibration was achieved with two calibrants, purine and hexakis (1H,1H,3Htetrafluoropropoxy) phosphazene (m/z 121.0509 and m/z 922.0098), providing a high mass accuracy better than 10 ppm. The data-dependent MS/MS events were acquired for the five most intense ions detected by full-scan MS, from the 200-1000 m/z range, above an absolute threshold of 1000 counts. Selected precursor ions were fragmented at a fixed collision energy of 30 eV and with an isolation window of 1.3 amu. The mass range of the precursor and fragment ions was set as m/z 200-1000.

Data Processing and Analysis
The data files were converted from the .d standard data format (Agilent Technologies, Waldbronn, Germany) to .mzXML format using MSConvert software, part of the ProteoWizard package 3.0 [53]. All mzxml values were processed using MZmine2v51, as previously described [33,34]. Mass detection was realized with an MS1 noise level of 1000 and an MS/MS noise level of 0. The ADAP chromatogram builder was employed with a minimum group size of scans of 3, a group intensity threshold of 1000, a minimum highest intensity of 1000, and m/z tolerance of 0.008 (or 20 ppm). Deconvolution was performed with the ADAP wavelet algorithm according to the following settings: S/N threshold = 10, minimum feature height = 1000, coefficient/area threshold = 10, peak duration range 0.01-1.5 min, and t R wavelet range 0.00-0.04 min. MS/MS scans were paired using a m/z tolerance range of 0.05 Da and t R tolerance range of 0.5 min. Isotopologues were grouped using the isotopic peak grouper algorithm with a m/z tolerance of 0.008 (or 20 ppm) and a t R tolerance of 0.2 min. Peaks were filtered using a feature list row filter, keeping only peaks with MS/MS scans (GNPS). Adduct identification, i.e., sodium-or potassium-cationized species, was performed on the peak list with a retention time tolerance of 0.1 min, a m/z tolerance of 0.008 or 20 ppm, and a maximum relative peak height of 150%. A complex search, such as dimers, was performed with a retention time tolerance of 0.1 min, a m/z tolerance of 0.008 or 20 ppm, and a maximum relative peak height of 150%. Peak alignment was performed using the join aligner with a m/z tolerance of 0.008 (or 20 ppm), a weight for m/z at 20, a retention time tolerance of 0.2 min, and weight for t R at 50. The MGF file and the metadata were generated using the export/submit to GNPS option (https://doi.org/10.5281/zenodo.5733173, accessed on 10 January 2019).
Molecular networks were calculated and visualized using MetGem 1.34 software (https://metgem.github.io/, accessed on 1 March 2019) [34], and MS/MS spectra were window-filtered by choosing only the top 6 peaks in the ±50 Da window throughout the spectrum. The data were filtered by removing all peaks in the ±17 Da range around the precursor m/z. The m/z tolerance windows used to find the matching peaks were set to 0.02 Da, and cosine scores were kept in consideration for spectra sharing at least 2 matching peaks. The number of iterations, perplexity, learning rate, and early exaggeration parameters were set to 5000, 25, 200, and 12 for the t-SNE view.
3.8. Genomic Analyses 3.8.1. gDNA Extraction of P. boydii Strains SNB-CN71, -CN73, -CN81, and -CN85 For high molecular size gDNA extraction, the strains were cultivated on PDA medium on permeable cellophane membranes that were prepared as described by Fauchery et al. [54]. The cellophane membranes were trimmed to the size of the Petri dish, put into boiling distilled water containing EDTA (1 g·L −1 ) 20 min in order to permeabilize the membrane, rinsed four times in a big container with distilled water, and autoclaved. After 7 days at 26 • C microorganisms were removed from cellophane membrane, snap-frozen, and milled using mortar and pestle. gDNA was extracted from microorganisms using NucleoBond Buffer Set III and AXG 100 (Macherey-Nagel, Hoerdt, France) with slight modifications: 250 mg of milled mycelium was resuspended in 5 mL of buffer G3, 37 • C for 14 h; 1.2 mL of buffer G4 was added; and the solution was gently mixed then incubated at 50 • C for 4 h. Sample was clarified through centrifugation 5000× g, 5 min. After ethanol precipitation, samples were re-dissolved in 100 µL of mQ water. Purity control was carried out on a Nanodrop 2000 spectrophotometer, DNA quantity using a Qubit dsDNA BR assay kit according to the manufacture's recommendation and a Qubit fluorometer, and quality control on 0.7% agarose gel.
3.9. Biological Tests 3.9.1. Antimicrobial Assays Pure isolated compounds were tested on human pathogenic microorganisms such as Candida albicans (ATCC 10213) to evaluate their antifungal activities. The test was performed in conformance with a reference protocol from the European Committee on Antimicrobial Susceptibility Testing (EUCAST). The MIC value was obtained after 48 h for C. albicans, and fluconazole (for fungi) was used as a positive control [58].

Antiparasitic Activities Assay for In Vitro Inhibition of T. brucei Brucei Growth
Bloodstream forms of Trypanosoma brucei brucei strain 93 were cultured in HMI9 medium supplemented with 10% FCS at 37 • C under an atmosphere of 5% CO 2 . In all experiments, log-phage cell cultures were harvested by centrifugation at 3000× g and immediately used. Drug assays were based on the conversion of a redox-sensitive dye (resazurin) to a fluorescent product by viable cells and were performed according to the manufacturer recommendations (AlamarBlue ® Assay, Invitrogen Corporation, Carlsbad, CA, USA) [59,60]. Drug stock solutions were prepared in pure DMSO. T. brucei brucei bloodstream forms (3 × 10 4 cells/mL) were cultured as described above in 96-well plates (200 µL per well) either in the absence or in the presence of different concentrations of inhibitors, with a final DMSO concentration that did not exceed 1%. After a 72 h incubation, resazurin solution was added to each well at with a final concentration of 45 µM. Fluorescence was measured at 530 nm excitation and 590 nm emission wavelengths after a further 4 h incubation. The percentage of inhibition of the parasite growth rate was calculated by comparing the fluorescence of parasites maintained in the presence of drug with that in the absence of drug. DMSO was used as control. Pentamidine was used as an anti-trypanosomal drug control (IC 50 = 0.011 ± 0.0017 µM). IC 50 s were determined from the dose-response curves with drug concentrations ranging from 100 µM to 50 nM. The IC 50 value is the mean ± the standard deviation of three independent experiments. Pentamidine was used as an anti-trypanosomal drug control (IC 50 = 0.011 ± 0.0017 µM).

In Vitro Growth Inhibition of Chloroquine-Resistant P. falciparum
The FcB1/Colombian strain of Plasmodium falciparum was maintained in vitro in human erythrocytes in an RPMI 1640 reaction mixture supplemented with 8% (v/v) heatinactivated human serum at 37 • C under an atmosphere of 3% CO 2 , 6% O 2 , and 91% N 2 . The chloroquine-resistant strain FcB1/Colombia of Plasmodium falciparum was obtained from the National Museum Natural History collection, Paris, France (n • MNHN-CEU-224-PfFCB1). Parasites were maintained in vitro on human A+ erythrocytes in RPMI 1640 medium supplemented by 8% (v/v) heat-inactivated human serum at 37 • C under an atmosphere of 3% CO 2 , 6% O 2 , and 91% N 2 . Human red blood cells and serum were provided by the Etablissement Français du Sang, Ile de France (EFS, Paris, France) under the C-CPSL-UNT approval n • 13/EFS/126. In vitro drug susceptibility assays were measured by [ 3 H]-hypoxanthine incorporation as described previously [61]. Drug solutions were prepared in DMSO at a 10 mM concentration. Compounds were serially diluted 2-fold with 100 µL culture medium in 96-well plates. Asynchronous parasite cultures (100 µL, 1% parasitemia and 1% final hematocrit) were then added to each well and incubated for 24 h at 37 • C prior to the addition of 0.5 µCi of [ 3 H]-hypoxanthine (GE Healthcare, Tremblayen-France, France, 1 to 5 Ci·mmol/mL) per well. After a further incubation of 24 h, plates were frozen and thawed. Cell lysates were then collected onto fiberglass filters and counted in a liquid scintillation spectrometer. The growth inhibition for each drug concentration was determined by comparison of the radioactivity incorporated in the treated culture with that in the control culture maintained on the same plate. The concentration causing 50% growth inhibition (IC 50 ) was obtained from the drug concentration-response curve, and the results are expressed as the mean values ± standard deviations determined from several independent experiments. Chloroquine was used as an antimalarial drug control (IC 50 = 0.072 ± 0.0074 µM).

Cytotoxic Assay
Human umbilical vein endothelial cells (HUVECs) were obtained from Clonetics (Lonza; Walkersville, MD, USA) and cultured according to the supplier's instructions. Briefly, HUVECs from three to six passages were subcultured to confluence onto 0.2% gelatin-coated tissue culture flasks in endothelial cell growth medium (EGM2) containing growth factors and 2% FCS. All cell lines were maintained at 37 • C in a humidified atmosphere containing 5% CO 2 . Cell viability was assessed using Promega (Promega France, Charbonnière-les-Bains, France) CellTiter-BlueTM reagent according to the manufacturer's instructions. Cells were seeded in 96-well plates (5 × 103 cells per well) containing 50 mL growth medium. After 24 h of culture, the cells were supplemented with 50 mL of test compound dissolved in DMSO (<0.1% in each preparation). After incubation for 72 h, 20 mL resazurin was added for 2 h before recording fluorescence (λex = 560 nm, λem = 590 nm) using a Victor microtiter plate fluorimeter (PerkinElmer France, Villebon-sur-Yvette, France). IC 50 values correspond to the concentration of test compound that elicited a 50% decrease in fluorescence for drug-treated cells relative to untreated cells. Experiments were performed in triplicate. Fumagillin was used as a positive control and showed an IC 50 of 0.006 ± 0.002 µM.

Conclusions
This study combined genomic and metabolomic approaches to study metabolome produced by mutualistic fungal strains of termites from French Guiana of Pseudallescheria genus. It allowed the annotation of tyroscherin analogs previously described, as well as pseurotin A, but also new ovalicin analogs.
Our P. boydii strain collection is a source of diverse ovalicin derivatives, including new fatty acid esters of ovalicin and boyden C that have been isolated and structurally characterized. Compounds 2-4, similar to ovalicin, exhibited antiplasmodial and antitrypanozomal properties. Our work also confirms that P. boydii strains are associated with termite biosynthesized compounds with biological activities. We were able to demonstrate that antiparasitic compounds such as Compounds 1 to 4 are produced in different amounts by the different strains of P. boydii. These same compounds reduce cytotoxicities, with IC 50 s greater than 10 µM on HUVECs.
The molecular networking analysis allowed the annotation of all isolated compounds and many other ovalicin analogues, together with the isolation and identification of the new ovalicin analog boyden C. Finally, it was possible by genome mining to propose a common biosynthetic pathway for the ovalicin analogs compounds produced by Pseudallescheria sp. strains of our collection. The biosynthesis also permitted the deduction of absolute configurations of ovalicin analogs.
Fungal natural products possess biological activities that are of great value for medicine, agriculture, and manufacturing. Recent "omics" studies accentuated the immense taxonomic diversity of fungi, and the accompanying specialized metabolic diversity provides a significant and still largely untapped resource for natural product discovery. Particular attention will be given in the future to integrative pan-genomic approaches that use a combination of genomic and metabolomic data for parallelized natural product discovery across multiple strains. Such novel technologies will not only expedite the natural product discovery process but will also allow the assembly of a high-quality toolbox for the re-design or even de novo design of biosynthetic pathways using synthetic biology approaches, while continuing to focus on specific environments such as symbiotic microorganisms.