BACKGROUND
colletotrichum graminicola is a filamentous ascomycete that causes anthracnose leaf blight and stalk rot of maize leading to estimated annual yield losses in the us of about 6%.
in addition to its economical importance, c. graminicola is among the best characterized and most tractable fungi of the genus colletotrichum, a genus comprising a broad range of hemibiotrophic plant pathogens that represent a constant threat to fruit and vegetable production worldwide .
interesting aspects of colletotrichum pathology are the pressure-driven penetration strategy and the hemibiotrophic lifestyle that is featured by a fundamental switch in nutrition.
an initial short biotrophic phase, during which the host cell remains alive, is followed by highly destructive necrotrophic development characterized by extended areas of killed host tissue
<cit> .
following germination of conidia, specialized infection cells called appressoria are differentiated on the host epidermis.
accumulation of osmolytes in appressoria of c. graminicola generates a hydrostatic pressure corresponding to approx.
<dig>  mpa, which is among the highest turgor pressures known in living cells
<cit> .
this enormous pressure is utilized by the fungus to forcefully penetrate plant epidermis cells
<cit> .
upon successful penetration of the host epidermis, biotrophic hyphae of c. graminicola initially colonize living plant cells as evidenced by the integrity of the host plasma membrane
<cit> .
this stage of symptomless parasitic development is terminated at the onset of necrotrophy.
fast-growing hyphae breach the walls of adjacent cells, ramify within the leave and kill host cells, leading to massive destruction of maize tissue and visible necrosis.
finally, acervuli are formed on necrotic tissue harboring millions of conidia that serve as secondary inoculum.
these conidia are distributed by splashing and wind-dispersed raindrops which ensures efficient dispersal and hence spread of the disease
<cit> .
measures to combat fungal infections require detailed knowledge of the fungal infection biology.
although some aspects of the pathogenic development of c. graminicola are well-understood, the molecular basis underlying the fundamental morphological and physiological switches that are indispensable for virulence of c. graminicola is still elusive.
large-scale genomics and transcriptomics approaches have been proven to enable the identification of pathogenicity genes in fungi affecting humans or plants
<cit> , which makes them valuable tools towards a better understanding of fungal infection mechanisms.
a first step towards the application of such methods in c. graminicola research was the release of the annotated genome assembly in  <dig> by the colletotrichum sequencing project group .
the genome of c. graminicola was sequenced with a coverage of about 9-fold using a combination of  <dig> and sanger whole genome shotgun methodology.
paired-end reads from  <dig>  plasmids and  <dig>  fosmids were employed to improve the continuity of the assembly.
protein-coding genes were annotated using multiple lines of evidence from blast, pfam searches and est alignments, and gene structures were predicted using the broad institute automated gene-calling pipeline as described
<cit> .
further analysis revealed an exceptionally large and diverse inventory of cell wall-degrading enzymes and secondary metabolism enzymes.
moreover,  <dig> genus-specific genes were identified that potentially encode effector proteins.
based on data obtained from optical mapping,  <dig> % of the genome could not be matched to chromosomes due to extensive repetitive sequences that were estimated to comprise  <dig> % of the c. graminicola genome.
the size of the assembled genome of c. graminicola was found to be  <dig>  mb with a total of  <dig>  genes,  <dig> trnas and  <dig> rrnas.
the identity of  <dig>  genes has been confirmed by sequencing of  <dig>  ests; however, the majority of gene models still lack experimental support
<cit> .
predictions of gene models that are predominantly based on automated algorithms have been shown to be often incorrect or partially correct, especially for species-specific and non-conserved genes
<cit> .
accordingly, using traditional bioinformatics tools makes it difficult to identify novel genes, utrs, intron boundaries and alternative splicing.
since the quality of further comparative and functional genomics or transcriptomics studies depends largely on the correctness of predicted gene models, large-scale transcript information is required to improve the current annotation of the c. graminicola genome.
in this context, next generation sequencing technology has become a powerful tool with the potential to provide a global view not only of the genes present in c. graminicola, but also of their transcriptional and post-transcriptional regulation at single-nucleotide resolution.
in fact, application of transcript data derived from rna-seq approaches offers the opportunity to identify sophisticated mechanisms of gene regulation, e. g. alternative splicing, rna editing or the presence of regulatory rnas
<cit> .
recent advances in gene model prediction based on the application of rna-seq have been published in a number of research articles.
for example, in sordaria macrospora both utr regions of about  <dig> of the predicted genes could be identified, in addition to approx.
<dig> gene models that were improved or newly annotated
<cit> .
using rna-seq data, zhao and coworkers
<cit>  identified  <dig> novel transcriptionally active regions and revised  <dig> incorrectly predicted gene models in fusarium graminearum.
moreover,  <dig> genes were identified with two or more alternative splice variants.
a similar approach by wang et al.
<cit>  yielded approx.
<dig> novel protein-coding genes,  <dig> new candidate exons, utr regions of more than  <dig> genes and as many as  <dig> genes differentially spliced in aspergillus oryzae.
these numbers impressively show the additional benefit provided by a transcript-based revision of the genome annotation.
in order to improve the predicted gene models in c. graminicola, we used an extensive set of rna-seq data covering the transcriptome of this fungus during pathogenic development in maize leaves.
in this study, we revised predicted gene models, identified novel genes and ntars, and searched for alternative splicing and rna editing in c. graminicola.
based on the analysis of our rna-seq data, an improved annotation of the c. graminicola genome together with a comprehensive overview of alternatively spliced genes has been made available  in order to promote colletotrichum research.
RESULTS
quality analysis of rna-seq data
we totally obtained  <dig>  billion reads of a length of  <dig> nucleotides.
of these, approx.
<dig> million quality-processed reads with an average length of  <dig> nucleotides mapped to the genome of c. graminicola.
this significant difference can be explained by the relatively low amount of fungal biomass that is present in the host, especially during the early stages of infection.
since about 67% of the c. graminicola genome contains protein-coding sequence information, the average coverage of these regions was approx.
122-fold.
analysis of the localization and distribution of the obtained reads in the genomic context revealed that among the matched reads 93% were uniquely mapped to the genome with their full read length allowing two mismatches.
about 7% mapped to multiple locations in the genome, of which 96% matched  <dig> to  <dig> different locations, and 4% to more than  <dig> different locations.
rna-seq data showed that  <dig> % of the reads that matched unique regions could be mapped to known coding/utr regions,  <dig> % to intergenic/unknown utr sequences, and  <dig> % mapped to genomic regions encoding rrna and trna .
of all reads that mapped to multiple locations,  <dig> % mapped to genomic regions encoding rrna and trna,  <dig> % mapped to known utr/protein coding regions, and  <dig> % mapped to intergenic regions and yet unidentified utrs .figure  <dig>
statistics and analysis of the quality of rna-seq data.
a. distribution of mapped reads over exons, introns, intergenic and untranslated regions .
about 93% of all reads mapped to unique locations in the genome of c. graminicola, of which  <dig> % matched coding sequences  and identified utrs,  <dig> % to intergenic sequences and unidentified utrs, and  <dig> % to sequences encoding rrna/trna.
b. the remaining 7% of reads mapped to multiple locations, of which  <dig> % matched sequences encoding rrna/trna,  <dig> % to cds and identified utrs, and  <dig> % to intergenic sequences and unidentified utrs.
the read numbers of each category are given in brackets.
c. total coverage of cds by rna-seq reads.
the cds were divided into  <dig> equal windows.
d. gene coverage vs. sequencing depth.
approx.
87% of c. graminicola genes were detected at about  <dig>  million uniquely mapped rna-seq reads, and coverage reaches a plateau afterwards despite the increasing sequencing depth.
e-f. scatter plot analysis of rna-seq data of biological repeats  <dig> vs.  <dig>  and  <dig> vs.  <dig> .
each repeat consists of a combined data set of six individual rna samples taken at  <dig>   <dig>   <dig>   <dig>   <dig>   <dig> hours postinoculation.
each blue point represents a particular gene.
the color shading correlates with the density of these points.
an indispensable prerequisite for reliable rna-seq data is the sufficient coverage of transcripts, preferably by a large number of reads that are evenly distributed across the target sequences
<cit> .
although 3-prime ends of transcripts were found to be represented by a somewhat higher quantity of reads, the overall coverage shows a largely even distribution and sufficient coverage of the coding transcript sequences by rna-seq reads .
altogether, transcripts corresponding to  <dig>   of the  <dig>  genes predicted in the current annotation of the c. graminicola genome could be detected.
this number contrasts the  <dig>  genes that had been previously validated by est contigs
<cit> .
about  <dig>   of all predicted genes were found to be covered in full-length by rna-seq reads.
reads corresponding to about 87% of c. graminicola genes were detected at approx.
<dig>  million uniquely mapped rna-seq reads.
subsequently, gene coverage increased only moderately despite the increasing sequencing depth .
further quality evaluation was performed in order to assess the reproducibility of biological repeats.
merged data sets of biological repeats  <dig> vs.  <dig> and  <dig> vs.  <dig> show consistent distribution of data points in the scatter plots .
these results clearly demonstrate the good quality of the obtained rna-seq data as proven by the extensive coverage of hitherto identified genes, sufficient depth of the sequencing approach and the broad diversity of rnas transcribed during pathogenic development of c. graminicola, which is a prerequisite for both the identification of a large number of ntars and improvement of the annotated gene models.
identification of novel gene models and non-coding rnas
in c. graminicola, a total of  <dig>  genes, encoding proteins of more than  <dig> amino acids, have already been predicted by gene-finding software tools
<cit> .
however,  <dig> % of the mapped reads fell within intergenic/unidentified utr regions.
a substantial part of these reads correspond presumably to non-coding rnas and transcripts of as yet unrecognized genes or exons of incorrectly annotated genes.
altogether, we found  <dig> ntars, of which  <dig> harbor introns indicating that they could be genes.
in addition, about 60% of all ntars were longer than  <dig> bp .
conceptional 6-frame translations employed to search the non-redundant protein database of ncbi  and the protein domain database of the sanger-institute  revealed the presence of coding sequences with similarity to known proteins in  <dig> ntars , of which  <dig> displayed a conserved protein domain .
notably, 74% of these putative proteins were found to have orthologs in other colletotrichum species .figure  <dig>
identification of novel gene models and non-coding rnas.
a. of the  <dig> novel transcriptionally active regions  identified,  <dig> ntars display similarity with non-coding rnas and  <dig> with known proteins .
b. about 80% of these share significant identity with proteins from other colletotrichum species.
c. one example of a novel transcript that has not been identified by automated annotation of the genomic sequence.
blast analysis revealed a protein of  <dig> amino acids with similarity  to the putative effector chec <dig> from colletotrichum higginsianum.
the green bar highlights the cds of the chec <dig> ortholog, grey bars represent individual rna-seq reads.
d. transcript profile of the chec <dig> ortholog in c. graminicola.
transcriptional activity at 24- <dig> hours postinoculation  correlates with necrotrophic development of c. graminicola in maize leaves.
the majority of ntars could not be assigned to particular genes, probably due to the fact that they either lack conserved orthologous genes in other species or these genes have not been annotated so far.
alternatively, these tars could correspond to non-coding or regulatory rnas.
in order to identify non-coding/regulatory rnas we compared database entries of known non-coding rnas with nucleotide sequences of the c. graminicola ntars.
despite the enrichment of mrna sequences carrying terminal poly-tails during sample preparation we found  <dig> ntars with similarity to non-coding rnas, i. e.  <dig> transfer rnas,  <dig> spliceosomal rna, <dig> small nuclear rna and  <dig> ribosomal rna .
none of the ntars shared similarity with database entries of known regulatory rnas.
one reason for this could be that regulatory rnas are species-specific, which would prevent identification of similar rnas in public databases containing predominantly entries of mammals and drosophila melanogaster.
though the function of the majority of ntars remains elusive,  <dig> of the ntars that were found to have orthologs in other fungal species display intact full-length orfs .
these results convincingly show that a substantial proportion of ntars are most likely transcripts of protein-encoding genes rather than pseudogenes or regulatory rnas.
as a consequence, a significant number of genes, some of which with putative functions in virulence or pathogenicity, have not been recorded in the current version of the genome annotation.
for example, an ntar designated cogr_ <dig>   encodes a protein with assumed effector function as evidenced by the occurrence of an ortholog with putative effector function in c. higginsianum
<cit> .
transcripts of cogr_ <dig>  were detected exclusively between 24- <dig> hpi , which could be indicative of a role in the switch from the biotrophic to the necrotrophic lifestyle.
identification of incorrect gene models
following mapping and assembly of the rna-seq reads, we conducted a thorough manual inspection of the currently annotated genes using the integrative genome viewer
<cit>  that allowed simultaneous survey of mapped reads, genome annotation of the broad institute and a genome annotation based on the algorithms of augustus, a gene prediction software dedicated to the analysis of fungal genomes that has been trained using the supplied fusarium graminearum data set
<cit> .
based on the coverage of the genomic sequence by rna-seq reads the following categories of incorrectly annotated features have been identified:figure  <dig>
identification of incorrect gene models in the
c. graminicola
database.
a. statistical analysis of incorrect gene models.
in total,  <dig> gene models were improved with respect to incorrect introns or exons, incorrect splice sites, incorrectly split or joined genes and sequencing errors.
b-e. examples of the four categories of incorrectly annotated genes.
b. read coverage of the coding sequence of glrg_ <dig> convincingly shows that the first intron is somewhat longer than predicted.
c. exon five of glrg_ <dig> harbours an unidentified intron.
using rna-seq data, an additional short exon at the immediate 3′-terminus needed to be added to the gene structure.
d. gene glrg_ <dig> and glrg_ <dig> were incorrectly split due to the occurrence of two introns .
the correct annotation shows a single orf that consists of four exons .
e. insertion of an adenine into the genomic dna sequence of glrg_ <dig> led to frame-shift and premature stop of translation.
in the current annotation, the conflict was apparently solved by inserting an intron which could be identified as incorrect.
blue bars indicate incorrect cds predicted by automated gene annotation, green bars highlight the re-annotated gene structure based on rna-seq data, grey bars represent individual rna-seq reads and red arrowheads indicate the positions of errors in the current annotation of the c. graminicola genome.
i. incorrect intron-exon boundaries: this type of incorrect annotation was identified either by reads that mapped partially to annotated intron regions or partial missing coverage of annotated exon regions.
prerequisite for the correction of these errors was a sufficient read coverage of the genes to be re-annotated.
in comparison with the current annotation we found  <dig> introns that were too long and  <dig> that were too short .
ii.
incorrect exons/introns: the procedure that led to the identification of this type of error was similar to the approach chosen to identify errors of type i with the exception that only absent or redundant full-length introns were captured.
at least one intron and one exon needed to be added to the sequence of  <dig> genes whereas one intron and one exon was removed in  <dig> genes .
iii.
incorrectly split or joined genes: this category refers to adjacent genes that were incorrectly joined in the annotation of the broad institute or individual genes that were incorrectly split due to misinterpretation of the genomic sequence data.
incorrectly joined genes were identified by the lack of read coverage in the sequence of the putative gene which could not be explained by the occurrence of an intron.
vice versa, incorrectly split genes could be found by a largely uniform read coverage throughout the sequence of the incorrectly annotated gene.
in both cases, newly annotated genes were subjected to blast search in order to confirm the identity of the new gene by sequence comparisons with orthologous genes of closely-related fungal species.
in total,  <dig> incorrectly joined genes and  <dig> incorrectly split genes could be identified .
iv.
incorrect genomic sequence: by comparing the genomic nucleotide sequence provided by the broad institute with our rna-seq and genomic sequence data we were able to detect sequencing errors or point mutations that led to incorrect gene models .
among the  <dig> aberrations identified,  <dig> did not alter the genetic information of coding sequences, i. e. they were silent nucleotide substitutions, insertions or deletions within intergenic regions and nucleotide substitutions located either upstream or downstream of coding sequences or within tars .
the gene structure or sequence of the translation product was affected in  <dig> genes by nucleotide substitutions leading to missense, loss of the stop codon and insertions or deletions of nucleotides in coding sequences .
in some cases, these errors had fatal effects with respect to the quality of the algorithm-based gene annotation.
for example, an insertion of adenine at nucleotide position  <dig> of gene glrg_ <dig> encoding a cyclohexanone monooxygenase resulted in a frameshift and hence stop of translation.
automated annotation by the broad institute led to the incorrect insertion of a short intron in order to solve the conflict.
the correct orf was restored by removing the adenine  as proven by blastp homology search against the ncbi database.
altogether, incorrect annotations of  <dig> genes originating from substitutions, deletions or insertions of single nucleotides could be improved by manual revision of gene models.
using rna-seq data, a total of  <dig> incorrect gene models identified in the current genome annotation were successfully revised, which will significantly contribute to both increased quality and reliability of the genomic data of c. graminicola.
identification of alternative splicing
alternative splicing describes the regulated processing of mrnas that leads to the expression of multiple proteins encoded by a single gene, which is an inevitable prerequisite for diversity and functional complexity of eukaryotic proteomes
.
in order to assess alternative splicing in c. graminicola we inspected the re-annotated genome with special emphasis on reads that match intronic regions.
candidate genes for alternative splicing were identified by a significant read coverage within intron sequences as compared to the read coverage of the surrounding exons.
we found  <dig> putative as events in  <dig> genes, including exon skipping, intron retention and alternative 5′ or 3′ splice sites .figure  <dig>
identification of alternative splicing.
a. overview of types and frequency of putative alternative splicing events in c. graminicola.
b. classification of alternatively spliced genes according to their functions.
c. alternative 3′- or  5′- splice sites were identified by a significant number of reads covering intron sequences at their 3′- or 5′-ends whereas intron retention  could be identified by reads covering the entire intron sequence.
exon skipping  was identified by a significant number of reads that cover predicted intron sequences, but do not bridge the gap between the alternate exon and adjacent exons.
blue bars visualize the structure of gene models predicted in the current annotation of the c. graminicola genome.
manually improved alternate gene model are illustrated by green bars, grey bars represent individual rna-seq reads.
intron retention appears to be the most prevalent form of as in c. graminicola, accounting for about 56% of all as events.
however, analysis of the alternate orfs revealed premature termination of translation in most of the candidate genes irrespective of the as mechanism.
there was no evidence for a stage-specific exchange of entire domains, a characteristic effect of as in higher organisms that leads to functionally diverse proteins.
the majority of the as candidate genes that have been characterized are involved in transcription or translation.
further putative functions include maintenance of dna, signal perception or transmission and cellular structure or transport .
in total, about  <dig> % of all c. graminicola genes were estimated to undergo as, though it must be emphasized that delayed processing of premature mrnas may mimic alternative splicing
<cit> .
in line with this, detection of alternatively spliced rnas by rt-pcr, which has been frequently used, does not appear to be an appropriate technique providing sufficient evidence for alternative splicing.
identification of untranslated regions
we identified the 5′- and 3′-boundaries of transcripts by searching for a sharp decline of rna-seq reads signals at both ends of the annotated transcripts.
utrs of genes that overlap with utrs of adjacent genes were excluded from further analysis.
altogether, 5′- and 3′-utrs of  <dig> genes could be retrieved .
the average length of utrs in the genome of c. graminicola was found to be  <dig> nucleotides.
utrs of the 5′- and 3′-regions have an average length of  <dig> and  <dig> nucleotides, respectively.
among the identified utrs,  <dig> harbor intronic sequences, of which  <dig> are located in the 5′-utrs and  <dig> in the 3′-utrs.
screening for rna editing
comparing the genomic sequence and rna-seq data obtained from c. graminicola strain m <dig> offered the opportunity to identify nucleotide modifications that occurred due to rna editing.
a genome-wide search for snps, deletions and insertions  of single nucleotides yielded  <dig> of these modifications, however, careful manual examination revealed that the majority of putative rna editing-based substitutions or insertions were restricted to repetitive sequences of a particular nucleotide or located near intron splice sites.
therefore, we concluded that single nucleotides were erroneously introduced or deleted by dna polymerases used in the sequencing reactions rather than by rna editing mechanisms.
in addition, several snps appeared to be caused by misalignment of rna-seq reads to the genomic dna sequence.
in order to verify the nucleotide sequence of regions harboring putative snps we performed local re-sequencing of the genomic dna.
the obtained results strongly suggest that no rna editing occurs in c. graminicola.
discussion
in the present study we employed an extensive set of rna-seq data covering the infection-related transcriptome of the maize pathogen c. graminicola to evaluate the correctness and completeness of the predicted gene models that are present in the annotated genome database hosted by the broad institute.
moreover, detailed sequence information of transcripts obtained by rna-seq enabled the investigation of post-transcriptional processing and regulation in c. graminicola, such as rna editing or alternative splicing.
deep sequencing techniques provide considerable advantage over microarrays that have been frequently used in transcriptome analysis
.
among these advantages reduced costs, increased sensitivity and detection range as well as the absence of cross-hybridization are the most prevalent.
using rna-seq, transcript abundances can be determined with high accuracy over several orders of magnitude based essentially on the number of reads covering a particular cdna.
however, extensive pcr steps during sample preparation and sequencing may partially bias the obtained results.
this is particularly true for transcripts containing nucleotide repeats or gc-rich regions that interfere with polymerase activity
<cit> .
rna-seq has been successfully applied to improve genome annotations of several species
<cit> , though transcriptomics of sample material comprising the transcriptome of host and pathogen is still a challenge with respect to correct mapping of reads and the significant disproportion in the amount of reads derived from host and pathogen.
the latter applies especially to fungal pathogens with a hemibiotrophic lifestyle which is featured by little or no increase of fungal biomass during the early stages of infection.
despite the excision of inoculated leaf areas as a method to increase the percentage of fungal transcripts, we observed approx.
100-fold more plant than fungal transcripts until  <dig> hpi.
in contrast, later time points showed a significantly improved ratio due to progressive colonization of leaf tissue and spread of infection during necrotrophic development of c. graminicola
<cit> .
the relative low number of about  <dig> million reads that mapped to the fungal genome compared to the overall number of reads  clearly reflects this relationship.
nevertheless, the excellent coverage of transcriptional active regions of 122-fold as well as the overall distribution of reads across transcripts shows sufficient sequencing depth and consistent 5′ to 3′ coverage of transcripts.
analysis of the read distribution revealed that  <dig> % of the reads mapped to intergenic regions.
this percentage is somewhat higher compared to  <dig> % in f. graminearum
<cit> , 3% in homo sapiens
<cit>  and 5% in arabidopsis thaliana
<cit> .
zhao and associates
<cit>  suggested that a high percentage of intergenic rna-seq reads may at least partly reflect the lower quality of gene model predictions in organisms that have either been sequenced recently or their genomes had not been subjected to several rounds of genome annotation as, for example, the well-studied genomes of h. sapiens or a. thaliana.
about 7% of the rna-seq reads could be matched to multiple locations in the genome.
analysis of these reads revealed predominantly rrna sequences.
in addition,  <dig> % of all quality-processed reads matched to unique positions in rrna-encoding sequences indicating the presence of relative high amounts of reads covering rrna, even after enrichment of mrnas using oligo magnetic beads prior to cdna synthesis.
this may be partly due to the fact that hairpin structures in rrna molecules enabled cdna synthesis.
another reason could be the presence of poly tails that have been proposed to serve as degradation mark in rrnas
<cit> .
in contrast, only very few reads that mapped to multiple locations matched coding regions of genes which indicates excellent specificity of our rna-seq reads with respect to the corresponding target transcript.
therefore, the numeric read coverage of a particular transcript can be regarded as a reliable reflection of the gene expression level.
consequently, rna-seq data collected in the present project can be applied for further studies including analysis of the transcriptional changes during pathogenic development of c. graminicola or comparative transcriptomics using different maize cultivars as hosts.
in c. graminicola,  <dig> % of the rna-seq reads fell within intergenic and intronic regions.
this relatively large number can be explained by the presence of as yet unrecognized gene transcripts or utrs and non-coding rnas.
alternatively, these transcripts correspond to novel exons of known genes which had not been identified as integral part of the orf.
among the  <dig> novel tars,  <dig> were likely representing protein-coding genes as evidenced by the presence of full-length orfs, i. e. there was no stop codon present in the orf that yielded the blast result.
among the remaining ntars, only  <dig> could be identified as transfer, splicesosomal, small nuclear and ribosomal rnas leaving 82% of the ntars without a clear function.
a possible explanation for this discrepancy is that transcripts without a proposed function could be either species-specific or the corresponding transcripts have not been annotated yet in related species.
as a consequence, the role of the vast majority of ntars remains elusive even though many of those represent most likely non-coding rnas with regulatory function, pseudogenes or by-products
<cit> .
in addition to several hundreds of ntars we identified novel candidate exons in  <dig> annotated genes and revised intronic sequences in  <dig> genes.
altogether,  <dig> incorrect gene models were improved based on our rna-seq data, which defines a fraction of  <dig> % of incorrect gene models deposited in the published annotation of c. graminicola.
this relatively high percentage of incorrectly annotated genes convincingly shows both the need for manual revision of genome annotations that rely solely on software-assisted prediction of gene models as well as the potential of rna-seq as a powerful tool for significantly improved genome annotations.
rna-seq data of this study were generated from fungal infection structures sampled at six time points of the infection cycle.
this experimental approach was chosen not only to increase the diversity of fungal transcripts but also to enable the detection of infection stage-specific as.
evidence for as has been found in many organisms including h. sapiens, caenorhabditis elegans, arabidopsis thaliana, f. graminearum and aspergillus oryzae, however, experimental proof that includes sequence data of proteins is scarce for fungal organisms.
whereas about 95% of all genes in h. sapiens undergo alternative splicing
<cit> , the fraction of alternatively spliced genes is 42% in a. thaliana,  <dig> % in a. oryzae,  <dig> % in f. graminearum and  <dig> % in m. oryzae
<cit> .
in general, as events in fungi appear to be restricted to relative few genes compared with higher organisms even though a percentage of  <dig> % alternatively spliced genes in c. graminicola and almost 9% in a. oryzae indicate a surprising variability in related species.
this significant spread may at least partially reflect the dilemma of as analysis.
on the one hand, prediction of as based on transcript sequence data seems to be overestimated in several filamentous fungi as cdna-derived sequence data without supporting protein sequence data do not allow discrimination between as and unprocessed or partially processed transcripts.
on the other hand, the frequency of as tends to be underestimated in some studies due to inappropriate experimental conditions that do not reflect the challenges of the natural habitat.
as has been evolved as a sophisticated mechanism that increases the flexibility of gene expression in order to enable organisms to adapt to various environmental conditions.
experimental approaches dedicated to the analysis of as should therefore reflect this dependency as the percentage of as events will certainly increase under more different growth conditions or stresses.
to meet these requirements, we collected sample material during pathogenic development of c. graminicola.
however, only  <dig> putative as events in  <dig> genes could be detected, most of which reduced the size of the corresponding gene product significantly.
although we applied strict criteria in defining putative as events with respect to the proportion of alternate transcripts, it cannot be excluded that at least some of these candidates constitute unprocessed transcripts rather than isoforms.
moreover, there was no evidence for developmentally-regulated as thus questioning the relevance of as in pathogenic development of c. graminicola.
CONCLUSIONS
we have analyzed the transcriptome of the corn anthracnose fungus c. graminicola during infection of maize leaves by rna-seq.
combined trancriptomes of host and pathogen require increased sequencing depth compared to similar approaches with axenic fungal cultures, however, samples of infection-related rnas display a greater variety of transcripts, which is a precondition for the identification of ntars and the improvement of hitherto annotated gene models.
we obtained an extensive set of rna-seq data covering about 98% of the predicted genes which enabled both the revision of incorrect gene models and the identification of hundreds of ntars.
in addition, we obtained preliminary information on alternative splicing in c. graminicola thereby questioning the frequency and relevance of as in filamentous fungi.
unbiased genome annotations are an indispensable prerequisite for modern research approaches.
combining automated gene annotations and manual revisions of gene models based on sequence data from genomic dna and mrna, which has been shown to increase the quality of genome annotations significantly, should therefore become a standard procedure in genome research.

 colletotrichum graminicola is a filamentous ascomycete that causes anthracnose leaf blight and stalk rot of maize leading to estimated annual yield losses in the us of about 6%. in addition to its economical importance, c. graminicola is among the best characterized and most tractable fungi of the genus colletotrichum, a genus comprising a broad range of hemibiotrophic plant pathogens that represent a constant threat to fruit and vegetable production worldwide. interesting aspects of colletotrichum pathology are the pressure-driven penetration strategy and the hemibiotrophic lifestyle that is featured by a fundamental switch in nutrition. an initial short biotrophic phase, during which the host cell remains alive, is followed by highly destructive necrotrophic development characterized by extended areas of killed host tissue. following germination of conidia, specialized infection cells called appressoria are differentiated on the host epidermis. accumulation of osmolytes in appressoria of c. graminicola generates a hydrostatic pressure corresponding to approx. mpa, which is among the highest turgor pressures known in living cells. this enormous pressure is utilized by the fungus to forcefully penetrate plant epidermis cells. upon successful penetration of the host epidermis, biotrophic hyphae of c. graminicola initially colonize living plant cells as evidenced by the integrity of the host plasma membrane. this stage of symptomless parasitic development is terminated at the onset of necrotrophy. fast-growing hyphae breach the walls of adjacent cells, ramify within the leave and kill host cells, leading to massive destruction of maize tissue and visible necrosis. finally, acervuli are formed on necrotic tissue harboring millions of conidia that serve as secondary inoculum. these conidia are distributed by splashing and wind-dispersed raindrops which ensures efficient dispersal and hence spread of the disease. measures to combat fungal infections require detailed knowledge of the fungal infection biology. although some aspects of the pathogenic development of c. graminicola are well-understood, the molecular basis underlying the fundamental morphological and physiological switches that are indispensable for virulence of c. graminicola is still elusive. large-scale genomics and transcriptomics approaches have been proven to enable the identification of pathogenicity genes in fungi affecting humans or plants, which makes them valuable tools towards a better understanding of fungal infection mechanisms. a first step towards the application of such methods in c. graminicola research was the release of the annotated genome assembly in by the colletotrichum sequencing project group. the genome of c. graminicola was sequenced with a coverage of about 9-fold using a combination of and sanger whole genome shotgun methodology. paired-end reads from plasmids and fosmids were employed to improve the continuity of the assembly. protein-coding genes were annotated using multiple lines of evidence from blast, pfam searches and est alignments, and gene structures were predicted using the broad institute automated gene-calling pipeline as described. further analysis revealed an exceptionally large and diverse inventory of cell wall-degrading enzymes and secondary metabolism enzymes. moreover, genus-specific genes were identified that potentially encode effector proteins. based on data obtained from optical mapping, % of the genome could not be matched to chromosomes due to extensive repetitive sequences that were estimated to comprise % of the c. graminicola genome. the size of the assembled genome of c. graminicola was found to be mb with a total of genes, trnas and rrnas. the identity of genes has been confirmed by sequencing of ests; however, the majority of gene models still lack experimental support. predictions of gene models that are predominantly based on automated algorithms have been shown to be often incorrect or partially correct, especially for species-specific and non-conserved genes. accordingly, using traditional bioinformatics tools makes it difficult to identify novel genes, utrs, intron boundaries and alternative splicing. since the quality of further comparative and functional genomics or transcriptomics studies depends largely on the correctness of predicted gene models, large-scale transcript information is required to improve the current annotation of the c. graminicola genome. in this context, next generation sequencing technology has become a powerful tool with the potential to provide a global view not only of the genes present in c. graminicola, but also of their transcriptional and post-transcriptional regulation at single-nucleotide resolution. in fact, application of transcript data derived from rna-seq approaches offers the opportunity to identify sophisticated mechanisms of gene regulation, e. g. alternative splicing, rna editing or the presence of regulatory rnas. recent advances in gene model prediction based on the application of rna-seq have been published in a number of research articles. for example, in sordaria macrospora both utr regions of about of the predicted genes could be identified, in addition to approx. gene models that were improved or newly annotated. using rna-seq data, zhao and coworkers identified novel transcriptionally active regions and revised incorrectly predicted gene models in fusarium graminearum. moreover, genes were identified with two or more alternative splice variants. a similar approach by wang et al. yielded approx. novel protein-coding genes, new candidate exons, utr regions of more than genes and as many as genes differentially spliced in aspergillus oryzae. these numbers impressively show the additional benefit provided by a transcript-based revision of the genome annotation. in order to improve the predicted gene models in c. graminicola, we used an extensive set of rna-seq data covering the transcriptome of this fungus during pathogenic development in maize leaves. in this study, we revised predicted gene models, identified novel genes and ntars, and searched for alternative splicing and rna editing in c. graminicola. based on the analysis of our rna-seq data, an improved annotation of the c. graminicola genome together with a comprehensive overview of alternatively spliced genes has been made available in order to promote colletotrichum research.  quality analysis of rna-seq data we totally obtained billion reads of a length of nucleotides. of these, approx. million quality-processed reads with an average length of nucleotides mapped to the genome of c. graminicola. this significant difference can be explained by the relatively low amount of fungal biomass that is present in the host, especially during the early stages of infection. since about 67% of the c. graminicola genome contains protein-coding sequence information, the average coverage of these regions was approx. 122-fold. analysis of the localization and distribution of the obtained reads in the genomic context revealed that among the matched reads 93% were uniquely mapped to the genome with their full read length allowing two mismatches. about 7% mapped to multiple locations in the genome, of which 96% matched to different locations, and 4% to more than different locations. rna-seq data showed that % of the reads that matched unique regions could be mapped to known coding/utr regions, % to intergenic/unknown utr sequences, and % mapped to genomic regions encoding rrna and trna. of all reads that mapped to multiple locations, % mapped to genomic regions encoding rrna and trna, % mapped to known utr/protein coding regions, and % mapped to intergenic regions and yet unidentified utrs.figure statistics and analysis of the quality of rna-seq data. a. distribution of mapped reads over exons, introns, intergenic and untranslated regions. about 93% of all reads mapped to unique locations in the genome of c. graminicola, of which % matched coding sequences and identified utrs, % to intergenic sequences and unidentified utrs, and % to sequences encoding rrna/trna. b. the remaining 7% of reads mapped to multiple locations, of which % matched sequences encoding rrna/trna, % to cds and identified utrs, and % to intergenic sequences and unidentified utrs. the read numbers of each category are given in brackets. c. total coverage of cds by rna-seq reads. the cds were divided into equal windows. d. gene coverage vs. sequencing depth. approx. 87% of c. graminicola genes were detected at about million uniquely mapped rna-seq reads, and coverage reaches a plateau afterwards despite the increasing sequencing depth. e-f. scatter plot analysis of rna-seq data of biological repeats vs. and vs.. each repeat consists of a combined data set of six individual rna samples taken at hours postinoculation. each blue point represents a particular gene. the color shading correlates with the density of these points. an indispensable prerequisite for reliable rna-seq data is the sufficient coverage of transcripts, preferably by a large number of reads that are evenly distributed across the target sequences. although 3-prime ends of transcripts were found to be represented by a somewhat higher quantity of reads, the overall coverage shows a largely even distribution and sufficient coverage of the coding transcript sequences by rna-seq reads. altogether, transcripts corresponding to of the genes predicted in the current annotation of the c. graminicola genome could be detected. this number contrasts the genes that had been previously validated by est contigs. about of all predicted genes were found to be covered in full-length by rna-seq reads. reads corresponding to about 87% of c. graminicola genes were detected at approx. million uniquely mapped rna-seq reads. subsequently, gene coverage increased only moderately despite the increasing sequencing depth. further quality evaluation was performed in order to assess the reproducibility of biological repeats. merged data sets of biological repeats vs. and vs. show consistent distribution of data points in the scatter plots. these results clearly demonstrate the good quality of the obtained rna-seq data as proven by the extensive coverage of hitherto identified genes, sufficient depth of the sequencing approach and the broad diversity of rnas transcribed during pathogenic development of c. graminicola, which is a prerequisite for both the identification of a large number of ntars and improvement of the annotated gene models. identification of novel gene models and non-coding rnas in c. graminicola, a total of genes, encoding proteins of more than amino acids, have already been predicted by gene-finding software tools. however, % of the mapped reads fell within intergenic/unidentified utr regions. a substantial part of these reads correspond presumably to non-coding rnas and transcripts of as yet unrecognized genes or exons of incorrectly annotated genes. altogether, we found ntars, of which harbor introns indicating that they could be genes. in addition, about 60% of all ntars were longer than bp. conceptional 6-frame translations employed to search the non-redundant protein database of ncbi and the protein domain database of the sanger-institute revealed the presence of coding sequences with similarity to known proteins in ntars, of which displayed a conserved protein domain. notably, 74% of these putative proteins were found to have orthologs in other colletotrichum species.figure identification of novel gene models and non-coding rnas. a. of the novel transcriptionally active regions identified, ntars display similarity with non-coding rnas and with known proteins. b. about 80% of these share significant identity with proteins from other colletotrichum species. c. one example of a novel transcript that has not been identified by automated annotation of the genomic sequence. blast analysis revealed a protein of amino acids with similarity to the putative effector chec from colletotrichum higginsianum. the green bar highlights the cds of the chec ortholog, grey bars represent individual rna-seq reads. d. transcript profile of the chec ortholog in c. graminicola. transcriptional activity at 24- hours postinoculation correlates with necrotrophic development of c. graminicola in maize leaves. the majority of ntars could not be assigned to particular genes, probably due to the fact that they either lack conserved orthologous genes in other species or these genes have not been annotated so far. alternatively, these tars could correspond to non-coding or regulatory rnas. in order to identify non-coding/regulatory rnas we compared database entries of known non-coding rnas with nucleotide sequences of the c. graminicola ntars. despite the enrichment of mrna sequences carrying terminal poly-tails during sample preparation we found ntars with similarity to non-coding rnas, i. e. transfer rnas, spliceosomal rna, small nuclear rna and ribosomal rna. none of the ntars shared similarity with database entries of known regulatory rnas. one reason for this could be that regulatory rnas are species-specific, which would prevent identification of similar rnas in public databases containing predominantly entries of mammals and drosophila melanogaster. though the function of the majority of ntars remains elusive, of the ntars that were found to have orthologs in other fungal species display intact full-length orfs. these results convincingly show that a substantial proportion of ntars are most likely transcripts of protein-encoding genes rather than pseudogenes or regulatory rnas. as a consequence, a significant number of genes, some of which with putative functions in virulence or pathogenicity, have not been recorded in the current version of the genome annotation. for example, an ntar designated cogr_ encodes a protein with assumed effector function as evidenced by the occurrence of an ortholog with putative effector function in c. higginsianum. transcripts of cogr_ were detected exclusively between 24- hpi, which could be indicative of a role in the switch from the biotrophic to the necrotrophic lifestyle. identification of incorrect gene models following mapping and assembly of the rna-seq reads, we conducted a thorough manual inspection of the currently annotated genes using the integrative genome viewer that allowed simultaneous survey of mapped reads, genome annotation of the broad institute and a genome annotation based on the algorithms of augustus, a gene prediction software dedicated to the analysis of fungal genomes that has been trained using the supplied fusarium graminearum data set. based on the coverage of the genomic sequence by rna-seq reads the following categories of incorrectly annotated features have been identified:figure identification of incorrect gene models in the c. graminicola database. a. statistical analysis of incorrect gene models. in total, gene models were improved with respect to incorrect introns or exons, incorrect splice sites, incorrectly split or joined genes and sequencing errors. b-e. examples of the four categories of incorrectly annotated genes. b. read coverage of the coding sequence of glrg_ convincingly shows that the first intron is somewhat longer than predicted. c. exon five of glrg_ harbours an unidentified intron. using rna-seq data, an additional short exon at the immediate 3′-terminus needed to be added to the gene structure. d. gene glrg_ and glrg_ were incorrectly split due to the occurrence of two introns. the correct annotation shows a single orf that consists of four exons. e. insertion of an adenine into the genomic dna sequence of glrg_ led to frame-shift and premature stop of translation. in the current annotation, the conflict was apparently solved by inserting an intron which could be identified as incorrect. blue bars indicate incorrect cds predicted by automated gene annotation, green bars highlight the re-annotated gene structure based on rna-seq data, grey bars represent individual rna-seq reads and red arrowheads indicate the positions of errors in the current annotation of the c. graminicola genome. i. incorrect intron-exon boundaries: this type of incorrect annotation was identified either by reads that mapped partially to annotated intron regions or partial missing coverage of annotated exon regions. prerequisite for the correction of these errors was a sufficient read coverage of the genes to be re-annotated. in comparison with the current annotation we found introns that were too long and that were too short. ii. incorrect exons/introns: the procedure that led to the identification of this type of error was similar to the approach chosen to identify errors of type i with the exception that only absent or redundant full-length introns were captured. at least one intron and one exon needed to be added to the sequence of genes whereas one intron and one exon was removed in genes. iii. incorrectly split or joined genes: this category refers to adjacent genes that were incorrectly joined in the annotation of the broad institute or individual genes that were incorrectly split due to misinterpretation of the genomic sequence data. incorrectly joined genes were identified by the lack of read coverage in the sequence of the putative gene which could not be explained by the occurrence of an intron. vice versa, incorrectly split genes could be found by a largely uniform read coverage throughout the sequence of the incorrectly annotated gene. in both cases, newly annotated genes were subjected to blast search in order to confirm the identity of the new gene by sequence comparisons with orthologous genes of closely-related fungal species. in total, incorrectly joined genes and incorrectly split genes could be identified. iv. incorrect genomic sequence: by comparing the genomic nucleotide sequence provided by the broad institute with our rna-seq and genomic sequence data we were able to detect sequencing errors or point mutations that led to incorrect gene models. among the aberrations identified, did not alter the genetic information of coding sequences, i. e. they were silent nucleotide substitutions, insertions or deletions within intergenic regions and nucleotide substitutions located either upstream or downstream of coding sequences or within tars. the gene structure or sequence of the translation product was affected in genes by nucleotide substitutions leading to missense, loss of the stop codon and insertions or deletions of nucleotides in coding sequences. in some cases, these errors had fatal effects with respect to the quality of the algorithm-based gene annotation. for example, an insertion of adenine at nucleotide position of gene glrg_ encoding a cyclohexanone monooxygenase resulted in a frameshift and hence stop of translation. automated annotation by the broad institute led to the incorrect insertion of a short intron in order to solve the conflict. the correct orf was restored by removing the adenine as proven by blastp homology search against the ncbi database. altogether, incorrect annotations of genes originating from substitutions, deletions or insertions of single nucleotides could be improved by manual revision of gene models. using rna-seq data, a total of incorrect gene models identified in the current genome annotation were successfully revised, which will significantly contribute to both increased quality and reliability of the genomic data of c. graminicola. identification of alternative splicing alternative splicing describes the regulated processing of mrnas that leads to the expression of multiple proteins encoded by a single gene, which is an inevitable prerequisite for diversity and functional complexity of eukaryotic proteomes. in order to assess alternative splicing in c. graminicola we inspected the re-annotated genome with special emphasis on reads that match intronic regions. candidate genes for alternative splicing were identified by a significant read coverage within intron sequences as compared to the read coverage of the surrounding exons. we found putative as events in genes, including exon skipping, intron retention and alternative 5′ or 3′ splice sites.figure identification of alternative splicing. a. overview of types and frequency of putative alternative splicing events in c. graminicola. b. classification of alternatively spliced genes according to their functions. c. alternative 3′- or 5′- splice sites were identified by a significant number of reads covering intron sequences at their 3′- or 5′-ends whereas intron retention could be identified by reads covering the entire intron sequence. exon skipping was identified by a significant number of reads that cover predicted intron sequences, but do not bridge the gap between the alternate exon and adjacent exons. blue bars visualize the structure of gene models predicted in the current annotation of the c. graminicola genome. manually improved alternate gene model are illustrated by green bars, grey bars represent individual rna-seq reads. intron retention appears to be the most prevalent form of as in c. graminicola, accounting for about 56% of all as events. however, analysis of the alternate orfs revealed premature termination of translation in most of the candidate genes irrespective of the as mechanism. there was no evidence for a stage-specific exchange of entire domains, a characteristic effect of as in higher organisms that leads to functionally diverse proteins. the majority of the as candidate genes that have been characterized are involved in transcription or translation. further putative functions include maintenance of dna, signal perception or transmission and cellular structure or transport. in total, about % of all c. graminicola genes were estimated to undergo as, though it must be emphasized that delayed processing of premature mrnas may mimic alternative splicing. in line with this, detection of alternatively spliced rnas by rt-pcr, which has been frequently used, does not appear to be an appropriate technique providing sufficient evidence for alternative splicing. identification of untranslated regions we identified the 5′- and 3′-boundaries of transcripts by searching for a sharp decline of rna-seq reads signals at both ends of the annotated transcripts. utrs of genes that overlap with utrs of adjacent genes were excluded from further analysis. altogether, 5′- and 3′-utrs of genes could be retrieved. the average length of utrs in the genome of c. graminicola was found to be nucleotides. utrs of the 5′- and 3′-regions have an average length of and nucleotides, respectively. among the identified utrs, harbor intronic sequences, of which are located in the 5′-utrs and in the 3′-utrs. screening for rna editing comparing the genomic sequence and rna-seq data obtained from c. graminicola strain m offered the opportunity to identify nucleotide modifications that occurred due to rna editing. a genome-wide search for snps, deletions and insertions of single nucleotides yielded of these modifications, however, careful manual examination revealed that the majority of putative rna editing-based substitutions or insertions were restricted to repetitive sequences of a particular nucleotide or located near intron splice sites. therefore, we concluded that single nucleotides were erroneously introduced or deleted by dna polymerases used in the sequencing reactions rather than by rna editing mechanisms. in addition, several snps appeared to be caused by misalignment of rna-seq reads to the genomic dna sequence. in order to verify the nucleotide sequence of regions harboring putative snps we performed local re-sequencing of the genomic dna. the obtained results strongly suggest that no rna editing occurs in c. graminicola. discussion in the present study we employed an extensive set of rna-seq data covering the infection-related transcriptome of the maize pathogen c. graminicola to evaluate the correctness and completeness of the predicted gene models that are present in the annotated genome database hosted by the broad institute. moreover, detailed sequence information of transcripts obtained by rna-seq enabled the investigation of post-transcriptional processing and regulation in c. graminicola, such as rna editing or alternative splicing. deep sequencing techniques provide considerable advantage over microarrays that have been frequently used in transcriptome analysis. among these advantages reduced costs, increased sensitivity and detection range as well as the absence of cross-hybridization are the most prevalent. using rna-seq, transcript abundances can be determined with high accuracy over several orders of magnitude based essentially on the number of reads covering a particular cdna. however, extensive pcr steps during sample preparation and sequencing may partially bias the obtained results. this is particularly true for transcripts containing nucleotide repeats or gc-rich regions that interfere with polymerase activity. rna-seq has been successfully applied to improve genome annotations of several species, though transcriptomics of sample material comprising the transcriptome of host and pathogen is still a challenge with respect to correct mapping of reads and the significant disproportion in the amount of reads derived from host and pathogen. the latter applies especially to fungal pathogens with a hemibiotrophic lifestyle which is featured by little or no increase of fungal biomass during the early stages of infection. despite the excision of inoculated leaf areas as a method to increase the percentage of fungal transcripts, we observed approx. 100-fold more plant than fungal transcripts until hpi. in contrast, later time points showed a significantly improved ratio due to progressive colonization of leaf tissue and spread of infection during necrotrophic development of c. graminicola. the relative low number of about million reads that mapped to the fungal genome compared to the overall number of reads clearly reflects this relationship. nevertheless, the excellent coverage of transcriptional active regions of 122-fold as well as the overall distribution of reads across transcripts shows sufficient sequencing depth and consistent 5′ to 3′ coverage of transcripts. analysis of the read distribution revealed that % of the reads mapped to intergenic regions. this percentage is somewhat higher compared to % in f. graminearum, 3% in homo sapiens and 5% in arabidopsis thaliana. zhao and associates suggested that a high percentage of intergenic rna-seq reads may at least partly reflect the lower quality of gene model predictions in organisms that have either been sequenced recently or their genomes had not been subjected to several rounds of genome annotation as, for example, the well-studied genomes of h. sapiens or a. thaliana. about 7% of the rna-seq reads could be matched to multiple locations in the genome. analysis of these reads revealed predominantly rrna sequences. in addition, % of all quality-processed reads matched to unique positions in rrna-encoding sequences indicating the presence of relative high amounts of reads covering rrna, even after enrichment of mrnas using oligo magnetic beads prior to cdna synthesis. this may be partly due to the fact that hairpin structures in rrna molecules enabled cdna synthesis. another reason could be the presence of poly tails that have been proposed to serve as degradation mark in rrnas. in contrast, only very few reads that mapped to multiple locations matched coding regions of genes which indicates excellent specificity of our rna-seq reads with respect to the corresponding target transcript. therefore, the numeric read coverage of a particular transcript can be regarded as a reliable reflection of the gene expression level. consequently, rna-seq data collected in the present project can be applied for further studies including analysis of the transcriptional changes during pathogenic development of c. graminicola or comparative transcriptomics using different maize cultivars as hosts. in c. graminicola, % of the rna-seq reads fell within intergenic and intronic regions. this relatively large number can be explained by the presence of as yet unrecognized gene transcripts or utrs and non-coding rnas. alternatively, these transcripts correspond to novel exons of known genes which had not been identified as integral part of the orf. among the novel tars, were likely representing protein-coding genes as evidenced by the presence of full-length orfs, i. e. there was no stop codon present in the orf that yielded the blast result. among the remaining ntars, only could be identified as transfer, splicesosomal, small nuclear and ribosomal rnas leaving 82% of the ntars without a clear function. a possible explanation for this discrepancy is that transcripts without a proposed function could be either species-specific or the corresponding transcripts have not been annotated yet in related species. as a consequence, the role of the vast majority of ntars remains elusive even though many of those represent most likely non-coding rnas with regulatory function, pseudogenes or by-products. in addition to several hundreds of ntars we identified novel candidate exons in annotated genes and revised intronic sequences in genes. altogether, incorrect gene models were improved based on our rna-seq data, which defines a fraction of % of incorrect gene models deposited in the published annotation of c. graminicola. this relatively high percentage of incorrectly annotated genes convincingly shows both the need for manual revision of genome annotations that rely solely on software-assisted prediction of gene models as well as the potential of rna-seq as a powerful tool for significantly improved genome annotations. rna-seq data of this study were generated from fungal infection structures sampled at six time points of the infection cycle. this experimental approach was chosen not only to increase the diversity of fungal transcripts but also to enable the detection of infection stage-specific as. evidence for as has been found in many organisms including h. sapiens, caenorhabditis elegans, arabidopsis thaliana, f. graminearum and aspergillus oryzae, however, experimental proof that includes sequence data of proteins is scarce for fungal organisms. whereas about 95% of all genes in h. sapiens undergo alternative splicing, the fraction of alternatively spliced genes is 42% in a. thaliana, % in a. oryzae, % in f. graminearum and % in m. oryzae. in general, as events in fungi appear to be restricted to relative few genes compared with higher organisms even though a percentage of % alternatively spliced genes in c. graminicola and almost 9% in a. oryzae indicate a surprising variability in related species. this significant spread may at least partially reflect the dilemma of as analysis. on the one hand, prediction of as based on transcript sequence data seems to be overestimated in several filamentous fungi as cdna-derived sequence data without supporting protein sequence data do not allow discrimination between as and unprocessed or partially processed transcripts. on the other hand, the frequency of as tends to be underestimated in some studies due to inappropriate experimental conditions that do not reflect the challenges of the natural habitat. as has been evolved as a sophisticated mechanism that increases the flexibility of gene expression in order to enable organisms to adapt to various environmental conditions. experimental approaches dedicated to the analysis of as should therefore reflect this dependency as the percentage of as events will certainly increase under more different growth conditions or stresses. to meet these requirements, we collected sample material during pathogenic development of c. graminicola. however, only putative as events in genes could be detected, most of which reduced the size of the corresponding gene product significantly. although we applied strict criteria in defining putative as events with respect to the proportion of alternate transcripts, it cannot be excluded that at least some of these candidates constitute unprocessed transcripts rather than isoforms. moreover, there was no evidence for developmentally-regulated as thus questioning the relevance of as in pathogenic development of c. graminicola.  we have analyzed the transcriptome of the corn anthracnose fungus c. graminicola during infection of maize leaves by rna-seq. combined trancriptomes of host and pathogen require increased sequencing depth compared to similar approaches with axenic fungal cultures, however, samples of infection-related rnas display a greater variety of transcripts, which is a precondition for the identification of ntars and the improvement of hitherto annotated gene models. we obtained an extensive set of rna-seq data covering about 98% of the predicted genes which enabled both the revision of incorrect gene models and the identification of hundreds of ntars. in addition, we obtained preliminary information on alternative splicing in c. graminicola thereby questioning the frequency and relevance of as in filamentous fungi. unbiased genome annotations are an indispensable prerequisite for modern research approaches. combining automated gene annotations and manual revisions of gene models based on sequence data from genomic dna and mrna, which has been shown to increase the quality of genome annotations significantly, should therefore become a standard procedure in genome research. 
