Global Transcriptional Response to CRISPR/Cas9-AAV6-Based Genome Editing in CD34+ Hematopoietic Stem and Progenitor Cells.

Genome-editing technologies are currently being translated to the clinic. However, cellular effects of the editing machinery have yet to be fully elucidated. Here, we performed global microarray-based gene expression measurements on human CD34+ hematopoietic stem and progenitor cells that underwent editing. We probed effects of the entire editing process as well as each component individually, including electroporation, Cas9 (mRNA or protein) with chemically modified sgRNA, and AAV6 transduction. We identified differentially expressed genes relative to control treatments, which displayed enrichment for particular biological processes. All editing machinery components elicited immune, stress, and apoptotic responses. Cas9 mRNA invoked the greatest amount of transcriptional change, eliciting a distinct viral response and global transcriptional downregulation, particularly of metabolic and cell cycle processes. Electroporation also induced significant transcriptional change, with notable downregulation of metabolic processes. Surprisingly, AAV6 evoked no detectable viral response. We also found Cas9/sgRNA ribonucleoprotein treatment to be well tolerated, in spite of eliciting a DNA damage signature. Overall, this data establishes a benchmark for cellular tolerance of CRISPR/Cas9-AAV6-based genome editing, ensuring that the clinical protocol is as safe and efficient as possible.

Serum N-glycan analysis in breast cancer patients--Relation to tumour biology and clinical outcome.

Glycosylation and related processes play important roles in cancer development and progression, including metastasis. Several studies have shown that N-glycans have potential diagnostic value as cancer serum biomarkers. We have explored the significance of the abundance of particular serum N-glycan structures as important features of breast tumour biology by studying the serum glycome and tumour transcriptome (mRNA and miRNA) of 104 breast cancer patients. Integration of these types of molecular data allows us to study the relationship between serum glycans and transcripts representing functional pathways, such as metabolic pathways or DNA damage response. We identified tri antennary trigalactosylated trisialylated glycans in serum as being associated with lower levels of tumour transcripts involved in focal adhesion and integrin-mediated cell adhesion. These glycan structures were also linked to poor prognosis in patients with ER negative tumours. High abundance of simple monoantennary glycan structures were associated with increased survival, particularly in the basal-like subgroup. The presence of circulating tumour cells was found to be significantly associated with several serum glycome structures like bi and triantennary, di- and trigalactosylated, di- and trisialylated. The link between tumour miRNA expression levels and N-glycan production is also examined.

Chemically modified guide RNAs enhance CRISPR-Cas genome editing in human primary cells.

CRISPR-Cas-mediated genome editing relies on guide RNAs that direct site-specific DNA cleavage facilitated by the Cas endonuclease. Here we report that chemical alterations to synthesized single guide RNAs (sgRNAs) enhance genome editing efficiency in human primary T cells and CD34(+) hematopoietic stem and progenitor cells. Co-delivering chemically modified sgRNAs with Cas9 mRNA or protein is an efficient RNA- or ribonucleoprotein (RNP)-based delivery method for the CRISPR-Cas system, without the toxicity associated with DNA delivery. This approach is a simple and effective way to streamline the development of genome editing with the potential to accelerate a wide array of biotechnological and therapeutic applications of the CRISPR-Cas technology.

ENViz: a Cytoscape App for integrated statistical analysis and visualization of sample-matched data with multiple data types.

ENViz (Enrichment Analysis and Visualization) is a Cytoscape app that performs joint enrichment analysis of two types of sample matched datasets in the context of systematic annotations. Such datasets may be gene expression or any other high-throughput data collected in the same set of samples. The enrichment analysis is done in the context of pathway information, gene ontology or any custom annotation of the data. The results of the analysis consist of significant associations between profiled elements of one of the datasets to the annotation terms (e.g. miR-19 was associated to the cell-cycle process in breast cancer samples). The results of the enrichment analysis are visualized as an interactive Cytoscape network.

Specific genomic and transcriptomic aberrations in tumors induced by partial hepatectomy of a chronically inflamed murine liver.

