Organization and distribution pattern of MGLR-3, a novel retrotransposon in the rice blast fungus Magnaporthe grisea

Comparative Genomics Reveals the High Copy Number Variation of a Retro Transposon in Different Magnaporthe Isolates

Magnaporthe oryzae is one of the fungal pathogens of rice which results in heavy yield losses worldwide. Understanding the genomic structure of M. oryzae is essential for appropriate deployment of the blast resistance in rice crop improvement programs. In this study we sequenced two M. oryzae isolates, RML-29 (avirulent) and RP-2421 (highly virulent) and performed comparative study along with three publically available genomes of 70-15, P131, and Y34. We identified several candidate effectors (>600) and isolate specific sequences from RML-29 and RP-2421, while a core set of 10013 single copy orthologs were found among the isolates. Pan-genome analysis showed extensive presence and absence variations (PAVs). We identified isolate-specific genes across 12 isolates using the pan-genome information. Repeat analysis was separately performed for each of the 15 isolates. This analysis revealed ∼25 times higher copy number of short interspersed nuclear elements (SINE) in virulent than avirulent isolate. We conclude that the extensive PAVs and occurrence of SINE throughout the genome could be one of the major mechanisms by which pathogenic variability is emerging in M. oryzae isolates. The knowledge gained in this comparative genome study can provide understandings about the fungal genome variations in different hosts and environmental conditions, and it will provide resources to effectively manage this important disease of rice.

Genome-Wide Comparison of Magnaporthe Species Reveals a Host-Specific Pattern of Secretory Proteins and Transposable Elements

Blast disease caused by the Magnaporthe species is a major factor affecting the productivity of rice, wheat and millets. This study was aimed at generating genomic information for rice and non-rice Magnaporthe isolates to understand the extent of genetic variation. We have sequenced the whole genome of the Magnaporthe isolates, infecting rice (leaf and neck), finger millet (leaf and neck), foxtail millet (leaf) and buffel grass (leaf). Rice and finger millet isolates infecting both leaf and neck tissues were sequenced, since the damage and yield loss caused due to neck blast is much higher as compared to leaf blast. The genome-wide comparison was carried out to study the variability in gene content, candidate effectors, repeat element distribution, genes involved in carbohydrate metabolism and SNPs. The analysis of repeat element footprints revealed some genes such as naringenin, 2-oxoglutarate 3-dioxygenase being targeted by Pot2 and Occan, in isolates from different host species. Some repeat insertions were host-specific while other insertions were randomly shared between isolates. The distributions of repeat elements, secretory proteins, CAZymes and SNPs showed significant variation across host-specific lineages of Magnaporthe indicating an independent genome evolution orchestrated by multiple genomic factors.

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.

Molecular cloning and characterization of the AVR-Pia locus from a Japanese field isolate of Magnaporthe oryzae

In order to clone and analyse the avirulence gene AVR-Pia from Japanese field isolates of Magnaporthe oryzae, a mutant of the M. oryzae strain Ina168 was isolated. This mutant, which was named Ina168m95-1, gained virulence towards the rice cultivar Aichi-asahi, which contains the resistance gene Pia. A DNA fragment (named PM01) that was deleted in the mutant and that co-segregated with avirulence towards Aichi-asahi was isolated. Three cosmid clones that included the regions that flanked PM01 were isolated from a genomic DNA library. One of these clones (46F3) complemented the mutant phenotype, which indicated clearly that this clone contained the avirulence gene AVR-Pia. Clone 46F3 contained insertions of transposable elements. The 46F3 insert was divided into fragments I-VI, and these were cloned individually into a hygromycin-resistant vector for the transformation of the mutant Ina168m95-1. An inoculation assay of the transformants revealed that fragment V (3.5 kb) contained AVR-Pia. By deletion analysis of fragment V, AVR-Pia was localized to an 1199-bp DNA fragment, which included a 255-bp open reading frame with weak homology to a bacterial cytochrome-c-like protein. Restriction fragment length polymorphism analysis of this region revealed that this DNA sequence co-segregated with the AVR-Pia locus in a genetic map that was constructed using Chinese isolates.

Characterization of chromosome ends in the filamentous fungus Neurospora crassa

