Candidate Effector Proteins of the Rust Pathogen Melampsora larici-populina Target Diverse Plant Cell Compartments

Rust fungi are devastating crop pathogens that deliver effector proteins into infected tissues to modulate plant functions and promote parasitic growth. The genome of the poplar leaf rust fungus Melampsora larici-populina revealed a large catalog of secreted proteins, some of which have been considered candidate effectors. Unraveling how these proteins function in host cells is a key to understanding pathogenicity mechanisms and developing resistant plants. In this study, we used an effectoromics pipeline to select, clone, and express 20 candidate effectors in Nicotiana benthamiana leaf cells to determine their subcellular localization and identify the plant proteins they interact with. Confocal microscopy revealed that six candidate effectors target the nucleus, nucleoli, chloroplasts, mitochondria, and discrete cellular bodies. We also used coimmunoprecipitation (coIP) and mass spectrometry to identify 606 N. benthamiana proteins that associate with the candidate effectors. Five candidate effectors specifically associated with a small set of plant proteins that may represent biologically relevant interactors. We confirmed the interaction between the candidate effector MLP124017 and TOPLESS-related protein 4 from poplar by in planta coIP. Altogether, our data enable us to validate effector proteins from M. larici-populina and reveal that these proteins may target multiple compartments and processes in plant cells. It also shows that N. benthamiana can be a powerful heterologous system to study effectors of obligate biotrophic pathogens.

The General Amino Acid Permease FfGap1 of Fusarium fujikuroi Is Sorted to the Vacuole in a Nitrogen-Dependent, but Npr1 Kinase-Independent Manner

The rice pathogenic fungus Fusarium fujikuroi is well known for the production of a broad spectrum of secondary metabolites (SMs) such as gibberellic acids (GAs), mycotoxins and pigments. The biosynthesis of most of these SMs strictly depends on nitrogen availability and of the activity of permeases of nitrogen sources, e.g. the ammonium and amino acid permeases. One of the three ammonium permeases, MepB, was recently shown to act not only as a transporter but also as a nitrogen sensor affecting the production of nitrogen-repressed SMs. Here we describe the identification of a general amino acid permease, FfGap1, among the 99 putative amino acid permeases (AAPs) in the genome of F. fujikuroi. FfGap1 is able to fully restore growth of the yeast gap1∆ mutant on several amino acids including citrulline and tryptophane. In S. cerevisiae, Gap1 activity is regulated by shuttling between the plasma membrane (nitrogen limiting conditions) and the vacuole (nitrogen sufficiency), which we also show for FfGap1. In yeast, the Npr1 serine/threonine kinase stabilizes the Gap1 position at the plasma membrane. Here, we identified and characterized three NPR1-homologous genes, encoding the putative protein kinases FfNpr1-1, FfNpr1-2 and FfNpr1-3 with significant similarity to yeast Npr1. Complementation of the yeast npr1Δ mutant with each of the three F. fujikuroi NPR1 homologues, resulted in partial restoration of ammonium, arginine and proline uptake by FfNPR1-1 while none of the three kinases affect growth on different nitrogen sources and nitrogen-dependent sorting of FfGap1 in F. fujikuroi. However, exchange of the putative ubiquitin-target lysine 9 (K9A) and 15 (K15A) residues of FfGap1 resulted in extended localization to the plasma membrane and increased protein stability independently of nitrogen availability. These data suggest a similar regulation of FfGap1 by nitrogen-dependent ubiquitination, but differences regarding the role of Fusarium Npr1 homologues compared to yeast.

Global mapping of DNA conformational flexibility on Saccharomyces cerevisiae

