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

Transposable Elements as a Source of Novel Repetitive DNA in the Eukaryote Genome

The impact of transposable elements (TEs) on the evolution of the eukaryote genome has been observed in a number of biological processes, such as the recruitment of the host's gene expression network or the rearrangement of genome structure. However, TEs may also provide a substrate for the emergence of novel repetitive elements, which contribute to the generation of new genomic components during the course of the evolutionary process. In this review, we examine published descriptions of TEs that give rise to tandem sequences in an attempt to comprehend the relationship between TEs and the emergence of de novo satellite DNA families in eukaryotic organisms. We evaluated the intragenomic behavior of the TEs, the role of their molecular structure, and the chromosomal distribution of the paralogous copies that generate arrays of repeats as a substrate for the emergence of new repetitive elements in the genome. We highlight the involvement and importance of TEs in the eukaryote genome and its remodeling processes.

Dynamics of Verticillium dahliae race 1 population under managed agricultural ecosystems

BACKGROUND

Plant pathogens and their hosts undergo adaptive changes in managed agricultural ecosystems, by overcoming host resistance, but the underlying genetic adaptations are difficult to determine in natural settings. Verticillium dahliae is a fungal pathogen that causes Verticillium wilt on many economically important crops including lettuce. We assessed the dynamics of changes in the V. dahliae genome under selection in a long-term field experiment.

RESULTS

In this study, a field was fumigated before the Verticillium dahliae race 1 strain (VdLs.16) was introduced. A derivative 145-strain population was collected over a 6-year period from this field in which a seggregating population of lettuce derived from Vr1/vr1 parents were evaluated. We de novo sequenced the parental genome of VdLs.16 strain and resequenced the derivative strains to analyze the genetic variations that accumulate over time in the field cropped with lettuce. Population genomics analyses identified 2769 single-nucleotide polymorphisms (SNPs) and 750 insertion/deletions (In-Dels) in the 145 isolates compared with the parental genome. Sequence divergence was identified in the coding sequence regions of 378 genes and in the putative promoter regions of 604 genes. Five-hundred and nine SNPs/In-Dels were identified as fixed. The SNPs and In-Dels were significantly enriched in the transposon-rich, gene-sparse regions, and in those genes with functional roles in signaling and transcriptional regulation.

CONCLUSIONS

Under the managed ecosystem continuously cropped to lettuce, the local adaptation of V. dahliae evolves at a whole genome scale to accumulate SNPs/In-Dels nonrandomly in hypervariable regions that encode components of signal transduction and transcriptional regulation.

Genome wide analysis of the transition to pathogenic lifestyles in Magnaporthales fungi

The rice blast fungus Pyricularia oryzae (syn. Magnaporthe oryzae, Magnaporthe grisea), a member of the order Magnaporthales in the class Sordariomycetes, is an important plant pathogen and a model species for studying pathogen infection and plant-fungal interaction. In this study, we generated genome sequence data from five additional Magnaporthales fungi including non-pathogenic species, and performed comparative genome analysis of a total of 13 fungal species in the class Sordariomycetes to understand the evolutionary history of the Magnaporthales and of fungal pathogenesis. Our results suggest that the Magnaporthales diverged ca. 31 millon years ago from other Sordariomycetes, with the phytopathogenic blast clade diverging ca. 21 million years ago. Little evidence of inter-phylum horizontal gene transfer (HGT) was detected in Magnaporthales. In contrast, many genes underwent positive selection in this order and the majority of these sequences are clade-specific. The blast clade genomes contain more secretome and avirulence effector genes, which likely play key roles in the interaction between Pyricularia species and their plant hosts. Finally, analysis of transposable elements (TE) showed differing proportions of TE classes among Magnaporthales genomes, suggesting that species-specific patterns may hold clues to the history of host/environmental adaptation in these fungi.

Elevated temperature increases in planta expression levels of virulence related genes in Magnaporthe oryzae and compromises resistance in Oryza sativa cv. Nipponbare

