Skip to main content
  • Research article
  • Open access
  • Published:

Cell cycle control, DNA damage repair, and apoptosis-related pathways control pre-ameloblasts differentiation during tooth development

Abstract

Background

Ameloblast differentiation is the most critical stepwise process in amelogenesis, and it is controlled by precise molecular events. To better understand the mechanism controlling pre-ameloblasts (PABs) differentiation into secretory ameloblasts (SABs), a more precise identification of molecules and signaling networks will elucidate the mechanisms governing enamel formation and lay a foundation for enamel regeneration.

Results

We analyzed transcriptional profiles of human PABs and SABs. From a total of 28,869 analyzed transcripts, we identified 923 differentially expressed genes (DEGs) with p < 0.05 and Fold-change > 2. Among the DEGs, 647 genes showed elevated expression in PABs compared to SABs. Notably, 38 DEGs displayed greater than eight-fold changes. Comparative analysis revealed that highly expressed genes in PABs were involved in cell cycle control, DNA damage repair and apoptosis, while highly expressed genes in SABs were related to cell adhesion and extracellular matrix. Moreover, coexpression network analysis uncovered two highly conserved sub-networks contributing to differentiation, containing transcription regulators (RUNX2, ETV1 and ETV5), solute carrier family members (SLC15A1 and SLC7A11), enamel matrix protein (MMP20), and a polymodal excitatory ion channel (TRPA1).

Conclusions

By combining comparative analysis and coexpression networks, this study provides novel biomarkers and research targets for ameloblast differentiation and the potential for their application in enamel regeneration.

Background

Enamel is one of the hardest mineralized tissues in vertebrates [1], and it is formed through extracellular matrix deposition by secretory ameloblasts (SABs) and the subsequent mineralization of organic matrix by mature ameloblasts (MABs) [2]. SABs and MABs are terminally differentiated cells derived from pre-ameloblasts (PABs), which are undifferentiated ameloblasts that are specified during odontogenesis [3]. PABs are elongated columnar dental epithelial cells characterized by nuclei localized near the stellate reticulum side and cytoplasm filled with organelles required for enamel protein synthesis and secretion.

Amelogenesis is a progressive differentiation process characterized by enamel formation and ameloblast differentiation, and it is regulated by various molecular and morphogenetic events. The transition of PAB to SAB is considered a hallmark of amelogenesis. Signals from the underlying dental mesenchymal compartment initiate this differentiation process by inducing PABs to SABs after basement membrane degradation during the bell stage of tooth development [4]. Differentiated SABs are tall, columnar, and polarized cells, and they can synthesize and secrete various enamel-specific proteins, such as amelogenin, ameloblastin, amelotin and enamelin [5]. These extracellular matrix proteins are essential for enamel formation and mineralization [6], and perturbation in ameloblast protein synthesis in humans results in malformed enamel [7]. After depositing full-thickness enamel matrix, ameloblasts further differentiate into short MABs to regulate enamel mineralization. Following tooth eruption, MABs undergo apoptosis [8]. Because there are no readily available ameloblast-lineage cells in human adults, they cannot restore or regenerate enamel.

Ameloblast differentiation is critical for enamel formation and enamel regeneration, thus, recent studies have focused on the molecular mechanisms of ameloblast differentiation [3, 9]. Roles of growth factors, enamel matrix proteins and transcription factors have been identified; for example, bone morphogenetic protein 2 (BMP2) and transforming growth factor-β1 (TGFβ1) were shown to induce ameloblast differentiation in vitro [10]. Studies on cultured tooth explants have also indicated that ameloblast differentiation is regulated by antagonistic effects between bone morphogenetic protein 4 (BMP4) and activin A [11]. In addition to these non-cell autonomous molecules, some secreted proteins, such as ameloblastin (AMBN), amelogenin (AMEL), enamelin (ENAM), and enamel proteinases such as matrix metalloproteinase 20 (MMP20) play key roles in full ameloblast differentiation [8]. The stage-specific expression of these genes is controlled by specific transcription factors. For example, the overexpression of Runt-related transcription factor (Runx2) and Distal-less homeobox 3 (Dlx3) result in the up-regulation of AMEL and ENAM mRNA levels [12]. These studies identified fundamental molecular events that control ameloblast differentiation and enamel regeneration. However, the differentiation from PAB to SAB is likely controlled by a precise molecular network like most cell differentiation processes. Therefore, further study of the transcriptional profile during ameloblast differentiation is important to identify novel biomarkers and research targets for the future investigation of tooth development and regeneration.

Here, we performed a comprehensive analysis of PABs and SABs transcriptional profiles from individual gene transcriptional levels to functional categories of enriched genes. This strategy revealed that the transition from PABs to SABs coincided with diverse DEGs involved in multiple pathways. Particularly, genes encoding growth factors, signaling molecules, transcription factors and their corresponding receptors, as well as membrane transport proteins were specifically expressed, and these molecules likely govern the biological processes essential for ameloblast differentiation and amelogenesis, including cell cycle, cell proliferation, cell adhesion, apoptosis and biomineral tissue development. Finally, we used bioinformatic approaches to explore significant coexpression networks and identify novel coexpression associations of central hubs in these networks.

Results

PAB and SAB morphology and gene marker characterization

To isolate PABs and SABs from human tooth buds, we performed rapid H&E staining to identify each cell type by location and morphology, then collected the cells by laser capture microdissection (LCM) (Fig. 1). In addition, we performed quantitative real-time PCR (qPCR) on micro-dissected cells (PABs and SABs) from human incisors to examine the expression of four genes known to be involved in ameloblast differentiation. As expected, SABs more highly expressed AMLX and AMBN compared to PABs, and PABs expressed higher levels of PCNA compared to SABs. PITX2 expression was up-regulated in PABs compared to SABs (Fig. 2a).

Fig. 1
figure 1

Identification of ameloblasts from a developing human incisor for laser capture microdissection. Left panel: H&E staining of a developing human incisor. Right panel: a Magnification of (a) in left panel. Pre-ameloblasts (PAB) were polarized inner enamel epithelial cells that directly contacted the basement membrane and were adjacent to polarized odontoblasts. b Magnification of (b) in left panel. Secretory ameloblasts (SAB) were identified as polarized epithelial cells in direct contact with the enamel matrix. (PAB: preameloblast, SAB: secretory ameloblast, POB: preodontoblast, SOB: secretory odontoblast, D: dentin, E: enamel, SR: Stellate reticulum, SI: Stratum intermedium, P: dental pulp)

Fig. 2
figure 2

a. Characterization of different stages of micro-dissected human dental epithelial cells. AMLX and AMBN expression levels were up-regulated as cells differentiated, while PCNA expression decreased in SABs compared to PABs. **p < 0.01. PITX2 expression was up-regulated in PABs compared to SABs. b. PCA clustering showed a clear divergence of PABs and SABs

Distinct transcriptional patterns in PABs and SABs

