Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II.

The most polymorphic part of the human genome, the MHC, encodes over 160 proteins of diverse function. Half of them, including the HLA class I and II genes, are directly involved in immune responses. Consequently, the MHC region strongly associates with numerous diseases and clinical therapies. Notoriously, the MHC region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp MHC region from genomic DNA. For 95 MHC homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative MHC reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the MHC region shows the approach accurately determines the sequences of the highly polymorphic HLA class I and HLA class II genes and the complex structural diversity of complement factor C4A/C4B It has also uncovered extensive and unexpected diversity in other MHC genes; an example is MUC22, which encodes a lung mucin and exhibits more coding sequence alleles than any HLA class I or II gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference MHC haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

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.

In vitro, long-range sequence information for de novo genome assembly via transposase contiguity.

We describe a method that exploits contiguity preserving transposase sequencing (CPT-seq) to facilitate the scaffolding of de novo genome assemblies. CPT-seq is an entirely in vitro means of generating libraries comprised of 9216 indexed pools, each of which contains thousands of sparsely sequenced long fragments ranging from 5 kilobases to > 1 megabase. These pools are "subhaploid," in that the lengths of fragments contained in each pool sums to ∼5% to 10% of the full genome. The scaffolding approach described here, termed fragScaff, leverages coincidences between the content of different pools as a source of contiguity information. Specifically, CPT-seq data is mapped to a de novo genome assembly, followed by the identification of pairs of contigs or scaffolds whose ends disproportionately co-occur in the same indexed pools, consistent with true adjacency in the genome. Such candidate "joins" are used to construct a graph, which is then resolved by a minimum spanning tree. As a proof-of-concept, we apply CPT-seq and fragScaff to substantially boost the contiguity of de novo assemblies of the human, mouse, and fly genomes, increasing the scaffold N50 of de novo assemblies by eight- to 57-fold with high accuracy. We also demonstrate that fragScaff is complementary to Hi-C-based contact probability maps, providing midrange contiguity to support robust, accurate chromosome-scale de novo genome assemblies without the need for laborious in vivo cloning steps. Finally, we demonstrate CPT-seq as a means of anchoring unplaced novel human contigs to the reference genome as well as for detecting misassembled sequences.

Very long haplotype tracts characterized at high resolution from HLA homozygous cell lines.

The HLA region of chromosome 6 contains the most polymorphic genes in humans. Spanning ~5 Mbp the densely packed region encompasses approximately 175 expressed genes including the highly polymorphic HLA class I and II loci. Most of the other genes and functional elements are also polymorphic, and many of them are directly implicated in immune function or immune-related disease. For these reasons, this complex genomic region is subject to intense scrutiny by researchers with the common goal of aiding further understanding and diagnoses of multiple immune-related diseases and syndromes. To aid assay development and characterization of the classical loci, a panel of cell lines partially or fully homozygous for HLA class I and II was assembled over time by the International Histocompatibility Working Group (IHWG). Containing a minimum of 88 unique HLA haplotypes, we show that this panel represents a significant proportion of European HLA allelic and haplotype diversity (60-95 %). Using a high-density whole genome array that includes 13,331 HLA region SNPs, we analyzed 99 IHWG cells to map the coordinates of the homozygous tracts at a fine scale. The mean homozygous tract length within chromosome 6 from these individuals is 21 Mbp. Within HLA, the mean haplotype length is 4.3 Mbp, and 65 % of the cell lines were shown to be homozygous throughout the entire region. In addition, four cell lines are homozygous throughout the complex KIR region of chromosome 19 (~250 kbp). The data we describe will provide a valuable resource for characterizing haplotypes, designing and refining imputation algorithms and developing assay controls.

Highly parallel oligonucleotide purification and functionalization using reversible chemistry.

We have developed a cost-effective, highly parallel method for purification and functionalization of 5'-labeled oligonucleotides. The approach is based on 5'-hexa-His phase tag purification, followed by exchange of the hexa-His tag for a functional group using reversible reaction chemistry. These methods are suitable for large-scale (micromole to millimole) production of oligonucleotides and are amenable to highly parallel processing of many oligonucleotides individually or in high complexity pools. Examples of the preparation of 5'-biotin, 95-mer, oligonucleotide pools of >40K complexity at micromole scale are shown. These pools are prepared in up to ~16% yield and 90-99% purity. Approaches for using this method in other applications are also discussed.

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.

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.

Computational analysis of whole-genome differential allelic expression data in human.

