Clonostachys rosea 'omics profiling: identification of putative metabolite-gene associations mediating its in vitro antagonism against Fusarium graminearum.

BACKGROUND: Clonostachys rosea is an established biocontrol agent. Selected strains have either mycoparasitic activity against known pathogens (e.g. Fusarium species) and/or plant growth promoting activity on various crops. Here we report outcomes from a comparative 'omics analysis leveraging a temporal variation in the in vitro antagonistic activities of C. rosea strains ACM941 and 88-710, toward understanding the molecular mechanisms underpinning mycoparasitism. RESULTS: Transcriptomic data highlighted specialized metabolism and membrane transport related genes as being significantly upregulated in ACM941 compared to 88-710 at a time point when the ACM941 strain had higher in vitro antagonistic activity than 88-710. In addition, high molecular weight specialized metabolites were differentially secreted by ACM941, with accumulation patterns of some metabolites matching the growth inhibition differences displayed by the exometabolites of the two strains. In an attempt to identify statistically relevant relationships between upregulated genes and differentially secreted metabolites, transcript and metabolomic abundance data were associated using IntLIM (Integration through Linear Modeling). Of several testable candidate associations, a putative C. rosea epidithiodiketopiperazine (ETP) gene cluster was identified as a prime candidate based on both co-regulation analysis and transcriptomic-metabolomic data association. CONCLUSIONS: Although remaining to be validated functionally, these results suggest that a data integration approach may be useful for identification of potential biomarkers underlying functional divergence in C. rosea strains.

Draft Genome Resources for Plant-Beneficial Fungi Clonostachys rosea Strains ACM941 and 88-710.

Clonostachys rosea strains ACM941 and 88-710 are beneficial microbes recognized for their plant disease control and growth promotion properties, respectively, when applied to economically important crops. In addition to their geographical and functional overlap, the two strains also share a high degree of genetic similarity. In an effort to identify the subtleties that underlie their strain-specific applications, their genomic sequence is reported here. The genome size of ACM941 was estimated to be 56.9 Mb, encoding 17,585 putative genes, while strain 88-710 was estimated to have a 55.5 Mb genome size, containing 17,188 predicted genes. Overall, ACM941 and 88-710 share >96% of their encoded genomes, such that their strain-specific characteristics are likely encoded in either the remaining variable 4% or differentially regulated shared genes or both. These genomic sequences form a foundation for future studies aimed at identifying the genomic and metabolic machinery driving their respective beneficial properties.

Barley "uzu" and Wheat "uzu-like" Brassinosteroid Receptor BRI1 Kinase Domain Variations Modify Phosphorylation Activity In Vitro.

Brassinosteroid insensitive1 (BRI1), a leucine-rich repeat receptor kinase, is responsible for the perception of the brassinosteroid (BR) phytohormone in plants. While recent evidence has implicated a naturally occurring Hordeum vulgare V. (barley) HvBRI1 kinase domain (KD) variant (H857R; "uzu" variation) in increased fungal disease resistance, the impact of the variation on receptor function and thus the mechanism by which disease resistance might be imparted remain enigmatic. Here, the functional implications of the uzu variation as well as the effects of newly identified naturally occurring Triticum aestivum L. (wheat) TaBRI1-KD variants are investigated. Recombinantly produced KDs of wild-type (WT) and uzu HvBRI1 were assessed for phosphorylation activity in vitro, yielding WT KM and VMAX values similar to those of other reports, but the uzu variation delayed saturation and reduced turnover levels. In silico modeling of the H857R variation showed it to be surface-exposed and distal from the catalytic site. Further evaluation of three naturally occurring wheat TaBRI1 variants, A907T, A970V, and G1019R (barley numbering) identified in the A, B, and D subgenomic genes, respectively, highlighted a significant loss of activity for A907T. A907T is located on the same surface as the H857R variation and a negative regulatory phosphorylation site (T982) in Arabidopsis thaliana BRI1. A fourth variation, T1031A (barley numbering), unique to both subgenomic A proteins and localized to the BKI1 binding site, also decreased activity. The outcomes are discussed with respect to the predicted structural contexts of the variations and their implications with respect to mechanisms of action.

