mRNA stability and m6A are major determinants of subcellular mRNA localization in neurons.

For cells to perform their biological functions, they need to adopt specific shapes and form functionally distinct subcellular compartments. This is achieved in part via an asymmetric distribution of mRNAs within cells. Currently, the main model of mRNA localization involves specific sequences called "zipcodes" that direct mRNAs to their proper locations. However, while thousands of mRNAs localize within cells, only a few zipcodes have been identified, suggesting that additional mechanisms contribute to localization. Here, we assess the role of mRNA stability in localization by combining the isolation of the soma and neurites of mouse primary cortical and mESC-derived neurons, SLAM-seq, m6A-RIP-seq, the perturbation of mRNA destabilization mechanisms, and the analysis of multiple mRNA localization datasets. We show that depletion of mRNA destabilization elements, such as m6A, AU-rich elements, and suboptimal codons, functions as a mechanism that mediates the localization of mRNAs associated with housekeeping functions to neurites in several types of neurons.

Mime-seq 2.0: a method to sequence microRNAs from specific mouse cell types.

Many microRNAs (miRNAs) are expressed with high spatiotemporal specificity during organismal development, with some being limited to rare cell types, often embedded in complex tissues. Yet, most miRNA profiling efforts remain at the tissue and organ levels. To overcome challenges in accessing the microRNomes from tissue-embedded cells, we had previously developed mime-seq (miRNome by methylation-dependent sequencing), a technique in which cell-specific miRNA methylation in C. elegans and Drosophila enabled chemo-selective sequencing without the need for cell sorting or biochemical purification. Here, we present mime-seq 2.0 for profiling miRNAs from specific mouse cell types. We engineered a chimeric RNA methyltransferase that is tethered to Argonaute protein and efficiently methylates miRNAs at their 3'-terminal 2'-OH in mouse and human cell lines. We also generated a transgenic mouse for conditional expression of this methyltransferase, which can be used to direct methylation of miRNAs in a cell type of choice. We validated the use of this mouse model by profiling miRNAs from B cells and bone marrow plasma cells.

Natural RNA Polymerase Aptamers Regulate Transcription in E. coli.

In search for RNA signals that modulate transcription via direct interaction with RNA polymerase (RNAP), we deep sequenced an E. coli genomic library enriched for RNAP-binding RNAs. Many natural RNAP-binding aptamers, termed RAPs, were mapped to the genome. Over 60% of E. coli genes carry RAPs in their mRNA. Combining in vitro and in vivo approaches, we characterized a subset of inhibitory RAPs (iRAPs) that promote Rho-dependent transcription termination. A representative iRAP within the coding region of the essential gene, nadD, greatly reduces its transcriptional output in stationary phase and under oxidative stress, demonstrating that iRAPs control gene expression in response to changing environment. The mechanism of iRAPs involves active uncoupling of transcription and translation, making nascent RNA accessible to Rho. iRAPs encoded in the antisense strand also promote gene expression by reducing transcriptional interference. In essence, our work uncovers a broad class of cis-acting RNA signals that globally control bacterial transcription.

GREEN-DB: a framework for the annotation and prioritization of non-coding regulatory variants from whole-genome sequencing data.

Non-coding variants have long been recognized as important contributors to common disease risks, but with the expansion of clinical whole genome sequencing, examples of rare, high-impact non-coding variants are also accumulating. Despite recent advances in the study of regulatory elements and the availability of specialized data collections, the systematic annotation of non-coding variants from genome sequencing remains challenging. Here, we propose a new framework for the prioritization of non-coding regulatory variants that integrates information about regulatory regions with prediction scores and HPO-based prioritization. Firstly, we created a comprehensive collection of annotations for regulatory regions including a database of 2.4 million regulatory elements (GREEN-DB) annotated with controlled gene(s), tissue(s) and associated phenotype(s) where available. Secondly, we calculated a variation constraint metric and showed that constrained regulatory regions associate with disease-associated genes and essential genes from mouse knock-outs. Thirdly, we compared 19 non-coding impact prediction scores providing suggestions for variant prioritization. Finally, we developed a VCF annotation tool (GREEN-VARAN) that can integrate all these elements to annotate variants for their potential regulatory impact. In our evaluation, we show that GREEN-DB can capture previously published disease-associated non-coding variants as well as identify additional candidate disease genes in trio analyses.

Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia.