In this study we provide the first comprehensive map of DNA conformational flexibility in Saccharomyces cerevisiae complete genome. Flexibility plays a key role in DNA supercoiling and DNA/protein binding, regulating DNA transcription, replication or repair. Specific interest in flexibility analysis concerns its relationship with human genome instability. Enrichment in flexible sequences has been detected in unstable regions of human genome defined fragile sites, where genes map and carry frequent deletions and rearrangements in cancer. Flexible sequences have been suggested to be the determinants of fragile gene proneness to breakage; however, their actual role and properties remain elusive. Our in silico analysis carried out genome-wide via the StabFlex algorithm, shows the conserved presence of highly flexible regions in budding yeast genome as well as in genomes of other Saccharomyces sensu stricto species. Flexibile peaks in S. cerevisiae identify 175 ORFs mapping on their 3'UTR, a region affecting mRNA translation, localization and stability. (TA)n repeats of different extension shape the central structure of peaks and co-localize with polyadenylation efficiency element (EE) signals. ORFs with flexible peaks share common features. Transcripts are characterized by decreased half-life: this is considered peculiar of genes involved in regulatory systems with high turnover; consistently, their function affects biological processes such as cell cycle regulation or stress response. Our findings support the functional importance of flexibility peaks, suggesting that the flexible sequence may be derived by an expansion of canonical TAYRTA polyadenylation efficiency element. The flexible (TA)n repeat amplification could be the outcome of an evolutionary neofunctionalization leading to a differential 3'-end processing and expression regulation in genes with peculiar function. Our study provides a new support to the functional role of flexibility in genomes and a strategy for its characterization inside human fragile sites.

Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists

To elucidate the genetic bases of mycorrhizal lifestyle evolution, we sequenced new fungal genomes, including 13 ectomycorrhizal (ECM), orchid (ORM) and ericoid (ERM) species, and five saprotrophs, which we analyzed along with other fungal genomes. Ectomycorrhizal fungi have a reduced complement of genes encoding plant cell wall-degrading enzymes (PCWDEs), as compared to their ancestral wood decayers. Nevertheless, they have retained a unique array of PCWDEs, thus suggesting that they possess diverse abilities to decompose lignocellulose. Similar functional categories of nonorthologous genes are induced in symbiosis. Of induced genes, 7-38% are orphan genes, including genes that encode secreted effector-like proteins. Convergent evolution of the mycorrhizal habit in fungi occurred via the repeated evolution of a 'symbiosis toolkit', with reduced numbers of PCWDEs and lineage-specific suites of mycorrhiza-induced genes.

Horizontal gene transfer and gene dosage drives adaptation to wood colonization in a tree pathogen

Some of the most damaging tree pathogens can attack woody stems, causing lesions (cankers) that may be lethal. To identify the genomic determinants of wood colonization leading to canker formation, we sequenced the genomes of the poplar canker pathogen, Mycosphaerella populorum, and the closely related poplar leaf pathogen, M. populicola. A secondary metabolite cluster unique to M. populorum is fully activated following induction by poplar wood and leaves. In addition, genes encoding hemicellulose-degrading enzymes, peptidases, and metabolite transporters were more abundant and were up-regulated in M. populorum growing on poplar wood-chip medium compared with M. populicola. The secondary gene cluster and several of the carbohydrate degradation genes have the signature of horizontal transfer from ascomycete fungi associated with wood decay and from prokaryotes. Acquisition and maintenance of the gene battery necessary for growth in woody tissues and gene dosage resulting in gene expression reconfiguration appear to be responsible for the adaptation of M. populorum to infect, colonize, and cause mortality on poplar woody stems.

The Genome of the Saprophytic Fungus Verticillium tricorpus Reveals a Complex Effector Repertoire Resembling That of Its Pathogenic Relatives

Vascular wilts caused by Verticillium spp. are destructive plant diseases affecting hundreds of hosts. Only a few Verticillium spp. are causal agents of vascular wilt diseases, of which V. dahliae is the most notorious pathogen, and several V. dahliae genomes are available. In contrast, V. tricorpus is mainly known as a saprophyte and causal agent of opportunistic infections. Based on a hybrid approach that combines second and third generation sequencing, a near-gapless V. tricorpus genome assembly was obtained. With comparative genomics, we sought to identify genomic features in V. dahliae that confer the ability to cause vascular wilt disease. Unexpectedly, both species encode similar effector repertoires and share a genomic structure with genes encoding secreted proteins clustered in genomic islands. Intriguingly, V. tricorpus contains significantly fewer repetitive elements and an extended spectrum of secreted carbohydrate- active enzymes when compared with V. dahliae. In conclusion, we highlight the technical advances of a hybrid sequencing and assembly approach and show that the saprophyte V. tricorpus shares many hallmark features with the pathogen V. dahliae.