Transcriptomic and Exometabolomic Profiling Reveals Antagonistic and Defensive Modes of Clonostachys rosea Action Against Fusarium graminearum.

The mycoparasite Clonostachys rosea ACM941 is under development as a biocontrol organism against Fusarium graminearum, the causative agent of Fusarium head blight in cereals. To identify molecular factors associated with this interaction, the transcriptomic and exometabolomic profiles of C. rosea and F. graminearum GZ3639 were compared during coculture. Prior to physical contact, the antagonistic activity of C. rosea correlated with a response heavily dominated by upregulation of polyketide synthase gene clusters, consistent with the detected accumulation of corresponding secondary metabolite products. Similarly, prior to contact, trichothecene gene clusters were upregulated in F. graminearum, while those responsible for fusarielin and fusarin biosynthesis were downregulated, correlating with an accumulation of trichothecene products in the interaction zone over time. A concomitant increase in 15-acetyl deoxynivalenol-3-glucoside in the interaction zone was also detected, with C. rosea established as the source of this detoxified mycotoxin. After hyphal contact, C. rosea was found to predominantly transcribe genes encoding cell wall-degradation enzymes, major facilitator superfamily sugar transporters, anion:cation symporters, as well as alternative carbon source utilization pathways, together indicative of a transition to necrotropism at this stage. F. graminearum notably activated the transcription of phosphate starvation pathway signature genes at this time. Overall, a number of signature molecular mechanisms likely contributing to antagonistic activity by C. rosea against F. graminearum, as well as its mycotoxin tolerance, are identified in this report, yielding several new testable hypotheses toward understanding the basis of C. rosea as a biocontrol agent for continued agronomic development and application.

A lavender ABC transporter confers resistance to monoterpene toxicity in yeast.

MAIN CONCLUSION: Functional expression of a multidrug resistance-type ABC transporter from Lavandulaangustifolia improved yeast resistance to geraniol, a monoterpene constituent of lavender essential oil. Plant ATP-binding cassette (ABC) transporters are a large family of membrane proteins involved in active and selective transport of structurally diverse compounds. In this study, we functionally evaluated LaABCB1, a multidrug resistance (MDR)-type ABC transporter strongly expressed in the secretory cells of lavender glandular trichomes, where monoterpene essential oil constituents are synthesized and secreted. We used LaABCB1 to complement a yeast knockout mutant in which 16 ABC transporters were deleted. Expression of LaABCB1 enhanced tolerance of yeast mutants to geraniol, a key constituent of essential oils in lavenders and numerous other plants. Our findings suggest a role for the MDR-type ABC transporters in the toxicity tolerance of at least certain essential oil constituents in lavender oil glands.

Identification, validation and cross-species transferability of novel Lavandula EST-SSRs.

MAIN CONCLUSION: We identified and characterized EST-SSRs with strong discrimination power against Lavandula angustifolia and Lavandula x intermedia . The markers also showed considerable cross-species transferability rate into six related Lavandula species. Lavenders (Lavandula) are important economical crops grown around the globe for essential oil production. In an attempt to develop genetic markers for these plants, we analyzed over 13,000 unigenes developed from L. angustifolia and L. x intermedia EST databases, and identified 3,459 simple sequence repeats (SSR), which were dominated by trinucleotides (41.2 %) and dinucleotides (31.45 %). Approximately, 19 % of the unigenes contained at least one SSR marker, over 60 % of which were localized in the UTRs. Only 252 EST-SSRs were 18 bp or longer from which 31 loci were validated, and 24 amplified discrete fragments with 85 % polymorphism in L. x intermedia and L. angustifolia. The average number of alleles in L. x intermedia and L. angustifolia were 3.42 and 3.71 per marker with average PIC values of 0.47 and 0.52, respectively. These values suggest a moderate to strong level of informativeness for the markers, with some loci producing unique fingerprints. The cross-species transferability rate of the markers ranges 50-100 % across eight species. The utility of these markers was assessed in eight Lavandula species and 15 L. angustifolia and L. x intermedia cultivars, and the dendrogram deduced from their similarity indexes successfully delineated the species into their respective sections and the cultivars into their respective species. These markers have potential for application in fingerprinting, diversity studies and marker-assisted breeding of Lavandula.