The transformation of chronic lymphocytic leukemia (CLL) to high-grade B-cell lymphoma is known as Richter syndrome (RS), a rare event with dismal prognosis. In this study, we conducted whole-genome sequencing (WGS) of paired circulating CLL (PB-CLL) and RS biopsies (tissue-RS) from 17 patients recruited into a clinical trial (CHOP-O). We found that tissue-RS was enriched for mutations in poor-risk CLL drivers and genes in the DNA damage response (DDR) pathway. In addition, we identified genomic aberrations not previously implicated in RS, including the protein tyrosine phosphatase receptor (PTPRD) and tumor necrosis factor receptor-associated factor 3 (TRAF3). In the noncoding genome, we discovered activation-induced cytidine deaminase-related and unrelated kataegis in tissue-RS affecting regulatory regions of key immune-regulatory genes. These include BTG2, CXCR4, NFATC1, PAX5, NOTCH-1, SLC44A5, FCRL3, SELL, TNIP2, and TRIM13. Furthermore, differences between the global mutation signatures of pairs of PB-CLL and tissue-RS samples implicate DDR as the dominant mechanism driving transformation. Pathway-based clonal deconvolution analysis showed that genes in the MAPK and DDR pathways demonstrate high clonal-expansion probability. Direct comparison of nodal-CLL and tissue-RS pairs from an independent cohort confirmed differential expression of the same pathways by RNA expression profiling. Our integrated analysis of WGS and RNA expression data significantly extends previous targeted approaches, which were limited by the lack of germline samples, and it facilitates the identification of novel genomic correlates implicated in RS transformation, which could be targeted therapeutically. Our results inform the future selection of investigative agents for a UK clinical platform study. This trial was registered at www.clinicaltrials.gov as #NCT03899337.

An inhibitor of complement C5 provides structural insights into activation.

The complement system is a crucial part of innate immune defenses against invading pathogens. The blood-meal of the tick Rhipicephalus pulchellus lasts for days, and the tick must therefore rely on inhibitors to counter complement activation. We have identified a class of inhibitors from tick saliva, the CirpT family, and generated detailed structural data revealing their mechanism of action. We show direct binding of a CirpT to complement C5 and have determined the structure of the C5-CirpT complex by cryoelectron microscopy. This reveals an interaction with the peripheral macro globulin domain 4 (C5_MG4) of C5. To achieve higher resolution detail, the structure of the C5_MG4-CirpT complex was solved by X-ray crystallography (at 2.7 Å). We thus present the fold of the CirpT protein family, and provide detailed mechanistic insights into its inhibitory function. Analysis of the binding interface reveals a mechanism of C5 inhibition, and provides information to expand our biological understanding of the activation of C5, and thus the terminal complement pathway.

Nanopanel2 calls phased low-frequency variants in Nanopore panel sequencing data.

MOTIVATION: Clinical decision making is increasingly guided by accurate and recurrent determination of presence and frequency of (somatic) variants and their haplotype through panel sequencing of disease-relevant genomic regions. Haplotype calling (phasing), however, is difficult and error prone unless variants are located on the same read which limits the ability of short-read sequencing to detect, e.g. co-occurrence of drug-resistance variants. Long-read panel sequencing enables direct phasing of amplicon variants besides having multiple other benefits, however, high error rates of current technologies prevented their applicability in the past. RESULTS: We have developed Nanopanel2, a variant caller for Nanopore panel sequencing data. Nanopanel2 works directly on base-called FAST5 files and uses allele probability distributions and several other filters to robustly separate true from false positive (FP) calls. It effectively calls SNVs and INDELs with variant allele frequencies as low as 1% and 5%, respectively, and produces only few low-frequency false-positive calls (∼1 FP call with VAF<5% per kb amplicon). Haplotype compositions are then determined by direct phasing. Nanopanel2 is the first somatic variant caller for Nanopore data, enabling accurate, fast (turnaround <48 h) and cheap (sequencing costs ∼10$/sample) diagnostic workflows. AVAILABILITYAND IMPLEMENTATION: The data for this study have been deposited at zenodo.org under DOIs accession numbers 4110691 and 4110698. Nanopanel2 is open source and available at https://github.com/popitsch/nanopanel2. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VARAN-GIE: curation of genomic interval sets.

