Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

Delimitation and phylogeny of Dictyochaeta, and introduction of Achrochaeta and Tubulicolla, genera nova

Dictyochaeta (Chaetosphaeriaceae) is a phialidic dematiaceous hyphomycete with teleomorphs classified in Chaetosphaeria. It is associated with significant variability of asexual morphological traits, which led to its broad delimitation. In the present study, six loci: nuc rDNA internal transcribed spacer region ITS1-5.8S-ITS2 (ITS barcode), nuc 18S rDNA (18S), nuc 28S rDNA (28S), DNA-directed RNA polymerase II second largest subunit gene (RPB2), translation elongation factor 1-α (TEF1-α), and β-tubulin (TUB2), along with comparative morphological and cultivation studies, are used to reevaluate the concept of Dictyochaeta and establish species boundaries. Based on revised species, morphological characteristics of conidia (shape, septation, absence or presence of setulae), collarettes (shape), and setae (presence or absence) and an extension of the conidiogenous cell proved to be important at the generic level. The dual DNA barcoding using ITS and TEF1-α, together with TUB2, facilitated accurate identification of Dictyochaeta species. Thirteen species are accepted, of which seven are characterized in this study; an identification key is provided. It was revealed that D. fuegiana, the type species, is a complex of three distinct species including D. querna and the newly described D. stratosa. Besides, a new species, D. detriticola, and two new combinations, D. callimorpha and D. montana, are proposed. An epitype of D. montana is selected. Dictyochaeta includes saprobes on decaying wood, bark, woody fruits, and fallen leaves. Dictyochaeta is shown to be distantly related to the morphologically similar Codinaea, which is resolved as paraphyletic. Chaetosphaeria talbotii with a Dictyochaeta anamorph represents a novel lineage in the Chaetosphaeriaceae; it is segregated from Dictyochaeta, and a new genus Achrochaeta is proposed. Multigene phylogenetic analysis revealed that D. cylindrospora belongs to the Vermiculariopsiellales, and a new genus Tubulicolla is introduced.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

Transposon-mediated telomere destabilization: a driver of genome evolution in the blast fungus

The fungus Magnaporthe oryzae causes devastating diseases of crops, including rice and wheat, and in various grasses. Strains from ryegrasses have highly unstable chromosome ends that undergo frequent rearrangements, and this has been associated with the presence of retrotransposons (Magnaporthe oryzae Telomeric Retrotransposons-MoTeRs) inserted in the telomeres. The objective of the present study was to determine the mechanisms by which MoTeRs promote telomere instability. Targeted cloning, mapping, and sequencing of parental and novel telomeric restriction fragments (TRFs), along with MinION sequencing of genomic DNA allowed us to document the precise molecular alterations underlying 109 newly-formed TRFs. These included truncations of subterminal rDNA sequences; acquisition of MoTeR insertions by 'plain' telomeres; insertion of the MAGGY retrotransposons into MoTeR arrays; MoTeR-independent expansion and contraction of subtelomeric tandem repeats; and a variety of rearrangements initiated through breaks in interstitial telomere tracts that are generated during MoTeR integration. Overall, we estimate that alterations occurred in approximately sixty percent of chromosomes (one in three telomeres) analyzed. Most importantly, we describe an entirely new mechanism by which transposons can promote genomic alterations at exceptionally high frequencies, and in a manner that can promote genome evolution while minimizing collateral damage to overall chromosome architecture and function.

Telomere-to-telomere assembled and centromere annotated genomes of the two main subspecies of the button mushroom Agaricus bisporus reveal especially polymorphic chromosome ends

Agaricus bisporus, the most cultivated edible mushroom worldwide, is represented mainly by the subspecies var. bisporus and var. burnettii. var. bisporus has a secondarily homothallic life cycle with recombination restricted to chromosome ends, while var. burnettii is heterothallic with recombination seemingly equally distributed over the chromosomes. To better understand the relationship between genomic make-up and different lifestyles, we have de novo sequenced a burnettii homokaryon and synchronised gene annotations with updated versions of the published genomes of var. bisporus. The genomes were assembled into telomere-to-telomere chromosomes and a consistent set of gene predictions was generated. The genomes of both subspecies were largely co-linear, and especially the chromosome ends differed in gene model content between the two subspecies. A single large cluster of repeats was found on each chromosome at the same respective position in all strains, harbouring nearly 50% of all repeats and likely representing centromeres. Repeats were all heavily methylated. Finally, a mapping population of var. burnettii confirmed an even distribution of crossovers in meiosis, contrasting the recombination landscape of var. bisporus. The new findings using the exceptionally complete and well annotated genomes of this basidiomycete demonstrate the importance for unravelling genetic components underlying the different life cycles.

Novel Variation and Evolution of AvrPiz-t of Magnaporthe oryzae in Field Isolates

The product of the avirulence (Avr) gene of Magnaporthe oryzae can be detected by the product of the corresponding resistance (R) gene of rice and activates immunity to rice mediated by the R gene. The high degree of variability of M. oryzae isolates in pathogenicity makes the control of rice blast difficult. That resistance of the R gene in rice has been lost has been ascribed to the instability of the Avr gene in M. oryzae. Further study on the variation of the Avr genes in M. oryze field isolates may yield valuable information on the durable and effective deployment of R genes in rice production areas. AvrPiz-t and Piz-t are a pair of valuable genes in the Rice-Magnaporthe pathosystem. AvrPiz-t is detectable by Piz-t and determines the effectiveness of Piz-t. To effectively deploy the R gene Piz-t, the distribution, variation, and evolution of the corresponding Avr gene AvrPiz-t were found among 312 M. oryzae isolates collected from Yunnan rice production areas of China. PCR amplification and pathogenicity assays of AvrPiz-t showed that 202 isolates (64.7%) held AvrPiz-t alleles and were avirulent to IRBLzt-T (holding Piz-t). There were 42.3–83.3% avirulent isolates containing AvrPiz-t among seven regions in Yunnan Province. Meanwhile, 11 haplotypes of AvrPiz-t encoding three novel AvrPiz-t variants were identified among 100 isolates. A 198 bps insertion homologous to solo-LTR of the retrotransposon inago2 in the promoter region of AvrPiz-t in one isolate and a frameshift mutation of CDS in another isolate were identified among 100 isolates, and those two isolates had evolved to virulent from avirulent. Synonymous mutation and non-AUG-initiated N-terminal extensions keeps the AvrPiz-t gene avirulence function in M. oryzae field isolates in Yunnan. A haplotype network showed that H3 was an ancestral haplotype. Structure variance for absence (28.2%) or partial fragment loss (71.8%) of AvrPiz-t was found among 39 virulent isolates and may cause the AvrPiz-t avirulence function to be lost. Overall, AvrPiz-t evolved to virulent from avirulent forms via point mutation, retrotransposon, shift mutation, and structure variance under field conditions.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Analysis of Epichloë festucae small secreted proteins in the interaction with Lolium perenne

