Critical DNA binding interactions of the insulator protein CTCF: a small number of zinc fingers mediate strong binding, and a single finger-DNA interaction controls binding at imprinted loci

The DNA-binding protein CTCF (CCCTC binding factor) mediates enhancer blocking insulation at sites throughout the genome and plays an important role in regulating allele-specific expression at the Igf2/H19 locus and at other imprinted loci. Evidence is also accumulating that CTCF is involved in large scale organization of genomic chromatin. Although CTCF has 11 zinc fingers, we show here that only 4 of these are essential to strong binding and that they recognize a core 12-bp DNA sequence common to most CTCF sites. By deleting individual fingers and mutating individual sites, we determined the orientation of binding. Furthermore, we were able to identify the specific finger and its point of DNA interaction that are responsible for the loss of CTCF binding when CpG residues are methylated in the imprinted Igf2/H19 locus. This single interaction appears to be critical for allele-specific binding and insulation by CTCF.     

CTCF-dependent chromatin bias constitutes transient epigenetic memory of the mother at the H19-Igf2 imprinting control region in prospermatogonia

Genomic imprints-parental allele-specific DNA methylation marks at the differentially methylated regions (DMRs) of imprinted genes-are erased and reestablished in germ cells according to the individual's sex. Imprint establishment at paternally methylated germ line DMRs occurs in fetal male germ cells. In prospermatogonia, the two unmethylated alleles exhibit different rates of de novo methylation at the H19/Igf2 imprinting control region (ICR) depending on parental origin. We investigated the nature of this epigenetic memory using bisulfite sequencing and allele-specific ChIP-SNuPE assays. We found that the chromatin composition in fetal germ cells was biased at the ICR between the two alleles with the maternally inherited allele exhibiting more H3K4me3 and less H3K9me3 than the paternally inherited allele. We determined genetically that the chromatin bias, and also the delayed methylation establishment in the maternal allele, depended on functional CTCF insulator binding sites in the ICR. Our data suggest that, in primordial germ cells, maternally inherited allele-specific CTCF binding sets up allele-specific chromatin differences at the ICR. The erasure of these allele-specific chromatin marks is not complete before the process of de novo methylation imprint establishment begins. CTCF-dependent allele-specific chromatin composition imposes a maternal allele-specific delay on de novo methylation imprint establishment at the H19/Igf2 ICR in prospermatogonia.     

Genome-wide identification of allele-specific effects on gene expression for single and multiple individuals

The analysis of allele-specific gene expression (ASE) is essential for the mapping of genetic variants that affect gene regulation, and for the identification of alleles that modify disease risk. Although RNA sequencing offers the opportunity to measure expression at allele levels, the availability of powerful statistical methods for mapping ASE in single or multiple individuals is limited. We developed a maximum likelihood model to characterize ASE in the human genome. Approximately 17% of genes displayed an allele-specific effect on gene expression in a single individual. Simulations using our model gave a better performance and improved robustness when compared with the binomial test, with different coverage levels, allelic expression fractions and random noise. In addition, our method can identify ASE in multiple individuals, with enhanced performance. This is helpful in understanding the mechanism of genetic regulation leading to expression changes, alternative splicing variants and even disease susceptibility.     

Refinement of single-nucleotide polymorphism genotyping methods on human genomic DNA: amplifluor allele-specific polymerase chain reaction versus ligation detection reaction-TaqMan

Single-nucleotide polymorphisms (SNPs) have proven to be powerful genetic markers for a variety of genetic applications, e.g., association studies leading to dissection of both monogenetic and complex diseases. However, no single SNP genotyping method has been broadly accepted. In the present study, we compared and refined two promising methods with potential for research and for diagnostic SNP genotyping: Amplifluor allele-specific polymerase chain reaction (PCR) and ligation detection reaction (LDR)-TaqMan. The methods are based on allele-specific primer extension and allele-specific ligation, respectively. Since LDR-TaqMan had previously been tested on just Arabidopsis thaliana, we adjusted the method for the more complex human genome. Amplifluor allele-specific PCR has a single-step and closed-tube format, whereas the LDR-TaqMan assay comprises two simple steps. Contrary to the primer-extension-based method, the ligation-based method can be multiplexed. Refining the LDR-TaqMan technique, we successfully replaced a previously suggested three-step multiplexing procedure with a less laborious two-step approach. Comparing refined LDR-TaqMan with Amplifluor allele-specific PCR in a family-based study, both techniques appeared similar with respect to high robustness and accuracy. As both approaches utilize primers with common tails, all SNPs can be assayed with the same couple of fluorescence reporting reagents, ensuring low establishing and running expenses.     

Detection of single-nucleotide polymorphisms in the <Emphasis Type="Italic">p53</Emphasis> gene by LDR/RCA in hydrogel microarrays