To generate the transcriptional profiles of PABs and SABs, we performed Affymetrix Human Genome 1.0ST Arrays on RNA from laser-captured cells. First, to obtain the whole transcriptional changes related to the differentiation from PABs to SABs, we applied principal component analysis (PCA) to the normalized microarray data. As anticipated, when the group number was set as 2, the transcriptional landscapes of the PABs and SABs separated into independent clusters, indicating distinct gene expression patterns for the two cell types (Fig. 2b). By comparing the probe intensity of 28,869 transcripts we identified a total of 923 DEGs with a threshold fold change (FC) > 2 and p < 0.05 between PABs and SABs. Among these DEGs, 647 genes were upregulated in PABs compared to SABs and 276 genes were upregulated in SABs compared to PABs. (See Additional file 1 for a full list of identified genes). Furthermore, 38 genes were dramatically changed (FC > 8) between PABs and SABs (Table 1).

Table 1 Most significant DEGs between PABs and SABs

KEGG pathway analysis of differentially expressed genes

We annotated the DEGs identified in the comparison of PABs to SABs using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. We assigned DEGs to 280 pathways, and we identified a total of 23 significantly enriched pathways (threshold p-value < 0.05) (Fig. 3a). We applied size-dependent enrichment analysis to map each gene cluster into pathways. Cluster maps of DNA replication, mismatch repair, cell cycle, p53 signaling, Notch signaling, small cell lung cancer and ECM-receptor interaction pathways are shown in Fig. 3b.

Fig. 3
figure 3

KEGG pathway analysis. a. KEGG pathway enrichment analysis was performed by Fisher exact test. Significantly enriched KEGG pathways (p < 0.05) are presented. For each KEGG pathway, the bar shows the fold-enrichment of the pathway. b. The expression intensity of genes in the selected KEGG pathway. The italic values on the right of gene names were the fold changes of corresponding genes in PABs compared to SABs

Up-regulated genes in PABs control cell cycle, DNA damage and apoptosis

Our microarray results revealed that the up-regulated genes in PABs compared to SABs are primarily involved in cell cycling, cell proliferation, apoptosis and DNA mismatch repair. Cyclin-dependent kinases (CDKs) have been implicated as involved in the cell cycle control system, including chromosome condensation, nuclear envelop breakdown and spindle assembly, by phosphorylating intracellular proteins [13, 14]. Genes encoding CDK family members, including CDC2, CDKL2, CDK6, and CABLES1, showed a higher level of expression in PABs compared to SABs. Since it is known that CDKs tightly bind to cyclins to form cyclin/CDK complexes, which in turn activate CDKs and trigger cell-cycle progression [15], we next investigated the corresponding cyclin family members. In accordance with the detected increase of CDKs, genes encoding cyclins (for instance, CCNA2, CCNB1, CCNB2, and CCNE2) and the cyclins interacted protein CCNBLIP1 (cyclin B1 interacting protein) were highly expressed in PABs compared to SABs. Interestingly, several negative cyclin/CDK complex regulators, WEE1, cyclin-dependent kinase inhibitor 3 (CDKN3) and CDK inhibitor (CKI), were also up-regulated in PABs compared to SABs. These inhibitors, WEE1 for example, phosphorylate the kinase active site to inhibit cyclin/CDK complex activity (Additional file 2: Table S3) [16, 17].

Another noteworthy feature in PABs is the DNA damage, which is an important factor in cell cycle progression. DNA damage initiates a signaling pathway, and the cell cycle control system can detect signals to arrest the cell cycle at late G1 or G2/M checkpoints [18]. We found that MDM2 expression was increased in PABs compared to SABs. MDM2 is a critical regulator of the cell cycle, as it ubiquitinates P53, targeting it for destruction by the proteasome [19]. Moreover, members of the DNA mismatch repair pathway, including MSH2, PMS2l5, RMI1 and EXO1, were enriched in PABs compared to SABs (Additional file 2: Table S3) [20, 21].

Notably, CASP6 and CASP8AP2 expression was increased in PABs compared to SABs (Additional file 2: Table S3), and they encode the apoptosis-related executioner caspase (Caspase 6) and initiator caspase (Caspase 8), respectively [22].

Taken together, the significantly up-regulated genes in PABs will work coordinately to form a sophisticated network that tightly controls the cell cycle, DNA damage and apoptosis, altogether contributing to ensuring that PABs remain in a proper state to support differentiation.

SABs are enriched for genes associated with extracellular matrix and adhesion

Amelogenin, ameloblastin, enamelin and amelotin are four EMPs, which are highly expressed during the secretory and early maturation stages of amelogenesis [23]. They are subsequently degraded step-wise from the matrix by MMP20 and kallikrein 4 (KLK4) during the ameloblast secretory stage and maturation stage, respectively [2]. Among them, amelogenin is the most abundant protein in enamel matrix and is continuously expressed throughout amelogenesis. Ameloblastin is the second most abundant EMP, and it localizes between enamel rods. Enamelin is the largest enamel matrix protein of developing teeth and comprises approximately 5 % of total EMPs. As expected, genes encoding AMBN, ENAM, AMTN and MMP20 were more highly expressed in SAB cells compared to PAB cells [6] (Additional file 2: Table S3). However, AMLX and KLK4 expression did not significantly differ between PABs and SABs. Other genes encoding proteins involved in biomineral tissue development (DSPP, DMP1, PHEX, ALPL and MMP16) [2427] were also more highly expressed in SABs compared to PABs (Additional file 2: Table S3).

In addition to genes encoding secreted proteins, SABs were enriched in five genes encoding laminins, including LAMC1, LAMC2, LAMB3, LAMA3 and LAMA2 (Additional file 2: Table S3). Laminins are primary basal lamina components, and they play key roles in differentiation, migration and adhesion [28]. Laminin alpha5 (LAMA5)-null mice have been shown to display aberrant cusp formation and dental epithelium proliferation [29].

Apart from adhering to the extracellular matrix, cell-cell adhesions also play an indispensable role in stabilizing the cellular and tissue structure. Cell-cell adhesion sites are strengthened by the connected underlying cytoskeleton, which then forms a network across several cells. In vertebrates, cell-cell adhesion is primarily mediated by the cadherin superfamily, which are important for cell rearrangements during tissue morphogenesis [30]. We found that cadherin superfamily members CDH8 and CDH13 were enriched in SABs compared to PABs (Additional file 2: Table S3). In addition, adhesion molecule with Ig-like domain 2 (AMIGO2) and IGJ expression were both elevated in SABs compared to PABs (Additional file 2: Table S3) [31, 32].

Differential expression of transcription factors and signaling pathways between PABs and SABs

By comparing PAB and SAB transcriptional profiles, we found that several growth factors, signaling molecules and transcription factors were differentially expressed between the two cell types. The transforming growth factor family members transforming growth factor β2 (TGFβ2) was more highly expressed in PABs compared to SABs, while bone morphogenetic protein 8A (BMP8A) was more highly expressed in SABs compared to PABs. Several other growth factors were differentially expressed between the cell types, suggesting that they may be critical regulators in ameloblast differentiation. One fibroblast growth factor family member (FGF20) was increased in SABs, while fibroblast growth factor 9 (FGF9) and nuclear factor I b (NFIB) expression were decreased in SABs compared to PABs. Interestingly, several negative regulators of the fibroblast growth factor receptor signaling pathway (SPRY1, SPRY2, SPRY4, THBS1 and SULF1) [3336] were enriched in PABs compared to SABs (Additional file 2: Table S3).

