Purpose To identify novel therapeutic drug targets for p53 mutant head and neck squamous cell carcinoma (HNSCC). Experimental Design RNAi kinome viability screens were performed on HNSCC cells including autologous pairs from primary tumor and recurrent/metastatic lesions, and in parallel on murine squamous cell carcinoma (MSCC) cells derived from tumors of inbred mice bearing germline mutations in Trp53, and p53 regulatory genes: Atm, Prkdc, and p19Arf. Cross-species analysis of cell lines stratified by p53 mutational status and metastatic phenotype was utilized to select 38 kinase targets. Both primary and secondary RNAi validation assays were performed on additional HNSCC cell lines to credential these kinase targets utilizing multiple phenotypic endpoints. Kinase targets were also examined via chemical inhibition utilizing a panel of kinase inhibitors. A preclinical study was conducted on the WEE1 kinase inhibitor, MK-1775. Results Our functional kinomics approach identified novel survival kinases in HNSCC involved in G2/M cell cycle checkpoint, SFK, PI3K and FAK pathways. RNAi mediated knockdown and chemical inhibition of the WEE1 kinase with a specific inhibitor, MK-1775, had a significant effect on both viability and apoptosis. Sensitivity to the MK-1775 kinase inhibitor is in part determined by p53 mutational status, and due to unscheduled mitotic entry. MK-1775 displays single-agent activity and potentiates the efficacy of cisplatin in a p53 mutant HNSCC xenograft model. Conclusions WEE1 kinase is a potential therapeutic drug target for HNSCC. This study supports the application of a functional kinomics strategy to identify novel therapeutic targets for cancer.

BackgroundLymphotropism in oral squamous cell carcinoma (OSCC) is one of the most important prognostic factors of 5-year survival. In an effort to identify genes that may be responsible for the initiation of OSCC lymphotropism, we examined DNA copy number gains and losses and corresponding gene expression changes from tumor cells in metastatic lymph nodes of patients with OSCC.ResultsWe performed integrative analysis of DNA copy number alterations (CNA) and corresponding mRNA expression from OSCC cells isolated from metastatic lymph nodes of 20 patients using Affymetrix 250 K Nsp I SNP and U133 Plus 2.0 arrays, respectively. Overall, genome CNA accounted for expression changes in 31% of the transcripts studied. Genome region 11q13.2-11q13.3 shows the highest correlation between DNA CNA and expression. With a false discovery rate < 1%, 530 transcripts (461 genes) demonstrated a correlation between CNA and expression. Among these, we found two subsets that were significantly associated with OSCC (n = 122) when compared to controls, and with survival (n = 27), as tested using an independent dataset with genome-wide expression profiles for 148 primary OSCC and 45 normal oral mucosa. We fit Cox models to calculate a principal component analysis-derived risk-score for these two gene sets ('122-' or '27-transcript PC'). The models combining the 122- or 27-transcript PC with stage outperformed the model using stage alone in terms of the Area Under the Curve (AUC = 0.82 or 0.86 vs. 0.72, with p = 0.044 or 0.011, respectively).ConclusionsGenes exhibiting CNA-correlated expression may have biological impact on carcinogenesis and cancer progression in OSCC. Determination of copy number-associated transcripts associated with clinical outcomes in tumor cells with an aggressive phenotype (i.e., cells metastasized to the lymph nodes) can help prioritize candidate transcripts from high-throughput data for further studies.

