Novel Secreted Effectors Conserved Among Smut Fungi Contribute to the Virulence of Ustilago maydis

Fungal pathogens deploy a set of molecules (proteins, specialized metabolites, and sRNAs), so-called effectors, to aid the infection process. In comparison to other plant pathogens, smut fungi have small genomes and secretomes of 20 Mb and around 500 proteins, respectively. Previous comparative genomic studies have shown that many secreted effector proteins without known domains, i.e., novel, are conserved only in the Ustilaginaceae family. By analyzing the secretomes of 11 species within Ustilaginaceae, we identified 53 core homologous groups commonly present in this lineage. By collecting existing mutants and generating additional ones, we gathered 44 Ustilago maydis strains lacking single core effectors as well as 9 strains containing multiple deletions of core effector gene families. Pathogenicity assays revealed that 20 of these 53 mutant strains were affected in virulence. Among the 33 mutants that had no obvious phenotypic changes, 13 carried additional, sequence-divergent, structurally similar paralogs. We report a virulence contribution of seven previously uncharacterized single core effectors and of one effector family. Our results help to prioritize effectors for understanding U. maydis virulence and provide genetic resources for further characterization. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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

BACKGROUND

Fungi are important sources for bioactive compounds that find their applications in many important sectors like in the pharma-, food- or agricultural industries. In an environmental monitoring project for fungi involved in soil nitrogen cycling we also isolated Cephalotrichum gorgonifer (strain NG_p51). In the course of strain characterisation work we found that this strain is able to naturally produce high amounts of rasfonin, a polyketide inducing autophagy, apoptosis, necroptosis in human cell lines and showing anti-tumor activity in KRAS-dependent cancer cells.

RESULTS

In order to elucidate the biosynthetic pathway of rasfonin, the strain was genome sequenced, annotated, submitted to transcriptome analysis and genetic transformation was established. Biosynthetic gene cluster (BGC) prediction revealed the existence of 22 BGCs of which the majority was not expressed under our experimental conditions. In silico prediction revealed two BGCs with a suite of enzymes possibly involved in rasfonin biosynthesis. Experimental verification by gene-knock out of the key enzyme genes showed that one of the predicted BGCs is indeed responsible for rasfonin biosynthesis.

CONCLUSIONS

This study identified a biosynthetic gene cluster containing a key-gene responsible for rasfonin production. Additionally, molecular tools were established for the non-model fungus Cephalotrichum gorgonifer which allows strain engineering and heterologous expression of the BGC for high rasfonin producing strains and the biosynthesis of rasfonin derivates for diverse applications.

Machine Learning Identifies New Predictors on Restenosis Risk after Coronary Artery Stenting in 10,004 Patients with Surveillance Angiography

OBJECTIVE

Machine learning (ML) approaches have the potential to uncover regular patterns in multi-layered data. Here we applied self-organizing maps (SOMs) to detect such patterns with the aim to better predict in-stent restenosis (ISR) at surveillance angiography 6 to 8 months after percutaneous coronary intervention with stenting.

METHODS

In prospectively collected data from 10,004 patients receiving percutaneous coronary intervention (PCI) for 15,004 lesions, we applied SOMs to predict ISR angiographically 6-8 months after index procedure. SOM findings were compared with results of conventional uni- and multivariate analyses. The predictive value of both approaches was assessed after random splitting of patients into training and test sets (50:50).

RESULTS

Conventional multivariate analyses revealed 10, mostly known, predictors for restenosis after coronary stenting: balloon-to-vessel ratio, complex lesion morphology, diabetes mellitus, left main stenting, stent type (bare metal vs. first vs. second generation drug eluting stent), stent length, stenosis severity, vessel size reduction, and prior bypass surgery. The SOM approach identified all these and nine further predictors, including chronic vessel occlusion, lesion length, and prior PCI. Moreover, the SOM-based model performed well in predicting ISR (AUC under ROC: 0.728); however, there was no meaningful advantage in predicting ISR at surveillance angiography in comparison with the conventional multivariable model (0.726, p = 0.3).

CONCLUSIONS

The agnostic SOM-based approach identified-without clinical knowledge-even more contributors to restenosis risk. In fact, SOMs applied to a large prospectively sampled cohort identified several novel predictors of restenosis after PCI. However, as compared with established covariates, ML technologies did not improve identification of patients at high risk for restenosis after PCI in a clinically relevant fashion.

Use of big data from health insurance for assessment of cardiovascular outcomes

Outcome research that supports guideline recommendations for primary and secondary preventions largely depends on the data obtained from clinical trials or selected hospital populations. The exponentially growing amount of real-world medical data could enable fundamental improvements in cardiovascular disease (CVD) prediction, prevention, and care. In this review we summarize how data from health insurance claims (HIC) may improve our understanding of current health provision and identify challenges of patient care by implementing the perspective of patients (providing data and contributing to society), physicians (identifying at-risk patients, optimizing diagnosis and therapy), health insurers (preventive education and economic aspects), and policy makers (data-driven legislation). HIC data has the potential to inform relevant aspects of the healthcare systems. Although HIC data inherit limitations, large sample sizes and long-term follow-up provides enormous predictive power. Herein, we highlight the benefits and limitations of HIC data and provide examples from the cardiovascular field, i.e. how HIC data is supporting healthcare, focusing on the demographical and epidemiological differences, pharmacotherapy, healthcare utilization, cost-effectiveness and outcomes of different treatments. As an outlook we discuss the potential of using HIC-based big data and modern artificial intelligence (AI) algorithms to guide patient education and care, which could lead to the development of a learning healthcare system and support a medically relevant legislation in the future.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

Transcriptome-wide association study of coronary artery disease identifies novel susceptibility genes

The majority of risk loci identified by genome-wide association studies (GWAS) are in non-coding regions, hampering their functional interpretation. Instead, transcriptome-wide association studies (TWAS) identify gene-trait associations, which can be used to prioritize candidate genes in disease-relevant tissue(s). Here, we aimed to systematically identify susceptibility genes for coronary artery disease (CAD) by TWAS. We trained prediction models of nine CAD-relevant tissues using EpiXcan based on two genetics-of-gene-expression panels, the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and the Genotype-Tissue Expression (GTEx). Based on these prediction models, we imputed gene expression of respective tissues from individual-level genotype data on 37,997 CAD cases and 42,854 controls for the subsequent gene-trait association analysis. Transcriptome-wide significant association (i.e. P < 3.85e-6) was observed for 114 genes. Of these, 96 resided within previously identified GWAS risk loci and 18 were novel. Stepwise analyses were performed to study their plausibility, biological function, and pathogenicity in CAD, including analyses for colocalization, damaging mutations, pathway enrichment, phenome-wide associations with human data and expression-traits correlations using mouse data. Finally, CRISPR/Cas9-based gene knockdown of two newly identified TWAS genes, RGS19 and KPTN, in a human hepatocyte cell line resulted in reduced secretion of APOB100 and lipids in the cell culture medium. Our CAD TWAS work (i) prioritized candidate causal genes at known GWAS loci, (ii) identified 18 novel genes to be associated with CAD, and iii) suggested potential tissues and pathways of action for these TWAS CAD genes.

Polygenic risk for coronary artery disease in the Scottish and English population

Background

Epidemiological studies have repeatedly observed a markedly higher risk for coronary artery disease (CAD) in Scotland as compared to England. Up to now, it is unclear whether environmental or genetic factors might explain this phenomenon.

Methods

Using UK Biobank (UKB) data, we assessed CAD risk, based on the Framingham risk score (FRS) and common genetic variants, to explore the respective contribution to CAD prevalence in Scotland (n = 31,963) and England (n = 317,889). We calculated FRS based on sex, age, body mass index (BMI), total cholesterol (TC), high density lipoprotein cholesterol (HDL-C), systolic blood pressure (SBP), antihypertensive medication, smoking status, and diabetes. We determined the allele frequency of published genome-wide significant risk CAD alleles and a weighted genetic risk score (wGRS) for quantifying genetic CAD risk.

Results

