A new inducible CRISPR-Cas9 system useful for genome editing and study of double-strand break repair in Candida glabrata

In recent years, the CRISPR-Cas9 system has proven extremely useful for genome editing in many species, including the model yeast Saccharomyces cerevisiae and other yeast species such as Candida glabrata. Inducible CRISPR-Cas9 systems have the additional advantage of allowing to separate the transformation step of the organism by the CRISPR-Cas9 system, from the cutting and repair steps. This has indeed been developed in S. cerevisiae, where most inducible expression systems rely on the GAL promoters. Unfortunately, C. glabrata is gal- and lacks the GAL genes, like many other yeast species. We report here the use of a vector expressing cas9 under the control of the MET3 promoter, with the guide RNA cloned into the same plasmid. We show that it can be used efficiently in C. glabrata, for both described outcomes of CRISPR-Cas9-induced chromosome breaks; nonhomologous end joining in the absence of a homologous repair template; and homologous recombination in the presence of such a template. This system therefore allows easy editing of the genome of C. glabrata, and its inducibility may allow identification of essential genes in this asexual yeast, where spore lethality cannot be observed, as well as the study of double-strand break repair.

The High-Quality Complete Genome Sequence of the Opportunistic Fungal Pathogen Candida vulturna CBS 14366T

Candida vulturna is a new member of the Candida haemulonii species complex that recently received much attention as it includes the emerging multidrug-resistant pathogen Candida auris. Here, we describe the high-quality genome sequence of C. vulturna type strain CBS 14366^T to cover all genomes of pathogenic C. haemulonii species complex members.

Whole Genome Sequence Resource of the Asian Pear Scab Pathogen Venturia nashicola

Venturia nashicola, the cause of scab disease of Asian pears, is a host-specific, biotrophic fungus. It is restricted to Asia and is regarded as a quarantine threat outside this region. European pear displays nonhost resistance (NHR) to V. nashicola and Asian pears are nonhosts of V. pyrina (the cause of European pear scab disease). The host specificity of these two fungi is likely governed by differences in their effector arsenals, with a subset hypothesized to activate NHR. The Pyrus-Venturia pathosystem provides an opportunity to dissect the underlying genetics of nonhost interactions in this potentially more durable form of resistance. The V. nashicola genome will enable comparisons to other Venturia spp. genomes to identify effectors that potentially activate NHR in the pear scab pathosystem.

Genome Sequence Resource for Ilyonectria mors-panacis, Causing Rusty Root Rot of Panax notoginseng

Ilyonectria mors-panacis is the cause of a serious disease hampering the production of Panax notoginseng, an important Chinese medicinal herb, widely used for its anti-inflammatory, antifatigue, hepato-protective, and coronary heart disease prevention effects. Here, we report the first Illumina-Pacbio hybrid sequenced draft genome assembly of I. mors-panacis strain G3B and its annotation. The availability of this genome sequence not only represents an important tool toward understanding the genetics behind the infection mechanism of I. mors-panacis strain G3B but also will help illuminate the complexities of the taxonomy of this species.

Genome Data of Fusarium oxysporum f. sp. cubense Race 1 and Tropical Race 4 Isolates Using Long-Read Sequencing

Fusarium wilt of banana is caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. cubense. We generated two chromosome-level assemblies of F. oxysporum f. sp. cubense race 1 and tropical race 4 strains using single-molecule real-time sequencing. The F. oxysporum f. sp. cubense race 1 and tropical race 4 assemblies had 35 and 29 contigs with contig N50 lengths of 2.08 and 4.28 Mb, respectively. These two new references genomes represent a greater than 100-fold improvement over the contig N50 statistics of the previous short-read-based F. oxysporum f. sp. cubense assemblies. The two high-quality assemblies reported here will be a valuable resource for the comparative analysis of F. oxysporum f. sp. cubense races at the pathogenic level.

Genomic features and evolution of the conditionally dispensable chromosome in the tangerine pathotype of Alternaria alternata

The tangerine pathotype of the ascomycete fungus Alternaria alternata is the causal agent of citrus brown spot, which can result in significant losses of both yield and marketability for tangerines worldwide. A conditionally dispensable chromosome (CDC), which harbours the host-selective ACT toxin gene cluster, is required for tangerine pathogenicity of A. alternata. To understand the genetic makeup and evolution of the tangerine pathotype CDC, we isolated and sequenced the CDCs of the A. alternata Z7 strain and analysed the function and evolution of their genes. The A. alternata Z7 strain has two CDCs (~1.1 and ~0.8 Mb, respectively), and the longer Z7 CDC contains all but one contig of the shorter one. Z7 CDCs contain 254 predicted protein-coding genes, which are enriched in functional categories associated with 'metabolic process' (55 genes, P = 0.037). Relatively few of the CDC genes can be classified as carbohydrate-active enzymes (CAZymes) (4) and transporters (19) and none as kinases. Evolutionary analysis of the 254 CDC proteins showed that their evolutionary conservation tends to be restricted within the genus Alternaria and that the CDC genes evolve faster than genes in the essential chromosomes, likely due to fewer selective constraints. Interestingly, phylogenetic analysis suggested that four of the 25 genes responsible for the ACT toxin production were likely transferred from Colletotrichum (Sordariomycetes). Functional experiments showed that two of them are essential for the virulence of the tangerine pathotype of A. alternata. These results provide new insights into the function and evolution of CDC genes in Alternaria.

Whole Genome Sequences of the Tea Leaf Spot Pathogen Didymella segeticola

The fungal pathogen Didymella segeticola (basionym Phoma segeticola) causes leaf spot on tea (Camellia sinensis), which leads to a loss in tea leaf production in Guizhou Province, China. D. segeticola isolate GZSQ-4 was sequenced using Illumina HiSeq and Pacific Biosciences (PacBio) RS technologies, and then assembled to approximately 33.4 Mbp with a scaffold N50 value of approximately 2.3 Mbp. In total, 10,893 genes were predicted using the Nonredundant, Gene Ontology, Clusters of Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes, and SWISS-PROT databases. The whole-genome sequence of D. segeticola will provide a resource for future research on host-pathogen interactions, determination of trait-specific genes, pathogen evolution, and plant-host adaptation mechanisms.

Liberibacter crescens Is a Cultured Surrogate for Functional Genomics of Uncultured Pathogenic 'Candidatus Liberibacter' spp. and Is Naturally Competent for Transformation

