Autosomal monoallelic expression: genetics of epigenetic diversity?

The Spectrum of the Heterozygous Effect in Biallelic Mendelian Diseases-The Symptomatic Heterozygote Issue

Heterozygous carriers of pathogenic/likely pathogenic variants in autosomal recessive disorders seem to be asymptomatic. However, in recent years, an increasing number of case reports have suggested that mild and unspecific symptoms can occur in some heterozygotes, as symptomatic heterozygotes have been identified across different disease types, including neurological, neuromuscular, hematological, and pulmonary diseases. The symptoms are usually milder in heterozygotes than in biallelic variants and occur "later in life". The status of symptomatic heterozygotes as separate entities is often disputed, and alternative diagnoses are considered. Indeed, often only a thin line exists between dual, dominant, and recessive modes of inheritance and symptomatic heterozygosity. Interestingly, recent population studies have found global disease effects in heterozygous carriers of some genetic variants. What makes the few heterozygotes symptomatic, while the majority show no symptoms? The molecular basis of this phenomenon is still unknown. Possible explanations include undiscovered deep-splicing variants, genetic and environmental modifiers, digenic/oligogenic inheritance, skewed methylation patterns, and mutational burden. Symptomatic heterozygotes are rarely reported in the literature, mainly because most did not undergo the complete diagnostic procedure, so alternative diagnoses could not be conclusively excluded. However, despite the increasing accessibility to high-throughput technologies, there still seems to be a small group of patients with mild symptoms and just one variant of autosomes in biallelic diseases. Here, we present some examples, the current state of knowledge, and possible explanations for this phenomenon, and thus argue against the existing dominant/recessive classification.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

Thyroid cancer harboring PTEN and TP53 mutations: A peculiar molecular and clinical case report

To date, the molecular mechanisms that underline aggressiveness and resistance to tyrosine kinase inhibitors in some thyroid carcinomas (TCs) are not known yet. We report the case of a young patient with a metastatic poorly differentiated (PDTC) and follicular thyroid carcinoma (FTC) refractory to conventional therapies and to Sorafenib. The patient, despite an initial partial response, died of progressive disease 21 months after diagnosis. The genetic analysis performed on the primary tumor and on lymph nodes and distant metastases allowed to identify a frameshift mutation (p.P248Tfs*5) in the PTEN gene, never described in TC. This mutation was present in the primary tumor and, with a lower allelic frequency, in metastases diagnosed after treatment with Sorafenib. Mutations in TP53 (p.C135Y and c.920-2A>G previously detected in anaplastic carcinomas and p.M133R never found in TC) were also detected in the primary tissue together with a mono-allelic expression of the p.C135Y mutant at RNA level. At metastatic sites level, we found only the TP53 splicing mutation c.920-2A>G. The presence of defects in mismatch repair (MMR) proteins and genomic instability was also evaluated. The primary tumor showed a partial expression of MMR proteins together with a strong genomic instability. In conclusion, we demonstrated that the rare combination of somatic PTEN and TP53 mutations in a patient with a metastatic FTC, together with the presence of tumor heterogeneity and genomic instability, might be associated with a high tumor aggressiveness and resistance to treatments.

Human GPR17 missense variants identified in metabolic disease patients have distinct downstream signaling profiles