Genome, Transcriptome, and Functional Analyses of Penicillium expansum Provide New Insights Into Secondary Metabolism and Pathogenicity

The relationship between secondary metabolism and infection in pathogenic fungi has remained largely elusive. The genus Penicillium comprises a group of plant pathogens with varying host specificities and with the ability to produce a wide array of secondary metabolites. The genomes of three Penicillium expansum strains, the main postharvest pathogen of pome fruit, and one Pencillium italicum strain, a postharvest pathogen of citrus fruit, were sequenced and compared with 24 other fungal species. A genomic analysis of gene clusters responsible for the production of secondary metabolites was performed. Putative virulence factors in P. expansum were identified by means of a transcriptomic analysis of apple fruits during the course of infection. Despite a major genome contraction, P. expansum is the Penicillium species with the largest potential for the production of secondary metabolites. Results using knockout mutants clearly demonstrated that neither patulin nor citrinin are required by P. expansum to successfully infect apples. Li et al. ( MPMI-12-14-0398-FI ) reported similar results and conclusions in their recently accepted paper.

High-throughput biochemical fingerprinting of Saccharomyces cerevisiae by Fourier transform infrared spectroscopy

Single-channel optical density measurements of population growth are the dominant large scale phenotyping methodology for bridging the gene-function gap in yeast. However, a substantial amount of the genetic variation induced by single allele, single gene or double gene knock-out technologies fail to manifest in detectable growth phenotypes under conditions readily testable in the laboratory. Thus, new high-throughput phenotyping technologies capable of providing information about molecular level consequences of genetic variation are sorely needed. Here we report a protocol for high-throughput Fourier transform infrared spectroscopy (FTIR) measuring biochemical fingerprints of yeast strains. It includes high-throughput cultivation for FTIR spectroscopy, FTIR measurements and spectral pre-treatment to increase measurement accuracy. We demonstrate its capacity to distinguish not only yeast genera, species and populations, but also strains that differ only by a single gene, its excellent signal-to-noise ratio and its relative robustness to measurement bias. Finally, we illustrated its applicability by determining the FTIR signatures of all viable Saccharomyces cerevisiae single gene knock-outs corresponding to lipid biosynthesis genes. Many of the examined knock-out strains showed distinct, highly reproducible FTIR phenotypes despite having no detectable growth phenotype. These phenotypes were confirmed by conventional lipid analysis and could be linked to specific changes in lipid composition. We conclude that the introduced protocol is robust to noise and bias, possible to apply on a very large scale, and capable of generating biologically meaningful biochemical fingerprints that are strain specific, even when strains lack detectable growth phenotypes. Thus, it has a substantial potential for application in the molecular functionalization of the yeast genome.

Single site suppressors of a fission yeast temperature-sensitive mutant in cdc48 identified by whole genome sequencing

The protein called p97 in mammals and Cdc48 in budding and fission yeast is a homo-hexameric, ring-shaped, ubiquitin-dependent ATPase complex involved in a range of cellular functions, including protein degradation, vesicle fusion, DNA repair, and cell division. The cdc48⁺ gene is essential for viability in fission yeast, and point mutations in the human orthologue have been linked to disease. To analyze the function of p97/Cdc48 further, we performed a screen for cold-sensitive suppressors of the temperature-sensitive cdc48-353 fission yeast strain. In total, 29 independent pseudo revertants that had lost the temperature-sensitive growth defect of the cdc48-353 strain were isolated. Of these, 28 had instead acquired a cold-sensitive phenotype. Since the suppressors were all spontaneous mutants, and not the result of mutagenesis induced by chemicals or UV irradiation, we reasoned that the genome sequences of the 29 independent cdc48-353 suppressors were most likely identical with the exception of the acquired suppressor mutations. This prompted us to test if a whole genome sequencing approach would allow us to map the mutations. Indeed genome sequencing unambiguously revealed that the cold-sensitive suppressors were all second site intragenic cdc48 mutants. Projecting these onto the Cdc48 structure revealed that while the original temperature-sensitive G338D mutation is positioned near the central pore in the hexameric ring, the suppressor mutations locate to subunit-subunit and inter-domain boundaries. This suggests that Cdc48-353 is structurally compromized at the restrictive temperature, but re-established in the suppressor mutants. The last suppressor was an extragenic frame shift mutation in the ufd1 gene, which encodes a known Cdc48 co-factor. In conclusion, we show, using a novel whole genome sequencing approach, that Cdc48-353 is structurally compromized at the restrictive temperature, but stabilized in the suppressors.

