Cell-free DNA from nail clippings as source of normal control for genomic studies in hematologic malignancies.

Comprehensive genomic sequencing is becoming a critical component in the assessment of hematologic malignancies, with broad implications for patient management. In this context, unequivocally discriminating somatic from germline events is challenging but greatly facilitated by matched analysis of tumor:normal pairs. In contrast to solid tumors, conventional sources of normal control (peripheral blood, buccal swabs, saliva) could be highly involved by the neoplastic process, rendering them unsuitable. In this work we describe our real-world experience using cell free DNA (cfDNA) isolated from nail clippings as an alternate source of normal control, through the dedicated review of 2,610 tumor:nail pairs comprehensively sequenced by MSK-IMPACT-heme. Overall, we find nail cfDNA is a robust source of germline control for paired genomic studies. In a subset of patients, nail DNA may have tumor DNA contamination, reflecting unique attributes of the hematologic disease and transplant history. Contamination is generally low level, but significantly more common among patients with myeloid neoplasms (20.5%; 304/1482) compared to lymphoid diseases (5.4%; 61/1128) and particularly enriched in myeloproliferative neoplasms with marked myelofibrosis. When identified in patients with lymphoid and plasma-cell neoplasms, mutations commonly reflected a myeloid profile and correlated with a concurrent/evolving clonal myeloid neoplasm. For nails collected after allogeneic stem-cell transplantation, donor DNA was identified in 22% (11/50). In this cohort, an association with recent history of graft-vs-host disease was identified. These findings should be considered as a potential limitation for the use of nail as normal control but could also provide important diagnostic information regarding the disease process.

A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Enhanced clinical assessment of hematologic malignancies through routine paired tumor and normal sequencing.

Genomic profiling of hematologic malignancies has augmented our understanding of variants that contribute to disease pathogenesis and supported development of prognostic models that inform disease management in the clinic. Tumor only sequencing assays are limited in their ability to identify definitive somatic variants, which can lead to ambiguity in clinical reporting and patient management. Here, we describe the MSK-IMPACT Heme cohort, a comprehensive data set of somatic alterations from paired tumor and normal DNA using a hybridization capture-based next generation sequencing platform. We highlight patterns of mutations, copy number alterations, and mutation signatures in a broad set of myeloid and lymphoid neoplasms. We also demonstrate the power of appropriate matching to make definitive somatic calls, including in patients who have undergone allogeneic stem cell transplant. We expect that this resource will further spur research into the pathobiology and clinical utility of clinical sequencing for patients with hematologic neoplasms.

Prospective Clinical Genomic Profiling of Ewing Sarcoma: ERF and FGFR1 Mutations as Recurrent Secondary Alterations of Potential Biologic and Therapeutic Relevance.

PURPOSE: Ewing sarcoma (ES) is a primitive sarcoma defined by EWSR1-ETS fusions as the primary driver alteration. To better define the landscape of cooperating secondary genetic alterations in ES, we analyzed clinical genomic profiling data of 113 patients with ES, a cohort including more adult patients (> 18 years) and more patients with advanced stage at presentation than previous genomic cohorts. METHODS: The data set consisted of patients with ES prospectively tested with the US Food and Drug Administration-cleared Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets large panel, hybrid capture-based next-generation sequencing assay. To assess the functional significance of ERF loss, we generated ES cell lines with increased expression of ERF and lines with knockdown of ERF. We assessed cell viability, clonogenic growth, and motility in these ES lines and performed transcriptomic and epigenetic analyses. Finally, we validated our findings in vivo using cell line xenografts. RESULTS: Novel subsets were defined by recurrent secondary alterations in ERF, which encodes an ETS domain transcriptional repressor, in 7% of patients (five truncating mutations, one deep deletion, and two missense mutations) and in FGFR1 in another 2.7% (one amplification and two known activating mutations). ERF alterations were nonoverlapping with STAG2 alterations. In vitro, increased expression of ERF decreased tumor cell growth, colony formation, and motility in two ES cell lines, whereas ERF loss induced cellular proliferation and clonogenic growth. Transcriptomic analysis of cell lines with ERF loss revealed an increased expression of genes and pathways associated with aggressive tumor biology, and epigenetic, chromatin-based studies revealed that ERF competes with EWSR1-FLI1 at ETS-binding sites. CONCLUSION: Our findings open avenues to new insights into ES pathobiology and to novel therapeutic approaches in a subset of patients with ES.