Telomeres and subtelomere regions have vital roles in cellular homeostasis and can facilitate niche adaptation. However, information on telomere/subtelomere structure is still limited to a small number of organisms. Prior to initiation of this project, the Neurospora crassa genome assembly contained only 3 of the 14 telomeres. The missing telomeres were identified through bioinformatic mining of raw sequence data from the genome project and from clones in new cosmid and plasmid libraries. Their chromosomal locations were assigned on the basis of paired-end read information and/or by RFLP mapping. One telomere is attached to the ribosomal repeat array. The remaining chromosome ends have atypical structures in that they lack distinct subtelomere domains or other sequence features that are associated with telomeres in other organisms. Many of the chromosome ends terminate in highly AT-rich sequences that appear to be products of repeat-induced point mutation, although most are not currently repeated sequences. Several chromosome termini in the standard Oak Ridge wild-type strain were compared to their counterparts in an exotic wild type, Mauriceville. This revealed that the sequences immediately adjacent to the telomeres are usually genome specific. Finally, despite the absence of many features typically found in the telomere regions of other organisms, the Neurospora chromosome termini still retain the dynamic nature that is characteristic of chromosome ends.

Genome organization and evolution of the AVR-Pita avirulence gene family in the Magnaporthe grisea species complex

The avirulence (AVR) gene AVR-Pita in Magnaporthe oryzae prevents the fungus from infecting rice cultivars containing the resistance gene Pi-ta. A survey of isolates of the M. grisea species complex from diverse hosts showed that AVR-Pita is a member of a gene family, which led us to rename it to AVR-Pita1. Avirulence function, distribution, and genomic context of two other members, named AVR-Pita2 and AVR-Pita3, were characterized. AVR-Pita2, but not AVR-Pita3, was functional as an AVR gene corresponding to Pi-ta. The AVR-Pita1 and AVR-Pita2 genes were present in isolates of both M. oryzae and M. grisea, whereas the AVR-Pita3 gene was present only in isolates of M. oryzae. Orthologues of members of the AVR-Pita family could not be found in any fungal species sequenced to date, suggesting that the gene family may be unique to the M. grisea species complex. The genomic context of its members was analyzed in eight strains. The AVR-Pita1 and AVR-Pita2 genes in some isolates appeared to be located near telomeres and flanked by diverse repetitive DNA elements, suggesting that frequent deletion or amplification of these genes within the M. grisea species complex might have resulted from recombination mediated by repetitive DNA elements.

Genetics of avirulence genes in Blumeria graminis f.sp. hordei and physical mapping of AVR(a22) and AVR(a12)

Powdery mildew fungi are parasites that cause disease on a wide range of important crops. Plant resistance (R) genes, which induce host defences against powdery mildews, encode proteins that recognise avirulence (AVR) molecules from the parasite in a gene-for-gene manner. To gain insight into how virulence evolves in Blumeria graminis f.sp. hordei, associations between segregating AVR genes were established. As a prerequisite to the isolation of AVR genes, two loci were selected for further analysis. AVR(a22) is located in a tightly linked cluster comprising AVR(a10) and AVR(k1) as well as up to five other AVR genes. The ratio between physical and genetic distance in the cluster ranged between 0.7 and 35 kB/cM. The AVR(a22) locus was delimited by the previously isolated gene AVR(a10) and two cleaved amplified polymorphic sequence (CAPS) markers, 19H12R and 74E9L. By contrast, AVR(a12) was not linked to other AVR genes in two crosses. Bulk segregant analysis of over 100,000 AFLP fragments yielded two markers, ETAMTG-285 and PAAMACT-473, mapping 10 and 2cM from AVR(a12), respectively, thus delimiting AVR(a12) on one side. All markers obtained for AVR(a12) mapped proximal to it, indicating that the gene is located at the end of a chromosome. Three more AVR(a10) paralogues were identified at the locus interspersed among genes for metabolic enzymes and abundant repetitive elements, especially those homologous to the CgT1 class of retrotransposons. The flanking and close markers obtained will facilitate the isolation of AVR(a22) and AVR(a12) and provide useful tools for studies of the evolution of powdery mildew fungi in agriculture and nature.

Telomeres in the rice blast fungus Magnaporthe oryzae: the world of the end as we know it

