Although mechanisms of acquired resistance of epidermal growth factor receptor (EGFR)-mutant non-small-cell lung cancers to EGFR inhibitors have been identified, little is known about how resistant clones evolve during drug therapy. Here we observe that acquired resistance caused by the EGFR(T790M) gatekeeper mutation can occur either by selection of pre-existing EGFR(T790M)-positive clones or via genetic evolution of initially EGFR(T790M)-negative drug-tolerant cells. The path to resistance impacts the biology of the resistant clone, as those that evolved from drug-tolerant cells had a diminished apoptotic response to third-generation EGFR inhibitors that target EGFR(T790M); treatment with navitoclax, an inhibitor of the anti-apoptotic factors BCL-xL and BCL-2 restored sensitivity. We corroborated these findings using cultures derived directly from EGFR inhibitor-resistant patient tumors. These findings provide evidence that clinically relevant drug-resistant cancer cells can both pre-exist and evolve from drug-tolerant cells, and they point to therapeutic opportunities to prevent or overcome resistance in the clinic.

Expression of the p53-inducible antiproliferative gene BTG2 is suppressed in many cancers in the absence of inactivating gene mutations, suggesting alternative mechanisms of silencing. Using a shRNA screen targeting 43 histone lysine methyltransferases (KMTs), we show that SETD1A suppresses BTG2 expression through its induction of several BTG2-targeting miRNAs. This indirect but highly specific mechanism, by which a chromatin regulator that mediates transcriptional activating marks can lead to the downregulation of a critical effector gene, is shared with multiple genes in the p53 pathway. Through such miRNA-dependent effects, SETD1A regulates cell cycle progression in vitro and modulates tumorigenesis in mouse xenograft models. Together, these observations help explain the remarkably specific genetic consequences associated with alterations in generic chromatin modulators in cancer.

Proteins that function in the same pathways, protein complexes or the same environmental conditions can show similar patterns of sequence conservation across phylogenetic clades. In species that no longer require a specific protein complex or pathway, these proteins, as a group, tend to be lost or diverge. Analysis of the similarity in patterns of sequence conservation across a large set of eukaryotes can predict functional associations between different proteins, identify new pathway members and reveal the function of previously uncharacterized proteins. We used normalized phylogenetic profiling to predict protein function and identify new pathway members and disease genes. The phylogenetic profiles of tens of thousands conserved proteins in the human, mouse, Caenorhabditis elegans and Drosophila genomes can be queried on our new web server, PhyloGene. PhyloGene provides intuitive and user-friendly platform to query the patterns of conservation across 86 animal, fungal, plant and protist genomes. A protein query can be submitted either by selecting the name from whole-genome protein sets of the intensively studied species or by entering a protein sequence. The graphic output shows the profile of sequence conservation for the query and the most similar phylogenetic profiles for the proteins in the genome of choice. The user can also download this output in numerical form.

Auditory neuropathy (AN) is a recently recognized hearing disorder characterized by intact outer hair cell function, disrupted auditory nerve synchronization and poor speech perception and recognition. Cochlear implants (CIs) are currently the most promising intervention for improving hearing and speech in individuals with AN. Although previous studies have shown optimistic results, there was large variability concerning benefits of CIs among individuals with AN. The data indicate that different criteria are needed to evaluate the benefit of CIs in these children compared to those with sensorineural hearing loss. We hypothesized that a hierarchic assessment would be more appropriate to evaluate the benefits of cochlear implantation in AN individuals.

Acquired chromosomal DNA amplifications are features of many tumors. Although overexpression and stabilization of the histone H3 lysine 9/36 (H3K9/36) tri-demethylase KDM4A generates transient site-specific copy number gains (TSSGs), additional mechanisms directly controlling site-specific DNA copy gains are not well defined. In this study, we uncover a collection of H3K4-modifying chromatin regulators that function with H3K9 and H3K36 regulators to orchestrate TSSGs. Specifically, the H3K4 tri-demethylase KDM5A and specific COMPASS/KMT2 H3K4 methyltransferases modulate different TSSG loci through H3K4 methylation states and KDM4A recruitment. Furthermore, a distinct chromatin modifier network, MLL1-KDM4B-KDM5B, controls copy number regulation at a specific genomic locus in a KDM4A-independent manner. These pathways comprise an epigenetic addressing system for defining site-specific DNA rereplication and amplifications.

