Relationship between BCL2 mutations and follicular lymphoma outcome in the chemoimmunotherapy era.

How to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression is unclear. Building on a prior study demonstrating the early transformation of FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS, and MYC, in 199 newly diagnosed grade 1 and 2 FLs. BCL2 mutations with VAF ≥20% occurred in 52% of cases. Among 97 FL patients who did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF ≥20% were associated with increased transformation risk (HR 3.01, 95% CI 1.04-8.78, p = 0.043) and a trend toward shorter event-free survival (EFS, median 20 months with mutations versus 54 months without, p = 0.052). Other sequenced genes were less frequently mutated and did not increase the prognostic value of the panel. Across the entire population, nonsynonymous BCL2 mutations at VAF ≥20% were associated with decreased EFS (HR 1.55, 95% CI 1.02-2.35, p = 0.043 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.82, 95% CI 1.05-3.17, p = 0.034). Thus, high VAF nonsynonymous BCL2 mutations remain prognostic even in the chemoimmunotherapy era.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair.

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics.

PURPOSE: Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies. METHODS: Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations. RESULTS: Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases. CONCLUSION: NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

NUP155 insufficiency recalibrates a pluripotent transcriptome with network remodeling of a cardiogenic signaling module.

BACKGROUND: Atrial fibrillation is a cardiac disease driven by numerous idiopathic etiologies. NUP155 is a nuclear pore complex protein that has been identified as a clinical driver of atrial fibrillation, yet the precise mechanism is unknown. The present study employs a systems biology algorithm to identify effects of NUP155 disruption on cardiogenicity in a model of stem cell-derived differentiation. METHODS: Embryonic stem (ES) cell lines (n = 5) with truncated NUP155 were cultured in parallel with wild type (WT) ES cells (n = 5), and then harvested for RNAseq. Samples were run on an Illumina HiSeq 2000. Reads were analyzed using Strand NGS, Cytoscape, DAVID and Ingenuity Pathways Analysis to deconvolute the NUP155-disrupted transcriptome. Network topological analysis identified key features that controlled framework architecture and functional enrichment. RESULTS: In NUP155 truncated ES cells, significant expression changes were detected in 326 genes compared to WT. These genes segregated into clusters that enriched for specific gene ontologies. Deconvolution of the collective framework into discrete sub-networks identified a module with the highest score that enriched for Cardiovascular System Development, and revealed NTRK1/TRKA and SRSF2/SC35 as critical hubs within this cardiogenic module. CONCLUSIONS: The strategy of pluripotent transcriptome deconvolution used in the current study identified a novel association of NUP155 with potential drivers of arrhythmogenic AF. Here, NUP155 regulates cardioplasticity of a sub-network embedded within a larger framework of genome integrity, and exemplifies how transcriptome cardiogenicity in an embryonic stem cell genome is recalibrated by nucleoporin dysfunction.

Targeted sequencing of 36 known or putative colorectal cancer susceptibility genes.

BACKGROUND: Mutations in several genes predispose to colorectal cancer. Genetic testing for hereditary colorectal cancer syndromes was previously limited to single gene tests; thus, only a very limited number of genes were tested, and rarely those infrequently mutated in colorectal cancer. Next-generation sequencing technologies have made it possible to sequencing panels of genes known and suspected to influence colorectal cancer susceptibility. METHODS: Targeted sequencing of 36 known or putative CRC susceptibility genes was conducted for 1231 CRC cases from five subsets: (1) Familial Colorectal Cancer Type X (n = 153); (2) CRC unselected by tumor immunohistochemical or microsatellite stability testing (n = 548); (3) young onset (age <50 years) (n = 333); (4) proficient mismatch repair (MMR) in cases diagnosed at ≥50 years (n = 68); and (5) deficient MMR CRCs with no germline mutations in MLH1, MSH2, MSH6, or PMS2 (n = 129). Ninety-three unaffected controls were also sequenced. RESULTS: Overall, 29 nonsense, 43 frame-shift, 13 splice site, six initiator codon variants, one stop codon, 12 exonic deletions, 658 missense, and 17 indels were identified. Missense variants were reviewed by genetic counselors to determine pathogenicity; 13 were pathogenic, 61 were not pathogenic, and 584 were variants of uncertain significance. Overall, we identified 92 cases with pathogenic mutations in APC,MLH1,MSH2,MSH6, or multiple pathogenic MUTYH mutations (7.5%). Four cases with intact MMR protein expression by immunohistochemistry carried pathogenic MMR mutations. CONCLUSIONS: Results across case subsets may help prioritize genes for inclusion in clinical gene panel tests and underscore the issue of variants of uncertain significance both in well-characterized genes and those for which limited experience has accumulated.

