Different disordered structures either in the nacreous or prismatic shell layers, suggesting potential roles in the regulation of shell formation. Our results have described the temporal expression and relative levels of RNA accumulation during larval shell development, increased the understanding of the molecular mechanisms, and the knowledge of biomineralization-related genes.ResultsGlobal gene expression analysisThe Raw and normalized fluorescence microarray data have been deposited in the GEO database under Accession Number GSE63824. A principal components analysis (PCA) on the entire probe set separated all 15 sample pools into five groups (Figure 1), which were relevant to the five larval developmental stages. The first and second components explained 70.93 of the variation in the entire data set. Stages of larval development was explained by the first principle components (PC1), which explained 48.99 of the variation, with low component loading for earlier stages and high component loading for later stages. A curved “horseshoe effect” along the second axis was apparent, which placed the trochophore and D-shaped stage larvae on the opposite side from fertilized eggs and the juvenile stage. Comparable results were obtained according to similarities in the expression patterns after a hierarchical clustering analysis (Additional file 1). These results indicate different gene expression patterns corresponding to different developmental stages. Meanwhile, the three biological replicate pools at each developmental stage shared a similar pattern, which would also increase the reliability and repeatability of this work.Transcriptional changes across larval stage transitionsA comparative analysis was performed using SAM software to identify changes in expression profiles between two consecutive developmental stages [30]. Among all 58,749 unigenes, the greatest number of differentially expressed genes was found between fertilized eggs and trochophore stage larvae, with 23,300 buy Ixazomib citrate significant genes (12,279 up-regulated and 11,021 down-regulated). However, the fewest genes changed between the D-shaped stage and trochophore stage, with 14,197 significant PubMed ID:https://www.ncbi.nlm.nih.gov/pubmed/25768400 genes (7,885 up-regulated and 6,312 down-regulated). The expression profiles of most genes changed significantly during larval development (Figure 2). Cellular Component Gen Ontology (GO) (Additional file 2) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (Additional file 3) terms were representedLiu et al. BMC Genomics (2015) 16:Page 3 ofFigure 1 3D score plot using the first 3PCs identified by principal component analysis on the entire larval gene probe set. Fifteen pools of samples were separated into five groups, relevant to the five larval developmental stages. Each stage point included a triplicate.significantly (q-value < 0.05) between the two consecutive stages.Comparison of fertilized eggs and trochophore stage larvaeThe genes that changed significantly between these two stages were related to GO:0005634//nuclear and GO:0005622//intracellular processes in a Cellular Components GO term analysis.Comparison of trochophore stage and D-shaped stage larvaeSeventeen KEGG pathways (q-value < 0.05) were significantly enriched between fertilized eggs and trochophore stage larvae. The majority was related to cellular progression (e.g., ko03030://DNA replication, ko04110://cell cycle, ko03420://nucleotide excision repair, ko03430:// mismatch repair, ko03440://homologous recombination).As.