Developmental dynamics and disease potential of random monoallelic gene expression

Detection of allele-specific expression in spatial transcriptomics with spASE

Spatial transcriptomics technologies permit the study of the spatial distribution of RNA at near-single-cell resolution genome-wide. However, the feasibility of studying spatial allele-specific expression (ASE) from these data remains uncharacterized. Here, we introduce spASE, a computational framework for detecting and estimating spatial ASE. To tackle the challenges presented by cell type mixtures and a low signal to noise ratio, we implement a hierarchical model involving additive mixtures of spatial smoothing splines. We apply our method to allele-resolved Visium and Slide-seq from the mouse cerebellum and hippocampus and report new insight into the landscape of spatial and cell type-specific ASE therein.

Gene-specific somatic epigenetic mosaicism of FDFT1 underlies a non-hereditary localized form of porokeratosis

Porokeratosis is a clonal keratinization disorder characterized by solitary, linearly arranged, or generally distributed multiple skin lesions. Previous studies showed that genetic alterations in MVK, PMVK, MVD, or FDPS-genes in the mevalonate pathway-cause hereditary porokeratosis, with skin lesions harboring germline and lesion-specific somatic variants on opposite alleles. Here, we identified non-hereditary porokeratosis associated with epigenetic silencing of FDFT1, another gene in the mevalonate pathway. Skin lesions of the generalized form had germline and lesion-specific somatic variants on opposite alleles in FDFT1, representing FDFT1-associated hereditary porokeratosis identified in this study. Conversely, lesions of the solitary or linearly arranged localized form had somatic bi-allelic promoter hypermethylation or mono-allelic promoter hypermethylation with somatic genetic alterations on opposite alleles in FDFT1, indicating non-hereditary porokeratosis. FDFT1 localization was uniformly diminished within the lesions, and lesion-derived keratinocytes showed cholesterol dependence for cell growth and altered expression of genes related to cell-cycle and epidermal development, confirming that lesions form by clonal expansion of FDFT1-deficient keratinocytes. In some individuals with the localized form, gene-specific promoter hypermethylation of FDFT1 was detected in morphologically normal epidermis adjacent to methylation-related lesions but not distal to these lesions, suggesting that asymptomatic somatic epigenetic mosaicism of FDFT1 predisposes certain skin areas to the disease. Finally, consistent with its genetic etiology, topical statin treatment ameliorated lesions in FDFT1-deficient porokeratosis. In conclusion, we identified bi-allelic genetic and/or epigenetic alterations of FDFT1 as a cause of porokeratosis and shed light on the pathogenesis of skin mosaicism involving clonal expansion of epigenetically altered cells.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA

In placental females, one copy of the two X chromosomes is largely silenced during a narrow developmental time window, in a process mediated by the non-coding RNA Xist1. Here, we demonstrate that Xist can initiate X-chromosome inactivation (XCI) well beyond early embryogenesis. By modifying its endogenous level, we show that Xist has the capacity to actively silence genes that escape XCI both in neuronal progenitor cells (NPCs) and in vivo, in mouse embryos. We also show that Xist plays a direct role in eliminating TAD-like structures associated with clusters of escapee genes on the inactive X chromosome, and that this is dependent on Xist's XCI initiation partner, SPEN2. We further demonstrate that Xist's function in suppressing gene expression of escapees and topological domain formation is reversible for up to seven days post-induction, but that sustained Xist up-regulation leads to progressively irreversible silencing and CpG island DNA methylation of facultative escapees. Thus, the distinctive transcriptional and regulatory topologies of the silenced X chromosome is actively, directly - and reversibly - controlled by Xist RNA throughout life.

Maternal histone methyltransferases antagonistically regulate monoallelic expression in C. elegans

Undefined epigenetic programs act to probabilistically silence individual autosomal alleles, generating unique individuals, even from genetic clones. This sort of random monoallelic expression can explain variation in traits and diseases that differences in genes and environments cannot. Here, we developed the nematode Caenorhabditis elegans to study monoallelic expression in whole tissues, and defined a developmental genetic regulation pathway. We found maternal H3K9 histone methyltransferase (HMT) SET-25/SUV39/G9a works with HPL-2/HP1 and LIN-61/L3MBTL2 to randomly silence alleles in the intestinal progenitor E-cell of 8-cell embryos to cause monoallelic expression. SET-25 was antagonized by another maternal H3K9 HMT, MET-2/SETDB1, which works with LIN-65/ATF7ZIP and ARLE-14/ARL14EP to prevent monoallelic expression. The HMT-catalytic SET domains of both MET-2 and SET-25 were required for regulating monoallelic expression. Our data support a model wherein SET-25 and MET-2 regulate histones during development to generate patterns of somatic monoallelic expression that are persistent but not heritable.

Mechanisms of pathogenicity and the quest for genetic modifiers of kidney disease in branchiootorenal syndrome

Backgound

Branchiootorenal (BOR) syndrome is an autosomal dominant disorder caused by pathogenic EYA1 variants and clinically characterized by auricular malformations with hearing loss, branchial arch anomalies, and congenital anomalies of the kidney and urinary tract. BOR phenotypes are highly variable and heterogenous. While random monoallelic expression is assumed to explain this phenotypic heterogeneity, the potential role of modifier genes has not yet been explored.

Methods

Through thorough phenotyping and exome sequencing, we studied one family with disease presentation in at least four generations in both clinical and genetic terms. Functional investigation of the single associated EYA1 variant c.1698+1G>A included splice site analysis and assessment of EYA1 distribution in patient-derived fibroblasts. The candidate modifier gene CYP51A1 was evaluated by histopathological analysis of murine Cyp51+/- and Cyp51-/- kidneys. As the gene encodes the enzyme lanosterol 14α-demethylase, we assessed sterol intermediates in patient blood samples as well.

Results

The EYA1 variant c.1698+1G>A resulted in functional deletion of the EYA domain by exon skipping. The EYA domain mediates protein-protein interactions between EYA1 and co-regulators of transcription. EYA1 abundance was reduced in the nuclear compartment of patient-derived fibroblasts, suggesting impaired nuclear translocation of these protein complexes. Within the affected family, renal phenotypes spanned from normal kidney function in adulthood to chronic kidney failure in infancy. By analyzing exome sequencing data for variants that potentially play roles as genetic modifiers, we identified a canonical splice site alteration in CYP51A1 as the strongest candidate variant.