Androgen Receptor Variant AR-V9 Is Coexpressed with AR-V7 in Prostate Cancer Metastases and Predicts Abiraterone Resistance.

Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies.Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate.Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3' terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31-12.2; P = 0.02).Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704-15. ©2017 AACR.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Target-enrichment sequencing and copy number evaluation in inherited polyneuropathy.

OBJECTIVE: To assess the efficiency of target-enrichment next-generation sequencing (NGS) with copy number assessment in inherited neuropathy diagnosis. METHODS: A 197 polyneuropathy gene panel was designed to assess for mutations in 93 patients with inherited or idiopathic neuropathy without known genetic cause. We applied our novel copy number variation algorithm on NGS data, and validated the identified copy number mutations using CytoScan (Affymetrix). Cost and efficacy of this targeted NGS approach was compared to earlier evaluations. RESULTS: Average coverage depth was ∼760× (median = 600, 99.4% > 100×). Among 93 patients, 18 mutations were identified in 17 cases (18%), including 3 copy number mutations: 2 PMP22 duplications and 1 MPZ duplication. The 2 patients with PMP22 duplication presented with bulbar and respiratory involvement and had absent extremity nerve conductions, leading to axonal diagnosis. Average onset age of these 17 patients was 25 years (2-61 years), vs 45 years for those without genetic discovery. Among those with onset age less than 40 years, the diagnostic yield of targeted NGS approach is high (27%) and cost savings is significant (∼20%). However, the cost savings for patients with late onset age and without family history is not demonstrated. CONCLUSIONS: Incorporating copy number analysis in target-enrichment NGS approach improved the efficiency of mutation discovery for chronic, inherited, progressive length-dependent polyneuropathy diagnosis. The new technology is facilitating a simplified genetic diagnostic algorithm utilizing targeted NGS, clinical phenotypes, age at onset, and family history to improve diagnosis efficiency. Our findings prompt a need for updating the current practice parameters and payer guidelines.

Prevalence of Pathogenic Mutations in Cancer Predisposition Genes among Pancreatic Cancer Patients.

The prevalence of germline pathogenic mutations in a comprehensive panel of cancer predisposition genes is not well-defined for patients with pancreatic ductal adenocarcinoma (PDAC). To estimate the frequency of mutations in a panel of 22 cancer predisposition genes, 96 patients unselected for a family history of cancer who were recruited to the Mayo Clinic Pancreatic Cancer patient registry over a 12-month period were screened by next-generation sequencing. Fourteen pathogenic mutations in 13 patients (13.5%) were identified in eight genes: four in ATM, two in BRCA2, CHEK2, and MSH6, and one in BARD1, BRCA1, FANCM, and NBN. These included nine mutations (9.4%) in established pancreatic cancer genes. Three mutations were found in patients with a first-degree relative with PDAC, and 10 mutations were found in patients with first- or second-degree relatives with breast, pancreas, colorectal, ovarian, or endometrial cancers. These results suggest that a substantial proportion of patients with PDAC carry germline mutations in predisposition genes associated with other cancers and that a better understanding of pancreatic cancer risk will depend on evaluation of families with broad constellations of tumors. These findings highlight the need for recommendations governing germline gene-panel testing of patients with pancreatic cancer.

