Regions of microsynteny in Magnaporthe grisea and Neurospora crassa

Analysis of Magnaporthe oryzae Genome Reveals a Fungal Effector, Which Is Able to Induce Resistance Response in Transgenic Rice Line Containing Resistance Gene, Pi54

Rice blast caused by Magnaporthe oryzae is one of the most important diseases of rice. Pi54, a rice gene that imparts resistance to M. oryzae isolates prevalent in India, was already cloned but its avirulent counterpart in the pathogen was not known. After decoding the whole genome of an avirulent isolate of M. oryzae, we predicted 11440 protein coding genes and then identified four candidate effector proteins which are exclusively expressed in the infectious structure, appresoria. In silico protein modeling followed by interaction analysis between Pi54 protein model and selected four candidate effector proteins models revealed that Mo-01947_9 protein model encoded by a gene located at chromosome 4 of M. oryzae, interacted best at the Leucine Rich Repeat domain of Pi54 protein model. Yeast-two-hybrid analysis showed that Mo-01947_9 protein physically interacts with Pi54 protein. Nicotiana benthamiana leaf infiltration assay confirmed induction of hypersensitive response in the presence of Pi54 gene in a heterologous system. Genetic complementation test also proved that Mo-01947_9 protein induces avirulence response in the pathogen in presence of Pi54 gene. Here, we report identification and cloning of a new fungal effector gene which interacts with blast resistance gene Pi54 in rice.

Evolution of a subtilisin-like protease gene family in the grass endophytic fungus Epichloë festucae

BACKGROUND

Subtilisin-like proteases (SLPs) form a superfamily of enzymes that act to degrade protein substrates. In fungi, SLPs can play either a general nutritive role, or may play specific roles in cell metabolism, or as pathogenicity or virulence factors.

RESULTS

Fifteen different genes encoding SLPs were identified in the genome of the grass endophytic fungus Epichloë festucae. Phylogenetic analysis indicated that these SLPs belong to four different subtilisin families: proteinase K, kexin, pyrolysin and subtilisin. The pattern of intron loss and gain is consistent with this phylogeny. E. festucae is exceptional in that it contains two kexin-like genes. Phylogenetic analysis in Hypocreales fungi revealed an extensive history of gene loss and duplication.

CONCLUSION

This study provides new insights into the evolution of the SLP superfamily in filamentous fungi.

Novel bacterial artificial chromosome vector pUvBBAC for use in studies of the functional genomics of Listeria spp

Bacterial artificial chromosome (BAC) vectors are important tools for microbial genome research. We constructed a novel BAC vector, pUvBBAC, for replication in both gram-negative and gram-positive bacterial hosts. The pUvBBAC vector was used to generate a BAC library for the facultative intracellular pathogen Listeria monocytogenes EGD-e. The library had insert sizes ranging from 68 to 178 kb. We identified two recombinant BACs from the L. monocytogenes pUvBBAC library that each contained the entire virulence gene cluster (vgc) of L. monocytogenes and transferred them to a nonpathogenic Listeria innocua strain. Recombinant L. innocua strains harboring pUvBBAC+vgc1 and pUvBBAC+vgc2 produced the vgc-specific listeriolysin (LLO) and actin assembly protein ActA and represent the first reported cloning of the vgc locus in its entirety. The use of the novel broad-host-range BAC vector pUvBBAC extends the versatility of this technology and provides a powerful platform for detailed functional genomics of gram-positive bacteria as well as its use in explorative functional metagenomics.

Functional analysis of a fungal endophyte stress-activated MAP kinase

The ability of fungi to sense and respond rapidly to environmental stress is crucial for their survival in the wild. One of the most important pathways involved in this response is the stress-activated MAP (mitogen-activated protein) kinase pathway. We report here on the isolation of the stress-activated MAP kinase, sakA, from the fungal endophyte Epichloë festucae. Complementation of the stress sensitivity and cell cycle defects of an Schizosaccharomyces pombe sty1Delta mutant with sakA confirmed it encodes a functional MAP kinase. Analysis of an E. festucae DeltasakA mutant revealed sakA is essential for growth under conditions of temperature and osmotic stress in culture, and for sensitivity to the fungicide fludioxonil. However, the DeltasakA mutant shows no increased sensitivity to hydrogen peroxide. Given sakA can rescue the sty1Delta mutant from sensitivity to oxidative stress, SakA has the potential to sense and transduce oxidative stress signals. The DeltasakA mutant is also defective in conidia formation, suggesting a role for SakA in asexual development of E. festucae. The detection of elevated hydrogen peroxide production in the DeltasakA mutant suggests there may be a link between MAP kinase and ROS (reactive oxygen species) signalling pathways in E. festucae.

