Genomic mutation landscape of skin cancers from DNA repair-deficient xeroderma pigmentosum patients.

Xeroderma pigmentosum (XP) is a genetic disorder caused by mutations in genes of the Nucleotide Excision Repair (NER) pathway (groups A-G) or in Translesion Synthesis DNA polymerase η (V). XP is associated with an increased skin cancer risk, reaching, for some groups, several thousand-fold compared to the general population. Here, we analyze 38 skin cancer genomes from five XP groups. We find that the activity of NER determines heterogeneity of the mutation rates across skin cancer genomes and that transcription-coupled NER extends beyond the gene boundaries reducing the intergenic mutation rate. Mutational profile in XP-V tumors and experiments with POLH knockout cell line reveal the role of polymerase η in the error-free bypass of (i) rare TpG and TpA DNA lesions, (ii) 3' nucleotides in pyrimidine dimers, and (iii) TpT photodimers. Our study unravels the genetic basis of skin cancer risk in XP and provides insights into the mechanisms reducing UV-induced mutagenesis in the general population.

Integrative Pan-Cancer Genomic and Transcriptomic Analyses of Refractory Metastatic Cancer.

Metastatic relapse after treatment is the leading cause of cancer mortality, and known resistance mechanisms are missing for most treatments administered to patients. To bridge this gap, we analyze a pan-cancer cohort (META-PRISM) of 1,031 refractory metastatic tumors profiled via whole-exome and transcriptome sequencing. META-PRISM tumors, particularly prostate, bladder, and pancreatic types, displayed the most transformed genomes compared with primary untreated tumors. Standard-of-care resistance biomarkers were identified only in lung and colon cancers-9.6% of META-PRISM tumors, indicating that too few resistance mechanisms have received clinical validation. In contrast, we verified the enrichment of multiple investigational and hypothetical resistance mechanisms in treated compared with nontreated patients, thereby confirming their putative role in treatment resistance. Additionally, we demonstrated that molecular markers improve 6-month survival prediction, particularly in patients with advanced breast cancer. Our analysis establishes the utility of the META-PRISM cohort for investigating resistance mechanisms and performing predictive analyses in cancer. SIGNIFICANCE: This study highlights the paucity of standard-of-care markers that explain treatment resistance and the promise of investigational and hypothetical markers awaiting further validation. It also demonstrates the utility of molecular profiling in advanced-stage cancers, particularly breast cancer, to improve the survival prediction and assess eligibility to phase I clinical trials. This article is highlighted in the In This Issue feature, p. 1027.

Frequency and Genomic Aspects of Intrinsic Resistance to Vismodegib in Locally Advanced Basal Cell Carcinoma.

PURPOSE: Vismodegib is approved for the treatment of locally advanced basal cell carcinoma (laBCC), but some cases demonstrate intrinsic resistance (IR) to the drug. We sought to assess the frequency of IR to vismodegib in laBCC and its underlying genomic mechanisms. EXPERIMENTAL DESIGN: Response to vismodegib was evaluated in a cohort of 148 laBCC patients. Comprehensive genomic and transcriptomic profiling was performed in a subset of five intrinsically resistant BCC (IR-BCC). RESULTS: We identified that IR-BCC represents 6.1% of laBCC in the studied cohort. Prior treatment with chemotherapy was associated with IR. Genetic events that were previously associated with acquired resistance (AR) in BCC or medulloblastoma were observed in three out of five IR-BCC. However, IR-BCCs were distinct by highly rearranged polyploid genomes. Functional analyses identified hyperactivation of the HIPPO-YAP and WNT pathways at RNA and protein levels in IR-BCC. In vitro assay on the BCC cell line further confirmed that YAP1 overexpression increases the cell proliferation rate. CONCLUSIONS: IR to vismodegib is a rare event in laBCC. IR-BCCs frequently harbor resistance mutations in the Hh pathway, but also are characterized by hyperactivation of the HIPPO-YAP and WNT pathways.

RAS mutations drive proliferative chronic myelomonocytic leukemia via a KMT2A-PLK1 axis.

Proliferative chronic myelomonocytic leukemia (pCMML), an aggressive CMML subtype, is associated with dismal outcomes. RAS pathway mutations, mainly NRASG12D, define the pCMML phenotype as demonstrated by our exome sequencing, progenitor colony assays and a Vav-Cre-NrasG12D mouse model. Further, these mutations promote CMML transformation to acute myeloid leukemia. Using a multiomics platform and biochemical and molecular studies we show that in pCMML RAS pathway mutations are associated with a unique gene expression profile enriched in mitotic kinases such as polo-like kinase 1 (PLK1). PLK1 transcript levels are shown to be regulated by an unmutated lysine methyl-transferase (KMT2A) resulting in increased promoter monomethylation of lysine 4 of histone 3. Pharmacologic inhibition of PLK1 in RAS mutant patient-derived xenografts, demonstrates the utility of personalized biomarker-driven therapeutics in pCMML.