Clinical, genomic, and transcriptomic correlates of response to immune checkpoint blockade-based therapy in a cohort of patients with angiosarcoma treated at a single center.

BACKGROUND: Angiosarcoma is a histologically and molecularly heterogeneous vascular neoplasm with aggressive clinical behavior. Emerging data suggests that immune checkpoint blockade (ICB) is efficacious against some angiosarcomas, particularly cutaneous angiosarcoma of the head and neck (CHN). METHODS: Patients with histologically confirmed angiosarcoma treated with ICB-based therapy at a comprehensive cancer center were retrospectively identified. Clinical characteristics and the results of targeted exome sequencing, transcriptome sequencing, and immunohistochemistry analyses were examined for correlation with clinical benefit. Durable clinical benefit was defined as a progression-free survival (PFS) of ≥16 weeks. RESULTS: For the 35 patients included in the analyses, median PFS and median overall survival (OS) from the time of first ICB-based treatment were 11.9 (95% CI 7.4 to 31.9) and 42.5 (95% CI 19.6 to 114.2) weeks, respectively. Thirteen patients (37%) had PFS ≥16 weeks. Clinical factors associated with longer PFS and longer OS in multivariate analyses were ICB plus other therapy regimens, CHN disease, and white race. Three of 10 patients with CHN angiosarcoma evaluable for tumor mutational burden (TMB) had a TMB ≥10. Five of six patients with CHN angiosarcoma evaluable for mutational signature analysis had a dominant mutational signature associated with ultraviolet (UV) light. No individual gene or genomic pathway was significantly associated with PFS or OS; neither were TMB or UV signature status. Analyses of whole transcriptomes from nine patient tumor samples found upregulation of angiogenesis, inflammatory response, and KRAS signaling pathways, among others, in patients with PFS ≥16 weeks, as well as higher levels of cytotoxic T cells, dendritic cells, and natural killer cells. Patients with PFS <16 weeks had higher numbers of cancer-associated fibroblasts. Immunohistochemistry findings for 12 patients with baseline samples available suggest that neither PD-L1 expression nor presence of tumor-infiltrating lymphocytes at baseline appears necessary for a response to ICB-based therapy. CONCLUSIONS: ICB-based therapy benefits only a subset of angiosarcoma patients. Patients with CHN angiosarcoma are more likely to have PFS ≥16 weeks, a dominant UV mutational signature, and higher TMB than angiosarcomas arising from other primary sites. However, clinical benefit was seen in other angiosarcomas also and was not restricted to tumors with a high TMB, a dominant UV signature, PD-L1 expression, or presence of tumor infiltrating lymphocytes at baseline.

Defining Novel DNA Virus-Tumor Associations and Genomic Correlates Using Prospective Clinical Tumor/Normal Matched Sequencing Data.

This study is the largest analysis of DNA viruses in solid tumors with associated genomics. To achieve this, a novel method for discovery of DNA viruses from matched tumor/normal next-generation sequencing samples was developed and validated. This method performed comparably to reference methods for the detection of high-risk (HR) human papilloma virus (HPV) (area under the receiver operating characteristic curve = 0.953). After virus identification in 48,148 consecutives samples from 42,846 unique patients, novel virus tumor associations were established by segregating tumor types to determine whether each DNA virus was enriched in each of the tumor types compared with the remaining cohort. All firmly established solid tumor-virus associations (eg, HR HPV in cervical cancer) were confirmed, and the novel associations discovered included: human herpes virus 6 in neuroblastoma, human herpes virus 7 in esophagogastric cancer, and HPV42 in digital papillary adenocarcinoma. These associations were confirmed in an independent validation cohort. HR HPV- and Epstein-Barr virus-associated tumors showed newly discovered genomic associations, including a lower tumor mutation burden. The study demonstrated the ability to study the role of DNA viruses in human cancer from clinical genomics data and established the largest cohort that can be utilized as a validation set for future discovery efforts.