Temperature changes have the potential to alter the incidence and severity of plant disease epidemics and pressures, as well as to reshape the co-evolutionary relationships between plants and pathogens. However, the molecular basis of temperature modulation of pathogenicity of plant pathogens is still unclear. Here, we studied the effect of temperature on biomass of Magnaporthe oryzae in planta using qPCR. Additionally, the transcriptomes of M. oryzae and rice were analysed using RNA-seq. Rice seedlings were exposed to 35°C and 28°C for 7 days before pathogen inoculation. Inoculated plants were kept in the dark at 28°C for 24h and later re-exposed to 35°C and 28°C for an additional 24h before sample collection. Plants grown and predisposed to 35°C prior to inoculation exhibited accelerated tissue necrosis compared with plants grown and inoculated at 28°C. In accordance with the disease severity observed on infected leaves, in planta fungal biomass was significantly higher at 35°C than 28°C. Moreover, M. oryzae exhibited increased expression levels of putative fungal effector genes in plants exposed to 35°C compared with plants exposed to 28°C. Collectively, this study revealed that temperature elevation could favour M. oryzae infection by compromising plant resistance and accelerating plant tissue colonisation with the pathogen.

A BAC based physical map and genome survey of the rice false smut fungus Villosiclava virens

Background

Rice false smut caused by Villosiclava virens is a devastating fungal disease that spreads in major rice-growing regions throughout the world. However, the genomic information for this fungal pathogen is limited and the pathogenic mechanism of this disease is still not clear. To facilitate genetic, molecular and genomic studies of this fungal pathogen, we constructed the first BAC-based physical map and performed the first genome survey for this species.

Results

High molecular weight genomic DNA was isolated from young mycelia of the Villosiclava virens strain UV-8b and a high-quality, large-insert and deep-coverage Bacterial Artificial Chromosome (BAC) library was constructed with the restriction enzyme HindIII. The BAC library consisted of 5,760 clones, which covers 22.7-fold of the UV-8b genome, with an average insert size of 140 kb and an empty clone rate of lower than 1%. BAC fingerprinting generated successful fingerprints for 2,290 BAC clones. Using the fingerprints, a whole genome-wide BAC physical map was constructed that contained 194 contigs (2,035 clones) spanning 51.2 Mb in physical length. Bidirectional-end sequencing of 4,512 BAC clones generated 6,560 high quality BAC end sequences (BESs), with a total length of 3,030,658 bp, representing 8.54% of the genome sequence. Analysis of the BESs revealed general genome information, including 51.52% GC content, 22.51% repetitive sequences, 376.12/Mb simple sequence repeat (SSR) density and approximately 36.01% coding regions. Sequence comparisons to other available fungal genome sequences through BESs showed high similarities to Metarhizium anisopliae, Trichoderma reesei, Nectria haematococca and Cordyceps militaris, which were generally in agreement with the 18S rRNA gene analysis results.

Conclusion