Head and neck squamous cell carcinoma (HNSCC) is characterized by significant genomic instability that could lead to clonal diversity. Intratumor clonal heterogeneity has been proposed as a major attribute underlying tumor evolution, progression, and resistance to chemotherapy and radiation. Understanding genetic heterogeneity could lead to treatments specific to resistant and metastatic tumor cells. To characterize the degree of intratumor genetic heterogeneity within a single tumor, we performed whole-genome sequencing on three separate regions of an human papillomavirus (HPV)-positive oropharyngeal squamous cell carcinoma and two separate regions from one corresponding cervical lymph node metastasis. This approach achieved coverage of approximately 97.9% of the genome across all samples. In total, 5701 somatic point mutations (SPMs) and 4347 small somatic insertions and deletions (indels)were detected in at least one sample. Ninety-two percent of SPMs and 77% of indels were validated in a second set of samples adjacent to the discovery set. All five tumor samples shared 41% of SPMs, 57% of the 1805 genes with SPMs, and 34 of 55 cancer genes. The distribution of SPMs allowed phylogenetic reconstruction of this tumor's evolutionary pathway and showed that the metastatic samples arose as a late event. The degree of intratumor heterogeneity showed that a single biopsy may not represent the entire mutational landscape of HNSCC tumors. This approach may be used to further characterize intratumor heterogeneity in more patients, and their sample-to-sample variations could reveal the evolutionary process of cancer cells, facilitate our understanding of tumorigenesis, and enable the development of novel targeted therapies.

Background Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance. Results In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5–100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels. Conclusion These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

The Apollo butterfly, Parnassius glacialis, is one of the most charming members of its genus and includes two subspecies locally distributed in montane areas of south-central China and Japan. In this study, we investigated the genetic structure and demographic history of P. glacialis by analyzing partial sequences of four mitochondrial genes and nuclear single nucleotide polymorphisms (SNPs) via genotyping-by-sequencing (GBS) of samples from nearly the entire known distributional range in China. The mitochondrial DNA (mtDNA) data demonstrated that a total of 39 haplotypes were present, and the species was estimated to have diverged about 0.95 million years ago during the middle Pleistocene transition into two main clades that likely formed during the Kunlun-Huanghe tectonic movement. The two clades then dispersed independently in distinct geographic areas alongside the mountainous routes in central and southern China, most likely driven by the Pleistocene glacial-interglacial cycles. Nuclear SNP analysis was generally congruent with mtDNA results at the individual level. A minor incongruence of genetic structures that was detected between mtDNA and nuclear SNP data from the Laojunshan and Tiantangzhai populations was likely due to secondary contact and male-biased dispersal. Our work demonstrates that complicated dispersal-vicariance evolutionary processes likely led to the current geographic distribution of P. glacialis in China, particularly the uplift of the Qinghai-Tibet Plateau and related climatic oscillations during the Quaternary period.

No abstract

Background Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing. Results All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5–20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden. Conclusion This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

In oral squamous cell carcinoma (OSCC), metastasis to lymph nodes is associated with a 50% reduction in 5-year survival. To identify a metastatic gene set based on DNA copy number abnormalities (CNAs) of differentially expressed genes, we compared DNA and RNA of OSCC cells laser-microdissected from non-metastatic primary tumors (n = 17) with those from lymph node metastases (n = 20), using Affymetrix 250K Nsp single-nucleotide polymorphism (SNP) arrays and U133 Plus 2.0 arrays, respectively. With a false discovery rate (FDR)<5%, 1988 transcripts were found to be differentially expressed between primary and metastatic OSCC. Of these, 114 were found to have a significant correlation between DNA copy number and gene expression (FDR<0.01). Among these 114 correlated transcripts, the corresponding genomic regions of each of 95 transcripts had CNAs differences between primary and metastatic OSCC (FDR<0.01). Using an independent dataset of 133 patients, multivariable analysis showed that the OSCC–specific and overall mortality hazards ratio (HR) for patients carrying the 95-transcript signature were 4.75 (95% CI: 2.03–11.11) and 3.45 (95% CI: 1.84–6.50), respectively. To determine the degree by which these genes impact cell survival, we compared the growth of five OSCC cell lines before and after knockdown of over-amplified transcripts via a high-throughput siRNA–mediated screen. The expression-knockdown of 18 of the 26 genes tested showed a growth suppression ≥30% in at least one cell line (P<0.01). In particular, cell lines derived from late-stage OSCC were more sensitive to the knockdown of G3BP1 than cell lines derived from early-stage OSCC, and the growth suppression was likely caused by increase in apoptosis. Further investigation is warranted to examine the biological role of these genes in OSCC progression and their therapeutic potentials.