Transcriptome of Protoplasts Reprogrammed into Stem Cells in Physcomitrella patens

Abstract

Background

Differentiated plant cells can retain the capacity to be reprogrammed into pluripotent stem cells during regeneration. This capacity is associated with both cell cycle reactivation and acquisition of specific cellular characters. However, the molecular mechanisms underlying the reprogramming of protoplasts into stem cells remain largely unknown. Protoplasts of the moss Physcomitrella patens easily regenerate into protonema and therefore provide an ideal system to explore how differentiated cells can be reprogrammed to produce stem cells.

Principal findings

We obtained genome-wide digital gene expression tag profiles within the first three days of P. patens protoplast reprogramming. At four time-points during protoplast reprogramming, the transcript levels of 4827 genes changed more than four-fold and their expression correlated with the reprogramming phase. Gene ontology (GO) and pathway enrichment analysis of differentially expressed genes (DEGs) identified a set of significantly enriched GO terms and pathways, most of which were associated with photosynthesis, protein synthesis and stress responses. DEGs were grouped into six clusters that showed specific expression patterns using a K-means clustering algorithm. An investigation of function and expression patterns of genes identified a number of key candidate genes and pathways in early stages of protoplast reprogramming, which provided important clues to reveal the molecular mechanisms responsible for protoplast reprogramming.

Conclusions

We identified genes that show highly dynamic changes in expression during protoplast reprogramming into stem cells in P. patens. These genes are potential targets for further functional characterization and should be valuable for exploration of the mechanisms of stem cell reprogramming. In particular, our data provides evidence that protoplasts of P. patens are an ideal model system for elucidation of the molecular mechanisms underlying differentiated plant cell reprogramming.

About the Authors:

Lihong Xiao

Affiliation: College of Life Science, Capital Normal University, Beijing, People's Republic of China

Liechi Zhang

Affiliation: College of Life Science, Capital Normal University, Beijing, People's Republic of China

Ge Yang

Affiliation: College of Life Science, Capital Normal University, Beijing, People's Republic of China

Honglin Zhu

Affiliation: Beijing Computing Center, Beijing, People's Republic of China

Yikun He

* E-mail: [email protected]

Affiliation: College of Life Science, Capital Normal University, Beijing, People's Republic of China

Introduction

Cell reprogramming is an important biological phenomenon whereby cells regress from a specialized, differentiated state to a simple, undifferentiated cell type reminiscent of stem cells with the capability for both self-renewal and to give rise to most other cell types in multicellular organisms [1], [2]. Asymmetric division of a stem cell generates two different daughter cells: a self-renewed stem cell daughter that retains the stem cell's pluripotent characteristics, and a differentiated non-stem cell daughter [3]. Stem cells in plant shoot and root meristems are maintained throughout the life of plants and produce somatic daughter cells that make up the body of plants [4]. Commonly, specialized cells are formed by a one-way process as a fertilized egg develops into an adult, and the cells become increasingly, and normally irreversibly, committed to their fate. However, certain experimental procedures can reverse the cell differentiation process and cause cells to acquire a new fate by reprogramming, a term that describes a switch in nuclear gene expression in one kind of cell to induce it to differentiate into a different cell type. A distinctive feature of cell reprogramming is the withdrawal from a given differentiated state into a stem cell-like state that confers pluripotentiality, a process that precedes the switch in cell fate [5]. This process underlies the totipotency, regeneration and formation of new stem cells. Elucidation of how cell reprogramming takes place is important to help us understand the mechanisms by which cell division and differentiation occur.

Reprogramming of a differentiated cell to become a pluripotent stem cell is frequently observed in plants and is more easily induced in plants than in animals. Differentiated plant cells, in contrast to those of animals, hold multiple developmental potentialities during development and retain a plasticity that enables dedifferentiation [6]. However, the genetic and molecular bases of this difference between plant and animal cells are mostly unknown.

Recently, artificial expression of two transcription factors, Oct4 and Sox2, together with other factors made it possible to reprogram differentiated somatic cells into pluripotent stem cells in mice and humans [7]. The study of cellular reprogramming in animals is limited because of the lack of a suitable, convenient experimental system [8]. Plant protoplasts (plant cells devoid of cell walls) provide an outstanding experimental tool for the study of the biochemical and molecular basis of cellular reprogramming into stem cells [6], [9], [10].

Application of phytohormones, such as auxins and cytokinins, stimulate protoplasts from different tissues to reenter the cell cycle, proliferate, and undergo regenerative processes to give rise to new plantlets [11]–[13]. Using a plant protoplast system, Zhao et al. [9] demonstrated that protoplasts can be isolated easily from fully differentiated, non-dividing mesophyll cells of tobacco leaves, and reenter the cell cycle and proliferate following treatment with auxin and cytokinin. These authors also found that the reprogramming of tobacco mesophyll cells proceeds by two functionally distinct phases of chromatin decondensation: the first is a transitory phase that confers competence for a switch in cell fate followed, under appropriate hormonal conditions, by a second phase that represents commitment to the mitotic cycle. Subsequent studies of protoplast reprogramming concentrated on the relationship between transcription and the structure of chromatin that is involved in heterochromatin decondensation and histone modification [9], [10], [14], [15].

Compared to vascular plants, bryophytes and ferns have a single stem cell in the protonema tip and leafy shoot apex, lack a stem cell niche and organization center, and possess an accessible haploid and relatively simple structure [16], [17], and thus represent a simpler experimental system. Although they diverged hundreds of millions of years ago, bryophytes share similar fundamental genetic and physiological features with seed plants [18] and are phylogenetically intermediate between algae and seed plants [19]. Among bryophytes, the moss Physcomitrella patens has emerged recently as the bryophyte model of choice for studies of development, genetics and stress responses [20]–[23].

The P. patens apical stem cell system has received much attention following its description in 2007 (). Specifically, differentiated cells from any part of the gametophyte or sporophyte, including an excised leaf, protonemal cells or freshly isolated protoplasts, can be easily reprogrammed into protonemal apical stem cells without exogenous phytohormone treatment within a few days [24]. A recent study on reprogramming of excised leaf cells indicated that a cell cycle protein kinase A (CDKA) regulates cell division and acquisition of new cell characteristics in the reprogramming of differentiated cells to become stem cells in plants [25]. However, the factors that coordinate cell cycle reactivation with acquisition of other cellular characteristics during protoplast reprogramming into stem cells have not been determined.

In this study, we utilized the P. patens protoplast system to explore the mechanisms and key candidate regulators involved in stem cell reprogramming. The objectives of the present study were to: (1) characterize changes in gene expression associated with protoplast reprogramming, (2) investigate the molecular mechanisms responsible for protoplast reprogramming, and (3) identify candidate genes and key factors involved in protoplast reprogramming into stem cells. By combining the P. patens protoplast system with a digital gene expression tag profiling (DGEP) strategy, we obtained spatiotemporal-specific gene regulation models for protoplast reprogramming. These results provided a comprehensive catalogue of gene expression changes during protoplast reprogramming, from which potentially key regulatory factors can be mined.

Results

Morphogenesis of protoplasts reprogrammed into stem cells

To investigate how protoplasts were reprogrammed into stem cells, six-day-old subcultured P. patens protonemata (Fig. 1B) were used to establish an efficient and reproducible ‘protoplast system’ in P. patens. When protonemal tissue was treated to form protoplasts, the cells changed their nuclear program and lost their differentiated state (Fig. 1C). Protoplasts freshly isolated from protonemata were round and green because of the presence of chloroplasts. Upon culture in the dark, more than 90% of the protoplasts developed a new cell wall within 24 h of culture (Fig. 1D). The new polar axes were re-established within 48 h and up to 87% of the cells were pyriform (Fig. 1E), which indicated that protoplasts are reprogrammed to acquire at least some stem cell characteristics before mitosis. After 72 h, 85% of the originally plated protoplasts divided asymmetrically and yielded chloronema that contained two to three cells, of which the apical cell was a self-renewing stem cell and other cell(s) were differentiated non-stem cell(s) (Fig. 1F). Subsequently, the cultures were transferred to a propagation and regeneration medium (BCDAG medium) for formation of protonema clones, the apical cells of which were stem cells. Using this protoplast system, we followed alterations in the cell cycle and DGEP analysis at subsequent time intervals up to 72 h after protoplast isolation.

[Figure omitted. See PDF.]

Figure 1. Morphology and evaluation of cell cycle by flow cytometry.

(A) Two-month-old leafy shoots. (B) Six-day-old protonemata. (C) Freshly prepared protoplast. (D) Cells cultured for 24 hours. (E) Cells cultured for 48 hours. (F) Cells cultured for 72 hours. The stem cells are indicated by arrows in each stage. The S phase of cell cycle is represented by the area between the two red broken lines, and the G2/M phase is indicated by a blue star in each stage.

https://doi.org/10.1371/journal.pone.0035961.g001

Protoplast reprogramming into stem cells occurs within 48 hours