Epichloë festucae is an endophyte of the agriculturally important perennial ryegrass. This species systemically colonises the aerial tissues of this host where its growth is tightly regulated thereby maintaining a mutualistic symbiotic interaction. Recent studies have suggested that small secreted proteins, termed effectors, play a vital role in the suppression of host defence responses. To date only a few effectors with important roles in mutualistic interactions have been described. Here we make use of the fully assembled E. festucae genome and EffectorP to generate a suite of 141 effector candidates. These were analysed with respect to their genome location and expression profiles in planta and in several symbiosis-defective mutants. We found an association between effector candidates and a class of transposable elements known as MITEs, but no correlation with other dynamic features of the E. festucae genome, such as transposable element-rich regions. Three effector candidates and a small GPI-anchored protein were chosen for functional analysis based on their high expression in planta compared to in culture and their differential regulation in symbiosis defective E. festucae mutants. All three candidate effector proteins were shown to possess a functional signal peptide and two could be detected in the extracellular medium by western blotting. Localization of the effector candidates in planta suggests that they are not translocated into the plant cell, but rather, are localized in the apoplastic space or are attached to the cell wall. Deletion and overexpression of the effector candidates, as well as the putative GPI-anchored protein, did not affect the plant growth phenotype or restrict growth of E. festucae mutants in planta. These results indicate that these proteins are either not required for the interaction at the observed life stages or that there is redundancy between effectors expressed by E. festucae.

The Structure of a Conserved Telomeric Region Associated with Variant Antigen Loci in the Blood Parasite Trypanosoma congolense

African trypanosomiasis is a vector-borne disease of humans and livestock caused by African trypanosomes (Trypanosoma spp.). Survival in the vertebrate bloodstream depends on antigenic variation of Variant Surface Glycoproteins (VSGs) coating the parasite surface. In T. brucei, a model for antigenic variation, monoallelic VSG expression originates from dedicated VSG expression sites (VES). Trypanosoma brucei VES have a conserved structure consisting of a telomeric VSG locus downstream of unique, repeat sequences, and an independent promoter. Additional protein-coding sequences, known as “Expression Site Associated Genes (ESAGs)”, are also often present and are implicated in diverse, bloodstream-stage functions. Trypanosoma congolense is a related veterinary pathogen, also displaying VSG-mediated antigenic variation. A T. congolense VES has not been described, making it unclear if regulation of VSG expression is conserved between species. Here, we describe a conserved telomeric region associated with VSG loci from long-read DNA sequencing of two T. congolense strains, which consists of a distal repeat, conserved noncoding elements and other genes besides the VSG; although these are not orthologous to T. brucei ESAGs. Most conserved telomeric regions are associated with accessory minichromosomes, but the same structure may also be associated with megabase chromosomes. We propose that this region represents the T. congolense VES, and through comparison with T. brucei, we discuss the parallel evolution of antigenic switching mechanisms, and unique adaptation of the T. brucei VES for developmental regulation of bloodstream-stage genes. Hence, we provide a basis for understanding antigenic switching in T. congolense and the origins of the African trypanosome VES.

Phylogenetic classification and generic delineation of Calyptosphaeria gen. nov., Lentomitella, Spadicoides and Torrentispora (Sordariomycetes)

The genus Ceratostomella has a long history of taxonomic confusion. While species with evanescent asci have been transferred to the Microascales and Ophiostomatales, the taxonomic status of species with persistent asci has not been completely resolved. In previous studies using DNA sequence data, cultures and morphology, several Ceratostomella spp. were allocated in 13 genera in the Eurotiomycetes and Sordariomycetes. In our study, the systematics of the remaining Ceratostomella spp. with persistent asci is revisited with new collection data, cultures and phylogeny based on novel DNA sequences from six nuclear loci. Bayesian inference and Maximum Likelihood analyses support the monophyly of several wood-inhabiting species formerly classified in Ceratostomella and other unknown morphologically similar taxa and their division into four genera, i.e. Lentomitella, Spadicoides, Torrentispora and the newly described Calyptosphaeria. This robust clade represents the order Xenospadicoidales in the Sordariomycetidae. Comparative analysis of the ITS2 secondary structure revealed a genetic variation among Lentomitella isolates; 11 species were recognised, of which five are newly introduced and two are new combinations. Other taxonomic novelties include four new species and eight new combinations in Calyptosphaeria, Spadicoides, and Torrentispora. Molecular data suggest that Spadicoides is polyphyletic. The core of the genus is positioned in the Xenospadicoidales; Spadicoides s. str. is experimentally linked with sexual morphs for the first time. Based on DNA sequence data, the monotypic genera Xenospadicoides and Pseudodiplococcium are reduced to synonymy under Spadicoides, while Fusoidispora and Pseudoannulatascus are synonymised with Torrentispora. Members of the Xenospadicoidales inhabit decaying wood in terrestrial and freshwater environments and share a few morphological characters such as the absence of stromatic tissue, ascomata with a cylindrical or rostrate neck, similar anatomies of the ascomatal walls, thin-walled unitunicate asci with a non-amyloid apical annulus, disintegrating paraphyses, usually ellipsoidal to fusiform ascospores and holoblastic-denticulate or tretic conidiogenesis. Revised Ceratostomella spp. with persistent asci are listed and the taxonomic status of each species is re-evaluated based on revision of the holotype and other representative material, published details and available phylogenetic data.

Identification and characterization of suppressors of plant cell death (SPD) effectors from Magnaporthe oryzae

