DiNAMIC.Duo: detecting somatic DNA copy number differences without a normal reference.

MOTIVATION: Somatic DNA copy number alterations (CNAs) arise in tumor tissue because of underlying genomic instability. Recurrent CNAs that occur in the same genomic region across multiple independent samples are of interest to researchers because they may contain genes that contribute to the cancer phenotype. However, differences in copy number states between cancers are also commonly of interest, for example when comparing tumors with distinct morphologies in the same anatomic location. Current methodologies are limited by their inability to perform direct comparisons of CNAs between tumor cohorts, and thus they cannot formally assess the statistical significance of observed copy number differences or identify regions of the genome where these differences occur. RESULTS: We introduce the DiNAMIC.Duo R package that can be used to identify recurrent CNAs in a single cohort or recurrent copy number differences between two cohorts, including when neither cohort is copy neutral. The package utilizes Python scripts for computational efficiency and provides functionality for producing figures and summary output files. AVAILABILITY AND IMPLEMENTATION: The DiNAMIC.Duo R package is available from CRAN at https://cran.r-project.org/web/packages/DiNAMIC.Duo/index.html. This article uses publicly available data from the Broad Institute TCGA Genome Data Analysis Center, https://doi.org/10.7908/C11G0KM9. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Synthesis of α-Diimine Complex Enabling Rapidly Covalent Attachment to Silica Supports and Application of Homo-/Heterogeneous Catalysts in Ethylene Polymerization.

For covalent attachment-supported α-diimine catalysts, on the basis of ensuring the thermal stability and activity of the catalysts, the important problem is that the active group on the catalyst can quickly react with the support, anchoring it firmly on the support, shortening the loading time, reducing the negative impact of the support on the active centers, and further improving the polymer morphology, which makes them suitable for use in industrial polymerization temperatures. Herein, we synthesized a α-diimine nickel(II) catalyst bearing four hydroxyl substituents. The hydroxyl substituents enable the catalyst to be immobilized firmly on silica support by covalent linkage in 5-10 min. Compared with the toluene solvent system, the homogeneous catalysts show high activity and thermal stability in hexane solvent at the same conditions. Compared with homogeneous catalysts, heterogeneous catalysis leads to improvements in catalyst lifetime, polymer morphology control, catalytic activity, and the molecular weight of polyethylene (up to 679 kg/mol). The silica-supported catalysts resulted in higher melting temperatures as well as lower branching densities in polyethylenes. Even at 70 °C, the polyethylene prepared by S-CatA-2 still exhibits dispersed particle morphology, and there is no phenomenon of reactor fouling, which is suitable for industrial polymerization processes.

BRPF3-HBO1 regulates replication origin activation and histone H3K14 acetylation.

During DNA replication, thousands of replication origins are activated across the genome. Chromatin architecture contributes to origin specification and usage, yet it remains unclear which chromatin features impact on DNA replication. Here, we perform a RNAi screen for chromatin regulators implicated in replication control by measuring RPA accumulation upon replication stress. We identify six factors required for normal rates of DNA replication and characterize a function of the bromodomain and PHD finger-containing protein 3 (BRPF3) in replication initiation. BRPF3 forms a complex with HBO1 that specifically acetylates histone H3K14, and genomewide analysis shows high enrichment of BRPF3, HBO1 and H3K14ac at ORC1-binding sites and replication origins found in the vicinity of TSSs. Consistent with this, BRPF3 is necessary for H3K14ac at selected origins and efficient origin activation. CDC45 recruitment, but not MCM2-7 loading, is impaired in BRPF3-depleted cells, identifying a BRPF3-dependent function of HBO1 in origin activation that is complementary to its role in licencing. We thus propose that BRPF3-HBO1 acetylation of histone H3K14 around TSS facilitates efficient activation of nearby replication origins.

An atlas of active enhancers across human cell types and tissues.

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

Swedish Nerve Growth Factor Mutation (NGFR100W) Defines a Role for TrkA and p75NTR in Nociception.