To find single-nucleotide polymorphisms (SNPs) in the human genome, three modern technologies of molecular genetic analysis were combined: the ligase detection reaction (LDR), rolling circle amplification (RCA), and immobilized microarray of gel elements (IMAGE). SNPs were detected in target DNA by selective ligation of allele-specific nucleotides in microarrays. The ligation product was assayed in microarray gel pads by RCA. Two variants of microarray analysis were compared. One included selective ligation of short oligonu-cleotides immobilized in a microarray with subsequent amplification with a preformed circular probe (a common circle). The probe was especially designed for human genome research. The other variant employed immobilized allele-specific padlock probes, which could be circularized as a result of selective ligation. Codon 72 SNP of the human p53 gene was used as a model. RCA in microarrays proved to be a quantitative assay and, in combination with LDR, allowed efficient discrimination of alleles. The principles and prospects of LDR/RCA in microarrays are discussed.Translated from Molekulyarnaya Biologiya, Vol. 39, No. 1, 2005, pp. 30–39.Original Russian Text Copyright © 2005 by Kashkin, Strizhkov, Gryadunov, Surzhikov, Grechishnikova, Kreindlin, Chupeeva, Evseev, Turygin, Mirzabekov.     

Analysis of allelic differential expression in the human genome using allele-specific serial analysis of gene expression tags

Recent reports have demonstrated that a significant proportion of human genes display allelic differential expression (ADE). ADE is associated with phenotypic variability and may contribute to complex genetic diseases. Here, we present a computational analysis of ADE using allele-specific serial analysis of gene expression (SAGE) tags representing 1295 human genes. We identified 472 genes for which unequal representation (>3-fold) of allele-specific SAGE tags was observed in at least one SAGE library, suggesting the occurrence of ADE. For 235 out of these 472 genes, the difference in the expression level between both allele-specific SAGE tags was statistically significant (p < 0.05). Eleven candidate genes were then subjected to experimental validation and ADE was confirmed for 8 out of these 11 genes. Our results suggest that at least 25% of the human genes display ADE and that allele-specific SAGE tags can be efficiently used for the identification of such genes.     

Detection of single base polymorphism in p53 gene by ligase detection reaction and rolling circle amplification on microarrays

We combined three modern technologies of single base polymorphism detection in human genome: ligase detection reaction, rolling circle amplification and IMAGE hydro-gel microarrays. Polymorphism in target DNA was tested by selective ligation on microarray. Product of the ligase reaction was determined in microarray gel pads by rolling circle amplification. Two different methods were compared. In first, selective ligation of short oligonucleotides immobilized on microarray was used with subsequent amplification on preformed circle probe ("common circle"). The circle probe was designed especially for human genome research. In second variant, allele-specific padlock probes that may be circularized by selective ligation were immobilized on microarray. Polymorphism of codon 72 in human p53 gene was used as a biological model. It was shown that LDR/RCA on microarray is a quantitative reaction and gives high discrimination of alleles. Principles and perspectives of selective ligation and rolling circle amplification are being discussed.     

ASHIC: hierarchical Bayesian modeling of diploid chromatin contacts and structures

The recently developed Hi-C technique has been widely applied to map genome-wide chromatin interactions. However, current methods for analyzing diploid Hi-C data cannot fully distinguish between homologous chromosomes. Consequently, the existing diploid Hi-C analyses are based on sparse and inaccurate allele-specific contact matrices, which might lead to incorrect modeling of diploid genome architecture. Here we present ASHIC, a hierarchical Bayesian framework to model allele-specific chromatin organizations in diploid genomes. We developed two models under the Bayesian framework: the Poisson-multinomial (ASHIC-PM) model and the zero-inflated Poisson-multinomial (ASHIC-ZIPM) model. The proposed ASHIC methods impute allele-specific contact maps from diploid Hi-C data and simultaneously infer allelic 3D structures. Through simulation studies, we demonstrated that ASHIC methods outperformed existing approaches, especially under low coverage and low SNP density conditions. Additionally, in the analyses of diploid Hi-C datasets in mouse and human, our ASHIC-ZIPM method produced fine-resolution diploid chromatin maps and 3D structures and provided insights into the allelic chromatin organizations and functions. To summarize, our work provides a statistically rigorous framework for investigating fine-scale allele-specific chromatin conformations. The ASHIC software is publicly available at https://github.com/wmalab/ASHIC.     

Allele-specific amplification in cancer revealed by SNP array analysis

