Extensive variation between tissues in allele specific expression in an outbred mammal.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations.

Maternally derived Dnmt1o cytosine methyltransferase maintains parental, sex-specific genomic imprints during preimplantation development, an event that is postulated to occur at the 8-cell stage, the only time when Dnmt1o is normally found in the nucleus. Midgestation embryos developing in the absence of Dnmt1o have significant delays in development and a wide range of morphological abnormalities. The goal of this study was to determine the molecular imprinting defects in embryos and placentae of Dnmt1o-deficient mothers and better define the timing of loss of imprinting during preimplantation development. Females lacking oocyte-derived Dnmt1o as well as control females were crossed with CAST-7 males to generate Dnmt1o-deficient and control embryos, so that the parental alleles could be distinguished. Embryos and placentae were collected from Dnmt1o-deficient and control mice on 9.5 dpc and morulae were collected on 3.5 dpc. Bisulfite sequencing was used to determine the methylation of CpG sites in the 5′ regions of H19 and Snrpn and allele specific gene expression was assessed by fluorescent hybridization probe analysis. At 9.5 dpc, control embryos showed normal methylation and expression: H19 was methylated on the paternal allele and expressed exclusively from the maternal allele while Snrpn was methylated on the maternal allele and expressed from the paternal allele. In contrast, control placentae showed a loss of methylation imprints (loss of 39% for H19 and 35% for Snrpn) while maintaining normal appropriately mono-allelic expression for the two genes. Embryos from Dnmt1o-deficient mothers showed three different patterns of methylation for both genes: completely normal, 50% loss and complete loss of methylation, with variable levels of biallelic expression. Methylation patterns in placentae of conceptuses from Dnmt1o deficient mice were similar to those in controls, even though the expression became biallelic for both genes. Together, these results provide evidence that methylation is not required to maintain monoallelic expression of H19 and Snrpn in the placentae and thus was not informative in assessing disruption of imprinting in the Dnmt1o model. Methylation status was also investigated in pooled morulae (n=8) from control and Dnmt1o-deficient mothers, in order to determine the timing of the initial methylation defect. Consistent with DNA methylation being lost after the 8-cell stage, methylation of H19 and Snrpn were reduced by about 50% in morulae from Dnmt1o deficient mothers. The variation observed in the phenotype and imprinted gene methylation/expression appears to be due to the epigenetic mosaic make-up of embryos created by the lack of Dnmt1o protein in the preimplantation period. These findings have important implications for understanding the epigenetic controls of imprinted gene expression and the types of embryonic phenotypes related to the disruption of inherited imprints. (Supported by the NIH and CIHR) (poster)

Genome-Wide Analysis of Allele-Specific Expression Genes in Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia

Chickens are a valuable livestock species, providing an affordable and nutritious food source worldwide through eggs and meat. Additionally, chickens are extensively used in scientific research as model organisms in virology, immunology, epigenetics, development, and conservation biology. Considering this, there are increasing efforts to deepen our understanding of their genome complexity, such as efforts by the Functional Annotation of Animal Genomes (FAANG) Consortium. Toward this goal, we are working to contribute to the annotation of the chicken transcriptome. Allele-specific expression (ASE) is an imbalance of allelic gene expression often driven by genetic and epigenetic changes in cis-regulatory regions. We, therefore, aimed to determine allelic imbalances in gene expression across 20 tissues and cells to shed light on mechanisms regulating gene expression. Using replicate paired-end (PE) RNA-seq samples, we identified around 7 million variants across 20 tissues and cells (e.g., reproductive tissues, intestine tissues, muscle, and immune tissues and cells). We applied quality control measures and filtered out monoallelic and homozygous variants. After that, we calculated read counts per allele to determine allelic expression. We found 365,894 significant ASE SNPs and 11,530 ASE genes in at least one tissue or cell type (FDR ≤ 0.05). We discovered that ASE SNPs are widely distributed throughout the chicken genome, and ASE SNPs and ASE genes showed a varied distribution across tissues and cells. Primary macrophage exhibited the most abundant, while the pectoralis major (breast muscle) had the lowest ASE genes and ASE SNPs, ranging from 7,592 to 32,505 ASE SNPs and 773 to 2,387 ASE genes. Our findings reveal several important pathways affected by allelic imbalance of expression. These pathways included fatty acid biosynthesis and regulation, cell proliferation and differentiation, cell metabolism, adipocytokine signaling, proteolysis and autophagy, and focal adhesion. In summary, our study provides insight into relevant biological processes critical to chickens and can contribute to improving genome annotation, helping to expand the transcriptomic reference source.

