Targeting NHE6 gene expression identifies lysosome and neurodevelopmental mechanisms in a haploid in vitro cell model.

Christianson syndrome (CS) is an X-linked disorder resulting from loss-of-function (LoF) mutations in SLC9A6 encoding the endosomal Na+/H+ exchanger 6 (NHE6). CS presents with developmental delay, seizures, intellectual disability, nonverbal status, postnatal microcephaly, and ataxia. To define transcriptome signatures of NHE6 LoF, we conducted in-depth RNA-sequencing (RNA-seq) analysis on a haploid NHE6 null cell model. CRIPSR/Cas9 genome editing introduced multiple LoF mutations into SLC9A6 in the near haploid human cell line Hap1. Isogenic, paired parental controls were also studied. NHE6 mutant cell lines were confirmed to have intra-endosomal over-acidification as was seen in other NHE6 null cells. RNA-seq analysis was performed by two widely used pipelines: HISAT2-StringTie-DEseq2 and STAR-HTseq-DEseq2. We identified 1056 differentially expressed genes in mutant NHE6 lines, including genes associated with neurodevelopment, synapse function, voltage-dependent calcium channels, and neuronal signaling. Weighted gene co-expression network analysis was then applied and identified a critical module enriched for genes governing lysosome function. By identifying significantly changed gene expression that is associated with lysosomal mechanisms in NHE6-null cells, our analyses suggest that loss of NHE6 function may converge on mechanisms implicated in lysosome-related neurologic disease. Further, this haploid cell model will serve as an important tool for translational science in CS.

Salvage of ribose from uridine or RNA supports glycolysis in nutrient-limited conditions.

Glucose is vital for life, serving as both a source of energy and carbon building block for growth. When glucose is limiting, alternative nutrients must be harnessed. To identify mechanisms by which cells can tolerate complete loss of glucose, we performed nutrient-sensitized genome-wide genetic screens and a PRISM growth assay across 482 cancer cell lines. We report that catabolism of uridine from the medium enables the growth of cells in the complete absence of glucose. While previous studies have shown that uridine can be salvaged to support pyrimidine synthesis in the setting of mitochondrial oxidative phosphorylation deficiency1, our work demonstrates that the ribose moiety of uridine or RNA can be salvaged to fulfil energy requirements via a pathway based on: (1) the phosphorylytic cleavage of uridine by uridine phosphorylase UPP1/UPP2 into uracil and ribose-1-phosphate (R1P), (2) the conversion of uridine-derived R1P into fructose-6-P and glyceraldehyde-3-P by the non-oxidative branch of the pentose phosphate pathway and (3) their glycolytic utilization to fuel ATP production, biosynthesis and gluconeogenesis. Capacity for glycolysis from uridine-derived ribose appears widespread, and we confirm its activity in cancer lineages, primary macrophages and mice in vivo. An interesting property of this pathway is that R1P enters downstream of the initial, highly regulated steps of glucose transport and upper glycolysis. We anticipate that 'uridine bypass' of upper glycolysis could be important in the context of disease and even exploited for therapeutic purposes.

Copper induces cell death by targeting lipoylated TCA cycle proteins.

Copper is an essential cofactor for all organisms, and yet it becomes toxic if concentrations exceed a threshold maintained by evolutionarily conserved homeostatic mechanisms. How excess copper induces cell death, however, is unknown. Here, we show in human cells that copper-dependent, regulated cell death is distinct from known death mechanisms and is dependent on mitochondrial respiration. We show that copper-dependent death occurs by means of direct binding of copper to lipoylated components of the tricarboxylic acid (TCA) cycle. This results in lipoylated protein aggregation and subsequent iron-sulfur cluster protein loss, which leads to proteotoxic stress and ultimately cell death. These findings may explain the need for ancient copper homeostatic mechanisms.

Methionine synthase is essential for cancer cell proliferation in physiological folate environments.

Folate metabolism can be an effective target for cancer treatment. However, standard cell culture conditions utilize folic acid, a non-physiological folate source for most tissues. We find that the enzyme that couples folate and methionine metabolic cycles, methionine synthase, is required for cancer cell proliferation and tumour growth when 5-methyl tetrahydrofolate (THF), the major folate found in circulation, is the extracellular folate source. In such physiological conditions, methionine synthase incorporates 5-methyl THF into the folate cycle to maintain intracellular levels of the folates needed for nucleotide production. 5-methyl THF can sustain intracellular folate metabolism in the absence of folic acid. Therefore, cells exposed to 5-methyl THF are more resistant to methotrexate, an antifolate drug that specifically blocks folic acid incorporation into the folate cycle. Together, these data argue that the environmental folate source has a profound effect on folate metabolism, determining how both folate cycle enzymes and antifolate drugs impact proliferation.

Functional Assessment In Vivo of the Mouse Homolog of the Human Ala-9-Ser NHE6 Variant.