Discovery of a novel small secreted protein family with conserved N-terminal IGY motif in Dikarya fungi

BACKGROUND

Small secreted proteins (SSPs) are employed by plant pathogenic fungi as essential strategic tools for their successful colonization. SSPs are often species-specific and so far only a few widely phylogenetically distributed SSPs have been identified.

RESULTS

A novel fungal SSP family consisting of 107 members was identified in the poplar tree fungal pathogen Marssonina brunnea, which accounts for over 17% of its secretome. We named these proteins IGY proteins (IGYPs) based on the conserved three amino acids at the N-terminus. In spite of overall low sequence similarity among IGYPs; they showed conserved N- and C-terminal motifs and a unified gene structure. By RT-PCR-seq, we analyzed the IGYP gene models and validated their expressions as active genes during infection. IGYP homologues were also found in 25 other Dikarya fungal species, all of which shared conserved motifs and the same gene structure. Furthermore, 18 IGYPs from 11 fungi also shared similar genomic contexts. Real-time RT-PCR showed that 8 MbIGYPs were highly expressed in the biotrophic stage. Interestingly, transient assay of 12 MbIGYPs showed that the MbIGYP13 protein induced cell death in resistant poplar clones.

CONCLUSIONS

In total, 154 IGYPs in 26 fungi of the Dikarya subkingdom were discovered. Gene structure and genomic context analyses indicated that IGYPs originated from a common ancestor. In M. brunnea, the expansion of highly divergent MbIGYPs possibly is associated with plant-pathogen arms race.

Formation of hydrogen sulfide from cysteine in Saccharomyces cerevisiae BY4742: genome wide screen reveals a central role of the vacuole

Discoveries on the toxic effects of cysteine accumulation and, particularly, recent findings on the many physiological roles of one of the products of cysteine catabolism, hydrogen sulfide (H₂S), are highlighting the importance of this amino acid and sulfur metabolism in a range of cellular activities. It is also highlighting how little we know about this critical part of cellular metabolism. In the work described here, a genome-wide screen using a deletion collection of Saccharomyces cerevisiae revealed a surprising set of genes associated with this process. In addition, the yeast vacuole, not previously associated with cysteine catabolism, emerged as an important compartment for cysteine degradation. Most prominent among the vacuole-related mutants were those involved in vacuole acidification; we identified each of the eight subunits of a vacuole acidification sub-complex (V₁ of the yeast V-ATPase) as essential for cysteine degradation. Other functions identified included translation, RNA processing, folate-derived one-carbon metabolism, and mitochondrial iron-sulfur homeostasis. This work identified for the first time cellular factors affecting the fundamental process of cysteine catabolism. Results obtained significantly contribute to the understanding of this process and may provide insight into the underlying cause of cysteine accumulation and H₂S generation in eukaryotes.

Evolutionary Perspectives on Human Fungal Pathogens

Owing to their small size and paucity of phenotypic characters, progress in the evolutionary biology of microbes in general, and human pathogenic fungi in particular, has been linked to a series of advances in DNA sequencing over the past quarter century. Phylogenetics was the first area to benefit, with the achievement of a basic understanding of fungal phylogeny. Population genetics was the next advance, finding cryptic species everywhere, and recombination in species previously thought to be asexual. Comparative genomics saw the next advance, in which variation in gene content and changes in gene family size were found to be important sources of variation. Fungal population genomics is showing that gene flow among closely related populations and species provides yet another source of adaptive, genetic variation. Now, two means to associate genetic variation with phenotypic variation, "reverse ecology" for adaptive phenotypes, and genome-wide association of any phenotype, are letting evolutionary biology make a profound contribution to molecular developmental biology of pathogenic fungi.

Identification and functional analysis of AG1-IA specific genes of Rhizoctonia solani