GPR17 is a G-protein-coupled receptor (GPCR) implicated in the regulation of glucose metabolism and energy homeostasis. Such evidence is primarily drawn from mouse knockout studies and suggests GPR17 as a potential novel therapeutic target for the treatment of metabolic diseases. However, links between human GPR17 genetic variants, downstream cellular signaling, and metabolic diseases have yet to be reported. Here, we analyzed GPR17 coding sequences from control and disease cohorts consisting of individuals with adverse clinical metabolic deficits including severe insulin resistance, hypercholesterolemia, and obesity. We identified 18 nonsynonymous GPR17 variants, including eight variants that were exclusive to the disease cohort. We characterized the protein expression levels, membrane localization, and downstream signaling profiles of nine GPR17 variants (F43L, V96M, V103M, D105N, A131T, G136S, R248Q, R301H, and G354V). These nine GPR17 variants had similar protein expression and subcellular localization as wild-type GPR17; however, they showed diverse downstream signaling profiles. GPR17-G136S lost the capacity for agonist-mediated cAMP, Ca²⁺, and β-arrestin signaling. GPR17-V96M retained cAMP inhibition similar to GPR17-WT, but showed impaired Ca²⁺ and β-arrestin signaling. GPR17-D105N displayed impaired cAMP and Ca²⁺ signaling, but unaffected agonist-stimulated β-arrestin recruitment. The identification and functional profiling of naturally occurring human GPR17 variants from individuals with metabolic diseases revealed receptor variants with diverse signaling profiles, including differential signaling perturbations that resulted in GPCR signaling bias. Our findings provide a framework for structure–function relationship studies of GPR17 signaling and metabolic disease.

Trappc9 deficiency causes parent-of-origin dependent microcephaly and obesity

Some imprinted genes exhibit parental origin specific expression bias rather than being transcribed exclusively from one copy. The physiological relevance of this remains poorly understood. In an analysis of brain-specific allele-biased expression, we identified that Trappc9, a cellular trafficking factor, was expressed predominantly (~70%) from the maternally inherited allele. Loss-of-function mutations in human TRAPPC9 cause a rare neurodevelopmental syndrome characterized by microcephaly and obesity. By studying Trappc9 null mice we discovered that homozygous mutant mice showed a reduction in brain size, exploratory activity and social memory, as well as a marked increase in body weight. A role for Trappc9 in energy balance was further supported by increased ad libitum food intake in a child with TRAPPC9 deficiency. Strikingly, heterozygous mice lacking the maternal allele (70% reduced expression) had pathology similar to homozygous mutants, whereas mice lacking the paternal allele (30% reduction) were phenotypically normal. Taken together, we conclude that Trappc9 deficient mice recapitulate key pathological features of TRAPPC9 mutations in humans and identify a role for Trappc9 and its imprinting in controlling brain development and metabolism.

Monoallelic, antisense and total RNA transcription in an in vitro neural differentiation system based on F1 hybrid mice

We developed an in vitro system to differentiate embryonic stem cells (ESCs) derived from reciprocally crossed F1 hybrid mice into neurons, and used it to investigate poly(A)+ and total RNA transcription at different stages of cell differentiation. By comparing expression profiles of transcripts assembled from 20 RNA sequencing datasets [2 alleles×(2 cell lines×4 time-points+2 mouse brains)], the relative influence of strain, cell and parent specificities to overall expression could be assessed. Divergent expression profiles of ESCs converged tightly at neural progenitor stage. Patterns of temporal variation of monoallelically expressed transcripts and antisense transcripts were quantified. Comparison of sense and antisense transcript pairs within the poly(A)+ sample, within the total RNA sample, and across poly(A)+ and total RNA samples revealed distinct rates of pairs showing anti-correlated expression variation. Unique patterns of sharing of poly(A)+ and poly(A)- transcription were identified in distinct RNA species. Regulation and functionality of monoallelic expression, antisense transcripts and poly(A)- transcription remain elusive. We demonstrated the effectiveness of our approach to capture these transcriptional activities, and provided new resources to elucidate the mammalian developmental transcriptome.

Molecular defects in thyroid dysgenesis