This study provides the first BAC-based physical map and genome information for the important rice fungal pathogen Villosiclava virens. The BAC clones, physical map and genome information will serve as fundamental resources to accelerate the genetic, molecular and genomic studies of this pathogen, including positional cloning, comparative genomic analysis and whole genome sequencing. The BAC library and physical map have been opened to researchers as public genomic resources (http://gresource.hzau.edu.cn/resource/resource.html).

Genomic organization and sequence dynamics of the AvrPiz-t locus in Magnaporthe oryzae

Plants utilize multiple layers of defense mechanisms to fight against the invasion of diverse pathogens. The R gene mediates resistance, in most cases, dependent on the co-existence of its cognate pathogen-derived avirulence (Avr) gene. The rice blast R gene Piz-t corresponds in gene-for-gene fashion to the Magnaporthe oryzae Avr gene AvrPiz-t. In this study, we determined and compared the genomic sequences surrounding the AvrPiz-t gene in both avirulent and virulent isolates, designating as AvrPiz-t-ZB15 and avrPiz-t-70-15 regions, respectively. The sequence of the AvrPiz-t-ZB15 region is 120966 bp whereas avrPiz-t-70-15 is 146292 bp in length. The extreme sequence similarity and good synteny in gene order and content along with the absence of two predicted genes in the avrPiz-t-70-15 region were observed in the predicted protein-coding regions in the AvrPiz-t locus. Nevertheless, frequent presence/absence and highly dynamic organization of transposable elements (TEs) were identified, representing the major variation of the AvrPiz-t locus between different isolates. Moreover, TEs constitute 27.3% and 43.2% of the genomic contents of the AvrPiz-t-ZB15 and avrPiz-t-70-15 regions, respectively, indicating that TEs contribute largely to the organization and evolution of AvrPiz-t locus. The findings of this study suggest that M. oryzae could benefit in an evolutionary sense from the presence of active TEs in genes conferring avirulence and provide an ability to rapidly change and thus to overcome host R genes.

A physical map of Brassica oleracea shows complexity of chromosomal changes following recursive paleopolyploidizations

Background

Evolution of the Brassica species has been recursively affected by polyploidy events, and comparison to their relative, Arabidopsis thaliana, provides means to explore their genomic complexity.

Results

A genome-wide physical map of a rapid-cycling strain of B. oleracea was constructed by integrating high-information-content fingerprinting (HICF) of Bacterial Artificial Chromosome (BAC) clones with hybridization to sequence-tagged probes. Using 2907 contigs of two or more BACs, we performed several lines of comparative genomic analysis. Interspecific DNA synteny is much better preserved in euchromatin than heterochromatin, showing the qualitative difference in evolution of these respective genomic domains. About 67% of contigs can be aligned to the Arabidopsis genome, with 96.5% corresponding to euchromatic regions, and 3.5% (shown to contain repetitive sequences) to pericentromeric regions. Overgo probe hybridization data showed that contigs aligned to Arabidopsis euchromatin contain ~80% of low-copy-number genes, while genes with high copy number are much more frequently associated with pericentromeric regions. We identified 39 interchromosomal breakpoints during the diversification of B. oleracea and Arabidopsis thaliana, a relatively high level of genomic change since their divergence. Comparison of the B. oleracea physical map with Arabidopsis and other available eudicot genomes showed appreciable 'shadowing' produced by more ancient polyploidies, resulting in a web of relatedness among contigs which increased genomic complexity.

Conclusions

A high-resolution genetically-anchored physical map sheds light on Brassica genome organization and advances positional cloning of specific genes, and may help to validate genome sequence assembly and alignment to chromosomes.

All the physical mapping data is freely shared at a WebFPC site (http://lulu.pgml.uga.edu/fpc/WebAGCoL/brassica/WebFPC/; Temporarily password-protected: account: pgml; password: 123qwe123.

Repeat induced point mutation in two asexual fungi, Aspergillus niger and Penicillium chrysogenum

Repeat induced point mutation (RIP) is a gene silencing mechanism present in fungal genomes. During RIP, duplicated sequences are efficiently and irreversibly mutated by transitions from C:G to T:A. For the first time, we have identified traces of RIP in transposable elements of Aspergillus niger and Penicillium chrysogenum, two biotechnologically relevant fungi. We found that RIP in P. chrysogenum has affected a large set of sequences, which also contain other mutations. On the other hand, RIP in A. niger is limited to only few sequences, but literally all mutations are RIP-like. Surprisingly, RIP occurred only in transposon sequences that have disrupted open reading frames in A. niger, a phenomenon not yet reported for other fungi. In both fungal species, we identified two sequences with strong sequence similarity to Neurospora crassa RID. RID is a putative DNA methyltransferase and the only known enzyme involved in the RIP process. Our findings suggest that both A. niger and P. chrysogenum either had a sexual past or have a sexual potential. These findings have important implications for future strain development of these fungi.

Chromodomains direct integration of retrotransposons to heterochromatin

The enrichment of mobile genetic elements in heterochromatin may be due, in part, to targeted integration. The chromoviruses are Ty3/gypsy retrotransposons with chromodomains at their integrase C termini. Chromodomains are logical determinants for targeting to heterochromatin, because the chromodomain of heterochromatin protein 1 (HP1) typically recognizes histone H3 K9 methylation, an epigenetic mark characteristic of heterochromatin. We describe three groups of chromoviruses based on amino acid sequence relationships of their integrase C termini. Genome sequence analysis indicates that representative chromoviruses from each group are enriched in gene-poor regions of the genome relative to other retrotransposons, and when fused to fluorescent marker proteins, the chromodomains target proteins to specific subnuclear foci coincident with heterochromatin. The chromodomain of the fungal element, MAGGY, interacts with histone H3 dimethyl- and trimethyl-K9, and when the MAGGY chromodomain is fused to integrase of the Schizosaccharomyces pombe Tf1 retrotransposon, new Tf1 insertions are directed to sites of H3 K9 methylation. Repetitive sequences such as transposable elements trigger the RNAi pathway resulting in their epigenetic modification. Our results suggest a dynamic interplay between retrotransposons and heterochromatin, wherein mobile elements recognize heterochromatin at the time of integration and then perpetuate the heterochromatic mark by triggering epigenetic modification.

Bacterial artificial chromosome-based physical map of Gibberella zeae (Fusarium graminearum)

Fusarium graminearum is the primary causal pathogen of Fusarium head blight of wheat and barley. To accelerate genomic analysis of F. graminearum, we developed a bacterial artificial chromosome (BAC)-based physical map and integrated it with the genome sequence and genetic map. One BAC library, developed in the HindIII restriction enzyme site, consists of 4608 clones with an insert size of approximately 107 kb and covers about 13.5 genome equivalents. The other library, developed in the BamHI restriction enzyme site, consists of 3072 clones with an insert size of approximately 95 kb and covers about 8.0 genome equivalents. We fingerprinted 2688 clones from the HindIII library and 1536 clones from the BamHI library and developed a physical map of F. graminearum consisting of 26 contigs covering 39.2 Mb. Comparison of our map with the F. graminearum genome sequence showed that the size of our physical map is equivalent to the 36.1 Mb of the genome sequence. We used 31 sequence-based genetic markers, randomly spaced throughout the genome, to integrate the physical map with the genetic map. We also end-sequenced 17 BamHI BAC clones and identified 4 clones that spanned gaps in the genome sequence. Our new integrated map is highly reliable and useful for a variety of genomics studies.

Genome-wide analysis of T-DNA integration into the chromosomes of Magnaporthe oryzae

Agrobacterium tumefaciens-mediated transformation (ATMT) has become a prevalent tool for functional genomics of fungi, but our understanding of T-DNA integration into the fungal genome remains limited relative to that in plants. Using a model plant-pathogenic fungus, Magnaporthe oryzae, here we report the most comprehensive analysis of T-DNA integration events in fungi and the development of an informatics infrastructure, termed a T-DNA analysis platform (TAP). We identified a total of 1110 T-DNA-tagged locations (TTLs) and processed the resulting data via TAP. Analysis of the TTLs showed that T-DNA integration was biased among chromosomes and preferred the promoter region of genes. In addition, irregular patterns of T-DNA integration, such as chromosomal rearrangement and readthrough of plasmid vectors, were also observed, showing that T-DNA integration patterns into the fungal genome are as diverse as those of their plant counterparts. However, overall the observed junction structures between T-DNA borders and flanking genomic DNA sequences revealed that T-DNA integration into the fungal genome was more canonical than those observed in plants. Our results support the potential of ATMT as a tool for functional genomics of fungi and show that the TAP is an effective informatics platform for handling data from large-scale insertional mutagenesis.

The evolutionary history of cytochrome P450 genes in four filamentous Ascomycetes

Background

The Cytochrome P450 system is important in fungal evolution for adapting to novel ecological niches. To elucidate the evolutionary process of cytochrome P450 genes in fungi with different life styles, we studied the patterns of gene gains and losses in the genomes of four filamentous Ascomycetes, including two saprotrophs (Aspergillus nidulans (AN) and Neurospora crassa (NC)) and two plant pathogens (Fusarium graminearum (FG) and Magnaporthe grisea (MG)).

Results

A total of 376 P450 genes were assigned to 168 families according to standard nomenclature. On average, only 1 to 2 genes per family were in each genome. To resolve conflicting results between different clustering analyses and standard family designation, a higher order relationship was formulated. 376 genes were clustered into 115 clans. Subsequently a novel approach based on parsimony was developed to build the evolutionary models. Based on these analyses, a core of 30 distinct clans of P450s was defined. The core clans experienced contraction in all four fungal lineages while new clans expanded in all with exception of NC. MG experienced more genes and clans gains compared to the other fungi. Parsimonious analyses unanimously supported one species topology for the four fungi.

Conclusion

The four studied fungi exhibit unprecedented diversity in their P450omes in terms of coding sequence, intron-exon structures and genome locations, suggesting a complicated evolutionary history of P450s in filamentous Ascomycetes. Clan classification and a novel strategy were developed to study evolutionary history. Contraction of core clans and expansion of novel clans were identified. The exception was the NC lineage, which exhibited pure P450 gene loss.

Magnaporthe as a model for understanding host-pathogen interactions

The rice blast pathosystem has been the subject of intense interest in part because of the importance of the disease to world agriculture, but also because both Magnaporthe oryzae and its host are amenable to advanced experimental approaches. The goal of this review is to provide an overview of the system and to point out recent significant studies that update our understanding of the biology of M. oryzae. The genome sequence of M. oryzae has provided insight into how genome structure and pathogen population genetic variability has been shaped by transposable elements. The sequence allows systematic approaches to long-standing areas of investigation, including pathogen development and the molecular basis of compatible and incompatible interactions with its host. Rice blast provides an integrated system to illustrate most of the important concepts governing fungal/plant interactions and serves as an excellent starting point for gaining a broad perspective of issues in plant pathology.

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.

The role of transposable element clusters in genome evolution and loss of synteny in the rice blast fungus Magnaporthe oryzae

Analysis of the Magnaporthe oryzae chromosome 7 and comparison with syntenic regions in other fungal genomes suggests that transposable elements create localized segments with increased rates of chromosomal rearrangements, gene duplications and gene evolution.

Background

Transposable elements are abundant in the genomes of many filamentous fungi, and have been implicated as major contributors to genome rearrangements and as sources of genetic variation. Analyses of fungal genomes have also revealed that transposable elements are largely confined to distinct clusters within the genome. Their impact on fungal genome evolution is not well understood. Using the recently available genome sequence of the plant pathogenic fungus Magnaporthe oryzae, combined with additional bacterial artificial chromosome clone sequences, we performed a detailed analysis of the distribution of transposable elements, syntenic blocks, and other features of chromosome 7.

Results

We found significant levels of conserved synteny between chromosome 7 and the genomes of other filamentous fungi, despite more than 200 million years of divergent evolution. Transposable elements are largely restricted to three clusters located in chromosomal segments that lack conserved synteny. In contradiction to popular evolutionary models and observations from other model organism genomes, we found a positive correlation between recombination rate and the distribution of transposable element clusters on chromosome 7. In addition, the transposable element clusters are marked by more frequent gene duplications, and genes within the clusters have greater sequence diversity to orthologous genes from other fungi.

Conclusion

Together, these data suggest that transposable elements have a profound impact on the M. oryzae genome by creating localized segments with increased rates of chromosomal rearrangements, gene duplications and gene evolution.

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.

Comparative physical mapping links conservation of microsynteny to chromosome structure and recombination in grasses

Nearly finished sequences for model organisms provide a foundation from which to explore genomic diversity among other taxonomic groups. We explore genome-wide microsynteny patterns between the rice sequence and two sorghum physical maps that integrate genetic markers, bacterial artificial chromosome (BAC) fingerprints, and BAC hybridization data. The sorghum maps largely tile a genomic component containing 41% of BACs but 80% of single-copy genes that shows conserved microsynteny with rice and partially tile a nonsyntenic component containing 46% of BACs but only 13% of single-copy genes. The remaining BACs are centromeric (4%) or unassigned (8%). The two genomic components correspond to cytologically discernible "euchromatin" and "heterochromatin." Gene and repetitive DNA distributions support this classification. Greater microcolinearity in recombinogenic (euchromatic) than nonrecombinogenic (heterochromatic) regions is consistent with the hypothesis that genomic rearrangements are usually deleterious, thus more likely to persist in nonrecombinogenic regions by virtue of Muller's ratchet. Interchromosomal centromeric rearrangements may have fostered diploidization of a polyploid cereal progenitor. Model plant sequences better guide studies of related genomes in recombinogenic than nonrecombinogenic regions. Bridging of 35 physical gaps in the rice sequence by sorghum BAC contigs illustrates reciprocal benefits of comparative approaches that extend at least across the cereals and perhaps beyond.

The genome sequence of the rice blast fungus Magnaporthe grisea

Magnaporthe grisea is the most destructive pathogen of rice worldwide and the principal model organism for elucidating the molecular basis of fungal disease of plants. Here, we report the draft sequence of the M. grisea genome. Analysis of the gene set provides an insight into the adaptations required by a fungus to cause disease. The genome encodes a large and diverse set of secreted proteins, including those defined by unusual carbohydrate-binding domains. This fungus also possesses an expanded family of G-protein-coupled receptors, several new virulence-associated genes and large suites of enzymes involved in secondary metabolism. Consistent with a role in fungal pathogenesis, the expression of several of these genes is upregulated during the early stages of infection-related development. The M. grisea genome has been subject to invasion and proliferation of active transposable elements, reflecting the clonal nature of this fungus imposed by widespread rice cultivation.

The C-terminal chromodomain-like module in the integrase domain is crucial for high transposition efficiency of the retrotransposon MAGGY

MAGGY is a Ty3/Gypsy retrotransposon, which was identified in the rice blast fungus Magnaporthe oryzae. Some Ty3/Gypsy retrotransposons, including MAGGY, contain a chromodomain-like module (CLM) in the C-terminus of the integrase domain. We have made a series of MAGGY mutants to examine the role of the CLM in the transposition activity of the element. Introduction of a mutation at different positions in the MAGGY integrase revealed that a loss or alteration of the CLM resulted in a drastic decrease in the transposition activity of the element. Our results indicate that the CLM may confer high transposition activity to the element.