XPC deficiency increases risk of hematologic malignancies through mutator phenotype and characteristic mutational signature.

Recent studies demonstrated a dramatically increased risk of leukemia in patients with a rare genetic disorder, Xeroderma Pigmentosum group C (XP-C), characterized by constitutive deficiency of global genome nucleotide excision repair (GG-NER). The genetic mechanisms of non-skin cancers in XP-C patients remain unexplored. In this study, we analyze a unique collection of internal XP-C tumor genomes including 6 leukemias and 2 sarcomas. We observe a specific mutational pattern and an average of 25-fold increase of mutation rates in XP-C versus sporadic leukemia which we presume leads to its elevated incidence and early appearance. We describe a strong mutational asymmetry with respect to transcription and the direction of replication in XP-C tumors suggesting association of mutagenesis with bulky purine DNA lesions of probably endogenous origin. These findings suggest existence of a balance between formation and repair of bulky DNA lesions by GG-NER in human body cells which is disrupted in XP-C patients.

Topoisomerase 1 prevents replication stress at R-loop-enriched transcription termination sites.

R-loops have both positive and negative impacts on chromosome functions. To identify toxic R-loops in the human genome, here, we map RNA:DNA hybrids, replication stress markers and DNA double-strand breaks (DSBs) in cells depleted for Topoisomerase I (Top1), an enzyme that relaxes DNA supercoiling and prevents R-loop formation. RNA:DNA hybrids are found at both promoters (TSS) and terminators (TTS) of highly expressed genes. In contrast, the phosphorylation of RPA by ATR is only detected at TTS, which are preferentially replicated in a head-on orientation relative to the direction of transcription. In Top1-depleted cells, DSBs also accumulate at TTS, leading to persistent checkpoint activation, spreading of γ-H2AX on chromatin and global replication fork slowdown. These data indicate that fork pausing at the TTS of highly expressed genes containing R-loops prevents head-on conflicts between replication and transcription and maintains genome integrity in a Top1-dependent manner.

Mec1 Is Activated at the Onset of Normal S Phase by Low-dNTP Pools Impeding DNA Replication.

The Mec1 and Rad53 kinases play a central role during acute replication stress in budding yeast. They are also essential for viability in normal growth conditions, but the signal that activates the Mec1-Rad53 pathway in the absence of exogenous insults is currently unknown. Here, we show that this pathway is active at the onset of normal S phase because deoxyribonucleotide triphosphate (dNTP) levels present in G1 phase may not be sufficient to support processive DNA synthesis and impede DNA replication. This activation can be suppressed experimentally by increasing dNTP levels in G1 phase. Moreover, we show that unchallenged cells entering S phase in the absence of Rad53 undergo irreversible fork collapse and mitotic catastrophe. Together, these data indicate that cells use suboptimal dNTP pools to detect the onset of DNA replication and activate the Mec1-Rad53 pathway, which in turn maintains functional forks and triggers dNTP synthesis, allowing the completion of DNA replication.

MRX Increases Chromatin Accessibility at Stalled Replication Forks to Promote Nascent DNA Resection and Cohesin Loading.

The recovery of stalled replication forks depends on the controlled resection of nascent DNA and on the loading of cohesin. These processes operate in the context of nascent chromatin, but the impact of nucleosome structure on a fork restart remains poorly understood. Here, we show that the Mre11-Rad50-Xrs2 (MRX) complex acts together with the chromatin modifiers Gcn5 and Set1 and the histone remodelers RSC, Chd1, and Isw1 to promote chromatin remodeling at stalled forks. Increased chromatin accessibility facilitates the resection of nascent DNA by the Exo1 nuclease and the Sgs1 and Chl1 DNA helicases. Importantly, increased ssDNA promotes the recruitment of cohesin to arrested forks in a Scc2-Scc4-dependent manner. Altogether, these results indicate that MRX cooperates with chromatin modifiers to orchestrate the action of remodelers, nucleases, and DNA helicases, promoting the resection of nascent DNA and the loading of cohesin, two key processes involved in the recovery of arrested forks.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

Biased allelic expression in human primary fibroblast single cells.

The study of gene expression in mammalian single cells via genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblasts over 133,633 unique heterozygous single-nucleotide variants (hetSNVs). We observed that at the snapshot of analyses, each cell contained mostly transcripts from one allele from the majority of genes; indeed, 76.4% of the hetSNVs displayed stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibited a haplotype-consistent allelic ratio; in contrast, distant sites located in two different genes were independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in the majority of the single cells at a given time point were rare and enriched with highly expressed genes. The relative abundance of each allele in a cell was controlled by some regulatory mechanisms given that we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing.

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans' lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

