Search

Abstract Transition of cytosine to thymine in CpG dinucleotides is the most frequent type of mutation in cancer. This increased mutability is commonly attributed to the spontaneous deamination of 5-methylcytosine (5mC), which is normally repaired by the base-excision repair (BER) pathway. However, the contribution of 5mC deamination in the increasing diversity of cancer mutational signatures remains poorly explored. We integrate mutational signatures analysis in a large series of tumor whole genomes with lineage-specific epigenomic data to draw a detailed view of 5mC deamination in cancer. We uncover tumor type-specific patterns of 5mC deamination signatures in CpG and non-CpG contexts. We demonstrate that the BER glycosylase MBD4 preferentially binds to active chromatin and early replicating DNA, which correlates with lower mutational burden in these domains. We validate our findings by modeling BER deficiencies in isogenic cell models. Here, we establish MBD4 as the main actor responsible for 5mC deamination repair in humans.

Abstract Aim Chemoresistance is a major cause of treatment failure in colorectal cancer (CRC) therapy. In this study, the impact of the IGF2BP family of RNA-binding proteins on CRC chemoresistance was investigated using in silico, in vitro, and in vivo approaches. Methods Gene expression data from a well-characterized cohort and publicly available cross-linking immunoprecipitation sequencing (CLIP-Seq) data were collected. Resistance to chemotherapeutics was assessed in patient-derived xenografts (PDXs) and patient-derived organoids (PDOs). Functional studies were performed in 2D and 3D cell culture models, including proliferation, spheroid growth, and mitochondrial respiration analyses. Results We identified IGF2BP2 as the most abundant IGF2BP in primary and metastastatic CRC, correlating with tumor stage in patient samples and tumor growth in PDXs. IGF2BP2 expression in primary tumor tissue was significantly associated with resistance to selumetinib, gefitinib, and regorafenib in PDOs and to 5-fluorouracil and oxaliplatin in PDX in vivo. IGF2BP2 knockout (KO) HCT116 cells were more susceptible to regorafenib in 2D and to oxaliplatin, selumitinib, and nintedanib in 3D cell culture. Further, a bioinformatic analysis using CLIP data suggested stabilization of target transcripts in primary and metastatic tumors. Measurement of oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) revealed a decreased basal OCR and an increase in glycolytic ATP production rate in IGF2BP2 KO. In addition, real-time reverse transcriptase polymerase chain reaction (qPCR) analysis confirmed decreased expression of genes of the respiratory chain complex I, complex IV, and the outer mitochondrial membrane in IGF2BP2 KO cells. Conclusions IGF2BP2 correlates with CRC tumor growth in vivo and promotes chemoresistance by altering mitochondrial respiratory chain metabolism. As a druggable target, IGF2BP2 could be used in future CRC therapy to overcome CRC chemoresistance.

Summary: The noncoding genome presents a largely untapped source of new biological insights, including thousands of long noncoding RNA (lncRNA) loci. While lncRNA dysregulation has been reported in myeloid malignancies, their functional relevance remains to be systematically interrogated. We performed CRISPRi screens of lncRNA signatures from normal and malignant hematopoietic cells and identified MYNRL15 as a myeloid leukemia dependency. Functional dissection suggests an RNA-independent mechanism mediated by two regulatory elements embedded in the locus. Genetic perturbation of these elements triggered a long-range chromatin interaction and downregulation of leukemia dependency genes near the gained interaction sites, as well as overall suppression of cancer dependency pathways. Thus, this study describes a new noncoding myeloid leukemia vulnerability and mechanistic concept for myeloid leukemia. Importantly, MYNRL15 perturbation caused strong and selective impairment of leukemia cells of various genetic backgrounds over normal hematopoietic stem and progenitor cells in vitro, and depletion of patient-derived xenografts in vivo.

5-Fluorouracil (5-FU) is a commonly used chemo drug to treat cancer patients. Using in vitro 5-FU treated organoids and in vivo 5-FU treated human cancer data, the authors describe a 5-FU mutational signature that is similar to COSMIC signature 17 frequently observed in untreated cancers

Transcriptional profiling of hematopoietic cell subpopulations has helped to characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases. Previously, only the genes targeted by expression microarrays could be profiled genome-wide. High-throughput RNA sequencing, however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analyzed the BLUEPRINT consortium RNA-sequencing data for mature hematopoietic cell types. The data comprised 90 total RNA-sequencing samples, each composed of one of 27 cell types, and 32 small RNA-sequencing samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type-dependent expression. We also characterized the expression of circular RNA and found that these are also cell type-specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematologic development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

Abstract Existing methods for paired antibody heavy- and light-chain repertoire sequencing rely on specialized equipment and are limited by their commercial availability and high costs. Here, we report a novel simple and cost-effective emulsion-based single-cell paired antibody repertoire sequencing method that employs only basic laboratory equipment. We performed a proof-of-concept using mixed mouse hybridoma cells and we also showed that our method can be used for discovery of novel antigen-specific monoclonal antibodies by sequencing human CD19+ B cell IgM and IgG repertoires isolated from peripheral whole blood before and seven days after Td (Tetanus toxoid/Diphtheria toxoid) booster immunization. We anticipate broad applicability of our method for providing insights into adaptive immune responses associated with various diseases, vaccinations, and cancer immunotherapies.

