Long-Read Genome Sequence Resource of Ascochyta versabilis Causing Leaf Spot Disease in Pseudostellaria heterophylla

Ascochyta versabilis is the fungal pathogen that causes the severe leaf spot disease of Pseudostellaria heterophylla (Miq.) Pax, a vital Chinese herbal plant. Here, we deployed PacBio single-molecule real-time long-read sequencing technology to generate a near-complete genome assembly for the A. versabilis KC1 strain and obtained a total of 9.80 Gb raw reads. These reads were processed into a 41.05 Mb genome assembly containing 95 contigs with N50 of 1.70 Mb and a maximum length of 3.93 Mb. A total of 10,457 gene models, of which 1,004 encode putatively secreted proteins, were identified in the genome. This high-quality genome assembly and gene annotation resource will facilitate the institution of functional genetic studies aimed at providing a better insight into the infection mechanisms of A. versabilis to support the development of effective control strategies for leaf spot disease of P. heterophylla.

Identifying candidate Aspergillus pathogenicity factors by annotation frequency

BACKGROUND

Members of the genus Aspergillus display a variety of lifestyles, ranging from saprobic to pathogenic on plants and/or animals. Increased genome sequencing of economically important members of the genus permits effective use of "-omics" comparisons between closely related species and strains to identify candidate genes that may contribute to phenotypes of interest, especially relating to pathogenicity. Protein-coding genes were predicted from 216 genomes of 12 Aspergillus species, and the frequencies of various structural aspects (exon count and length, intron count and length, GC content, and codon usage) and functional annotations (InterPro, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes terms) were compared.

RESULTS

Using principal component analyses, the three sets of functional annotations for each strain were clustered by species. The species clusters appeared to separate by pathogenicity on plants along the first dimensions, which accounted for over 20% of the variance. More annotations for genes encoding pectinases and secondary metabolite biosynthetic enzymes were assigned to phytopathogenic strains from species such as Aspergillus flavus. In contrast, Aspergillus fumigatus strains, which are pathogenic to animals but not plants, were assigned relatively more terms related to phosphate transferases, and carbohydrate and amino-sugar metabolism. Analyses of publicly available RNA-Seq data indicated that one A. fumigatus protein among 17 amino-sugar processing candidates, a hexokinase, was up-regulated during co-culturing with human immune system cells.

CONCLUSION

Genes encoding hexokinases and other proteins of interest may be subject to future manipulations to further refine understanding of Aspergillus pathogenicity factors.

A haplotype-phased genome of wheat stripe rust pathogen Puccinia striiformis f. sp. tritici, race PST-130 from the Western USA

More virulent and aggressive races of Puccinia striiformis f. sp. tritici (Pst), the pathogen causing wheat stripe rust, have been spreading around the world since 2000 causing large grain yield losses. A better understanding of the genome and genetic diversity of these new Pst races will be useful to develop new strategies to ameliorate these losses. In this study, we generated an improved genome assembly of a post-2000 virulent race from the Western USA designated as PST-130. We implemented a haplotype phasing strategy using the diploid-aware assembler, Falcon-Unzip and new long-read technology from PacBio to phase the two genomes of this dikaryotic organism. The combination of these new technologies resulted in an improved PST-130 assembly with only 151 contigs (85.4 Mb, N50 of 1.44 Mb), and a complementary assembly (haplotigs) with 458 contigs (65.9 Mb, N50 of 0.235 Mb, PRJNA650506). This new assembly improved gene predictions resulting in 228 more predicted complete genes than in the initial Illumina assembly (29,178 contigs, N50 of 5 kb). The alignment of the non-repetitive primary and haplotig contigs revealed and average of 5.22 SNP/kb, with 39.1% showing <2 SNP/kb and 15.9% >10 SNP/kb. This large divergent regions may represent introgressions of chromosome segments from more divergent Pst races in regions where a complete sexual cycle and recombination are possible. We hypothesize that some of the divergent regions in PST-130 may be related to the European "Warrior" race PST-DK0911 because this genome is more similar to PST-130 (3.18 SNP/kb) than to the older European race PST-104E (3.75 SNP/kb). Complete phasing of additional Pst genomes or sequencing individual nuclei will facilitate the tracing of the haploid genomes introduced by the new Pst races into local populations.

An intergenic "safe haven" region in Cryptococcus neoformans serotype D genomes

Cryptococcus neoformans is an opportunistic human fungal pathogen and serves as a model organism for studies of eukaryotic microbiology and microbial pathogenesis. C. neoformans species complex is classified into serotype A, serotype D, and AD hybrids, which are currently considered different subspecies. Different serotype strains display varied phenotypes, virulence, and gene regulation. Genetic investigation of important pathways is often performed in both serotype A and D reference strains in order to identify diversification or conservation of the interrogated signaling network. Many genetic tools have been developed for C. neoformans serotype A reference strain H99, including the gene free "safe haven" (SH) regions for DNA integration identified based on genomic features. However, no such a genomic safe haven region has been identified in serotype D strains. Here, capitalizing on the available genomic, transcriptomic, and chromatin data, we identified an intergenic region named as SH3 for the serotype D reference strains JEC21 and XL280. We also designed a sgRNA and a vector facilitating any alien gene integration into SH3 through a CRISPR-Cas9 system. We found that gene inserted in this region complemented the corresponding gene deletion mutant. Fluorescent reporter gene inserted in SH3 can also be expressed efficiently. Insertion in SH3 itself did not alter the expression of adjacent genes and did not affect the growth or mating of C. neoformans. Thus, SH3 provides a resource for genetic manipulations in serotype D strains and will facilitate comparative analyses of gene functions in this species complex. In addition, the incorporation of the multi-omic data in our selection of the safe haven region could help similar studies in other organisms.

Penicillium expansum: biology, omics, and management tools for a global postharvest pathogen causing blue mould of pome fruit

Blue mould, caused primarily by Penicillium expansum, is a major threat to the global pome fruit industry, causing multimillion-dollar losses annually. The blue mould fungus negatively affects fruit quality, thereby reducing fresh fruit consumption, and significantly contributes to food loss. P. expansum also produces an array of mycotoxins that are detrimental to human health. Management options are limited and the emergence of fungicide-resistant Penicillium spp. makes disease management difficult, therefore new approaches and tools are needed to combat blue mould in storage. This species profile comprises a comprehensive literature review of this aggressive pathogen associated with pomes (apple, pear, quince), focusing on biology, mechanisms of disease, control, genomics, and the newest developments in disease management.

TAXONOMY

Penicillium expansum Link 1809. Domain Eukaryota, Kingdom Fungi, Phylum Ascomycota, Subphylum Pezizomycotina, Class Eurotiomycetes, Subclass: Eurotiomycetidae, Order Eurotiales; Family Trichocomaceae, Genus Penicillium, Species expansum.

BIOLOGY