Putative cis-regulatory drivers in colorectal cancer.

The cis-regulatory effects responsible for cancer development have not been as extensively studied as the perturbations of the protein coding genome in tumorigenesis. To better characterize colorectal cancer (CRC) development we conducted an RNA-sequencing experiment of 103 matched tumour and normal colon mucosa samples from Danish CRC patients, 90 of which were germline-genotyped. By investigating allele-specific expression (ASE) we show that the germline genotypes remain important determinants of allelic gene expression in tumours. Using the changes in ASE in matched pairs of samples we discover 71 genes with excess of somatic cis-regulatory effects in CRC, suggesting a cancer driver role. We correlate genotypes and gene expression to identify expression quantitative trait loci (eQTLs) and find 1,693 and 948 eQTLs in normal samples and tumours, respectively. We estimate that 36% of the tumour eQTLs are exclusive to CRC and show that this specificity is partially driven by increased expression of specific transcription factors and changes in methylation patterns. We show that tumour-specific eQTLs are more enriched for low CRC genome-wide association study (GWAS) P values than shared eQTLs, which suggests that some of the GWAS variants are tumour specific regulatory variants. Importantly, tumour-specific eQTL genes also accumulate more somatic mutations when compared to the shared eQTL genes, raising the possibility that they constitute germline-derived cancer regulatory drivers. Collectively the integration of genome and the transcriptome reveals a substantial number of putative somatic and germline cis-regulatory cancer changes that may have a role in tumorigenesis.

Comprehensive human virus screening using high-throughput sequencing with a user-friendly representation of bioinformatics analysis: a pilot study.

High-throughput sequencing (HTS) provides the means to analyze clinical specimens in unprecedented molecular detail. While this technology has been successfully applied to virus discovery and other related areas of research, HTS methodology has yet to be exploited for use in a clinical setting for routine diagnostics. Here, a bioinformatics pipeline (ezVIR) was designed to process HTS data from any of the standard platforms and to evaluate the entire spectrum of known human viruses at once, providing results that are easy to interpret and customizable. The pipeline works by identifying the most likely viruses present in the specimen given the sequencing data. Additionally, ezVIR can generate optional reports for strain typing, can create genome coverage histograms, and can perform cross-contamination analysis for specimens prepared in series. In this pilot study, the pipeline was challenged using HTS data from 20 clinical specimens representative of those most often collected and analyzed in daily practice. The specimens (5 cerebrospinal fluid, 7 bronchoalveolar lavage fluid, 5 plasma, 2 serum, and 1 nasopharyngeal aspirate) were originally found to be positive for a diverse range of DNA or RNA viruses by routine molecular diagnostics. The ezVIR pipeline correctly identified 14 of 14 specimens containing viruses with genomes of <40,000 bp, and 4 of 6 specimens positive for large-genome viruses. Although further validation is needed to evaluate sensitivity and to define detection cutoffs, results obtained in this pilot study indicate that the overall detection success rate, coupled with the ease of interpreting the analysis reports, makes it worth considering using HTS for clinical diagnostics.

The histone deacetylases sir2 and rpd3 act on ribosomal DNA to control the replication program in budding yeast.

In S. cerevisiae, replication timing is controlled by epigenetic mechanisms restricting the accessibility of origins to limiting initiation factors. About 30% of these origins are located within repetitive DNA sequences such as the ribosomal DNA (rDNA) array, but their regulation is poorly understood. Here, we have investigated how histone deacetylases (HDACs) control the replication program in budding yeast. This analysis revealed that two HDACs, Rpd3 and Sir2, control replication timing in an opposite manner. Whereas Rpd3 delays initiation at late origins, Sir2 is required for the timely activation of early origins. Moreover, Sir2 represses initiation at rDNA origins, whereas Rpd3 counteracts this effect. Remarkably, deletion of SIR2 restored normal replication in rpd3Δ cells by reactivating rDNA origins. Together, these data indicate that HDACs control the replication timing program in budding yeast by modulating the ability of repeated origins to compete with single-copy origins for limiting initiation factors.

Comprehensive metagenomic analysis of glioblastoma reveals absence of known virus despite antiviral-like type I interferon gene response.