To evaluate the time course of protoplast reprogramming into stem cells, we measured the DNA content of nuclei isolated from freshly prepared (t-0) protoplasts and protoplast-derived cells cultured for 24, 48 or 72 hours (t-24, t-48 and t-72), respectively, by flow cytometry (FCM). In order to ascertain the cell cycle phase of protoplasts, we used DAPI-stained nuclei from two-month-old leafy shoots and six-day-old protonemata of P. patens as a control (Fig. 1A–B). Under the experimental conditions, six-day-old protonemata mainly contained chloronema cells and rarely caulonemal cells (Fig. 1B). Nuclei from two-month-old leafy shoots exhibited three peaks that corresponded to G1- and G2-phase as well as polyploidization, which may be the results of endoduplication (Fig. 1A). Six-day-old P. patens protonemal nuclei displayed double peaks that corresponded to G1- and G2/M-phase DNA content, respectively (Fig. 1B). Comparison of the FCM results for the controls, we can speculate that most nuclei were in the G1 phase of the cell cycle with a only small number of cells in the S and G2/M phases, which differed from previous reports that chloronema cells were arrested in the G2 phase [26], [27].

Comparison of the flow cytometry histograms of DAPI-stained nuclei isolated from t-0 protoplasts and protoplast-derived cells (t-24, t-48 and t-72 cells) with nuclei prepared from two-month-old leafy shoots and six-day-old protonemata indicated that nuclei prepared from both t-0 protoplasts (Fig. 1C; G1 nuclei) and t-24 cells (Fig. 1D; G1 nuclei) reproducibly showed a single peak. These results indicated that t-0 protoplasts and t-24 protoplast-derived cells underwent synchronization. The t-24 protoplast cultures were highly synchronized and almost all cells were in the G1 phase. Few t-48 cells showed an increase in fluorescence intensity compared with t-24 protoplasts (Fig. 1E). A portion of the t-48 cells were in the S phase and the other cells were in the G2/M phase (Fig. 1). This result implied that the protoplasts had been reprogrammed and re-entered the cell cycle. Features of t-72 cells were similar to those of t-48 cells (Fig. 1). From integration of FCM analysis with morphogenetic observations during protoplast culture, we concluded that protoplast reprogramming into stem cells occurred within 48 h.

Tag identification and quantification and depth of sequencing

To obtain global patterns of gene expression during protoplast reprogramming, RNA extracted from fresh protoplasts and cells cultured for 24, 48 and 72 h was used for DGEP analyses. More than 3.2 million raw tags (Table 1) were sequenced using the cDNA library derived from fresh protoplasts and cells cultured for 24, 48 and 72 h (Table 1). Custom Perl scripts were used for adaptor trimming and read parsing. Before mapping these tag sequences to the reference sequence, low-quality tags (tags containing ‘N’ and adaptor sequences) were filtered. To increase the robustness of the approach, single-copy tags in the four libraries were excluded [28]. The distribution of distinct clean tag counts over different tag abundance categories showed very similar tendencies for all four libraries (Fig. S1) [29], [30]. Common and specific tags within and among samples are shown in Fig. S1.

[Figure omitted. See PDF.]

Table 1. Major characteristics of DGEP libraries and mapping information.

https://doi.org/10.1371/journal.pone.0035961.t001

Saturation of the library was determined by identification of unique tags. Sequencing reached saturation when no new unique tags were detected. The results shown in Fig. S2 indicated that all four sampling libraries were sequenced to saturation, and thus a full representation of the transcripts in the experimental conditions was obtained. In the four libraries fewer unique tags were identified as the number of sequencing tags increased, and reached a plateau shortly after 2 million tags were sequenced and a negligible increase in the number of genes detected in the four libraries was observed.

Mapping of short reads to the reference genome and detection of differentially expressed genes

Bowtie 0.12.7 was used to map unique consensus sequence tags (a total of two or more reads from all libraries) to the P. patens reference genome. Bowtie is an ultrafast, memory-efficient short-read aligner [31]. Bowtie indexes the genome with a Burrows-Wheeler index to keep its memory footprint small. This method performs effectively with DGEP data sets, which are reduced in size and complexity since reads are collapsed to unique tags before mapping. Finally, a preprocessed database of all possible CATG+17-nt tag sequences was created using reference gene sequences. All clean tags were mapped to the reference sequences and allowed no more than 1-nt mismatches. Clean tags mapped to reference sequences from multiple genes were filtered (Table 1). For genes that have multitags found in Solexa tags, the sum of all tags was considered as the gene expression value.

To compare gene expression profiles, we employed the TMM method from edgeR (empirical analysis of digital gene expression in R) to normalize the tag distribution per library and determine significance values for differentially expressed genes based on their relative abundance, which reflected the difference in number of tags between each two libraries. The edgeR algorithm uses an empirical Bayes approach to improve power in small sample sizes [32]–[34]. This approach accounts for biological and technical variation and has been implemented for tag-based data sets where small numbers of replicates are tested and standard errors disperse further from the mean at low versus high levels of expression [33], [35]. We used P values≤0.01, false discovery rates (FDR)≤0.01, and |logFC|≥2 as the threshold to determine differentially expressed genes (DEGs). The transcripts detected in the four libraries are shown in Fig. S3. The red dots were deemed to be differentially expressed transcripts. The black dots represented transcripts that were arbitrarily designated as ‘no difference in expression’ between the two comparative libraries (Fig. S3). As a result, 4827 differences were identified (Fig. 2; Table S1). The full transcriptomic data set was deposited in the GEO database (accession no. GSE36117).

[Figure omitted. See PDF.]

Figure 2. Genes differentially expressed at different time-points during protoplast reprogramming into stem cells.

(A) Genes differentially expressed during specific phases of protoplast reprogramming into stem cells were separated into two groups on the basis of whether they were significantly up-regulated or down-regulated. (B) Venn diagrams showing the number of differentially expressed genes during specific time-points of protoplast reprogramming into stem cells.

https://doi.org/10.1371/journal.pone.0035961.g002

Quantitative real-time PCR analysis

To test the reliability, accuracy and reproducibility of the DGEP data, quantitative real-time PCR analysis was applied to validate the expression pattern of eight randomly selected genes (Table S2). Among these genes, a homologue of the Arabidopsis stem cell maintenance gene WOX13 (Pp1s224_106V6, which we designated PpWOX13a) was included, which exhibited a relatively high expression level in t-0 protoplasts and 24-hour cultures and subsequently was down-regulated. A putative heat shock cognate protein-encoding gene (Pp1s97_279V6), for which expression was relatively stable in all four libraries, was chosen as a reference gene for data normalization. For each gene, the trends in quantitative real-time PCR (qRT-PCR) expression profile were in agreement (i.e., up- or down-regulation) with that obtained by DGEP analysis (Fig. 3, Table S3).

[Figure omitted. See PDF.]

Figure 3. Digital gene expression tag profiling and quantitative real-time PCR analysis of the expression of five randomly selected genes.

All real-time PCR reactions were repeated three times and the data are presented as the mean ± SD. The x-axis indicates the sampling time-points and cell types. The y-axis shows the expression levels: the left bar (red color) shows tag number per million tags by DGEP and the right (blue color) shows the relative expression level by qRT-PCR.

https://doi.org/10.1371/journal.pone.0035961.g003

Global gene expression profiles during protoplast reprogramming

Among all DEGs, 1095 genes showed a change in expression level in t-24 cells compared to t-0 protoplasts, of which 383 genes were up-regulated and 712 genes were down-regulated. A total of 2059 DEGs were detected between culture for 24 h to 48 h and 1074 genes showed transcript accumulation. Of the 1870 DEGs detected between 48 h to 72 h culture, 969 genes were up-regulated and 901 genes were down-regulated (Fig. 2A; Table S1). Among the DEGs, 4630 genes were preferentially expressed in a single comparison, 197 genes were expressed in two comparisons, and none were shared in all three comparisons (Fig. 2B).

To understand further details of the DEGs, significant enriched gene ontology (GO) terms were identified according to their P value and enrichment factor. The top 15 significant enriched GO terms in each group are summarized in Table S4 and Fig. S4. The DEGs were significantly enriched in the processes of a variety of biotic and abiotic responses (including responses to salt stress, cold, cadmium ions and bacterial infection), photosynthesis-related processes, glycolysis, ribosome biosynthesis and translation. However, the number, expression level and description of DEGs differed at different stages during protoplast reprogramming. These results indicated that these biological processes played important roles in protoplast reprogramming. The expression changes of DEGs involved in photosystem II repair and responses to blue light, high light intensity and heat occurred within the first 24 h of protoplast reprogramming. In addition, alteration in expression of DEGs involved in thylakoid membrane organization, protein importation into chloroplast stroma, photorespiration, chloroplast organization, and response to abscisic acid (ABA) might also be necessary for stem cells to re-enter the cell cycle. Similarly, a requirement for stem-cell division and differentiation was reflected in expression of DEGs involved in purine nucleotide biosynthesis, embryo development, fatty acid biosynthesis, protein folding, photoinhibition and response to wounding. Significantly enriched molecular functions of DEGs were 3-chloroallylaldehyde dehydrogenase activity, poly (U) RNA binding, chlorophyll binding, structural constituent of ribosome, copper ion binding and catalytic activity. Some DEGs that functioned in glyceraldehydes-3-phosphate dehydrogenase activity, p-p-bond-hydrolysis-driven protein factor activity, translation initiation factor activity binding and rRNA binding were significantly enriched from 24 and 72 h. These results indicated DEGs that functioned in binding, protein synthesis and catalytic activity were important during protoplast reprogramming. Furthermore, DEGs that functioned in phosphoglycerate kinase activity, aminomethyltransferase activity, acetyl-coA carboxylase activity, glutamate-ammonia ligase activity, oxidoreductase activity that acted on paired donors, glutathione binding, ATPase activity coupled to transmembrane movement of substances and calmodulin binding (from 0 to 24 h), alpha-amylase activity, NADH dehydrogenase (ubiquinone) activity and GTP binding (from 24 to 48 h), and phosphoglucomutase activity, ubiquitin binding, protein histidine kinase activity and cobalt ion binding (from 48 to 72 h), might also be vital for the transition from protoplasts to stem cells and stem-cell division and differentiation.

