Derivation of pre-X inactivation human embryonic stem cells under physiological oxygen concentrations.

The presence of two active X chromosomes (XaXa) is a hallmark of the ground state of pluripotency specific to murine embryonic stem cells (ESCs). Human ESCs (hESCs) invariably exhibit signs of X chromosome inactivation (XCI) and are considered developmentally more advanced than their murine counterparts. We describe the establishment of XaXa hESCs derived under physiological oxygen concentrations. Using these cell lines, we demonstrate that (1) differentiation of hESCs induces random XCI in a manner similar to murine ESCs, (2) chronic exposure to atmospheric oxygen is sufficient to induce irreversible XCI with minor changes of the transcriptome, (3) the Xa exhibits heavy methylation of the XIST promoter region, and (4) XCI is associated with demethylation and transcriptional activation of XIST along with H3K27-me3 deposition across the Xi. These findings indicate that the human blastocyst contains pre-X-inactivation cells and that this state is preserved in vitro through culture under physiological oxygen.

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).

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.

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.

Myc Dysregulation in Activated Macrophages Initiates Iron-Mediated Lipid Peroxidation that Fuels Type I Interferon and Compromises TB Resistance.

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop clinical, mostly pulmonary, TB. To dissect mechanisms of susceptibility in immunocompetent individuals, we developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: development of necrotic lung lesions after infection with virulent Mtb. In this study, we explored the connectivity of the sst1-regulated pathways during prolonged macrophage activation with TNF. We determined that the aberrant response of the sst1-susceptible macrophages to TNF was primarily driven by conflicting Myc and antioxidant response pathways that resulted in a coordinated failure to properly sequester intracellular iron and activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNß superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. Our findings suggest a novel link between metabolic dysregulation in macrophages and susceptibility to TB, offering insights into potential therapeutic targets aimed at modulating macrophage function and improving TB control.

Coordination of the random asynchronous replication of autosomal loci.

Random monoallelic expression and asynchronous replication define an unusual class of autosomal mammalian genes. We show that every cell has randomly chosen either the maternal or paternal copy of each given autosome pair, such that alleles of these genes scattered across the chosen chromosome replicate earlier than the alleles on the homologous chromosome. Thus, chromosome-pair non-equivalence, rather than being limited to X-chromosome inactivation, is a fundamental property of mouse chromosomes.

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.

Cell state transition analysis identifies interventions that improve control of M. tuberculosis infection by susceptible macrophages.

Understanding cell state transitions and purposefully controlling them to improve therapies is a longstanding challenge in biological research and medicine. Here, we identify a transcriptional signature that distinguishes activated macrophages from TB-susceptible and TB-resistant mice. We then apply the cSTAR (cell State Transition Assessment and Regulation) approach to data from screening-by-RNA sequencing to identify chemical perturbations that shift the. transcriptional state of the TB-susceptible macrophages towards that of TB-resistant cells. Finally, we demonstrate that the compounds identified with this approach enhance resistance of the TB-susceptible mouse macrophages to virulent M. tuberculosis .

Cell state transition analysis identifies interventions that improve control of Mycobacterium tuberculosis infection by susceptible macrophages.

Understanding cell state transitions and purposefully controlling them to improve therapies is a longstanding challenge in biological research and medicine. Here, we identify a transcriptional signature that distinguishes activated macrophages from the tuberculosis (TB) susceptible and resistant mice. We then apply the cSTAR (cell state transition assessment and regulation) approach to data from screening-by-RNA sequencing to identify chemical perturbations that shift the transcriptional state of tumor necrosis factor (TNF)-activated TB-susceptible macrophages toward that of TB-resistant cells, i.e., prevents their aberrant activation without suppressing beneficial TNF responses. Last, we demonstrate that the compounds identified with this approach enhance the resistance of the TB-susceptible mouse macrophages to virulent Mycobacterium tuberculosis.

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.

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.

Diversity of developing peripheral glia revealed by single-cell RNA sequencing.

The peripheral nervous system responds to a wide variety of sensory stimuli, a process that requires great neuronal diversity. These diverse neurons are closely associated with glial cells originating from the neural crest. However, the molecular nature and diversity among peripheral glia are not understood. Here, we used single-cell RNA sequencing to profile developing and mature glia from somatosensory dorsal root ganglia and auditory spiral ganglia. We found that glial precursors (GPs) in these two systems differ in their transcriptional profiles. Despite their unique features, somatosensory and auditory GPs undergo convergent differentiation to generate molecularly uniform myelinating and non-myelinating Schwann cells. By contrast, somatosensory and auditory satellite glial cells retain system-specific features. Lastly, we identified a glial signature gene set, providing new insights into commonalities among glia across the nervous system. This survey of gene expression in peripheral glia constitutes a resource for understanding functions of glia across different sensory modalities.

KDM5 Histone Demethylase Activity Links Cellular Transcriptomic Heterogeneity to Therapeutic Resistance.