The subtelomeres of many microbial eukaryotes are highly enriched in genes with roles in niche adaptation. Host and cultivar specificity genes in the rice blast fungus Magnaporthe oryzae also tend to be located near telomeres. In addition, the M. oryzae telomeres are highly variable chromosome regions. These observations suggested that plant pathogenic fungi might also use subtelomere regions for the amplification of genes with adaptive significance. Targeted sequencing of the M. oryzae telomeres provided an opportunity to test this hypothesis, and has yielded valuable insights into the organization and dynamics of these important chromosome regions.

Structure and transcriptional impact of divergent repetitive elements inserted within Phanerochaete chrysosporium strain RP-78 genes

We describe the structure, organization, and transcriptional impact of repetitive elements within the lignin-degrading basidiomycete, Phanerochaete chrysosporium. Searches of the P. chrysosporium genome revealed five copies of pce1, a ~1,750-nt non-autonomous, class II element. Alleles encoding a putative glucosyltransferase and a cytochrome P450 harbor pce insertions and produce incomplete transcripts. Class I elements included pcret1, an intact 8.14-kb gypsy-like retrotransposon inserted within a member of the multicopper oxidase gene family. Additionally, we describe a complex insertion of nested transposons within another putative cytochrome P450 gene. The disrupted allele lies within a cluster of >14 genes, all of which encode family 64 cytochrome P450s. Components of the insertion include a disjoint copia-like element, pcret3, the pol domain of a second retroelement, pcret2, and a duplication of an extended ORF of unknown function. As in the case of the pce elements, pcret1 and pcret2/3 insertions are confined to single alleles, transcripts of which are truncated. The corresponding wild-type alleles are apparently unaffected. In aggregate, P. chrysosporium harbors a complex array of repetitive elements, at least five of which directly influence expression of genes within families of structurally related sequences.

Transposition of MINE, a composite retrotransposon, in the avirulence gene ACE1 of the rice blast fungus Magnaporthe grisea

The ACE1 avirulence gene allele from the rice blast fungus Magnaporthe grisea was characterized in virulent isolate 2/0/3, revealing the insertion of a 1.9 kb MINE retrotransposon in the last ACE1 exon. MINE is a novel chimeric element composed of a transcribed non-coding sequence of 1.1 kb (WEIRD) fused to a 5'-truncated MGL retrotransposon. MINEs were found in high copy number in M. grisea isolates from rice (68 copies) and as a single copy in isolate CD156 from Eleusine. MINEs vary in size (1.3-6.7 kb) with conserved 5' WEIRD sequences and variable 3' MGL sequences. MGLs fused to WEIRDs correspond to different 5'-truncated MGLs with conserved 3' ends. The organization and diversity of MINEs suggest that these retrotransposons result from independent fusions between WEIRD and 5'-truncated MGLs. Such chimera could be formed during MGL reverse transcription as proposed for human U6-LINE1 chimeric retrotransposons and integrated into M. grisea genome using MGL machinery.

BAC end sequences and a physical map reveal transposable element content and clustering patterns in the genome of Magnaporthe grisea

Transposable elements (TEs) are viewed as major contributors to the evolution of fungal genomes. Genomic resources such as BAC libraries are an underutilized resource for studying genome-wide TE distribution. Using the BAC end sequences and physical map that are available for the rice blast fungus, Magnaporthe grisea, we describe a likelihood ratio test designed to identify clustering of TEs in the genome. A significant variation in the distribution of three TEs, MAGGY, MGL, and Pot2 was observed among the fingerprint contigs of the physical map. We utilized a draft sequence of M. grisea chromosome 7 to validate our results and found a similar pattern of clustering. By examining individual BAC end sequences, we found evidence for 11 unique integrations of MAGGY or MGL into Pot2 but no evidence for the reciprocal integration of Pot2 into another TE. This suggests that: (a) the presence of Pot2 in the genome predates that of the other TEs, (b) Pot2 was less transpositionally active than other TEs, or (c) that MAGGY and MGL have integration site preference for Pot2. High transition/transversion mutation ratios as well as bias in transition site context was observed in MAGGY and MGL elements, but not in Pot2 elements. These features are consistent with the effects of a Repeat-Induced Point (RIP) mutation-like process occurring in MAGGY and MGL elements. This study illustrates the general utility of a physical map and BAC end sequences for the study of genome-wide repetitive DNA content and organization.