Allelic imbalance (AI) is a phenomenon where the two alleles of a given gene are expressed at different levels in a given cell, either because of epigenetic inactivation of one of the two alleles, or because of genetic variation in regulatory regions. Recently, Bing et al. have described the use of genotyping arrays to assay AI at a high resolution (approximately 750,000 SNPs across the autosomes). In this paper, we investigate computational approaches to analyze this data and identify genomic regions with AI in an unbiased and robust statistical manner. We propose two families of approaches: (i) a statistical approach based on z-score computations, and (ii) a family of machine learning approaches based on Hidden Markov Models. Each method is evaluated using previously published experimental data sets as well as with permutation testing. When applied to whole genome data from 53 HapMap samples, our approaches reveal that allelic imbalance is widespread (most expressed genes show evidence of AI in at least one of our 53 samples) and that most AI regions in a given individual are also found in at least a few other individuals. While many AI regions identified in the genome correspond to known protein-coding transcripts, others overlap with recently discovered long non-coding RNAs. We also observe that genomic regions with AI not only include complete transcripts with consistent differential expression levels, but also more complex patterns of allelic expression such as alternative promoters and alternative 3' end. The approaches developed not only shed light on the incidence and mechanisms of allelic expression, but will also help towards mapping the genetic causes of allelic expression and identify cases where this variation may be linked to diseases.

Power to detect risk alleles using genome-wide tag SNP panels.

Advances in high-throughput genotyping and the International HapMap Project have enabled association studies at the whole-genome level. We have constructed whole-genome genotyping panels of over 550,000 (HumanHap550) and 650,000 (HumanHap650Y) SNP loci by choosing tag SNPs from all populations genotyped by the International HapMap Project. These panels also contain additional SNP content in regions that have historically been overrepresented in diseases, such as nonsynonymous sites, the MHC region, copy number variant regions and mitochondrial DNA. We estimate that the tag SNP loci in these panels cover the majority of all common variation in the genome as measured by coverage of both all common HapMap SNPs and an independent set of SNPs derived from complete resequencing of genes obtained from SeattleSNPs. We also estimate that, given a sample size of 1,000 cases and 1,000 controls, these panels have the power to detect single disease loci of moderate risk (lambda approximately 1.8-2.0). Relative risks as low as lambda approximately 1.1-1.3 can be detected using 10,000 cases and 10,000 controls depending on the sample population and disease model. If multiple loci are involved, the power increases significantly to detect at least one locus such that relative risks 20%-35% lower can be detected with 80% power if between two and four independent loci are involved. Although our SNP selection was based on HapMap data, which is a subset of all common SNPs, these panels effectively capture the majority of all common variation and provide high power to detect risk alleles that are not represented in the HapMap data.

Whole-genome genotyping on bead arrays.

In this review, we describe the laboratory implementation of Infinium whole genome genotyping (WGG) technology for whole genome association studies and copy number studies. Briefly, the Infinium WGG assay employs a single tube whole genome amplification reaction to amplify the entire genome; genomic loci of interest are captured on an array by specific hybridization of picomolar concentrations amplified gDNA. After target capture, single nucleotide polymorphisms (SNPs) are genotyped on the array by a primer extension reaction using hapten-labeled nucleotides. The resultant hapten signal is amplified by immunhistochemical sandwich staining and the array is read out on a high resolution confocal scanner. We have combined this Infinium assay with high-density BeadChips to create the first array platform capable of genotyping over 1 million SNPs per slide. Additionally, the complete Infinium assay is automated using Tecan GenePaint slide processing system. Hybridization, washing, array-based primer extension and staining are performed directly in the Tecan capillary gap Te-Flow Through chambers. This automation process greatly increases assay robustness and throughput while enabling Laboratory Information Management System (LIMS) control of sample tracking. Finally, we give several examples of how this advance in genotyping technology is being applied in whole genome association and copy number studies.

Design of tag SNP whole genome genotyping arrays.

Whole genome association studies have recently been enabled by combining tag SNP information derived from the International HapMap project with novel whole genome genotyping array technologies. In particular, Infinium whole genome genotyping (WGG) technology now has the power to genotype over 1 million SNPs on a single array. Additionally, this assay provides access to virtually any SNP in the genome enabling selection of optimized SNP content . In this chapter, we provide an overview of the tag SNP-based selection strategy for Infinium whole-genome genotyping BeadChips, including the Human 1 M BeadChip. These advances in both SNP content and technology have enabled both large-scale whole-genome disease association (WGAS) and copy number variation (CNV) studies with the ultimate goal of identifying common genetic variants, disease-associated loci, proteins, and biomarkers.

Delineation of the proximal 3q microdeletion syndrome.

Interstitial deletions of the proximal long arm of chromosome 3 are very rare and a defined clinical phenotype is not established yet. We report on the clinical, cytogenetic and molecular findings of a 20-month-old Hispanic male with a 2.5 Mb de novo deletion on q13.11q13.12. Up to now, this is the smallest deletion reported among patients with the proximal 3q microdeletion syndrome. The patient has distinct facial features including brachycephaly, broad and prominent forehead, flat nasal bridge, prominent ears, anteverted nose, tetralogy of Fallot, bilateral cryptorchidism, and peripheral skeletal abnormalities. To further delineate the proximal 3q deletion syndrome, the phenotype of our patient was compared with 10 other patients previously described. We found that ALCAM and CBLB are the only genes deleted in our patient and based on previously published data, we propose that the CBLB gene is responsible for the craniofacial phenotype in patients with deletions of proximal 3q region.