Christianson syndrome (CS) is an X-linked neurogenetic disorder resulting from loss-of-function (LoF) mutations in SLC9A6, which encodes the endosomal Na+/H+ exchanger 6 (NHE6). NHE6 regulates proton efflux from endosomes and, thus, participates in regulating cargo processing and trafficking. LoF mutations in NHE6 cause aberrant acidification of endosomes. While CS arises in males generally due to clear LoF mutations, other potentially hypomorphic variants have emerged, yet most of these variants have not been evaluated for functional effects, particularly in vivo Here we characterize an SLC9A6 variant that has been previously reported in patients, yet now also appears in exome datasets of largely control individuals-c.25G>T, p.A9S. By heterologous expression in cell lines, we show that human NHE6A9S is expressed and localizes in a manner comparable to control NHE6. By genome editing, we generated the equivalent NHE6 mutation in mouse-p.A11S-and determined that male NHE6A11S mice have normal brain size at 6 months of age and do not show cerebellar degeneration or defective neuronal arborization. Neurons from male NHE6A11S mice also did not demonstrate an abnormality in intraendosomal pH compared with controls. These findings are in contrast to findings in NHE6-null mice previously reported and indicate that the NHE6A11S variant functions at a level equivalent to control NHE6 for many of the assays performed. These data stand in support of the population genetic data, which are also evaluated here, indicating that the A9S variant is unlikely to confer disease susceptibility with high penetrance.

GPT2 mutations in autosomal recessive developmental disability: extending the clinical phenotype and population prevalence estimates.

The glutamate pyruvate transaminase 2 (GPT2) gene produces a nuclear-encoded mitochondrial enzyme that catalyzes the reversible transfer of an amino group from glutamate to pyruvate, generating alanine and alpha-ketoglutarate. Recessive mutations in GPT2 have been recently identified in a new syndrome involving intellectual and developmental disability (IDD), postnatal microcephaly, and spastic paraplegia. We have identified additional families with recessive GPT2 mutations and expanded the phenotype to include small stature. GPT2 loss-of-function mutations were identified in four families, nine patients total, including: a homozygous mutation in one child [c.775T>C (p.C259R)]; compound heterozygous mutations in two siblings [c.812A>C (p.N271T)/c.1432_1433delGT (p.V478Rfs*73)]; a novel homozygous, putative splicing mutation [c.1035C>T (p.G345=)]; and finally, a recurrent mutation, previously identified in a distinct family [c.1210C>T (p.R404*)]. All patients were diagnosed with IDD. A majority of patients had remarkably small stature throughout development, many < 1st percentile for height and weight. Given the potential biological function of GPT2 in cellular growth, this phenotype is strongly suggestive of a newly identified clinical susceptibility. Further, homozygous GPT2 mutations manifested in at least 2 of 176 families with IDD (approximately 1.1%) in a Pakistani cohort, thereby representing a relatively common cause of recessive IDD in this population, with recurrence of the p.R404* mutation in this population. Based on variants in the ExAC database, we estimated that approximately 1 in 248 individuals are carriers of moderately or severely deleterious variants in GPT2.

Novel metrics for quantifying bacterial genome composition skews.

BACKGROUND: Bacterial genomes have characteristic compositional skews, which are differences in nucleotide frequency between the leading and lagging DNA strands across a segment of a genome. It is thought that these strand asymmetries arise as a result of mutational biases and selective constraints, particularly for energy efficiency. Analysis of compositional skews in a diverse set of bacteria provides a comparative context in which mutational and selective environmental constraints can be studied. These analyses typically require finished and well-annotated genomic sequences. RESULTS: We present three novel metrics for examining genome composition skews; all three metrics can be computed for unfinished or partially-annotated genomes. The first two metrics, (dot-skew and cross-skew) depend on sequence and gene annotation of a single genome, while the third metric (residual skew) highlights unusual genomes by subtracting a GC content-based model of a library of genome sequences. We applied these metrics to 7738 available bacterial genomes, including partial drafts, and identified outlier species. A phylogenetically diverse set of these outliers (i.e., Borrelia, Ehrlichia, Kinetoplastibacterium, and Phytoplasma) display similar skew patterns but share lifestyle characteristics, such as intracellularity and biosynthetic dependence on their hosts. CONCLUSIONS: Our novel metrics appear to reflect the effects of biosynthetic constraints and adaptations to life within one or more hosts on genome composition. We provide results for each analyzed genome, software and interactive visualizations at http://db.systemsbiology.net/gestalt/ skew_metrics .

Whole genome sequence and comparative analysis of Borrelia burgdorferi MM1.

Lyme disease is caused by spirochaetes of the Borrelia burgdorferi sensu lato genospecies. Complete genome assemblies are available for fewer than ten strains of Borrelia burgdorferi sensu stricto, the primary cause of Lyme disease in North America. MM1 is a sensu stricto strain originally isolated in the midwestern United States. Aside from a small number of genes, the complete genome sequence of this strain has not been reported. Here we present the complete genome sequence of MM1 in relation to other sensu stricto strains and in terms of its Multi Locus Sequence Typing. Our results indicate that MM1 is a new sequence type which contains a conserved main chromosome and 15 plasmids. Our results include the first contiguous 28.5 kb assembly of lp28-8, a linear plasmid carrying the vls antigenic variation system, from a Borrelia burgdorferi sensu stricto strain.

Differences between the genomes of lymphoblastoid cell lines and blood-derived samples.

Lymphoblastoid cell lines (LCLs) represent a convenient research tool for expanding the amount of biologic material available from an individual. LCLs are commonly used as reference materials, most notably from the Genome in a Bottle Consortium. However, the question remains how faithfully LCL-derived genome assemblies represent the germline genome of the donor individual as compared to the genome assemblies derived from peripheral blood mononuclear cells. We present an in-depth comparison of a large collection of LCL- and peripheral blood mononuclear cell-derived genomes in terms of distributions of coverage and copy number alterations. We found significant differences in the depth of coverage and copy number calls, which may be driven by differential replication timing. Importantly, these copy number changes preferentially affect regions closer to genes and with higher GC content. This suggests that genomic studies based on LCLs may display locus-specific biases, and that conclusions based on analysis of depth of coverage and copy number variation may require further scrutiny.