Draft Genome Sequences of Nine Pseudomonas aeruginosa Strains, Including Eight Clinical Isolates.

We report on nine draft genomes of Pseudomonas aeruginosa isolates, assembled using a hybrid paired-end and Nextera mate-pair library approach. Eight are of clinical origin, and one is the ATCC 27853 strain. We also report their multilocus sequence types.

Recessive MYH6 Mutations in Hypoplastic Left Heart With Reduced Ejection Fraction.

BACKGROUND: The molecular underpinnings of hypoplastic left heart are poorly understood. Staged surgical palliation has dramatically improved survival, yet eventual failure of the systemic right ventricle necessitates cardiac transplantation in a subset of patients. We sought to identify genetic determinants of hypoplastic left heart with latent right ventricular dysfunction in individuals with a Fontan circulation. METHODS AND RESULTS: Evaluation of cardiac structure and function by echocardiography in patients with hypoplastic left heart and their first-degree relatives identified 5 individuals with right ventricular ejection fraction ≤40% after Fontan operation. Whole genome sequencing was performed on DNA from 21 family members, filtering for genetic variants with allele frequency <1% predicted to alter protein structure or expression. Secondary family-based filtering for de novo and recessive variants revealed rare inherited missense mutations on both paternal and maternal alleles of MYH6, encoding myosin heavy chain 6, in 2 patients who developed right ventricular dysfunction 3 to 11 years postoperatively. Parents and siblings who were heterozygous carriers had normal echocardiograms. Protein modeling of the 4 highly conserved amino acid substitutions, residing in both head and tail domains, predicted perturbation of protein structure and function. CONCLUSIONS: In contrast to dominant MYH6 mutations with variable penetrance identified in other congenital heart defects and dilated cardiomyopathy, this study reveals compound heterozygosity for recessive MYH6 mutations in patients with hypoplastic left heart and reduced systemic right ventricular ejection fraction. These findings implicate a shared molecular basis for the developmental arrest and latent myopathy of left and right ventricles, respectively.

Homozygous/Compound Heterozygous Triadin Mutations Associated With Autosomal-Recessive Long-QT Syndrome and Pediatric Sudden Cardiac Arrest: Elucidation of the Triadin Knockout Syndrome.

BACKGROUND: Long-QT syndrome (LQTS) may result in syncope, seizures, or sudden cardiac arrest. Although 16 LQTS-susceptibility genes have been discovered, 20% to 25% of LQTS remains genetically elusive. METHODS AND RESULTS: We performed whole-exome sequencing child-parent trio analysis followed by recessive and sporadic inheritance modeling and disease-network candidate analysis gene ranking to identify a novel underlying genetic mechanism for LQTS. Subsequent mutational analysis of the candidate gene was performed with polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing on a cohort of 33 additional unrelated patients with genetically elusive LQTS. After whole-exome sequencing and variant filtration, a homozygous p.D18fs*13 TRDN-encoded triadin frameshift mutation was discovered in a 10-year-old female patient with LQTS with a QTc of 500 milliseconds who experienced recurrent exertion-induced syncope/cardiac arrest beginning at 1 year of age. Subsequent mutational analysis of TRDN revealed either homozygous or compound heterozygous frameshift mutations in 4 of 33 unrelated cases of LQTS (12%). All 5 TRDN-null patients displayed extensive T-wave inversions in precordial leads V1 through V4, with either persistent or transient QT prolongation and severe disease expression of exercise-induced cardiac arrest in early childhood (≤3 years of age) and required aggressive therapy. The overall yield of TRDN mutations was significantly greater in patients ≤10 years of age (5 of 10, 50%) compared with older patients (0 of 24, 0%; P=0.0009). CONCLUSIONS: We identified TRDN as a novel underlying genetic basis for recessively inherited LQTS. All TRDN-null patients had strikingly similar phenotypes. Given the recurrent nature of potential lethal arrhythmias, patients fitting this phenotypic profile should undergo cardiac TRDN genetic testing.