Phytopathogenic microorganisms, including the fungal pathogen Magnaporthe oryzae, secrete a myriad of effector proteins to facilitate infection. Utilizing the transient expression of candidate effectors in the leaves of the model plant Nicotiana benthamiana, we identified 11 suppressors of plant cell death (SPD) effectors from M. oryzae that were able to block the host cell death reaction induced by Nep1. Ten of these 11 were also able to suppress BAX-mediated plant cell death. Five of the 11 SPD genes have been identified previously as either essential for the pathogenicity of M. oryzae, secreted into the plant during disease development, or as suppressors or homologues of other characterized suppressors. In addition, of the remaining six, we showed that SPD8 (previously identified as BAS162) was localized to the rice cytoplasm in invaded and surrounding uninvaded cells during biotrophic invasion. Sequence analysis of the 11 SPD genes across 43 re-sequenced M. oryzae genomes revealed that SPD2, SPD4 and SPD7 have nucleotide polymorphisms amongst the isolates. SPD4 exhibited the highest level of nucleotide diversity of any currently known effector from M. oryzae in addition to the presence/absence polymorphisms, suggesting that this gene is potentially undergoing selection to avoid recognition by the host. Taken together, we have identified a series of effectors, some of which were previously unknown or whose function was unknown, that probably act at different stages of the infection process and contribute to the virulence of M. oryzae.

Small RNA Functions Are Required for Growth and Development of Magnaporthe oryzae

RNA interference (RNAi) is conserved in eukaryotic organisms, and it has been well studied in many animal and plant species and some fungal species, yet it is not well studied in fungal plant pathogens. In the rice blast fungus Magnaporthe oryzae, we examined small RNA (sRNA) and their biogenesis in the context of growth and pathogenicity. Through genetic and genomic analyses, we demonstrate that loss of a single gene encoding Dicer, RNA-dependent RNA polymerase, or Argonaute reduces sRNA levels. These three proteins are required for the biogenesis of sRNA-matching genome-wide regions (coding regions, repeats, and intergenic regions). The loss of one Argonaute reduced both sRNA and fungal virulence on barley leaves. Transcriptome analysis of multiple mutants revealed that sRNA play an important role in transcriptional regulation of repeats and intergenic regions in M. oryzae. Together, these data support that M. oryzae sRNA regulate developmental processes including, fungal growth and virulence.

A Matter of Scale and Dimensions: Chromatin of Chromosome Landmarks in the Fungi

Chromatin and chromosomes of fungi are highly diverse and dynamic, even within species. Much of what we know about histone modification enzymes, RNA interference, DNA methylation, and cell cycle control was first addressed in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, and Neurospora crassa. Here, we examine the three landmark regions that are required for maintenance of stable chromosomes and their faithful inheritance, namely, origins of DNA replication, telomeres and centromeres. We summarize the state of recent chromatin research that explains what is required for normal function of these specialized chromosomal regions in different fungi, with an emphasis on the silencing mechanism associated with subtelomeric regions, initiated by sirtuin histone deacetylases and histone H3 lysine 27 (H3K27) methyltransferases. We explore mechanisms for the appearance of "accessory" or "conditionally dispensable" chromosomes and contrast what has been learned from studies on genome-wide chromosome conformation capture in S. cerevisiae, S. pombe, N. crassa, and Trichoderma reesei. While most of the current knowledge is based on work in a handful of genetically and biochemically tractable model organisms, we suggest where major knowledge gaps remain to be closed. Fungi will continue to serve as facile organisms to uncover the basic processes of life because they make excellent model organisms for genetics, biochemistry, cell biology, and evolutionary biology.

Directional Selection from Host Plants Is a Major Force Driving Host Specificity in Magnaporthe Species

One major threat to global food security that requires immediate attention, is the increasing incidence of host shift and host expansion in growing number of pathogenic fungi and emergence of new pathogens. The threat is more alarming because, yield quality and quantity improvement efforts are encouraging the cultivation of uniform plants with low genetic diversity that are increasingly susceptible to emerging pathogens. However, the influence of host genome differentiation on pathogen genome differentiation and its contribution to emergence and adaptability is still obscure. Here, we compared genome sequence of 6 isolates of Magnaporthe species obtained from three different host plants. We demonstrated the evolutionary relationship between Magnaporthe species and the influence of host differentiation on pathogens. Phylogenetic analysis showed that evolution of pathogen directly corresponds with host divergence, suggesting that host-pathogen interaction has led to co-evolution. Furthermore, we identified an asymmetric selection pressure on Magnaporthe species. Oryza sativa-infecting isolates showed higher directional selection from host and subsequently tends to lower the genetic diversity in its genome. We concluded that, frequent gene loss or gain, new transposon acquisition and sequence divergence are host adaptability mechanisms for Magnaporthe species, and this coevolution processes is greatly driven by directional selection from host plants.

Subtelomere organization in the genome of the microsporidian Encephalitozoon cuniculi: patterns of repeated sequences and physicochemical signatures

Background

The microsporidian Encephalitozoon cuniculi is an obligate intracellular eukaryotic pathogen with a small nuclear genome (2.9 Mbp) consisting of 11 chromosomes. Although each chromosome end is known to contain a single rDNA unit, the incomplete assembly of subtelomeric regions following sequencing of the genome identified only 3 of the 22 expected rDNA units. While chromosome end assembly remains a difficult process in most eukaryotic genomes, it is of significant importance for pathogens because these regions encode factors important for virulence and host evasion.

Results

Here we report the first complete assembly of E. cuniculi chromosome ends, and describe a novel mosaic structure of segmental duplications (EXT repeats) in these regions. EXT repeats range in size between 3.5 and 23.8 kbp and contain four multigene families encoding membrane associated proteins. Twenty-one recombination sites were identified in the sub-terminal region of E. cuniculi chromosomes. Our analysis suggests that these sites contribute to the diversity of chromosome ends organization through Double Strand Break repair mechanisms. The region containing EXT repeats at chromosome extremities can be differentiated based on gene composition, GC content, recombination sites density and chromosome landscape.

Conclusion

Together this study provides the complete structure of the chromosome ends of E. cuniculi GB-M1, and identifies important factors, which could play a major role in parasite diversity and host-parasite interactions. Comparison with other eukaryotic genomes suggests that terminal regions could be distinguished precisely based on gene content, genetic instability and base composition biais. The diversity of processes assciated with chromosome extremities and their biological consequences, as they are presented in the present study, emphasize the fact that great effort will be necessary in the future to characterize more carefully these regions during whole genome sequencing efforts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1920-7) contains supplementary material, which is available to authorized users.

A complete annotation of the chromosomes of the cellulase producer Trichoderma reesei provides insights in gene clusters, their expression and reveals genes required for fitness

BACKGROUND

Investigations on a few eukaryotic model organisms showed that many genes are non-randomly distributed on chromosomes. In addition, chromosome ends frequently possess genes that are important for the fitness of the organisms. Trichoderma reesei is an industrial producer of enzymes for food, feed and biorefinery production. Its seven chromosomes have recently been assembled, thus making an investigation of its chromosome architecture possible.