Same-Cell Co-Occurrence of RAS Hotspot and BRAF V600E Mutations in Treatment-Naive Colorectal Cancer.

PURPOSE: Mitogen-activated protein kinase pathway-activating mutations occur in the majority of colorectal cancer (CRC) cases and show mutual exclusivity. We identified 47 epidermal growth factor receptor/BRAF inhibitor-naive CRC patients with dual RAS hotspot/BRAF V600E mutations (CRC-DD) from a cohort of 4,561 CRC patients with clinical next-generation sequencing results. We aimed to define the molecular phenotypes of the CRC-DD and to test if the dual RAS hotspot/BRAF V600E mutations coexist within the same cell. MATERIALS AND METHODS: We developed a single-cell genotyping method with a mutation detection rate of 96.3% and a genotype prediction accuracy of 92.1%. Mutations in the CRC-DD cohort were analyzed for clonality, allelic imbalance, copy number, and overall survival. RESULTS: Application of single-cell genotyping to four CRC-DD revealed the co-occurrence of both mutations in the following percentages of cells per case: NRAS G13D/KRAS G12C, 95%; KRAS G12D/NRAS G12V, 48%; BRAF V600E/KRAS G12D, 44%; and KRAS G12D/NRAS G13V, 14%, respectively. Allelic imbalance favoring the oncogenic allele was less frequent in CRC-DD (24 of 76, 31.5%, somatic mutations) compared with a curated cohort of CRC with a single-driver mutation (CRC-SD; 119 of 232 mutations, 51.3%; P = .013). Microsatellite instability-high status was enriched in CRC-DD compared with CRC-SD (23% v 11.4%, P = .028). Of the seven CRC-DD cases with multiregional sequencing, five retained both driver mutations throughout all sequenced tumor sites. Both CRC-DD cases with discordant multiregional sequencing were microsatellite instability-high. CONCLUSION: Our findings indicate that dual-driver mutations occur in a rare subset of CRC, often within the same tumor cells and across multiple tumor sites. Their presence and a lower rate of allelic imbalance may be related to dose-dependent signaling within the mitogen-activated protein kinase pathway.

Next-generation sequencing analysis suggests varied multistep mutational pathogenesis for endocrine mucin-producing sweat gland carcinoma with comments on INSM1 and MUC2 suggesting a conjunctival origin.

Endocrine mucin-producing sweat gland carcinoma is a low-grade eyelid tumor. Small biopsies and insensitive immunohistochemistry predispose to misdiagnosis. We aimed to identify clarifying immunohistochemical markers, molecular markers, or both. Clinicopathologic data (22 cases) were reviewed. Immunohistochemistry (insulinoma-associated protein 1, BCL-2, mucin 2 [MUC2], mucin 4, androgen receptor, ß-catenin, and Merkel cell polyomavirus) and next-generation sequencing (Memorial Sloan Kettering integrated mutation profiling of actionable cancer targets, 468 genes) were performed (3 cases). Female patients (n = 15) and male patients (n = 7) (mean age 71.8 years; range 53-88 years) had eyelid or periorbital tumors (>90%) with mucin-containing solid or cystic neuroendocrine pathology. Immunohistochemistry (insulinoma-associated protein 1, BCL2, androgen receptor, retinoblastoma-associated protein 1, and ß-catenin) was diffusely positive (5/5), MUC2 partial, mucin 4 focal, and Merkel cell polyomavirus negative. Memorial Sloan Kettering integrated mutation profiling of actionable cancer targets identified 12 single-nucleotide variants and 1 in-frame deletion in 3 cases, each with DNA damage response or repair (BRD4, PPP4R2, and RTEL1) and tumor-suppressor pathway (BRD4, TP53, TSC1, and LATS2) mutations. Microsatellite instability, copy number alterations, and structural alterations were absent. Insulinoma-associated protein 1 and MUC2 are positive in endocrine mucin-producing sweat gland carcinoma. MUC2 positivity suggests conjunctival origin. Multistep pathogenesis involving DNA damage repair and tumor-suppressor pathways may be implicated.

