Kevin L. Gunderson

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R2 of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.

A genome-wide scalable SNP genotyping assay using microarray technology

Oligonucleotide probe arrays have enabled massively parallel analysis of gene expression levels from a single cDNA sample. Application of microarray technology to analyzing genomic DNA has been stymied by the sequence complexity of the entire human genome. A robust, single base–resolution direct genomic assay would extend the reach of microarray technology. We developed an array-based whole-genome genotyping assay that does not require PCR and enables effectively unlimited multiplexing. The assay achieves a high signal-to-noise ratio by combining specific hybridization of picomolar concentrations of whole genome–amplified DNA to arrayed probes with allele-specific primer extension and signal amplification. As proof of principle, we genotyped several hundred previously characterized SNPs. The conversion rate, call rate and accuracy were comparable to those of high-performance PCR-based genotyping assays.

Highly parallel genomic assays

Recent developments in highly parallel genome-wide assays are transforming the study of human health and disease. High-resolution whole-genome association studies of complex diseases are finally being undertaken after much hypothesizing about their merit for finding disease loci. The availability of inexpensive high-density SNP-genotyping arrays has made this feasible. Cancer biology will also be transformed by high-resolution genomic and epigenomic analysis. In the future, most cancers might be staged by high-resolution molecular profiling rather than by gross cytological analysis. Here, we describe the key developments that enable highly parallel genomic assays.

Whole-genome genotyping with the single-base extension assay.

We describe an efficient, accurate and robust whole-genome genotyping (WGG) assay based on a two-color, single-base extension (SBE), single-nucleotide polymorphism (SNP)-scoring step. We report genotyping results for biallelic International HapMap quality control (QC) SNPs using a single probe per locus. We show scalability, throughput and accuracy of the system by resequencing homozygous loci from our 100k Human-1 Genotyping BeadChip.

Illumina universal bead arrays.

This chapter describes an accurate, scalable, and flexible microarray technology. It includes a miniaturized array platform where each individual feature is quality controlled and a versatile assay that can be adapted for various genetic analyses, such as single nucleotide polymorphism genotyping, DNA methylation detection, and gene expression profiling. This chapter describes the concept of the BeadArray technology, two different Array of Arrays formats, the assay scheme and protocol, the performance of the system, and its use in large-scale genetic, epigenetic, and expression studies.

Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing

Chromatin state and the single cell Identifying the chromatin state of any single cell, which may or may not have a different function or represent different stages relative to others collected within any single culture, experiment, or tissue, has been challenging. Cusanovitch et al. skirted previously identified technological limitations to identify regions of accessible chromatin at single-cell resolution. Combinatorial cellular indexing, a strategy for multiplex barcoding of thousands of single cells per experiment, was successfully used to investigate the genome-wide chromatin accessibility landscape in each of over 15,000 single cells. Science, this issue p. 910 Combinatorial indexing can identify chromatin states at single-cell resolution. Technical advances have enabled the collection of genome and transcriptome data sets with single-cell resolution. However, single-cell characterization of the epigenome has remained challenging. Furthermore, because cells must be physically separated before biochemical processing, conventional single-cell preparatory methods scale linearly. We applied combinatorial cellular indexing to measure chromatin accessibility in thousands of single cells per assay, circumventing the need for compartmentalization of individual cells. We report chromatin accessibility profiles from more than 15,000 single cells and use these data to cluster cells on the basis of chromatin accessibility landscapes. We identify modules of coordinately regulated chromatin accessibility at the level of single cells both between and within cell types, with a scalable method that may accelerate progress toward a human cell atlas.

Global patterns of cis variation in human cells revealed by high-density allelic expression analysis.

Cis-acting variants altering gene expression are a source of phenotypic differences. The cis-acting components of expression variation can be identified through the mapping of differences in allelic expression (AE), which is the measure of relative expression between two allelic transcripts. We generated a map of AE associated SNPs using quantitative measurements of AE on Illumina Human1M BeadChips. In 53 lymphoblastoid cell lines derived from donors of European descent, we identified common cis variants affecting 30% (2935/9751) of the measured RefSeq transcripts at 0.001 permutation significance. The pervasive influence of cis-regulatory variants, which explain 50% of population variation in AE, extend to full-length transcripts and their isoforms as well as to unannotated transcripts. These strong effects facilitate fine mapping of cis-regulatory SNPs, as demonstrated by dissection of heritable control of transcripts in the systemic lupus erythematosus–associated C8orf13-BLK region in chromosome 8. The dense collection of associations will facilitate large-scale isolation of cis-regulatory SNPs.

Speech delay and autism spectrum behaviors are frequently associated with duplication of the 7q11.23 Williams-Beuren syndrome region

Purpose: Williams-Beuren syndrome is among the most well-characterized microdeletion syndromes, caused by recurrent de novo microdeletions at 7q11.23 mediated by nonallelic homologous recombination between low copy repeats flanking this critical region. However, the clinical phenotype associated with reciprocal microduplication of this genomic region is less well described. We investigated the molecular, clinical, neurodevelopmental, and behavioral features of seven patients with dup(7)(q11.23), including two children who inherited the microduplication from one of their parents, to more fully characterize this emerging microduplication syndrome.Methods: Patients were identified by array-based comparative genomic hybridization. Clinical examinations were performed on seven affected probands, and detailed cognitive and behavioral evaluations were carried out on four of the affected probands.Results: Our findings confirm initial reports of speech delay seen in patients with dup(7)(q11.23) and further delineate and expand the phenotypic spectrum of this condition to include communication, social interactions, and repetitive interests that are often observed in individuals diagnosed with autism spectrum disorders.Conclusions: Array-based comparative genomic hybridization is a powerful means of detecting genomic imbalances and identifying molecular etiologies in the clinic setting, including genomic disorders such as Williams-Beuren syndrome and dup(7)(q11.23). We propose that dup(7)(q11.23) syndrome may be as frequent as Williams-Beuren syndrome and a previously unrecognized cause of language delay and behavioral abnormalities. Indeed, these individuals may first be referred for evaluation of autism, even if they do not ultimately meet diagnostic criteria for an autism spectrum disorder.

Whole‐Genome Genotyping

We have developed an array-based whole-genome genotyping (WGG) assay (Infinium) using our BeadChip platform that effectively enables unlimited multiplexing and unconstrained single nucleotide polymorphism (SNP) selection. A single tube whole-genome amplification reaction is used to amplify the genome, and loci of interest are captured by specific hybridization of amplified gDNA to 50-mer probe arrays. After target capture, SNPs are genotyped on the array by a primer extension reaction in the presence of hapten-labeled nucleotides. The resultant signal is amplified during staining and the array is read out on a high-resolution confocal scanner. We have employed our high-density BeadChips supporting up to 288,000 bead types to create an array that can query over 100,000 SNPs using the Infinium assay. In addition, we have developed an automated BeadChip processing platform using Tecans GenePaint slide processing system. Hybridization, washing, array-based primer extension, and staining are performed directly in Tecans capillary gap Te-Flow chambers. This automation process increases assay robustness and throughput greatly while enabling laboratory information management system control of sample tracking.

Massively multiplex single-cell Hi-C

We present single-cell combinatorial indexed Hi-C (sciHi-C), a method that applies combinatorial cellular indexing to chromosome conformation capture. In this proof of concept, we generate and sequence six sciHi-C libraries comprising a total of 10,696 single cells. We use sciHi-C data to separate cells by karyotypic and cell-cycle state differences and identify cell-to-cell heterogeneity in mammalian chromosomal conformation. Our results demonstrate that combinatorial indexing is a generalizable strategy for single-cell genomics.