Genome-wide identification of autosomal genes with allelic imbalance of chromatin state

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single-nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, and with Polycomb occupancy. We show that histone genes are expressed with allelic imbalance in mESCs, and are involved in haplotype-specific chromatin contacts marked by H3K27me3. Conditional knockouts of Polycomb enzymatic subunits, Ezh2 or Ring1, show that one-third of ASE genes, including histone genes, is regulated through Polycomb repression. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

Detection of DNA copy number alterations by matrix-assisted laser desorption/ionization time-of-flight mass spectrometric analysis of single nucleotide polymorphisms

OBJECTIVES

Copy number alterations (CNAs) are frequently found in malignant tissues. Different approaches have been used for CNA detection. However, it is not easy to detect a large panel of CNA targets in heterogenous tumors.

METHODS

We have developed a CNAs detection approach through quantitatively analyzed allelic imbalance by allelotyping single nucleotide polymorphisms (SNPs) by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Furthermore, the copy number changes were quantified by real-competitive PCR (rcPCR) to distinguish loss of heterozygosity (LOH) and genomic amplification. The approach was used to validate the CNA regions detected by next generation sequencing (NGS) in early-stage lung carcinoma.

RESULTS

CNAs were detected in heterogeneous DNA samples where tumor DNA is present at only 10% through the SNP based allelotyping. In addition, two different types of CNAs (loss of heterozygosity and chromosome amplification) were able to be distinguished quantitatively by rcPCR. Validation on a total of 41 SNPs from the selected CNA regions showed that copy number changes did occur, and the tissues from early-stage lung carcinoma were distinguished from normal.

CONCLUSIONS

CNA detection by MALDI-TOF MS can be used for validating potentially interesting genomic regions identified from next generation sequencing, and for detecting CNAs in tumor tissues consisting of a mixture of neoplastic and normal cells.

Asynchronous Replication Timing: A Mechanism for Monoallelic Choice During Development

Developmental programming is carried out by a sequence of molecular choices that epigenetically mark the genome to generate the stable cell types which make up the total organism. A number of important processes, such as genomic imprinting, selection of immune or olfactory receptors, and X-chromosome inactivation in females are dependent on the ability to stably choose one single allele in each cell. In this perspective, we propose that asynchronous replication timing (ASRT) serves as the basis for a sophisticated universal mechanism for mediating and maintaining these decisions.

Locus specific epigenetic modalities of random allelic expression imbalance

Most autosomal genes are thought to be expressed from both alleles, with some notable exceptions, including imprinted genes and genes showing random monoallelic expression (RME). The extent and nature of RME has been the subject of debate. Here we investigate the expression of several candidate RME genes in F1 hybrid mouse cells before and after differentiation, to define how they become persistently, monoallelically expressed. Clonal monoallelic expression is not present in embryonic stem cells, but we observe high frequencies of monoallelism in neuronal progenitor cells by assessing expression status in more than 200 clones. We uncover unforeseen modes of allelic expression that appear to be gene-specific and epigenetically regulated. This non-canonical allelic regulation has important implications for development and disease, including autosomal dominant disorders and opens up therapeutic perspectives.

Induced pluripotent stem cell reprogramming-associated methylation at the GABRA2 promoter and chr4p12 GABAA subunit gene expression in the context of alcohol use disorder

Twin studies indicate that there is a significant genetic contribution to the risk of developing alcohol use disorder (AUD). With the exception of coding variants in ADH1B and ALDH2, little is known about the molecular effects of AUD-associated loci. We previously reported that the AUD-associated synonymous polymorphism rs279858 within the GABAA α2 receptor subunit gene, GABRA2, was associated with gene expression of the chr4p12 GABAA subunit gene cluster in induced pluripotent stem cell (iPSC)-derived neural cultures. Based on this and other studies that showed changes in GABRA2 DNA methylation associated with schizophrenia and aging, we examined methylation in GABRA2. Specifically, using 69 iPSC lines and neural cultures derived from 47 of them, we examined whether GABRA2 rs279858 genotype predicted methylation levels and whether methylation was related to GABAA receptor subunit gene expression. We found that the GABRA2 CpG island undergoes random stochastic methylation during reprogramming and that methylation is associated with decreased GABRA2 gene expression, an effect that extends to the GABRB1 gene over 600 kb distal to GABRA2. Further, we identified additive effects of GABRA2 CpG methylation and GABRA2 rs279858 genotype on expression of the GABRB1 subunit gene in iPSC-derived neural cultures.

Maternal 5mCpG Imprints at the PARD6G-AS1 and GCSAML Differentially Methylated Regions Are Decoupled From Parent-of-Origin Expression Effects in Multiple Human Tissues

A hallmark of imprinted genes in mammals is the occurrence of parent-of-origin-dependent asymmetry of DNA cytosine methylation (5^mC) of alleles at CpG islands (CGIs) in their promoter regions. This 5^mCpG asymmetry between the parental alleles creates allele-specific imprinted differentially methylated regions (iDMRs). iDMRs are often coupled to the transcriptional repression of the methylated allele and the activation of the unmethylated allele in a tissue-specific, developmental-stage-specific and/or isoform-specific fashion. iDMRs function as regulatory platforms, built through the recruitment of chemical modifications to histones to achieve differential, parent-of-origin-dependent chromatin segmentation states. Here, we used a comparative computational data mining approach to identify 125 novel constitutive candidate iDMRs that integrate the maximal number of allele-specific methylation region records overlapping CGIs in human methylomes. Twenty-nine candidate iDMRs display gametic 5^mCpG asymmetry, and another 96 are candidate secondary iDMRs. We established the maternal origin of the 5^mCpG imprints of one gametic (PARD6G-AS1) and one secondary (GCSAML) iDMRs. We also found a constitutively hemimethylated, nonimprinted domain at the PWWP2AP1 promoter CGI with oocyte-derived methylation asymmetry. Given that the 5^mCpG level at the iDMRs is not a sufficient criterion to predict active or silent locus states and that iDMRs can regulate genes from a distance of more than 1 Mb, we used RNA-Seq experiments from the Genotype-Tissue Expression project and public archives to assess the transcriptional expression profiles of SNPs across 4.6 Mb spans around the novel maternal iDMRs. We showed that PARD6G-AS1 and GCSAML are expressed biallelically in multiple tissues. We found evidence of tissue-specific monoallelic expression of ZNF124 and OR2L13, located 363 kb upstream and 419 kb downstream, respectively, of the GCSAML iDMR. We hypothesize that the GCSAML iDMR regulates the tissue-specific, monoallelic expression of ZNF124 but not of OR2L13. We annotated the non-coding epigenomic marks in the two maternal iDMRs using data from the Roadmap Epigenomics project and showed that the PARD6G-AS1 and GCSAML iDMRs achieve contrasting activation and repression chromatin segmentations. Lastly, we found that the maternal 5^mCpG imprints are perturbed in several hematopoietic cancers. We conclude that the maternal 5^mCpG imprints at PARD6G-AS1 and GCSAML iDMRs are decoupled from parent-of-origin transcriptional expression effects in multiple tissues.