Pan-squamous genomic profiling stratified by anatomic tumor site and viral association

Background: Squamous cell carcinomas (SCC) have diverse anatomic etiologies but may share common genomic biomarkers. We profiled 7,871 unique SCCs across nine anatomic sites to investigate commonality in genomic alterations (GA), tumor mutational burden (TMB), human papillomavirus (HPV) association, and mutational signatures. Methods: Tissue from over 8,100 unique SCC samples originating from nine anatomic sites (anogenital (anus, cervix, penis, vagina, vulva), esophagus, head and neck, lung, and skin) were sequenced by hybrid capture-based comprehensive genomic profiling to evaluate GA and TMB. About 3% of non-cutaneous SCC samples had UV signatures, indicative of potential primary site misdiagnoses, and were filtered from the analysis. Detection of HPV, including high-risk strains 16, 18, 31, 33, and 45, was implemented through de novo assembly of non-human sequencing reads and BLASTn comparison against all viral nucleotide sequences in the NCBI database. Results: The proportion of HPV+ patients by anatomic site varied, with the highest being anal (91%) and cervical (83%). The mutational landscape of each cohort was similar, regardless of anatomic origin, but clustered based on HPV status. The largest differences in GA frequency as stratified by HPV- vs. HPV+ were TP53 (87% vs. 12%), CDKN2A (45% vs. 6%), and PIK3CA (22% vs. 33%). The median TMB in cases originating from HPV-associated sites was similar, regardless of HPV status. Higher median TMB was observed in lung and skin cases, which exhibited significant enrichment of mutational signatures indicative of tobacco- and UV-induced DNA damage, respectively. Conclusions: HPV+ and HPV- SCC populations have distinct genomic profiles and, for the latter, anatomic site is correlated with TMB distribution, secondary to associated carcinogen exposure. As such, biomarkers such as TMB and UV signature can provide unexpected insight into site of origin misdiagnoses and may correlate with benefit from immune checkpoint inhibitors

CYLD mutation characterizes a subset of HPV-positive head and neck squamous cell carcinomas with distinctive genomics and frequent cylindroma-like histologic features

Mutations in the tumor suppressor CYLD, known to be causative of cylindromas, were recently described in a subset of high-risk (hr) HPV-positive head and neck squamous cell carcinomas (HNSCC). Pathologic and genetic characterization of these CYLD-mutant carcinomas, however, remains limited. Here, we investigated whether CYLD mutations characterize a histopathologically and genomically distinct subset of hrHPV-positive HNSCC. Comprehensive genomic profiling via hybrid capture-based DNA sequencing was performed on 703 consecutive head and neck carcinomas with hrHPV sequences, identifying 148 unique cases (21%) harboring CYLD mutations. Clinical data, pathology reports, and histopathology were reviewed. CYLD mutations included homozygous deletions (n = 61/148; 41%), truncations (n = 52; 35%), missense (n = 26; 18%) and splice-site (n = 9; 6%) mutations, and in-frame deletion (n = 1; 1%). Among hrHPV-positive HNSCC, the CYLD-mutant cohort showed substantially lower tumor mutational burden than CYLD-wildtype cases (n = 555) (median 2.6 vs. 4.4 mut/Mb, p \u3c 0.00001) and less frequent alterations in PIK3CA (11% vs. 34%, p \u3c 0.0001), KMT2D (1% vs. 16%, p \u3c 0.0001), and FBXW7 (3% vs. 11%, p = 0.0018). Male predominance (94% vs. 87%), median age (58 vs. 60 years), and detection of HPV16 (95% vs. 89%) were similar. On available histopathology, 70% of CYLD-mutant HNSCC (98/141 cases) contained hyalinized material, consistent with basement membrane inclusions, within crowded aggregates of tumor cells. Only 7% of CYLD-wildtype cases demonstrated this distinctive pattern (p \u3c 0.0001). Histopathologic patterns of CYLD-mutant HNSCC lacking basement membrane inclusions included nonkeratinizing (n = 22, 16%), predominantly nonkeratinizing (nonkeratinizing SCC with focal maturation; n = 10, 7%), and keratinizing (n = 11, 8%) patterns. The latter two groups showed significantly higher frequency of PTEN alterations compared with other CYLD-mutant cases (38% [8/21] vs. 7% [8/120], p = 0.0004). Within our cohort of hrHPV-positive HNSCCs, CYLD mutations were frequent (21%) and demonstrated distinctive clinical, histopathologic, and genomic features that may inform future study of prognosis and treatment

Differential Expression of HERV-K (HML-2) Proviruses in Cells and Virions of the Teratocarcinoma Cell Line Tera-1

