Search

DNA methylation is an important epigenetic mark that is involved in the regulation of many cellular processes such as gene expression, genomic imprinting and silencing of repetitive elements. Because of their ability to cause and capture phenotypic plasticity, epigenetic marks such as DNA methylation represent potential biomarkers to distinguish between different types of tissues and stages of differentiation. Here, we have identified differential DNA methylation in the gene body of the nitric oxide inhibitor Ddah2 that discriminates embryonic stem cells from neural stem cells and is positively correlated with differential gene expression.

The magnetic solids Cu(terpy)Mo2O7 (terpy = terpyridine) and Cu(OH)(p-pyc)H2O (p-pyc = p-pyridinecarboxylate) have a spin gap and possess chains of Cu2+ ions in which two different Cu...Cu distances alternate. On the basis of their reported crystal structures, the spin-exchange interactions of these compounds were examined by performing spin dimer analysis to determine whether an antiferromagnetic dimer or an alternating antiferromagnetic chain model is appropriate for their magnetic properties. Our analysis shows that an antiferromagnetic dimer model is correct for both compounds because of the anisotropic overlap between the magnetic orbitals of their Cu2+ sites.

Other than an association with HPV infection, little is known about the genetic alterations determining the development of penile cancer. Although penile cancer is rare in the developed world, it presents a significant burden in developing countries. Here, we report the findings of whole-exome sequencing (WES) to determine the somatic mutational landscape of penile cancer. WES was performed on penile cancer and matched germline DNA from 27 patients undergoing surgical resection. Targeted resequencing of candidate genes was performed in an independent 70 patient cohort. Mutation data were also integrated with DNA methylation and copy-number information from the same patients. We identified an HPV-associated APOBEC mutation signature and an NpCpG signature in HPV-negative disease. We also identified recurrent mutations in the novel penile cancer tumor suppressor genes CSN1(GPS1) and FAT1. Expression of CSN1 mutants in cells resulted in colocalization with AGO2 in cytoplasmic P-bodies, ultimately leading to the loss of miRNA-mediated gene silencing, which may contribute to disease etiology. Our findings represent the first comprehensive analysis of somatic alterations in penile cancer, highlighting the complex landscape of alterations in this malignancy. Cancer Res; 76(16); 4720–7. ©2016 AACR.

Summary Glioblastoma (GBM) is an aggressive brain tumor whose growth is driven by stem cell-like cells. BMP signaling triggers cell-cycle exit and differentiation of GBM stem cells (GSCs) and, therefore, might have therapeutic value. However, the epigenetic mechanisms that accompany differentiation remain poorly defined. It is also unclear whether cell-cycle arrest is terminal. Here we find only a subset of GSC cultures exhibit astrocyte differentiation in response to BMP. Although overtly differentiated non-cycling astrocytes are generated, they remain vulnerable to cell-cycle re-entry and fail to appropriately reconfigure DNA methylation patterns. Chromatin accessibility mapping identified loci that failed to alter in response to BMP and these were enriched in SOX transcription factor-binding motifs. SOX transcription factors, therefore, may limit differentiation commitment. A similar propensity for cell-cycle re-entry and de-differentiation was observed in GSC-derived oligodendrocyte-like cells. These findings highlight significant obstacles to BMP-induced differentiation as therapy for GBM., Graphical Abstract, Highlights • Genome-wide profiling shows DNA methylation patterns during glioblastoma (GBM) differentiation • Delayed and incomplete epigenetic changes appear in GBM stem cells in response to BMP • SOX transcription factors may explain the lack of terminal differentiation • Lack of differentiation commitment limits the effectiveness of BMP-based therapies, BMP induces differentiation of glioblastoma stem cells (GSCs), but it remains unclear if differentiation commitment and permanent cell-cycle arrest occurs. Pollard, Beck, and colleagues report that differentiated progeny of GSCs fail to reconfigure DNA methylation patterns and are vulnerable to de-differentiation. Failure to suppress the activity of SOX transcription factors may explain this deficit.