Members of the KDM5 histone H3 lysine 4 demethylase family are associated with therapeutic resistance, including endocrine resistance in breast cancer, but the underlying mechanism is poorly defined. Here we show that genetic deletion of KDM5A/B or inhibition of KDM5 activity increases sensitivity to anti-estrogens by modulating estrogen receptor (ER) signaling and by decreasing cellular transcriptomic heterogeneity. Higher KDM5B expression levels are associated with higher transcriptomic heterogeneity and poor prognosis in ER+ breast tumors. Single-cell RNA sequencing, cellular barcoding, and mathematical modeling demonstrate that endocrine resistance is due to selection for pre-existing genetically distinct cells, while KDM5 inhibitor resistance is acquired. Our findings highlight the importance of cellular phenotypic heterogeneity in therapeutic resistance and identify KDM5A/B as key regulators of this process.

Tumor-suppressor genes that escape from X-inactivation contribute to cancer sex bias.

There is a striking and unexplained male predominance across many cancer types. A subset of X-chromosome genes can escape X-inactivation, which would protect females from complete functional loss by a single mutation. To identify putative 'escape from X-inactivation tumor-suppressor' (EXITS) genes, we examined somatic alterations from >4,100 cancers across 21 tumor types for sex bias. Six of 783 non-pseudoautosomal region (PAR) X-chromosome genes (ATRX, CNKSR2, DDX3X, KDM5C, KDM6A, and MAGEC3) harbored loss-of-function mutations more frequently in males (based on a false discovery rate < 0.1), in comparison to zero of 18,055 autosomal and PAR genes (Fisher's exact P < 0.0001). Male-biased mutations in genes that escape X-inactivation were observed in combined analysis across many cancers and in several individual tumor types, suggesting a generalized phenomenon. We conclude that biallelic expression of EXITS genes in females explains a portion of the reduced cancer incidence in females as compared to males across a variety of tumor types.

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.

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.

Olfactory-enriched transcripts are cell-specific markers in the lobster olfactory organ.

Genes expressed specifically in a tissue are often involved in the defining functions of that tissue. We used representational difference analysis of cDNA to amplify 20 cDNA fragments representing transcripts that were more abundant in the lobster olfactory organ than in brain, eye/eyestalk, dactyl, pereiopod, or second antenna. We then independently confirmed that the transcripts represented by these clones were enriched in the olfactory organ. The 20 cDNA fragments represent between 6 and 15 different genes. Six of the cDNAs contained sequences highly similar to known gene families. We performed in situ hybridization with these six and found that all were expressed in subsets of cells associated with the aesthetasc sensilla in the olfactory organ. Clones OET-07, an ionotropic receptor, and OET-10, an alpha tubulin, were specific to the olfactory receptor neurons. OET-02, a monooxygenase, was expressed only in the outer auxiliary cells. OET-03, a serine protease, was specific to the collar cells. OET-11, an alpha(2) macroglobulin, was expressed by the receptor neurons and the collar cells. OET-17, a calcyphosine, was expressed in the receptor neurons, inner auxiliary cells, and collar cells. The identities and expression patterns of these six transcripts predict involvement in both known and novel properties of the lobster olfactory organ.

Autosomal monoallelic expression: genetics of epigenetic diversity?

In mammals, relative expression of the two parental alleles of many genes is controlled by one of three major epigenetic phenomena: X chromosome inactivation, imprinting, and mitotically stable autosomal monoallelic expression (MAE). MAE affects a large fraction of human autosomal genes and introduces enormous epigenetic heterogeneity in otherwise similar cell populations. Despite its prevalence, many functional and mechanistic aspects of MAE biology remain unknown. Several lines of evidence imply that MAE establishment and maintenance are controlled by a variety of genetic elements. Based on known genomic features regulating X-inactivation and imprinting, we outline likely features of MAE-controlling elements. We also assess implications of MAE for genotype-phenotype relationship, with a focus on haploinsufficiency.

Chromatin signature of widespread monoallelic expression.

In mammals, numerous autosomal genes are subject to mitotically stable monoallelic expression (MAE), including genes that play critical roles in a variety of human diseases. Due to challenges posed by the clonal nature of MAE, very little is known about its regulation; in particular, no molecular features have been specifically linked to MAE. In this study, we report an approach that distinguishes MAE genes in human cells with great accuracy: a chromatin signature consisting of chromatin marks associated with active transcription (H3K36me3) and silencing (H3K27me3) simultaneously occurring in the gene body. The MAE signature is present in ∼20% of ubiquitously expressed genes and over 30% of tissue-specific genes across cell types. Notably, it is enriched among key developmental genes that have bivalent chromatin structure in pluripotent cells. Our results open a new approach to the study of MAE that is independent of polymorphisms, and suggest that MAE is linked to cell differentiation. DOI: http://dx.doi.org/10.7554/eLife.01256.001.