RESULTS

We manually annotated and mapped 9194 ORFs on their respective chromosomes and investigated the clustering of the major gene categories and of genes encoding carbohydrate-active enzymes (CAZymes), and the relationship between clustering and expression. Genes responsible for RNA processing and modification, amino acid metabolism, transcription, translation and ribosomal structure and biogenesis indeed showed loose clustering, but this had no impact on their expression. A third of the genes encoding CAZymes also occurred in loose clusters that also contained a high number of genes encoding small secreted cysteine-rich proteins. Five CAZyme clusters were located less than 50 kb apart from the chromosome ends. These genes exhibited the lowest basal (but not induced) expression level, which correlated with an enrichment of H3K9 methylation in the terminal 50 kb areas indicating gene silencing. No differences were found in the expression of CAZyme genes present in other parts of the chromosomes. The putative subtelomeric areas were also enriched in genes encoding secreted proteases, amino acid permeases, enzyme clusters for polyketide synthases (PKS)-non-ribosomal peptide synthase (NRPS) fusion proteins (PKS-NRPS) and proteins involved in iron scavenging. They were strongly upregulated during conidiation and interaction with other fungi.

CONCLUSIONS

Our findings suggest that gene clustering on the T. reesei chromosomes occurs but generally has no impact on their expression. CAZyme genes, located in subtelomers, however, exhibited a much lower basal expression level. The gene inventory of the subtelomers suggests a major role of competition for nitrogen and iron supported by antibiosis for the fitness of T. reesei. The availability of fully annotated chromosomes will facilitate the use of genetic crossings in identifying still unknown genes responsible for specific traits of T. reesei.

Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes

BACKGROUND

Supernumerary chromosomes have been found in many organisms. In fungi, these "accessory" or "dispensable" chromosomes are present at different frequencies in populations and are usually characterized by higher repetitive DNA content and lower gene density when compared to the core chromosomes. In the reference strain of the wheat pathogen, Zymoseptoria tritici, eight discrete accessory chromosomes have been found. So far, no functional role has been assigned to these chromosomes; however, they have existed as separate entities in the karyotypes of Zymoseptoria species over evolutionary time. In this study, we addressed what-if anything-distinguishes the chromatin of accessory chromosomes from core chromosomes. We used chromatin immunoprecipitation combined with high-throughput sequencing ("ChIP-seq") of DNA associated with the centromere-specific histone H3, CENP-A (CenH3), to identify centromeric DNA, and ChIP-seq with antibodies against dimethylated H3K4, trimethylated H3K9 and trimethylated H3K27 to determine the relative distribution and proportion of euchromatin, obligate and facultative heterochromatin, respectively.

RESULTS

Centromeres of the eight accessory chromosomes have the same sequence composition and structure as centromeres of the 13 core chromosomes and they are of similar length. Unlike those of most other fungi, Z. tritici centromeres are not composed entirely of repetitive DNA; some centromeres contain only unique DNA sequences, and bona fide expressed genes are located in regions enriched with CenH3. By fluorescence microscopy, we showed that centromeres of Z. tritici do not cluster into a single chromocenter during interphase. We found dramatically higher enrichment of H3K9me3 and H3K27me3 on the accessory chromosomes, consistent with the twofold higher proportion of repetitive DNA and poorly transcribed genes. In contrast, no single histone modification tested here correlated with the distribution of centromeric nucleosomes.

CONCLUSIONS

All centromeres are similar in length and composed of a mixture of unique and repeat DNA, and most contain actively transcribed genes. Centromeres, subtelomeric regions or telomere repeat length cannot account for the differences in transfer fidelity between core and accessory chromosomes, but accessory chromosomes are greatly enriched in nucleosomes with H3K27 trimethylation. Genes on accessory chromosomes appear to be silenced by trimethylation of H3K9 and H3K27.

Genome Sequences of Three Phytopathogenic Species of the Magnaporthaceae Family of Fungi

Magnaporthaceae is a family of ascomycetes that includes three fungi of great economic importance: Magnaporthe oryzae, Gaeumannomyces graminis var. tritici, and Magnaporthe poae. These three fungi cause widespread disease and loss in cereal and grass crops, including rice blast disease (M. oryzae), take-all disease in wheat and other grasses (G. graminis), and summer patch disease in turf grasses (M. poae). Here, we present the finished genome sequence for M. oryzae and draft sequences for M. poae and G. graminis var. tritici. We used multiple technologies to sequence and annotate the genomes of M. oryzae, M. poae, and G. graminis var. tritici. The M. oryzae genome is now finished to seven chromosomes whereas M. poae and G. graminis var. tritici are sequenced to 40.0× and 25.0× coverage respectively. Gene models were developed by the use of multiple computational techniques and further supported by RNAseq data. In addition, we performed preliminary analysis of genome architecture and repetitive element DNA.

Complete Genome Sequence of Sporisorium scitamineum and Biotrophic Interaction Transcriptome with Sugarcane

Sporisorium scitamineum is a biotrophic fungus responsible for the sugarcane smut, a worldwide spread disease. This study provides the complete sequence of individual chromosomes of S. scitamineum from telomere to telomere achieved by a combination of PacBio long reads and Illumina short reads sequence data, as well as a draft sequence of a second fungal strain. Comparative analysis to previous available sequences of another strain detected few polymorphisms among the three genomes. The novel complete sequence described herein allowed us to identify and annotate extended subtelomeric regions, repetitive elements and the mitochondrial DNA sequence. The genome comprises 19,979,571 bases, 6,677 genes encoding proteins, 111 tRNAs and 3 assembled copies of rDNA, out of our estimated number of copies as 130. Chromosomal reorganizations were detected when comparing to sequences of S. reilianum, the closest smut relative, potentially influenced by repeats of transposable elements. Repetitive elements may have also directed the linkage of the two mating-type loci. The fungal transcriptome profiling from in vitro and from interaction with sugarcane at two time points (early infection and whip emergence) revealed that 13.5% of the genes were differentially expressed in planta and particular to each developmental stage. Among them are plant cell wall degrading enzymes, proteases, lipases, chitin modification and lignin degradation enzymes, sugar transporters and transcriptional factors. The fungus also modulates transcription of genes related to surviving against reactive oxygen species and other toxic metabolites produced by the plant. Previously described effectors in smut/plant interactions were detected but some new candidates are proposed. Ten genomic islands harboring some of the candidate genes unique to S. scitamineum were expressed only in planta. RNAseq data was also used to reassure gene predictions.