Congenital hypothyroidism (CH) is a neonatal endocrine disorder that might occur as itself or be associated to congenital extra-thyroidal defects. About 85% of affected subjects experience thyroid dysgenesis (TD), characterized by defect in thyroid gland development. In vivo experiments on null mice paved the way for the identification of genes involved thyroid morphogenesis and development, whose mutation has been strongly associated to TD. Most of them are thyroid-specific transcription factors expressed during early thyroid development. Despite the arduous effort in unraveling the genetics of TD in animal models, up to now these data have been discontinuously confirmed in humans and only 5% of TD have associated with known null mice-related mutations (mainly PAX8 and TSHR). Notwithstanding, the advance in genetic testing represented by the next-generation sequencing (NGS) approach is steadily increasing the list of genes whose highly penetrant mutation predisposes to TD. In this review we intend to outline the molecular bases of TD, summarizing the current knowledge on thyroid development in both mice and humans and delineating the genetic features of its monogenetic forms. We will also highlight current strategies to enhance the insight into the non-Mendelian mechanisms of abnormal thyroid development.

MaGIC: a machine learning tool set and web application for monoallelic gene inference from chromatin

Background

A large fraction of human and mouse autosomal genes are subject to random monoallelic expression (MAE), an epigenetic mechanism characterized by allele-specific gene expression that varies between clonal cell lineages. MAE is highly cell-type specific and mapping it in a large number of cell and tissue types can provide insight into its biological function. Its detection, however, remains challenging.

Results

We previously reported that a sequence-independent chromatin signature identifies, with high sensitivity and specificity, genes subject to MAE in multiple tissue types using readily available ChIP-seq data. Here we present an implementation of this method as a user-friendly, open-source software pipeline for monoallelic gene inference from chromatin (MaGIC). The source code for the MaGIC pipeline and the Shiny app is available at https://github.com/gimelbrantlab/magic.

Conclusion

The pipeline can be used by researchers to map monoallelic expression in a variety of cell types using existing models and to train new models with additional sets of chromatin marks.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2679-7) contains supplementary material, which is available to authorized users.

Risk alleles of genes with monoallelic expression are enriched in gain-of-function variants and depleted in loss-of-function variants for neurodevelopmental disorders

Over 3000 human genes can be expressed from a single allele in one cell, and from the other allele-or both-in neighboring cells. Little is known about the consequences of this epigenetic phenomenon, monoallelic expression (MAE). We hypothesized that MAE increases expression variability, with a potential impact on human disease. Here, we use a chromatin signature to infer MAE for genes in lymphoblastoid cell lines and human fetal brain tissue. We confirm that across clones MAE status correlates with expression level, and that in human tissue data sets, MAE genes show increased expression variability. We then compare mono- and biallelic genes at three distinct scales. In the human population, we observe that genes with polymorphisms influencing expression variance are more likely to be MAE (P<1.1 × 10-6). At the trans-species level, we find gene expression differences and directional selection between humans and chimpanzees more common among MAE genes (P<0.05). Extending to human disease, we show that MAE genes are under-represented in neurodevelopmental copy number variants (CNVs) (P<2.2 × 10-10), suggesting that pathogenic variants acting via expression level are less likely to involve MAE genes. Using neuropsychiatric single-nucleotide polymorphism (SNP) and single-nucleotide variant (SNV) data, we see that genes with pathogenic expression-altering or loss-of-function variants are less likely MAE (P<7.5 × 10-11) and genes with only missense or gain-of-function variants are more likely MAE (P<1.4 × 10-6). Together, our results suggest that MAE genes tolerate a greater range of expression level than biallelic expression (BAE) genes, and this information may be useful in prediction of pathogenicity.

Lysine-specific demethylase-1 (LSD1) depletion disrupts monogenic and monoallelic odorant receptor (OR) expression in an olfactory neuronal cell line

Function of the mammalian olfactory system depends on specialized olfactory sensory neurons (OSNs) that each express only one allele ("monoallelic") of one odorant receptor (OR) gene ("monogenic"). The lysine-specific demethylase-1 (LSD1) protein removes activating H3K4 or silencing H3K9 methylation marks in a variety of developmental contexts, and is thought to be important for proper OR regulation. Most of the focus in the field has been on a potential "activating" function for LSD1; e.g., in the demethylation of H3K9 associated with the expressed OR allele. Here we show that depletion of LSD1 in an immortalized olfactory-placode-derived cell line (OP6) results in multigenic and multiallelic OR transcription per cell, while not seemingly disrupting the ability of these cells to activate new OR genes during clonal expansion. These results are consistent with LSD1 having a role in silencing additional OR alleles, as opposed to being required for the activation of OR alleles, within the OP6 cellular context.