Quantitative Off-Target Detection of Epstein-Barr Virus-Derived DNA in Routine Molecular Profiling of Hematopoietic Neoplasms by Panel-Based Hybrid-Capture Next-Generation Sequencing.

Epstein-Barr virus (EBV) is associated with hematologic and solid tumors. We utilized a hybridization capture-based next-generation sequencing (NGS) platform targeting 400 genes associated with hematological malignancies to detect and quantify nontargeted viral-derived EBV reads that aligned to the EBV reference contig (NC_007605). We evaluated 5234 samples from 3636 unique patients with hematological neoplasms and found that 100 samples (1.9%) in 93 unique patients had ≥6 EBV reads (range, 6 to 32,325; mean, 827.5; median, 54). Most (n = 73, 73%) represented known EBV-associated conditions, and the most common was post-transplant lymphoproliferative disorders (n = 21, 29%). Documented EBV viremia was found in 4 of 27 samples with a moderate quantity of EBV reads and conditions not known to be EBV associated, whereas suspected viremia or low-level activation was likely in the remaining 23 samples. A good correlation (Spearman r = 0.8; 95% CI, 0.74-0.85) was found between EBV reads by NGS and systematic semiquantitative EBV-encoded RNA in situ hybridization in 162 available samples, particularly at greater EBV involvement. An optimal threshold for significant morphologic EBV involvement was found to be ≥10 reads by the receiver operating characteristic analysis (area under the curve, 0.990; 95% CI, 0.9974%-1.000%). Thus, in addition to mutational analysis, hybrid-capture-based NGS panels can detect and quantitate off-target EBV-derived viral DNA, which correlates well with morphology.

Enhanced specificity of clinical high-sensitivity tumor mutation profiling in cell-free DNA via paired normal sequencing using MSK-ACCESS.

Circulating cell-free DNA from blood plasma of cancer patients can be used to non-invasively interrogate somatic tumor alterations. Here we develop MSK-ACCESS (Memorial Sloan Kettering - Analysis of Circulating cfDNA to Examine Somatic Status), an NGS assay for detection of very low frequency somatic alterations in 129 genes. Analytical validation demonstrated 92% sensitivity in de-novo mutation calling down to 0.5% allele frequency and 99% for a priori mutation profiling. To evaluate the performance of MSK-ACCESS, we report results from 681 prospective blood samples that underwent clinical analysis to guide patient management. Somatic alterations are detected in 73% of the samples, 56% of which have clinically actionable alterations. The utilization of matched normal sequencing allows retention of somatic alterations while removing over 10,000 germline and clonal hematopoiesis variants. Our experience illustrates the importance of analyzing matched normal samples when interpreting cfDNA results and highlights the importance of cfDNA as a genomic profiling source for cancer patients.

Comprehensive Molecular Profiling of Desmoplastic Small Round Cell Tumor.

Desmoplastic small round cell tumor (DSRCT) is characterized by the EWSR1-WT1 t(11;22) (p13:q12) translocation. Few additional putative drivers have been identified, and research has suffered from a lack of model systems. Next-generation sequencing (NGS) data from 68 matched tumor-normal samples, whole-genome sequencing data from 10 samples, transcriptomic and affymetrix array data, and a bank of DSRCT patient-derived xenograft (PDX) are presented. EWSR1-WT1 fusions were noted to be simple, balanced events. Recurrent mutations were uncommon, but were noted in TERT (3%), ARID1A (6%), HRAS (5%), and TP53 (3%), and recurrent loss of heterozygosity (LOH) at 11p, 11q, and 16q was identified in 18%, 22%, and 34% of samples, respectively. Comparison of tumor-normal matched versus unmatched analysis suggests overcalling of somatic mutations in prior publications of DSRCT NGS data. Alterations in fibroblast growth factor receptor 4 (FGFR4) were identified in 5 of 68 (7%) of tumor samples, whereas differential overexpression of FGFR4 was confirmed orthogonally using 2 platforms. PDX models harbored the pathognomic EWSR1-WT1 fusion and were highly representative of corresponding tumors. Our analyses confirm DSRCT as a genomically quiet cancer defined by the balanced translocation, t(11;22)(p13:q12), characterized by a paucity of secondary mutations but a significant number of copy number alterations. Against this genomically quiet background, recurrent activating alterations of FGFR4 stood out, and suggest that this receptor tyrosine kinase, also noted to be highly expressed in DSRCT, should be further investigated. Future studies of DSRCT biology and preclinical therapeutic strategies should benefit from the PDX models characterized in this study. IMPLICATIONS: These data describe the general quiescence of the desmoplastic small round cell tumor (DSRCT) genome, present the first available bank of DSRCT model systems, and nominate FGFR4 as a key receptor tyrosine kinase in DSRCT, based on high expression, recurrent amplification, and recurrent activating mutations.