Most loci that are regulated by genomic imprinting have differentially methylated regions (DMRs). Previously, we showed that the DMRs of the mouse Snrpn and U2af1-rs1 genes have paternal allele-specific patterns of acetylation on histones H3 and H4. To investigate the maintenance of acetylation at these DMRs, we performed chromatin immunoprecipitation on trichostatin-A (TSA)-treated and control cells. In embryonic stem (ES) cells and fibroblasts, brief (6-h) TSA treatment induces global hyperacetylation of H3 and H4. In ES cells only, TSA led to a selective increase in maternal acetylation at U2af1-rs1, at lysine 5 of H4 and at lysine 14 of H3. TSA treatment of ES cells did not affect DNA methylation or expression of U2af1-rs1, but was sufficient to increase DNase I sensitivity along the maternal allele to a level comparable with that of the paternal allele. In fibroblasts, TSA did not alter U2af1-rs1 acetylation, and the parental alleles retained their differential DNase I sensitivity. At Snrpn, no changes in acetylation were observed in the TSA-treated cells. Our data suggest that the mechanisms regulating histone acetylation at DMRs are locus and developmental stage-specific and are distinct from those effecting global levels of acetylation. Furthermore, it seems that the allelic U2af1-rs1 acetylation determines DNase I sensitivity/chromatin conformation.

Cytokinins play important roles in the growth and development of plants. Physiological and photosynthetic characteristics are common indicators to measure the growth and development in plants. However, few reports have described the molecular mechanisms of physiological and photosynthetic changes in response to cytokinin, particularly in woody plants. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression in response to the external environment. In this study, we examined genome-wide DNA methylation variation and transcriptional variation in poplar (Populus tomentosa) after short-term treatment with the synthetic cytokinin 6-benzylaminopurine (6-BA). We identified 460 significantly differentially methylated regions (DMRs) in response to 6-BA treatment. Transcriptome analysis showed that 339 protein-coding genes, 262 long non-coding RNAs (lncRNAs), and 15,793 24-nt small interfering RNAs (siRNAs) were differentially expressed under 6-BA treatment. Among these, 79% were differentially expressed between alleles in P. tomentosa, and 102,819 allele-specific expression (ASE) loci in 19,200 genes were detected showing differences in ASE levels after 6-BA treatment. Combined DNA methylation and gene expression analysis demonstrated that DNA methylation plays an important role in regulating allele-specific gene expression. To further investigate the relationship between these 6-BA-responsive genes and phenotypic variation, we performed SNP analysis of 460 6-BA-responsive DMRs via re-sequencing using a natural population of P. tomentosa, and we identified 206 SNPs that were significantly associated with growth and wood properties. Association analysis indicated that 53% of loci with allele-specific expression had primarily dominant effects on poplar traits. Our comprehensive analyses of P. tomentosa DNA methylation and the regulation of allele-specific gene expression suggest that DNA methylation is an important regulator of imbalanced expression between allelic loci.

Comprehensive analyses of 435 goat transcriptomes provides insight into male reproduction