The cellular components of DEGs were mainly enriched in the chloroplast stroma, chloroplast thylakoid, and chloroplast envelope, followed by the vacuolar membrane, ribosome, stromule and apoplast, which indicated that the molecular biological reactions mainly occurred in the chloroplast, vacuolar membrane, ribosome, stromule and apoplast. In addition, some DEGs functioned in mitochondria from 24 to 72 h. Few DEGs that functioned in the cell wall were specifically expressed.

To investigate the functions of DEGs during protoplast reprogramming, significant enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified according to the P value and enrichment factor. The top 15 significant enriched KEGG pathways in each comparison are summarized in Table S5 and Fig. S5. Hierarchical clustering of significant pathways showed that the photosynthesis, butanoate metabolism, ribosome, glycolysis/gluconeogenesis and pyruvate metabolism pathways were involved in all processes of protoplast reprogramming (Fig. 4), which indicated that these pathways are essential for cell survival, division and differentiation. Some genes involved in the nitrogen metabolism, alanine, aspartate and glutamate metabolism and selenoamino acid metabolism pathways were specifically enriched within the first 24 h of protoplast reprogramming, early stage of stem cell reprogramming. Porphyrin and chlorophyll metabolism, lysine degradation, fatty acid metabolism, citrate cycle (TCA cycle), C5-branched dibasic acid metabolism, glyoxylate and dicarboxylate metabolism, oxidative phosphorylation, and pentose phosphate pathways and valine, leucine and isoleucine biosynthesis were specifically enriched from 24 to 48 h, a stage of stem cell re-entering cell cycle. Thus these pathways were indicated to be closely associated with the cell fate transition during protoplast reprogramming into stem cells. Meanwhile, spliceosome, peroxisome, glutathione metabolism and mRNA surveillance pathways were specifically enriched from 48 to 72 h (Figs. 4–5), which indicated that these pathways were closely related with the protonemal regeneration after stem cell reprogrammed.

[Figure omitted. See PDF.]

Figure 4. Heatmap of significantly enriched pathways.

The yellow and red color shows the p-value of significantly enriched pathways.

https://doi.org/10.1371/journal.pone.0035961.g004

[Figure omitted. See PDF.]

Figure 5. Specific enriched pathways at specific stages during protoplast reprogramming into stem cells.

https://doi.org/10.1371/journal.pone.0035961.g005

Key candidate factors involved in protoplast reprogramming

To elucidate the molecular controls of protoplast reprogramming into stem cells in P. patens, we further transformed the comparison of two sequential time-points into a comparison with the expression level in t-0 protoplasts as a common reference. Based on the transformed comparison data and on the basis of mRNA accumulation trends, we performed a K-means clustering analysis using MEV, which partitioned the 4827 DEGs into six clusters (Fig. 6) [36]. The gene number in each pattern ranged from 119 (pattern B) to 1817 (pattern F). Patterns A to D contained genes preferentially expressed in fresh protoplasts (pattern A), cells that had regenerated a new cell wall (pattern B), reprogrammed stem cells (pattern C) and regenerated 2- to 3-celled chloronemata (pattern D) (Fig. 6).

[Figure omitted. See PDF.]

Figure 6. Patterns of gene expression by K-means cluster analysis in the developing gametophyte of P. patens.

Differentially expressed genes across all four time-points were grouped into six patterns using the K-means clustering algorithm. The y-axis gives the tag count (on a log₂ scale) of differentially expressed genes. Each line represents a different gene.

https://doi.org/10.1371/journal.pone.0035961.g006

Compared with cells of six-day-old protonemata, freshly isolated protoplasts differed in shape and gene expression because of removal of cell walls. Most of the 270 preferentially expressed genes in fresh protoplasts functioned in material synthesis and degradation, energy supply and abiotic stress responses (Figs. 1C and 6; Table S6). In addition, more than 20 transcription factors, 16 kinase activity-related genes, a number of phytohormone metabolism and signaling-related genes, and three methyltransferase-encoding genes also showed peak expression level in freshly isolated protoplasts. It is interesting that several homologues of stem cell-associated transcription factors, such as NAC2, CUC2, RD26 and WOX13, in vascular plants demonstrated maximal expression levels.

One hundred and nineteen genes in pattern B showed peak expression levels in t-24 cells, which represented a key stage in the reprogramming of P. patens protoplasts (Figs. 1D and 6; Table S5). Among the annotated genes, most were related to material and energy metabolism, including a few genes involved in amino acid metabolism and starch synthesis and genes encoding organic materials degradation (Fig. 4; Table S5). Less than 15 genes that encoded transcription factors were specifically expressed. In addition, five methylation-associated genes and two chromatin remodeling genes also exhibited significantly high expression levels.

Genes in pattern C were specifically expressed at 48 h, at which time the reprogrammed protoplasts re-entered the cell cycle (Figs. 1E and 6; Table S5). In contrast to fresh protoplasts and t-24 cells, genes involved in photosynthesis and structural constituents of the ribosome were the most abundant categories at this time-point, followed by energy metabolism and fatty acid biosynthesis genes. Some DEGs were related to phytohormone signaling transduction. As for regulatory factors, four transcription factors were maximally expressed, including a putative bZIP1, KAN2 and HRD.

Pattern D contained the genes preferentially expressed in cells cultured for 72 h, at which point stem cell reprogramming was complete and chloronemata regeneration had been initiated (Figs. 1F and 4; Table S5). Of the 196 genes preferentially expressed, the majority functioned in metabolism and protein synthesis (Table S5). A few genes that encoded putative transcription factors and phytohormone signal transduction-associated factors were also included in this pattern. In addition, a putative gene (root hair specific 8; RHS8) that involved in root hair development in vascular plants [37] was specifically expressed. According to previous reports, protonemal apical cells in mosses, like root hairs in flowering plants, are tip-growing cells, and some genes involved in plant cell tip growth were also functionally conserved between Physcomitrella and flowering plants [38], [39] Thus, we speculated that RHS8 might play an important role during protonema/chloronema development after stem cell reprogrammed in P. patens. In addition, a putative cyclin-dependent protein kinase, CAK1AT, which encoded a homologue of CDK-activating kinase 1AT and functioned in maintenance of root meristem identity in Arabidopsis [40], was maximally expressed.

It is well known that protein methylation and phytohormones play vital roles throughout all aspects of plant growth and development, including stress responses. In the present study, several genes that encode homologues of Arabidopsis methyltransferase showed maximal expression within 24 h of protoplast culture (i.e. Patterns C and D; Fig. 4, Table S5). More than 30 phytohormone-related genes, which included almost all phytohormone categories, showed maximal expression in fresh protoplasts. However, a few genes that associated with three main phytohormones, namely auxin, ABA and jasmonic acid, showed dynamic expression levels throughout the protoplast reprogramming process. From comparison of the gene number and category at each stage, we suggest that phytohormones might play important roles in the early stage of protoplast reprogramming into stem cells.

Discussion

Physcomitrella patens is an excellent system for exploration of the mechanisms of protoplast reprogramming into stem cells

At recent annual moss meetings (2007) and the 21^st International Conference on Arabidopsis Research in 2010, the P. patens apical stem cell system received considerable attention. In addition to the advantages of the system already outlined, cell divisions can be followed easily in time and space because nearly all gametophytic cells are in direct contact with the environment (i.e., there are no ‘hidden’ cell layers). In addition, the control of cell fate over time and space can be studied both in the protonemal cell lineage and in the more complex tissue of the gametophore. Orthologues for many genes are present in angiosperms and close paralogs are present in gymnosperms. The genome of P. patens encodes homologues of stem cell-related genes, such as three WOX paralogs and NAC transcription factors [16], [41]. However, how these genes function during development or protoplast regeneration is still unclear. The morphological and structural characteristics allow us to track the mechanisms both at the single-cell level and in living plants. Compared to other protoplast systems, P. patens protoplasts are readily isolated from young protonemal tissue [2]. A major advantage of the P. patens system is the ease with which cells can be reprogrammed into pluripotent stem cells. Another advantage of the P. patens ‘protoplast–stem cell system’ is that protoplast regeneration occurs at a high frequency, does not require supplementation of the medium with hormones, and results directly in the formation of chloronemal filaments or even a whole gametophyte and does not undergo a callus stage [24]. In these respects the process thus resembles spore germination [24]. Flow cytometry analysis of the cell cycle also showed high synchronization of the process after protoplast isolation.

Because of the higher degree of synchronization, the ease with which reprogramming to stem cells is induced, and the shorter times for regeneration, the P. patens protoplast–stem cell system is a more favorable system with which to explore molecular mechanisms for stem cell reprogramming than other protoplast systems, such as tobacco, cucumber and clover blossom [6], [9], [42].

Altered metabolism and stress responses are necessary for protoplast survival and reprogramming into stem cells