The heterogeneity of colorectal cancer has important clinical and therapeutic implications. Here the authors analysed the responses of a large biobank of organoids and xenografts derived from colorectal patients to a panel of clinically relevant therapeutic agents to identify genes signatures associated with drug response.

Congenital Cytomegalovirus infection (cCMV) is the leading infection in determining permanent long-term impairments (LTI), and its pathogenesis is largely unknown due to the complex interplay between viral, maternal, placental, and child factors. The cellular activity, considered to be the result of the response to exogenous and endogenous factors, is captured by the determination of gene expression profiles. In this study, we determined whole blood transcriptomes in relation to cCMV, CMV viral load and LTI development at 6 years of age by using RNA isolated from neonatal dried blood spots (DBS) stored at room temperature for 8 years. As DBS were assumed to mainly reflect the neonatal immune system, particular attention was given to the immune pathways using the global test. Additionally, differential expression of individual genes was performed using the voom/limma function packages. We demonstrated feasibility of RNA sequencing from archived neonatal DBS of children with cCMV, and non-infected controls, in relation to LTI and CMV viral load. Despite the lack of statistical power to detect individual genes differences, pathway analysis suggested the involvement of innate immune response with higher CMV viral loads, and of anti-inflammatory markers in infected children that did not develop LTI. Finally, the T cell exhaustion observed in infected neonates, in particular with higher viral load, did not correlate with LTI, therefore other mechanisms are likely to be involved in the long-term immune dysfunction. Despite these data demonstrate limitation in determining prognostic markers for LTI by means of transcriptome analysis, this exploratory study represents a first step in unraveling the pathogenesis of cCMV, and the aforementioned pathways certainly merit further evaluation.

The Prep1 (Pknox1) homeodomain transcription factor is essential at multiple stages of embryo development. In the E11.5 embryo trunk, we previously estimated that Prep1 binds about 3,300 genomic sites at a highly specific decameric consensus sequence, mainly governing basal cellular functions. We now show that in embryonic stem (ES) cells Prep1 binding pattern only partly overlaps that of the embryo trunk, with about 2,000 novel sites. Moreover, in ES cells Prep1 still binds mostly to promoters, as in total embryo trunk but, among the peaks bound exclusively in ES cells, the percentage of enhancers was three-fold higher. RNA-seq identifies about 1800 genes down-regulated in Prep1-/- ES cells which belong to gene ontology categories not enriched in the E11.5 Prep1i/i differentiated embryo, including in particular essential components of the Wnt and Fgf pathways. These data agree with aberrant Wnt and Fgf expression levels in the Prep1-/- ES cells with a deficient embryoid bodies (EBs) formation and differentiation. Re-establishment of the Prep1 level rescues the phenotype.

Large-scale genomic analyses of patient cohorts have revealed extensive heterogeneity between individual tumors, contributing to treatment failure and drug resistance. In malignant melanoma, heterogeneity is thought to arise as a consequence of the differentiation of melanoma-initiating cells that are defined by cell-surface markers like CD271 or CD133.Here we confirmed that the nerve growth factor receptor (CD271) is a crucial determinant of tumorigenicity, stem-like properties, heterogeneity and plasticity in melanoma cells. Stable shRNA mediated knock-down of CD271 in patient-derived melanoma cells abrogated their tumor-initiating and colony-forming capacity. A genome-wide expression profiling and gene-set enrichment analysis revealed novel connections of CD271 with melanoma-associated genes like CD133 and points to a neural crest stem cell (NCSC) signature lost upon CD271 knock-down. In a meta-analysis we have determined a shared set of 271 differentially regulated genes, linking CD271 to SOX10, a marker that specifies the neural crest. To dissect the connection of CD271 and CD133 we have analyzed 10 patient-derived melanoma-cell strains for cell-surface expression of both markers compared to established cell lines MeWo and A375. We found CD271+ cells in the majority of cell strains analyzed as well as in a set of 16 different patient-derived melanoma metastases. Strikingly, only 2/12 cell strains harbored a CD133+ sub-set that in addition comprised a fraction of cells of a CD271+/CD133+ phenotype. Those cells were found in the label-retaining fraction and in vitro deduced from CD271+ but not CD271 knock-down cells.Our present study provides a deeper insight into the regulation of melanoma cell properties and points CD271 out as a regulator of several melanoma-associated genes. Further, our data strongly suggest that CD271 is a crucial determinant of stem-like properties of melanoma cells like colony-formation and tumorigenicity.