G protein αs (GNAS) mediates receptor-stimulated cAMP signalling, which integrates diverse environmental cues with intracellular responses. GNAS is mutationally activated in multiple tumour types, although its oncogenic mechanisms remain elusive. We explored this question in pancreatic tumourigenesis where concurrent GNAS and KRAS mutations characterize pancreatic ductal adenocarcinomas (PDAs) arising from intraductal papillary mucinous neoplasms (IPMNs). By developing genetically engineered mouse models, we show that GnasR201C cooperates with KrasG12D to promote initiation of IPMN, which progress to invasive PDA following Tp53 loss. Mutant Gnas remains critical for tumour maintenance in vivo. This is driven by protein-kinase-A-mediated suppression of salt-inducible kinases (Sik1-3), associated with induction of lipid remodelling and fatty acid oxidation. Comparison of Kras-mutant pancreatic cancer cells with and without Gnas mutations reveals striking differences in the functions of this network. Thus, we uncover Gnas-driven oncogenic mechanisms, identify Siks as potent tumour suppressors, and demonstrate unanticipated metabolic heterogeneity among Kras-mutant pancreatic neoplasms.

The gene encoding the Insulin-like Growth Factor 2 mRNA binding protein 2/IMP2 is amplified and overexpressed in many human cancers, accompanied by a poorer prognosis. Mice lacking IMP2 exhibit a longer lifespan and a reduced tumor burden at old age. Herein we show in a diverse array of human cancer cells that IMP2 overexpression stimulates and IMP2 elimination diminishes proliferation by 50-80%. In addition to its known ability to promote the abundance of Insulin-like Growth Factor 2/IGF2, we find that IMP2 strongly promotes IGF action, by binding and stabilizing the mRNA encoding the DNA binding protein HMGA1, a known oncogene. HMGA1 suppresses the abundance of IGF binding protein 2/IGFBP2 and Grb14, inhibitors of IGF action. IMP2 stabilization of HMGA1 mRNA plus IMP2 stimulated IGF2 production synergistically drive cancer cell proliferation and account for IMP2's tumor promoting action. IMP2's ability to promote proliferation and IGF action requires IMP2 phosphorylation by mTOR.

Chromatin-modifying enzymes play a fundamental role in regulating chromatin structure so that DNA replication is spatially and temporally coordinated. For example, the lysine demethylase 4A/Jumonji domain-containing 2A (KDM4A/JMJD2A) is tightly regulated during the cell cycle. Overexpression of JMJD2A leads to altered replication timing and faster S phase progression. In this study, we demonstrate that degradation of JMJD2A is regulated by the proteasome. JMJD2A turnover is coordinated through the SKP1-Cul1-F-box ubiquitin ligase complex that contains cullin 1 and the F-box and leucine-rich repeat protein 4 (FbxL4). This complex interacted with JMJD2A. Ubiquitin overexpression restored turnover and blocked the JMJD2A-dependent faster S phase progression in a cullin 1-dependent manner. Furthermore, increased ubiquitin levels decreased JMJD2A occupancy and BrdU incorporation at target sites. This study highlights a finely tuned mechanism for regulating histone demethylase levels and emphasizes the need to tightly regulate chromatin modifiers so that the cell cycle occurs properly.

The KDM4/JMJD2 family of histone demethylases is amplified in human cancers. However, little is known about their physiologic or tumorigenic roles. We have identified a conserved and unappreciated role for the JMJD2A/KDM4A H3K9/36 tridemethylase in cell cycle progression. We demonstrate that JMJD2A protein levels are regulated in a cell cycle-dependent manner and that JMJD2A overexpression increased chromatin accessibility, S phase progression, and altered replication timing of specific genomic loci. These phenotypes depended on JMJD2A enzymatic activity. Strikingly, depletion of the only C. elegans homolog, JMJD-2, slowed DNA replication and increased ATR/p53-dependent apoptosis. Importantly, overexpression of HP1γ antagonized JMJD2A-dependent progression through S phase, and depletion of HPL-2 rescued the DNA replication-related phenotypes in jmjd-2(-/-) animals. Our findings describe a highly conserved model whereby JMJD2A regulates DNA replication by antagonizing HP1γ and controlling chromatin accessibility.

DNA double-strand break (DSB) repair involves complex interactions between chromatin and repair proteins, including Tip60, a tumour suppressor. Tip60 is an acetyltransferase that acetylates both histones and ATM (ataxia telangiectasia mutated) kinase. Inactivation of Tip60 leads to defective DNA repair and increased cancer risk. However, how DNA damage activates the acetyltransferase activity of Tip60 is not known. Here, we show that direct interaction between the chromodomain of Tip60 and histone H3 trimethylated on lysine 9 (H3K9me3) at DSBs activates the acetyltransferase activity of Tip60. Depletion of intracellular H3K9me3 blocks activation of the acetyltransferase activity of Tip60, resulting in defective ATM activation and widespread defects in DSB repair. In addition, the ability of Tip60 to access H3K9me3 is dependent on the DNA damage-induced displacement of HP1beta (heterochromatin protein 1beta) from H3K9me3. Finally, we demonstrate that the Mre11-Rad50-Nbs1 (MRN) complex targets Tip60 to H3K9me3, and is required to activate the acetyltransferase activity of Tip60. These results reveal a new function for H3K9me3 in coordinating activation of Tip60-dependent DNA repair pathways, and imply that aberrant patterns of histone methylation may contribute to cancer by altering the efficiency of DSB repair.