Comprehensive assessment of genetic variants within TCF4 in Fuchs' endothelial corneal dystrophy.

PURPOSE: The single nucleotide variant (SNV), rs613872, in the transcription factor 4 (TCF4) gene was previously found to be strongly associated (P = 6 × 10(-26)) with Fuchs' endothelial corneal dystrophy (FECD). Subsequently, an intronic expansion of the repeating trinucleotides, TGC, was found to be even more predictive of disease. We performed comprehensive sequencing of the TCF4 gene region in order to identify the best marker for FECD within TCF4 and to identify other novel variants that may be associated with FECD. METHODS: Leukocyte DNA was isolated from 68 subjects with FECD and 16 unaffected individuals. A custom capture panel was used to isolate the region surrounding the two previously validated markers of FECD. Sequencing of the TCF4 coding region, introns and flanking sequence, spanning 465 kb was performed at >1000× average coverage using the Illumina HiSequation 2000. RESULTS: TGC expansion (>50 repeats) was present in 46 (68%) FECD-affected subjects and one (6%) normal subject. A total of 1866 variants, including 1540 SNVs, were identified. Only two previously reported SNVs resided in the TCF4 coding region, neither of which segregated with disease. No variant, including TGC expansion, correlated perfectly with disease status. Trinucleotide repeat expansion was a better predictor of disease than any other variant. CONCLUSIONS: Complete sequencing of the TCF4 genomic region revealed no single causative variant for FECD. The intronic trinucleotide repeat expansion within TCF4 continues to be more strongly associated with FECD than any other genetic variant.

Clinical biomarkers of pulmonary carcinoid tumors in never smokers via profiling miRNA and target mRNA.

BACKGROUND: miRNAs play key regulatory roles in cellular pathological processes. We aimed to identify clinically meaningful biomarkers in pulmonary carcinoid tumors (PCTs), a member of neuroendocrine neoplasms, via profiling miRNAs and mRNAs. RESULTS: From the total of 1145 miRNAs, we obtained 16 and 17 miRNAs that showed positive and negative fold changes (FCs, tumors vs. normal tissues) in the top 1% differentially expressed miRNAs, respectively. We uncovered the target genes that were predicted by at least two prediction tools and overlapped by at least one-half of the top miRNAs, which yielded 44 genes (FC<-2) and 56 genes (FC>2), respectively. Higher expressions of CREB5, PTPRB and COL4A3 predicted favorable disease free survival (Hazard ratio: 0.03, 0.19 and 0.36; P value: 0.03, 0.03 and 0.08). Additionally, 79 mutated genes have been found in nine PCTs where TP53 was the only repeated mutation. CONCLUSION: We identified that the expressions of three genes have clinical implications in PCTs. The biological functions of these biomarkers warrant further studies.

Conserved recurrent gene mutations correlate with pathway deregulation and clinical outcomes of lung adenocarcinoma in never-smokers.

BACKGROUND: Novel and targetable mutations are needed for improved understanding and treatment of lung cancer in never-smokers. METHODS: Twenty-seven lung adenocarcinomas from never-smokers were sequenced by both exome and mRNA-seq with respective normal tissues. Somatic mutations were detected and compared with pathway deregulation, tumor phenotypes and clinical outcomes. RESULTS: Although somatic mutations in DNA or mRNA ranged from hundreds to thousands in each tumor, the overlap mutations between the two were only a few to a couple of hundreds. The number of somatic mutations from either DNA or mRNA was not significantly associated with clinical variables; however, the number of overlap mutations was associated with cancer subtype. These overlap mutants were preferentially expressed in mRNA with consistently higher allele frequency in mRNA than in DNA. Ten genes (EGFR, TP53, KRAS, RPS6KB2, ATXN2, DHX9, PTPN13, SP1, SPTAN1 and MYOF) had recurrent mutations and these mutations were highly correlated with pathway deregulation and patient survival. CONCLUSIONS: The recurrent mutations present in both DNA and RNA are likely the driver for tumor biology, pathway deregulation and clinical outcomes. The information may be used for patient stratification and therapeutic target development.