Resection of hepatocellular carcinoma (HCC) tumors by partial hepatectomy (PHx) is associated with promoting hepatocarcinogenesis. We have previously reported that PHx promotes hepatocarcinogenesis in the Mdr2-knockout (Mdr2-KO) mouse, a model for inflammation-mediated HCC. Now, to explore the molecular mechanisms underlying the tumor-promoting effect of PHx, we compared genomic and transcriptomic profiles of HCC tumors developing in the Mdr2-KO mice either spontaneously or following PHx. PHx accelerated HCC development in these mice by four months. PHx-induced tumors had major chromosomal aberrations: all were amplifications affecting multiple chromosomes. Most of these amplifications were located near the acrocentric centromeres of murine chromosomes. Four different chromosomal regions were amplified each in at least three tumors. The human orthologs of these common amplified regions are known to be amplified in HCC. All tumors of untreated mice had chromosomal aberrations, including both deletions and amplifications. Amplifications in spontaneous tumors affected fewer chromosomes and were not located preferentially at the chromosomal edges. Comparison of gene expression profiles revealed a significantly enriched expression of oncogenes, chromosomal instability markers and E2F1 targets in the post-PHx compared to spontaneous tumors. Both tumor groups shared the same frequent amplification at chromosome 18. Here, we revealed that one of the regulatory genes encoded by this amplified region, Crem, was over-expressed in the nuclei of murine and human HCC cells in vivo, and that it stimulated proliferation of human HCC cells in vitro. Our results demonstrate that PHx of a chronically inflamed liver directed tumor development to a discrete pathway characterized by amplification of specific chromosomal regions and expression of specific tumor-promoting genes. Crem is a new candidate HCC oncogene frequently amplified in this model and frequently over-expressed in human HCC.

Association of N-glycosylation with breast carcinoma and systemic features using high-resolution quantitative UPLC.

An improved separation of the human serum N-glycome using hydrophilic interaction chromatography technology with UPLC is described, where more than 140 N-glycans were assigned. Using this technique, serum samples from 107 healthy controls and 62 newly diagnosed breast cancer patients were profiled. The most statistically significant alterations were observed in cancer patients compared with healthy controls: an increase in sialylation, branching, and outer-arm fucosylation and a decrease in high-mannosylated and biantennary core-fucosylated glycans. In the controls and cases combined systemic features were analyzed; serum estradiol was associated with increase in digalactosylated glycans, and higher mammographic density was associated with increase in biantennary digalactosylated glycans and with decrease in trisialylated and in outer-arm fucosylated glycans. Furthermore, particular glycans were altered in some features of the breast carcinomas; bisected biantennary nonfucosylated glycans were decreased in patients with progesterone receptor positive tumors, and core-fucosylated biantennary bisected monogalactosylated glycans were decreased in patients with the TP53 mutation. Systemic features show more significant associations with the serum N-glycome than do the features of the breast carcinomas. In conclusion, the UPLC-based glycan analysis technique described here reveals highly significant differences between healthy women and breast cancer patients. Significant associations with breast carcinoma and systemic features are described.

Molecular rules governing de novo methylation in cancer.

De novo methylation of CpG islands is seen in many cancers, but the general rules governing this process are not known. By analyzing DNA from tumors, as well as normal tissues, and by utilizing a range of published data, we have identified a universal set of tumor targets, each with its own "coefficient" of methylation that is largely correlated with its inherent relative ability to recruit polycomb. This pattern is initially formed by a slow process of de novo methylation that occurs during aging and then undergoes expansion early in tumorigenesis, where we hypothesize that it may act as an inhibitor of development-associated gene activation.

Deregulation of cancer-related miRNAs is a common event in both benign and malignant human breast tumors.

MicroRNAs (miRNAs) are endogenous non-coding RNAs, which play an essential role in the regulation of gene expression during carcinogenesis. The role of miRNAs in breast cancer has been thoroughly investigated, and although many miRNAs are identified as cancer related, little is known about their involvement in benign tumors. In this study, we investigated miRNA expression profiles in the two most common types of human benign tumors (fibroadenoma/fibroadenomatosis) and in malignant breast tumors and explored their role as oncomirs and tumor suppressor miRNAs. Here, we identified 33 miRNAs with similar deregulated expression in both benign and malignant tumors compared with the expression levels of those in normal tissue, including breast cancer-related miRNAs such as let-7, miR-21 and miR-155. Additionally, messenger RNA (mRNA) expression profiles were obtained for some of the same samples. Using integrated mRNA/miRNA expression analysis, we observed that overexpression of certain miRNAs co-occurred with a significant downregulation of their candidate target mRNAs in both benign and malignant tumors. In support of these findings, in vitro functional screening of the downregulated miRNAs in non-malignant and breast cancer cell lines identified several possible tumor suppressor miRNAs, including miR-193b, miR-193a-3p, miR-126, miR-134, miR-132, miR-486-5p, miR-886-3p, miR-195 and miR-497, showing reduced growth when re-expressed in cancer cells. The finding of deregulated expression of oncomirs and tumor suppressor miRNAs in benign breast tumors is intriguing, indicating that they may play a role in proliferation. A role of cancer-related miRNAs in the early phases of carcinogenesis and malignant transformation can, therefore, not be ruled out.

