- Research article
- Open Access
Polyketide synthases of Diaporthe helianthi and involvement of DhPKS1 in virulence on sunflower
BMC Genomics volume 19, Article number: 27 (2018)
The early phases of Diaporthe helianthi pathogenesis on sunflower are characterized by the production of phytotoxins that may play a role in host colonisation. In previous studies, phytotoxins of a polyketidic nature were isolated and purified from culture filtrates of virulent strains of D. helianthi isolated from sunflower. A highly aggressive isolate (7/96) from France contained a gene fragment of a putative nonaketide synthase (lovB) which was conserved in a virulent D. helianthi population.
In order to investigate the role of polyketide synthases in D. helianthi 7/96, a draft genome of this isolate was examined. We were able to find and phylogenetically analyse 40 genes putatively coding for polyketide synthases (PKSs). Analysis of their domains revealed that most PKS genes of D. helianthi are reducing PKSs, whereas only eight lacked reducing domains. Most of the identified PKSs have orthologs shown to be virulence factors or genetic determinants for toxin production in other pathogenic fungi. One of the genes (DhPKS1) corresponded to the previously cloned D. helianthi lovB gene fragment and clustered with a nonribosomal peptide synthetase (NRPS) -PKS hybrid/lovastatin nonaketide like A. nidulans LovB. We used DhPKS1 as a case study and carried out its disruption through Agrobacterium-mediated transformation in the isolate 7/96. D. helianthi DhPKS1 deleted mutants were less virulent to sunflower compared to the wild type, indicating a role for this gene in the pathogenesis of the fungus.
The PKS sequences analysed and reported here constitute a new genomic resource that will be useful for further research on the biology, ecology and evolution of D. helianthi and generally of fungal plant pathogens.
Polyketides are a large and diverse group of secondary metabolites with different biological activities, including pathogenicity, such as T-toxin produced by Cochliobolus heterostrophus  and melanin, a pigment essential for plant pathogenesis by many fungi [2, 3]. Biosynthesis of these metabolites is accomplished by polyketide synthases . Polyketide synthases (PKSs) are modular enzymes classified on the basis of their molecular architecture and operating mechanisms into types I, II and III [5,6,7,8]. Fungal PKs show a structural diversity that can vary from simple aromatics to highly modified complex reduced-type compounds . Many fungal PKSs have an architecture based on single modular iterative type I polyketide synthases (iPKSs), containing ketosynthase (KS), acyltransferase (AT) and acyl carrier protein (ACP) domains. In addition to these, other functional domains, such as ketoreductase (KR), dehydratase (DH), enoylreductase (ER), methyltransferase (MeT) and thioesterase (TE), may be present in PKSs .
Diaporthe helianthi Munt.-Cvetk., Mihaljč. & M. Petrov (syn. Phomopsis helianthi Munt.-Cvetk., Mihaljč. & M. Petrov) is a phytopathogenic fungus which causes stem canker and leaf shedding in sunflower (Helianthus annuus L.). The fungus invades and spreads through the leaves, progresses towards the petioles along foliar veins, and finally enters the stem, where cankers form in the advanced stage of pathogenesis . D. helianthi is an important pathogen with a worldwide distribution. First reported in the former Yugoslavia , it subsequently spread to several other countries [13,14,15]. It can cause significant losses in yield and reduction of oil content when environmental conditions are favourable for disease development . In Italy this disease has been observed since 1987, but even when climatic conditions are favourable to the disease [17, 18], its appearance is sporadic . Epidemiological differences of this widespread disease could be explained by a high genetic diversity occurring in the D. helianthi populations depending on their geographic origin. Intraspecific variability has been previously displayed by isolates representative of diverse sunflower-growing areas [20,21,22,23]. All the isolates collected in France and in the former Yugoslavia, where epidemics of sunflower stem canker are severe, formed a monophyletic clade clearly distinct from all other isolates, while all the Italian isolates were phylogenetically distant from this cluster, evidencing a clear link between genetic biotype and pathogenic behaviour .
It has been nowadays accepted that sunflower stem canker is associated with a complex of Diaporthe species with different levels of pathogenicity. The main causal agent of the disease has been identified as D. gulyae sp. nov. in association with two less virulent species, D. kochmanii sp. nov. and D. kongii sp. nov. in Australia , whereas in the United States, D. gulyae and D. helianthi were both identified as causal agents of the disease having similar levels of aggressiveness .
The highly virulent French isolate 7/96 can be referred as D. helianthi sensu stricto. A DNA region of 532 bp from this isolate, shared with all highly virulent strains but not with the Italian mildly virulent isolates, was identified and sequenced . This sequence, named lovB (accession number: AJ512137) showed a high similarity to genes encoding polyketide synthases (PKSs) from several species of filamentous fungi, including mlcA and mlcB of Penicillium citrinum, lovB of Aspergillus terreus, fum5 of Gibberella moniliformis, pks1 of Cochliobolus heterostrophus and pks1 of G. fujikuroi.
The mechanisms of pathogenicity and symptom induction are poorly understood in D. helianthi. Mazars et al. [26, 27] have demonstrated the production of a polyketidic phytotoxin, named phomozin, during pathogenesis on sunflower leaves and in culture filtrates of a French D. helianthi strain. The purified toxin produced symptoms comparable to those caused by pathogen infection. Subsequently, Avantaggiato et al.  purified other two phytotoxic metabolites, identified as cis- and trans-4,6-dihydroxymellein, from cultures of French and Italian isolates with different degrees of virulence. These toxins show a structure similar to phomozin, sharing with it the same precursor, known as orsellinic acid .
In the present study we generated a draft genome of the highly virulent isolate D. helianthi 7/96 and annotated 40 genes coding for putative PKSs . By using the Pathogen-Host Interaction database (PHI-base) (http://www.phi-base.org) , orthologs of genes known to be involved in biosynthesizing PKs, which are virulence factors in other fungal species, were identified. Furthermore, the role in virulence of the gene DhPHS1 (= D. helianthi polyketide synthase 1) was evaluated through a gene disruption approach.
Fungal isolate and media
D. helianthi highly virulent French isolate 7/96 belongs to the fungal collection of Department of Agriculture, Food and Environment, University of di Pisa (Italy)  and was maintained on slants of PDA (potato dextrose agar, Difco) under mineral oil at 4 °C.
Diaporthe helianthi strain 7/96 draft genome sequence v2
Based on raw data available from a previous project , we generated a new D. helianthi strain 7/96 genome assembly. Paired end reads of 90 bp (1.80 Gbp) were assembled using SPAdes 3.11.0 . The genome of D. helianthi consists of 7376 sequence scaffolds with a total assembly length of 63.67 Mbp (N50 = 20,184 and L50 = 860), 43.99% GC-content, and a maximum scaffold size of 151,286 bp. The completeness of the assembly was assessed using BUSCO v12 , which estimated the genome sequence to be 99.65% complete. The genome was annotated using the MAKER2 pipeline . Overall, 13,139 protein-coding gene models were predicted.
The new genome assembly of D. helianthi strain 7/96 is present in GenBank with accession number: MAVT00000000.2.
Genomic characterization of putative PKSs genes
Putative PKS genes were identified according to Klarsson et al.  and manually inspected for conserved domain (acyl transferase [AT] – InterPro domain IPR014043, acyl carrier protein or phosphopantetheine attachment site [ACP or PP] – IPR009081, beta-ketoacyl synthase N-terminal domain [KS-N] – IPR014030, beta-ketoacyl synthase C-terminal domain [KS-C] – IPR014031, Ketoreductase [KR] – IPR013968, polyketide synthase dehydratase [DH] – IPR020807, polyketide synthase, enoylreductase domain [ER] – IPR020843, methyltransferase [MT] – IPR013217, thioesterase [TE] – IPR001031) using InterProScan . Moreover, putative D. helianthi PKSs sequences were aligned with references of other ascomycetes using MAFFT 7.310  and a phylogenetic analysis was performed with PhyML 3.0 .
DNA molecular techniques
Total DNA was obtained from D. helianthi grown on PDA plates overlaid with a cellophane membrane. Plates were inoculated with mycelial plugs and incubated at 24 °C for 5 days. Mycelium mats were peeled from membranes, freeze-dried overnight, and used for DNA extraction according to Raeder and Broda .
Plasmids were purified from E. coli DH5α cultures, grown on LB (yeast extract 5 g/L, NaCl 5 g/L, tryptone 10 g/L) by using QIAprep Spin Miniprep Kit (QIAGEN).
All PCR amplifications were performed as follows: initial denaturation at 94 °C for 3 min, 35 cycles of denaturation (45 s at 94 °C), annealing (45 s at 60 °C) and extension (1 min 50 s at 72 °C), and a final extension at 72 °C for 10 min. PCR reactions were carried out in PCR buffer (Promega, Madison, WI, USA), 0.2 mM dNTPs (Roche Applied Science, Mannheim, Germany), 0.2 M each primer, 0.4 unit of Taq polymerase (Promega), and ca. 5 ng of template DNA.
To amplify the complete sequence of DhPKS1 gene, PCR were performed by using Platinum® Taq DNA Polymerase High Fidelity (ThermoFisher) with the primers DhPKS1for (ATGTCCAAGGCAATTTGTACTAC) and DhPKS1rev (CCTATCGCTAACAATCTTGT). The PCR cycles were as follows: initial denaturation at 95 °C for 5 min, 35 cycles of denaturation (45 s at 95 °C), annealing (45 s at 60 °C) and extension (8 min at 68 °C), and a final extension at 68 °C for 10 min.
Sequence editing and analysis were carried out using BioEdit 184.108.40.206 (http://www.mbio.ncsu.edu/bioedit/bioedit.html) software and online tools available at European Bioinformatics Institute (http://www.ebi.ac.uk/).
DNA restriction, elution and ligation, and Southern blot analysis were carried out as described by , performing hybridization at 60 °C for 16 h and washing of blots in 2× SSC and 0.5 SSC at room temperature.
Construction of plasmid pUR5750∙ΔDhPKS1
For D. helianthi transformation, the binary vector pUR5750, described by De Groot et al.  was used. This vector contains the neomycin phosphotransferase gene under the control of the nopaline synthase promoter and trpC terminator that confers resistance to kanamycin; it also carries, between the HindIII and KpnI restriction sites, the E. coli hygromycin B (hph)-resistance cassette, coding for the hygromycin B phosphotransferase enzyme (hph) under the control of the A. nidulans gpdA promoter and trpC terminator. The two flanking regions of the previously isolated sequence LovB  were identified and cloned with GenomeWalker kit (Clontech Laboratories, Palo Alto, CA) as follows: Separate fungal DNA aliquots were digested with four different restriction enzymes (EcoRV, DraI, PvuII, StuI) leaving blunt ends and ligated to adaptors. For each fragment library, two primary PCR amplifications were carried out using an adaptor primer provided with the kit and an outer, gene-specific primer for downstream and upstream walking, 5′-AAG GTG GAC ACG GCA TAC CAC TCA TT-3′ and 5′-CCA AGT CTT CAG CAG GAA TAT CAA CCA C-3′, respectively. The primary PCR product was then diluted and used as a template for a secondary PCR amplification using a nested adaptor primer and nested gene-specific primers (5′-AGC TGC AAG TGC CTT ACC ACG GAT TAC-3′ for downstream walking and 5′-ATG AGT GGT ATG CCG TGT CCA CCT TC-3′ for upstream walking). The resulting DNA, flanking LovB, of 548 and 1256 bp were singly cloned in p-GEM-T Easy vector system (Promega) to form clones and DL2, respectively. To generate pUR5750∙ΔDhPKS1, the two previously cloned regions were excised from p-GEM-T Easy vector and fragments of about 548 bp (5′ flanking gene fragment, called hereafter “A”) and 944 bp (3′ flanking gene fragment, called hereafter “B”), were inserted in pUR5750 KpnI and HindIII restriction sites, respectively, at the sides of (hph)-resistance cassette (Fig. 1).
In detail, the A fragment was excised from the plasmid by EcoRI digestion, and PCR-amplified with specific primers (5′AGGTACCATTCGATTACTATAGGGCACG 3′ and 5′AGGTACC GTACTCAGGCATGGAGCAAA 3′), carrying at ends the KpnI cutting site sequence. The B fragment, was excised from the plasmid pGEM-T easy by EcoRI digestion, and inserted in the EcoRI site of pBluescript KS; from this construct a smaller fragment of about 944 bp was excised, with HindIII restricion enzyme. The two DhPKS1 gene fragments (KpnI 548 bp and HindIII 944 bp) flanking the designed knock-out site, were inserted in the corresponding cutting-site in the plasmid pUR5750, upstream and downstream of the hph resistance cassette, to form a new plasmid named pUR5750∙ΔDhpks1. Finally, pUR5750∙ΔDhPKS1 was transferred into A. tumefaciens LBA1100 by electroporation ; electroporation conditions were 25 μF, 200 Ω, 2.5 kV (0.2 cm cuvettes) in Gene Pulser® electroporator (Bio-Rad, USA), and transformants were selected on LB agar with and kanamycin (100 μg/ml). Agrobacterium strains containing the binary vectors were identified by PCR.
Agrobacterium tumefaciens-mediated gene disruption
The LBA1100-derived strain transformed with pUR5750∙ΔDhPKS1 was grown at 28 °C for 3 days in Petri dishes containing LB medium supplemented with kanamycin (100 μg/ml). One hundred milliliter of liquid LB supplemented with kanamycin (100 μg ml−1) were inoculated with a single bacterial colony and incubated at 28 °C overnight on an orbital shaker at 150 rpm. Bacterial cells were harvested by centrifugation at 12000 g at 4 °C and re-suspended in 5 ml of IM (1 mM KH2PO4 pH 4.8; 2.4 mM MgSO4•7H2O; 5 mM NaCl; 0.068 mM CaCl2; 0.003 mM FeSO4•7H2O; 0.0015 mM ZnSO4•7H2O; 0.002 mM CuSO4 5H2O; 0.008 mM H3BO3; 0.003 mM MnSO4•H2O; 0.002 mM Na2MoO4•2H2O; 6.25 mM NH4NO3; 54.2 mM glycerol; 40 mM 2-(N-morpholino)ethanesulfonic acid; 0.05 mM glucose). The bacterial suspension was distributed in aliquots of 1 ml and further incubated at 30 °C for 6 h under stirring. Cultures of the fungal isolate 7/96b were grown on IM-agar supplemented with acetosyringone (AS, 0.2 mM) and incubated at 24 °C. After 4–5 days, mycelial plugs of these colonies were used to inoculate 10 ml of liquid IM + AS. The cultures were then incubated for 4 days at 24 °C under gentle stirring and periodically vortexed for 5 min. The mycelium was recovered by centrifugation at 4000 g for 15 min and re-suspended in 5 ml of IM + AS.
For co-cultivation, 1 ml of mycelial suspension was mixed with an equal volume of bacterial suspension at 24 °C and stirred for about 20 min. The mix was plated onto nitrocellulose filters on a co-cultivation medium (IM-agar + AS) and incubated at 24 °C for 3 days. After growth on co-cultivation medium, the filters were transferred to PDA amended with hygromycin (25 μg/ml) as selective medium for fungal transformants and cefotaxime (200 μM) to inhibit growth of A. tumefaciens, and incubated at 24 °C for 20 days. Hygromicin resistant colonies were purified by three successive hyphal transfers on PDA supplemented with 100 μg/ml hygromicin to select hyphae containing only transformed nuclei. The analysis of transformants was performed by PCR using specific primers designed to amplify both hygromicin resistance cassette (hph2for: 5′-ATGGCAACAAATGTTGGACTG-3′; gpdArev: 5′-CAAGGAGGAGTAAGCTCCTT-3′), and DhPKS1 gene replacement site (Dhpks537for: 5′-GTACTCCAAGGCTTTATCGC-3′; Dhpks1326rev: 5′-TGATGTAGAACTGGGCCACA-3′) of D. helianthi. As control, D. helianthi mutants containing pUR5750 empty vector insertion were obtained.
Single copy insertion was verified by Southern blot analysis of genomic DNA digested with SpeI and BglII restriction enzymes, which do not cut in the Agrobacterium-transferred construct. Hybridization was carried out with digoxigenin-labelled probe obtained by PCR-amplification of binary vector pUR5750ΔDhPKS1 plasmidic DNA using hph2for - gpdArev primers. D. helianthi wild type DNA and PCR-amplified construct ΔDhPKS1 containing the hygromycin B (hph)-resistance cassette, were used as negative and positive controls, respectively.
Phytotoxicity and virulence assays
For phytotoxicity tests, mycelial plugs from actively growing colonies of D. helianthi isolate 7/96 and transformants were grown, in static condition, on liquid medium in 2 L volume Erlenmeyer flasks containing 1 L of substrate made with: 3 g/L L-asparagine, 15 g/L sucrose, 1 g/L K2HPO4, 0.5 g/L MgSO47 H2O, 0.5 g/L KCl, 0.018 g/L FeSO47H2O, 5 g/L fresh sunflower tissue. The same liquid medium without fungal inoculation was used as a control. After 28 days, the liquid cultures were vacuum filtered, sterilised through Millipore filters of 0.2 μm and used in phytotoxicity tests, performed with cuttings of 15-day-old sunflower seedlings according to Avantaggiato et al. . Ten sunflower cuttings were used per experiment and the experiment was replicated three times.
For virulence tests inoculum was prepared by growing D. helianthi strain 7/96 and transformants in Petri dishes containing different parts of sunflower (stems, leaves and flowers). Fresh plant tissues, surface sterilized with sodium hypochlorite (2.5%/vol.) and washed in sterile water, were dried, chopped (250 g/L) and mixed with agar (15 g/L). The fungus inoculated on the sunflower-agar substrate was left to grow in the dark for 1 week at 25 °C. Seventy-eight-day-old sunflower plants (very susceptible cv Ala) at the flower bud stage were inoculated by placing a mycelial plug both on the upper surface of leaves and at the insertion point of the leaf petiole (sticking with a wooden toothpick). The diseased area was measured recording two diameters of the necrotic spot, by using a digital caliper. Inoculated plants were kept in a moist chamber for 7 days at 24 °C. The bioassay was carried out on five plants for each fungal strain and was repeated twice. Data were analysed by ANOVA test with Welsch modification, because data were normally distributed but not homoscedastic. Different means were separated by T3 Dunnett test. Statistical analyses were performed by SPSS 20.0 software.
PKS genes in Diaporthe helianthi genome
By searching the draft nuclear genome of D. helianthi isolate 7/96 , we found a large number of putative PKS homologues which have been deposited in GenBank (Additional file 1). Most PKS genes coded by D. helianthi were reducing PKSs, whereas only eight PKSs lacked reducing domains and clustered with non-reducing PKSs (Additional file 1). DhPKS8, DhPKS18, DhPKS27 proteins did not cluster into any of the clades indicated in Fig. 2. DhPKS1 was placed in nonribosomal peptide synthetase (NRPS)-PKS clade, with the closest sequence being ATEG_00325 protein from A. terreus involved in isoflavipucine biosynthesis.
Homologies of DhPKSs with other genes involved in host-pathogen interaction
By whole genome sequencing of D. helianthi, we were able to decipher the complete sequence of all DhPKS genes and in silico deduced amino acid sequences (Additional file 1). Through a PHI-base interrogation we identified homologs for all DhPKS genes, except for two (DhPKS11 and DhPKS35), with experimentally verified pathogenicity, virulence and effector genes from fungal, oomycete and bacterial pathogens, infecting animal, plant, fungal and insect hosts. Results of this search are reported in Table 1.
Replacement of DhPKS1
Using oligonucleotides designed at 5′ and 3′ ends of genomic sequence of DhPKS1 we sequenced the complete gene (7877 bp) confirming its in silico prediction. In silico analysis showed that the coded protein corresponded to a highly reducing (HR) type I iPKS, containing the full set of domains, ketoacyl synthase (KS), acyl transferase, (AT), ketoreductase (KR), dehydratase (DH) and acyl carrier protein (ACP) (Additional file 1). Due to the identity of DhPKS1 with lovB (AN AJ512137) fragment identified by Vergara et al.  we decided to test its potential involvement in pathogenicity of D. helianthi.
Transformation of D. helianthi 7/96 with the gene disruption cassette in DhPKS1 yielded 160 hygromicin-resistant colonies. These putative transformants were purified by successive transfers of mono-hyphal and single protoplasts to selective medium amended with hygromycin 100 μg/L. Hygromycin resistance of transformants was not lost upon successive transfers; 74 of the hygromicin resistant colonies were screened by two PCR analyses. The first PCR was performed on the genomic DNA of each transformant, with the hph2for – trpc2rev primers specific to amplify the hygromicin resistance cassette (Fig. 3). Hph-containing transformants were then screened with a second PCR, to verify the homologous integration of the construct. In this case, we used Dhpks537for - Dhpks1326rev primers, designed on sequences internal to the replacement site of DhPKS1 (Fig. 4). This second PCR showed that 13 out of the 74 (18%) transformants integrated the ΔDhPKS1 construct at the homologous site. The number of integrated copies of the construct in the genome of D. helianthi transformants was ascertained by Southern blot analysis with a probe specific for the hygromicin-resistance cassette (Fig. 5). Single copy integration of the construct occurred in five transformants designated as Tr1-5, as the enzymes used for genomic DNA digestion do not cut the Agrobacterium-transferred construct.
Phytotoxicity and virulence bioassays
ΔDhPKS1 mutants showed regular in vitro development and their growth rates were not significantly different from that of the wild type. Sunflower cuttings steeped in liquid cultures of the D. helianthi isolate 7/96 showed, after 5 days, brown marginal necrotic lesions and leaf chlorosis attributable to filtrate phytotoxicity. Cuttings steeped in culture extracts of DhPKS1 knock-out mutants, did not show clear evidence of phytotoxicity (Fig. 6).
Inoculation of sunflower plants with D. helianthi strain 7/96 wild type resulted in typical symptoms of stem canker disease when observed 20 days after inoculation (Fig. 7a). Plants inoculated with transformants showed less intense symptoms (Fig. 7b). The disease severity in terms of necrotic area was measured in cm2 on leaf blade and stem, in sunflower plants inoculated with the parental strain and transformants. For each time point, infection was markedly less severe with ΔDhPKS1 mutant than with the other strains. In particular the average necrotic area (cm2) at 84 h was 1.3 ± 0.28, 8.0 ± 1.65 and 8.5 ± 1.56, for ΔDhPKS1, wild type and empty vector transformant strains, respectively (Fig. 7c).
The management of sunflower stem canker requires the identification of diverse species of Diaporthe associated with it, the determination of their aggressiveness and the study of pathogenicity mechanisms. In this paper we analysed the genome of the highly aggressive isolate 7/96 of D. helianthi, regarded as the main causal agent of sunflower stem canker epidemics in Europe , on the basis of its putative polyketide synthase genes. We found in its genome at least 40 PKS genes, more than those reported in other ascomycetes: 15 in G. moniliformis, 16 in G. zeae, 20 in B. cinerea, 25 in C. heterostrophus  and 27 in Aspergillus nidulans . Such a repertoire of PKSs makes possible the synthesis of almost all known types of polyketide compounds by D. helianthi 7/96.
The PHI-base search demonstrated that all the putative 40 PKS proteins found in D. helianthi, with the exception of DhPKS11 and DhPKS35, have at least one ortholog in other pathogenic fungi or bacteria (Table 1). The majority of the reported orthologous genes have been demonstrated to be important factors in pathogenicity or virulence by gene deletion experiments.
A fragment of LovB-like coding sequence was isolated by Vergara et al.  from the genome of the aggressive isolate of D. helianthi 7/96, and proved to be present in the genome of the most aggressive isolates of D. helianthi, such as the French and Yugoslavian ones. In the present work this gene, named DhPKS1, was fully characterized. In silico sequence analysis revealed that DhPKS1 belongs to subclade II of type I reducing PKS , having the typical conserved domain of this protein family KS-AT-DH-(ME)KR-PP-(CON)-(AMP-PP). Reducing PKS subclade II is characterized by enzymes missing the ER domain; PKs synthesized by PKSs of this subclade are predicted to either lack reduced alkyl groups or to contain alkyl groups whose reduction is completed by the product of an external ER domain-containing gene, as A. terreus lovC  and P. citrinum mlcG. The PKSs of this clade were also found to have either a condensation (CON) domain typical of nonribosomal peptide synthetases (NPSs)  or an entire NPS module consisting of a CON domain, and an acyl carrier or phosphopantetheine attachment site (ACP or PP) domain. The absence of a DhPKS1 homologous gene in the less aggressive Italian isolates prompted us to further investigate the possible significance of this gene as a virulence determinant. Results from PHI-base search revealed that DHPKS1 is very similar to PKSs from filamentous fungi known to be involved in the production of important pathogenicity factors. Among these, it showed 35% identity with BcPKS6 gene of B. cinerea which encodes a key enzyme for botcinic acid biosynthesis, a phytotoxin involved in virulence of the fungus on tomato . Moreover, BcPKS6 belongs to a PKS cluster co-regulated by the Gα subunit BCG1, which, in turn, is essential for pathogenicity of B. cinerea on bean leaves  and is part of a clade including several PKSs responsible for the synthesis of cyclic polyketides like: LovB (LNKS, nonaketide part of lovastatin ); MlcA (nonaketide part of citrinin; ), EQS (equisetin; ) and FusS (fusarin; ). DHPKS1 has also a relatively high homology (Table 1) with the Magnaporthe grisea gene ACE1 encoding a putative PKS expressed exclusively during penetration of M. grisea into leaves and involved in the recognition of the fungus by resistant rice varieties carrying the resistance gene Pi33, thus revealing a role for avirulence [50, 51]. Orthologous to DhPKS1 is also AbPKS9 (DEP5) gene involved in biosynthesis of depudecin, an 11 linear polyketide inhibitor of histone deacetylase (HDAC) virulence factor of the fungus A. brassicicola .
Among the functionally characterized PKSs DhPKS1 orthologs, there is also lovB of A. terreus, which synthesizes the nonaketide chain of lovastatin, a 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase inhibitor . A. terreus has been known to carry an unusual PKS gene cluster for the polyketides in which two PKS genes (lovB and lovF) are closely linked in the cluster and are required for the biosynthesis of the nonaketide and the diketide moieties of the compound, respectively . Interestingly, we also found that LovF of A. terreus was homologous to DhPKS24, suggesting the presence in D. helianthi 7/96 of a similar PKS gene cluster. The same hypothesis has been formulated for the biosynthesis of both zearalenone in Gibberella zeae  and the linear polyketide T-toxin in race T of C. heterostrophus. The synthesis of T-toxin, essential for fungal virulence on male sterile cytoplasm corn , seems to involve a PKS1 together with PKS2 .
In our work, as a case study, knock-out mutants of DhPKS1 were produced by targeted DNA integration through Agrobacterium-mediated transformation. This method has been successfully applied to study gene functions in other phytopathogenic and toxigenic fungi [2, 56]. The presence of a single copy of the T-DNA made the analysis of the transformants straightforward. Unfortunately, it was not possible to confirm the above results by performing experiments with complementation mutants. The reason why no double mutants were obtained is unknown, but it seems that this fungal strain cannot be doubly transformed. As a result of DhPKS1 gene inactivation, reduced fungal virulence on a susceptible sunflower cultivar was observed. Indeed, in contrast with the wild type, mutants caused less visible symptoms after artificial inoculation of fungal mycelia on stem and leaves. This indicates that DhPKS1 could be involved in virulence of the sunflower stem canker agent. The DhPKS1 gene disruption also seems to affect in vitro production of toxic secondary metabolites by D. helianthi. Liquid filtrates οf ΔDhPKS1 mutants showed a reduced phytotoxicity on sunflower seedlings. The highly reducing synthase DhPKS1 probably produces a metabolite, not yet identified, which could be involved in the D. helianthi 7/96 strain virulence on sunflower. However, for a better understanding of the possible role of DhPKS1, more studies should be carried out considering also the role of the other PKSs genes that cluster with it.
In previous work  the polyketidic metabolite phomozin, an ester of orsellinic acid, was isolated both from culture filtrates of D. helianthi and from infected plants, and its possible role in development of symptoms was suggested. Orsellinic acid synthase is the simplest tetraketide synthase and is grouped with the nonreducing PKSs (NR-PKSs). Our results indicate that, in the case of D. helianthi 7/96 orsellinc acid synthase could be coded by DhPKS7 due to its high homology with the A. nidulans gene EAA59563 coding for orsellinc acid synthase in A. nidulans . The toxinic theory is supported by the evidence that other phytopathogenic species related to the genus Phomopsis produce toxic metabolites involved in pathogenesis [58,59,60]. For a disease caused by a Phomopsis species affecting soybean, similar to stem canker of sunflower, the possible involvement of a phytotoxin was also implied . Polyketide synthases have been reported as important virulence factors in other several phytopathogenic fungi such as C. heterostrophus , Colletotrichum lageniarum  and A. alternata . In C. nicotianae, the genes CTB1 and CTB3 encode two polyketide synthases, involved in the biosynthesis of cercosporin, a photoactivated perylenequinone toxin, which play a key role in fungal pathogenesis [64, 65].
The PKS sequences reported here are a new important resource that will be useful for further research in the biology, ecology and evolution of D. helianthi and in general of fungal plant pathogens. Further investigation is necessary to fully understand the role of D. helianthi PKS genes. Moreover, it will be very important to verify if the PKSs identified in D. helianthi sensu stricto are also present in other species of Diaporthe, which have been associated with sunflower stem canker.
Acyl carrier protein
3-hydroxy-3-methylglutaryl coenzyme A
Nonaketide part of lovastatin
Nonribosomal peptide synthetase
Baker SE, Kroken S, Inderbitzin P, Asvarak T, Li BY, Shi L, Yoder OC, Turgeon BG. Two polyketide synthase-encoding genes are required for biosynthesis of the polyketide virulence factor, T-toxin, by Cochliobolus heterostrophus. Mol Plant-Microbe Interact. 2006;19(2):139–49.
Zhang A, Lu P, Dahl-Roshak AM, Paress PS, Kennedy S, Tkacz JS, An Z. Efficient disruption of a polyketide synthase gene (pks1) required for melanin synthesis through Agrobacterium-mediated transformation of Glarea lozoyeasis. Mol Gen Genomics. 2003;268(5):645–55.
Yang G, Rose MS, Turgeon BG, Yoder OC. A polyketide synthase is required for fungal virulence and production of the polyketide T-toxin. Plant Cell. 1996;8(11):2139–50.
Khosla C, Gokhale RS, Jacobsen JR, Cane DE. Tolerance and specificity of polyketide synthases. Annu Rev Biochem. 1999;68:219–53.
Shen B. Polyketide biosynthesis beyond the type I, II and III polyketide synthase paradigms. Curr Opin Chem Biol. 2003;7(2):285–95.
Crawford JM, Townsend CA. New insights into the formation of fungal aromatic polyketides. Nat Rev Microbiol. 2010;8(12):879–89.
Simpson TJ. Fungal polyketide biosynthesis - a personal perspective. Nat Prod Rep. 2014;31(10):1247–52.
Vederas JC. Explorations of fungal biosynthesis of reduced polyketides - a personal viewpoint. Nat Prod Rep. 2014;31(10):1253–9.
Fujii I. Functional analysis of fungal polyketide biosynthesis genes. J Antibiot (Tokyo). 2010;63(5):207–18.
Fujii I, Watanabe A, Ebizuka Y. More functions for multifunctional Polyketide Synthases. In: Tkacz JS, Lange L, editors. Advances in fungal biotechnology for industry, agriculture, and medicine. Boston: Springer US; 2004. p. 97–125.
Muntanola-Cvetkovic M, Vukojevic J, Mihaljcevic M. The systemic nature of the sunfluower disease caused by Diaporthe helianthi. Can J Bot. 1991;69(7):1552–6.
Acimovic M, Straser N. Phomopsis sp. - novi parazit suncokreta. (Phomopsis sp. a new sunflower parasite). Zastita Bilja. 1982;33(2)
Herr LJ, Lipps PE, Watters BL. Diaporthe stem canker of sunflower. Plant Dis. 1983;67:911–3.
Delos M, Moinard J. Évolution du Phomopsis du tournesol en France: un bref historique. Phytoma. 1995;473:22–4.
Gulya T, Rashid KY, Masirevic SM. Sunflower diseases. In: AAE S, editor. Sunflower technology and production. Madison: ASA-CSSA-SSSA; 1997. p. 263–379.
Diaz Franco A, Ortegon Morales A. Influence of sunflower stem canker (Diaporthe helianthi) on seed quality and yield during seed development. Agronomie. 2003;23(7):581–92. https://doi.org/10.1051/agro:2003032.
Zazzerini A, Tosi L, Losavio N. Rilievi fitopatoligici su varietà di girasole a confronto nel 1987. L’informatore Agrario. 1988;13:85–88.
Zazzerini A, Tosi L. Situazione delle fitopatie del girasole in Italia (1970-88). Informatore Fitopatologico. 1990;2:25–8.
Battilani P, Rossi V, Girometta B, Delos M, Rouzet J, André N, Esposito S. Estimating the potential development of Diaporthe helianthi epidemics in Italy. EPPO Bull. 2003;33:427–31.
Pecchia S, Mercatelli E, Vannacci G. Intraspecific diversity within Diaporthe helianthi: evidence from rDNA intergenic spacer (IGS) sequence analysis. Mycopathologia. 2004;157(3):317–26.
Rekab D, Del Sorbo G, Reggio C, Zoina A, Firrao G. Polymorphisms in nuclear rDNA and mtDNA reveal the polyphyletic nature of isolates of Phomopsis pathogenic to sunflower and a tight monophyletic clade of defined geographic origin. Mycol Res. 2004;108:393–402.
Vergara M, Capasso T, Gobbi E, Vannacci G. Plasmid distribution in European Diaporthe helianthi isolates. Mycopathologia. 2005;159(4):591–9.
Thompson SM, Tan YP, Young AJ, Neate SM, Aitken EAB, Shivas RG. Stem cankers on sunflower (Helianthus annuus) in Australia reveal a complex of pathogenic Diaporthe (Phomopsis) species. Persoonia. 2011;27:80–9.
Mathew FM, Alananbeh KM, Jordahl JG, Meyer SM, Castlebury LA, Gulya TJ, Markell SG. Phomopsis stem canker: a reemerging threat to sunflower (Helianthus annuus) in the United States. Phytopathology. 2015;105(7):990–7.
Vergara M, Cristani C, Regis C, Vannacci G. A coding region in Diaporthe helianthi reveals genetic variability among isolates of different geographic origin. Mycopathologia. 2004;158(1):123–30.
Mazars C, Rossignol M, Auriol P, Klaebe A. Phomozin, a phytotoxin from Phomopsis helianthi, the causal agent of stem canker of sunflower. Phytochemistry. 1990;29(11):3441–4.
Mazars C, Canivenc E, Rossignol M, Auriol P. Production of phomozin in sunflower following artificial inoculation with Phomopsis helianthi. Plant Sci. 1991;75(2):155–60.
Avantaggiato G, Solfrizzo M, Tosi L, Zazzerini A, Fanizzi FP, Visconti A. Isolation and characterization of phytotoxic compounds produced by Phomopsis helianthi. Nat Toxins. 1999;7(3):119–27.
Baroncelli R, Scala F, Vergara M, Thon MR, Ruocco M. Draft whole-genome sequence of the Diaporthe helianthi 7/96 strain, causal agent of sunflower stem canker. Genomics Data. 2016;10:151–2.
Urban M, Cuzick A, Rutherford K, Irvine A, Pedro H, Pant R, Sadanadan V, Khamari L, Billal S, Mohanty S, et al. PHI-base: a new interface and further additions for the multi-species pathogen-host interactions database. Nucleic Acids Res. 2017;45:D604–10.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.
Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12:491.
Karlsson M, Durling MB, Choi J, Kosawang C, Lackner G, Tzelepis GD, Nygren K, Dubey MK, Kamou N, Levasseur A, et al. Insights on the evolution of mycoparasitism from the genome of Clonostachys rosea. Genome Biol Evol. 2015;7(2):465–80.
Goujon M, McWilliam H, Li WZ, Valentin F, Squizzato S, Paern J, Lopez R. A new bioinformatics analysis tools framework at EMBL-EBI. Nucleic Acids Res. 2010;38:W695–9.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2017;1–7. https://doi.org/10.1093/bib/bbx108.
Raeder U, Broda P. Rapid preparation of DNA from filamentous fungi. Lett Appl Microbiol. 1985;1:17–20.
Pane C, Rekab D, Firrao G, Ruocco M, Scala F. A novel gene coding for an ABC transporter in Botrytis cinerea (Botryotinia fuckeliana) is involved in resistance to H2O2. Journal of Plant Pathology. 2008;90:(3):453–462.
de Groot MJ, Bundock P, Hooykaas PJ, Beijersbergen AG. Agrobacterium tumefaciens-mediated transformation of filamentous fungi. Nat Biotechnol. 1998;16(9):839–42.
Mozo T, Hooykaas PJJ. Electroporation of megaplasmids into Agrobacterium. Plant Mol Biol. 1991;16(5):917–8.
Kroken S, Glass NL, Taylor JW, Yoder OC, Turgeon BG. Phylogenomic analysis of type I polyketide synthase genes in pathogenic and saprobic ascomycetes. Proc Natl Acad Sci U S A. 2003;100(26):15670–5.
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422(6934):859–68.
Kennedy J, Auclair K, Kendrew SG, Park C, Vederas JC, Hutchinson CR. Modulation of polyketide synthase activity by accessory proteins during lovastatin biosynthesis. Science. 1999;284(5418):1368–72.
Dalmais B, Schumacher J, Moraga J, LEP P, Tudzynski B, Collado IG, Viaud M. The Botrytis cinerea phytotoxin botcinic acid requires two polyketide synthases for production and has a redundant role in virulence with botrydial. Mol Plant Pathol. 2011;12(6):564–79.
Schumacher J, Viaud M, Simon A, Tudzynski B. The Gα subunit BCG1, the phospholipase C (BcPLC1) and the calcineurin phosphatise co-ordinately regulate gene expression in the grey mould fungus Botrytis cinerea. Mol Microbiol. 2008;67:1027–50.
Sutherland A, Auclair K, Vederas JC. Recent advances in the biosynthetic studies of lovastatin. Curr Opin Drug Discov Dev. 2001;4(2):229–36.
Abe Y, Suzuki T, Ono C, Iwamoto K, Hosobuchi M, Yoshikawa H. Molecular cloning and characterization of an ML-236B (compactin) biosynthetic gene cluster in Penicillium citrinum. Mol Gen Genomics. 2002;267(5):636–46.
Sims JW, Fillmore JP, Warner DD, Schmidt EW. Equisetin biosynthesis in Fusarium heterosporum. Chem Commun (Camb). 2005;2:186–8.
Song Z, Cox RJ, Lazarus CM, Simpson TT. Fusarin C biosynthesis in Fusarium moniliforme and Fusarium venenatum. Chembiochem. 2004;5(9):1196–203.
Bohnert HU, Fudal I, Dioh W, Tharreau D, Notteghem JL, Lebrun MH. A putative polyketide synthase/peptide synthetase from Magnaporthe grisea signals pathogen attack to resistant rice. Plant Cell. 2004;16(9):2499–513.
Fudal I, Collemare J, Bohnert HU, Melayah D, Lebrun MH. Expression of Magnaporthe grisea avirulence gene ACE1 is connected to the initiation of appressorium-mediated penetration. Eukaryot Cell. 2007;6(3):546–54.
Wight WD, Kim KH, Lawrence CB, Walton JD. Biosynthesis and role in virulence of the histone deacetylase inhibitor depudecin from Alternaria brassicicola. Mol Plant-Microbe Interact. 2009;22(10):1258–67.
Hendrickson L, Davis CR, Roach C, Nguyen DK, Aldrich T, McAda PC, Reeves CD. Lovastatin biosynthesis in Aspergillus terreus: characterization of blocked mutants, enzyme activities and a multifunctional polyketide synthase gene. Chem Biol. 1999;6(7):429–39.
Ma SM, Tang Y. Biochemical characterization of the minimal polyketide synthase domains in the lovastatin nonaketide synthase LovB. FEBS J. 2007;274(11):2854–64.
Kim YT, Lee YR, Jin J, Han KH, Kim H, Kim JC, Lee T, Yun SH, Lee YW. Two different polyketide synthase genes are required for synthesis of zearalenone in Gibberella zeae. Mol Microbiol. 2005;58(4):1102–13.
Schumann J, Hertweck C. Advances in cloning, functional analysis and heterologous expression of fungal polyketide synthase genes. J Biotechnol. 2006;124(4):690–703.
Schroeckh V, Scherlach K, Nutzmann HW, Shelest E, Schmidt-Heck W, Schuemann J, Martin K, Hertweck C, Brakhage AA. Intimate bacterial-fungal interaction triggers biosynthesis of archetypal polyketides in Aspergillus nidulans. Proc Natl Acad Sci U S A. 2009;106(34):14558–63.
Evidente A, Rodeva R, Andolfi A, Stoyanova Z, Perrone C, Motta A. Phytotoxic polyketides produced by Phomopsis foeniculi, a strain isolated from diseased Bulgarian fennel. Eur J Plant Pathol. 2011;130(2):173–82.
Patwardhan SA, Pandey RC, Dev S, Pendse GS. Toxic cytochalasins of Phomopsis paspalli, a pathogen of kodo millet. Phytochemistry. 1974;13:1985–8.
Wheeler MM, Wheeler DMS, Peterson GW. Anthraquinone pigments from the phytopathogen Phomopsis juniperovora hahn. Phytochemistry. 1974;14:288–9.
Lalitha B, Snow JP, Berggren GT. Phytotoxin production by Diaporthe phaseolorum var. caulivora, the causal organism of stem canker of soybean. Phytopathology. 1989;79:499–504.
Takano Y, Kubo Y, Shimizu K, Mise K, Okuno T, Furusawa I. Structural analysis of PKS1, a polyketide synthase gene involved in melanin biosynthesis in Colletotrichum lagenarium. Mol Gen Genet. 1995;249(2):162–7.
Izumi Y, Ohtani K, Miyamoto Y, Masunaka A, Fukumoto T, Gomi K, Tada Y, Ichimura K, Peever TL, Akimitsu K. A polyketide synthase gene, ACRTS2, is responsible for biosynthesis of host-selective ACR-toxin in the rough lemon pathotype of Alternaria alternata. Mol Plant-Microbe Interact. 2012;25(11):1419–29.
Choquer M, Dekkers KL, Chen HQ, Cao L, Ueng PP, Daub ME, Chung KR. The CTB1 gene encoding a fungal polyketide synthase is required for cercosporin biosynthesis and fungal virulence of Cercospora nicotianae. Mol Plant-Microbe Interact. 2005;18(5):468–76.
Dekkers KL, You BJ, Gowda VS, Liao HL, Lee MH, Bau HJ, Ueng PP, Chung KR. The Cercospora nicotianae gene encoding dual O-methyltransferase and FAD-dependent monooxygenase domains mediates cercosporin toxin biosynthesis. Fungal Genet Biol. 2007;44(5):444–54.
Desjardins AE, Munkvold GP, Plattner RD, Proctor RH. FUM1--a gene required for fumonisin biosynthesis but not for maize ear rot and ear infection by Gibberella moniliformis in field tests. Mol Plant-Microbe Interact. 2002;15(11):1157–64.
Zhao X, Wang Y, Zhao Y, Huang Y, Zhang KQ, Yang J. Malate synthase gene AoMls in the nematode-trapping fungus Arthrobotrys oligospora contributes to conidiation, trap formation, and pathogenicity. Appl Microbiol Biotechnol. 2014;98(6):2555–63.
Brown DWL, Lee S-H, Kim L-H, Ryu J-G, Lee S, Seo Y, Kim YH, Busman M, Yun S-H, Proctor RH, Lee T. Identification of a 12-gene fusaric acid biosynthetic gene cluster in Fusarium species through comparative and functional genomics. Mol Plant-Microbe Interact. 2015;28(3):319–32.
Gaffoor I, Trail F. Characterization of two polyketide synthase genes involved in zearalenone biosynthesis in Gibberella zeae. Appl Environ Microbiol. 2006;72(3):1793–9.
Ludwig N, Lohrer M, Hempel M, Mathea S, Schliebner I, Menzel M, Kiesow A, Schaffrath U, Deising HB, Horbach R. Melanin is not required for turgor generation but enhances cell-wall rigidity in appressoria of the corn pathogen Colletotrichum graminicola. Mol Plant-Microbe Interact. 2014;27(4):315–27.
Gaffoor I, Brown DW, Plattner R, Proctor RH, Qi W, Trail F. Functional analysis of the polyketide synthase genes in the filamentous fungus Gibberella zeae (anamorph Fusarium graminearum). Eukaryot Cell. 2005;4(11):1926–33.
Garcia-Rodas R, Trevijano-Contador N, Roman E, Janbon G, Moyrand F, Pla J, Casadevall A, Zaragoza O. Role of Cln1 during melanization of Cryptococcus neoformans. Front Microbiol. 2015;6:798.
Watanabe A, Fujii I, Tsai H, Chang YC, Kwon-Chung KJ, Ebizuka Y. Aspergillus fumigatus alb1 encodes naphthopyrone synthase when expressed in Aspergillus oryzae. FEMS Microbiol Lett. 2000;192(1):39–44.
Langfelder K, Jahn B, Gehringer H, Schmidt A, Wanner G, Brakhage AA. Identification of a polyketide synthase gene (pksP) of Aspergillus fumigatus involved in conidial pigment biosynthesis and virulence. Med Microbiol Immunol. 1998;187(2):79–89.
Feng B, Wang X, Hauser M, Kaufmann S, Jentsch S, Haase G, Becker JM, Szaniszlo PJ. Molecular cloning and characterization of WdPKS1, a gene involved in dihydroxynaphthalene melanin biosynthesis and virulence in Wangiella (Exophiala) dermatitidis. Infect Immun. 2001;69(3):1781–94.
Leng Y, Zhong S. Sfp-type 4′-phosphopantetheinyl transferase is required for lysine synthesis, tolerance to oxidative stress and virulence in the plant pathogenic fungus Cochliobolus sativus. Mol Plant Pathol. 2012;13(4):375–87.
Chooi YH, Muria-Gonzalez MJ, Mead OL, Solomon PS. SnPKS19 encodes the Polyketide Synthase for Alternariol Mycotoxin Biosynthesis in the wheat pathogen Parastagonospora nodorum. Appl Environ Microbiol. 2015;81(16):5309–17.
The Authors take this opportunity to remember prof. Giovanni Del Sorbo who started this study and to thanks dr. Paolo Alfonso Pedata who helped with statistical analysis.
This project was financially supported by the projects: BIP BioIndustrial Processes, CUP: B25C13000290007 and CARINA, CUP: B25B09000080007 and by the PhD program of the University of Palermo.
Availability of data and materials
The reassembled WGS genome of D. helianthi strain 7/96 is present in NCBI GenBank with accession number: MAVT00000000.2.
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
List of the 40 putative PKS in Diaporthe helianthi isolate 7/96 and related information. Domain analyses revealed that most PKS genes coded by D. helianthi are highly reducing PKSs, whereas only eight PKSs lack reducing domains and cluster with non-reducing PKSs. In grey are highlighted the genes partially sequenced. (XLSX 76 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Ruocco, M., Baroncelli, R., Cacciola, S.O. et al. Polyketide synthases of Diaporthe helianthi and involvement of DhPKS1 in virulence on sunflower. BMC Genomics 19, 27 (2018). https://doi.org/10.1186/s12864-017-4405-z
- Polyketide synthases
- Diaporthe helianthi
- Plant pathogen
- Pathogen virulence, toxins