Gene loss rate: a probabilistic measure for the conservation of eukaryotic genes

The rate of conservation of a gene in evolution is believed to be correlated with its biological importance. Recent studies have devised various conservation measures for genes and have shown that they are correlated with several biological characteristics of functional importance. Specifically, the state-of-the-art propensity for gene loss (PGL) measure was shown to be strongly correlated with gene essentiality and its number of protein-protein interactions (PPIs). The observed correlation between conservation and functional importance varies however between conservation measures, underscoring the need for accurate and general measures for the rate of gene conservation. Here we develop a novel maximum-likelihood approach to computing the rate in which a gene is lost in evolution, motivated by the same principles as those underlying PGL. However, in difference to PGL which considers only the most parsimonious ancestral states of the internal nodes of the phylogenetic tree relating the species, our approach weighs in a probabilistic manner all possible ancestral states, and includes the branch length information as part of the probabilistic model. In application to data of 16 eukaryotic genomes, our approach shows higher correlations with experimental data than PGL, including data on gene lethality, level of connectivity in a PPI network and coherence within functionally related genes.

Comparative analysis of Phytophthora genes encoding secreted proteins reveals conserved synteny and lineage-specific gene duplications and deletions

Comparative analysis of two Phytophthora genomes revealed overall colinearity in four genomic regions consisting of a 1.5-Mb sequence of Phytophthora sojae and a 0.9-Mb sequence of P. ramorum. In these regions with conserved synteny, the gene order is largely similar; however, genome rearrangements also have occurred. Deletions and duplications often were found in association with genes encoding secreted proteins, including effectors that are important for interaction with host plants. Among secreted protein genes, different evolutionary patterns were found. Elicitin genes that code for a complex family of highly conserved Phytophthora-specific elicitors show conservation in gene number and order, and often are clustered. In contrast, the race-specific elicitor gene Avrlb-1 appeared to be missing from the region with conserved synteny, as were its five homologs that are scattered over the four genomic regions. Some gene families encoding secreted proteins were found to be expanded in one species compared with the other. This could be the result of either repeated gene duplications in one species or specific deletions in the other. These different evolutionary patterns may shed light on the functions of these secreted proteins in the biology and pathology of the two Phytophthora spp.

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.

Transcription factor families have much higher expansion rates in plants than in animals

Transcription factors (TFs), which are central to the regulation of gene expression, are usually members of multigene families. In plants, they are involved in diverse processes such as developmental control and elicitation of defense and stress responses. To investigate if differences exist in the expansion patterns of TF gene families between plants and other eukaryotes, we first used Arabidopsis (Arabidopsis thaliana) TFs to identify TF DNA-binding domains. These DNA-binding domains were then used to identify related sequences in 25 other eukaryotic genomes. Interestingly, among 19 families that are shared between animals and plants, more than 14 are larger in plants than in animals. After examining the lineage-specific expansion of TF families in two plants, eight animals, and two fungi, we found that TF families shared among these organisms have undergone much more dramatic expansion in plants than in other eukaryotes. Moreover, this elevated expansion rate of plant TF is not simply due to higher duplication rates of plant genomes but also to a higher degree of expansion compared to other plant genes. Further, in many Arabidopsis-rice (Oryza sativa) TF orthologous groups, the degree of lineage-specific expansion in Arabidopsis is correlated with that in rice. This pattern of parallel expansion is much more pronounced than the whole-genome trend in rice and Arabidopsis. The high rate of expansion among plant TF genes and their propensity for parallel expansion suggest frequent adaptive responses to selection pressure common among higher plants.

Insight into the genome of Aspergillus fumigatus: analysis of a 922 kb region encompassing the nitrate assimilation gene cluster

Aspergillus fumigatus is the most ubiquitous opportunistic filamentous fungal pathogen of human. As an initial step toward sequencing the entire genome of A. fumigatus, which is estimated to be approximately 30 Mb in size, we have sequenced a 922 kb region, contained within 16 overlapping bacterial artificial chromosome (BAC) clones. Fifty-four percent of the DNA is predicted to be coding with 341 putative protein coding genes. Functional classification of the proteins showed the presence of a higher proportion of enzymes and membrane transporters when compared to those of Saccharomyces cerevisiae. In addition to the nitrate assimilation gene cluster, the quinate utilisation gene cluster is also present on this 922 kb genomic sequence. We observed large scale synteny between A. fumigatus and Aspergillus nidulans by comparing this sequence to the A. nidulans genetic map of linkage group VIII.