Conclusion

In this study, we demonstrate pathogenicity of EYA1 c.1698+1G>A, propose a mechanism for dysfunction of mutant EYA1, and conjecture CYP51A1 as a potential genetic modifier of renal involvement in BOR syndrome.

Genomic Imprinting and Random Monoallelic Expression

The review discusses the mechanisms of monoallelic expression, such as genomic imprinting, in which gene transcription depends on the parental origin of the allele, and random monoallelic transcription. Data on the regulation of gene activity in the imprinted regions are summarized with a particular focus on the areas controlling imprinting and factors influencing the variability of the imprintome. The prospects of studies of the monoallelic expression are discussed.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Foreign RNA spike-ins enable accurate allele-specific expression analysis at scale

MOTIVATION

Analysis of allele-specific expression is strongly affected by the technical noise present in RNA-seq experiments. Previously, we showed that technical replicates can be used for precise estimates of this noise, and we provided a tool for correction of technical noise in allele-specific expression analysis. This approach is very accurate but costly due to the need for two or more replicates of each library. Here, we develop a spike-in approach which is highly accurate at only a small fraction of the cost.

RESULTS

We show that a distinct RNA added as a spike-in before library preparation reflects technical noise of the whole library and can be used in large batches of samples. We experimentally demonstrate the effectiveness of this approach using combinations of RNA from species distinguishable by alignment, namely, mouse, human, and Caenorhabditis elegans. Our new approach, controlFreq, enables highly accurate and computationally efficient analysis of allele-specific expression in (and between) arbitrarily large studies at an overall cost increase of ∼5%.

AVAILABILITY AND IMPLEMENTATION

Analysis pipeline for this approach is available at GitHub as R package controlFreq (github.com/gimelbrantlab/controlFreq).

Rethinking (again) Hardy-Weinberg and genetic drift in undergraduate biology

Designing effective curricula is challenging. Content decisions can impact both learning outcomes and student engagement. As an example consider the place of Hardy-Weinberg equilibria (HWE) and genetic drift calculations in introductory biology courses, as discussed by Masel (2012). Given that population genetics, "a fairly arcane speciality", can be difficult to grasp, there is little justification for introducing introductory students to HWE calculations. It is more useful to introduce them to the behavior of alleles in terms of basic features of biological systems, and that in the absence of selection recessive alleles are no "weaker" or preferentially lost from a population than are dominant alleles. On the other hand, stochastic behaviors, such as genetic drift, are ubiquitous in biological systems and often play functionally significant roles; they can be introduced to introductory students in mechanistic and probabilistic terms. Specifically, genetic drift emerges from the stochastic processes involved in meiotic chromosome segregation and recombination. A focus on stochastic processes may help counteract naive bio-deterministic thinking and can reinforce, for students, the value of thinking quantitatively about biological processes.

Foreign RNA spike-ins enable accurate allele-specific expression analysis at scale

Motivation

Analysis of allele-specific expression is strongly affected by the technical noise present in RNA-seq experiments. Previously, we showed that technical replicates can be used for precise estimates of this noise, and we provided a tool for correction of technical noise in allele-specific expression analysis. This approach is very accurate but costly due to the need for two or more replicates of each library. Here, we develop a spike-in approach that is highly accurate at only a small fraction of the cost.

Results

We show that a distinct RNA added as a spike-in before library preparation reflects technical noise of the whole library and can be used in large batches of samples. We experimentally demonstrate the effectiveness of this approach using combinations of RNA from species distinguishable by alignment, namely, mouse, human, and C.elegans . Our new approach, controlFreq , enables highly accurate and computationally efficient analysis of allele-specific expression in (and between) arbitrarily large studies at an overall cost increase of ~ 5%.

Availability

Analysis pipeline for this approach is available at GitHub as R package controlFreq ( github.com/gimelbrantlab/controlFreq ).

Contact

agimelbrant@altius.org.

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.

Allele-specific differential regulation of monoallelically expressed autosomal genes in the cardiac lineage

Each mammalian autosomal gene is represented by two alleles in diploid cells. To our knowledge, no insights have been made in regard to allele-specific regulatory mechanisms of autosomes. Here we use allele-specific single cell transcriptomic analysis to elucidate the establishment of monoallelic gene expression in the cardiac lineage. We find that monoallelically expressed autosomal genes in mESCs and mouse blastocyst cells are differentially regulated based on the genetic background of the parental alleles. However, the genetic background of the allele does not affect the establishment of monoallelic genes in differentiated cardiomyocytes. Additionally, we observe epigenetic differences between deterministic and random autosomal monoallelic genes. Moreover, we also find a greater contribution of the maternal versus paternal allele to the development and homeostasis of cardiac tissue and in cardiac health, highlighting the importance of maternal influence in male cardiac tissue homeostasis. Our findings emphasize the significance of allele-specific insights into gene regulation in development, homeostasis and disease.

In Vivo Clonal Analysis Reveals Random Monoallelic Expression in Lymphocytes That Traces Back to Hematopoietic Stem Cells

Evaluating the epigenetic landscape in the stem cell compartment at the single-cell level is essential to assess the cells’ heterogeneity and predict their fate. Here, using a genome-wide transcriptomics approach in vivo, we evaluated the allelic expression imbalance in the progeny of single hematopoietic cells (HSCs) as a read-out of epigenetic marking. After 4 months of extensive proliferation and differentiation, we found that X-chromosome inactivation (XCI) is tightly maintained in all single-HSC derived hematopoietic cells. In contrast, the vast majority of the autosomal genes did not show clonal patterns of random monoallelic expression (RME). However, a persistent allele-specific autosomal transcription in HSCs and their progeny was found in a rare number of cases, none of which has been previously reported. These data show that: 1) XCI and RME in the autosomal chromosomes are driven by different mechanisms; 2) the previously reported high frequency of genes under RME in clones expanded in vitro (up to 15%) is not found in clones undergoing multiple differentiation steps in vivo; 3) prior to differentiation, HSCs have stable patterns of autosomal RME. We propose that most RME patterns in autosomal chromosomes are erased and established de novo during cell lineage differentiation.

