First genome-wide analysis of the endangered, endemic lichen Cetradonia linearis reveals isolation by distance and strong population structure

PREMISE OF THE STUDY

Lichenized fungi are evolutionarily diverse and ecologically important, but little is known about the processes that drive their diversification and genetic differentiation. Distributions are often assumed to be wholly shaped by ecological requirements rather than dispersal limitations. Furthermore, although asexual and sexual reproductive structures are observable, the lack of information about recombination rates makes inferences about reproductive strategies difficult. We investigated the population genomics of Cetradonia linearis, a federally endangered lichen in the southern Appalachians of eastern North America, to test the relative contributions of environmental and geographic distance in shaping genetic structure, and to characterize the mating system and genome-wide recombination.

METHODS

Whole-genome shotgun sequencing was conducted to generate data for 32 individuals of C. linearis. A reference genome was assembled, and reads from all samples were aligned to generate a set of single-nucleotide polymorphisms for further analyses.

KEY RESULTS

We found evidence for low rates of recombination and for isolation by distance, but not for isolation by environment. The species is putatively unisexual, given that only one mating-type locus was found. Hindcast species distribution models and the distribution of genetic diversity support C. linearis having a larger range during the Last Glacial Maximum in the southern portion of its current extent.

CONCLUSIONS

Our findings contribute to the understanding of factors that shape genetic diversity in C. linearis and in fungi more broadly. Because all populations are highly genetically differentiated, the extirpation of any population would mean the loss of unique genetic diversity; therefore, our results support the continued conservation of this species.

Genome wide analysis of meiotic recombination in yeast: For a few SNPs more

Diploid organisms undergo meiosis to produce haploid germ cells. Crossover events during meiosis promote genetic diversity and facilitate accurate chromosome segregation. The baker's yeast Saccharomyces cerevisiae is extensively used as a model for analysis of meiotic recombination. Conventional methods for measuring recombination events in S. cerevisiae have been limited by the number and density of genetic markers. Next generation sequencing (NGS)-based analysis of hybrid yeast genomes bearing thousands of heterozygous single nucleotide polymorphism (SNP) markers has revolutionized analysis of meiotic recombination. By facilitating analysis of marker segregation in the whole genome with unprecedented resolution, this method has resulted in the generation of high-resolution recombination maps in wild-type and meiotic mutants. These studies have provided novel insights into the mechanism of meiotic recombination. In this review, we discuss the methodology, challenges, insights and future prospects of using NGS-based methods for whole genome analysis of meiotic recombination. The objective is to facilitate the use of these high through-put sequencing methods for the analysis of meiotic recombination given their power to provide significant new insights into the process. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology, 70(8):743-752, 2018.

Accessories Make the Outfit: Accessory Chromosomes and Other Dispensable DNA Regions in Plant-Pathogenic Fungi

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

Annotation resource of tandem repeat-containing secretory proteins in sixty fungi

Fungal secretory proteins that interact with host plants are regarded as effectors. Because fungal effectors rarely contain conserved sequence features, identification and annotation of fungal effectors from predicted secretory proteins are difficult using outward comparison methods such as BLAST or hidden Markov model. In desire of more sequence features to prioritize research interests of fungal secretory proteins, this study developed a pipeline to identify tandem repeat (TR) domain within putative secretory proteins and tested a hypothesis that at least one type of TR domain in non-orthologous secretory proteins has emerged from convergent evolution for plant pathogenicity. There were 2804 types of TR domains and a total of 2925 TR-containing secretory proteins found from 60 fungi. There was no conserved type of TR domain shared only by plant pathogens, indicating functional divergence for different types of TR domain and TR-containing secretory proteins. The annotation resource of putative fungal TR-containing secretory proteins provides new sequence features that will be useful for the community interested in fungal effector biology.

Genome-Wide Characterization and Expression Analyses of Pleurotus ostreatus MYB Transcription Factors during Developmental Stages and under Heat Stress Based on de novo Sequenced Genome

Pleurotus ostreatus is a commercially grown mushroom species in China. However, studies on the mechanisms of the fruiting body development and stress response of P. ostreatus are still at a primary stage. In this study, we report the entire genome sequence of P. ostreatus CCMSSC03989. Then, we performed comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during a series of developmental stages and under the condition of heat stress. A 34.76 Mb genome was obtained through next-generation sequencing (NGS) and Bionano optical mapping approaches. The genome has a scaffold N50 of 1.1 Mb and contains 10.11% repeats, and 10,936 gene models were predicted. A total of 20 MYB genes (PoMYB) were identified across the genome, and the full-length open reading frames were isolated. The PoMYBs were classified into 1 repeat (1R), 2R, and 3R-MYB groups according to their MYB domain repeat numbers, and 3R-MYBs possessed relatively more introns than 1R and 2R-MYBs. Based on phylogenetic analysis, the PoMYBs were divided into four groups and showed close relationships with the MYB genes of plants and fungi. RNA-sequencing (RNA-Seq) and quantitative PCR (qPCR) analyses revealed that PoMYB expression showed stage-specific patterns in reproductive stages and could be induced by heat stress. The P. ostreatus draft genome will promote genome-wide analysis, and our study of PoMYBs will promote further functional analysis of MYB genes in mushrooms.

The origin and adaptive evolution of domesticated populations of yeast from Far East Asia

The yeast Saccharomyces cerevisiae has been an essential component of human civilization because of its long global history of use in food and beverage fermentation. However, the diversity and evolutionary history of the domesticated populations of the yeast remain elusive. We show here that China/Far East Asia is likely the center of origin of the domesticated populations of the species. The domesticated populations form two major groups associated with solid- and liquid-state fermentation and appear to have originated from heterozygous ancestors, which were likely formed by outcrossing between diverse wild isolates primitively for adaptation to maltose-rich niches. We found consistent gene expansion and contraction in the whole domesticated population, as well as lineage-specific genome variations leading to adaptation to different environments. We show a nearly panoramic view of the diversity and life history of S. cerevisiae and provide new insights into the origin and evolution of the species.

Whole genome comparison of Aspergillus flavus L-morphotype strain NRRL 3357 (type) and S-morphotype strain AF70

Aspergillus flavus is a saprophytic fungus that infects corn, peanuts, tree nuts and other agriculturally important crops. Once the crop is infected the fungus has the potential to secrete one or more mycotoxins, the most carcinogenic of which is aflatoxin. Aflatoxin contaminated crops are deemed unfit for human or animal consumption, which results in both food and economic losses. Within A. flavus, two morphotypes exist: the S strains (small sclerotia) and L strains (large sclerotia). Significant morphological and physiological differences exist between the two morphotypes. For example, the S-morphotypes produces sclerotia that are smaller (< 400 μm), greater in quantity, and contain higher concentrations of aflatoxin than the L-morphotypes (>400 μm). The morphotypes also differ in pigmentation, pH homeostasis in culture and the number of spores produced. Here we report the first full genome sequence of an A. flavus S morphotype, strain AF70. We provide a comprehensive comparison of the A. flavus S-morphotype genome sequence with a previously sequenced genome of an L-morphotype strain (NRRL 3357), including an in-depth analysis of secondary metabolic clusters and the identification SNPs within their aflatoxin gene clusters.

dbCAN2: a meta server for automated carbohydrate-active enzyme annotation