Notch signaling and MAPK signaling contribute to cell differentiation during tissue development [37, 38]. In this study, we found that Notch 2 and MAPK6 expression was increased in SABs compared to PABs. However, two other MAPK members (MAPK10 and MAPK14) were enriched in PABs compared to SABs. Moreover, the distalless related homeodomain proteins (Dlx5) and transcription factor AP-2 alpha (TFAP2A) were highly expressed in PABs compared to SABs (Additional file 2: Table S3). Two G protein-coupled receptor pathway members (GPR110 and GPSM2) were more highly expressed in PABs. In addition, three interleukin genes were differentially expressed between PABs and SABs. IL17 receptor (IL17RD) expression increased, while IL1 receptor (IL1R1) and IL33 expression decreased in PABs compared to SABs (Additional file 2: Table S3).

Transporters and ion channels are differentially expressed between PABs and SABs

Transcellular Ca2+ transport plays an important part in the formation of calcium hydroxyapatite during amelogenesis. Previous studies have shown that stromal interaction molecule 1 (STIM1) and solute carrier family 24, member 4 (SLC24A4) were critical ion transport systems in enamel maturation [39]. We found that in addition to SLC24A4, nine other solute carrier family members (SLC15A1, SLC7A11, SLC41A2, SLC13A5, SLC4A4, SLC17A5, SLC25A33 and SLC45A4) were enriched in SABs compared to PABs. In contrast, six family members (SLC26A4, SLC4A7, SLC35A3, SLC27A2, SLC16A10 and SLC35F1) were increased in PABs compared to SABs (Additional file 2: Table S3).

Channels are another major class of membrane transport proteins [40]. We found that three genes encoding transient receptor potential (TRP) channels, including TRPM7, TRPA1 and TRPC1, were differentially expressed between the two cell types. TRPM7 is essential for many cellular processes, including proliferation, survival, differentiation, growth, and migration [41]. Its expression increased in SABs compared to PABs, while TRPA1 and TRPC1 were more highly expressed in PABs. Strikingly, TRPA1 expression was more than 16 fold higher in PABs than SABs (Additional file 2: Table S3).

Other ion transport-associated genes that were differentially expressed include members of Na-K-ATPase (ATP1B1) and V-ATPase (ATP6V0D1 and ATP6V0A1), which were all increased in SABs compared to PABs (Additional file 2: Table S3).

Genome-wide coexpression network construction identified ameloblast development-associated modules

Different cell types synthesize different sets of proteins, so, to a large extent, gene expression defines the cell differentiation program. However, rather than functioning as completely independent components, genes interface in intricate ways to create a network. In particular, signaling networks precisely regulate diverse cellular events, including cell cycle, morphogenesis, proliferation, cell adhesion, and programmed cell death [42].

To further explore specific networks that orchestrate ameloblast differentiation and identify the gene coexpression associations of DEGs in these networks, we used weighted gene coexpression network analysis (WGCNA) to identify modules containing highly coexpressed genes. The coexpression network was based on the correlation values between a pair of genes across microarray data. As shown in Fig. 4a, we identified a total of 33 conserved modules and assigned them with unique colors. To determine how modules correlated with the differentiation of preameloblasts to secretory ameloblasts, we defined the gene signature measure (GS) as a transformation from p-values for testing differential gene expression between PABs and SABs (details in Methods). In particular, brown and blue modules (named by color) were enriched with high GS values (Fig. 4b). After mapping the DEGs to the coexpression modules, we found that brown and blue modules contained a significant number of DEGs (413 in brown and 252 in blue), as shown in Fig. 4a (second color bar). These data indicate that these two modules were highly related to differential gene expression in ameloblast differentiation. We also visualized the global coexpression network in brown and blue modules with a proper adjacency threshold (Fig. 5), and the detailed interaction list is in Additional file 3: Table S4 and S5.

Fig. 4
figure 4

Gene coexpression modules and assigned module colors. a. Clustering dendrogram of genes with the x-axis corresponding to each gene and the y-axis representing the dissimilarity based on topological overlap of the whole genome throughout ameloblast differentiation. Top color bar: modules were assigned colors corresponding to different dendrogram branches detected by a cluster algorithm. Second color bar: each gene was annotated with DEG information such that significantly up-regulated genes in PABs are marked with red bands and genes down-regulated in PABs marked with green bands. b. Boxplot depicting gene significance defined as -log10(p-value) of each module. The left legend indicates module colors and the number of genes (stars with a number sign) in the corresponding module (for example, 1527 genes in the blue module)

Fig. 5
figure 5

Gene coexpression networks of modules. Global network visualization of the top two modules “brown” (a) and “blue” (b) that are highly correlated with gene significance. Each node corresponds to a gene, and each edge joining two nodes indicates the connection determined by the adjacency. The nodes in boxes represent DEGs and circles indicate genes without significant differential expression between PABs and SABs. Nodes are also marked with colors to indicate their representative GO biological processes, and grey nodes are those not included in either of the GO terms. The left legend in each network shows colors and corresponding GO terms

Brown and blue module hub genes

To investigate the biological functions of the WGCNA modules, we selected the highly connected “hub genes” within each module. Due to their high degrees of connectivity, intramodular hub genes were centrally located in coexpression networks and may play an important role in their module’s cellular function and signaling pathways [43, 44]. We defined the module membership measure (MM) as the correlation between the gene expression and module eigengene. We selected the top 20 genes with high MM (>0.9) and GS values as reprehensive hub genes for the brown and blue modules (Fig. 6a and b) [45].

Fig. 6
figure 6

Hub gene selection and gene ontology analysis of brown and blue modules. a and b are the scatterplots between module membership measure (x-axis) and the gene significance of the “brown” and “blue” module. The top 20 genes with both the highest gene significance and membership measure larger than 0.9 were chosen as “hub genes”, and they are annotated with their gene symbols. The grey dashed line in the plot is the threshold for choosing significantly expressed genes, and the threshold value is –log10(0.05). c and d are the hub gene network visualizations for the “brown” and “blue” modules. The colors of the hub genes and their directly connected genes indicate the GO terms from Fig. 5 (a) and (b). e. Barplots represent enriched GO terms of all genes included in the “brown” (left with brown color) and “blue” (right with blue color) networks, and –log2(p-value) represents the relative enrichment of each GO term

In the brown module, there was a significant correlation (cor = 0.50, p-value = 1.50 × 10−92) between the intramodular connectivity and the gene significance, indicating the potential gene expression pattern in the brown module (Fig. 6a). Interestingly, RUNX2 was one of the highly connected genes in the brown module (MM = 0.93), which is consistent with a previous study that showed RUNX2 was a key regulator in ameloblast differentiation by regulating amelogenin and enamelin expression [12]. More importantly, we detected RUNX2 directly connected to 18 genes, including EXO1 (adjacency = 0.420), SULF2 (adjacency = 0.422) and ETV1 (adjacency = 0.423) (Fig. 6c and Additional file 3: Table S4). These genes may be RUNX2 partners during ameloblast differentiation. In addition, we identified other genes encoding RUNX proteins (RUNX1 and RUNX3) in the brown module. The hub genes RUNX1 (MM = 0.99) and RUNX3 (MM = 0.95) directly connected to 21 and 8 genes in the brown module, respectively (Fig. 6c and Additional file 3: Table S4). This was consistent with previous reports that they act as master developmental regulators [46] and indicated their function and partners in ameloblast differentiation. Another highly connected gene MMP20 (MM = 0.94) was significantly increased (FC = 25) in SABs and was associated with anatomical structure development. This result was consistent with a recent study that MMP20 facilitates ameloblast movement by cleaving ameloblast cell-cell contacts [47]. FGF9 (MM = 0.97) was also located in the center of the coexpression network (Fig. 6c), which suggests that they also play an important role in ameloblast differentiation. This was in line with their role in transcription regulation throughout development [48].

