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.

Cell culture-induced aberrant methylation of the imprinted IG DMR in human lymphoblastoid cell lines.

DNA methylation patterns are often poorly conserved through cell culturing. To determine the effect of cell immortalization and culture on DNA methylation profiles, we analyzed methylation in the differentially methylated regions (DMR) of five imprinted domains: the intergenic (IG) DMR on chromosome 14q32; potassium voltage-gated channel, KQT-like subfamily, member 1, (KCNQ1); small nuclear ribonucleoprotein polypeptide N (SNRPN), mesoderm specific transcript homolog (MEST); and H19 in lymphoblastoid cell lines (LCLs). In the IG DMR we found an aberrant methylation pattern that was consistent through all the cell lines tested and significantly different from that of noncultured peripheral blood cells. Using a generalized linear mixed model to compare methylation profiles, we show that recently derived LCLs significantly differ from the CEPH LCLs. This implies a gradual cell-culture related deterioration of DNA methylation in the IG DMR with at least two steps that may be identified: loss of methylation at CG sites 1 and 8; and loss of allelic differences in DNA methylation. The IG DMR methylation profile also confirms the high level of clonality of the CEPH LCLs. We conclude that non-transformed primary cells may be less susceptible to epigenetic anomalies and therefore may provide a more accurate reflection of gene expression in vivo.

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.