Comparative sequence analysis of Sordaria macrospora and Neurospora crassa as a means to improve genome annotation

One of the most challenging parts of large scale sequencing projects is the identification of functional elements encoded in a genome. Recently, studies of genomes of up to six different Saccharomyces species have demonstrated that a comparative analysis of genome sequences from closely related species is a powerful approach to identify open reading frames and other functional regions within genomes [Science 301 (2003) 71, Nature 423 (2003) 241]. Here, we present a comparison of selected sequences from Sordaria macrospora to their corresponding Neurospora crassa orthologous regions. Our analysis indicates that due to the high degree of sequence similarity and conservation of overall genomic organization, S. macrospora sequence information can be used to simplify the annotation of the N. crassa genome.

Contig assembly and microsynteny analysis using a bacterial artificial chromosome library for Epichloë festucae, a mutualistic fungal endophyte of grasses

We constructed and characterized a bacterial artificial chromosome (BAC) library for Epichloë festucae, a genetically tractable fungal plant mutualist. The 6144 clone library with an average insert size of 87kb represents at least 18-fold coverage of the 29 Mb genome. We used the library to assemble a 110kb contig spanning the putative ornithine decarboxylase (odc) ortholog and subsequently expanded it to 228kb with a single walking step in each direction. Furthermore, we evaluated conservation of microsynteny between E. festucae and some model filamentous fungi by comparing sequence available from a 43kb region at the end of one BAC to publicly available fungal genome sequences. Orthologs to the 13 contiguous open reading frames (ORFs) identified in E. festucae are syntenic in Neurospora crassa and Magnaporthe grisea occurring in small sets of two, three or four colinear ORFs. This library is a valuable resource for research into traits important for the development and maintenance of a plant-fungus mutualistic symbiosis.

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.

An ordered collection of expressed sequences from Cryphonectria parasitica and evidence of genomic microsynteny with Neurospora crassa and Magnaporthe grisea

Cryphonectria parasitica, the causative agent of chestnut blight, has proven to be a tractable experimental system for studying fungal pathogenesis. Moreover, the development of infectious cDNA clones of C. parasitica hypoviruses, capable of attenuating fungal virulence, has provided the opportunity to examine molecular aspects of fungal plant pathogenesis in the context of biological control. In order to establish a genomic base for future studies of C. parasitica, the authors have analysed a collection of expressed sequences. A mixed cDNA library was prepared from RNA isolated from wild-type (virus-free) and hypovirus-infected C. parasitica strains. Plasmid DNA was recovered from individual transformants and sequenced from the 5' end of the insert. Contig analysis of the collected sequences revealed that they represented approximately 2200 individual ORFs. An assessment of functional diversity present in this collection was achieved by using the BLAST software utilities and the NCBI protein database. Candidate genes were identified with significant potential relevance to C. parasitica growth, development, pathogenesis and vegetative incompatibility. Additional investigations of a 12.9 kbp genomic region revealed microsynteny between C. parasitica and both Neurospora crassa and Magnaporthe grisea, two closely related fungi. These data represent the largest collection of sequence information currently available for C. parasitica and are now forming the basis of further studies using microarray analyses to determine global changes in transcription that occur in response to hypovirus infection.

Characterization of the genomic organization of the region bordering the centromere of chromosome V of Podospora anserina by direct sequencing

A Podospora anserina BAC library of 4800 clones has been constructed in the vector pBHYG allowing direct selection in fungi. Screening of the BAC collection for centromeric sequences of chromosome V allowed the recovery of clones localized on either sides of the centromere, but no BAC clone was found to contain the centromere. Seven BAC clones containing 322,195 and 156,244bp from either sides of the centromeric region were sequenced and annotated. One 5S rRNA gene, 5 tRNA genes, and 163 putative coding sequences (CDS) were identified. Among these, only six CDS seem specific to P. anserina. The gene density in the centromeric region is approximately one gene every 2.8kb. Extrapolation of this gene density to the whole genome of P. anserina suggests that the genome contains about 11,000 genes. Synteny analyses between P. anserina and Neurospora crassa show that co-linearity extends at the most to a few genes, suggesting rapid genome rearrangements between these two species.

The rice blast pathosystem as a case study for the development of new tools and raw materials for genome analysis of fungal plant pathogens