BACKGROUND:Non-alcoholic fatty liver disease (NAFLD) has a broad spectrum of disease states ranging from mild steatosis characterized by an abnormal retention of lipids within liver cells to steatohepatitis (NASH) showing fat accumulation, inflammation, ballooning and degradation of hepatocytes, and fibrosis. Ultimately, steatohepatitis can result in liver cirrhosis and hepatocellular carcinoma. METHODOLOGY AND RESULTS:In this study we have analyzed three different mouse strains, A/J, C57BL/6J, and PWD/PhJ, that show different degrees of steatohepatitis when administered a 3,5-diethoxycarbonyl-1,4-dihydrocollidine (DDC) containing diet. RNA-Seq gene expression analysis, protein analysis and metabolic profiling were applied to identify differentially expressed genes/proteins and perturbed metabolite levels of mouse liver samples upon DDC-treatment. Pathway analysis revealed alteration of arachidonic acid (AA) and S-adenosylmethionine (SAMe) metabolism upon other pathways. To understand metabolic changes of arachidonic acid metabolism in the light of disease expression profiles a kinetic model of this pathway was developed and optimized according to metabolite levels. Subsequently, the model was used to study in silico effects of potential drug targets for steatohepatitis. CONCLUSIONS:We identified AA/eicosanoid metabolism as highly perturbed in DDC-induced mice using a combination of an experimental and in silico approach. Our analysis of the AA/eicosanoid metabolic pathway suggests that 5-hydroxyeicosatetraenoic acid (5-HETE), 15-hydroxyeicosatetraenoic acid (15-HETE) and prostaglandin D2 (PGD2) are perturbed in DDC mice. We further demonstrate that a dynamic model can be used for qualitative prediction of metabolic changes based on transcriptomics data in a disease-related context. Furthermore, SAMe metabolism was identified as being perturbed due to DDC treatment. Several genes as well as some metabolites of this module show differences between A/J and C57BL/6J on the one hand and PWD/PhJ on the other.

Ascertaining when and where genes are expressed is of crucial importance to understanding or predicting the physiological role of genes and proteins and how they interact to form the complex networks that underlie organ development and function. It is, therefore, crucial to determine on a genome-wide level, the spatio-temporal gene expression profiles at cellular resolution. This information is provided by colorimetric RNA in situ hybridization that can elucidate expression of genes in their native context and does so at cellular resolution. We generated what is to our knowledge the first genome-wide transcriptome atlas by RNA in situ hybridization of an entire mammalian organism, the developing mouse at embryonic day 14.5. This digital transcriptome atlas, the Eurexpress atlas (http://www.eurexpress.org), consists of a searchable database of annotated images that can be interactively viewed. We generated anatomy-based expression profiles for over 18,000 coding genes and over 400 microRNAs. We identified 1,002 tissue-specific genes that are a source of novel tissue-specific markers for 37 different anatomical structures. The quality and the resolution of the data revealed novel molecular domains for several developing structures, such as the telencephalon, a novel organization for the hypothalamus, and insight on the Wnt network involved in renal epithelial differentiation during kidney development. The digital transcriptome atlas is a powerful resource to determine co-expression of genes, to identify cell populations and lineages, and to identify functional associations between genes relevant to development and disease.

Matrin 3 (MATR3) is a highly conserved, inner nuclear matrix protein with two zinc finger domains and two RNA recognition motifs (RRM), whose function is largely unknown. Recently we found MATR3 to be phosphorylated by the protein kinase ATM, which activates the cellular response to double strand breaks in the DNA. Here, we show that MATR3 interacts in an RNA-dependent manner with several proteins with established roles in RNA processing, and maintains its interaction with RNA via its RRM2 domain. Deep sequencing of the bound RNA (RIP-seq) identified several small noncoding RNA species. Using microarray analysis to explore MATR3's role in transcription, we identified 77 transcripts whose amounts depended on the presence of MATR3. We validated this finding with nine transcripts which were also bound to the MATR3 complex. Finally, we demonstrated the importance of MATR3 for maintaining the stability of several of these mRNA species and conclude that it has a role in mRNA stabilization. The data suggest that the cellular level of MATR3, known to be highly regulated, modulates the stability of a group of gene transcripts.

Summary: Neutrophils are short-lived blood cells that play a critical role in host defense against infections. To better comprehend neutrophil functions and their regulation, we provide a complete epigenetic overview, assessing important functional features of their differentiation stages from bone marrow-residing progenitors to mature circulating cells. Integration of chromatin modifications, methylation, and transcriptome dynamics reveals an enforced regulation of differentiation, for cellular functions such as release of proteases, respiratory burst, cell cycle regulation, and apoptosis. We observe an early establishment of the cytotoxic capability, while the signaling components that activate these antimicrobial mechanisms are transcribed at later stages, outside the bone marrow, thus preventing toxic effects in the bone marrow niche. Altogether, these data reveal how the developmental dynamics of the chromatin landscape orchestrate the daily production of a large number of neutrophils required for innate host defense and provide a comprehensive overview of differentiating human neutrophils. : Grassi et al. report that the establishment of transcriptional enhancers drives neutrophil differentiation. Coordinated waves of gene expression establish the cytotoxic capability of these cells at early stages of maturation. A set of super-enhancers is specifically opened at the end of the differentiation process to control neutrophil activation. Keywords: neutrophil, epigenome, transcriptome, myeloid differentiation