We also detected another group of hub genes including TRPA1, ETV1, ASPM, SPC25 and NUSAP1 (Additional file 3: Table S6). Although these genes have not been previously associated with ameloblast differentiation, their functions have been demonstrated in previous studies. TRPA1 is a stress sensor [49], and ETV1 is an ETS transcription factor required for retinal neuron differentiation [50]. Their appearance as hub genes in the brown module and up-regulation (FC > 8) in PABs highly suggest that they play critical functions in ameloblast differentiation. Notably, TRPA1 and ETV1 directly interacted with each other (adjacency = 0.428), and both of them had direct interactions with three other mitosis-associated hub genes (ASPM, SPC25 and NUSAP1) (Fig. 6c and Additional file 3: Table S4 and S6) [5153]. It is therefore reasonable to infer that TRPA1 and ETV1 influence ameloblast differentiation by regulating cell division. Strikingly, ETV1 also directly interacted with RUNX2, suggesting that ETV1 plays pleiotropic roles in ameloblast differentiation.

Moreover, ETV5 was a hub gene with the highest MM (0.97), but it was not directly connected with ETV1. Rather, it was surrounded by another hub gene in the brown module (FAM111B) with an adjacency of 0.42 (Fig. 6c and Additional file 3: Table S4). These data suggest that the two ETS transcription factors are not regulated in the same fashion and may therefore perform distinct physiological functions. This is consistent with the recent discovery that ETV4 and ETV5 often perform similar functions during morphogenesis, whereas ETV1 has distinct roles [54].

Similar to the brown module, we observed that the intramodular connectivity and GS were highly correlated in the blue module (cor = 0.49, p = 1.00 × 10−90), as shown in Fig. 6b. SLC15A1, SLC7A11 and SPHKAP were hubs in the blue module with a module membership of 0.96, 0.98 and 0.97, respectively. We found that SLC15A1 directly interacted with 27 genes and SLC7A11 directly connected with 24 genes (Fig. 6d and Additional file 3: Table S5). Interestingly, SPHKAP showed a coexpression relationship with both SLC15A1 (adjacency = 0.414) and SLC7A11 (adjacency = 0.414) (Fig. 6d and Additional file 3: Table S5). SPHKAP is annotated as an SPHK1-interacting protein. However, its role is still unclear. The interaction between SPHKAP and solute carrier family members (SLC15A1 and SLC7A11) suggests that it may function in solute transport during ameloblast differentiation.

LAMC2 showed a high MM of 0.97 in the blue module, suggesting that it has an indispensable function (Additional file 3: Table S7). By performing GO analysis, we found that it was involved in tissue development (Fig. 6d), which could suggest that laminin plays an important role in ameloblast differentiation. Though the exact function of laminin in ameloblast differentiation has not been reported, its function as an extracellular matrix component suggests that it is involved in cell proliferation and differentiation during tooth morphogenesis [55].

As an epithelial-to-mesenchymal (EMT) transcription factor, zinc finger E-box binding protein 1 (ZEB1) has been implicated in neural crest cell biology [56]. Though ZEB1 was not detected as a hub gene in the blue module, it directly connected with the SLC15A1 hub gene with an adjacency of 0.41, and it was surrounded by DEGs in the blue module hub gene network (Fig. 6d), indicating its potential role in ameloblast differentiation.

Interestingly, we also identified two BMP family members in the two coexpression modules. However, BMP6 (MM = 0.95) localized to the brown module (Additional file 3: Table S4 and S6), and BMP8A (MM = 0.95) localized to the blue module (Additional file 3: Table S5 and S7), indicating an independent regulatory function in ameloblast differentiation.

Gene Ontology (GO) analysis reveals gene association and biological processes essential to ameloblast differentiation

To functionally characterize genes in the brown and blue modules, we performed GO analysis. Genes associated with cellular component organization or biogenesis, cell division and proliferation were significantly enriched in the brown module, whereas genes related to ion binding, cell adhesion, and tissue development were enriched in the blue module (Fig. 6e), which indicated a striking difference in the functions of the two modules. The brown module closely associated with cell division and proliferation, and almost all of the hub genes in the brown module were up-regulated in PABs compared to SABs. While the blue module was likely related to cell adhesion and tissue development, its hub genes were upregulated in SABs compared to PABs. Collectively, these results suggested that the blue and brown modules play critical and distinct roles in ameloblast differentiation. We also identified several genes with high MM but low GS in both modules (Fig. 6a and b), suggesting that they were closely related to certain functions of the corresponding modules but that they not serve as DEGs. By separately plotting the genes within each GO term (Additional file 4), we found that the genes with high MM and low GS were functionally annotated. For example, those genes in the brown module were annotated as being involved in response to stimulus, while those in the blue module were annotated as functioning in cell differentiation. We speculated that these genes may act by interacting with hub genes or as indispensable components in biological pathways, and they should be examined as potential players in ameloblast differentiation.

qPCR validation of microarray results

To validate microarray results, we confirmed the transcriptional levels of 12 genes in PABs and SABs by qPCR, including seven DEGs and AMLX, KRT14, RUNX2, KLK4 and STAT2 (Table 2). Most of the genes exhibited differential expression in qPCR that was consistent with the microarray data, indicating good concordance of both methods.

Table 2 Results of microarray and qPCR experiments

Discussion

Ameloblasts are a tooth-specific cell type, and they secrete extracellular matrix proteins and deposit teeth enamel [6]. Full ameloblast differentiation occurs by three distinct processes: the presecretory, secretory and maturation stages [4]. The sequential differentiation of ameloblast lineage cells requires signals from the underlying mesenchymal compartment and dentinal matrix. Ameloblasts differentiate in a niche composed of secretory proteins, enzymes, signaling molecules and other factors. Ameloblast differentiation is defined by the shift from preameloblast to secretory ameloblast, and it can significantly influence enamel matrix protein secretion and mineralization initiation. This process is regulated by a complex signaling network composed of signaling molecules, their receptors and transcriptional regulators [3]. However, the field lacks a full understanding of the transcriptional program underlying ameloblast differentiation. In this study, we aimed to identify molecular events governing the shift from PABs to SABs and gene coexpression networks that control ameloblast maturation.