In the course of protoplast isolation from protonemata, a cell rapidly (within seconds) changes its form and is converted into a geometrically ideal sphere. The communication between cells is interrupted. In order to survive, protoplasts must first cope with the alteration of the surrounding environment and renew the cell configuration. We propose that a cascade of stress signal reactions is generated when a protoplast is completely separated from a cell wall. Such reactions might be similar to those accompanied by plasmolysis of plant cells caused by drought or salt. Changes in the expression levels of stress-response genes in the protoplasts might be partially associated with the process of protoplast separation per se. This process is likely to be stressful for a protonemal cell because the cell wall is enzymatically degraded and the protoplast is separated from its mechanically rigid, protective polysaccharide envelope into a hypertonic solution. The process of protoplast separation from a cell wall is a type of cell wounding. To survive, the protoplast must initiate the repair mechanism to confront the damage induced by cell wall removal, which represent alteration of metabolism of other materials and energy supply. Our results show that the significantly enriched pathways, which were concentrated in decreased photosynthesis, increased fatty acid, O-glycan, flavonoid, flavone and flavonol biosynthesis, and propanoate metabolism, contribute to the survival and rejuvenation of the protoplast.

Results from morphogenetic observation and FCM (Fig. 1) show that after culture for 48 h most protoplast-derived cells have been successfully reprogrammed into stem cells and re-enter the cell cycle. In order to satisfy the need for rapid division and growth, the cells show enhanced photosynthesis and protein synthesis as well as metabolism of other materials and energy supply.

After 72 h culture, the reprogrammed stem cells had divided once or twice, and chloronema development was initiated. The characteristics of metabolic processes were similar to those at 48 h. The increased expression levels of several embryogenesis-related genes indicated that protoplast-derived chloronema development showed similarities to embryo development.

Key candidate factors and regulatory mechanisms involved in protoplast reprogramming into stem cells in P. patens

To maintain the stem cell state, it is important that multiple regulation factors are coordinated, such as transcription factors, phytohormone synthesis and signaling transduction, and DNA and protein modification, except for alterations in metabolism.

1. Stem cell-related genes and their roles

For survival, the protoplast must also cope with induction of reprogramming into stem cells, therefore activation of genes associated with stem cell sustenance and meristem identity is essential. Our results showed that a number of homologues of meristem identity-related transcription factors, including NAC2, CUC2, RD26, WOX13 and BAM2, and one putative protein kinase exhibit maximal expression levels in freshly isolated protoplasts (Table S6). NAC2, CUC2 and RD26 encode transcriptional activators of the NAC gene family in Arabidopsis [43]. NAC2 is known to control age-dependent senescence and salt-promoted senescence [44]. RD26 is induced in response to desiccation, is localized to the nucleus and acts as a transcriptional activator in the ABA-mediated dehydration response [45]. The higher expression level of these genes may be correlated to the alteration of the culture medium, which contained a higher concentration of mannitol after removal of the cell wall. CUC2 expression in the leaf sinus region is required for serration and the extent of serration is modulated by mir164A-mediated repression of CUC2 [46], [47].

It was interesting that a homologue of WOX13, a member of the Wuschel-related homeobox (WOX) gene family in vascular plants, showed peak expression in freshly isolated protoplasts and a relatively high level of transcription was maintained after 24 h culture (Fig. 3, Tables S1 and S5). The WOX gene family belongs to the homeodomain-containing transcription factors, which are key regulators implicated in the determination of cell fate and cell differentiation in plants [48], [49]. WOX genes are specifically expressed in different plant organs and cell types [49]–[54]. WUS and WOX5 genes are two vital members that function in stem cell maintenance in a restricted region of the shoot apical meristem and root apical meristem, respectively [55]–[57]. WOX13 is among the most highly conserved WOX genes and affects root and flower development in Arabidopsis [16]. In P. patens, no homologues of WUS and WOX5 are known; of three WOX genes, two are homologues of WOX13 and one is a homologue of WOX14 [16]. Therefore, we propose that WOX genes in P. patens possess broader functions. Our results indicate that the homologue of WOX13 may play a key role in stem cell identity during protoplast reprogramming into stem cells. Down-regulation of PpWOX13a in cells after culture for 24 h may be associated with the smaller apical stem number that led to cell division or the other two WOX homologues may have complementary functions. This finding provides an important clue for further functional study of WOX genes in P. patens and might also contribute to elucidation of the mechanisms involved in protoplast reprogramming and maintenance of a single apical stem cell in the protonema and leafy shoot apex.

Pp1s352_22V6 encodes a homologue of an Arabidopsis CLAVATA1-related receptor kinase-like protein, Barely Any Meristem 2 (BAM2), which is a member of the leucine-rich receptor-like protein kinase family that is required for both shoot and flower meristem function in Arabidopsis [58]. This gene shows a broad expression pattern and is involved in multiple developmental processes, such as vascular strand development in the leaf, control of leaf shape, size and symmetry, male gametophyte (especially anther) development, and ovule specification and function [55], [56]. In addition, the BAM2 expression pattern supports both an early role in promoting somatic cell fates and a subsequent function in pollen mother cells.

In Arabidopsis, stem cell identity maintenance is dependent on the WUS-CLV feedback regulation loop in the shoot apical meristem or on WOX5-SCR/SHR feedback regulation in root apical meristems, which are both dependent on stem niche [59]. P. patens lacks a stem cell niche and only one apical stem cell is present in the protonema tip and leafy shoot apex [16]. In the present study, the higher expression levels of homologues of WOX13 and BAM2 indicate that the two genes may play pivotal roles in protoplast transformation into stem cells. We further hypothesize that protoplast reprogramming into stem cells and maintenance of the single stem cell identity may be partially determined by interaction of WOX and BAM proteins. However, further experiments are needed to test this hypothesis.

2. Epigenetic/methylation modification of protein during stem cell reprogramming

Genetic evidence indicates that, similar to animals, stem cells in plants possess a specialized chromatin structure. This may reflect their capacity for a variety of gene-expression programs, as well as their ability to divide repeatedly without either differentiation or senescence. In recent years, epigenetic mechanisms that control chromatin structure and function, including DNA methylation and histone modification, have emerged as key factors in the regulation of cell growth and differentiation and, thereby, the nuclear reprogramming necessary for dedifferentiation [60]. Protein methylation is one type of protein post-translational modification and has been most studied for histones, which can act epigenetically to repress or activate gene expression [61], [62].

Several genes that encode homologues of methyltransferase in Arabidopsis, including Pp1s233_104V6 (protein arginine methyltransferase 4B, PRMT4B), Pp1s271_68V6 (cobalamin-independent synthase family protein, ATMS1), Pp1s31_85V6 (protein-l-isoaspartate methyltransferase 1, PIMT1), Pp1s35_262V6 (protein-l-isoaspartate methyltransferase 2, PIMT2) and Pp1s8_46V6 (homocysteine methyltransferase 2, HMT2), are preferentially expressed within 24 h of protoplast isolation, a stage of stem cell reformation. These results indicate that protein methylation may be an additional important epigenetic mechanism for establishment and/or maintenance of the undifferentiated state during protoplast reprogramming into stem cells.

3. Endogenous phytohormones and phytohormone-responsive genes

How the regulators described above actually control the behavior of protoplast reprogramming into stem cells in P. patens is still largely unknown. Phytohormones are the most likely candidates as regulators of developmental switches, and it has been proposed that hormones play a central role in mediating the signal transduction cascade that leads to the reprogramming of gene expression. Our results indicate that the requirement for a wide variety of endogenous phytohormones during protoplast reprogramming into stem cells is largely determined by the stage of the cultures.

Cytokinins have many roles in plant development, one of which is to stimulate cell division in the shoot in vascular plants [63]. Cytokinins induce bud formation in mosses in a concentration-dependent manner [64]. Cytokinin biosynthesis and signal transduction might be necessary for the activation of protoplast reprogramming processes, but not for subsequent stem cell reprogramming processes. For example, isopentenyl transferase 9 (IPT9), which catalyzes the first step in isopentenyl-type cytokinin (a major cytokinin in moss) de novo biosynthesis [65], shows maximal expression in freshly isolated protoplasts.

Moreover, several ethylene-response genes show preferential expression during both stages of initiation of reprogramming and chloronemal cell division and differentiation after stem cell formation. However, some genes related to the synthesis and signaling of auxin, ABA and jasmonic acid are also altered the transcriptional levels throughout the reprogramming process after removal of the cell wall. Thus, we speculate that a dynamic balance in the interactions among phytohormones, especially of auxin, ABA and jasmonic acid may be pivotal in distinct stages of protoplast reprogramming and switching cell fate during protoplast reprogramming into stem cells, through coordinated interactions with many metabolic pathways (for example, photosynthesis, cell respiration, and protein synthesis and degradation). Nevertheless, further experiments are necessary to verify this conclusion.

Conclusion

In this study, we analyzed the transcriptome during P. patens protoplast reprogramming into stem cells. The transcript levels of 4827 genes were significantly increased or reduced, of which the majority changed only during a specific stage. Our results provide an extensive catalogue of regulatory factors and related genes involved in protoplast survival, reprogramming, and cell division during protoplast reprogramming into stem cells in P. patens. Potential applications of these data include identification of candidate genes and as targets for reverse genetic studies of stem cell maintenance and its evolution, and as tools for exploration of the molecular mechanisms underlying stem cell reprogramming in plants.

Materials and Methods

Preparation of protoplasts and culture conditions

