DNA methylation profiles of human active and inactive X chromosomes.

X-chromosome inactivation (XCI) is a dosage compensation mechanism that silences the majority of genes on one X chromosome in each female cell. To characterize epigenetic changes that accompany this process, we measured DNA methylation levels in 45,X patients carrying a single active X chromosome (X(a)), and in normal females, who carry one X(a) and one inactive X (X(i)). Methylated DNA was immunoprecipitated and hybridized to high-density oligonucleotide arrays covering the X chromosome, generating epigenetic profiles of active and inactive X chromosomes. We observed that XCI is accompanied by changes in DNA methylation specifically at CpG islands (CGIs). While the majority of CGIs show increased methylation levels on the X(i), XCI actually results in significant reductions in methylation at 7% of CGIs. Both intra- and inter-genic CGIs undergo epigenetic modification, with the biggest increase in methylation occurring at the promoters of genes silenced by XCI. In contrast, genes escaping XCI generally have low levels of promoter methylation, while genes that show inter-individual variation in silencing show intermediate increases in methylation. Thus, promoter methylation and susceptibility to XCI are correlated. We also observed a global correlation between CGI methylation and the evolutionary age of X-chromosome strata, and that genes escaping XCI show increased methylation within gene bodies. We used our epigenetic map to predict 26 novel genes escaping XCI, and searched for parent-of-origin-specific methylation differences, but found no evidence to support imprinting on the human X chromosome. Our study provides a detailed analysis of the epigenetic profile of active and inactive X chromosomes.

Identification of cis- and trans-regulatory variation modulating microRNA expression levels in human fibroblasts.

MicroRNAs (miRNAs) are regulatory noncoding RNAs that affect the production of a significant fraction of human mRNAs via post-transcriptional regulation. Interindividual variation of the miRNA expression levels is likely to influence the expression of miRNA target genes and may therefore contribute to phenotypic differences in humans, including susceptibility to common disorders. The extent to which miRNA levels are genetically controlled is largely unknown. In this report, we assayed the expression levels of miRNAs in primary fibroblasts from 180 European newborns of the GenCord project and performed association analysis to identify eQTLs (expression quantitative traits loci). We detected robust expression for 121 miRNAs out of 365 interrogated. We have identified significant cis- (10%) and trans- (11%) eQTLs. Furthermore, we detected one genomic locus (rs1522653) that influences the expression levels of five miRNAs, thus unraveling a novel mechanism for coregulation of miRNA expression.

Methylation profiling in individuals with uniparental disomy identifies novel differentially methylated regions on chromosome 15.

The maternal and paternal genomes possess distinct epigenetic marks that distinguish them at imprinted loci. In order to identify imprinted loci, we used a novel method, taking advantage of the fact that uniparental disomy (UPD) provides a system that allows the two parental chromosomes to be studied independently. We profiled the paternal and maternal methylation on chromosome 15 using immunoprecipitation of methylated DNA and hybridization to tiling oligonucleotide arrays. Comparison of six individuals with maternal versus paternal UPD15 revealed 12 differentially methylated regions (DMRs). Putative DMRs were validated by bisulfite sequencing, confirming the presence of parent-of-origin-specific methylation marks. We detected DMRs associated with known imprinted genes within the Prader-Willi/Angelman syndrome region, such as SNRPN and MAGEL2, validating this as a method of detecting imprinted loci. Of the 12 DMRs identified, eight were novel, some of which are associated with genes not previously thought to be imprinted. These include a site within intron 2 of IGF1R at 15q26.3, a gene that plays a fundamental role in growth, and an intergenic site upstream of GABRG3 that lies within a previously defined candidate region conferring an increased maternal risk of psychosis. These data provide a map of parent-of-origin-specific epigenetic modifications on chromosome 15, identifying DNA elements that may play a functional role in the imprinting process. Application of this methodology to other chromosomes for which UPD has been reported will allow the systematic identification of imprinted sites throughout the genome.

Extensive natural variation for cellular hydrogen peroxide release is genetically controlled.