SUMMARY: Genomic interval sets are fundamental elements of genome annotation and are the output of countless bioinformatics applications. Nevertheless, tool support for the manual curation of these data is currently limited. We developed VARAN-GIE, an extension of the popular Integrative Genomics Viewer (IGV) that adds functionality to edit, annotate and merge genomic interval sets. Data can easily be shared with other users and imported/exported from/to multiple common data formats. AVAILABILITY AND IMPLEMENTATION: VARAN-GIE binary releases, source-code, user guides and tutorials are available at https://github.com/popitsch/varan-gie/.

The complete costs of genome sequencing: a microcosting study in cancer and rare diseases from a single center in the United Kingdom.

PURPOSE: The translation of genome sequencing into routine health care has been slow, partly because of concerns about affordability. The aspirational cost of sequencing a genome is $1000, but there is little evidence to support this estimate. We estimate the cost of using genome sequencing in routine clinical care in patients with cancer or rare diseases. METHODS: We performed a microcosting study of Illumina-based genome sequencing in a UK National Health Service laboratory processing 399 samples/year. Cost data were collected for all steps in the sequencing pathway, including bioinformatics analysis and reporting of results. Sensitivity analysis identified key cost drivers. RESULTS: Genome sequencing costs £6841 per cancer case (comprising matched tumor and germline samples) and £7050 per rare disease case (three samples). The consumables used during sequencing are the most expensive component of testing (68-72% of the total cost). Equipment costs are higher for rare disease cases, whereas consumable and staff costs are slightly higher for cancer cases. CONCLUSION: The cost of genome sequencing is underestimated if only sequencing costs are considered, and likely surpasses $1000/genome in a single laboratory. This aspirational sequencing cost will likely only be achieved if consumable costs are considerably reduced and sequencing is performed at scale.

Shep interacts with posttranscriptional regulators to control dendrite morphogenesis in sensory neurons.

RNA binding proteins (RBPs) mediate posttranscriptional gene regulatory events throughout development. During neurogenesis, many RBPs are required for proper dendrite morphogenesis within Drosophila sensory neurons. Despite their fundamental role in neuronal morphogenesis, little is known about the molecular mechanisms in which most RBPs participate during neurogenesis. In Drosophila, alan shepard (shep) encodes a highly conserved RBP that regulates dendrite morphogenesis in sensory neurons. Moreover, the C. elegans ortholog sup-26 has also been implicated in sensory neuron dendrite morphogenesis. Nonetheless, the molecular mechanism by which Shep/SUP-26 regulate dendrite development is not understood. Here we show that Shep interacts with the RBPs Trailer Hitch (Tral), Ypsilon schachtel (Yps), Belle (Bel), and Poly(A)-Binding Protein (PABP), to direct dendrite morphogenesis in Drosophila sensory neurons. Moreover, we identify a conserved set of Shep/SUP-26 target RNAs that include regulators of cell signaling, posttranscriptional gene regulators, and known regulators of dendrite development.

Impaired plasticity of macrophages in X-linked adrenoleukodystrophy.

X-linked adrenoleukodystrophy is caused by ATP-binding cassette transporter D1 (ABCD1) mutations and manifests by default as slowly progressive spinal cord axonopathy with associated demyelination (adrenomyloneuropathy). In 60% of male cases, however, X-linked adrenoleukodystrophy converts to devastating cerebral inflammation and demyelination (cerebral adrenoleukodystrophy) with infiltrating blood-derived monocytes and macrophages and cytotoxic T cells that can only be stopped by allogeneic haematopoietic stem cell transplantation or gene therapy at an early stage of the disease. Recently, we identified monocytes/macrophages but not T cells to be severely affected metabolically by ABCD1 deficiency. Here we found by whole transcriptome analysis that, although monocytes of patients with X-linked adrenoleukodystrophy have normal capacity for macrophage differentiation and phagocytosis, they are pro-inflammatory skewed also in patients with adrenomyloneuropathy in the absence of cerebral inflammation. Following lipopolysaccharide activation, the ingestion of myelin debris, normally triggering anti-inflammatory polarization, did not fully reverse the pro-inflammatory status of X-linked adrenoleukodystrophy macrophages. Immunohistochemistry on post-mortem cerebral adrenoleukodystrophy lesions reflected the activation pattern by prominent presence of enlarged lipid-laden macrophages strongly positive for the pro-inflammatory marker co-stimulatory molecule CD86. Comparative analyses of lesions with matching macrophage density in cases of cerebral adrenoleukodystrophy and acute multiple sclerosis showed a similar extent of pro-inflammatory activation but a striking reduction of anti-inflammatory mannose receptor (CD206) and haemoglobin-haptoglobin receptor (CD163) expression on cerebral adrenoleukodystrophy macrophages. Accordingly, ABCD1-deficiency leads to an impaired plasticity of macrophages that is reflected in incomplete establishment of anti-inflammatory responses, thus possibly contributing to the devastating rapidly progressive demyelination in cerebral adrenoleukodystrophy that only in rare cases arrests spontaneously. These findings emphasize monocytes/macrophages as crucial therapeutic targets for preventing or stopping myelin destruction in patients with X-linked adrenoleukodystrophy.