Protonemal tissue of the Gransden 2004 wild-type strain of Physcomitrella patens subsp. patens was subcultured at weekly intervals on solid BCDAG medium overlaid on cellophane and containing 5 mM ammonium tartrate [66]. Protoplasts were isolated from six-day-old protonemal tissues following the modified protocol of Grimsley et al. [67]. Ten Petri dishes contained preplasmolyzed protonemal filaments pretreated with 0.48 M mannitol supplemented with 10 ml of 1% solution of Driselase in 0.48 M mannitol (Sigma-Aldrich, USA). After agitated incubation for 30 min in the dark at 25°C, the culture was successively passed through sieves of pore size 100 µm and 50 µm. The filtrate was incubated without agitation under the same conditions for 15 min. Then the protoplast suspension was precipitated by centrifugation for 5 min at 600 rpm and the pellet was washed twice with 0.48 M mannitol by centrifugation at the same rate. The number of protoplasts was determined with a hemocytometer. A portion of the freshly isolated protoplasts were used immediately for FCM analysis and RNA isolation. The remainder of the protoplasts were cultured in BCDAG liquid medium without hormones, but supplemented with 5 mM ammonium tartrate and 0.48 M mannitol, and cultured in the dark for 24 h at 25±1°C. Then the cultures were incubated under a 16/8 h (light/dark) photoperiod regime with light intensity at the surface of the vessels of 55 µmol m⁻² s⁻¹ provided by Philips TLD25 fluorescent lamps.

Flow cytometry analysis of DNA content

To analyze the DNA content and cell cycle alteration during protoplast culture, FCM analysis was performed according to the method of Ulrich and Ulrich [68]. Suspensions of intact nuclei were prepared from six-day-old protonema (control), freshly isolated protoplasts (t-0) and cultures incubated for 24 h (t-24), 48 h (t-48) and 72 h (t-72) after protoplast isolation by chopping the material with a razor blade in a glass Petri dish that contained 500 µl specific buffer (45 mM MgCl₂•6H₂O, 30 mM sodium citrate, 20 mM MOPS, 0.2 mg/ml Triton X-100) [69]. Nuclei from mature Arabidopsis leaves and rice root tips were prepared for determination of the cell cycle stage of the protoplast and its derived cultures. The solution was filtered twice through a sieve of 30 µm pore size. After filtration, a final concentration of 2 µg/ml DAPI (Sigma-Aldrich, USA) was added and incubated for at least 5 min in the dark to stain the nuclei. FCM analysis was conducted with a Becton Dickinson FACSort system (Becton Dickinson, Mountain View, CA, USA). Data were analyzed using Quanta SC software.

Preparation of digital expression libraries

To achieve comprehensive gene expression profiling during protoplast reprogramming in P. patens, samples from t-0 protoplasts and t-24, t-48 and t-72 cultured cells were pooled for RNA isolation and library construction. Total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions. RNase-free DNase (RQ1; Promega Corporation, USA) was used to remove genomic DNA contamination. RNA was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, USA). RNA quality was tested at an absorbance of A₂₆₀/A₂₈₀>2.0 with an Agilent 2100 Bioanalyzer using the RNA 6000 LabChip Kit (Agilent Technologies, Palo Alto, CA). Only RNA with RNA integrity values greater than 8.0 were used for digital gene expression tag profiling (DGEP).

Tag libraries for each mRNA 3′ terminus were prepared with the Digital Gene Expression Tag Profiling Kit version 2.1B (Illumina, San Diego, CA, USA) in accordance with the manufacturer's protocol. Six micrograms of total RNA was extracted and mRNA was purified using biotin-Oligo (dT) magnetic bead adsorption. Messenger RNA enriched on the beads by oligo (dT) was converted into double-stranded cDNA through reverse transcription. While on the beads, double-strand cDNA was ligated with NlaIII endonuclease to produce a bead-bound cDNA fragment that comprised the sequence from the 3′-most CATG to the poly (A)-tail. Magnetic bead precipitation was used to purify the cDNA fragments with 3′ ends and add the Illumina adapter 1 to their new 5′ ends. The junction of the Illumina adapter 1 and the CATG site was recognized by MmeI, which is a type I endonuclease with separate recognition sites and digestion sites. MmeI cuts 17 bp downstream of the CATG site to produce 17 bp cDNA sequence tags with adapter 1. After removing 3′ fragments by magnetic bead precipitation, the Illumina adapter 2 was ligated to the 3′ end of the cDNA tags, which generated tags with different adapters at both ends to form a tag library.

Digital gene expression tag profiling

Linear PCR amplification with 15 cycles was performed with primers complementary to the adapter sequences to enrich the samples for the desired fragments. The 85 bp fragments were purified by 6% TBE PAGE gel electrophoresis. These fragments were then digested and the single-chain molecules fixed onto a Solexa Sequencing Chip (flowcell) for Solexa sequencing by an Illumina Cluster Station and Illumina Genome Analyzer II System (version 1.0). Each molecule grows into a single-molecule cluster sequencing template by in situ amplification. Four types of nucleotides, which were labeled by different colors, were added and sequencing by synthesis was performed. Each tunnel generated millions of raw reads with a sequence length of 35 bp.

Sequence annotation and identification of differentially expressed genes