Background A substantial proportion of cancer cases present with a metastatic tumor and require further testing to determine the primary site; many of these are never fully diagnosed and remain cancer of unknown primary origin (CUP). It has been previously demonstrated that the somatic point mutations detected in a tumor can be used to identify its site of origin with limited accuracy. We hypothesized that higher accuracy could be achieved by a classification algorithm based on the following feature sets: 1) the number of nonsynonymous point mutations in a set of 232 specific cancer-associated genes, 2) frequencies of the 96 classes of single-nucleotide substitution determined by the flanking bases, and 3) copy number profiles, if available. Methods We used publicly available somatic mutation data from the COSMIC database to train random forest classifiers to distinguish among those tissues of origin for which sufficient data was available. We selected feature sets using cross-validation and then derived two final classifiers (with or without copy number profiles) using 80 % of the available tumors. We evaluated the accuracy using the remaining 20 %. For further validation, we assessed accuracy of the without-copy-number classifier on three independent data sets: 1669 newly available public tumors of various types, a cohort of 91 breast metastases, and a set of 24 specimens from 9 lung cancer patients subjected to multiregion sequencing. Results The cross-validation accuracy was highest when all three types of information were used. On the left-out COSMIC data not used for training, we achieved a classification accuracy of 85 % across 6 primary sites (with copy numbers), and 69 % across 10 primary sites (without copy numbers). Importantly, a derived confidence score could distinguish tumors that could be identified with 95 % accuracy (32 %/75 % of tumors with/without copy numbers) from those that were less certain. Accuracy in the independent data sets was 46 %, 53 % and 89 % respectively, similar to the accuracy expected from the training data. Conclusions Identification of primary site from point mutation and/or copy number data may be accurate enough to aid clinical diagnosis of cancers of unknown primary origin. Electronic supplementary material The online version of this article (doi:10.1186/s12920-015-0130-0) contains supplementary material, which is available to authorized users.

Esophageal adenocarcinomas are associated with a dismal prognosis. Deciphering the evolutionary history of this disease may shed light on therapeutically tractable targets and reveal dynamic mutational processes during the disease course and following neoadjuvant chemotherapy (NAC). We exome sequenced 40 tumor regions from 8 patients with operable esophageal adenocarcinomas, before and after platinum-containing NAC. This revealed the evolutionary genomic landscape of esophageal adenocarcinomas with the presence of heterogeneous driver mutations, parallel evolution, early genome-doubling events, and an association between high intratumor heterogeneity and poor response to NAC. Multiregion sequencing demonstrated a significant reduction in thymine to guanine mutations within a CpTpT context when comparing early and late mutational processes and the presence of a platinum signature with enrichment of cytosine to adenine mutations within a CpC context following NAC. Esophageal adenocarcinomas are characterized by early chromosomal instability leading to amplifications containing targetable oncogenes persisting through chemotherapy, providing a rationale for future therapeutic approaches. Significance: This work illustrates dynamic mutational processes occurring during esophageal adenocarcinoma evolution and following selective pressures of platinum exposure, emphasizing the iatrogenic impact of therapy on cancer evolution. Identification of amplifications encoding targetable oncogenes maintained through NAC suggests the presence of stable vulnerabilities, unimpeded by cytotoxics, suitable for therapeutic intervention. Cancer Discov; 5(8); 821–31. ©2015 AACR. See related commentary by Devarakonda and Govindan, p. 796. This article is highlighted in the In This Issue feature, p. 783

The genomic revolution is fully upon us in 2006 and the pace of discovery is set to accelerate with the emergence of ultra-high-throughput sequencing technologies. Our complete genome collection of bacteria and archaea continues to grow in number and diversity, as genome sequencing is applied to an array of new problems, from the characterization of the pan-genome to the detection of mutation after experimentation and the exploration of microbial communities in unprecedented detail. The benefits of large-scale comparative genomic analyses are driving the community to think about how to manage our public collections of genomes in novel ways.

The explosion in the number of complete genomes over the past decade has spawned a new and exciting discipline, that of comparative genomics. To exploit the full potential of this approach requires the development of novel algorithms, databases and software which are sophisticated enough to draw meaningful comparisons between complete genome sequences and are widely accessible to the scientific community at large. This article reviews progress towards the development of computational tools and databases for organizing and extracting biological meaning from the comparison of large collections of genomes.