The biosynthetic origin of irregular monoterpenes in Lavandula: isolation and biochemical characterization of a novel cis-prenyl diphosphate synthase gene, lavandulyl diphosphate synthase.

Lavender essential oils are constituted predominantly of regular monoterpenes, for example linalool, 1,8-cineole, and camphor. However, they also contain irregular monoterpenes including lavandulol and lavandulyl acetate. Although the majority of genes responsible for the production of regular monoterpenes in lavenders are now known, enzymes (including lavandulyl diphosphate synthase (LPPS)) catalyzing the biosynthesis of irregular monoterpenes in these plants have not been described. Here, we report the isolation and functional characterization of a novel cis-prenyl diphosphate synthase cDNA, termed Lavandula x intermedia lavandulyl diphosphate synthase (LiLPPS), through a homology-based cloning strategy. The LiLPPS ORF, encoding for a 305-amino acid long protein, was expressed in Escherichia coli, and the recombinant protein was purified by nickel-nitrilotriacetic acid affinity chromatography. The approximately 34.5-kDa bacterially produced protein specifically catalyzed the head-to-middle condensation of two dimethylallyl diphosphate units to LPP in vitro with apparent Km and kcat values of 208 ± 12 µm and 0.1 s(-1), respectively. LiLPPS is a homodimeric enzyme with a sigmoidal saturation curve and Hill coefficient of 2.7, suggesting a positive co-operative interaction among its catalytic sites. LiLPPS could be used to modulate the production of lavandulol and its derivatives in plants through metabolic engineering.

Cloning of a sesquiterpene synthase from Lavandula x intermedia glandular trichomes.

The essential oil (EO) of Lavandula is dominated by monoterpenes, but can also contain small amounts of sesquiterpenes, depending on species and environmental conditions. For example, the sesquiterpene 9-epi-caryophyllene can make up to 8 % of the EO in a few species, including those commercially propagated for EO production. Here, we report the cloning and functional characterization of 9-epi-caryophyllene synthase (LiCPS) from the glandular trichomes of Lavandula x intermedia, cv. Grosso. The 1,617 bp open reading frame of LiCPS, which did not encode a transit peptide, was expressed in Escherichia coli and the recombinant protein purified by Ni-NTA agarose affinity chromatography. The ca. 60 kDa recombinant protein specifically converted farnesyl diphosphate to 9-epi-caryophyllene. LiCPS also produced a few monoterpenes when assayed with the monoterpene precursor geranyl diphosphate (GPP), but--unlike most monoterpene synthases--was not able to derive detectable amounts of any products from the cis isomer of GPP, neryl diphosphate. The LiCPS transcripts accumulated in developing L. x intermedia flowers and were highly enriched in glandular trichomes, but were not detected in leaves suggesting that the transcriptional expression of this gene is spatially and developmentally regulated.

Multiple Clonostachys rosea UDP-Glycosyltransferases Contribute to the Production of 15-Acetyl-Deoxynivalenol-3-O-Glycoside When Confronted with Fusarium graminearum.

Mycotoxins, derived from toxigenic fungi such as Fusarium, Aspergillus, and Penicillium species have impacted the human food chain for thousands of years. Deoxynivalenol (DON), is a tetracyclic sesquiterpenoid type B trichothecene mycotoxin predominantly produced by F. culmorum and F. graminearum during the infection of corn, wheat, oats, barley, and rice. Glycosylation of DON is a protective detoxification mechanism employed by plants. More recently, DON glycosylating activity has also been detected in fungal microparasitic (biocontrol) fungal organisms. Here we follow up on the reported conversion of 15-acetyl-DON (15-ADON) into 15-ADON-3-O-glycoside (15-ADON-3G) in Clonostachys rosea. Based on the hypothesis that the reaction is likely being carried out by a uridine diphosphate glycosyl transferase (UDP-GTase), we applied a protein structural comparison strategy, leveraging the availability of the crystal structure of rice Os70 to identify a subset of potential C. rosea UDP-GTases that might have activity against 15-ADON. Using CRISPR/Cas9 technology, we knocked out several of the selected UDP-GTases in the C. rosea strain ACM941. Evaluation of the impact of knockouts on the production of 15-ADON-3G in confrontation assays with F. graminearum revealed multiple UDP-GTase enzymes, each contributing partial activities. The relationship between these positive hits and other UDP-GTases in fungal and plant species is discussed.

