Transcriptome diversity is a systematic source of variation in RNA-sequencing data.

RNA sequencing has been widely used as an essential tool to probe gene expression. While standard practices have been established to analyze RNA-seq data, it is still challenging to interpret and remove artifactual signals. Several biological and technical factors such as sex, age, batches, and sequencing technology have been found to bias these estimates. Probabilistic estimation of expression residuals (PEER), which infers broad variance components in gene expression measurements, has been used to account for some systematic effects, but it has remained challenging to interpret these PEER factors. Here we show that transcriptome diversity-a simple metric based on Shannon entropy-explains a large portion of variability in gene expression and is the strongest known factor encoded in PEER factors. We then show that transcriptome diversity has significant associations with multiple technical and biological variables across diverse organisms and datasets. In sum, transcriptome diversity provides a simple explanation for a major source of variation in both gene expression estimates and PEER covariates.

Carcinogen susceptibility is regulated by genome architecture and predicts cancer mutagenesis.

The development of many sporadic cancers is directly initiated by carcinogen exposure. Carcinogens induce malignancies by creating DNA lesions (i.e., adducts) that can result in mutations if left unrepaired. Despite this knowledge, there has been remarkably little investigation into the regulation of susceptibility to acquire DNA lesions. In this study, we present the first quantitative human genome-wide map of DNA lesions induced by ultraviolet (UV) radiation, the ubiquitous carcinogen in sunlight that causes skin cancer. Remarkably, the pattern of carcinogen susceptibility across the genome of primary cells significantly reflects mutation frequency in malignant melanoma. Surprisingly, DNase-accessible euchromatin is protected from UV, while lamina-associated heterochromatin at the nuclear periphery is vulnerable. Many cancer driver genes have an intrinsic increase in carcinogen susceptibility, including the BRAF oncogene that has the highest mutation frequency in melanoma. These findings provide a genome-wide snapshot of DNA injuries at the earliest stage of carcinogenesis. Furthermore, they identify carcinogen susceptibility as an origin of genome instability that is regulated by nuclear architecture and mirrors mutagenesis in cancer.

The INO80 chromatin remodeler sustains metabolic stability by promoting TOR signaling and regulating histone acetylation.

Chromatin remodeling complexes are essential for gene expression programs that coordinate cell function with metabolic status. However, how these remodelers are integrated in metabolic stability pathways is not well known. Here, we report an expansive genetic screen with chromatin remodelers and metabolic regulators in Saccharomyces cerevisiae. We found that, unlike the SWR1 remodeler, the INO80 chromatin remodeling complex is composed of multiple distinct functional subunit modules. We identified a strikingly divergent genetic signature for the Ies6 subunit module that links the INO80 complex to metabolic homeostasis. In particular, mitochondrial maintenance is disrupted in ies6 mutants. INO80 is also needed to communicate TORC1-mediated signaling to chromatin, as ino80 mutants exhibit defective transcriptional profiles and altered histone acetylation of TORC1-responsive genes. Furthermore, comparative analysis reveals subunits of INO80 and mTORC1 have high co-occurrence of alterations in human cancers. Collectively, these results demonstrate that the INO80 complex is a central component of metabolic homeostasis that influences histone acetylation and may contribute to disease when disrupted.

Genome-wide profiles of UV lesion susceptibility, repair, and mutagenic potential in melanoma.

Exposure to the ultraviolet (UV) radiation in sunlight creates DNA lesions, which if left unrepaired can induce mutations and contribute to skin cancer. The two most common UV-induced DNA lesions are the cis-syn cyclobutane pyrimidine dimers (CPDs) and pyrimidine (6-4) pyrimidone photoproducts (6-4PPs), both of which can initiate mutations. Interestingly, mutation frequency across the genomes of many cancers is heterogenous with significant increases in heterochromatin. Corresponding increases in UV lesion susceptibility and decreases in repair are observed in heterochromatin versus euchromatin. However, the individual contributions of CPDs and 6-4PPs to mutagenesis have not been systematically examined in specific genomic and epigenomic contexts. In this study, we compared genome-wide maps of 6-4PP and CPD lesion abundances in primary cells and conducted comprehensive analyses to determine the genetic and epigenetic features associated with susceptibility. Overall, we found a high degree of similarity between 6-4PP and CPD formation, with an enrichment of both in heterochromatin regions. However, when examining the relative levels of the two UV lesions, we found that bivalent and Polycomb-repressed chromatin states were uniquely more susceptible to 6-4PPs. Interestingly, when comparing UV susceptibility and repair with melanoma mutation frequency in these regions, disparate patterns were observed in that susceptibility was not always inversely associated with repair and mutation frequency. Functional enrichment analysis hint at mechanisms of negative selection for these regions that are essential for cell viability, immune function and induce cell death when mutated. Ultimately, these results reveal both the similarities and differences between UV-induced lesions that contribute to melanoma.