Nerve growth factor (NGF) exerts multiple functions on target neurons throughout development. The recent discovery of a point mutation leading to a change from arginine to tryptophan at residue 100 in the mature NGFß sequence (NGFR100W) in patients with hereditary sensory and autonomic neuropathy type V (HSAN V) made it possible to distinguish the signaling mechanisms that lead to two functionally different outcomes of NGF: trophic versus nociceptive. We performed extensive biochemical, cellular, and live-imaging experiments to examine the binding and signaling properties of NGFR100W Our results show that, similar to the wild-type NGF (wtNGF), the naturally occurring NGFR100W mutant was capable of binding to and activating the TrkA receptor and its downstream signaling pathways to support neuronal survival and differentiation. However, NGFR100W failed to bind and stimulate the 75 kDa neurotrophic factor receptor (p75NTR)-mediated signaling cascades (i.e., the RhoA-Cofilin pathway). Intraplantar injection of NGFR100W into adult rats induced neither TrkA-mediated thermal nor mechanical acute hyperalgesia, but retained the ability to induce chronic hyperalgesia based on agonism for TrkA signaling. Together, our studies provide evidence that NGFR100W retains trophic support capability through TrkA and one aspect of its nociceptive signaling, but fails to engage p75NTR signaling pathways. Our findings suggest that wtNGF acts via TrkA to regulate the delayed priming of nociceptive responses. The integration of both TrkA and p75NTR signaling thus appears to regulate neuroplastic effects of NGF in peripheral nociception.SIGNIFICANCE STATEMENT In the present study, we characterized the naturally occurring nerve growth factor NGFR100W mutant that is associated with hereditary sensory and autonomic neuropathy type V. We have demonstrated for the first time that NGFR100W retains trophic support capability through TrkA, but fails to engage p75NTR signaling pathways. Furthermore, after intraplantar injection into adult rats, NGFR100W induced neither thermal nor mechanical acute hyperalgesia, but retained the ability to induce chronic hyperalgesia. We have also provided evidence that the integration of both TrkA- and p75NTR-mediated signaling appears to regulate neuroplastic effects of NGF in peripheral nociception. Our study with NGFR100W suggests that it is possible to uncouple trophic effect from nociceptive function, both induced by wild-type NGF.

TRiC subunits enhance BDNF axonal transport and rescue striatal atrophy in Huntington's disease.

Corticostriatal atrophy is a cardinal manifestation of Huntington's disease (HD). However, the mechanism(s) by which mutant huntingtin (mHTT) protein contributes to the degeneration of the corticostriatal circuit is not well understood. We recreated the corticostriatal circuit in microfluidic chambers, pairing cortical and striatal neurons from the BACHD model of HD and its WT control. There were reduced synaptic connectivity and atrophy of striatal neurons in cultures in which BACHD cortical and striatal neurons were paired. However, these changes were prevented if WT cortical neurons were paired with BACHD striatal neurons; synthesis and release of brain-derived neurotrophic factor (BDNF) from WT cortical axons were responsible. Consistent with these findings, there was a marked reduction in anterograde transport of BDNF in BACHD cortical neurons. Subunits of the cytosolic chaperonin T-complex 1 (TCP-1) ring complex (TRiC or CCT for chaperonin containing TCP-1) have been shown to reduce mHTT levels. Both CCT3 and the apical domain of CCT1 (ApiCCT1) decreased the level of mHTT in BACHD cortical neurons. In cortical axons, they normalized anterograde BDNF transport, restored retrograde BDNF transport, and normalized lysosomal transport. Importantly, treating BACHD cortical neurons with ApiCCT1 prevented BACHD striatal neuronal atrophy by enhancing release of BDNF that subsequently acts through tyrosine receptor kinase B (TrkB) receptor on striatal neurons. Our findings are evidence that TRiC reagent-mediated reductions in mHTT enhanced BDNF delivery to restore the trophic status of BACHD striatal neurons.

JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net) is the largest open-access database of matrix-based nucleotide profiles describing the binding preference of transcription factors from multiple species. The fifth major release greatly expands the heart of JASPAR-the JASPAR CORE subcollection, which contains curated, non-redundant profiles-with 135 new curated profiles (74 in vertebrates, 8 in Drosophila melanogaster, 10 in Caenorhabditis elegans and 43 in Arabidopsis thaliana; a 30% increase in total) and 43 older updated profiles (36 in vertebrates, 3 in D. melanogaster and 4 in A. thaliana; a 9% update in total). The new and updated profiles are mainly derived from published chromatin immunoprecipitation-seq experimental datasets. In addition, the web interface has been enhanced with advanced capabilities in browsing, searching and subsetting. Finally, the new JASPAR release is accompanied by a new BioPython package, a new R tool package and a new R/Bioconductor data package to facilitate access for both manual and automated methods.

Defective axonal transport of Rab7 GTPase results in dysregulated trophic signaling.

Retrograde trophic signaling of nerve growth factor (NGF) supports neuronal survival and differentiation. Dysregulated trophic signaling could lead to various neurological disorders. Charcot-Marie-Tooth type 2B (CMT2B) is one of the most common inherited peripheral neuropathies characterized by severe terminal axonal loss. Genetic analysis of human CMT2B patients has revealed four missense point mutations in Rab7, a small GTPase that regulates late endosomal/lysosomal pathways, but the exact pathological mechanism remains poorly understood. Here, we show that these Rab7 mutants dysregulated axonal transport and diminished the retrograde signaling of NGF and its TrkA receptor. We found that all CMT2B Rab7 mutants were transported significantly faster than Rab7(wt) in the anterograde direction, accompanied with an increased percentile of anterograde Rab7-vesicles within axons of rat E15.5 dorsal root ganglion (DRG) neurons. In PC12M cells, the CMT2B Rab7 mutants drastically reduced the level of surface TrkA and NGF binding, presumably by premature degradation of TrkA. On the other hand, siRNA knock-down of endogenous Rab7 led to the appearance of large TrkA puncta in enlarged Rab5-early endosomes within the cytoplasm, suggesting delayed TrkA degradation. We also show that CMT2B Rab7 mutants markedly impaired NGF-induced Erk1/2 activation and differentiation in PC12M cells. Further analysis revealed that CMT2B Rab7 mutants caused axonal degeneration in rat E15.5 DRG neurons. We propose that Rab7 mutants induce premature degradation of retrograde NGF-TrkA trophic signaling, which may potentially contribute to the CMT2B disease.

Identification of driver genes in hepatocellular carcinoma by exome sequencing.

UNLABELLED: Genetic alterations in specific driver genes lead to disruption of cellular pathways and are critical events in the instigation and progression of hepatocellular carcinoma (HCC). As a prerequisite for individualized cancer treatment, we sought to characterize the landscape of recurrent somatic mutations in HCC. We performed whole-exome sequencing on 87 HCCs and matched normal adjacent tissues to an average coverage of 59×. The overall mutation rate was roughly two mutations per Mb, with a median of 45 nonsynonymous mutations that altered the amino acid sequence (range, 2-381). We found recurrent mutations in several genes with high transcript levels: TP53 (18%); CTNNB1 (10%); KEAP1 (8%); C16orf62 (8%); MLL4 (7%); and RAC2 (5%). Significantly affected gene families include the nucleotide-binding domain and leucine-rich repeat-containing family, calcium channel subunits, and histone methyltransferases. In particular, the MLL family of methyltransferases for histone H3 lysine 4 were mutated in 20% of tumors. CONCLUSION: The NFE2L2-KEAP1 and MLL pathways are recurrently mutated in multiple cohorts of HCC.

Interaction of antidepressants with the serotonin and norepinephrine transporters: mutational studies of the S1 substrate binding pocket.