Real-Time Evolution of a Subtelomeric Gene Family in Candida albicans

Subtelomeric regions of the genome are notable for high rates of sequence evolution and rapid gene turnover. Evidence of subtelomeric evolution has relied heavily on comparisons of historical evolutionary patterns to infer trends and frequencies of these events. Here, we describe evolution of the subtelomeric TLO gene family in Candida albicans during laboratory passaging for over 4000 generations. C. albicans is a commensal and opportunistic pathogen of humans and the TLO gene family encodes a subunit of the Mediator complex that regulates transcription and affects a range of virulence factors. We identified 16 distinct subtelomeric recombination events that altered the TLO repertoire. Ectopic recombination between subtelomeres on different chromosome ends occurred approximately once per 5000 generations and was often followed by loss of heterozygosity, resulting in the complete loss of one TLO gene sequence with expansion of another. In one case, recombination within TLO genes produced a novel TLO gene sequence. TLO copy number changes were biased, with some TLOs preferentially being copied to novel chromosome arms and other TLO genes being frequently lost. The majority of these nonreciprocal recombination events occurred either within the 3′ end of the TLO coding sequence or within a conserved 50-bp sequence element centromere-proximal to TLO coding sequence. Thus, subtelomeric recombination is a rapid mechanism of generating genotypic diversity through alterations in the number and sequence of related gene family members.

Motif-independent de novo detection of secondary metabolite gene clusters-toward identification from filamentous fungi

Secondary metabolites are produced mostly by clustered genes that are essential to their biosynthesis. The transcriptional expression of these genes is often cooperatively regulated by a transcription factor located inside or close to a cluster. Most of the secondary metabolism biosynthesis (SMB) gene clusters identified to date contain so-called core genes with distinctive sequence features, such as polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS). Recent efforts in sequencing fungal genomes have revealed far more SMB gene clusters than expected based on the number of core genes in the genomes. Several bioinformatics tools have been developed to survey SMB gene clusters using the sequence motif information of the core genes, including SMURF and antiSMASH. More recently, accompanied by the development of sequencing techniques allowing to obtain large-scale genomic and transcriptomic data, motif-independent prediction methods of SMB gene clusters, including MIDDAS-M, have been developed. Most these methods detect the clusters in which the genes are cooperatively regulated at transcriptional levels, thus allowing the identification of novel SMB gene clusters regardless of the presence of the core genes. Another type of the method, MIPS-CG, uses the characteristics of SMB genes, which are highly enriched in non-syntenic blocks (NSBs), enabling the prediction even without transcriptome data although the results have not been evaluated in detail. Considering that large portion of SMB gene clusters might be sufficiently expressed only in limited uncommon conditions, it seems that prediction of SMB gene clusters by bioinformatics and successive experimental validation is an only way to efficiently uncover hidden SMB gene clusters. Here, we describe and discuss possible novel approaches for the determination of SMB gene clusters that have not been identified using conventional methods.

Biodiversity of genes encoding anti-microbial traits within plant associated microbes

The plant is an attractive versatile home for diverse associated microbes. A subset of these microbes produces a diversity of anti-microbial natural products including polyketides, non-ribosomal peptides, terpenoids, heterocylic nitrogenous compounds, volatile compounds, bacteriocins, and lytic enzymes. In recent years, detailed molecular analysis has led to a better understanding of the underlying genetic mechanisms. New genomic and bioinformatic tools have permitted comparisons of orthologous genes between species, leading to predictions of the associated evolutionary mechanisms responsible for diversification at the genetic and corresponding biochemical levels. The purpose of this review is to describe the biodiversity of biosynthetic genes of plant-associated bacteria and fungi that encode selected examples of antimicrobial natural products. For each compound, the target pathogen and biochemical mode of action are described, in order to draw attention to the complexity of these phenomena. We review recent information of the underlying molecular diversity and draw lessons through comparative genomic analysis of the orthologous coding sequences (CDS). We conclude by discussing emerging themes and gaps, discuss the metabolic pathways in the context of the phylogeny and ecology of their microbial hosts, and discuss potential evolutionary mechanisms that led to the diversification of biosynthetic gene clusters.

Effectiveness and durability of the rice pi-ta gene in Yunnan province of China

Rice blast is one of the most damaging diseases of rice worldwide. In the present study, we analyzed DNA sequence variation of avirulence (AVR) genes of AVR-Pita1 in field isolates of Magnaporthe oryzae in order to understand the effectiveness of the resistance gene Pi-ta in China. Genomic DNA of 366 isolates of M. oryzae collected from Yunnan province of China were used for polymerase chain reaction (PCR) amplification to examine the existence of AVR-Pita1 using gene-specific PCR markers. Results of PCR products revealed that 218 isolates of M. oryzae carry AVR-Pita1. Among of them, 62.5, 56.3, 58.5, 46.7, 72.4, and 57.4% of M. oryzae carry AVR-Pita1 from northeastern, southeast, western, northwest, southwestern, and central Yunnan province, respectively. The detection rate of AVR-Pita1 was, in order: southwestern > northeastern > western > central > southeastern > northwestern Yunnan province. Moreover, in total, 18 AVR-Pita1 haplotypes encoding 13 novel AVR-Pita1 variants were identified among 60 isolates. Most DNA sequence variation was found to occur in the exon region, resulting in amino acid substitution. Six virulent haplotypes of AVR-Pita1 to Pita were identified among 60 field isolates. The AVR-Pita1 has evolved to virulence from avirulent origins via base substitution. These findings demonstrate that AVR-Pita1 is under positive selection and mutations of AVR-Pita1 are responsible for defeating race-specific resistance in nature.

Novel telomere-anchored PCR approach for studying sexual stage telomeres in Aspergillus nidulans

Telomere length varies between germline and somatic cells of the same organism, leading to the hypothesis that telomeres are lengthened during meiosis. However, little is known about the meiotic telomere length in many organisms. In the filamentous fungus Aspergillus nidulans, the telomere lengths in hyphae and asexual spores are invariant. No study using existing techniques has determined the telomere length of the sexual ascospores due to the relatively low abundance of pure meiotic cells in A. nidulans and the small quantity of DNA present. To address this, we developed a simple and sensitive PCR strategy to measure the telomere length of A. nidulans meiotic cells. This novel technique, termed “telomere-anchored PCR,” measures the length of the telomere on chromosome II-L using a small fraction of the DNA required for the traditional terminal restriction fragment (TRF) Southern analysis. Using this approach, we determined that the A. nidulans ascospore telomere length is virtually identical to telomeres of other cell types from this organism, approximately 110 bp, indicating that a surprisingly strict telomere length regulation exists in the major cell types of A. nidulans. When the hyphal telomeres were measured in a telomerase reverse transcriptase (TERT) knockout strain, small decreases in length were readily detected. Thus, this technique can detect telomeres in relatively rare cell types and is particularly sensitive in measuring exceptionally short telomeres. This rapid and inexpensive telomere-anchored PCR method potentially can be utilized in other filamentous fungi and types of organisms.