Next-generation assessment of human epidermal growth factor receptor 2 gene (ERBB2) amplification status in invasive breast carcinoma: a focus on Group 4 by use of the 2018 American Society of Clinical Oncology/College of American Pathologists HER2 testing guideline.

AIMS: The American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) updated the testing guideline in 2018 to address issues arising from uncommon human epidermal growth factor receptor 2 (HER2) fluorescence in-situ hybridisation (FISH) results according to the 2013 guideline. Next-generation sequencing (NGS) may be used to better classify patients. The aim of this study was to assess the ERBB2 amplification status of invasive breast carcinoma with equivocal HER2 immunohistochemistry (IHC) results by using NGS, focusing on Group 4 (HER2/CEP17 ratio of <2.0; average HER2 signals/cell of ≥4.0 and <6.0). METHODS AND RESULTS: We retrospectively reviewed HER2 FISH and NGS data of HER2 IHC-equivocal breast carcinomas at our centre between January 2009 and September 2019, wherein all three assays were performed on the same tissue block, and compared HER2 FISH results, according to the 2018 ASCO/CAP guideline, and the ERBB2 amplification status determined with NGS. A total of 52 HER2 FISH and NGS results from 51 patients with HER2 IHC-equivocal breast carcinomas were reviewed. The cohort included eight cases classified as 2018 ASCO/CAP in-situ hybridisation Group 1, three classified as Group 2, three classified as Group 3, 14 classified as Group 4, and 24 classified as Group 5. Thirteen of 14 (92.9%) Group 4 (HER2-negative) cases were classified as ERBB2-non-amplified by the use of NGS; the discordant case was later classified as Group 1 with alternative sample FISH testing. NGS revealed no significant difference in somatic mutations or copy number alterations between Groups 4 and 5. CONCLUSIONS: Our NGS findings support the reclassification of HER2 FISH-equivocal cases as HER2-negative under the 2018 ASCO/CAP guideline.

Therapeutic Potential of NTRK3 Inhibition in Desmoplastic Small Round Cell Tumor.

PURPOSE: Desmoplastic small round cell tumor (DSRCT) is a highly lethal intra-abdominal sarcoma of adolescents and young adults. DSRCT harbors a t(11;22)(p13:q12) that generates the EWSR1-WT1 chimeric transcription factor, the key oncogenic driver of DSRCT. EWSR1-WT1 rewires global gene expression networks and activates aberrant expression of targets that together mediate oncogenesis. EWSR1-WT1 also activates a neural gene expression program. EXPERIMENTAL DESIGN: Among these neural markers, we found prominent expression of neurotrophic tyrosine kinase receptor 3 (NTRK3), a druggable receptor tyrosine kinase. We investigated the regulation of NTRK3 by EWSR1-WT1 and its potential as a therapeutic target in vitro and in vivo, the latter using novel patient-derived models of DSRCT. RESULTS: We found that EWSR1-WT1 binds upstream of NTRK3 and activates its transcription. NTRK3 mRNA is highly expressed in DSRCT compared with other major chimeric transcription factor-driven sarcomas and most DSRCTs are strongly immunoreactive for NTRK3 protein. Remarkably, expression of NTRK3 kinase domain mRNA in DSRCT is also higher than in cancers with NTRK3 fusions. Abrogation of NTRK3 expression by RNAi silencing reduces growth of DSRCT cells and pharmacologic targeting of NTRK3 with entrectinib is effective in both in vitro and in vivo models of DSRCT. CONCLUSIONS: Our results indicate that EWSR1-WT1 directly activates NTRK3 expression in DSRCT cells, which are dependent on its expression and activity for growth. Pharmacologic inhibition of NTRK3 by entrectinib significantly reduces growth of DSRCT cells both in vitro and in vivo, providing a rationale for clinical evaluation of NTRK3 as a therapeutic target in DSRCT.