Binary outcomes of enhancer activity underlie stable random monoallelic expression

Mitotically stable random monoallelic gene expression (RME) is documented for a small percentage of autosomal genes. We developed an in vivo genetic model to study the role of enhancers in RME using high-resolution single-cell analysis of natural killer (NK) cell receptor gene expression and enhancer deletions in the mouse germline. Enhancers of the RME NK receptor genes were accessible and enriched in H3K27ac on silent and active alleles alike in cells sorted according to allelic expression status, suggesting enhancer activation and gene expression status can be decoupled. In genes with multiple enhancers, enhancer deletion reduced gene expression frequency, in one instance converting the universally expressed gene encoding NKG2D into an RME gene, recapitulating all aspects of natural RME including mitotic stability of both the active and silent states. The results support the binary model of enhancer action, and suggest that RME is a consequence of general properties of gene regulation by enhancers rather than an RME-specific epigenetic program. Therefore, many and perhaps all genes may be subject to some degree of RME. Surprisingly, this was borne out by analysis of several genes that define different major hematopoietic lineages, that were previously thought to be universally expressed within those lineages: the genes encoding NKG2D, CD45, CD8α, and Thy-1. We propose that intrinsically probabilistic gene allele regulation is a general property of enhancer-controlled gene expression, with previously documented RME representing an extreme on a broad continuum.

RNA sequencing-based screen for reactivation of silenced alleles of autosomal genes

In mammalian cells, maternal and paternal alleles usually have similar transcriptional activity. Epigenetic mechanisms such as X-chromosome inactivation (XCI) and imprinting were historically viewed as rare exceptions to this rule. Discovery of autosomal monoallelic autosomal expression (MAE) a decade ago revealed an additional allele-specific mode regulating thousands of mammalian genes. Despite MAE prevalence, its mechanistic basis remains unknown. Using an RNA sequencing-based screen for reactivation of silenced alleles, we identified DNA methylation as key mechanism of MAE mitotic maintenance. In contrast with the all-or-nothing allelic choice in XCI, allele-specific expression in MAE loci is tunable, with exact allelic imbalance dependent on the extent of DNA methylation. In a subset of MAE genes, allelic imbalance was insensitive to DNA demethylation, implicating additional mechanisms in MAE maintenance in these loci. Our findings identify a key mechanism of MAE maintenance and provide basis for understanding the biological role of MAE.

A cis-acting mechanism mediates transcriptional memory at Polycomb target genes in mammals

Epigenetic inheritance of gene expression states enables a single genome to maintain distinct cellular identities. How histone modifications contribute to this process remains unclear. Using global chromatin perturbations and local, time-controlled modulation of transcription, we establish the existence of epigenetic memory of transcriptional activation for genes that can be silenced by the Polycomb group. This property emerges during cell differentiation and allows genes to be stably switched after a transient transcriptional stimulus. This transcriptional memory state at Polycomb targets operates in cis; however, rather than relying solely on read-and-write propagation of histone modifications, the memory is also linked to the strength of activating inputs opposing Polycomb proteins, and therefore varies with the cellular context. Our data and computational simulations suggest a model whereby transcriptional memory arises from double-negative feedback between Polycomb-mediated silencing and active transcription. Transcriptional memory at Polycomb targets thus depends not only on histone modifications but also on the gene-regulatory network and underlying identity of a cell.

Heterogeneity is not always a source of noise: Stochastic gene expression in regulatory heterozygote

Zygosity of diploid genome (i.e., degree to which two parental alleles of a gene have varied genetic sequences) adds another dimension to stochastic gene expression. The allelic imbalance in chromatin accessibility or divergence in regulatory sequences leads to fitness effects but the quantitative aspects thereof are largely left unexplored. We investigate diploid gene expression systems with homozygous (the same) and heterozygous (varied) combination of alleles in cis-regulatory sequences, not in structural gene loci, and characterize the zygosity-associated stochastic fluctuations in protein abundance. An emerging feat of heterozygosity is its counterintuitive capacity for genetic noise control. Especially when the noise is dominantly contributed to by the fluctuations in duty cycle ("reliability") rather than in process speed ("productivity") of gene expression machinery, its interallelic discrepancy acts to reduce the gene expression noise. These findings offer a novel insight into the rich repertoire of balancing selection enriched in the regulatory elements of immune response genes.

Introns control stochastic allele expression bias

Monoallelic expression (MAE) or extreme allele bias can account for incomplete penetrance, missing heritability and non-Mendelian diseases. In cancer, MAE is associated with shorter patient survival times and higher tumor grade. Prior studies showed that stochastic MAE is caused by stochastic epigenetic silencing, in a gene and tissue-specific manner. Here, we used C. elegans to study stochastic MAE in vivo. We found allele bias/MAE to be widespread within C. elegans tissues, presenting as a continuum from fully biallelic to MAE. We discovered that the presence of introns within alleles robustly decreases MAE. We determined that introns control MAE at distinct loci, in distinct cell types, with distinct promoters, and within distinct coding sequences, using a 5'-intron position-dependent mechanism. Bioinformatic analysis showed human intronless genes are significantly enriched for MAE. Our experimental evidence demonstrates a role for introns in regulating MAE, possibly explaining why some mutations within introns result in disease.

A Comprehensive Characterization of Monoallelic Expression During Hematopoiesis and Leukemogenesis via Single-Cell RNA-Sequencing

Single-cell RNA-sequencing (scRNA-seq) is becoming a powerful tool to investigate monoallelic expression (MAE) in various developmental and pathological processes. However, our knowledge of MAE during hematopoiesis and leukemogenesis is limited. In this study, we conducted a systematic interrogation of MAEs in bone marrow mononuclear cells (BMMCs) at single-cell resolution to construct a MAE atlas of BMMCs. We identified 1,020 constitutive MAEs in BMMCs, which included imprinted genes such as MEG8, NAP1L5, and IRAIN. We classified the BMMCs into six cell types and identified 74 cell type specific MAEs including MTSS1, MOB1A, and TCF12. We further identified 114 random MAEs (rMAEs) at single-cell level, with 78.1% single-allele rMAE and 21.9% biallelic mosaic rMAE. Many MAEs identified in BMMCs have not been reported and are potentially hematopoietic specific, supporting MAEs are functional relevance. Comparison of BMMC samples from a leukemia patient with multiple clinical stages showed the fractions of constitutive MAE were correlated with fractions of leukemia cells in BMMCs. Further separation of the BMMCs into leukemia cells and normal cells showed that leukemia cells have much higher constitutive MAE and rMAEs than normal cells. We identified the leukemia cell-specific MAEs and relapsed leukemia cell-specific MAEs, which were enriched in immune-related functions. These results indicate MAE is prevalent and is an important gene regulation mechanism during hematopoiesis and leukemogenesis. As the first systematical interrogation of constitutive MAEs, cell type specific MAEs, and rMAEs during hematopoiesis and leukemogenesis, the study significantly increased our knowledge about the features and functions of MAEs.