The serotonin transporter (SERT) and the norepinephrine transporter (NET) are sodium-dependent neurotransmitter transporters responsible for reuptake of released serotonin and norepinephrine, respectively, into nerve terminals in the brain. A wide range of inhibitors of SERT and NET are used as treatment of depression and anxiety disorders or as psychostimulant drugs of abuse. Despite their clinical importance, the molecular mechanisms by which various types of antidepressant drugs bind and inhibit SERT and NET are still elusive for the majority of the inhibitors, including the molecular basis for SERT/NET selectivity. Mutational analyses have suggested that a central substrate binding site (denoted the S1 pocket) also harbors an inhibitor binding site. In this study, we determine the effect of mutating six key S1 residues in human SERT (hSERT) and NET (hNET) on the potency of 15 prototypical SERT/NET inhibitors belonging to different drug classes. Analysis of the resulting drug sensitivity profiles provides novel information on drug binding modes in hSERT and hNET and identifies specific S1 residues as important molecular determinants for inhibitor potency and hSERT/hNET selectivity.

GMD: measuring the distance between histograms with applications on high-throughput sequencing reads.

UNLABELLED: GMD (generalized minimum distance of distributions) is an R package to assess the similarity between spatial distributions of read-based sequencing data such as ChIP-seq and RNA-seq. GMD calculates the optimal distance between pairs of normalized signal distributions, optionally sliding one distribution over the other to 'align' the distributions. GMD also provides graphical and downstream clustering tools. AVAILABILITY: The R package GMD source code is available at http://cran.r-project.org/web/packages/GMD/ under GPL license.

The mouse ortholog of NEIL3 is a functional DNA glycosylase in vitro and in vivo.

To protect cells from oxidative DNA damage and mutagenesis, organisms possess multiple glycosylases to recognize the damaged bases and to initiate the Base Excision Repair pathway. Three DNA glycosylases have been identified in mammals that are homologous to the Escherichia coli Fpg and Nei proteins, Neil1, Neil2, and Neil3. Neil1 and Neil2 in human and mouse have been well characterized while the properties of the Neil3 protein remain to be elucidated. In this study, we report the characterization of Mus musculus (house mouse) Neil3 (MmuNeil3) as an active DNA glycosylase both in vitro and in vivo. In duplex DNA, MmuNeil3 recognizes the oxidized purines, spiroiminodihydantoin (Sp), guanidinohydantoin (Gh), 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapyG) and 4,6-diamino- 5-formamidopyrimidine (FapyA), but not 8-oxo-7,8-dihydroguanine (8-oxoG). Interestingly, MmuNeil3 prefers lesions in single-stranded DNA and in bubble structures. In contrast to other members of the family that use the N-terminal proline as the nucleophile, MmuNeil3 forms a Schiff base intermediate via its N-terminal valine. We expressed the glycosylase domain of MmuNeil3 (MmuNeil3Delta324) in an Escherichia coli triple mutant lacking Fpg, Nei, and MutY glycosylase activities and showed that MmuNeil3 greatly reduced both the spontaneous mutation frequency and the level of FapyG in the DNA, suggesting that Neil3 plays a role in repairing FapyG in vivo.

JASPAR 2010: the greatly expanded open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net) is the leading open-access database of matrix profiles describing the DNA-binding patterns of transcription factors (TFs) and other proteins interacting with DNA in a sequence-specific manner. Its fourth major release is the largest expansion of the core database to date: the database now holds 457 non-redundant, curated profiles. The new entries include the first batch of profiles derived from ChIP-seq and ChIP-chip whole-genome binding experiments, and 177 yeast TF binding profiles. The introduction of a yeast division brings the convenience of JASPAR to an active research community. As binding models are refined by newer data, the JASPAR database now uses versioning of matrices: in this release, 12% of the older models were updated to improved versions. Classification of TF families has been improved by adopting a new DNA-binding domain nomenclature. A curated catalog of mammalian TFs is provided, extending the use of the JASPAR profiles to additional TFs belonging to the same structural family. The changes in the database set the system ready for more rapid acquisition of new high-throughput data sources. Additionally, three new special collections provide matrix profile data produced by recent alternative high-throughput approaches.

Slurry Homopolymerization of Ethylene Using Thermostable α-Diimine Nickel Catalysts Covalently Linked to Silica Supports via Substituents on Acenaphthequinone-Backbone.