Molecular profiling and analysis of genetic aberrations aimed at identifying potential therapeutic targets in fibrolamellar carcinoma of the liver.

BACKGROUND: Fibrolamellar carcinoma (FLC) is a rare primary liver cancer of young adults. A functional chimeric transcript resulting from the in-frame fusion of the DNAJ homolog, subfamily B, member 1 (DNAJB1), and the catalytic subunit of protein kinase A (PRKACA) genes on chromosome 19 is believed to be unique in FLC, with a possible role in pathogenesis, yet with no established therapeutic value. The objective of the current study was to understand the molecular landscape of FLC and to identify potential novel therapeutic targets. METHODS: Archival fresh, formalin-fixed, paraffin-embedded samples from patients with FLC who prospectively consented to an institutional review board-approved protocol were analyzed using Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT), a next-generation sequencing assay encompassing up to 468 key cancer genes. Custom targeted RNA-Seq was performed in selected patients. Demographics, treatment, and outcome data were collected prospectively. Survival outcomes were estimated and correlated with mutation and/or copy number alterations. RESULTS: A total of 33 tumor samples from 31 patients with FLC were analyzed. The median age of the patients at the time of diagnosis was 18 years and approximately 53% were women. The DNAJB1-PRKACA fusion transcript was detected in 100% of patients. In 10 of 31 patients in which MSK-IMPACT did not detect the fusion, its presence was confirmed by targeted RNA-Seq. TERT promoter mutation was the second most common, and was detected in 7 patients. The median follow up was 30 months (range, 6-153 months). The 3-year overall survival rate was 84% (95% CI, 61%-93%). CONCLUSIONS: The DNAJB1-PRKACA fusion transcript is nonspecific and nonsensitive to FLC. Its potential therapeutic value currently is under evaluation. Opportunities currently are under development for therapy that may be driven or related to the DNAJB1-PRKACA fusion transcript or any therapeutic target identified from next-generation sequencing in patients with FLC.

Identification of prognostic molecular biomarkers in 157 HPV-positive and HPV-negative squamous cell carcinomas of the oropharynx.

The incidence of oropharyngeal squamous cell carcinoma (OPSCC) has been increasing due to high-risk HPV infection. We explored the significance of genetic alterations in HPV-positive (HPV-P) and HPV-negative (HPV-N) OPSCC patients on long-term outcome. A total of 157 cases of primary resected OPSCC diagnosed from 1978 to 2005 were subjected to a targeted exome sequencing by MSK-IMPACT™ interrogating somatic mutations in 410 cancer-related genes. Mutational profiles were correlated to recurrence and survival outcomes. OPSCC included 47% HPV-positive (HPV-P) and 53% HPV-negative (HPV-N) tumors arising in the base of tongue (BOT, 43%), palatine tonsil (30%) and soft palate (SP, 27%). HPV negative status, SP location and smoking were associated with poorer outcome. Poorer overall survival was found in NOTCH1-mutated HPV-P (p = 0.039), and in SOX2-amplified HPV-N cases (p = 0.036). Chromosomal arm gains in 8p and 8q, and 16q loss were more common in HPV-P (p = 0.005, 0.04 and 0.01, respectively), while 9p, 18q and 21q losses were more frequent in HPV-N OPSCC (p = 0.006, 0.002 and 0.01, respectively). Novel, potentially functional JAK3, MYC and EP300 intragenic deletions were found in HPV-P, and FOXP1, CDKN2A, CCND1 and RUNX1 intragenic deletions and one FGFR3 inversion were detected in HPV-N tumors. HPV-N/TP53-wild-type OPSCC harbored recurrent mutations in NOTCH1/3/4 (39%), PIK3CA, FAT1 and TERT. In comparison to their oral and laryngeal counterparts, HPV-N OPSCC were genetically distinct. In OPSCC, HPV status, tumor subsite and smoking determine outcome. Risk-stratification can be further refined based on the mutational signature, namely, NOTCH1 and SOX2 mutation status.