Allelic specific gene expression (ASGE) appears to be an important factor in human phenotypic variability and as a consequence, for the development of complex traits and diseases. In order to study ASGE across the human genome, we have performed a study in which genotyping was coupled with an analysis of ASGE by screening 11,500 SNPs using the Mapping 10 K Array to identify differential allelic expression. We found that from the 5,133 SNPs that were suitable for analysis (heterozygous in our sample and expressed in peripheral blood mononuclear cells), 2,934 (57%) SNPs had differential allelic expression. Such SNPs were equally distributed along human chromosomes and biological processes. We validated the presence or absence of ASGE in 18 out 20 SNPs (90%) randomly selected by real time PCR in 48 human subjects. In addition, we observed that SNPs close to -but not included in- segmental duplications had increased levels of ASGE. Finally, we found that transcripts of unknown function or non-coding RNAs, also display ASGE: from a total of 2,308 intronic SNPs, 1510 (65%) SNPs underwent differential allelic expression. In summary, ASGE is a widespread mechanism in the human genome whose regulation seems to be far more complex than expected.

Different sheep breeds show distinct phenotypic plasticity in fat deposition in the tails. The genetic background underlying fat deposition in the tail of sheep is complex, multifactorial, and may involve allele-specific expression (ASE) mechanism to modulate allelic expression. ASE is a common phenomenon in mammals and refers to allelic imbalanced expression modified by cis-regulatory genetic variants that can be observed at heterozygous loci. Therefore, regulatory processes behind the fat-tail formation in sheep may be to some extent explained by cis- regulatory variants, through ASE mechanism, which was investigated in the present study. An RNA-Seq-based variant calling was applied to perform genome-wide survey of ASE genes using 45 samples from seven independent studies comparing the transcriptome of fat-tail tissue between fat- and thin-tailed sheep breeds. Using a rigorous computational pipeline, 115 differential ASE genes were identified, which were narrowed down to four genes (LPL, SOD3, TCP1 and LRPAP1) for being detected in at least two studies. Functional analysis revealed that the ASE genes were mainly involved in fat metabolism. Of these, LPL was of greater importance, as 1) observed in five studies, 2) reported as ASE gene in the previous studies and 3) with a known role in fat deposition. Our findings implied that complex physiological traits, like fat-tail formation, can be better explained by considering various genetic mechanisms, which can be more finely mapped through ASE analyses. The insights gained in this study indicate that biallelic expression may not be a common mechanism in sheep fat-tail development. Hence, allelic imbalance of the fat deposition-related genes can be considered a novel layer of information for future research on genetic improvement and increased efficiency in sheep breeding programs.

In mammals, perception of pheromones is based on the expression in each vomeronasal sensory neuron of a limited set of receptor genes, chosen among a large repertoire. Here, we report an extremely tight control of the monogenic and monoallelic transcription of the V1rb2 receptor gene. Combining genetic and electrophysiological approaches, we show that the transcription of a non-functional V1r allele leads to the coexpression of another, functional V1r gene. The choice of this coexpressed gene surprisingly includes genes located on the cluster homologous to the one from which the mutant allele is transcribed. However, V1r genes located in cis relative to the transcribed mutant allele are excluded from the coexpression choice. Our observations strongly suggest a monogenic regulatory mechanism acting (a) at a general level, via the expression of the V1r receptor itself, and (b) at a more local level, defined by the V1r gene cluster.

SNP-based large-scale identification of allele-specific gene expression in human B cells

Recent large-scale genomic studies established the occurrence of multiple DNA sequence variants in genomes of healthy individuals that differ from the reference sequence. Among these variants mostly represented by germline single nucleotide polymorphisms disease-related alleles are detected including alleles which are associated with monogenic disorders, and putative deleterious genetic variants. Apart from functional significance of a particular variant and of a gene harboring it, the penetrance of these allelic variants depends on their expression level and can be determined by preferential expression of a particular allele, or allele-specific expression. It is estimated that 20–30 % of genes present in the human genome display allelic bias in a tissue-specific manner. Allele-specific expression is defined by a range of genetic and epigenetic mechanisms including cis-regulatory polymorphisms, allele-specific binding of transcription factors, allele-specific DNA methylation and regulation through non-coding RNA. Although the data on the issue are scarce, allele-specific expression has been reported to be implicated in several hereditary disorders including benign and malignant tumors of the large intestine. Recent studies that estimate allele-specific expression incidence in tumors and identify wide range of genes displaying allelic imbalance indicate that allele-specific expression might play a significant role in carcinogenesis. Eventually, estimation of transcriptional rate of allelic variants which cause dysfunction of oncogenes and tumor suppressors may prove to be essential for rational choice of antitumor therapeutic strategy. In this review, we outline the main concepts and mechanisms of allele-specific expression and the data on allelic imbalance in tumors.

