The most striking characteristic of CHO cells is their adaptability, which enables efficient production of proteins as well as growth under a variety of culture conditions, but also results in genomic and phenotypic instability. To investigate the relative contribution of genomic and epigenetic modifications towards phenotype evolution, comprehensive genome and epigenome data are presented for six related CHO cell lines, both in response to perturbations (different culture conditions and media as well as selection of a specific phenotype with increased transient productivity) and in steady state (prolonged time in culture under constant conditions). Clear transitions were observed in DNA-methylation patterns upon each perturbation, while few changes occurred over time under constant conditions. Only minor DNA-methylation changes were observed between exponential and stationary growth phase; however, throughout a batch culture the histone modification pattern underwent continuous adaptation. Variation in genome sequence between the six cell lines on the level of SNPs, InDels, and structural variants is high, both upon perturbation and under constant conditions over time. The here presented comprehensive resource may open the door to improved control and manipulation of gene expression during industrial bioprocesses based on epigenetic mechanisms. Biotechnol. Bioeng. 2016;113: 2241-2253. © 2016 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

The Ensembl gene annotation system has been used to annotate over 70 different vertebrate species across a wide range of genome projects. Furthermore, it generates the automatic alignment-based annotation for the human and mouse GENCODE gene sets. The system is based on the alignment of biological sequences, including cDNAs, proteins and RNA-seq reads, to the target genome in order to construct candidate transcript models. Careful assessment and filtering of these candidate transcripts ultimately leads to the final gene set, which is made available on the Ensembl website. Here, we describe the annotation process in detail.Database URL: http://www.ensembl.org/index.html.

Large scale transcript analysis of human glomerular microvascular endothelial cells (HGMEC) has never been accomplished. We designed this study to define the transcriptome of HGMEC and facilitate a better characterization of these endothelial cells with unique features. Serial analysis of gene expression (SAGE) was used for its unbiased approach to quantitative acquisition of transcripts.

Immune responses are tightly regulated to ensure efficient pathogen clearance while avoiding tissue damage. Here we report that Setdb2 was the only protein lysine methyltransferase induced during infection with influenza virus. Setdb2 expression depended on signaling via type I interferons, and Setdb2 repressed expression of the gene encoding the neutrophil attractant CXCL1 and other genes that are targets of the transcription factor NF-κB. This coincided with occupancy by Setdb2 at the Cxcl1 promoter, which in the absence of Setdb2 displayed diminished trimethylation of histone H3 Lys9 (H3K9me3). Mice with a hypomorphic gene-trap construct of Setdb2 exhibited increased infiltration of neutrophils during sterile lung inflammation and were less sensitive to bacterial superinfection after infection with influenza virus. This suggested that a Setdb2-mediated regulatory crosstalk between the type I interferons and NF-κB pathways represents an important mechanism for virus-induced susceptibility to bacterial superinfection.

We provide a catalog for the effects of the human kinome on cell survival in response to DNA-damaging agents, covering all major DNA repair pathways. By treating 313 kinase-deficient cell lines with ten diverse DNA-damaging agents, including seven commonly used chemotherapeutics, we identified examples of vulnerability and resistance that are kinase specific. To investigate synthetic lethal interactions, we tested the response to carmustine for 25 cell lines by establishing a phenotypic fluorescence-activated cell sorting (FACS) assay designed to validate gene-drug interactions. We show apoptosis, cell cycle changes, and DNA damage and proliferation after alkylation- or crosslink-induced damage. In addition, we reconstitute the cellular sensitivity of DYRK4, EPHB6, MARK3, and PNCK as a proof of principle for our study. Furthermore, using global phosphoproteomics on cells lacking MARK3, we provide evidence for its role in the DNA damage response. Our data suggest that cancers with inactivating mutations in kinases, including MARK3, are particularly vulnerable to alkylating chemotherapeutic agents.

Patients with refractory or relapsed haematological malignancies have few treatment options and short survival times. Identification of effective therapies with genomic-based precision medicine is hampered by intratumour heterogeneity and incomplete understanding of the contribution of various mutations within specific cancer phenotypes. Ex-vivo drug-response profiling in patient biopsies might aid effective treatment identification; however, proof of its clinical utility is limited.

STAT5B is often mutated in hematopoietic malignancies. The most frequent STAT5B mutation, Asp642His (N642H), has been found in over 90 leukemia and lymphoma patients. Here, we used the Vav1 promoter to generate transgenic mouse models that expressed either human STAT5B or STAT5BN642H in the hematopoietic compartment. While STAT5B-expressing mice lacked a hematopoietic phenotype, the STAT5BN642H-expressing mice rapidly developed T cell neoplasms. Neoplasia manifested as transplantable CD8+ lymphoma or leukemia, indicating that the STAT5BN642H mutation drives cancer development. Persistent and enhanced levels of STAT5BN642H tyrosine phosphorylation in transformed CD8+ T cells led to profound changes in gene expression that were accompanied by alterations in DNA methylation at potential histone methyltransferase EZH2-binding sites. Aurora kinase genes were enriched in STAT5BN642H-expressing CD8+ T cells, which were exquisitely sensitive to JAK and Aurora kinase inhibitors. Together, our data suggest that JAK and Aurora kinase inhibitors should be further explored as potential therapeutics for lymphoma and leukemia patients with the STAT5BN642H mutation who respond poorly to conventional chemotherapy.

