A transcriptional blueprint for a spiral-cleaving embryo
© The Author(s). 2016
Received: 27 January 2016
Accepted: 29 June 2016
Published: 5 August 2016
The spiral cleavage mode of early development is utilized in over one-third of all animal phyla and generates embryonic cells of different size, position, and fate through a conserved set of stereotypic and invariant asymmetric cell divisions. Despite the widespread use of spiral cleavage, regulatory and molecular features for any spiral-cleaving embryo are largely uncharted. To address this gap we use RNA-sequencing on the spiralian model Platynereis dumerilii to capture and quantify the first complete genome-wide transcriptional landscape of early spiral cleavage.
RNA-sequencing datasets from seven stages in early Platynereis development, from the zygote to the protrochophore, are described here including the de novo assembly and annotation of ~17,200 Platynereis genes. Depth and quality of the RNA-sequencing datasets allow the identification of the temporal onset and level of transcription for each annotated gene, even if the expression is restricted to a single cell. Over 4000 transcripts are maternally contributed and cleared by the end of the early spiral cleavage phase. Small early waves of zygotic expression are followed by major waves of thousands of genes, demarcating the maternal to zygotic transition shortly after the completion of spiral cleavages in this annelid species.
Our comprehensive stage-specific transcriptional analysis of early embryonic stages in Platynereis elucidates the regulatory genome during early spiral embryogenesis and defines the maternal to zygotic transition in Platynereis embryos. This transcriptome assembly provides the first systems-level view of the transcriptional and regulatory landscape for a spiral-cleaving embryo.
KeywordsSpiral cleavage Lophotrochozoan Annelid Transcriptome Maternal contribution Zygotic transcription
Over the last decade the spiralian Platynereis dumerilii has emerged as an excellent model organism for the study of development, evolution and marine biology (Fig. 1a) [6, 7]. One of the main advantages of this marine annelid is that its body plan and genome has maintained many ancient features [8, 9]. For example, the Platynereis genome has retained ancestral complements for gene families and gene structure, like the wnt family, with 12 of 13 ancient wnt genes conserved in Platynereis [10, 11]. Furthermore, another comparative study analyzing exon-intron structure shows that Platynereis genes are more similar to human genes than to genes from insects and nematodes . This suggests more conserved genomic features between this annelid and vertebrates and an increase in evolutionary changes in insect and nematode lineages (Fig. 1b). Platynereis also exhibits many common features with vertebrates including similar signatures of developmental gene expression during the formation of the brain, central nervous system and eye development, and a similar neuropeptide complement [13–17], features that were lost or strongly modified in the evolutionarily closer model systems Caenorhabditis elegans and Drosophila melanogaster. Although more recent molecular analyses have placed nereids, like Platynereis, into a more derived phylogenetic position within annelids [18, 19], Platynereis has emerged as a prominent model for comparative studies to infer early bilaterian characteristics.
Another important advantage of Platynereis is its accessibility and amenability for experimental analyses. The entire life-cycle of Platynereis can be recreated under laboratory conditions, and its lunar synchronized mating behavior makes it possible to collect thousands of synchronously developing embryos at distinct embryonic stages . Additionally, several experimental avenues have been pioneered in recent years in Platynereis including zygote microinjection, transient and stable transgenesis, and various genome-modifying technologies that allow functional studies [20, 21].
The first 14 h of Platynereis development comprises early embryogenesis, from fertilization to an early protrochophore stage (~330 cells) that has hatched from the vitelline membrane formed shortly after fertilization (Fig. 1c and d) [22–24]. At fertilization and triggered by sperm contact, the fertilized egg completes the meiotic divisions and generates two polar bodies at the animal pole before the zygote enters the first mitotic cell division shortly after 2 h post fertilization (hpf). The first two cell divisions are highly unequal giving rise to four large embryonic founder cells of different size called A, B, C, and D (Fig. 1c). The spiral cleavage mode of embryogenesis is mainly confined to the next four rounds of cell divisions. Each of the four founder cells and their progeny exhibit a similar series of asymmetric cell divisions oriented along the animal-vegetal axis of the embryo, generating animal pole and vegetal pole daughter cells. During consecutive divisions, the orientation of mitotic spindles alternate and generate quartets of daughter cells that are positioned clockwise and counterclockwise when viewed from the animal pole, assuming a spiral arrangement. Each embryonic cell can be identified by its size and position, and has been shown to have a distinct cell fate [24, 25]. The embryonic cells begin to transition to a bilaterally symmetric mode of division after 7hpf, though alternating spindle orientations can still be observed. The bilaterally symmetric pattern of cell divisions generates a pattern of similarly sized cells on the left and right side of the embryo. The larger D-quadrant bears special significance in this transition by creating two large cells, the somatoblast named ‘2d112’, and the mesentoblast named ‘4d’ [22, 24]. Both cells divide symmetrically and perpendicular to the animal-vegetal axis, segregating two founder cells each whose progeny will form the left and right side of the trunk ectoderm and trunk mesoderm [22, 24–26]. Thus, by 7hpf, the founder cells for each germ layer have been segregated, and after 9hpf, the germ cells segregate from the mesodermal cell lineages [27, 28]. Twenty-four embryonic cells have stopped dividing by 12hpf and begin to differentiate into multi-ciliated cell types forming a ciliated ring, the prototroch. By 14hpf, an early protrochophore stage of ~330 cells, Platynereis has hatched and rotates freely [23, 24].
Although the early development of Platynereis has been described and some molecular data for early embryonic stages is emerging [11, 23, 29, 30], a comprehensive transcriptome study for early embryonic stages is lacking. Current transcriptome studies are confined to genes of interests or are focused on later stages [31–33]. Comprehensive sequencing and functional annotation of a developmental Platynereis transcriptome is the first crucial step for understanding the molecular dynamics during early development, including maternal and zygotic contributions, and the complex regulatory networks that lead to the different cell fates. The advent of next-generation RNA-sequencing technologies allows us to not only reconstruct gene models, but also to obtain the precise expression level of every transcript throughout early development.
Here we describe the first comprehensive transcriptome draft during early development in Platynereis using a de novo assembly strategy. We performed mRNA deep sequencing of seven different early developmental stages using Illumina HiSeq sequencing with read lengths ~75–100 bp. 273,087 transcripts were assembled and 51,260 of the transcripts have potential protein-coding regions larger than 100 amino acids. The assembled genes were annotated by comparison of various known protein and pathway databases such as Swiss-Prot , Pfam , Gene Ontology , and KEGG . We identify around 1000 to 2520 differentially expressed genes between adjacent stages with a FDR < 0.001, and our analysis allows for the temporal onset and level of transcription to be described for each of 28,500 genes. Importantly, our data analyses identify maternal contributions and several waves of zygotic transcription during and after the completion of spiral cleavages. This Platynereis transcriptome provides the first in-depth view of the transcriptional and regulatory landscapes for a spiral-cleaving embryo, generates a molecular platform for comparison with vertebrates and other model organisms to identify conserved core biological processes necessary for early animal embryogenesis, and enables insights into plasticity and evolution of early animal development.
De novo assembly of the early Platynereis dumerilii transcriptome
Since there is no quality reference genome for Platynereis available yet, we used the software Trinity  to assemble the RNA-sequencing libraries in a genome-independent manner into transcripts. To create unified gene models, we combined all biological replicates together (~785 million paired-end reads), and assembled them in a single pass into 273,087 non-redundant transcripts belonging to 193,310 genes (N50 size: 1466 bp) (see Methods for details). We focused our subsequent analyses on sequences with potential open reading frames (ORFs) larger than 100aa (51,260 transcripts corresponding to 32,257 genes; N50: 3230 bp; N90: 1011 bp). Similar to results from previous de novo assemblies [38–40], the vast majority of the remaining 221,826 transcripts consist of shorter sequences (N50: 528 bp; N90: 235 bp) with 93 % not expressed at any stage (FPKM (fragments per kilobase per million mapped reads) < 1: 205,693 transcripts). To evaluate the quality of our transcriptome assembly without a reference genome or transcriptome, we aligned the assembled transcripts to (1) known homologous proteins (Swiss-Prot), (2) highly conserved sets of eukaryotic genes (CEGMA, BUSCO), (3) available Platynereis sequences at NCBI, and (4) a well-curated set of Platynereis wnt genes.
At the time of our analysis, there were 1775 Platynereis mRNA sequences submitted to the NCBI nucleotide database, including full-length cDNA and fragmented EST sequences. Aligning our assembly to these sequences finds 1763 out of 1775 sequences have at least one hit using BLASTN with an E-value cutoff of < 10−10. Most of these sequences (1447) can be almost completely reconstructed (90 % coverage), indicating the consistency between our assembly and the known Platynereis mRNA sequences (Fig. 3c).
The BLASTN alignment of curated wnt genes
Aligned Length (bp)
0.1 (6 h)
9.0 (10 h)
15.7 (12 h)
3.8 (14 h)
0.1 (8 h)
0.2 (6 h)
3.5 (10 h)
0.2 (10 h)
0.1 (14 h)
2.3 (14 h)
0.2 (10 h)
12.7 (8 h)
In summary, given that our de novo Platynereis transcriptome was generated exclusively from early embryonic sources it contains a remarkably high fraction of complete gene models. It recovers most core eukaryotic genes (Fig. 3b and Additional file 1: Figure S1), and a major portion of previously submitted Platynereis genes that were mostly isolated from later larval stages (Fig. 3c). Even several wnt genes that have very low expression levels (Table 1), below the sensitivity of detection by in situ hybridization in early stages , can be partially reconstructed. Therefore, we conclude that the de novo Platynereis transcriptome is of high quality, and may represent every gene that is transcribed at early embryonic stages.
Next, we used databases to identify potential protein domains (Pfam and Swiss-Prot) encoded by Platynereis transcripts, and pathways (Kyoto Encyclopedia of Genes and Genomes, KEGG) that these gene models may be involved in to better understand the genes and transcripts with potential ORFs. A total of 431,701 Pfam domains were identified in 18,146 genes (32,464 transcripts) (Fig. 4b). Our previous Swiss-Prot BLASTP search identified 17,213 genes (31,806 transcripts), with 15,690 of these genes (28,326 transcripts) overlapping with Pfam identified genes and transcripts (Fig. 4b). Furthermore, we found that 10,132 genes (18,532 transcripts) are associated with known KEGG pathways and processes. Thus, almost 50 % of all automatically annotated genes (9642) and transcripts (17,498) were found with all three databases, and 2377 (3988), 1211 (2773), and 99 (177) genes (transcripts) were uniquely identified in Pfam, Swiss-Prot, and KEGG, respectively. In summary, automated methods using Swiss-Prot, Pfam, and KEGG databases enabled us to annotate 69.2 % of all genes, and 70.5 % of all transcripts that contain a predicted ORF larger than 100 amino acids.
To further evaluate how many of the 28,580 Platynereis gene models with ORFs might be legitimate gene models, we investigated and compared supporting evidence using the previous Swiss-Prot analysis and two additional comparative analyses, finding evidence for 21,870 genes (Fig. 4c). Compared to the Swiss-Prot analysis that found support for 17,213 genes, an OrthoMCL analysis to identify orthologous genes among 18 selected species (see Methods for details) found support for 13,295 genes. Comparison to a Platynereis transcriptome assembly recently published by Achim and colleagues  found support for 13,559 genes. This is only 62 % of the 21,870 Platynereis genes with supporting evidence (47.4 % of gene models with ORFs), suggesting that our focused effort on early stages has significantly increased the number of supported gene models in Platynereis. 6710 gene models of the 28,580 gene models with ORFs in our assembly do not have supporting evidence and may constitute either novel Platynereis genes or assembly artifacts. However, 4592 genes among these 6710 are unique paralogous genes within Platynereis, and may therefore represent novel Platynereis gene families.
By comparative analysis we asked how many of the 21,870 supported gene models are conserved in other species from different metazoan taxa of various evolutionary distance (Additional file 2: Figure S2; see Fig. 1b). For genes to be considered conserved, the genes had to satisfy two main criteria: 1) shared genes must represent the best reciprocal blast hits between each species, and 2) the ORFs must share 50 % identity at the protein level. By these stringent criteria the closely related annelids Capitella telata and the leech Helobdella robusta shared 10,351 genes and 6931 genes with Platynereis, respectively. The lower number in leech is likely the result of a derived leech genome . Interestingly, Platynereis shared 8900 genes with Lottia gigantea, and 8249 with Crassostrea gigas, two mollusks more distantly related to Platynereis than the leech, suggesting slower genome evolution in both mollusks. Comparing our Platynereis gene models with ecdysozoans and deuterostomes we found higher shared numbers with invertebrate deuterostomes such as the echinoderm Strongylocentrotus purpuratus (8276), the hemichordate Saccoglossus kowalevskii (8482), and the cephalochordate Branchiostoma lanceolatum (8144), and lower shared numbers with the ecdysozoans Daphnia pulex (6434) and Drosophila melanogaster (5669), suggesting much slower genomic changes within the invertebrate deuterostome lineages than in the ecdysozoan lineages consistent with previous studies [12, 47]. This comparative analysis also found 7240 genes, and 7322 genes shared with the vertebrates Homo sapiens, and the teleost fish Danio rerio, respectively. The lower number of shared genes is likely due to faster genome evolution in the deuterostome lineage leading to the vertebrates, and also to the many additional paralogous genes in vertebrates  that obscure/interfere with the best reciprocal blast hit analysis. Comparison with the cnidarian Nematostella vectensis found 7717 genes shared with Platynereis indicating the deep evolutionary conservation of many protein-coding genes within metazoans . Due to the stringency of this analysis, the number of shared genes between the species represented here is strongly underestimated. This type of systematic analysis with fixed criteria defines the conserved ancestral gene sets present within the early Platynereis transcriptome. Although these comparisons point to species that are less and more derived within annelid, spiralian, protostome, and deuterostome model systems, they also suggest that Platynereis is one metazoan species that contains a more ancestral gene set.
Mapping and expression level estimation
To investigate the dynamics of gene expression in early embryonic stages, the normalized expression abundances were estimated for the samples from 2, 4, 6, 8, 10, 12 and 14 hpf. The total number of raw, filtered and mapped reads for each stage is depicted in Fig. 2. 93 % of reads were kept after trimming the adapter sequences and the low quality sections of the reads. Among the remaining reads, 88 % were successfully mapped to the assembled transcripts. The FPKM (fragments per kilobase per million mapped reads) for each sample was obtained by normalization of total mappable reads and transcript length. We also used a scaling normalization method called TMM (trimmed mean of M values)  to determine the TMM-normalized FPKM, and a gene or transcript was considered as “expressed” if its TMM–normalized FPKM was greater than 1. An FPKM > 1 was empirically chosen because the onset of gene expression in a single embryonic cell as determined by in situ hybridization corresponded generally to an FPKM higher than 5 (see below and previous studies [11, 29], and therefore an FPKM of 1 seemed to be an inclusionary though arbitrary cutoff. Thus, we consider an FPKM > 1 as a reasonable operational definition for a gene being expressed within our datasets. Within the first 14 h of embryogenesis 20,977 genes (34,944 transcripts) are expressed at least once during one of the seven time points sampled. In the one cell zygote at 2hpf, 11,906 genes (16,521 transcripts) are expressed, and these numbers increase until 10hpf (14,961 genes, 20,363 transcripts), and then remain constant at 12hpf and 14hpf. All expressed genes at each stage are listed in Additional file 6: Table S4.
Gene expression profiling of early stages in Platynereis development
Next, we examined the gene expression profiles of early developmental stages and focused on the expression of genes with predicted ORFs. Among the assembled gene models with predicted ORFs, there are 13,160 genes and 18,940 transcripts that are expressed in at least one of the seven stages. Figure 6c shows the correlation of gene expression between different time points including the correlation between the technical and biological replicates. In general, each individual stage shows the highest correlation with directly adjacent stages. This is significant because it indicates that the selected developmental stages are close enough to share some or most of their expression characteristics with neighboring stages. To determine whether the developmental stages are spaced close enough to capture most transcripts during early embryogenesis, we determined which genes are exclusively expressed in one stage but not in the two neighboring stages. If our sampled time points are close enough, the frequency of highly and exclusively expressed transcripts in any stage should be low, as genes with high levels of expression should be detected in neighboring stages as well. Therefore, the distribution of the highest exclusively expressed transcripts may serve as an indicator of whether the sampling chosen for the transcriptome analysis was insufficient. Indeed, the 30 highest exclusively expressed genes have low levels of expression, range from 1 to 10 FPKM (Additional file 7: Table S5), and are slightly higher at the first and last stages included in the analysis, 2hpf and 14hpf. This might indicate that some genes would continue to be expressed at the next adjacent stages 0hpf and 16hpf. Thus, these low FPKMs of the highest exclusively expressed transcripts in each stage indicate that transcriptional profiling of more intermediate stages eg 7hpf and 9hpf would not capture significantly more novel transcripts. Together, the high correlations between adjacent stages, and the low FPKM of exclusively expressed genes in each stage suggest that our dataset captures important transcriptional dynamics during this early developmental time window.
To gain insights into the stage-specific transcriptional landscape during early embryogenesis we examined the enrichment of KEGG pathways and GO terms for each stage (Additional file 8: Figure S3, Additional file 9: Figure S4; see Methods and Additional files for details). The pathways highly represented include metabolic pathways and cellular functions that generate and shape the transcriptional and translational landscape in embryos, such as the spliceosome, RNA degradation, protein synthesis, and protein degradation. Other enriched functions are associated with the fast embryonic cell divisions during early development, thereby providing insights into the transcriptional regulation of the cell cycle and DNA replication during spiral cleavages. Interestingly, several signaling pathways including Wnt, Hippo, TGF-beta, and Notch/Delta are present in this high frequency group, pointing to potentially important roles for non-cell autonomous mechanisms through these pathways during spiralian embryogenesis. High transcriptional input for key cellular functions implicated in asymmetric cell division, cell polarity and cell adhesion is also detected. Most importantly, enrichments of processes that constitute the early gene regulatory networks including basic and specific transcription factors and the signal transduction machinery potentially involved in cell lineage and cell fate specification in spiral cleaving embryos are found (see also Additional file 9: Figure S4A). Each of these processes or pathways warrants a more in-depth description of specific components and these will be targeted in follow-up studies.
Overall, the enrichment of pathways and processes appears to be similar between different stages with the largest observed differences between the zygote (2hpf) and the other 6 stages. This was also indicated in a correlation analysis (Additional file 8: Figure S3D, Additional file 9: Figure S4D). These differences may be attributed to fast drops in the level of maternally provided transcripts between 2hpf and 4hpf by decay of distinct maternal mRNA species, and to early mRNA synthesis of others. A second, less prominent cluster of processes are mostly enriched in stages after 10hpf and may be involved in the early onset of cell differentiation in some cells, such as the cells that will form the first differentiated cell types involved in motility and light perception [22–24] (Additional file 8: Figure S3B, C, Additional file 9: Figure S4B, C).
With the assembled Platynereis early transcriptome, it is possible to dissect the dynamics of the transcriptional landscape between different embryonic stages. We assessed this by (1) determining the number of differentially expressed genes between stages, (2) identifying every gene whose expression is up-regulated or down-regulated between adjacent stages, and (3) defining clusters of genes whose transcriptional developmental profile during the first 14 h of development was similar. Each of these analyses offers unique opportunities to identify overall trends during early spiralian embryogenesis and to detect scores of individual candidate genes whose temporal expression appears to be tightly regulated during development.
The RNA-sequencing approach used here can provide a developmental profile for individual genes that show both the transcript level and the timing of the expression. Comparing all expression profiles can reveal individual and global embryonic regulatory mechanisms. This type of analysis begins by identifying all significantly up-regulated and down-regulated transcripts between adjacent stages, including only those genes that increase to a level that is above 10 FPKM in the later stage or decrease from a level that is above 10 FPKM in the earlier stage, respectively. Using these criteria, hundreds of genes have been identified with significant gene expression differences between each adjacent stage (Additional file 10: Tables S6, Additional file 11: Tables S7). While the identity of specific transcripts is useful, this analysis is also informative for discovering larger trends during early spiralian embryogenesis.
During early stages, there is a dominant trend of transcript down-regulation, with 858, 618, 663, and 1041 down-regulated genes between 2–4hpf, 4–6hpf, 6–8hpf, 8–10hpf, respectively. Only 152 and 478 are down regulated between 10–12hpf and 12–14hpf, respectively (Additional file 10: Table S6). Examples for down-regulated genes include genes with previously known early embryonic functions such as Regulators of G-protein signaling (RGS) proteins that are implicated in regulating asymmetric cell division, regulators of the cell cycle like Cyclins and CDK20, and replacing ‘inhibitory’ histone marks necessary for zygotic genome activation. Interestingly, there is a dramatic drop in transcript level for ribosomal proteins between 8 and 10hpf that seems to coincide with a major onset of zygotic gene expression in the embryo.
In an opposing trend, up-regulation of genes is lowest between the first stages (125 genes between 2 and 4hpf), and highest between the last stages (1131 genes between 12 and 14hpf) (Additional file 11: Table S7). Interestingly, the list of genes in the later stages (10hpf to 12hpf) exhibit the onset of gene expression for differentiation genes such as those involved in ciliogenesis, like tektins and axonemal components, as well as hatching enzymes. Between 12 and 14hpf new waves of synthesis of ribosomal and cytoskeleton components may indicate replenishing and remodeling of the machinery for protein synthesis and the cytoskeleton.
In this study, we generated the first transcriptomic map of Platynereis dumerilii from more than one and a half billion reads sequenced during seven early developmental stages. Such deep sequencing provides an opportunity to explore gene expression profiles with extremely high resolution. As the reference genome is not available yet, a de novo transcriptome assembly strategy was employed. We evaluated and demonstrated the high quality of our assembly in terms of the number of full-length reconstruction for the known Platynereis dumerilii cDNA sequences as well as homologous genes. The correlation analysis also indicates a high level of reproducibility in our replicates. Our data provides a detailed blueprint of spiralian embryogenesis that will impact future research endeavors in many realms of biology including: (1) fundamental processes of development such as the maternal to zygotic transition and embryonic gene regulatory networks, (2) core cellular processes such as the regulation of cell cycle and asymmetric cell division, (3) comparative analyses of evolution and development such as ancestral reconstruction of maternal contributions to eggs, germ layer specification, and the deciphering of gene regulatory networks that specify lineages and cell type, as well as (4) genome evolution illuminating the content of the ancestral bilaterian transcriptome and genome.
The early transcriptome: expanding the Platynereis dumerilii toolbox
The presented early transcriptome for Platynereis dumerilii is an additional valuable molecular resource, adding to several recent efforts to define the broader transcriptome in Platynereis . Previously published Platynereis transcriptomes were more focused, such as the immune-related transcriptome , the neuropeptide complement throughout the life cycle , and developmental processes and cell type specification in later larval stages [33, 52]. Our transcriptome analysis is the first to center on early embryogenesis and is therefore the first to focus on spiral cleaving stages, the maternal to zygotic transition, early cell lineage and germ layer specification, and the emergence of the first differentiated cell types. Although the depths and scope of the presented early transcriptome are unprecedented so far, the previous works with different sequencing technologies and a focus on later stages are advantageous to synergistically define the broader Platynereis transcriptome. Together, the transcriptomes add to the growing molecular toolbox for this lophotrochozoan spiralian model system by allowing a systems level analysis of any biological process.
Other spiralian transcriptome resources
In recent years several spiralian/lophotrochozoan transcriptome and genome resources were established. These include the first reported spiralian/lophotrochozoan genomes and/or transcriptomes for two plathyhelminths Schistosoma species (Schistosoma mansoni  and Schistosoma japonicum ), the annelids Capitella teleta and Helobdella robusta , the mollusks Lottia gigantea , Aplysia californica , Ilyanassa obsoleta , Crepidula fornicate , the polyclad flatworm Maritigrella crozierie , and the annelids Pristina leidyi and Hermodice carunuculata [59, 60], with most sequencing derived from mixed stage cDNA-libraries. In addition, several recent extensive phylogenomic studies contributed transcriptomes from dozens of annelid species [4, 19, 61, 62]. Although not focused on early stages, these transcriptomes provide an enormous resource to investigate the conservation and evolution of protein coding genes within annelid phylogeny. Stage-specific transcriptome sequencing has been done for three selected stages (1 cell, 2 cell and 32 cells) for the pond snail Lymnaea stagnalis , and stages throughout the life cycle of the pacific oyster Crassostrea gigas  including a similar set of early embryonic stages to the presented Platynereis dataset. Although the oyster sequencing data is not as comprehensive as our Platynereis data, it will be an exciting source to unravel conserved and derived molecular mechanisms between the spiral cleaving programs of mollusks and annelids.
Identifying a comprehensive “Parts” list for early spiralian embryogenesis
The depth and scope of the presented sequencing effort was designed to capture every potentially expressed transcript during the first 14 h of Platynereis development. The correlation analysis between adjacent stages (Fig. 6c) and the low frequency of exclusively expressed, stage-specific genes (Additional file 7: Table S5) suggest that only a few transcripts below 2–3 FPKM might have been missed in this dataset. Thus, the > 20,000 transcribed genes with an expression level of > 1FPKM (Additional file 6: Table S4) likely represents a comprehensive early transcriptome, yielding an extensive parts list for every process that comprises spiralian embryogenesis. To subdivide the extensive parts list, 60 GO term classes were used to classify each gene further and provide new parts lists for each process or function within the annotated transcriptome. Thus, these identified and categorized genes comprise a blueprint for early Platynereis development, will be instrumental in dissecting spiralian embryogenesis in Platynereis, and will facilitate cross species comparisons among spiralians and other metazoans.
Dynamics of the early Platynereis transcriptomes during spiral cleavages (I to V)
The Platynereis transcriptome described here was assembled from deep sequencing to allow for the opportunity to (1) explore gene expression profiles with extremely high resolution, and (2) define broader trends in overall gene transcription. In general, the transcriptional landscape is very dynamic with ~ 200–2500 genes being differentially expressed between each stage. Further, we found that more genes are down-regulated between each adjacent stage during the first four stages sampled (2hpf to 8hpf), while up-regulation of genes between adjacent stages is more prominent between the later stages (10hpf to 14hpf). This suggests a major switch in the transcriptional landscape between 8 and 10hpf, and this timing is corroborated by the cluster analysis of differentially expressed genes (Fig. 8). We believe this marks the transition of maternal to zygotic gene transcription in Platynereis embryogenesis, the Maternal to Zygotic Transition (MZT) .
Maternal contribution to a spiralian egg
The first stage included in our transcriptome was 2hpf, the Platynereis zygote. There are four maternal clusters of genes (Clusters 1–4) that share the trend of moderate to high enrichment in the zygote and then a decrease in level of expression during consecutive stages (Fig. 8). While the early stages, 2hpf to 8hpf, share a high number of expressed genes (7107 genes), the KEGG and GO term enrichment analysis between stages found that the 2hpf zygote differs substantially from the other early stages (Additional file 9: Figure S4, Additional file 10: Tables S6, Additional file 12: Tables S8, Additional file 13: Table S9, Additional file 14: Table S10). For example, in the Platynereis zygote there are enrichments for functions associated with nucleotide binding, protein degradation and synthesis, and the cell cycle, and this is consistent with a previous analysis that compared maternal contributions of a mollusk egg, the snail Lymnaea stagnalis to different bilaterians . Although considered mostly housekeeping functions, the dynamics of these gene cohorts facilitates new insights into the orchestration of genes in early development and warrants follow-up studies. For example, individual genes in Cluster 1 (136 genes) that show an extreme decrease in expression in the early stages, like the RHOGTP-activating enzymes and several Regulator of G-protein signaling molecules implicated in regulation of early cell division in other early embryos [66–68], will be especially interesting to study (Fig. 8; Additional file 12: Table S8). Future analysis with this transcriptome data for the Platynereis zygote can begin to answer fundamental questions such as ‘What is inside the egg?’ or ‘How similar is the egg composition among different metazoans?’ or ‘What are conserved and novel maternally provided components between different egg types?’. This can open research avenues to further define maternal contributions and enable cross species comparisons, promising general insights into early embryogenesis, and evolution and development.
Clearance of maternal transcripts in a spiralian embryo
Clearance of maternal transcripts is one of the two major mechanisms for the transfer of developmental control during the maternal to zygotic transition in embryos , and this clearance is visible in the Platynereis transcriptome presented here. The four clusters with maternally provided genes (Clusters 1–4, Fig. 8) each show different degrees and timing in the down-regulation of maternal transcripts. In Cluster 1, the maternal transcripts are cleared rapidly between 2 and 4hpf. In contrast, in Cluster 4, the maternal transcripts are cleared later and much more slowly. This suggests the possible existence of several different mechanisms for maternal clearance in Platynereis.
Early clearance mechanisms often target regulators of the cell cycle . In this regard it is interesting that transcripts for a cdc25/string-like gene are dramatically decreased from above 250 FPKM to below 15 FPKM between 2 and 4hpf. Our cluster analysis also suggests a final clearance step for distinct transcript species between 8 and 10hpf (Fig. 8). Interestingly, this includes a very pronounced decay of many ribosomal transcripts (Additional file 10: Table S6) that may constitute a defining signature for the end of maternal clearance in Platynereis. Thus, the transcriptional landscape for early Platynereis development shown here may be used for the discovery and regulation of clearance mechanisms, and for understanding the shift from maternal to zygotic control of the cell cycle in early embryos.
Zygotic activation of gene expression
The second major mechanism for the transfer of developmental control during the maternal to zygotic transition in embryos constitutes the zygotic activation of gene expression . Similar to several recent studies that utilize RNA-sequencing or microarray approaches to investigate the MZT in various other species, Platynereis embryos seem to exhibit early waves of zygotic gene expression: 47 genes show peak expression at 4hpf, 225 genes at 6hpf, and 370 genes at 8hpf (Fig. 8; Additional file 12: Table S8). Concurrent with the later peak, a first larger cohort of genes (Clusters 9, 11 and 12; Fig. 8) exhibits a moderate increase in overall gene expression level between 6 and 8hpf, followed by a more pronounced increase of a second major cohort of genes between 8 and 10hpf (Clusters 13 and 14, Fig. 8). Thus, the major wave of zygotic transcription appears to coincide with the end of the 6th or 7th cycle of cell division after 8hpf, and at the transition to bilaterally symmetric cell divisions in some cell lineages in Platynereis [22, 24]. Although more experiments are necessary to identify if these changes in transcript level correspond to active zygotic gene expression , the presented developmental transcriptional landscape provides hints for future studies.
Early small waves of zygotic expression
Our study identified several early small waves of zygotic expression comprising three clusters (Clusters 6, 7, and 8, Fig. 8; Additional file 12: Table S8). The three clusters contain 642 genes including many transcription factors such as members of the glia cell missing (GCM)-, forkhead (Fox)-, paired box (PRD)-, hes-like, and zinc finger-like gene families that will be further described in subsequent studies. These small, early waves of transcription have been found in several species (for review see ), and it will be interesting to see whether common molecular mechanisms will be found. However, many of the early expressed genes appear to be more derived genes within their particular conserved gene families. For example, preliminary phylogenetic analyses of early expressed transcription factors identified at least seven derived Fox-like, two derived GCM-like, three derived Hes-like, and three derived PRD-like genes that defy easy classification within their conserved gene families (Additional file 12: Table S8; data not shown). This observation is in agreement with recent findings that novel, lophotrochozoan-specific homeobox genes are preferentially recruited to and utilized in early development in the oyster Crassostrea gigas . Paps and colleagues (2015) suggest the observations with Crassostrea support an hourglass or phylotypic stage model of developmental evolution, predicting more novelty and evolutionary lability in the regulation of early developmental stages. Rapid divergence and novelty in transcriptional regulation of early embryogenesis is further supported by the many divergent and derived GATA transcription factors that exist in nematodes like C. elegans, with many crucial functions in the early endo- and mesodermal gene regulatory networks [70–72]. Consistent with these observations is also a recent comparative study that found that early zygotic genes in zebrafish exhibit little evolutionary conservation and are composed of mostly novel genes . Thus, these observations in Platynereis (this study), Crassostrea , C. elegans , and zebrafish , may lend support to the phylotypic stage model of developmental evolution (see also Levin and colleagues  for alternative observations), and may also reflect recent notions that the conserved pattern of spiral cleavage might be regulated differently within and between different spiralian clades .
Later major waves of zygotic expression
In Platynereis, the major waves of zygotic transcription begin between 6 and 8hpf with a modest increase in expression in ~3500 genes (Clusters 11 and 12, Fig. 8). This major wave of zygotic transcription is shortly followed by another wave of zygotic expression between 8 and 10hpf (Clusters 10, 13, and 14). The second wave of zygotic expression includes ~2100 genes, most of which have a higher increase in expression level compared to the waves between 6 and 8hpf. The gene clusters from the two major waves of zygotic expression include many components of gene regulatory networks including conserved transcription factors and signaling components (see below). Many critical components of major pathways are expressed like receptors and ligands of Wnt pathways such as frizzled 1/2/7, sFRP1/2/5, wnt-4 and wnt-5 [11, 29]. These clusters also harbor many genes that support the formation of the first differentiating cell types in the Platynereis embryo such as early onset ciliogenesis genes (Additional file 12: Table S8). Thus, these clusters comprise of many crucial developmental regulators and the first differentiating gene cascades, all of which reflect the hybrid nature of the developing protrochophore that is composed of developing cell lineages as well as the first functional cell types [22, 24]. Due to the higher amplitude of the major wave of zygotic transcription between 8 and 10hpf, we suggest that this period represents the major transition to zygotic control in the Platynereis embryo.
Comparison of the maternal to zygotic transition in Platynereis to other species
The fine-scale transcriptional landscape of the MZT has been reported recently for several animal model systems (reviewed in ). Detailed comparative analyses have just begun , but are complicated by the many unique early developmental modes that need to be considered. Different developmental modes greatly change the early transcriptional landscape, especially the time point and scale of zygotic gene activation, which can differ depending on size of the embryo, the speed of early cell divisions, and the time when gastrulation occurs. Early zygotic gene expression has been reported as early as the formation of the male and female pronuclei in the nematode Ascaris suum and in the 1-cell mouse embryo (both exhibit very slow development) [75–77], but also in the 1-cell stage of the much faster developing sea urchin embryo . Many animal embryos exhibit early and small waves of zygotic expression followed later by larger zygotic waves of gene expression [73, 79–85] similar to what is seen with our Platynereis dataset. The scale of early and late zygotic transcription involving hundreds to several thousands of genes, respectively, appears to be similar between Platynereis and other animal model systems , and in all animal embryos, including Platynereis, the major zygotic activation is well underway before the beginning of gastrulation.
Systems-level insights into the early embryogenesis of Platynereis
Remarkably, the transcript level for individual genes in biological replicates at the same stage was very similar (Figs. 6b, c and 8), highlighting three important aspects of this study: (1) the accuracy of the RNA-sequencing approach to capture the level of thousands of transcripts reliably within a sample, (2) the suitability, amenability and strength of the Platynereis model system that enables the collection of thousands of synchronized developmental stages  and subsequent sequencing library preparations without necessary amplification steps, and (3) the apparent tight regulation within developing Platynereis embryos that determines the exact level of each transcript. The observed level of accuracy seen in the dataset suggests each mechanism in early development is tightly controlled on the transcriptional level.
Thus, this study accurately captures and dissects the changing genome-wide transcriptional landscape between a one-cell zygote to a hatched and rotating 330-cell stage protrochophore that exhibits the first differentiated cell types. Further, the Platynereis dataset presented here reveals the complete transcriptional input for any process related to early spiralian embryogenesis, whether considered ‘house keeping’ or ‘developmental regulative’ functions (Fig. 5, Additional file 3: Table S1, Additional file 4: Table S2, Additional file 5: Table S3, Additional file 8: Figure S3, Additional file 9: Figure S4). For example, there is a prominent and coordinated up-regulation of ribosomal transcripts and components of microtubules between 12 and 14hpf, and expectedly, up-regulation of ciliogenesis genes between 8 and 12hpf when the cells that will form a ciliated ring stop dividing and start to differentiate [22, 24, 28]. Furthermore, our dataset identified the evident expression of several key developmental pathways such as Wnt, Hippo, Notch, and TGF-beta during early embryonic stages. These findings will enable future studies to begin systems-level approaches to understand the coordinated regulation of hundreds of components that constitute these individual pathways throughout early development. We demonstrate that by determining the transcriptional profiles between close enough stages one can capture the outline of the underlying gene regulatory networks in embryos comprehensively.
Defining the gene regulatory networks in early spiralian development
One of the motivations for this study was to identify all components of each gene regulatory network (GRN) that operates during early spiralian embryogenesis in the annelid Platynereis dumerilii (Additional file 6: Table S4). The functional annotation of our dataset allowed us to identify hundreds of transcription factors and several signaling pathways that may constitute the embryonic GRNs. Furthermore, our transcriptional developmental profiles for each gene expressed identifies the exact time during development when a gene is turned on and when its transcripts are depleted. The developmental profiles and the analysis to determine each gene up-regulated or down-regulated between adjacent stages helped to pinpoint the onset of expression for many transcription factors, signaling ligands, receptors, and intracellular components of signaling pathways during early spiralian development (Additional file 10: Tables S6, Additional file 11: Tables S7, Additional file 12: Tables S8).
To validate the data from the RNA-sequencing approach, specific genes that are important developmental regulators were cloned and their spatial and temporal expression was determined in early Platynereis embryos (Fig. 9). The showcased genes have low to moderate levels of expression, with peak values between 9 and 30 FPKM, to highlight the ability and sensitivity of the data to capture lowly expressed transcripts. Each biological replicate has similar FPKM values, and we believe that any observed differences are due to slightly accelerated development in some batches of our embryos caused by minute local fluctuations in our temperature-controlled culture room. Using whole mount in situ hybridization to localize the transcripts, the onset of expression can be seen in a single cell (nodal, bmp2/4; Fig. 9a, c), or a few cells (hes-like 1, fz5/8; Fig. 9b, d). The expression can even be visualized at the two chromosomal loci in individual nuclei (hes-like 1, bmp2/4; Fig. 9b, c). These results demonstrate the transcriptional landscape of early Platynereis development produced by RNA-sequencing is reproducible, sensitive and predictive for low to moderate level expression of developmental regulators, and therefore can provide a comprehensive map outlining embryonic GRNs in this spiralian embryo.
In summary, we have provided the first comprehensive transcriptome study of spiralian embryogenesis for the marine annelid Platynereis dumerilii, reporting the sequence, annotation, and stage-specific expression level for 28,580 genes. The unprecedented depth in sequencing of seven developmental stages identified the dynamic and invariant transcriptional landscape from the one cell zygote to ~330 cell stage, including the maternal contribution of 11,904 genes, the maternal clearance of over 4000 expressed genes, and both the minor and major waves of zygotic expression. Quantifying the expression levels for each expressed gene at each developmental stage analyzed captures the transcriptional input into every biological process during spiralian embryogenesis and outlines the gene regulatory networks specifying cell lineages and germ layers in the Platynereis embryo. Our gene models and functional annotation in this spiralian/lophotrochozoan model system may serve as a valuable resource to further decipher the gene regulatory networks patterning cell lineages and cell fates during early spiralian development, providing crucial systems-level data to infer conserved and novel molecular features of early bilaterian development and evolution.
Platynereis dumerilii culture, embryo collection, cDNA preparation, and sequencing
All Platynereis embryos were obtained from a single breeding culture maintained at Iowa State University according to standards and protocols established at www.platynereis.de. In order to ensure consistent developmental timing and reproducibility, embryos were kept in an 18 °C incubator where the temperature is continuously monitored. Newly fertilized eggs were monitored for quality including formation of jelly coat and normal cleavage patterns before selection for RNA collection. Biological replicates for each time point were obtained from independently fertilized batches of embryos developing synchronously under identical conditions. Embryos were collected at 2, 4, 6, 8, 10, 12, and 14hpf and homogenized in Trizol (Ambion). In each case ~100–200 embryos were left behind as quality controls to be monitored for normal development for 48 h before proceeding to RNA extraction from Trizol in accordance with the manufacturer’s protocol. Extracted RNA was treated with the RNase-free DNase set (QIAGEN), purified with the RNeasy Mini Kit (QIAGEN), and checked for RNA degradation on a 1 % agarose gel. Additional total RNA quality controls using the Bioanalyzer system (Agilent), preparation of each barcoded Illumina mRNA-seq library, and Illumina deep sequencing with 75 bp–100 bp paired-end reads (4 samples per lane) were performed by the Genome Sequencing & Analysis Core Resource at Duke Institute for Genome Sciences and Policy using an Illumina HiSeq sequencing system.
Read processing and de novo assembly
The raw reads were preprocessed with Trimmomatic  to remove the adapter sequences and low quality regions. Bases were removed from both ends if the Phred score was lower than 20. A sliding window of 4 was used to continue removing bases with the threshold of an average score of less than 20. The quality of the preprocessed reads were analyzed and visualized by FASTX-Toolkit . To create unified gene models, all samples were combined together and assembled in a single pass. Trinity’s in silico normalization utility was used, with a k-mer size of 25, to reduce the memory requirement and improve the running time. Trinity assembled the normalized reads into 357,961 transcripts with an N50 size of 2331 bp. The digital normalization and assembly were performed on the Stampede high performance computing system, which is a part of Texas Advanced Computing Center (TACC) . With the final assembly, transcripts with more than 95 % identity were further clustered using CD-HIT’s “cd-hit-est” program . From this 273,087 transcripts were obtained with an N50 size of 1466 bp. These transcripts belong to 193,310 genes and 50,237 of the genes are longer than 1000 bp.
Functional annotation of the Platynereis transcriptome
The putative coding regions of the final transcripts were predicted using TransDecoder , which can extract the long ORFs and estimate the posterior probability of all six possible open reading frames using the Markov Model. The likely protein sequences were aligned to the Swiss-Prot database  using BLASTP search with options “-evalue 1e-10 -max_target_seqs 1”. We also made use of the curated gene ontology and eggNOG information provided by Swiss-Prot FTP site (ftp://ftp.uniprot.org/pub/databases/uniprot). The Pfam  domains of the protein sequences were predicted using default parameters in the hmmscan program in HMMER . For the KEGG pathway analysis , KEGG’s REST APIs were used to download the KEGG curated protein sequences. A BLASTP search (E-value < 10−10) was employed to identify the homologous KEGG genes and the associated pathways. An R package, VennDiagram, was used to count the number of annotated genes and plot the Venn diagram. All annotation information was uploaded to a MySQL database and the data used in enrichment analyses were filtered by custom SQL commands.
To perform whole transcriptome comparison to other metazoan species, potential orthologous groups were identified using the OrthoMCL pipeline . Except for our Platynereis dumerilii gene models, protein sequences were collected from 17 species: Capitella teleta , Helobdella robusta , Lottia gigantea , Crassostrea gigas , Daphnia pulex , Tribolium castaneum , Drosophila melanogaster , Strongylocentrotus purpuratus , Saccoglossus kowalevskii , Branchiostoma floridae , Danio rerio , Xenopus tropicalis , Homo sapiens , Nematostella vectensis , Amphimedon queenslandica , Trichoplax adherens  and Monosiga brevicollis . Only the longest splicing isoform for each gene was included in the comparison. An all-against-all BLASTP search was performed using 1e-05 cutoff. Any hits with less than 50 % identity match were removed. The orthologous relationships were confirmed if genes are the reciprocal best hits for any two species. 40,206 groups were identified using the MCL clustering algorithm , and 32,482 groups have at least two species.
Our models were also compared with another transcriptome assembly from Platynereis preferentially generated by sequencing mixed and later larval stages  that includes 102,433 transcripts. BLASTN was used to search our gene models at the nucleotide level against the other assembly with an E-value cutoff 1e-20.
Assessment of the Platynereis transcriptome assembly
To evaluate the quality of our assembly, the assembled transcripts were aligned to known orthologous models from the Swiss-Prot database, the NCBI Platynereis database, the CEGMA database, and the BUSCO datasets  The completeness was measured as the percentage of target sequences aligned to the Platynereis gene models in the transcriptome presented here. BLAST searches and an in-house Perl script was implemented to parse the BLAST results and to calculate the coverage. Two CEGMA datasets (Caenorhabditis elegans and Drosophila melanogaster) were downloaded from CEGMA website (http://korflab.ucdavis.edu/datasets/cegma/#SCT5). For BUSCO analysis, three BUSCO datasets (Arthropods, Metazoans and Eukaryotes) were downloaded from BUSCO website (http://busco.ezlab.org/). The BUSCO main script (BUSCO_v1.1b.py) was used to generate assessment results with option “-m trans”.
Mapping and expression level estimation
The read count of each sample was calculated by RSEM, which used Bowtie  as an aligner to map the raw reads back to the assembled transcripts. The number of aligned reads was calculated using SAMtools “flagstat” program . An R package “edgeR”  was used to incorporate the TMM scaling factor and to calculate the differentially expressed genes. The differentially expressed genes were filtered by False Discovery Rate (FDR) adjusted p-value < = 0.001. For the clustering analysis, only genes with an FPKM > 1 in at least one stage were selected and the genes were clustered into 32 groups using median-centered log2-transformed FPKM with the R function “hclust” (hierarchical clustering). The clusters were further manually combined into 15 clusters and were ordered by the expression patterns using a custom Perl script. The heatmap was generated by the R function “heatmap.2”. For the correlation analysis, the pairwise Spearman correlation between all biological and technical replicates was calculated using the R function “cor” and the correlation scatter plots were visualized using the R function “plot”.
Cloning genes and whole-mount in situ hybridization
All cloning, anti-sense RNA probe preparation, and whole-mount in situ hybridizations (WMISH) were performed as previously published [11, 29]. In short, gene sequences were found using the assembled transcriptome, primers designed with Primer3 , stage-specific cDNA was used with the gene-specific primers to amplify the target gene. Standard cloning procedures were used to clone the PCR-amplified gene target into a vector that allows for anti-sense RNA transcription. Ages of embryos used for WMISH was based on the expression profile for the specific target. WMISH procedure was performed as previously described [11, 29].
BLAST, Basic Local Alignment Search Tool; bp, base pair; CEGMA, Core Eukaryotic Genes Mapping Approach; EST, Expressed sequence tag; FPKM, Fragments Per Kilobase per Million mapped reads; GO, gene ontology; GRN, gene regulatory networks; hpf, hours post fertilization; KEGG, Kyoto Encyclopedia of Genes and Genomes; MCL, Markov Cluster Algorithm; MZT, maternal to zygotic transition; nt, nucleotide; ORF, open reading frame; TMM, Trimmed Mean of M-values; WMISH, whole-mount in situ hybridization
We appreciate Dr. Volker Brendel’s essential support in the initial phase of the project. We would like to thank Edward Letcher, Roy Holmes, Ben Wu, and Kali Levsen for handling the polychaete culture. We would also like to thank the anonymous reviewers for comments on the manuscript.
Funding for this work was provided by the Roy J. Carver Charitable Trust and the National Science Foundation (Award ID 1455185) to SQS.
Availability of data and material
The data sets supporting the results of this article are available under the URL:
HCC and SQS conceived of the study. HCC performed the bioinformatics tasks and drafted the manuscript. MMP cloned genes, carried out in situ hybridizations in early stages, and helped to write the manuscript. BRB collected the RNA for sequencing, cloned genes, and performed the phylogenetic analyses. SQS designed and coordinated the study, and wrote the manuscript. All authors read, corrected and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Henry JQ. Spiralian model systems. Int J Dev Biol. 2014;58(6–8):389–401.View ArticlePubMedGoogle Scholar
- Dunn CW, Giribet G, Edgecombe GD, Hejnol A. Animal Phylogeny and Its Evolutionary Implications*. Annu Rev Ecol Evol Syst. 2014;45(1):371–95.View ArticleGoogle Scholar
- Hejnol A. A twist in time--the evolution of spiral cleavage in the light of animal phylogeny. Integr Comp Biol. 2010;50(5):695–706.View ArticlePubMedGoogle Scholar
- Laumer CE, Bekkouche N, Kerbl A, Goetz F, Neves RC, Sorensen MV, Kristensen RM, Hejnol A, Dunn CW, Giribet G, et al. Spiralian phylogeny informs the evolution of microscopic lineages. Curr Biol. 2015;25(15):2000–6.View ArticlePubMedGoogle Scholar
- Telford MJ, Budd GE, Philippe H. Phylogenomic Insights into Animal Evolution. Curr Biol. 2015;25(19):R876–87.View ArticlePubMedGoogle Scholar
- Fischer A, Dorresteijn A. The polychaete Platynereis dumerilii (Annelida): a laboratory animal with spiralian cleavage, lifelong segment proliferation and a mixed benthic/pelagic life cycle. Bioessays. 2004;26(3):314–25.View ArticlePubMedGoogle Scholar
- Raible F, Tessmar-Raible K. Platynereis dumerilii. Curr Biol. 2014;24(15):R676–7.View ArticlePubMedGoogle Scholar
- Hui JH, Raible F, Korchagina N, Dray N, Samain S, Magdelenat G, Jubin C, Segurens B, Balavoine G, Arendt D, et al. Features of the ancestral bilaterian inferred from Platynereis dumerilii ParaHox genes. BMC Biol. 2009;7:43.View ArticlePubMedPubMed CentralGoogle Scholar
- Raible F, Arendt D. Metazoan Evolution: Some Animals Are More Equal than Others. Curr Biol. 2004;14(3):R106–8.View ArticlePubMedGoogle Scholar
- Cho SJ, Valles Y, Giani Jr VC, Seaver EC, Weisblat DA. Evolutionary dynamics of the wnt gene family: a lophotrochozoan perspective. Mol Biol Evol. 2010;27(7):1645–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Pruitt MM, Letcher EJ, Chou HC, Bastin BR, Schneider SQ. Expression of the wnt gene complement in a spiral-cleaving embryo and trochophore larva. Int J Dev Biol. 2014;58(6–8):563–73.View ArticlePubMedGoogle Scholar
- Raible F, Tessmar-Raible K, Osoegawa K, Wincker P, Jubin C, Balavoine G, Ferrier DE, Benes V, de Jong P, Weissenbach J, et al. Vertebrate-Type Intron-Rich Genes in the Marine Annelid Platynereis dumerilii. Science. 2005;310(25):2.Google Scholar
- Arendt D, Tessmar-Raible K, Snyman H, Dorresteijn AW, Wittbrodt J. Ciliary Photoreceptors with a Vertebrate-Type Opsin in an Invertebrate Brain. Science. 2004;306:3.View ArticleGoogle Scholar
- Conzelmann M, Williams EA, Krug K, Franz-Wachtel M, Macek B, Jekely G. The neuropeptide complement of the marine annelid Platynereis dumerilii. BMC Genomics. 2013;14(906):15.Google Scholar
- Denes AS, Jekely G, Steinmetz PR, Raible F, Snyman H, Prud’homme B, Ferrier DE, Balavoine G, Arendt D. Molecular architecture of annelid nerve cord supports common origin of nervous system centralization in bilateria. Cell. 2007;129(2):277–88.View ArticlePubMedGoogle Scholar
- Tessmar-Raible K, Raible F, Christodoulou F, Guy K, Rembold M, Hausen H, Arendt D. Conserved sensory-neurosecretory cell types in annelid and fish forebrain: insights into hypothalamus evolution. Cell. 2007;129(7):1389–400.View ArticlePubMedGoogle Scholar
- Tomer R, Denes AS, Tessmar-Raible K, Arendt D. Profiling by image registration reveals common origin of annelid mushroom bodies and vertebrate pallium. Cell. 2010;142(5):800–9.View ArticlePubMedGoogle Scholar
- Weigert A, Bleidorn C. Current status of annelid phylogeny. Organ Divers Evol. 2016;16(2):345–62.View ArticleGoogle Scholar
- Weigert A, Helm C, Meyer M, Nickel B, Arendt D, Hausdorf B, Santos SR, Halanych KM, Purschke G, Bleidorn C, et al. Illuminating the base of the annelid tree using transcriptomics. Mol Biol Evol. 2014;31(6):1391–401.View ArticlePubMedGoogle Scholar
- Backfish B, Rajan VBV, Fischer RM, Lohs C, Arboleda E, Tessmar-Raible K, Raible F. Stable transgenesis in the marine annelid Platynereis dumerilii sheds new light on photoreceptor evolution. Proc Natl Acad Sci U S A. 2013;110(1):6.Google Scholar
- Zantke J, Bannister S, Rajan VB, Raible F, Tessmar-Raible K. Genetic and genomic tools for the marine annelid Platynereis dumerilii. Genetics. 2014;197(1):19–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Dorresteijn AWC. Quantitative analysis of cellular differentiation during early embryogenesis of Platynereis dumerilii. Roux’s Arch Dev Biol. 1990;199:17.View ArticleGoogle Scholar
- Fischer AHL, Henrich T, Arendt D. The normal development of Platynereis dumerilii (Nereididae, Annelida). Front Zool. 2010;7(31):39.Google Scholar
- Wilson EB. The Cell-Lineage of Nereis. J Morphol. 1892;6(3):119.View ArticleGoogle Scholar
- Ackermann C, Dorresteijn A, Fischer A. Clonal domains in postlarval Platynereis dumerilii (Annelida: Polychaeta). J Morphol. 2005;266(3):258–80.View ArticlePubMedGoogle Scholar
- Fischer AH, Arendt D. Mesoteloblast-like mesodermal stem cells in the polychaete annelid Platynereis dumerilii (Nereididae). J Exp Zool B Mol Dev Evol. 2013;320(2):94–104.View ArticlePubMedGoogle Scholar
- Rebscher N, Zelada-Gonzalez F, Banisch TU, Raible F, Arendt D. Vasa unveils a common origin of germ cells and of somatic stem cells from the posterior growth zone in the polychaete Platynereis dumerilii. Dev Biol. 2007;306(2):599–611.View ArticlePubMedGoogle Scholar
- Schneider SQ, Bowerman B. beta-Catenin asymmetries after all animal/vegetal- oriented cell divisions in Platynereis dumerilii embryos mediate binary cell-fate specification. Dev Cell. 2007;13(1):73–86.View ArticlePubMedGoogle Scholar
- Bastin BR, Chou HC, Pruitt MM, Schneider SQ. Structure, phylogeny, and expression of the frizzled-related gene family in the lophotrochozoan annelid Platynereis dumerilii. Evodevo. 2015;6:37.View ArticlePubMedPubMed CentralGoogle Scholar
- Pfeifer K, Schaub C, Domsch K, Dorresteijn A, Wolfstetter G. Maternal inheritance of twist and analysis of MAPK activation in embryos of the polychaete annelid Platynereis dumerilii. PLoS One. 2014;9(5):e96702.View ArticlePubMedPubMed CentralGoogle Scholar
- Altincicek B, Vilcinskas A. Analysis of the immune-related transcriptome of a lophotrochozoan model, the marine annelid Platynereis dumerilii. Front Zool. 2007;4:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Kulakova M, Bakalenko N, Novikova E, Cook CE, Eliseeva E, Steinmetz PR, Kostyuchenko RP, Dondua A, Arendt D, Akam M, et al. Hox gene expression in larval development of the polychaetes Nereis virens and Platynereis dumerilii (Annelida, Lophotrochozoa). Dev Genes Evol. 2007;217(1):39–54.View ArticlePubMedGoogle Scholar
- Achim K, Pettit JB, Saraiva LR, Gavriouchkina D, Larsson T, Arendt D, Marioni JC. High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin. Nat Biotechnol. 2015;33(5):503–9.View ArticlePubMedGoogle Scholar
- Boeckmann B. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003;31(1):365–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40(Database issue):D290–301.View ArticlePubMedGoogle Scholar
- Consortium G. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25:5.Google Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Jäger M, Ott C-E, Grünhagen J, Hecht J, Schell H, Mundlos S, Duda GN, Robinson PN, Lienau J. Composite transcriptome assembly of RNA-seq data in a sheep model for delayed bone healing. BMC Genomics. 2011;12(1):1–12.View ArticleGoogle Scholar
- Theissinger K, Falckenhayn C, Blande D, Toljamo A, Gutekunst J, Makkonen J, Jussila J, Lyko F, Schrimpf A, Schulz R, et al. De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome. Mar Genomics. 2016.Google Scholar
- Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23(9):1061–7.View ArticlePubMedGoogle Scholar
- Cornman RS, Bennett AK, Murray KD, Evans JD, Elsik CG, Aronstein K. Transcriptome analysis of the honey bee fungal pathogen, Ascosphaera apis: implications for host pathogenesis. BMC Genomics. 2012;13(285):13.Google Scholar
- Li YI, Copley RR. Scaffolding low quality genomes using orthologous protein sequences. Bioinformatics. 2013;29(2):160–5.View ArticlePubMedGoogle Scholar
- Parra G, Bradnam K, Ning Z, Keane T, Korf I. Assessing the gene space in draft genomes. Nucleic Acids Res. 2009;37(1):289–97.View ArticlePubMedGoogle Scholar
- Simao FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3.View ArticleGoogle Scholar
- Janssen R, Gouar ML, Pechmann M, Poulin F, Bolognesi R, Schwager EE, Hopfen C, Colbourne JK, Budd GE, Brown SJ, et al. Conservation, loss, and redeployment of Wnt ligands in protostomes: implications for understanding the evolution of segment formation. BMC Evol Biol. 2010;10(374):21.Google Scholar
- Simakov O, Marletaz F, Cho SJ, Edsinger-Gonzales E, Havlak P, Hellsten U, Kuo DH, Larsson T, Lv J, Arendt D, et al. Insights into bilaterian evolution from three spiralian genomes. Nature. 2013;493(7433):526–31.View ArticlePubMedGoogle Scholar
- Dehal P, Boore JL. Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005;3(10):e314.View ArticlePubMedPubMed CentralGoogle Scholar
- Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman JA, Salamov A, Terry AY, Shapiro H, Lindquist E, Kapitonov VV, et al. Sea Anemone Genome Reveals Ancestral Eumetazoan Gene Repertoire and Genomic Organization. Science. 2007;317:9.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Levin M, Anavy L, Cole AG, Winter E, Mostov N, Khair S, Senderovich N, Kovalev E, Silver DH, Feder M, et al. The mid-developmental transition and the evolution of animal body plans. Nature. 2016;531(7596):637–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Berriman M, Haas BJ, LoVerde PT, Wilson RA, Dillon GP, Cerqueira GC, Mashiyama ST, Al-Lazikani B, Andrade LF, Ashton PD, et al. The genome of the blood fluke Schistosoma mansoni. Nature. 2009;460(7253):352–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Schistosoma japonicum Genome S, Functional Analysis C. The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature. 2009;460(7253):345–51.View ArticleGoogle Scholar
- Moroz LL, Edwards JR, Puthanveettil SV, Kohn AB, Ha T, Heyland A, Knudsen B, Sahni A, Yu F, Liu L, et al. Neuronal transcriptome of Aplysia: neuronal compartments and circuitry. Cell. 2006;127(7):1453–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Lambert JD, Chan XY, Spiecker B, Sweet HC. Characterizing the embryonic transcriptome of the snail Ilyanassa. Integr Comp Biol. 2010;50(5):768–77.View ArticlePubMedGoogle Scholar
- Henry JJ, Perry KJ, Fukui L, Alvi N. Differential localization of mRNAs during early development in the mollusc, Crepidula fornicata. Integr Comp Biol. 2010;50(5):720–33.View ArticlePubMedGoogle Scholar
- Lapraz F, Rawlinson KA, Girstmair J, Tomiczek B, Berger J, Jekely G, Telford MJ, Egger B. Put a tiger in your tank: the polyclad flatworm Maritigrella crozieri as a proposed model for evo-devo. Evodevo. 2013;4(29):16.Google Scholar
- Mehr S, Verdes A, DeSalle R, Sparks J, Pieribone V, Gruber DF. Transcriptome sequencing and annotation of the polychaete Hermodice carunculata (Annelida, Amphinomidae). BMC Genomics. 2015;16:445.View ArticlePubMedPubMed CentralGoogle Scholar
- Nyberg KG, Conte MA, Kostyun JL, Forde A, Bely AE. Transcriptome characterization via 454 pyrosequencing of the annelid Pristina leidyi, an emerging model for studying the evolution of regeneration. BMC Genomics. 2012;13(287):14.Google Scholar
- Andrade SCS, Novo M, Kawauchi GY, Worsaae K, Pleijel F, Giribet G, Rouse GW. Articulating “Archiannelids”: Phylogenomics and Annelid Relationships, with Emphasis on Meiofaunal Taxa. Mol Biol Evol. 2015;32(11):16.View ArticleGoogle Scholar
- Struck Torsten H, Golombek A, Weigert A, Franke Franziska A, Westheide W, Purschke G, Bleidorn C, Halanych Kenneth M. The Evolution of Annelids Reveals Two Adaptive Routes to the Interstitial Realm. Curr Biol. 2015;25(15):1993–9.View ArticlePubMedGoogle Scholar
- Liu MM, Davey JW, Jackson DJ, Blaxter ML, Davison A. A conserved set of maternal genes? Insights from a molluscan transcriptome. Int J Dev Biol. 2014;58(6–8):501–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang G, Fang X, Guo X, Li L, Luo R, Xu F, Yang P, Zhang L, Wang X, Qi H, et al. The oyster genome reveals stress adaptation and complexity of shell formation. Nature. 2012;490(7418):49–54.View ArticlePubMedGoogle Scholar
- Lee MT, Bonneau AR, Giraldez AJ. Zygotic genome activation during the maternal-to-zygotic transition. Annu Rev Cell Dev Biol. 2014;30:581–613.View ArticlePubMedPubMed CentralGoogle Scholar
- Hess HA, Roper JC, Grill SW, Koelle MR. RGS-7 completes a receptor-independent heterotrimeric G protein cycle to asymmetrically regulate mitotic spindle positioning in C. elegans. Cell. 2004;119(2):209–18.View ArticlePubMedGoogle Scholar
- Martin-McCaffrey L, Willard FS, Oliveira-dos-Santos AJ, Natale DR, Snow BE, Kimple RJ, Pajak A, Watson AJ, Dagnino L, Penninger JM, et al. RGS14 is a mitotic spindle protein essential from the first division of the mammalian zygote. Dev Cell. 2004;7(5):763–9.View ArticlePubMedGoogle Scholar
- Porter MY, Koelle MR. Insights into RGS protein function from studies in Caenorhabditis elegans. Prog Mol Biol Transl Sci. 2009;86:15–47.View ArticlePubMedGoogle Scholar
- Paps J, Xu F, Zhang G, Holland PW. Reinforcing the egg-timer: recruitment of novel lophotrochozoa homeobox genes to early and late development in the pacific oyster. Genome Biol Evol. 2015;7(3):677–88.View ArticlePubMedGoogle Scholar
- Gillis WQ, Bowerman BA, Schneider SQ. The evolution of protostome GATA factors: molecular phylogenetics, synteny, and intron/exon structure reveal orthologous relationships. BMC Evol Biol. 2008;8:112.View ArticlePubMedPubMed CentralGoogle Scholar
- Maduro MF. Endomesoderm specification in Caenorhabditis elegans and other nematodes. Bioessays. 2006;28(10):1010–22.View ArticlePubMedGoogle Scholar
- Maduro MF, Broitman-Maduro G, Choi H, Carranza F, Wu AC, Rifkin SA. MED GATA factors promote robust development of the C. elegans endoderm. Dev Biol. 2015;404(1):66–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Heyn P, Kircher M, Dahl A, Kelso J, Tomancak P, Kalinka AT, Neugebauer KM. The earliest transcribed zygotic genes are short, newly evolved, and different across species. Cell Rep. 2014;6(2):285–92.View ArticlePubMedGoogle Scholar
- Seaver EC. Variation in spiralian development: insights from polychaetes. Int J Dev Biol. 2014;58(6–8):457–67.View ArticlePubMedGoogle Scholar
- Park SJ, Komata M, Inoue F, Yamada K, Nakai K, Ohsugi M, Shirahige K. Inferring the choreography of parental genomes during fertilization from ultralarge-scale whole-transcriptome analysis. Genes Dev. 2013;27(24):2736–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Garrey J, Davis RE. Transcription in pronuclei and one- to four-cell embryos drives early development in a nematode. Curr Biol. 2014;24(2):124–33.View ArticlePubMedGoogle Scholar
- Xue Z, Huang K, Cai C, Cai L, Jiang CY, Feng Y, Liu Z, Zeng Q, Cheng L, Sun YE, et al. Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013;500(7464):593–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Poccia D, Wolff R, Kragh S, Williamson P. RNA synthesis in male pronuclei of the sea urchin. Biochim Biophys Acta. 1985;824(4):349–56.View ArticlePubMedGoogle Scholar
- Lee MT, Bonneau AR, Takacs CM, Bazzini AA, DiVito KR, Fleming ES, Giraldez AJ. Nanog, Pou5f1 and SoxB1 activate zygotic gene expression during the maternal-to-zygotic transition. Nature. 2013;503(7476):360–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Baugh LR, Hill AA, Slonim DK, Brown EL, Hunter CP. Composition and dynamics of the Caenorhabditis elegans early embryonic transcriptome. Development. 2003;130(5):889–900.View ArticlePubMedGoogle Scholar
- Ali-Murthy Z, Lott SE, Eisen MB, Kornberg TB. An essential role for zygotic expression in the pre-cellular Drosophila embryo. PLoS Genet. 2013;9(4):e1003428.View ArticlePubMedPubMed CentralGoogle Scholar
- Collart C, Owens ND, Bhaw-Rosun L, Cooper B, De Domenico E, Patrushev I, Sesay AK, Smith JN, Smith JC, Gilchrist MJ. High-resolution analysis of gene activity during the Xenopus mid-blastula transition. Development. 2014;141(9):1927–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Paranjpe SS, Jacobi UG, van Heeringen SJ, Veenstra GJ. A genome-wide survey of maternal and embryonic transcripts during Xenopus tropicalis development. BMC Genomics. 2013;14:762.View ArticlePubMedPubMed CentralGoogle Scholar
- Aanes H, Winata CL, Lin CH, Chen JP, Srinivasan KG, Lee SG, Lim AY, Hajan HS, Collas P, Bourque G, et al. Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition. Genome Res. 2011;21(8):1328–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Harvey SA, Sealy I, Kettleborough R, Fenyes F, White R, Stemple D, Smith JC. Identification of the zebrafish maternal and paternal transcriptomes. Development. 2013;140(13):2703–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolger A, Giorgi F. Trimmomatic: a flexible read trimming tool for illumina NGS data. http://www usadellab org/cms/index php 2014.Google Scholar
- Hannon G. fastx-Toolkit. In.: Cold Spring Harbor Laboratory. http://hannonlab.cshl.edu/fastx_toolkit//(27 Feb. 2014); 2010.
- Texas Advanced Computing Center (TACC), The University of Texas at Austin. https://portal.tacc.utexas.edu/tacc-citation.
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.View ArticlePubMedGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.View ArticlePubMedGoogle Scholar
- UniProt C. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43(Database issue):D204–12.Google Scholar
- Zhang Z, Wood WI. A profile hidden Markov model for signal peptides generated by HMMER. Bioinformatics. 2003;19(2):2.View ArticleGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42(Database issue):D199–205.View ArticlePubMedGoogle Scholar
- Li L, Stoeckert Jr CJ, Roos DS. OrthoMCL: Identification of Ortholog Groups for Eukaryotic Genomes. Genome Res. 2003;13:12.Google Scholar
- Nordberg H, Cantor M, Dusheyko S, Hua S, Poliakov A, Shabalov I, Smirnova T, Grigoriev IV, Dubchak I. The genome portal of the Department of Energy Joint Genome Institute: 2014 updates. Nucleic Acids Res. 2014;42(Database issue):D26–31.View ArticlePubMedGoogle Scholar
- Tribolium Genome Sequencing C, Richards S, Gibbs RA, Weinstock GM, Brown SJ, Denell R, Beeman RW, Gibbs R, Beeman RW, Brown SJ, et al. The genome of the model beetle and pest Tribolium castaneum. Nature. 2008;452(7190):949–55.View ArticleGoogle Scholar
- Attrill H, Falls K, Goodman JL, Millburn GH, Antonazzo G, Rey AJ, et al. FlyBase: establishing a Gene Group resource for Drosophila melanogaster. Nucleic Acids Res. 2016;44(D1):D786-92.Google Scholar
- Cameron RA, Samanta M, Yuan A, He D, Davidson E. SpBase: the sea urchin genome database and web site. Nucleic Acids Res. 2009;37(Database issue):D750–4.View ArticlePubMedGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Billis K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fitzgerald S, et al. Ensembl 2014. Nucleic Acids Res. 2014;42(Database issue):D749–55.View ArticlePubMedGoogle Scholar
- Van Dongen SM. Graph clustering by flow simulation. 2001.Google Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26(4):493–500.View ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky H. Primer3 on the WWW for General Users and for Biologist Programmers. In: Bioinformatics Methods and Protocols, vol. 132. 1999.View ArticleGoogle Scholar