Differential Allelic Expression among Long Non-Coding RNAs

Long non-coding RNAs (lncRNA) comprise a diverse group of non-protein-coding RNAs >200 bp in length that are involved in various normal cellular processes and disease states, and can affect coding gene expression through mechanisms in cis or in trans. Since the discovery of the first functional lncRNAs transcribed by RNA Polymerase II, H19 and Xist, many others have been identified and noted for their unusual transcriptional pattern, whereby expression from one chromosome homolog is strongly favored over the other, also known as mono-allelic or differential allelic expression. lncRNAs with differential allelic expression have been observed to play critical roles in developmental gene regulation, chromosome structure, and disease. Here, we will focus on known examples of differential allelic expression of lncRNAs and highlight recent research describing functional lncRNAs expressed from both imprinted and random mono-allelic expression domains.

X-Chromosome Inactivation and Autosomal Random Monoallelic Expression as "Faux Amis"

X-chromosome inactivation (XCI) and random monoallelic expression of autosomal genes (RMAE) are two paradigms of gene expression regulation where, at the single cell level, genes can be expressed from either the maternal or paternal alleles. X-chromosome inactivation takes place in female marsupial and placental mammals, while RMAE has been described in mammals and also other species. Although the outcome of both processes results in random monoallelic expression and mosaicism at the cellular level, there are many important differences. We provide here a brief sketch of the history behind the discovery of XCI and RMAE. Moreover, we review some of the distinctive features of these two phenomena, with respect to when in development they are established, their roles in dosage compensation and cellular phenotypic diversity, and the molecular mechanisms underlying their initiation and stability.

Semicoordinated allelic-bursting shape dynamic random monoallelic expression in pregastrulation embryos

Recently, allele-specific single-cell RNA-seq analysis has demonstrated widespread dynamic random monoallelic expression of autosomal genes (aRME) in different cell types. However, the prevalence of dynamic aRME during pregastrulation remains unknown. Here, we show that dynamic aRME is widespread in different lineages of pregastrulation embryos. Additionally, the origin of dynamic aRME remains elusive. It is believed that independent transcriptional bursting from each allele leads to dynamic aRME. Here, we show that allelic burst is not perfectly independent; instead it happens in a semicoordinated fashion. Importantly, we show that semicoordinated allelic bursting of genes, particularly with low burst frequency, leads to frequent asynchronous allelic bursting, thereby contributing to dynamic aRME. Furthermore, we found that coordination of allelic bursting is lineage specific and genes regulating the development have a higher degree of coordination. Altogether, our study provides significant insights into the prevalence and origin of dynamic aRME and their developmental relevance during early development.

•

Dynamic aRME is widespread in different lineages of pregastrulation embryos

•

Semicoordinated bursting of genes with low burst frequency leads to dynamic aRME

•

Degree of coordination of allelic bursting is lineage specific

•

Developmental genes have higher degree of coordination of allelic bursting

Systematic Analysis of Monoallelic Gene Expression and Chromatin Accessibility Across Multiple Tissues in Hybrid Mice

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.

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.

Starvation and Climate Change—How to Constrain Cancer Cell Epigenetic Diversity and Adaptability to Enhance Treatment Efficacy

Advanced metastatic cancer is currently not curable and the major barrier to eliminating the disease in patients is the resistance of subpopulations of tumor cells to drug treatments. These resistant subpopulations can arise stochastically among the billions of tumor cells in a patient or emerge over time during therapy due to adaptive mechanisms and the selective pressures of drug therapies. Epigenetic mechanisms play important roles in tumor cell diversity and adaptability, and are regulated by metabolic pathways. Here, I discuss knowledge from ecology, evolution, infectious disease, species extinction, metabolism and epigenetics to synthesize a roadmap to a clinically feasible approach to help homogenize tumor cells and, in combination with drug treatments, drive their extinction. Specifically, cycles of starvation and hyperthermia could help synchronize tumor cells and constrain epigenetic diversity and adaptability by limiting substrates and impairing the activity of chromatin modifying enzymes. Hyperthermia could also help prevent cancer cells from entering dangerous hibernation-like states. I propose steps to a treatment paradigm to help drive cancer extinction that builds on the successes of fasting, hyperthermia and immunotherapy and is achievable in patients. Finally, I highlight the many unknowns, opportunities for discovery and that stochastic gene and allele level epigenetic mechanisms pose a major barrier to cancer extinction that warrants deeper investigation.

Olfactory expression of trace amine-associated receptors requires cooperative cis-acting enhancers

Olfactory sensory neurons express a large family of odorant receptors (ORs) and a small family of trace amine-associated receptors (TAARs). While both families are subject to so-called singular expression (expression of one allele of one gene), the mechanisms underlying TAAR gene choice remain obscure. Here, we report the identification of two conserved sequence elements in the mouse TAAR cluster (T-elements) that are required for TAAR gene expression. We observed that cell-type-specific expression of a TAAR-derived transgene required either T-element. Moreover, deleting either element reduced or abolished expression of a subset of TAAR genes, while deleting both elements abolished olfactory expression of all TAARs in cis with the mutation. The T-elements exhibit several features of known OR enhancers but also contain highly conserved, unique sequence motifs. Our data demonstrate that TAAR gene expression requires two cooperative cis-acting enhancers and suggest that ORs and TAARs share similar mechanisms of singular expression.

Replicate sequencing libraries are important for quantification of allelic imbalance