Maintaining sufficient levels of Pdx1 activity is a prerequisite for proper regulation of blood glucose homeostasis and beta cell function. Mice that are haploinsufficient for Pdx1 display impaired glucose tolerance and lack the ability to increase beta cell mass in response to decreased insulin signaling. Several studies have shown that post-translational modifications are regulating Pdx1 activity through intracellular localization and binding to co-factors. Understanding the signaling cues converging on Pdx1 and modulating its activity is therefore an attractive approach in diabetes treatment. We employed a novel technique called Nanofluidic Proteomic Immunoassay to characterize the post-translational profile of Pdx1. Following isoelectric focusing in nano-capillaries, this technology relies on a pan specific antibody for detection and it therefore allows the relative abundance of differently charged protein species to be examined simultaneously. In all eukaryotic cells tested we find that the Pdx1 protein separates into four distinct peaks whereas Pdx1 protein from bacteria only produces one peak. Of the four peaks in eukaryotic cells we correlate one of them to a phosphorylation Using alanine scanning and mass spectrometry we map this phosphorylation to serine 61 in both Min6 cells and in exogenous Pdx1 over-expressed in HEK293 cells. A single phosphorylation is also present in cultured islets but it remains unaffected by changes in glucose levels. It is present during embryogenesis but is not required for pancreas development.

In chronic myeloid leukemia (CML) resistance against one or more BCR-ABL1 tyrosine kinase inhibitors (TKI) remains a clinical challenge. Preclinical data suggest that TKI combinations may overcome resistance. We report on a heavily pre-treated 78 year-old female patient with CML who developed multi-resistant blast crisis with bone marrow fibrosis and a Ph- clone. Treatment with ponatinib resulted in blast cell clearance, decrease in fibrosis, and disappearance of BCR-ABL1, but also in severe thrombocytopenia with bleedings requiring platelet transfusions. We therefore switched from ponatinib to bosutinib. During bosutinib, platelet counts recovered. However, after 6 months, BCR-ABL1 mRNA levels increased to > 1%. Therefore, we ´switched back´ to ponatinib, and this was again followed by disappearance of BCR-ABL1 and a decrease in platelets. During the next 2 years, we applied ponatinib and bosutinib in continuous rotation-cycles and added hydroxyurea in order to suppress all sub-clones and to balance between efficacy and potential side effects following the principle of personalized medicine. With this approach the patient remained in complete molecular response and reached normal blood counts and a normal quality of life without vascular or other side effects. In conclusion, TKI rotation is a novel potent approach to suppress multiple resistant sub-clones and to balance between clinical efficacy and side effects in patients with advanced CML. Clinical trials are now warranted to show that TKI-rotation is in general safe and effective in these patients.