Comparison of transcriptional profiles of PAB and SAB together with KEGG pathway analysis revealed that PAB cells express genes that allow them to control cell cycle, repair DNA mismatch and eliminate abnormal, misplaced, or nonfunctional cells by apoptosis, which are necessary quality-control processes for PAB differentiation. In contrast, pathways enriched in SAB were primarily related to extracellular matrix and cell adhesion. It is widely accepted that SABs secrete EMPs into the enamel matrix, which is primarily composed of three structural proteins (AMELX, AMBN and ENAM). These proteins are degraded step-wise from the matrix by MMP20 and KLK4 at ameloblast secretory stage and maturation stage, respectively [2]. In this study, we identified the increased expression of four secretory ameloblast-specific genes (AMBN, AMTN, ENAM and MMP20) in SABs by microarray, which is consistent with previous studies and further validated by our LCM and microarray approach. Although microarray failed to detect AMELX enrichment in SABs, we confirmed its increased expression in SABs by real time PCR. This discrepancy may be due to the detection limit of microarrays, as fold-changes observed by microarray are usually lower than those observed by real-time PCR [57]. In addition, we identified five other genes (DSPP, DMP1, PHEX, ALPL and MMP16) related to biomineral tissue development enriched in SABs, suggesting that these genes may regulate ameloblast maturation and enamel formation.

In addition to genes encoding extracellular matrix proteins, we also noted the enrichment of cadherin and laminin genes in SABs. Cadherins control the selective cell sorting, and the upregulation and downregulation of specific cadherins correlate with specific steps in embryonic development [58]. Previous studies have shown increased expression of E-cadherin and N-cadherin in PABs and SABs, respectively [47]. We found that cadherin superfamily genes (CDH8, CDH13 and CDHR3) and laminin genes (LAMC1, LAMC2, LAMB3, LAMA3 and LAMA2) were more highly expressed in SABs compared to PABs, suggesting that they may regulate secretory ameloblast cell rearrangement and morphogenesis. However, this is the first implication that they may function in amelogenesis; thus, their precise function in ameloblast differentiation require further investigation.

Some signaling pathways are critical for animal development, including receptor tyrosine kinase (RTK), transforming growth factor family, Wnt, Hedgehog and Notch. In addition, JAK/STAT, nuclear hormone receptor, and G-protein-coupled receptor pathways are important in some developmental processes [42, 59]. Terminal ameloblast differentiation is regulated by the interactions between the epithelium and mesenchyme. PAB differentiation can be influenced by odontoblasts signals, such as TGFβ1, BMP2 and BMP4 [3, 60]. Our study provides proof that members of growth factors (BMP8A, TGFβ2, NFIB, FGF9 and FGF20), signaling pathways (MAPK6, MAPK10, MAPK14 and Notch2), G protein coupled receptor pathways (GPR110 and GPSM2) and transcription factors (Dlx5 and TFAP2A) participate in the differentiation of PABs to SABs.

Differentially expressed genes and KEGG pathway analysis revealed key genes and metabolic characteristics that may regulate PAB differentiation to SABs. Further coexpression network analysis allowed us to identify modules of highly correlated genes from whole transcriptome profiles and identify hub genes and their expression associations during ameloblast differentiation. Moreover, a functional enrichment analysis of the brown and blue modules revealed basic functional pathways preserved in these modules.

Conclusions

In this study, we found that the PAB to SAB differentiation relies on a highly regulated network of interactions between conserved signal transduction pathways, including members of the BMP/TGF-β, Notch, MAPK pathways to coordinate all aspects of ameloblasts in intracellular processes and their social contexts. Specifically, expression of genes associated with cell cycle control, DNA damage repair, and apoptosis pathways regulate pre-ameloblast maturation during tooth development. SAB cells are regulated by several signaling pathways that control enamel matrix protein secretion and cell adhesion, which are critical for enamel formation and cell-cell interactions. Furthermore, the application of bioinformatic analysis allowed us to explore potential key genes and gene-associations involved in ameloblast differentiation. These findings will aid in the design of new strategies to promote ameloblast differentiation in tooth regeneration and tissue engineering. However, further experiments are required to confirm the roles of these genes and their interacting genes. Although transcriptional controls are paramount for most gene expression, these genes and pathways can also be regulated by a complex array of genetic (e.g., RNA splicing and selecting), epigenetic (e.g., histone modification), non-coding RNA (e.g., microRNA and long non-coding RNA) and exogenous signaling factors that serve to guide cell fate and behavior during development and differentiation.

Methods

Ethics statement

All human tissues were collected from legally aborted fetuses at West China Women and Children’s Hospital under the approved guidelines of Sichuan University. The written informed consent of all human subjects was obtained. The study and the consent procedure were approved by the Ethical Committees of West China School of Stomatology, Sichuan University and State Key Laboratory of Oral Diseases.

Sample collection

Human tooth buds (18-22 weeks) were obtained from fetal cadaver tissue within three hours after legal abortion under the guidelines of the West China School of Stomatology, Sichuan University Committee on Human Research. Teeth were dissected from the mandibles under a laminar flow hood, embedded in OCT compound (Sakura Finetek, Torrance, CA, USA), and cryosectioned at 10-μm thickness. These sections were used for LCM.

Laser capture microdissection (LCM)

A PALM MicroBeam system (Carl Zeiss Micro Imaging, Inc. Thornwood, NY, USA) laser capture microscope was used for cell dissection [61, 62]. All reagents used for LCM and RNA purification were prepared in RNase-free water. Cryosectioned human tooth organs were mounted on polyethylene naphthalate (PEN) foil glass slides. After rapid hematoxylin-eosin (H&E) staining, slides were air-dried, and dental epithelial cells were identified by location and morphology. PABs and SABs were separately dissected and catapulted to Adhesive Cap (Carl Zeiss MicroImaging Gmbh, München, Germany) containing ß-Mercaptoethanol-supplemented RLT buffer (RNeasy Plus Micro Kit, RNeasy Lysis Buffer, Qiagen, Valencia, CA, USA). The caps containing buffer and cells were incubated upside-down at room temperature (RT) for 30 min and then briefly centrifuged to collect the buffer and cell complexes. RNA purification was performed with an RNeasy Plus Micro Kit (Qiagen), and reverse transcription was performed with SuperScript™ III First-Strand Synthesis System (Invitrogen Corporation, Carlsbad, CA, USA) following the manufacturer’s instructions.

Microarray procedures

Affymetrix Human Genome 1.0ST Arrays (Affymetrix, Santa Clara, CA, USA) covering 28,869 transcripts in the human genome were employed for transcriptional profiling and performed by Affymetrix. Three to four pre-ameloblast and secretory ameloblast samples collected by LCM were pooled for each array. RNA samples were extracted and purified by RNeasy Mini kits (Qiagen and Rnase-Free DNase Set) as described previously [63]. A Nanodrop ND 1000 spectrophotometer (Thermo Fisher Scientific, Pittsburgh, PA, USA) was used to determine total RNA concentrations, and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used to evaluate RNA quality [64] Both RNA samples were of acceptable quality (Average RIN of 9.0 and 8.0 for PABs and SABs, respectively), and there was no significant difference between their qualities. For each replicate, 100 ng of total RNA was amplified and labeled using the Affymetrix Whole-Transcript (WT) Sense Target Labeling Protocol without rRNA reduction. Array hybridization, washing, and scanning were performed according to the protocol described in WT Sense Target Labeling Assay Manual.