'Candidatus Liberibacter' spp. are uncultured insect endosymbionts and phloem-limited bacterial plant pathogens associated with diseases ranging from severe to nearly asymptomatic. 'Ca. L. asiaticus', causal agent of Huanglongbing or citrus "greening," and 'Ca. L. solanacearum', causal agent of potato zebra chip disease, respectively threaten citrus and potato production worldwide. Research on both pathogens has been stymied by the inability to culture these agents and to reinoculate into any host. Only a single isolate of a single species of Liberibacter, Liberibacter crescens, has been axenically cultured. L. crescens strain BT-1 is genetically tractable to standard molecular manipulation techniques and has been developed as a surrogate model for functional studies of genes, regulatory elements, promoters, and secreted effectors derived from the uncultured pathogenic Liberibacters. Detailed, step-by-step, and highly reproducible protocols for axenic culture, transformation, and targeted gene knockouts of L. crescens are described. In the course of developing these protocols, we found that L. crescens is also naturally competent for direct uptake and homology-guided chromosomal integration of both linear and circular plasmid DNA. The efficiency of natural transformation was about an order of magnitude higher using circular plasmid DNA compared with linearized fragments. Natural transformation using a replicative plasmid was obtained at a rate of approximately 900 transformants per microgram of plasmid, whereas electroporation using the same plasmid resulted in 6 × 10⁴ transformants. Homology-guided marker interruptions using either natural uptake or electroporation of nonreplicative plasmids yielded 10 to 12 transformation events per microgram of DNA, whereas similar interruptions using linear fragments via natural uptake yielded up to 34 transformation events per microgram of DNA.

Whole-genome resequencing of Trichophyton rubrum provides insights into population differentiation and drug resistance

Trichophyton rubrum (T. rubrum) is anthropophilic fungus and thus a very common cause of dermatophyte infections around the world. Infection of T. rubrum could result in conditions such as tinea capitis, tinea corporis, tinea inguinalis, tinea manus, tinea unguium, or tinea pedis. Because of this, the resistance of T. rubrum to antifungal therapies has drawn extensive research interest. However, the pathogenic characteristics of T. rubrum, such as site of infections, geographic location and host groups, have yet to be explored. In this study, the whole genome of 48 strains from different regions is resequenced and the population structure and association of single nucleotide polymorphism with resistance to six widely used antifungal drugs are analyzed. A total of 23,394 genomic variations are detected, which cover 2165 genes with only 15.14% of the variations located in exons. The population structure of T. rubrum is monomorphic, and genetic diversity is very low. Population structure analysis shows that the 48 sampled strains can be divided into two sub-populations. The gene TERG_08771 harboring the highest SNPs density is found to be associated with resistance to voriconazole. Although many proteins have yet to be identified and explored, association studies could still be useful to identify drug resistance or drug-susceptible loci, which would warrant further insightful investigations.

Complex genetic and epigenetic regulation deviates gene expression from a unifying global transcriptional program

Environmental or genetic perturbations lead to gene expression changes. While most analyses of these changes emphasize the presence of qualitative differences on just a few genes, we now know that changes are widespread. This large-scale variation has been linked to the exclusive influence of a global transcriptional program determined by the new physiological state of the cell. However, given the sophistication of eukaryotic regulation, we expect to have a complex architecture of specific control affecting this program. Here, we examine this architecture. Using data of Saccharomyces cerevisiae expression in different nutrient conditions, we first propose a five-sector genome partition, which integrates earlier models of resource allocation, as a framework to examine the deviations from the global control. In this scheme, we recognize invariant genes, whose regulation is dominated by physiology, specific genes, which substantially depart from it, and two additional classes that contain the frequently assumed growth-dependent genes. Whereas the invariant class shows a considerable absence of specific regulation, the rest is enriched by regulation at the level of transcription factors (TFs) and epigenetic modulators. We nevertheless find markedly different strategies in how these classes deviate. On the one hand, there are TFs that act in a unique way between partition constituents, and on the other, the action of chromatin modifiers is significantly diverse. The balance between regulatory strategies ultimately modulates the action of the general transcription machinery and therefore limits the possibility of establishing a unifying program of expression change at a genomic scale.

How can we understand expression changes observed as a result of environmental or genetic perturbations? This issue has been conventionally answered by examining small groups of genes whose expression becomes qualitatively altered after these perturbations. But this approach is too simplistic, as we now know that extensive variation is typically observed. To explain this variation, recent works proposed a model in which genome-wide changes are explained by the action of a general program of transcription. Our manuscript reasons that given the complexity of eukaryotic transcriptional control, a unifying program of regulation cannot be achievable. Instead, we propose within an integrated framework of resource allocation that a rich structure of deviations from it exists and that by characterizing these deviations we can fully appreciate large-scale expression change.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

The emerging blast disease on wheat is proving even harder to control than the ancient, still-problematic rice blast disease. Potential wheat resistance identified using strains isolated soon after disease emergence are no longer effective in controlling recent aggressive field isolates from wheat in South America and South Asia. We construct a high-quality assembly of an aggressive, recently-isolated wheat blast fungal strain and the first assembled mini-chromosome genome sequence of wheat and rice blast pathogens. We report that recent wheat pathogens can contain one or two highly-variable dispensable mini-chromosomes, each with an amalgamation of fungal effector genes and other sequences that are duplicated or absent from indispensable core chromosome ends. Well-studied effectors found on different core chromosomes in rice pathogens appear side-by-side in wheat pathogen mini-chromosomes. The rice pathogen often overcomes deployed resistance genes by deleting triggering effector genes. We propose that the fast-evolving effector-rich compartment of the wheat blast fungus is a combination of core chromosome ends and mobile mini-chromosomes that are easily lost from individual strains. Localization of effectors on mini-chromosomes would therefore accelerate pathogen adaptation in the field.

High Contiguity Whole Genome Sequence and Gene Annotation Resource for Two Venturia nashicola Isolates

Venturia nashicola is a fungal pathogen that causes Asian pear scab disease. This pathogen is of particular importance in Northeast Asian countries, where Asian pears are grown industrially. Scab disease in Asian pear is currently controlled by fungicide spraying and this situation calls for developing scab resistant cultivars. High-quality genome data are therefore required for in-depth comparative genome analysis of different isolates of V. nashicola and V. pyrina, a closely related species, which only infects European pear plants. Here, we report the high-contiguity whole genome assembly of two V. nashicola isolates, which is expected to enable genome comparisons for identification of the genes involved in host range determination of V. nashicola.