Raw reads were transformed into clean tags by trimming adaptor-only tags, filtering low-quality tags (those that contained ambiguous bases), tags that were too long (>21 nt) or too short (<21 nt), and tags that consisted of a single copy (probably the result of sequencing error). The remaining tags were designated distinct clean tags. The genomic sequences corresponding to the tags were retrieved from the Physcomitrella genome sequence assembly (ftp://ftp.jgi-psf.org/pub/JGI_data/phytozome/v7.0/Ppatens/annotation/). All distinct clean tags were aligned to the reference P. patens database with Bowtie 0.12.7 (http://bowtie-bio.sourceforge.net/) and tolerances were set to allow no more than 1 nt mismatch. Tags that mapped multigenes were removed. The number of distinct clean tags for each gene was calculated and then normalized using the TMM method in edgeR (http://www.bioconductor.org/packages/2.3/bioc/html/). Saturation analysis was performed to check whether the number of detected genes continued to increase as the total tag number increases. Gene models for all of the up- and down-regulated tags in this study can be obtained from the Physcomitrella Genome Browser using the gene ID number assigned to each tag feature.

A rigorous algorithm to identify DEGs between two samples was developed for significance testing [70]. The P value corresponds to the differential gene expression test. The FDR was applied to determine the threshold P value in multiple tests and analyses [71]. The number of tags mapped to a reference gene was considered as the expression abundance of the gene. Expression abundances of a gene from two samples were compared to determine differences in expression. A gene was considered to show a statistically significant change in expression between samples when the expression difference was more than four-fold the cutoff value (|log₂Ratio|≥2) with P value≤0.01 and FDR<0.01.

Quantitative real-time PCR assay

Samples were prepared using the method described above and total RNA was isolated from three biological repeats of t-0 protoplasts and t-24, t-48 and t-72 cultured cells. DNA contamination was removed from each RNA sample using RQ1 RNase-free DNase (RQ1; Promega Corporation, USA). cDNA synthesis was performed with SuperScript™ III reverse transcriptase (Invitrogen) and 5 µg total RNA for each sample using oligo (dT18) primers. For real-time PCR, gene-specific primers were designed using Primer 5 software (Primer 5, Applied Biosystems, USA), assessed by Oligo 6 and synthesized by Shenggong Cooperation (Shenggong, Shanghai, China). Primers used in the real-time PCR assay are listed in Table S2. Real-time PCR was performed using a Corbett Research Rotor-Gene 3000 under the following conditions: 94°C for 5 min (1 cycle); 94°C for 20 s, 50–60°C for 20 s, and 72°C for 20 s (50 cycles). Transcript abundances were identified using the SYBR Green PCR Master Mix (TaKaRa). Each reaction contained 1× mix buffer, 0.25 µM each primer, and about 2 ng cDNA in a final volume of 20 µl. qRT-PCR for each gene was performed on three biological replicates with three technical replicates per biological replicate. Melting curves were performed on the product to test if only a single product was amplified without primer–dimers and other bands. The resulting products with all primer combinations were visualized on a 2% agarose gel to confirm the generation of a single product of the correct size. Pp1s97_279V6 expression was used as an internal control to normalize all data. Relative quantitative analysis was performed using comparative quantitation with Rotor-Gene Real-Time Analysis Software 6.1. Significant variation from the internal control was determined using Student's t-test where p≤0.05 was considered to be differentially expressed.

Gene Ontology functional and pathway enrichment analysis of differentially expressed genes

Gene Ontology (GO, http://www.geneontology.org/) is an international standardized gene functional classification system that offers a dynamic-updated controlled vocabulary and a strictly defined concept to comprehensively describe properties of genes and their products in any organism. GO covers three domains: cellular component, molecular function, and biological process. GO functional enrichment analysis was used to identify significant enriched GO terms in DEGs. Statistical significance of GO terms was calculated with the following formula [36]:where N is the number of all genes with GO annotation, n is the number of DEGs in N, M is the number of all genes that are annotated to certain GO terms, and m is the number of DEGs in M. GO terms with a P value<0.05 were significantly enriched in DEGs. The top 15 significantly enriched GO terms in each DEG set were considered for further analysis on the basis of the enrichment factor.

Pathway enrichment analysis was used to identify significantly enriched metabolic pathways or signal transduction pathways in DEGs compared with the whole genome background. Pathway analysis of DEGs was performed using the KEGG database and Fisher's exact test algorithm to calculate the statistical significance of each pathway [36]. The formula used was the same as that for GO analysis. In this instance, N is the number of all genes with KEGG annotation, n is the number of DEGs in N, M is the number of all genes annotated to specific pathways, and m is number of DEGs in M. Since a large number of pathways were involved, we implemented FDR correction to control the overall Type I error rate of multiple hypotheses. The Q value was calculated to determine the threshold P value in multiple tests and analyses [72]. Pathways with a Q value<0.05 were significantly enriched in DEGs. In addition, a cluster analysis was performed using the statistical significance pathway data of the three sample groups with Cluster 3.0 (http://smd.stanford.edu/resources/restech.shtml). The x-axis was the pathway's p value, and the y-axis was the different sample groups [73].

Identification of key candidate factors during protoplast reprogramming

To identify the key candidate factors during protoplast reprogramming, the comparison of two consecutive time-points was transformed into a comparison using the t-0 time-point as a common reference point. We obtained the union of the three DEG groups, and performed a K-means clustering analysis with MEV (http://www.tm4.org/mev/) using the normalized gene expression level of each DEG.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Distribution of distinct clean tags in each sample.

https://doi.org/10.1371/journal.pone.0035961.s001

(TIF)

Figure S2.

Saturation evaluation of expression of unique tags in each sample.

https://doi.org/10.1371/journal.pone.0035961.s002

(TIF)

Figure S3.

Smear plots from the edgeR-based analysis of gene expression. Genes are plotted based on their log-fold change of transcript abundance between two compared samples on the y-axis and log concentration on the x-axis for raw tag libraries separately. Differentially expressed genes are shown in red.

https://doi.org/10.1371/journal.pone.0035961.s003

(TIF)

Figure S4.

Top 15 significantly enriched GO terms in each sample.

https://doi.org/10.1371/journal.pone.0035961.s004

(TIF)

Figure S5.

Top 15 significantly enriched pathways in each sample.

https://doi.org/10.1371/journal.pone.0035961.s005

(TIF)

Table S1.

List of differentially expressed genes during protoplast reprogramming into stem cells.

https://doi.org/10.1371/journal.pone.0035961.s006

(XLS)

Table S2.

Primer sequences of the selected and internal genes used for real-time PCR analysis.

https://doi.org/10.1371/journal.pone.0035961.s007

(DOC)

Table S3.

List of the R² value between DGEP and qRT-PCR.

https://doi.org/10.1371/journal.pone.0035961.s008

(XLS)

Table S4.

Significantly enriched GO terms of differentially expressed genes.

https://doi.org/10.1371/journal.pone.0035961.s009

(XLS)

Table S5.

Significantly enriched pathways.

https://doi.org/10.1371/journal.pone.0035961.s010

(XLS)

Table S6.

Complete gene list of each K-means cluster during protoplast reprogramming in P. patens .

https://doi.org/10.1371/journal.pone.0035961.s011

(XLS)

Acknowledgments

We thank Fang Bao at the College of Life Sciences, Capital Normal University, Beijng, China, for her aid on the qRT-PCR validation of DGEP.

Author Contributions

Conceived and designed the experiments: YH. Performed the experiments: LX LZ GY. Analyzed the data: LX HZ. Contributed reagents/materials/analysis tools: YH LX. Wrote the paper: LX.

Citation: Xiao L, Zhang L, Yang G, Zhu H, He Y (2012) Transcriptome of Protoplasts Reprogrammed into Stem Cells in Physcomitrella patens. PLoS ONE 7(4): e35961. https://doi.org/10.1371/journal.pone.0035961

References

1. Gilbert SF (2006) Developmental Biology. Sunderland, MA: Sinauer Associates. SF Gilbert2006Developmental BiologySunderland, MASinauer Associates

2. Slack JM (2008) Origin of stem cells in organogenesis. Science 322: 1498–1501.JM Slack2008Origin of stem cells in organogenesis.Science32214981501

3. McKay R (2000) Stem cells – hype and houpe. Nature 406: 361–364.R. McKay2000Stem cells – hype and houpe.Nature406361364

4. Byrne ME, Kidner CA, Martienssen RA (2003) Plant stem cells: divergent pathways and common themes in shoots and roots. Current Opinion in Genetics & Development 13: 551–557.ME ByrneCA KidnerRA Martienssen2003Plant stem cells: divergent pathways and common themes in shoots and roots.Current Opinion in Genetics & Development13551557

5. Grafi G, Avivi Y (2004) Stem cells: a lesson from dedifferentiation. TRENDS in Biotechnology 22(8): 388–389.G. GrafiY. Avivi2004Stem cells: a lesson from dedifferentiation.TRENDS in Biotechnology228388389

6. Ondřej V, Kitner M, Doležalov I, Nádvorník P, Navrátilová B, et al. (2009) Chromatin Structural Rearrangement during Dedifferentiation of Protoplasts of Cucumis sativus L. Mol Cells 27: 443–447.V. OndřejM. KitnerI. DoležalovP. NádvorníkB. Navrátilová2009Chromatin Structural Rearrangement during Dedifferentiation of Protoplasts of Cucumis sativus L.Mol Cells27443447

7. Masip M, Veiga A, Izpisua JC, Simon C (2010) Reprogramming with defined factors: From induced pluripotency to induced transdifferentiation. Mol Hum Reprod 16: 856–868.M. MasipA. VeigaJC IzpisuaC. Simon2010Reprogramming with defined factors: From induced pluripotency to induced transdifferentiation.Mol Hum Reprod16856868

8. Takebe I, Labib G, Melchers G (1971) Regeneration of whole plants from isolated mesophyll protoplasts of tobacco. Naturwissenschaften 58: 318–320.I. TakebeG. LabibG. Melchers1971Regeneration of whole plants from isolated mesophyll protoplasts of tobacco.Naturwissenschaften58318320

9. Zhao J, Morozova N, Williams L, Libs L, Avivi Y, et al. (2001) Two phases of chromatin decondensation during dedifferentiation of plant cells: distinction between competence for cell fate switch and a commitment for S phase. J Biol Chem 276: 22772–22778.J. ZhaoN. MorozovaL. WilliamsL. LibsY. Avivi2001Two phases of chromatin decondensation during dedifferentiation of plant cells: distinction between competence for cell fate switch and a commitment for S phase.J Biol Chem2762277222778

10. Tessadori F, Chupeau MC, Chupeau Y, Knip M, Germann S, et al. (2007) Large-scale dissociation and sequential reassembly of pericentric heterochromatin in dedifferentiated Arabidopsis cells. J Cell Sci 120: 1200–1208.F. TessadoriMC ChupeauY. ChupeauM. KnipS. Germann2007Large-scale dissociation and sequential reassembly of pericentric heterochromatin in dedifferentiated Arabidopsis cells.J Cell Sci12012001208

11. Cocking EC (1960) A method for the isolation of plant protoplasts and vacuoles. Nature 187: 927–929.EC Cocking1960A method for the isolation of plant protoplasts and vacuoles.Nature187927929

12. Debeaujon I, Branchard M (1992) Induction of somatic embryogenesis and caulogenesis from cotyledon and leaf protoplast-derived colonies of melon (Cucumis melo L.). Plant Cell Rep 12: 37–40.I. DebeaujonM. Branchard1992Induction of somatic embryogenesis and caulogenesis from cotyledon and leaf protoplast-derived colonies of melon (Cucumis melo L.).Plant Cell Rep123740

13. Takebe I, Otsuki Y (1969) Infection of tobacco mesophyll protoplasts by tobacco mosaic virus. Proc Natl Acad Sci US 64: 843–848.I. TakebeY. Otsuki1969Infection of tobacco mesophyll protoplasts by tobacco mosaic virus.Proc Natl Acad Sci US64843848

14. van Driel R, Fransz P (2004) Nuclear architecture and genome functioning in plants and animals: what can we learn from both? Exp Cell Res 296: 86–90.R. van DrielP. Fransz2004Nuclear architecture and genome functioning in plants and animals: what can we learn from both?Exp Cell Res2968690

15. Tessadori F, van Driel R, Fransz P (2004) Cytogenetics as a tool to study gene regulation. Trends Plant Sci 9: 147–153.F. TessadoriR. van DrielP. Fransz2004Cytogenetics as a tool to study gene regulation.Trends Plant Sci9147153

16. Deveaux Y, Toffano-Nioche C, Claisse G, Thareau V, Morin H, et al. (2008) Genes of the most conserved WOX clade in plants affect root and flower development in Arabidopsis. BMC Evolutionary Biology 8(291): 1–19.Y. DeveauxC. Toffano-NiocheG. ClaisseV. ThareauH. Morin2008Genes of the most conserved WOX clade in plants affect root and flower development in Arabidopsis.BMC Evolutionary Biology8291119

17. Harrison CJ, Roeder AH, Meyerowitz EM, Langdale JA (2009) Local cues and asymmetric cell divisions underpin body plan transitions in the moss Physcomitrella patens. Curr Biol 19: 461–471.CJ HarrisonAH RoederEM MeyerowitzJA Langdale2009Local cues and asymmetric cell divisions underpin body plan transitions in the moss Physcomitrella patens.Curr Biol19461471

18. Rensing , et al. (2008) The Physcomitrella Genome reveals evolutionary insights into the conquest of land by plants. Science 319(5859): 64–69.Rensing2008The Physcomitrella Genome reveals evolutionary insights into the conquest of land by plants.Science31958596469

19. Reski R (1998) Development, genetics and molecular biology of mosses. Bot Acta 111: 1–15.R. Reski1998Development, genetics and molecular biology of mosses.Bot Acta111115

20. Frank W, Decker EL, Reski R (2005) Molecular tools to study Physcomitrella patens. Plant Biol 7: 220–227.W. FrankEL DeckerR. Reski2005Molecular tools to study Physcomitrella patens.Plant Biol7220227

21. Wang XQ, Yang PF, Gao Q, Liu X, Kuang TY, et al. (2008) Proteomic analysis of the response to high-salinity stress in Physcomitrella patens. Planta 228: 167–177.XQ WangPF YangQ. GaoX. LiuTY Kuang2008Proteomic analysis of the response to high-salinity stress in Physcomitrella patens.Planta228167177

22. Cumming AC, Cho SH, Kamisugi Y, Graham H, Quatrano RS (2007) Microarray analysis of transcriptional responses to abscisic acid and osmotic, salt, and drought stress in the moss, Physcomitrella patens. New Phytologist 176: 275–287.AC CummingSH ChoY. KamisugiH. GrahamRS Quatrano2007Microarray analysis of transcriptional responses to abscisic acid and osmotic, salt, and drought stress in the moss, Physcomitrella patens.New Phytologist176275287

23. Wang XQ, Yang PF, Liu Z, Liu WZ, Hu Y, et al. (2009) Exploring the Mechanism of Physcomitrella patens Desiccation Tolerance through a Proteomic Strategy. Plant Physiology 149: 1739–1750.XQ WangPF YangZ. LiuWZ LiuY. Hu2009Exploring the Mechanism of Physcomitrella patens Desiccation Tolerance through a Proteomic Strategy.Plant Physiology14917391750

24. Cove DJ (2000) The Moss, Physcomitrella patens. J Plant Growth Regul 19: 275–283.DJ Cove2000The Moss, Physcomitrella patens.J Plant Growth Regul19275283

25. Ishikawa M, Murata T, Sato Y, Nishiyama T, Hiwatashi Y, et al. (2011) Physcomitrella Cyclin-Dependent Kinase A Links Cell Cycle Reactivation to Other Cellular Changes during Reprogramming of Leaf Cells. The Plant Cell. M. IshikawaT. MurataY. SatoT. NishiyamaY. Hiwatashi2011Physcomitrella Cyclin-Dependent Kinase A Links Cell Cycle Reactivation to Other Cellular Changes during Reprogramming of Leaf Cells.The Plant Cell

26. den Boer BG, Murray JA (2000) Triggering the cell cycle in plants. Trends Cell Biology 10: 245–250.BG den BoerJA Murray2000Triggering the cell cycle in plants.Trends Cell Biology10245250

27. Schween G, Gorr G, Hohe A, Reski R (2003) Unique tissue-specific delection of DNA synthesis in vivo. Proc Natl Acad Sci U S A 105: 2415–2420.G. SchweenG. GorrA. HoheR. Reski2003Unique tissue-specific delection of DNA synthesis in vivo.Proc Natl Acad Sci U S A10524152420

28. Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, et al. (2010) Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant. Genomics 96(6): 369–376.QQ WangF. LiuXS ChenXJ MaHQ Zeng2010Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant.Genomics966369376

29. Xiang LX, He D, Dong WR, et al. (2010) Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics 11: 472.LX XiangD. HeWR Dong2010Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish.BMC Genomics11472DOI:10.1186/1471-2164-11-472. DOI:10.1186/1471-2164-11-472.

30. Xiao S, Mo D, Wang Q, Jia J, Qin L, et al. (2010) Aberrant host immune response induced by highly virulent PRRSV identified by digital gene expression tag profiling. BMC Genomics 11: 544.S. XiaoD. MoQ. WangJ. JiaL. Qin2010Aberrant host immune response induced by highly virulent PRRSV identified by digital gene expression tag profiling.BMC Genomics11544DOI:10.1186/1471-2164-11-544. DOI:10.1186/1471-2164-11-544.

31. Langmead B, Trapnell C, Pop M, Salzberg SL (2010) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.B. LangmeadC. TrapnellM. PopSL Salzberg2010Ultrafast and memory-efficient alignment of short DNA sequences to the human genome.Genome Biol10R25

32. Robinson MD, Smyth GK (2007) Moderated statistical tests for assessing differences in tag abundance. Bioinformatics 23: 2881–2887.MD RobinsonGK Smyth2007Moderated statistical tests for assessing differences in tag abundance.Bioinformatics2328812887

33. Robinson MD, Smyth GK (2008) Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics 9: 321–332.MD RobinsonGK Smyth2008Small-sample estimation of negative binomial dispersion, with applications to SAGE data.Biostatistics9321332

34. Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.MD RobinsonDJ McCarthyGK Smyth2010edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.Bioinformatics26139140

35. Babbitt CC, Fedrigo O, Pfefferle AD, Boyle AP, Horvath JE, et al. (2010) Both non-coding and protein-coding RNAs contribute to gene expression evolution in the primate brain. Genome Biol and Evol 18(2): 67–79.CC BabbittO. FedrigoAD PfefferleAP BoyleJE Horvath2010Both non-coding and protein-coding RNAs contribute to gene expression evolution in the primate brain.Genome Biol and Evol1826779

36. Trajkovski I, Lavrac N, Tolar J (2008) SEGS: search for enriched gene sets in microarray data. J Biomed Inform 41(4): 588–601.I. TrajkovskiN. LavracJ. Tolar2008SEGS: search for enriched gene sets in microarray data.J Biomed Inform414588601

37. Won S, Lee Y, Lee H, Heo Y, Cho M, et al. (2009) cis-Element- and Transcriptome-Based Screening of Root Hair-Specific Genes and Their Functional Characterization in Arabidopsis. Plant Physiology 150: 1459–1473.S. WonY. LeeH. LeeY. HeoM. Cho2009cis-Element- and Transcriptome-Based Screening of Root Hair-Specific Genes and Their Functional Characterization in Arabidopsis.Plant Physiology15014591473

38. Menand B, Yi K, Jouannic S, Hoffmann L, Ryan E, et al. (2007) An ancient mechanism controls the development of cells with a rooting function in land plants. Science 316(5830): 1477–80.B. MenandK. YiS. JouannicL. HoffmannE. Ryan2007An ancient mechanism controls the development of cells with a rooting function in land plants.Science3165830147780

39. Jang G, Yi K, Pires ND, Menand B, Dolan L (2011) RSL genes are sufficient for rhizoid system development in early diverging land plants. Development 138(11): 2273–81.G. JangK. YiND PiresB. MenandL. Dolan2011RSL genes are sufficient for rhizoid system development in early diverging land plants.Development13811227381

40. Umeda M, Shimotohno A, Yamaguchi M (2005) Control of Cell Division and Transcription by Cyclin-dependent Kinaseactivating Kinases in Plants. Plant Cell Physiol 46(9): 1437–1442.M. UmedaA. ShimotohnoM. Yamaguchi2005Control of Cell Division and Transcription by Cyclin-dependent Kinaseactivating Kinases in Plants.Plant Cell Physiol46914371442

41. Richardt S, Lang D, Frank W, Reski R, Rensing SA (2007) PlanTAPDB: A phylogeny-based resource of plant transcription associated proteins. Plant Physiology 143: 1452–1466.S. RichardtD. LangW. FrankR. ReskiSA Rensing2007PlanTAPDB: A phylogeny-based resource of plant transcription associated proteins.Plant Physiology14314521466

42. Fehér A, Pasternak T, Ötvös K, Miskolczi P, Dudits D (2002) Induction of embryogenic competence in somatic plant cells: a review. Biologia 57: 5–12.A. FehérT. PasternakK. ÖtvösP. MiskolcziD. Dudits2002Induction of embryogenic competence in somatic plant cells: a review.Biologia57512

43. Ooka H, Satoh K, Doi K, Nagata T, Otomo Y, et al. (2003) Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Research 10(6): 239–247.H. OokaK. SatohK. DoiT. NagataY. Otomo2003Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana.DNA Research106239247

44. Balazadeh S, Siddiqui H, Allu AD, Matallana-Ramirez LP, Caldana C, Mehrnia M, et al. (2010) A gene regulatory network controlled by the NAC transcription factor ANAC092/AtNAC2/ORE1 during salt-promoted senescence. The Plant Journal 62: 250–264.S. BalazadehH. SiddiquiAD AlluLP Matallana-RamirezC. CaldanaM. Mehrnia2010A gene regulatory network controlled by the NAC transcription factor ANAC092/AtNAC2/ORE1 during salt-promoted senescence.The Plant Journal62250264

45. Bray EA (2004) Genes commonly regulated by water-deficit stress in Arabidopsis thaliana. Journal of Experimental Botany 55(407): 2331–2341.EA Bray2004Genes commonly regulated by water-deficit stress in Arabidopsis thaliana.Journal of Experimental Botany5540723312341

46. Aida M, Ishida T, Tasaka M (1999) Shoot apical meristem and cotyledon formation during Arabidopsis embryogenesis: interaction among the CUP-SHAPED COTYLEDON and SHOOT MERISTEMLESS genes. Development 126: 1563–1570.M. AidaT. IshidaM. Tasaka1999Shoot apical meristem and cotyledon formation during Arabidopsis embryogenesis: interaction among the CUP-SHAPED COTYLEDON and SHOOT MERISTEMLESS genes.Development12615631570

47. Nikovics K, Blein T, Peaucelle A, Ishida T, Morin H, et al. (2006) The Balance between the MIR164A and CUC2 Genes Controls Leaf Margin Serration in Arabidopsis. The Plant Cell 18: 2929–2945.K. NikovicsT. BleinA. PeaucelleT. IshidaH. Morin2006The Balance between the MIR164A and CUC2 Genes Controls Leaf Margin Serration in Arabidopsis.The Plant Cell1829292945

48. Haecker A, Gross-Hardt R, Geiges B, Sarkar A, Breuninger H, et al. (2004) Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arabidopsis thaliana. Development 131(3): 657–668.A. HaeckerR. Gross-HardtB. GeigesA. SarkarH. Breuninger2004Expression dynamics of WOX genes mark cell fate decisions during early embryonic patterning in Arabidopsis thaliana.Development1313657668

49. Nardmann J, Werr W (2006) The Shoot Stem Cell Niche in Angiosperms: Expression Patterns of WUS Orthologues in Rice and Maize Imply Major Modifications in the Course of Mono- and Dicot Evolution. Mol Biol Evol 23(12): 2492–2504.J. NardmannW. Werr2006The Shoot Stem Cell Niche in Angiosperms: Expression Patterns of WUS Orthologues in Rice and Maize Imply Major Modifications in the Course of Mono- and Dicot Evolution.Mol Biol Evol231224922504

50. Gallois JL, Nora FR, Mizukami Y, Sablowski R (2004) WUSCHEL induces shoot stem cell activity and developmental plasticity in the root meristem. Genes Dev 18(4): 375–380.JL GalloisFR NoraY. MizukamiR. Sablowski2004WUSCHEL induces shoot stem cell activity and developmental plasticity in the root meristem.Genes Dev184375380

51. Park SO, Zheng Z, Oppenheimer DG, Hauser BA (2005) The PRETTY FEW SEEDS2 gene encodes an Arabidopsis homeodomain protein that regulates ovule development. Development 132(4): 841–849.SO ParkZ. ZhengDG OppenheimerBA Hauser2005The PRETTY FEW SEEDS2 gene encodes an Arabidopsis homeodomain protein that regulates ovule development.Development1324841849

52. Kamiya N, Nagasaki H, Morikami A, Sato Y, Matsuoka M (2003) Isolation and characterization of a rice WUSCHEL-type homeobox gene that is specifically expressed in the central cells of a quiescent center in the root apical meristem. Plant J 35(4): 429–441.N. KamiyaH. NagasakiA. MorikamiY. SatoM. Matsuoka2003Isolation and characterization of a rice WUSCHEL-type homeobox gene that is specifically expressed in the central cells of a quiescent center in the root apical meristem.Plant J354429441

53. Wurschum T, Gross-Hardt R, Laux T (2006) APETALA2 regulates the stem cell niche in the Arabidopsis shoot meristem. Plant Cell 18(2): 295–307.T. WurschumR. Gross-HardtT. Laux2006APETALA2 regulates the stem cell niche in the Arabidopsis shoot meristem.Plant Cell182295307

54. Zuo J, Niu QW, Frugis G, Chua NH (2002) The WUSCHEL gene promotes vegetative-to-embryonic transition in Arabidopsis. Plant J 30(3): 49–359.J. ZuoQW NiuG. FrugisNH Chua2002The WUSCHEL gene promotes vegetative-to-embryonic transition in Arabidopsis.Plant J30349359

55. Mayer KF, Schoof H, Haecker A, Lenhard M, Jurgens G, et al. (1998) Role of WUSCHEL in regulating stem cell fate in the Arabidopsis shoot meristem. Cell 95(6): 805–815.KF MayerH. SchoofA. HaeckerM. LenhardG. Jurgens1998Role of WUSCHEL in regulating stem cell fate in the Arabidopsis shoot meristem.Cell956805815

56. Sarkar AK, Luijten M, Miyashima S, Lenhard M, Hashimoto T, et al. (2007) Conserved factors regulate signalling in Arabidopsis thaliana shoot and root stem cell organizers. Nature 446(7137): 811–814.AK SarkarM. LuijtenS. MiyashimaM. LenhardT. Hashimoto2007Conserved factors regulate signalling in Arabidopsis thaliana shoot and root stem cell organizers.Nature4467137811814

57. DeYoung BJ, Bickle KL, Schrage KJ, Muskett P, Patel K (2006) The CLAVATA1-related BAM1, BAM2 and BAM3 receptor kinase-like proteins are required for meristem function in Arabidopsis. Plant Journal 45: 1–16.BJ DeYoungKL BickleKJ SchrageP. MuskettK. Patel2006The CLAVATA1-related BAM1, BAM2 and BAM3 receptor kinase-like proteins are required for meristem function in Arabidopsis.Plant Journal45116

58. Hord CLH, Chen C, DeYoung BJ, Clark SE, Ma H (2006) The BAM1/BAM2 Receptor-Like Kinases Are Important Regulators of Arabidopsis Early Anther Development. The Plant Cell 18: 1667–1680.CLH HordC. ChenBJ DeYoungSE ClarkH. Ma2006The BAM1/BAM2 Receptor-Like Kinases Are Important Regulators of Arabidopsis Early Anther Development.The Plant Cell1816671680

59. Sablowski R (2007) The dynamic plant stem cell niches. Current Opinion in Plant Biology 10: 639–644.R. Sablowski2007The dynamic plant stem cell niches.Current Opinion in Plant Biology10639644

60. Costa S, Shaw P (2007) ‘Open minded’ cells: how cells can change fate. Trends Cell Biol 17: 101–106.S. CostaP. Shaw2007‘Open minded’ cells: how cells can change fate.Trends Cell Biol17101106

61. Grewal SI, Rice JC (2004) Regulation of heterochromatin by histone methylation and small RNAs. Curr Opin Cell Biol 16(3): 230–8.SI GrewalJC Rice2004Regulation of heterochromatin by histone methylation and small RNAs.Curr Opin Cell Biol1632308

62. Nakayama J, Rice JC, Strahl BD, Allis CD, Grewal SI (2001) Role of histone H3 lysine 9 methylation in epigenetic control of heterochromatin assembly. Science 292(5514): 110–3.J. NakayamaJC RiceBD StrahlCD AllisSI Grewal2001Role of histone H3 lysine 9 methylation in epigenetic control of heterochromatin assembly.Science29255141103

63. Ferreira FJ, Kieber JJ (2005) Cytokinin signaling. Curr Opin Plant Biol 8: 518–525.FJ FerreiraJJ Kieber2005Cytokinin signaling.Curr Opin Plant Biol8518525

64. Decker EL, Frank W, Sarnighausen E, Reski R (2006) Moss Systems Biology en Route: Phytohormones in Physcomitrella Development. Plant Biol 8: 397–406.EL DeckerW. FrankE. SarnighausenR. Reski2006Moss Systems Biology en Route: Phytohormones in Physcomitrella Development.Plant Biol8397406

65. Schulz P, Reski R, Maldiney R, Laloue M, v. Schwartzenberg K (2000) Kinetics of cytokinin production and bud formation in Physcomitrella: analysis of a wild type, a developmental mutant and two of its ipt transgenics. Journal of Plant Physiology 156: 768–774.P. SchulzR. ReskiR. MaldineyM. LaloueK. v. Schwartzenberg2000Kinetics of cytokinin production and bud formation in Physcomitrella: analysis of a wild type, a developmental mutant and two of its ipt transgenics.Journal of Plant Physiology156768774

66. Ashton NW, Cove DJ (1997) The isolation and preliminary characterisation of auxotrophic and analogue resistant mutants in the moss Physcomitrella patens. Molecular and General Genetics 154: 87–95.NW AshtonDJ Cove1997The isolation and preliminary characterisation of auxotrophic and analogue resistant mutants in the moss Physcomitrella patens.Molecular and General Genetics1548795

67. Grimsley NH, Ashton NW, Cove DJ (1977) The production of somatic hybrids by protoplast fusion in the moss, Physcomitrella patens. MOI Gen Genet 154: 97–100.NH GrimsleyNW AshtonDJ Cove1977The production of somatic hybrids by protoplast fusion in the moss, Physcomitrella patens.MOI Gen Genet15497100

68. Ulrich I, Ulrich W (1991) High-resolution flow cytometry of nuclear DNA in plants. Protoplasma 165: 212–215.I. UlrichW. Ulrich1991High-resolution flow cytometry of nuclear DNA in plants.Protoplasma165212215

69. Johnston JS, Bennett M, Rayburn AL, Galbraith DW, Price HJ (1999) Reference standards for determination of DNA content of plant nuclei. American Journal of Botany 86: 609–613.JS JohnstonM. BennettAL RayburnDW GalbraithHJ Price1999Reference standards for determination of DNA content of plant nuclei.American Journal of Botany86609613

70. Audic S, Claverie JM (1997) The significance of digital gene expression profiles. Genome Res 7(10): 986–95.S. AudicJM Claverie1997The significance of digital gene expression profiles.Genome Res71098695

71. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125: 279–284.Y. BenjaminiD. DraiG. ElmerN. KafkafiI. Golani2001Controlling the false discovery rate in behavior genetics research.Behav Brain Res125279284

72. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, Series B (Methodological) 57: 289–300.Y. BenjaminiY. Hochberg1995Controlling the false discovery rate: a practical and powerful approach to multiple testing.Journal of the Royal Statistical Society, Series B (Methodological)57289300

73. Zeeberg BR, Qin H, Narasimhan S, Sunshine M, Cao H, Kane DW, et al. (2005) High-Throughput GoMiner, an ‘industrial-strength’ integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of Common Variable Immune Deficiency (CVID). BMC Bioinformatics 6: 168.BR ZeebergH. QinS. NarasimhanM. SunshineH. CaoDW Kane2005High-Throughput GoMiner, an ‘industrial-strength’ integrative gene ontology tool for interpretation of multiple-microarray experiments, with application to studies of Common Variable Immune Deficiency (CVID).BMC Bioinformatics6168