Rhizoctonia solani is an important necrotrophic fungal pathogen which causes disease on diverse plant species. It has been classified into 14 genetically distinct anastomosis groups (AGs), however, very little is known about their genomic diversity. AG1-IA causes sheath blight disease in rice and controlling this disease remains a challenge for sustainable rice cultivation. Recently the draft genome sequences of AG1-IA (rice isolate) and AG1-IB (lettuce isolate) had become publicly available. In this study, using comparative genomics, we report identification of 3,942 R. solani genes that are uniquely present in AG1-IA. Many of these genes encode important biological, molecular functions and exhibit dynamic expression during in-planta growth of the pathogen in rice. Based upon sequence similarity with genes that are required for plant and human/zoonotic diseases, we identified several putative virulence/pathogenicity determinants amongst AG1-IA specific genes. While studying the expression of 19 randomly selected genes, we identified three genes highly up-regulated during in-planta growth. The detailed in silico characterization of these genes and extent of their up-regulation in different rice genotypes, having variable degree of disease susceptibility, suggests their importance in rice-Rhizoctonia interactions. In summary, the present study reports identification, functional characterization of AG1-IA specific genes and predicts important virulence determinants that might enable the pathogen to grow inside hostile plant environment. Further characterization of these genes would shed useful insights about the pathogenicity mechanism of AG1-IA on rice.

Evolution of a morphological novelty occurred before genome compaction in a lineage of extreme parasites

Intracellular parasitism results in extreme adaptations, whose evolutionary history is difficult to understand, because the parasites and their known free-living relatives are so divergent from one another. Microsporidia are intracellular parasites of humans and other animals, which evolved highly specialized morphological structures, but also extreme physiologic and genomic simplification. They are suggested to be an early-diverging branch on the fungal tree, but comparisons to other species are difficult because their rates of molecular evolution are exceptionally high. Mitochondria in microsporidia have degenerated into organelles called mitosomes, which have lost a genome and the ability to produce ATP. Here we describe a gut parasite of the crustacean Daphnia that despite having remarkable morphological similarity to the microsporidia, has retained genomic features of its fungal ancestors. This parasite, which we name Mitosporidium daphniae gen. et sp. nov., possesses a mitochondrial genome including genes for oxidative phosphorylation, yet a spore stage with a highly specialized infection apparatus--the polar tube--uniquely known only from microsporidia. Phylogenomics places M. daphniae at the root of the microsporidia. A comparative genomic analysis suggests that the reduction in energy metabolism, a prominent feature of microsporidian evolution, was preceded by a reduction in the machinery controlling cell cycle, DNA recombination, repair, and gene expression. These data show that the morphological features unique to M. daphniae and other microsporidia were already present before the lineage evolved the extreme host metabolic dependence and loss of mitochondrial respiration for which microsporidia are well known.

Nothing special in the specialist? Draft genome sequence of Cryomyces antarcticus, the most extremophilic fungus from Antarctica

The draft genome of the Antarctic endemic fungus Cryomyces antarcticus is presented. This rock inhabiting, microcolonial fungus is extremely stress tolerant and it is a model organism for exobiology and studies on stress resistance in Eukaryots. Since this fungus is a specialist in the most extreme environment of the Earth, the analysis of its genome is of important value for the understanding of fungal genome evolution and stress adaptation. A comparison with Neurospora crassa as well as with other microcolonial fungi shows that the fungus has a genome size of 24 Mbp, which is the average in the fungal kingdom. Although sexual reproduction was never observed in this fungus, 34 mating genes are present with protein homologs in the classes Eurotiomycetes, Sordariomycetes and Dothideomycetes. The first analysis of the draft genome did not reveal any significant deviations of this genome from comparative species and mesophilic hyphomycetes.

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

BACKGROUND

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