Ensembl (http://www.ensembl.org) integrates genomic information for a comprehensive set of chordate genomes with a particular focus on resources for human, mouse, rat, zebrafish and other high-value sequenced genomes. We provide complete gene annotations for all supported species in addition to specific resources that target genome variation, function and evolution. Ensembl data is accessible in a variety of formats including via our genome browser, API and BioMart. This year marks the tenth anniversary of Ensembl and in that time the project has grown with advances in genome technology. As of release 56 (September 2009), Ensembl supports 51 species including marmoset, pig, zebra finch, lizard, gorilla and wallaby, which were added in the past year. Major additions and improvements to Ensembl since our previous report include the incorporation of the human GRCh37 assembly, enhanced visualisation and data-mining options for the Ensembl regulatory features and continued development of our software infrastructure.

Human respiratory syncytial virus (HRSV) is a major cause of pediatric lower respiratory tract disease to which there is no vaccine or efficacious chemotherapeutic strategy. Although RNA synthesis and virus assembly occur in the cytoplasm, HRSV is known to induce nuclear responses in the host cell as replication alters global gene expression. Quantitative proteomics was used to take an unbiased overview of the protein changes in transformed human alveolar basal epithelial cells infected with HRSV. Underpinning this was the use of stable isotope labeling with amino acids in cell culture coupled to LC-MS/MS, which allowed the direct and simultaneous identification and quantification of both cellular and viral proteins. To reduce sample complexity and increase data return on potential protein localization, cells were fractionated into nuclear and cytoplasmic extracts. This resulted in the identification of 1,140 cellular proteins and six viral proteins. The proteomics data were analyzed using Ingenuity Pathways Analysis to identify defined canonical pathways and functional groupings. Selected data were validated using Western blot, direct and indirect immunofluorescence confocal microscopy, and functional assays. The study served to validate and expand upon known HRSV-host cell interactions, including those associated with the antiviral response and alterations in subnuclear structures such as the nucleolus and ND10 (promyelocytic leukemia bodies). In addition, novel changes were observed in mitochondrial proteins and functions, cell cycle regulatory molecules, nuclear pore complex proteins and nucleocytoplasmic trafficking proteins. These data shed light into how the cell is potentially altered to create conditions more favorable for infection. Additionally, the study highlights the application and advantage of stable isotope labeling with amino acids in cell culture coupled to LC-MS/MS for the analysis of virus-host interactions.

SMNDC1 is a Tudor domain protein that recognizes di-methylated arginines and controls gene expression as an essential splicing factor. Here, we study the specific contributions of the SMNDC1 Tudor domain to protein-protein interactions, subcellular localization, and molecular function. To perturb the protein function in cells, we develop small molecule inhibitors targeting the dimethylarginine binding pocket of the SMNDC1 Tudor domain. We find that SMNDC1 localizes to phase-separated membraneless organelles that partially overlap with nuclear speckles. This condensation behavior is driven by the unstructured C-terminal region of SMNDC1, depends on RNA interaction and can be recapitulated in vitro. Inhibitors of the protein's Tudor domain drastically alter protein-protein interactions and subcellular localization, causing splicing changes for SMNDC1-dependent genes. These compounds will enable further pharmacological studies on the role of SMNDC1 in the regulation of nuclear condensates, gene regulation and cell identity.

Pathways modulating glucose homeostasis independently of insulin would open new avenues to combat insulin resistance and diabetes. Here, we report the establishment, characterization, and use of a vertebrate 'insulin-free' model to identify insulin-independent modulators of glucose metabolism. insulin knockout zebrafish recapitulate core characteristics of diabetes and survive only up to larval stages. Utilizing a highly efficient endoderm transplant technique, we generated viable chimeric adults that provide the large numbers of insulin mutant larvae required for our screening platform. Using glucose as a disease-relevant readout, we screened 2233 molecules and identified three that consistently reduced glucose levels in insulin mutants. Most significantly, we uncovered an insulin-independent beneficial role for androgen receptor antagonism in hyperglycemia, mostly by reducing fasting glucose levels. Our study proposes therapeutic roles for androgen signaling in diabetes and, more broadly, offers a novel in vivo model for rapid screening and decoupling of insulin-dependent and -independent mechanisms.

Ca2+-sensor proteins are generally implicated in insulin release through SNARE interactions. Here, secretagogin, whose expression in human pancreatic islets correlates with their insulin content and the incidence of type 2 diabetes, is shown to orchestrate an unexpectedly distinct mechanism. Single-cell RNA-seq reveals retained expression of the TRP family members in β-cells from diabetic donors. Amongst these, pharmacological probing identifies Ca2+-permeable transient receptor potential vanilloid type 1 channels (TRPV1) as potent inducers of secretagogin expression through recruitment of Sp1 transcription factors. Accordingly, agonist stimulation of TRPV1s fails to rescue insulin release from pancreatic islets of glucose intolerant secretagogin knock-out(-/-) mice. However, instead of merely impinging on the SNARE machinery, reduced insulin availability in secretagogin-/- mice is due to β-cell loss, which is underpinned by the collapse of protein folding and deregulation of secretagogin-dependent USP9X deubiquitinase activity. Therefore, and considering the desensitization of TRPV1s in diabetic pancreata, a TRPV1-to-secretagogin regulatory axis seems critical to maintain the structural integrity and signal competence of β-cells.

Transcription factor fusion proteins can transform cells by inducing global changes of the transcriptome, often creating a state of oncogene addiction. Here, we investigate the role of epigenetic mechanisms in this process, focusing on Ewing sarcoma cells that are dependent on the EWS-FLI1 fusion protein. We established reference epigenome maps comprising DNA methylation, seven histone marks, open chromatin states, and RNA levels, and we analyzed the epigenome dynamics upon downregulation of the driving oncogene. Reduced EWS-FLI1 expression led to widespread epigenetic changes in promoters, enhancers, and super-enhancers, and we identified histone H3K27 acetylation as the most strongly affected mark. Clustering of epigenetic promoter signatures defined classes of EWS-FLI1-regulated genes that responded differently to low-dose treatment with histone deacetylase inhibitors. Furthermore, we observed strong and opposing enrichment patterns for E2F and AP-1 among EWS-FLI1-correlated and anticorrelated genes. Our data describe extensive genome-wide rewiring of epigenetic cell states driven by an oncogenic fusion protein.

The Ensembl project (http://www.ensembl.org) provides genome information for sequenced chordate genomes with a particular focus on human, mouse, zebrafish and rat. Our resources include evidenced-based gene sets for all supported species; large-scale whole genome multiple species alignments across vertebrates and clade-specific alignments for eutherian mammals, primates, birds and fish; variation data resources for 17 species and regulation annotations based on ENCODE and other data sets. Ensembl data are accessible through the genome browser at http://www.ensembl.org and through other tools and programmatic interfaces.

Mice lacking Cdk6 kinase activity suffer from mild anemia accompanied by elevated numbers of Ter119+ cells in the bone marrow. The animals show hardly any alterations in erythroid development, indicating that Cdk6 is not required for proliferation and maturation of erythroid cells. There is also no difference in stress erythropoiesis following hemolysis in vivo However, Cdk6-/- erythrocytes have a shortened lifespan and are more sensitive to mechanical stress in vitro, suggesting differences in cytoskeletal architecture. Erythroblasts contain both Cdk4 and Cdk6, while mature erythrocytes apparently lack Cdk4 and their Cdk6 is partly associated with the cytoskeleton. We used mass spectrometry to show that Cdk6 interacts with a number of proteins involved in cytoskeleton organization. Cdk6-/- erythroblasts show impaired F-actin formation and lower levels of gelsolin, which interacts with Cdk6. We also found that Cdk6 regulates the transcription of a panel of genes involved in actin (de-)polymerization. Cdk6-deficient cells are sensitive to drugs that interfere with the cytoskeleton, suggesting that our findings are relevant to the treatment of patients with anemia - and may be relevant to cancer patients treated with the new generation of CDK6 inhibitors.

Superspreading events shaped the coronavirus disease 2019 (COVID-19) pandemic, and their rapid identification and containment are essential for disease control. Here, we provide a national-scale analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) superspreading during the first wave of infections in Austria, a country that played a major role in initial virus transmissions in Europe. Capitalizing on Austria's well-developed epidemiological surveillance system, we identified major SARS-CoV-2 clusters during the first wave of infections and performed deep whole-genome sequencing of more than 500 virus samples. Phylogenetic-epidemiological analysis enabled the reconstruction of superspreading events and charts a map of tourism-related viral spread originating from Austria in spring 2020. Moreover, we exploited epidemiologically well-defined clusters to quantify SARS-CoV-2 mutational dynamics, including the observation of low-frequency mutations that progressed to fixation within the infection chain. Time-resolved virus sequencing unveiled viral mutation dynamics within individuals with COVID-19, and epidemiologically validated infector-infectee pairs enabled us to determine an average transmission bottleneck size of 103 SARS-CoV-2 particles. In conclusion, this study illustrates the power of combining epidemiological analysis with deep viral genome sequencing to unravel the spread of SARS-CoV-2 and to gain fundamental insights into mutational dynamics and transmission properties.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Antiapoptotic Bcl-2 family members have recently (re)emerged as key drug targets in cancer, with a tissue- and tumor-specific activity profile of available BH3 mimetics. In multiple myeloma, MCL-1 has been described as a major gatekeeper of apoptosis. This discovery has led to the rapid establishment of clinical trials evaluating the impact of various MCL-1 inhibitors. However, our understanding about the clinical impact and optimal use of MCL-1 inhibitors is still limited. We therefore explored mechanisms of acquired MCL-1 inhibitor resistance and optimization strategies in myeloma. Our findings indicated heterogeneous paths to resistance involving baseline Bcl-2 family alterations of proapoptotic (BAK, BAX, and BIM) and antiapoptotic (Bcl-2 and MCL-1) proteins. These manifestations depend on the BH3 profile of parental cells that guide the enhanced formation of Bcl-2:BIM and/or the dynamic (ie, treatment-induced) formation of Bcl-xL:BIM and Bcl-xL:BAK complexes. Accordingly, an unbiased high-throughput drug-screening approach (n = 528) indicated alternative BH3 mimetics as top combination partners for MCL-1 inhibitors in sensitive and resistant cells (Bcl-xL>Bcl-2 inhibition), whereas established drug classes were mainly antagonistic (eg, antimitotic agents). We also revealed reduced activity of MCL-1 inhibitors in the presence of stromal support as a drug-class effect that was overcome by concurrent Bcl-xL or Bcl-2 inhibition. Finally, we demonstrated heterogeneous Bcl-2 family deregulation and MCL-1 inhibitor cross-resistance in carfilzomib-resistant cells, a phenomenon linked to the MDR1-driven drug efflux of MCL-1 inhibitors. The implications of our findings for clinical practice emphasize the need for patient-adapted treatment protocols, with the tracking of tumor- and/or clone-specific adaptations in response to MCL-1 inhibition.