Four supported α-diimine nickel(II) catalysts covalently linked to silica via hydroxyl functionality on α-diimine acenaphthequinone-backbone were prepared and used in slurry polymerizations of ethylene to produce branched polyethylenes. The catalytic activities of these still reached 106 g/molNi·h at 70 °C. The life of the supported catalyst is prolonged, as can be seen from the kinetic profile. The molecular weight of the polyethylene obtained by the 955 silica gel supported catalyst was higher than that obtained by the 2408D silica gel supported catalyst. The melting points of polyethylene obtained by the supported catalysts S-C1-a/b are all above 110 °C. Compared with the homogeneous catalyst, the branching numbers of the polyethylenes obtained by the supported catalysts S-C1-a/b is significantly lower. The polyethylenes obtained by supported catalyst S-C1-a/b at 30-50 °C are free-flowing particles, which is obviously better than the rubber-like cluster polymer obtained from homogeneous catalyst.

IL-2K35C-moFA, a Long-Acting Engineered Cytokine with Decreased Interleukin 2 Receptor α Binding, Improved the Cellular Selectivity Profile and Antitumor Efficacy in a Mouse Tumor Model.

Human interleukin 2 (IL-2) has shown impressive results as a therapeutic agent for cancer. However, IL-2-based cancer therapy is limited by strong Treg amplification owing to its high binding affinity to IL-2 receptor α (IL-2Rα) and its short half-life owing to its small molecular size. In this study, we solved these problems using a covalent modification strategy of the IL-2 variant, i.e., substituting cysteine (C) for lysine (K) at position 35, using octadecanedicarboxylic acid through maleimide chemistry, creating IL-2K35C-moFA. IL-2K35C-moFA was equipotent to human IL-2 wild type (IL-2WT) in activating tumor-killing CD8+ memory effector T cells (CD8+ T) and NK cells bearing the intermediate affinity IL-2 receptors, and less potent than IL-2WT on CTLL-2 cells bearing the high-affinity IL-2 receptors. Moreover, it was shown to support the preferential activation of IL-2 receptor ß (IL-2Rß) over IL-2Rα because of the mutation and fatty acid conjugation. In a B16F10 murine tumor model, IL-2K35C-moFA showed efficacy as a single dose and provided durable immunity for 1 week. Our results support the further evaluation of IL-2K35C-moFA as a novel cancer immunotherapy.

Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data.

BACKGROUND: Initiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models. RESULTS: Here, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling. CONCLUSIONS: In this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.

Design, synthesis and structure-activity relationship (SAR) studies of 2,4-disubstituted pyrimidine derivatives: dual activity as cholinesterase and Aß-aggregation inhibitors.

A novel class of 2,4-disubstituted pyrimidines (7a-u, 8a-f, 9a-e) that possess substituents with varying steric and electronic properties at the C-2 and C-4 positions, were designed, synthesized and evaluated as dual cholinesterase and amyloid-ß (Aß)-aggregation inhibitors. In vitro screening identified N-(naphth-1-ylmethyl)-2-(pyrrolidin-1-yl)pyrimidin-4-amine (9a) as the most potent AChE inhibitor (IC(50)=5.5 µM). Among this class of compounds, 2-(4-methylpiperidin-1-yl)-N-(naphth-1-ylmethyl)pyrimidin-4-amine (9e) was identified as the most potent and selective BuChE inhibitor (IC(50)=2.2 µM, selectivity index=11.7) and was about 5.7-fold more potent compared to the commercial, approved reference drug galanthamine (BuChE IC(50)=12.6 µM). In addition, the selective AChE inhibitor N-benzyl-2-(4-methylpiperazin-1-yl)pyrimidin-4-amine (7d), exhibited good inhibition of hAChE-induced aggregation of Aß(1-40) fibrils (59% inhibition). Furthermore, molecular modeling studies indicate that a central pyrimidine ring serves as a suitable template to develop dual inhibitors of cholinesterase and AChE-induced Aß aggregation thereby targeting multiple pathological routes in AD.

SCISSOR: a framework for identifying structural changes in RNA transcripts.