Mosaic tetrasomy 12p with triplication of 12p detected by array-based comparative genomic hybridization of peripheral blood DNA.

A patient whose dysmorphism at birth was not diagnostic for Pallister-Killian syndrome (PKS) was found to have mosaic tetrasomy 12p by an array-based comparative genomic hybridization of peripheral blood DNA. He was determined to be mosaic for 46,XY,trp(12)(p11.2 --> p13) in cultured skin fibroblasts. His appearance was typical for PKS at 4 months of age.

Contiguity-Preserving Transposition Sequencing (CPT-Seq) for Genome-Wide Haplotyping, Assembly, and Single-Cell ATAC-Seq.

Most genomes to date have been sequenced without taking into account the diploid nature of the genome. However, the distribution of variants on each individual chromosome can (1) significantly impact gene regulation and protein function, (2) have important implications for analyses of population history and medical genetics, and (3) be of great value for accurate interpretation of medically relevant genetic variation. Here, we describe a comprehensive and detailed protocol for an ultra fast (<3 h library preparation), cost-effective, and scalable haplotyping method, named Contiguity Preserving Transposition sequencing or CPT-seq (Amini et al., Nat Genet 46(12):1343-1349, 2014). CPT-seq accurately phases >95 % of the whole human genome in Mb-scale phasing blocks. Additionally, the same workflow can be used to aid de novo assembly (Adey et al., Genome Res 24(12):2041-2049, 2014), detect structural variants, and perform single cell ATAC-seq analysis (Cusanovich et al., Science 348(6237):910-914, 2015).

Haplotype phasing of whole human genomes using bead-based barcode partitioning in a single tube.

Haplotype-resolved genome sequencing promises to unlock a wealth of information in population and medical genetics. However, for the vast majority of genomes sequenced to date, haplotypes have not been determined because of cumbersome haplotyping workflows that require fractions of the genome to be sequenced in a large number of compartments. Here we demonstrate barcode partitioning of long DNA molecules in a single compartment using "on-bead" barcoded tagmentation. The key to the method that we call "contiguity preserving transposition" sequencing on beads (CPTv2-seq) is transposon-mediated transfer of homogenous populations of barcodes from beads to individual long DNA molecules that get fragmented at the same time (tagmentation). These are then processed to sequencing libraries wherein all sequencing reads originating from each long DNA molecule share a common barcode. Single-tube, bulk processing of long DNA molecules with ∼150,000 different barcoded bead types provides a barcode-linked read structure that reveals long-range molecular contiguity. This technology provides a simple, rapid, plate-scalable and automatable route to accurate, haplotype-resolved sequencing, and phasing of structural variants of the genome.

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.

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 Tecan's GenePaint slide processing system. Hybridization, washing, array-based primer extension, and staining are performed directly in Tecan's capillary gap Te-Flow chambers. This automation process increases assay robustness and throughput greatly while enabling laboratory information management system control of sample tracking.

Whole-genome genotyping of haplotype tag single nucleotide polymorphisms.

The International HapMap Consortium recently completed genotyping over 3.8 million single nucleotide polymorphisms (SNPs) in three major populations, and the results of studying patterns of linkage disequilibrium indicate that characterization of 300,000-500,000 tag SNPs is sufficient to provide good genomic coverage for linkage-disequilibrium-based association studies in many populations. These whole-genome association studies will be used to dissect the genetics of complex diseases and pharmacogenomic drug responses. As such, the development of a cost-effective genotyping platform that can assay hundred of thousands of SNPs across thousands of samples is essential. In this review, we describe the development of a whole-genome genotyping (WGG) assay that enables unconstrained SNP selection and effectively unlimited multiplexing from a single sample preparation. The development of WGG in concert with high-density BeadChips has enabled the creation of three different high-density SNP genotyping BeadChips: the Sentrix Human-1 Genotyping BeadChip containing over 109,000 exon-centric SNPs; the HumanHap300 BeadChip containing over 317,000 tag SNPs, and the HumanHap550 Beadchip containing over 550,000 tag SNPs.

Illumina, Inc.

Illumina, Inc., based in San Diego (CA, USA), is a genomics tool company that develops and markets integrated array-based systems and assays for a broad range of applications including genotyping, gene expression and epigenetics. Product offerings range from focused assay sets (up to 1,536 multiplexed assays) to whole-genome analysis (>100,000 assays/sample). Illumina's two microarray platforms, the Sentrix Array Matrix and the Sentrix BeadChip, are characterized by small (3 microm) feature size, dense feature packing (over 10 million features can be deployed on a single microarray) and the ability to analyze, in parallel, multiple samples on the same device. Illumina has developed a spectrum of proprietary assays (GoldenGate, Infinium, and DASL) for application on our microarray platform. These assays have been successfully employed for a number of applications, including collection of the majority of the Phase I genotyping data for the International HapMap Project. Illumina are focusing our activities on extending the applications of the above assays and developing new assays for future products. Illumina's goal is to deliver high-performance, high-throughput solutions that enable researchers to expand experimental scale while reducing the cost of large-scale research.

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

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.