RESULTS

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG's annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii's metabolism in comparison to the one of S. cerevisiae (post-WGD - whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

CONCLUSIONS

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii's metabolism.

Cytochrome P450 monooxygenase CYP53 family in fungi: comparative structural and evolutionary analysis and its role as a common alternative anti-fungal drug target

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins whose role as a drug target against pathogenic microbes has been explored because of their stereo- and regio-specific oxidation activity. We aimed to assess the CYP53 family's role as a common alternative drug target against animal (including human) and plant pathogenic fungi and its role in fungal-mediated wood degradation. Genome-wide analysis of fungal species revealed the presence of CYP53 members in ascomycetes and basidiomycetes. Basidiomycetes had a higher number of CYP53 members in their genomes than ascomycetes. Only two CYP53 subfamilies were found in ascomycetes and six subfamilies in basidiomycetes, suggesting that during the divergence of phyla ascomycetes lost CYP53 P450s. According to phylogenetic and gene-structure analysis, enrichment of CYP53 P450s in basidiomycetes occurred due to the extensive duplication of CYP53 P450s in their genomes. Numerous amino acids (103) were found to be conserved in the ascomycetes CYP53 P450s, against only seven in basidiomycetes CYP53 P450s. 3D-modelling and active-site cavity mapping data revealed that the ascomycetes CYP53 P450s have a highly conserved protein structure whereby 78% amino acids in the active-site cavity were found to be conserved. Because of this rigid nature of ascomycetes CYP53 P450s' active site cavity, any inhibitor directed against this P450 family can serve as a common anti-fungal drug target, particularly toward pathogenic ascomycetes. The dynamic nature of basidiomycetes CYP53 P450s at a gene and protein level indicates that these P450s are destined to acquire novel functions. Functional analysis of CYP53 P450s strongly supported our hypothesis that the ascomycetes CYP53 P450s ability is limited for detoxification of toxic molecules, whereas basidiomycetes CYP53 P450s play an additional role, i.e. involvement in degradation of wood and its derived components. This study is the first report on genome-wide comparative structural (gene and protein structure-level) and evolutionary analysis of a fungal P450 family.

De-novo learning of genome-scale regulatory networks in S. cerevisiae

De-novo reverse-engineering of genome-scale regulatory networks is a fundamental problem of biological and translational research. One of the major obstacles in developing and evaluating approaches for de-novo gene network reconstruction is the absence of high-quality genome-scale gold-standard networks of direct regulatory interactions. To establish a foundation for assessing the accuracy of de-novo gene network reverse-engineering, we constructed high-quality genome-scale gold-standard networks of direct regulatory interactions in Saccharomyces cerevisiae that incorporate binding and gene knockout data. Then we used 7 performance metrics to assess accuracy of 18 statistical association-based approaches for de-novo network reverse-engineering in 13 different datasets spanning over 4 data types. We found that most reconstructed networks had statistically significant accuracies. We also determined which statistical approaches and datasets/data types lead to networks with better reconstruction accuracies. While we found that de-novo reverse-engineering of the entire network is a challenging problem, it is possible to reconstruct sub-networks around some transcription factors with good accuracy. The latter transcription factors can be identified by assessing their connectivity in the inferred networks. Overall, this study provides the gene network reverse-engineering community with a rigorous assessment of the accuracy of S. cerevisiae gene network reconstruction and variability in performance of various approaches for learning both the entire network and sub-networks around transcription factors.

The genomic landscape of compensatory evolution

The Genomic Landscape of Compensatory Evolution Laboratory selection experiment explains how organisms compensate for the loss of genes during evolution, and reveals the deleterious side-effects of this process when adapting to novel environments.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.

While core cellular processes are generally conserved during evolution, the constituent genes differ somewhat between related species with similar lifestyles. Why should this be so? In this work, we propose that gene loss may initially be deleterious, but organisms can recover fitness by the accumulation of compensatory mutations elsewhere in the genome. To investigate this process in the laboratory, we investigated 180 haploid yeast strains, each of which initially displayed slow growth owing to the deletion of a single gene. Laboratory evolutionary experiments revealed that defects in a broad range of molecular processes can readily be compensated during evolution. Genomic analyses and functional assays demonstrated that compensatory evolution generates hidden genetic and physiological variation across parallel evolving lines, which can be revealed when the environment changes. Strikingly, despite nearly full recovery of fitness, the wild-type genomic expression pattern is generally not restored. Based on these results, we argue that genomes undergo major changes not simply to adapt to external conditions but also to compensate for previously accumulated deleterious mutations.

Population genomics of the fission yeast Schizosaccharomyces pombe

The fission yeast Schizosaccharomyces pombe has been widely used as a model eukaryote to study a diverse range of biological processes. However, population genetic studies of this species have been limited to date, and we know very little about the evolutionary processes and selective pressures that are shaping its genome. Here, we sequenced the genomes of 32 worldwide S. pombe strains and examined the pattern of polymorphisms across their genomes. In addition to introns and untranslated regions (UTRs), intergenic regions also exhibited lower levels of nucleotide diversity than synonymous sites, suggesting that a considerable amount of noncoding DNA is under selective constraint and thus likely to be functional. A number of genomic regions showed a reduction of nucleotide diversity probably caused by selective sweeps. We also identified a region close to the end of chromosome 3 where an extremely high level of divergence was observed between 5 of the 32 strains and the remain 27, possibly due to introgression, strong positive selection, or that region being responsible for reproductive isolation. Our study should serve as an important starting point in using a population genomics approach to further elucidate the biology of this important model organism.

Fungi on the skin: dermatophytes and Malassezia

Several human skin diseases and disorders are associated with two groups of fungi, the dermatophytes and Malassezia. Although these skin-related problems are not generally life threatening, they are among the most common diseases and disorders of mankind. These fungi are phylogenetically divergent, with the dermatophytes within the Ascomycota and Malassezia within Basidiomycota. Genome analysis indicates that the adaptations to the skin environment are different in these two groups of fungi. Malassezia are dependent on host lipids and secrete lipases and phospholipases that likely release host fatty acids. The dermatophytes encode multiple enzymes with potential roles in modulating host interactions: polyketide synthases, nonribosomal peptide synthetases, LysM, proteases, kinases, and pseudokinases. These two fungal groups have maximized their interactions with the host using two very different mechanisms.

The endo-arabinanase BcAra1 is a novel host-specific virulence factor of the necrotic fungal phytopathogen Botrytis cinerea

The plant cell wall is one of the first physical interfaces encountered by plant pathogens and consists of polysaccharides, of which arabinan is an important constituent. During infection, the necrotrophic plant pathogen Botrytis cinerea secretes a cocktail of plant cell-wall-degrading enzymes, including endo-arabinanase activity, which carries out the breakdown of arabinan. The roles of arabinan and endo-arabinanases during microbial infection were thus far elusive. In this study, the gene Bcara1 encoding for a novel α-1,5-L-endo-arabinanase was identified and the heterologously expressed BcAra1 protein was shown to hydrolyze linear arabinan with high efficiency whereas little or no activity was observed against the other oligo- and polysaccharides tested. The Bcara1 knockout mutants displayed reduced arabinanase activity in vitro and severe retardation in secondary lesion formation during infection of Arabidopsis leaves. These results indicate that BcAra1 is a novel endo-arabinanase and plays an important role during the infection of Arabidopsis. Interestingly, the level of Bcara1 transcript was considerably lower during the infection of Nicotiana benthamiana compared with Arabidopsis and, consequently, the ΔBcara1 mutants showed the wild-type level of virulence on N. benthamiana leaves. These results support the conclusion that the expression of Bcara1 is host dependent and is a key determinant of the disease outcome.

Genome sequencing of four Aureobasidium pullulans varieties: biotechnological potential, stress tolerance, and description of new species

BACKGROUND

Aureobasidium pullulans is a black-yeast-like fungus used for production of the polysaccharide pullulan and the antimycotic aureobasidin A, and as a biocontrol agent in agriculture. It can cause opportunistic human infections, and it inhabits various extreme environments. To promote the understanding of these traits, we performed de-novo genome sequencing of the four varieties of A. pullulans.

RESULTS

The 25.43-29.62 Mb genomes of these four varieties of A. pullulans encode between 10266 and 11866 predicted proteins. Their genomes encode most of the enzyme families involved in degradation of plant material and many sugar transporters, and they have genes possibly associated with degradation of plastic and aromatic compounds. Proteins believed to be involved in the synthesis of pullulan and siderophores, but not of aureobasidin A, are predicted. Putative stress-tolerance genes include several aquaporins and aquaglyceroporins, large numbers of alkali-metal cation transporters, genes for the synthesis of compatible solutes and melanin, all of the components of the high-osmolarity glycerol pathway, and bacteriorhodopsin-like proteins. All of these genomes contain a homothallic mating-type locus.

CONCLUSIONS

The differences between these four varieties of A. pullulans are large enough to justify their redefinition as separate species: A. pullulans, A. melanogenum, A. subglaciale and A. namibiae. The redundancy observed in several gene families can be linked to the nutritional versatility of these species and their particular stress tolerance. The availability of the genome sequences of the four Aureobasidium species should improve their biotechnological exploitation and promote our understanding of their stress-tolerance mechanisms, diverse lifestyles, and pathogenic potential.

Verticillium dahliae race 2-specific PCR reveals a high frequency of race 2 strains in commercial spinach seed lots and delineates race structure

Two pathogenic races of Verticillium dahliae have been described on lettuce and tomato. Host resistance to race 1 is governed by plant immune receptors that recognize the race 1-specific fungal effector Ave1. Only partial resistance to race 2 exists in lettuce. Although polymerase chain reaction (PCR) assays are available to identify race 1, no complementary test exists to positively identify race 2, except for lengthy pathogenicity assays on host differentials. Using the genome sequences of two isolates of V. dahliae, one each from races 1 and 2, we identified potential markers and PCR primers to distinguish the two races. Several primer pairs based on polymorphisms between the races were designed and tested on reference isolates of known race. One primer pair, VdR2F-VdR2R, consistently yielded a 256-bp amplicon in all race 2 isolates exclusively. We screened DNA from 677 V. dahliae isolates, including 340 from spinach seedlots, with the above primer pair and a previously published race 1-specific primer pair. DNA from isolates that did not amplify with race 1-specific PCRs amplified with the race 2-specific primers. To validate this, two differential lines of lettuce were inoculated with 53 arbitrarily selected isolates from spinach seed and their pathogenicity and virulence were assessed in a greenhouse. The reactions of the differential cultivars strongly supported the PCR data. V. dahliae race structure was investigated in crops in coastal California and elsewhere using primers specific to the two races. All artichoke isolates from California were race 1, whereas nearly all tomato isolates were race 2. Isolates from lettuce, pepper, and strawberry from California as well as isolates from spinach seed from two of four countries comprised both races, whereas only race 2 was observed in cotton, mint, olive, and potato. This highlights the importance of identifying resistance against race 2 in different hosts. The technique developed in this study will benefit studies in ecology, population biology, disease surveillance, and epidemiology at local and global scales, and resistance breeding against race 2 in lettuce and other crops.

Rif1 regulates initiation timing of late replication origins throughout the S. cerevisiae genome

Chromosomal DNA replication involves the coordinated activity of hundreds to thousands of replication origins. Individual replication origins are subject to epigenetic regulation of their activity during S-phase, resulting in differential efficiencies and timings of replication initiation during S-phase. This regulation is thought to involve chromatin structure and organization into timing domains with differential ability to recruit limiting replication factors. Rif1 has recently been identified as a genome-wide regulator of replication timing in fission yeast and in mammalian cells. However, previous studies in budding yeast have suggested that Rif1’s role in controlling replication timing may be limited to subtelomeric domains and derives from its established role in telomere length regulation. We have analyzed replication timing by analyzing BrdU incorporation genome-wide, and report that Rif1 regulates the timing of late/dormant replication origins throughout the S. cerevisiae genome. Analysis of pfa4Δ cells, which are defective in palmitoylation and membrane association of Rif1, suggests that replication timing regulation by Rif1 is independent of its role in localizing telomeres to the nuclear periphery. Intra-S checkpoint signaling is intact in rif1Δ cells, and checkpoint-defective mec1Δ cells do not comparably deregulate replication timing, together indicating that Rif1 regulates replication timing through a mechanism independent of this checkpoint. Our results indicate that the Rif1 mechanism regulates origin timing irrespective of proximity to a chromosome end, and suggest instead that telomere sequences merely provide abundant binding sites for proteins that recruit Rif1. Still, the abundance of Rif1 binding in telomeric domains may facilitate Rif1-mediated repression of non-telomeric origins that are more distal from centromeres.