Making choices-how stochastic decisions determine disease progression

Park and Pfeifer discuss Ginart et al.’s use of RNA FISH to measure allele-specific expression in single cells derived from a mouse model of a human imprinting disorder and their finding that differential methylation at the imprinting control region underlies the observed cell-to-cell variability.

In this issue of Genes & Development, Ginart and colleagues (pp. 567–578) study a mouse model for Russell-Silver syndrome (RSS) and show that similar cells within one individual can display distinct gene expression patterns because of epigenetic marks that are established stochastically during early development. Their results provide an excellent explanation for phenotypes seen in RSS and other imprinting disorders and especially help us understand how patients with similar or even identical genetic mutations can display distinct disease profiles.

Visualizing allele-specific expression in single cells reveals epigenetic mosaicism in an H19 loss-of-imprinting mutant

Imprinting defects can lead to inappropriate expression from the normally silenced allele, but it remains unclear whether every cell in a mutant organism follows the population average. Ginart et al. applied a new fluorescence in situ hybridization method that measures allele-specific expression in single cells to address this question in mutants exhibiting aberrant H19/Igf2 imprinting.

Imprinting is a classic mammalian epigenetic phenomenon that results in expression from a single parental allele. Imprinting defects can lead to inappropriate expression from the normally silenced allele, but it remains unclear whether every cell in a mutant organism follows the population average, which would have profound implications for human imprinting disorders. Here, we apply a new fluorescence in situ hybridization method that measures allele-specific expression in single cells to address this question in mutants exhibiting aberrant H19/Igf2 (insulin-like growth factor 2) imprinting. We show that mutant primary embryonic mouse fibroblasts are comprised of two subpopulations: one expressing both H19 alleles and another expressing only the maternal copy. Only in the latter cell population is Igf2 expression detected. Furthermore, the two subpopulations are stable in that cells do not interconvert between the two expression patterns. Combined small input methylation analysis and transcriptional imaging revealed that these two mutant subpopulations exhibit distinct methylation patterns at their imprinting control regions. Consistently, pharmacological inhibition of DNA methylation reduced the proportion of monoallelic cells. Importantly, we observed that the same two subpopulations are also present in vivo within murine cardiac tissue. Our results establish that imprinting disorders can display striking single-cell heterogeneity in their molecular phenotypes and suggest that such heterogeneity may underlie epigenetic mosaicism in human imprinting disorders.

Demonstration of Autosomal Monoallelic Expression in Thyroid Tissue Assessed by Whole-Exome and Bulk RNA Sequencing

BACKGROUND

Congenital hypothyroidism due to thyroid dysgenesis (CHTD) is a disorder with a prevalence of 1/4000 live births, the cause of which remains unknown. The most common diagnostic category is thyroid ectopy, which occurs in up to 80% of CHTD cases. CHTD is predominantly not inherited and has a high discordance rate (>92%) between monozygotic (MZ) twins. The sporadic nature of CHTD might be explained by somatic events such as autosomal monoallelic expression (AME), given that genes expressed in a monoallelic way are more vulnerable to otherwise benign monoallelict genetic or epigenetic mutations.

OBJECTIVE

The aim of this study was to search for complete (90%) AME in normal and dysgenetic thyroid tissues.

METHODS

Aggregated analysis of whole-exome and bulk RNA sequencing was performed on two ectopic thyroids, four normal thyroids, and the human thyroid cell line Nthy-ori.

RESULTS