Complete Transcriptome Profiling of Normal and Age-Related Macular Degeneration Eye Tissues Reveals Dysregulation of Anti-Sense Transcription.

Age-related macular degeneration (AMD) predominantly affects the retina and retinal pigment epithelium in the posterior eye. While there are numerous studies investigating the non-coding transcriptome of retina and RPE, few significant differences between AMD and normal tissues have been reported. Strand specific RNA sequencing of both peripheral retina (PR) and RPE-Choroid-Sclera (PRCS), in both AMD and matched normal controls were generated. The transcriptome analysis reveals a highly significant and consistent impact on anti-sense transcription as well as moderate changes in the regulation of non-coding (sense) RNA. Hundreds of genes that do not express anti-sense transcripts in normal PR and PRCS demonstrate significant anti-sense expression in AMD in all patient samples. Several pathways are highly enriched in the upregulated anti-sense transcripts-in particular the EIF2 signaling pathway. These results call for a deeper exploration into anti-sense and noncoding RNA regulation in AMD and their potential as therapeutic targets.

Association of primary open-angle glaucoma with mitochondrial variants and haplogroups common in African Americans.

PURPOSE: To estimate the population frequencies of all common mitochondrial variants and ancestral haplogroups among 1,999 subjects recruited for the Primary Open-Angle African American Glaucoma Genetics (POAAGG) Study, including 1,217 primary open-angle glaucoma (POAG) cases and 782 controls, and to identify ancestral subpopulations and mitochondrial mutations as potential risk factors for POAG susceptibility. METHODS: Subject classification by characteristic glaucomatous optic nerve findings and corresponding visual field defects, as defined by enrolling glaucoma specialists, stereo disc photography, phlebotomy, extraction of total DNA from peripheral blood or saliva, DNA quantification and normalization, PCR amplification of whole mitochondrial genomes, Ion Torrent deep semiconductor DNA sequencing on DNA pools ("Pool-seq"), Sanger sequencing of 3,479 individual mitochondrial DNAs, and bioinformatic analysis. RESULTS: The distribution of common African haplogroups within the POAAGG study population was broadly similar to prior surveys of African Americans. However, the POAG case population was found to be enriched in L1c2 haplogroups, which are defined in part by missense mutations m.6150G>A (Val83Ile, odds ratio [OR] 1.8, p=0.01), m.6253C>T (Met117Thr, rs200165736, OR 1.6, p=0.04), and m.6480G>A (Val193Ile, rs199476128, OR 4.6, p=0.04) in the cytochrome c oxidase subunit 1 (MT-CO1) gene and by a variant, m.2220A>G (OR 2.0, p=0.01), in MT-RNR2, which encodes the mitochondrial ribosomal 16s RNA gene. L2 haplogroups were predicted to be overrepresented in the POAG case population by Pool-seq, and the difference was confirmed to be significant with Sanger sequencing, that targeted the L2-associated variants m.2416T>C (rs28358580, OR 1.2, p=0.02) and m.2332C>T (OR 1.2, p=.02) in MT-RNR2. Another variant within MT-RNR2, m.3010G>A (rs3928306), previously implicated in sensitivity to the optic neuropathy-associated antibiotic linezolid, and arising on D4 and J1 lineages, associated with Leber hereditary optic neuropathy (LHON) severity, was confirmed to be common (>5%) but was not significantly enriched in the POAG cases. Two variants linked to the composition of the gut microbiome, m.15784T>C (rs527236194, haplogroup L2a1) and m.16390G>A (rs41378955, L2 haplogroups), were also enriched in the case DNA pools. CONCLUSIONS: These results implicate African mtDNA haplogroups L1c2, L1c2b, and L2 as risk factors for POAG. Approximately one in four African Americans have these mitochondrial ancestries, which may contribute to their elevated glaucoma risk. These haplogroups are defined in part by ancestral variants in the MT-RNR2 and/or MT-CO1 genes, several of which have prior disease associations, such as MT-CO1 missense variants that have been implicated in prostate cancer.