Conservation, Divergence, and Genome-Wide Distribution of PAL and POX A Gene Families in Plants

Genome-wide identification and phylogenetic and syntenic comparison were performed for the genes responsible for phenylalanine ammonia lyase (PAL) and peroxidase A (POX A) enzymes in nine plant species representing very diverse groups like legumes (Glycine max and Medicago truncatula), fruits (Vitis vinifera), cereals (Sorghum bicolor, Zea mays, and Oryza sativa), trees (Populus trichocarpa), and model dicot (Arabidopsis thaliana) and monocot (Brachypodium distachyon) species. A total of 87 and 1045 genes in PAL and POX A gene families, respectively, have been identified in these species. The phylogenetic and syntenic comparison along with motif distributions shows a high degree of conservation of PAL genes, suggesting that these genes may predate monocot/eudicot divergence. The POX A family genes, present in clusters at the subtelomeric regions of chromosomes, might be evolving and expanding with higher rate than the PAL gene family. Our analysis showed that during the expansion of POX A gene family, many groups and subgroups have evolved, resulting in a high level of functional divergence among monocots and dicots. These results will act as a first step toward the understanding of monocot/eudicot evolution and functional characterization of these gene families in the future.

Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice worldwide. The fungal pathogen is notorious for its ability to overcome host resistance. To better understand its genetic variation in nature, we sequenced the genomes of two field isolates, Y34 and P131. In comparison with the previously sequenced laboratory strain 70-15, both field isolates had a similar genome size but slightly more genes. Sequences from the field isolates were used to improve genome assembly and gene prediction of 70-15. Although the overall genome structure is similar, a number of gene families that are likely involved in plant-fungal interactions are expanded in the field isolates. Genome-wide analysis on asynonymous to synonymous nucleotide substitution rates revealed that many infection-related genes underwent diversifying selection. The field isolates also have hundreds of isolate-specific genes and a number of isolate-specific gene duplication events. Functional characterization of randomly selected isolate-specific genes revealed that they play diverse roles, some of which affect virulence. Furthermore, each genome contains thousands of loci of transposon-like elements, but less than 30% of them are conserved among different isolates, suggesting active transposition events in M. oryzae. A total of approximately 200 genes were disrupted in these three strains by transposable elements. Interestingly, transposon-like elements tend to be associated with isolate-specific or duplicated sequences. Overall, our results indicate that gain or loss of unique genes, DNA duplication, gene family expansion, and frequent translocation of transposon-like elements are important factors in genome variation of the rice blast fungus.

Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae is one of the most destructive diseases of rice worldwide. The fungal pathogen is notorious for its ability to overcome host resistance. To better understand its genetic variation in nature, we sequenced the genomes of two field isolates, Y34 and P131. In comparison with the previously sequenced laboratory strain 70-15, both field isolates had a similar genome size but slightly more genes. Sequences from the field isolates were used to improve genome assembly and gene prediction of 70-15. Although the overall genome structure is similar, a number of gene families that are likely involved in plant-fungal interactions are expanded in the field isolates. Genome-wide analysis on asynonymous to synonymous nucleotide substitution rates revealed that many infection-related genes underwent diversifying selection. The field isolates also have hundreds of isolate-specific genes and a number of isolate-specific gene duplication events. Functional characterization of randomly selected isolate-specific genes revealed that they play diverse roles, some of which affect virulence. Furthermore, each genome contains thousands of loci of transposon-like elements, but less than 30% of them are conserved among different isolates, suggesting active transposition events in M. oryzae. A total of approximately 200 genes were disrupted in these three strains by transposable elements. Interestingly, transposon-like elements tend to be associated with isolate-specific or duplicated sequences. Overall, our results indicate that gain or loss of unique genes, DNA duplication, gene family expansion, and frequent translocation of transposon-like elements are important factors in genome variation of the rice blast fungus.

Magnaporthe oryzae is the causal agent of rice blast that is mainly controlled with resistance cultivars. However, genetic variations in the pathogen often lead to overcoming R gene-mediated resistance in rice cultivars. In this study we sequenced two field isolates from China and Japan. In comparison with the laboratory strain that was previously sequenced, the field isolates have a similar genome size and overall genome structure. However, they have slightly more genes and contain a number of expanded gene families that are likely involved in plant-fungal interactions. Each of the isolates has specific genes, some of which affect virulence and some others are important for asexual development. The three strains differ noticeably in the distribution of transposon-like elements. Many of the transposable elements tend to be associated with isolate-specific or duplicated sequences. This study revealed genetic factors involved in genome variation of the rice blast fungus.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

Microevolution of Cryptococcus neoformans driven by massive tandem gene amplification

The subtelomeric regions of organisms ranging from protists to fungi undergo a much higher rate of rearrangement than is observed in the rest of the genome. While characterizing these ~40-kb regions of the human fungal pathogen Cryptococcus neoformans, we have identified a recent gene amplification event near the right telomere of chromosome 3 that involves a gene encoding an arsenite efflux transporter (ARR3). The 3,177-bp amplicon exists in a tandem array of 2-15 copies and is present exclusively in strains with the C. neoformans var. grubii subclade VNI A5 MLST profile. Strains bearing the amplification display dramatically enhanced resistance to arsenite that correlates with the copy number of the repeat; the origin of increased resistance was verified as transport-related by functional complementation of an arsenite transporter mutant of Saccharomyces cerevisiae. Subsequent experimental evolution in the presence of increasing concentrations of arsenite yielded highly resistant strains with the ARR3 amplicon further amplified to over 50 copies, accounting for up to ~1% of the whole genome and making the copy number of this repeat as high as that seen for the ribosomal DNA. The example described here therefore represents a rare evolutionary intermediate-an array that is currently in a state of dynamic flux, in dramatic contrast to relatively common, static relics of past tandem duplications that are unable to further amplify due to nucleotide divergence. Beyond identifying and engineering fungal isolates that are highly resistant to arsenite and describing the first reported instance of microevolution via massive gene amplification in C. neoformans, these results suggest that adaptation through gene amplification may be an important mechanism that C. neoformans employs in response to environmental stresses, perhaps including those encountered during infection. More importantly, the ARR3 array will serve as an ideal model for further molecular genetic analyses of how tandem gene duplications arise and expand.