A median of 5062 (range 2081-5270) genes per sample showed sufficient numbers of heterozygous single nucleotide polymorphisms to be informative. The median monoallelic expression represented 22 (range 16-32) of the informative genes for each thyroid sample. Examples of genes displaying AME are FCGBP, ZNF331, USP10, BCLAF1, and some HLA genes; these genes are involved in epithelial-mesenchymal transition, cell migration, cancer, and immunity.

CONCLUSIONS

AME may account for the high discordance rate observed between MZ twins and for the sporadic nature of CHTD. These findings also have implications for other pathologies, including cancers and autoimmune disorders of the thyroid.

Genes with monoallelic expression contribute disproportionately to genetic diversity in humans

An unexpectedly large number of human autosomal genes are subject to monoallelic expression (MAE). Our analysis of 4,227 such genes uncovers surprisingly high genetic variation across human populations. This increased diversity is unlikely to reflect relaxed purifying selection. Remarkably, MAE genes exhibit an elevated recombination rate and an increased density of hypermutable sequence contexts. However, these factors do not fully account for the increased diversity. We find that the elevated nucleotide diversity of MAE genes is also associated with greater allelic age: variants in these genes tend to be older and are enriched in polymorphisms shared by Neanderthals and chimpanzees. Both synonymous and nonsynonymous alleles of MAE genes have elevated average population frequencies. We also observed strong enrichment of the MAE signature among genes reported to evolve under balancing selection. We propose that an important biological function of widespread MAE might be the generation of cell-to-cell heterogeneity; the increased genetic variation contributes to this heterogeneity.

The 2016 database issue of Nucleic Acids Research and an updated molecular biology database collection

The 2016 Database Issue of Nucleic Acids Research starts with overviews of the resources provided by three major bioinformatics centers, the U.S. National Center for Biotechnology Information (NCBI), the European Bioinformatics Institute (EMBL-EBI) and Swiss Institute for Bioinformatics (SIB). Also included are descriptions of 62 new databases and updates on 95 databases that have been previously featured in NAR plus 17 previously described elsewhere. A number of papers in this issue deal with resources on nucleic acids, including various kinds of non-coding RNAs and their interactions, molecular dynamics simulations of nucleic acid structure, and two databases of super-enhancers. The protein database section features important updates on the EBI's Pfam, PDBe and PRIDE databases, as well as a variety of resources on pathways, metabolomics and metabolic modeling. This issue also includes updates on popular metagenomics resources, such as MG-RAST, EBI Metagenomics, and probeBASE, as well as a newly compiled Human Pan-Microbe Communities database. A significant fraction of the new and updated databases are dedicated to the genetic basis of disease, primarily cancer, and various aspects of drug research, including resources for patented drugs, their side effects, withdrawn drugs, and potential drug targets. A further six papers present updated databases of various antimicrobial and anticancer peptides. The entire Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/). The NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/c/, has been updated with the addition of 88 new resources and removal of 23 obsolete websites, which brought the current listing to 1685 databases.

dbMAE: the database of autosomal monoallelic expression

