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.

Mutations in PPCS, Encoding Phosphopantothenoylcysteine Synthetase, Cause Autosomal-Recessive Dilated Cardiomyopathy.

Coenzyme A (CoA) is an essential metabolic cofactor used by around 4% of cellular enzymes. Its role is to carry and transfer acetyl and acyl groups to other molecules. Cells can synthesize CoA de novo from vitamin B5 (pantothenate) through five consecutive enzymatic steps. Phosphopantothenoylcysteine synthetase (PPCS) catalyzes the second step of the pathway during which phosphopantothenate reacts with ATP and cysteine to form phosphopantothenoylcysteine. Inborn errors of CoA biosynthesis have been implicated in neurodegeneration with brain iron accumulation (NBIA), a group of rare neurological disorders characterized by accumulation of iron in the basal ganglia and progressive neurodegeneration. Exome sequencing in five individuals from two unrelated families presenting with dilated cardiomyopathy revealed biallelic mutations in PPCS, linking CoA synthesis with a cardiac phenotype. Studies in yeast and fruit flies confirmed the pathogenicity of identified mutations. Biochemical analysis revealed a decrease in CoA levels in fibroblasts of all affected individuals. CoA biosynthesis can occur with pantethine as a source independent from PPCS, suggesting pantethine as targeted treatment for the affected individuals still alive.

Correction to: ESR1 mutations are frequent in newly diagnosed metastatic and loco-regional recurrence of endocrine-treated breast cancer and carry worse prognosis.

After the publication of the original article [1], we were notified the upper panel of the Fig. 1, where the patients' codes are listed, was cropped by mistake so the patients 1-8 are repeated.

ESR1 mutations are frequent in newly diagnosed metastatic and loco-regional recurrence of endocrine-treated breast cancer and carry worse prognosis.

BACKGROUND: Emerging mutations in the ESR1 gene that encodes for the estrogen receptor (ER) are associated with resistance to endocrine therapy. ESR1 mutations rarely exist in primary tumors (~ 1%) but are relatively common (10-50%) in metastatic, endocrine therapy-resistant cancers and are associated with a shorter progression-free survival. Little is known about the incidence and clinical implication of these mutations in early recurrence events, such as local recurrences or newly diagnosed metastatic disease. METHODS: We collected 130 archival tumor samples from 103 breast cancer patients treated with endocrine therapy prior to their local/metastatic recurrence. The cohort consisted of 41 patients having at least 1 sample from local/loco-regional recurrence and 62 patients with metastatic disease (of whom 41 newly diagnosed and 28 with advanced disease). The 5 most common ESR1 hotspot mutations (D538G, L536R, Y537S/N/C) were analyzed either by targeted sequencing or by droplet digital PCR. Progression-free survival (PFS), disease-free survival (DFS), and distant recurrence-free survival (DRFS) were statistically tested by Kaplan-Meier analysis. RESULTS: The prevalence of ESR1 mutations was 5/41 (12%) in newly diagnosed metastatic patients and 5/28 (18%) for advanced metastases, detected at allele frequency > 1%. All mutations in advanced metastases were detected in patients previously treated with both tamoxifen (TAM) and aromatase inhibitors (AI). However, in newly diagnosed metastatic patients, 4/5 mutations occurred in patients treated with TAM alone. PFS on AI treatment in metastatic patients was significantly shorter for ESR1 mutation carriers (p = 0.017). In the local recurrence cohort, ESR1 mutations were identified in 15/41 (36%) patients but only 4/41 (10%) were detected at allele frequency > 1%. Again, most mutations (3/4) were detected under TAM monotherapy. Notably, 1 patient developed ESR1 mutation while on neoadjuvant endocrine therapy. DFS and DRFS were significantly shorter (p = 0.04 and p = 0.017, respectively) in patients that had ESR1 mutations (> 1%) in their loco-regional recurrence tumor. CONCLUSIONS: Clinically relevant ESR1 mutations are prevalent in newly diagnosed metastatic and local recurrence of endocrine-treated breast cancer. Since local recurrences are amenable to curative therapy, these mutations may inform the selection of subsequent endocrine therapies.

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.

MECR Mutations Cause Childhood-Onset Dystonia and Optic Atrophy, a Mitochondrial Fatty Acid Synthesis Disorder.