A wide host range necrotrophic postharvest pathogen that requires a wound (e.g., stem pull, punctures, bruises, shoulder cracks) or natural openings (e.g., lenticel, stem end, calyx sinus) to gain ingress and infect.

TOXINS

Patulin, citrinin, chaetoglobosins, communesins, roquefortine C, expansolides A and B, ochratoxin A, penitrem A, rubratoxin B, and penicillic acid.

HOST RANGE

Primarily apples, European pear, Asian pear, medlar, and quince. Blue mould has also been reported on stone fruits (cherry, plum, peach), small fruits (grape, strawberry, kiwi), and hazel nut.

DISEASE SYMPTOMS

Blue mould initially appears as light tan to dark brown circular lesions with a defined margin between the decayed and healthy tissues. The decayed tissue is soft and watery, and blue-green spore masses appear on the decayed area, starting at the infection site and radiating outward as the decayed area ages.

DISEASE CONTROL

Preharvest fungicides with postharvest activity and postharvest fungicides are primarily used to control decay. Orchard and packinghouse sanitation methods are also critical components of an integrated pest management strategy.

USEFUL WEBSITES

Penn State Tree Fruit Production Guide (https://extension.psu.edu/forage-and-food-crops/fruit), Washington State Comprehensive Tree Fruit (http://treefruit.wsu.edu/crop-protection/disease-management/blue-mold/), The Apple Rot Doctor (https://waynejurick.wixsite.com/applerotdr), penicillium expansum genome sequences and resources (https://www.ncbi.nlm.nih.gov/genome/browse/#!/eukaryotes/11336/).

Genome-scale analyses and characteristics of putative pathogenicity genes of Stagonosporopsis cucurbitacearum, a pumpkin gummy stem blight fungus

Outbreaks of gummy stem blight (GSB), an emerging seed pumpkin disease, have increased in number and have become more widespread in recent years. Previously we reported that Stagonosporopsis cucurbitacearum (Sc.) is the dominant fungal cause of pumpkin seedling GSB in Northeast China, where it has greatly reduced crop yields in that region. Here, high-throughput whole-genome sequencing and assembly of the Sc. genome were conducted toward revealing pathogenic molecular regulatory mechanisms involved in fungal growth and development. Zq-1 as representative Sc. strain, DNA of Zq-1was prepared for genomic sequencing, we obtained 5.24 Gb of high-quality genomic sequence data via PacBio RS II sequencing. After sequence data was processed to filter out low quality reads, a hierarchical genome-assembly process was employed that generated a genome sequence of 35.28 Mb in size. A total of 9844 genes were predicted, including 237 non-coding RNAs, 1024 genes encoding proteins with signal peptides, 2066 transmembrane proteins and 756 secretory proteins.Transcriptional identification revealed 54 differentially expressed secretory proteins. Concurrently, 605, 130 and 2869 proteins were matched in the proprietary databases Carbohydrate-Active EnZymes database (CAZyme), Transporter Classification Database (TCDB) and Pathogen-Host Interactions database (PHI), respectively. And 96 and 36 DEGs were identified form PHI database and CAZyme database, respectively. In addition, contig00011.93 was an up-regulated DEG involving ATP-binding cassette metabolism in the procession of infection. In order to test relevance of gene predictions to GSB, DEGs with potential pathogenic relevance were revealed through transcriptome data analysis of Sc. strains pre- and post-infection of pumpkin. Interestingly, Sc. and Leptosphaeria maculans (Lm.) exhibited relatively similar with genome lengths, numbers of protein-coding genes and other characteristics. This work provides a foundation for future exploration of additional Sc. gene functions toward the development of more effective GSB control strategies.

Trajectory and uniqueness of mutational signatures in yeast mutators

The acquisition of mutations plays critical roles in adaptation, evolution, senescence, and tumorigenesis. Massive genome sequencing has allowed extraction of specific features of many mutational landscapes but it remains difficult to retrospectively determine the mechanistic origin(s), selective forces, and trajectories of transient or persistent mutations and genome rearrangements. Here, we conducted a prospective reciprocal approach to inactivate 13 single or multiple evolutionary conserved genes involved in distinct genome maintenance processes and characterize de novo mutations in 274 diploid Saccharomyces cerevisiae mutation accumulation lines. This approach revealed the diversity, complexity, and ultimate uniqueness of mutational landscapes, differently composed of base substitutions, small insertions/deletions (InDels), structural variants, and/or ploidy variations. Several landscapes parallel the repertoire of mutational signatures in human cancers while others are either novel or composites of subsignatures resulting from distinct DNA damage lesions. Notably, the increase of base substitutions in the homologous recombination-deficient Rad51 mutant, specifically dependent on the Polζ translesion polymerase, yields COSMIC signature 3 observed in BRCA1/BRCA2-mutant breast cancer tumors. Furthermore, "mutome" analyses in highly polymorphic diploids and single-cell bottleneck lineages revealed a diverse spectrum of loss-of-heterozygosity (LOH) signatures characterized by interstitial and terminal chromosomal events resulting from interhomolog mitotic cross-overs. Following the appearance of heterozygous mutations, the strong stimulation of LOHs in the rad27/FEN1 and tsa1/PRDX1 backgrounds leads to fixation of homozygous mutations or their loss along the lineage. Overall, these mutomes and their trajectories provide a mechanistic framework to understand the origin and dynamics of genome variations that accumulate during clonal evolution.

Genomic instability in an interspecific hybrid of the genus Saccharomyces: a matter of adaptability

Ancient events of polyploidy have been linked to huge evolutionary leaps in the tree of life, while increasing evidence shows that newly established polyploids have adaptive advantages in certain stress conditions compared to their relatives with a lower ploidy. The genus Saccharomyces is a good model for studying such events, as it contains an ancient whole-genome duplication event and many sequenced Saccharomyces cerevisiae are, evolutionary speaking, newly formed polyploids. Many polyploids have unstable genomes and go through large genome erosions; however, it is still unknown what mechanisms govern this reduction. Here, we sequenced and studied the natural S. cerevisiae × Saccharomyces kudriavzevii hybrid strain, VIN7, which was selected for its commercial use in the wine industry. The most singular observation is that its nuclear genome is highly unstable and drastic genomic alterations were observed in only a few generations, leading to a widening of its phenotypic landscape. To better understand what leads to the loss of certain chromosomes in the VIN7 cell population, we looked for genetic features of the genes, such as physical interactions, complex formation, epistatic interactions and stress responding genes, which could have beneficial or detrimental effects on the cell if their dosage is altered by a chromosomal copy number variation. The three chromosomes lost in our VIN7 population showed different patterns, indicating that multiple factors could explain the mechanisms behind the chromosomal loss. However, one common feature for two out of the three chromosomes is that they are among the smallest ones. We hypothesize that small chromosomes alter their copy numbers more frequently as a low number of genes is affected, meaning that it is a by-product of genome instability, which might be the chief driving force of the adaptability and genome architecture of this hybrid.

Effect of H2A.Z deletion is rescued by compensatory mutations in Fusarium graminearum

Fusarium head blight is a destructive disease of grains resulting in reduced yields and contamination of grains with mycotoxins worldwide; Fusarium graminearum is its major causal agent. Chromatin structure changes play key roles in regulating mycotoxin biosynthesis in filamentous fungi. Using a split-marker approach in three F. graminearum strains INRA156, INRA349 and INRA812 (PH-1), we knocked out the gene encoding H2A.Z, a ubiquitous histone variant reported to be involved in a diverse range of biological processes in yeast, plants and animals, but rarely studied in filamentous fungi. All ΔH2A.Z mutants exhibit defects in development including radial growth, sporulation, germination and sexual reproduction, but with varying degrees of severity between them. Heterogeneity of osmotic and oxidative stress response as well as mycotoxin production was observed in ΔH2A.Z strains. Adding-back wild-type H2A.Z in INRA349ΔH2A.Z could not rescue the phenotypes. Whole genome sequencing revealed that, although H2A.Z has been removed from the genome and the deletion cassette is inserted at H2A.Z locus only, mutations occur at other loci in each mutant regardless of the genetic background. Genes affected by these mutations encode proteins involved in chromatin remodeling, such as the helicase Swr1p or an essential subunit of the histone deacetylase Rpd3S, and one protein of unknown function. These observations suggest that H2A.Z and the genes affected by such mutations are part or the same genetic interaction network. Our results underline the genetic plasticity of F. graminearum facing detrimental gene perturbation. These findings suggest that intergenic suppressions rescue deleterious phenotypes in ΔH2A.Z strains, and that H2A.Z may be essential in F. graminearum. This assumption is further supported by the fact that H2A.Z deletion failed in another Fusarium spp., i.e., the rice pathogen Fusarium fujikuroi.

Transcriptome-wide expression profiling of Sporothrix schenckii yeast and mycelial forms and the establishment of the Sporothrix Genome DataBase

Sporothrix schenckii is a dimorphic fungus existing as mould in the environment and as yeast in the host. The morphological shift between mycelial/yeast phases is crucial for its virulence, but the transcriptional networks implicated in dimorphic transition are still not fully understood. Here, we report the global transcriptomic differences occurring between mould and yeast phases of S. schenckii, including changes in gene expression profiles associated with these distinct cellular phenotypes. Moreover, we also propose a new genome annotation, which reveals a more complex transcriptional architecture than previously assumed. Using RNA-seq, we identified a total of 17 307 genes, of which 11 217 were classified as protein-encoding genes, whereas 6090 were designated as non-coding RNAs (ncRNAs). Approximately ~71 % of all annotated genes were found to overlap and the different-strand overlapping type was the most common. Gene expression analysis revealed that 8795 genes were differentially regulated among yeast and mould forms. Differential gene expression was also observed for antisense ncRNAs overlapping neighbouring protein-encoding genes. The release of transcriptome-wide data and the establishment of the Sporothrix Genome DataBase (http://sporothrixgenomedatabase.unime.it) represent an important milestone for Sporothrix research, because they provide a strong basis for future studies on the molecular pathways involved in numerous biological processes.

A Signature-Tagged Mutagenesis (STM)-based murine-infectivity assay for Cryptococcus neoformans

Signature-tagged mutagenesis (STM) is a high-throughput genetic technique that can be used to investigate the function of genes by constructing a large number of mutant strains with unique DNA identification tags, pooling them, and screening them for a particular phenotypic trait. STM was first designed for the identification of genes that contribute to the virulence or infectivity of a pathogen in its host. Recently, this method has also been applied for the identification of mutants with specific phenotypes, such as antifungal drug resistance and proliferation. In the present study, we describe an STM method for the identification of genes contributing to the infectivity of Cryptococcus neoformans using a mutant library, in which each strain was tagged with a unique DNA sequence.

Protocol for Genome-Scale Reconstruction and Melanogenesis Analysis of Exophiala dermatitidis

Exophiala dermatitidis is a polyextremotolerant fungus with a small genome, thus suitable as a model system for melanogenesis and carotenogensis. A genome-scale model, iEde2091, is reconstructed to increase metabolic understanding and used in a shadow price analysis of pigments, as detailed here. Important to this reconstruction is OptFill, a recently developed alternative gap-filling method useful in the holistic and conservative reconstruction of genome-scale models of metabolism, particularly for understudied organisms like E. dermatitidis where gaps in metabolic knowledge are abundant. For complete details on the use and execution of this protocol, please refer to Schroeder and Saha (2020) and Schroeder et al. (2020).

Dual-PKS Cluster for Biosynthesis of a Light-Induced Secondary Metabolite Found from Genome Sequencing of Hyphodiscus hymeniophilus Fungus

Filamentous fungi are known producers of important secondary metabolites. In spite of this, the majority of these organisms have not been studied at the genome level, leaving many of the bioactive molecules they produce undiscovered. In this study, we explore the secondary metabolite potential of an understudied fungus, Hyphodiscus hymeniophilus. By sequencing and assembling the first genome from this genus, we show that this fungus has genes for at least 20 natural products and that many of these products are likely novel. One of these metabolites is identified: a new, red-pigmented member of the azaphilone class, hyphodiscorubrin. We show that this metabolite is only produced when the fungus is grown in the light. Furthermore, the biosynthetic gene cluster of hyphodiscorubrin is identified though homology to other known azaphilone producing clusters.

High efficiency CRISPR/Cas9 genome editing system with an eliminable episomal sgRNA plasmid in Pichia pastoris

Pichia pastoris is a methylotrophic yeast in which host heterologous expression of proteins has been developed owing to the strong inducible alcohol oxidase promoter (PAOX1). However, it is difficult to manipulate the genome in P. pastoris. Based on previous attempts to apply the CRISPR/Cas9 system in P. pastoris, a CRISPR/Cas9 system with episomal sgRNA plasmid was developed and 100 % genome editing efficiency, high multicopy gene editing and stable multigene editing were obtained without a sharp decline caused by multi-sgRNA. And 28/34 (∼82 %) sgRNAs tested were effective. The CGG may have a slightly higher and more stable cleavage efficiency than the other three NGG motifs, and a low GC content may be preferable for higher cleavage efficiency. This provides researchers with a stable genome editing tool that shows a high editing efficiency, shortening the experimentation period. Furthermore, we introduced dCas9 into P. pastoris and achieved target gene interference, expanding the CRISPR/Cas9 toolbox in P. pastoris.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

Genomic diversity of the human pathogen Paracoccidioides across the South American continent

Paracoccidioidomycosis (PCM) is a life-threatening systemic mycosis widely reported in the Gran Chaco ecosystem. The disease is caused by different species from the genus Paracoccidioides, which are all endemic to South and Central America. Here, we sequenced and analyzed 31 isolates of Paracoccidioides across South America, with particular focus on isolates from Argentina and Paraguay. The de novo sequenced isolates were compared with publicly available genomes. Phylogenetics and population genomics revealed that PCM in Argentina and Paraguay is caused by three distinct Paracoccidioides genotypes, P. brasiliensis (S1a and S1b) and P. restrepiensis (PS3). P. brasiliensis S1a isolates from Argentina are frequently associated with chronic forms of the disease. Our results suggest the existence of extensive molecular polymorphism among Paracoccidioides species, and provide a framework to begin to dissect the connection between genotypic differences in the pathogen and the clinical outcomes of the disease.

A simple approach to mediate genome editing in the filamentous fungus Trichoderma reesei by CRISPR/Cas9-coupled in vivo gRNA transcription

OBJECTIVE

To simplify CRISPR/Cas9 genome editing in the industrial filamentous fungus Trichoderma reesei based on in vivo guide RNA (gRNA) transcription.

RESULTS

Two putative RNA polymerase III U6 snRNA genes were identified in the genome of T. reesei QM6a by BLASTN using Myceliophthora. thermophila U6 snRNA gene as the template. The regions approximately 500 bp upstream of two U6 genes were efficient promoters for the in vivo expression of gRNA. The CRISPR system consisting of Cas9 and in vivo synthesized gRNA under control of the T. reesei U6 snRNA promoters was sufficient to cause a frameshift mutation in the ura5 gene via non-homologous end-joining-mediated events.

CONCLUSIONS

We report a simple gene editing method using a CRISPR/Cas9-coupled in vivo gRNA transcription system in T. reesei.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.

Genome-wide mRNA and miRNA analysis in the early stages of germ tube outgrowth in Coprinopsis cinerea

Coprinopsis cinerea is a model mushroom-forming basidiomycete which produces basidiospores during sexual reproduction. This fungus is widely used to study fruiting body formation and development. Molecular mechanisms controlling its growth from vegetative mycelium to multicellular mature fruiting body have been studied extensively. However, little is known about the underlying biological processes during germ tube outgrowth or the transition from basidiospores to multinucleate hyphae. To better understand sexual spore germination in fungi, here we examined the time-dependent cellular events at resting, germinating and fully germinated basidiospores of C. cinerea by genome-wide transcriptional and post-transcriptional analyses and by carbohydrate composition analysis. Our results revealed a high demand of protein degradation, and biosynthesis of various compounds at the early stage of basidiospore gemination and dynamic changes of carbohydrate metabolism throughout the germination process. Seven microRNA-like RNAs (milRNAs) were identified in the resting basidiospores of C. cinerea, six of which were basidiospore-specific. Glycogen and trehalose were shown to be the carbon sources supporting the initiation of germ tube outgrowth. One basidiospore-specific milRNA, cci-milR-37, was found to be a potential regulator of glycogen metabolic pathways related to vegetative hyphal growth. Our results demonstrated the mRNA and miRNA-mediated regulation on energy production, protein and carbohydrate metabolisms at the early developmental stages of germ tube to form totipotent hyphae. To our knowledge, this is the first study to show the roles of miRNAs in mushroom basidiospore germination and out-growth.

Convergent genes shape budding yeast pericentromeres

The three-dimensional architecture of the genome governs its maintenance, expression and transmission. The cohesin protein complex organizes the genome by topologically linking distant loci, and is highly enriched in specialized chromosomal domains surrounding centromeres, called pericentromeres1-6. Here we report the three-dimensional structure of pericentromeres in budding yeast (Saccharomyces cerevisiae) and establish the relationship between genome organization and function. We find that convergent genes mark pericentromere borders and, together with core centromeres, define their structure and function by positioning cohesin. Centromeres load cohesin, and convergent genes at pericentromere borders trap it. Each side of the pericentromere is organized into a looped conformation, with border convergent genes at the base. Microtubule attachment extends a single pericentromere loop, size-limited by convergent genes at its borders. Reorienting genes at borders into a tandem configuration repositions cohesin, enlarges the pericentromere and impairs chromosome biorientation during mitosis. Thus, the linear arrangement of transcriptional units together with targeted cohesin loading shapes pericentromeres into a structure that is competent for chromosome segregation. Our results reveal the architecture of the chromosomal region within which kinetochores are embedded, as well as the restructuring caused by microtubule attachment. Furthermore, we establish a direct, causal relationship between the three-dimensional genome organization of a specific chromosomal domain and cellular function.

Genome mining of the citrus pathogen Elsinoë fawcettii; prediction and prioritisation of candidate effectors, cell wall degrading enzymes and secondary metabolite gene clusters

Elsinoë fawcettii, a necrotrophic fungal pathogen, causes citrus scab on numerous citrus varieties around the world. Known pathotypes of E. fawcettii are based on host range; additionally, cryptic pathotypes have been reported and more novel pathotypes are thought to exist. E. fawcettii produces elsinochrome, a non-host selective toxin which contributes to virulence. However, the mechanisms involved in potential pathogen-host interactions occurring prior to the production of elsinochrome are unknown, yet the host-specificity observed among pathotypes suggests a reliance upon such mechanisms. In this study we have generated a whole genome sequencing project for E. fawcettii, producing an annotated draft assembly 26.01 Mb in size, with 10,080 predicted gene models and low (0.37%) coverage of transposable elements. A small proportion of the assembly showed evidence of AT-rich regions, potentially indicating genomic regions with increased plasticity. Using a variety of computational tools, we mined the E. fawcettii genome for potential virulence genes as candidates for future investigation. A total of 1,280 secreted proteins and 276 candidate effectors were predicted and compared to those of other necrotrophic (Botrytis cinerea, Parastagonospora nodorum, Pyrenophora tritici-repentis, Sclerotinia sclerotiorum and Zymoseptoria tritici), hemibiotrophic (Leptosphaeria maculans, Magnaporthe oryzae, Rhynchosporium commune and Verticillium dahliae) and biotrophic (Ustilago maydis) plant pathogens. Genomic and proteomic features of known fungal effectors were analysed and used to guide the prioritisation of 120 candidate effectors of E. fawcettii. Additionally, 378 carbohydrate-active enzymes were predicted and analysed for likely secretion and sequence similarity with known virulence genes. Furthermore, secondary metabolite prediction indicated nine additional genes potentially involved in the elsinochrome biosynthesis gene cluster than previously described. A further 21 secondary metabolite clusters were predicted, some with similarity to known toxin producing gene clusters. The candidate virulence genes predicted in this study provide a comprehensive resource for future experimental investigation into the pathogenesis of E. fawcettii.

Spatial inter-centromeric interactions facilitated the emergence of evolutionary new centromeres

Centromeres of Candida albicans form on unique and different DNA sequences but a closely related species, Candida tropicalis, possesses homogenized inverted repeat (HIR)-associated centromeres. To investigate the mechanism of centromere type transition, we improved the fragmented genome assembly and constructed a chromosome-level genome assembly of C. tropicalis by employing PacBio sequencing, chromosome conformation capture sequencing (3C-seq), chromoblot, and genetic analysis of engineered aneuploid strains. Further, we analyzed the 3D genome organization using 3C-seq data, which revealed spatial proximity among the centromeres as well as telomeres of seven chromosomes in C. tropicalis. Intriguingly, we observed evidence of inter-centromeric translocations in the common ancestor of C. albicans and C. tropicalis. Identification of putative centromeres in closely related Candida sojae, Candida viswanathii and Candida parapsilosis indicates loss of ancestral HIR-associated centromeres and establishment of evolutionary new centromeres (ENCs) in C. albicans. We propose that spatial proximity of the homologous centromere DNA sequences facilitated karyotype rearrangements and centromere type transitions in human pathogenic yeasts of the CUG-Ser1 clade.

Lectins from the Edible Mushroom Agaricus bisporus and Their Therapeutic Potentials

The mushroom Agaricus bisporus secretes biologically active compounds and proteins with benefits for human health. Most reported proteins from A. bisporus are tyrosinases and lectins. Lectins are of therapeutic or pharmaceutical interest. To date, only limited information is available on A. bisporus lectins and lectin-like proteins. No therapeutic products derived from A. bisporus lectin (ABL) are available on the market despite its extensive exploration. Recently, A. bisporus mannose-binding protein (Abmb) was discovered. Its discovery enriches the information and increases the interest in proteins with therapeutic potential from this mushroom. Furthermore, the A. bisporus genome reveals the possible occurrence of other lectins in this mushroom that may also have therapeutic potential. Most of these putative lectins belong to the same lectin groups as ABL and Abmb. Their relationship is discussed. Particular attention is addressed to ABL and Abmb, which have been explored for their potential in medicinal or pharmaceutical applications. ABL and Abmb have anti-proliferative activities toward cancer cells and a stimulatory effect on the immune system. Possible scenarios for their use in therapy and modification are also presented.

Chemical, Bioactivity, and Biosynthetic Screening of Epiphytic Fungus Zasmidium pseudotsugae

We report the first secondary metabolite, 8,8′-bijuglone, obtained from pure cultures of the slow growing Douglas fir- (Pseudotsuga menziesii var. menziesii) foliage-associated fungus Zasmidium pseudotsugae. The quinone was characterized using extensive LC/MS and NMR-based spectroscopic methods. 8,8′-Bijuglone exhibited moderate antibiotic activity against Gram-positive pathogens and weak cytotoxic activity in the NCI-60 cell line panel and in our in-house human colon carcinoma (HCT-116) cell line. An analysis of the fungal genome sequence to assess its metabolic potential was implemented using the bioinformatic tool antiSMASH. In total, 36 putative biosynthetic gene clusters were found with a majority encoding for polyketides (17), followed by non-ribosomal peptides (14), terpenes (2), ribosomal peptides (1), and compounds with mixed biosynthetic origin (2). This study demonstrates that foliage associated fungi of conifers produce antimicrobial metabolites and suggests this guild of fungi may present a rich source of novel molecules.

Evolution of Aspergillus oryzae before and after domestication inferred by large-scale comparative genomic analysis

Aspergillus oryzae is an industrially useful species, of which various strains have been identified; however, their genetic relationships remain unclear. A. oryzae was previously thought to be asexual and unable to undergo crossbreeding. However, recent studies revealed the sexual reproduction of Aspergillus flavus, a species closely related to A. oryzae. To investigate potential sexual reproduction in A. oryzae and evolutionary history among A. oryzae and A. flavus strains, we assembled 82 draft genomes of A. oryzae strains used practically. The phylogenetic tree of concatenated genes confirmed that A. oryzae was monophyletic and nested in one of the clades of A. flavus but formed several clades with different genomic structures. Our results suggest that A. oryzae strains have undergone multiple inter-genomic recombination events between A. oryzae ancestors, although sexual recombination among domesticated species did not appear to have occurred during the domestication process, at least in the past few decades. Through inter- and intra-cladal comparative analysis, we found that evolutionary pressure induced by the domestication of A. oryzae appears to selectively cause non-synonymous and gap mutations in genes involved in fermentation characteristics, as well as intra-genomic rearrangements, with the conservation of industrially useful catalytic enzyme-encoding genes.

Natural selection drives population divergence for local adaptation in a wheat pathogen

Evolution favors the emergence of locally-adapted optimum phenotypes that are likely to differ across a wide array of environmental conditions. The emergence of favorable adaptive characteristics is accelerated in agricultural pathogens due to the unique properties of agro-ecosystems. We performed a Q_ST - F_ST comparison using 164 strains of Parastagonospora nodorum sampled from eight global field populations to disentangle the predominant evolutionary forces driving population divergence in a wheat pathogen. We used digital image analysis to obtain quantitative measurements of growth rate and melanization at different temperatures and under different fungicide concentrations in a common garden experiment. F_ST measures were based on complete genome sequences obtained for all 164 isolates. Our analyses indicated that all measured traits were under selection. Growth rates at 18 °C and 24 °C were under stabilizing selection (Q_ST < F_ST), while diversifying selection (Q_ST > F_ST) was the predominant evolutionary force affecting growth under fungicide and high temperature stress. Stabilizing selection (Q_ST < F_ST) was the predominant force affecting melanization across the different environments. Melanin production increased at 30 °C but was negatively correlated with higher growth rates, consistent with a trade-off under heat stress. Our results demonstrate that global populations of P. nodorum possess significant evolutionary potential to adapt to changing local conditions, including warmer temperatures and applications of fungicides.

Parasexual and Sexual Reproduction in Arbuscular Mycorrhizal Fungi: Room for Both

Arbuscular mycorrhizal fungi (AMF) harbor thousands of nuclei in a large syncytium at all times. Although mating processes have not been observed in AMF, their cells and genomes show many signatures of sexual reproduction. Here, we describe how some of these signatures could also arise from parasexual processes in these widespread plant symbionts. As such, parasexual and sexual evolution could both be at play in generating nuclear diversity in AMF.

Genomic analysis of a riboflavin-overproducing Ashbya gossypii mutant isolated by disparity mutagenesis

BACKGROUND

Ashbya gossypii naturally overproduces riboflavin and has been utilized for industrial riboflavin production. To improve riboflavin production, various approaches have been developed. In this study, to investigate the change in metabolism of a riboflavin-overproducing mutant, namely, the W122032 strain (MT strain) that was isolated by disparity mutagenesis, genomic analysis was carried out.

RESULTS

In the genomic analysis, 33 homozygous and 1377 heterozygous mutations in the coding sequences of the genome of MT strain were detected. Among these heterozygous mutations, the proportion of mutated reads in each gene was different, ranging from 21 to 75%. These results suggest that the MT strain may contain multiple nuclei containing different mutations. We tried to isolate haploid spores from the MT strain to prove its ploidy, but this strain did not sporulate under the conditions tested. Heterozygous mutations detected in genes which are important for sporulation likely contribute to the sporulation deficiency of the MT strain. Homozygous and heterozygous mutations were found in genes encoding enzymes involved in amino acid metabolism, the TCA cycle, purine and pyrimidine nucleotide metabolism and the DNA mismatch repair system. One homozygous mutation in AgILV2 gene encoding acetohydroxyacid synthase, which is also a flavoprotein in mitochondria, was found. Gene ontology (GO) enrichment analysis showed heterozygous mutations in all 22 DNA helicase genes and genes involved in oxidation-reduction process.

CONCLUSION

This study suggests that oxidative stress and the aging of cells were involved in the riboflavin over-production in A. gossypii riboflavin over-producing mutant and provides new insights into riboflavin production in A. gossypii and the usefulness of disparity mutagenesis for the creation of new types of mutants for metabolic engineering.

Guide RNA Engineering Enables Dual Purpose CRISPR-Cpf1 for Simultaneous Gene Editing and Gene Regulation in Yarrowia lipolytica

Yarrowia lipolytica has fast become a biotechnologically significant yeast for its ability to accumulate lipids to high levels. While there exists a suite of synthetic biology tools for genetic engineering in this yeast, there is a need for multipurposed tools for rapid strain generation. Here, we describe a dual purpose CRISPR-Cpf1 system that is capable of simultaneous gene disruption and gene regulation. Truncating guide RNA spacer length to 16 nt inhibited nuclease activity but not binding to the target loci, enabling gene activation and repression with Cpf1-fused transcriptional regulators. Gene repression was demonstrated using a Cpf1-Mxi1 fusion achieving a 7-fold reduction in mRNA, while CRISPR-activation with Cpf1-VPR increased hrGFP expression by 10-fold. High efficiency disruptions were achieved with gRNAs 23-25 bp in length, and efficiency and repression levels were maintained with multiplexed expression of truncated and full-length gRNAs. The developed CRISPR-Cpf1 system should prove useful in metabolic engineering, genome wide screening, and functional genomics studies.

Genome Resource for Two Stemphylium vesicarium Isolates Causing Stemphylium Leaf Blight of Onion in New York

Stemphylium leaf blight caused by Stemphylium vesicarium was recently identified as an emerging disease and dominant in the foliar disease complex affecting onion in New York. Here, we report the genomes of two isolates of S. vesicarium, On16-63 and On16-391. The availability of the genomes will accelerate genomic studies of S. vesicarium, including population biology, sexual reproduction, and fungicide resistance. Additionally, comparative genomics with the other published genome of S. vesicarium causing brown spot of pear will help understand pathogen biology and underpin the development of management strategies for this disease.

Genome Sequence of Fusarium oxysporum f. sp. matthiolae, a Brassicaceae Pathogen

The filamentous fungus Fusarium oxysporum is a soilborne pathogen of many cultivated species and an opportunistic pathogen of humans. F. oxysporum f. sp. matthiolae is one of three formae speciales that are pathogenic to crucifers, including Arabidopsis thaliana, a premier model for plant molecular biology and genetics. Here, we report a genome assembly of F. oxysporum f. sp. matthiolae strain PHW726, generated using a combination of PacBio and Illumina sequencing technologies. The genome assembly presented here should facilitate in-depth investigation of F. oxysporum-Arabidopsis interactions and shed light on the genetics of fungal pathogenesis and plant immunity.

Genome Sequence Resource for Elsinoë ampelina, the Causal Organism of Grapevine Anthracnose

Elsinoë ampelina is an ascomycetous fungus that causes grape anthracnose, a potentially devastating disease worldwide. Here, we report a 28.29 Mb high-quality genome sequence of E. ampelina YL-1 that encodes 8,057 predicted protein-coding genes and represents the first sequenced genome assembly of E. ampelina. This study adds to the current genomic resources for the genus Elsinoë and paves the way for research on comparative genomic studies, E. ampelina-grape interactions, and improvement of management strategies.

Genome Sequences of Verticillium dahliae Defoliating Strain XJ592 and Nondefoliating Strain XJ511

Verticillium dahliae is a widely distributed soilborne pathogen that causes vascular wilt in more than 200 plant species. Defoliating and nondefoliating symptoms caused by the disease that result in either the loss or retention of leaves in infected plants, respectively, in hosts such as cotton, olive, and okra, divide the causal agent into defoliating and nondefoliating pathotypes. Our goal in this current work was to generate genome resources for the defoliating strain XJ592 and the nondefoliating strain XJ511 of V. dahliae isolated from cotton in China.

Fungal secretome profile categorization of CAZymes by function and family corresponds to fungal phylogeny and taxonomy: Example Aspergillus and Penicillium

Fungi secrete an array of carbohydrate-active enzymes (CAZymes), reflecting their specialized habitat-related substrate utilization. Despite its importance for fitness, enzyme secretome composition is not used in fungal classification, since an overarching relationship between CAZyme profiles and fungal phylogeny/taxonomy has not been established. For 465 Ascomycota and Basidiomycota genomes, we predicted CAZyme-secretomes, using a new peptide-based annotation method, Conserved-Unique-Peptide-Patterns, enabling functional prediction directly from sequence. We categorized each enzyme according to CAZy-family and predicted molecular function, hereby obtaining a list of "EC-Function;CAZy-Family" observations. These "Function;Family"-based secretome profiles were compared, using a Yule-dissimilarity scoring algorithm, giving equal consideration to the presence and absence of individual observations. Assessment of "Function;Family" enzyme profile relatedness (EPR) across 465 genomes partitioned Ascomycota from Basidiomycota placing Aspergillus and Penicillium among the Ascomycota. Analogously, we calculated CAZyme "Function;Family" profile-similarities among 95 Aspergillus and Penicillium species to form an alignment-free, EPR-based dendrogram. This revealed a stunning congruence between EPR categorization and phylogenetic/taxonomic grouping of the Aspergilli and Penicillia. Our analysis suggests EPR grouping of fungi to be defined both by "shared presence" and "shared absence" of CAZyme "Function;Family" observations. This finding indicates that CAZymes-secretome evolution is an integral part of fungal speciation, supporting integration of cladogenesis and anagenesis.

Genome sequence of the fungus Pycnoporus sanguineus, which produces cinnabarinic acid and pH- and thermo- stable laccases

Pycnoporus sanguineus, an edible mushroom, produces antimicrobial and antitumor bioactive compounds and pH- and thermo- stable laccases that have multiple potential biotechnological applications. Here we reported the complete genome of the species Pycnoporus sanguineus ACCC 51,180 by using the combination of Illumina HiSeq X Ten and the PacBio sequencing technology. The represented genome is 36.6 Mb composed of 59 scaffolds with 12,086 functionally annotated protein-coding genes. The genome of Pycnoporus sanguineus encodes at least 19 biosynthetic gene clusters for secondary metabolites, including a terpene cluster for biosynthesis of the antitumor clavaric acid. Seven laccases were identified, while 22 genes were found to be involved in the kynurenine pathway in which the intermediate metabolite 3-hydroxyanthranilic acid were catalyzed by laccases into cinnabarinic acid. This study represented the third genome of the genus Pycnoporus, and wound facilitate the exploration of useful sources from Pycnoporus sanguineus for future industrial applications.

Impact of Chromosome Fusions on 3D Genome Organization and Gene Expression in Budding Yeast

The three-dimensional (3D) organization of chromosomes can influence transcription. However, the frequency and magnitude of these effects remain debated. To determine how changes in chromosome positioning affect transcription across thousands of genes with minimal perturbation, we characterized nuclear organization and global gene expression in budding yeast containing chromosome fusions. We used computational modeling and single-cell imaging to determine chromosome positions, and integrated these data with genome-wide transcriptional profiles from RNA sequencing. We find that chromosome fusions dramatically alter 3D nuclear organization without leading to strong genome-wide changes in transcription. However, we observe a mild but significant and reproducible increase in the expression of genes displaced away from the periphery. The increase in transcription is inversely proportional to the propensity of a given locus to be at the nuclear periphery; for example, a 10% decrease in the propensity of a gene to reside at the nuclear envelope is accompanied by a 10% increase in gene expression. Modeling suggests that this is due to both deletion of telomeres and to displacement of genes relative to the nuclear periphery. These data suggest that basal transcriptional activity is sensitive to radial changes in gene position, and provide insight into the functional relevance of budding yeast chromosome-level 3D organization in gene expression.

Homologous recombination repair intermediates promote efficient de novo telomere addition at DNA double-strand breaks

The healing of broken chromosomes by de novo telomere addition, while a normal developmental process in some organisms, has the potential to cause extensive loss of heterozygosity, genetic disease, or cell death. However, it is unclear how de novo telomere addition (dnTA) is regulated at DNA double-strand breaks (DSBs). Here, using a non-essential minichromosome in fission yeast, we identify roles for the HR factors Rqh1 helicase, in concert with Rad55, in suppressing dnTA at or near a DSB. We find the frequency of dnTA in rqh1Δ rad55Δ cells is reduced following loss of Exo1, Swi5 or Rad51. Strikingly, in the absence of the distal homologous chromosome arm dnTA is further increased, with nearly half of the breaks being healed in rqh1Δ rad55Δ or rqh1Δ exo1Δ cells. These findings provide new insights into the genetic context of highly efficient dnTA within HR intermediates, and how such events are normally suppressed to maintain genome stability.

Efficient targeted mutation of genomic essential genes in yeast Saccharomyces cerevisiae

Targeted gene mutation by allelic replacement is important for functional genomic analysis and metabolic engineering. However, it is challenging in mutating the essential genes with the traditional method by using a selection marker, since the first step of essential gene knockout will result in a lethal phenotype. Here, we developed a two-end selection marker (Two-ESM) method for site-directed mutation of essential genes in Saccharomyces cerevisiae with the aid of the CRISPR/Cas9 system. With this method, single and double mutations of the essential gene ERG20 (encoding farnesyl diphosphate synthase) in S. cerevisiae were successfully constructed with high efficiencies of 100%. In addition, the Two-ESM method significantly improved the mutation efficiency and simplified the genetic manipulation procedure compared with traditional methods. The genome integration and mutation efficiencies were further improved by dynamic regulation of mutant gene expression and optimization of the integration modules. This Two-ESM method will facilitate the construction of genomic mutations of essential genes for functional genomic analysis and metabolic flux regulation in yeasts. KEY POINTS: • A Two-ESM strategy achieves mutations of essential genes with high efficiency of 100%. • The optimized three-module method improves the integration efficiency by more than three times. • This method will facilitate the functional genomic analysis and metabolic flux regulation.

Comparative genome-scale analysis of Pichia pastoris variants informs selection of an optimal base strain

Komagataella phaffii, also known as Pichia pastoris, is a common host for the production of biologics and enzymes, due to fast growth, high productivity, and advancements in host engineering. Several K. phaffii variants are commonly used as interchangeable base strains, which confounds efforts to improve this host. In this study, genomic and transcriptomic analyses of Y-11430 (CBS7435), GS115, X-33, and eight other variants enabled a comparative assessment of the relative fitness of these hosts for recombinant protein expression. Cell wall integrity explained the majority of the variation among strains, impacting transformation efficiency, growth, methanol metabolism, and secretion of heterologous proteins. Y-11430 exhibited the highest activity of genes involved in methanol utilization, up to two-fold higher transcription of heterologous genes, and robust growth. With a more permeable cell wall, X-33 displayed a six-fold higher transformation efficiency and up to 1.2-fold higher titers than Y-11430. X-33 also shared nearly all mutations, and a defective variant of HIS4, with GS115, precluding robust growth. Transferring two beneficial mutations identified in X-33 into Y-11430 resulted in an optimized base strain that provided up to four-fold higher transformation efficiency and three-fold higher protein titers, while retaining robust growth. The approach employed here to assess unique banked variants in a species and then transfer key beneficial variants into a base strain should also facilitate rational assessment of a broad set of other recombinant hosts.

Draft Genome Sequences of Elsinoë fawcettii and Elsinoë australis Causing Scab Diseases on Citrus

Elsinoë fawcettii and E. australis (phylum Ascomycota) are phytopathogenic fungi causing scab diseases on citrus plants. We report here the high-quality draft genome sequences and ab initio gene predictions of two E. fawcettii strains and one E. australis strain, which differ in their host range. This genome sequence information will provide valuable resources to underpin genomic attributes for determining host range through comparative genomic analyses of citrus scab fungi.

Draft Genomic Resources for the Brown Rot Fungal Pathogen Monilinia laxa

Monilinia laxa is the causal agent of brown rot on stone fruit, and it can cause heavy yield losses during field production and postharvest storage. This article reports the draft genome assembly of the M. laxa Mlax316 strain, obtained using a hybrid genome assembly with both Illumina short-reads and PacBio long-reads sequencing technologies. The complete draft genome consists of 49 scaffolds with total size of 42.81 Mb, and scaffold N₅₀ of 2,449.4 kb. Annotation of the M. laxa assembly identified 11,163 genes and 12,424 proteins which were functionally annotated. This new genome draft improves current genomic resources available for M. laxa and represents a useful tool for further research into its interactions with host plants and into evolution in the Monilinia genus.

The Genome Sequence of Five Genotypes of Fusarium oxysporum f. sp. vasinfectum: A Resource for Studies on Fusarium Wilt of Cotton

Fusarium oxysporum f. sp. vasinfectum is an important plant pathogen responsible for vascular wilt disease on cotton. Members of this group are known to carry supernumerary chromosomes that encode virulence factors. We sequenced the genomes of five F. oxysporum f. sp. vasinfectum isolates, including the genome of a representative of the highly virulent genotype race 4, at a high coverage to assemble reference-quality genomes. These genomes provide a necessary resource for comparative genomic analyses to identify genes or genome features that are involved in pathogenicity on cotton and may ultimately be used to identify improved management strategies.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.

Colletotrichum: species complexes, lifestyle, and peculiarities of some sources of genetic variability

The genus Colletotrichum comprises species with different lifestyles but is mainly known for phytopathogenic species that infect crops of agronomic relevance causing considerable losses. The fungi of the genus Colletotrichum are distributed in species complexes and within each complex some species have particularities regarding their lifestyle. The most commonly found and described lifestyles in Colletotrichum are endophytic and hemibiotrophic phytopathogenic. Several of these phytopathogenic species show wide genetic variability, which makes long-term maintenance of resistance in plants difficult. Different mechanisms may play an important role in the emergence of genetic variants but are not yet fully understood in this genus. These mechanisms include heterokaryosis, a parasexual cycle, sexual cycle, transposable element activity, and repeat-induced point mutations. This review provides an overview of the genus Colletotrichum, the species complexes described so far and the most common lifestyles in the genus, with a special emphasis on the mechanisms that may be responsible, at least in part, for the emergence of new genotypes under field conditions.

SCRaMbLE-in: A Fast and Efficient Method to Diversify and Improve the Yields of Heterologous Pathways in Synthetic Yeast

The synthetic chromosome rearrangement and modification by LoxP-mediated evolution (SCRaMbLE) system is a key component of the synthetic yeast genome (Sc2.0) project, an international effort to construct an entire synthetic genome in yeast. SCRaMbLE involves the introduction of thousands of symmetrical LoxP (LoxPsym) recombination sites downstream of every nonessential gene in all 16 chromosomes, enabling numerous genome rearrangements in the form of deletions, inversions, duplications, and translocations by the Cre-LoxPsym recombination system. We highlight a two-step protocol for SCRaMbLE-in (Liu, Nat Commun 9(1):1936, 2018), a recombinase-based combinatorial method to expedite genetic engineering and exogenous pathway optimization, using a synthetic β-carotene pathway as an example. First, an in vitro phase uses a recombinase toolkit to diversify gene expression by integrating various regulatory elements into the target pathway. This combinatorial pathway library can be transformed directly into yeast for traditional screening. Once an optimized pathway which is flanked by LoxPsym sites is identified, it is transformed into Sc2.0 yeast for the in vivo SCRaMbLE phase, where LoxPsym sites in the synthetic yeast genome and Cre recombinase catalyze massive genome rearrangements. We describe all the conditions necessary to perform SCRaMbLE and post-SCRaMbLE experiments including screening, spot test analysis, and PCRTag analysis to elucidate genotype-phenotype relationships.

The Fml1-MHF complex suppresses inter-fork strand annealing in fission yeast

Previously we reported that a process called inter-fork strand annealing (IFSA) causes genomic deletions during the termination of DNA replication when an active replication fork converges on a collapsed fork (Morrow et al., 2017). We also identified the FANCM-related DNA helicase Fml1 as a potential suppressor of IFSA. Here, we confirm that Fml1 does indeed suppress IFSA, and show that this function depends on its catalytic activity and ability to interact with Mhf1-Mhf2 via its C-terminal domain. Finally, a plausible mechanism of IFSA suppression is demonstrated by the finding that Fml1 can catalyse regressed fork restoration in vitro.

Origin of Lactose Fermentation in Kluyveromyces lactis by Interspecies Transfer of a Neo-functionalized Gene Cluster during Domestication

Humans have used yeasts to make cheese and kefir for millennia, but the ability to ferment the milk sugar lactose is found in only a few yeast species, of which the foremost is Kluyveromyces lactis [1]. Two genes, LAC12 (lactose permease) and LAC4 (lactase), are sufficient for lactose uptake and hydrolysis to glucose and galactose [2]. Here, we show that these genes have a complex evolutionary history in the genus Kluyveromyces that is likely the result of human activity during domestication. We show that the ancestral Lac12 was bifunctional, able to import both lactose and cellobiose into the cell. These disaccharides were then hydrolyzed by Lac4 in the case of lactose or Cel2 in the case of cellobiose. A second cellobiose transporter, Cel1, was also present ancestrally. In the K. lactis lineage, the ancestral LAC12 and LAC4 were lost and a separate upheaval in the sister species K. marxianus resulted in loss of CEL1 and quadruplication of LAC12. One of these LAC12 genes became neofunctionalized to encode an efficient lactose transporter capable of supporting fermentation, specifically in dairy strains of K. marxianus, where it formed a LAC4-LAC12-CEL2 gene cluster, although another remained a cellobiose transporter. Then, the ability to ferment lactose was acquired very recently by K. lactis var. lactis by introgression of LAC12 and LAC4 on a 15-kb subtelomeric region from a dairy strain of K. marxianus. The genomic history of the LAC genes shows that strong selective pressures were imposed on yeasts by early dairy farmers.

The Zymoseptoria tritici ORFeome: A Functional Genomics Community Resource

Libraries of protein-encoding sequences can be generated by identification of open reading frames (ORFs) from a genome of choice that are then assembled into collections of plasmids termed ORFeome libraries. These represent powerful resources to facilitate functional genomic characterization of genes and their encoded products. Here, we report the generation of an ORFeome for Zymoseptoria tritici, which causes the most serious disease of wheat in temperate regions of the world. We screened the genome of strain IP0323 for high confidence gene models, identifying 4,075 candidates from 10,933 predicted genes. These were amplified from genomic DNA, were cloned into the Gateway entry vector pDONR207, and were sequenced, providing a total of 3,022 quality-controlled plasmids. The ORFeome includes genes predicted to encode effectors (n = 410) and secondary metabolite biosynthetic proteins (n = 171) in addition to genes residing at dispensable chromosomes (n = 122) or those that are preferentially expressed during plant infection (n = 527). The ORFeome plasmid library is compatible with our previously developed suite of Gateway destination vectors, which have various combinations of promoters, selection markers, and epitope tags. The Z. tritici ORFeome constitutes a powerful resource for functional genomics and offers unparalleled opportunities to understand the biology of Z. tritici.[Formula: see text] Copyright © 2019 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.