The m6A epitranscriptome: transcriptome plasticity in brain development and function.

The field of epitranscriptomics examines the recently deciphered form of gene expression regulation that is mediated by type- and site-specific RNA modifications. Similarly to the role played by epigenetic mechanisms - which operate via DNA and histone modifications - epitranscriptomic modifications are involved in the control of the delicate gene expression patterns that are needed for the development and activity of the nervous system and are essential for basic and higher brain functions. Here we describe the mechanisms that are involved in the writing, erasing and reading of N6-methyladenosine, the most prevalent internal mRNA modification, and the emerging roles played by N6-methyladenosine in the nervous system.

The dynamic N(1)-methyladenosine methylome in eukaryotic messenger RNA.

Gene expression can be regulated post-transcriptionally through dynamic and reversible RNA modifications. A recent noteworthy example is N(6)-methyladenosine (m(6)A), which affects messenger RNA (mRNA) localization, stability, translation and splicing. Here we report on a new mRNA modification, N(1)-methyladenosine (m(1)A), that occurs on thousands of different gene transcripts in eukaryotic cells, from yeast to mammals, at an estimated average transcript stoichiometry of 20% in humans. Employing newly developed sequencing approaches, we show that m(1)A is enriched around the start codon upstream of the first splice site: it preferentially decorates more structured regions around canonical and alternative translation initiation sites, is dynamic in response to physiological conditions, and correlates positively with protein production. These unique features are highly conserved in mouse and human cells, strongly indicating a functional role for m(1)A in promoting translation of methylated mRNA.

Whole exome sequencing (WES) approach for diagnosing primary immunodeficiencies (PIDs) in a highly consanguineous community.

Primary immunodeficiencies (PIDs) are a heterogeneous group of monogenic inborn errors of immunity. The genetic causes of these diseases can be identified using whole exome sequencing (WES). Here, DNA samples from 106 patients with a clinical suspicion of PID were subjected to WES in order to test the diagnostic yield of this test in a highly consanguineous community. A likely genetic diagnosis was achieved in 70% of patients. Several factors were considered to possibly influence the diagnostic rate of WES among our cohort including early age, presence of consanguinity, family history suggestive of PID, the number of family members who underwent WES and the clinical phenotype of the patient. The highest diagnostic rate was in patients with combined immunodeficiency or with a syndrome. Notably, WES findings altered the clinical management in 39% (41/106) of patients in our cohort. Our findings support the use of WES as an important diagnostic tool in patients with suspected PID, especially in highly consanguineous communities.

Nm-seq maps 2'-O-methylation sites in human mRNA with base precision.

The ribose of RNA nucleotides can be 2'-O-methylated (Nm). Despite advances in high-throughput detection, the inert chemical nature of Nm still limits sensitivity and precludes mapping in mRNA. We leveraged the differential reactivity of 2'-O-methylated and 2'-hydroxylated nucleosides to periodate oxidation to develop Nm-seq, a sensitive method for transcriptome-wide mapping of Nm with base precision. Nm-seq uncovered thousands of Nm sites in human mRNA with features suggesting functional roles.

Newly diagnosed multiple myeloma patients carrying monoallelic deletion of the whole locus of immunoglobulin heavy chain gene have a better prognosis compared to those with t(4;14) and t(14;16).

The current study evaluated the prognostic significance of the monoallelic deletion of the whole locus of the immunoglobulin heavy-chain (w_del(IGH)) gene compared to translocations t(4;14) and t(14;16) among newly diagnosed multiple myeloma (MM) patients. We retrospectively analyzed clinical (age, gender, and staging) and laboratory data at diagnosis and the overall survival (OS) of 255 newly diagnosed MM patients carrying w_del(IGH) or translocations t(4;14) or t(14;16). Bone marrow samples were examined by morphological and sequential interphase fluorescense in situ hybridization analyses. Among 255 patients, 117 (45.8%) had w_del(IGH), 99 (38.8%) had t(4;14), and 39 (15.3%) had t(14;16). Mean age was 61.6 ± 11.6 years. Groups did not differ significantly in age, gender, or lactate dehydrogenase levels. Patients in the w_del(IGH) group presented more frequently at International Staging System stage I than at stage II/III. Patients in the w_del(IGH) group had significantly fewer additional chromosomal aberrations (1.58) than the other two groups (2.3 and 2.13 in the del(IGH), t(14;16) and t(4;14) groups, respectively, P < 0.0001). Furthermore, the w_del(IGH) group had significantly longer estimated median OS (9.47 years) compared to those with translocations t(14;16) (3.02 years, P = 0.002) or t(4;14) (4.18 years, P = 0.001), respectively. These findings suggest a potential prognostic significance of monoallelic deletion of IGH among these patients. Additional studies are needed to better understand the nature and mechanism of this prognostic factor.