Mitochondrial fatty acid synthesis (mtFAS) is an evolutionarily conserved pathway essential for the function of the respiratory chain and several mitochondrial enzyme complexes. We report here a unique neurometabolic human disorder caused by defective mtFAS. Seven individuals from five unrelated families presented with childhood-onset dystonia, optic atrophy, and basal ganglia signal abnormalities on MRI. All affected individuals were found to harbor recessive mutations in MECR encoding the mitochondrial trans-2-enoyl-coenzyme A-reductase involved in human mtFAS. All six mutations are extremely rare in the general population, segregate with the disease in the families, and are predicted to be deleterious. The nonsense c.855T>G (p.Tyr285∗), c.247_250del (p.Asn83Hisfs∗4), and splice site c.830+2_830+3insT mutations lead to C-terminal truncation variants of MECR. The missense c.695G>A (p.Gly232Glu), c.854A>G (p.Tyr285Cys), and c.772C>T (p.Arg258Trp) mutations involve conserved amino acid residues, are located within the cofactor binding domain, and are predicted by structural analysis to have a destabilizing effect. Yeast modeling and complementation studies validated the pathogenicity of the MECR mutations. Fibroblast cell lines from affected individuals displayed reduced levels of both MECR and lipoylated proteins as well as defective respiration. These results suggest that mutations in MECR cause a distinct human disorder of the mtFAS pathway. The observation of decreased lipoylation raises the possibility of a potential therapeutic strategy.

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.

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.

Congenital Sucrase-isomaltase Deficiency: A Novel Compound Heterozygous Mutation Causing Aberrant Protein Localization.

OBJECTIVES: Congenital diarrheal disorders is a group of inherited enteropathies presenting in early life and requiring parenteral nutrition. In most cases, genetics may be the key for precise diagnosis. We present an infant girl with chronic congenital diarrhea that resolved after introduction of fructose-based formula but had no identified mutation in the SLC5A1 gene. Using whole exome sequencing (WES) we identified other mutations that better dictated dietary adjustments. METHODS: WES of the patient and her parents was performed. The analysis focused on recessive model including compound heterozygous mutations. Sanger sequencing was used to validate identified mutations and to screen the patient's newborn sister and grandparents. Expression and localization analysis were performed in the patient's duodenal biopsies using immunohistochemistry. RESULTS: Using WES we identified a new compound heterozygote mutation in sucrase-isomaltase (SI) gene; a maternal inherited known V577G mutation, and a novel paternal inherited C1531W mutation. Importantly, the newborn offspring carried similar compound heterozygous mutations. Computational predictions suggest that both mutations highly destabilize the protein. SI expression and localization studies determined that the mutated SI protein was not expressed on the brush border membrane in the patient's duodenal biopsies, verifying the diagnosis of congenital sucrase-isomaltase deficiency (CSID). CONCLUSIONS: The novel compound heterozygote V577G/C1531W SI mutations lead to lack of SI expression in the duodenal brush border, confirming the diagnosis of CSID. These cases of CSID extend the molecular spectrum of this condition, further directing a more adequate dietary intervention for the patient and newborn sibling.

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 Novel Mutation in a Critical Region for the Methyl Donor Binding in DNMT3B Causes Immunodeficiency, Centromeric Instability, and Facial Anomalies Syndrome (ICF).

PURPOSE: Immunodeficiency, centromeric instability, and facial anomalies (ICF) syndrome is an extremely rare autosomal recessive disease. The immune phenotype is characterized by hypogammaglobulinemia in the presence of B cells. T cell lymphopenia also develops in some patients. We sought to further investigate the immune defect in an ICF patient with a novel missense mutation in DNMT3B and a severe phenotype. METHODS: Patient lymphocytes were examined for subset counts, immunoglobulin levels, T and B cell de novo production (via excision circles) and receptor repertoire diversity. Mutated DNMT3B protein structure was modeled to assess the effect of a mutation located outside of the catalytic region on protein function. RESULTS: A novel homozygous missense mutation, Ala585Thr, was found in DNMT3B. The patient had decreased B cell counts with hypogammaglobulinemia, and normal T cell counts. CD4+ T cells decreased over time, leading to an inversion of the CD4+ to CD8+ ratio. Excision circle copy numbers were normal, signifying normal de novo lymphocyte production, but the ratio between naïve and total B cells was low, indicating decreased in vivo B cell replication. T and B cell receptor repertoires displayed normal diversity. Computerized modeling of the mutated Ala585 residue suggested reduced thermostability, possibly affecting the enzyme kinetics. CONCLUSIONS: Our results highlight the existence of a T cell defect that develops over time in ICF patient, in addition to the known B cell dysfunction. With intravenous immunoglobulin (IVIG) treatment ameliorating the B cell defect, the extent of CD4+ lymphopenia may determine the severity of ICF immunodeficiency.

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.).