Defining mechanisms that generate intratumour heterogeneity and branched evolution may inspire novel therapeutic approaches to limit tumour diversity and adaptation. SETD2 (Su(var), Enhancer of zeste, Trithorax-domain containing 2) trimethylates histone-3 lysine-36 (H3K36me3) at sites of active transcription and is mutated in diverse tumour types, including clear cell renal carcinomas (ccRCCs). Distinct SETD2 mutations have been identified in spatially separated regions in ccRCC, indicative of intratumour heterogeneity. In this study, we have addressed the consequences of SETD2 loss-of-function through an integrated bioinformatics and functional genomics approach. We find that bi-allelic SETD2 aberrations are not associated with microsatellite instability in ccRCC. SETD2 depletion in ccRCC cells revealed aberrant and reduced nucleosome compaction and chromatin association of the key replication proteins minichromosome maintenance complex component (MCM7) and DNA polymerase δ hindering replication fork progression, and failure to load lens epithelium-derived growth factor and the Rad51 homologous recombination repair factor at DNA breaks. Consistent with these data, we observe chromosomal breakpoint locations are biased away from H3K36me3 sites in SETD2 wild-type ccRCCs relative to tumours with bi-allelic SETD2 aberrations and that H3K36me3-negative ccRCCs display elevated DNA damage in vivo. These data suggest a role for SETD2 in maintaining genome integrity through nucleosome stabilization, suppression of replication stress and the coordination of DNA repair.

Background Common human diseases are caused by the complex interplay of genetic susceptibility as well as environmental factors. Due to the environment’s influence on the epigenome, and therefore genome function, as well as conversely the genome’s facilitative effect on the epigenome, analysis of this level of regulation may increase our knowledge of disease pathogenesis. Methods In order to identify human-specific epigenetic influences, we have performed a novel genome-wide DNA methylation analysis comparing human, chimpanzee and rhesus macaque. Results We have identified that the immunological Leukotriene B4 receptor (LTB4R, BLT1 receptor) is the most epigenetically divergent human gene in peripheral blood in comparison with other primates. This difference is due to the co-ordinated active state of human-specific hypomethylation in the promoter and human-specific increased gene body methylation. This gene is significant in innate immunity and the LTB4/LTB4R pathway is involved in the pathogenesis of the spectrum of human inflammatory diseases. This finding was confirmed by additional neutrophil-only DNA methylome and lymphoblastoid H3K4me3 chromatin comparative data. Additionally we show through functional analysis that this receptor has increased expression and a higher response to the LTB4 ligand in human versus rhesus macaque peripheral blood mononuclear cells. Genome-wide we also find human species-specific differentially methylated regions (human s-DMRs) are more prevalent in CpG island shores than within the islands themselves, and within the latter are associated with the CTCF motif. Conclusions This result further emphasises the exclusive nature of the human immunological system, its divergent adaptation even from very closely related primates, and the power of comparative epigenomics to identify and understand human uniqueness.

Isocitrate dehydrogenase (IDH) genes 1 and 2 are frequently mutated in acute myeloid leukaemia (AML), low-grade glioma, cholangiocarcinoma (CC) and chondrosarcoma (CS). For AML, low-grade glioma and CC, mutant IDH status is associated with a DNA hypermethylation phenotype, implicating altered epigenome dynamics in the aetiology of these cancers. Here we show that the IDH variants in CS are also associated with a hypermethylation phenotype and display increased production of the oncometabolite 2-hydroxyglutarate, supporting the role of mutant IDH-produced 2-hydroxyglutarate as an inhibitor of TET-mediated DNA demethylation. Meta-analysis of the acute myeloid leukaemia, low-grade glioma, cholangiocarcinoma and CS methylation data identifies cancer-specific effectors within the retinoic acid receptor activation pathway among the hypermethylated targets. By analysing sequence motifs surrounding hypermethylated sites across the four cancer types, and using chromatin immunoprecipitation and western blotting, we identify the transcription factor EBF1 (early B-cell factor 1) as an interaction partner for TET2, suggesting a sequence-specific mechanism for regulating DNA methylation., Cancer-associated mutations in isocitrate dehydrogenase are proposed to impair TET2-dependent DNA demethylation. By comparing the methylomes of IDH-mutant cancers, the authors identify the transcription factor EBF1 as a partner of TET2, suggesting a possible means for targeting TET2 to specific DNA sequences.