High-throughput sequencing protocols such as RNA-seq have made it possible to interrogate the sequence, structure and abundance of RNA transcripts at higher resolution than previous microarray and other molecular techniques. While many computational tools have been proposed for identifying mRNA variation through differential splicing/alternative exon usage, challenges in its analysis remain. Here, we propose a framework for unbiased and robust discovery of aberrant RNA transcript structures using short read sequencing data based on shape changes in an RNA-seq coverage profile. Shape changes in selecting sample outliers in RNA-seq, SCISSOR, is a series of procedures for transforming and normalizing base-level RNA sequencing coverage data in a transcript independent manner, followed by a statistical framework for its analysis ( https://github.com/hyochoi/SCISSOR ). The resulting high dimensional object is amenable to unsupervised screening of structural alterations across RNA-seq cohorts with nearly no assumption on the mutational mechanisms underlying abnormalities. This enables SCISSOR to independently recapture known variants such as splice site mutations in tumor suppressor genes as well as novel variants that are previously unrecognized or difficult to identify by any existing methods including recurrent alternate transcription start sites and recurrent complex deletions in 3' UTRs.

UNMASC: tumor-only variant calling with unmatched normal controls.

Despite years of progress, mutation detection in cancer samples continues to require significant manual review as a final step. Expert review is particularly challenging in cases where tumors are sequenced without matched normal control DNA. Attempts have been made to call somatic point mutations without a matched normal sample by removing well-known germline variants, utilizing unmatched normal controls, and constructing decision rules to classify sequencing errors and private germline variants. With budgetary constraints related to computational and sequencing costs, finding the appropriate number of controls is a crucial step to identifying somatic variants. Our approach utilizes public databases for canonical somatic variants as well as germline variants and leverages information gathered about nearby positions in the normal controls. Drawing from our cohort of targeted capture panel sequencing of tumor and normal samples with varying tumortypes and demographics, these served as a benchmark for our tumor-only variant calling pipeline to observe the relationship between our ability to correctly classify variants against a number of unmatched normals. With our benchmarked samples, approximately ten normal controls were needed to maintain 94% sensitivity, 99% specificity and 76% positive predictive value, far outperforming comparable methods. Our approach, called UNMASC, also serves as a supplement to traditional tumor with matched normal variant calling workflows and can potentially extend to other concerns arising from analyzing next generation sequencing data.

Integrative Analysis of miRNAs Identifies Clinically Relevant Epithelial and Stromal Subtypes of Head and Neck Squamous Cell Carcinoma.

PURPOSE: The objective of this study is to characterize the role of miRNAs in the classification of head and neck squamous cell carcinoma (HNSCC). EXPERIMENTAL DESIGN: Here, we analyzed 562 HNSCC samples, 88 from a novel cohort and 474 from The Cancer Genome Atlas, using miRNA microarray and miRNA sequencing, respectively. Using an integrative correlations method followed by miRNA expression-based hierarchical clustering, we validated miRNA clusters across cohorts. Evaluation of clusters by logistic regression and gene ontology approaches revealed subtype-based clinical and biological characteristics. RESULTS: We identified two independently validated and statistically significant (P < 0.01) tumor subtypes and named them "epithelial" and "stromal" based on associations with functional target gene ontology relating to differing stages of epithelial cell differentiation. miRNA-based subtypes were correlated with individual gene expression targets based on miRNA seed sequences, as well as with miRNA families and clusters including the miR-17 and miR-200 families. These correlated genes defined pathways relevant to normal squamous cell function and pathophysiology. miRNA clusters statistically associated with differential mutation patterns including higher proportions of TP53 mutations in the stromal class and higher NSD1 and HRAS mutation frequencies in the epithelial class. miRNA classes correlated with previously reported gene expression subtypes, clinical characteristics, and clinical outcomes in a multivariate Cox proportional hazards model with stromal patients demonstrating worse prognoses (HR, 1.5646; P = 0.006). CONCLUSIONS: We report a reproducible classification of HNSCC based on miRNA that associates with known pathologically altered pathways and mutations of squamous tumors and is clinically relevant.