Strand-Specific RNA-Seq Analyses of Fruiting Body Development in Coprinopsis cinerea

The basidiomycete fungus Coprinopsis cinerea is an important model system for multicellular development. Fruiting bodies of C. cinerea are typical mushrooms, which can be produced synchronously on defined media in the laboratory. To investigate the transcriptome in detail during fruiting body development, high-throughput sequencing (RNA-seq) was performed using cDNA libraries strand-specifically constructed from 13 points (stages/tissues) with two biological replicates. The reads were aligned to 14,245 predicted transcripts, and counted for forward and reverse transcripts. Differentially expressed genes (DEGs) between two adjacent points and between vegetative mycelium and each point were detected by Tag Count Comparison (TCC). To validate RNA-seq data, expression levels of selected genes were compared using RPKM values in RNA-seq data and qRT-PCR data, and DEGs detected in microarray data were examined in MA plots of RNA-seq data by TCC. We discuss events deduced from GO analysis of DEGs. In addition, we uncovered both transcription factor candidates and antisense transcripts that are likely to be involved in developmental regulation for fruiting

Strand-Specific RNA-Seq Analyses of Fruiting Body Development in Coprinopsis cinerea.

The basidiomycete fungus Coprinopsis cinerea is an important model system for multicellular development. Fruiting bodies of C. cinerea are typical mushrooms, which can be produced synchronously on defined media in the laboratory. To investigate the transcriptome in detail during fruiting body development, high-throughput sequencing (RNA-seq) was performed using cDNA libraries strand-specifically constructed from 13 points (stages/tissues) with two biological replicates. The reads were aligned to 14,245 predicted transcripts, and counted for forward and reverse transcripts. Differentially expressed genes (DEGs) between two adjacent points and between vegetative mycelium and each point were detected by Tag Count Comparison (TCC). To validate RNA-seq data, expression levels of selected genes were compared using RPKM values in RNA-seq data and qRT-PCR data, and DEGs detected in microarray data were examined in MA plots of RNA-seq data by TCC. We discuss events deduced from GO analysis of DEGs. In addition, we uncovered both transcription factor candidates and antisense transcripts that are likely to be involved in developmental regulation for fruiting

Strand-Specific RNA-Seq Analyses of Fruiting Body Development in <i>Coprinopsis cinerea</i>

<div><p>The basidiomycete fungus <i>Coprinopsis cinerea</i> is an important model system for multicellular development. Fruiting bodies of <i>C</i>. <i>cinerea</i> are typical mushrooms, which can be produced synchronously on defined media in the laboratory. To investigate the transcriptome in detail during fruiting body development, high-throughput sequencing (RNA-seq) was performed using cDNA libraries strand-specifically constructed from 13 points (stages/tissues) with two biological replicates. The reads were aligned to 14,245 predicted transcripts, and counted for forward and reverse transcripts. Differentially expressed genes (DEGs) between two adjacent points and between vegetative mycelium and each point were detected by Tag Count Comparison (TCC). To validate RNA-seq data, expression levels of selected genes were compared using RPKM values in RNA-seq data and qRT-PCR data, and DEGs detected in microarray data were examined in MA plots of RNA-seq data by TCC. We discuss events deduced from GO analysis of DEGs. In addition, we uncovered both transcription factor candidates and antisense transcripts that are likely to be involved in developmental regulation for fruiting.</p></div

Samples used for RNA-seq.

<p>Samples used for RNA-seq.</p

GO analysis of DEGs in each transition.

<p>GO analysis of DEGs in each transition.</p

A model for significant developmental transitions based on GO analysis.

<p>Notable events among top 3 categories are depicted in fruiting. Events detected by up-regulated DEGs and down-regulated DEGs are indicated in the upper part with upward red arrows and in the lower part with downward blue arrows, respectively. The asterisks and bidirectional light green arrows indicate that similar annotation terms are enriched in both up-regulated and down-regulated DEGs.</p

Schematic diagram of fruiting body development in <i>C</i>. <i>cinerea</i>.

<p>The 13 stages/tissues (numbers in circles) were selected to investigate the transcriptome by RNA-seq. Developmental and cellular events in each stage/tissue are shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0141586#pone.0141586.t001" target="_blank">Table 1</a>.</p

Transcription factor candidates in <i>C</i>. <i>cinerea</i> and changes in the expression levels.

<p>(A) Relative distribution of Pfam domains in the TFCs. A total of 848 TFCs were sorted on the basis of Pfam domains, and the number and percentage of each Pfam are indicated. (B) A histogram dividing 848 transcription factor candidates into bins based on their RPKM values. Red line shows a total number of the TFC genes, 848. Red arrows indicate notable changes.</p