Mapping and sequencing of structural variation from eight human genomes.

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

A method for rapid, targeted CNV genotyping identifies rare variants associated with neurocognitive disease.

Copy-number variants (CNVs) are substantial contributors to human disease. A central challenge in CNV-disease association studies is to characterize the pathogenicity of rare and possibly incompletely penetrant events, which requires the accurate detection of rare CNVs in large numbers of individuals. Cost and throughput issues limit our ability to perform these studies. We have adapted the Illumina BeadXpress SNP genotyping assay and developed an algorithm, SNP-Conditional OUTlier detection (SCOUT), to rapidly and accurately detect both rare and common CNVs in large cohorts. This approach is customizable, cost effective, highly parallelized, and largely automated. We applied this method to screen 69 loci in 1105 children with unexplained intellectual disability, identifying pathogenic variants in 3.1% of these individuals and potentially pathogenic variants in an additional 2.3%. We identified seven individuals (0.7%) with a deletion of 16p11.2, which has been previously associated with autism. Our results widen the phenotypic spectrum of these deletions to include intellectual disability without autism. We also detected 1.65-3.4 Mbp duplications at 16p13.11 in 1.1% of affected individuals and 350 kbp deletions at 15q11.2, near the Prader-Willi/Angelman syndrome critical region, in 0.8% of affected individuals. Compared to published CNVs in controls they are significantly (P = 4.7 x 10(-5) and 0.003, respectively) enriched in these children, supporting previously published hypotheses that they are neurocognitive disease risk factors. More generally, this approach offers a previously unavailable balance between customization, cost, and throughput for analysis of CNVs and should prove valuable for targeted CNV detection in both research and diagnostic settings.

Targeted interrogation of copy number variation using SCIMMkit.

SUMMARY: Copy number variants (CNVs) contribute substantially to human genomic diversity, and development of accurate and efficient methods for CNV genotyping is a central problem in exploring human genotype-phenotype associations. SCIMMkit provides a robust, integrated implementation of three previously validated algorithms [SCIMM (SNP-Conditional Mixture Modeling), SCIMM-Search and SCOUT (SNP-Conditional OUTlier detection)] for targeted interrogation of CNVs using Illumina Infinium II and GoldenGate SNP assays. SCIMMkit is applicable to standardized genome-wide SNP arrays and customized multiplexed SNP panels, providing economy, efficiency and flexibility in experimental design. AVAILABILITY: Source code and documentation are available for noncommercial use at http://droog.gs.washington.edu/scimmkit.

High-throughput discovery of rare human nucleotide polymorphisms by Ecotilling.

Human individuals differ from one another at only approximately 0.1% of nucleotide positions, but these single nucleotide differences account for most heritable phenotypic variation. Large-scale efforts to discover and genotype human variation have been limited to common polymorphisms. However, these efforts overlook rare nucleotide changes that may contribute to phenotypic diversity and genetic disorders, including cancer. Thus, there is an increasing need for high-throughput methods to robustly detect rare nucleotide differences. Toward this end, we have adapted the mismatch discovery method known as Ecotilling for the discovery of human single nucleotide polymorphisms. To increase throughput and reduce costs, we developed a universal primer strategy and implemented algorithms for automated band detection. Ecotilling was validated by screening 90 human DNA samples for nucleotide changes in 5 gene targets and by comparing results to public resequencing data. To increase throughput for discovery of rare alleles, we pooled samples 8-fold and found Ecotilling to be efficient relative to resequencing, with a false negative rate of 5% and a false discovery rate of 4%. We identified 28 new rare alleles, including some that are predicted to damage protein function. The detection of rare damaging mutations has implications for models of human disease.

Automated band mapping in electrophoretic gel images using background information.

Some popular methods for polymorphism and mutation discovery involve ascertainment of novel bands by the examination of electrophoretic gel images. Although existing strategies for mapping bands work well for specific applications, such as DNA sequencing, these strategies are not well suited for novel band detection. Here, we describe a general strategy for band mapping that uses background banding patterns to facilitate lane calling and size calibration. We have implemented this strategy in GelBuddy, a user-friendly Java-based program for PC and Macintosh computers, which includes several utilities to assist discovery of mutations and polymorphisms. We demonstrate the use of GelBuddy in applications based on single-base mismatch cleavage of heteroduplexed PCR products. Use of software designed to facilitate novel band detection can significantly shorten the time needed for image analysis and data entry in a high-throughput setting. Furthermore, the interactive strategy implemented in GelBuddy has been successfully applied to DNA fingerprinting applications, such as AFLP. GelBuddy promises to make electrophoretic gel analysis a viable alternative to DNA resequencing for discovery of mutations and polymorphisms.

Systematic assessment of copy number variant detection via genome-wide SNP genotyping.

SNP genotyping has emerged as a technology to incorporate copy number variants (CNVs) into genetic analyses of human traits. However, the extent to which SNP platforms accurately capture CNVs remains unclear. Using independent, sequence-based CNV maps, we find that commonly used SNP platforms have limited or no probe coverage for a large fraction of CNVs. Despite this, in 9 samples we inferred 368 CNVs using Illumina SNP genotyping data and experimentally validated over two-thirds of these. We also developed a method (SNP-Conditional Mixture Modeling, SCIMM) to robustly genotype deletions using as few as two SNP probes. We find that HapMap SNPs are strongly correlated with 82% of common deletions, but the newest SNP platforms effectively tag about 50%. We conclude that currently available genome-wide SNP assays can capture CNVs accurately, but improvements in array designs, particularly in duplicated sequences, are necessary to facilitate more comprehensive analyses of genomic variation.

A protocol for TILLING and Ecotilling in plants and animals.

We describe Targeting-Induced Local Lesions IN Genomes (TILLING), a reverse-genetic strategy for the discovery and mapping of induced mutations. TILLING is suitable for essentially any organism that can be mutagenized. The TILLING procedure has also been adapted for the discovery and cataloguing of natural polymorphisms, a method called Ecotilling. To discover nucleotide changes within a particular gene, PCR is performed with gene-specific primers that are end-labeled with fluorescent molecules. After PCR, samples are denatured and annealed to form heteroduplexes between polymorphic DNA strands. Mismatched base pairs in these heteroduplexes are cleaved by digestion with a single-strand specific nuclease. The resulting products are size-fractionated using denaturing polyacrylamide gel electrophoresis and visualized by fluorescence detection. The migration of cleaved products indicates the approximate location of nucleotide polymorphisms. Throughput is increased and costs are reduced by sample pooling, multi-well liquid handling and automated gel band mapping. Once genomic DNA samples have been obtained, pooled and arrayed, thousands of samples can be screened daily.