The somatic mutation landscape of the human body.

BACKGROUND: Somatic mutations in healthy tissues contribute to aging, neurodegeneration, and cancer initiation, yet they remain largely uncharacterized. RESULTS: To gain a better understanding of the genome-wide distribution and functional impact of somatic mutations, we leverage the genomic information contained in the transcriptome to uniformly call somatic mutations from over 7500 tissue samples, representing 36 distinct tissues. This catalog, containing over 280,000 mutations, reveals a wide diversity of tissue-specific mutation profiles associated with gene expression levels and chromatin states. For example, lung samples with low expression of the mismatch-repair gene MLH1 show a mutation signature of deficient mismatch repair. In addition, we find pervasive negative selection acting on missense and nonsense mutations, except for mutations previously observed in cancer samples, which are under positive selection and are highly enriched in many healthy tissues. CONCLUSIONS: These findings reveal fundamental patterns of tissue-specific somatic evolution and shed light on aging and the earliest stages of tumorigenesis.

Transcriptome profiling of Set5 and Set1 methyltransferases: Tools for visualization of gene expression.

Cells regulate transcription by coordinating the activities of multiple histone modifying complexes. We recently identified the yeast histone H4 methyltransferase Set5 and discovered functional overlap with the histone H3 methyltransferase Set1 in gene expression. Specifically, using next-generation RNA sequencing (RNA-Seq), we found that Set5 and Set1 function synergistically to regulate specific transcriptional programs at subtelomeres and transposable elements [1]. Here we provide a comprehensive description of the methodology and analysis tools corresponding to the data deposited in NCBI's Gene Expression Omnibus (GEO) under the accession number GSE52086. This data complements the experimental methods described in Mas Martín G et al., 2014, and provides the means to explore the cooperative functions of histone H3 and H4 methyltransferases in the regulation of transcription. Furthermore, a fully annotated R code is included to enable researchers to use the following computational tools: comparison of significant differential expression (SDE) profiles; gene ontology enrichment of SDE; and enrichment of SDE relative to chromosomal features, such as centromeres, telomeres, and transposable elements. Overall, we present a bioinformatics platform that can be generally implemented for similar analyses with different datasets and in different organisms.

Set5 and Set1 cooperate to repress gene expression at telomeres and retrotransposons.

A complex interplay between multiple chromatin modifiers is critical for cells to regulate chromatin structure and accessibility during essential DNA-templated processes such as transcription. However, the coordinated activities of these chromatin modifiers in the regulation of gene expression are not fully understood. We previously determined that the budding yeast histone H4 methyltransferase Set5 functions together with Set1, the H3K4 methyltransferase, in specific cellular contexts. Here, we sought to understand the relationship between these evolutionarily conserved enzymes in the regulation of gene expression. We generated a comprehensive genetic interaction map of the functionally uncharacterized Set5 methyltransferase and expanded the existing genetic interactome of the global chromatin modifier Set1, revealing functional overlap of the two enzymes in chromatin-related networks, such as transcription. Furthermore, gene expression profiling via RNA-Seq revealed an unexpected synergistic role of Set1 and Set5 in repressing transcription of Ty transposable elements and genes located in subtelomeric regions. This study uncovers novel pathways in which the methyltransferase Set5 participates and, more importantly, reveals a partnership between Set1 and Set5 in transcriptional repression near repetitive DNA elements in budding yeast. Together, our results define a new functional relationship between histone H3 and H4 methyltransferases, whose combined activity may be implicated in preserving genomic integrity.