CRISPR/Cas9 genome editing technology in filamentous fungi: progress and perspective

Filamentous fungi play an important role in human health and industrial/agricultural production. With the increasing number of full genomes available for fungal species, the study of filamentous fungi has brought about a wider range of genetic manipulation opportunities. However, the utilization of traditional methods to study fungi is time consuming and laborious. Recent rapid progress and wide application of a versatile genome editing technology, i.e., the CRISPR (clustered regularly interspaced short palindromic repeat)-Cas9 (CRISPR-related nuclease 9) system, has revolutionized biological research and has many innovative applications in a wide range of fields showing great promise in research and application of filamentous fungi. In this review, we introduce the CRISPR/Cas9 genome editing technology focusing on its application in research of filamentous fungi and we discuss the general considerations of genome editing using CRISPR/Cas9 system illustrating vector construction, multiple editing strategies, technical consideration of different sizes of homology arms on genome editing efficiency, off-target effects, and different transformation methodologies. In addition, we discuss the challenges encountered using CRISPR/Cas9 technology and give the perspectives of future applications of CRISPR/Cas9 technology for basic research and practical application of filamentous fungi.

Genome architecture and stability in the Saccharomyces cerevisiae knockout collection

Despite major progress in defining the functional roles of genes, a complete understanding of their influences is far from being realized, even in relatively simple organisms. A major milestone in this direction arose via the completion of the yeast Saccharomyces cerevisiae gene-knockout collection (YKOC), which has enabled high-throughput reverse genetics, phenotypic screenings and analyses of synthetic-genetic interactions1-3. Ensuing experimental work has also highlighted some inconsistencies and mistakes in the YKOC, or genome instability events that rebalance the effects of specific knockouts4-6, but a complete overview of these is lacking. The identification and analysis of genes that are required for maintaining genomic stability have traditionally relied on reporter assays and on the study of deletions of individual genes, but whole-genome-sequencing technologies now enable-in principle-the direct observation of genome instability globally and at scale. To exploit this opportunity, we sequenced the whole genomes of nearly all of the 4,732 strains comprising the homozygous diploid YKOC. Here, by extracting information on copy-number variation of tandem and interspersed repetitive DNA elements, we describe-for almost every single non-essential gene-the genomic alterations that are induced by its loss. Analysis of this dataset reveals genes that affect the maintenance of various genomic elements, highlights cross-talks between nuclear and mitochondrial genome stability, and shows how strains have genetically adapted to life in the absence of individual non-essential genes.

Temperature preference can bias parental genome retention during hybrid evolution

Interspecific hybridization can introduce genetic variation that aids in adaptation to new or changing environments. Here, we investigate how hybrid adaptation to temperature and nutrient limitation may alter parental genome representation over time. We evolved Saccharomyces cerevisiae x Saccharomyces uvarum hybrids in nutrient-limited continuous culture at 15°C for 200 generations. In comparison to previous evolution experiments at 30°C, we identified a number of responses only observed in the colder temperature regime, including the loss of the S. cerevisiae allele in favor of the cryotolerant S. uvarum allele for several portions of the hybrid genome. In particular, we discovered a genotype by environment interaction in the form of a loss of heterozygosity event on chromosome XIII; which species' haplotype is lost or maintained is dependent on the parental species' temperature preference and the temperature at which the hybrid was evolved. We show that a large contribution to this directionality is due to a temperature dependent fitness benefit at a single locus, the high affinity phosphate transporter gene PHO84. This work helps shape our understanding of what forces impact genome evolution after hybridization, and how environmental conditions may promote or disfavor the persistence of hybrids over time.

Effector gene reshuffling involves dispensable mini-chromosomes in the wheat blast fungus

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome.

Reshuffling yeast chromosomes with CRISPR/Cas9

Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions.

Chromosomes are highly dynamic objects that often undergo large structural variations such as reciprocal translocations. Such rearrangements can have dramatic functional consequences, as they can disrupt genes, change their regulation or create novel fusion genes at their breakpoints. For instance, 90–95% of patients diagnosed with chronic myeloid leukemia carry the Philadelphia chromosome characterized by a reciprocal translocation between chromosomes 9 and 22. In addition, translocations reorganize the genetic information along chromosomes, which in turn can modify the 3D architecture of the genome and potentially affect its functioning. Quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. Here, we report a novel CRISPR/Cas9-based technology allowing to generate with high efficiency and at a base-pair precision either uniquely targeted or multiple reciprocal translocations in yeast, without leaving any marker or scar in the genome. Engineering targeted reciprocal translocations allowed us for the first time to untangle the phenotypic impacts of large chromosomal rearrangements from that of point mutations. In addition, the generation of multiple translocations led to a large reorganization of the genetic information along the chromosomes, often including unanticipated large segmental duplications. We showed that reshuffling the genome resulted in the emergence of fitness advantage in stressful environmental conditions, even in strains where no gene was disrupted or amplified by the translocations.

Processing of eukaryotic Okazaki fragments by redundant nucleases can be uncoupled from ongoing DNA replication in vivo

Prior to ligation, each Okazaki fragment synthesized on the lagging strand in eukaryotes must be nucleolytically processed. Nuclease cleavage takes place in the context of 5′ flap structures generated via strand-displacement synthesis by DNA polymerase delta. At least three DNA nucleases: Rad27 (Fen1), Dna2 and Exo1, have been implicated in processing Okazaki fragment flaps. However, neither the contributions of individual nucleases to lagging-strand synthesis nor the structure of the DNA intermediates formed in their absence have been fully defined in vivo. By conditionally depleting lagging-strand nucleases and directly analyzing Okazaki fragments synthesized in vivo in Saccharomyces cerevisiae, we conduct a systematic evaluation of the impact of Rad27, Dna2 and Exo1 on lagging-strand synthesis. We find that Rad27 processes the majority of lagging-strand flaps, with a significant additional contribution from Exo1 but not from Dna2. When nuclease cleavage is impaired, we observe a reduction in strand-displacement synthesis as opposed to the widespread generation of long Okazaki fragment 5′ flaps, as predicted by some models. Further, using cell cycle-restricted constructs, we demonstrate that both the nucleolytic processing and the ligation of Okazaki fragments can be uncoupled from DNA replication and delayed until after synthesis of the majority of the genome is complete.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50–75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

Genome Sequence Resource of the Wide-Host-Range Anthracnose Pathogen Colletotrichum siamense

