Ameloblastoma is the most frequent odontogenic tumor. Various evidence has highlighted the role of somatic mutations, including recurrent mutation BRAF V600E, in the tumorigenesis of Ameloblastoma, but the intact genetic pathology remains unknown. We sequenced the whole exome of both tumor tissue and healthy bone tissue from four mandibular ameloblastoma patients. The identified somatic mutations were integrated into Weighted Gene Co-expression Network Analysis on publicly available expression data of odontoblast, ameloblast, and Ameloblastoma. We identified a total of 70 rare and severe somatic mutations. We found BRAF V600E on all four patients, supporting previous discovery. HSAP4 was also hit by two missense mutations on two different patients. By applying Weighted Gene Co-expression Network Analysis on expression data of odontoblast, ameloblast, and Ameloblastoma, we found a proliferation-associated gene module that was significantly disrupted in tumor tissues. Each patient carried at least two rare, severe somatic mutations affecting genes within this module, including HSPA4, GNAS, CLTC, NES, and KMT2D. All these mutations had a ratio of variant-support reads lower than BRAF V600E, indicating that they occurred later than BRAF V600E. We suggest that a severe somatic mutation on the gene network of cell proliferation other than BRAF V600E, namely second hit, may contribute to the tumorigenesis of Ameloblastoma.