The Duffy antigen receptor for chemokines, ACKR1,- 'Jeanne DARC' of benign neutropenia.

Benign neutropenia, observed in different ethnic groups, is the most common form of neutropenia worldwide. A specific single nucleotide polymorphism, rs2814778, located at the promoter of the ACKR1 (previously termed DARC) gene, which disrupts a binding site for the GATA1 erythroid transcription factor, resulting in a ACKR1-null phenotype, was found to serve as a predictor of low white blood cell and neutrophil counts in African-Americans and Yemenite Jews. Individuals with benign neutropenia due to the ACKR1-null allele have been found to have an increased susceptibility to human immunodeficiency virus infection and, on the other hand, a protective effect against malaria. The associated protective effect may explain the spread of the ACKR1-null allele by natural selection. The reviewed relationships between ACKR1 polymorphism and various pathological states may have important clinical implications to individuals with and without benign neutropenia. Potential mechanisms for ACKR1 (previously termed DARC) modulation during neutrophil recruitment to inflammation, and chemokine bioavailability in the circulation and in local tissue are reviewed and discussed.

Cerebral and portal vein thrombosis, macrocephaly and atypical absence seizures in Glycosylphosphatidyl inositol deficiency due to a PIGM promoter mutation.

Defects of the glycosylphosphatidylinositol (GPI) biosynthesis pathway constitute an emerging subgroup of congenital disorders of glycosylation with heterogeneous phenotypes. A mutation in the promoter of PIGM, resulting in a syndrome with portal vein thrombosis and persistent absence seizures, was previously described in three patients. We now report four additional patients in two unrelated families, with further clinical, biochemical and molecular delineation of this unique entity. We also describe the first prenatal diagnosis of PIGM deficiency, allowing characterization of the natural history of the disease from birth. The patients described herein expand the phenotypic spectrum of PIGM deficiency to include macrocephaly and infantile-onset cerebrovascular thrombotic events. Finally, we offer insights regarding targeted treatment of this rare disorder with sodium phenylbutyrate.

Disruption of Thrombocyte and T Lymphocyte Development by a Mutation in ARPC1B.

Regulation of the actin cytoskeleton is crucial for normal development and function of the immune system, as evidenced by the severe immune abnormalities exhibited by patients bearing inactivating mutations in the Wiskott-Aldrich syndrome protein (WASP), a key regulator of actin dynamics. WASP exerts its effects on actin dynamics through a multisubunit complex termed Arp2/3. Despite the critical role played by Arp2/3 as an effector of WASP-mediated control over actin polymerization, mutations in protein components of the Arp2/3 complex had not previously been identified as a cause of immunodeficiency. Here, we describe two brothers with hematopoietic and immunologic symptoms reminiscent of Wiskott-Aldrich syndrome (WAS). However, these patients lacked mutations in any of the genes previously associated with WAS. Whole-exome sequencing revealed a homozygous 2 bp deletion, n.c.G623DEL-TC (p.V208VfsX20), in Arp2/3 complex component ARPC1B that causes a frame shift resulting in premature termination. Modeling of the disease in zebrafish revealed that ARPC1B plays a critical role in supporting T cell and thrombocyte development. Moreover, the defects in development caused by ARPC1B loss could be rescued by the intact human ARPC1B ortholog, but not by the p.V208VfsX20 variant identified in the patients. Moreover, we found that the expression of ARPC1B is restricted to hematopoietic cells, potentially explaining why a mutation in ARPC1B has now been observed as a cause of WAS, whereas mutations in other, more widely expressed, components of the Arp2/3 complex have not been observed.

Two phases of mitogenic signaling unveil roles for p53 and EGR1 in elimination of inconsistent growth signals.

Normal cells require continuous exposure to growth factors in order to cross a restriction point and commit to cell-cycle progression. This can be replaced by two short, appropriately spaced pulses of growth factors, where the first pulse primes a process, which is completed by the second pulse, and enables restriction point crossing. Through integration of comprehensive proteomic and transcriptomic analyses of each pulse, we identified three processes that regulate restriction point crossing: (1) The first pulse induces essential metabolic enzymes and activates p53-dependent restraining processes. (2) The second pulse eliminates, via the PI3K/AKT pathway, the suppressive action of p53, as well as (3) sets an ERK-EGR1 threshold mechanism, which digitizes graded external signals into an all-or-none decision obligatory for S phase entry. Together, our findings uncover two gating mechanisms, which ensure that cells ignore fortuitous growth factors and undergo proliferation only in response to consistent mitogenic signals.