Correlation of gene expression and protein production rate - a system wide study

BACKGROUND

Growth rate is a major determinant of intracellular function. However its effects can only be properly dissected with technically demanding chemostat cultivations in which it can be controlled. Recent work on Saccharomyces cerevisiae chemostat cultivations provided the first analysis on genome wide effects of growth rate. In this work we study the filamentous fungus Trichoderma reesei (Hypocrea jecorina) that is an industrial protein production host known for its exceptional protein secretion capability. Interestingly, it exhibits a low growth rate protein production phenotype.

RESULTS

We have used transcriptomics and proteomics to study the effect of growth rate and cell density on protein production in chemostat cultivations of T. reesei. Use of chemostat allowed control of growth rate and exact estimation of the extracellular specific protein production rate (SPPR). We find that major biosynthetic activities are all negatively correlated with SPPR. We also find that expression of many genes of secreted proteins and secondary metabolism, as well as various lineage specific, mostly unknown genes are positively correlated with SPPR. Finally, we enumerate possible regulators and regulatory mechanisms, arising from the data, for this response.

CONCLUSIONS

Based on these results it appears that in low growth rate protein production energy is very efficiently used primarly for protein production. Also, we propose that flux through early glycolysis or the TCA cycle is a more fundamental determining factor than growth rate for low growth rate protein production and we propose a novel eukaryotic response to this i.e. the lineage specific response (LSR).

Multiple translocation of the AVR-Pita effector gene among chromosomes of the rice blast fungus Magnaporthe oryzae and related species

Magnaporthe oryzae is the causal agent of rice blast disease, a devastating problem worldwide. This fungus has caused breakdown of resistance conferred by newly developed commercial cultivars. To address how the rice blast fungus adapts itself to new resistance genes so quickly, we examined chromosomal locations of AVR-Pita, a subtelomeric gene family corresponding to the Pita resistance gene, in various isolates of M. oryzae (including wheat and millet pathogens) and its related species. We found that AVR-Pita (AVR-Pita1 and AVR-Pita2) is highly variable in its genome location, occurring in chromosomes 1, 3, 4, 5, 6, 7, and supernumerary chromosomes, particularly in rice-infecting isolates. When expressed in M. oryzae, most of the AVR-Pita homologs could elicit Pita-mediated resistance, even those from non-rice isolates. AVR-Pita was flanked by a retrotransposon, which presumably contributed to its multiple translocation across the genome. On the other hand, family member AVR-Pita3, which lacks avirulence activity, was stably located on chromosome 7 in a vast majority of isolates. These results suggest that the diversification in genome location of AVR-Pita in the rice isolates is a consequence of recognition by Pita in rice. We propose a model that the multiple translocation of AVR-Pita may be associated with its frequent loss and recovery mediated by its transfer among individuals in asexual populations. This model implies that the high mobility of AVR-Pita is a key mechanism accounting for the rapid adaptation toward Pita. Dynamic adaptation of some fungal plant pathogens may be achieved by deletion and recovery of avirulence genes using a population as a unit of adaptation.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

Targeted cloning of fungal telomeres

Telomeres are the sequences that form the ends of eukaryotic chromosomes and are essential structures that confer genome stability and guide chromosome behavior. In addition, the terminal regions of the chromosomes tend to house genes with predicted roles in ecological adaptation. Unfortunately, however, most fungal genome assemblies contain very few telomeres and, therefore, the identities of genes residing near the chromosome ends are often unknown. In an effort to develop a complete understanding of the organization and gene content of chromosome ends in a number of fungi, we developed efficient methods for the identification and targeted cloning of telomeres. This chapter describes the basic steps and shows exemplary results from the targeted cloning of Epichloë festucae telomeres.

Telomere position effect is regulated by heterochromatin-associated proteins and NkuA in Aspergillus nidulans

Gene-silencing mechanisms are being shown to be associated with an increasing number of fungal developmental processes. Telomere position effect (TPE) is a eukaryotic phenomenon resulting in gene repression in areas immediately adjacent to telomere caps. Here, TPE is shown to regulate expression of transgenes on the left arm of chromosome III and the right arm of chromosome VI in Aspergillus nidulans. Phenotypes found to be associated with transgene repression included reduction in radial growth and the absence of sexual spores; however, these pleiotropic phenotypes were remedied when cultures were grown on media with appropriate supplementation. Simple radial growth and ascosporogenesis assays provided insights into the mechanism of TPE, including a means to determine its extent. These experiments revealed that the KU70 homologue (NkuA) and the heterochromatin-associated proteins HepA, ClrD and HdaA were partially required for transgene silencing. This study indicates that TPE extends at least 30 kb on chromosome III, suggesting that this phenomenon may be important for gene regulation in subtelomeric regions of A. nidulans.

Detecting novel genes with sparse arrays

Species-specific genes play an important role in defining the phenotype of an organism. However, current gene prediction methods can only efficiently find genes that share features such as sequence similarity or general sequence characteristics with previously known genes. Novel sequencing methods and tiling arrays can be used to find genes without prior information and they have demonstrated that novel genes can still be found from extensively studied model organisms. Unfortunately, these methods are expensive and thus are not easily applicable, e.g., to finding genes that are expressed only in very specific conditions. We demonstrate a method for finding novel genes with sparse arrays, applying it on the 33.9 Mb genome of the filamentous fungus Trichoderma reesei. Our computational method does not require normalisations between arrays and it takes into account the multiple-testing problem typical for analysis of microarray data. In contrast to tiling arrays, that use overlapping probes, only one 25 mer microarray oligonucleotide probe was used for every 100b. Thus, only relatively little space on a microarray slide was required to cover the intergenic regions of a genome. The analysis was done as a by-product of a conventional microarray experiment with no additional costs. We found at least 23 good candidates for novel transcripts that could code for proteins and all of which were expressed at high levels. Candidate genes were found to neighbour ire1 and cre1 and many other regulatory genes. Our simple, low-cost method can easily be applied to finding novel species-specific genes without prior knowledge of their sequence properties.