Human endogenous retrovirus (HERV-K (HML-2)) proviruses are among the few endogenous retroviral elements in the human genome that retain coding sequence. HML-2 expression has been widely associated with human disease states, including different types of cancers as well as with HIV-1 infection. Understanding of the potential impact of this expression requires that it be annotated at the proviral level. Here, we utilized the high throughput capabilities of next-generation sequencing to profile HML-2 expression at the level of individual proviruses and secreted virions in the teratocarcinoma cell line Tera-1. We identified well-defined expression patterns, with transcripts emanating primarily from two proviruses located on chromosome 22, only one of which was efficiently packaged. Interestingly, there was a preference for transcripts of recently integrated proviruses, over those from other highly expressed but older elements, to be packaged into virions. We also assessed the promoter competence of the 5’ long terminal repeats (LTRs) of expressed proviruses via a luciferase assay following transfection of Tera-1 cells. Consistent with the RNASeq results, we found that the activity of most LTRs corresponded to their transcript levels

Promoter expression of HERV-K (HML-2) provirus-derived sequences is related to LTR sequence variation and polymorphic transcription factor binding sites

Abstract Background Increased transcription of the human endogenous retrovirus group HERV-K (HML-2) is often seen during disease. Although the mechanism of its tissue-specific activation is unclear, research shows that LTR CpG hypomethylation alone is not sufficient to induce its promoter activity and that the transcriptional milieu of a malignant cell contributes, at least partly, to differential HML-2 expression. Results We analyzed the relationship between LTR sequence variation and promoter expression patterns in human breast cancer cell lines, finding them to be positively correlated. In particular, two proviruses (3q12.3 and 11p15.4) displayed increased activity in almost all tumorigenic cell lines sampled. Using a transcription factor binding site prediction algorithm, we identified two unique binding sites in each 5′ LTR that appeared to be associated with inducing promoter activity during neoplasia. Genomic analysis of the homologous proviruses in several non-human primates indicated post-integration genetic drift in two transcription factor binding sites, away from the ancestral sequence and towards the active form. Based on the sequences of 2504 individuals from the 1000 Genomes Project, the active form of the 11p15.4 site was found to be polymorphic within the human population, with an allele frequency of 51%, whereas the activating mutation in the 3q12.3 provirus was fixed in humans but not present in the orthologous provirus in chimpanzees or gorillas. Conclusions These data suggest that stage-specific transcription factors at least partly contribute to LTR promoter activity during transformation and that, in some cases, transcription factor binding site polymorphisms may be responsible for the differential HML-2 expression often seen between individuals

Budget impact analysis of comprehensive genomic profiling in patients with advanced non-small-cell lung cancer

PURPOSE This study assessed the economic impact of increased use of comprehensive genomic profiling (CGP) versus conventional testing strategies among patients with advanced non-small-cell lung cancer (aNSCLC) from a US commercial health plan perspective. METHODS A decision analytic model was developed to estimate the incremental benefits and costs across testing methodologies (CGP v non-CGP), as well as across sample types (tissue-based and liquid-based), for patients with newly diagnosed aNSCLC. Model outcomes included total direct costs, testing costs, and per member per month budget impact. Secondary model outcomes included the number of patients needed to test with CGP to add 1 life-year, and the number of patients needed to test with CGP to treat one individual with a biomarkermatched therapy. RESULTS In a hypothetical 2,000,000-member health plan, 790 members were estimated to have incident aNSCLC; 609 underwent molecular diagnostic testing with 122 (20%) tested with CGP (109 tissue-based and 13 liquid) in the base-case. An increase in CGP from 20% to 30% (an additional 61 patients tested with CGP) was associated with 3.11 additional life-years gained and a $0.01 in US dollars per member per month budget impact. Approximately 19.6 patients would need to be tested with CGP versus non-CGP to add one life-year and 5.9 patients would need to be tested with CGP to treat at least one patient with a biomarker-matched therapy. CONCLUSION An increase in CGP from 20% to 30% among patients with aNSCLC undergoing molecular diagnostic testing was associated with modest budget impact, most of which was attributable to prolonged survival associated with increased use of more effective treatments. © 2021 American Society of Clinical Oncology. All rights reserved.12 month embargo; first published: 14 October 2021This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Impaired outcome of controlled ovarian hyperstimulation in women with thyroid autoimmune disease.

Background: Controlled ovarian hyperstimulation (COH) is a crucial step of assisted reproductive technology (ART). Thyroid dysfunction and autoimmune thyroid disease (ATD) may negatively affect the outcome of ART, but the underlying mechanisms are still poorly understood. Our aim was to evaluate the respective role of ATD and thyroid function, as assessed by serum thyrotropin (TSH), on the early outcome of COH. Methods: In total, 262 (202 ATD-negative and 60 ATD-positive) euthyroid subfertile women underwent ART. Before COH, serum follicle-stimulating hormone (FSH), luteinizing hormone, and estradiol (E2) were measured at cycle day 3, and progesterone at cycle day 21. At oocyte pickup and at embryo transfer, we evaluated the performance of recombinant FSH (r-FSH), as assessed by serum E2 concentration/total administered r-FSH units (E2/r-FSH) ratio and by oocyte quality. Results: At both oocyte pickup and embryo transfer, the performance of r-FSH was significantly poorer in ATD-positive than in ATD-negative women. In the ATD-positive group, women with a TSH <2.5 mIU/L displayed a higher serum E2 concentration at oocyte pickup, a higher E2/r-FSH ratio, and a greater number of mature metaphase II oocytes than women with a TSH >2.5 mIU/L. When ATD-positive women were divided into quartiles according to their serum TSH level, both the ovarian response to r-FSH and the number of mature metaphase II oocytes significantly increased from the lowest to the highest quartiles of serum TSH concentration. Conclusions: ATD has a negative effect on the early outcome of COH, but this negative influence may be avoided with adequate levothyroxine therapy aimed at keeping TSH <2.5 mU/L. Thyroid antibodies and serum TSH should be checked in any woman undergoing ART

