Genomic insights into the rice blast fungus through estimation of gene emergence time in phylogenetic context

The rice blast fungus, Magnaporthe oryzae, is an important pathogen of rice plants. It is well known that genes encoded in the genome have different evolutionary histories that are related to their functions. Phylostratigraphy is a method that correlates the evolutionary origin of genes with evolutionary transitions. Here we applied phylostratigraphy to partition total gene content of M. oryzae into distinct classes (phylostrata), which we designated PS1 to PS7, based on estimation of their emergence time. Genes in individual phylostrata did not show significant biases in their global distribution among seven chromosomes, but at the local level, clustering of genes belonging to the same phylostratum was observed. Our phylostrata-wide analysis of genes revealed that genes in the same phylostratum tend to be similar in many physical and functional characteristics such as gene length and structure, GC contents, codon adaptation index, and level of transcription, which correlates with biological functions in evolutionary context. We also found that a significant proportion of genes in the genome are orphans, for which no orthologs can be detected in the database. Among them, we narrowed down to seven orphan genes having transcriptional and translational evidences, and showed that one of them is implicated in asexual reproduction and virulence, suggesting ongoing evolution in this fungus through lineage-specific genes. Our results provide genomic basis for linking functions of pathogenicity factors and gene emergence time.