Prevalence of CAD was 16% higher in Scotland as compared to England (8.98% vs. 7.68%, P < 0.001). However, the FRS only predicted a marginally higher CAD risk (less than 1%) in Scotland (12.5 ± 10.5 vs.12.6 ± 10.6, P = 0.03). Likewise, the overall number of genome-wide significant variants affecting CAD risk (157.6 ± 7.7 and 157.5 ± 7.7; P = 0.12) and a wGRS for CAD (2.49 ± 0.25 in both populations, P = 0.14) were remarkably similar in the English and Scottish population. Interestingly, we observed substantial differences in the allele frequencies of individual risk variants. Of the previously described 163 genome-wide significant variants studied here, 35 variants had higher frequencies in Scotland, whereas 37 had higher frequencies in England (P < 0.001 each).

Conclusions

Neither the traditional risk factors included in the FRS nor a genetic risk score (GRS) based on established common risk alleles explained the higher CAD prevalence in Scotland. However, we observed marked differences in the distribution of individual risk alleles, which emphasizes that even geographically and ethnically closely related populations may display relevant differences in the genetic architecture of a common disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12872-021-02398-4.

Vascular Tissue Specific miRNA Profiles Reveal Novel Correlations with Risk Factors in Coronary Artery Disease

Cardiovascular disease (CVD) is the leading cause of morbidity and mortality worldwide. Non-coding RNAs have already been linked to CVD development and progression. While microRNAs (miRs) have been well studied in blood samples, there is little data on tissue-specific miRs in cardiovascular relevant tissues and their relation to cardiovascular risk factors. Tissue-specific miRs derived from Arteria mammaria interna (IMA) from 192 coronary artery disease (CAD) patients undergoing coronary artery bypass grafting (CABG) were analyzed. The aims of the study were 1) to establish a reference atlas which can be utilized for identification of novel diagnostic biomarkers and potential therapeutic targets, and 2) to relate these miRs to cardiovascular risk factors. Overall, 393 individual miRs showed sufficient expression levels and passed quality control for further analysis. We identified 17 miRs-miR-10b-3p, miR-10-5p, miR-17-3p, miR-21-5p, miR-151a-5p, miR-181a-5p, miR-185-5p, miR-194-5p, miR-199a-3p, miR-199b-3p, miR-212-3p, miR-363-3p, miR-548d-5p, miR-744-5p, miR-3117-3p, miR-5683 and miR-5701-significantly correlated with cardiovascular risk factors (correlation coefficient >0.2 in both directions, p-value (p < 0.006, false discovery rate (FDR) <0.05). Of particular interest, miR-5701 was positively correlated with hypertension, hypercholesterolemia, and diabetes. In addition, we found that miR-629-5p and miR-98-5p were significantly correlated with acute myocardial infarction. We provide a first atlas of miR profiles in IMA samples from CAD patients. In perspective, these miRs might play an important role in improved risk assessment, mechanistic disease understanding and local therapy of CAD.

Genome sequencing, annotation and exploration of the SO2-tolerant non-conventional yeast Saccharomycodes ludwigii

BACKGROUND

Saccharomycodes ludwigii belongs to the poorly characterized Saccharomycodeacea family and is known by its ability to spoil wines, a trait mostly attributable to its high tolerance to sulfur dioxide (SO₂). To improve knowledge about Saccharomycodeacea our group determined whole-genome sequences of Hanseniaspora guilliermondii (UTAD222) and S. ludwigii (UTAD17), two members of this family. While in the case of H. guilliermondii the genomic information elucidated crucial aspects concerning the physiology of this species in the context of wine fermentation, the draft sequence obtained for S. ludwigii was distributed by more than 1000 contigs complicating extraction of biologically relevant information. In this work we describe the results obtained upon resequencing of S. ludwigii UTAD17 genome using PacBio as well as the insights gathered from the exploration of the annotation performed over the assembled genome.

RESULTS

Resequencing of S. ludwigii UTAD17 genome with PacBio resulted in 20 contigs totaling 13 Mb of assembled DNA and corresponding to 95% of the DNA harbored by this strain. Annotation of the assembled UTAD17 genome predicts 4644 protein-encoding genes. Comparative analysis of the predicted S. ludwigii ORFeome with those encoded by other Saccharomycodeacea led to the identification of 213 proteins only found in this species. Among these were six enzymes required for catabolism of N-acetylglucosamine, four cell wall β-mannosyltransferases, several flocculins and three acetoin reductases. Different from its sister Hanseniaspora species, neoglucogenesis, glyoxylate cycle and thiamine biosynthetic pathways are functional in S. ludwigii. Four efflux pumps similar to the Ssu1 sulfite exporter, as well as robust orthologues for 65% of the S. cerevisiae SO₂-tolerance genes, were identified in S. ludwigii genome.

CONCLUSIONS

This work provides the first genome-wide picture of a S. ludwigii strain representing a step forward for a better understanding of the physiology and genetics of this species and of the Saccharomycodeacea family. The release of this genomic sequence and of the information extracted from it can contribute to guide the design of better wine preservation strategies to counteract spoilage prompted by S. ludwigii. It will also accelerate the exploration of this species as a cell factory, specially in production of fermented beverages where the use of Non-Saccharomyces species (including spoilage species) is booming.

Functional investigation of the coronary artery disease gene SVEP1

A missense variant of the sushi, von Willebrand factor type A, EGF and pentraxin domain containing protein 1 (SVEP1) is genome-wide significantly associated with coronary artery disease. The mechanisms how SVEP1 impacts atherosclerosis are not known. We found endothelial cells (EC) and vascular smooth muscle cells to represent the major cellular source of SVEP1 in plaques. Plaques were larger in atherosclerosis-prone Svep1 haploinsufficient (ApoE^−/−Svep1^+/−) compared to Svep1 wild-type mice (ApoE^−/−Svep1^+/+) and ApoE^−/−Svep1^+/− mice displayed elevated plaque neutrophil, Ly6C^high monocyte, and macrophage numbers. We assessed how leukocytes accumulated more inside plaques in ApoE^−/−Svep1^+/− mice and found enhanced leukocyte recruitment from blood into plaques. In vitro, we examined how SVEP1 deficiency promotes leukocyte recruitment and found elevated expression of the leukocyte attractant chemokine (C-X-C motif) ligand 1 (CXCL1) in EC after incubation with missense compared to wild-type SVEP1. Increasing wild-type SVEP1 levels silenced endothelial CXCL1 release. In line, plasma Cxcl1 levels were elevated in ApoE^−/−Svep1^+/− mice. Our studies reveal an atheroprotective role of SVEP1. Deficiency of wild-type Svep1 increased endothelial CXCL1 expression leading to enhanced recruitment of proinflammatory leukocytes from blood to plaque. Consequently, elevated vascular inflammation resulted in enhanced plaque progression in Svep1 deficiency.

Genetically determined intelligence and coronary artery disease risk

BACKGROUND

Epidemiological studies have shown inverse association between intelligence and coronary artery disease (CAD) risk, but the underlying mechanisms remain unclear.

METHODS

Based on 242 SNPs independently associated with intelligence, we calculated the genetic intelligence score (gIQ) for participants from 10 CAD case-control studies (n = 34,083) and UK Biobank (n = 427,306). From UK Biobank, we extracted phenotypes including body mass index (BMI), type 2 diabetes (T2D), smoking, hypertension, HDL cholesterol, LDL cholesterol, measured intelligence score, and education attainment. To estimate the effects of gIQ on CAD and its related risk factors, regression analyses was applied. Next, we studied the mediatory roles of measured intelligence and educational attainment. Lastly, Mendelian randomization was performed to validate the findings.

RESULTS

In CAD case-control studies, one standard deviation (SD) increase of gIQ was related to a 5% decrease of CAD risk (odds ratio [OR] of 0.95; 95% confidence interval [CI] 0.93 to 0.98; P = 4.93e-5), which was validated in UK Biobank (OR = 0.97; 95% CI 0.96 to 0.99; P = 6.4e-4). In UK Biobank, we also found significant inverse correlations between gIQ and risk factors of CAD including smoking, BMI, T2D, hypertension, and a positive correlation with HDL cholesterol. The association signals between gIQ and CAD as well as its risk factors got largely attenuated after the adjustment of measured intelligence and educational attainment. The causal role of intelligence in mediating CAD risk was confirmed by Mendelian randomization analyses.

CONCLUSION

Genetic components of intelligence affect measured intelligence and educational attainment, which subsequently affect the prevalence of CAD via a series of unfavorable risk factor profiles.

Infection cushions of Fusarium graminearum are fungal arsenals for wheat infection

Fusarium graminearum is one of the most destructive plant pathogens worldwide, causing fusarium head blight (FHB) on cereals. F. graminearum colonizes wheat plant surfaces with specialized unbranched hyphae called runner hyphae (RH), which develop multicelled complex appressoria called infection cushions (IC). IC generate multiple penetration sites, allowing the fungus to enter the plant cuticle. Complex infection structures are typical for several economically important plant pathogens, yet with unknown molecular basis. In this study, RH and IC formed on the surface of wheat paleae were isolated by laser capture microdissection. RNA-Seq-based transcriptomic analyses were performed on RH and IC and compared to mycelium grown in complete medium (MY). Both RH and IC displayed a high number of infection up-regulated genes (982), encoding, among others, carbohydrate-active enzymes (CAZymes: 140), putative effectors (PE: 88), or secondary metabolism gene clusters (SMC: 12 of 67 clusters). RH specifically up-regulated one SMC corresponding to aurofusarin biosynthesis, a broad activity antibiotic. IC specifically up-regulated 248 genes encoding mostly putative virulence factors such as 7 SMC, including the mycotoxin deoxynivalenol and the newly identified fusaoctaxin A, 33 PE, and 42 CAZymes. Furthermore, we studied selected candidate virulence factors using cellular biology and reverse genetics. Hence, our results demonstrate that IC accumulate an arsenal of proven and putative virulence factors to facilitate the invasion of epidermal cells.

Two Is Better Than One: Studying Ustilago bromivora-Brachypodium Compatibility by Using a Hybrid Pathogen

Pathogenic fungi can have devastating effects on agriculture and health. One potential challenge in dealing with pathogens is the possibility of a host jump (i.e., when a pathogen infects a new host species). This can lead to the emergence of new diseases or complicate the management of existing threats. We studied host specificity by using a hybrid fungus formed by mating two closely related fungi: Ustilago bromivora, which normally infects Brachypodium spp., and U. hordei, which normally infects barley. Although U. hordei was unable to infect Brachypodium spp., the hybrid could. These hybrids also displayed the same mating-type bias that had been observed in U. bromivora and provide evidence of a dominant spore-killer-like system on the sex chromosome of U. bromivora. By analyzing the genomic composition of 109 hybrid strains, backcrossed with U. hordei over four generations, we identified three regions associated with infection on Brachypodium spp. and 75 potential virulence candidates. The most strongly associated region was located on chromosome 8, where seven genes encoding predicted secreted proteins were identified. The fact that we identified several regions relevant for pathogenicity on Brachypodium spp. but that none were essential suggests that host specificity, in the case of U. bromivora, is a multifactorial trait which can be achieved through different subsets of virulence factors.

Biochemical Characterization of the Fusarium graminearum Candidate ACC-Deaminases and Virulence Testing of Knockout Mutant Strains

Fusarium graminearum is a plant pathogenic fungus which is able to infect wheat and other economically important cereal crop species. The role of ethylene in the interaction with host plants is unclear and controversial. We have analyzed the inventory of genes with a putative function in ethylene production or degradation of the ethylene precursor 1-aminocyclopropane carboxylic acid (ACC). F. graminearum, in contrast to other species, does not contain a candidate gene encoding ethylene-forming enzyme. Three genes with similarity to ACC synthases exist; heterologous expression of these did not reveal enzymatic activity. The F. graminearum genome contains in addition two ACC deaminase candidate genes. We have expressed both genes in E. coli and characterized the enzymatic properties of the affinity-purified products. One of the proteins had indeed ACC deaminase activity, with kinetic properties similar to ethylene-stress reducing enzymes of plant growth promoting bacteria. The other candidate was inactive with ACC but turned out to be a d-cysteine desulfhydrase. Since it had been reported that ethylene insensitivity in transgenic wheat increased Fusarium resistance and reduced the content of the mycotoxin deoxynivalenol (DON) in infected wheat, we generated single and double knockout mutants of both genes in the F. graminearum strain PH-1. No statistically significant effect of the gene disruptions on fungal spread or mycotoxin content was detected, indicating that the ability of the fungus to manipulate the production of the gaseous plant hormones ethylene and H₂S is dispensable for full virulence.

Genome sequence of the non-conventional wine yeast Hanseniaspora guilliermondii UTAD222 unveils relevant traits of this species and of the Hanseniaspora genus in the context of wine fermentation

Hanseanispora species, including H. guilliermondii, are long known to be abundant in wine grape-musts and to play a critical role in vinification by modulating, among other aspects, the wine sensory profile. Despite this, the genetics and physiology of Hanseniaspora species remains poorly understood. The first genomic sequence of a H. guilliermondii strain (UTAD222) and the discussion of its potential significance are presented in this work. Metabolic reconstruction revealed that H. guilliermondii is not equipped with a functional gluconeogenesis or glyoxylate cycle, nor does it harbours key enzymes for glycerol or galactose catabolism or for biosynthesis of biotin and thiamine. Also, no fructose-specific transporter could also be predicted from the analysis of H. guilliermondii genome leaving open the mechanisms underlying the fructophilic character of this yeast. Comparative analysis involving H. guilliermondii, H. uvarum, H. opuntiae and S. cerevisiae revealed 14 H. guilliermondii-specific genes (including five viral proteins and one β-glucosidase). Furthermore, 870 proteins were only found within the Hanseniaspora proteomes including several β-glucosidases and decarboxylases required for catabolism of biogenic amines. The release of H. guilliermondii genomic sequence and the comparative genomics/proteomics analyses performed, is expected to accelerate research focused on Hanseniaspora species and to broaden their application in the wine industry and in other bio-industries in which they could be explored as cell factories.

Comparative genomic and transcriptomic analyses unveil novel features of azole resistance and adaptation to the human host in Candida glabrata

The frequent emergence of azole resistance among Candida glabrata strains contributes to increase the incidence of infections caused by this species. Whole-genome sequencing of a fluconazole and voriconazole-resistant clinical isolate (FFUL887) and subsequent comparison with the genome of the susceptible strain CBS138 revealed prominent differences in several genes documented to promote azole resistance in C. glabrata. Among these was the transcriptional regulator CgPdr1. The CgPdr1 FFUL887 allele included a K274Q modification not documented in other azole-resistant strains. Transcriptomic profiling evidenced the upregulation of 92 documented targets of CgPdr1 in the FFUL887 strain, supporting the idea that the K274Q substitution originates a CgPdr1 gain-of-function mutant. The expression of CgPDR1K274Q in the FFUL887 background sensitised the cells against high concentrations of organic acids at a low pH (4.5), but had no detectable effect in tolerance towards other environmental stressors. Comparison of the genome of FFUL887 and CBS138 also revealed prominent differences in the sequence of adhesin-encoding genes, while comparison of the transcriptome of the two strains showed a significant remodelling of the expression of genes involved in metabolism of carbohydrates, nitrogen and sulphur in the FFUL887 strain; these responses likely reflecting adaptive responses evolved by the clinical strain during colonisation of the host.

Production and Role of Hormones During Interaction of Fusarium Species With Maize (Zea mays L.) Seedlings

It has long been known that hormones affect the interaction of a phytopathogen with its host plant. The pathogen can cause changes in plant hormone homeostasis directly by affecting biosynthesis or metabolism in the plant or by synthesizing and secreting the hormone itself. We previously demonstrated that pathogenic fungi of the Fusarium species complex are able to produce three major types of hormones: auxins, cytokinins, and gibberellins. In this work, we explore changes in the levels of these hormones in maize and mango plant tissues infected with Fusarium. The ability to produce individual phytohormones varies significantly across Fusarium species and such differences likely impact host specificity inducing the unique responses noted in planta during infection. For example, the production of gibberellins by F. fujikuroi leads to elongated rice stalks and the suppression of gibberellin biosynthesis in plant tissue. Although all Fusarium species are able to synthesize auxin, sometimes by multiple pathways, the ratio of its free form and conjugates in infected tissue is affected more than the total amount produced. The recently characterized unique pathway for cytokinin de novo synthesis in Fusarium appears silenced or non-functional in all studied species during plant infection. Despite this, a large increase in cytokinin levels was detected in F. mangiferae infected plants, caused likely by the up-regulation of plant genes responsible for their biosynthesis. Thus, the accumulation of active cytokinins may contribute to mango malformation of the reproductive organs upon infection of mango trees. Together, our findings provide insight into the complex role fungal and plant derived hormones play in the fungal-plant interactions.

A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

Ramularia leaf spot has recently emerged as a major threat to barley production world-wide, causing 25% yield loss in many barley growing regions. Here, we provide a new reference genome of the causal agent, the Dothideomycete Ramularia collo-cygni. The assembly of 32 Mb consists of 78 scaffolds. We used RNA-seq to identify 11,622 genes of which 1,303 and 282 are coding for predicted secreted proteins and putative effectors respectively.

The pathogen separated from its nearest sequenced relative, Zymoseptoria tritici ∼27 Ma. We calculated the divergence of the two species on protein level and see remarkably high synonymous and nonsynonymous divergence. Unlike in many other plant pathogens, the comparisons of transposable elements and gene distributions, show a very homogeneous genome for R. collo-cygni. We see no evidence for higher selective pressure on putative effectors or other secreted proteins and repetitive sequences are spread evenly across the scaffolds. These findings could be associated to the predominantly endophytic life-style of the pathogen. We hypothesize that R. collo-cygni only recently became pathogenic and that therefore its genome does not yet show the typical pathogen characteristics. Because of its high scaffold length and improved CDS annotations, our new reference sequence provides a valuable resource for the community for future comparative genomics and population genetics studies.

Set1 and Kdm5 are antagonists for H3K4 methylation and regulators of the major conidiation-specific transcription factor gene ABA1 in Fusarium fujikuroi

Here we present the identification and characterization of the H3K4‐specific histone methyltransferase Set1 and its counterpart, the Jumonji C demethylase Kdm5, in the rice pathogen Fusarium fujikuroi. While Set1 is responsible for all detectable H3K4me2/me3 in this fungus, Kdm5 antagonizes the H3K4me3 mark. Notably, deletion of both SET1 and KDM5 mainly resulted in the upregulation of genome‐wide transcription, also affecting a large set of secondary metabolite (SM) key genes. Although H3K4 methylation is a hallmark of actively transcribed euchromatin, several SM gene clusters located in subtelomeric regions were affected by Set1 and Kdm5. While the regulation of many of them is likely indirect, H3K4me2 levels at gibberellic acid (GA) genes correlated with GA biosynthesis in the wild type, Δkdm5 and OE::KDM5 under inducing conditions. Whereas Δset1 showed an abolished GA₃ production in axenic culture, phytohormone biosynthesis was induced in planta, so that residual amounts of GA₃ were detected during rice infection. Accordingly, Δset1 exhibited a strongly attenuated, though not abolished, virulence on rice. Apart from regulating secondary metabolism, Set1 and Kdm5 function as activator and repressor of conidiation respectively. They antagonistically regulate H3K4me3 levels and expression of the major conidiation‐specific transcription factor gene ABA1 in F. fujikuroi.

The 'PhenoBox', a flexible, automated, open-source plant phenotyping solution

There is a need for flexible and affordable plant phenotyping solutions for basic research and plant breeding. We demonstrate our open source plant imaging and processing solution ('PhenoBox'/'PhenoPipe') and provide construction plans, source code and documentation to rebuild the system. Use of the PhenoBox is exemplified by studying infection of the model grass Brachypodium distachyon by the head smut fungus Ustilago bromivora, comparing phenotypic responses of maize to infection with a solopathogenic Ustilago maydis (corn smut) strain and effector deletion strains, and studying salt stress response in Nicotiana benthamiana. In U. bromivora-infected grass, phenotypic differences between infected and uninfected plants were detectable weeks before qualitative head smut symptoms. Based on this, we could predict the infection outcome for individual plants with high accuracy. Using a PhenoPipe module for calculation of multi-dimensional distances from phenotyping data, we observe a time after infection-dependent impact of U. maydis effector deletion strains on phenotypic response in maize. The PhenoBox/PhenoPipe system is able to detect established salt stress responses in N. benthamiana. We have developed an affordable, automated, open source imaging and data processing solution that can be adapted to various phenotyping applications in plant biology and beyond.

Analysis of the global regulator Lae1 uncovers a connection between Lae1 and the histone acetyltransferase HAT1 in Fusarium fujikuroi

The fungus Fusarium fujikuroi causes "bakanae" disease of rice due to its ability to produce gibberellins (GAs), a family of plant hormones. Recent genome sequencing revealed the genetic capacity for the biosynthesis of 46 additional secondary metabolites besides the industrially produced GAs. Among them are the pigments bikaverin and fusarubins, as well as mycotoxins, such as fumonisins, fusarin C, beauvericin, and fusaric acid. However, half of the potential secondary metabolite gene clusters are silent. In recent years, it has been shown that the fungal specific velvet complex is involved in global regulation of secondary metabolism in several filamentous fungi. We have previously shown that deletion of the three components of the F. fujikuroi velvet complex, vel1, vel2, and lae1, almost totally abolished biosynthesis of GAs, fumonisins and fusarin C. Here, we present a deeper insight into the genome-wide regulatory impact of Lae1 on secondary metabolism. Over-expression of lae1 resulted in de-repression of GA biosynthetic genes under otherwise repressing high nitrogen conditions demonstrating that the nitrogen repression is overcome. In addition, over-expression of one of five tested histone acetyltransferase genes, HAT1, was capable of returning GA gene expression and GA production to the GA-deficient Δlae1 mutant. Deletion and over-expression of HAT1 in the wild type resulted in downregulation and upregulation of GA gene expression, respectively, indicating that HAT1 together with Lae1 plays an essential role in the regulation of GA biosynthesis.

Comparative genomics of geographically distant Fusarium fujikuroi isolates revealed two distinct pathotypes correlating with secondary metabolite profiles

Fusarium fujikuroi causes bakanae ("foolish seedling") disease of rice which is characterized by hyper-elongation of seedlings resulting from production of gibberellic acids (GAs) by the fungus. This plant pathogen is also known for production of harmful mycotoxins, such as fusarins, fusaric acid, apicidin F and beauvericin. Recently, we generated the first de novo genome sequence of F. fujikuroi strain IMI 58289 combined with extensive transcriptional, epigenetic, proteomic and chemical product analyses. GA production was shown to provide a selective advantage during infection of the preferred host plant rice. Here, we provide genome sequences of eight additional F. fujikuroi isolates from distant geographic regions. The isolates differ in the size of chromosomes, most likely due to variability of subtelomeric regions, the type of asexual spores (microconidia and/or macroconidia), and the number and expression of secondary metabolite gene clusters. Whilst most of the isolates caused the typical bakanae symptoms, one isolate, B14, caused stunting and early withering of infected seedlings. In contrast to the other isolates, B14 produced no GAs but high amounts of fumonisins during infection on rice. Furthermore, it differed from the other isolates by the presence of three additional polyketide synthase (PKS) genes (PKS40, PKS43, PKS51) and the absence of the F. fujikuroi-specific apicidin F (NRPS31) gene cluster. Analysis of additional field isolates confirmed the strong correlation between the pathotype (bakanae or stunting/withering), and the ability to produce either GAs or fumonisins. Deletion of the fumonisin and fusaric acid-specific PKS genes in B14 reduced the stunting/withering symptoms, whereas deletion of the PKS51 gene resulted in elevated symptom development. Phylogenetic analyses revealed two subclades of F. fujikuroi strains according to their pathotype and secondary metabolite profiles.

The GATA-Type Transcription Factor Csm1 Regulates Conidiation and Secondary Metabolism in Fusarium fujikuroi

GATA-type transcription factors (TFs) such as the nitrogen regulators AreA and AreB, or the light-responsive TFs WC-1 and WC-2, play global roles in fungal growth and development. The conserved GATA TF NsdD is known as an activator of sexual development and key repressor of conidiation in Aspergillus nidulans, and as light-regulated repressor of macroconidia formation in Botrytis cinerea. In the present study, we functionally characterized the NsdD ortholog in Fusarium fujikuroi, named Csm1. Deletion of this gene resulted in elevated microconidia formation in the wild-type (WT) and restoration of conidiation in the non-sporulating velvet mutant Δvel1 demonstrating that Csm1 also plays a role as repressor of conidiation in F. fujikuroi. Furthermore, biosynthesis of the PKS-derived red pigments, bikaverin and fusarubins, is de-regulated under otherwise repressing conditions. Cross-species complementation of the Δcsm1 mutant with the B. cinerea ortholog LTF1 led to full restoration of WT-like growth, conidiation and pigment formation. In contrast, the F. fujikuroi CSM1 rescued only the defects in growth, the tolerance to H₂O₂ and virulence, but did not restore the light-dependent differentiation when expressed in the B. cinerea Δltf1 mutant. Microarray analysis comparing the expression profiles of the F. fujikuroi WT and the Δcsm1 mutant under different nitrogen conditions revealed a strong impact of this GATA TF on 19 of the 47 gene clusters in the genome of F. fujikuroi. One of the up-regulated silent gene clusters is the one containing the sesquiterpene cyclase-encoding key gene STC1. Heterologous expression of STC1 in Escherichia coli enabled us to identify the product as the volatile bioactive compound (-)-germacrene D.

Genome sequence of the highly weak-acid-tolerant Zygosaccharomyces bailii IST302, amenable to genetic manipulations and physiological studies

Zygosaccharomyces bailii is one of the most problematic spoilage yeast species found in the food and beverage industry particularly in acidic products, due to its exceptional resistance to weak acid stress. This article describes the annotation of the genome sequence of Z. bailii IST302, a strain recently proven to be amenable to genetic manipulations and physiological studies. The work was based on the annotated genomes of strain ISA1307, an interspecies hybrid between Z. bailii and a closely related species, and the Z. bailii reference strain CLIB 213T. The resulting genome sequence of Z. bailii IST302 is distributed through 105 scaffolds, comprising a total of 5142 genes and a size of 10.8 Mb. Contrasting with CLIB 213T, strain IST302 does not form cell aggregates, allowing its manipulation in the laboratory for genetic and physiological studies. Comparative cell cycle analysis with the haploid and diploid Saccharomyces cerevisiae strains BY4741 and BY4743, respectively, suggests that Z. bailii IST302 is haploid. This is an additional trait that makes this strain attractive for the functional analysis of non-essential genes envisaging the elucidation of mechanisms underlying its high tolerance to weak acid food preservatives, or the investigation and exploitation of the potential of this resilient yeast species as cell factory.

Comparative transcriptome and proteome analysis reveals a global impact of the nitrogen regulators AreA and AreB on secondary metabolism in Fusarium fujikuroi

The biosynthesis of multiple secondary metabolites in the phytopathogenic ascomycete Fusarium fujikuroi is strongly affected by nitrogen availability. Here, we present the first genome-wide transcriptome and proteome analysis that compared the wild type and deletion mutants of the two major nitrogen regulators AreA and AreB. We show that AreB acts not simply as an antagonist of AreA counteracting the expression of AreA target genes as suggested based on the yeast model. Both GATA transcription factors affect a large and diverse set of common as well as specific target genes and proteins, acting as activators and repressors. We demonstrate that AreA and AreB are not only involved in fungal nitrogen metabolism, but also in the control of several complex cellular processes like carbon metabolism, transport and secondary metabolism. We show that both GATA transcription factors can be considered as master regulators of secondary metabolism as they affect the expression of more than half of the 47 putative secondary metabolite clusters identified in the genome of F. fujikuroi. While AreA acts as a positive regulator of many clusters under nitrogen-limiting conditions, AreB is able to activate and repress gene clusters (e.g. bikaverin) under nitrogen limitation and sufficiency. In addition, ChIP analyses revealed that loss of AreA or AreB causes histone modifications at some of the regulated gene clusters.

Comparative "Omics" of the Fusarium fujikuroi Species Complex Highlights Differences in Genetic Potential and Metabolite Synthesis

Species of the Fusarium fujikuroi species complex (FFC) cause a wide spectrum of often devastating diseases on diverse agricultural crops, including coffee, fig, mango, maize, rice, and sugarcane. Although species within the FFC are difficult to distinguish by morphology, and their genes often share 90% sequence similarity, they can differ in host plant specificity and life style. FFC species can also produce structurally diverse secondary metabolites (SMs), including the mycotoxins fumonisins, fusarins, fusaric acid, and beauvericin, and the phytohormones gibberellins, auxins, and cytokinins. The spectrum of SMs produced can differ among closely related species, suggesting that SMs might be determinants of host specificity. To date, genomes of only a limited number of FFC species have been sequenced. Here, we provide draft genome sequences of three more members of the FFC: a single isolate of F. mangiferae, the cause of mango malformation, and two isolates of F. proliferatum, one a pathogen of maize and the other an orchid endophyte. We compared these genomes to publicly available genome sequences of three other FFC species. The comparisons revealed species-specific and isolate-specific differences in the composition and expression (in vitro and in planta) of genes involved in SM production including those for phytohormome biosynthesis. Such differences have the potential to impact host specificity and, as in the case of F. proliferatum, the pathogenic versus endophytic life style.

Globally distributed root endophyte Phialocephala subalpina links pathogenic and saprophytic lifestyles

BACKGROUND

Whereas an increasing number of pathogenic and mutualistic ascomycetous species were sequenced in the past decade, species showing a seemingly neutral association such as root endophytes received less attention. In the present study, the genome of Phialocephala subalpina, the most frequent species of the Phialocephala fortinii s.l. - Acephala applanata species complex, was sequenced for insight in the genome structure and gene inventory of these wide-spread root endophytes.

RESULTS

The genome of P. subalpina was sequenced using Roche/454 GS FLX technology and a whole genome shotgun strategy. The assembly resulted in 205 scaffolds and a genome size of 69.7 Mb. The expanded genome size in P. subalpina was not due to the proliferation of transposable elements or other repeats, as is the case with other ascomycetous genomes. Instead, P. subalpina revealed an expanded gene inventory that includes 20,173 gene models. Comparative genome analysis of P. subalpina with 13 ascomycetes shows that P. subalpina uses a versatile gene inventory including genes specific for pathogens and saprophytes. Moreover, the gene inventory for carbohydrate active enzymes (CAZymes) was expanded including genes involved in degradation of biopolymers, such as pectin, hemicellulose, cellulose and lignin.

CONCLUSIONS

The analysis of a globally distributed root endophyte allowed detailed insights in the gene inventory and genome organization of a yet largely neglected group of organisms. We showed that the ubiquitous root endophyte P. subalpina has a broad gene inventory that links pathogenic and saprophytic lifestyles.

Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses

BACKGROUND

The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts.

RESULTS

The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type.

CONCLUSION

Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

A complete toolset for the study of Ustilago bromivora and Brachypodium sp. as a fungal-temperate grass pathosystem

Due to their economic relevance, the study of plant pathogen interactions is of importance. However, elucidating these interactions and their underlying molecular mechanisms remains challenging since both host and pathogen need to be fully genetically accessible organisms. Here we present milestones in the establishment of a new biotrophic model pathosystem: Ustilago bromivora and Brachypodium sp. We provide a complete toolset, including an annotated fungal genome and methods for genetic manipulation of the fungus and its host plant. This toolset will enable researchers to easily study biotrophic interactions at the molecular level on both the pathogen and the host side. Moreover, our research on the fungal life cycle revealed a mating type bias phenomenon. U. bromivora harbors a haplo-lethal allele that is linked to one mating type region. As a result, the identified mating type bias strongly promotes inbreeding, which we consider to be a potential speciation driver.

Comparison of Fusarium graminearum Transcriptomes on Living or Dead Wheat Differentiates Substrate-Responsive and Defense-Responsive Genes

Fusarium graminearum is an opportunistic pathogen of cereals where it causes severe yield losses and concomitant mycotoxin contamination of the grains. The pathogen has mixed biotrophic and necrotrophic (saprophytic) growth phases during infection and the regulatory networks associated with these phases have so far always been analyzed together. In this study we compared the transcriptomes of fungal cells infecting a living, actively defending plant representing the mixed live style (pathogenic growth on living flowering wheat heads) to the response of the fungus infecting identical, but dead plant tissues (cold-killed flowering wheat heads) representing strictly saprophytic conditions. We found that the living plant actively suppressed fungal growth and promoted much higher toxin production in comparison to the identical plant tissue without metabolism suggesting that molecules signaling secondary metabolite induction are not pre-existing or not stable in the plant in sufficient amounts before infection. Differential gene expression analysis was used to define gene sets responding to the active or the passive plant as main impact factor and driver for gene expression. We correlated our results to the published F. graminearum transcriptomes, proteomes, and secretomes and found that only a limited number of in planta- expressed genes require the living plant for induction but the majority uses simply the plant tissue as signal. Many secondary metabolite (SM) gene clusters show a heterogeneous expression pattern within the cluster indicating that different genetic or epigenetic signals govern the expression of individual genes within a physically linked cluster. Our bioinformatic approach also identified fungal genes which were actively repressed by signals derived from the active plant and may thus represent direct targets of the plant defense against the invading pathogen.

A Tale of Genome Compartmentalization: The Evolution of Virulence Clusters in Smut Fungi

Smut fungi are plant pathogens mostly parasitizing wild species of grasses as well as domesticated cereal crops. Genome analysis of several smut fungi including Ustilago maydis revealed a singular clustered organization of genes encoding secreted effectors. In U. maydis, many of these clusters have a role in virulence. Reconstructing the evolutionary history of clusters of effector genes is difficult because of their intrinsically fast evolution, which erodes the phylogenetic signal and homology relationships. Here, we describe the use of comparative evolutionary analyses of quality draft assemblies of genomes to study the mechanisms of this evolution. We report the genome sequence of a South African isolate of Sporisorium scitamineum, a smut fungus parasitizing sugar cane with a phylogenetic position intermediate to the two previously sequenced species U. maydis and Sporisorium reilianum. We show that the genome of S. scitamineum contains more and larger gene clusters encoding secreted effectors than any previously described species in this group. We trace back the origin of the clusters and find that their evolution is mainly driven by tandem gene duplication. In addition, transposable elements play a major role in the evolution of the clustered genes. Transposable elements are significantly associated with clusters of genes encoding fast evolving secreted effectors. This suggests that such clusters represent a case of genome compartmentalization that restrains the activity of transposable elements on genes under diversifying selection for which this activity is potentially beneficial, while protecting the rest of the genome from its deleterious effect.

The global regulator FfSge1 is required for expression of secondary metabolite gene clusters but not for pathogenicity in Fusarium fujikuroi

The plant pathogenic fungus Fusarium fujikuroi is the causal agent of bakanae disease on rice due to its ability to produce gibberellins. Besides these phytohormones, F. fujikuroi is able to produce several other secondary metabolites (SMs). Although much progress has been made in the field of secondary metabolism, the transcriptional regulation of SM biosynthesis is complex and still incompletely understood. Environmental conditions, global as well as pathway-specific regulators and chromatin remodelling have been shown to play major roles. Here, the role of FfSge1, a homologue of the morphological switch regulators Wor1 and Ryp1 in Candida albicans and Histoplasma capsulatum, respectively, is explored with emphasis on secondary metabolism. FfSge1 is not required for formation of conidia and pathogenicity but is involved in vegetative growth. Transcriptome analysis of the mutant Δffsge1 compared with the wild type, as well as comparative chemical analysis between the wild type, Δffsge1 and OE:FfSGE1, revealed that FfSge1 functions as a global activator of secondary metabolism in F. fujikuroi. Double mutants of FfSGE1 and other SM regulatory genes brought insights into the hierarchical regulation of secondary metabolism. In addition, FfSge1 is also required for expression of a yet uncharacterized SM gene cluster containing a non-canonical non-ribosomal peptide synthetase.

The genome sequence of the highly acetic acid-tolerant Zygosaccharomyces bailii-derived interspecies hybrid strain ISA1307, isolated from a sparkling wine plant

In this work, it is described the sequencing and annotation of the genome of the yeast strain ISA1307, isolated from a sparkling wine continuous production plant. This strain, formerly considered of the Zygosaccharomyces bailii species, has been used to study Z. bailii physiology, in particular, its extreme tolerance to acetic acid stress at low pH. The analysis of the genome sequence described in this work indicates that strain ISA1307 is an interspecies hybrid between Z. bailii and a closely related species. The genome sequence of ISA1307 is distributed through 154 scaffolds and has a size of around 21.2 Mb, corresponding to 96% of the genome size estimated by flow cytometry. Annotation of ISA1307 genome includes 4385 duplicated genes (∼90% of the total number of predicted genes) and 1155 predicted single-copy genes. The functional categories including a higher number of genes are ‘Metabolism and generation of energy’, ‘Protein folding, modification and targeting’ and ‘Biogenesis of cellular components’. The knowledge of the genome sequence of the ISA1307 strain is expected to contribute to accelerate systems-level understanding of stress resistance mechanisms in Z. bailii and to inspire and guide novel biotechnological applications of this yeast species/strain in fermentation processes, given its high resilience to acidic stress. The availability of the ISA1307 genome sequence also paves the way to a better understanding of the genetic mechanisms underlying the generation and selection of more robust hybrid yeast strains in the stressful environment of wine fermentations.

Genetic manipulation of the Fusarium fujikuroi fusarin gene cluster yields insight into the complex regulation and fusarin biosynthetic pathway

In this work, the biosynthesis and regulation of the polyketide synthase/nonribosomal peptide synthetase (PKS/NRPS)-derived mutagenic mycotoxin fusarin C was studied in the fungus Fusarium fujikuroi. The fusarin gene cluster consists of nine genes (fus1-fus9) that are coexpressed under high-nitrogen and acidic pH conditions. Chromatin immunoprecipitation revealed a correlation between high expression and enrichment of activating H3K9-acetylation marks under inducing conditions. We provide evidence that only four genes are sufficient for the biosynthesis. The combination of genetic engineering with nuclear magnetic resonance and mass-spectrometry-based structure elucidation allowed the discovery of the putative fusarin biosynthetic pathway. Surprisingly, we indicate that PKS/NRPS releases its product with an open ring structure, probably as an alcohol. Our data indicate that 2-pyrrolidone ring closure, oxidation at C-20, and, finally, methylation at C-20 are catalyzed by Fus2, Fus8, and Fus9, respectively.

Deciphering the cryptic genome: genome-wide analyses of the rice pathogen Fusarium fujikuroi reveal complex regulation of secondary metabolism and novel metabolites

The fungus Fusarium fujikuroi causes "bakanae" disease of rice due to its ability to produce gibberellins (GAs), but it is also known for producing harmful mycotoxins. However, the genetic capacity for the whole arsenal of natural compounds and their role in the fungus' interaction with rice remained unknown. Here, we present a high-quality genome sequence of F. fujikuroi that was assembled into 12 scaffolds corresponding to the 12 chromosomes described for the fungus. We used the genome sequence along with ChIP-seq, transcriptome, proteome, and HPLC-FTMS-based metabolome analyses to identify the potential secondary metabolite biosynthetic gene clusters and to examine their regulation in response to nitrogen availability and plant signals. The results indicate that expression of most but not all gene clusters correlate with proteome and ChIP-seq data. Comparison of the F. fujikuroi genome to those of six other fusaria revealed that only a small number of gene clusters are conserved among these species, thus providing new insights into the divergence of secondary metabolism in the genus Fusarium. Noteworthy, GA biosynthetic genes are present in some related species, but GA biosynthesis is limited to F. fujikuroi, suggesting that this provides a selective advantage during infection of the preferred host plant rice. Among the genome sequences analyzed, one cluster that includes a polyketide synthase gene (PKS19) and another that includes a non-ribosomal peptide synthetase gene (NRPS31) are unique to F. fujikuroi. The metabolites derived from these clusters were identified by HPLC-FTMS-based analyses of engineered F. fujikuroi strains overexpressing cluster genes. In planta expression studies suggest a specific role for the PKS19-derived product during rice infection. Thus, our results indicate that combined comparative genomics and genome-wide experimental analyses identified novel genes and secondary metabolites that contribute to the evolutionary success of F. fujikuroi as a rice pathogen.

Plant-symbiotic fungi as chemical engineers: multi-genome analysis of the clavicipitaceae reveals dynamics of alkaloid loci

The fungal family Clavicipitaceae includes plant symbionts and parasites that produce several psychoactive and bioprotective alkaloids. The family includes grass symbionts in the epichloae clade (Epichloë and Neotyphodium species), which are extraordinarily diverse both in their host interactions and in their alkaloid profiles. Epichloae produce alkaloids of four distinct classes, all of which deter insects, and some-including the infamous ergot alkaloids-have potent effects on mammals. The exceptional chemotypic diversity of the epichloae may relate to their broad range of host interactions, whereby some are pathogenic and contagious, others are mutualistic and vertically transmitted (seed-borne), and still others vary in pathogenic or mutualistic behavior. We profiled the alkaloids and sequenced the genomes of 10 epichloae, three ergot fungi (Claviceps species), a morning-glory symbiont (Periglandula ipomoeae), and a bamboo pathogen (Aciculosporium take), and compared the gene clusters for four classes of alkaloids. Results indicated a strong tendency for alkaloid loci to have conserved cores that specify the skeleton structures and peripheral genes that determine chemical variations that are known to affect their pharmacological specificities. Generally, gene locations in cluster peripheries positioned them near to transposon-derived, AT-rich repeat blocks, which were probably involved in gene losses, duplications, and neofunctionalizations. The alkaloid loci in the epichloae had unusual structures riddled with large, complex, and dynamic repeat blocks. This feature was not reflective of overall differences in repeat contents in the genomes, nor was it characteristic of most other specialized metabolism loci. The organization and dynamics of alkaloid loci and abundant repeat blocks in the epichloae suggested that these fungi are under selection for alkaloid diversification. We suggest that such selection is related to the variable life histories of the epichloae, their protective roles as symbionts, and their associations with the highly speciose and ecologically diverse cool-season grasses.

Genome comparison of barley and maize smut fungi reveals targeted loss of RNA silencing components and species-specific presence of transposable elements

Ustilago hordei is a biotrophic parasite of barley (Hordeum vulgare). After seedling infection, the fungus persists in the plant until head emergence when fungal spores develop and are released from sori formed at kernel positions. The 26.1-Mb U. hordei genome contains 7113 protein encoding genes with high synteny to the smaller genomes of the related, maize-infecting smut fungi Ustilago maydis and Sporisorium reilianum but has a larger repeat content that affected genome evolution at important loci, including mating-type and effector loci. The U. hordei genome encodes components involved in RNA interference and heterochromatin formation, normally involved in genome defense, that are lacking in the U. maydis genome due to clean excision events. These excision events were possibly a result of former presence of repetitive DNA and of an efficient homologous recombination system in U. maydis. We found evidence of repeat-induced point mutations in the genome of U. hordei, indicating that smut fungi use different strategies to counteract the deleterious effects of repetitive DNA. The complement of U. hordei effector genes is comparable to the other two smuts but reveals differences in family expansion and clustering. The availability of the genome sequence will facilitate the identification of genes responsible for virulence and evolution of smut fungi on their respective hosts.

Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica

Recent sequencing projects have provided deep insight into fungal lifestyle-associated genomic adaptations. Here we report on the 25 Mb genome of the mutualistic root symbiont Piriformospora indica (Sebacinales, Basidiomycota) and provide a global characterization of fungal transcriptional responses associated with the colonization of living and dead barley roots. Extensive comparative analysis of the P. indica genome with other Basidiomycota and Ascomycota fungi that have diverse lifestyle strategies identified features typically associated with both, biotrophism and saprotrophism. The tightly controlled expression of the lifestyle-associated gene sets during the onset of the symbiosis, revealed by microarray analysis, argues for a biphasic root colonization strategy of P. indica. This is supported by a cytological study that shows an early biotrophic growth followed by a cell death-associated phase. About 10% of the fungal genes induced during the biotrophic colonization encoded putative small secreted proteins (SSP), including several lectin-like proteins and members of a P. indica-specific gene family (DELD) with a conserved novel seven-amino acids motif at the C-terminus. Similar to effectors found in other filamentous organisms, the occurrence of the DELDs correlated with the presence of transposable elements in gene-poor repeat-rich regions of the genome. This is the first in depth genomic study describing a mutualistic symbiont with a biphasic lifestyle. Our findings provide a significant advance in understanding development of biotrophic plant symbionts and suggest a series of incremental shifts along the continuum from saprotrophy towards biotrophy in the evolution of mycorrhizal association from decomposer fungi.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Pathogenicity determinants in smut fungi revealed by genome comparison

Biotrophic pathogens, such as the related maize pathogenic fungi Ustilago maydis and Sporisorium reilianum, establish an intimate relationship with their hosts by secreting protein effectors. Because secreted effectors interacting with plant proteins should rapidly evolve, we identified variable genomic regions by sequencing the genome of S. reilianum and comparing it with the U. maydis genome. We detected 43 regions of low sequence conservation in otherwise well-conserved syntenic genomes. These regions primarily encode secreted effectors and include previously identified virulence clusters. By deletion analysis in U. maydis, we demonstrate a role in virulence for four previously unknown diversity regions. This highlights the power of comparative genomics of closely related species for identification of virulence determinants.

FGDB: revisiting the genome annotation of the plant pathogen Fusarium graminearum

The MIPS Fusarium graminearum Genome Database (FGDB) was established as a comprehensive genome database on one of the most devastating fungal plant pathogens of wheat, barley and maize. The current version of FGDB v3.1 provides information on the full manually revised gene set based on the Broad Institute assembly FG3 genome sequence. The results of gene prediction tools were integrated with the help of comparative data on related species to result in a set of 13.718 annotated protein coding genes. This rigorous approach involved adding or modifying gene models and represents a coding sequence gold standard for the genus Fusarium. The gene loci improvements results in 2461 genes which either are new or have different structures compared to the Broad Institute assembly 3 gene set. Moreover the database serves as a convenient entry point to explore expression data results and to obtain information on the Affymetrix GeneChip probe sets. The resource is accessible on http://mips.gsf.de/genre/proj/FGDB/.

Transcriptome analysis of nitrate assimilation in Aspergillus nidulans reveals connections to nitric oxide metabolism

Nitrate is a dominant form of inorganic nitrogen (N) in soils and can be efficiently assimilated by bacteria, fungi and plants. We studied here the transcriptome of the short-term nitrate response using assimilating and non-assimilating strains of the model ascomycete Aspergillus nidulans. Among the 72 genes positively responding to nitrate, only 18 genes carry binding sites for the pathway-specific activator NirA. Forty-five genes were repressed by nitrate metabolism. Because nirA(-) strains are N-starved at nitrate induction conditions, we also compared the nitrate transcriptome with N-deprived conditions and found a partial overlap of differentially regulated genes between these conditions. Nitric oxide (NO)-metabolizing flavohaemoglobins were found to be co-regulated with nitrate assimilatory genes. Subsequent molecular characterization revealed that the strongly inducible FhbA is required for full activity of nitrate and nitrite reductase enzymes. The co-regulation of NO-detoxifying and nitrate/nitrite assimilating systems may represent a conserved mechanism, which serves to neutralize nitrosative stress imposed by an external NO source in saprophytic and pathogenic fungi. Our analysis using membrane-permeable NO donors suggests that signalling for NirA activation only indirectly depends on the nitrate transporters NrtA (CrnA) and NrtB (CrnB).

PEDANT covers all complete RefSeq genomes

The PEDANT genome database provides exhaustive annotation of nearly 3000 publicly available eukaryotic, eubacterial, archaeal and viral genomes with more than 4.5 million proteins by a broad set of bioinformatics algorithms. In particular, all completely sequenced genomes from the NCBI's Reference Sequence collection (RefSeq) are covered. The PEDANT processing pipeline has been sped up by an order of magnitude through the utilization of precalculated similarity information stored in the similarity matrix of proteins (SIMAP) database, making it possible to process newly sequenced genomes immediately as they become available. PEDANT is freely accessible to academic users at http://pedant.gsf.de. For programmatic access Web Services are available at http://pedant.gsf.de/webservices.jsp.

The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization

We sequenced and annotated the genome of the filamentous fungus Fusarium graminearum, a major pathogen of cultivated cereals. Very few repetitive sequences were detected, and the process of repeat-induced point mutation, in which duplicated sequences are subject to extensive mutation, may partially account for the reduced repeat content and apparent low number of paralogous (ancestrally duplicated) genes. A second strain of F. graminearum contained more than 10,000 single-nucleotide polymorphisms, which were frequently located near telomeres and within other discrete chromosomal segments. Many highly polymorphic regions contained sets of genes implicated in plant-fungus interactions and were unusually divergent, with higher rates of recombination. These regions of genome innovation may result from selection due to interactions of F. graminearum with its plant hosts.

Gene expression shifts during perithecium development in Gibberella zeae (anamorph Fusarium graminearum), with particular emphasis on ion transport proteins

Gibberella zeae, the causal agent of Fusarium head blight, is a devastating pathogen of small grains worldwide. The sexual cycle is a crucial component of head blight epidemiology, as forcibly discharged ascospores serve as the primary inoculum. The recent development of an Affymetrix GeneChip containing probesets representative of all predicted genes of G. zeae has opened the door to studies of differential gene expression during sexual development. Using GeneChips, a developmental time course was performed in culture, from vegetative hyphae to mature perithecia with multiseptate ascospores. Time-points represent the development of the major cell types comprising the mature perithecium. The majority of the 17,830 G. zeae probesets, 78%, were expressed during at least one of the developmental stages; 12% of these appear to be specific to sexual development. Analysis of the 162 predicted ion transporter genes is reported in detail, due to their association with perithecium function. Expression patterns of the MirA-type siderophores, chloride channels, P-type ATPases and potassium transporters show some specialization in regard to developmental stage. This is the first whole-genome analysis of differential transcript accumulation during sexual development in a filamentous fungus.

Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis

Ustilago maydis is a ubiquitous pathogen of maize and a well-established model organism for the study of plant-microbe interactions. This basidiomycete fungus does not use aggressive virulence strategies to kill its host. U. maydis belongs to the group of biotrophic parasites (the smuts) that depend on living tissue for proliferation and development. Here we report the genome sequence for a member of this economically important group of biotrophic fungi. The 20.5-million-base U. maydis genome assembly contains 6,902 predicted protein-encoding genes and lacks pathogenicity signatures found in the genomes of aggressive pathogenic fungi, for example a battery of cell-wall-degrading enzymes. However, we detected unexpected genomic features responsible for the pathogenicity of this organism. Specifically, we found 12 clusters of genes encoding small secreted proteins with unknown function. A significant fraction of these genes exists in small gene families. Expression analysis showed that most of the genes contained in these clusters are regulated together and induced in infected tissue. Deletion of individual clusters altered the virulence of U. maydis in five cases, ranging from a complete lack of symptoms to hypervirulence. Despite years of research into the mechanism of pathogenicity in U. maydis, no 'true' virulence factors had been previously identified. Thus, the discovery of the secreted protein gene clusters and the functional demonstration of their decisive role in the infection process illuminate previously unknown mechanisms of pathogenicity operating in biotrophic fungi. Genomic analysis is, similarly, likely to open up new avenues for the discovery of virulence determinants in other pathogens.

Development of a Fusarium graminearum Affymetrix GeneChip for profiling fungal gene expression in vitro and in planta

Recently the genome sequences of several filamentous fungi have become available, providing the opportunity for large-scale functional analysis including genome-wide expression analysis. We report the design and validation of the first Affymetrix GeneChip microarray based on the entire genome of a filamentous fungus, the ascomycetous plant pathogen Fusarium graminearum. To maximize the likelihood of representing all putative genes (approximately 14,000) on the array, two distinct sets of automatically predicted gene calls were used and integrated into the online F. graminearum Genome DataBase. From these gene sets, a subset of calls was manually annotated and a non-redundant extract of all calls together with additional EST sequences and controls were submitted for GeneChip design. Experiments were conducted to test the performance of the F. graminearum GeneChip. Hybridization experiments using genomic DNA demonstrated the usefulness of the array for experimentation with F. graminearum and at least four additional pathogenic Fusarium species. Differential transcript accumulation was detected using the F. graminearum GeneChip with treatments derived from the fungus grown in culture under three nutritional regimes and in comparison with fungal growth in infected barley. The ability to detect fungal genes in planta is surprisingly sensitive even without efforts to enrich for fungal transcripts. The Plant Expression Database (PLEXdb, http://www.plexdb.org) will be used as a public repository for raw and normalized expression data from the F. graminearum GeneChip. The F. graminearum GeneChip will help to accelerate exploration of the pathogen-host pathways that may involve interactions between pathogenicity genes in the fungus and disease response in the plant.

FGDB: a comprehensive fungal genome resource on the plant pathogen Fusarium graminearum

The MIPS Fusarium graminearum Genome Database (FGDB) is a comprehensive genome database on one of the most devastating fungal plant pathogens of wheat and barley. FGDB provides information on two gene sets independently derived by automated annotation of the F.graminearum genome sequence. A complete manually revised gene set will be completed within the near future. The initial results of systematic manual correction of gene calls are already part of the current gene set. The database can be accessed to retrieve information from bioinformatics analyses and functional classifications of the proteins. The data are also organized in the well established MIPS catalogs and novel query techniques are available to search the data. The comprehensive set of gene calls was also used for the design of an Affymetrix GeneChip. The resource is accessible on .

MPact: the MIPS protein interaction resource on yeast

In recent years, the Munich Information Center for Protein Sequences (MIPS) yeast protein–protein interaction (PPI) dataset has been used in numerous analyses of protein networks and has been called a gold standard because of its quality and comprehensiveness [H. Yu, N. M. Luscombe, H. X. Lu, X. Zhu, Y. Xia, J. D. Han, N. Bertin, S. Chung, M. Vidal and M. Gerstein (2004) Genome Res., 14, 1107–1118]. MPact and the yeast protein localization catalog provide information related to the proximity of proteins in yeast. Beside the integration of high-throughput data, information about experimental evidence for PPIs in the literature was compiled by experts adding up to 4300 distinct PPIs connecting 1500 proteins in yeast. As the interaction data is a complementary part of CYGD, interactive mapping of data on other integrated data types such as the functional classification catalog [A. Ruepp, A. Zollner, D. Maier, K. Albermann, J. Hani, M. Mokrejs, I. Tetko, U. Güldener, G. Mannhaupt, M. Münsterkötter and H. W. Mewes (2004) Nucleic Acids Res., 32, 5539–5545] is possible. A survey of signaling proteins and comparison with pathway data from KEGG demonstrates that based on these manually annotated data only an extensive overview of the complexity of this functional network can be obtained in yeast. The implementation of a web-based PPI-analysis tool allows analysis and visualization of protein interaction networks and facilitates integration of our curated data with high-throughput datasets. The complete dataset as well as user-defined sub-networks can be retrieved easily in the standardized PSI-MI format. The resource can be accessed through .