To identify novel genes and pathways associated with AMD, we performed microarray gene expression and linkage analysis which implicated the candidate gene, retinoic acid receptor-related orphan receptor alpha (RORA, 15q). Subsequent genotyping of 159 RORA single nucleotide polymorphisms (SNPs) in a family-based cohort, followed by replication in an unrelated case-control cohort, demonstrated that SNPs and haplotypes located in intron 1 were significantly associated with neovascular AMD risk in both cohorts. This is the first report demonstrating a possible role for RORA, a receptor for cholesterol, in the pathophysiology of AMD. Moreover, we found a significant interaction between RORA and the ARMS2/HTRA1 locus suggesting a novel pathway underlying AMD pathophysiology.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

Despite continued widespread use, the genomic effects of cisplatin-based chemotherapy and implications for subsequent treatment are incompletely characterized. Here, we analyze whole exome sequencing of matched pre- and post-neoadjuvant cisplatin-based chemotherapy primary bladder tumor samples from 30 muscle-invasive bladder cancer patients. We observe no overall increase in tumor mutational burden post-chemotherapy, though a significant proportion of subclonal mutations are unique to the matched pre- or post-treatment tumor, suggesting chemotherapy-induced and/or spatial heterogeneity. We subsequently identify and validate a novel mutational signature in post-treatment tumors consistent with known characteristics of cisplatin damage and repair. We find that post-treatment tumor heterogeneity predicts worse overall survival, and further observe alterations in cell-cycle and immune checkpoint regulation genes in post-treatment tumors. These results provide insight into the clinical and genomic dynamics of tumor evolution with cisplatin-based chemotherapy, suggest mechanisms of clinical resistance, and inform development of clinically relevant biomarkers and trials of combination therapies.

A sedentary lifestyle, chronic inflammation and leukocytosis increase atherosclerosis; however, it remains unclear whether regular physical activity influences leukocyte production. Here we show that voluntary running decreases hematopoietic activity in mice. Exercise protects mice and humans with atherosclerosis from chronic leukocytosis but does not compromise emergency hematopoiesis in mice. Mechanistically, exercise diminishes leptin production in adipose tissue, augmenting quiescence-promoting hematopoietic niche factors in leptin-receptor-positive stromal bone marrow cells. Induced deletion of the leptin receptor in Prrx1-creERT2; Leprfl/fl mice reveals that leptin's effect on bone marrow niche cells regulates hematopoietic stem and progenitor cell (HSPC) proliferation and leukocyte production, as well as cardiovascular inflammation and outcomes. Whereas running wheel withdrawal quickly reverses leptin levels, the impact of exercise on leukocyte production and on the HSPC epigenome and transcriptome persists for several weeks. Together, these data show that physical activity alters HSPCs via modulation of their niche, reducing hematopoietic output of inflammatory leukocytes.

Current genomics methods are designed to handle tens to thousands of samples but will need to scale to millions to match the pace of data and hypothesis generation in biomedical science. Here, we show that high efficiency at low cost can be achieved by leveraging general-purpose libraries for computing using graphics processing units (GPUs), such as PyTorch and TensorFlow. We demonstrate > 200-fold decreases in runtime and ~ 5-10-fold reductions in cost relative to CPUs. We anticipate that the accessibility of these libraries will lead to a widespread adoption of GPUs in computational genomics.

Chromatin modulation provides a key checkpoint for controlling cell cycle regulated gene networks. The replicative canonical histone genes are one such gene family under tight regulation during cell division. These genes are most highly expressed during S phase when histones are needed to chromatinize the new DNA template. While this fact has been known for a while, limited knowledge exists about the specific chromatin regulators controlling their temporal expression during cell cycle. Since histones and their associated mutations are emerging as major players in diseases such as cancer, identifying the chromatin factors modulating their expression is critical. The histone lysine tri-demethylase KDM4A is regulated over cell cycle and plays a direct role in DNA replication timing, site-specific rereplication, and DNA amplifications during S phase. Here, we establish an unappreciated role for the catalytically active KDM4A in directly regulating canonical replicative histone gene networks during cell cycle. Of interest, we further demonstrate that KDM4A interacts with proteins controlling histone expression and RNA processing (i.e., hnRNPUL1 and FUS/TLS). Together, this study provides a new function for KDM4A in modulating canonical histone gene expression.

