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.

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.

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).

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.

Subangstrom single-molecule measurements of motor proteins using a nanopore.

Techniques for measuring the motion of single motor proteins, such as FRET and optical tweezers, are limited to a resolution of ∼300 pm. We use ion current modulation through the protein nanopore MspA to observe translocation of helicase Hel308 on DNA with up to ∼40 pm sensitivity. This approach should be applicable to any protein that translocates on DNA or RNA, including helicases, polymerases, recombinases and DNA repair enzymes.

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.

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.

Haplotype-resolved whole-genome sequencing by contiguity-preserving transposition and combinatorial indexing.

Haplotype-resolved genome sequencing enables the accurate interpretation of medically relevant genetic variation, deep inferences regarding population history and non-invasive prediction of fetal genomes. We describe an approach for genome-wide haplotyping based on contiguity-preserving transposition (CPT-seq) and combinatorial indexing. Tn5 transposition is used to modify DNA with adaptor and index sequences while preserving contiguity. After DNA dilution and compartmentalization, the transposase is removed, resolving the DNA into individually indexed libraries. The libraries in each compartment, enriched for neighboring genomic elements, are further indexed via PCR. Combinatorial 96-plex indexing at both the transposition and PCR stage enables the construction of phased synthetic reads from each of the nearly 10,000 'virtual compartments'. We demonstrate the feasibility of this method by assembling >95% of the heterozygous variants in a human genome into long, accurate haplotype blocks (N50 = 1.4-2.3 Mb). The rapid, scalable and cost-effective workflow could enable haplotype resolution to become routine in human genome sequencing.

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.

Mobile genetic elements and genome evolution 2014.

The Mobile Genetic Elements and Genome Evolution conference was hosted by Keystone Symposia in Santa Fe, NM USA, 9 March through 14 March 2014. The goal of this conference was to bring together scientists from around the world who study transposable elements in diverse organisms and researchers who study the impact these elements have on genome evolution. The meeting included over 200 scientists who participated through poster presentations, short talks selected from abstracts, and invited speakers. The talks were organized into eight sessions and two workshops. The topics varied from diverse mechanisms of mobilization to the evolution of genomes and their defense strategies against transposable elements.

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.

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.

Genome-wide assessment of imprinted expression in human cells.

BACKGROUND: Parent-of-origin-dependent expression of alleles, imprinting, has been suggested to impact a substantial proportion of mammalian genes. Its discovery requires allele-specific detection of expressed transcripts, but in some cases detected allelic expression bias has been interpreted as imprinting without demonstrating compatible transmission patterns and excluding heritable variation. Therefore, we utilized a genome-wide tool exploiting high density genotyping arrays in parallel measurements of genotypes in RNA and DNA to determine allelic expression across the transcriptome in lymphoblastoid cell lines (LCLs) and skin fibroblasts derived from families. RESULTS: We were able to validate 43% of imprinted genes with previous demonstration of compatible transmission patterns in LCLs and fibroblasts. In contrast, we only validated 8% of genes suggested to be imprinted in the literature, but without clear evidence of parent-of-origin-determined expression. We also detected five novel imprinted genes and delineated regions of imprinted expression surrounding annotated imprinted genes. More subtle parent-of-origin-dependent expression, or partial imprinting, could be verified in four genes. Despite higher prevalence of monoallelic expression, immortalized LCLs showed consistent imprinting in fewer loci than primary cells. Random monoallelic expression has previously been observed in LCLs and we show that random monoallelic expression in LCLs can be partly explained by aberrant methylation in the genome. CONCLUSIONS: Our results indicate that widespread parent-of-origin-dependent expression observed recently in rodents is unlikely to be captured by assessment of human cells derived from adult tissues where genome-wide assessment of both primary and immortalized cells yields few new imprinted loci.

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.

Deletion at fragile sites is a common and early event in Barrett's esophagus.

Barrett's esophagus (BE) is a premalignant intermediate to esophageal adenocarcinoma, which develops in the context of chronic inflammation and exposure to bile and acid. We asked whether there might be common genomic alterations that could be identified as potential clinical biomarker(s) for BE by whole genome profiling. We detected copy number alterations and/or loss of heterozygosity at 56 fragile sites in 20 patients with premalignant BE. Chromosomal fragile sites are particularly sensitive to DNA breaks and are frequent sites of rearrangement or loss in many human cancers. Seventy-eight percent of all genomic alterations detected by array-CGH were associated with fragile sites. Copy number losses in early BE were observed at particularly high frequency at FRA3B (81%), FRA9A/C (71.4%), FRA5E (52.4%), and FRA 4D (52.4%), and at lower frequencies in other fragile sites, including FRA1K (42.9%), FRAXC (42.9%), FRA 12B (33.3%), and FRA16D (33.3%). Due to the consistency of the region of copy number loss, we were able to verify these results by quantitative PCR, which detected the loss of FRA3B and FRA16D, in 83% and 40% of early molecular stage BE patients, respectively. Loss of heterozygosity in these cases was confirmed through pyrosequencing at FRA3B and FRA16D (75% and 70%, respectively). Deletion and genomic instability at FRA3B and other fragile sites could thus be a biomarker of genetic damage in BE patients and a potential biomarker of cancer risk.

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.

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.

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.

Increased LIS1 expression affects human and mouse brain development.

Deletions of the PAFAH1B1 gene (encoding LIS1) in 17p13.3 result in isolated lissencephaly sequence, and extended deletions including the YWHAE gene (encoding 14-3-3epsilon) cause Miller-Dieker syndrome. We identified seven unrelated individuals with submicroscopic duplication in 17p13.3 involving the PAFAH1B1 and/or YWHAE genes, and using a 'reverse genomics' approach, characterized the clinical consequences of these duplications. Increased PAFAH1B1 dosage causes mild brain structural abnormalities, moderate to severe developmental delay and failure to thrive. Duplication of YWHAE and surrounding genes increases the risk for macrosomia, mild developmental delay and pervasive developmental disorder, and results in shared facial dysmorphologies. Transgenic mice conditionally overexpressing LIS1 in the developing brain showed a decrease in brain size, an increase in apoptotic cells and a distorted cellular organization in the ventricular zone, including reduced cellular polarity but preserved cortical cell layer identity. Collectively, our results show that an increase in LIS1 expression in the developing brain results in brain abnormalities in mice and humans.

Allele-specific gene expression patterns in primary leukemic cells reveal regulation of gene expression by CpG site methylation.

To identify genes that are regulated by cis-acting functional elements in acute lymphoblastic leukemia (ALL) we determined the allele-specific expression (ASE) levels of 2, 529 genes by genotyping a genome-wide panel of single nucleotide polymorphisms in RNA and DNA from bone marrow and blood samples of 197 children with ALL. Using a reproducible, quantitative genotyping method and stringent criteria for scoring ASE, we found that 16% of the analyzed genes display ASE in multiple ALL cell samples. For most of the genes, the level of ASE varied largely between the samples, from 1.4-fold overexpression of one allele to apparent monoallelic expression. For genes exhibiting ASE, 55% displayed bidirectional ASE in which overexpression of either of the two SNP alleles occurred. For bidirectional ASE we also observed overall higher levels of ASE and correlation with the methylation level of these sites. Our results demonstrate that CpG site methylation is one of the factors that regulates gene expression in ALL cells.