In vivo expression technology and 5' end mapping of the Borrelia burgdorferi transcriptome identify novel RNAs expressed during mammalian infection.

Borrelia burgdorferi, the bacterial pathogen responsible for Lyme disease, modulates its gene expression profile in response to the environments encountered throughout its tick-mammal infectious cycle. To begin to characterize the B. burgdorferi transcriptome during murine infection, we previously employed an in vivo expression technology-based approach (BbIVET). This identified 233 putative promoters, many of which mapped to un-annotated regions of the complex, segmented genome. Herein, we globally identify the 5' end transcriptome of B. burgdorferi grown in culture as a means to validate non-ORF associated promoters discovered through BbIVET. We demonstrate that 119 BbIVET promoters are associated with transcription start sites (TSSs) and validate novel RNA transcripts using Northern blots and luciferase promoter fusions. Strikingly, 49% of BbIVET promoters were not found to associate with TSSs. This finding suggests that these sequences may be primarily active in the mammalian host. Furthermore, characterization of the 6042 B. burgdorferi TSSs reveals a variety of RNAs including numerous antisense and intragenic transcripts, leaderless RNAs, long untranslated regions and a unique nucleotide frequency for initiating intragenic transcription. Collectively, this is the first comprehensive map of TSSs in B. burgdorferi and characterization of previously un-annotated RNA transcripts expressed by the spirochete during murine infection.

ReliableGenome: annotation of genomic regions with high/low variant calling concordance.

MOTIVATION: The increasing adoption of clinical whole-genome resequencing (WGS) demands for highly accurate and reproducible variant calling (VC) methods. The observed discordance between state-of-the-art VC pipelines, however, indicates that the current practice still suffers from non-negligible numbers of false positive and negative SNV and INDEL calls that were shown to be enriched among discordant calls but also in genomic regions with low sequence complexity. RESULTS: Here, we describe our method ReliableGenome (RG) for partitioning genomes into high and low concordance regions with respect to a set of surveyed VC pipelines. Our method combines call sets derived by multiple pipelines from arbitrary numbers of datasets and interpolates expected concordance for genomic regions without data. By applying RG to 219 deep human WGS datasets, we demonstrate that VC concordance depends predominantly on genomic context rather than the actual sequencing data which manifests in high recurrence of regions that can/cannot be reliably genotyped by a single method. This enables the application of pre-computed regions to other data created with comparable sequencing technology and software. RG outperforms comparable efforts in predicting VC concordance and false positive calls in low-concordance regions which underlines its usefulness for variant filtering, annotation and prioritization. RG allows focusing resource-intensive algorithms (e.g. consensus calling methods) on the smaller, discordant share of the genome (20-30%) which might result in increased overall accuracy at reasonable costs. Our method and analysis of discordant calls may further be useful for development, benchmarking and optimization of VC algorithms and for the relative comparison of call sets between different studies/pipelines. AVAILABILITY AND IMPLEMENTATION: RG was implemented in Java, source code and binaries are freely available for non-commercial use at https://github.com/popitsch/wtchg-rg/ CONTACT: niko@wtchg.ox.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: pilot study for the 100,000 Genomes Project.

