Combining laser microdissection and RNA-seq to chart the transcriptional landscape of fungal development
© Teichert et al.; licensee BioMed Central Ltd. 2012
Received: 28 May 2012
Accepted: 26 September 2012
Published: 27 September 2012
During sexual development, filamentous ascomycetes form complex, three-dimensional fruiting bodies for the protection and dispersal of sexual spores. Fruiting bodies contain a number of cell types not found in vegetative mycelium, and these morphological differences are thought to be mediated by changes in gene expression. However, little is known about the spatial distribution of gene expression in fungal development. Here, we used laser microdissection (LM) and RNA-seq to determine gene expression patterns in young fruiting bodies (protoperithecia) and non-reproductive mycelia of the ascomycete Sordaria macrospora.
Quantitative analysis showed major differences in the gene expression patterns between protoperithecia and total mycelium. Among the genes strongly up-regulated in protoperithecia were the pheromone precursor genes ppg1 and ppg2. The up-regulation was confirmed by fluorescence microscopy of egfp expression under the control of ppg1 regulatory sequences. RNA-seq analysis of protoperithecia from the sterile mutant pro1 showed that many genes that are differentially regulated in these structures are under the genetic control of transcription factor PRO1.
We have generated transcriptional profiles of young fungal sexual structures using a combination of LM and RNA-seq. This allowed a high spatial resolution and sensitivity, and yielded a detailed picture of gene expression during development. Our data revealed significant differences in gene expression between protoperithecia and non-reproductive mycelia, and showed that the transcription factor PRO1 is involved in the regulation of many genes expressed specifically in sexual structures. The LM/RNA-seq approach will also be relevant to other eukaryotic systems in which multicellular development is investigated.
Fungi are a large group of eukaryotes consisting of a great number of species with a worldwide distribution and great impact on ecology and human society[1, 2]. Fungi comprise both unicellular and multicellular species (yeasts and filamentous fungi, respectively), as well as species capable of both growth forms (dimorphic fungi). All filamentous fungi form a network of vegetative hyphae, called mycelium, that usually grows within or on substrates to acquire nutrients. In addition, many filamentous fungi are capable of developing complex, three-dimensional structures for the generation, protection, and dispersal of spores. Examples are conidiophores for the production of vegetative spores, and fruiting bodies for the production of sexual spores. Fruiting bodies are produced by many ascomycetes and basidiomycetes, and contain a number of specialized cell types that are not present in the vegetative mycelium[3, 4]. The differentiation of these cell types is thought to be orchestrated by spatio-temporal changes in gene expression under the control of regulatory genetic networks. To address the question of developmental regulation of gene expression on a larger scale, several expression studies have been performed with the ascomycetes Gibberella zeae (anamorph Fusarium graminearum), Neurospora crassa, and Sordaria macrospora, all of which belong to the Sordariomycetes and form flask-like fruiting bodies called perithecia. Expression analyses were carried out using high-throughput methods, such as EST sequencing and microarray hybridization[5–13]. In most of the analyses, either time courses of developing mycelia and fruiting bodies were analyzed, or wild-type strains were compared to developmental mutants at certain time points during development. However, the tissues used in these studies usually contained cells from fruiting bodies and vegetative mycelium in varying proportions. One reason for this is that fruiting bodies in ascomycetes are often surrounded by or embedded in vegetative mycelium, from which they are difficult to separate; another reason is that especially the early stages of fruiting body development are quite small (< 50 μm), and even if collected would yield low amounts of material for RNA preparation and subsequent detection methods, such as microarray hybridizations. Therefore, little information is available regarding the spatial control of gene expression patterns, especially for the early stages of fruiting body development.
Laser microdissection (LM) can be used to isolate specific structures consisting of a few cells from samples mounted on microscope slides. LM has been used to isolate cells from animal and plant tissues, and in the case of fungi to study the growth of phytopathogenic or symbiotic species in planta and for the analysis of gene expression differences in single, neighboring hyphae[14–23]. Here, we have established an LM protocol for isolating fruiting body precursors called protoperithecia (young fruiting bodies that are more-or-less spherical without a differentiated neck) from the filamentous fungus S. macrospora. This ascomycete is a model system for the analysis of fungal sexual development and cell differentiation[24, 25]. The genome was sequenced recently using next-generation sequencing techniques, and a number of developmental mutants have already been characterized by classical complementation analyses or by the sequencing of mutant genomes[27, 28]. Prior to the availability of the S. macrospora genome sequence, we had already conducted large-scale expression analyses using cross-species microarray hybridizations with microarrays based on N. crassa cDNAs or oligonucleotides to study gene expression during development in the wild-type and several sterile mutants[8–11]. However, these analyses were limited in sensitivity because less conserved genes give low signal-to-noise ratios in the cross-species array hybridizations. With the genome sequence available, RNA-seq is now the method of choice for large-scale expression analysis. The unprecedented sequencing depths that can be achieved using next-generation sequencing techniques to sequence cDNAs allows much higher sensitivity than microarray hybridization, and the RNA-seq data can also be used for annotation purposes[29–31].
RNA-seq has been used in combination with LM to study gene expression in apical meristems and female gametophytes of Arabidopsis thaliana, in ripening tomato fruits, and in nucleus accumbens neurons in rats[32–35]; however, the combination of LM and RNA-seq has not yet been applied to the analysis of fungal organ-specific transcriptomes. Therefore, in the present study we established an LM protocol for isolating protoperithecia of S. macrospora, and used amplified RNA from the microdissected samples in subsequent RNA-seq analysis. Based on the RNA-seq data, we modeled untranslated regions (UTRs) for more than 50% of the predicted S. macrospora genes, and improved the annotation of roughly 1000 genes. We then compared gene expression patterns in wild-type protoperithecia to those of non-reproductive mycelium from the wild-type, as well as to protoperithecia from the developmental mutant pro1. The sterility of the pro1 mutant is caused by deletion of the transcription factor gene pro1; and one aim of the study was to identify genes that are differentially regulated in protoperithecia, depending on or independent of the PRO1 transcription factor.
Laser microdissection of protoperithecia and RNA-seq analysis
Summary of sequence reads generated in this study
read length in bases
no. of reads
no. of reads mappedto reference genome
% of reads mappedto reference genome
Improving the S. macrospora genome annotation based on RNA-seq data
Mapped RNA-seq reads were visualized in the genome browser Artemis. Though the reads from vegetative and sexual mycelium were evenly distributed along the transcripts, the reads derived from microdissected samples were clustered towards the 3’ end of most transcripts, an example is shown for pro41 in Additional file1 Figure S1. This 3’ bias was expected, because each round of RNA amplification leads to some loss at the 5’ end. In addition, any UV damage during microdissection that might lead to strand cleavage would cause loss of the 5’ portion of the corresponding RNA because reverse transcription was based on polyA tails. In general, a 3’ bias does not hinder quantitative analysis as long as the relative amount of RNA in each sample is preserved, as it does not matter for general quantitation whether the reads from a transcript are evenly spread out or clustered towards one end. However, due to the 3’ bias, many of the reads from the protoperithecial LM samples mapped to the 3’ untranslated regions (UTRs), and therefore would only be taken into account when the 3’ UTR is annotated in the genome sequence. When we started the study, UTRs were annotated for only two S. macrospora genes. Therefore, we used the RNA-seq data to model UTRs and improve the exon-intron structures of the predicted S. macrospora genes in order to make full use of the RNA-seq data in subsequent quantitative expression analyses.
UTR predictions for annotated genes
no. of genes
% of genes
both UTRs found
only 5’ UTR found
only 3’ UTR found
no UTR found
Overview of gene expression across the S. macrospora genome
To obtain an overview of the genome-wide expression as represented by the RNA-seq data, reads that mapped to exons of predicted mRNAs as well as reads that mapped to intergenic regions or introns were counted based on the improved annotation using custom-made Perl scripts ( Additional file1 Figure S4). The majority of reads mapped to annotated exons of protein-coding genes ( Additional file1 Figure S5). This percentage was even higher for the protoperithecial samples, most likely because the two rounds of RNA amplification were based on polyA-dependent reverse transcription constituting an even stronger selection for protein-coding mRNAs. The true percentage of reads that map to intergenic regions might be even lower, because some genes are most likely still missing in the current annotation and not all UTRs for all predicted genes have been annotated yet (see previous section). Therefore, reads mapping to those regions would erroneously be counted as mapping to intergenic regions. A small percentage of reads mapped to introns and might represent incompletely processed transcripts or alternative splicing events, although an analysis of splice sites that were predicted by Tophat showed that, for the majority of genes, no significant alternative splicing could be detected in the different samples (data not shown).
Overall, only 764 predicted genes had no reads map to them in at least one condition ( Additional file1 Figure S6); thus, the majority of genes were expressed in at least one of the four conditions investigated (sexual mycelium, vegetative mycelium, wild-type protoperithecia and pro1 protoperithecia). The majority (295 genes) of the 764 genes were not expressed in wild-type protoperithecia, but were expressed in the other three conditions including protoperithecia from mutant pro1. This finding might indicate that pro1 protoperithecia retain, at least to some degree, properties of non-reproductive mycelia, including the expression of genes that are not required in wild-type protoperithecia. In summary, expression was detected for more than 90% of all annotated genes (version 02) in at least one condition, confirming the high sensitivity of this deep-sequencing approach.
Gene expression in protoperithecia and non-reproductive mycelia
For a quantitative analysis of gene expression in the different samples (vegetative and sexual mycelium, and protoperithecia from wild-type and pro1), sequence reads that mapped to predicted genes were counted using custom-made Perl scripts ( Additional file1 Figure S4) and used for quantitative analysis. Results from LOX and “classical analysis” agreed best with the results from other methods, therefore these approaches were chosen for the final analysis (see Methods and Additional file2). qRT-PCR was used to determine the expression of 17 genes in microdissected samples of wild-type protoperithecia without RNA amplification, and the results were compared to the RNA-seq data. In addition, qRT-PCR results and RNA-seq results were compared for gene expression in vegetative mycelium versus sexual mycelium ( Additional file1 Figure S7). The overall results agree well, and tendencies (up- or down-regulation) are conserved with both methods.
Comparison of gene expression in different samples
with 3’ UTR1
with 3’ UTR1
veg / sex
wt proto / sex
pro1 proto / sex
pro1 proto / wt proto
Analysis of pheromone gene expression in protoperithecia
Analysis of pro1-dependent gene expression in protoperithecia
Genes up-regulated in wt proto/sex and pro1 proto/wt proto are predicted to encode secreted proteins
best Blast hit withknown function
loosenin (Bjerkandera adusta)
riboflavin aldehyde-forming enzyme
MoFLP1 (M. grisea)
We also looked at genes that are physically clustered within the genome, differentially regulated in wild-type protoperithecia, and dependent on pro1 for correct expression in protoperithecia. Physical clustering of co-regulated genes is often found in fungi for genes involved in secondary metabolism, and it can be used as a tool for identifying novel secondary metabolism pathways[50–53]. We found only two instances of clusters that were differentially regulated in wild-type and pro1 protoperithecia. One cluster comprised putative polyketide synthase genes that were described previously and were most likely acquired by horizontal gene transfer. This cluster was down-regulated in wild-type protoperithecia compared to sexual mycelium and up-regulated in pro1 protoperithecia (data not shown). The second group of clustered genes (SMAC_09002 to SMAC_09009) has the opposite expression pattern, namely up-regulated in wild-type protoperithecia compared to sexual mycelium and down-regulated in pro1 protoperithecia. This cluster does not contain a polyketide or non-ribosomal peptide synthase typical for the corresponding gene clusters; however, one unifying theme of this cluster is that three of its genes encode proteins with the predicted domain of unknown function DUF3328 ( Additional file1 Figure S9). Whether genes from this family play a role in fruiting body formation remains to be elucidated.
Fungal fruiting body formation is a complex process that requires coordinated patterns of gene expression in time and space. Even though a number of genes that are essential for this process have been isolated from several model organisms, no unifying theory yet explains the spatio-temporal succession of developmental events leading to the mature fruiting body. One way to learn more about the genes that are active during this process is to look at genome-wide expression patterns at different developmental stages, but this is difficult in many ascomycetes, because fruiting bodies are often rather small (< 500 μm for the mature fruiting body) and difficult to separate from surrounding, non-reproductive hyphae. In the present study, we used LM and RNA-seq to analyze gene expression in protoperithecia from the model organism S. macrospora. To the best of our knowledge, this study is the first time that a combination of these methods has been used for the analysis of fungal gene expression. Our study demonstrates that young fruiting bodies of S. macrospora can be isolated using LM, and RNA extracted from these samples in sufficient amounts for RNA-seq after two rounds of linear amplification. The amplification process largely conserves expression ratios, as demonstrated by comparing the expression of selected genes prior to amplification with the RNA-seq results, which is consistent with results from other organisms where linear RNA amplification was used to prepare samples for microarray hybridization[19, 20]. We also found some overlap with prior microarray experiments in which we compared gene expression in vegetative and sexual mycelia; however, in these experiments we used mycelia grown in defined medium, whereas for RNA-seq analysis, RNA from mycelia grown in defined medium and cornmeal medium were pooled. Therefore, some differences between these experiments might be due to different growth conditions ( Additional file1 Table S2). Furthermore, we demonstrated that the regulatory regions of ppg1 can drive expression of an egfp reporter gene in protoperithecia, as predicted by the RNA-seq analysis. In addition, fluorescence microscopy analysis revealed distinct expression patterns for ppg1 in the outer layers of the protoperithecium. This finding might be consistent with a hypothesis that has been put forward for N. crassa that predicts that pheromones are not only signaling molecules that enable the recognition of mating partners, but that they also play a role in the attachment (“conglutination”) of hyphae forming the rigid outer perithecial wall. A role for pheromones as “molecular glue” might explain the expression of ppg1 in cells of the protoperithecial outer layers in S. macrospora.
In a study of gene expression in several tissues of different metazoans, Hebenstreit et al. found that genes can be grouped into two classes, namely genes with high and low expression, independent of tissue type, species or type of experiment (microarray analysis or RNA-seq). This classification resulted in two distinct peaks when plotting the distribution of gene expression levels. We wondered whether this distribution might also be found in fungi, but plots of the distribution of gene expression levels showed different patterns for our data ( Additional file1 Figure S10, Additional file4). We observed a single main peak in both vegetative and sexual mycelium, whereas the frequency distribution in wild-type and pro1 protoperithecia could be dissected into three peaks. This difference indicates that, in contrast to metazoans, fungal genes might not generally fall into two main classes of expression. One reason might be that in the case of sexual and vegetative mycelium, pooled RNA samples were used from mycelia grown in different types of media. These mycelia might express different sets of genes at high and low levels, and such a mixture would drive overall expression frequencies towards intermediate values, resulting in a single peak as observed. However, the multiple peaks for protoperithecia cannot be explained by a mixture of different samples. The analysis by Hebenstreit indicated that the genes from the high expression group constitute the active and functional transcriptome of the cell, whereas the genes from the low expression group show “leaky” expression. Our data indicate that the situation in fungi might be different, but further analysis will be needed to clarify this point.
In previous studies, we used cross-species microarray hybridization to hybridize S. macrospora targets on N. crassa cDNA or oligonucleotide microarrays[8–11]; however, in these analyses, less than 50% of all genes on the arrays gave a significant signal. The use of RNA-seq dramatically improves detection levels, with more than 90% of all genes being detected in at least one of the sequenced samples. Also, the comparison of gene expression revealed that the overall expression in sexual mycelium is more similar to that of vegetative mycelium than protoperithecia. This finding indicates that gene expression in the sexual mycelium is most likely driven to a large extent by genes expressed in the non-reproductive hyphae making up the bulk of the mycelium and that, in order to study genes specifically expressed in developing fruiting bodies, the microdissection method applied here provides a much better spatial resolution and a much more detailed and specific picture of gene expression during development. Especially weakly expressed, fruiting body-specific genes would most likely not be detected as differentially expressed (or at all) in an expression study using only sexual mycelium. Previous approaches for isolating fruiting bodies for gene expression studies were performed in N. crassa and F. graminearum, using EST sequencing with RNA from mature fruiting bodies, or by analyzing different stages of fruiting bodies by microarray hybridization[6, 7, 13]. The analysis by Hallen et al. was performed with Affymetrix GeneChips for F. graminearum, and signals were detected for nearly 80% of all transcripts, whereas the EST analysis was limited by a comparatively low sequencing depth, and in the other two microarray studies[6, 13], only 10% of all genes gave signals or were represented on the arrays. In all studies, fruiting bodies were harvested by scraping developing structures from a plate, and these preparations might contain an undetermined amount of non-fruiting body mycelia, especially in the early stages of development when fruiting body precursors are small. Therefore, the present analysis of microdissected protoperithecia allowed the analysis of gene expression solely in these structures for the first time. An additional advantage of RNA-seq is that the data can be used also for annotation purposes and, in the case of S. macrospora, allowed the modeling of more than 50% of the UTRs, and the improvement of exon-intron structures for about 1,000 genes (~ 10% of the predicted genes in the genome).
The analysis of gene expression ratios and the 500 genes with the highest number of reads in each of the four sequenced samples showed that expression in protoperithecia from the wild-type and mutant pro1 is more similar to each other than to either vegetative or sexual mycelium, indicating that the transcriptional landscape of protoperithecia is distinct from that of non-reproductive mycelium. However, there are also significant differences between protoperithecia from the wild-type and the sterile mutant pro1 that can form protoperithecia, but not mature fruiting bodies. More than 400 genes were significantly up- or downregulated in pro1 protoperithecia compared to wild-type protoperithecia, and therefore might be direct or indirect targets of PRO1. Among the genes that are dependent on pro1 for correct expression in protoperithecia are the pheromone precursor genes, several genes that might be involved in perithecial wall morphogenesis, and a number of transcription factors. Previous analyses identified several mutants in which the pheromone precursor genes are differentially regulated in sexual mycelium compared to the wild-type[8, 10, 11, 60]; however, no information was available about the spatial regulation of the expression of developmental genes prior to this study. One might hypothesize that pro1 is involved in balancing the expression of genes involved in the formation of the rigid perithecial wall because the pheromone precursor genes and several other genes predicted to be involved in cell-wall biosynthesis are up-regulated in pro1 protoperithecia (Table4).
We have established a combined LM/RNA-seq approach to analyze gene expression in developing fungal sexual structures, and used it to analyze the transcriptome of young fruiting bodies in the wild-type and the sterile mutant pro1 of S. macrospora. pro1, which encodes a transcription factor, is essential for sexual development. To the best of our knowledge, this is the first genome-wide analysis of genes that are dependent on a development-specific transcription factor for correct expression in a defined developmental structure in fungi. Genes that are differentially expressed in protoperithecia are prime candidates for further functional analysis to unravel the spatio-temporal sequence of events leading to the mature fungal fruiting body. Together with three recent studies in Arabidopsis and tomato, as well as rat neurons[32–34], our analysis of fungal development demonstrates the power of a combined approach of LM and RNA-seq to analyze cell type-specific or tissue-specific gene expression in complex, multicellular structures.
Strains and culture conditions
S. macrospora strains used in this study were the wild-type (FGSC 10222,) and the sterile mutant pro1 from the culture collection of the Department of General and Molecular Botany. For propagation, strains were grown on cornmeal medium as described previously. For RNA extraction from vegetative or sexual mycelium, strains were grown in liquid medium as surface cultures (sexual mycelium) or submerged (vegetative mycelium) in cornmeal medium or defined medium as described previously[10, 41].
Laser microdissection was performed with a CellCut Plus system (MMI, Molecular Machines and Industries, Zürich, Switzerland) comprising an Olympus IX81 inverted microscope equipped with a UV laser (355 nm), microscope stage and isolation cup holder. For microdissection, strains were first grown on cellophane-coated cornmeal agar plates for 3 days. Hyphae from these cultures were used to inoculate MMI membrane slides coated with a thin layer of medium (150–200 μl of cornmeal medium with 1% agar). Slides were incubated in a glass Petri dish with approximately 5 ml of water to prevent the samples from drying. Strains were grown for 4–6 days at 25°C in constant light. For fixation, acetone and ethanol gave similar results in RNA extraction and qRT-PCR analysis (data not shown); in the experiments described here, slides were fixed in ethanol at 4°C over night. Prior to microdissection, the slides were air-dried for 30 min, then the mycelium-bearing side was covered with a glass microscope slide and the sample was inserted into the microscope stage with the membrane slide on top (Figure1). Protoperithecia (~20 μm in diameter) were labeled manually using the MMI CellTools software and cut with the UV laser in a distance of ~2-5 μm from the edge of the protoperithecium. This distance minimizes laser damage to the protoperithecium, and ensures that only minimal amounts of unrelated hyphae were isolated. Approximately 100–300 protoperithecia were collected from each slide using an MMI isolation cup (Figure1).
RNA preparation, RNA amplification, and qRT-PCR
RNA was isolated from microdissected protoperithecia using the Arcturus PicoPure kit (Applied Biosystems, Carlsbad, CA, USA) according to the manufacturer’s protocol with the following modifications: 50 μl extraction buffer was pipetted into the isolation cup which was then inverted to cover the cap with the attached microdissected samples and incubated at 42°C for 30 min. Isolation cups were centrifuged (5 min, 800 g) to collect the samples, and then the following procedure was performed three times: cups were frozen in liquid nitrogen for 30 s and then vortexed at room temperature for 30 s. Afterwards, the solution was thawed at 42°C; then extraction was performed using the PicoPure columns according to the manufacturer’s protocol, including a 15 min DNase incubation step (RNase-free DNase Set, Qiagen, Hilden, Germany) after the first washing step. RNA was eluted in 20 μl water. Reverse transcription and qRT-PCR of RNAs isolated from microdissected samples was performed as described previously[10, 41], but instead of 1 μg of RNA, the complete 20 μl of eluted RNA was used for reverse transcription. Primers for qRT-PCR are given in Additional file1 Table S3. Amplification of RNA from microdissected samples was performed with the TargetAmp 2-round aRNA amplification kit 2.0 (Epicentre Biotechnologies, Madison, WI, USA) according to the manufacturer’s protocol with the following modifications. In the second amplification round, first strand cDNA synthesis was performed using primer T7N9 and second strand synthesis using primer oligo-dT(24)-anchored ( Additional file1 Table S3) instead of the primers supplied with the kit in order to obtain polyA-tailed aRNA that could be used directly for Illumina/Solexa library preparation with the same protocol used for total RNA from non-microdissected samples. For vegetative and sexual mycelium, RNA preparation, reverse transcription and qRT-PCR were performed as described previously[10, 41, 62].
Illumina/Solexa sequencing by synthesis and cleaning of primary sequence data
Five micrograms of amplified RNA from microdissected sample was used for Illumina/Solexa sequencing. For sexual and vegetative mycelium, RNA concentrations were quantified photometrically, and equal amounts (50 μg per condition) of RNA from the following growth conditions (independent biological replicates for each growth condition) were pooled: growth on cornmeal medium for 3d, 4d, and 5d, growth on defined medium for 3d, 4d, and 5d; surface cultures for sexual mycelium and shaken cultures for vegetative mycelium. cDNA preparation and Illumina/Solexa sequencing were performed at GATC Biotech (Konstanz, Germany). cDNA preparation was performed using the SMART cDNA library construction kit (Clontech, Mountain View, CA, USA) based on oligo-dT priming for first strand synthesis. For each sample, two independent biological replicates were analyzed (single reads of 35 to 101 bases), each was sequenced in one lane of the GAII (samples SM1 and SM2) or HiSeq 2000 (all other samples, for overview of read numbers, see Table1). Raw sequence data were analyzed and trimmed with custom-made Perl scripts (available athttp://c4-1-8.serverhosting.rub.de/public/software.html). Sequence reads that contained undetermined bases (“N”) were removed. The remaining reads were checked for base quality from the 3’ end, and bases with a quality score of less than 10 (standard Sanger phred scores,) were removed consecutively; reads longer than 20 bases after 3’ trimming were kept for mapping.
Mapping of RNA-seq reads, UTR predictions, and improvement of annotations
The cleaned sequence reads were mapped to the S. macrospora reference genome using Tophat. UTRs were predicted according to the principle shown in Additional file1 Figure S2 with search algorithms described in Additional file1 Method S1. The annotation of predicted open reading frames was checked and improved based on the RNA-seq data using custom-made Perl scripts to implement the algorithm shown in Additional file1 Figure S3 and Additional file1 Method S2. Novel open reading frames were annotated manually based on confirmed splice reads outside of predicted genes. The improved annotation of the S. macrospora genome is available from the ENA database under accession numbers CABT02000001-CABT02001583 and fromhttp://c4-1-8.serverhosting.rub.de/public/.
Quantitative analysis of gene expression based on RNA-seq data
Custom-made Perl scripts were used to determine the number of reads that mapped to each annotated protein-coding gene based on the SAM files with the mapping information (output from Tophat) and using the algorithm shown in Additional file1 Figure S4. Reads were counted stringently in that only reads (> 34 bases) were counted where both ends mapped to the same annotated feature. This approach leads to some loss of reads that map with one end to the UTR of a gene where UTRs are not yet annotated; however, the more stringent counting also prevents spuriously mapped reads from being counted. Raw read counts were used for quantitative analysis with four different methods. The first two approaches were with the Bioconductor packages DESeq and baySeq in the R computing environment (version 2.12.1). However, differential expression estimated by these methods was not in good agreement with previous results from other techniques (microarray, qRT-PCR, and Northern blot, data not shown), most likely due to the fact that the statistical models upon which these methods are based are only valid if the majority of genes (~ 90%) are not differentially expressed. However, in some of our samples, considerably more genes are differentially expressed. Therefore, we employed two other methods to calculate gene expression ratios. One method was based on the LOX program that calculates expression ratios and Bayesian credible intervals and P-values for differential expression. The other method (called “classical analysis”) consists of the calculation of expression ratios, standard deviation, and coefficient of variance from read counts normalized to the total number of read counts for the sample, similar to what was described previously for microarray analyses. In the classical analysis, genes were sorted into five groups (0–4) according to the following criteria: genes in group 4 have ratios of ≤ 0.25 or ≥ 4 in all independent biological replicates, genes in group 3 have a mean ratio of ≤ 0.25 or ≥ 4 and a coefficient of variance < 0.5, genes in group 2 have ratios of ≤ 0.5 or ≥ 2 in all independent biological replicates, genes in group 1 have a mean ratio of ≤ 0.5 or ≥ 2 and a coefficient of variance < 0.5, and group 0 contains all other genes (with the exception of genes for which no ratios could be calculated due to a lack of read coverage, these were not included in the analysis). Overall, results from both methods agreed better with previous results and additional qRT-PCR analyses than results from DESeq or baySeq. To classify genes as differentially expressed, a consensus was determined for each gene based on the results from both the classical and LOX analysis; a gene was labeled as up-regulated (1), down-regulated (−1) or not differentially expressed under the conditions that were compared, when the following criteria were met: (a) for a gene to be classified as differentially regulated in the comparison veg/sex (vegetative mycelium versus sexual mycelium) or pro1 proto/wt proto (pro1 protoperithecia versus wild-type protoperithecia), expression ratios from both classical and LOX analysis had to be > 4 and < 0.25, LOX Bayesian probability for differential expression = 1, and the gene had to be in groups 1–4 in the classical analysis; (b) for a gene to be classified as differentially regulated in the comparison wt proto/sex (wild-type protoperithecia versus sexual mycelium) or pro1 proto/sex (pro1 protoperithecia versus sexual mycelium), expression ratios from both classical and LOX analysis had to be > 8 and < 0.125, LOX Bayesian probability for differential expression = 1, and the gene had to be in groups 1–4 in the classical analysis. We used two different thresholds for differential expression depending on the conditions that were compared, because our analyses showed that lower cutoffs in comparisons of RNA-seq data from microdissected samples and non-microdissected samples resulted in relatively large numbers of false-positives, therefore a more stringent cutoff was used. Results from the classical and LOX analysis as well as the consensus analysis are given in Additional file2. For the analysis of reads that mapped to different genomic regions (e.g., exons, introns, intergenic regions), reads were counted based on the SAM files with the mapping information using custom-made Perl scripts as described above. To determine the distribution of expression frequencies, the coverage for locus tags of protein-coding gene was determined as the average coverage for the bases of the predicted mRNA (normalized to coverage per kilobase per million counted bases in the sample, Additional file4). Curve fitting and clustering of the data by expectation-maximization was performed on the log2-transformed RNA-seq data using the R package mclust.
Multiple alignments were created in CLUSTALX and trimmed with Jalview, and the same alignment was used for analysis by neighbor joining (NJ) and maximum parsimony (MP). Phylogenetic analyses were performed with PAUP version 4.0b10 for Windows (D.L. Swofford, distributed by Sinauer Associates, copyright 2001 Smithsonian Institution) for NJ and MP analyses using 10,000 bootstrap replicates. Consensus trees were graphically displayed with TREEVIEW.
Cloning procedures and transformation of S. macrospora
Plasmid pSNPTppg1 containing egfp flanked by the upstream and downstream sequences (1 kb each) of ppg1 was generated by homologous recombination in yeast as described by Colot et al. based on plasmid pRSnat which contains a nourseothricin resistance cassette for selection in S. macrospora. S. macrospora was transformed as described previously using a combination of Glucanex 100 G (Novozymes A/S, Bagsvaerd, Denmark) and 1.4 U/ml chitinase (ASA Spezialenzyme GmbH) for protoplast generation[28, 71].
For fluorescence microscopy of transformants carrying egfp expression plasmids, strains were grown on glass slides with a thin layer of cornmeal medium and analyzed with an AxioImager fluorescence microscope as described previously[72, 73]. EGFP fluorescence under the control of ppg1 regulatory regions in protoperithecia was observed through a 6% neutral density filter because of very strong fluorescence. All other observations were made without neutral density filters.
The improved version 2 of the S. macrospora genome generated during this work was submitted to the ENA database and is available under the accession numbers CABT02000001-CABT02001583. The RNA-seq reads and results of expression quantification generated in this project were submitted to the GEO database (accession number GSE33668).
The authors would like to thank Swenja Ellßel, Ingeborg Godehardt, Silke Nimtz, Regina Ricke, and Susanne Schlewinski for excellent technical assistance, Dr. Carsten Balczun (Ruhr-Universität Bochum) for help with the BioRad Experion system, and Dr. Jeffrey Townsend (Yale University) for suggesting the use of LOX. This research was supported by the Protein Research Department of the Ruhr-Universität Bochum and the Deutsche Forschungsgemeinschaft (DFG, grants FOR1334, PAK489 KU517/11-1 and NO407/4-1).
- Stajich JE, Berbee ML, Blackwell M, Hibbett DS, James TY, Spatafora JW, Taylor JW: The fungi. Curr Biol. 2009, 19: R840-R845.PubMed CentralView ArticlePubMed
- Hawksworth DL: The magnitude of fungal diversity: the 1.5 million species estimate revisited. Mycol Res. 2001, 105: 1422-1432.View Article
- Pöggeler S, Nowrousian M, Kück U: Fruiting-body development in ascomycetes. The Mycota I. Edited by: Kües U, Fischer R. 2006, Springer, Berlin, Heidelberg
- Bistis GN, Perkins DD, Read ND: Different cell types in Neurospora crassa. Fungal Genet Newsl. 2003, 50: 17-19.
- Nelson MA, Kang S, Braun EL, Crawford ME, Dolan PL, Leonard PM, Mitchell J, Armijo AM, Bean L, Blueyes E, Cushing T, Errett A, Fleharty M, Gorman M, Judson K, Miller R, Ortega J, Pavlova I, Perea J, Todisco S, Trujillo R, Valentine J, Wells A, Werner-Washburn M, Yazzie S, Natvig DO: Expressed sequences from conidial, mycelial, and sexual stages of Neurospora crassa. Fungal Genet Biol. 1997, 21: 348-363.View ArticlePubMed
- Qi W, Kwon C, Trail F: Microarray analysis of transcript accumulation during perithecium development in the filamentous fungus Gibberella zeae (anamorph Fusarium graminearum). Mol Genet Genomics. 2006, 276: 87-100.View ArticlePubMed
- Hallen HE, Huebner M, Shiu SH, Güldener U, Trail F: Gene expression shifts during perithecium development in Gibberella zeae (anamorph Fusarium graminearum), with particular emphasis on ion transport proteins. Fungal Genet Biol. 2007, 44: 1146-1156.View ArticlePubMed
- Klix V, Nowrousian M, Ringelberg C, Loros JJ, Dunlap JC, Pöggeler S: Functional characterization of MAT1-1-specific mating-type genes in the homothallic ascomycete Sordaria macrospora provides new insights into essential and non-essential sexual regulators. Eukaryot Cell. 2010, 9: 894-905.PubMed CentralView ArticlePubMed
- Nowrousian M, Frank S, Koers S, Strauch P, Weitner T, Ringelberg C, Dunlap JC, Loros JJ, Kück U: The novel ER membrane protein PRO41 is essential for sexual development in the filamentous fungus Sordaria macrospora. Mol Microbiol. 2007, 64: 923-937.PubMed CentralView ArticlePubMed
- Nowrousian M, Ringelberg C, Dunlap JC, Loros JJ, Kück U: Cross-species microarray hybridization to identify developmentally regulated genes in the filamentous fungus Sordaria macrospora. Mol Genet Genomics. 2005, 273: 137-149.View ArticlePubMed
- Pöggeler S, Nowrousian M, Ringelberg C, Loros JJ, Dunlap JC, Kück U: Microarray and real time PCR analyses reveal mating type-dependent gene expression in a homothallic fungus. Mol Genet Genomics. 2006, 275: 492-503.View ArticlePubMed
- Lee SH, Lee S, Choi D, Lee YW, Yun SH: Identifcation of the down-regulated genes in a mat1-2- deleted strain of Gibberella zeae, using cDNA subtraction and microarray analysis. Fungal Genet Biol. 2006, 43: 295-310.View ArticlePubMed
- Wang Z, Lehr N, Trail F, Townsend JP: Differential impact of nutrition on developmental and metabolic gene expression during fruiting body development in Neurospora crassa. Fungal Genet Biol. 2012, 49: 405-413.PubMed CentralView ArticlePubMed
- Jiao Y, Tausta SL, Gandotra N, Sun N, Liu T, Clay NK, Ceserani T, Chen M, Ma L, Holford M, Zhang HY, Zhao H, Deng XW, Nelson T: A transcriptome atlas of rice cell types uncovers cellular, funtional and developmental hierarchies. Nat Genet. 2009, 41: 258-263.View ArticlePubMed
- Tremblay A, Li S, Sheffler BE, Matthews BF: Laser capture microdissection and expressed sequence tag analysis of uredinia formed by Phakopsora pachyrhizi, the causal agent of Asian soybean rust. Physiol Mol Plant Pathol. 2009, 73: 163-174.View Article
- Tang W, Coughlan S, Crane E, Beatty M, Duvick J: The application of laser microdissection to in planta gene expression profiling of the maize anthracnose stalk rot fungus Colletotrichum graminicola. Mol Plant Microbe Interact. 2006, 19: 1240-1250.View ArticlePubMed
- Emrich SJ, Barbazuk WB, Li L, Schnable PS: Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Res. 2007, 17: 69-73.PubMed CentralView ArticlePubMed
- Nakazono M, Qiu F, Borsuk L, Schnable PS: Laser-capture microdissection, a tool for the global analysis of gene expression in specific plant cell types: identification of genes expressed differentially in epidermal cells or vascular tissue of maize. Plant Cell. 2003, 15: 583-596.PubMed CentralView ArticlePubMed
- Gomez SK, Harrison MJ: Laser microdissection and its application to analyze gene expression in arbuscular mycorrhizal symbiosis. Pest Manag Sci. 2009, 65: 504-511.View ArticlePubMed
- Fosu-Nyarko J, Jones MG, Wang Z: Application of laser microdissection to study plant-fungal pathogen interactions. Methods Mol Biol. 2010, 638: 153-163.View ArticlePubMed
- Chandran D, Inada N, Hather G, Kleindt CK, Wildermuth MC: Laser microdissection of Arabidopsis cells at the powdery mildew infection site reveals site-specific processes and regulators. Proc Nat Acad Sci USA. 2010, 107: 460-465.PubMed CentralView ArticlePubMed
- Espina V, Wulfkuhle JD, Calvert VS, VanMeter A, Zhou W, Coukos G, Geho DH, Petricoin EF, Liotta LA: Laser-capture microdissection. Nat Prot. 2006, 1: 586-603.View Article
- de Bekker C, Bruning O, Jonker MJ, Breit TM, Wösten HA: Single cell transcriptomics of neighboring hyphae of Aspergillus niger. Genome Biol. 2011, 12: R71-PubMed CentralView ArticlePubMed
- Zickler D: Observing meiosis in filamentous fungi: Sordaria and Neurospora. Methods Mol Biol. 2009, 558: 91-114.View ArticlePubMed
- Engh I, Nowrousian M, Kück U: Sordaria macrospora, a model organism to study fungal cellular development. Eur J Cell Biol. 2010, 89: 864-872.View ArticlePubMed
- Nowrousian M, Stajich JE, Chu M, Engh I, Espagne E, Halliday K, Kamerewerd J, Kempken F, Knab B, Kuo HC, Osiewacz HD, Pöggeler S, Read ND, Seiler S, Smith KM, Zickler D, Kück U, Freitag M: De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sordaria macrospora, a model organism for fungal morphogenesis. PLoS Genet. 2010, 6: e1000891-PubMed CentralView ArticlePubMed
- Kück U, Pöggeler S, Nowrousian M, Nolting N, Engh I: Sordaria macrospora, a model system for fungal development. The Mycota XV, Physiology and Genetics. Edited by: Anke T, Weber D. 2009, Springer, Berlin, Heidelberg, 17-39. 1View Article
- Nowrousian M, Teichert I, Masloff S, Kück U: Whole-genome sequencing of Sordaria macrospora mutants identifies developmental genes. G3 (Bethesda). 2012, 2: 261-270.View Article
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320: 1344-1349.PubMed CentralView ArticlePubMed
- Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, Penkett CJ, Rogers J, Bähler J: Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature. 2008, 453: 1239-1243.View ArticlePubMed
- Nowrousian M: Next-generation sequencing techniques for eukaryotic microorganisms: sequencing-based solutions to biological problems. Eukaryot Cell. 2010, 9: 1300-1310.PubMed CentralView ArticlePubMed
- Matas AJ, Yeats TH, Buda GJ, Zheng Y, Chatterjee S, Tohge T, Ponnala L, Adato A, Aharoni A, Stark R, Fernie AR, Fei Z, Giovannoni JJ, Rose JKC: Tissue- and cell-type specific transcriptome profiling of expanding tomato fruit provides insights into metabolic and regulatory specialization and cuticle formation. Plant Cell. 2011, 23: 3893-3910.PubMed CentralView ArticlePubMed
- Schmid MW, Schmidt A, Klostermeier UC, Barann M, Rosenstiel P, Grossniklaus U: A powerful method for transcriptional profiling of specific cell-types in eukaryotes: laser-assisted microdissection and RNA sequencing. PLoS One. 2012, 7: e29685-PubMed CentralView ArticlePubMed
- Chen H, Liu Z, Gong S, Wu X, Taylor WL, Williams RW, Matta SG, Sharp BM: Genome-wide gene expression profiling of nucleus accumbens neurons projecting to ventral pallidum using both microarray and transcriptome sequencing. Front Neurosci. 2011, 5: 98-PubMed CentralView ArticlePubMed
- Torti S, Fornara F, Vincent C, Andrés F, Nordström K, Göbel U, Knoll D, Schoof H, Coupland G: Analysis of the Arabidopsis shoot meristem transcriptome during floral transition identifies distinct regulatory patterns and a leucine-rich repeat protein that promotes flowering. Plant Cell. 2012, 24: 444-462.PubMed CentralView ArticlePubMed
- Masloff S, Pöggeler S, Kück U: The pro1+ gene from Sordaria macrospora encodes a C6 zinc finger transcription factor required for fruiting body development. Genetics. 1999, 152: 191-199.PubMed CentralPubMed
- Van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH: Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Nat Acad Sci USA. 1990, 87: 1663-1667.PubMed CentralView ArticlePubMed
- Eberwine JH, Yeh H, Miyashiro K, Cao Y, Nair S, Finnell R, Zettel M, Coleman P: Analysis of gene expression in single live neurons. Proc Nat Acad Sci USA. 1992, 89: 3010-3014.PubMed CentralView ArticlePubMed
- Wang E: RNA amplification for successful gene profiling analysis. J Transl Med. 2005, 3: 28-PubMed CentralView ArticlePubMed
- Masloff S, Jacobsen S, Pöggeler S, Kück U: Functional analysis of the C6 zinc finger gene pro1 involved in fungal sexual development. Fungal Genet Biol. 2002, 36: 107-116.View ArticlePubMed
- Nowrousian M, Cebula P: The gene for a lectin-like protein is transcriptionally activated during sexual development, but is not essential for fruiting body formation in the filamentous fungus Sordaria macrospora. BMC Microbiol. 2005, 5: 64-PubMed CentralView ArticlePubMed
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotech. 2010, 28: 511-515.View Article
- Carver T, Harris SR, Berriman M, Parkhill J, McQuillan JA: Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinf. 2012, 28: 464-469.View Article
- Wang B, Guo G, Wang C, Lin Y, Wang X, Zhao M, Guo Y, He M, Zhang Y, Pan L: Survey of the transcriptome of Aspergillus oryzae via massively parallel mRNA sequencing. Nucl Acids Res. 2010, 38: 5075-5087.PubMed CentralView ArticlePubMed
- Zhang Z, López-Giráldez F, Townsend JP: LOX: inferring Level Of eXpression from diverse methods of census sequencing. Bioinf. 2010, 26: 1918-1919.View Article
- Wessels JGH: Fruiting in the higher fungi. Adv Microb Physiol. 1993, 34: 147-202.View ArticlePubMed
- Mayrhofer S, Weber JM, Pöggeler S: Pheromones and pheromone receptors are required for proper sexual development in the homothallic ascomycete Sordaria macrospora. Genetics. 2006, 172: 1521-1533.PubMed CentralView ArticlePubMed
- Quiroz-Castañeda RE, Martínez-Anaya C, Cuervo-Soto LI, Segovia L, Folch-Mallol JL: Loosenin, a novel protein with cellulose-disrupting activity from Bjerkandera adusta. Microb Cell Fact. 2011, 10: 8-PubMed CentralView ArticlePubMed
- Liu TB, Chen GQ, Min H, Lin FC: MoFLP1, encoding a novel fungal fasciclin-like protein, is involved in conidiation and pathogenicity in Magnaporthe oryzae. J Zhejiang Univ Sci B. 2009, 10: 434-444.PubMed CentralView ArticlePubMed
- Hoffmeister D, Keller NP: Natural products of filamentous fungi: enzymes, genes, and their regulation. Nat Prod Rep. 2007, 24: 393-416.View ArticlePubMed
- Brakhage AA, Schroeckh V: Fungal secondary metabolites - strategies to activate silent gene clusters. Fungal Genet Biol. 2011, 48: 15-22.View ArticlePubMed
- Teichert I, Nowrousian M: Evolution of genes for secondary metabolism in fungi. Evolution of fungi and fungal-like organisms, The Mycota XIV. Edited by: Pöggeler S, Wöstemeyer J. 2011, Springer-Verlag, Berlin, Heidelberg, 231-255.View Article
- Nowrousian M: A novel polyketide biosynthesis gene cluster is involved in fruiting body morphogenesis in the filamentous fungi Sordaria macrospora and Neurospora crassa. Curr Genet. 2009, 55: 185-198.View ArticlePubMed
- Nolting N, Pöggeler S: A MADS box protein interacts with a mating-type protein and is required for fruiting body development in the homothallic ascomycete Sordaria macrospora. Eukaryot Cell. 2006, 5: 1043-1056.PubMed CentralView ArticlePubMed
- Colot HV, Park G, Turner GE, Ringelberg C, Crew CM, Litvinkova L, Weiss RL, Borkovich KA, Dunlap JC: A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc Nat Acad Sci USA. 2006, 103: 10352-10357.PubMed CentralView ArticlePubMed
- Son H, Seo YS, Min K, Park AR, Lee J, Jin JM, Lin Y, Cao P, Hong SY, Kim EK, Lee SH, Cho A, Lee S, Kim MG, Kim Y, Kim JE, Kim JC, Choi G, Yun SH, Lim JY, Kim M, Lee YH, Choi YD, Lee YW: A phenome-based functional analysis of transcription factors in the cereal head blight fungus. Fusarium graminearum. PLoS Pathog. 2011, 7: e1002310-View ArticlePubMed
- Han KH, Han KY, Yu JH, Chae KS, Jahng KY, Han DM: The nsdD gene encodes a putative GATA-type transcription factor necessary for sexual development of Aspergillus nidulans. Mol Microbiol. 2001, 41: 299-309.View ArticlePubMed
- Kim H, Metzenberg RL, Nelson MA: Multiple functions of mfa-1, a putative pheromone precursor gene of Neurospora crassa. Eukaryot Cell. 2002, 1: 987-999.PubMed CentralView ArticlePubMed
- Hebenstreit D, Fang M, Gu M, Charoensawan V, van Oudenaarden A, Teichmann SA: RNA sequencing reveals two major classes of gene expression levels in metazoan cells. Mol Syst Biol. 2011, 7: 497-PubMed CentralView ArticlePubMed
- Gesing S, Schindler D, Fränzel B, Wolters D, Nowrousian M: The histone chaperone ASF1 is essential for sexual development in the filamentous fungus Sordaria macrospora. Mol Microbiol. 2012, 84: 748-765.View ArticlePubMed
- Esser K: Cryptogams - Cyanobacteria, Algae, Fungi, Lichens. 1982, Cambridge University Press, London
- Engh I, Nowrousian M, Kück U: Regulation of melain biosynthesis via the dihydroxynaphtalene pathway is dependent on sexual development in the ascomycete Sordaria macrospora. FEMS Microbiol Lett. 2007, 275: 62-70.View ArticlePubMed
- Cock PJA, Fields C, Goto N, Heuer ML, Rice PM: The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucl Acids Res. 2010, 38: 1767-1771.PubMed CentralView ArticlePubMed
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-PubMed CentralView ArticlePubMed
- Hardcastle TJ, Kelly KA: baySeq: empirical Bayesian methods for identifying differential expression in sequence count data. BMC Bioinf. 2010, 11: 422-View Article
- Fraley C, Raftery AE: Model-based clustering, discriminant analysis, and density estimation. J Amer Stat Assoc. 2002, 97: 611-631.View Article
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 24: 4876-4882.View Article
- Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview Version 2 - a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25: 1189-1191.PubMed CentralView ArticlePubMed
- Hall BG: Phylogenetic trees made easy. 2004, Sinauer Associates, Sunderland, 2
- Page R: TREEVIEW: an application to display phylogenetic trees on personal computers. Appl Biosci. 1996, 12: 357-358.
- Nowrousian M, Masloff S, Pöggeler S, Kück U: Cell differentiation during sexual development of the fungus Sordaria macrospora requires ATP citrate lyase activity. Mol Cell Biol. 1999, 19: 450-460.PubMed CentralView ArticlePubMed
- Nowrousian M, Piotrowski M, Kück U: Multiple layers of temporal and spatial control regulate accumulation of the fruiting body-specific protein APP in Sordaria macrospora and Neurospora crassa. Fungal Genet Biol. 2007, 44: 602-614.View ArticlePubMed
- Engh I, Würtz C, Witzel-Schlömp K, Zhang HY, Hoff B, Nowrousian M, Rottensteiner H, Kück U: The WW domain protein PRO40 is required for fungal fertility and associates with Woronin bodies. Eukaryot Cell. 2007, 6: 831-843.PubMed CentralView ArticlePubMed
- Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz H-R, Ceric G, Forslund K, Eddy SR, Sonnhammer ELL, Bateman A: The Pfam protein families database. Nucl Acids Res. 2008, 36 (Suppl 1): D281-288.PubMed CentralPubMed
- Eddy SR: Profile Hidden Markov Models. Bioinformatics. 1998, 14: 755-763.View ArticlePubMed
- 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.View ArticlePubMed
- Horton P, Park KJ, Obayashi T, Nakai K: Proceedings of the 4th Annual Asia Pacific Bioinformatics Conference APBC06. Protein subcellular localization prediction with WoLF PSORT. 2006, Taipei, Taiwan, 39-48.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.