A sensitive approach to quantitative analysis of transcriptional regulation in diploid organisms is analysis of allelic imbalance (AI) in RNA sequencing (RNA-seq) data. A near-universal practice in such studies is to prepare and sequence only one library per RNA sample. We present theoretical and experimental evidence that data from a single RNA-seq library is insufficient for reliable quantification of the contribution of technical noise to the observed AI signal; consequently, reliance on one-replicate experimental design can lead to unaccounted-for variation in error rates in allele-specific analysis. We develop a computational approach, Qllelic, that accurately accounts for technical noise by making use of replicate RNA-seq libraries. Testing on new and existing datasets shows that application of Qllelic greatly decreases false positive rate in allele-specific analysis while conserving appropriate signal, and thus greatly improves reproducibility of AI estimates. We explore sources of technical overdispersion in observed AI signal and conclude by discussing design of RNA-seq studies addressing two biologically important questions: quantification of transcriptome-wide AI in one sample, and differential analysis of allele-specific expression between samples.

Parental bias in expression and interaction of genes in the equine placenta

Most autosomal genes in the placenta show a biallelic expression pattern. However, some genes exhibit allele-specific transcription depending on the parental origin of the chromosomes on which the copy of the gene resides. Parentally expressed genes are involved in the reciprocal interaction between maternal and paternal genes, coordinating the allocation of resources between fetus and mother. One of the main challenges of studying parental-specific allelic expression (allele-specific expression [ASE]) in the placenta is the maternal cellular remnant at the fetomaternal interface. Horses (Equus caballus) have an epitheliochorial placenta in which both the endometrial epithelium and the epithelium of the chorionic villi are juxtaposed with minimal extension into the uterine mucosa, yet there is no information available on the allelic gene expression of equine chorioallantois (CA). In the current study, we present a dataset of 1,336 genes showing ASE in the equine CA (https://pouya-dini.github.io/equine-gene-db/) along with a workflow for analyzing ASE genes. We further identified 254 potentially imprinted genes among the parentally expressed genes in the equine CA and evaluated the expression pattern of these genes throughout gestation. Our gene ontology analysis implies that maternally expressed genes tend to decrease the length of gestation, while paternally expressed genes extend the length of gestation. This study provides fundamental information regarding parental gene expression during equine pregnancy, a species with a negligible amount of maternal cellular remnant in its placenta. This information will provide the basis for a better understanding of the role of parental gene expression in the placenta during gestation.

Transcriptional bursts explain autosomal random monoallelic expression and affect allelic imbalance

Transcriptional bursts render substantial biological noise in cellular transcriptomes. Here, we investigated the theoretical extent of allelic expression resulting from transcriptional bursting and how it compared to the amount biallelic, monoallelic and allele-biased expression observed in single-cell RNA-sequencing (scRNA-seq) data. We found that transcriptional bursting can explain the allelic expression patterns observed in single cells, including the frequent observations of autosomal monoallelic gene expression. Importantly, we identified that the burst frequency largely determined the fraction of cells with monoallelic expression, whereas the burst size had little effect on monoallelic observations. The high consistency between the bursting model predictions and scRNA-seq observations made it possible to assess the heterogeneity of a group of cells as their deviation in allelic observations from the expected. Finally, both burst frequency and size contributed to allelic imbalance observations and reinforced that studies of allelic imbalance can be confounded from the inherent noise in transcriptional bursting. Altogether, we demonstrate that allele-level transcriptional bursting renders widespread, although predictable, amounts of monoallelic and biallelic expression in single cells and cell populations.

Genes are transcribed into RNA and further translated into proteins. The maternal and paternal copy of each gene are typically transcribed independently, and transcription itself occur in discrete stochastic bursts (transcriptional bursts). Pioneering single-cell analysis of RNA across cells revealed abundant fluctuations in the amounts of maternal and paternal RNA in cells, with frequent observations of RNA from only the maternal or paternal gene copy (monoallelic expression). In this study, we investigated to which extent the observed monoallelic expression across single cells can be explained by transcriptional bursting. We demonstrate that the process of transcriptional bursting is sufficient to explain the amount of monoallelic expression, and we further demonstrate that the frequency of bursts mainly determines the frequency of monoallelic observations. Furthermore, we show that transcriptional bursts may lead to false positive observations of monoallelic expression across cell populations. Therefore, stochastic transcription renders large fluctuations in allelic origin of RNA in cells over time, including frequent monoallelic observations when profiling single cells.

Genome-wide analysis of allele-specific expression of genes in the model diatom Phaeodactylum tricornutum

Recent advances in next generation sequencing technologies have allowed the discovery of widespread autosomal allele-specific expression (aASE) in mammals and plants with potential phenotypic effects. Extensive numbers of genes with allele-specific expression have been described in the diatom Fragilariopsis cylindrus in association with adaptation to external cues, as well as in Fistulifera solaris in the context of natural hybridization. However, the role of aASE and its extent in diatoms remain elusive. In this study, we investigate allele-specific expression in the model diatom Phaeodactylum tricornutum by the re-analysis of previously published whole genome RNA sequencing data and polymorphism calling. We found that 22% of P. tricornutum genes show moderate bias in allelic expression while 1% show nearly complete monoallelic expression. Biallelic expression associates with genes encoding components of protein metabolism while moderately biased genes associate with functions in catabolism and protein transport. We validated candidate genes by pyrosequencing and found that moderate biases in allelic expression were less stable than monoallelically expressed genes that showed consistent bias upon experimental validations at the population level and in subcloning experiments. Our approach provides the basis for the analysis of aASE in P. tricornutum and could be routinely implemented to test for variations in allele expression under different environmental conditions.

Cannabinoid Receptor Interacting Protein 1a (CRIP1a) in Health and Disease

Endocannabinoid signaling depends upon the CB₁ and CB₂ cannabinoid receptors, their endogenous ligands anandamide and 2-arachidonoylglycerol, and intracellular proteins that mediate responses via the C-terminal and other intracellular receptor domains. The CB₁ receptor regulates and is regulated by associated G proteins predominantly of the Gi/o subtypes, β-arrestins 1 and 2, and the cannabinoid receptor-interacting protein 1a (CRIP1a). Evidence for a physiological role for CRIP1a is emerging as data regarding the cellular localization and function of CRIP1a are generated. Here we summarize the neuronal distribution and role of CRIP1a in endocannabinoid signaling, as well as discuss investigations linking CRIP1a to development, vision and hearing sensory systems, hippocampus and seizure regulation, and psychiatric disorders including schizophrenia. We also examine the genetic and epigenetic association of CRIP1a within a variety of cancer subtypes. This review provides evidence upon which to base future investigations on the function of CRIP1a in health and disease.

CTCF-Mediated Genome Architecture Regulates the Dosage of Mitotically Stable Mono-allelic Expression of Autosomal Genes

The mechanisms that guide the clonally stable random mono-allelic expression of autosomal genes remain enigmatic. We show that (1) mono-allelically expressed (MAE) genes are assorted and insulated from bi-allelically expressed (BAE) genes through CTCF-mediated chromatin loops; (2) the cell-type-specific dynamics of mono-allelic expression coincides with the gain and loss of chromatin insulator sites; (3) dosage of MAE genes is more sensitive to the loss of chromatin insulation than that of BAE genes; and (4) inactive alleles of MAE genes are significantly more insulated than active alleles and are de-repressed upon CTCF depletion. This alludes to a topology wherein the inactive alleles of MAE genes are insulated from the spatial interference of transcriptional states from the neighboring bi-allelic domains via CTCF-mediated loops. We propose that CTCF functions as a typical insulator on inactive alleles, but facilitates transcription through enhancer-linking on active allele of MAE genes, indicating widespread allele-specific regulatory roles of CTCF.

Functional effects of protein variants

Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.

Quantifying Genomic Imprinting at Tissue and Cell Resolution in the Brain

Imprinted genes are a group of ~150 genes that are preferentially expressed from one parental allele owing to epigenetic marks asymmetrically distributed on inherited maternal and paternal chromosomes. Altered imprinted gene expression causes human brain disorders such as Prader-Willi and Angelman syndromes and additional rare brain diseases. Research data principally obtained from the mouse model revealed how imprinted genes act in the normal and pathological brain. However, a better understanding of imprinted gene functions calls for building detailed maps of their parent-of-origin-dependent expression and of associated epigenetic signatures. Here we review current methods for quantifying genomic imprinting at tissue and cell resolutions, with a special emphasis on methods to detect parent-of-origin dependent expression and their applications to the brain. We first focus on bulk RNA-sequencing, the main method to detect parent-of-origin-dependent expression transcriptome-wide. We discuss the benefits and caveats of bulk RNA-sequencing and provide a guideline to use it on F1 hybrid mice. We then review methods for detecting parent-of-origin-dependent expression at cell resolution, including single-cell RNA-seq, genetic reporters, and molecular probes. Finally, we provide an overview of single-cell epigenomics technologies that profile additional features of genomic imprinting, including DNA methylation, histone modifications and chromatin conformation and their combination into sc-multimodal omics approaches, which are expected to yield important insights into genomic imprinting in individual brain cells.

The function of a heterozygous p53 mutation in a Li-Fraumeni syndrome patient

p53 is one of the most extensively studied proteins in cancer research. Mutations in p53 generally abolish normal p53 function, and some mutants can gain new oncogenic functions. However, the mechanisms underlying p53 mutation-driven cancer remains to be elucidated. Our study investigated the function of a heterozygous p53 mutation (p.Asn268Glufs*4) in a Li-Fraumeni syndrome (LFS) patient. We used episomal technology to perform somatic reprogramming, and used molecular and cell biology methods to determine the p53 mutation levels in patient-originated induced pluripotent stem (iPS) cells at the RNA and protein levels. We found that p53 protein expression was not increased in this patient’s somatic cells compared with those of a healthy control. p53 mutation facilitates the proliferation of tumor cells by inhibiting apoptosis and promoting cell division. It can inhibit the efficiency of somatic reprogramming by inhibiting OCT4 expression during reprogramming stage. Moreover, not all p53 mutant iPS cell lines have mutant p53 RNA sequences. A small percentage of mutant p53 mRNA is present in the somatic cells from the patient and his mother. In summary, this p53 mutation can promote tumor cell proliferation, inhibit somatic reprogramming, and exhibit random p53 allelic expression of heterozygous mutations in the patient and iPS cells which may be one of the reasons why the people with p53 mutations develop cancer at random. This finding suggested that mutant p53 allelic expression should be added to the risk forecasting of cancer.

The influence of DNA methylation on monoallelic expression

Monoallelic gene expression occurs in diploid cells when only one of the two alleles of a gene is active. There are three main classes of genes that display monoallelic expression in mammalian genomes: (1) imprinted genes that are monoallelically expressed in a parent-of-origin dependent manner; (2) X-linked genes that undergo random X-chromosome inactivation in female cells; (3) random monoallelically expressed single and clustered genes located on autosomes. The heritability of monoallelic expression patterns during cell divisions implies that epigenetic mechanisms are involved in the cellular memory of these expression states. Among these, methylation of CpG sites on DNA is one of the best described modification to explain somatic inheritance. Here, we discuss the relevance of DNA methylation for the establishment and maintenance of monoallelic expression patterns among these three groups of genes, and how this is intrinsically linked to development and cellular states.

Reciprocal monoallelic expression of ASAR lncRNA genes controls replication timing of human chromosome 6

DNA replication occurs on mammalian chromosomes in a cell-type distinctive temporal order known as the replication timing program. We previously found that disruption of the noncanonical lncRNA genes ASAR6 and ASAR15 results in delayed replication timing and delayed mitotic chromosome condensation of human chromosomes 6 and 15, respectively. ASAR6 and ASAR15 display random monoallelic expression and display asynchronous replication between alleles that is coordinated with other random monoallelic genes on their respective chromosomes. Disruption of the expressed allele, but not the silent allele, of ASAR6 leads to delayed replication, activation of the previously silent alleles of linked monoallelic genes, and structural instability of human chromosome 6. In this report, we describe a second lncRNA gene (ASAR6-141) on human chromosome 6 that when disrupted results in delayed replication timing in cisASAR6-141 is subject to random monoallelic expression and asynchronous replication and is expressed from the opposite chromosome 6 homolog as ASAR6 ASAR6-141 RNA, like ASAR6 and ASAR15 RNAs, contains a high L1 content and remains associated with the chromosome territory where it is transcribed. Three classes of cis-acting elements control proper chromosome function in mammals: origins of replication, centromeres, and telomeres, which are responsible for replication, segregation, and stability of all chromosomes. Our work supports a fourth type of essential chromosomal element, the "Inactivation/Stability Center," which expresses ASAR lncRNAs responsible for proper replication timing, monoallelic expression, and structural stability of each chromosome.

Incomplete penetrance in primary immunodeficiency: a skeleton in the closet

Primary immunodeficiencies (PIDs) comprise a diverse group of over 400 genetic disorders that result in clinically apparent immune dysfunction. Although PIDs are classically considered as Mendelian disorders with complete penetrance, we now understand that absent or partial clinical disease is often noted in individuals harboring disease-causing genotypes. Despite the frequency of incomplete penetrance in PID, no conceptual framework exists to categorize and explain these occurrences. Here, by reviewing decades of reports on incomplete penetrance in PID we identify four recurrent themes of incomplete penetrance, namely genotype quality, (epi)genetic modification, environmental influence, and mosaicism. For each of these principles, we review what is known, underscore what remains unknown, and propose future experimental approaches to fill the gaps in our understanding. Although the content herein relates specifically to inborn errors of immunity, the concepts are generalizable across genetic diseases.

New subtypes of allele-specific epigenetic effects: implications for brain development, function and disease

Typically, it is assumed that the maternal and paternal alleles for most genes are equally expressed. Known exceptions include canonical imprinted genes, random X-chromosome inactivation, olfactory receptors and clustered protocadherins. Here, we highlight recent studies showing that allele-specific expression is frequent in the genome and involves subtypes of epigenetic allelic effects that differ in terms of heritability, clonality and stability over time. Different forms of epigenetic allele regulation could have different roles in brain development, function, and disease. An emerging area involves understanding allelic effects in a cell-type and developmental stage-specific manner and determining how these effects influence the impact of genetic variants and mutations on the brain. A deeper understanding of epigenetics at the allele and cellular level in the brain could help clarify the mechanisms underlying phenotypic variance.

Characterization of genetic-origin-dependent monoallelic expression in mouse embryonic stem cells

Monoallelic gene expression occurs in various mammalian cells and can be regulated genetically, epigenetically and/or stochastically. We identified 145 monoallelically expressed genes (MoEGs), including seven known imprinted genes, in mouse embryonic stem cells (ESCs) derived from reciprocal F1 hybrid blastocysts and cultured in 2i/LIF. As all MoEGs except for the imprinted genes were expressed in a genetic-origin-dependent manner, we focused on this class of MoEGs for mechanistic studies. We showed that a majority of the genetic-origin-dependent MoEGs identified in 2i/LIF ESCs remain monoallelically expressed in serum/LIF ESCs, but become more relaxed or even biallelically expressed upon differentiation. These MoEGs and their regulatory regions were highly enriched for single nucleotide polymorphisms. In addition, some MoEGs were associated with retrotransposon insertions/deletions, consistent with the fact that certain retrotransposons act as regulatory elements in pluripotent stem cells. Interestingly, most MoEGs showed allelic differences in enrichment of histone H3K27me and H3K4me marks, linking allelic epigenetic differences and monoallelic expression. In contrast, there was little or no allelic difference in CpG methylation or H3K9me. Taken together, our study highlights the impact of genetic variation including single nucleotide polymorphisms and retrotransposon insertions/deletions on monoallelic epigenetic marks and expression in ESCs.

Role of the Chromosome Architectural Factor SMCHD1 in X-Chromosome Inactivation, Gene Regulation, and Disease in Humans

Structural maintenance of chromosomes flexible hinge domain-containing 1 (SMCHD1) is an architectural factor critical for X-chromosome inactivation (XCI) and the repression of select autosomal gene clusters. In mice, homozygous nonsense mutations in Smchd1 cause female-specific embryonic lethality due to an XCI defect. However, although human mutations in SMCHD1 are associated with congenital arhinia and facioscapulohumeral muscular dystrophy type 2 (FSHD2), the diseases do not show a sex-specific bias, despite the essential nature of XCI in humans. To investigate whether there is a dosage imbalance for the sex chromosomes, we here analyze transcriptomic data from arhinia and FSHD2 patient blood and muscle cells. We find that X-linked dosage compensation is maintained in these patients. In mice, SMCHD1 controls not only protocadherin (Pcdh) gene clusters, but also Hox genes critical for craniofacial development. Ablating Smchd1 results in aberrant expression of these genes, coinciding with altered chromatin states and three-dimensional (3D) topological organization. In a subset of FSHD2 and arhinia patients, we also found dysregulation of clustered PCDH, but not HOX genes. Overall, our study demonstrates preservation of XCI in arhinia and FSHD2, and implicates SMCHD1 in the regulation of the 3D organization of select autosomal gene clusters.

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.

Epigenetic Control of Gene Expression in the Normal and Malignant Human Prostate: A Rapid Response Which Promotes Therapeutic Resistance

A successful prostate cancer must be capable of changing its phenotype in response to a variety of microenvironmental influences, such as adaptation to treatment or successful proliferation at a particular metastatic site. New cell phenotypes emerge by selection from the large, genotypically heterogeneous pool of candidate cells present within any tumor mass, including a distinct stem cell-like population. In such a multicellular model of human prostate cancer, flexible responses are primarily governed not only by de novo mutations but appear to be dominated by a combination of epigenetic controls, whose application results in treatment resistance and tumor relapse. Detailed studies of these individual cell populations have resulted in an epigenetic model for epithelial cell differentiation, which is also instructive in explaining the reported high and inevitable relapse rates of human prostate cancers to a multitude of treatment types.

Functional evaluation of transposable elements as enhancers in mouse embryonic and trophoblast stem cells

Transposable elements (TEs) are thought to have helped establish gene regulatory networks. Both the embryonic and extraembryonic lineages of the early mouse embryo have seemingly co-opted TEs as enhancers, but there is little evidence that they play significant roles in gene regulation. Here we tested a set of long terminal repeat TE families for roles as enhancers in mouse embryonic and trophoblast stem cells. Epigenomic and transcriptomic data suggested that a large number of TEs helped to establish tissue-specific gene expression programmes. Genetic editing of individual TEs confirmed a subset of these regulatory relationships. However, a wider survey via CRISPR interference of RLTR13D6 elements in embryonic stem cells revealed that only a minority play significant roles in gene regulation. Our results suggest that a subset of TEs are important for gene regulation in early mouse development, and highlight the importance of functional experiments when evaluating gene regulatory roles of TEs.

Monoallelic expression in melanoma

BACKGROUND

Monoallelic expression (MAE) is a frequent genomic phenomenon in normal tissues, however its role in cancer is yet to be fully understood. MAE is defined as the expression of a gene that is restricted to one allele in the presence of a diploid heterozygous genome. Constitutive MAE occurs for imprinted genes, odorant receptors and random X inactivation. Several studies in normal tissues have showed MAE in approximately 5-20% of the cases. However, little information exists on the MAE rate in cancer. In this study we assessed the presence and rate of MAE in melanoma. The genetic basis of melanoma has been studied in depth over the past decades, leading to the identification of mutations/genetic alterations responsible for melanoma development.

METHODS

To examine the role of MAE in melanoma we used 15 melanoma cell lines and compared their RNA-seq data with genotyping data obtained by the parental TIL (tumor infiltrating lymphocytes). Genotyping was performed using the Illumina HumanOmni1 beadchip. The RNA-seq library preparation and sequencing was performed using the Illumina TruSeq Stranded Total RNA Human Kit and subsequently sequenced using a HiSeq 2500 according to manufacturer's guidelines. By comparing genotyping data with RNA-seq data, we identified SNPs in which DNA genotypes were heterozygous and corresponding RNA genotypes were homozygous. All homozygous DNA genotypes were removed prior to the analysis. To confirm the validity to detect MAE, we examined heterozygous DNA genotypes from X chromosome of female samples as well as for imprinted and olfactory receptor genes and confirmed MAE.

RESULTS

MAE was detected in all 15 cell lines although to a different rate. When looking at the B-allele frequencies we found a preferential pattern of complete monoallelic expression rather then differential monoallelic expression across the 15 melanoma cell lines. As some samples showed high differences in the homozygous and heterozygous call rate, we looked at the single chromosomes and showed that MAE may be explained by underlying large copy number imbalances in some instances. Interestingly these regions included genes known to play a role in melanoma initiation and progression. Nevertheless, some chromosome regions showed MAE without CN imbalances suggesting that additional mechanisms (including epigenetic silencing) may explain MAE in melanoma.

CONCLUSION

The biological implications of MAE are yet to be realized. Nevertheless, our findings suggest that MAE is a common phenomenon in melanoma cell lines. Further analyses are currently being undertaken to evaluate whether MAE is gene/pathway specific and to understand whether MAE can be employed by cancers to achieve a more aggressive phenotype.

Allele-specific RNA imaging shows that allelic imbalances can arise in tissues through transcriptional bursting

Extensive cell-to-cell variation exists even among putatively identical cells, and there is great interest in understanding how the properties of transcription relate to this heterogeneity. Differential expression from the two gene copies in diploid cells could potentially contribute, yet our ability to measure from which gene copy individual RNAs originated remains limited, particularly in the context of tissues. Here, we demonstrate quantitative, single molecule allele-specific RNA FISH adapted for use on tissue sections, allowing us to determine the chromosome of origin of individual RNA molecules in formaldehyde-fixed tissues. We used this method to visualize the allele-specific expression of Xist and multiple autosomal genes in mouse kidney. By combining these data with mathematical modeling, we evaluated models for allele-specific heterogeneity, in particular demonstrating that apparent expression from only one of the alleles in single cells can arise as a consequence of low-level mRNA abundance and transcriptional bursting.

In mammals, most cells of the body contain two genetic datasets: one from the mother and one from the father, and—in theory—these two sets of information could contribute equally to produce the molecules in a given cell. In practice, however, this is not always the case, which can have dramatic implications for many traits, including visible features (such as fur color) and even disease outcomes. However, it remains difficult to measure the parental origin of individual molecules in a given cell and thus to assess what processes lead to an imbalance of the maternal and paternal contribution. We adapted a microscopy technique—called allele-specific single molecule RNA FISH—that uses a combination of fluorescent tags to specifically label one type of molecule, RNA, depending on its origin, for use in mouse kidney sections. Focusing on RNAs that were previously reported to show imbalance, we performed measurements and combined these with mathematical modeling to quantify imbalance in tissues and explain how these can arise. We found that we could recapitulate the observed imbalances using models that only take into account the random processes that produce RNA, without the need to invoke special regulatory mechanisms to create unequal contributions.

Characterization of human mosaic Rett syndrome brain tissue by single-nucleus RNA sequencing

In females with X-linked genetic disorders, wild-type and mutant cells coexist within brain tissue because of X-chromosome inactivation, posing challenges for interpreting the effects of X-linked mutant alleles on gene expression. We present a single-nucleus RNA sequencing approach that resolves mosaicism by using single-nucleotide polymorphisms in genes expressed in cis with the X-linked mutation to determine which nuclei express the mutant allele even when the mutant gene is not detected. This approach enables gene expression comparisons between mutant and wild-type cells within the same individual, eliminating variability introduced by comparisons to controls with different genetic backgrounds. We apply this approach to mosaic female mouse models and humans with Rett syndrome, an X-linked neurodevelopmental disorder caused by mutations in the gene encoding the methyl-DNA-binding protein MECP2, and observe that cell-type-specific DNA methylation predicts the degree of gene upregulation in MECP2-mutant neurons. This approach can be broadly applied to study gene expression in mosaic X-linked disorders.

Parallels between Mammalian Mechanisms of Monoallelic Gene Expression

Different types of monoallelic gene expression are present in mammals, some of which are highly flexible, whereas others are more rigid. These include allelic exclusion at antigen receptor loci, the expression of olfactory receptor genes, genomic imprinting, X-chromosome inactivation, and random monoallelic expression (MAE). Although these processes play diverse biological roles, and arose through different selective pressures, the underlying epigenetic mechanisms show striking resemblances. Regulatory transcriptional events are important in all systems, particularly in the specification of MAE. Combined with comparative studies between species, this suggests that the different MAE systems found in mammals may have evolved from analogous ancestral processes.