Microarray expression analysis

The R/Bioconductor Package “Affy” was used to perform gene expression quartile normalization to adjust the marginal distribution of each sample. We used the “limma” package to identify DEGs [65]. Genes with fold-change greater than 2 and p < 0.05 between PABs and SABs were considered statistically significant.

KEGG pathway analysis

The complete list of human pathways was retrieved from the KEGG Pathway Database (http://www.genome.jp/kegg/). A functional gene enrichment analysis was performed based on the KEGG pathways [66]. We applied Fisher’s exact test to determine whether a set of DEGs was selectively enriched in a pathway. KEGG pathways containing at least two DEGs with p-value less than 0.05 were considered significant.

Weighted gene coexpression network analysis (WGCNA)

The WGCNA package in R was used to construct unsigned coexpression networks [45, 67]. After filtering raw data without gene function annotation, 22,247 probes were used to construct the pair-wise Pearson’s correlation coefficient matrix. Next, the correlation matrix was transformed to an adjacency matrix by a soft power adjacency function provided by Bin Zhang and Steven Horvath [45]. The scale-free topology criterion method was used to choose the power parameter β. Specifically, β = 16 was used as a balance between the maximization of scale-free topology model fit (topological scale R square was 0.73) and the number of connections. The adjacency value for each pair of gene was defined as raising the absolute value of their Pearson’s correlation coefficient to the power β, and the adjacency values ranged from 0 to 1. By using the adjacency matrix, the topological overlap measure, which is a similarity measurement, was calculated for each pair of genes. The topological overlap-based dissimilarity was then used for clustering and dynamic tree cutting with the minimum module size of 30 and merged tree-cut height of 0.25 [68]. A group of genes that were tightly gathered into the same cluster was considered a coexpression module. The blockwiseModules function in WGCNA R package was used to automatically select modules in PABs and SABs, and the maximum block size was set to 10,000 due to the limitation of computer memory.

The T-statistics-based p-value yielded from “limma” R package [65] was used for each gene as a criterion to choose proper modules associated with ameloblast differentiation. The gene significance (GS) should be a positive value calculated by the formula: GS = −log10 (p-value). The top two modules with the highest GS values (equal to the enrichment of more DEGs) were selected and identified as the brown and blue modules. To make a representative visualization of the correlation network, only the strongest linkages were drawn (correlation was greater than 0.95). Correlation networks of the brown module and blue module were generated in Cytoscape with adjacency thresholds of 0.42 and 0.41, respectively [69].

Hub gene selection and gene ontology analysis

The most highly connected genes, which were also well-known hub genes, were used to represent the expression profiles and biological characters of the entire module. The module membership measure (MM) for a gene was defined as the correlation between the eigengene of the module and the gene expression; therefore, the MM was also named the module eigengene based connectivity (kME) [45]. As we used the p-value based GS for module selection, the top 20 genes with both highest GS and MM (0.9 or greater) were defined as intramodular hub genes. Subnetworks consisting of intramodular hub genes and their directly connected genes were also plotted.

Functional enrichment analysis for each selected module was conducted by the R/Bioconductor package clusterProfiler [70]. The background gene set was the total list of genes used for coexpression network construction. The threshold level to determine the significance of all the GO enrichment was set as p < 0.05.

Microarray data accession

All raw microarray data were submitted to the NCBI Gene Expression Omnibus database with the accession number GSE59214, which can now be viewed through the link http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE59214

Quantitative real-time PCR

For cell-type confirmation after LCM, AMLX, AMBN, PITX2, and PCNA expression was examined by quantitative real-time PCR in an ABI 7500 system (Applied Biosystems, Foster City, CA, USA). For microarray data validation, we examined the transcriptional levels of 12 genes in PABs and SABs by qPCR. RNA was extracted with the same method as described for microarray. The relative abundance of mRNA transcripts was quantified relative to GAPDH levels using the ∆∆CT method [71]. Primers and probe sets including an endogenous GAPDH control were purchased from Applied Biosystems (Applied Biosystems, Foster City, CA, USA). The corresponding arithmetic formulas used are as follows: ΔCTcondition = CTtarget gene − CTendogeneous control and ΔΔCT = ΔCTcondition1 − ΔCTcondition2. For example, we compared the gene expression fold change between different cell types. Take amelogenin expression in PABs and SABs as an example, ΔCTPAB = CTamelogenin at PAB-CTGAPDH at PAB, ΔCTSAB = CTamelogenin at SAB-CTGAPDH at SAB, ΔΔCT = ΔCTSAB-ΔCTPAB. The fold change of SAB to PAB is calculated as FoldSAB to PAB = 2(−ΔΔCT) (amelogenin expression at PAB was used as a baseline value of 1). Data were compared by one-way ANOVA followed by the post hoc Tukey’s test.

References

  1. Wald T, Osickova A, Sulc M, Benada O, Semeradtova A, Rezabkova L, et al. Intrinsically disordered enamel matrix protein ameloblastin forms ribbon-like supramolecular structures via an N-terminal segment encoded by exon 5. J Biol Chem. 2013;288(31):22333–45.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  2. Gibson CW. The Amelogenin proteins and enamel development in humans and mice. J Oral Biosci. 2011;53(3):248–56.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Bei M. Molecular genetics of ameloblast cell lineage. J Exp Zool B Mol Dev Evol. 2009;312B(5):437–44.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Zeichner-David M, Diekwisch T, Fincham A, Lau E, MacDougall M, Moradian-Oldak J, et al. Control of ameloblast differentiation. Int J Dev Biol. 1995;39(1):69–92.

    CAS  PubMed  Google Scholar 

  5. Bartlett JD, Smith CE. Modulation of cell-cell junctional complexes by matrix metalloproteinases. J Dent Res. 2013;92(1):10–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Bartlett JD. Dental enamel development: proteinases and their enamel matrix substrates. ISRN Dent. 2013;2013:684607.

    PubMed Central  PubMed  Google Scholar 

  7. Wright JT. The molecular etiologies and associated phenotypes of amelogenesis imperfecta. Am J Med Genet A. 2006;140(23):2547–55.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Sierant ML, Bartlett JD. Stress response pathways in ameloblasts: implications for amelogenesis and dental fluorosis. Cells. 2012;B1(3):631–45.

    Article  Google Scholar 

  9. Li Z, Yu M, Tian W. An inductive signalling network regulates mammalian tooth morphogenesis with implications for tooth regeneration. Cell Prolif. 2013;46(5):501–8.

    CAS  PubMed  Google Scholar 

  10. Coin R, Haikel Y, Ruch JV. Effects of apatite, transforming growth factor beta-1, bone morphogenetic protein-2 and interleukin-7 on ameloblast differentiation in vitro. Eur J Oral Sci. 1999;107(6):487–95.

    Article  CAS  PubMed  Google Scholar 

  11. Wang XP, Suomalainen M, Jorgez CJ, Matzuk MM, Werner S, Thesleff I. Follistatin regulates enamel patterning in mouse incisors by asymmetrically inhibiting BMP signaling and ameloblast differentiation. Dev Cell. 2004;7(5):719–30.

    Article  CAS  PubMed  Google Scholar 

  12. Athanassiou-Papaefthymiou M, Kim D, Harbron L, Papagerakis S, Schnell S, Harada H, et al. Molecular and circadian controls of ameloblasts. Eur J Oral Sci. 2011;119 Suppl 1:35–40.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Murray AW, Kirschner MW. Cyclin synthesis drives the early embryonic cell cycle. Nature. 1989;339(6222):275–80.

    Article  CAS  PubMed  Google Scholar 

  14. Bruyere C, Meijer L. Targeting cyclin-dependent kinases in anti-neoplastic therapy. Curr Opin Cell Biol. 2013;25(6):772–9.

    Article  CAS  PubMed  Google Scholar 

  15. Ear PH, Booth MJ, Abd-Rabbo D, Kowarzyk Moreno J, Hall C, Chen D, et al. Dissection of Cdk1-cyclin complexes in vivo. Proc Natl Acad Sci U S A. 2013;110(39):15716–21.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Do K, Doroshow JH, Kummar S. Wee1 kinase as a target for cancer therapy. Cell cycle. 2013;12(19):3159–64.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Lim S, Kaldis P. Cdks, cyclins and CKIs: roles beyond cell cycle regulation. Development. 2013;140(15):3079–93.

    Article  CAS  PubMed  Google Scholar 

  18. Harrison JC, Haber JE. Surviving the breakup: the DNA damage checkpoint. Annu Rev Genet. 2006;40:209–35.

    Article  CAS  PubMed  Google Scholar 

  19. Tollini LA, Jin A, Park J, Zhang Y. Regulation of p53 by Mdm2 E3 ligase function is dispensable in embryogenesis and development, but essential in response to DNA damage. Cancer cell. 2014;26(2):235–47.

    Article  CAS  PubMed  Google Scholar 

  20. Hsieh P, Yamane K. DNA mismatch repair: molecular mechanism, cancer, and ageing. Mech Ageing Dev. 2008;129(7–8):391–407.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Bocquet N, Bizard AH, Abdulrahman W, Larsen NB, Faty M, Cavadini S, et al. Structural and mechanistic insight into Holliday-junction dissolution by topoisomerase IIIalpha and RMI1. Nat Struct Mol Biol. 2014;21(3):261–8.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. van Raam BJ, Ehrnhoefer DE, Hayden MR, Salvesen GS. Intrinsic cleavage of receptor-interacting protein kinase-1 by caspase-6. Cell Death Differ. 2013;20(1):86–96.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Hu JC, Chun YH, Al Hazzazzi T, Simmer JP. Enamel formation and amelogenesis imperfecta. Cells Tissues Organs. 2007;186(1):78–85.

    Article  PubMed  Google Scholar 

  24. Guo S, Lim D, Dong Z, Saunders TL, Ma PX, Marcelo CL, et al. DSPP: A regulatory protein for dental pulp stem cell identity and fate. Stem Cells Dev. 2014;23(23):2883–94.

    Article  CAS  PubMed  Google Scholar 

  25. Lin SX, Zhang Q, Zhang H, Yan K, Ward L, Lu YB, et al. Nucleus-targeted Dmp1 transgene fails to rescue dental defects in Dmp1 null mice. Int J Oral Sci. 2014;6(3):133–41.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Uchihashi K, Aoki S, Matsunobu A, Toda S. Osteoblast migration into type I collagen gel and differentiation to osteocyte-like cells within a self-produced mineralized matrix: a novel system for analyzing differentiation from osteoblast to osteocyte. Bone. 2013;56(1):227.

    Article  Google Scholar 

  27. Liu J, Nam HK, Campbell C, Gasque KC, Millan JL, Hatch NE. Tissue-nonspecific alkaline phosphatase deficiency causes abnormal craniofacial bone development in the Alpl mouse model of infantile hypophosphatasia. Bone. 2014;67C:81–94.

    Article  Google Scholar 

  28. Timpl R, Rohde H, Robey PG, Rennard SI, Foidart JM, Martin GR. Laminin--a glycoprotein from basement membranes. J Biol Chem. 1979;254(19):9933–7.

    CAS  PubMed  Google Scholar 

  29. Fukumoto S, Miner JH, Ida H, Fukumoto E, Yuasa K, Miyazaki H, et al. Laminin alpha 5 is required for dental epithelium growth and polarity and the development of tooth bud and shape. J Biol Chem. 2006;281(8):5008–16.

    Article  CAS  PubMed  Google Scholar 

  30. Gumbiner BM. Regulation of cadherin-mediated adhesion in morphogenesis. Nat Rev Mol Cell Biol. 2005;6(8):622–34.

    Article  CAS  PubMed  Google Scholar 

  31. Stoeckli ET, Kilinc D, Kunz B, Kunz S, Lee GU, Martines E, et al. Analysis of cell-cell contact mediated by Ig superfamily cell adhesion molecules. Curr Protoc Cell Biol. 2013;61:9.5.1–9.5.85.

    Article  Google Scholar 

  32. Marie PJ, Hay E, Saidak Z. Integrin and cadherin signaling in bone: role and potential therapeutic targets. Trends Endocrinol Metab. 2014;25(11):567–75.

    Article  CAS  PubMed  Google Scholar 

  33. Felfly H, Klein OD. Sprouty genes regulate proliferation and survival of human embryonic stem cells. Sci Rep. 2013;3:2277.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Edwards AK, van den Heuvel MJ, Wessels JM, Lamarre J, Croy BA, Tayade C. Expression of angiogenic basic fibroblast growth factor, platelet derived growth factor, thrombospondin-1 and their receptors at the porcine maternal-fetal interface. Reprod Biol Endocrinol. 2011;9:5.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Otsuki S, Hanson SR, Miyaki S, Grogan SP, Kinoshita M, Asahara H, et al. Extracellular sulfatases support cartilage homeostasis by regulating BMP and FGF signaling pathways. Proc Natl Acad Sci U S A. 2010;107(22):10202–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Nan Su MJ, Chen L. Role of FGF/FGFR signaling in skeletal development and homeostasis: learning from mouse models. Bone Research. 2014;2:14003.

    Article  Google Scholar 

  37. Sun F, Wan M, Xu X, Gao B, Zhou Y, Sun J, et al. Crosstalk between miR-34a and Notch signaling promotes differentiation in Apical Papilla Stem Cells (SCAPs). J Dent Res. 2014;93(6):589–95.

    Article  CAS  PubMed  Google Scholar 

  38. Boucherat O, Nadeau V, Berube-Simard FA, Charron J, Jeannotte L. Crucial requirement of ERK/MAPK signaling in respiratory tract development. Development. 2014;141(16):3197–211.

    Article  CAS  PubMed  Google Scholar 

  39. Wang S, Choi M, Richardson AS, Reid BM, Seymen F, Yildirim M, et al. STIM1 and SLC24A4 are critical for enamel maturation. J Dent Res. 2014;93(7 suppl):94S–100S.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Duan X. Ion channels, channelopathies, and tooth formation. J Dent Res. 2014;93(2):117–25.

    Article  CAS  PubMed  Google Scholar 

  41. Yee NS, Kazi AA, Yee RK. Cellular and developmental biology of TRPM7 Channel-Kinase: implicated roles in cancer. Cells. 2014;3(3):751–77.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Molecular biology of the cell. 5th ed. New York: Garland Science; 2008.

    Google Scholar 

  43. Han JD, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, et al. Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature. 2004;430(6995):88–93.

    Article  CAS  PubMed  Google Scholar 

  44. Carlson MR, Zhang B, Fang Z, Mischel PS, Horvath S, Nelson SF. Gene connectivity, function, and sequence conservation: predictions from modular yeast co-expression networks. BMC Genomics. 2006;7:40.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.

    PubMed  Google Scholar 

  46. Appleford PJ, Woollard A. RUNX genes find a niche in stem cell biology. J Cell Biochem. 2009;108(1):14–21.

    Article  CAS  PubMed  Google Scholar 

  47. Guan X, Bartlett JD. MMP20 modulates cadherin expression in ameloblasts as enamel develops. J Dent Res. 2013;92(12):1123–8.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Yin Y, Wang F, Ornitz DM. Mesothelial- and epithelial-derived FGF9 have distinct functions in the regulation of lung development. Development. 2011;138(15):3169–77.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Laursen WJ, Bagriantsev SN, Gracheva EO. TRPA1 channels: chemical and temperature sensitivity. Curr Top Membr. 2014;74:89–112.

    Article  PubMed  Google Scholar 

  50. Willardsen M, Hutcheson DA, Moore KB, Vetter ML. The ETS transcription factor Etv1 mediates FGF signaling to initiate proneural gene expression during Xenopus laevis retinal development. Mech Dev. 2014;131:57–67.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Bond J, Roberts E, Mochida GH, Hampshire DJ, Scott S, Askham JM, et al. ASPM is a major determinant of cerebral cortical size. Nat Genet. 2002;32(2):316–20.

    Article  CAS  PubMed  Google Scholar 

  52. Sun SC, Lee SE, Xu YN, Kim NH. Perturbation of Spc25 expression affects meiotic spindle organization, chromosome alignment and spindle assembly checkpoint in mouse oocytes. Cell cycle. 2010;9(22):4552–9.

    Article  CAS  PubMed  Google Scholar 

  53. Nie J, Wang H, He F, Huang H. Nusap1 is essential for neural crest cell migration in zebrafish. Protein Cell. 2010;1(3):259–66.

    Article  CAS  PubMed  Google Scholar 

  54. Oh S, Shin S, Janknecht R. ETV1, 4 and 5: an oncogenic subfamily of ETS transcription factors. Biochim Biophys Acta. 2012;1826(1):1–12.

    CAS  PubMed Central  PubMed  Google Scholar 

  55. Fukumoto S, Yamada Y. Review: extracellular matrix regulates tooth morphogenesis. Connect Tissue Res. 2005;46(4–5):220–6.

    Article  CAS  PubMed  Google Scholar 

  56. Denecker G, Vandamme N, Akay O, Koludrovic D, Taminau J, Lemeire K, et al. Identification of a ZEB2-MITF-ZEB1 transcriptional network that controls melanogenesis and melanoma progression. Cell Death Differ. 2014;21(8):1250–61.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Gao H, Wang Y, Liu X, Yan T, Wu L, Alm E, et al. Global transcriptome analysis of the heat shock response of Shewanella oneidensis. J Bacteriol. 2004;186(22):7796–803.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Pokutta S, Weis WI. Structure and mechanism of cadherins and catenins in cell-cell contacts. Annu Rev Cell Dev Biol. 2007;23:237–61.

    Article  CAS  PubMed  Google Scholar 

  59. Maupin KA, Scha CJ, Williams BO. A comprehensive overview of skeletal phenotypes associated with alterations in Wnt/β-catenin signaling in humans and mice. Bone Research. 2013;1(1):27–71.

    Article  CAS  Google Scholar 

  60. Vainio S, Karavanova I, Jowett A, Thesleff I. Identification of BMP-4 as a signal mediating secondary induction between epithelial and mesenchymal tissues during early tooth development. Cell. 1993;75(1):45–58.

    Article  CAS  PubMed  Google Scholar 

  61. Espina V, Milia J, Wu G, Cowherd S, Liotta LA. Laser capture microdissection. Methods Mol Biol. 2006;319:213–29.

    Article  CAS  PubMed  Google Scholar 

  62. Espina V, Wulfkuhle JD, Calvert VS, VanMeter A, Zhou W, Coukos G, et al. Laser-capture microdissection. Nat Protoc. 2006;1(2):586–603.

    Article  CAS  PubMed  Google Scholar 

  63. Wan M, Gao B, Sun F, Tang Y, Ye L, Fan Y, et al. microRNA miR-34a regulates cytodifferentiation and targets multi-signaling pathways in human dental papilla cells. PLoS One. 2012;7(11):e50090.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Sun JX, Horst OV, Bumgarner R, Lakely B, Somerman MJ, Zhang H. Laser capture microdissection enables cellular and molecular studies of tooth root development. Int J Oral Sci. 2012;4(1):7–13.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Wettenhall JM, Smyth GK. limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004;20(18):3705–6.

    Article  CAS  PubMed  Google Scholar 

  66. Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP. GSEA-P: a desktop application for Gene Set Enrichment Analysis. Bioinformatics. 2007;23(23):3251–3.

    Article  CAS  PubMed  Google Scholar 

  67. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Langfelder P, Zhang B, Horvath S. Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics. 2008;24(5):719–20.

    Article  CAS  PubMed  Google Scholar 

  69. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  70. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grant from the National Natural Science Foundation of China (grant number: 81371136 to XZ, 81470711 and 81200760 to LZ). YN was supported by China Scholarship Council.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Yan Zhang or Liwei Zheng.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CL and YN performed the majority of the experiments, carried data analysis and drafted the manuscript. XZ helped to design the study and draft the manuscript. XX and YY provided suggestions for the project and critically reviewed the manuscript. YZ and LZ supervised the project and wrote most of the manuscript. All authors read and approved the final manuscript.

Chengcheng Liu and Yulong Niu contributed equally to this work.

Additional files

Additional file 1:

Differentially expressed genes between PABs and SABs. Table S1. Genes significantly up-regulated in PABs. Table S2. Genes significantly up-regulated in SABs. (XLS 194 kb)

Additional file 2: Table S3.

Expression of selected genes in PABs compared to SABs. (XLS 40 kb)

Additional file 3: Table S4 and Table S5.

Selected connections, including gene symbols of nodes and adjacency, in the brown and blue modules. Table S6 and Table S7: Selected connections with membership measure and gene significance in the brown and blue modules. (XLS 933 kb)

Additional file 4:

Scatterplots between module membership measure (x-axis) and gene significance of GO terms in the brown module (A-E) and blue module (F-J). The grey dashed line in the plot is the threshold for choosing significantly expressed genes, and the threshold value is –log10(0.05). (TIFF 1712 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, C., Niu, Y., Zhou, X. et al. Cell cycle control, DNA damage repair, and apoptosis-related pathways control pre-ameloblasts differentiation during tooth development. BMC Genomics 16, 592 (2015). https://doi.org/10.1186/s12864-015-1783-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-015-1783-y

Keywords