Identifying in-trans process associated genes in breast cancer by integrated analysis of copy number and expression data.

Genomic copy number alterations are common in cancer. Finding the genes causally implicated in oncogenesis is challenging because the gain or loss of a chromosomal region may affect a few key driver genes and many passengers. Integrative analyses have opened new vistas for addressing this issue. One approach is to identify genes with frequent copy number alterations and corresponding changes in expression. Several methods also analyse effects of transcriptional changes on known pathways. Here, we propose a method that analyses in-cis correlated genes for evidence of in-trans association to biological processes, with no bias towards processes of a particular type or function. The method aims to identify cis-regulated genes for which the expression correlation to other genes provides further evidence of a network-perturbing role in cancer. The proposed unsupervised approach involves a sequence of statistical tests to systematically narrow down the list of relevant genes, based on integrative analysis of copy number and gene expression data. A novel adjustment method handles confounding effects of co-occurring copy number aberrations, potentially a large source of false positives in such studies. Applying the method to whole-genome copy number and expression data from 100 primary breast carcinomas, 6373 genes were identified as commonly aberrant, 578 were highly in-cis correlated, and 56 were in addition associated in-trans to biological processes. Among these in-trans process associated and cis-correlated (iPAC) genes, 28% have previously been reported as breast cancer associated, and 64% as cancer associated. By combining statistical evidence from three separate subanalyses that focus respectively on copy number, gene expression and the combination of the two, the proposed method identifies several known and novel cancer driver candidates. Validation in an independent data set supports the conclusion that the method identifies genes implicated in cancer.

miRNA target enrichment analysis reveals directly active miRNAs in health and disease.

microRNAs (miRNAs) are short non-coding regulatory RNA molecules. The activity of a miRNA in a biological process can often be reflected in the expression program that characterizes the outcome of the activity. We introduce a computational approach that infers such activity from high-throughput data using a novel statistical methodology, called minimum-mHG (mmHG), that examines mutual enrichment in two ranked lists. Based on this methodology, we provide a user-friendly web application that supports the statistical assessment of miRNA target enrichment analysis (miTEA) in the top of a ranked list of genes or proteins. Using miTEA, we analyze several target prediction tools by examining performance on public miRNA constitutive expression data. We also apply miTEA to analyze several integrative biology data sets, including a novel matched miRNA/mRNA data set covering nine human tissue types. Our novel findings include proposed direct activity of miR-519 in placenta, a direct activity of the oncogenic miR-15 in different healthy tissue types and a direct activity of the poorly characterized miR-768 in both healthy tissue types and cancer cell lines. The miTEA web application is available at http://cbl-gorilla.cs.technion.ac.il/miTEA/.

Global methylation patterns in idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic Pulmonary Fibrosis (IPF) is characterized by profound changes in the lung phenotype including excessive extracellular matrix deposition, myofibroblast foci, alveolar epithelial cell hyperplasia and extensive remodeling. The role of epigenetic changes in determining the lung phenotype in IPF is unknown. In this study we determine whether IPF lungs exhibit an altered global methylation profile. METHODOLOGY/PRINCIPAL FINDINGS: Immunoprecipitated methylated DNA from 12 IPF lungs, 10 lung adenocarcinomas and 10 normal histology lungs was hybridized to Agilent human CpG Islands Microarrays and data analysis was performed using BRB-Array Tools and DAVID Bioinformatics Resources software packages. Array results were validated using the EpiTYPER MassARRAY platform for 3 CpG islands. 625 CpG islands were differentially methylated between IPF and control lungs with an estimated False Discovery Rate less than 5%. The genes associated with the differentially methylated CpG islands are involved in regulation of apoptosis, morphogenesis and cellular biosynthetic processes. The expression of three genes (STK17B, STK3 and HIST1H2AH) with hypomethylated promoters was increased in IPF lungs. Comparison of IPF methylation patterns to lung cancer or control samples, revealed that IPF lungs display an intermediate methylation profile, partly similar to lung cancer and partly similar to control with 402 differentially methylated CpG islands overlapping between IPF and cancer. Despite their similarity to cancer, IPF lungs did not exhibit hypomethylation of long interspersed nuclear element 1 (LINE-1) retrotransposon while lung cancer samples did, suggesting that the global hypomethylation observed in cancer was not typical of IPF. CONCLUSIONS/SIGNIFICANCE: Our results provide evidence that epigenetic changes in IPF are widespread and potentially important. The partial similarity to cancer may signify similar pathogenetic mechanisms while the differences constitute IPF or cancer specific changes. Elucidating the role of these specific changes will potentially allow better understanding of the pathogenesis of IPF.

miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors.

INTRODUCTION: Few studies have performed expression profiling of both miRNA and mRNA from the same primary breast carcinomas. In this study we present and analyze data derived from expression profiling of 799 miRNAs in 101 primary human breast tumors, along with genome-wide mRNA profiles and extensive clinical information. METHODS: We investigate the relationship between these molecular components, in terms of their correlation with each other and with clinical characteristics. We use a systems biology approach to examine the correlative relationship between miRNA and mRNAs using statistical enrichment methods. RESULTS: We identify statistical significant differential expression of miRNAs between molecular intrinsic subtypes, and between samples with different levels of proliferation. Specifically, we point to miRNAs significantly associated with TP53 and ER status. We also show that several cellular processes, such as proliferation, cell adhesion and immune response, are strongly associated with certain miRNAs. We validate the role of miRNAs in regulating proliferation using high-throughput lysate-microarrays on cell lines and point to potential drivers of this process. CONCLUSION: This study provides a comprehensive dataset as well as methods and system-level results that jointly form a basis for further work on understanding the role of miRNA in primary breast cancer.

EGF decreases the abundance of microRNAs that restrain oncogenic transcription factors.

Epidermal growth factor (EGF) stimulates cells by launching gene expression programs that are frequently deregulated in cancer. MicroRNAs, which attenuate gene expression by binding complementary regions in messenger RNAs, are broadly implicated in cancer. Using genome-wide approaches, we showed that EGF stimulation initiates a coordinated transcriptional program of microRNAs and transcription factors. The earliest event involved a decrease in the abundance of a subset of 23 microRNAs. This step permitted rapid induction of oncogenic transcription factors, such as c-FOS, encoded by immediate early genes. In line with roles as suppressors of EGF receptor (EGFR) signaling, we report that the abundance of this early subset of microRNAs is decreased in breast and in brain tumors driven by the EGFR or the closely related HER2. These findings identify specific microRNAs as attenuators of growth factor signaling and oncogenesis.

Novel rank-based statistical methods reveal microRNAs with differential expression in multiple cancer types.

BACKGROUND: MicroRNAs (miRNAs) regulate target genes at the post-transcriptional level and play important roles in cancer pathogenesis and development. Variation amongst individuals is a significant confounding factor in miRNA (or other) expression studies. The true character of biologically or clinically meaningful differential expression can be obscured by inter-patient variation. In this study we aim to identify miRNAs with consistent differential expression in multiple tumor types using a novel data analysis approach. METHODS: Using microarrays we profiled the expression of more than 700 miRNAs in 28 matched tumor/normal samples from 8 different tumor types (breast, colon, liver, lung, lymphoma, ovary, prostate and testis). This set is unique in putting emphasis on minimizing tissue type and patient related variability using normal and tumor samples from the same patient. We develop scores for comparing miRNA expression in the above matched sample data based on a rigorous characterization of the distribution of order statistics over a discrete state set, including exact p-values. Specifically, we compute a Rank Consistency Score (RCoS) for every miRNA measured in our data. Our methods are also applicable in various other contexts. We compare our methods, as applied to matched samples, to paired t-test and to the Wilcoxon Signed Rank test. RESULTS: We identify consistent (across the cancer types measured) differentially expressed miRNAs. 41 miRNAs are under-expressed in cancer compared to normal, at FDR (False Discovery Rate) of 0.05 and 17 are over-expressed at the same FDR level. Differentially expressed miRNAs include known oncomiRs (e.g miR-96) as well as miRNAs that were not previously universally associated with cancer. Specific examples include miR-133b and miR-486-5p, which are consistently down regulated and mir-629* which is consistently up regulated in cancer, in the context of our cohort. Data is available in GEO. Software is available at: http://bioinfo.cs.technion.ac.il/people/zohar/RCoS/