Fungi have an astounding and diverse impact on this planet. While they are agents of human diseases and the cause of allergic reactions, factories for the conversion of carbon in environmental and industrially adapted systems, and potential biological weapons, their importance as plant pathogens is unparalleled. In plants alone, fungi cause tens of thousands of different diseases and are responsible for massive losses of food, fiber and forestry at an estimated annual cost of hundreds of billions of dollars. These losses are not only realized in the incomes of individual farmers and state economies, but contribute significantly to world hunger problems and issues relating to safeguarding a global food supply. Our collective understanding of how fungi, particularly plant pathogens, grow, reproduce, identify a host and cause disease is still at a formative stage. There is an equal lack of detailed knowledge about how a plant recognizes that it is being attacked and then mounts an adequate defense response. The advent of genomic technologies has given researchers an unprecedented opportunity to address these mysteries in a powerful and more holistic manner. Where the genetic revolution of only a few years ago allowed for the characterization of single genes, today's genomic technologies are facilitating the evaluation of the entire complement of genes in an organism and the discovery of the suites of genes that act during any one time or particular condition. This review will describe the recent development of tools for whole or partial genome analysis and multigenome comparisons. Th discussion focuses on the rice blast pathosystem as a case study.

A genomics approach to crop pest and disease research

Genome-wide analyses of gene function and gene expression are beginning to yield valuable information in many areas of biological research, and these genomic tools are now being applied to crop pest and disease research. DNA sequencing of cDNA libraries to generate sets of expressed sequence tags (ESTs) are allowing gene compendiums for crop diseases to be compiled. Annotation of such data collections is also providing a wealth of functional information about gene products through similarities to proteins with known function. The next phase of the functional genomics era will be to employ large-scale techniques to knock out or silence genes in order to synthesize gene-specific mutants for phenotypic analysis and to use micro-array methodology to analyze global gene expression, protein turnover and protein processing during the processes of parasitism and colonization. Application of these technologies promises to accelerate the pace that biological information relevant to crop protection accrues. The ability of researchers to assimilate this information into complex models and workable hypotheses is, thus, set to revolutionize the way we study pests and diseases of crop plants.

Trichothecene nonproducer Gibberella species have both functional and nonfunctional 3-O-acetyltransferase genes

The trichothecene 3-O-acetyltransferase gene (FgTri101) required for trichothecene production by Fusarium graminearum is located between the phosphate permease gene (pho5) and the UTP-ammonia ligase gene (ura7). We have cloned and sequenced the pho5-to-ura7 regions from three trichothecene nonproducing Fusarium (i.e., F. oxysporum, F. moniliforme, and Fusarium species IFO 7772) that belong to the teleomorph genus Gibberella. BLASTX analysis of these sequences revealed portions of predicted polypeptides with high similarities to the TRI101 polypeptide. While FspTri101 (Fusarium species Tri101) coded for a functional 3-O-acetyltransferase, FoTri101 (F. oxysporum Tri101) and FmTri101 (F. moniliforme Tri101) were pseudogenes. Nevertheless, F. oxysporum and F. moniliforme were able to acetylate C-3 of trichothecenes, indicating that these nonproducers possess another as yet unidentified 3-O-acetyltransferase gene. By means of cDNA expression cloning using fission yeast, we isolated the responsible FoTri201 gene from F. oxysporum; on the basis of this sequence, FmTri201 has been cloned from F. moniliforme by PCR techniques. Both Tri201 showed only a limited level of nucleotide sequence similarity to FgTri101 and FspTri101. The existence of Tri101 in a trichothecene nonproducer suggests that this gene existed in the fungal genome before the divergence of producers from nonproducers in the evolution of Fusarium species.