Somatic NRAS mutation in patient with generalized lymphatic anomaly.

Generalized lymphatic anomaly (GLA or lymphangiomatosis) is a rare disease characterized by a diffuse proliferation of lymphatic vessels in skin and internal organs. It often leads to progressive respiratory failure and death, but its etiology is unknown. Here, we isolated lymphangiomatosis endothelial cells from GLA tissue. These cells were characterized by high proliferation and survival rates, but displayed impaired capacities for migration and tube formation. We employed whole exome sequencing to search for disease-causing genes and identified a somatic mutation in NRAS. We used mouse and zebrafish model systems to initially evaluate the role of this mutation in the development of the lymphatic system, and we studied the effect of drugs blocking the downstream effectors, mTOR and ERK, on this disease.

Deleterious variants in TRAK1 disrupt mitochondrial movement and cause fatal encephalopathy.

Cellular distribution and dynamics of mitochondria are regulated by several motor proteins and a microtubule network. In neurons, mitochondrial trafficking is crucial because of high energy needs and calcium ion buffering along axons to synapses during neurotransmission. The trafficking kinesin proteins (TRAKs) are well characterized for their role in lysosomal and mitochondrial trafficking in cells, especially neurons. Using whole exome sequencing, we identified homozygous truncating variants in TRAK1 (NM_001042646:c.287-2A > C), in six lethal encephalopathic patients from three unrelated families. The pathogenic variant results in aberrant splicing and significantly reduced gene expression at the RNA and protein levels. In comparison with normal cells, TRAK1-deficient fibroblasts showed irregular mitochondrial distribution, altered mitochondrial motility, reduced mitochondrial membrane potential, and diminished mitochondrial respiration. This study confirms the role of TRAK1 in mitochondrial dynamics and constitutes the first report of this gene in association with a severe neurodevelopmental disorder.

Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq.

An extensive repertoire of modifications is known to underlie the versatile coding, structural and catalytic functions of RNA, but it remains largely uncharted territory. Although biochemical studies indicate that N(6)-methyladenosine (m(6)A) is the most prevalent internal modification in messenger RNA, an in-depth study of its distribution and functions has been impeded by a lack of robust analytical methods. Here we present the human and mouse m(6)A modification landscape in a transcriptome-wide manner, using a novel approach, m(6)A-seq, based on antibody-mediated capture and massively parallel sequencing. We identify over 12,000 m(6)A sites characterized by a typical consensus in the transcripts of more than 7,000 human genes. Sites preferentially appear in two distinct landmarks--around stop codons and within long internal exons--and are highly conserved between human and mouse. Although most sites are well preserved across normal and cancerous tissues and in response to various stimuli, a subset of stimulus-dependent, dynamically modulated sites is identified. Silencing the m(6)A methyltransferase significantly affects gene expression and alternative splicing patterns, resulting in modulation of the p53 (also known as TP53) signalling pathway and apoptosis. Our findings therefore suggest that RNA decoration by m(6)A has a fundamental role in regulation of gene expression.

e23D: database and visualization of A-to-I RNA editing sites mapped to 3D protein structures.

UNLABELLED: e23D, a database of A-to-I RNA editing sites from human, mouse and fly mapped to evolutionary related protein 3D structures, is presented. Genomic coordinates of A-to-I RNA editing sites are converted to protein coordinates and mapped onto 3D structures from PDB or theoretical models from ModBase. e23D allows visualization of the protein structure, modeling of recoding events and orientation of the editing with respect to nearby genomic functional sites from databases of disease causing mutations and genomic polymorphism. AVAILABILITY AND IMPLEMENTATION: http://www.sheba-cancer.org.il/e23D CONTACT: oz.solomon@live.biu.ac.il or Eran.Eyal@sheba.health.gov.il.

ADAR1 deletion induces NFκB and interferon signaling dependent liver inflammation and fibrosis.