The anisotropic network model web server at 2015 (ANM 2.0).

SUMMARY: The anisotropic network model (ANM) is one of the simplest yet powerful tools for exploring protein dynamics. Its main utility is to predict and visualize the collective motions of large complexes and assemblies near their equilibrium structures. The ANM server, introduced by us in 2006 helped making this tool more accessible to non-sophisticated users. We now provide a new version (ANM 2.0), which allows inclusion of nucleic acids and ligands in the network model and thus enables the investigation of the collective motions of protein-DNA/RNA and -ligand systems. The new version offers the flexibility of defining the system nodes and the interaction types and cutoffs. It also includes extensive improvements in hardware, software and graphical interfaces. AVAILABILITY AND IMPLEMENTATION: ANM 2.0 is available at http://anm.csb.pitt.edu CONTACT: eran.eyal@sheba.health.gov.il, eyal.eran@gmail.com.

Expanding the molecular diversity and phenotypic spectrum of glycerol 3-phosphate dehydrogenase 1 deficiency.

Transient infantile hypertriglyceridemia (HTGT1; OMIM #614480) is a rare autosomal recessive disorder, which manifests in early infancy with transient hypertriglyceridemia, hepatomegaly, elevated liver enzymes, persistent fatty liver and hepatic fibrosis. This rare clinical entity is caused by inactivating mutations in the GPD1 gene, which encodes the cytosolic isoform of glycerol-3-phosphate dehydrogenase. Here we report on four patients from three unrelated families of diverse ethnic origins, who presented with hepatomegaly, liver steatosis, hypertriglyceridemia, with or without fasting ketotic hypoglycemia. Whole exome sequencing revealed the affected individuals to harbor deleterious biallelic mutations in the GPD1 gene, including the previously undescribed c.806G > A (p.Arg269Gln) and c.640T > C (p.Cys214Arg) mutations. The clinical features in three of our patients showed several differences compared to the original reports. One subject presented with recurrent episodes of fasting hypoglycemia along with hepatomegaly, hypetriglyceridemia, and elevated liver enzymes; the second showed a severe liver disease, with intrahepatic cholestasis associated with kidney involvement; finally, the third presented persistent hypertriglyceridemia at the age of 30 years. These findings expand the current knowledge of this rare disorder, both with regard to the phenotype and molecular basis. The enlarged phenotypic spectrum of glycerol-3-phosphate dehydrogenase 1 deficiency can mimic other inborn errors of metabolism with liver involvement and should alert clinicians to recognize this entity by considering GPD1 mutations in appropriate clinical settings.

Loss of function of NaPiIIa causes nephrocalcinosis and possibly kidney insufficiency.

BACKGROUND: Inherited metabolic disorders associated with nephrocalcinosis are rare conditions. The aim of this study was to identify the genetic cause of an Israeli-Arab boy from a consanguineous family with severe nephrocalcinosis and kidney insufficiency. METHODS: Clinical and biochemical data of the proband and family members were obtained from both previous and recent medical charts. Genomic DNA was isolated from peripheral blood cells. The coding sequence and splice sites of candidate genes (CYP24A1, CYP27B1, FGF23, KLOTHO, SLC34A3 and SLC34A1) were sequenced directly. Functional studies were performed in Xenopus laevis oocytes and in transfected opossum kidney (OK) cells. RESULTS: Our patient was identified as having nephrocalcinosis in utero, and at the age of 16.5 years, he had kidney insufficiency but no bone disease. Genetic analysis revealed a novel homozygous missense mutation, Arg215Gln, in SLC34A1, which encodes the renal sodium phosphate cotransporter NaPiIIa. Functional studies of the Arg215Gln mutant revealed reduced transport activity in Xenopus laevis oocytes and increased intracellular cytoplasmic accumulation in OK cells. CONCLUSIONS: Our findings show that dysfunction of the human NaPiIIa causes severe renal calcification that may eventually lead to reduced kidney function, rather than complications of phosphate loss.

Global regulation of alternative splicing by adenosine deaminase acting on RNA (ADAR).

Alternative mRNA splicing is a major mechanism for gene regulation and transcriptome diversity. Despite the extent of the phenomenon, the regulation and specificity of the splicing machinery are only partially understood. Adenosine-to-inosine (A-to-I) RNA editing of pre-mRNA by ADAR enzymes has been linked to splicing regulation in several cases. Here we used bioinformatics approaches, RNA-seq and exon-specific microarray of ADAR knockdown cells to globally examine how ADAR and its A-to-I RNA editing activity influence alternative mRNA splicing. Although A-to-I RNA editing only rarely targets canonical splicing acceptor, donor, and branch sites, it was found to affect splicing regulatory elements (SREs) within exons. Cassette exons were found to be significantly enriched with A-to-I RNA editing sites compared with constitutive exons. RNA-seq and exon-specific microarray revealed that ADAR knockdown in hepatocarcinoma and myelogenous leukemia cell lines leads to global changes in gene expression, with hundreds of genes changing their splicing patterns in both cell lines. This global change in splicing pattern cannot be explained by putative editing sites alone. Genes showing significant changes in their splicing pattern are frequently involved in RNA processing and splicing activity. Analysis of recently published RNA-seq data from glioblastoma cell lines showed similar results. Our global analysis reveals that ADAR plays a major role in splicing regulation. Although direct editing of the splicing motifs does occur, we suggest it is not likely to be the primary mechanism for ADAR-mediated regulation of alternative splicing. Rather, this regulation is achieved by modulating trans-acting factors involved in the splicing machinery.

Characterizing of functional human coding RNA editing from evolutionary, structural, and dynamic perspectives.

A-to-I RNA editing has been recently shown to be a widespread phenomenon with millions of sites spread in the human transcriptome. However, only few are known to be located in coding sequences and modify the amino acid sequence of the protein product. Here, we used high-throughput data, variant prediction tools, and protein structural information in order to find structural and functional preferences for coding RNA editing. We show that RNA editing has a unique pattern of amino acid changes characterized by enriched stop-to-tryptophan changes, positive-to-neutral and neutral-to-positive charge changes. RNA editing tends to have stronger structural effect than equivalent A-to-G SNPs but weaker effect than random A-to-G mutagenesis events. Sites edited at low level tend to be located at conserved positions with stronger predicted deleterious effect on proteins comparing to sites edited at high frequencies. Lowly edited sites tend to destabilize the protein structure and affect amino acids with larger number of intra-molecular contacts. Still, some highly edited sites are predicted also to prominently affect structure and tend to be located at critical positions of the protein matrix and are likely to be functionally important. Using our pipeline, we identify and discuss several novel putative functional coding changing editing sites in the genes COPA (I164V), GIPC1 (T62A), ZN358 (K382R), and CCNI (R75G).

Detection of BCR-ABL1 mutations in chronic myeloid leukaemia by massive parallel sequencing.

The ability to sequence nucleic acids at an unprecedented pace and decreased costs using massive parallel sequencing (MPS) strongly affects biomedical research. Here we applied MPS for the detection of rare, clinically relevant mutations in a chronic myeloid leukaemia (CML) patient. Tyrosine kinase inhibitors revolutionized CML therapy but in some patients the disease progresses due to resistance-conferring mutations. MPS was applied herein to monitor such mutations in BCR-ABL1 transcripts at different time points. The large volume of sequencing data increases sensitivity compared to direct sequencing and allows detection of marginally represented and previously uncharacterized mutations. We detected changes in the frequency of mutated clones including the emergence and disappearance of the resistance-associated ABL1 T315I mutation. We also observed correlation in appearance of adjacent mutations, and exploited this observation to demonstrate the existence of mutated clones at the time of diagnosis. A tool is provided for detection of low frequency single nucleotide variants/mutations from deep coverage MPS data, applicable to clinical translation of advanced sequencing technologies.

A comparative analysis of total serum miRNA profiles identifies novel signature that is highly indicative of metastatic melanoma: a pilot study.

CONTEXT: Quantification of circulating microRNAs (miRNAs) has recently become feasible and reliable, with most efforts focusing on miRNAs overexpressed by cancer cells. OBJECTIVE: Identification of a characteristic circulating miRNAs profile in melanoma patients. METHODS: We conducted a pilot study comprised of unbiased qPCR comparison of serum miRNA profiles between metastatic melanoma patients and healthy donors. RESULTS: Loss of two normal serum-miRNAs, miR-29c and miR-324-3p, is highly indicative of metastatic melanoma. Hierarchical clustering analysis supported the results and clearly distinguished melanoma patients from healthy donors, metastatic colon and renal cancer patients. DISCUSSION AND CONCLUSIONS: This approach is independent of tumor heterogeneity and is expected to have superior biomarker performances.