Recently, data on ‘random’ autosomal monoallelic expression has become available for the entire genome in multiple human and mouse tissues and cell types, creating a need for better access and dissemination. The database of autosomal monoallelic expression (dbMAE; https://mae.hms.harvard.edu) incorporates data from multiple recent reports of genome-wide analyses. These include transcriptome-wide analyses of allelic imbalance in clonal cell populations based on sequence polymorphisms, as well as indirect identification, based on a specific chromatin signature present in MAE gene bodies. Currently, dbMAE contains transcriptome-wide chromatin identification calls for 8 human and 21 mouse tissues, and describes over 16 000 murine and ∼700 human cases of directly measured biased expression, compiled from allele-specific RNA-seq and genotyping array data. All data are manually curated. To ensure cross-publication uniformity, we performed re-analysis of transcriptome-wide RNA-seq data using the same pipeline. Data are accessed through an interface that allows for basic and advanced searches; all source references, including raw data, are clearly described and hyperlinked. This ensures the utility of the resource as an initial screening tool for those interested in investigating the role of monoallelic expression in their specific genes and tissues of interest.

Chromatin Signature Identifies Monoallelic Gene Expression Across Mammalian Cell Types

Monoallelic expression of autosomal genes (MAE) is a widespread epigenetic phenomenon which is poorly understood, due in part to current limitations of genome-wide approaches for assessing it. Recently, we reported that a specific histone modification signature is strongly associated with MAE and demonstrated that it can serve as a proxy of MAE in human lymphoblastoid cells. Here, we use murine cells to establish that this chromatin signature is conserved between mouse and human and is associated with MAE in multiple cell types. Our analyses reveal extensive conservation in the identity of MAE genes between the two species. By analyzing MAE chromatin signature in a large number of cell and tissue types, we show that it remains consistent during terminal cell differentiation and is predominant among cell-type specific genes, suggesting a link between MAE and specification of cell identity.

Allelic Imbalance Is a Prevalent and Tissue-Specific Feature of the Mouse Transcriptome

In mammals, several classes of monoallelic genes have been identified, including those subject to X-chromosome inactivation (XCI), genomic imprinting, and random monoallelic expression (RMAE). However, the extent to which these epigenetic phenomena are influenced by underlying genetic variation is unknown. Here we perform a systematic classification of allelic imbalance in mouse hybrids derived from reciprocal crosses of divergent strains. We observe that deviation from balanced biallelic expression is common, occurring in ∼20% of the mouse transcriptome in a given tissue. Allelic imbalance attributed to genotypic variation is by far the most prevalent class and typically is tissue-specific. However, some genotype-based imbalance is maintained across tissues and is associated with greater genetic variation, especially in 5' and 3' termini of transcripts. We further identify novel random monoallelic and imprinted genes and find that genotype can modify penetrance of parental origin even in the setting of large imprinted regions. Examination of nascent transcripts in single cells from inbred parental strains reveals that genes showing genotype-based imbalance in hybrids can also exhibit monoallelic expression in isogenic backgrounds. This surprising observation may suggest a competition between alleles and/or reflect the combined impact of cis- and trans-acting variation on expression of a given gene. Our findings provide novel insights into gene regulation and may be relevant to human genetic variation and disease.

RNA sequencing: current and prospective uses in metabolic research

Previous global RNA analysis was restricted to known transcripts in species with a defined transcriptome. Next generation sequencing has transformed transcriptomics by making it possible to analyse expressed genes with an exon level resolution from any tissue in any species without any a priori knowledge of which genes that are being expressed, splice patterns or their nucleotide sequence. In addition, RNA sequencing is a more sensitive technique compared with microarrays with a larger dynamic range, and it also allows for investigation of imprinting and allele-specific expression. This can be done for a cost that is able to compete with that of a microarray, making RNA sequencing a technique available to most researchers. Therefore RNA sequencing has recently become the state of the art with regards to large-scale RNA investigations and has to a large extent replaced microarrays. The only drawback is the large data amounts produced, which together with the complexity of the data can make a researcher spend far more time on analysis than performing the actual experiment.

Random monoallelic expression: regulating gene expression one allele at a time

Monoallelic gene expression is a remarkable process in which transcription occurs from only one of two homologous alleles in a diploid cell. Interestingly, between 0.5% and 15% of autosomal genes exhibit random monoallelic gene expression, in which different cells express only one allele independently of the underlying genomic sequence, in a cell type-specific manner. Recently, genome-wide studies have increased our understanding of the cell type-specific incidence of random monoallelic gene expression, and how the imbalance in allelic expression is distinguished within the cell and potentially maintained across cell generations. Monoallelic gene expression is likely generated through stochastic independent regulation of the two alleles upon differentiation, and has varied implications for the cell and organism, in particular with respect to disease.