PURPOSE: Fresh-frozen (FF) tissue is the optimal source of DNA for whole-genome sequencing (WGS) of cancer patients. However, it is not always available, limiting the widespread application of WGS in clinical practice. We explored the viability of using formalin-fixed, paraffin-embedded (FFPE) tissues, available routinely for cancer patients, as a source of DNA for clinical WGS. METHODS: We conducted a prospective study using DNAs from matched FF, FFPE, and peripheral blood germ-line specimens collected from 52 cancer patients (156 samples) following routine diagnostic protocols. We compared somatic variants detected in FFPE and matching FF samples. RESULTS: We found the single-nucleotide variant agreement reached 71% across the genome and somatic copy-number alterations (CNAs) detection from FFPE samples was suboptimal (0.44 median correlation with FF) due to nonuniform coverage. CNA detection was improved significantly with lower reverse crosslinking temperature in FFPE DNA extraction (80 °C or 65 °C depending on the methods). Our final data showed somatic variant detection from FFPE for clinical decision making is possible. We detected 98% of clinically actionable variants (including 30/31 CNAs). CONCLUSION: We present the first prospective WGS study of cancer patients using FFPE specimens collected in a routine clinical environment proving WGS can be applied in the clinic.

Temperature-dependent sRNA transcriptome of the Lyme disease spirochete.

BACKGROUND: Transmission of Borrelia burgdorferi from its tick vector to a vertebrate host requires extensive reprogramming of gene expression. Small regulatory RNAs (sRNA) have emerged in the last decade as important regulators of bacterial gene expression. Despite the widespread observation of sRNA-mediated gene regulation, only one sRNA has been characterized in the Lyme disease spirochete B. burgdorferi. We employed an sRNA-specific deep-sequencing approach to identify the small RNA transcriptome of B. burgdorferi at both 23 °C and 37 °C, which mimics in vitro the transmission from the tick vector to the mammalian host. RESULTS: We identified over 1000 sRNAs in B. burgdorferi revealing large amounts of antisense and intragenic sRNAs, as well as characteristic intergenic and 5' UTR-associated sRNAs. A large fraction of the novel sRNAs (43%) are temperature-dependent and differentially expressed at the two temperatures, suggesting a role in gene regulation for adaptation during transmission. In addition, many genes important for maintenance of Borrelia during its enzootic cycle are associated with antisense RNAs or 5' UTR sRNAs. RNA-seq data were validated for twenty-two of the sRNAs via Northern blot analyses. CONCLUSIONS: Our study demonstrates that sRNAs are abundant and differentially expressed by environmental conditions suggesting that gene regulation via sRNAs is a common mechanism utilized in B. burgdorferi. In addition, the identification of antisense and intragenic sRNAs impacts the broadly used loss-of-function genetic approach used to study gene function and increases the coding potential of a small genome. To facilitate access to the analyzed RNA-seq data we have set-up a website at http://www.cibiv.at/~niko/bbdb/ that includes a UCSC browser track hub. By clicking on the respective link, researchers can interactively inspect the data in the UCSC genome browser (Kent et al., Genome Res 12:996-1006, 2002).

Germline recessive mutations in PI4KA are associated with perisylvian polymicrogyria, cerebellar hypoplasia and arthrogryposis.

Polymicrogyria (PMG) is a structural brain abnormality involving the cerebral cortex that results from impaired neuronal migration and although several genes have been implicated, many cases remain unsolved. In this study, exome sequencing in a family where three fetuses had all been diagnosed with PMG and cerebellar hypoplasia allowed us to identify regions of the genome for which both chromosomes were shared identical-by-descent, reducing the search space for causative variants to 8.6% of the genome. In these regions, the only plausibly pathogenic mutations were compound heterozygous variants in PI4KA, which Sanger sequencing confirmed segregated consistent with autosomal recessive inheritance. The paternally transmitted variant predicted a premature stop mutation (c.2386C>T; p.R796X), whereas the maternally transmitted variant predicted a missense substitution (c.5560G>A; p.D1854N) at a conserved residue within the catalytic domain. Functional studies using expressed wild-type or mutant PI4KA enzyme confirmed the importance of p.D1854 for kinase activity. Our results emphasize the importance of phosphoinositide signalling in early brain development.

Mutations in PIGY: expanding the phenotype of inherited glycosylphosphatidylinositol deficiencies.