Complex carbohydrates of plants are the main food sources of animals and microbes, and serve as promising renewable feedstock for biofuel and biomaterial production. Carbohydrate active enzymes (CAZymes) are the most important enzymes for complex carbohydrate metabolism. With an increasing number of plant and plant-associated microbial genomes and metagenomes being sequenced, there is an urgent need of automatic tools for genomic data mining of CAZymes. We developed the dbCAN web server in 2012 to provide a public service for automated CAZyme annotation for newly sequenced genomes. Here, dbCAN2 (http://cys.bios.niu.edu/dbCAN2) is presented as an updated meta server, which integrates three state-of-the-art tools for CAZome (all CAZymes of a genome) annotation: (i) HMMER search against the dbCAN HMM (hidden Markov model) database; (ii) DIAMOND search against the CAZy pre-annotated CAZyme sequence database and (iii) Hotpep search against the conserved CAZyme short peptide database. Combining the three outputs and removing CAZymes found by only one tool can significantly improve the CAZome annotation accuracy. In addition, dbCAN2 now also accepts nucleotide sequence submission, and offers the service to predict physically linked CAZyme gene clusters (CGCs), which will be a very useful online tool for identifying putative polysaccharide utilization loci (PULs) in microbial genomes or metagenomes.

Pathogenic adaptations of Colletotrichum fungi revealed by genome wide gene family evolutionary analyses

The fungal genus Colletotrichum contains hemibiotrophic phytopathogens being highly variable in host and tissue specificities. We sequenced a C. fructicola genome (1104–7) derived from an isolate of apple in China and compared it with the reference genome (Nara_gc5) derived from an isolate of strawberry in Japan. Mauve alignment and BlastN search identified 0.62 Mb lineage-specific (LS) genomic regions in 1104–7 with a length criterion of 10 kb. Genes located within LS regions evolved more dynamically, and a strongly elevated proportion of genes were closely related to non-Colletotrichum sequences. Two LS regions, containing nine genes in total, showed features of fungus-to-fungus horizontal transfer supported by both gene order collinearity and gene phylogeny patterns. We further compared the gene content variations among 13 Colletotrichum and 11 non-Colletotrichum genomes by gene function annotation, OrthoMCL grouping and CAFE analysis. The results provided a global evolutionary picture of Colletotrichum gene families, and identified a number of strong duplication/loss events at key phylogenetic nodes, such as the contraction of the detoxification-related RTA1 family in the monocot-specializing graminicola complex and the expansions of several ammonia production-related families in the fruit-infecting gloeosporioides complex. We have also identified the acquirement of a RbsD/FucU fucose transporter from bacterium by the Colletotrichum ancestor. In sum, this study summarized the pathogenic evolutionary features of Colletotrichum fungi at multiple taxonomic levels and highlights the concept that the pathogenic successes of Colletotrichum fungi require shared as well as lineage-specific virulence factors.

Complete genomic and transcriptional landscape analysis using third-generation sequencing: a case study of Saccharomyces cerevisiae CEN.PK113-7D

Completion of eukaryal genomes can be difficult task with the highly repetitive sequences along the chromosomes and short read lengths of second-generation sequencing. Saccharomyces cerevisiae strain CEN.PK113-7D, widely used as a model organism and a cell factory, was selected for this study to demonstrate the superior capability of very long sequence reads for de novo genome assembly. We generated long reads using two common third-generation sequencing technologies (Oxford Nanopore Technology (ONT) and Pacific Biosciences (PacBio)) and used short reads obtained using Illumina sequencing for error correction. Assembly of the reads derived from all three technologies resulted in complete sequences for all 16 yeast chromosomes, as well as the mitochondrial chromosome, in one step. Further, we identified three types of DNA methylation (5mC, 4mC and 6mA). Comparison between the reference strain S288C and strain CEN.PK113-7D identified chromosomal rearrangements against a background of similar gene content between the two strains. We identified full-length transcripts through ONT direct RNA sequencing technology. This allows for the identification of transcriptional landscapes, including untranslated regions (UTRs) (5' UTR and 3' UTR) as well as differential gene expression quantification. About 91% of the predicted transcripts could be consistently detected across biological replicates grown either on glucose or ethanol. Direct RNA sequencing identified many polyadenylated non-coding RNAs, rRNAs, telomere-RNA, long non-coding RNA and antisense RNA. This work demonstrates a strategy to obtain complete genome sequences and transcriptional landscapes that can be applied to other eukaryal organisms.

The yeast Ty1 retrotransposon requires components of the nuclear pore complex for transcription and genomic integration

Nuclear pore complexes (NPCs) orchestrate cargo between the cytoplasm and nucleus and regulate chromatin organization. NPC proteins, or nucleoporins (Nups), are required for human immunodeficiency virus type 1 (HIV-1) gene expression and genomic integration of viral DNA. We utilize the Ty1 retrotransposon of Saccharomyces cerevisiae (S. cerevisiae) to study retroviral integration because retrotransposons are the progenitors of retroviruses and have conserved integrase (IN) enzymes. Ty1-IN targets Ty1 elements into the genome upstream of RNA polymerase (Pol) III transcribed genes such as transfer RNA (tRNA) genes. Evidence that S. cerevisiae tRNA genes are recruited to NPCs prompted our investigation of a functional role for the NPC in Ty1 targeting into the genome. We find that Ty1 mobility is reduced in multiple Nup mutants that cannot be accounted for by defects in Ty1 gene expression, cDNA production or Ty1-IN nuclear entry. Instead, we find that Ty1 insertion upstream of tRNA genes is impaired. We also identify Nup mutants with wild type Ty1 mobility but impaired Ty1 targeting. The NPC nuclear basket, which interacts with chromatin, is required for both Ty1 expression and nucleosome targeting. Deletion of components of the NPC nuclear basket causes mis-targeting of Ty1 elements to the ends of chromosomes.

Control of yeast retrotransposons mediated through nucleoporin evolution

Yeasts serve as hosts to several types of genetic parasites. Few studies have addressed the evolutionary trajectory of yeast genes that control the stable co-existence of these parasites with their host cell. In Saccharomyces yeasts, the retrovirus-like Ty retrotransposons must access the nucleus. We show that several genes encoding components of the yeast nuclear pore complex have experienced natural selection for substitutions that change the encoded protein sequence. By replacing these S. cerevisiae genes with orthologs from other Saccharomyces species, we discovered that natural sequence changes have affected the mobility of Ty retrotransposons. Specifically, changing the genetic sequence of NUP84 or NUP82 to match that of other Saccharomyces species alters the mobility of S. cerevisiae Ty1 and Ty3. Importantly, all tested housekeeping functions of NUP84 and NUP82 remained equivalent across species. Signatures of natural selection, resulting in altered interactions with viruses and parasitic genetic elements, are common in host defense proteins. Yet, few instances have been documented in essential housekeeping proteins. The nuclear pore complex is the gatekeeper of the nucleus. This study shows how the evolution of this large, ubiquitous eukaryotic complex can alter the replication of a molecular parasite, but concurrently maintain essential host functionalities regarding nucleocytoplasmic trafficking.

Molecular Characterization and Functional Analysis of PR-1-Like Proteins Identified from the Wheat Head Blight Fungus Fusarium graminearum

The group 1 pathogenesis-related (PR-1) proteins originally identified from plants and their homologs are also found in other eukaryotic kingdoms. Studies on nonplant PR-1-like (PR-1L) proteins have been pursued widely in humans and animals but rarely in filamentous ascomycetes. Here, we report the characterization of four PR-1L proteins identified from the ascomycete fungus Fusarium graminearum, the primary cause of Fusarium head blight of wheat and barley (designated FgPR-1L). Molecular cloning revealed that the four FgPR-1L proteins are all encoded by small open reading frames (612 to 909 bp) that are often interrupted by introns, in contrast to plant PR-1 genes that lack introns. Sequence analysis indicated that all FgPR-1L proteins contain the PR-1-specific three-dimensional structure, and one of them features a C-terminal transmembrane (TM) domain that has not been reported for any stand-alone PR-1 proteins. Transcriptional analysis revealed that the four FgPR-1L genes are expressed in axenic cultures and in planta with different spatial or temporal expression patterns. Phylogenetic analysis indicated that fungal PR-1L proteins fall into three major groups, one of which harbors FgPR-1L-2-related TM-containing proteins from both phytopathogenic and human-pathogenic ascomycetes. Low-temperature sodium dodecyl sulfate polyacrylamide gel electrophoresis and proteolytic assays indicated that the recombinant FgPR-1L-4 protein exists as a monomer and is resistant to subtilisin of the serine protease family. Functional analysis confirmed that deletion of the FgPR-1L-4 gene from the fungal genome results in significantly reduced virulence on susceptible wheat. This study provides the first example that the F. graminearum-wheat interaction involves a pathogen-derived PR-1L protein that affects fungal virulence on the host.

Comparative genomic analysis of fungal TPP-riboswitches

Riboswitches are conserved RNA structures located in non-coding regions of mRNA and able to bind small molecules (e.g. metabolites) changing conformation upon binding. This feature enables them to function as regulators of gene expression. The thiamin pyrophosphate (TPP) riboswitch is the only type of riboswitches found not only in bacteria, but also in eukaryotes - in plants, green algae, protists, and fungi. Two main mechanisms of fungal TPP riboswitch action, involving alternative splicing, have been established so far. Here, we report a large-scale bioinformatic study of riboswitch structural features, action mechanisms, and distribution along the fungal taxonomy groups. For each putatively regulated gene, we reconstruct the riboswitch structure, identify other components of the regulation machinery, and establish mechanisms of riboswitch-mediated regulation. In addition to three genes known to be regulated by TPP riboswitches, thiazole synthase THI4, hydroxymethilpyrimidine-syntase NMT1, and putative transporter NCU01977, we identify two new genes, a putative thiamin transporter THI9 and a transporter of unknown specificity. While the riboswitch sequence and structure remain highly conserved in all species and genes, the mode of riboswitch-mediated regulation varies between regulated genes. The riboswitch usage varies strongly between fungal taxa, with the largest number of riboswitch-regulated genes found in Pezizomycotina and no riboswitch-mediated regulation established in Saccaromycotina.

10KP: A phylodiverse genome sequencing plan

Understanding plant evolution and diversity in a phylogenomic context is an enormous challenge due, in part, to limited availability of genome-scale data across phylodiverse species. The 10KP (10,000 Plants) Genome Sequencing Project will sequence and characterize representative genomes from every major clade of embryophytes, green algae, and protists (excluding fungi) within the next 5 years. By implementing and continuously improving leading-edge sequencing technologies and bioinformatics tools, 10KP will catalogue the genome content of plant and protist diversity and make these data freely available as an enduring foundation for future scientific discoveries and applications. 10KP is structured as an international consortium, open to the global community, including botanical gardens, plant research institutes, universities, and private industry. Our immediate goal is to establish a policy framework for this endeavor, the principles of which are outlined here.

The fungal snoRNAome

Small nucleolar RNAs (snoRNAs) are essential players in the rRNA biogenesis due to their involvement in the nucleolytic processing of the precursor and the subsequent guidance of nucleoside modifications. Within the kingdom Fungi, merely a few species-specific surveys have explored their snoRNA repertoire. However, the wide range of the snoRNA landscape spanning all major fungal lineages has not been mapped so far, mainly because of missing tools for automatized snoRNA detection and functional analysis. For the first time, we report here a comprehensive inventory of fungal snoRNAs together with a functional analysis and an in-depth investigation of their evolutionary history including innovations, deletions, and target switches. This large-scale analysis, incorporating more than 120 snoRNA families with more than 7700 individual snoRNA sequences, catalogs and clarifies the landscape of fungal snoRNA families, assigns functions to previously orphan snoRNAs, and increases the number of sequences by 450%. We also show that the snoRNAome is subject to ongoing rearrangements and adaptations, e.g., through lineage-specific targets and redundant guiding functions.

CDK activity provides temporal and quantitative cues for organizing genome duplication

In eukaryotes, the spatial and temporal organization of genome duplication gives rise to distinctive profiles of replication origin usage along the chromosomes. While it has become increasingly clear that these programs are important for cellular physiology, the mechanisms by which they are determined and modulated remain elusive. Replication initiation requires the function of cyclin-dependent kinases (CDKs), which associate with various cyclin partners to drive cell proliferation. Surprisingly, although we possess detailed knowledge of the CDK regulators and targets that are crucial for origin activation, little is known about whether CDKs play a critical role in establishing the genome-wide pattern of origin selection. We have addressed this question in the fission yeast, taking advantage of a simplified cell cycle network in which cell proliferation is driven by a single cyclin-CDK module. This system allows us to precisely control CDK activity in vivo using chemical genetics. First, in contrast to previous reports, our results clearly show that distinct cyclin-CDK pairs are not essential for regulating specific subsets of origins and for establishing a normal replication program. Importantly, we then demonstrate that the timing at which CDK activity reaches the S phase threshold is critical for the organization of replication in distinct efficiency domains, while the level of CDK activity at the onset of S phase is a dose-dependent modulator of overall origin efficiencies. Our study therefore implicates these different aspects of CDK regulation as versatile mechanisms for shaping the architecture of DNA replication across the genome.

The duplication of the genetic material is a highly conserved and tightly regulated process that is essential for cell proliferation. DNA synthesis initiates at sites called origins that are distributed throughout the genome. Replication origins are not all used equivalently, and their patterns of activation along the chromosomes give rise to specific profiles, or programs, of DNA replication. These programs change during development and in response to external stimuli, and we have previously shown that they have important consequences for cellular function. However, we still do not understand the mechanisms by which cells establish different replication patterns. Here we investigate the role of the cyclin-dependent kinase (CDK) family of proteins, whose activities are critical for cell cycle progression, in regulating the organization of genome duplication. Taking advantage of a system that allows us to precisely modulate CDK activity levels in living cells, we demonstrate that both the temporal and quantitative controls of CDK function are crucial for determining distinct programs of DNA replication. Our work therefore uncovers a fundamental link between CDK activity, a central input in a variety of cellular and developmental processes, and the architecture of genome duplication.

Reference Quality Genome Assemblies of Three Parastagonospora nodorum Isolates Differing in Virulence on Wheat

Parastagonospora nodorum, the causal agent of Septoria nodorum blotch in wheat, has emerged as a model necrotrophic fungal organism for the study of host-microbe interactions. To date, three necrotrophic effectors have been identified and characterized from this pathogen, including SnToxA, SnTox1, and SnTox3. Necrotrophic effector identification was greatly aided by the development of a draft genome of Australian isolate SN15 via Sanger sequencing, yet it remained largely fragmented. This research presents the development of nearly finished genomes of P. nodorum isolates Sn4, Sn2000, and Sn79-1087 using long-read sequencing technology. RNAseq analysis of isolate Sn4, consisting of eight time points covering various developmental and infection stages, mediated the annotation of 13,379 genes. Analysis of these genomes revealed large-scale polymorphism between the three isolates, including the complete absence of contig 23 from isolate Sn79-1087, and a region of genome expansion on contig 10 in isolates Sn4 and Sn2000. Additionally, these genomes exhibit the hallmark characteristics of a "two-speed" genome, being partitioned into two distinct GC-equilibrated and AT-rich compartments. Interestingly, isolate Sn79-1087 contains a lower proportion of AT-rich segments, indicating a potential lack of evolutionary hotspots. These newly sequenced genomes, consisting of telomere-to-telomere assemblies of nearly all 23 P. nodorum chromosomes, provide a robust foundation for the further examination of effector biology and genome evolution.

High-resolution mapping of the recombination landscape of the phytopathogen Fusarium graminearum suggests two-speed genome evolution

Recombination is a major evolutionary force, increasing genetic diversity and permitting efficient coevolution of fungal pathogen(s) with their host(s). The ascomycete Fusarium graminearum is a devastating pathogen of cereal crops, and can contaminate food and feed with harmful mycotoxins. Previous studies have suggested a high adaptive potential of this pathogen, illustrated by an increase in pathogenicity and resistance to fungicides. In this study, we provide the first detailed picture of the crossover events occurring during meiosis and discuss the role of recombination in pathogen evolution. An experimental recombinant population (n = 88) was created and genotyped using 1306 polymorphic markers obtained from restriction site-associated DNA sequencing (RAD-seq) and aligned to the reference genome. The construction of a high-density linkage map, anchoring 99% of the total length of the reference genome, allowed the identification of 1451 putative crossovers, positioned at a median resolution of 24 kb. The majority of crossovers (87.2%) occurred in a relatively small portion of the genome (30%). All chromosomes demonstrated recombination-active sections, which had a near 15-fold higher crossover rate than non-active recombinant sections. The recombination rate showed a strong positive correlation with nucleotide diversity, and recombination-active regions were enriched for genes with a putative role in host-pathogen interaction, as well as putative diversifying genes. Our results confirm the preliminary analysis observed in other F. graminearum strains and suggest a conserved 'two-speed' recombination landscape. The consequences with regard to the evolutionary potential of this major fungal pathogen are also discussed.

Evaluating Programs for Predicting Genes and Transcripts with RNA-Seq Support in Fungal Genomes

The steps needed to computationally predict genes and transcripts in fungal genomes with support from RNA-Seq data are described in detail for three prediction programs: CodingQuarry, BRAKER1, and Harfang. These programs predicted from 86% to 92% (Harfang) of the genes in a manually curated reference set for Aspergillus niger strain NRRL3. Genes with little or no RNA-Seq read coverage were predicted less successfully than genes with adequate coverage.

Efficient Extraction Method for High Quality Fungal RNA from Complex Lignocellulosic Substrates

Here we describe an efficient and reproducible method for the extraction of fungal RNA from complex lignocellulose containing materials. The fungal cells are snap-frozen and disrupted in chaotropic guanidinium thiocyanate buffer, after which the extracted RNA is isolated by using CsCl gradient ultracentrifugation. By lowering the pH of the extraction buffer, the procedure is also suitable for sample materials rich in humic acids. The method results in high quantity and quality RNA that is separated from endogenous contaminants (e.g., RNases) and substances derived from plant biomass (e.g., colored aromatic compounds). In addition, no further steps such as DNase treatment are needed. The extracted RNA is highly suitable for downstream gene expression analyses such as RNA sequencing.

CRISPR-Cas9-Mediated Genome Editing and Transcriptional Control in Yarrowia lipolytica

The discovery and adaptation of RNA-guided nucleases has resulted in the rapid development of efficient, scalable, and easily accessible synthetic biology tools for targeted genome editing and transcriptional control. In these systems, for example CRISPR-Cas9 from Streptococcus pyogenes, a protein with nuclease activity is targeted to a specific nucleotide sequence by a short RNA molecule, whereupon binding it cleaves the targeted nucleotide strand. To extend this genome-editing ability to the industrially important oleaginous yeast Yarrowia lipolytica, we developed a set of easily usable and effective CRISPR-Cas9 episomal vectors. In this protocols chapter, we first present a method by which arbitrary protein-coding genes can be disrupted via indel formation after CRISPR-Cas9 targeting. A second method demonstrates how the same CRISPR-Cas9 system can be used to induce markerless gene cassette integration into the genome by inducing homologous recombination after DNA cleavage by Cas9. Finally, we describe how a catalytically inactive form of Cas9 fused to a transcriptional repressor can be used to control transcription of native genes in Y. lipolytica. The CRISPR-Cas9 tools and strategies described here greatly increase the types of genome editing and transcriptional control that can be achieved in Y. lipolytica, and promise to facilitate more advanced engineering of this important oleaginous host.

ChIP-Seq Analysis in Neurospora crassa

Chromatin immunoprecipitation paired with next-generation sequencing (ChIP-seq) can be used to determine genome-wide distribution of transcriptions factors, transcriptional machinery, or histone modifications. DNA-protein interactions are covalently cross-linked with the addition of formaldehyde. Chromatin is prepared and sheared, then immunoprecipitated with the appropriate antibody. After reversal of cross-linking and treating with protease, the resulting DNA fragments are sequenced and mapped to the reference genome to determine overall enrichment. Here we describe a method of ChIP-seq for investigating protein-DNA interactions in the filamentous fungus Neurospora crassa.

Methods for Genomic Characterization and Maintenance of Anaerobic Fungi

The rapid development of molecular biology and bioinformatics has fueled renewed interests in anaerobic fungi from the phylum Neocallimastigomycota. This chapter presents well-established methods for isolation, routine cultivation, and cryopreservation of anaerobic fungi. Moreover, detailed nucleic acid extraction protocols are provided, which should enable readers to isolate high-quality DNA and RNA from a variety of anaerobic fungal culture media for downstream applications such as next-generation sequencing.

Evolutionary Adaptation to Generate Mutants

In this chapter we describe a method to generate mutants of filamentous fungi using their genomic plasticity and rapid adaptability to their environment. This method is based on spontaneous mutations occurring in relation to improved growth of fungi on media by repeated inoculation resulting in adaptation of the strain to the condition. The critical aspect of this method is the design of the selective media, which will depend strongly on the phenomenon that will be studied. This method is advantageous over UV or chemical random mutagenesis as it results in a lower frequency of undesired mutations and can result in strains that combined with (post)genomic approaches can enhance our understanding of the mechanisms driving various biological processes. In addition, it can be used to obtain better strains for various industrial applications. The method described here is specific for sporulating fungi and has so far not yet been tested for nonsporulating fungi.

The genome sequence of Bipolaris cookei reveals mechanisms of pathogenesis underlying target leaf spot of sorghum

Bipolaris cookei (=Bipolaris sorghicola) causes target leaf spot, one of the most prevalent foliar diseases of sorghum. Little is known about the molecular basis of pathogenesis in B. cookei, in large part due to a paucity of resources for molecular genetics, such as a reference genome. Here, a draft genome sequence of B. cookei was obtained and analyzed. A hybrid assembly strategy utilizing Illumina and Pacific Biosciences sequencing technologies produced a draft nuclear genome of 36.1 Mb, organized into 321 scaffolds with L50 of 31 and N50 of 378 kb, from which 11,189 genes were predicted. Additionally, a finished mitochondrial genome sequence of 135,790 bp was obtained, which contained 75 predicted genes. Comparative genomics revealed that B. cookei possessed substantially fewer carbohydrate-active enzymes and secreted proteins than closely related Bipolaris species. Novel genes involved in secondary metabolism, including genes implicated in ophiobolin biosynthesis, were identified. Among 37 B. cookei genes induced during sorghum infection, one encodes a putative effector with a limited taxonomic distribution among plant pathogenic fungi. The draft genome sequence of B. cookei provided novel insights into target leaf spot of sorghum and is an important resource for future investigation.

Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species

Filamentous fungi produce a diverse array of secondary metabolites (SMs) critical for defense, virulence, and communication. The metabolic pathways that produce SMs are found in contiguous gene clusters in fungal genomes, an atypical arrangement for metabolic pathways in other eukaryotes. Comparative studies of filamentous fungal species have shown that SM gene clusters are often either highly divergent or uniquely present in one or a handful of species, hampering efforts to determine the genetic basis and evolutionary drivers of SM gene cluster divergence. Here, we examined SM variation in 66 cosmopolitan strains of a single species, the opportunistic human pathogen Aspergillus fumigatus. Investigation of genome-wide within-species variation revealed 5 general types of variation in SM gene clusters: nonfunctional gene polymorphisms; gene gain and loss polymorphisms; whole cluster gain and loss polymorphisms; allelic polymorphisms, in which different alleles corresponded to distinct, nonhomologous clusters; and location polymorphisms, in which a cluster was found to differ in its genomic location across strains. These polymorphisms affect the function of representative A. fumigatus SM gene clusters, such as those involved in the production of gliotoxin, fumigaclavine, and helvolic acid as well as the function of clusters with undefined products. In addition to enabling the identification of polymorphisms, the detection of which requires extensive genome-wide synteny conservation (e.g., mobile gene clusters and nonhomologous cluster alleles), our approach also implicated multiple underlying genetic drivers, including point mutations, recombination, and genomic deletion and insertion events as well as horizontal gene transfer from distant fungi. Finally, most of the variants that we uncover within A. fumigatus have been previously hypothesized to contribute to SM gene cluster diversity across entire fungal classes and phyla. We suggest that the drivers of genetic diversity operating within a fungal species shown here are sufficient to explain SM cluster macroevolutionary patterns.

All organisms produce metabolites, which are small molecules important for growth, reproduction, and other essential functions. Some organisms, including fungi, plants, and bacteria, make specialized forms of metabolites known as “secondary” metabolites that are ecologically important and improve their producers’ chances of survival and reproduction. In fungi, the genes in pathways that synthesize secondary metabolites are typically located next to each other in the genome and organized in contiguous gene clusters. These gene clusters, along with the metabolites they produce, are highly distinct, even between otherwise similar fungi, and it is often difficult to reconstruct how these differences evolved. To understand how secondary metabolic pathways evolve in fungi, we compared secondary metabolic gene clusters in 66 strains of one species of filamentous fungus, the human pathogen Aspergillus fumigatus. We show that these gene clusters vary extensively within this species, and describe the genetic processes that cause these differences. We identify 5 types of variants: single nucleotide changes, gene and gene cluster gain and loss, different gene clusters at the same genomic position, and mobile gene clusters that “jump” around the genome. These results provide a road map to the types and frequencies of genomic changes underlying the extensive diversity of fungal secondary metabolites.

Lignin degradation: microorganisms, enzymes involved, genomes analysis and evolution

Extensive research efforts have been dedicated to describing degradation of wood, which is a complex process; hence, microorganisms have evolved different enzymatic and non-enzymatic strategies to utilize this plentiful plant material. This review describes a number of fungal and bacterial organisms which have developed both competitive and mutualistic strategies for the decomposition of wood and to thrive in different ecological niches. Through the analysis of the enzymatic machinery engaged in wood degradation, it was possible to elucidate different strategies of wood decomposition which often depend on ecological niches inhabited by given organism. Moreover, a detailed description of low molecular weight compounds is presented, which gives these organisms not only an advantage in wood degradation processes, but seems rather to be a new evolutionatory alternative to enzymatic combustion. Through analysis of genomics and secretomic data, it was possible to underline the probable importance of certain wood-degrading enzymes produced by different fungal organisms, potentially giving them advantage in their ecological niches. The paper highlights different fungal strategies of wood degradation, which possibly correlates to the number of genes coding for secretory enzymes. Furthermore, investigation of the evolution of wood-degrading organisms has been described.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Genome sequence of the ectophytic fungus Ramichloridium luteum reveals unique evolutionary adaptations to plant surface niche

BACKGROUND

Ectophytic fungi occupy the waxy plant surface, an extreme environment characterized by prolonged desiccation, nutrient limitation, and exposure to solar radiation. The nature of mechanisms that facilitate adaptation to this environment remains unclear. In this study, we sequenced the complete genome of an ectophytic fungus, Ramichloridium luteum, which colonizes the surface of apple fruit, and carried out comparative genomic and transcriptome analysis.

RESULTS

The R. luteum genome was 28.18 Mb and encoded 9466 genes containing 1.85% repetitive elements. Compared with cell-penetrating pathogens, genes encoding plant cell wall degrading enzymes (PCWDEs), PTH11-like G protein-coupled receptors (GPCRs) and effectors were drastically reduced. In contrast, genes encoding cutinases and secretory lipases were strikingly expanded, and four of nine secretory lipases were probably acquired by horizontal gene transfer from Basidiomycota. Transcriptomic analysis revealed elevated expression of genes involved in cuticle degradation (cutinase, secretory lipase) and stress responses (melanin biosynthesis, aquaporins, lysozymes and HOG pathway).

CONCLUSIONS

Taken together, our results highlight genomic features associated with evolution of surface niche adaptation by the ectophytic fungus R. luteum, namely the contraction of PCWDEs, PTH11-like GPCRs and effectors, and the expansion of cuticle degradation and stress tolerance.

Intra and Interspecific Variations of Gene Expression Levels in Yeast Are Largely Neutral: (Nei Lecture, SMBE 2016, Gold Coast)

It is commonly, although not universally, accepted that most intra and interspecific genome sequence variations are more or less neutral, whereas a large fraction of organism-level phenotypic variations are adaptive. Gene expression levels are molecular phenotypes that bridge the gap between genotypes and corresponding organism-level phenotypes. Yet, it is unknown whether natural variations in gene expression levels are mostly neutral or adaptive. Here we address this fundamental question by genome-wide profiling and comparison of gene expression levels in nine yeast strains belonging to three closely related Saccharomyces species and originating from five different ecological environments. We find that the transcriptome-based clustering of the nine strains approximates the genome sequence-based phylogeny irrespective of their ecological environments. Remarkably, only ∼0.5% of genes exhibit similar expression levels among strains from a common ecological environment, no greater than that among strains with comparable phylogenetic relationships but different environments. These and other observations strongly suggest that most intra and interspecific variations in yeast gene expression levels result from the accumulation of random mutations rather than environmental adaptations. This finding has profound implications for understanding the driving force of gene expression evolution, genetic basis of phenotypic adaptation, and general role of stochasticity in evolution.

Confirmation of independent introductions of an exotic plant pathogen of Cornus species, Discula destructiva, on the east and west coasts of North America

Cornus florida (flowering dogwood) and C. nuttallii (Pacific dogwood) are North American native tree species that belong to the big-bracted group of dogwoods. Cornus species are highly valued for their ornamental characteristics, and have fruits that contain high fat content for animals. Also, they are an important understory tree in natural forests. Dogwood anthracnose, caused by Discula destructiva, was observed in the late 1970s on the east and west coasts of the United States and by 1991 had quickly spread throughout most of the native ranges of C. florida and C. nuttalli. We investigated the genetic diversity and population structure of 93 D. destructiva isolates using 47 microsatellite loci developed from the sequenced genome of the type strain of D. destructiva. Clone-corrected data indicated low genetic diversity and the presence of four genetic clusters that corresponded to two major geographic areas, the eastern United States and the Pacific Northwest, and to the two collection time periods when the isolates were collected (pre- and post-1993). Linkage disequilibrium was present in five out of six subpopulations, suggesting that the fungus only reproduced asexually. Evidence of population bottlenecks was indicated across four identified genetic clusters, and was probably the result of the limited number of founding individuals on both coasts. These results support the hypothesis that D. destructiva is an exotic pathogen with independent introductions on the east and west coasts of North America. We also tested the cross-amplification of these microsatellite primers to other Discula species. Genomic DNA from 17 isolates of four other Discula species and two isolates of Juglanconis species (formerly Melanconis species) were amplified by 17 of 47 primer pairs. These primers may be useful for investigating the genetic diversity and population structure of these Discula species.

Spatial organization of the budding yeast genome in the cell nucleus and identification of specific chromatin interactions from multi-chromosome constrained chromatin model

Nuclear landmarks and biochemical factors play important roles in the organization of the yeast genome. The interaction pattern of budding yeast as measured from genome-wide 3C studies are largely recapitulated by model polymer genomes subject to landmark constraints. However, the origin of inter-chromosomal interactions, specific roles of individual landmarks, and the roles of biochemical factors in yeast genome organization remain unclear. Here we describe a multi-chromosome constrained self-avoiding chromatin model (mC-SAC) to gain understanding of the budding yeast genome organization. With significantly improved sampling of genome structures, both intra- and inter-chromosomal interaction patterns from genome-wide 3C studies are accurately captured in our model at higher resolution than previous studies. We show that nuclear confinement is a key determinant of the intra-chromosomal interactions, and centromere tethering is responsible for the inter-chromosomal interactions. In addition, important genomic elements such as fragile sites and tRNA genes are found to be clustered spatially, largely due to centromere tethering. We uncovered previously unknown interactions that were not captured by genome-wide 3C studies, which are found to be enriched with tRNA genes, RNAPIII and TFIIS binding. Moreover, we identified specific high-frequency genome-wide 3C interactions that are unaccounted for by polymer effects under landmark constraints. These interactions are enriched with important genes and likely play biological roles.

The architecture of the cell nucleus and the spatial organization of the genome are important in determining nuclear functions. Single-cell imaging techniques and chromosome conformation capture (3C) based methods have provided a wealth of information on the spatial organization of chromosomes. Here we describe a multi-chromosome ensemble model of chromatin chains for understanding the folding principles of budding yeast genome. By overcoming severe challenges in sampling self-avoiding chromatin chains in nuclear confinement, we succeed in generating a large number of model genomes of budding yeast. Our model predicts chromatin interactions that have good correlation with experimental measurements. Our results showed that the spatial confinement of cell nucleus and excluded-volume effect are key determinants of the folding behavior of yeast chromosomes, and largely account for the observed intra-chromosomal interactions. Furthermore, we determined the specific roles of individual nuclear landmarks and biochemical factors, and our analysis showed that centromere tethering largely determines inter-chromosomal interactions. In addition, we were able to infer biological properties from the organization of modeled genomes. We found that the spatial locations of important elements such as fragile sites and tRNA genes are largely determined by the tethering of centromeres to the Spindle Pole Body. We further showed that many of these spatial locations can be predicted by using the genomic distances to the centromeres. Overall, our results revealed important insight into the organizational principles of the budding yeast genome and predicted a number of important biological findings that are fully experimentally testable.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

Large-scale molecular genetic analysis in plant-pathogenic fungi: a decade of genome-wide functional analysis

Plant-pathogenic fungi cause diseases to all major crop plants world-wide and threaten global food security. Underpinning fungal diseases are virulence genes facilitating plant host colonization that often marks pathogenesis and crop failures, as well as an increase in staple food prices. Fungal molecular genetics is therefore the cornerstone to the sustainable prevention of disease outbreaks. Pathogenicity studies using mutant collections provide immense function-based information regarding virulence genes of economically relevant fungi. These collections are rich in potential targets for existing and new biological control agents. They contribute to host resistance breeding against fungal pathogens and are instrumental in searching for novel resistance genes through the identification of fungal effectors. Therefore, functional analyses of mutant collections propel gene discovery and characterization, and may be incorporated into disease management strategies. In the light of these attributes, mutant collections enhance the development of practical solutions to confront modern agricultural constraints. Here, a critical review of mutant collections constructed by various laboratories during the past decade is provided. We used Magnaporthe oryzae and Fusarium graminearum studies to show how mutant screens contribute to bridge existing knowledge gaps in pathogenicity and fungal-host interactions.

Network-assisted target identification for haploinsufficiency and homozygous profiling screens

Chemical genomic screens have recently emerged as a systematic approach to drug discovery on a genome-wide scale. Drug target identification and elucidation of the mechanism of action (MoA) of hits from these noisy high-throughput screens remain difficult. Here, we present GIT (Genetic Interaction Network-Assisted Target Identification), a network analysis method for drug target identification in haploinsufficiency profiling (HIP) and homozygous profiling (HOP) screens. With the drug-induced phenotypic fitness defect of the deletion of a gene, GIT also incorporates the fitness defects of the gene's neighbors in the genetic interaction network. On three genome-scale yeast chemical genomic screens, GIT substantially outperforms previous scoring methods on target identification on HIP and HOP assays, respectively. Finally, we showed that by combining HIP and HOP assays, GIT further boosts target identification and reveals potential drug's mechanism of action.

Dynamic genome-scale metabolic modeling of the yeast Pichia pastoris

BACKGROUND

Pichia pastoris shows physiological advantages in producing recombinant proteins, compared to other commonly used cell factories. This yeast is mostly grown in dynamic cultivation systems, where the cell's environment is continuously changing and many variables influence process productivity. In this context, a model capable of explaining and predicting cell behavior for the rational design of bioprocesses is highly desirable. Currently, there are five genome-scale metabolic reconstructions of P. pastoris which have been used to predict extracellular cell behavior in stationary conditions.

RESULTS

In this work, we assembled a dynamic genome-scale metabolic model for glucose-limited, aerobic cultivations of Pichia pastoris. Starting from an initial model structure for batch and fed-batch cultures, we performed pre/post regression diagnostics to ensure that model parameters were identifiable, significant and sensitive. Once identified, the non-relevant ones were iteratively fixed until a priori robust modeling structures were found for each type of cultivation. Next, the robustness of these reduced structures was confirmed by calibrating the model with new datasets, where no sensitivity, identifiability or significance problems appeared in their parameters. Afterwards, the model was validated for the prediction of batch and fed-batch dynamics in the studied conditions. Lastly, the model was employed as a case study to analyze the metabolic flux distribution of a fed-batch culture and to unravel genetic and process engineering strategies to improve the production of recombinant Human Serum Albumin (HSA). Simulation of single knock-outs indicated that deviation of carbon towards cysteine and tryptophan formation improves HSA production. The deletion of methylene tetrahydrofolate dehydrogenase could increase the HSA volumetric productivity by 630%. Moreover, given specific bioprocess limitations and strain characteristics, the model suggests that implementation of a decreasing specific growth rate during the feed phase of a fed-batch culture results in a 25% increase of the volumetric productivity of the protein.

CONCLUSION

In this work, we formulated a dynamic genome scale metabolic model of Pichia pastoris that yields realistic metabolic flux distributions throughout dynamic cultivations. The model can be calibrated with experimental data to rationally propose genetic and process engineering strategies to improve the performance of a P. pastoris strain of interest.

Draft genome sequence of Karnal bunt pathogen (Tilletia indica) of wheat provides insights into the pathogenic mechanisms of quarantined fungus

Karnal bunt disease in wheat is caused by hemibiotrophic fungus, Tilletia indica that has been placed as quarantine pest in more than 70 countries. Despite its economic importance, little knowledge about the molecular components of fungal pathogenesis is known. In this study, first time the genome sequence of T. indica has been deciphered for unraveling the effectors' functions of molecular pathogenesis of Karnal bunt disease. The T. indica genome was sequenced employing hybrid approach of PacBio Single Molecule Real Time (SMRT) and Illumina HiSEQ 2000 sequencing platforms. The genome was assembled into 10,957 contigs (N50 contig length 3 kb) with total size of 26.7 Mb and GC content of 53.99%. The number of predicted putative genes were 11,535, which were annotated with Gene Ontology databases. Functional annotation of Karnal bunt pathogen genome and classification of identified effectors into protein families revealed interesting functions related to pathogenesis. Search for effectors' genes using pathogen host interaction database identified 135 genes. The T. indica genome sequence and putative genes involved in molecular pathogenesis would further help in devising novel and effective disease management strategies including development of resistant wheat genotypes, novel biomarkers for pathogen detection and new targets for fungicide development.

Genome-wide exonic small interference RNA-mediated gene silencing regulates sexual reproduction in the homothallic fungus Fusarium graminearum

Various ascomycete fungi possess sex-specific molecular mechanisms, such as repeat-induced point mutations, meiotic silencing by unpaired DNA, and unusual adenosine-to-inosine RNA editing, for genome defense or gene regulation. Using a combined analysis of functional genetics and deep sequencing of small noncoding RNA (sRNA), mRNA, and the degradome, we found that the sex-specifically induced exonic small interference RNA (ex-siRNA)-mediated RNA interference (RNAi) mechanism has an important role in fine-tuning the transcriptome during ascospore formation in the head blight fungus Fusarium graminearum. Approximately one-third of the total sRNAs were produced from the gene region, and sRNAs with an antisense direction or 5'-U were involved in post-transcriptional gene regulation by reducing the stability of the corresponding gene transcripts. Although both Dicers and Argonautes partially share their functions, the sex-specific RNAi pathway is primarily mediated by FgDicer1 and FgAgo2, while the constitutively expressed RNAi components FgDicer2 and FgAgo1 are responsible for hairpin-induced RNAi. Based on our results, we concluded that F. graminearum primarily utilizes ex-siRNA-mediated RNAi for ascosporogenesis but not for genome defenses and other developmental stages. Each fungal species appears to have evolved RNAi-based gene regulation for specific developmental stages or stress responses. This study provides new insights into the regulatory role of sRNAs in fungi and other lower eukaryotes.

Biology, dynamics, and applications of transposable elements in basidiomycete fungi

The phylum Basidiomycota includes filamentous fungi and yeast species with different ecological and genomic characteristics. Transposable elements (TEs) are abundant components of most eukaryotic genomes, and their transition from being genomic parasites to key drivers of genomic architecture, functionality, and evolution is a subject receiving much attention. In light of the abundant genomic information released during the last decade, the aims of this mini-review are to discuss the dynamics and impact of TEs in basidiomycete fungi. To do this, we surveyed and explored data from 75 genomes, which encompass the phylogenetic diversity of the phylum Basidiomycota. We describe annotation approaches and analyze TE distribution in the context of species phylogeny and genome size. Further, we review the most relevant literature about the role of TEs in species lifestyle, their impact on genome architecture and functionality, and the defense mechanisms evolved to control their proliferation. Finally, we discuss potential applications of TEs that can drive future innovations in fungal research.

Exploring Protein Function Using the Saccharomyces Genome Database

Elucidating the function of individual proteins will help to create a comprehensive picture of cell biology, as well as shed light on human disease mechanisms, possible treatments, and cures. Due to its compact genome, and extensive history of experimentation and annotation, the budding yeast Saccharomyces cerevisiae is an ideal model organism in which to determine protein function. This information can then be leveraged to infer functions of human homologs. Despite the large amount of research and biological data about S. cerevisiae, many proteins' functions remain unknown. Here, we explore ways to use the Saccharomyces Genome Database (SGD; http://www.yeastgenome.org ) to predict the function of proteins and gain insight into their roles in various cellular processes.

Haplotyping a Non-meiotic Diploid Fungal Pathogen Using Induced Aneuploidies and SNP/CGH Microarray Analysis

The generation of haplotype information has recently become very attractive due to its utility for identifying mutations associated with human disease and for the development of personalized medicine. Haplotype information also is crucial for studying recombination mechanisms and genetic diversity, and for analyzing allele-specific gene expression. Classic haplotyping methods require the analysis of hundreds of meiotic progeny. To facilitate haplotyping in the non-meiotic human fungal pathogen Candida albicans, we exploited trisomic heterozygous chromosomes generated via the UAU1 selection strategy. Using this system, we obtained phasing information from allelic biases, detected by SNP/CGH microarray analysis. This strategy has the potential to be applicable to other diploid, asexual Candida species that are important causes of human disease.

Comparative study on fermentation performance in the genome shuffled Candida versatilis and wild-type salt tolerant yeast strain

BACKGROUND

The fermentation performance of a genome-shuffled strain of Candida versatilis S3-5, isolated for improved tolerance to salt, and wild-type (WT) strain were analysed. The fermentation parameters, such as growth, reducing sugar, ethanol, organic acids and volatile compounds, were detected during soy sauce fermentation process.

RESULTS

The results showed that ethanol produced by the genome shuffled strain S3-5 was increasing at a faster rate and to a greater extent than WT. At the end of the fermentation, malic acid, citric acid and succinic acid formed in tricarboxylic acid cycle after S3-5 treatment elevated by 39.20%, 6.85% and 17.09% compared to WT, respectively. Moreover, flavour compounds such as phenethyl acetate, ethyl vanillate, ethyl acetate, isoamyl acetate, ethyl myristate, ethyl pentadecanoate, ethyl palmitate and phenylacetaldehyde produced by S3-5 were 2.26, 2.12, 2.87, 34.41, 6.32, 13.64, 2.23 and 78.85 times as compared to WT.

CONCLUSIONS

S3-5 exhibited enhanced metabolic ability as compared to the wild-type strain, improved conversion of sugars to ethanol, metabolism of organic acid and formation of volatile compounds, especially esters, Moreover, S3-5 might be an ester-flavour type salt-tolerant yeast. © 2016 Society of Chemical Industry.

Refined Pichia pastoris reference genome sequence

Strains of the species Komagataella phaffii are the most frequently used "Pichia pastoris" strains employed for recombinant protein production as well as studies on peroxisome biogenesis, autophagy and secretory pathway analyses. Genome sequencing of several different P. pastoris strains has provided the foundation for understanding these cellular functions in recent genomics, transcriptomics and proteomics experiments. This experimentation has identified mistakes, gaps and incorrectly annotated open reading frames in the previously published draft genome sequences. Here, a refined reference genome is presented, generated with genome and transcriptome sequencing data from multiple P. pastoris strains. Twelve major sequence gaps from 20 to 6000 base pairs were closed and 5111 out of 5256 putative open reading frames were manually curated and confirmed by RNA-seq and published LC-MS/MS data, including the addition of new open reading frames (ORFs) and a reduction in the number of spliced genes from 797 to 571. One chromosomal fragment of 76kbp between two previous gaps on chromosome 1 and another 134kbp fragment at the end of chromosome 4, as well as several shorter fragments needed re-orientation. In total more than 500 positions in the genome have been corrected. This reference genome is presented with new chromosomal numbering, positioning ribosomal repeats at the distal ends of the four chromosomes, and includes predicted chromosomal centromeres as well as the sequence of two linear cytoplasmic plasmids of 13.1 and 9.5kbp found in some strains of P. pastoris.

Cercospora zeina from Maize in South Africa Exhibits High Genetic Diversity and Lack of Regional Population Differentiation

South Africa is one of the leading maize-producing countries in sub-Saharan Africa. Since the 1980s, Cercospora zeina, a causal agent of gray leaf spot of maize, has become endemic in South Africa, and is responsible for substantial yield reductions. To assess genetic diversity and population structure of C. zeina in South Africa, 369 isolates were collected from commercial maize farms in three provinces (KwaZulu-Natal, Mpumalanga, and North West). These isolates were evaluated with 14 microsatellite markers and species-specific mating type markers that were designed from draft genome sequences of C. zeina isolates from Africa (CMW 25467) and the United States (USPA-4). Sixty alleles were identified across 14 loci, and gene diversity values within each province ranged from 0.18 to 0.35. High levels of gene flow were observed (Nm = 5.51), and in a few cases, identical multilocus haplotypes were found in different provinces. Overall, 242 unique multilocus haplotypes were identified with a low clonal fraction of 34%. No distinct population clusters were identified using STRUCTURE, principal coordinate analysis, or Weir's theta θ statistic. The lack of population differentiation was supported by analysis of molecular variance tests, which indicated that only 2% of the variation was attributed to variability between populations from each province. Mating type ratios of MAT1-1 and MAT1-2 idiomorphs from 335 isolates were not significantly different from a 1:1 ratio in all provinces, which provided evidence for sexual reproduction. The draft genome of C. zeina CMW 25467 exhibited a complete genomic copy of the MAT1-1 idiomorph as well as exonic fragments of MAT genes from both idiomorphs. The high level of gene diversity, shared haplotypes at different geographical locations within South Africa, and presence of both MAT idiomorphs at all sites indicates widespread dispersal of C. zeina between maize fields in the country as well as evidence for sexual recombination. The outcomes of this genome-enabled study are important for disease management since the high diversity has implications for dispersal of fungicide resistance should it emerge and the need for diversified resistance breeding.

Molecular diagnostic toolkit for Rhizophagus irregularis isolate DAOM-197198 using quantitative PCR assay targeting the mitochondrial genome

Rhizophagus irregularis (previously named Glomus irregulare) is one of the most widespread and common arbuscular mycorrhizal fungal (AMF) species. It has been recovered worldwide in agricultural and natural soils, and the isolate DAOM-197198 has been utilized as a commercial inoculant for two decades. Despite the ecological and economical importance of this taxon, specific markers for quantification of propagules by quantitative real-time PCR (qPCR) are extremely limited and none have been rigorously validated for quality control of manufactured products such as biofertilizers. From the sequencing of 14 complete AMF mitochondrial (mt) genomes, a qPCR assay using a hydrolysis probe designed in the single copy cox3-rnl intergenic region was tested and validated to specifically and accurately quantify the spores of R. irregularis isolate DAOM-197198. Specificity tests were performed using standard PCR and qPCR, and results clearly showed that the primers specifically amplified the isolate DAOM-197198, yielding a PCR product of 106 bp. According to the qPCR analyses on spores produced in vitro, the average copy number of mt genomes per spore was 3172 ± 304 SE (n = 6). Quantification assays were successfully undertaken on known and unknown samples in liquid suspensions and commercial dry formulations to show the accuracy, precision, robustness, and reproducibility of the qPCR assay. This study provides a powerful molecular toolkit specifically designed to quantify spores of the model AMF isolate DAOM-197198. The approach of molecular toolkit used in our study could be applied to other AMF taxa and will be useful to research institutions and governmental and industrial laboratories running routine quality control of AMF-based products.

System-level genome editing in microbes

The release of the first complete microbial genome sequences at the end of the past century opened the way for functional genomics and systems-biology to uncover the genetic basis of various phenotypes. The surge of available sequence data facilitated the development of novel genome editing techniques for system-level analytical studies. Recombineering allowed unprecedented throughput and efficiency in microbial genome editing and the recent discovery and widespread use of RNA-guided endonucleases offered several further perspectives: (i) previously recalcitrant species became editable, (ii) the efficiency of recombineering could be elevated, and as a result (iii) diverse genomic libraries could be generated more effectively. Supporting recombineering by RNA-guided endonucleases has led to success stories in metabolic engineering, but their use for system-level analysis is mostly unexplored. For the full exploitation of opportunities that are offered by the genome editing proficiency, future development of large scale analytical procedures is also vitally needed.

The large (134.9 kb) mitochondrial genome of the glomeromycete Funneliformis mosseae

Funneliformis mosseae is among the most ecologically and economically important glomeromycete species and occurs both in natural and disturbed areas in a wide range of habitats and climates. In this study, we report the sequencing of the complete mitochondrial (mt) genome of F. mosseae isolate FL299 using 454 pyrosequencing and Illumina HiSeq technologies. This mt genome is a full-length circular chromosome of 134,925 bp, placing it among the largest mitochondrial DNAs (mtDNAs) in the fungal kingdom. A comparative analysis with publically available arbuscular mycorrhizal fungal mtDNAs revealed that the mtDNA of F. mosseae FL299 contained a very large number of insertions contributing to its expansion. The gene synteny was completely reshuffled compared to previously published glomeromycotan mtDNAs and several genes were oriented in an anti-sense direction. Furthermore, the presence of different types of introns and insertions in rnl (14 introns) made this gene very distinctive in Glomeromycota. The presence of alternative genetic codes in both initiation (GUG) and termination (UGA) codons was another new feature in this mtDNA compared to previously published glomeromycotan mt genomes. The phylogenetic analysis inferred from the analysis of 14 protein mt genes confirmed the position of the Glomeromycota clade as a sister group of Mortierellomycotina. This mt genome is the largest observed so far in Glomeromycota and the first mt genome within the Funneliformis clade, providing new opportunities to better understand their evolution and to develop molecular markers.

Genome Anatomy of Pyrenochaeta unguis-hominis UM 256, a Multidrug Resistant Strain Isolated from Skin Scraping

Pyrenochaeta unguis-hominis is a rare human pathogen that causes infection in human skin and nail. P. unguis-hominis has received little attention, and thus, the basic biology and pathogenicity of this fungus is not fully understood. In this study, we performed in-depth analysis of the P. unguis-hominis UM 256 genome that was isolated from the skin scraping of a dermatitis patient. The isolate was identified to species level using a comprehensive multilocus phylogenetic analysis of the genus Pyrenochaeta. The assembled UM 256 genome has a size of 35.5 Mb and encodes 12,545 putative genes, and 0.34% of the assembled genome is predicted transposable elements. Its genomic features propose that the fungus is a heterothallic fungus that encodes a wide array of plant cell wall degrading enzymes, peptidases, and secondary metabolite biosynthetic enzymes. Antifungal drug resistance genes including MDR, CDR, and ERG11/CYP51 were identified in P. unguis-hominis UM 256, which may confer resistance to this fungus. The genome analysis of P. unguis-hominis provides an insight into molecular and genetic basis of the fungal lifestyles, understanding the unrevealed biology of antifungal resistance in this fungus.

Genomic Analyses of Cladophialophora bantiana, a Major Cause of Cerebral Phaeohyphomycosis Provides Insight into Its Lifestyle, Virulence and Adaption in Host

Cladophialophora bantiana is a dematiaceous fungus with a predilection for causing central nervous system (CNS) infection manifesting as brain abscess in both immunocompetent and immunocompromised patients. In this paper, we report comprehensive genomic analyses of C. bantiana isolated from the brain abscess of an immunocompetent man, the first reported case in Malaysia and Southeast Asia. The identity of the fungus was determined using combined morphological analysis and multilocus phylogeny. The draft genome sequence of a neurotrophic fungus, C. bantiana UM 956 was generated using Illumina sequencing technology to dissect its genetic fundamental and basic biology. The assembled 37.1 Mb genome encodes 12,155 putative coding genes, of which, 1.01% are predicted transposable elements. Its genomic features support its saprophytic lifestyle, renowned for its versatility in decomposing hemicellulose and pectin components. The C. bantiana UM 956 was also found to carry some important putative genes that engaged in pathogenicity, iron uptake and homeostasis as well as adaptation to various stresses to enable the organism to survive in hostile microenvironment. This wealth of resource will further catalyse more downstream functional studies to provide better understanding on how this fungus can be a successful and persistent pathogen in human.