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.

/pdf/the-genome-analysis-toolkit-a-mapreduce-framework-for-4aygfa9b66.pdf

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

Summary: Research over the last few years has revealed significant haplotype structure in the human genome. The characterization of these patterns, particularly in the context of medical genetic association studies, is becoming a routine research activity. Haploview is a software package that provides computation of linkage disequilibrium statistics and population haplotype patterns from primary genotype data in a visually appealing and interactive interface.

Availability: http://www.broad.mit.edu/mpg/haploview/

Contact: jcbarret@broad.mit.edu

/pdf/haploview-analysis-and-visualization-of-ld-and-haplotype-1p9iyir38a.pdf

Haploview: analysis and visualization of LD and haplotype maps

Motivation: DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser.

Availability: Freely available to academic users from: http://www.ub.edu/dnasp

Contact: [email protected]

DnaSP v5

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies.

https://cdr.lib.unc.edu/downloads/v692t8920

A global reference for human genetic variation.

Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

Genomic disorders and rare copy number abnormalities are identified in 15–25% of patients with syndromic conditions, but their prevalence in individuals with isolated birth defects is less clear. A...

https://royalsocietypublishing.org/doi/pdf/10.1098/rstb.2015.0406

Copy number variation as a genetic basis for heterotaxy and heterotaxy-spectrum congenital heart defects.

Gene therapy is one of several approaches that are being tested in the search for an effective anti-human immunodeficiency virus (HIV) treatment. In this strategy, a “protective” gene would be introduced into target cells, rendering them relatively resistance to the virus-induced cytopathicity. Tat and Rev are viral proteins essential for HIV gene expression. Tat increases viral gene transcription and Rev is responsible for the nuclear export of mRNA encoding structural viral proteins. A fusion protein (Trev) was constructed, joining Tat and Rev transdominant mutant gene sequences. Previously, we showed that Trev inhibits both Tat and Rev activities in Jurkat T cells. To determine whether Trev could inhibit HIV replication in primary cells, we transferred the trev gene to peripheral blood lymphocytes and challenged them with different HIV strains. Levels of HIV p24 antigen (Ag) were reduced 4- to 15-fold in cultures of Trev-CD4+ T cells infected with two HIV primary clinical isolates and were not...

Protection of Primary Human T Cells from HIV Infection by Trev: A Transdominant Fusion Gene

ZIC3, a GLI superfamily transcription factor, is involved in establishing normal embryonic left–right patterning. Multiple abnormalities in the central nervous system (CNS) and axial skeleton have also been observed in mice bearing a Zic3 null allele, mice with a Zic3 overexpression allele, and the majority of patients carrying ZIC3 mutations. Previous studies indicate that ZIC3 protein can bind to the GLI consensus binding site (GLIBS) and physically interact with GLI3, a transcription factor involved in multiple aspects of neural and skeletal development. We investigated in vitro interactions of ZIC3 with GLI3 and the effect of ZIC3 mutations identified in patients with either heterotaxy or isolated cardiovascular malformations. Electrophoresis mobility shift assay (EMSA) revealed that all five intact zinc finger (ZF) domains were necessary for binding of ZIC3 to GLIBS. Inclusion of GLIBS upstream of a basal TK promoter had no effect on the activation of the promoter by ZIC3 alone, but it enhanced the synergistic activation of ZIC3 and GLI3. Wild-type (WT) ZIC3 showed specific binding to GLI3 in GST-pull-down assays. Nonsense and frameshift ZIC3 mutants lacking one or more of the zinc finger domains did not physically interact with GST-GLI3; however, two missense mutants c.1213A>G (p.K405E, fifth ZF domain), and c.649C>G (p.P217A, conserved N-terminal domain) retained binding. Luciferase reporter assays indicated that both p.P217A and p.K405E mutants also retained coactivation with GLI3 of reporter gene expression activity, while all the GLI3-nonbinding ZIC3 mutants lacked this activity. Interestingly, no CNS or skeletal abnormalities were observed in patients bearing the p.P217A or p.K405E mutations. Hum Mutat 29(1), 99–105, 2008. Published 2007 Wiley-Liss, Inc.

Characterization of the interactions of human ZIC3 mutants with GLI3.

Biotinidase deficiency leads to a biotin-deficient state, with cardinal symptoms of ataxia, alopecia, and skin rash presenting in infancy. Previous reports of head CTs in patients with biotinidase deficiency did not note basal ganglia calcifications. We report the first case of biotinidase deficiency with basal ganglia calcifications. There were no symptoms referable to basal ganglia dysfunction.

John W. Belmont

Papers

Copy number variation as a genetic basis for heterotaxy and heterotaxy-spectrum congenital heart defects.

Protection of Primary Human T Cells from HIV Infection by Trev: A Transdominant Fusion Gene

Characterization of the interactions of human ZIC3 mutants with GLI3.

Basal ganglia calcifications in a case of biotinidase deficiency.

Identification and cloning of differentially expressed genes by long-distance differential display.