Background Methylated DNA immunoprecipitation (MeDIP) is a popular enrichment based method and can be combined with sequencing (termed MeDIP-seq) to interrogate the methylation status of cytosines across entire genomes. However, quality control and analysis of MeDIP-seq data have remained to be a challenge. Results We report genome-wide DNA methylation profiles of wild type (wt) and mutant mouse cells, comprising 3 biological replicates of Thymine DNA glycosylase (Tdg) knockout (KO) embryonic stem cells (ESCs), in vitro differentiated neural precursor cells (NPCs) and embryonic fibroblasts (MEFs). The resulting 18 methylomes were analysed with MeDUSA (Methylated DNA Utility for Sequence Analysis), a novel MeDIP-seq computational analysis pipeline for the identification of differentially methylated regions (DMRs). The observed increase of hypermethylation in MEF promoter-associated CpG islands supports a previously proposed role for Tdg in the protection of regulatory regions from epigenetic silencing. Further analysis of genes and regions associated with the DMRs by gene ontology, pathway, and ChIP analyses revealed further insights into Tdg function, including an association of TDG with low-methylated distal regulatory regions. Conclusions We demonstrate that MeDUSA is able to detect both large-scale changes between cells from different stages of differentiation and also small but significant changes between the methylomes of cells that only differ in the KO of a single gene. These changes were validated utilising publicly available datasets and confirm TDG's function in the protection of regulatory regions from epigenetic silencing.

Objective: A significant proportion of HNSCC is caused by human papillomavirus. Patients diagnosed with this subtype have a significantly better outcome. It is the aim of the UCL Head and Neck Cancer Epigenome Project to generate and analyze comprehensive, genome-wide profiles of epigenetic aberrations in this specific type of human cancer.Method: Whole-genome DNA methylation (methylome) analysis of 3 HPV positive and 3 HPV negative frozen HNSCC samples (obtained from the UCL Head and Neck Cancer Tissue Bank as part of a larger discovery set) was conducted using methylated DNA immunoprecipitation coupled with next-generation sequencing (MeDIP-Seq).Results: On average, each sample in the pilot phase generated 2 x 107 sequencing reads and 70% coverage of all CpG sites within the host genome. HPV type 16 was observed in all HPV positive samples and, interestingly, methylated regions of this virus were identified—most notably at genes E1 and L1 (previously shown to have been methylated in cervical cancer). We...

Aberrant DNA methylation (DNAm) was first linked to cancer over 25 yr ago. Since then, many studies have associated hypermethylation of tumor suppressor genes and hypomethylation of oncogenes to the tumorigenic process. However, most of these studies have been limited to the analysis of promoters and CpG islands (CGIs). Recently, new technologies for whole-genome DNAm (methylome) analysis have been developed, enabling unbiased analysis of cancer methylomes. By using MeDIP-seq, we report a sequencing-based comparative methylome analysis of malignant peripheral nerve sheath tumors (MPNSTs), benign neurofibromas, and normal Schwann cells. Analysis of these methylomes revealed a complex landscape of DNAm alterations. In contrast to what has been reported for other tumor types, no significant global hypomethylation was observed in MPNSTs using methylome analysis by MeDIP-seq. However, a highly significant (P < 10−100) directional difference in DNAm was found in satellite repeats, suggesting these repeats to be the main target for hypomethylation in MPNSTs. Comparative analysis of the MPNST and Schwann cell methylomes identified 101,466 cancer-associated differentially methylated regions (cDMRs). Analysis showed these cDMRs to be significantly enriched for two satellite repeat types (SATR1 and ARLα) and suggests an association between aberrant DNAm of these sequences and transition from healthy cells to malignant disease. Significant enrichment of hypermethylated cDMRs in CGI shores (P < 10−60), non–CGI-associated promoters (P < 10−4) and hypomethylated cDMRs in SINE repeats (P < 10−100) was also identified. Integration of DNAm and gene expression data showed that the expression pattern of genes associated with CGI shore cDMRs was able to discriminate between disease phenotypes. This study establishes MeDIP-seq as an effective method to analyze cancer methylomes.

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10-4, permutation p = 1.0×10-3). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10-7). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within diseaseassociated loci, thus providing a novel route to aid unravelling common complex diseases. © 2010 Bell et al.