Transposable elements in filamentous fungi

The past 10 years have been productive in the characterization of fungal transposable elements (TEs). All eukaryotic TEs described are found including an extraordinary prevalence of active members of the pogo family. The role of TEs in mutation and genome organization is well documented, leading to significant advances in our perception of the mechanisms underlying genetic changes in these organisms. TE-mediated changes, associated with transposition and recombination, provide a broad range of genetic variation, which is useful for natural populations in their adaptation to environmental constraints, especially for those lacking the sexual stage. Interestingly, some fungal species have evolved distinct silencing mechanisms that are regarded as host defense systems against TEs. The examination of forces acting on the evolutionary dynamics of TEs should provide important insights into the interactions between TEs and the fungal genome. Another issue of major significance is the practical applications of TEs in gene tagging and population analysis, which will undoubtedly facilitate research in systematic biology and functional genomics.

On the Trail of a Cereal Killer: Exploring the Biology ofMagnaporthe grisea

The blast fungus Magnaporthe grisea causes a serious disease on a wide variety of grasses including rice, wheat, and barley. Rice blast is the most serious disease of cultivated rice and therefore poses a threat to the world's most important food security crop. Here, I review recent progress toward understanding the molecular biology of plant infection by M. grisea, which involves development of a specialized cell, the appressorium. This dome-shaped cell generates enormous turgor pressure and physical force, allowing the fungus to breach the host cuticle and invade plant tissue. The review also considers the role of avirulence genes in M. grisea and the mechanisms by which resistant rice cultivars are able to perceive the fungus and defend themselves. Finally, the likely mechanisms that promote genetic diversity in M. grisea and our current understanding of the population structure of the blast fungus are evaluated.

Occan, a novel transposon in the Fot1 family, is ubiquitously found in several Magnaporthe grisea isolates

We investigated a DNA fragment and its flanking region deleted in the spontaneous Pi-a virulent mutant of Magnaporthe grisea Ina168. A new transposon-like sequence was identified from a region adjacent to the deleted fragment and was named Occan. Occan contained a 2,259-bp ORF interrupted by one 63-bp intron and had both a TA dinucleotide and 77 bp of perfect inverted repeats at both termini, without direct repeats. These features indicated that Occan is a member of the Fot1 family. RT-PCR analysis confirmed the expression of the putative transposase and the presence of an intron. Southern analysis of pulse-field gel electrophoresis-separated chromosomes indicated that Occan was dispersed in all chromosomes of the rice pathogen, Ina168. Copy numbers of Occan were also preserved in a host-specific manner amongst M. grisea isolates. In particular, rice pathogens contained a large number of the element inserted into their genome. Phylogenetic analysis with other known members of the Fot1 family revealed that Occan was dissimilar to any other known elements and it is thus proposed that Occan be separated to a new subfamily.

Impala, a transposon from Fusarium oxysporum, is active in the genome of Penicillium griseoroseum

An autonomous impala transposon trapped in Fusarium oxysporum by insertion within the niaD gene encoding nitrate reductase was introduced in the genome of the fungus Penicillium griseoroseum, a producer of pectinase enzymes. Through a phenotypic assay, we demonstrate that this element is able to excise from the niaD gene and to reinsert at new genomic positions. As in the original host, impala inserts into a TA site and footprints left by impala excisions are generally 5 bp. The fact that impala is able to transpose in P. griseoroseum offers the opportunity to develop a gene-tagging system based on this element with the objective to detect and clone genes related in pectinase production.

Evolution and Organization of a Highly Dynamic, Subtelomeric Helicase Gene Family in the Rice Blast Fungus Magnaporthe grisea