Saliva DNA quality and genotyping efficiency in a predominantly elderly population.

BACKGROUND: The question of whether DNA obtained from saliva is an acceptable alternative to DNA from blood is a topic of considerable interest for large genetics studies. We compared the yields, quality and performance of DNAs from saliva and blood from a mostly elderly study population. METHODS: Two thousand nine hundred ten DNAs from primarily elderly subjects (mean age ± standard deviation (SD): 65 ± 12 years), collected for the Primary Open-Angle African-American Glaucoma Genetics (POAAGG) study, were evaluated by fluorometry and/or spectroscopy. These included 566 DNAs from blood and 2344 from saliva. Subsets of these were evaluated by Sanger sequencing (n = 1555), and by microarray SNP genotyping (n = 94) on an Illumina OmniExpress bead chip platform. RESULTS: The mean age of subjects was 65, and 68 % were female in both the blood and saliva groups. The mean ± SD of DNA yield per ml of requested specimen was significantly higher for saliva (17.6 ± 17.8 µg/ml) than blood (13.2 ± 8.5 µg/ml), but the mean ± SD of total DNA yield obtained per saliva specimen (35 ± 36 µg from 2 ml maximum specimen volume) was approximately three-fold lower than from blood (106 ± 68 µg from 8 ml maximum specimen volume). The average genotyping call rates were >99 % for 43 of 44 saliva DNAs and >99 % for 50 of 50 for blood DNAs. For 22 of 23 paired blood and saliva samples from the same individuals, the average genotyping concordance rate was 99.996 %. High quality PCR Sanger sequencing was obtained from ≥ 98 % of blood (n = 297) and saliva (n = 1258) DNAs. DNA concentrations ≥10 ng/µl, corresponding to total yields ≥ 2 µg, were obtained for 94 % of the saliva specimens (n = 2344). CONCLUSIONS: In spite of inferior purity, the performance of saliva DNAs for microarray genotyping was excellent. Our results agree with other studies concluding that saliva collection is a viable alternative to blood. The potential to boost study enrollments and reduce subject discomfort is not necessarily offset by a reduction in genotyping efficiency. Saliva DNAs performed comparably to blood DNAs for PCR Sanger sequencing.

Transcriptome of the human retina, retinal pigmented epithelium and choroid.

The retina and its adjacent supporting tissues - retinal pigmented epithelium (RPE) and choroid - are critical structures in human eyes required for normal visual perception. Abnormal changes in these layers have been implicated in diseases such as age-related macular degeneration and glaucoma. With the advent of high-throughput methods, such as serial analysis of gene expression, cDNA microarray, and RNA sequencing, there is unprecedented opportunity to facilitate our understanding of the normal retina, RPE, and choroid. This information can be used to identify dysfunction in age-related macular degeneration and glaucoma. In this review, we describe the current status in our understanding of these transcriptomes through the use of high-throughput techniques.

Comprehensive analysis of gene expression in human retina and supporting tissues.

Understanding the influence of gene expression on the molecular mechanisms underpinning human phenotypic diversity is fundamental to being able to predict health outcomes and treat disease. We have carried out whole transcriptome expression analysis on a series of eight normal human postmortem eyes by RNA sequencing. Here we present data showing that ∼80% of the transcriptome is expressed in the posterior layers of the eye and that there is significant differential expression not only between the layers of the posterior part of the eye but also between locations of a tissue layer. These differences in expression also extend to alternative splicing and splicing factors. Differentially expressed genes are enriched for genes associated with psychiatric, immune and cardiovascular disorders. Enrichment categories for gene ontology included ion transport, synaptic transmission and visual and sensory perception. Lastly, allele-specific expression was found to be significant for CFH, C3 and CFB, which are known risk genes for age-related macular degeneration. These expression differences should be useful in determining the underlying biology of associations with common diseases of the human retina, retinal pigment epithelium and choroid and in guiding the analysis of the genomic regions involved in the control of normal gene expression.