A pan-sarcoma landscape of telomeric content shows that alterations in RAD51B and GID4 are associated with higher telomeric content

Abstract Tumor cells need to activate a telomere maintenance mechanism, enabling limitless replication. The bulk of evidence supports that sarcomas predominantly use alternative lengthening of telomeres (ALT) mechanism, commonly associated with alterations in ATRX and DAXX. In our dataset, only 12.3% of sarcomas harbored alterations in these genes. Thus, we checked for the presence of other genomic determinants of high telomeric content in sarcomas. Our dataset consisted of 13555 sarcoma samples, sequenced as a part of routine clinical care on the FoundationOne®Heme platform. We observed a median telomeric content of 622.3 telomeric reads per GC-matched million reads (TRPM) across all samples. In agreement with previous studies, telomeric content was significantly higher in ATRX altered and POT1 altered sarcomas. We further observed that sarcomas with alterations in RAD51B or GID4 were enriched in samples with high telomeric content, specifically within uterus leiomyosarcoma for RAD51B and soft tissue sarcoma (not otherwise specified, nos) for GID4, Furthermore, RAD51B and POT1 alterations were mutually exclusive with ATRX and DAXX alterations, suggestive of functional redundancy. Our results propose a role played by RAD51B and GID4 in telomere elongation in sarcomas and open research opportunities for agents aimed at targeting this critical pathway in tumorigenesis

Phenotypic and Genomic Determinants of Immunotherapy Response Associated with Squamousness

Advanced and metastatic squamous cell carcinomas (SCC) are common and difficult-to-treat malignancies. We assessed 75 immunotherapy-treated patients with SCC from a clinically annotated database of 2,651 patients, as well as 9,407 patients from a deidentified database for molecular features that might influence checkpoint blockade response. SCCs had higher tumor mutational burdens (TMB) than non-SCCs (P &lt; 0.0001). Cutaneous SCCs had the highest TMB (P &lt; 0.0001), with 41.3% demonstrating a very high TMB (≥50 mutations/Mb). In immunotherapy-treated patients with SCC, higher TMB (≥12 mutations/Mb) correlated with a trend to higher clinical benefit rate [stable disease ≥ 6 months or partial/complete remission; 60% vs. 29%; (high vs. low TMB); P = 0.06] and significantly longer median time-to-treatment failure (TTF; 9.9 vs. 4.4 months; P = 0.0058). Cutaneous SCCs had the highest clinical benefit [11/15 patients (73%) vs. 20/60 (33%) non-cutaneous (P = 0.008)], TTF (P = 0.0015), and overall survival (P = 0.06) with immunotherapy treatment. In conclusion, among a diverse set of SCCs, higher TMB and cutaneous disease associated with better immunotherapy outcome

A computational approach to distinguish somatic vs. germline origin of genomic alterations from deep sequencing of cancer specimens without a matched normal.

A key constraint in genomic testing in oncology is that matched normal specimens are not commonly obtained in clinical practice. Thus, while well-characterized genomic alterations do not require normal tissue for interpretation, a significant number of alterations will be unknown in whether they are germline or somatic, in the absence of a matched normal control. We introduce SGZ (somatic-germline-zygosity), a computational method for predicting somatic vs. germline origin and homozygous vs. heterozygous or sub-clonal state of variants identified from deep massively parallel sequencing (MPS) of cancer specimens. The method does not require a patient matched normal control, enabling broad application in clinical research. SGZ predicts the somatic vs. germline status of each alteration identified by modeling the alteration's allele frequency (AF), taking into account the tumor content, tumor ploidy, and the local copy number. Accuracy of the prediction depends on the depth of sequencing and copy number model fit, which are achieved in our clinical assay by sequencing to high depth (>500x) using MPS, covering 394 cancer-related genes and over 3,500 genome-wide single nucleotide polymorphisms (SNPs). Calls are made using a statistic based on read depth and local variability of SNP AF. To validate the method, we first evaluated performance on samples from 30 lung and colon cancer patients, where we sequenced tumors and matched normal tissue. We examined predictions for 17 somatic hotspot mutations and 20 common germline SNPs in 20,182 clinical cancer specimens. To assess the impact of stromal admixture, we examined three cell lines, which were titrated with their matched normal to six levels (10-75%). Overall, predictions were made in 85% of cases, with 95-99% of variants predicted correctly, a significantly superior performance compared to a basic approach based on AF alone. We then applied the SGZ method to the COSMIC database of known somatic variants in cancer and found >50 that are in fact more likely to be germline