Colletotrichum siamense causes fruit or foliar disease called anthracnose on a variety of plant hosts such as vegetables, fruits, ornamental plants, and others, including chili pepper, apple, American cranberry, mango, orange, papaya, guava, rubber plant, jasmine, coffee berry, and tea plants. Here, we report the first Illumina-sequenced draft genome assembly of C. siamense strain ICMP 18578 and its annotation. This genome sequence provides a unique resource that will be useful for future research on the evolution of Colletotrichum spp. and improvement of anthracnose management strategies.

Enzymes of early-diverging, zoosporic fungi

The secretome, the complement of extracellular proteins, is a reflection of the interaction of an organism with its host or substrate, thus a determining factor for the organism’s fitness and competitiveness. Hence, the secretome impacts speciation and organismal evolution. The zoosporic Chytridiomycota, Blastocladiomycota, Neocallimastigomycota, and Cryptomycota represent the earliest diverging lineages of the Fungal Kingdom. The review describes the enzyme compositions of these zoosporic fungi, underscoring the enzymes involved in biomass degradation. The review connects the lifestyle and substrate affinities of the zoosporic fungi to the secretome composition by examining both classical phenotypic investigations and molecular/genomic-based studies. The carbohydrate-active enzyme profiles of 19 genome-sequenced species are summarized. Emphasis is given to recent advances in understanding the functional role of rumen fungi, the basis for the devastating chytridiomycosis, and the structure of fungal cellulosome. The approach taken by the review enables comparison of the secretome enzyme composition of anaerobic versus aerobic early-diverging fungi and comparison of enzyme portfolio of specialized parasites, pathogens, and saprotrophs. Early-diverging fungi digest most major types of biopolymers: cellulose, hemicellulose, pectin, chitin, and keratin. It is thus to be expected that early-diverging fungi in its entirety represents a rich and diverse pool of secreted, metabolic enzymes. The review presents the methods used for enzyme discovery, the diversity of enzymes found, the status and outlook for recombinant production, and the potential for applications. Comparative studies on the composition of secretome enzymes of early-diverging fungi would contribute to unraveling the basal lineages of fungi.

antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline

Secondary metabolites produced by bacteria and fungi are an important source of antimicrobials and other bioactive compounds. In recent years, genome mining has seen broad applications in identifying and characterizing new compounds as well as in metabolic engineering. Since 2011, the 'antibiotics and secondary metabolite analysis shell-antiSMASH' (https://antismash.secondarymetabolites.org) has assisted researchers in this, both as a web server and a standalone tool. It has established itself as the most widely used tool for identifying and analysing biosynthetic gene clusters (BGCs) in bacterial and fungal genome sequences. Here, we present an entirely redesigned and extended version 5 of antiSMASH. antiSMASH 5 adds detection rules for clusters encoding the biosynthesis of acyl-amino acids, β-lactones, fungal RiPPs, RaS-RiPPs, polybrominated diphenyl ethers, C-nucleosides, PPY-like ketones and lipolanthines. For type II polyketide synthase-encoding gene clusters, antiSMASH 5 now offers more detailed predictions. The HTML output visualization has been redesigned to improve the navigation and visual representation of annotations. We have again improved the runtime of analysis steps, making it possible to deliver comprehensive annotations for bacterial genomes within a few minutes. A new output file in the standard JavaScript object notation (JSON) format is aimed at downstream tools that process antiSMASH results programmatically.

ParaDB: A manually curated database containing genomic annotation for the human pathogenic fungi Paracoccidioides spp

BACKGROUND

The genus Paracoccidioides consists of thermodymorphic fungi responsible for Paracoccidioidomycosis (PCM), a systemic mycosis that has been registered to affect ~10 million people in Latin America. Biogeographical data subdivided the genus Paracoccidioides in five divergent subgroups, which have been recently classified as different species. Genomic sequencing of five Paracoccidioides isolates, representing each of these subgroups/species provided an important framework for the development of post-genomic studies with these fungi. However, functional annotations of these genomes have not been submitted to manual curation and, as a result, ~60-90% of the Paracoccidioides protein-coding genes (depending on isolate/annotation) are currently described as responsible for hypothetical proteins, without any further functional/structural description.

PRINCIPAL FINDINGS

The present work reviews the functional assignment of Paracoccidioides genes, reducing the number of hypothetical proteins to ~25-28%. These results were compiled in a relational database called ParaDB, dedicated to the main representatives of Paracoccidioides spp. ParaDB can be accessed through a friendly graphical interface, which offers search tools based on keywords or protein/DNA sequences. All data contained in ParaDB can be partially or completely downloaded through spreadsheet, multi-fasta and GFF3-formatted files, which can be subsequently used in a variety of downstream functional analyses. Moreover, the entire ParaDB environment has been configured in a Docker service, which has been submitted to the GitHub repository, ensuring long-term data availability to researchers. This service can be downloaded and used to perform fully functional local installations of the database in alternative computing ecosystems, allowing users to conduct their data mining and analyses in a personal and stable working environment.

CONCLUSIONS

These new annotations greatly reduce the number of genes identified solely as hypothetical proteins and are integrated into a dedicated database, providing resources to assist researchers in this field to conduct post-genomic studies with this group of human pathogenic fungi.

Genome-wide sequencing and metabolic annotation of Pythium irregulare CBS 494.86: understanding Eicosapentaenoic acid production

BACKGROUND

Pythium irregulare is an oleaginous Oomycete able to accumulate large amounts of lipids, including Eicosapentaenoic acid (EPA). EPA is an important and expensive dietary supplement with a promising and very competitive market, which is dependent on fish-oil extraction. This has prompted several research groups to study biotechnological routes to obtain specific fatty acids rather than a mixture of various lipids. Moreover, microorganisms can use low cost carbon sources for lipid production, thus reducing production costs. Previous studies have highlighted the production of EPA by P. irregulare, exploiting diverse low cost carbon sources that are produced in large amounts, such as vinasse, glycerol, and food wastewater. However, there is still a lack of knowledge about its biosynthetic pathways, because no functional annotation of any Pythium sp. exists yet. The goal of this work was to identify key genes and pathways related to EPA biosynthesis, in P. irregulare CBS 494.86, by sequencing and performing an unprecedented annotation of its genome, considering the possibility of using wastewater as a carbon source.

RESULTS

Genome sequencing provided 17,727 candidate genes, with 3809 of them associated with enzyme code and 945 with membrane transporter proteins. The functional annotation was compared with curated information of oleaginous organisms, understanding amino acids and fatty acids production, and consumption of carbon and nitrogen sources, present in the wastewater. The main features include the presence of genes related to the consumption of several sugars and candidate genes of unsaturated fatty acids production.

CONCLUSIONS

The whole metabolic genome presented, which is an unprecedented reconstruction of P. irregulare CBS 494.86, shows its potential to produce value-added products, in special EPA, for food and pharmaceutical industries, moreover it infers metabolic capabilities of the microorganism by incorporating information obtained from literature and genomic data, supplying information of great importance to future work.

A Flp-SUMO hybrid recombinase reveals multi-layered copy number control of a selfish DNA element through post-translational modification

Mechanisms for highly efficient chromosome-associated equal segregation, and for maintenance of steady state copy number, are at the heart of the evolutionary success of the 2-micron plasmid as a stable multi-copy extra-chromosomal selfish DNA element present in the yeast nucleus. The Flp site-specific recombination system housed by the plasmid, which is central to plasmid copy number maintenance, is regulated at multiple levels. Transcription of the FLP gene is fine-tuned by the repressor function of the plasmid-coded partitioning proteins Rep1 and Rep2 and their antagonist Raf1, which is also plasmid-coded. In addition, the Flp protein is regulated by the host’s post-translational modification machinery. Utilizing a Flp-SUMO fusion protein, which functionally mimics naturally sumoylated Flp, we demonstrate that the modification signals ubiquitination of Flp, followed by its proteasome-mediated degradation. Furthermore, reduced binding affinity and cooperativity of the modified Flp decrease its association with the plasmid FRT (Flp recombination target) sites, and/or increase its dissociation from them. The resulting attenuation of strand cleavage and recombination events safeguards against runaway increase in plasmid copy number, which is deleterious to the host—and indirectly—to the plasmid. These results have broader relevance to potential mechanisms by which selfish genomes minimize fitness conflicts with host genomes by holding in check the extra genetic load they pose.

Plasmids of budding yeasts, exemplified by the 2-micron plasmid of Saccharomyces cerevisiae, and mammalian papilloma and gammaherpes viruses typify eukaryotic extra-chromosomal selfish DNA elements. The plasmid and the viral episomes, despite the long evolutionary divergence of their hosts, share striking similarities in lifestyles. These include the ability to segregate to daughter cells by hitchhiking on chromosomes and to switch from cell cycle regulated replication to iterative replication for copy number maintenance. While selfish elements, including those integrated into chromosomes, rely on their hosts’ genetic potential for long-term survival, their genetic load is carefully regulated to minimize fitness conflicts with the hosts. Our study focuses on the Flp site-specific recombinase, which is central to the copy number control of the 2-micron plasmid and whose steady state levels are regulated through transcriptional control by plasmid coded proteins and through post-translational modification by the host’s sumoylation machinery. We demonstrate that sumoylation, in addition, attenuates the catalytic activity of Flp by diminishing its DNA binding affinity and inter-monomer cooperativity, providing another layer of protection against runaway increase in plasmid copy number. Population control by self-imposed and host-mediated mechanisms is likely a general strategy among selfish elements to ensure nearly conflict-free coexistence with host genomes.

Seasons of change: Mechanisms of genome evolution in human fungal pathogens

Fungi are a diverse kingdom of organisms capable of thriving in various niches across the world including those in close association with multicellular eukaryotes. Fungal pathogens that contribute to human disease reside both within the host as commensal organisms of the microbiota and the environment. Their niche of origin dictates how infection initiates but also places specific selective pressures on the fungal pathogen that contributes to its genome organization and genetic repertoire. Recent efforts to catalogue genomic variation among major human fungal pathogens have unveiled evolutionary themes that shape the fungal genome. Mechanisms ranging from large scale changes such as aneuploidy and ploidy cycling as well as more targeted mutations like base substitutions and gene copy number variations contribute to the evolution of these species, which are often under multiple competing selective pressures with their host, environment, and other microbes. Here, we provide an overview of the major selective pressures and mechanisms acting to evolve the genome of clinically important fungal pathogens of humans.

Coccinonectria pachysandricola, Causal Agent of a New Foliar Blight Disease of Sarcococca hookeriana

Woody plants of the Buxaceae, including species of Buxus, Pachysandra, and Sarcococca, are widely grown evergreen shrubs and groundcovers. Severe leaf spot symptoms were observed on S. hookeriana at the U.S. National Arboretum in Washington, DC, in 2016. Affected plants were growing adjacent to P. terminalis exhibiting Volutella blight symptoms. Fungi isolated from both hosts were identical based on morphology and multilocus phylogenetic analysis and were identified as Coccinonectria pachysandricola (Nectriaceae, Hypocreales), causal agent of Volutella blight of Pachysandra species. Pathogenicity tests established that Co. pachysandricola isolated from both hosts caused disease symptoms on P. terminalis and S. hookeriana, but not on B. sempervirens. Artificial inoculations with Pseudonectria foliicola, causal agent of Volutella blight of B. sempervirens, did not result in disease on P. terminalis or S. hookeriana. Wounding enhanced infection by Co. pachysandricola and Ps. foliicola on all hosts tested but was not required for disease development. Genome assemblies were generated for the Buxaceae pathogens that cause Volutella diseases: Co. pachysandricola, Ps. buxi, and Ps. foliicola; these ranged in size from 25.7 to 28.5 Mb. To our knowledge, this foliar blight of S. hookeriana represents a new disease for this host and is capable of causing considerable damage to infected plants.

Timing the Origin of Cryptococcus gattii sensu stricto, Southeastern United States

We conducted molecular clock analysis of whole-genome sequences from a set of autochthonous isolates of Cryptococcus gattii sensu stricto from the southeastern United States. Our analysis indicates that C. gattii arrived in the southeastern United States approximately 9,000–19,000 years ago, long before its arrival in the Pacific Northwest.

qDSB-Seq is a general method for genome-wide quantification of DNA double-strand breaks using sequencing

DNA double-strand breaks (DSBs) are among the most lethal types of DNA damage and frequently cause genome instability. Sequencing-based methods for mapping DSBs have been developed but they allow measurement only of relative frequencies of DSBs between loci, which limits our understanding of the physiological relevance of detected DSBs. Here we propose quantitative DSB sequencing (qDSB-Seq), a method providing both DSB frequencies per cell and their precise genomic coordinates. We induce spike-in DSBs by a site-specific endonuclease and use them to quantify detected DSBs (labeled, e.g., using i-BLESS). Utilizing qDSB-Seq, we determine numbers of DSBs induced by a radiomimetic drug and replication stress, and reveal two orders of magnitude differences in DSB frequencies. We also measure absolute frequencies of Top1-dependent DSBs at natural replication fork barriers. qDSB-Seq is compatible with various DSB labeling methods in different organisms and allows accurate comparisons of absolute DSB frequencies across samples.

Three phylogenetic groups have driven the recent population expansion of Cryptococcus neoformans

Cryptococcus neoformans (C. neoformans var. grubii) is an environmentally acquired pathogen causing 181,000 HIV-associated deaths each year. We sequenced 699 isolates, primarily C. neoformans from HIV-infected patients, from 5 countries in Asia and Africa. The phylogeny of C. neoformans reveals a recent exponential population expansion, consistent with the increase in the number of susceptible hosts. In our study population, this expansion has been driven by three sub-clades of the C. neoformans VNIa lineage; VNIa-4, VNIa-5 and VNIa-93. These three sub-clades account for 91% of clinical isolates sequenced in our study. Combining the genome data with clinical information, we find that the VNIa-93 sub-clade, the most common sub-clade in Uganda and Malawi, was associated with better outcomes than VNIa-4 and VNIa-5, which predominate in Southeast Asia. This study lays the foundation for further work investigating the dominance of VNIa-4, VNIa-5 and VNIa-93 and the association between lineage and clinical phenotype.

Prediction and identification of recurrent genomic rearrangements that generate chimeric chromosomes in Saccharomyces cerevisiae

Genomes are dynamic structures. Different mechanisms participate in the generation of genomic rearrangements. One of them is nonallelic homologous recombination (NAHR). This rearrangement is generated by recombination between pairs of repeated sequences with high identity. We analyzed rearrangements mediated by repeated sequences located in different chromosomes. Such rearrangements generate chimeric chromosomes. Potential rearrangements were predicted by localizing interchromosomal identical repeated sequences along the nuclear genome of the Saccharomyces cerevisiae S288C strain. Rearrangements were identified by a PCR-based experimental strategy. PCR primers are located in the unique regions bordering each repeated region of interest. When the PCR is performed using forward primers from one chromosome and reverse primers from another chromosome, the break point of the chimeric chromosome structure is revealed. In all cases analyzed, the corresponding chimeric structures were found. Furthermore, the nucleotide sequence of chimeric structures was obtained, and the origin of the unique regions bordering the repeated sequence was located in the expected chromosomes, using the perfect-match genomic landscape strategy (PMGL). Several chimeric structures were searched in colonies derived from single cells. All of the structures were found in DNA isolated from each of the colonies. Our findings indicate that interchromosomal rearrangements that generate chimeric chromosomes are recurrent and occur, at a relatively high frequency, in cell populations of S. cerevisiae.

Targeted metagenomics approach to capture the biodiversity of Saccharomyces genus in wild environments

The species of the genus Saccharomyces are commonly inhabiting tree bark and the surrounding soil, but their abundance have likely been underestimated due to biases in culturing methods. Metagenomic studies have so far been unable to detect Saccharomyces species in wild environments. Here, we sequenced the mycobiome of soils surrounding different trees at various altitudes in the Italian Alps. To survey for yeasts species belonging to Saccharomyces genus rather than other fungal species, we performed a selectivity step involving the isolation of the internal transcribed spacer (ITS) region that is specific to this yeast group. Reads mapping to Saccharomyces species were detected in all soil samples, including reads for S. mikatae and for S. eubayanus. ITS1 alignment of the S. cerevisiae, S. paradoxus and S. kudriavzevii sequences showed up to three base pair polymorphisms with other known strains, indicating possible new lineages. Basidiomycetous fungi were still the dominant species, compared to the Ascomycota, but the selectivity step allowed for the first time the detection and study of the biodiversity of the Saccharomyces species in their natural environment.

Genome resources for the stem and bark canker pathogens Corinectria fuckeliana, Neonectria hederae and N. punicea

Corinectria fuckeliana, Neonectria hederae, and N. punicea are fungi in the family Nectriaceae that cause canker diseases of numerous hardwood trees, conifers, and woody perennials, often leading to plant mortality. Here, we report draft genome sequences for these three phytopathogenic fungal species. The genome sizes are consistent with those reported for other members of the Nectriaceae (28 to 43 Mb). These are the first genome resources available for C. fuckeliana, N. hederae, and N. punicea. These genome sequences may provide insights into the mechanisms of virulence and pathogenicity employed by these three destructive plant pathogens, and are resources suitable for the development of molecular markers that could be used for species identification, diagnostic tools and barcodes, and population studies.

Colletotrichum orbiculare MTF4 Is a Key Transcription Factor Downstream of MOR Essential for Plant Signal-Dependent Appressorium Development and Pathogenesis

The cucumber anthracnose fungus Colletotrichum orbiculare forms a specialized infection structure, called an appressorium. Appressorium differentiation relies on fungal perception of physical and biochemical signals at the plant surface. Our previous report showed that the morphogenesis-related NDR (nuclear Dbf2-related) kinase pathway (MOR) is crucial for translating plant-derived signals for appressorium development. Here, we focused on identifying transcriptional regulators downstream of MOR that are involved in plant signal sensing and transduction for appressorium development. Based on whole-genome transcript profiling, we identified a Zn(II)₂Cys₆ transcription factor, CoMTF4, as a potential downstream factor of MOR. CoMTF4 was expressed in planta rather than in vitro under the control of the NDR kinase CoCbk1. Phenotypes of comtf4 mutants, strains with constitutively active CoCbk1 and strains with constitutive overexpression of CoMTF4 suggested that CoMtf4 acts downstream of MOR. Furthermore, nuclear localization of CoMtf4 was dependent on the MOR and responsive to plant-derived signals that lead to appressorium morphogenesis. Thus, we conclude that CoMtf4 is a transcription factor downstream of MOR that is essential for appressorium morphogenesis and pathogenesis and is regulated in response to plant-derived signals. This study provides insights into fungal sensing of plant signals and subsequent responses critical for appressorium formation.

Evolutionary Analysis of the F-Box Gene Family in Saccharomycetaceae

F-box proteins are a core component of Skp1-Cul1-F-box (SCF) ubiquitin/ligase complexes and are involved in a lot of cellular processes in yeasts. However, the current knowledge of the molecular evolution of the F-box gene family in yeasts remains unclear. In this study, 136 F-box genes were identified in 10 yeast species of the Saccharomycetaceae. In addition to the F-box domain, the other six domains were identified in these F-box proteins. The evolutionary history of F-box gene numbers in 10 Saccharomycetaceae yeasts was reconstructed. Whole-genome duplication, interspersed repeats, and gene loss events were inferred. These events contributed to F-box gene number variation in the 10 yeast species. Eighty-seven and 33 positively selected sites were detected in program Selecton and Datamonkey web-server, respectively. Three of them were considered the significant positively selected sites, and 23 of them had changed radically in amino acid properties by using TreeSAAP. We investigated F-box gene number variation and underlying mechanisms, and selection patterns, all of which were beneficial to deeply understand genome evolution and figure out the function of the F-box proteins.

Genome Sequencing

Strategies for sequencing fungal genomes on next-generation sequencing (NGS) platforms depend on the characteristics of the genome of the targeted species, quantity and quality of the genomic DNA, and cost considerations. Massively parallel sequencing with sequencing by synthesis (SBS) approach by Illumina produces terabases of short read sequences (i.e., ~300 bp) in a time and cost-effective manner, though the read length can limit the assembly particularly in repetitive regions. The single molecule, real-time (SMRT) sequencing approach by Pacific Biosciences (PacBio) produces longer reads (i.e., ~12,500 bp) which can facilitate de novo assembly of genomes that contain long repetitive sequences, though due to the lower-throughput of this platform achieving the coverage needed for assembly is more expensive than by SBS. Additionally, the Illumina SBS platforms can handle low quantity/quality of genomic DNA materials, while the SMRT system requires undamaged long DNA fragments as input to ensure that high-quality data is produced. Both platforms are discussed in this chapter including key decision-making points.

Genome-enhanced detection and identification of fungal pathogens responsible for pine and poplar rust diseases

Biosurveillance is a proactive approach that may help to limit the spread of invasive fungal pathogens of trees, such as rust fungi which have caused some of the world's most damaging diseases of pines and poplars. Most of these fungi have a complex life cycle, with up to five spore stages, which is completed on two different hosts. They have a biotrophic lifestyle and may be propagated by asymptomatic plant material, complicating their detection and identification. A bioinformatics approach, based on whole genome comparison, was used to identify genome regions that are unique to the white pine blister rust fungus, Cronartium ribicola, the poplar leaf rust fungi Melampsora medusae and Melampsora larici-populina or to members of either the Cronartium and Melampsora genera. Species- and genus-specific real-time PCR assays, targeting these unique regions, were designed with the aim of detecting each of these five taxonomic groups. In total, twelve assays were developed and tested over a wide range of samples, including different spore types, different infected plant parts on the pycnio-aecial or uredinio-telial host, and captured insect vectors. One hundred percent detection accuracy was achieved for the three targeted species and two genera with either a single assay or a combination of two assays. This proof of concept experiment on pine and poplar leaf rust fungi demonstrates that the genome-enhanced detection and identification approach can be translated into effective real-time PCR assays to monitor tree fungal pathogens.

A new approach to Cas9-based genome editing in Aspergillus niger that is precise, efficient and selectable

Aspergillus niger and other filamentous fungi are widely used in industry, but efficient genetic engineering of these hosts remains nascent. For example, while molecular genetic tools have been developed, including CRISPR/Cas9, facile genome engineering of A. niger remains challenging. To address these challenges, we have developed a simple Cas9-based gene targeting method that provides selectable, iterative, and ultimately marker-free generation of genomic deletions and insertions. This method leverages locus-specific “pop-out” recombination to suppress off-target integrations. We demonstrated the effectiveness of this method by targeting the phenotypic marker albA and validated it by targeting the glaA and mstC loci. After two selection steps, we observed 100% gene editing efficiency across all three loci. This method greatly reduces the effort required to engineer the A. niger genome and overcomes low Cas9 transformations efficiency by eliminating the need for extensive screening. This method represents a significant addition to the A. niger genome engineering toolbox and could be adapted for use in other organisms. It is expected that this method will impact several areas of industrial biotechnology, such as the development of new strains for the secretion of heterologous enzymes and the discovery and optimization of metabolic pathways.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

Ultrahigh-Density Screens for Genome-Wide Yeast EMAPs in a Single Plate

Systematic measurements of genetic interactions have been used to classify gene functions and to categorize genes into protein complexes, functional pathways and biological processes. This protocol describes how to perform a high-throughput genetic interaction screen in S. cerevisiae using a variant of epistatic miniarray profiles (E-MAP) in which the fitnesses of 6144 colonies are measured simultaneously. We also describe the computational methods to analyze the resulting data.

An Exact Algorithm for Sorting by Weighted Preserving Genome Rearrangements

The preserving Genome Sorting Problem (pGSP) asks for a shortest sequence of rearrangement operations that transforms a given gene order into another given gene order by using rearrangement operations that preserve common intervals, i.e., groups of genes that form an interval in both given gene orders. The wpGSP is the weighted version of the problem were each type of rearrangement operation has a weight and a minimum weight sequence of rearrangement operations is sought. An exact algorithm - called CREx2 - is presented, which solves the wpGSP for arbitrary gene orders and the following types of rearrangement operations: inversions, transpositions, inverse transpositions, and tandem duplication random loss operations. CREx2 has a (worst case) exponential runtime, but a linear runtime for problem instances where the common intervals are organized in a linear structure. The efficiency of CREx2 and its usefulness for phylogenetic analysis is shown empirically for gene orders of fungal mitochondrial genomes.

A resource of variant effect predictions of single nucleotide variants in model organisms

The effect of single nucleotide variants (SNVs) in coding and noncoding regions is of great interest in genetics. Although many computational methods aim to elucidate the effects of SNVs on cellular mechanisms, it is not straightforward to comprehensively cover different molecular effects. To address this, we compiled and benchmarked sequence and structure-based variant effect predictors and we computed the impact of nearly all possible amino acid and nucleotide variants in the reference genomes of Homo sapiens, Saccharomyces cerevisiae and Escherichia coli Studied mechanisms include protein stability, interaction interfaces, post-translational modifications and transcription factor binding sites. We apply this resource to the study of natural and disease coding variants. We also show how variant effects can be aggregated to generate protein complex burden scores that uncover protein complex to phenotype associations based on a set of newly generated growth profiles of 93 sequenced S. cerevisiae strains in 43 conditions. This resource is available through mutfunc (www.mutfunc.com), a tool by which users can query precomputed predictions by providing amino acid or nucleotide-level variants.

Genome Sequence Resources for the Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) and the Barley Stripe Rust Pathogen (Puccinia striiformis f. sp. hordei)