Large-scale profiling and identification of potential regulatory mechanisms for allelic gene expression in colorectal cancer cells

A large number of genome-wide association studies have been performed during the past five years to identify associations between SNPs and human complex diseases and traits. The assignment of a functional role for the identified disease-associated SNP is not straight-forward. Genome-wide expression quantitative trait locus (eQTL) analysis is frequently used as the initial step to define a function while allele-specific gene expression (ASE) analysis has not yet gained a wide-spread use in disease mapping studies. We compared the power to identify cis-acting regulatory SNPs (cis-rSNPs) by genome-wide allele-specific gene expression (ASE) analysis with that of traditional expression quantitative trait locus (eQTL) mapping. Our study included 395 healthy blood donors for whom global gene expression profiles in circulating monocytes were determined by Illumina BeadArrays. ASE was assessed in a subset of these monocytes from 188 donors by quantitative genotyping of mRNA using a genome-wide panel of SNP markers. The performance of the two methods for detecting cis-rSNPs was evaluated by comparing associations between SNP genotypes and gene expression levels in sample sets of varying size. We found that up to 8-fold more samples are required for eQTL mapping to reach the same statistical power as that obtained by ASE analysis for the same rSNPs. The performance of ASE is insensitive to SNPs with low minor allele frequencies and detects a larger number of significantly associated rSNPs using the same sample size as eQTL mapping. An unequivocal conclusion from our comparison is that ASE analysis is more sensitive for detecting cis-rSNPs than standard eQTL mapping. Our study shows the potential of ASE mapping in tissue samples and primary cells which are difficult to obtain in large numbers.

In diploid eukaryotic organisms, both alleles of each autosomal gene are usually assumed to be simultaneously expressed at similar levels. However, some genes can be expressed preferentially or strictly from a single allele, a process known as monoallelic expression. Classic monoallelic expression of X-chromosome-linked genes, olfactory receptor genes and developmentally imprinted genes is the result of epigenetic modifications. Genetic-origin-dependent monoallelic expression, however, is caused by cis-regulatory differences between the alleles. There is a paucity of systematic study to investigate these phenomena across multiple tissues, and the mechanisms underlying such monoallelic expression are not yet fully understood. Here we provide a detailed portrait of monoallelic gene expression across multiple tissues/cell lines in a hybrid mouse cross between the Mus musculus strain C57BL/6J and the Mus spretus strain SPRET/EiJ. We observed pervasive tissue-dependent allele-specific gene expression: in total, 1,839 genes exhibited monoallelic expression in at least one tissue, and 410 genes in at least two tissues. Among these 88 are monoallelic genes with different active alleles between tissues, probably representing genetic-origin-dependent monoallelic expression. We also identified six autosomal monoallelic genes with the active allele being identical in all eight tissues, which are likely novel candidates of imprinted genes. To depict the underlying regulatory mechanisms at the chromatin layer, we performed ATAC-seq in two different cell lines derived from the F1 mouse. Consistent with the global expression pattern, cell-type dependent monoallelic peaks were found, and a higher proportion of C57BL/6J-active peaks were observed in both cell types, implying possible species-specific regulation. Finally, only a small part of monoallelic gene expression could be explained by allelic differences in chromatin organization in promoter regions, suggesting that other distal elements may play important roles in shaping the patterns of allelic gene expression across tissues.