Cloning, functional characterization and genomic organization of 1,8-cineole synthases from Lavandula.

Several members of the genus Lavandula produce valuable essential oils (EOs) that are primarily constituted of the low molecular weight isoprenoids, particularly monoterpenes. We isolated over 8,000 ESTs from the glandular trichomes of L. x intermedia flowers (where bulk of the EO is synthesized) to facilitate the discovery of genes that control the biosynthesis of EO constituents. The expression profile of these ESTs in L. x intermedia and its parents L. angustifolia and L. latifolia was established using microarrays. The resulting data highlighted a differentially expressed, previously uncharacterized cDNA with strong homology to known 1,8-cineole synthase (CINS) genes. The ORF, excluding the transit peptide, of this cDNA was expressed in E. coli, purified by Ni-NTA agarose affinity chromatography and functionally characterized in vitro. The ca. 63 kDa bacterially produced recombinant protein, designated L. x intermedia CINS (LiCINS), converted geranyl diphosphate (the linear monoterpene precursor) primarily to 1,8-cineole with K ( m ) and k ( cat ) values of 5.75 µM and 8.8 × 10(-3) s(-1), respectively. The genomic DNA of CINS in the studied Lavandula species had identical exon-intron architecture and coding sequences, except for a single polymorphic nucleotide in the L. angustifolia ortholog which did not alter protein function. Additional nucleotide variations restricted to L. angustifolia introns were also observed, suggesting that LiCINS was most likely inherited from L. latifolia. The LiCINS mRNA levels paralleled the 1,8-cineole content in mature flowers of the three lavender species, and in developmental stages of L. x intermedia inflorescence indicating that the production of 1,8 cineole in Lavandula is most likely controlled through transcriptional regulation of LiCINS.

Molecular cloning and functional characterization of borneol dehydrogenase from the glandular trichomes of Lavandula x intermedia.

Several varieties of Lavandula x intermedia (lavandins) are cultivated for their essential oils (EOs) for use in cosmetic, hygiene and personal care products. These EOs are mainly constituted of monoterpenes including camphor, which contributes an off odor reducing the olfactory appeal of the oil. We have recently constructed a cDNA library from the glandular trichomes (the sites of EO synthesis) of L. x intermedia plants. Here, we describe the cloning of a borneol dehydrogenase cDNA (LiBDH) from this library. The 780 bp open reading frame of the cDNA encoded a 259 amino acid short chain alcohol dehydrogenase with a predicted molecular mass of ca. 27.5 kDa. The recombinant LiBDH was expressed in Escherichia coli, purified by Ni-NTA agarose affinity chromatography, and functionally characterized in vitro. The bacterially produced enzyme specifically converted borneol to camphor as the only product with K(m) and k(cat) values of 53 µM and 4.0 × 10(-4) s(-1), respectively. The LiBDH transcripts were specifically expressed in glandular trichomes of mature flowers indicating that like other Lavandula monoterpene synthases the expression of this gene is regulated in a tissue-specific manner. The cloning of LiBDH has far reaching implications in improving the quality of Lavandula EOs through metabolic engineering.

Cloning and functional characterization of ß-phellandrene synthase from Lavandula angustifolia.