Puccinia striiformis f. sp. tritici causes devastating stripe (yellow) rust on wheat and P. striiformis f. sp. hordei causes stripe rust on barley. Several P. striiformis f. sp. tritici genomes are available, but no P. striiformis f. sp. hordei genome is available. More genomes of P. striiformis f. sp. tritici and P. striiformis f. sp. hordei are needed to understand the genome evolution and molecular mechanisms of their pathogenicity. We sequenced P. striiformis f. sp. tritici isolate 93-210 and P. striiformis f. sp. hordei isolate 93TX-2, using PacBio and Illumina technologies and RNA sequencing. Their genomic sequences were assembled to contigs with high continuity and showed significant structural differences. The circular mitochondria genomes of both were complete. These genomes provide high-quality resources for deciphering the genomic basis of rapid evolution and host adaptation, identifying genes for avirulence and other important traits, and studying host-pathogen interactions.

A protocol for introduction of multiple genetic modifications in Saccharomyces cerevisiae using CRISPR/Cas9

Here, two methods are described for efficient genetic modification of Saccharomyces cerevisiae using CRISPR/Cas9. The first method enables the modification of a single genetic locus using in vivo assembly of a guide RNA (gRNA) expression plasmid without the need for prior cloning. A second method using in vitro assembled plasmids that could contain up to two gRNAs was used to simultaneously introduce up to six genetic modifications (e.g. six gene deletions) in a single transformation step by transforming up to three gRNA expression plasmids simultaneously. The method is not only suitable for gene deletion but is also applicable for in vivo site-directed mutagenesis and integration of multiple DNA fragments in a single locus. In all cases, the strain transformed with the gRNA expression plasmids was equipped with a genomic integration of Spcas9, leading to strong and constitutive expression of SpCas9. The protocols detailed here have been streamlined to be executed by virtually any yeast molecular geneticist.

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.