Adenosine deaminase acting on RNA (ADAR) 1 binds and edits double-stranded (ds) RNA secondary structures found mainly within untranslated regions of many transcripts. In the current research, our aim was to study the role of ADAR1 in liver homeostasis. As previous studies show a conserved immunoregulatory function for ADAR1 in mammalians, we focused on its role in preventing chronic hepatic inflammation and the associated activation of hepatic stellate cells to produce extracellular matrix and promote fibrosis. We show that hepatocytes specific ADAR1 knock out (KO) mice display massive liver damage with multifocal inflammation and fibrogenesis. The bioinformatics analysis of the microarray gene-expression datasets of ADAR1 KO livers reveled a type-I interferons signature and an enrichment for immune response genes compared to control littermate livers. Furthermore, we found that in vitro silencing of ADAR1 expression in HepG2 cells leads to enhanced transcription of NFκB target genes, foremost of the pro-inflammatory cytokines IL6 and IL8. We also discovered immune cell-independent paracrine signaling among ADAR1-depleted HepG2 cells and hepatic stellate cells, leading to the activation of the latter cell type to adopt a profibrogenic phenotype. This paracrine communication dependent mainly on the production and secretion of the cytokine IL6 induced by ADAR1 silencing in hepatocytes. Thus, our findings shed a new light on the vital regulatory role of ADAR1 in hepatic immune homeostasis, chiefly its inhibitory function on the crosstalk between the NFκB and type-I interferons signaling cascades, restraining the development of liver inflammation and fibrosis.

Translocation t(11;14) in newly diagnosed patients with multiple myeloma: Is it always favorable?

The most common translocation in multiple myeloma (MM) is t(11;14)(q13;q32). According to several studies, this translocation represents a unique subset of patients with relatively favorable outcomes. Using combined analyses of morphology and fluorescence in situ hybridization (I-FISH), we examined the co-occurrence rates of t(11;14) with seven chromosomal aberrations (CAs), del(13q), del(17p), del(1p), gain(1q), multiple gains(1q), del(16q), and del(IGH), and assessed the effect of the different combinations on patient outcomes, with overall survival (OS) as the main outcome measure. Bone marrow samples and clinical data from 212 patients with MM with t(11;14) were analyzed. At least two additional CAs were found in 35% (75/205) of patients and a strong correlation between specific CAs. The occurrence of three CAs [multiple gains of (1q) (HR = 6.94, P = 0.001), del(1p) (HR = 4.47, P = 0.008), and del(IGH) (HR = 2.38, P = 0.002)] exerted a profoundly deleterious effect on median OS when compared with patients with t(11;14) only. Del(17p) and del(13q) have also exerted a deleterious effect albeit to a lesser extent (HR = 2.05, P = 0.07 and HR = 1.81, P = 0.03, respectively). When compared with t(11;14) alone, the addition of certain CAs lead to worse outcomes. These findings may have important clinical and biological implications. Patients with coexisting adverse lesions and t(11;14) may be considered at high risk and managed accordingly. © 2016 Wiley Periodicals, Inc.

G23D: Online tool for mapping and visualization of genomic variants on 3D protein structures.

BACKGROUND: Evaluation of the possible implications of genomic variants is an increasingly important task in the current high throughput sequencing era. Structural information however is still not routinely exploited during this evaluation process. The main reasons can be attributed to the partial structural coverage of the human proteome and the lack of tools which conveniently convert genomic positions, which are the frequent output of genomic pipelines, to proteins and structure coordinates. RESULTS: We present G23D, a tool for conversion of human genomic coordinates to protein coordinates and protein structures. G23D allows mapping of genomic positions/variants on evolutionary related (and not only identical) protein three dimensional (3D) structures as well as on theoretical models. By doing so it significantly extends the space of variants for which structural insight is feasible. To facilitate interpretation of the variant consequence, pathogenic variants, functional sites and polymorphism sites are displayed on protein sequence and structure diagrams alongside the input variants. G23D also provides modeling of the mutant structure, analysis of intra-protein contacts and instant access to functional predictions and predictions of thermo-stability changes. G23D is available at http://www.sheba-cancer.org.il/G23D . CONCLUSIONS: G23D extends the fraction of variants for which structural analysis is applicable and provides better and faster accessibility for structural data to biologists and geneticists who routinely work with genomic information.

A congenital neutrophil defect syndrome associated with mutations in VPS45.

