Global analyses of Ceratocystis cacaofunesta mitochondria: from genome to proteome
- Alinne Batista Ambrosio1,
- Leandro Costa do Nascimento1,
- Bruno V Oliveira1,
- Paulo José P L Teixeira1,
- Ricardo A Tiburcio1,
- Daniela P Toledo Thomazella1,
- Adriana F P Leme2,
- Marcelo F Carazzolle1, 3,
- Ramon O Vidal1, 2,
- Piotr Mieczkowski4,
- Lyndel W Meinhardt5,
- Gonçalo A G Pereira1Email author and
- Odalys G Cabrera1
© Ambrosio et al.; licensee BioMed Central Ltd. 2013
Received: 9 August 2012
Accepted: 27 January 2013
Published: 11 February 2013
The ascomycete fungus Ceratocystis cacaofunesta is the causal agent of wilt disease in cacao, which results in significant economic losses in the affected producing areas. Despite the economic importance of the Ceratocystis complex of species, no genomic data are available for any of its members. Given that mitochondria play important roles in fungal virulence and the susceptibility/resistance of fungi to fungicides, we performed the first functional analysis of this organelle in Ceratocystis using integrated “omics” approaches.
The C. cacaofunesta mitochondrial genome (mtDNA) consists of a single, 103,147-bp circular molecule, making this the second largest mtDNA among the Sordariomycetes. Bioinformatics analysis revealed the presence of 15 conserved genes and 37 intronic open reading frames in C. cacaofunesta mtDNA. Here, we predicted the mitochondrial proteome (mtProt) of C. cacaofunesta, which is comprised of 1,124 polypeptides - 52 proteins that are mitochondrially encoded and 1,072 that are nuclearly encoded. Transcriptome analysis revealed 33 probable novel genes. Comparisons among the Gene Ontology results of the predicted mtProt of C. cacaofunesta, Neurospora crassa and Saccharomyces cerevisiae revealed no significant differences. Moreover, C. cacaofunesta mitochondria were isolated, and the mtProt was subjected to mass spectrometric analysis. The experimental proteome validated 27% of the predicted mtProt. Our results confirmed the existence of 110 hypothetical proteins and 7 novel proteins of which 83 and 1, respectively, had putative mitochondrial localization.
The present study provides the first partial genomic analysis of a species of the Ceratocystis genus and the first predicted mitochondrial protein inventory of a phytopathogenic fungus. In addition to the known mitochondrial role in pathogenicity, our results demonstrated that the global function analysis of this organelle is similar in pathogenic and non-pathogenic fungi, suggesting that its relevance in the lifestyle of these organisms should be based on a small number of specific proteins and/or with respect to differential gene regulation. In this regard, particular interest should be directed towards mitochondrial proteins with unknown function and the novel protein that might be specific to this species. Further functional characterization of these proteins could enhance our understanding of the role of mitochondria in phytopathogenicity.
KeywordsCacao wilt disease Ceratocystis cacaofunesta Mitogenomics Mitochondrial proteome Fungal virulence
Ceratocystis cacaofunesta is an ascomycete fungus (Class Sordariomycetes, Order Microascales) that causes wilt disease in cacao (Theobroma cacao). This fungus is a member of the Latin American clade of the Ceratocystis fimbriata complex [1, 2], a taxonomic group that includes species with high genetic variability and wide host ranges [3, 4]. These species cause canker and wilt diseases in many economically important crops, such as Coffea arabica and Eucalyptus spp.. Engelbrecht and Harrington  proposed that host specialization may influence speciation in this group.
C. cacaofunesta is indigenous to Central and South America [1, 2]. In Brazil, this pathogen was first reported in the Amazon region . In the 1990s, C. cacaofunesta was introduced to the southern region of Bahia, which is the largest Brazilian cacao-producing state . This fungus is able to penetrate cacao trees through stem wounds that are caused either by insects (natural vectors) or through infected cutting tools . Unlike other cacao diseases that primarily affect branches and fruits, such as Witches’ broom disease and Frosty pod, wilt disease is a systemic infection that damages the entire plant. The fungus enters its host through the secondary xylem, resulting in the formation of deep spots and leading to the obstruction of water and nutrient transport . Consequently, the plant turns yellow and then brown; the infection culminates in the wilting and sudden death of the tree. This disease is responsible for reductions in the cacao population in plantation areas, which has resulted in great economic losses in the affected regions. Extensive damages have also been reported in Trinidad , Venezuela and Colombia .
Several studies have attempted to characterize the genetic variation, aggressiveness and host specialization of different populations of C. fimbriata complex including C cacaofunesta[6, 13–15]. However, little is known at the molecular level regarding this fungus and its interaction with cacao. Notably, no genomic data are available for any members of the Ceratocystis genus. In 2011, our group initiated the C. cacaofunesta Genome Project (http://www.lge.ibi.unicamp.br/ceratocystis), with the goal of understanding the mechanisms that underlie the interactions between C. cacaofunesta and Theobroma cacao in the development of wilt disease. We initially focused on the study of C. cacaofunesta mitochondria for two reasons: (i) their potential role in fungal pathogenesis and (ii) the relevance of this organelle as a target for fungicides .
In addition to their canonical function as an energy-producing compartment, mitochondria are involved in multiple cellular processes . These organelles play important roles in calcium homeostasis , the biosynthesis of iron-sulfur clusters , lipid and amino acid metabolism, aging and the signaling of programmed cell death . Due to the importance of these processes, mitochondrial dysfunctions cause serious consequences for the cell and, ultimately, the entire organism. Functionally compromised mitochondria are associated with senescence in non-pathogenic fungi, such as Podospora anserina and N. crassa.
Mitogenomics has become a useful tool for evolutionary studies , and the continuous advances in this field have contributed to the understanding of the diverse topology, organization and structure of mtDNA in fungi . However, as mtDNA encodes approximately 1% of the mitochondrial proteome, scant information concerning the roles of mitochondria in the metabolism, development and lifestyle of organisms has been gleaned from exclusive analyses of the mitochondrial genome. Fungal mtDNA generally contains 14 genes that encode hydrophobic subunits of the respiratory chain complexes, two genes that encode ribosomal RNAs (small and large subunits) and genes that encode a full set of tRNAs . The other 99% of the mitochondrial proteome is encoded by nuclear DNA and imported into the mitochondria . Nuclear-encoded mitochondrial proteins (NMP) are produced by cytosolic ribosomes and are targeted to the proper mitochondrial subcompartment [27–29].
To understand mitochondrial function, studies using genomic, transcriptomic and proteomic approaches have been performed in several types of organisms, such as fungi , plants  and humans . The first mitochondrial analysis using integrative proteomic and genomic approaches in fungi was conducted in S. cerevisiae. By overcoming the individual limitations of each technique, this integrative study produced a powerful tool for the prediction of mitochondrial processes in yeast. The S. cerevisiae mitochondrial proteome contains an estimated 1,000 proteins, of which 851 were identified using proteomic assays . In N. crassa, proteomic studies have led to the identification of 438 mitochondrial proteins . Proteomics has previously been used to illuminate central processes in phytopathogenic fungi [36, 37]. Certain more recent studies have analyzed of mycelial, conidiospores and haustoria proteomes [37–41], fungal secretomes and proteomes that are associated with fungal virulence . However, the mitochondrial proteomes of phytopathogenic fungi have not been explored.
In the present study, we performed a global analysis of C. cacaofunesta mitochondria using an integrative approach (mitogenomics, transcriptomics and proteomics). We predicted the C. cacaofunesta mitochondrial proteome (mtProt), including mitochondrial-encoded proteins and mitochondrial proteins encoded in the nuclear genome (NMP). Moreover, 27% of predicted mtProt was validated by experimental proteome analysis (LC-MS/MS). We focused on comparisons with available mitochondrial proteomes from non-pathogenic model fungi to increase our knowledge on the role of mitochondria in pathogenicity. The present study is the first partial genomic analysis of a Ceratocystis species to be published. Additionally, to the best of our knowledge, this study is pioneering in that it presents a global analysis of the mitochondrial proteome of a phytopathogenic fungus using this integrative approach.
Results and discussion
C. cacaofunesta mitochondrial genome structure
The annotation of the C. cacaofunesta mt genome revealed that introns form 48.7% of this DNA molecule, explaining its large size. The invasion of mitochondrial genes by group I introns is a primary reason for the wide variation in fungal mitochondrial genome sizes (Additional file 1) . A total of 37 mitochondrial group I introns, with an average size of 1,535 bp, were identified in the conserved coding and rRNA genes. These introns are distributed within conserved genes throughout the C. cacaofunesta mitochondrial genome, with cox1 harboring the largest number of these elements (10 introns) (Figure 1, Additional file 1). Cox 1 has been described as a reservoir for mitochondrial group I introns in fungi, harboring as many as 18 of these elements in Agaricus bisporus. Of the total of mitochondrial group I interns, 36 encoded homing endonuclease genes (HEG) (Figure 1, Additional file 1). The dynamic properties of these introns suggest that they can be moved into other regions of the genome and between the genomes of phylogenetically distant species .
Although the presence of mitochondrial plasmids is a common feature of Sordariomycetes, such as G. zeae, N. crassaand P. anserine, no integrated or free linear mitochondrial plasmids were identified in the C. cacaofunesta mtDNA. However, a pseudogene with similarity to the RNA polymerase (rpo) of the P. anserina mitochondrial plasmid pAL2-1(e-value: 4e-92, 100% of query coverage and 57% of positives matches) was identified downstream of the cox1 gene (Figure 1, labeled, “pseudo rpo”). The presence of this pseudogene with plasmidial origin suggests an ancestral plasmid insertion in this region. We found no evidence of inverted repeats flanking this pseudogene.
The analysis of the C. cacaofunesta mtDNA also revealed an average of 26.8% G+C residues, which were uniformly distributed throughout the sequence (Figure 1, blue histogram). The intronic ORFs exhibited a slightly lower GC content (25.2%) than did the conserved genes (26.5%). Moreover, a GC skew analysis was performed, and the shift points of the GC skew graphs were consistent with the loci of the origin (org) and termination (ter) of replication in bacteria  and certain fungi .
Predicted mitochondrial proteome and annotation
Considering that approximately 99% of the mitochondrial proteins are encoded in the nucleus  and that we have the complete genome sequences of C. cacaofunesta, we estimated the total mitochondrial proteome of C. cacaofunesta based on the numbers of mitochondrial- and nuclear-encoded proteins.
Summary of MtProt annotation
A global analysis of the C. cacaofunesta transcriptome was performed using large-scale mRNA sequencing (RNA-Seq), as described in the corresponding section of Methods. This technique is widely used to analyze gene expression and to validate gene predictions [65–68]. This methodology was recently used to analyze mitochondrial gene expression [69, 70]. Approximately 55 million of reads were generated in two biological replicates (CER1 and CER2). As expected, only a very small number of reads mapped to mitochondrial genes (400 from CER1 and 600 from CER2) (Additional file 2A). The RNAseq methodology used here requires a mRNA polyA tail, but RNAs that are transcribed from fungal mtDNAs generally lack this feature. Polyadenylation of mitochondrial mRNA appears to be restricted to higher eukaryotes . Therefore, the 1,000 RNAseq reads aligned with C. cacaofunesta mtDNA most likely represents artifacts.
Experimental mitochondrial proteome
The isolated mitochondrial proteome was electrophoretically separated on a 12.5% SDS-PAGE gel, and 19 bands were excised from the gel and subjected to LC-MS/MS analysis in the corresponding section (Additional file 3). A total of 12,181 MS/MS spectra were generated and confidently assigned to 476 unique proteins with statistically significant (p < 0.05) Mascot-based scores and a 9.1% FDR (Additional file 4, Sheet 4A). Among the 476 proteins that were identified using LC-MS/MS, the 5 most abundant proteins presented more than 140 spectral counts in both samples (Additional file 4, Sheet 4A). The size of these proteins ranged from 34 to 138 kDa, and there was not a molecular weight bias in their abundance.
Among the total proteins that were identified using LC-MS/MS, 64% (304) were contained within the predicted mtProt (Figure 2). This result indicates that approximately 27% of the predicted mtProt was experimentally validated. Accordingly, the mitochondrial proteome analysis by LC-MS/MS is capable of identifying 23%-40% of the known mitochondrial proteins but does not capture low-abundance proteins , except for proteins that are produced under specific circumstances . Experimental mitochondrial proteome analysis of N. crassa has identified 473 proteins that correspond to 169 unique genes; these represent 20% of the total mtProt that has been described for this species .
Proteins identified by LC-MS/MS involved in the respiration metabolic process
Pyruvate dehydrogenase complex
Pyruvate Dehydrogenase (lipoamide) alpha
Pyruvate Dehydrogenase (lipoamide) beta
Pyruvate Dehydrogenase Kinase
Tricarboxylic acid cycle
Isocitrate Dehydrogenase 2
Remarkably, we identified 83 conserved hypothetical proteins and 1 novel protein in the experimental proteome that exhibited putative mitochondrial localization (Additional file 4, Sheet 4B). We believe that these results are an important contribution because they give the first experimental evidence of the existence of these putative mitochondrial proteins. The functional characterization of these proteins and their association with particular mitochondrial pathways is a great challenge but could certainly improve our understanding of this organelle.
We also identified 172 proteins that were not predicted to be mitochondrial based on our analysis (Figure 2, Additional file 4, Sheet 4C). Manual annotation was used to elucidate the composition of this group of proteins. Total number of mitochondrial proteins identified using LC-MS/MS could reach 388 (304 from prediction, and 84 from manual annotation), representing 81% of the total experimental proteome. This result indicates that the enriched-mitochondrial preparation of C. cacaofunesta was performed with high efficiency. However, the standard methodology chosen here was automatic annotation due to the high stringency of the parameters; hence, the 84 proteins that were classified as putative mitochondrial based on manual annotation were not included in the final catalogue of mitochondrial proteins that are proposed here. The experimental proteome was used to validate a fraction of predicted proteome, focusing on the identification of hypothetical and novel proteins.
One of the main contributions of this study is the experimental identification of two groups of predicted mitochondrial proteins with no assigned function: (i) proteins for which there was no experimental evidence (hypothetical proteins) and (ii) novel proteins. Table 1 shows the total of proteins of groups (i) and (ii) based both on the predicted and the experimental proteomes. Of the 349 predicted hypothetical proteins, 83 were identified by mass spectrometry, as mentioned above. With respect to group (ii), the predicted mtProt identified 38 novel putative proteins. However, RNAseq confirmed the expression of 33 of them, and one protein was identified experimentally. The 84 hypothetical proteins that were identified by mass spectrometry were reclassified as conserved unknown function (Additional file 4, Sheet 4B). The S. cerevisiae mitochondrion is the best-understood and characterized molecularly [24, 29, 33, 72]. However, approximately 19% of the identified proteins remain have no known function . The predicted C. cacaofunesta mtProt contains approximately 34% of proteins with unknown function. These data suggest that we are far from fully understanding the function of this organelle, a fact that reflects its functional plasticity. Further functional characterization of these proteins may enable a better understanding of mitochondrial function. It is important to highlight that there are more than 65 known Sordariomycetes genomes published, and we suggest that the 38 predicted novel genes (33 expressed) may be specific to this species. Moreover, it is very likely that these proteins are shared with other species from the Ceratocystis fimbriata complex, a hypothesis that will investigated with the upcoming availability of other Ceratocystis genomes. Our results provide a framework for examining the involvement of novel proteins in mitochondrial pathways.
The functional annotation of the C. cacaofunesta mitochondrial proteome and a comparison with those of S. cerevisiae and N. crassa
We used the same GO categorization criteria as above to annotate the S. cerevisiae and N. crassa mitochondrial proteomes [72, 73] to compare their global functional profiles to that from C. cacaofunesta (Figure 4 and Additional file 6, Sheet 6B and 6C). Due to the redundancy of GO categories associated with a same protein, it is difficult to perform a statistical analysis of these data. However, Figure 4 clearly shows that there is a similar pattern to the distribution of mitochondrial proteins in different GO functional categories among the three species. Considering the important role of mitochondrial function to basic cell metabolism, the conservation of core mitochondrial functions is expected . However, it is known that the expansion or reduction in the size of individual protein families, and hence individual functions, plays a specific role in lifecycle of the organism. In phytopathogenic fungi, differences in protein families have been associated with the importance of the related functions in plant-pathogen interactions . Considering that C. cacaofunesta is a phytopathogenic fungus and that S. cerevisiae and N. crassa are non-pathogenic, it would be plausible to expect more differences between these species regarding the distribution of proteins in functional categories that are related to transport of metabolites, signaling, and defense . However, no such major differences were identified. According the most recent publication on the subject, certain categories that were represented in this analysis, such as signaling, energy metabolism, lipid metabolism, protein transport and protein folding, have similar protein percentages as those that have been described for the S. cerevisiae mtProt . These results suggest that the involvement of the mitochondria in virulence and pathogenicity may rely on specific proteins and/or in the convergence in time and space of different components.
Notably, the global analysis of the C. cacaofunesta mitochondrial function that was based on functional classification (GO) of the experimentally identified proteins (Figure 5) reflects the functional grouping of the full set of predicted mitochondrial proteins (Figure 4). Considering that we performed mass spectrometry only for fungi in cultured conditions, which appears to be far from mimicking fungal growth in nature, the similarity of functions between the predicted and experimental proteomes reinforces the importance of the core mitochondrial functions. Additionally, this similarity suggests that the mitochondrial role in the adaptation of the fungus to certain growth conditions, including the environment of its host, depends of specific sets of proteins, pathways and/or differential gene regulation rather than on major changes in mitochondrial function.
C. cacaofunesta global proteome survey
It is well established that mitochondrial function is important for the virulence and pathogenicity of fungi [16, 77]. Here, we demonstrated that the global proteome compositions of S. cerevisiae, N. crassa and C. cacaofunesta do not contain major differences. However, of these three fungi, C. cacaofunesta is the only pathogenic fungus, which prompted us to perform a more detailed investigation of its proteome composition.
Early stages of necrotrophic infections, like those caused by C. cacaofunesta, are associated with host cell death, the production of secondary metabolites and the accumulation of reactive oxygen species (ROS) [78, 79]. Mitochondrial proteins that are associated with detoxification may play an important role in the success of pathogen colonization. We therefore analyzed the predicted mtProt of C. cacaofunesta to identify proteins that may play important roles in pathogenicity and sensitivity/resistance to fungicides.
The C. cacaofunesta predicted mtProt contains 5 proteins classified as ABC transporters. All of the genes codifying these proteins were expressed under the tested conditions, with average RPKM values of 140.6 and 1,463.7. Mass spectrometry identified 2 of this class of transporters (Additional file 2, Sheet 2B; Additional file 4, Sheet 4A). The ABC transporter of a Nectria haematococca (NhABC1) proved to be an important virulence factor that contributes to the N. haematococca tolerance of phytoalexin during its interaction with the host . Furthermore, ABC transporters have been associated with fungicide resistance .
Aldehyde dehydrogenase (ADH) has also been implicated in the pathogenicity of fungi. ADH is regarded as a detoxification enzyme due to its role in the metabolism of intermediates and exogenous aldehydes . An analysis of adh expression in Clasdosporium fulvum, both in vitro and in vivo, during its interactions with tomato indicated that the expression of this gene was associated with infection and starvation conditions . Interestingly, the predicted C. cacaofunesta mtProt survey identified 12 ADH-encoding genes with confirmed expression (RPKMs > 10), and 6 ADH proteins were identified by LC-MS/MS. Moreover, glutathione S-transferases (GSTs) are a large family of proteins that are involved in cellular detoxification processes during oxidative stress responses  and may be involved in fungicide resistance . Mao and coworkers (2012) suggested a link between the number of GST-encoding genes in the fungal genome and its lifecycle . The predicted C. cacaofunesta mtProt contains 4 proteins that are homologous to GST. The RNA-Seq analysis indicated varying average expression levels, with RPKM values of as high as 150; one of those proteins was identified in the experimental proteome (Additional file 2, Sheet 2B; Additional file 4, Sheet 4B).
Moreover, components of the alternative respiratory chain have been related to fungal pathogenicity, specifically alternative NADH dehydrogenase  and alternative oxidase (AOX) . Studies in which alternative NADH dehydrogenase was blocked in Aspergillus niger demonstrated the role of this protein in protection against oxidative stress . C. cacaofunesta has one encoded alternative NADH dehydrogenase gene with RPKM values > 60; this gene was also identified in the experimental proteome (Additional file 2, Sheet 2B; Additional file 4, Sheet 4B). Moreover, the predicted mtProt included a putative AOX. This enzyme has been implicated in fungal protection against oxidative stress, pathogenicity and fungicide resistance [89, 90]. C. cacaofunesta AOX was not identified in the experimental proteome, but the gene aox exhibited detectable expression, with a RPKM value of approximately 88 (Additional file 2, Sheet 2B; Additional file 4, Sheet 4A). Recently, Thomazella and coworkers (2012) provided strong evidence for the involvement of AOX in maintaining the infective phase of the basidiomycete Moniliophthora perniciosa.
The deletion of a Fusarium oxysporum mitochondrial carrier-encoding gene (fow 1) caused the loss of its ability to infect plant tissues . A putative FOW1 protein was identified in the C. cacaofunesta predicted mtProt. The RNA-Seq data indicated that the FOW1 homologue in C. cacaofunesta was highly expressed (with an RPKM value of 536), and its product was identified in the experimental proteome (Additional file 2, Sheet 2B; Additional file 4, Sheet 4A).
It has been reported that the ability of certain fungi to withstand different levels of hypoxia is a critical virulence factor in terms of fungal pathogenicity in mammals . Hypoxia-related genes have been identified as important virulence factors for Aspergillus fumigatus and Candida albicans given that they are exposed to low oxygen concentrations during infection [93, 94]. Oxygen levels in the secondary xylem, the plant tissue that is infected by C. cacaofunesta, are as low as are observed in human tissues that are infected by these fungal pathogens, reaching levels less than 5% [93, 95]. This observation suggests C. cacaofunesta should possess mechanism with which to adapt to hypoxia. In this regard, the predicted C. cacaofunesta mtProt contains one protein that has a conserved hypoxia response domain. This gene was expressed in in vitro conditions, with an RPKM value as high as 259 and its product was also identified in the experimental proteome (Additional file 2, Sheet 2B; Additional file 4, Sheet 4A). Although hypoxia response mechanisms are more complex and involve a series of regulatory pathways , these results open the way for the exploration of the role of C. cacaofunesta mitochondria in fungal phytopathogens responses to hypoxia.
Due to inherent limitations for each of the techniques that were used to predict the subcellular localization of proteins, we chose to perform a stringent integrative approach to predict the C. cacaofunesta mtProt. We therefore must state that the list of mitochondrial proteins proposed here as the predicted mitochondrial proteome of C. cacaofunesta may not correspond to the full set of mitochondrial polypeptides. This mtProt survey is a preliminary analysis due to the relatively small number of proteins discussed but reveals how these data can be exploited in future studies. The presented genome and estimated proteome data provide a starting point for the improved characterization of C. cacaofunesta mitochondrial function in two ways: (i) facilitating the identification of proteins that may be relevant to the lifestyle of this fungus and (ii) expanding the literature that relates to studies of mitochondrial function and evolution to the discovery of new mitochondrial genes. In this regard, the identification of both conserved proteins with unknown function and novel proteins is particularly relevant. The comprehensive characterization of these proteins opens up exciting new opportunities for understanding the role of mitochondria in phytopathogenicity. Due to the important role of this organelle in resistance/susceptibility to fungicides, this study could aid in the identification of more effective strategies by which to the control this fungus.
In the present study, we performed a global analysis of C. cacaofunesta mitochondria using an integrative approach (mitogenomics, transcriptomics and proteomics). We predicted the C. cacaofunesta mitochondrial proteome, including 52 mitochondrial encoded proteins and 1,072 proteins encoded in the nuclear genome (NMP). All of the NMP (except for 10 genes) were transcribed in mycelia that were grown in vitro; approximately 27% of these genes were also identified by mass spectrometry. RNAseq analysis allowed to detected expression of 33 novel putative genes. The comparison of the global analysis of mitochondrial function between C. cacaofunesta, S. cerevisiae and N. crassa revealed that, in addition to the known involvement of mitochondria in the pathogenesis of fungi, the differences in the lifestyle of these organisms (pathogenic and no-pathogenic) are not accompanied by major differences in the functional composition of their mitochondrial proteomes. Therefore, the relevance of this organelle in different fungi lifestyle could lie in a particular set of proteins and/or in a different pattern of gene regulation. In this regard, LC-MS analysis validated the existence of 84 proteins with unknown function, of which one is probably specific to this species and likely to the Ceratocystis species complex. Further functional characterization of these putative mitochondrial proteins could improve our understanding of the mitochondria’s role in fungal pathogenesis.
Biological material and nucleic acid isolation
Dr. Tomas Harrington, from the Iowa State University Department of Plant Pathology, generously donated Ceratocystis cacaofunesta strain C1593. This strain was isolated in 1999 from infected cacao trees located on a farm in the district of Uruçuca, Bahia, Brazil. In our laboratory, C. cacaofunesta cultures were maintained on 2% malt, yeast extract and agar (MYEA) plates inside a BOD chamber at 28°C. For liquid cultures, 20 to 30 agar blocks (1 mm2) that were removed from the solid cultures were inoculated into flasks that contained 100 mL of malt media and were cultivated at 28°C for 7 days, under constant agitation at 150 rpm.
The mycelia were separated from the media using filtration, were washed twice with sterile distilled water and were frozen with liquid nitrogen. The samples were ground to a fine powder in a mortar and were processed for DNA isolation as previously described . The RNA isolation was performed using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s protocol. All of the isolated genetic materials were qualitatively analyzed using denaturing formaldehyde/agarose gel electrophoresis and were quantified using a NanoDropTM 1000 spectrophotometer (Thermo Scientific).
Genome sequencing and assembly
The C. cacaofunesta mitochondrial genome was sequenced as part of the C. cacaofunesta Genome Project (http://www.lge.ibi.unicamp.br/ceratocystis). The DNA was sequenced on a Genome Analyzer IIx platform (Illumina) at the University of North Carolina, Chapel Hill High-Throughput Sequencing Facility. The whole-genome shotgun strategy was used to produce 76-bp paired-end reads (400-bp insert size) and 36-bp mate-pair reads (3-kb insert size).
The paired-end reads were assembled into longer scaffolds using de novo assembler VELVET 1.0.12 software  with a k-mer parameter of 69, which maximizes the length of the mitochondrial DNA (mtDNA) contig. The mtDNA was identified through comparisons with C. cacaofunesta scaffolds and Gibberella zeae mtDNA (NC_009493), resulting in the identification of a single mtDNA contig. To verify the topology of the mtDNA, we used the paired-end and mate-pair reads that had been aligned with the mtDNA using the SOAP2 aligner . The reads that mapped to the borders of the contig were located at an expected distance from their respective pairs.
The complete mtDNA sequence was deposited into the GenBank database and is available under the accession number JX185564.
C. cacaofunesta mtDNA was inspected by tBLASTn searches to identify 14 known conserved coding genes using their Gibberella zeae mtDNA orthologous as the query (Accession NC_009493). The rps3 gene was predicted using the C. fimbriata rps3 sequence available at GenBank (Accession FJ895616.1). Open Reading Frame Finder (ORF Finder) was used to search for intronic and hypothetical ORFs using Genetic code 4 (mold mitochondrial). Alternative start codons were allowed for predicted intronic ORFs when the initial ATG site resulted in an incomplete LAGLIDADG or GIY-YIG domain as indicated by a BLASTp search. The set of tRNAs was identified using tRNAscan-SE . Genes encoding both small and large rRNA subunits were identified through comparisons with orthologous Gibberella zeae sequences using BLAST2seq (BLASTn).
mtDNA GC content and GC skew
The GC content, local GC content and cumulative GC skew of the C. cacaofunesta mtDNA were calculated using customized Perl scripts (available upon request). For the local GC content and cumulative GC skew, a sliding window of over 5,000 bp with a 500-bp range was used.
Predicted mitochondrial proteome
The putative genes encoding NMPs were identified using two strategies: ab initio prediction and sequence comparison. For the ab initio prediction, SignalP 3.0 software  was used to analyze all of the predicted nuclear proteins to identify those with a low probability (≤ 50%) of containing a signal peptide. Using the SignalP results as the input, the WoLF PSORT  and TargetP  programs were used to classify putative mitochondrial proteins. Within the TargetP output, only the proteins that were identified with a high probability of encoding mitochondrial genes (mTP ≥ 50% and other ≤ 50%) were considered. The proteins that were identified by both programs (TargetP and WoLF PSORT) constituted the final ab initio protein dataset. For the sequence comparison approach, we built a database that contained 1,583 known fungal mitochondrial proteins (Additional file 7; 741 proteins from S. cerevisiae and 842 from N. crassa). After the C. cacaofunesta nuclear encoded proteins were compared against this database using BLASTp (with an e-value cutoff of 1e-5), the comparative protein dataset was obtained from the BLAST results by applying query coverage (≥ 70%) and similarity (≥ 70%) filters. The ab initio and comparative protein datasets were merged to generate the final NMP dataset. The sequences of the nuclear-encoded mitochondrial proteins from C. cacaofunesta are provided in Additional file 8.
The automatic annotation of the NMPs was performed using the NR/NCBI, KEGG  and UniRef90 databases (BLASTp with an e-value cutoff of 1e-5)  and summarized using the AutoFACT program . The CDD/Pfam database was used to identify the conserved domains. The Blast2GO program  was used to perform the gene ontology classification (BLASTp, with an e-value cutoff of 1e-5 on Generic GO Slim; the first 50 hits were considered). The proteins with identified biological functions according to GO Slim were manually grouped into 12 classes: amino acid metabolism, carbohydrate metabolism, cell organization, energy metabolism, genome maintenance and transcription, lipid metabolism, protein metabolism, protein transport and folding, response to stimulus, signaling, transport of metabolites and other. The proteins that were classified into multiple categories were considered to be members of both categories, with the exception of proteins that were classified as both “other” and another category; these were excluded from the “other” category and remained in the specific category. Thus, the “other” category includes proteins of known function that were not included in any of the above-described categories. The GO Slim biological processes that clustered in each category are shown in Additional file 5. Similar analyses were performed for the S. cerevisiae and N. crassa mitochondrial proteomes to compare them with the C. cacaofunesta data (Additional file 6, Sheets 6A, 6B and 6C).
Preparation of RNA-seq libraries and transcriptome analysis
RNA from C. cacaofunesta liquid cultures was extracted using the RNeasy Plant mini kit (Qiagen). Approximately 5 μg total RNA were used to prepare RNA-seq libraries following the procedures described by the manufacturer (Illumina). Briefly, the mRNA was purified using Sera-mag oligo (dT) beads and subsequently used for cDNA synthesis. Double-stranded cDNA was subjected to end-repair, A-tailing, adapter ligation and PCR amplification. The libraries were quantified using a Qubit fluorometer (Invitrogen), and quality control was performed using the Experion automated electrophoresis system (Bio-Rad). Two independent biological replicates were used for transcriptome sequencing. Each sample was sequenced in one lane of an Illumina Genome Analyzer IIx sequencer. Approximately 27 and 29 million single-end reads of 36 bp were produced for the CER1 and CER2 libraries, respectively.
The SOAP2 aligner  was used to align the RNA-Seq reads from mycelia that were grown in vitro with predicted genes that encode nuclear-encoded mitochondrial proteins. The program was configured to allow for as many as two mismatches, to discard sequences that contained “N”s and to return all optimal alignments. The expression level of each gene was estimated as an RPKM (reads per kilobase of exon per million reads mapped) value .
Mitochondrial isolation and proteomic assays
The mitochondrial isolation was performed according to the protocol that is described by Sorensen and coworkers , with minor modifications. C. cacaofunesta was cultured for seven days in standard MYEA growth medium. Subsequently, the fungal mycelia were washed in cold water and homogenized using a bead beater in cold extraction buffer that contained 330 mM sucrose, 10 mM Tris–HCl, 1 mM EDTA, 1% PVP, 0.1% BSA and 0.3 mM PMSF (pH 7.5). The homogenates were filtered and centrifuged at 1,500 g for 15 min to remove the cellular debris. The resulting supernatants were centrifuged at 15,500 g for 20 min. The samples were rinsed twice in wash buffer (330 mM sucrose, 10 mM Tris–HCl, 1 mM EDTA, pH 7.2), and the final mitochondrial preparation was used to perform the proteomic assays.
The identification of proteins using LC-MS/MS
Enzymatic in-gel digestion for mass spectrometry analysis
A total of 20 μg of the concentrated mitochondrial proteome was obtained as described above and was separated using one-dimensional SDS-PAGE electrophoresis. Nineteen bands were excised from the gel and subjected to in-gel trypsin digestion, as described previously , with modifications.
The resulting peptide solution was dried in a SpeedVac concentrator and resuspended in 100 μL of 0.1% formic acid. An aliquot of 4.5 μL was separated using C18 (75 μm × 100 mm) RP-nanoUPLC (nanoACQUITY, Waters) coupled with a Q-Tof Ultima mass spectrometer (Waters) with a nano-electrospray source. The flow rate was 600 nL/min, and the gradient was 2-90% acetonitrile in 0.1% formic acid over 45 min. The instrument was operated in the “top three” mode, in which one MS spectrum was acquired, followed by an MS/MS analysis of the three most intense peaks .
The spectra were acquired using MassLynx v.4.1 software, and the raw data files were converted to a peak list format (mgf) using the Mascot Distiller v.220.127.116.11, 2009 software (Matrix Science Ltd.) without summing the scans, allowing for a label-free analysis (Additional file 9). The files were then searched against the Ceratocystis cacaofunesta database (7,321 entries – 7,269 nuclear proteins, and 52 mitochondrial proteins) using Mascot engine v.2.3.01 (Matrix Science Ltd.). Carbamidomethylation was used as a fixed modification, methionine oxidation as a variable modification, one missed trypsin cleavage and a tolerance of 0.1 Da for both precursor and fragment ions. For the protein quantitation, the .dat files from the Mascot output were analyzed using Scaffold Q+ (version 3_00_03, Proteome Software), and quantitative values (normalized spectral counts) were obtained [106, 107]. For the endogenous peptide identification, methionine oxidation was set as a variable modification, with a tolerance of 0.1 Da for both the precursor and fragment ions. Only peptides with a minimum of 5 amino acid residues and significance (p < 0.05) based on the Mascot-based scores were considered in the results.
- ABC tranporter:
ATP binding cassette (family of membrane transport proteins)
Basic local alignment search tool
Bovine serum albumin
Conserved domains database
Cytochrome b oxidade
Cytochrome c oxidase
Fusarium oxysporum mitochondrial carrier
Isocitrate dehydrogenases (IDH1 and IDH2)
Kyoto encyclopedia of genes and genomes
Nicotinamide adenine dinucleotide (electron donor)
Reduced B-nicotinamide adenine dinucleotide
National Center for Biotechnology Information
Nuclear-encoded mitochondrial proteins
Non-redundant GenBank database
Open reading frame
Proteins of oxidative phosphorylation machinery
Proteins of the pyruvate dehydrogenase complex (
Large RNA subunit
Small RNA subunit
Reads per kilobase of exons per million
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Tricarboxylic acid cycle
Tris (hydroxymethyl) aminomethane hydrochloride
Unknown function protein
The authors would like to acknowledge to Dr. Tomas Harrington, from the Iowa State University Department of Plant Pathology, who generously donated Ceratocystis cacaofunesta strain C1593. This work was supported by research funding of Fundação de Amparo à Pesquisa do Estado de São Paulo, grant number 2009/50119-9.
- Baker CJ, Harrington TC, Krauss U, Alfenas AC: Genetic variability and host specialization in the latin american clade of ceratocystis fimbriata. Phytopathology. 2003, 93 (10): 1274-1284. 10.1094/PHYTO.2003.93.10.1274.PubMedGoogle Scholar
- Johnson JA, Harrington TC, Engelbrecht CJ: Phylogeny and taxonomy of the North American clade of the ceratocystis fimbriata complex. Mycologia. 2005, 97 (5): 1067-1092. 10.3852/mycologia.97.5.1067.PubMedGoogle Scholar
- Barnes I, Gaur A, Burgess T, Roux J, Wingfield BD, Wingfield MJ: Microsatellite markers reflect intra-specific relationships between isolates of the vascular wilt pathogen Ceratocystis fimbriata. Mol Plant Pathol. 2001, 2 (6): 319-325. 10.1046/j.1464-6722.2001.00080.x.PubMedGoogle Scholar
- Engelbrecht CJ, Harrington TC, Steimel J, Capretti P: Genetic variation in eastern North American and putatively introduced populations of Ceratocystis fimbriata f. platani. Mol Ecol. 2004, 13 (10): 2995-3005. 10.1111/j.1365-294X.2004.02312.x.PubMedGoogle Scholar
- Soanes DM, Alam I, Cornell M, Wong HM, Hedeler C, Paton NW, Rattray M, Hubbard SJ, Oliver SG, Talbot NJ: Comparative genome analysis of filamentous fungi reveals gene family expansions associated with fungal pathogenesis. PLoS One. 2008, 3 (6): e2300-10.1371/journal.pone.0002300.PubMed CentralPubMedGoogle Scholar
- Engelbrecht CJ, Harrington TC: Intersterility, morphology and taxonomy of Ceratocystis fimbriata on sweet potato, cacao and sycamore. Mycologia. 2005, 97 (1): 57-69. 10.3852/mycologia.97.1.57.PubMedGoogle Scholar
- Bastos CNE HC: Ocorrencia de Ceratocystis fimbriata. Acta Amazonica. 1978, 8 (4): 543-544.Google Scholar
- Bezerra JL: Ceratocystis fimbriata causing death of budded cocoa seedlings in Bahia, Brazil. Incoped Newsletter. 1997, 1:Google Scholar
- Malaguti G: Ceratostomella fimbriata en el cacao de Venezuela. Acta cientif venezol. 1952, 3 (3): 94-97.Google Scholar
- Harrington TC: Host specialization and speciation in the American wilt pathogen Ceratocystis fimbriata. Fitopatol Bras. 2000, 25S: 262-263.Google Scholar
- Spence JL: Preliminary observations on a wilt condition of cocoa. J Agric Soc Trinidad and Tobago. 1958, 58 (3): 349-359.Google Scholar
- Thorold CA: Diseseas of cocoa. 1975, Oxford: Clarendon Press, 423-XIIGoogle Scholar
- Engelbrecht CJ, Harrington TC, Alfenas A: Ceratocystis wilt of cacao-a disease of increasing importance. Phytopathology. 2007, 97 (12): 1648-1649. 10.1094/PHYTO-97-12-1648.PubMedGoogle Scholar
- Ferreira MA, Harrington TC, Alfenas AC, Mizubuti ES: Movement of genotypes of Ceratocystis fimbriata within and among Eucalyptus plantations in Brazil. Phytopathology. 2011, 101 (8): 1005-1012. 10.1094/PHYTO-01-11-0015.PubMedGoogle Scholar
- Harrington TC, Thorpe DJ, Alfenas AC: Genetic variation and variation in aggressiveness to native and exotic hosts among Brazilian populations of Ceratocystis fimbriata. Phytopathology. 2011, 101 (5): 555-566. 10.1094/PHYTO-08-10-0228.PubMedGoogle Scholar
- Shingu-Vazquez M, Traven A: Mitochondria and fungal pathogenesis: drug tolerance, virulence, and potential for antifungal therapy. Eukaryot Cell. 2011, 10 (11): 1376-1383. 10.1128/EC.05184-11.PubMed CentralPubMedGoogle Scholar
- Schatz G: Mitochondria: beyond oxidative phosphorylation. Biochim Biophys Acta. 1995, 1271 (1): 123-126. 10.1016/0925-4439(95)00018-Y.PubMedGoogle Scholar
- Gunter TE, Yule DI, Gunter KK, Eliseev RA, Salter JD: Calcium and mitochondria. FEBS Lett. 2004, 567 (1): 96-102. 10.1016/j.febslet.2004.03.071.PubMedGoogle Scholar
- Muhlenhoff U, Richhardt N, Gerber J, Lill R: Characterization of iron-sulfur protein assembly in isolated mitochondria. A requirement for ATP, NADH, and reduced iron. J Biol Chem. 2002, 277 (33): 29810-29816. 10.1074/jbc.M204675200.PubMedGoogle Scholar
- Eisenberg T, Buttner S, Kroemer G, Madeo F: The mitochondrial pathway in yeast apoptosis. Apoptosis. 2007, 12: 1011-1023. 10.1007/s10495-007-0758-0.PubMedGoogle Scholar
- Lorin S, Dufour E, Sainsard-Chanet A: Mitochondrial metabolism and aging in the filamentous fungus Podospora anserina. Biochimica et Biophysica Acta (BBA) - Bioenergetics. 1757, 5–6: 604-610.Google Scholar
- Court DA, Griffiths AJF, Kraus SR, Russell PJ, Bertrand H: A new senescence-inducing mitochondrial linear plasmid in field-isolated <i>Neurospora crassa strains from India. Curr Genet. 1991, 19 (2): 129-137. 10.1007/BF00326294.PubMedGoogle Scholar
- Torriani SF, Goodwin SB, Kema GH, Pangilinan JL, McDonald BA: Intraspecific comparison and annotation of two complete mitochondrial genome sequences from the plant pathogenic fungus Mycosphaerella graminicola. FG & B. 2008, 45 (5): 628-637.Google Scholar
- Solieri L: Mitochondrial inheritance in budding yeasts: towards an integrated understanding. Trends Microbiol. 2010, 18 (11): 521-530. 10.1016/j.tim.2010.08.001.PubMedGoogle Scholar
- Wiedemann N, Frazier AE, Pfanner N: The protein import machinery of mitochondria. J Biol Chem. 2004, 279 (15): 14473-14476. 10.1074/jbc.R400003200.PubMedGoogle Scholar
- Gray MW, Lang BF: Transcription in chloroplasts and mitochondria: a tale of two polymerases. Trends Microbiol. 1998, 6 (1): 1-3. 10.1016/S0966-842X(97)01182-7.PubMedGoogle Scholar
- Dolezal P, Likic V, Tachezy J, Lithgow T: Evolution of the molecular machines for protein import into mitochondria. Science. 2006, 313 (5785): 314-318. 10.1126/science.1127895.PubMedGoogle Scholar
- Neupert W, Herrmann JM: Translocation of proteins into mitochondria. Annu Rev Biochem. 2007, 76: 723-749. 10.1146/annurev.biochem.76.052705.163409.PubMedGoogle Scholar
- Schmidt O, Pfanner N, Meisinger C: Mitochondrial protein import: from proteomics to functional mechanisms. Nat Rev Mol Cell Biol. 2010, 11 (9): 655-667. 10.1038/nrm2959.PubMedGoogle Scholar
- Costa GG, Cabrera OG, Tiburcio RA, Medrano FJ, Carazzolle MF, Thomazella DP, Schuster SC, Carlson JE, Guiltinan MJ, Bailey BA: The mitochondrial genome of Moniliophthora roreri, the frosty pod rot pathogen of cacao. Fungal Biol. 2012, 116 (5): 551-562. 10.1016/j.funbio.2012.01.008.PubMedGoogle Scholar
- David Day A, Heazlewood J, Julian Tonti-Filippini S, Alexander Gout M, David Day A, James W, Harvey Millar A: Experimental analysis of the arabidopsis mitochondrial proteome highlights signaling and regulatory components, provides assessment of targeting prediction programs, and indicates plant-specific mitochondrial proteins. Plant Cell. 2004, 16 (1): 241-256. 10.1105/tpc.016055.PubMed CentralPubMedGoogle Scholar
- Mercer Tim R, Neph S, Dinger Marcel E, Crawford J, Smith Martin A, Shearwood A-Marie J, Haugen E, Bracken Cameron P, Rackham O, Stamatoyannopoulos John A: The human mitochondrial transcriptome. Cell. 2011, 146 (4): 645-658. 10.1016/j.cell.2011.06.051.PubMed CentralPubMedGoogle Scholar
- Prokisch H, Scharfe C, Camp DG, Xiao W, David L, Andreoli C, Monroe ME, Moore RJ, Gritsenko MA, Kozany C: Integrative analysis of the mitochondrial proteome in yeast. PLoS Biol. 2004, 2 (6): e160-10.1371/journal.pbio.0020160.PubMed CentralPubMedGoogle Scholar
- Reinders J, Zahedi RP, Pfanner N, Meisinger C, Sickmann A: Toward the complete yeast mitochondrial proteome: multidimensional separation techniques for mitochondrial proteomics. J Proteome Res. 2006, 5 (7): 1543-1554. 10.1021/pr050477f.PubMedGoogle Scholar
- Doyle C, Donaldson M, Morrison E, Saville B: Ustilago maydis transcript features identified through full-length cDNA analysis. Mol Genet Genomics. 2011, 286 (2): 143-159. 10.1007/s00438-011-0634-z.PubMedGoogle Scholar
- Tan K-C, Ipcho SVS, Trengove RD, Oliver RP, Solomon PS: Assessing the impact of transcriptomics, proteomics and metabolomics on fungal phytopathology. Mol Plant Pathol. 2009, 10 (5): 703-715. 10.1111/j.1364-3703.2009.00565.x.PubMedGoogle Scholar
- Gonzalez-Fernandez R, Prats E, Jorren-Novo JV: Proteomics of plant pathogenic fungi. J Biomed Biotechnol. 2010, 2010: 1-36.Google Scholar
- Yajima W, Kav NNV: The proteome of the phytopathogenic fungus Sclerotinia sclerotiorum. Proteomics. 2006, 6 (22): 5995-6007. 10.1002/pmic.200600424.PubMedGoogle Scholar
- Fernández-Acero FJ, Jorge I, Calvo E, Vallejo I, Carbú M, Camafeita E, López JA, Cantoral JM, Jorrín J: Two-dimensional electrophoresis protein profile of the phytopathogenic fungus Botrytis cinerea. Proteomics. 2006, 6 (S1): S88-S96. 10.1002/pmic.200500436.PubMedGoogle Scholar
- Medina M: Genomes, phylogeny, and evolutionary systems biology. Proc Natl Acad Sci U S A. 2005, 102 (Suppl 1): 6630-6635.PubMed CentralPubMedGoogle Scholar
- Godfrey D, Zhang Z, Saalbach G, Thordal-Christensen H: A proteomics study of barley powdery mildew haustoria. Proteomics. 2009, 9 (12): 3222-3232. 10.1002/pmic.200800645.PubMedGoogle Scholar
- Xu J, Saunders CW, Hu P, Grant RA, Boekhout T, Kuramae EE, Kronstad JW, DeAngelis YM, Reeder NL, Johnstone KR: Dandruff-associated Malassezia genomes reveal convergent and divergent virulence traits shared with plant and human fungal pathogens. Proc Natl Acad Sci. 2007, 104 (47): 18730-18735. 10.1073/pnas.0706756104.PubMed CentralPubMedGoogle Scholar
- Amlacher S, Sarges P, Flemming D, van Noort V, Kunze R, Devos Damien P, Arumugam M, Bork P, Hurt E: Insight into structure and assembly of the nuclear pore complex by utilizing the genome of a eukaryotic thermophile. Cell. 2011, 146 (2): 277-289. 10.1016/j.cell.2011.06.039.PubMedGoogle Scholar
- Paquin B, Franz Lang B: The mitochondrial DNA of allomyces macrogynus: the complete genomic sequence from an ancestral fungus. J Mol Biol. 1996, 255 (5): 688-701. 10.1006/jmbi.1996.0056.PubMedGoogle Scholar
- Turmel M, Lemieux C, Burger G, Lang BF, Otis C, Plante I, Gray MW: The complete mitochondrial DNA sequences of nephroselmis olivacea and pedinomonas minor: Two radically different evolutionary patterns within green algae. The Plant Cell Online. 1999, 11 (9): 1717-1729.Google Scholar
- Sethuraman J, Majer A, Friedrich NC, Edgell DR, Hausner G: Genes within genes: multiple LAGLIDADG homing endonucleases target the ribosomal protein S3 gene encoded within an rnl group I intron of ophiostoma and related taxa. Mol Biol Evol. 2009, 26 (10): 2299-2315. 10.1093/molbev/msp145.PubMedGoogle Scholar
- Sethuraman J, Majer A, Iranpour M, Hausner G: Molecular evolution of the mtDNA encoded rps3 gene among filamentous ascomycetes fungi with an emphasis on the ophiostomatoid fungi. J Mol Evol. 2009, 69 (4): 372-385. 10.1007/s00239-009-9291-9.PubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25 (5): 955-964.PubMed CentralPubMedGoogle Scholar
- Pantou MP, Kouvelis VN, Typas MA: The complete mitochondrial genome of the vascular wilt fungus Verticillium dahliae: a novel gene order for Verticillium and a diagnostic tool for species identification. Curr Genet. 2006, 50 (2): 125-136. 10.1007/s00294-006-0079-9.PubMedGoogle Scholar
- Pantou MP, Kouvelis VN, Typas MA: The complete mitochondrial genome of fusarium oxysporum: insights into fungal mitochondrial evolution. Gene. 2008, 419 (1–2): 7-15.PubMedGoogle Scholar
- Ferandon C, Moukha S, Callac P, Benedetto JP, Castroviejo M, Barroso G: The Agaricus bisporus cox1 gene: the longest mitochondrial gene and the largest reservoir of mitochondrial group i introns. PLoS One. 2010, 5 (11): e14048-10.1371/journal.pone.0014048.PubMed CentralPubMedGoogle Scholar
- Chevalier BS, Stoddard BL: Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility. Nucleic Acids Res. 2001, 29 (18): 3757-3774. 10.1093/nar/29.18.3757.PubMed CentralPubMedGoogle Scholar
- Jurgenson JE, Bowden RL, Zeller KA, Leslie JF, Alexander NJ, Plattner RD: A genetic map of Gibberella zeae (Fusarium graminearum). Genetics. 2002, 160 (4): 1451-1460.PubMed CentralPubMedGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S: The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003, 422 (6934): 859-868. 10.1038/nature01554.PubMedGoogle Scholar
- Osiewacz HD, Hermanns J, Marcou D, Triffi M, Esser K: Mitochondrial DNA rearrangements are correlated with a delayed amplification of the mobile intron (plDNA) in a long-lived mutant of Podospora anserina. Mutat Res. 1989, 219 (1): 9-15. 10.1016/0921-8734(89)90036-2.PubMedGoogle Scholar
- Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277 (5331): 1453-1462. 10.1126/science.277.5331.1453.PubMedGoogle Scholar
- Andersson SG, Kurland CG: Ancient and recent horizontal transfer events: the origins of mitochondria. APMIS Suppl. 1998, 84: 5-14.PubMedGoogle Scholar
- Schmitt S, Prokisch H, Schlunck T, Camp DG, Ahting U, Waizenegger T, Scharfe C, Meitinger T, Imhof A, Neupert W: Proteome analysis of mitochondrial outer membrane from Neurospora crassa. Proteomics. 2006, 6 (1): 72-80.PubMedGoogle Scholar
- Pagliarini DJ, Calvo SE, Chang B, Sheth SA, Vafai SB, Ong S-E, Walford GA, Sugiana C, Boneh A, Chen WK:A Mitochondrial Protein Compendium Elucidates Complex I Disease Biology. Cell. 2008, 134 (1): 112-123. 10.1016/j.cell.2008.06.016.PubMed CentralPubMedGoogle Scholar
- Koski LB, Gray MW, Lang BF, Burger G: AutoFACT: an automatic functional annotation and classification tool. BMC Bioinforma. 2005, 6: 151-10.1186/1471-2105-6-151.Google Scholar
- McGinnis S, Madden TL: BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004, 32 (suppl 2): W20-W25.PubMed CentralPubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralPubMedGoogle Scholar
- Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res. 2000, 28 (1): 263-266. 10.1093/nar/28.1.263.PubMed CentralPubMedGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedGoogle Scholar
- Wilhelm BT, Landry J-R: RNA-Seq—quantitative measurement of expression through massively parallel RNA-sequencing. Methods. 2009, 48 (3): 249-257. 10.1016/j.ymeth.2009.03.016.PubMedGoogle Scholar
- Au KF, Jiang H, Lin L, Xing Y, Wong WH: Detection of splice junctions from paired-end RNA-seq data by SpliceMap. Nucleic Acids Res. 2010, 38 (14): 4570-4578. 10.1093/nar/gkq211.PubMed CentralPubMedGoogle Scholar
- Richard H, Schulz MH, Sultan M, Nürnberger A, Schrinner S, Balzereit D, Dagand E, Rasche A, Lehrach H, Vingron M: Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments. Nucleic Acids Res. 2010, 38 (10): e112-10.1093/nar/gkq041.PubMed CentralPubMedGoogle Scholar
- Rogers M, Thomas J, Reddy A, Ben-Hur A: SpliceGrapher: detecting patterns of alternative splicing from RNA-Seq data in the context of gene models and EST data. Genome Biol. 2012, 13 (1): R4-10.1186/gb-2012-13-1-r4.PubMed CentralPubMedGoogle Scholar
- Torres TT, Dolezal M, Schlötterer C, Ottenwälder B: Expression profiling of Drosophila mitochondrial genes via deep mRNA sequencing. Nucleic Acids Res. 2009, 37 (22): 7509-7518. 10.1093/nar/gkp856.PubMed CentralPubMedGoogle Scholar
- Neira-Oviedo M, Tsyganov-Bodounov A, Lycett GJ, Kokoza V, Raikhel AS, Krzywinski J: The RNA-Seq approach to studying the expression of mosquito mitochondrial genes. Insect Mol Biol. 2011, 20 (2): 141-152. 10.1111/j.1365-2583.2010.01053.x.PubMedGoogle Scholar
- Chang JH, Tong L: Mitochondrial poly(A) polymerase and polyadenylation. Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms. 1819, 9–10: 992-997.Google Scholar
- Sickmann A, Reinders J, Wagner Y, Joppich C, Zahedi R, Meyer HE, Schonfisch B, Perschil I, Chacinska A, Guiard B: The proteome of Saccharomyces cerevisiae mitochondria. Proc Natl Acad Sci U S A. 2003, 100 (23): 13207-13212. 10.1073/pnas.2135385100.PubMed CentralPubMedGoogle Scholar
- Keeping A, DeAbreu D, DiBernardo M, Collins RA: Gel-based mass spectrometric and computational approaches to the mitochondrial proteome of Neurospora. Fungal Genet Biol. 2011, 48 (5): 526-536. 10.1016/j.fgb.2010.11.011.PubMedGoogle Scholar
- Scheffler IE: Mitochondria make a come back. Adv Drug Deliv Rev. 2001, 49 (1–2): 3-26.PubMedGoogle Scholar
- Idnurm A, Howlett BJ: Pathogenicity genes of phytopathogenic fungi. Mol Plant Pathol. 2001, 2 (4): 241-255. 10.1046/j.1464-6722.2001.00070.x.PubMedGoogle Scholar
- Vödisch M, Scherlach K, Winkler R, Hertweck C, Braun H-P, Roth M, Haas H, Werner ER, Brakhage AA, Kniemeyer O: Analysis of the aspergillus fumigatus proteome reveals metabolic changes and the activation of the pseurotin a biosynthesis gene cluster in response to hypoxia. J Proteome Res. 2011, 10 (5): 2508-2524. 10.1021/pr1012812.PubMed CentralPubMedGoogle Scholar
- Calvo S, Jain M, Xie X, Sheth SA, Chang B, Goldberger OA, Spinazzola A, Zeviani M, Carr SA, Mootha VK: Systematic identification of human mitochondrial disease genes through integrative genomics. Nat Genet. 2006, 38 (5): 576-582. 10.1038/ng1776.PubMedGoogle Scholar
- Mengiste T: Plant immunity to necrotrophs. Annu Rev Phytopathol. 2012, 50 (1): 267-294. 10.1146/annurev-phyto-081211-172955.PubMedGoogle Scholar
- Coleman JJ, White GJ, Rodriguez-Carres M, VanEtten HD: An ABC transporter and a cytochrome P450 of nectria haematococca MPVI Are virulence factors on Pea and Are the major tolerance mechanisms to the phytoalexin pisatin. Mol Plant Microbe Interact. 2010, 24 (3): 368-376.Google Scholar
- de Waard MA: Significance of ABC transporters in fungicide sensitivity and resistance. Pestic Sci. 1997, 51 (3): 271-275. 10.1002/(SICI)1096-9063(199711)51:3<271::AID-PS642>3.0.CO;2-#.Google Scholar
- Jacoby WBaZ DM: The enzymes of detoxification. J Biol Chem. 1990, 256: 20715-20718.Google Scholar
- Coleman M, Henricot B, Arnau J, Oliver RP: Starvation-induced genes of the tomato pathogen cladosporium fulvum are also induced during growth in planta. Mol Plant Microbe Interact. 1997, 10 (9): 1106-1109. 10.1094/MPMI.1918.104.22.1686.PubMedGoogle Scholar
- Mariani D, Mathias C, da Silva C, Herdeiro RS, Pereira R, Panek A, Eleutherio E, Pereira MD: Involvement of glutathione transferases, Gtt1and Gtt2, with oxidative stress response generated by H2O2 during growth of Saccharomyces cerevisiae. Redox Rep. 2008, 13 (6): 246-254. 10.1179/135100008X309028.PubMedGoogle Scholar
- Cohen E, Gamliel A, Katan J: Glutathione and glutathione-S-transferase in fungi: Effect of pentachloronitrobenzene and 1-chloro-2,4-dinitrobenzene; Purification and characterization of the transferase from Fusarium. Pestic Biochem Physiol. 1986, 26 (1): 1-9. 10.1016/0048-3575(86)90056-8.Google Scholar
- Mao P, Meshul CK, Thuillier P, Goldberg NRS, Reddy PH: CART peptide is a potential endogenous antioxidant and preferentially localized in mitochondria. PLoS One. 2012, 7 (1): e29343-10.1371/journal.pone.0029343.PubMed CentralPubMedGoogle Scholar
- Joseph-Horne T, Hollomon DW: Functional diversity within the mitochondrial electron transport chain of plant pathogenic fungi. Pest Manag Sci. 2000, 56 (1): 24-30. 10.1002/(SICI)1526-4998(200001)56:1<24::AID-PS71>3.0.CO;2-Y.Google Scholar
- Avila-Adame C, Köller W: Disruption of the alternative oxidase gene in magnaporthe grisea and its impact on host infection. Mol Plant Microbe Interact. 2002, 15 (5): 493-500. 10.1094/MPMI.2002.15.5.493.PubMedGoogle Scholar
- Voulgaris I, O'Donnell A, Harvey LM, McNeil B: Inactivating alternative NADH dehydrogenases: enhancing fungal bioprocesses by improving growth and biomass yield?. Sci Rep. 2012, 2 (322): 1-9.Google Scholar
- Miguez M, Reeve C, Wood PM, Hollomon DW: Alternative oxidase reduces the sensitivity of Mycosphaerella graminicola to QOI fungicides. Pest Manag Sci. 2004, 60 (1): 3-7. 10.1002/ps.837.PubMedGoogle Scholar
- Kaneko I, Ishii H: Effect of azoxystrobin on activities of antioxidant enzymes and alternative oxidase in wheat head blight pathogens Fusarium graminearum and Microdochium nivale. J Gen Plant Pathol. 2009, 75 (5): 388-398. 10.1007/s10327-009-0178-9.Google Scholar
- Thomazella DPT, Teixeira PJPL, Oliveira HC, Saviani EE, Rincones J, Toni IM, Reis O, Garcia O, Meinhardt LW, Salgado I: The hemibiotrophic cacao pathogen Moniliophthora perniciosa depends on a mitochondrial alternative oxidase for biotrophic development. New Phytol. 2012, 194 (4): 1025-1034. 10.1111/j.1469-8137.2012.04119.x.PubMed CentralPubMedGoogle Scholar
- Inoue I, Namiki F, Tsuge T: Plant colonization by the vascular wilt fungus fusarium oxysporum requires FOW1, a gene encoding a mitochondrial protein. The Plant Cell Online. 2002, 14 (8): 1869-1883. 10.1105/tpc.002576.Google Scholar
- Warn PA, Sharp A, Guinea J, Denning DW: Effect of hypoxic conditions on in vitro susceptibility testing of amphotericin B, itraconazole and micafungin against Aspergillus and Candida. J Antimicrob Chemother. 2004, 53 (5): 743-749. 10.1093/jac/dkh153.PubMedGoogle Scholar
- Adachi H, Shimizu K, Hattori H, Tanaka R, Chibana H, Takagi Y, Tomita Y, Kanbe T: Genotyping of candida albicans by fragment analysis of microsatellites combined with 25S rDNA and RPS-based strategies. Nippon Ishinkin Gakkai Zasshi. 2009, 50 (3): 167-174. 10.3314/jjmm.50.167.Google Scholar
- Spicer R, Holbrook NM: Within-stem oxygen concentration and sap flow in four temperate tree species: does long-lived xylem parenchyma experience hypoxia?. Plant Cell Environ. 2005, 28 (2): 192-201. 10.1111/j.1365-3040.2004.01262.x.Google Scholar
- Grahl N, Cramer RA: Regulation of hypoxia adaptation: an overlooked virulence attribute of pathogenic fungi?. Medical Mycology. 2010, 48 (1): 1-15. 10.3109/13693780902947342.PubMed CentralPubMedGoogle Scholar
- Mondego J, Carazzolle M, Costa G, Formighieri E, Parizzi L, Rincones J, Cotomacci C, Carraro D, Cunha A, Carrer H: A genome survey of moniliophthora perniciosa gives new insights into Witches' broom disease of cacao. BMC Genomics. 2008, 9 (1): 548-10.1186/1471-2164-9-548.PubMed CentralPubMedGoogle Scholar
- Zerbino DR, McEwen GK, Margulies EH, Birney E: Pebble and rock band: heuristic resolution of repeats and scaffolding in the velvet short-read de novo assembler. PLoS One. 2009, 4 (12): e8407-10.1371/journal.pone.0008407.PubMed CentralPubMedGoogle Scholar
- Li R, Yu C, Li Y, Lam T-W, Yiu S-M, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.PubMedGoogle Scholar
- Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protocols. 2007, 2 (4): 953-971. 10.1038/nprot.2007.131.PubMedGoogle Scholar
- Horton P, Park K-J, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007, 35 (suppl 2): W585-W587.PubMed CentralPubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.Google Scholar
- Sorensen M, Sanz A, Gómez J, Pamplona R, Portero-Otín M, Gredilla R, Barja G: Effects of fasting on oxidative stress in rat liver mitochondria. Free Radic Res. 2006, 40 (4): 339-347. 10.1080/10715760500250182.PubMedGoogle Scholar
- Hanna SL, Sherman NE, Kinter MT, Goldberg JB: Comparison of proteins expressed by Pseudomonas aeruginosa strains representing initial and chronic isolates from a cystic fibrosis patient: an analysis by 2-D gel electrophoresis and capillary column liquid chromatography–tandem mass spectrometry. Microbiology. 2000, 146 (10): 2495-2508.PubMedGoogle Scholar
- Paes Leme AF, Sherman NE, Smalley DM, Sizukusa LO, Oliveira AK, Menezes MC, Fox JW, Serrano SMT: Hemorrhagic activity of HF3, a snake venom metalloproteinase: insights from the proteomic analysis of mouse skin and blood plasma. J Proteome Res. 2011, 11 (1): 279-291.PubMedGoogle Scholar
- Escalante T, Rucavado A, Pinto AFM, Terra RMS, Gutiérrez JM, Fox JW: Wound exudate as a proteomic window to reveal different mechanisms of tissue damage by snake venom toxins. J Proteome Res. 2009, 8 (11): 5120-5131. 10.1021/pr900489m.PubMedGoogle Scholar
- Eming SA, Koch M, Krieger A, Brachvogel B, Kreft S, Bruckner-Tuderman L, Krieg T, Shannon JD, Fox JW: Differential proteomic analysis distinguishes tissue repair biomarker signatures in wound exudates obtained from normal healing and chronic wounds. J Proteome Res. 2010, 9 (9): 4758-4766. 10.1021/pr100456d.PubMedGoogle Scholar