Evolved resistance to tyrosine kinase inhibitor (TKI)-targeted therapies remains a major clinical challenge. In epidermal growth factor receptor (EGFR) mutant non-small-cell lung cancer (NSCLC), failure of EGFR TKIs can result from both genetic and epigenetic mechanisms of acquired drug resistance. Widespread reports of histologic and gene expression changes consistent with an epithelial-to-mesenchymal transition (EMT) have been associated with initially surviving drug-tolerant persister cells, which can seed bona fide genetic mechanisms of resistance to EGFR TKIs. While therapeutic approaches targeting fully resistant cells, such as those harboring an EGFRT790M mutation, have been developed, a clinical strategy for preventing the emergence of persister cells remains elusive. Using mesenchymal cell lines derived from biopsies of patients who progressed on EGFR TKI as surrogates for persister populations, we performed whole-genome CRISPR screening and identified fibroblast growth factor receptor 1 (FGFR1) as the top target promoting survival of mesenchymal EGFR mutant cancers. Although numerous previous reports of FGFR signaling contributing to EGFR TKI resistance in vitro exist, the data have not yet been sufficiently compelling to instigate a clinical trial testing this hypothesis, nor has the role of FGFR in promoting the survival of persister cells been elucidated. In this study, we find that combining EGFR and FGFR inhibitors inhibited the survival and expansion of EGFR mutant drug-tolerant cells over long time periods, preventing the development of fully resistant cancers in multiple vitro models and in vivo. These results suggest that dual EGFR and FGFR blockade may be a promising clinical strategy for both preventing and overcoming EMT-associated acquired drug resistance and provide motivation for the clinical study of combined EGFR and FGFR inhibition in EGFR-mutated NSCLCs.

Totipotency emerges in early embryogenesis, but its molecular underpinnings remain poorly characterized. In the present study, we employed DNA fiber analysis to investigate how pluripotent stem cells are reprogrammed into totipotent-like 2-cell-like cells (2CLCs). We show that totipotent cells of the early mouse embryo have slow DNA replication fork speed and that 2CLCs recapitulate this feature, suggesting that fork speed underlies the transition to a totipotent-like state. 2CLCs emerge concomitant with DNA replication and display changes in replication timing (RT), particularly during the early S-phase. RT changes occur prior to 2CLC emergence, suggesting that RT may predispose to gene expression changes and consequent reprogramming of cell fate. Slowing down replication fork speed experimentally induces 2CLCs. In vivo, slowing fork speed improves the reprogramming efficiency of somatic cell nuclear transfer. Our data suggest that fork speed regulates cellular plasticity and that remodeling of replication features leads to changes in cell fate and reprogramming.

The benefit of adjuvant cyclin-dependent kinase 4 and 6 (CDK4/6) inhibitors with endocrine therapy (ET) in hormone receptor-positive, human epidermal growth factor 2 receptor-negative (HR+/HER2-) early breast cancer (EBC) is uncertain. Hence, we performed a meta-analysis to determine the efficacy and safety of adjuvant CDK4/6 inhibitors plus ET and to identify potential preferred subpopulations for this regimen.

Circular RNAs (circRNAs) are group of noncoding RNAs derived from back-splicing events. Accumulating evidence certifies the critical roles of circRNAs in human tumorigenesis. However, the role and biogenesis of circRNAs in cervical cancer are still unclear. Here, a novel identified circRNA, circSLC26A4, was found to be upregulated in cervical cancer tissue and cells. Clinically, the high expression of circSLC26A4 was related to the poor survival of cervical cancer patients. Functionally, cellular experiments indicated that circSLC26A4 knockdown repressed the proliferation, invasion, and tumor growth in vitro and in vivo. Furthermore, circSLC26A4 acted as the sponge of miR-1287-5p; moreover, miR-1287-5p targeted the 3' UTR of HOXA7 mRNA. Mechanistically, RNA binding protein (RBP) quaking (QKI) was identified to interact with the QKI response elements (QREs) in SLC26A4 gene introns, thereby promoting circSLC26A4 biogenesis. In conclusion, these findings demonstrate that circSLC26A4 facilitates cervical cancer progression through the QKI/circSLC26A4/miR-1287-5p/HOXA7 axis, which might bring novel therapeutic strategies for cervical cancer.