With the quantity of genomic data increasing at an exponential rate, it is imperative that these data be captured electronically, in a standard format. Standardization activities must proceed within the auspices of open-access and international working bodies. To tackle the issues surrounding the development of better descriptions of genomic investigations, we have formed the Genomic Standards Consortium (GSC). Here, we introduce the minimum information about a genome sequence (MIGS) specification with the intent of promoting participation in its development and discussing the resources that will be required to develop improved mechanisms of metadata capture and exchange. As part of its wider goals, the GSC also supports improving the 'transparency' of the information contained in existing genomic databases.

Density functional calculations are used to obtain the electronic structure of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{SrRh}}_{2}{\mathrm{O}}_{4}$ in comparison with spinel ${\mathrm{ZnRh}}_{2}{\mathrm{O}}_{4}$. Both materials are band insulators, with substantial crystal field induced band gaps, reflecting strong transition-metal--O hybridization. However, due to the bonding topology in these materials, the valence bands are very narrow. This leads to high thermopowers within standard Boltzmann transport theory, and indicates that they can be the basis of good thermoelectric materials provided that they can be doped into metallic states with reasonable carrier mobility. In the case of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{SrRh}}_{2}{\mathrm{O}}_{4}$, scattering due to Sr disorder is important. Also, again in $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{SrRh}}_{2}{\mathrm{O}}_{4}$, the band gap may be large enough to be of interest for photoelectrochemical ${\mathrm{H}}_{2}$ production.

Background Lineage-specific, or taxonomically restricted genes (TRGs), especially those that are species and strain-specific, are of special interest because they are expected to play a role in defining exclusive ecological adaptations to particular niches. Despite this, they are relatively poorly studied and little understood, in large part because many are still orphans or only have homologues in very closely related isolates. This lack of homology confounds attempts to establish the likelihood that a hypothetical gene is expressed and, if so, to determine the putative function of the protein. Methodology/Principal Findings We have developed “QIPP” (“Quality Index for Predicted Proteins”), an index that scores the “quality” of a protein based on non-homology-based criteria. QIPP can be used to assign a value between zero and one to any protein based on comparing its features to other proteins in a given genome. We have used QIPP to rank the predicted proteins in the proteomes of Bacteria and Archaea. This ranking reveals that there is a large amount of variation in QIPP scores, and identifies many high-scoring orphans as potentially “authentic” (expressed) orphans. There are significant differences in the distributions of QIPP scores between orphan and non-orphan genes for many genomes and a trend for less well-conserved genes to have lower QIPP scores. Conclusions The implication of this work is that QIPP scores can be used to further annotate predicted proteins with information that is independent of homology. Such information can be used to prioritize candidates for further analysis. Data generated for this study can be found in the OrphanMine at http://www.genomics.ceh.ac.uk/orphan_min​e.