Glioblastoma is a deadly malignant brain tumor and one of the most incurable forms of cancer in need of new therapeutic targets. As some cancers are known to be caused by a virus, the discovery of viruses could open the possibility to treat, and perhaps prevent, such a disease. Although an association with viruses such as cytomegalovirus or Simian virus 40 has been strongly suggested, involvement of these and other viruses in the initiation and/or propagation of glioblastoma remains vague, controversial and warrants elucidation. To exhaustively address the association of virus and glioblastoma, we developed and validated a robust metagenomic approach to analyze patient biopsies via high-throughput sequencing, a sensitive tool for virus screening. In addition to traditional clinical diagnostics, glioblastoma biopsies were deep-sequenced and analyzed with a multistage computational pipeline to identify known or potentially discover unknown viruses. In contrast to the studies reporting the presence of viral signatures in glioblastoma, no common or recurring active viruses were detected, despite finding an antiviral-like type I interferon response in some specimens. Our findings highlight a discrete and non-specific viral signature and uncharacterized short RNA sequences in glioblastoma. This study provides new insights into glioblastoma pathogenesis and defines a general methodology that can be used for high-resolution virus screening and discovery in human cancers.

Coordinated effects of sequence variation on DNA binding, chromatin structure, and transcription.

DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Passive and active DNA methylation and the interplay with genetic variation in gene regulation.

DNA methylation is an essential epigenetic mark whose role in gene regulation and its dependency on genomic sequence and environment are not fully understood. In this study we provide novel insights into the mechanistic relationships between genetic variation, DNA methylation and transcriptome sequencing data in three different cell-types of the GenCord human population cohort. We find that the association between DNA methylation and gene expression variation among individuals are likely due to different mechanisms from those establishing methylation-expression patterns during differentiation. Furthermore, cell-type differential DNA methylation may delineate a platform in which local inter-individual changes may respond to or act in gene regulation. We show that unlike genetic regulatory variation, DNA methylation alone does not significantly drive allele specific expression. Finally, inferred mechanistic relationships using genetic variation as well as correlations with TF abundance reveal both a passive and active role of DNA methylation to regulatory interactions influencing gene expression. DOI:http://dx.doi.org/10.7554/eLife.00523.001.

easyRNASeq: a bioconductor package for processing RNA-Seq data.

MOTIVATION: RNA sequencing is becoming a standard for expression profiling experiments and many tools have been developed in the past few years to analyze RNA-Seq data. Numerous 'Bioconductor' packages are available for next-generation sequencing data loading in R, e.g. ShortRead and Rsamtools as well as to perform differential gene expression analyses, e.g. DESeq and edgeR. However, the processing tasks lying in between these require the precise interplay of many Bioconductor packages, e.g. Biostrings, IRanges or external solutions are to be sought. RESULTS: We developed 'easyRNASeq', an R package that simplifies the processing of RNA sequencing data, hiding the complex interplay of the required packages behind a single functionality. AVAILABILITY: The package is implemented in R (as of version 2.15) and is available from Bioconductor (as of version 2.10) at the URL: http://bioconductor.org/packages/release/bioc/html/easyRNASeq.html, where installation and usage instructions can be found. CONTACT: delhomme@embl.de.

Celiac disease T-cell epitopes from gamma-gliadins: immunoreactivity depends on the genome of origin, transcript frequency, and flanking protein variation.

BACKGROUND: Celiac disease (CD) is caused by an uncontrolled immune response to gluten, a heterogeneous mixture of wheat storage proteins. The CD-toxicity of these proteins and their derived peptides is depending on the presence of specific T-cell epitopes (9-mer peptides; CD epitopes) that mediate the stimulation of HLA-DQ2/8 restricted T-cells. Next to the thoroughly characterized major T-cell epitopes derived from the α-gliadin fraction of gluten, γ-gliadin peptides are also known to stimulate T-cells of celiac disease patients. To pinpoint CD-toxic γ-gliadins in hexaploid bread wheat, we examined the variation of T-cell epitopes involved in CD in γ-gliadin transcripts of developing bread wheat grains. RESULTS: A detailed analysis of the genetic variation present in γ-gliadin transcripts of bread wheat (T. aestivum, allo-hexaploid, carrying the A, B and D genome), together with genomic γ-gliadin sequences from ancestrally related diploid wheat species, enabled the assignment of sequence variants to one of the three genomic γ-gliadin loci, Gli-A1, Gli-B1 or Gli-D1. Almost half of the γ-gliadin transcripts of bread wheat (49%) was assigned to locus Gli-D1. Transcripts from each locus differed in CD epitope content and composition. The Gli-D1 transcripts contained the highest frequency of canonical CD epitope cores (on average 10.1 per transcript) followed by the Gli-A1 transcripts (8.6) and the Gli-B1 transcripts (5.4). The natural variants of the major CD epitope from γ-gliadins, DQ2-γ-I, showed variation in their capacity to induce in vitro proliferation of a DQ2-γ-I specific and HLA-DQ2 restricted T-cell clone. CONCLUSIONS: Evaluating the CD epitopes derived from γ-gliadins in their natural context of flanking protein variation, genome specificity and transcript frequency is a significant step towards accurate quantification of the CD toxicity of bread wheat. This approach can be used to predict relative levels of CD toxicity of individual wheat cultivars directly from their transcripts (cDNAs).