En route to building genomics resources for Lavandula, we have obtained over 14,000 ESTs for leaves and flowers of L. angustifolia, a major essential oil crop, and identified a number of previously uncharacterized terpene synthase (TPS) genes. Here we report the cloning, expression in E. coli, and functional characterization of ß-phellandrene synthase, LaßPHLS. The ORF--excluding the transit peptide--for this gene encoded a 62.3 kDa protein that contained all conserved motifs present in plant TPSs. Expression in bacteria resulted in the production of a soluble protein that was purified by Ni-NTA agarose affinity chromatography. While the recombinant LaßPHLS did not utilize FPP as a substrate, it converted GPP (the preferred substrate) and NPP into ß-phellandrene as the major product, with K (m) and k (cat) of 6.55 µM and 1.75 × 10(-2) s(-1), respectively, for GPP. The LaßPHLS transcripts were highly abundant in young leaves where ß-phellandrene is produced, but were barely detectable in flowers and older leaves, where ß-phellandrene is not synthesized in significant quantities. This data indicate that ß-phellandrene biosynthesis is transcriptionally and developmentally regulated. We also cloned and expressed in E. coli a second TPS-like protein, LaTPS-I, that lacks an internal stretch of 73 amino acids, including the signature DDxxD divalent metal binding motif, compared to other plant TPSs. The recombinant LaTPS-I did not produce detectable products in vitro when assayed with GPP, NPP or FPP as substrates. The lack of activity is most likely due to the absence of catalytically important amino acid residues within the missing region.

A Universally Primed-Polymerase Chain Reaction (UP-PCR) Marker to Discriminate Clonostachys rosea ACM941 from Related Strains.

Clonostachys rosea strain ACM941 is an effective biocontrol agent against several crop diseases including Fusarium head blight. In anticipation of its increased relevance going forward, the development of a reliable DNA-based molecular marker to track it is essential. Universally primed-PCR (UP-PCR) has been used successfully to differentiate other C. rosea strains. Herein, the development of a UP-PCR marker for ACM941 is described. A combination of two primers (AS15 and L45) produced a ~450 bp fragment that was unique to ACM941 compared to other commercial biocontrol agents. Primers subsequently designed based on the obtained fragment also produced a similarly unique band from ACM941 alone. BLAST analysis of the amplified sequence did not yield any homologous sequence in available online databases or within the closely related C. rosea IK726 and CBS125111 strains' genomes. The specificity of this marker for ACM941 was validated against ten additional C. rosea strains isolated from Canada, with ACM941 producing the brightest band. Taken together, these results imply that the UP-PCR primers AS15 and L45 and the amplified fragment can be used to detect and monitor the ACM941 strain after its release into the environment.

Profiling of the Transcriptomic Responses of Clonostachys rosea Upon Treatment With Fusarium graminearum Secretome.

Clonostachys rosea strain ACM941 is a fungal bio-control agent patented against the causative agent of Fusarium Head Blight, Fusarium graminearum. Although the molecular details remain enigmatic, previous studies have suggested that C. rosea may secrete F. graminearum growth inhibitors. Further toward this, experiments described herein show that induction of C. rosea cultures by the addition of an aliquot of F. graminearum(Fg)-spent media (including macroconidia), yield C. rosea (Cr)-spent media that elicited higher anti-F. graminearum activity than either control or deoxynivalenol (DON)-induced Cr-spent media. To gain additional insight into the genetic and metabolic factors modulating this interaction, transcriptomic (RNAseq) profiles of C. rosea in response to DON and Fg-spent media treatment, were developed. This analysis revealed 24,112 C. rosea unigenes, of which 5,605 and 6,285 were differentially regulated by DON and F-spent media, respectively. More than half of these unigenes were up-regulated, with annotations, most notably in the Fg-spent media treatment data, suggesting enhancement of polyketide (PK) and non-ribosomal peptide (NRP) secondary metabolite precursor synthesis, and PK/NRP-like synthases. Four ABC transporters were also up-regulated in response to Fg-spent media. Further analysis showed that the PK and NRP-like synthases belong to three gene clusters that also include ABC transporters, and other genes known to tailor secondary metabolite biosynthesis. The RNAseq data was further validated using quantitative RT-qPCR. Taken together, these results show that C. rosea responds to the presence of Fg-spent media (and to a lesser extent, DON-alone) by up-regulating unique aspects of its secondary metabolism-related genetic repertoire. The identities and roles of C. rosea secondary metabolites produced by the targeted gene clusters are now under investigation.