The Hypocrea jecorina gal10 (uridine 5'-diphosphate-glucose 4-epimerase-encoding) gene differs from yeast homologues in structure, genomic organization and expression

As part of a comprehensive study on lactose metabolism in Hypocrea jecorina (anamorph: Trichoderma reesei), a genomic clone of the gal10 gene encoding H. jecorina uridine 5'-diphosphate (UDP)-glucose 4-epimerase has been cloned and sequenced. It contains an open reading frame of 1548-base pair, interrupted by three introns, and encoding a 370-amino acids protein with similarity to pro- and eukaryotic UDP-glucose-4-epimerases. H. jecorina Gal10 does not contain the C-terminal mutarotase domain which is present in yeast Gal10 proteins but is able to functionally complement a corresponding Saccharomyces cerevisiae gal10 mutant. gal10 is not clustered with other H. jecorina gal genes (gal7, gene encoding galactose-1-phosphate uridylyltransferase and gal1, gene encoding galactokinase). The genomic location of H. jecorina gal10 and gal7 was syntenic with that in Neurospora crassa and colinear over an area of 6 and 3.5-kilobase. gal10 is constitutively expressed, and--unlike H. jecorina gal7--not further stimulated by D-galactose or L-arabinose or its corresponding polyols.

Time for a blast: genomics of Magnaporthe grisea

Summary The rice blast fungus Magnaporthe grisea causes one of the most destructive diseases of rice. Genetic studies of this important pathogen during the past decade have made it an excellent system for investigating fungal-plant interactions. Recently, the well coordinated efforts by the rice blast community have also made significant progresses in genomics studies of M. grisea. BAC contigs were assembled and integrated with the high density genetic map, and 12 674 BAC-ends were sequenced as genome survey sequences. Over ten BAC clones have been sequenced by the shot-gun approach. Preliminary analysis with one completely sequenced BAC clone indicated that M. grisea may have a gene density of 4.2 kb/gene and contain approximately 9000 genes. In addition, > 10 000 ESTs have been sequenced from several cDNA libraries representing different stages of fungal growth, differentiation and plant infection. Currently, there are ongoing projects to sequence additional ESTs, telomeres, and 6X coverage of the M. grisea genome. A large-scale functional genomics project on M. grisea and its interaction with rice is also underway. It is an exciting time for genomics studies in M. grisea, and the resources generated in these studies will certainly enhance our understanding of fungal pathogenicity.

Genomics of phytopathogenic fungi and the development of bioinformatic resources

Genomic resources available to researchers studying phytopathogenic fungi are limited. Here, we briefly review the genomic and bioinformatic resources available and the current status of fungal genomics. We also describe a relational database containing sequences of expressed sequence tags (ESTs) from three phytopathogenic fungi, Blumeria graminis, Magnaporthe grisea, and Mycosphaerella graminicola, and the methods and underlying principles required for its construction. The database contains significant annotation for each EST sequence and is accessible at http://cogeme.ex.ac.uk. An easy-to-use interface allows the user to identify gene sequences by using simple text queries or homology searches. New querying functions and large sequence sets from a variety of phytopathogenic species will be incorporated in due course.

Molecular characterization of a cystathionine beta-synthase gene, CBS1, in Magnaporthe grisea

CBS1 from Magnaporthe grisea is a structural and functional homolog of the cystathionine beta-synthase (CBS) gene from Saccharomyces cerevisiae. Our studies indicated that M. grisea can utilize homocysteine and methionine through a CBS-independent pathway. The results also revealed responses of M. grisea to homocysteine that are reminiscent of human homocystinuria.

A genetic map of Blumeria graminis based on functional genes, avirulence genes, and molecular markers

A genetic map of the powdery mildew fungus, Blumeria graminis f. sp. hordei, an obligate biotrophic pathogen of barley, is presented. The linkage analysis was conducted on 81 segregating haploid progeny isolates from a cross between 2 isolates differing in seven avirulence genes. A total of 359 loci were mapped, comprising 182 amplified fragment length polymorphism markers, 168 restriction fragment length polymorphism markers including 42 LTR-retrotransposon loci and 99 expressed sequence tags (ESTs), all the seven avirulence genes, and a marker closely linked to the mating type gene. The markers are distributed over 34 linkage groups covering a total of 2114 cM. Five avirulence genes were found to be linked and mapped in clusters of three and two, and two were unlinked. The Avr(a6) gene was found to be closely linked to markers suitable for a map-based cloning approach. A linkage between ESTs allowed us to demonstrate examples of synteny between genes in B. graminis and Neurospora crassa.

Conservation of long-range synteny and microsynteny between the genomes of two distantly related nematodes

BACKGROUND

Comparisons between the genomes of the closely related nematodes Caenorhabditis elegans and Caenorhabditis briggsae reveal high rates of rearrangement, with a bias towards within-chromosome events. To assess whether this pattern is true of nematodes in general, we have used genome sequence to compare two nematode species that last shared a common ancestor approximately 300 million years ago: the model C. elegans and the filarial parasite Brugia malayi.

RESULTS

An 83 kb region flanking the gene for Bm-mif-1 (macrophage migration inhibitory factor, a B. malayi homolog of a human cytokine) was sequenced. When compared to the complete genome of C. elegans, evidence for conservation of long-range synteny and microsynteny was found. Potential C. elegans orthologs for II of the 12 protein-coding genes predicted in the B. malayi sequence were identified. Ten of these orthologs were located on chromosome I, with eight clustered in a 2.3 Mb region. While several, relatively local, intrachromosomal rearrangements have occurred, the order, composition, and configuration of two gene clusters, each containing three genes, was conserved. Comparison of B. malayi BAC-end genome survey sequence to C. elegans also revealed a bias towards intrachromosome rearrangements.

CONCLUSIONS

We suggest that intrachromosomal rearrangement is a major force driving chromosomal organization in nematodes, but is constrained by the interdigitation of functional elements of neighboring genes.