Mate pair sequencing of oropharyngeal squamous cell carcinomas reveals that HPV integration occurs much less frequently than in cervical cancer.

BACKGROUND: Human papillomavirus (HPV) is now recognized to be very important in the pathogenesis of oropharyngeal squamous cell carcinoma (OPSCC). It is not clear yet whether the physical status of HPV in OPSCC is similar to what is found in cervical cancer. STUDY DESIGN: We performed genome-wide mate pair next generation sequencing from 20 OPSCCs patients, thirteen of which were positive for HPV16 to determine the HPV physical status and its relationship to HPV oncogene E6 and E7 expression. RESULTS: This high throughput approach detected HPV integration events and also determined the bridged HPV coverage in each sample. Two of the HPV16-positive OPSCCs had HPV integration and one of the HPV16-negative OPSCCs had an HPV26 integration. We mapped the site of integration in the HPV genome in all integration events and the integrations were located at E1, E5, E6 and L2 region respectively. One HPV positive OPSCC had two integration events but also had a very high bridged HPV coverage, while the other two just had HPV integrated into the human genome. CONCLUSION: Our results are thus different from what is routinely observed in cervical cancer where HPV is almost always integrated into the human genome with loss of episomal HPV sequences. Thus more investigation should be carried out to study how episomal HPV alone can contribute to the development of most OPSCCs.

Gene expression, single nucleotide variant and fusion transcript discovery in archival material from breast tumors.

Advantages of RNA-Seq over array based platforms are quantitative gene expression and discovery of expressed single nucleotide variants (eSNVs) and fusion transcripts from a single platform, but the sensitivity for each of these characteristics is unknown. We measured gene expression in a set of manually degraded RNAs, nine pairs of matched fresh-frozen, and FFPE RNA isolated from breast tumor with the hybridization based, NanoString nCounter (226 gene panel) and with whole transcriptome RNA-Seq using RiboZeroGold ScriptSeq V2 library preparation kits. We performed correlation analyses of gene expression between samples and across platforms. We then specifically assessed whole transcriptome expression of lincRNA and discovery of eSNVs and fusion transcripts in the FFPE RNA-Seq data. For gene expression in the manually degraded samples, we observed Pearson correlations of >0.94 and >0.80 with NanoString and ScriptSeq protocols, respectively. Gene expression data for matched fresh-frozen and FFPE samples yielded mean Pearson correlations of 0.874 and 0.783 for NanoString (226 genes) and ScriptSeq whole transcriptome protocols respectively, p<2x10(-16). Specifically for lincRNAs, we observed superb Pearson correlation (0.988) between matched fresh-frozen and FFPE pairs. FFPE samples across NanoString and RNA-Seq platforms gave a mean Pearson correlation of 0.838. In FFPE libraries, we detected 53.4% of high confidence SNVs and 24% of high confidence fusion transcripts. Sensitivity of fusion transcript detection was not overcome by an increase in depth of sequencing up to 3-fold (increase from ~56 to ~159 million reads). Both NanoString and ScriptSeq RNA-Seq technologies yield reliable gene expression data for degraded and FFPE material. The high degree of correlation between NanoString and RNA-Seq platforms suggests discovery based whole transcriptome studies from FFPE material will produce reliable expression data. The RiboZeroGold ScriptSeq protocol performed particularly well for lincRNA expression from FFPE libraries, but detection of eSNV and fusion transcripts was less sensitive.