Diploid human cells contain two non-identical genomes, and differences in their regulation underlie human development and disease. We present Fiber-seq Inferred Regulatory Elements (FIRE) and show that FIRE provides a more comprehensive and quantitative snapshot of the accessible chromatin landscape across the 6 Gbp diploid human genome, overcoming previously unrecognized biases in existing regulatory element catalogs. FIRE enables comprehensive detection of haplotype-selective chromatin accessibility (HSCA), exposing novel imprinted elements lacking underlying parent-of-origin CpG methylation differences, and gene regulatory modules that permit genes to escape X chromosome inactivation. We uncover that the human leukocyte antigen (HLA) locus harbors the most HSCA in immune cells, where we resolve specific transcription factor (TF) binding events disrupted by disease-associated variants. Finally, we demonstrate that the regulatory landscape of a cell is littered with autosomal somatic chromatin epimutations that are propagated by clonal expansions to create mitotically stable and non-genetically deterministic chromatin alterations.
Performance of Amplified DNA in an Illumina GoldenGate BeadArray Assay
Abstract Whole genome amplification (WGA) offers a means to enrich DNA quantities for epidemiologic studies. We used an ovarian cancer study of 1,536 single nucleotide polymorphisms (SNPs) and 2,368 samples to assess performance of multiple displacement amplification (MDA) WGA using an Illumina GoldenGate BeadArray. Initial screening revealed successful genotyping for 93.4% of WGA samples and 99.3% of genomic samples, and 93.2% of SNPs for WGA samples and 96.3% of SNPs for genomic samples. SNP failure was predicted by Illumina-provided designability rank, %GC (P ≤ 0.002), and for WGA only, distance to telomere and Illumina-provided SNP score (P ≤ 0.002). Distance to telomere and %GC were highly correlated; adjustment for %GC removed the association between distance to telomere and SNP failure. Although universally high, per-SNP call rates were related to designability rank, SNP score, %GC, minor allele frequency, distance to telomere (P ≤ 0.01), and, for WGA only, Illumina-provided validation class (P < 0.001). We found excellent concordance generally (>99.0%) among 124 WGA:genomic replicates, 15 WGA replicates, 88 replicate aliquots of the same WGA preparation, and 25 genomic replicates. Where there was discordance, it was across WGA:genomic replicates but limited to only a few samples among other replicates suggesting the introduction of error. Designability rank and SNP score correlated with WGA:genomic concordance (P < 0.001). In summary, use of MDA WGA DNA is feasible; however, caution is warranted regarding SNP selection and analysis. We recommend that biological SNP characteristics, notably distance to telomere and GC content (<50% GC recommended), as well as Illumina-provided metrics be considered in the creation of GoldenGate assays using MDA WGA DNA. (Cancer Epidemiol Biomarkers Prev 2008;17(7):1781–9)