Natural variation in DNA sequence contributes to individual differences in quantitative traits. While multiple studies have shown genetic control over gene expression variation, few additional cellular traits have been investigated. Here, we investigated the natural variation of NADPH oxidase-dependent hydrogen peroxide (H(2)O(2) release), which is the joint effect of reactive oxygen species (ROS) production, superoxide metabolism and degradation, and is related to a number of human disorders. We assessed the normal variation of H(2)O(2) release in lymphoblastoid cell lines (LCL) in a family-based 3-generation cohort (CEPH-HapMap), and in 3 population-based cohorts (KORA, GenCord, HapMap). Substantial individual variation was observed, 45% of which were associated with heritability in the CEPH-HapMap cohort. We identified 2 genome-wide significant loci of Hsa12 and Hsa15 in genome-wide linkage analysis. Next, we performed genome-wide association study (GWAS) for the combined KORA-GenCord cohorts (n = 279) using enhanced marker resolution by imputation (>1.4 million SNPs). We found 5 significant associations (p<5.00×10-8) and 54 suggestive associations (p<1.00×10-5), one of which confirmed the linked region on Hsa15. To replicate our findings, we performed GWAS using 58 HapMap individuals and ∼2.1 million SNPs. We identified 40 genome-wide significant and 302 suggestive SNPs, and confirmed genome signals on Hsa1, Hsa12, and Hsa15. Genetic loci within 900 kb from the known candidate gene p67phox on Hsa1 were identified in GWAS in both cohorts. We did not find replication of SNPs across all cohorts, but replication within the same genomic region. Finally, a highly significant decrease in H(2)O(2) release was observed in Down Syndrome (DS) individuals (p<2.88×10-12). Taken together, our results show strong evidence of genetic control of H(2)O(2) in LCL of healthy and DS cohorts and suggest that cellular phenotypes, which themselves are also complex, may be used as proxies for dissection of complex disorders.

Transcriptional and post-transcriptional profile of human chromosome 21.

Recent studies have demonstrated extensive transcriptional activity across the human genome, a substantial fraction of which is not associated with any functional annotation. However, very little is known regarding the post-transcriptional processes that operate within the different classes of RNA molecules. To characterize the post-transcriptional properties of expressed sequences from human chromosome 21 (HSA21), we separated RNA molecules from three cell lines (GM06990, HeLa S3, and SK-N-AS) according to their ribosome content by sucrose gradient fractionation. Polyribosomal-associated RNA and total RNA were subsequently hybridized to genomic tiling arrays. We found that approximately 50% of the transcriptional signals were located outside of annotated exons and were considered as TARs (transcriptionally active regions). Although TARs were observed among polysome-associated RNAs, RT-PCR and RACE experiments revealed that approximately 40% were likely to represent nonspecific cross-hybridization artifacts. Bioinformatics discrimination of TARs according to conservation and sequence complexity allowed us to identify a set of high-confidence TARs. This set of TARs was significantly depleted in the polysomes, suggesting that it was not likely to be involved in translation. Analysis of polysome representation of RefSeq exons showed that at least 15% of RefSeq transcripts undergo significant post-transcriptional regulation in at least two of the three cell lines tested. Among the regulated transcripts, enrichment analysis revealed an over-representation of genes involved in Alzheimer's disease (AD), including APP and the BACE1 protease that cleaves APP to produce the pathogenic beta 42 peptide. We demonstrate that the combination of RNA fractionation and tiling arrays is a powerful method to assess the transcriptional and post-transcriptional properties of genomic regions.

Promoter polymorphisms and allelic imbalance in ABCB1 expression.

OBJECTIVE: The ABCB1 (MDR1) gene, encoding the transporter P-glycoprotein, is known to act on a broad range of prescription medicines. For this reason a large number of studies have assessed the functional consequences of variation in this gene. Particular attention has focused on the ABCB1_3435C>T polymorphism, an exonic variant resulting in a synonymous change. This variant has been associated with mRNA, protein and serum levels, and with responses to a number of medicines. The results of association studies have, however, been variable and it is not currently clear whether this polymorphism is functional or is in linkage disequilibrium with functionally distinct alleles. RESULTS: To identify functional variation in the ABCB1 gene we assessed allelic imbalance by pyrosequencing cDNA from 80 lymphoblastoid B cell lines from the Centre d'Etude du Polymorphisme Humain (CEPH) collection, including 74 individuals heterozygous for 3435C>T. We found that the degree of ABCB1 allelic imbalance differed among B-cell lines. In an effort to fine-map variants that influence the proportion of the two allelic mRNA species we genotyped representative common variations near the 3435C>T polymorphism by using a tagging single nucleotide polymorphism (SNP) approach. In one approach, we assessed in segregating families the impact of cis-acting variants on mRNA levels by using allelic imbalance as the phenotype in a regression analysis that distinguishes the coupling arrangements (phase) among alleles. In a second approach, we assessed allelic imbalance levels in lymphoblastoid B-cell lines from unrelated HapMap individuals, and performed an association using tagSNPs in a 5-Mb region surrounding the gene. Two potential cis-acting variants, a promoter rs28656907/rs28373093 dinucleotide polymorphism (P=0.007) and the rs10245483 SNP (P=0.0003) located 2 Mb upstream from the gene, were predictors of ABCB1 expression. CONCLUSIONS: The study outlines a general experimental approach for fine mapping gene variants that influence mRNA expression by using cultured cell lines.