Glycosylphosphatidylinositol (GPI)-anchored proteins are ubiquitously expressed in the human body and are important for various functions at the cell surface. Mutations in many GPI biosynthesis genes have been described to date in patients with multi-system disease and together these constitute a subtype of congenital disorders of glycosylation. We used whole exome sequencing in two families to investigate the genetic basis of disease and used RNA and cellular studies to investigate the functional consequences of sequence variants in the PIGY gene. Two families with different phenotypes had homozygous recessive sequence variants in the GPI biosynthesis gene PIGY. Two sisters with c.137T>C (p.Leu46Pro) PIGY variants had multi-system disease including dysmorphism, seizures, severe developmental delay, cataracts and early death. There were significantly reduced levels of GPI-anchored proteins (CD55 and CD59) on the surface of patient-derived skin fibroblasts (∼20-50% compared with controls). In a second, consanguineous family, two siblings had moderate development delay and microcephaly. A homozygous PIGY promoter variant (c.-540G>A) was detected within a 7.7 Mb region of autozygosity. This variant was predicted to disrupt a SP1 consensus binding site and was shown to be associated with reduced gene expression. Mutations in PIGY can occur in coding and non-coding regions of the gene and cause variable phenotypes. This article contributes to understanding of the range of disease phenotypes and disease genes associated with deficiencies of the GPI-anchor biosynthesis pathway and also serves to highlight the potential importance of analysing variants detected in 5'-UTR regions despite their typically low coverage in exome data.

The Borrelia burgdorferi RelA/SpoT Homolog and Stringent Response Regulate Survival in the Tick Vector and Global Gene Expression during Starvation.

As the Lyme disease bacterium Borrelia burgdorferi traverses its enzootic cycle, alternating between a tick vector and a vertebrate host, the spirochete must adapt and persist in the tick midgut under prolonged nutrient stress between blood meals. In this study, we examined the role of the stringent response in tick persistence and in regulation of gene expression during nutrient limitation. Nutritionally starving B. burgdorferi in vitro increased the levels of guanosine tetraphosphate (ppGpp) and guanosine pentaphosphate (pppGpp), collectively referred to as (p)ppGpp, products of the bifunctional synthetase/hydrolase RelBbu (RelA/SpoT homolog). Conversely, returning B. burgdorferi to a nutrient-rich medium decreased (p)ppGpp levels. B. burgdorferi survival in ticks between the larval and nymph blood meals, and during starvation in vitro, was dependent on RelBbu. Furthermore, normal morphological conversion from a flat-wave shape to a condensed round body (RB) form during starvation was dependent on RelBbu; relBbu mutants more frequently formed RBs, but their membranes were compromised. By differential RNA sequencing analyses, we found that RelBbu regulates an extensive transcriptome, both dependent and independent of nutrient stress. The RelBbu regulon includes the glp operon, which is important for glycerol utilization and persistence in the tick, virulence factors and the late phage operon of the 32-kb circular plasmid (cp32) family. In summary, our data suggest that RelBbu globally modulates transcription in response to nutrient stress by increasing (p)ppGpp levels to facilitate B. burgdorferi persistence in the tick.

CODOC: efficient access, analysis and compression of depth of coverage signals.

UNLABELLED: Current data formats for the representation of depth of coverage data (DOC), a central resource for interpreting, filtering or detecting novel features in high-throughput sequencing datasets, were primarily designed for visualization purposes. This limits their applicability in stand-alone analyses of these data, mainly owing to inaccurate representation or mediocre data compression. CODOC is a novel data format and comprehensive application programming interface for efficient representation, access and analysis of DOC data. CODOC compresses these data ∼ 4-32× better than the best current comparable method by exploiting specific data characteristics while at the same time enabling more-exact signal recovery for lossy compression and very fast query answering times. AVAILABILITY AND IMPLEMENTATION: Java source code and binaries are freely available for non-commercial use at http://purl.org/bgraph/codoc.

NGC: lossless and lossy compression of aligned high-throughput sequencing data.

A major challenge of current high-throughput sequencing experiments is not only the generation of the sequencing data itself but also their processing, storage and transmission. The enormous size of these data motivates the development of data compression algorithms usable for the implementation of the various storage policies that are applied to the produced intermediate and final result files. In this article, we present NGC, a tool for the compression of mapped short read data stored in the wide-spread SAM format. NGC enables lossless and lossy compression and introduces the following two novel ideas: first, we present a way to reduce the number of required code words by exploiting common features of reads mapped to the same genomic positions; second, we present a highly configurable way for the quantization of per-base quality values, which takes their influence on downstream analyses into account. NGC, evaluated with several real-world data sets, saves 33-66% of disc space using lossless and up to 98% disc space using lossy compression. By applying two popular variant and genotype prediction tools to the decompressed data, we could show that the lossy compression modes preserve >99% of all called variants while outperforming comparable methods in some configurations.