Introduction Increasing evidence supports the heterogeneous clonal composition and branched evolution of non-small cell lung cancer (NSCLC) tumors. Circulating biomarkers, such as circulating cell-free tumor DNA (cfDNA), may potentially be used to track the evolution of such tumors in order to detect minimal residual disease, monitor treatment response, detect the emergence of resistance, and subsequent disease relapse. Experimental procedures In a cohort of 14 primary NSCLC tumors, DNA was extracted from fresh frozen tissue and multi-region whole-exome sequencing (mean depth of 219x) was performed on a total of 44 tumor regions. Single nucleotide variants, indels and copy number data were identified using VarScan2 (v2.3.6). Non-silent mutations were validated using ultra-deep amplicon sequencing (mean depth 1220x). The subclonal composition of each tumor was used to construct phylogenetic trees. Candidate driver genes were selected by referencing large-scale lung cancer or pan-cancer sequencing studies, and the COSMIC gene census. To validate the identified driver mutations and to further screen all tumor regions for additional low-level heterogeneity, we developed patient-specific UltraSEEK panels. The UltraSEEK technology utilizes mutation-specific single base extension and MassARRAY® MALDI-TOF Mass Spectrometry to identify low frequency mutated alleles down to 0.1%. This method was used to determine if potential driver mutations could be detected in cfDNA extracted at the time of diagnosis in a subset of patients. Summary of data Intratumor heterogeneity was evident in all tumors, with a median of 25% heterogeneous mutations (range 5-61%). Phylogenetic tree analysis revealed branched tumor evolution, with mutations in potential driver genes occurring both early (trunk) and late (branch). The UltraSEEK method identified truncal and branch mutations in tumor DNA, as well as truncal mutations in the EGFR, TP53, BRAF, RB1 and MLL2 genes in two out of six patients. These mutations are predicted to have deleterious functional consequences and may therefore play a role in disease pathogenesis. Conclusions Spatial and temporal heterogeneity of driver mutations was evident in all NSCLC tumors. The use of liquid biopsies, such as cfDNA, has the potential to aid early detection and guide treatment initiation in cancer patients. They may also avoid the need for repeated tissue sampling, and be of value in cases where disease sites, such as brain or bone, are not easily accessible to biopsy. Citation Format: Mariam Jamal-Hanjani, Elza De Bruin, Gareth Wilson, Nicholas McGranahan, Andrew Rowan, Michael Mosko, Heath Metzler, Anders Nygren, Charles Swanton. Intratumor heterogeneity and the detection of potential driver mutations in cfDNA in a NSCLC cohort using UltraSEEKTM. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 5249. doi:10.1158/1538-7445.AM2015-5249

Aim: Few studies have investigated the clonal architecture of non-small cell lung cancer (NSCLC) and its evolution over time. With increasing evidence for the potential impact of intratumour heterogeneity (ITH) on clinical outcome, and its recognition as a driver of cancer progression, there is an unmet need to elicit the implications of ITH for biomarker validation, therapeutic response, and therefore survival. Methods: In a cohort of 11 primary NSCLC tumours, multi-region whole-exome sequencing at a mean depth of 214x was performed on a total of 37 tumour regions. Single nucleotide variants, indels and copy number data were identified using VarScan2 (v2.3.6). Non-silent mutations were validated using ultra-deep amplicon sequencing (mean depth 1189x). The subclonal composition of each tumor was used to construct phylogenetic trees. Candidate driver genes were selected by referencing large-scale lung cancer or pan-cancer sequencing studies, and the Catalogue of Somatic Mutations in Cancer (COSMIC) gene census. Results: ITH was evident in all tumours, with a median of 33% heterogeneous mutations (range 5–61%). Phylogenetic tree analysis revealed branched tumour evolution, with mutations in potential driver genes occurring both early (trunk) and late (branch). Thirty-five percent (range 8-67%) of mutations in potential driver genes were heterogeneous, and included TP53, EGFR, KRAS, BRAF, ALK and ROS1. Overall, heterogeneous non-silent somatic mutations and copy number aberrations occurred frequently. Conclusions: This study demonstrates spatial and temporal heterogeneity of both driver mutations and copy number events, suggesting that NSCLCs may follow multiple distinct evolutionary pathways simultaneously. Knowledge of these pathways may help identify novel therapeutic targets to improve clinical outcomes, and large-scale longitudinal genomic studies, such as TRACERx (TRAcking NSCLC Evolution through therapy (Rx), NCT01888601), may become a central component to the delivery of precision cancer medicine. Disclosure: All authors have declared no conflicts of interest.

The integration of genomic and epigenomic data is an increasingly popular approach for studying the complex mechanisms driving cancer development. We have developed a method for evaluating both methylation and copy number from high-density DNA methylation arrays. Comparing copy number data from Infinium HumanMethylation450 BeadChips and SNP arrays, we demonstrate that Infinium arrays detect copy number alterations with the sensitivity of SNP platforms. These results show that high-density methylation arrays provide a robust and economic platform for detecting copy number and methylation changes in a single experiment. Our method is available in the ChAMP Bioconductor package: http://www.bioconductor.org/packages/2.13/bioc/html/ChAMP.html.