Sequence analysis of a 13-kb telomeric region in O-137, a rice pathogenic isolate of Magnaporthe grisea, uncovered a novel gene, designated TLH1 (telomere-linked helicase 1). The TLH1 gene is a member of a gene family, and the sequences flanking this gene family have also been amplified. Genetic mapping showed that most members of the TLH gene family are tightly linked to the telomeres. A physical mapping technique, termed RecA-mediated Achilles’ heel cleavage, and cloning and sequencing of two additional telomeres of O-137 associated with the TLH gene family confirmed that most members of the TLH gene family are located within 10 kb from the telomeric repeat. A survey of M. grisea strains from diverse hosts revealed that the gene family is ubiquitously present among rice pathogens, but is absent from almost all isolates of hosts other than rice. The gene family appears to be highly dynamic, undergoing frequent deletion/amplification events. Given the presence of similar helicase gene families in chromosome ends of Saccharomyces cerevisiae and Ustilago maydis, the initial association of helicase genes with fungal telomeres might date back to very early stages of the fungal evolution.

Pyret, a Ty3/Gypsy retrotransposon in Magnaporthe grisea contains an extra domain between the nucleocapsid and protease domains

A novel Ty3/Gypsy retrotransposon, named Pyret, was identified in the plant pathogenic fungus Magnaporthe grisea (anamorph Pyricularia oryzae). Pyret-related elements were distributed in a wide range of Pyricularia isolates from various gramineous plants. The Pyret element is 7250 bp in length with a 475 bp LTR and one conceptual ORF. The ORF contains seven nonsense mutations in the reading frame, indicating that the Pyret clone is lightly degenerate. Comparative domain analysis among retroelements revealed that Pyret exhibits an extra domain (WCCH domain) beyond the basic components of LTR retrotransposons. The WCCH domain consists of approximately 300 amino acids and is located downstream of the nucleocapsid domain. The WCCH domain is so named because it contains two repeats of a characteristic amino acid sequence, W-X(2)-C-X(4)-C-X(2)-H-X(3)-K. A WCCH motif-like sequence is found in the precoat protein of some geminiviruses, viral RNA-dependent RNA polymerase and also in an Arabidopsis protein of unknown function. Interestingly, detailed sequence analysis of the gag protein revealed that Pyret, as well as some other chromodomain-containing LTR retrotransposons, displays significant sequence homology with members of the gammaretroviruses (MLV-related retroviruses) in the capsid and nucleocapsid domains. This suggests that chromodomain-containing LTR retrotransposons and gammaretroviruses may share a common ancestor with the gag protein.

Pyricularia grisea Causing Gray Leaf Spot of Perennial Ryegrass Turf: Population Structure and Host Specificity

Gray leaf spot is a serious disease of perennial ryegrass (Lolium perenne) turf in the United States. Isolates of Pyricularia grisea causing the disease in perennial ryegrass were characterized using molecular markers and pathogenicity assays on various gramineous hosts. Genetic relationships among perennial ryegrass isolates were determined using different types of trans-posons as probes. Phylogenetic analysis using Pot2 and MGR586 probes, analyzed with AMOVA (analysis of molecular variance), showed that these isolates from perennial ryegrass consist of three closely related lineages. All the isolates belonged to a single mating type, MAT1-2. Among 20 isolates from 16 host species other than perennial ryegrass, only the isolates from wheat (Triticum aestivum) and triticale (× Triticosecale), showed notable similarity to the perennial ryegrass isolates based on their Pot2 fingerprints. The copy number and fingerprints of Pot2 and MGR586 in isolates of P. grisea from perennial ryegrass indicate that they are genetically distinct from the isolates derived from rice (Oryza sativa) in the United States. The perennial ryegrass isolates also had the same sequence in the internal transcribed spacer (ITS) region of the genes encoding ribosomal RNA as that of the wheat and triticale isolates, and exhibited rice isolate sequence polymorphisms. In pathogenicity assays, all the isolates of P. grisea from Legacy II perennial ryegrass caused characteristic blast symptoms on Marilee soft white winter wheat, Bennett hard red winter wheat, Era soft white spring wheat, and Presto triticale, and they were highly virulent on these hosts. An isolate from wheat and one from triticale (from Brazil) were also highly virulent on perennial ryegrass and Rebel III tall fescue (Festuca arundinacea). None of the isolates from perennial ryegrass caused the disease on Lagrue rice, and vice versa. Understanding the population structure of P. grisea isolates infecting perennial ryegrass and their relatedness to isolates from other gramineous hosts may aid in identifying alternate hosts for this pathogen.