Molecular epidemiology of environmental Cryptococcus species isolates based on amplified fragment length polymorphism

OBJECTIVE

Cryptococcosis is a major opportunistic fungal infection caused by members of the genus Cryptococcus, mainly those belonging to the Cryptococcus neoformans/Cryptococcus gattii species complexes. Here, we report a comprehensive molecular epidemiological study of the environmental distribution of Cryptococcus isolates in Shiraz, Iran with review of litreature.

METHOD

A total of 406 samples were obtained from Eucalyptus trees and 139 samples from pigeon droppings. Cryptococcus species identification and genotyping were performed by amplified fragment length polymorphism (AFLP) fingerprinting sequencing and sequencing of the ITS rDNA region.

RESULTS

Majority of the isolates belonged to the Naganishia taxon (n=69) including N. albida (formerly C. albidus, n=62), N. globosa (formerly C. saitoi, n=4), N. adeliensis (formerly C. adeliensis, n=2), N. diffluens (formerly C. diffluens, n=1), and the identified C. neoformans isolates (n=25) belonged to genotype AFLP1/VNI (n=22) and AFLP1B/VNII (n=3).

CONCLUSION

More research efforts should be employed to isolate C. gattii species complex from environmental niches in Iran and provide additional evidence related to novel molecular types.

