Mapping Copy Number Variation by Population Scale Genome Sequencing

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.Organismic and Evolutionary Biolog

The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data

Next-generation DNA sequencing (NGS) projects, such as the 1000 Genomes Project, are already revolutionizing our understanding of genetic variation among individuals. However, the massive data sets generated by NGS—the 1000 Genome pilot alone includes nearly five terabases—make writing feature-rich, efficient, and robust analysis tools difficult for even computationally sophisticated individuals. Indeed, many professionals are limited in the scope and the ease with which they can answer scientific questions by the complexity of accessing and manipulating the data produced by these machines. Here, we discuss our Genome Analysis Toolkit (GATK), a structured programming framework designed to ease the development of efficient and robust analysis tools for next-generation DNA sequencers using the functional programming philosophy of MapReduce. The GATK provides a small but rich set of data access patterns that encompass the majority of analysis tool needs. Separating specific analysis calculations from common data management infrastructure enables us to optimize the GATK framework for correctness, stability, and CPU and memory efficiency and to enable distributed and shared memory parallelization. We highlight the capabilities of the GATK by describing the implementation and application of robust, scale-tolerant tools like coverage calculators and single nucleotide polymorphism (SNP) calling. We conclude that the GATK programming framework enables developers and analysts to quickly and easily write efficient and robust NGS tools, many of which have already been incorporated into large-scale sequencing projects like the 1000 Genomes Project and The Cancer Genome Atlas

MTAP Deletions in Cancer Create Vulnerability to Targeting of the MAT2A/PRMT5/RIOK1 Axis

Homozygous deletions of p16/CDKN2A are prevalent in cancer, and these mutations commonly involve co-deletion of adjacent genes, including methylthioadenosine phosphorylase (MTAP). Here, we used shRNA screening and identified the metabolic enzyme, methionine adenosyltransferase II alpha (MAT2A), and the arginine methyltransferase, PRMT5, as vulnerable enzymes in cells with MTAP deletion. Metabolomic and biochemical studies revealed a mechanistic basis for this synthetic lethality. The MTAP substrate methylthioadenosine (MTA) accumulates upon MTAP loss. Biochemical profiling of a methyltransferase enzyme panel revealed that MTA is a potent and selective inhibitor of PRMT5. MTAP-deleted cells have reduced PRMT5 methylation activity and increased sensitivity to PRMT5 depletion. MAT2A produces the PRMT5 substrate S-adenosylmethionine (SAM), and MAT2A depletion reduces growth and PRMT5 methylation activity selectively in MTAP-deleted cells. Furthermore, this vulnerability extends to PRMT5 co-complex proteins such as RIOK1. Thus, the unique biochemical features of PRMT5 create an axis of targets vulnerable in CDKN2A/MTAP-deleted cancers

A small molecule inhibitor of mutant IDH2 rescues cardiomyopathy in a D-2-hydroxyglutaric aciduria type II mouse model

D-2-hydroxyglutaric aciduria (D2HGA) type II is a rare neurometabolic disorder caused by germline gain-of-function mutations in isocitrate dehydrogenase 2 (IDH2), resulting in accumulation of D-2-hydroxyglutarate (D2HG). Patients exhibit a wide spectrum of symptoms including cardiomyopathy, epilepsy, developmental delay and limited life span. Currently, there are no effective therapeutic interventions. We generated a D2HGA type II mouse model by introducing the Idh2R140Q mutation at the native chromosomal locus. Idh2R140Q mice displayed significantly elevated 2HG levels and recapitulated multiple defects seen in patients. AGI-026, a potent, selective inhibitor of the human IDH2R140Q-mutant enzyme, suppressed 2HG production, rescued cardiomyopathy, and provided a survival benefit in Idh2R140Q mice; treatment withdrawal resulted in deterioration of cardiac function. We observed differential expression of multiple genes and metabolites that are associated with cardiomyopathy, which were largely reversed by AGI-026. These findings demonstrate the potential therapeutic benefit of an IDH2R140Q inhibitor in patients with D2HGA type II