The incidence of papillary thyroid carcinomas (PTCs) has increased rapidly during the past several decades. Until now, the mechanisms underlying the tumorigenesis of PTCs have remained largely unknown. Next-generation-sequencing (NGS) provides new ways to investigate the molecular pathogenesis of PTCs. To characterize the somatic alterations associated with PTCs, we performed whole-exome sequencing (WES) of PTCs from 23 Chinese patients. This study revealed somatic mutations in genes with relevant functions for tumorigenesis, such as BRAF, BCR, CREB3L2, DNMT1, IRS2, MSH6, and TP53. We also identified novel somatic gene alterations which may be potentially involved in PTC progression. Gene set enrichment analysis revealed that the cellular response to hormone stimulus, epigenetic modifications, such as protein/histone methylation and protein alkylation, as well as MAPK, PI3K-AKT, and FoxO/mTOR signaling pathways, were significantly altered in the PTCs studied here. Moreover, Protein-Protein Interaction (PPI) network analysis of our mutated gene selection highlighted EP300, KRAS, PTEN, and TP53 as major core genes. The correlation between gene mutations and clinicopathologic features of the PTCs defined by conventional ultrasonography (US) and contrast-enhanced ultrasonography (CEUS) were assessed. These analyses established significant associations between subgroups of mutations and respectively taller-than-wide, calcified, and peak time iso- or hypo-enhanced and metastatic PTCs. In conclusion, our study supplements the genomic landscape of PTCs and identifies new actionable target candidates and clinicopathology-associated mutations. Extension of this study to larger cohorts will help define comprehensive genomic aberrations in PTCs and validate target candidates. These new targets may open methods of individualized treatments adapted to the clinicopathologic specifics of the patients.