our results significantly contribute to both routine laboratory experiments and large-scale genomics or transcriptomic studies in c. graminicola. an annotated genomic sequence of the corn anthracnose fungus colletotrichum graminicola has been published previously, but correct identification of gene models by means of automated gene annotation remains a challenge. among the ntars, share significant similarity with proteins that have been identified in other species suggesting that they are hitherto unidentified genes in c. graminicola. comparative analysis of rna-seq data and the genome sequence did not provide evidence for rna editing in c. graminicola. we successfully employed deep sequencing rna-seq data in combination with an elaborate bioinformatics strategy in order to identify novel genes, incorrect gene models and mechanisms of transcript processing in the corn anthracnose fungus c. graminicola. based on the nucleotide sequence information of transcripts, we identified novel transcriptionally active regions and revised incorrectly predicted gene models, including revisions of exon-intron structure, gene orientation and sequencing errors. sequence data of the revised genome annotation including several hundreds of novel transcripts, improved gene models and candidate genes for alternative splicing have been made accessible in a comprehensive database. colletotrichum graminicolaanthracnose of cornrna-seqgenome annotationissue-copyright-statement© the author 2014 rna-seq offers the potential for substantially improved gene annotations and for the identification of posttranscriptional rna modifications, such as alternative splicing and rna editing.
 an annotated genomic sequence of the corn anthracnose fungus colletotrichum graminicola has been published previously, but correct identification of gene models by means of automated gene annotation remains a challenge. rna-seq offers the potential for substantially improved gene annotations and for the identification of posttranscriptional rna modifications, such as alternative splicing and rna editing.  based on the nucleotide sequence information of transcripts, we identified novel transcriptionally active regions and revised incorrectly predicted gene models, including revisions of exon-intron structure, gene orientation and sequencing errors. among the ntars, share significant similarity with proteins that have been identified in other species suggesting that they are hitherto unidentified genes in c. graminicola. moreover, 5′- and 3′-utr sequences of genes have been retrieved and alternatively spliced variants of genes have been identified. comparative analysis of rna-seq data and the genome sequence did not provide evidence for rna editing in c. graminicola.  we successfully employed deep sequencing rna-seq data in combination with an elaborate bioinformatics strategy in order to identify novel genes, incorrect gene models and mechanisms of transcript processing in the corn anthracnose fungus c. graminicola. sequence data of the revised genome annotation including several hundreds of novel transcripts, improved gene models and candidate genes for alternative splicing have been made accessible in a comprehensive database. our results significantly contribute to both routine laboratory experiments and large-scale genomics or transcriptomic studies in c. graminicola. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords colletotrichum graminicolaanthracnose of cornrna-seqgenome annotationissue-copyright-statement© the author 2014 