BACKGROUND: Neutrophils are the predominant phagocytes that provide protection against bacterial and fungal infections. Genetically determined neutrophil disorders confer a predisposition to severe infections and reveal novel mechanisms that control vesicular trafficking, hematopoiesis, and innate immunity. METHODS: We clinically evaluated seven children from five families who had neutropenia, neutrophil dysfunction, bone marrow fibrosis, and nephromegaly. To identify the causative gene, we performed homozygosity mapping using single-nucleotide polymorphism arrays, whole-exome sequencing, immunoblotting, immunofluorescence, electron microscopy, a real-time quantitative polymerase-chain-reaction assay, immunohistochemistry, flow cytometry, fibroblast motility assays, measurements of apoptosis, and zebrafish models. Correction experiments were performed by transfecting mutant fibroblasts with the nonmutated gene. RESULTS: All seven affected children had homozygous mutations (Thr224Asn or Glu238Lys, depending on the child's ethnic origin) in VPS45, which encodes a protein that regulates membrane trafficking through the endosomal system. The level of VPS45 protein was reduced, as were the VPS45 binding partners rabenosyn-5 and syntaxin-16. The level of ß1 integrin was reduced on the surface of VPS45-deficient neutrophils and fibroblasts. VPS45-deficient fibroblasts were characterized by impaired motility and increased apoptosis. A zebrafish model of vps45 deficiency showed a marked paucity of myeloperoxidase-positive cells (i.e., neutrophils). Transfection of patient cells with nonmutated VPS45 corrected the migration defect and decreased apoptosis. CONCLUSIONS: Defective endosomal intracellular protein trafficking due to biallelic mutations in VPS45 underlies a new immunodeficiency syndrome involving impaired neutrophil function. (Funded by the National Human Genome Research Institute and others.).

Fluorescence lifetime imaging of DAPI-stained nuclei as a novel diagnostic tool for the detection and classification of B-cell chronic lymphocytic leukemia.

B-cell chronic lymphocytic leukaemia (B-CLL) and B-cell precursor acute lymphoblastic leukaemia (B-ALL) are the most common type of leukaemia in adults and children, respectively. Today, fluorescence in situ hybridization (FISH) is the standard for detecting chromosomal aberrations that reflect adverse and favorable outcome. This study revealed a new, simple, and fast diagnostic tool to detect pathological cells by measuring and imaging the fluorescence lifetime (FLT) using FLT imaging microscopy (FLIM) of the peripheral blood (PB) cells of B-CLL samples that were labeled with the DNA binder, DAPI. The FLT of DAPI in healthy individuals was found to be 2.66 ± 0.12 ns. In contrast, PB cells of B-CLL and BM cells of B-ALL patients were characterized by a specific group distribution of the FLT values. The FLT of DAPI was divided into four subgroups, relative to 2.66 ns: short+, normal, prolonged, and prolonged+. These alterations could be related to different chromatin arrangements of B-CLL and B-ALL interphase nuclei. Notably, extremely long FLT of nuclear DAPI correlate with the presence of extra chromosome 12, while moderate increases compared to normal characterize the deletion of p53. Such correlations potentially enable a FLT-based rapid automatic diagnosis and classification of B-CLL even when the frequency of genetic and chromosomal abnormalities is low. © 2016 International Society for Advancement of Cytometry.

A module of negative feedback regulators defines growth factor signaling.

Signaling pathways invoke interplays between forward signaling and feedback to drive robust cellular response. In this study, we address the dynamics of growth factor signaling through profiling of protein phosphorylation and gene expression, demonstrating the presence of a kinetically defined cluster of delayed early genes that function to attenuate the early events of growth factor signaling. Using epidermal growth factor receptor signaling as the major model system and concentrating on regulation of transcription and mRNA stability, we demonstrate that a number of genes within the delayed early gene cluster function as feedback regulators of immediate early genes. Consistent with their role in negative regulation of cell signaling, genes within this cluster are downregulated in diverse tumor types, in correlation with clinical outcome. More generally, our study proposes a mechanistic description of the cellular response to growth factors by defining architectural motifs that underlie the function of signaling networks.

Sensitive Replicate Real-Time Quantitative PCR of BCR-ABL Shows Deep Molecular Responses in Long-Term Post-Allogeneic Stem Cell Transplantation Chronic Myeloid Leukemia Patients.

Real-time quantitative PCR (RT-qPCR) is commonly used for follow-up of chronic myeloid leukemia (CML) patients treated with tyrosine kinase inhibitors, but its current sensitivity does not allow detection of very low BCR-ABL levels. Therefore RT-qPCR negativity is not synonymous with complete molecular response. Replicate RT-qPCR had shown increased sensitivity in tyrosine kinase inhibitor-treated patients and was, therefore, used here to evaluate whether RT-qPCR-negative post-allogeneic stem cell transplantation (SCT) patients harbor detectable disease. Samples from 12 patients were tested at 2 time points using 82 replicates of BCR-ABL RT-qPCR. One patient (38 months after SCT) had detectable transcripts at baseline and none at the follow-up test, done at a median of 107 months after SCT. This suggests cure from CML in the majority of allogeneic SCT patients who have no transcripts detectable by replicate RT-qPCR for BCR-ABL.