Array comparative genomic hybridization analysis of Trichoderma reesei strains with enhanced cellulase production properties

Background

Trichoderma reesei is the main industrial producer of cellulases and hemicellulases that are used to depolymerize biomass in a variety of biotechnical applications. Many of the production strains currently in use have been generated by classical mutagenesis. In this study we characterized genomic alterations in high-producing mutants of T. reesei by high-resolution array comparative genomic hybridization (aCGH). Our aim was to obtain genome-wide information which could be utilized for better understanding of the mechanisms underlying efficient cellulase production, and would enable targeted genetic engineering for improved production of proteins in general.

Results

We carried out an aCGH analysis of four high-producing strains (QM9123, QM9414, NG14 and Rut-C30) using the natural isolate QM6a as a reference. In QM9123 and QM9414 we detected a total of 44 previously undocumented mutation sites including deletions, chromosomal translocation breakpoints and single nucleotide mutations. In NG14 and Rut-C30 we detected 126 mutations of which 17 were new mutations not documented previously. Among these new mutations are the first chromosomal translocation breakpoints identified in NG14 and Rut-C30. We studied the effects of two deletions identified in Rut-C30 (a deletion of 85 kb in the scaffold 15 and a deletion in a gene encoding a transcription factor) on cellulase production by constructing knock-out strains in the QM6a background. Neither the 85 kb deletion nor the deletion of the transcription factor affected cellulase production.

Conclusions

aCGH analysis identified dozens of mutations in each strain analyzed. The resolution was at the level of single nucleotide mutation. High-density aCGH is a powerful tool for genome-wide analysis of organisms with small genomes e.g. fungi, especially in studies where a large set of interesting strains is analyzed.

PWT1, an avirulence gene of Magnaporthe oryzae tightly linked to the rDNA Locus, is recognized by two staple crops, common wheat and barley

The pathogenicity to wheat (Pwt1) locus conditions host species specificity of Magnaporthe oryzae on wheat. GFSI1-7-2 (Setaria isolate) carries the avirulence allele (PWT1) at this locus while Br48 (Triticum isolate) carries the virulence allele (pwt1). An F(1) culture derived from a cross between GFSI1-7-2 and Br48 was backcrossed with Br48 to produce a tester population in which PWT1 alone segregated. When hexaploid wheat lines were inoculated with the BC(1)F(1) testers, they were all resistant to all PWT1 carriers and susceptible to all pwt1 carriers, suggesting that they recognize PWT1. When barley cultivars were inoculated with the testers, they showed the same pattern of reactions as the hexaploid lines, suggesting that the barley cultivars also recognize PWT1. These results suggest that PWT1 is a fundamental gene that universally conditions the avirulence of Setaria isolates on two staple crops, hexaploid wheat and barley. Interestingly, tetraploid wheat lines did not recognize PWT1. Molecular mapping using the F(1) and BC(1)F(1) populations revealed that the Pwt1 locus is located on chromosome 2 and tightly linked to the ribosomal DNA locus and a telomere.

Systematic overrepresentation of DNA termini and underrepresentation of subterminal regions among sequencing templates prepared from hydrodynamically sheared linear DNA molecules

Background

Analysis of fungal genome sequence assemblies reveals that telomeres are poorly represented even though telomeric reads tend to be superabundant. We surmised that the problem might lie in the DNA shearing conditions used to create clone libraries for genome sequencing.

Results

A shotgun strategy was used to sequence and assemble circular and linear cosmid DNAs sheared using conditions typical for a genome project. The DNA sheared in circular form assembled into a single sequence contig. However, the linearized cosmid produced an incomplete assembly because the two DNA termini, though greatly overrepresented in the clone library used for sequencing, were separated from neighboring sequences by gaps of ~1.4 and 1.8 kb. These gap sizes were reduced, but not eliminated, by shearing the linear cosmid into smaller fragments. Mapping of shearing breakpoints revealed a paucity of breaks in the subterminal regions of the linearized cosmid and also near chromosome ends of the fungus Neurospora crassa.

Conclusion

Together, our data indicate that the ends of linear DNA molecules are recalcitrant to hydrodynamic shearing. We propose that this causes DNA termini to be overrepresented in the resulting fragment population but ultimately prevents their incorporation into sequence assemblies.

The ergot alkaloid gene cluster: functional analyses and evolutionary aspects

Ergot alkaloids and their derivatives have been traditionally used as therapeutic agents in migraine, blood pressure regulation and help in childbirth and abortion. Their production in submerse culture is a long established biotechnological process. Ergot alkaloids are produced mainly by members of the genus Claviceps, with Claviceps purpurea as best investigated species concerning the biochemistry of ergot alkaloid synthesis (EAS). Genes encoding enzymes involved in EAS have been shown to be clustered; functional analyses of EAS cluster genes have allowed to assign specific functions to several gene products. Various Claviceps species differ with respect to their host specificity and their alkaloid content; comparison of the ergot alkaloid clusters in these species (and of clavine alkaloid clusters in other genera) yields interesting insights into the evolution of cluster structure. This review focuses on recently published and also yet unpublished data on the structure and evolution of the EAS gene cluster and on the function and regulation of cluster genes. These analyses have also significant biotechnological implications: the characterization of non-ribosomal peptide synthetases (NRPS) involved in the synthesis of the peptide moiety of ergopeptines opened interesting perspectives for the synthesis of ergot alkaloids; on the other hand, defined mutants could be generated producing interesting intermediates or only single peptide alkaloids (instead of the alkaloid mixtures usually produced by industrial strains).

The choanoflagellate Monosiga brevicollis karyotype revealed by the genome sequence: telomere-linked helicase genes resemble those of some fungi

The approximately 42 Mbp assembled genome sequence for the choanoflagellate Monosiga brevicollis reveals that most of the large scaffolds of 300-2,600 kb represent entire chromosomes or chromosome arms. Telomeres are partially assembled at the termini of 37 scaffolds, while another 43 scaffolds end in telomere-associated regions containing distinctive gene sets. Potential centromeric regions were identified on 39 scaffolds. Together, these observations suggest a karyotype of approximately 40 metacentric and submetacentric chromosomes averaging 1 Mbp in size. Genes encoding RecQ family DNA helicases, along with ankyrin-domain proteins and serine/threonine kinases, are associated with most telomeres, a feature shared with some fungi. This telomere-linked helicase gene arrangement might be ancestral to both fungi and choanoflagellates in the super-kingdom Opisthokonta; however, the great lability of telomere architecture suggests that it could also be a convergent feature.