SummaryCharacterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

CblX (MIM309541) is an X-linked recessive disorder characterized by defects in cobalamin (vitamin B12) metabolism and other developmental defects. Mutations in HCFC1, a transcriptional co-regulator which interacts with multiple transcription factors, have been associated with cblX. HCFC1 regulates cobalamin metabolism via the regulation of MMACHC expression through its interaction with THAP11, a THAP domain-containing transcription factor. The HCFC1/THAP11 complex potentially regulates genes involved in diverse cellular functions including cell cycle, proliferation, and transcription. Thus, it is likely that mutation of THAP11 also results in biochemical and other phenotypes similar to those observed in patients with cblX. We report a patient who presented with clinical and biochemical phenotypic features that overlap cblX, but who does not have any mutations in either MMACHC or HCFC1. We sequenced THAP11 by Sanger sequencing and discovered a potentially pathogenic, homozygous variant, c.240C > G (p.Phe80Leu). Functional analysis in the developing zebrafish embryo demonstrated that both THAP11 and HCFC1 regulate the proliferation and differentiation of neural precursors, suggesting important roles in normal brain development. The loss of THAP11 in zebrafish embryos results in craniofacial abnormalities including the complete loss of Meckel's cartilage, the ceratohyal, and all of the ceratobranchial cartilages. These data are consistent with our previous work that demonstrated a role for HCFC1 in vertebrate craniofacial development. High throughput RNA-sequencing analysis reveals several overlapping gene targets of HCFC1 and THAP11. Thus, both HCFC1 and THAP11 play important roles in the regulation of cobalamin metabolism as well as other pathways involved in early vertebrate development.

B-cell precursor acute lymphoblastic leukemia (pre-B ALL) is the most common pediatric cancer. Although the genetic determinants underlying disease onset remain unclear, epigenetic modifications including DNA methylation are suggested to contribute significantly to leukemogenesis. Using the Illumina 450K array, we assessed DNA methylation in matched tumor-normal samples of 46 childhood patients with pre-B ALL, extending single CpG-site resolution analysis of the pre-B ALL methylome beyond CpG-islands (CGI). Unsupervised hierarchical clustering of CpG-site neighborhood, gene, or microRNA (miRNA) gene-associated methylation levels separated the tumor cohort according to major pre-B ALL subtypes, and methylation in CGIs, CGI shores, and in regions around the transcription start site was found to significantly correlate with transcript expression. Focusing on samples carrying the t(12;21) ETV6-RUNX1 fusion, we identified 119 subtype-specific high-confidence marker CpG-loci. Pathway analyses linked the CpG-loci-associated genes with hematopoiesis and cancer. Further integration with whole-transcriptome data showed the effects of methylation on expression of 17 potential drivers of leukemogenesis. Independent validation of array methylation and sequencing-derived transcript expression with Sequenom Epityper technology and realtime quantitative reverse transcriptase PCR, respectively, indicates more than 80% empirical accuracy of our genome-wide findings. In summary, genome-wide DNA methylation profiling enabled us to separate pre-B ALL according to major subtypes, to map epigenetic biomarkers specific for the t(12;21) subtype, and through a combined methylome and transcriptome approach to identify downstream effects on candidate drivers of leukemogenesis. Cancer Res; 73(14); 4323-36. Ó2013 AACR.