Convergent evolution of complex genomic rearrangements in two fungal meiotic drive elements

Meiotic drive is widespread in nature. The conflict it generates is expected to be an important motor for evolutionary change and innovation. In this study, we investigated the genomic consequences of two large multi-gene meiotic drive elements, Sk-2 and Sk-3, found in the filamentous ascomycete Neurospora intermedia. Using long-read sequencing, we generated the first complete and well-annotated genome assemblies of large, highly diverged, non-recombining regions associated with meiotic drive elements. Phylogenetic analysis shows that, even though Sk-2 and Sk-3 are located in the same chromosomal region, they do not form sister clades, suggesting independent origins or at least a long evolutionary separation. We conclude that they have in a convergent manner accumulated similar patterns of tandem inversions and dense repeat clusters, presumably in response to similar needs to create linkage between genes causing drive and resistance.

Bioinformatic characterisation of the effector repertoire of the strawberry pathogen Phytophthora cactorum

The oomycete pathogen Phytophthora cactorum causes crown rot, a major disease of cultivated strawberry. We report the draft genome of P. cactorum isolate 10300, isolated from symptomatic Fragaria x ananassa tissue. Our analysis revealed that there are a large number of genes encoding putative secreted effectors in the genome, including nearly 200 RxLR domain containing effectors, 77 Crinklers (CRN) grouped into 38 families, and numerous apoplastic effectors, such as phytotoxins (PcF proteins) and necrosis inducing proteins. As in other Phytophthora species, the genomic environment of many RxLR and CRN genes differed from core eukaryotic genes, a hallmark of the two-speed genome. We found genes homologous to known Phytophthora infestans avirulence genes including Avr1, Avr3b, Avr4, Avrblb1 and AvrSmira2 indicating effector sequence conservation between Phytophthora species of clade 1a and clade 1c. The reported P. cactorum genome sequence and associated annotations represent a comprehensive resource for avirulence gene discovery in other Phytophthora species from clade 1 and, will facilitate effector informed breeding strategies in other crops.

Characterisation of pathogen-specific regions and novel effector candidates in Fusarium oxysporum f. sp. cepae

A reference-quality assembly of Fusarium oxysporum f. sp. cepae (Foc), the causative agent of onion basal rot has been generated along with genomes of additional pathogenic and non-pathogenic isolates of onion. Phylogenetic analysis confirmed a single origin of the Foc pathogenic lineage. Genome alignments with other F. oxysporum ff. spp. and non pathogens revealed high levels of syntenic conservation of core chromosomes but little synteny between lineage specific (LS) chromosomes. Four LS contigs in Foc totaling 3.9 Mb were designated as pathogen-specific (PS). A two-fold increase in segmental duplication events was observed between LS regions of the genome compared to within core regions or from LS regions to the core. RNA-seq expression studies identified candidate effectors expressed in planta, consisting of both known effector homologs and novel candidates. FTF1 and a subset of other transcription factors implicated in regulation of effector expression were found to be expressed in planta.