Background Human papillomavirus-positive (HPV+) head and neck squamous cell carcinoma (HNSCC) represents a distinct clinical and epidemiological condition compared with HPV-negative (HPV-) HNSCC. To test the possible involvement of epigenetic modulation by HPV in HNSCC, we conducted a genome-wide DNA-methylation analysis. Methods Using laser-capture microdissection of 42 formalin-fixed paraffin wax-embedded (FFPE) HNSCCs, we generated DNA-methylation profiles of 18 HPV+ and 14 HPV- samples, using Infinium 450 k BeadArray technology. Methylation data were validated in two sets of independent HPV+/HPV- HNSCC samples (fresh-frozen samples and cell lines) using two independent methods (Infinium 450 k and whole-genome methylated DNA immunoprecipitation sequencing (MeDIP-seq)). For the functional analysis, an HPV- HNSCC cell line was transduced with lentiviral constructs containing the two HPV oncogenes (E6 and E7), and effects on methylation were assayed using the Infinium 450 k technology. Results and discussion Unsupervised clustering over the methylation variable positions (MVPs) with greatest variation showed that samples segregated in accordance with HPV status, but also that HPV+ tumors are heterogeneous. MVPs were significantly enriched at transcriptional start sites, leading to the identification of a candidate CpG island methylator phenotype in a sub-group of the HPV+ tumors. Supervised analysis identified a strong preponderance (87%) of MVPs towards hypermethylation in HPV+ HNSCC. Meta-analysis of our HNSCC and publicly available methylation data in cervical and lung cancers confirmed the observed DNA-methylation signature to be HPV-specific and tissue-independent. Grouping of MVPs into functionally more significant differentially methylated regions identified 43 hypermethylated promoter DMRs, including for three cadherins of the Polycomb group target genes. Integration with independent expression data showed strong negative correlation, especially for the cadherin gene-family members. Combinatorial ectopic expression of the two HPV oncogenes (E6 and E7) in an HPV- HNSCC cell line partially phenocopied the hypermethylation signature seen in HPV+ HNSCC tumors, and established E6 as the main viral effector gene. Conclusions Our data establish that archival FFPE tissue is very suitable for this type of methylome analysis, and suggest that HPV modulates the HNSCC epigenome through hypermethylation of Polycomb repressive complex 2 target genes such as cadherins, which are implicated in tumor progression and metastasis.

Type A and B influenza viruses can cause a wide spectrum of illness, and these viruses are responsible for considerable mortality and morbidity. Rapid typing of isolates is desirable when amantadine treatment or prophylaxis of contacts of type A influenza virus carriers is considered, but the available rapid techniques lack sensitivity and standard diagnostic methods require expansion of virus in tissue culture or embryonated hens' eggs. We developed a series of oligonucleotide primers able to detect, type, and subtype type A influenza viruses in a single reverse transcription-PCR. RNA was isolated from clinical specimens, and cDNA was generated with random primers. PCR was carried out with a mixture of primers specific for influenza viruses of types B, A/H1 and A/H3, and subtyping of the neuraminidase was carried out on the same cDNA template under identical conditions. Amplified products were detected by ethidium bromide staining of amplified products after agarose gel electrophoresis. When it was used to test 98 clinical specimens, this method was comparable to standard culture techniques in the detection, typing, and subtyping of influenza viruses.

Tumour mutational burden (TMB) predicts immunotherapy outcome in non-small cell lung cancer (NSCLC), consistent with immune recognition of tumour neoantigens. However, persistent antigen exposure is detrimental for T cell function. How TMB affects CD4 and CD8 T cell differentiation in untreated tumours, and whether this affects patient outcomes is unknown. Here we paired high-dimensional flow cytometry, exome, single-cell and bulk RNA sequencing from patients with resected, untreated NSCLC to examine these relationships. TMB was associated with compartment-wide T cell differentiation skewing, characterized by loss of TCF7-expressing progenitor-like CD4 T cells, and an increased abundance of dysfunctional CD8 and CD4 T cell subsets, with significant phenotypic and transcriptional similarity to neoantigen-reactive CD8 T cells. A gene signature of redistribution from progenitor-like to dysfunctional states associated with poor survival in lung and other cancer cohorts. Single-cell characterization of these populations informs potential strategies for therapeutic manipulation in NSCLC.