Amplification, deletion, and loss of heterozygosity of genomic DNA are hallmarks of cancer. In recent years a variety of studies have emerged measuring total chromosomal copy number at increasingly high resolution. Similarly, loss-of-heterozygosity events have been finely mapped using high-throughput genotyping technologies. We have developed a probe-level allele-specific quantitation procedure that extracts both copy number and allelotype information from single nucleotide polymorphism (SNP) array data to arrive at allele-specific copy number across the genome. Our approach applies an expectation-maximization algorithm to a model derived from a novel classification of SNP array probes. This method is the first to our knowledge that is able to (a) determine the generalized genotype of aberrant samples at each SNP site (e.g., CCCCT at an amplified site), and (b) infer the copy number of each parental chromosome across the genome. With this method, we are able to determine not just where amplifications and deletions occur, but also the haplotype of the region being amplified or deleted. The merit of our model and general approach is demonstrated by very precise genotyping of normal samples, and our allele-specific copy number inferences are validated using PCR experiments. Applying our method to a collection of lung cancer samples, we are able to conclude that amplification is essentially monoallelic, as would be expected under the mechanisms currently believed responsible for gene amplification. This suggests that a specific parental chromosome may be targeted for amplification, whether because of germ line or somatic variation. An R software package containing the methods described in this paper is freely available at http://genome.dfci.harvard.edu/~tlaframb/PLASQ.     

CTCF elements direct allele-specific undermethylation at the imprinted H19 locus

The H19 imprinted gene locus is regulated by an upstream 2 kb imprinting control region (ICR) that influences allele-specific expression, DNA methylation, and replication timing. This ICR becomes de novo methylated during late spermatogenesis in the male but emerges from oogenesis in an unmethylated form, and this allele-specific pattern is then maintained throughout early development and in all tissues of the mouse. We have used a genetic approach involving transfection into embryonic stem (ES) cells in order to decipher how the maternal allele is protected from de novo methylation at the time of implantation. Our studies show that CCCTC binding factor (CTCF) boundary elements within the ICR have the ability to prevent de novo methylation on the maternal allele. Since CTCF does not recognize its binding sequence when methylated, this reaction does not occur on the paternal allele, thus preserving the gamete-derived, allele-specific pattern. These results suggest that CTCF may play a general role in the maintenance of differential methylation patterns in vivo.     

Identification of apolipoprotein E polymorphism by using synthetic oligonucleotides

A procedure based on selective hybridization with allele-specific oligonucleotides was developed for typing apolipoprotein E variants from human genomic DNAs. Two sets of oligonucleotides were synthesized and used to discriminate either between epsilon 3 and epsilon 4 alleles or between epsilon 3 and epsilon 2 alleles. Combination of the allele-specific oligonucleotide hybridization with the method for in vitro DNA amplification (Polymerase Chain Reaction) Saiki, R. K. et al. 1985. Science, 230: 1350-1354) (1) dramatically improved the sensitivity and the reliability of the procedure. Adaptation of a simple strategy involving direct cloning and DNA sequencing of in vitro amplified DNA enables rapid identification of any mutation within the apoE gene area encoding the receptor binding domain.     

Comparative Anatomy of Chromosomal Domains with Imprinted and Non-Imprinted Allele-Specific DNA Methylation

Allele-specific DNA methylation (ASM) is well studied in imprinted domains, but this type of epigenetic asymmetry is actually found more commonly at non-imprinted loci, where the ASM is dictated not by parent-of-origin but instead by the local haplotype. We identified loci with strong ASM in human tissues from methylation-sensitive SNP array data. Two index regions (bisulfite PCR amplicons), one between the C3orf27 and RPN1 genes in chromosome band 3q21 and the other near the VTRNA2-1 vault RNA in band 5q31, proved to be new examples of imprinted DMRs (maternal alleles methylated) while a third, between STEAP3 and C2orf76 in chromosome band 2q14, showed non-imprinted haplotype-dependent ASM. Using long-read bisulfite sequencing (bis-seq) in 8 human tissues we found that in all 3 domains the ASM is restricted to single differentially methylated regions (DMRs), each less than 2kb. The ASM in the C3orf27-RPN1 intergenic region was placenta-specific and associated with allele-specific expression of a long non-coding RNA. Strikingly, the discrete DMRs in all 3 regions overlap with binding sites for the insulator protein CTCF, which we found selectively bound to the unmethylated allele of the STEAP3-C2orf76 DMR. Methylation mapping in two additional genes with non-imprinted haplotype-dependent ASM, ELK3 and CYP2A7, showed that the CYP2A7 DMR also overlaps a CTCF site. Thus, two features of imprinted domains, highly localized DMRs and allele-specific insulator occupancy by CTCF, can also be found in chromosomal domains with non-imprinted ASM. Arguing for biological importance, our analysis of published whole genome bis-seq data from hES cells revealed multiple genome-wide association study (GWAS) peaks near CTCF binding sites with ASM.