Although in vitro studies of embryonic stem cells have identified polycomb repressor complexes (PRCs) as key regulators of differentiation, it remains unclear as to how PRC-mediated mechanisms ...control fates of multipotent progenitors in developing tissues. Here, we show that an essential PRC component, Ezh2, is expressed in epidermal progenitors but diminishes concomitant with embryonic differentiation and with postnatal decline in proliferative activity. We show that Ezh2 controls proliferative potential of basal progenitors by repressing the
Ink4A-Ink4B locus and tempers the developmental rate of differentiation by preventing premature recruitment of AP1 transcriptional activator to the structural genes that are required for epidermal differentiation. Together, our studies reveal that PRCs control epigenetic modifications temporally and spatially in tissue-restricted stem cells. They maintain their proliferative potential and globally repressing undesirable differentiation programs while selectively establishing a specific terminal differentiation program in a stepwise fashion.
Genomic variant detection from next-generation sequencing has become established as an extremely important component of research and clinical diagnoses in both cancer and Mendelian disorders. ...Insertions and deletions (indels) are a common source of variation and can frequently impact functionality, thus making their detection vitally important. While substantial effort has gone into detecting indels from DNA, there is still opportunity for improvement. Further, detection of indels from RNA-Seq data has largely been an afterthought and offers another critical area for variant detection.
We present here ABRA2, a redesign of the original ABRA implementation that offers support for realignment of both RNA and DNA short reads. The process results in improved accuracy and scalability including support for human whole genomes. Results demonstrate substantial improvement in indel detection for a variety of data types, including those that were not previously supported by ABRA. Further, ABRA2 results in broad improvements to variant calling accuracy across a wide range of post-processing workflows including whole genomes, targeted exomes and transcriptome sequencing.
ABRA2 is implemented in a combination of Java and C/C++ and is freely available to all from: https://github.com/mozack/abra2.
Supplementary data are available at Bioinformatics online.
Polymerase theta (Pol θ, gene name Polq) is a widely conserved DNA polymerase that mediates a microhomology-mediated, error-prone, double strand break (DSB) repair pathway, referred to as Theta ...Mediated End Joining (TMEJ). Cells with homologous recombination deficiency are reliant on TMEJ for DSB repair. It is unknown whether deficiencies in other components of the DNA damage response (DDR) also result in Pol θ addiction. Here we use a CRISPR genetic screen to uncover 140 Polq synthetic lethal (PolqSL) genes, the majority of which were previously unknown. Functional analyses indicate that Pol θ/TMEJ addiction is associated with increased levels of replication-associated DSBs, regardless of the initial source of damage. We further demonstrate that approximately 30% of TCGA breast cancers have genetic alterations in PolqSL genes and exhibit genomic scars of Pol θ/TMEJ hyperactivity, thereby substantially expanding the subset of human cancers for which Pol θ inhibition represents a promising therapeutic strategy.
B cell lymphoma and melanoma harbor recurrent mutations in the gene encoding the EZH2 histone methyltransferase (EZH2), but the carcinogenic role of these mutations is unclear. Here we describe a ...mouse model in which the most common somatic Ezh2 gain-of-function mutation (EZH2(Y646F) in human; Ezh2(Y641F) in mouse) is conditionally expressed. Expression of Ezh2(Y641F) in mouse B cells or melanocytes caused high-penetrance lymphoma or melanoma, respectively. Overexpression of the anti-apoptotic protein Bcl2, but not the oncoprotein Myc, or loss of the tumor suppressor protein p53 (encoded by Trp53 in mice) further accelerated lymphoma progression. Expression of the mutant Braf but not the mutant Nras oncoprotein further accelerated melanoma progression. Although expression of Ezh2(Y641F) globally increased the abundance of trimethylated Lys27 of histone H3 (H3K27me3), it also caused a widespread redistribution of this repressive mark, including a loss of H3K27me3 that was associated with increased transcription at many loci. These results suggest that Ezh2(Y641F) induces lymphoma and melanoma through a vast reorganization of chromatin structure, inducing both repression and activation of polycomb-regulated loci.
The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (
I)-exposed children. We analyzed ...genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood
I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.
BACKGROUNDThe KRAS proto-oncogene is among the most frequently mutated genes in cancer, yet for 40 years it remained an elusive therapeutic target. Recently, allosteric inhibitors that covalently ...bind to KRAS G12C mutations have been approved for use in lung adenocarcinomas. Although responses are observed, they are often short-lived, thus making in-depth characterization of the mechanisms of resistance of paramount importance.METHODSHere, we present a rapid-autopsy case of a patient who had a KRASG12C-mutant lung adenocarcinoma who initially responded to a KRAS G12C inhibitor but then rapidly developed resistance. Using deep-RNA and whole-exome sequencing comparing pretreatment, posttreatment, and matched normal tissues, we uncover numerous mechanisms of resistance to direct KRAS inhibition.RESULTSIn addition to decreased KRAS G12C-mutant allele frequency in refractory tumors, we also found reactivation of the MAPK pathway despite no new mutations in KRAS or its downstream mediators. Tumor cell-intrinsic and non-cell autonomous mechanisms included increased complement activation, coagulation, and tumor angiogenesis, and several lines of evidence of immunologic evasion.CONCLUSIONTogether, our findings reveal numerous mechanisms of resistance to current KRAS G12C inhibitors through enrichment of clonal populations, KRAS-independent downstream signaling, and diverse remodeling of the tumor microenvironment.FUNDINGRichard and Fran Duley, Jimmy and Kay Mann, the NIH, and the North Carolina Biotechnology Center.
This study identifies mechanisms mediating responses to immune checkpoint inhibitors using mouse models of triple-negative breast cancer. By creating new mammary tumor models, we find that tumor ...mutation burden and specific immune cells are associated with response. Further, we developed a rich resource of single-cell RNA-seq and bulk mRNA-seq data of immunotherapy-treated and non-treated tumors from sensitive and resistant murine models. Using this, we uncover that immune checkpoint therapy induces T follicular helper cell activation of B cells to facilitate the anti-tumor response in these models. We also show that B cell activation of T cells and the generation of antibody are key to immunotherapy response and propose a new biomarker for immune checkpoint therapy. In total, this work presents resources of new preclinical models of breast cancer with large mRNA-seq and single-cell RNA-seq datasets annotated for sensitivity to therapy and uncovers new components of response to immune checkpoint inhibitors.
Display omitted
•New TNBC murine models with high mutation burden and immune cell activity•A genomics resource of immune checkpoint treated tumors from TNBC murine models•Immune checkpoint blockade activates Tfh and B cells in the anti-tumor response•B cells impact immunotherapy response by secreting antibody and activating T cells
Mouse models of triple-negative breast cancer provide insights into how T follicular helper cell activation of B cells facilitates the effects of immune checkpoint inhibitors.
RNA sequencing (RNA-Seq) is often used for transcriptome profiling as well as the identification of novel transcripts and alternative splicing events. Typically, RNA-Seq libraries are prepared from ...total RNA using poly(A) enrichment of the mRNA (mRNA-Seq) to remove ribosomal RNA (rRNA), however, this method fails to capture non-poly(A) transcripts or partially degraded mRNAs. Hence, a mRNA-Seq protocol will not be compatible for use with RNAs coming from Formalin-Fixed and Paraffin-Embedded (FFPE) samples.
To address the desire to perform RNA-Seq on FFPE materials, we evaluated two different library preparation protocols that could be compatible for use with small RNA fragments. We obtained paired Fresh Frozen (FF) and FFPE RNAs from multiple tumors and subjected these to different gene expression profiling methods. We tested 11 human breast tumor samples using: (a) FF RNAs by microarray, mRNA-Seq, Ribo-Zero-Seq and DSN-Seq (Duplex-Specific Nuclease) and (b) FFPE RNAs by Ribo-Zero-Seq and DSN-Seq. We also performed these different RNA-Seq protocols using 10 TCGA tumors as a validation set.The data from paired RNA samples showed high concordance in transcript quantification across all protocols and between FF and FFPE RNAs. In both FF and FFPE, Ribo-Zero-Seq removed rRNA with comparable efficiency as mRNA-Seq, and it provided an equivalent or less biased coverage on gene 3' ends. Compared to mRNA-Seq where 69% of bases were mapped to the transcriptome, DSN-Seq and Ribo-Zero-Seq contained significantly fewer reads mapping to the transcriptome (20-30%); in these RNA-Seq protocols, many if not most reads mapped to intronic regions. Approximately 14 million reads in mRNA-Seq and 45-65 million reads in Ribo-Zero-Seq or DSN-Seq were required to achieve the same gene detection levels as a standard Agilent DNA microarray.
Our results demonstrate that compared to mRNA-Seq and microarrays, Ribo-Zero-Seq provides equivalent rRNA removal efficiency, coverage uniformity, genome-based mapped reads, and consistently high quality quantification of transcripts. Moreover, Ribo-Zero-Seq and DSN-Seq have consistent transcript quantification using FFPE RNAs, suggesting that RNA-Seq can be used with FFPE-derived RNAs for gene expression profiling.
Topoisomerases are expressed throughout the developing and adult brain and are mutated in some individuals with autism spectrum disorder (ASD). However, how topoisomerases are mechanistically ...connected to ASD is unknown. Here we find that topotecan, a topoisomerase 1 (TOP1) inhibitor, dose-dependently reduces the expression of extremely long genes in mouse and human neurons, including nearly all genes that are longer than 200 kilobases. Expression of long genes is also reduced after knockdown of Top1 or Top2b in neurons, highlighting that both enzymes are required for full expression of long genes. By mapping RNA polymerase II density genome-wide in neurons, we found that this length-dependent effect on gene expression was due to impaired transcription elongation. Interestingly, many high-confidence ASD candidate genes are exceptionally long and were reduced in expression after TOP1 inhibition. Our findings suggest that chemicals and genetic mutations that impair topoisomerases could commonly contribute to ASD and other neurodevelopmental disorders.
Multiple breast cancer gene expression profiles have been developed that appear to provide similar abilities to predict outcome and may outperform clinical-pathologic criteria; however, the extent to ...which seemingly disparate profiles provide additive prognostic information is not known, nor do we know whether prognostic profiles perform equally across clinically defined breast cancer subtypes. We evaluated whether combining the prognostic powers of standard breast cancer clinical variables with a large set of gene expression signatures could improve on our ability to predict patient outcomes.
Using clinical-pathological variables and a collection of 323 gene expression "modules", including 115 previously published signatures, we build multivariate Cox proportional hazards models using a dataset of 550 node-negative systemically untreated breast cancer patients. Models predictive of pathological complete response (pCR) to neoadjuvant chemotherapy were also built using this approach.
We identified statistically significant prognostic models for relapse-free survival (RFS) at 7 years for the entire population, and for the subgroups of patients with ER-positive, or Luminal tumors. Furthermore, we found that combined models that included both clinical and genomic parameters improved prognostication compared with models with either clinical or genomic variables alone. Finally, we were able to build statistically significant combined models for pathological complete response (pCR) predictions for the entire population.
Integration of gene expression signatures and clinical-pathological factors is an improved method over either variable type alone. Highly prognostic models could be created when using all patients, and for the subset of patients with lymph node-negative and ER-positive breast cancers. Other variables beyond gene expression and clinical-pathological variables, like gene mutation status or DNA copy number changes, will be needed to build robust prognostic models for ER-negative breast cancer patients. This combined clinical and genomics model approach can also be used to build predictors of therapy responsiveness, and could ultimately be applied to other tumor types.